The content of the invention
The present invention proposes can be according to preset time value, the optimal movement of the foetus meter digital of discontinuity detection in a kind of movement of the foetus counting process
The detection method put, in more movement of the foetus sensor data acquisition, the optimal sensing station of selected section is as next
The movement of the foetus enumeration data collection source of period, to reduce movement of the foetus number system workload and improve count accuracy.
The technical scheme is that:
A kind of detection method of the optimal movement of the foetus counting position based on acceleration transducer, it is characterized in that it uses 8 three
Axle acceleration sensor carries out movement of the foetus detection to the belly of pregnant woman, the detection method of described optimal movement of the foetus counting position include with
Lower step:
The first step:8 3-axis acceleration sensors are close to the belly of pregnant woman, setting movement of the foetus detection in two horizontally-arranged arrangements
Time interval T, starting loop sensor detect to movement of the foetus, after time T, obtain m detection data, and record each time
Detect time point corresponding to data;
Second step:For the 1st time to the m times detection data, 8 3-axis acceleration sensors are chosen respectively in X-axis, Y-axis
4 maximum position P1-P4 with amplitude in each axle of Z axis;
3rd step:For any detection data of foregoing any 3-axis acceleration sensing, 4 in X-axis, Y-axis and Z axis
Descending arranges by size for position, and it is as follows to carry out computing:
It is P1-P4 for 4 maximum positions of amplitude in each axle of X, Y and Z, peak swing corresponding to each position is first
Amplitude is a1-a4;The fractional threshold parameter n1 of second amplitude and the first amplitude, the 3rd amplitude and the second amplitude ratio threshold are set
The fractional threshold parameter n3, n1 of value parameter n2, the 4th amplitude and the 3rd amplitude>n2>N3, S1, S2 and S3 are represented to assign respectively and passed
The score value of sensor, S1>S2>S3;
(1), for any detection data of any axle, it is corresponding that amplitude maximum a1, position P1 corresponding to imparting are selected
S1 point;
(2), for position P2 corresponding to a2, the ratio η of the second amplitude and the first amplitude is calculated1,If η1It is full
Foot:
η1≥n1, then s1 assign P2;
η1<n1Then s2 assigns P2
(3), for position P3 corresponding to a3, the ratio η of the 3rd amplitude and the second amplitude is calculated2,If η2It is full
Foot:
η2≥n2, then s2 assign P3;
η2<n2Then s3 assigns P3
(4), for position P4 corresponding to a4, the ratio η of the 4th amplitude and the 3rd amplitude is calculated3,If η3It is full
Foot:
η3≥n3, then s3 assign P4;
η3<n3Then s4 assigns P4
(5) score of remaining 4 position is calculated as 0
X, Y of 8 3-axis acceleration sensors and m detection data of Z axis are handled according to (1)-(5), obtained
The X-axis each time of each 3-axis acceleration sensor, Y-axis, 4 amplitude maximum positions corresponding to Z axis into m detection data
Score;
4th step:For m detection data of this detection, 8 3-axis acceleration sensors are calculated respectively in position X-axis, Y
The score of axle and each axle of Z axis and, calculate 8 3-axis acceleration sensors score summation, four higher positions of score are
4 optimal movement of the foetus counting positions that movement of the foetus counts in next period T.
In 8 3-axis acceleration sensors of the present invention, lower row is followed successively by position 1, position 2, position 3 and position from left to right
4 are put, upper row is followed successively by position 5, position 6, position 7 and position 8 from left to right, and setting position 1 and position 3 are one group, that is, organize A,
Position 2 and position 4 are one group, that is, organize B, position 5 and position 7 are one group, that is, organize C;Position 6 and position 8 are one group, that is, organize D;
After 4th step, compare the score of A groups and B groups, C groups and D groups, when the position corresponding to higher two groups of score is next
Between movement of the foetus counts in section T 4 optimal movement of the foetus counting positions.
Z-direction is paid the utmost attention in 8 3-axis acceleration sensors of the present invention, during selection, higher priority, X are set
Axle and Y-axis set lower priority, i.e., when seeking 8 position Synthesis scores, setting Z axis coefficient is higher than the coefficient of X-axis and Y-axis,
X, the comprehensive score of 8 positions is calculated between 0 to 1 in the coefficient of Y and Z axis.
The present invention's sets Z axis coefficient as 0.5, and X-axis and Y-axis coefficient are 0.25, so as to which the synthesis of 8 positions be calculated
Score.
8 3-axis acceleration sensors of the present invention are embedded in pregnant woman's intelligence clothing or abdominal belt.
In the first step of the present invention:Data are detected for m times, detect the interval time of data before and after calculating twice, if
Interval time is less than or equal to Δ t, then after a testing number evidence and it is preceding once belong to a movement of the foetus data, will before a testing number
According to time of the time as this movement of the foetus data, detect amplitude is big in data amplitude twice as this movement of the foetus number using front and rear
According to amplitude;Aforementioned processing is carried out for m detection data, obtains M movement of the foetus enumeration data, for subsequent treatment, wherein,
3min < Δ t < 6min.
The s1=1, s2=0.5, s3=0.2, s4=0 of the present invention;N1=0.8, n2=0.5, n3=0.4.
In the first step of the present invention, four groups of initial optimum detection positions are set, are position 1, position 3, the and of position 5 respectively
Position 7.
In the first step of the present invention, time interval T is 2-8 hours.
Beneficial effects of the present invention:
The present invention proposes a kind of detection method of optimal movement of the foetus counting position in movement of the foetus counting process, is passed in more movement of the foetus
In sensor data acquisition, collection source of the maximally effective sensing station of selected section as movement of the foetus enumeration data, catch to greatest extent
Catch Fetal Movement Signal, optimization Fetal Movement Signal gathered data source, the workload of operation of mitigation movement of the foetus number system software and data processing
Difficulty, while decrease the data space of system, real-time and accurately show that movement of the foetus calculates, and necessarily reduce system it is soft,
Hardware cost of manufacture;The individual difference that the present invention can also change according to pregnant woman's movement of the foetus, sets different test interval section T to enter
Row selection, the frequency of optimum position adjustment is increased or reduces, further to improve the accuracy of optimum position detection.
Embodiment
The present invention is further illustrated with reference to the accompanying drawings and examples.
As shown in figure 1, in order to obtain more comprehensively accurate Fetal Movement Signal, present invention introduces 83 axle acceleration sensors,
Be embedded in pregnant woman's intelligence clothing or abdominal belt, be close to the belly of pregnant woman, respectively by this 8 positions be defined as position 1,2 ... 8.Often
Start optimum position detection program, the test position counted as movement of the foetus in next T time after interval time T.
During each optimum position detection, sensor can detect that movement of the foetus several times changes, and be calculated as m times,.M movement of the foetus changes
Data are detected, detect the interval time of data before and after calculating twice, if interval time is less than or equal to Δ t, rear one-time detection
Data and it is preceding once belong to a movement of the foetus data, the time using the time of a preceding testing number evidence as this movement of the foetus data,
Amplitude using the big amplitude of amplitude in the front and rear data of detection twice as this movement of the foetus data;Before being carried out for m detection data
Processing is stated, obtains M movement of the foetus enumeration data, for subsequent treatment, wherein, 3min < Δ t < 6min.
Such as the present embodiment, start optimum position detection program after counting T hours away from last movement of the foetus.Because experiment condition limits
Make, in embodiment, only collect the counting of the first time movement of the foetus behind pregnant woman's T time interval, it is assumed here that T interval times are shorter and locate
In the fetus peace and quiet period, therefore movement of the foetus is counted only there occurs once in T time, and next T is used as using this movement of the foetus delta data
The optimum position detection data basis of period.
Comprise the following steps that:
The first step:Gather the testing number of 83 axle acceleration sensors during the M times movement of the foetus after preset time interval T
According to altogether there occurs movement of the foetus counting, and this movement of the foetus sensor for countering detects that movement of the foetus changes m times.Using digital coding+X,
Y, Z axis represents the data that some axle of some position receives.Such as 3x represents the data that the X-axis of position 3 receives, 4Y tables
Show the data that the Y-axis of position 4 receives.Obtained data are as shown in Table 2-4.
Second step:This movement of the foetus is that sensor obtains 14 detection data, i.e. m1=14 there occurs 14 movement of the foetus changes.
Time point corresponding to the change of movement of the foetus each time can also be recorded.Choose respectively in the 1st time to the m times movement of the foetus each time change
4 maximum positions of amplitude in X-axis, Y-axis, Z axis corresponding to 83 axle sensor positions.As shown in Table 2-4, the 1st movement of the foetus
Change, 4 maximum positions of amplitude are 1,2,3,6 in X-axis;4 maximum positions of amplitude are 2,3,5,7 in Y-axis;On Z axis
4 maximum positions of amplitude are 2,3,4,5.
3rd step:4 maximum positions of X-axis, Y-axis, Z axis amplitude will be pressed respectively in continuous movement of the foetus change procedure each time
Size descending arranges, and carries out score computing.Computing is as follows:
This 4 positions are arranged to position p1, p2, p3, p4, and peak swing corresponding to each position is a1, a2, a3, a4.
Threshold parameter n1, n2, n3 are set, the score of this 4 positions is obtained according to these threshold parameters.
1) the maximum a1 of amplitude is selected, assigns corresponding to the p1 of position s1 points
2) position p2 corresponding to a2 is investigated, if η (η=a2/a1) meets
η>Then s1 assigns p2 to=n1
η<Then s2 assigns p2 to n1
3) position p3 corresponding to a3 is investigated, if η (η=a3/a2) meets
η>Then s2 assigns p3 to=n2
η<Then s3 assigns p3 to n2
4) position p4 corresponding to a4 is investigated, if η (η=a4/a3) meets
η>Then s3 assigns p4 to=n3
η<Then s4 assigns p4 to n3
5) score of remaining 4 position is calculated as 0
Thus movement of the foetus X-axis, Y-axis, 4 amplitude maximum positions corresponding to Z axis each time can be calculated in m detection data
Score.
Such as, arrange parameter n1=0.8, n2=0.5, n3=0.4;S1=1, s2=0.5, s3=0.2, s4=0.6th time
4 maximum positions of the amplitude of Z axis are 2,3,7,8 respectively in continuous movement of the foetus change, and its corresponding amplitude is 0.025g, 0.08g,
0.01g, 0.03g.3,8,2,7 are adjusted to according to descending arrangement, the position for being set to answer is p1, p2, p3, and its corresponding amplitude of p4 is
A1, a2, a3, a4.
1) p1 responses are most strong, therefore p1=s1=1
2) p2 is investigated, because η=a2/a1=0.03/0.08=0.375<N1 (0.8), so p2=s2=0.5
3) p3 is investigated, because η=a3/a2=0.025/0.03=0.833>N2 (0.5), so p3=s2=0.5
4) p4 is investigated, because η=a4/a3=0.01/0.025=0.4>=n3 (0.4), so p4=s3=0.2
Therefore the grade form 1 of eight positions of Z-direction in the 6th continuous movement of the foetus change can be obtained:
Position |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
Scores |
0 |
0.5 |
1 |
0 |
0 |
0 |
0.2 |
0.5 |
The scoring of eight positions of Z-direction in the 6th continuous movement of the foetus change of table 1
The score for 14 continuous movement of the foetus change X, Y, Z axis data that this movement of the foetus collects is handled successively, obtains 8 three axles
Acceleration transducer position the score in tri- directions of Y, Z, is included in table 2, table 3, table 4 in X.
The left side of table represents that the sequence number of continuous movement of the foetus change frequency, such as FetMov 1 represent this movement of the foetus sensor for countering
The movement of the foetus delta data detected, Time are arranged apart from the initial time time at the time of representing to occur movement of the foetus, and e.g., the 1st movement of the foetus becomes
It is 5 seconds to change apart from the time of program triggering, and the 14th movement of the foetus change is 261 seconds apart from the time of program triggering.Selected
Times represents to choose as the number of one of 4 most strong amplitude in m detection data.Scores represents each position
Score accumulates.
Continuous movement of the foetus each time changes X-axis peak swing position score during this time movement of the foetus of table 2
Continuous movement of the foetus each time changes Y-axis peak swing position score during this time movement of the foetus of table 3
Continuous movement of the foetus each time changes Z axis peak swing position score during this time movement of the foetus of table 4
4th step:Add up the m continuous movement of the foetus change of this movement of the foetus process, 8 3-axis acceleration sensor position X-axis, Y
The score summation of axle, Z axis.As shown in table 4.
Position |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
X |
7.6 |
5.8 |
5 |
5.3 |
1.7 |
1.4 |
2 |
1 |
Y |
4.1 |
7 |
8.4 |
0 |
4.3 |
1.6 |
3.1 |
1.7 |
Z |
4.4 |
10.2 |
10.2 |
2.7 |
1.6 |
0.2 |
0.8 |
1 |
Score of each position in X, Y, Z axis during this time movement of the foetus of table 4
5th step:Accuracy in view of test and comprehensive, when being selected in position, setting position 1 and position 3 are one
Group, that is, organize A;Position 2 and position 4 are one group, that is, organize B;Select one group of seizure fetus dynamic below her abdominal in two groups
Make;Similarly, position 5 and position 7 are one group, that is, organize C;Position 6 and position 8 are one group, that is, organize D.One group of seizure fetus is selected to exist
Action above her abdominal.Group sets higher priority, X-axis and Y-axis with paying the utmost attention to Z-direction during the selection principle organized
Lower priority is set.Therefore when seeking 8 position Synthesis scores, it may be assumed that Z axis coefficient is 0.5, and X-axis and Y-axis coefficient are
0.25.So as to which the comprehensive score of 8 positions be calculated.
Such as, the final score of position 1 is:
P1=1Z*0.5+1X*0.25+1Y*0.25=4.4*0.5+4.1*0.25+7.6*0.25=5.12 5;
Similarly, finally giving for other 7 positions is obtained, as shown in table 5.
Position |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
Scores |
5.125 |
8.3 |
8.45 |
2.675 |
5.125 |
0.85 |
1.675 |
1.175 |
The final score of each position during this time movement of the foetus of table 5
6th step:Compare the score of A groups and B groups, C groups and D groups, corresponding position is as next in higher two groups of score
4 optimum detection positions that individual T time section movement of the foetus counts.
Such as:A group scores:5.125+8.45=13.575;Similarly obtain B, C, D group score be respectively 10.975,6.8,
2.025.The response intensity of A groups and C groups is more than its corresponding group.
Therefore 4 optimum detection positions that next T time section movement of the foetus counts are 1,3,5,7.
Part that the present invention does not relate to is same as the prior art or can be realized using prior art.