CN109567867A - A kind of processing method of Fetal Heart Rate amplitude-variation - Google Patents
A kind of processing method of Fetal Heart Rate amplitude-variation Download PDFInfo
- Publication number
- CN109567867A CN109567867A CN201811193762.2A CN201811193762A CN109567867A CN 109567867 A CN109567867 A CN 109567867A CN 201811193762 A CN201811193762 A CN 201811193762A CN 109567867 A CN109567867 A CN 109567867A
- Authority
- CN
- China
- Prior art keywords
- variation
- heart rate
- fetal heart
- peak
- amplitude
- 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.)
- Pending
Links
- 210000002458 fetal heart Anatomy 0.000 title claims abstract description 170
- 238000003672 processing method Methods 0.000 title claims abstract description 18
- 230000000694 effects Effects 0.000 claims abstract description 22
- 238000004458 analytical method Methods 0.000 claims abstract description 19
- 230000001605 fetal effect Effects 0.000 claims abstract description 18
- 230000008859 change Effects 0.000 claims abstract description 15
- 230000001133 acceleration Effects 0.000 claims abstract description 14
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 238000011084 recovery Methods 0.000 claims description 16
- 230000033764 rhythmic process Effects 0.000 claims description 9
- 230000033001 locomotion Effects 0.000 claims description 8
- 238000005311 autocorrelation function Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 208000031427 Foetal heart rate deceleration Diseases 0.000 claims description 3
- 238000012806 monitoring device Methods 0.000 claims description 3
- 230000035935 pregnancy Effects 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 19
- 210000004291 uterus Anatomy 0.000 abstract description 7
- 238000004364 calculation method Methods 0.000 abstract description 2
- 210000003754 fetus Anatomy 0.000 description 11
- 238000012544 monitoring process Methods 0.000 description 11
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 239000000523 sample Substances 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 208000006893 Fetal Hypoxia Diseases 0.000 description 1
- 208000036029 Uterine contractions during pregnancy Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000009602 intrauterine growth Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/02—Measuring pulse or heart rate
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Cardiology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Technical solution of the present invention includes a kind of processing method of Fetal Heart Rate amplitude-variation, includes the following steps: S1, acquires the Fetal Heart Rate of preset duration and be depicted as initial fetal heart rate curve；S2, data prediction is carried out to initial fetal heart rate curve, obtains effective fetal heart rate curve；S3, it distinguishes the thick variation in effective fetal heart rate curve and carefully makes a variation and extract thin variation data；S4, amplitude-variation curve is obtained to the thin variation data progress Sinusoid Analysis of extraction；Highest heart rate and the difference of HR min obtain amplitude-variation result in S5, calculated amplitude variation curve；Wherein, the thick variation includes the acceleration/deceleration activity change of Fetal Heart Rate, and the thin variation includes period variation and the amplitude-variation of fetal heart rate curve.The invention has the benefit that the accuracy rate of the analytical calculation to the amplitude-variation of fetal heart rate curve is improved without excessive artificial participation, so that the clinically accuracy to the judgement of fetal in utero situation is improved, and this method is simple, it is easy to accomplish.
Description
Technical field
The present invention relates to a kind of processing methods of Fetal Heart Rate amplitude-variation, belong to field of medical technology.
Background technique
With concern and attention of the people to fetal in utero health and safety, fetal monitoring in gynecotokology and antenatal faces
In occupation of critical role on bed.Before electronic fatal monitoring birth, people carry out the heart to fetus using traditional receiver mode
Hop-information obtains, and judges the situation in intrauterine growth with this.There is certain limitations for the mode of traditional acquisition Fetal Heart Rate
And uncertainty, especially in fetal anoxia state, some variations can occur for Fetal Heart Rate, be possible to that this can be missed with receiver mode
A little important informations, moreover different people is different to the identification capability of heart of fetus sound, influences whether the acquisition of information.With electricity
The birth and development of sub- fetal heart monitoring technology, relevant many technological means have evolved into it is ripe, especially the monitoring of fetal heart frequency,
Detection and calculating aspect, it is more superior than traditional means much in real time, accurately to reflect the variation of fetal heart frequency,
Improve and eliminate the limitation and uncertainty of traditional approach.
The thin variation variation of fetal heart rate data curve includes period variation and amplitude-variation, meanwhile, period and amplitude are also
Two component parts of one wavy curve.Amplitude-variation is also referred to as wobble amplitude or amplitude variation, according to " clinical detection "
In method, the measurement of wobble amplitude be by whole section of Fetal Heart Rate describe selected on paper one section not only without accelerating again without deceleration and wave
It is dynamic to compare as sinusoidal one minute curved section, see that the difference of its highest heart rate and HR min obtains.Judge Fetal Heart Rate song
The macroscopic wave swung up and down on line, plays a multiplier role to the analysis and calculating of amplitude-variation.
Existing electronic fatal monitoring technology, by the analysis and calculating to fetal heart rate curve amplitude-variation, by result with
Quantifiable form shows that this is clinically to the research important role of medical staff's Fetal Heart Rate.Current amplitude
There are two types of modes for the interpretation and calculating of variation: the first is the judgment mode of artificial experience；For this kind of mode, due in reality
In clinical data, the swing wave of similar sinusoidal variations is not obvious on fetal heart rate data curve, some are even difficult with the naked eye straight
The number for judging fluctuation is connect, difficulty is brought to the calculating of amplitude-variation, also increases the difficulty and complexity of algorithm accordingly
Degree.Second is the algorithm judgement based on coherent signal processing technique, Fetal Heart Rate is digitized, at corresponding algorithm
Reason, analyzes amplitude-variation and is calculated, such mode is widely used in a variety of electronic fatal monitoring products；But it is existing
Fetal heart monitoring product, there is different limitations, algorithms for processing and computing capability to fetal heart rate curve amplitude-variation
The excellent difference of performance is different, and the accuracy of the calculating level and judgement that result in amplitude-variation is irregular.
Summary of the invention
To solve the above problems, the purpose of the present invention is to provide a kind of processing methods of Fetal Heart Rate amplitude-variation, to prison
The fetal heart rate data of 20 minutes durations of shield is analyzed, and using certain algorithm process process, analyzes fetal heart rate curve
Situation of change judges the variation tendency for swinging wave on curve, calculates the numerical value of amplitude-variation；This method knows amplitude-variation
Other accuracy rate is high and algorithm is simple.
Technical solution used by the present invention solves the problems, such as it is:
A kind of processing method of Fetal Heart Rate amplitude-variation, includes the following steps:
S1, the Fetal Heart Rate for acquiring preset duration are simultaneously depicted as initial fetal heart rate curve；
S2, data prediction is carried out to initial fetal heart rate curve, obtains effective fetal heart rate curve；
S3, it distinguishes the thick variation in effective fetal heart rate curve and carefully makes a variation and extract thin variation data；
S4, amplitude-variation curve is obtained to the thin variation data progress Sinusoid Analysis of extraction；
Highest heart rate and the difference of HR min obtain amplitude-variation result in S5, calculated amplitude variation curve；
Wherein, thick variation includes the acceleration/deceleration activity change of Fetal Heart Rate, and thin variation includes becoming in the period of fetal heart rate curve
Different and amplitude-variation.
Further, step S1 further include: acquire the intrauterine fetal heartbeat of preset duration using fetal rhythm monitoring device
Voice signal；The voice signal of analysis identification acquisition, and handled to obtain initial fetal heart rate data by auto-correlation function；It will
Initial fetal heart rate data is depicted as initial fetal heart rate curve.
Further, step S2 further include: obtained in initial fetal heart rate curve according to preset effective Fetal Heart Rate value range
Effective fetal heart rate data；Low-pass filtering is carried out to effective fetal heart rate data and obtains effective fetal heart rate curve.
Further, step S3 further include: the first peak in effective fetal heart rate curve is obtained according to default first amplitude threshold
Valley relevant information；The first peak-to-valley value relevant information is analyzed and determined using thick variation features, obtains the data segment that slightly makes a variation；
The data segment other than the thick variation data segment in effective fetal heart rate curve is extracted, the data that carefully make a variation are obtained.
Specifically, the first peak-to-valley value relevant information corresponds to starting point and end point, first including the first peak value, the first peak value
Peak rise-time, the first peak value recovery time, the first peak value ascensional range, the first valley, the first valley correspond to starting point with
End point, the first valley fall time, the first valley recovery time and the first valley fall.
Further, the step of the first peak-to-valley value relevant information being analyzed and determined using thick variation features, comprising: by
One peak rise-time, the first peak value recovery time, the first peak value ascensional range add with corresponding Fetal Heart Rate in thick variation features
Fast active characteristics compare, and judge whether the fetal heart rate variation at first peak value belongs to the Fetal Heart Rate slightly to make a variation and accelerate to live
It is dynamic；First valley fall time, the first valley recovery time and the first valley fall is corresponding with thick variation features
Fetal Heart Rate deceleration active characteristics compare, and judge whether the fetal heart rate variation at first valley belongs to the Fetal Heart Rate slightly to make a variation
Deceleration activity；Wherein, thick variation features are depending on pregnancy time.
Further, step S4 further include: thin variance in effective fetal heart rate curve is searched according to default second amplitude threshold
According to the second peak value；Calculate adjacent second peak value peak it is peak-to-peak away from；According to the second peak value and peak it is peak-to-peak away from carry out periodicity analysis,
Judgement identification sine wave.
Further, step S4 further include: thin variance in effective fetal heart rate curve is searched according to default second amplitude threshold
According to the neighbouring local valley of the second peak value and each second peak value；To the data segment between the second peak value and local valley
It carries out curve fitting, judges that its variation tendency is curvilinear motion or linear change；It is the data of curvilinear motion to variation tendency
Section calculate adjacent second peak value peak it is peak-to-peak away from；Judge that obtained peak is peak-to-peak away from whether meeting cyclically-varying feature；Extract peak
It is peak-to-peak away from periodically variable data segment is met, obtain meeting the amplitude-variation curve of sine wave variation.
Further, step S5 further include: each section of obtained amplitude-variation curve is analyzed and processed, is respectively obtained every
The corresponding Magnitude Difference of section amplitude-variation curve, Magnitude Difference are the difference of maximum Fetal Heart Rate value and minimum Fetal Heart Rate value；According to
Probability statistics screen the Magnitude Difference closest to desired value, obtain amplitude-variation result.
The beneficial effects of the present invention are: a kind of processing method for Fetal Heart Rate amplitude-variation that the present invention uses, by using
Electronic equipment acquisition data are simultaneously automatically analyzed and processed the data of acquisition, without excessive artificial participation, improve to tire
The accuracy rate of the analytical calculation of the amplitude-variation of heart rate curve, to improve clinically to the accurate of fetal in utero situation judgement
Degree, and this method is simple, it is easy to accomplish.
Detailed description of the invention
Fig. 1 is total method flow diagram of the invention；
Fig. 2 is the method flow diagram of the specific embodiment of the invention one；
Fig. 3 is the partial process view of the method for the specific embodiment of the invention two；
Fig. 4 A is the exemplary diagram of the initial fetal heart rate curve C of the embodiment of the present invention；
Fig. 4 B is the exemplary diagram of the fetal heart rate curve F-1 of the embodiment of the present invention；
Fig. 4 C is the exemplary diagram of the fetal heart rate curve F-2 of the embodiment of the present invention；
Fig. 4 D is the exemplary diagram that the fetal heart rate curve of the embodiment of the present invention slightly makes a variation with thin variation data segment；
Fig. 4 E is the second peak value and the exemplary diagram of corresponding local valley of the fetal heart rate curve of the embodiment of the present invention.
Specific embodiment
In current gynecotokology and fetus clinical monitoring, monitoring and analysis to fetal heart frequency have been antenatal important
One of project.Using electronic fatal monitoring technology, the acquisition of fetal heart rate values generallys use ultrasonic Doppler technique, using super
The principle of reflection of sound wave acquires heart of fetus bounce by fetal rhythm probe placement in mother's abdomen close to heart of fetus position
Voice signal, and certain algorithm is used, signal is analyzed, heart rate value is calculated.
The clinically fetal heart rate data of general monitoring at least 10 minutes, is drawn in fetal heart rate curve for heart rate data and shows
Show, then fetal heart rate curve is analyzed.Fetal heart rate curve is to reflect fetus in the significant data information of in utero situation, the country
All it is using and by the Fischer, the improvement fetal rhythms methods of marking such as Fischer, Krebs, NST that are widely recognized as in the industry generally outside
Corresponding Fetal Heart Rate appraisal result is obtained by several important informations on analysis fetal heart rate curve.
There are the judgement and analysis to amplitude-variation in the methods of marking such as Fischer, improvement Fischer, Krebs, and
And also there is the analysis to baseline variation in " the three classes points-scoring system " generally used in recent years.The baseline variation point of fetal heart rate curve
Carefully to make a variation and slightly making a variation, thick variation refers to acceleration, the deceleration activity change of Fetal Heart Rate, and thin variation refers to the week of fetal heart rate curve
Phase variation and amplitude-variation.As it can be seen that amplitude-variation as baseline make a variation in one of the form of expression that carefully makes a variation, to amplitude-variation
Analysis and calculating already take up important position.Clinically, the amplitude-variation as directly reflection fetal heart rate curve variation,
Its situation that makes a variation is intuitive performance of the fetus in utero situation, reflects the reserve capabillity of fetus to a certain extent: if long-term
Without or there is shorter amplitude-variation, then illustrate that fetus is weak in reserve capabillity in utero, should draw attention, take when necessary
Intervening measure.
To make the objectives, technical solutions, and advantages of the present invention clearer, right in the following with reference to the drawings and specific embodiments
The present invention is described in detail.A kind of processing method of Fetal Heart Rate amplitude-variation of the invention is suitable for clinically to fetal heart frequency
Thin variation diagnosis.
As shown in Figure 1, a kind of processing method of Fetal Heart Rate amplitude-variation of the invention, includes the following steps:
S1, the Fetal Heart Rate for acquiring preset duration are simultaneously depicted as initial fetal heart rate curve；
S2, data prediction is carried out to initial fetal heart rate curve, obtains effective fetal heart rate curve；
S3, it distinguishes the thick variation in effective fetal heart rate curve and carefully makes a variation and extract thin variation data；
S4, amplitude-variation curve is obtained to the thin variation data progress Sinusoid Analysis of extraction；
Highest heart rate and the difference of HR min obtain amplitude-variation result in S5, calculated amplitude variation curve；
Wherein, thick variation includes the acceleration/deceleration activity change of Fetal Heart Rate, and thin variation includes becoming in the period of fetal heart rate curve
Different and amplitude-variation.
Referring to Fig. 2, further, step S1 further include:
S11, using fetal rhythm monitoring device acquisition preset duration intrauterine fetal heartbeat voice signal；
The voice signal that S12, analysis identification acquire, and handled to obtain initial Fetal Heart Rate number by auto-correlation function
According to；
S13, initial fetal heart rate data is depicted as initial fetal heart rate curve.
Since amplitude-variation is a kind of form of expression that fetal heart rate curve carefully makes a variation, the analysis to amplitude-variation is to be based on working as
Preceding fetal heart rate curve carries out, therefore the acquisition of fetal heart rate curve is the important prerequisite condition of calculated amplitude variation；The present embodiment
In, the voice signal of the intrauterine fetal heartbeat of 20 minutes durations is acquired using foetus ECG monitor, foetus ECG monitor is automatically right
The voice signal carries out analysis identification and is converted into digital signal, then by auto-correlation function to these digital signals at
Reason, obtains the initial fetal heart rate data of 20 minutes durations, these initial fetal heart rate data described point lines can be drawn out initially
Fetal heart rate curve, one section of exemplary initial fetal heart rate curve is as shown in Figure 4 A, here remembers the initial fetal heart rate curve in Fig. 4 A
For curve C, horizontal axis represents time shaft (or sampled point axis), and the longitudinal axis represents Fetal Heart Rate value.
Further, step S2 further include:
S21, effective fetal heart rate data in initial fetal heart rate curve is obtained according to preset effective Fetal Heart Rate value range；
S22, effective fetal heart rate curve is obtained to effective fetal heart rate data progress low-pass filtering.
Medical information regulation, Fetal Heart Rate value in a certain range are considered as virtual value, and being considered as except going beyond the scope is invalid
Heart rate value；Due to foetus ECG monitor when being acquired there may be it is some not heart of fetus bounce sound interference signal,
Such as the abnormal sudden change signal or noise signal acquired under fetal rhythm probe contacts are bad or the movable situation of lie, by this
The fetal heart rate data mutation that a little interference signals are calculated is often invalid, therefore, here first to initial fetal heart rate curve into
Line number Data preprocess rejects apparent invalid fetal heart rate data, to improve data accuracy.In the present embodiment, on curve C,
Exceed the invalid fetal heart rate data of the range according to preset normal tire heart rate value range-based searching and it is excluded, then uses
Low-pass filtering treatment filters out the invalid fetal heart rate data that noise signal obtains, then draws remaining data again described point line
Effective fetal heart rate curve is produced, curve F is denoted as.
Further, step S3 further include:
S31, the first peak-to-valley value relevant information in effective fetal heart rate curve is obtained according to default first amplitude threshold；
S32, the first peak-to-valley value relevant information is analyzed and determined using thick variation features, obtains the data segment that slightly makes a variation；
The data segment other than thick variation data segment in S33, the effective fetal heart rate curve of extraction, obtains the data that carefully make a variation.
Sometimes, when fetal in utero is by movable influences such as movement of the foetus, uterine contractions, certain variation can occur for heart rate value.
For example Fetal Heart Rate value rises, and shows on fetal heart rate curve, it may appear that the variation of some " arch " shapes, when amplitude changes and continues
Between be more than certain numerical value when, the referred to as acceleration variation activity of Fetal Heart Rate.Sometimes, fetus is influenced by anoxic, heart rate value
Also certain variation, such as the decline of Fetal Heart Rate value can occur, show on fetal heart rate curve, it may appear that some to drop to downwards centainly
The decline of valley changes, when amplitude changes and the duration is more than certain numerical value, the referred to as deceleration variation activity of Fetal Heart Rate.
Activity change and deceleration activity change is accelerated to belong to the thick variation version of fetal heart rate curve, and amplitude-variation is
The thin variation changing content of fetal heart rate curve.The two include and the data information that embodies be it is entirely different, need in fetal rhythm
It is effectively distinguished on curve.
Wherein, the first peak-to-valley value relevant information corresponds to starting point and end point, first peak including the first peak value, the first peak value
Value rise time, the first peak value recovery time, the first peak value ascensional range, the first valley, the first valley correspond to starting point and knot
Beam spot, the first valley fall time, the first valley recovery time and the first valley fall.
In the present embodiment, the first peak-to-valley value relevant information, is first arranged the first amplitude threshold in order to obtain, the first width here
The parameter setting of degree threshold value can refer to the range value range of acceleration activity, deceleration campaign definitions, according to the actual situation depending on；In song
Peak, the valley point that amplitude is more than the first amplitude threshold are searched on line F, and the size and location of these peak values is recorded in array PEAK
In, the size and location of these valleies is recorded in array VALE；And these peaks, valley point are used into circle on curve F
It is marked, such as Fig. 4 B example, obtains the first peak value and 4 the first valleies that 8 meet condition here, the curve note after label
For curve F-1；Its corresponding starting point and end point are searched near these first peak values and the first valley；Wherein, first peak
It is worth the time difference between corresponding starting point and the first peak point and amplitude difference, then respectively the first peak rise-time and first peak
It is worth ascensional range, it was then the first peak value recovery time that the first peak value, which corresponds to the time difference between end point and the first peak point,；Together
Reason, the first valley correspond to time difference and amplitude difference between starting point and the first valley point, then when respectively the first valley declines
Between and the first valley fall, it is then the first valley that the first valley, which corresponds to the time difference between end point and the first valley point,
Recovery time.
Further, the step of the first peak-to-valley value relevant information being analyzed and determined using thick variation features, comprising:
By the first peak rise-time, the first peak value recovery time, the first peak value ascensional range with it is right in thick variation features
The Fetal Heart Rate answered accelerates active characteristics to compare, and judges whether the fetal heart rate variation at first peak value belongs to the tire slightly to make a variation
Heart rate acceleration activity；
By the first valley fall time, the first valley recovery time and the first valley fall with it is right in thick variation features
The Fetal Heart Rate deceleration active characteristics answered compare, and judge whether the fetal heart rate variation at first valley belongs to the tire slightly to make a variation
Heart rate decelerations activity；
Wherein, thick variation features are depending on pregnancy time.
The first peak value and the corresponding fetal heart rate variation of the first valley only filtered out with the first amplitude threshold may be not all
Real acceleration activity or deceleration activity, because the baseline variation of some fetal heart rate curves also accelerates or subtracts there may be similar
Speed variation, i.e., can also generate peak, valley, it is therefore desirable to and the first peak value and the first valley obtain to these is screened again,
To obtain more accurately thick variation data；It is used as again here according to the actual characteristic of the acceleration or deceleration activity change slightly to make a variation
The foundation of screening, wherein thick variation features specifically:
A, pregnant 32 weeks hereinafter, the first peak rise-time > 10 second, the first peak value ascensional range > 10BPM, the first peak value are extensive
Multiple time > 10 second, for primary really acceleration activity；
B, pregnant 32 weeks or more, the first peak rise-time > 15 second, the first peak value ascensional range > 15BPM, the first peak value
Recovery time > 15 second, for primary really acceleration activity；
C, the first valley fall time > 15 second, the first valley fall > 15BPM, the first valley recovery time > 15 second,
It is referred to as primary really deceleration activity.
In the present embodiment, such as Fig. 4 C example, passed through again in curve F-1 by the first peak value of circles mark and the first valley
After secondary screening, the first peak value and 0 the first valley that 6 really meet thick variation features are obtained, and be marked with arrow, marked
Curve after note is denoted as curve F-2；It, will be where the first peak value of arrow mark in curve F-2 and the first valley such as Fig. 4 D example
Data segment is denoted as thick variation data segment, marked with EXE (1), EXE (2) ..., here because several meet thick variation features
First peak value is relatively close to therefore merging into a thick variation data segment for data segment where these first peak values, such as Fig. 4 D
In EXE (2), and the data other than the data segment that slightly makes a variation then are thin variation data, are marked with F (1), F (2) ....
Further, in embodiment one, step S4 further include:
S401, the second peak value that the data that carefully make a variation in effective fetal heart rate curve are searched according to default second amplitude threshold；
S402, the peak for calculating adjacent second peak value it is peak-to-peak away from；
S403, according to the second peak value and peak it is peak-to-peak away from carry out periodicity analysis, judge identify sine wave.
The measuring method of amplitude-variation is by the fetal heart rate curve of a period of time, selecting one section not only without acceleration but also nothing
Slow down and the curved section as sinusoidal one minute duration compared in fluctuation, then calculate its highest heart rate and HR min it
Difference obtains.Not only the data portion without accelerating but also without deceleration had obtained in above-mentioned steps, such as F (1), the F marked in Fig. 4 D
(2), how the data segment of F (3), F (4), F (5) judges that similar sinusoidal swing wave then becomes the pass of calculated amplitude variation
Key link.
Judge whether similar sinusoidal swing wave, a kind of relatively coarse method is directly by curve waveform
The peak-to-peak cyclically-varying away from judgment curves waveform in peak, be considered sine wave if having cyclically-varying feature.
Further, in embodiment two, step S4 further include:
S411, the second peak value that the data that carefully make a variation in effective fetal heart rate curve are searched according to default second amplitude threshold, with
And the local valley that each second peak value is neighbouring；
S412, it carries out curve fitting to the data segment between the second peak value and local valley, judges that its variation tendency is bent
Line variation or linear change；
S413, be to variation tendency curvilinear motion data segment calculate adjacent second peak value peak it is peak-to-peak away from；
S414, the obtained peak of judgement are peak-to-peak away from whether meeting cyclically-varying feature；
S415, extraction peak are peak-to-peak away from periodically variable data segment is met, and obtain the amplitude-variation for meeting sine wave variation
Curve.
As shown in figure 3, in addition to specific steps in embodiment two except step S4 are different, remaining step with one phase of embodiment
Together, i.e., step S11-S32 and S52 in Fig. 2 are omitted in Fig. 3；It is only that peak value and peak is peak-to-peak away from referring to as sine wave judgement
Mark, in certain curvilinear motion situations, such as the peak of sawtooth wave variation is peak-to-peak away from also having periodicity；Therefore in order to make judgement more
It is accurate to add, here peak value, peak it is peak-to-peak away from the basis of increase valley as judge index, specific judgment method is as follows:
Second amplitude threshold is set, and it is more than the second amplitude threshold that amplitude is searched in carefully variation data segment F (1)-F (5)
These peak value size and locations are kept and are recorded in array PEAKTC by peak value, such as Fig. 4 E example, find 20 here completely
Second peak value of sufficient condition is marked in Fig. 4 E with number 1-20.When handling more gentle fetal heart rate curve,
Lookup of one slightly larger amplitude threshold for peak value can be first set, if the peak value number found out is less than preset number,
Amplitude threshold is adjusted, searches peak value again, until peak value number is no less than preset number.
On one section of fetal heart rate curve that above-mentioned steps each of find near the second peak value, neighbouring minimum tire is searched
The size and location of these local valleies are kept records of in array VALETC, are being schemed as local valley by heart rate value
Corresponding second peak value is marked with 1-20 in 4E, then in Fig. 4 E, digital representation above curve is the second peak value, under curve
The digital representation of side is the corresponding local valley of the second peak value.
The data segment between the second peak value and corresponding local valley position is analyzed, judges the fetal rhythm rate score of the data segment
Variation tendency, by being fitted to data, judgement is curvilinear trend or linear trend.If plots changes, saves and remember
Second peak value and corresponding local valley are recorded, if linear trends of change, illustrates that peak-to-valley value variation herein is not sine swing
Wave trend excludes second peak value and corresponding local valley.
The peak that adjacent second peak value is calculated to the data segment that variation tendency is curvilinear motion is peak-to-peak away from being stored in digital PPD
In, it is peak-to-peak away from judging whether fetal heart rate curve waveform has cyclically-varying feature, have here by the peak of adjacent second peak value
The periodical deterministic process of body are as follows:
If the peak of adjacent second peak value is peak-to-peak away from more than disP, previous peak position is excluded；
If the peak of adjacent second peak value is peak-to-peak away from disP is not above, two peak positions are saved；Wherein disP is
Pre-set peak peak distance values.
The fetal heart rate curve at peak position preserved, the as sine wave wobble with cyclically-varying feature are bent
Wire shaped.See the example in Fig. 4 E, the fetal rhythm curve data section of similar sine wave variation is at the second peak value 3-4-5-6-7-8
Set, the second position peak value 10-11-12-13-14-15-16-17, at the second position peak value 19-20.
Further, step S5 further include:
S51, each section of obtained amplitude-variation curve is analyzed and processed, respectively obtains every section of amplitude-variation curve pair
The Magnitude Difference answered, Magnitude Difference are the difference of maximum Fetal Heart Rate value and minimum Fetal Heart Rate value；
S52, the Magnitude Difference closest to desired value is screened according to probability statistics, obtains amplitude-variation result.
In the present embodiment, in Fig. 4 E between the second peak position 3-8, between the second peak position 10-17, the second peak value
Sine at the 19-20 of position puts wave change curve part and carries out amplitude calculating, calculates separately out maximum Fetal Heart Rate value and minimum
The difference of Fetal Heart Rate value, referred to as amplitude difference keep records of as amP (1), amP (2) ...
In the present embodiment, there are several sections of sine wave delta datas of separation (similar on one section of 20 minutes fetal heart rate curve
Sample data situation in Fig. 4 E), the amplitude difference for the part that thin to entirety can make a variation here is for statistical analysis to calculate multiple width
The desired value of difference is spent, chooses the Magnitude Difference closest to the desired value as final amplitude-variation as a result, be recorded as Amp, and
This is exported as the result is shown.
The above, only presently preferred embodiments of the present invention, the invention is not limited to above embodiment, as long as
It reaches technical effect of the invention with identical means, all should belong to protection scope of the present invention.In protection model of the invention
Its technical solution and/or embodiment can have a variety of different modifications and variations in enclosing.
Claims (9)
1. a kind of processing method of Fetal Heart Rate amplitude-variation, which comprises the steps of:
S1, the Fetal Heart Rate for acquiring preset duration are simultaneously depicted as initial fetal heart rate curve；
S2, data prediction is carried out to initial fetal heart rate curve, obtains effective fetal heart rate curve；
S3, it distinguishes the thick variation in effective fetal heart rate curve and carefully makes a variation and extract thin variation data；
S4, amplitude-variation curve is obtained to the thin variation data progress Sinusoid Analysis of extraction；
Highest heart rate and the difference of HR min obtain amplitude-variation result in S5, calculated amplitude variation curve；
Wherein, the thick variation includes the acceleration/deceleration activity change of Fetal Heart Rate, and the thin variation includes the week of fetal heart rate curve
Phase variation and amplitude-variation.
2. the processing method of Fetal Heart Rate amplitude-variation according to claim 1, which is characterized in that the step S1 is also wrapped
It includes:
Use the voice signal of the intrauterine fetal heartbeat of fetal rhythm monitoring device acquisition preset duration；
The voice signal of analysis identification acquisition, and handled to obtain initial fetal heart rate data by auto-correlation function；
Initial fetal heart rate data is depicted as initial fetal heart rate curve.
3. the processing method of Fetal Heart Rate amplitude-variation according to claim 1, which is characterized in that the step S2 is also wrapped
It includes:
Effective fetal heart rate data in initial fetal heart rate curve is obtained according to preset effective Fetal Heart Rate value range；
Low-pass filtering is carried out to effective fetal heart rate data and obtains effective fetal heart rate curve.
4. the processing method of Fetal Heart Rate amplitude-variation according to claim 1, which is characterized in that the step S3 is also wrapped
It includes:
The first peak-to-valley value relevant information in effective fetal heart rate curve is obtained according to default first amplitude threshold；
The first peak-to-valley value relevant information is analyzed and determined using thick variation features, obtains the data segment that slightly makes a variation；
The data segment other than the thick variation data segment in effective fetal heart rate curve is extracted, the data that carefully make a variation are obtained.
5. the processing method of Fetal Heart Rate amplitude-variation according to claim 4, it is characterised in that: the first peak-to-valley value phase
When to close information include that first peak value, the first peak value correspond to starting point and end point, the first peak rise-time, the recovery of the first peak value
Between, the first peak value ascensional range, the first valley, the first valley correspond to starting point and end point, the first valley fall time, first
Valley recovery time and the first valley fall.
6. the processing method of Fetal Heart Rate amplitude-variation according to claim 5, which is characterized in that described special using thick variation
The step of sign analyzes and determines the first peak-to-valley value relevant information, comprising:
First peak rise-time, the first peak value recovery time, the first peak value ascensional range is corresponding with thick variation features
Fetal Heart Rate accelerates active characteristics to compare, and judges whether the fetal heart rate variation at first peak value belongs to the Fetal Heart Rate slightly to make a variation
Acceleration activity；
First valley fall time, the first valley recovery time and the first valley fall is corresponding with thick variation features
Fetal Heart Rate deceleration active characteristics compare, and judge whether the fetal heart rate variation at first valley belongs to the Fetal Heart Rate slightly to make a variation
Deceleration activity；
Wherein, the thick variation features are depending on pregnancy time.
7. the processing method of Fetal Heart Rate amplitude-variation according to claim 1, which is characterized in that the step S4 is also wrapped
It includes:
The second peak value of the data that carefully make a variation in effective fetal heart rate curve is searched according to default second amplitude threshold；
Calculate adjacent second peak value peak it is peak-to-peak away from；
It is peak-to-peak away from periodicity analysis is carried out according to the second peak value and peak, judge to identify sine wave.
8. the processing method of Fetal Heart Rate amplitude-variation according to claim 1, which is characterized in that the step S4 is also wrapped
It includes:
The second peak value and each second for the data that carefully make a variation in effective fetal heart rate curve is searched according to default second amplitude threshold
The neighbouring local valley of peak value；
Carry out curve fitting to the data segment between the second peak value and local valley, judge its variation tendency be curvilinear motion or
Linear change；
The peak that adjacent second peak value is calculated to the data segment that variation tendency is curvilinear motion it is peak-to-peak away from；
Judge that obtained peak is peak-to-peak away from whether meeting cyclically-varying feature；
Extraction peak is peak-to-peak away from periodically variable data segment is met, and obtains the amplitude-variation curve for meeting sine wave variation.
9. the processing method of Fetal Heart Rate amplitude-variation according to claim 1, which is characterized in that the step S5 is also wrapped
It includes:
Each section of obtained amplitude-variation curve is analyzed and processed, the corresponding amplitude of every section of amplitude-variation curve is respectively obtained
Difference, the Magnitude Difference are the difference of maximum Fetal Heart Rate value and minimum Fetal Heart Rate value；
According to probability statistics screening closest to the Magnitude Difference of desired value, amplitude-variation result is obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811193762.2A CN109567867A (en) | 2018-10-15 | 2018-10-15 | A kind of processing method of Fetal Heart Rate amplitude-variation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811193762.2A CN109567867A (en) | 2018-10-15 | 2018-10-15 | A kind of processing method of Fetal Heart Rate amplitude-variation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109567867A true CN109567867A (en) | 2019-04-05 |
Family
ID=65920073
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811193762.2A Pending CN109567867A (en) | 2018-10-15 | 2018-10-15 | A kind of processing method of Fetal Heart Rate amplitude-variation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109567867A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113425278A (en) * | 2021-08-02 | 2021-09-24 | 北京雪扬科技有限公司 | Method for monitoring fetal heart through data acquisition of wearable equipment |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060229518A1 (en) * | 2005-03-31 | 2006-10-12 | Physical Logic Ag | Method for continuous, non-invasive, non-radiating detection of fetal heart arrhythmia |
| CN102319063A (en) * | 2011-07-26 | 2012-01-18 | 深圳市理邦精密仪器股份有限公司 | Device and method for improving acceleration identification accuracy of fetal heart rate data |
| CN103565433A (en) * | 2013-10-30 | 2014-02-12 | 深圳市理邦精密仪器股份有限公司 | Method and device for improving fetal monitoring efficiency |
| CN103784130A (en) * | 2014-01-24 | 2014-05-14 | 深圳市理邦精密仪器股份有限公司 | Method and system for improving fetal heart rate curve efficiency |
| CN104287711A (en) * | 2014-09-24 | 2015-01-21 | 广州三瑞医疗器械有限公司 | Methods for calculating non-baseline part and baseline of fetal heart rate curve |
| CN104586378A (en) * | 2015-01-21 | 2015-05-06 | 深圳市理邦精密仪器股份有限公司 | Method and device for outputting fetal heart rate curve |
| CN108261176A (en) * | 2017-12-28 | 2018-07-10 | 深圳京柏医疗科技股份有限公司 | Fetal heart monitoring data processing method, device, system, storage medium and computer equipment |
-
2018
- 2018-10-15 CN CN201811193762.2A patent/CN109567867A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060229518A1 (en) * | 2005-03-31 | 2006-10-12 | Physical Logic Ag | Method for continuous, non-invasive, non-radiating detection of fetal heart arrhythmia |
| CN102319063A (en) * | 2011-07-26 | 2012-01-18 | 深圳市理邦精密仪器股份有限公司 | Device and method for improving acceleration identification accuracy of fetal heart rate data |
| CN103565433A (en) * | 2013-10-30 | 2014-02-12 | 深圳市理邦精密仪器股份有限公司 | Method and device for improving fetal monitoring efficiency |
| CN103784130A (en) * | 2014-01-24 | 2014-05-14 | 深圳市理邦精密仪器股份有限公司 | Method and system for improving fetal heart rate curve efficiency |
| CN104287711A (en) * | 2014-09-24 | 2015-01-21 | 广州三瑞医疗器械有限公司 | Methods for calculating non-baseline part and baseline of fetal heart rate curve |
| CN104586378A (en) * | 2015-01-21 | 2015-05-06 | 深圳市理邦精密仪器股份有限公司 | Method and device for outputting fetal heart rate curve |
| CN108261176A (en) * | 2017-12-28 | 2018-07-10 | 深圳京柏医疗科技股份有限公司 | Fetal heart monitoring data processing method, device, system, storage medium and computer equipment |
Non-Patent Citations (2)
| Title |
|---|
| 刘晶影等: "《现代临床妇产科学》", 30 September 2018, 天津出版传媒集团天津科学技术出版社 * |
| 周郅隆等: "《高危妊娠的监护与处理》", 31 May 1999, 上海科技教育出版社 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113425278A (en) * | 2021-08-02 | 2021-09-24 | 北京雪扬科技有限公司 | Method for monitoring fetal heart through data acquisition of wearable equipment |
| CN113425278B (en) * | 2021-08-02 | 2023-10-20 | 北京雪扬科技有限公司 | A method of collecting data to monitor fetal heart rate through wearable devices |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5271718B2 (en) | How to identify fetal and maternal ECGs across multiple time segments | |
| AU2004266555B2 (en) | EEG seizure analysis | |
| US8116855B2 (en) | Fetal surveillance | |
| EP0645035B1 (en) | Medical signal analyzer | |
| CN109567868A (en) | A kind of CTG fetal rhythm methods of marking and system | |
| JP7293778B2 (en) | NOISE JUDGMENT METHOD, NOISE JUDGMENT PROGRAM AND NOISE JUDGMENT DEVICE | |
| WO2023226223A1 (en) | Ppg signal quality evaluation method and apparatus and ppg signal processing method and system | |
| CN109567867A (en) | A kind of processing method of Fetal Heart Rate amplitude-variation | |
| CN110236594B (en) | A fetal movement signal identification and marking method, system, medium and equipment | |
| CN109589137B (en) | Fetal movement identification method, fetal movement identification device, terminal and computer-readable storage medium | |
| Nidhal et al. | A novel cardiotocography fetal heart rate baseline estimation algorithm | |
| Chang et al. | Validation of automated arrhythmia detection for Holter ECG | |
| CN109567866A (en) | A kind of processing method of Fetal Heart Rate period variation | |
| Wei et al. | Fetal heart rate analysis using a non-linear baseline and variability estimation method | |
| Alyousif et al. | Rule-based algorithm for intrapartum cardiotocograph pattern features extraction and classification | |
| CN116831589A (en) | An optimal lead signal selection method for long-term wearable ECG signals | |
| CN110292388A (en) | A kind of measurement method and terminal of cognitive load and psychological pressure | |
| JPH0622915A (en) | Heart arrhythmia detecting device based on ecg signal | |
| Pani et al. | A real-time algorithm for tracking of foetal ECG sources obtained by block-on-line BSS techniques | |
| CN114601439B (en) | Method for detecting fetal state | |
| CN120585301A (en) | Heart rate stability assessment method | |
| Li et al. | R-peak detection for ECG signal based on local maximums of signal magnitude and correlation | |
| Lweesy et al. | Extraction of fetal heart rate and fetal heart rate variability from mother's ECG signal | |
| Wang et al. | Research on R-wave localization algorithm based on improved Pan-Thompkins | |
| CN121003759A (en) | A smart sleep device and its control method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Yan Jinyuan Inventor after: Zhang Daoguo Inventor after: Liang Jin Inventor after: Kong Fengqiang Inventor before: Kong Fengqiang |
|
| CB03 | Change of inventor or designer information | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190405 |
|
| RJ01 | Rejection of invention patent application after publication |