CN105030231B - A kind of fetal heart monitoring accelerates the extracting method of characteristic parameter - Google Patents
A kind of fetal heart monitoring accelerates the extracting method of characteristic parameter Download PDFInfo
- Publication number
- CN105030231B CN105030231B CN201510420283.XA CN201510420283A CN105030231B CN 105030231 B CN105030231 B CN 105030231B CN 201510420283 A CN201510420283 A CN 201510420283A CN 105030231 B CN105030231 B CN 105030231B
- Authority
- CN
- China
- Prior art keywords
- fetal heart
- heart rate
- signal
- acceleration
- mrow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
A kind of acceleration parameter acquisition methods based on Fetal Heart Rate signal in fetal heart monitoring, including signal input detection module, signal for faster identification extraction module, unit rate of acceleration computing module, and acceleration parameter output module.Wherein：Signal input detection module includes Signal Pretreatment and useful signal is extracted, signal for faster identification extraction module includes the calculating of baseline, the identification of signal for faster, the extraction of signal for faster and the pretreatment to signal for faster, unit rate of acceleration computing module includes unit of account rate of acceleration and the unit interval accelerates lattice number.This method can be reduced the main error brought by visual method at present, effectively improved the deciphering efficiency to fetal heart monitoring figure using characteristic parameter of the Fetal Heart Rate signal acquisition about acceleration.
Description
Technical field
The present invention relates to processing of biomedical signals field, the feature ginseng that Fetal Heart Rate accelerates in specifically a kind of fetal heart monitoring
Number extracting method.
Background technology
Electronic fatal monitoring as a kind of real-time Intrauterine Monitoring technology for obtaining fetal heart rate signal and parent uterine contraction signal,
At present in clinical obstetricses extensive use.Because fetus is located in parent abdominal uterine, fetal heart sound signal is highly susceptible to come from
The interference of the noise of parent, sensor and external environment.Current fetal heart monitoring system mostly applies principle of Doppler simultaneously
Fetal Instantaneous Heart Rate is calculated using auto-correlation algorithm.Fetal Heart Rate accelerates to be used as the important composition portion in fetal electronic monitoring figure
Point, it is a kind of zooming phenomenon of Fetal Heart Rate, time duration is shown as on figure and heart rate amplitude all reaches necessarily
The projection of threshold value.
The analysis accelerated at present to Fetal Heart Rate is mainly made of the method for range estimation to acceleration to be classified and counts, and Fetal Heart Rate is believed
Number acceleration contained by information extraction it is not abundant enough.Using Digital Signal Processing, numeral is subject to Fetal Heart Rate signal
Change, more objective parameter value can be obtained, help to obtain more comprehensively information.
The content of the invention
Existing technical scheme, which is mainly used, accelerates definition in the fetal heart rate signal collected and electronic fatal monitoring
Fetal Heart Rate is accelerated to do to classify and count.Monitoring to prenatal foetal is a comprehensive complicated process, is obtained in the method for range estimation
The Fetal Heart Rate information taken is not objective enough, and otherness is larger, and only relies on category division and the counting of heart rate acceleration, it is difficult to obtain
Obtain the more information of Fetal Heart Rate signal.
To solve the above problems, the present invention proposes a kind of acceleration Eigenvalue Extraction Method based on Fetal Heart Rate signal, carry
The acceleration characteristic value name of taking-up is unit rate of acceleration, and the parameter utilizes the prolonged fetal heart rate data of fetal heart monitoring and fetal rhythm
Rate accelerates definition to calculate acquisition, can be from fetal rhythm with reference to the classification that Fetal Heart Rate accelerates and counting with clearer and more definite physiological significance
More comprehensively information is obtained in rate signal.It is possible thereby to improve acquisition more piece of the prior art to information in Fetal Heart Rate signal
The situation in face.
To reach above-mentioned purpose, the technical solution used in the present invention is：A kind of acceleration parameter based on Fetal Heart Rate signal
Acquisition methods, it is characterised in that：Detection module is inputted including signal, signal for faster identification extraction module, unit rate of acceleration is calculated
Module, acceleration parameter output module.By pretreatment and the extraction of useful signal after the input of Fetal Heart Rate signal, know in signal for faster
Baseline value is calculated in other extraction module, recognizes and pre-processes signal for faster, is calculated afterwards in unit rate of acceleration computing module
Unit rate of acceleration and unit interval accelerate lattice number, finally in the acceleration parameter output module output result.
The signal input detection module receives the fetal heart monitoring time domain Fetal Heart Rate signal f (t) of actual measurement.
The signal inputs detection module, the time domain Fetal Heart Rate signal f (t) of input is pre-processed, i.e., using existing
Technology is classified as tire in the Fetal Heart Rate signal sequence f (k) of N sections of certain time length, each section of Fetal Heart Rate signal sequence f (k)
Heart rate value is 0 lossing signal f (k0), its length is l (f (k0The ratio between)), with this section of Fetal Heart Rate signal sequence length l (f (k))
For dropout ratePick out and meet the Fetal Heart Rate signal that dropout rate lost is less than specific threshold low
Sequence g (k), i.e.,Wherein g (k0) be Fetal Heart Rate signal sequence g (k) Fetal Heart Rate value be 0 loss
Signal.Effective Fetal Heart Rate signal time total length is
The signal for faster identification extraction module receives the Fetal Heart Rate signal for meeting and requiringAnd root
The definition accelerated according to different Fetal Heart Rates, to accelerating identification, is below specific recognition methods.
For continually enter satisfaction require Fetal Heart Rate signal bpm (k), if meet heart rate value bpm (k) >=b (k) into
It is vertical, then start counting up, as bpm (k+aj) ＜ b (k+aj) set up when then stop count (count not include ajIt is individual, therefore have 0,1,
2...(aj- 1) common ajIndividual numerical value), if now count value aj≥m1And meet bpm (kr)≥b(kr)+m2(wherein kr∈[k,k+
aj- 1], m1With m2For specific threshold), then retain numerical value ajAnd wj=k+aj-1；If count value aj＜ m1, then by numerical value aj
Zero setting.Then Fetal Heart Rate acceleration time total length expression formula is：
Wherein, bpm (k) is the Fetal Heart Rate value at k moment, and b (k) is the baseline value at k moment, m1With m2It is threshold constant, can
Depending on pregnant woman's characteristic change.ajFor+1 count value of jth, wjTo retain numerical value ajK moment when -1.Baseline refers to remove Fetal Heart Rate
Average Fetal Heart Rate of 10 minutes inner tube of a tyre heart rate volatility scopes in 5bpm after acceleration, deceleration and significant variation fragment.
Fetal Heart Rate accelerate gross area expression formula be：
The unit rate of acceleration computing module, is added using Fetal Heart Rate acceleration time total length L obtained above with Fetal Heart Rate
Fast gross area S is calculated, and draws unit rate of acceleration.
Unit rate of acceleration expression formula is：
Wherein PAR is unit rate of acceleration, and the Fetal Heart Rate being characterized as within the unit interval during fetal heart monitoring accelerates accelerates
Amplitude, unit rate of acceleration PAR can obtain lattice number shared by unit interval acceleration area in fetal heart monitoring figure with reference to special parameter.
Unit interval accelerate area lattice number expression formula be：
Wherein SΔtFor the acceleration area lattice number in unit time Δ t, LtFor Fetal Heart Rate effective time total length, Δ s is tire
The area of the single acceleration grid of heart monitoring, a length of specific duration of the grid, a width of specific heart rate value.
The acceleration parameter output module, the above-mentioned fetal heart monitoring acceleration parameter result calculated is exported.
Fetal heart monitoring accelerates the beneficial effect of characteristic parameter extraction to be：
This method is only needed to input the Fetal Heart Rate primary signal of fetal heart monitoring, you can multiple relevant tires are obtained under certain precision
The parameter of cadion-acceleration.Step is simple and convenient, improves the less acquisition Fetal Heart Rate that is difficult to of the index accelerated at present on Fetal Heart Rate and believes
The situation of number more information.The error brought by visual method can be reduced in actual application, is effectively improved on tire
The heart guards the deciphering efficiency of figure.
Brief description of the drawings
Fig. 1 is that fetal heart monitoring accelerates characteristic parameter extraction structured flowchart.
Fig. 2 is that fetal heart monitoring accelerates characteristic parameter extraction operational flowchart.
Fig. 3 is Fetal Heart Rate signal for faster and relevant parameter schematic diagram.
Embodiment
A kind of more typical embodiment of the present invention is described in detail below in conjunction with accompanying drawing.
A kind of typical application scenarios of the present invention are to obtain using fetal heart monitoring heart rate signal joins about the feature accelerated
Number.It is possible thereby to obtain more objective complete using the mature technology and abundant acceleration characteristic parameter of existing Fetal Heart Rate signal acquisition
The information on Fetal Heart Rate signal in face.
As shown in Fig. 2 after extraction process starts, selecting corresponding given threshold according to the pregnant week of measured first.With one
The position age was 31 one full year of life, and exemplified by pregnant week is the measured of 39 weeks, pregnant week is set in step T1.According to setting in step T2,
System is the given threshold that correspondence crowd is automatically selected out according to the definition accelerated to different pregnant weeks, now m1For 15s, m2For
15bpm。
System starts to receive actual measurement Fetal Heart Rate signal in step T3, and step T4 is divided the Fetal Heart Rate signal f (t) of input
Section, calculates the Loss Rate lost of every segment signal sequence f (k).
Signal Effective judgement is carried out according to Loss Rate in step T5, specific threshold low's is higher than to Loss Rate lost
Original signal sequence is given up.If signal sequence f (k) occur to give up, next segment signal sequence f (k+1) is counted again
Calculate with judging.Measured's overall signal Loss Rate is 4.25%.
Step T6 carries out accelerating part identification to the useful signal bpm (k) that requirement is met in step T5, calculates single hop and accelerates
Partial duration ajWith area Sj, according to formula：
By the processing to whole signals in step T7, calculate Fetal Heart Rate acceleration time total length L and accelerate with Fetal Heart Rate
Gross area S, unit rate of acceleration PAR, the unit interval accelerate area lattice number to take unit interval Δ t to be 20 minutes, and expression formula is：
It is computed obtaining：
L=108658 Lt=258728 S=929378 PAR=8.56 SΔt=14.37
Step T8 exports above-mentioned parameter according to specified format.
Claims (1)
1. a kind of fetal heart monitoring accelerates the extracting method of characteristic parameter, it is characterised in that：Detection module is inputted including signal, is accelerated
Signal identification extraction module, unit rate of acceleration computing module, acceleration parameter output module；By locating in advance after the input of Fetal Heart Rate signal
The extraction of reason and useful signal, calculates baseline value in signal for faster identification extraction module, recognizes and pre-process signal for faster,
Unit rate of acceleration is calculated in unit rate of acceleration computing module afterwards and the unit interval accelerates lattice number, finally in the acceleration parameter
Output module output result；
The signal input detection module receives the fetal heart monitoring time domain Fetal Heart Rate signal f (t) of actual measurement；
The signal inputs detection module, and the time domain Fetal Heart Rate signal f (t) of input is pre-processed, that is, utilizes prior art
It is classified as Fetal Heart Rate in the Fetal Heart Rate signal sequence f (k) of N sections of certain time length, each section of Fetal Heart Rate signal sequence f (k)
It is worth the lossing signal f (k for 00), its length is l (f (k0)), it is letter with the ratio between this section of Fetal Heart Rate signal sequence length l (f (k))
Number Loss RatePick out and meet the Fetal Heart Rate signal sequence g that dropout rate lost is less than specific threshold low
(k), i.e.,Wherein g (k0) be Fetal Heart Rate signal sequence g (k) Fetal Heart Rate value be 0 lossing signal；
Effective Fetal Heart Rate signal time total length is
The signal for faster identification extraction module receives the Fetal Heart Rate signal for meeting and requiringAnd according to not
Definition with Fetal Heart Rate acceleration, to accelerating identification, is below specific recognition methods；
The Fetal Heart Rate signal bpm (k) required for the satisfaction continually entered, sets up if meeting heart rate value bpm (k) >=b (k),
Then start counting up, as bpm (k+aj) ＜ b (k+aj) set up when then stop count (count not include ajIt is individual, therefore have 0,1,2...
(aj- 1) common ajIndividual numerical value), if now count value aj≥m1And meet bpm (kr)≥b(kr)+m2(wherein kr∈[k,k+aj-
1], m1With m2For specific threshold), then retain numerical value ajAnd wj=k+aj-1；If count value aj＜ m1, then by numerical value ajPut
Zero；Then Fetal Heart Rate acceleration time total length expression formula is：
<mrow>
<mi>L</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
</mrow>
Wherein, bpm (k) is the Fetal Heart Rate value at k moment, and b (k) is the baseline value at k moment, m1With m2It is threshold constant, it is visual pregnant
Woman's characteristic change；ajFor+1 count value of jth, wjTo retain numerical value ajK moment when -1；Baseline refer to remove Fetal Heart Rate accelerate,
Slow down and average Fetal Heart Rate of 10 minutes inner tube of a tyre heart rate volatility scopes in 5bpm after significant variation fragment；
Fetal Heart Rate accelerate gross area expression formula be：
<mrow>
<mi>S</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<msub>
<mi>w</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>w</mi>
<mi>j</mi>
</msub>
</munderover>
<mo>&lsqb;</mo>
<mi>b</mi>
<mi>p</mi>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
The unit rate of acceleration computing module, is accelerated total using Fetal Heart Rate acceleration time total length L obtained above and Fetal Heart Rate
Area S is calculated, and draws unit rate of acceleration；
Unit rate of acceleration expression formula is：
<mrow>
<mi>P</mi>
<mi>A</mi>
<mi>R</mi>
<mo>=</mo>
<mfrac>
<mi>S</mi>
<mi>L</mi>
</mfrac>
</mrow>
Wherein PAR is unit rate of acceleration, and the Fetal Heart Rate being characterized as within the unit interval during fetal heart monitoring accelerates accelerates amplitude,
Unit rate of acceleration PAR can obtain lattice number shared by unit interval acceleration area in fetal heart monitoring figure with reference to special parameter；
Unit interval accelerate area lattice number expression formula be：
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>P</mi>
<mi>A</mi>
<mi>R</mi>
<mo>&times;</mo>
<mfrac>
<mi>L</mi>
<msub>
<mi>L</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>&times;</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>s</mi>
</mrow>
</mfrac>
</mrow>
Wherein SΔtFor the acceleration area lattice number in unit time Δ t, LtFor Fetal Heart Rate effective time total length, Δ s supervises for fetal rhythm
The area of the single acceleration grid of shield, a length of specific duration of the grid, a width of specific heart rate value；
The acceleration parameter output module, the above-mentioned fetal heart monitoring acceleration parameter result calculated is exported.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510420283.XA CN105030231B (en) | 2015-07-16 | 2015-07-16 | A kind of fetal heart monitoring accelerates the extracting method of characteristic parameter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510420283.XA CN105030231B (en) | 2015-07-16 | 2015-07-16 | A kind of fetal heart monitoring accelerates the extracting method of characteristic parameter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105030231A CN105030231A (en) | 2015-11-11 |
| CN105030231B true CN105030231B (en) | 2017-10-13 |
Family
ID=54437638
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510420283.XA Active CN105030231B (en) | 2015-07-16 | 2015-07-16 | A kind of fetal heart monitoring accelerates the extracting method of characteristic parameter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105030231B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106073746B (en) * | 2016-06-15 | 2018-10-19 | 北京工业大学 | Fetus quiet sleep periodicity extraction method based on fetal heart frequency curve |
| CN108261176B (en) * | 2017-12-28 | 2021-06-04 | 深圳京柏医疗科技股份有限公司 | Fetal heart monitoring data processing method, device, system, storage medium and computer equipment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102302363A (en) * | 2011-06-29 | 2012-01-04 | 深圳市理邦精密仪器股份有限公司 | Device and method for automatically identifying fetal heart rate baseline |
| CN102319063A (en) * | 2011-07-26 | 2012-01-18 | 深圳市理邦精密仪器股份有限公司 | Device and method for improving acceleration identification accuracy of fetal heart rate data |
| CN102405014A (en) * | 2008-12-29 | 2012-04-04 | 马克.埃文斯 | Identifying the level of fetal risk during labor |
| CN104287711A (en) * | 2014-09-24 | 2015-01-21 | 广州三瑞医疗器械有限公司 | Methods for calculating non-baseline part and baseline of fetal heart rate curve |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2447861A1 (en) * | 2002-11-01 | 2004-05-01 | Lms Medical Systems Ltd. | Method and apparatus for identifying heart rate feature events |
| US20060229518A1 (en) * | 2005-03-31 | 2006-10-12 | Physical Logic Ag | Method for continuous, non-invasive, non-radiating detection of fetal heart arrhythmia |
-
2015
- 2015-07-16 CN CN201510420283.XA patent/CN105030231B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102405014A (en) * | 2008-12-29 | 2012-04-04 | 马克.埃文斯 | Identifying the level of fetal risk during labor |
| CN102302363A (en) * | 2011-06-29 | 2012-01-04 | 深圳市理邦精密仪器股份有限公司 | Device and method for automatically identifying fetal heart rate baseline |
| CN102319063A (en) * | 2011-07-26 | 2012-01-18 | 深圳市理邦精密仪器股份有限公司 | Device and method for improving acceleration identification accuracy of fetal heart rate data |
| CN104287711A (en) * | 2014-09-24 | 2015-01-21 | 广州三瑞医疗器械有限公司 | Methods for calculating non-baseline part and baseline of fetal heart rate curve |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105030231A (en) | 2015-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102302363B (en) | Device and method for automatically identifying fetal heart rate baseline | |
| CN102319063B (en) | Device and method for improving acceleration identification accuracy of fetal heart rate data | |
| CN102073059B (en) | Digital pileup waveform processing method and system | |
| CN103565433B (en) | Method and device for improving fetal monitoring efficiency | |
| CN103323615B (en) | A kind of mobile terminal and method being calculated walking speed by acceleration transducer | |
| CN109858522A (en) | A kind of management line loss abnormality recognition method based on data mining | |
| CN104287711B (en) | The computational methods of the non-baseline of fetal heart rate curve and baseline | |
| CN108261176A (en) | Fetal heart monitoring data processing method, device, system, storage medium and computer equipment | |
| CN105030231B (en) | A kind of fetal heart monitoring accelerates the extracting method of characteristic parameter | |
| CN109977885A (en) | A kind of people's vehicle automatic identifying method and device based on Doppler Feature | |
| CN102090883A (en) | Automatic identification method and device for fetal movement | |
| CN111898443B (en) | Flow monitoring method for wire feeding mechanism of FDM type 3D printer | |
| CN106073746A (en) | Fetus quiet sleep periodicity extraction method based on fetal heart frequency curve | |
| CN104546007A (en) | Anti-interference processing method and device for fetal movement detection | |
| CN101995437A (en) | Method for extracting features of crack acoustic emission signal of drawing part | |
| JP2017073021A (en) | Driving support device | |
| CN107688820B (en) | Elevator fault diagnosis method based on BCSA optimized support vector machine | |
| CN110322870A (en) | A kind of Chinese speech signal cutting method and device | |
| CN109567869B (en) | Method and system for processing accelerated activities on fetal heart rate curve | |
| CN107436363B (en) | A kind of rail traffic vehicles speed dynamic measurement method | |
| CN105193451B (en) | A kind of recognition methods of fetal heart frequency curve baseline | |
| CN106033643A (en) | Data processing method and data processing device | |
| CN113633296A (en) | Reaction time prediction model construction method, device, equipment and readable storage medium | |
| CN107527022B (en) | Signal period calculation method | |
| Lu et al. | Nonlinear baseline estimation of FHR signal using empirical mode decomposition |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20170731 Address after: 100124 Chaoyang District, Beijing Ping Park, No. 100 Applicant after: Beijing University of Technology Applicant after: The Third Affiliated Hospital of Peking University Address before: 100124 Chaoyang District, Beijing Ping Park, No. 100 Applicant before: Beijing University of Technology |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zhang Song Inventor after: Hao Dongmei Inventor after: Yang Yimin Inventor after: Li Xuwen Inventor after: Li Guangfei Inventor after: Wang Yan Inventor after: Zhao Yangyu Inventor after: Yang Lin Inventor after: Li Shufang Inventor after: Chen Lian Inventor after: Wei Bin Inventor after: Zhu Liping Inventor before: Zhang Song Inventor before: Li Guangfei Inventor before: Yang Lin Inventor before: Wang Yan Inventor before: Li Shufang Inventor before: Yang Yimin Inventor before: Li Xuwen |