WO2018023697A1 - Fetal electrocardiosignal separation method and device - Google Patents
Fetal electrocardiosignal separation method and device Download PDFInfo
- Publication number
- WO2018023697A1 WO2018023697A1 PCT/CN2016/093525 CN2016093525W WO2018023697A1 WO 2018023697 A1 WO2018023697 A1 WO 2018023697A1 CN 2016093525 W CN2016093525 W CN 2016093525W WO 2018023697 A1 WO2018023697 A1 WO 2018023697A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- signal
- ecg
- abdominal wall
- qrs complex
- fetal
- Prior art date
Links
- 230000001605 fetal effect Effects 0.000 title claims abstract description 111
- 238000000926 separation method Methods 0.000 title claims abstract description 45
- 238000000718 qrs complex Methods 0.000 claims abstract description 124
- 210000003815 abdominal wall Anatomy 0.000 claims abstract description 100
- 230000008774 maternal effect Effects 0.000 claims abstract description 65
- 230000003044 adaptive effect Effects 0.000 claims abstract description 47
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 40
- 238000012545 processing Methods 0.000 claims abstract description 17
- 238000012935 Averaging Methods 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims description 34
- 238000009499 grossing Methods 0.000 claims description 26
- 238000012937 correction Methods 0.000 claims description 22
- 238000001514 detection method Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 description 12
- 230000000875 corresponding effect Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 10
- 238000004590 computer program Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 238000000605 extraction Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000013528 artificial neural network Methods 0.000 description 4
- 210000003754 fetus Anatomy 0.000 description 4
- 230000009984 peri-natal effect Effects 0.000 description 4
- 238000012549 training Methods 0.000 description 4
- 230000003187 abdominal effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 210000002458 fetal heart Anatomy 0.000 description 3
- 210000001015 abdomen Anatomy 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 210000004381 amniotic fluid Anatomy 0.000 description 2
- 239000004020 conductor Substances 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012880 independent component analysis Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 208000002330 Congenital Heart Defects Diseases 0.000 description 1
- 206010016855 Foetal distress syndrome Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 208000028831 congenital heart disease Diseases 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000009429 distress Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000003183 myoelectrical effect Effects 0.000 description 1
- 206010033675 panniculitis Diseases 0.000 description 1
- 230000001766 physiological effect Effects 0.000 description 1
- 210000002826 placenta Anatomy 0.000 description 1
- 230000035935 pregnancy Effects 0.000 description 1
- 230000005654 stationary process Effects 0.000 description 1
- 210000004003 subcutaneous fat Anatomy 0.000 description 1
- 210000004304 subcutaneous tissue Anatomy 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 210000004291 uterus Anatomy 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/344—Foetal cardiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
Definitions
- the invention relates to the field of biomedical engineering technology, in particular to a method and a device for separating fetal ECG, in particular to a method and a device for separating fetal ECG based on a smooth window.
- Fetal electrocardiogram is one of the most important objective indicators of intrauterine physiological activity. By analyzing fetal ECG signals, it can diagnose fetal intrauterine hypoxia, intrauterine distress and congenital heart disease. Effectively reduce the morbidity and mortality of various types of diseases in newborns. Therefore, the clinical application of the perinatal monitor will be one of the powerful measures to ensure the safety of the fetus at all stages of the perinatal period and to ensure the eugenics of the fetus.
- the fetus is surrounded by several layers of different conductivity in the mother's abdomen, such as the placenta, amniotic fluid, maternal uterus and the subcutaneous tissue of the mother's abdomen.
- amniotic fluid has the highest electrical conductivity
- fetal fat has the lowest electrical conductivity.
- the skin and subcutaneous fat of pregnant women also have weak electrical conductivity, but the electrical conductivity is about ten times smaller than that of muscle tissue.
- the fetal ECG signal is extracted from the maternal abdominal wall mixed signal to realize the monitoring of perinatal fetal life information. Is currently a widely studied and concerned issue.
- Matched filtering mainly uses statistical characteristics such as autocorrelation and cross-correlation of signals, but the situation to be faced in the actual environment is much more complicated than the assumption, and the fetal ECG signal obtained by matched filtering is used. Not very ideal.
- ICA Independent Component Analysis
- Wavelet threshold denoising fetal ECG signal extraction method filtering Using wavelet transform combined with prior signal knowledge can obtain a good fetal ECG signal, but the extraction process is susceptible to amplitude and fetal ECG signal-like interference. It often leads to the occurrence of R wave miss detection and misdetection in the fetal ECG signal. It will cause errors when calculating the fetal heart rate, which is not conducive to clinical diagnosis.
- Neural network method This method can detect the QRS complex of ECG signals according to the prediction error of the neural network.
- the PT segment signal is flat and can be predicted better.
- the QRS complex the difference between adjacent points is large, and the prediction error will be generated.
- the QRS complex is predicted, and the artificial neural network is predicted.
- the training sample is required to have a certain representativeness, which is difficult to realize in the actual abdominal wall ECG signal denoising application.
- the adaptive filtering method includes training the adaptive filtering to eliminate the mother electrocardiogram using one or several reference channels.
- the existing adaptive filtering method removes the parent ECG component and requires a parent ECG reference channel that is morphologically similar to the submerged FECG signal, but the existing adaptive cancellation method requires an additional reference signal, that is, in addition to the abdominal wall electrode An additional maternal chest lead electrode is also required, and an adaptive method is used to eliminate the maternal component associated with the chest lead and the abdominal wall lead, which means that the maternal abdominal wall ECG signal and the maternal chest lead ECG signal need to be correlated, Only the better pair can be eliminated, and the conduction path of the maternal ECG signal is very complicated.
- the maternal ECG signal in the abdominal wall is poorly correlated with the maternal ECG signal in the chest lead. Therefore, using the existing adaptive filtering algorithm to cancel the extraction of fetal ECG signals will have a large number of maternal ECG leakage, which will cause great interference to the correct identification of fetal ECG components. Moreover, the existing adaptive filtering cancellation cannot achieve single-lead fetal ECG signal extraction.
- the invention provides a method and a device for separating fetal electrocardiogram, which are used to solve the problem that the existing fetal electrocardiographic separation algorithm is complicated in calculation and cannot realize single lead.
- an embodiment of the present invention provides a method for separating a fetal electrocardiogram, comprising: performing QRS complex detection on a received abdominal wall electrocardiographic signal, determining a starting position of a QRS complex and a QRS complex; a start position; smoothing a signal segment between the QRS complexes by using a smoothing window; performing linear averaging on the QRS complex to generate an average template; and according to the parent QRS complex and the current center in the average template Electrical cycle correspondence Correlation of the parent QRS complex, correcting the dynamic amplitude of the parent QRS complex in the average template, generating a maternal abdominal wall ECG estimation signal; using the maternal abdominal wall ECG estimation signal as a reference signal, The maternal abdominal wall ECG signal is used as a desired signal, and the fetal ECG signal is separated from the maternal abdominal wall ECG signal by an adaptive algorithm.
- an embodiment of the present invention further provides a fetal electrocardiographic separation device, comprising: a wave group detecting module, configured to perform QRS complex detection on a received parent abdominal wall electrocardiographic signal, and determine a start of a QRS complex; a starting position between the position and the QRS complex; a smoothing module for smoothing the signal segments between the QRS complexes using a smoothing window; and an averaging processing module for linearly averaging the QRS complexes, Generating an average template; a group correction module, configured to: according to the correlation between the parent QRS complex in the average template and the parent QRS complex corresponding to the current ECG cycle, the parent QRS complex in the average template The dynamic amplitude correction is performed to generate a maternal abdominal wall ECG estimation signal; the separation module is configured to use the maternal abdominal wall ECG estimation signal as a reference signal, and use the maternal abdominal wall ECG signal as a desired signal through an adaptive algorithm.
- the fetal ECG signal is separated from the maternal abdominal wall
- the fetal electrocardiographic separation method and device proposed by the invention can eliminate the fetal component of the signal segment between the parent QRS groups by smoothing the signal between the parent QRS groups, and correct the amplitude of the parent QRS group to estimate the maternal abdominal wall electrocardiogram.
- this method avoids the incomplete influence of the maternal electrocardiogram caused by the uncorrelated maternal component in the chest lead and the abdominal wall lead, and has better separation effect; the realization of the present invention only needs one parent abdominal wall electrocardiogram signal Therefore, the single-lead fetal ECG signal separation can be completed; and the fetal separation algorithm design method of the present invention is simple and intuitive, and it is easy to find the cause and add the correction module for the abnormal clinical data.
- FIG. 1 is a schematic structural diagram of a system for a fetal electrocardiographic separation method according to an embodiment of the present invention
- FIG. 2 is a schematic block diagram of a mother window electrocardiographic estimation mechanism based on a smooth window according to an embodiment of the present invention
- FIG. 3 is a flowchart of processing of a fetal electrocardiogram separation method according to an embodiment of the present invention
- FIG. 4 is a schematic structural view of a fetal electrocardiographic separation device according to an embodiment of the present invention.
- FIG. 5 is a schematic structural diagram of a group correction module 104 of the embodiment shown in FIG. 4;
- FIG. 6 is a waveform diagram of a 5-channel mother abdominal wall mixed ECG signal used in a specific embodiment of the present invention.
- Figure 7 is an analysis of experimental results of a specific embodiment of the present invention.
- embodiments of the present invention can be implemented as a system, apparatus, device, method, or computer program product. Accordingly, the present disclosure may be embodied in the form of full hardware, complete software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.
- the invention provides an adaptive single-lead fetal electrocardiographic separation method and device based on smooth window, which uses a smoothing window to smooth the QRS inter-group signals in the maternal abdominal wall signal to smooth out the QRS inter-group signal segment.
- the fetal heartbeat beat point is extracted from the parent QRS complex in the abdominal wall ECG signal.
- the extracted N consecutive parent QRS complexes are linearly averaged by the parent R wave peak position, and the average QRS wave is generated.
- the group template intercepts the parent QRS complex in the input abdominal wall signal and the parent QRS complex in the average template, and uses the correlation coefficients of the two to copy and correct the parent QRS complex in the average template to obtain a dynamic matrix.
- Abdominal wall ECG estimation signal is used to obtain a dynamic matrix.
- the maternal component in the abdominal wall signal is eliminated, and the scheme avoids the influence of the maternal component leakage caused by the maternal component of the chest lead and the non-correlation of the maternal component in the abdominal wall lead, and can realize the single Lead fetal ECG separation.
- FIG. 1 is a schematic structural diagram of a system of a fetal electrocardiographic separation method according to an embodiment of the present invention.
- y(n) represents the acquired abdominal wall ECG signal
- f(n) is the fetal ECG signal in the abdominal wall signal
- x(n) is the maternal ECG signal in the abdominal wall signal.
- the desired signal of the adaptive filter, x'(n) is the abdominal wall of the abdominal wall estimated by the abdominal wall ECG estimation mechanism, and the adaptive filter continuously "learns" and "tracks” the input signal to make the output signal Maximum approximation of the abdominal wall ECG signal, thus making the error signal Minimal, This is the separated fetal ECG signal.
- two input signals are needed, namely: a parent abdominal wall ECG signal (main input signal), and an estimated abdominal wall ECG signal (reference input signal).
- a parent abdominal wall ECG signal main input signal
- an estimated abdominal wall ECG signal reference input signal
- the core of the invention is how to calculate the abdominal wall mother ECG estimation signal as a reference input signal to the adaptive algorithm.
- 2 is a schematic block diagram of the estimation mechanism of the abdominal wall mother body of the present invention.
- the QRS complex signal and the QRS complex signal are separately processed, and the smoothing window based method is used to perform the signal between the parent QRS complexes.
- the smoothing process eliminates the fetal ECG beat point of the signal segment between the parent QRS complexes, and uses the average correlation method to correct the amplitude of the parent QRS complex to estimate the abdominal wall signal waveform.
- FIG. 3 is a flowchart of processing of a fetal electrocardiographic separation method according to an embodiment of the present invention. As shown in Figure 3, it includes:
- Step S101 performing QRS complex detection on the received abdominal wall ECG signal to determine a starting position between the QRS complex and a starting position between the QRS complexes;
- Step S102 Smoothing a signal segment between the QRS groups by using a smoothing window
- Step S103 performing linear averaging on the QRS complex to generate an average template
- Step S104 according to the correlation between the parent QRS complex in the average template and the parent QRS complex corresponding to the current ECG cycle, correct the dynamic amplitude of the parent QRS complex in the average template to generate a matrix.
- Step S105 using the mother abdominal wall ECG estimation signal as a reference signal, using the mother abdominal wall ECG signal as a desired signal, and separating the fetal ECG signal from the mother abdominal wall ECG signal by an adaptive algorithm.
- step S102 a signal segment between the QRS complexes is smoothed by using a fixed length window, and the value of the fixed length window is determined according to a signal sampling rate and a fetal ECG group length.
- the signal segment between QRS groups is smoothed by a fixed length window M.
- the length of the fixed window is about the length of the fetal ECG QRS wave, and the length of the fetal ECG group is generally 47ms to 85ms, so the value of M is It can be (0.047*d, 0.085*d), where the parameter d is the signal sampling rate.
- the QRS inter-group signal can be expressed as:
- the maternal ECG amplitude is generally more than ten times the fetal ECG amplitude. After smoothing, the maternal QRS complex almost no longer contains the fetal QRS complex.
- step S103 linearly averaging the QRS complex to generate an average segment, comprising: extracting a parent ECG R wave peak position from the parent abdominal wall ECG signal; and using the R wave peak position as A central anchor point that linearly averages the QRS complex to generate the average template.
- the segment is subjected to a maximum search of the R wave peak position to determine the final R wave peak position point.
- step S104 of the present invention performs dynamic amplitude correction on the parent QRS complex in the average template based on the correlation between the latest parent QRS complex and the parent QRS complex in the average segment.
- the method is as follows:
- the QRS complex corresponding to the current ECG cycle is used to perform dynamic amplitude correction on the parent QRS complex in the average template, and the expression used is as follows:
- the average template generated by the average segment has the signal of the corresponding parent ECG waveform in the current abdominal wall electrocardiogram, and the fetal electrocardiogram waveform is substantially eliminated, the averaged parent ECG QRS cluster segment is obtained by the correlation coefficient. Dynamic amplitude correction can be used to estimate the maternal abdominal wall ECG signal to the maximum extent.
- step S105 the abdominal electrode obtains the main input signal, and the estimated mother abdominal wall ECG signal is used as an input signal, and the signal substantially does not contain the fetal ECG signal, and only contains the noise to be eliminated and the mother ECG signal;
- the electrical signal is the desired signal.
- the adaptive algorithm eliminates the correlation between the input signal and the desired signal, and preserves the uncorrelated portion. Therefore, the fetal ECG signal separation in the abdominal wall ECG signal can be maximized.
- the adaptive algorithm may be an adaptive least mean square (LMS) algorithm, an adaptive recursive least squares (RLS) algorithm, a least squares lattice (LSL) adaptive algorithm, or the like.
- LMS adaptive least mean square
- RLS adaptive recursive least squares
- LSL least squares lattice
- FIG. 4 is a schematic structural view of a fetal electrocardiographic separation device according to an embodiment of the present invention. As shown, it includes:
- the wave group detecting module 101 is configured to perform QRS complex detection on the received parent abdominal wall electrocardiographic signal to determine a starting position between the QRS complex and the QRS complex; and a smoothing module 102 for utilizing smoothing
- the window smoothes the signal segments between the QRS groups;
- the averaging processing module 103 is configured to linearly average the QRS groups to generate an average template;
- the group correction module 104 is configured to be used according to the average template Correlation of the parent QRS complex with the parent QRS complex corresponding to the current ECG cycle, correcting the dynamic amplitude of the parent QRS complex in the average template, and generating a maternal abdominal wall ECG estimation signal;
- the separation module 105 The parental abdominal wall ECG estimation signal is used as a reference signal, and the mother abdominal wall ECG signal is used as a desired signal, and the fetal ECG signal is separated from the mother abdominal wall ECG signal by an adaptive algorithm.
- the smoothing processing module 102 performs smoothing on the signal segments between the QRS groups by using a smoothing window, and specifically includes: using a fixed length window, the signal between the QRS groups The segment is smoothed, and the value of the fixed length window is determined according to the signal sampling rate and the length of the fetal ECG group.
- the averaging processing module 103 linearly averages the QRS complex to generate an average template, including: extracting a parent ECG R wave peak from the mother abdominal wall ECG signal Position: positioning the point of the R wave peak as a center, and linearly averaging the QRS group to generate the average template.
- the group correction module 104 compares the correlation between the parent QRS group in the average template and the parent QRS group corresponding to the current ECG cycle.
- the dynamic amplitude correction is performed to generate the mother abdominal wall ECG estimation signal, as shown in FIG. 5, which specifically includes:
- the correlation coefficient calculation unit 1041 is configured to calculate a correlation between the QRS group segment corresponding to the current ECG cycle and the QRS group segment in the average template, and generate a correlation coefficient;
- the dynamic correction unit 1042 is configured to perform dynamic amplitude correction on the parent QRS group in the average template by using the QRS complex corresponding to the current ECG cycle according to the correlation coefficient to generate a maternal abdominal wall ECG estimation signal. .
- the separation module 105 separates a fetal ECG signal from the mother abdominal wall ECG signal by an adaptive algorithm, where the adaptive algorithm includes an adaptive least mean square algorithm, Adaptive recursive least squares algorithm, least squares lattice adaptive algorithm.
- FIG. 6 a waveform diagram of a mixed 5-channel maternal abdominal wall ECG signal used in this embodiment is shown.
- the data length of each signal is 2500 sample points, and the sampling rate is 250 Hz.
- the first road abdominal wall signal is the one signal with the largest telecommunication noise ratio of the fetal heart, that is, the signal of this road is the best, and the fourth and fifth abdominal wall ECG signals can not recognize the fetal heart with the naked eye.
- the position of the electric R peak which is not used by the algorithm of the present invention, is not considered by the algorithm of the present invention.
- the second path signal is selected as the algorithm verification data, and the adaptive filtering adopts the adaptive recursive least squares algorithm, and the experimental result analysis is shown in FIG. 7.
- the fetal QRS complex is far from the parent QRS complex:
- the fetal QRS complex is far from the parent QRS complex, as labeled a in Figure 7, as can be seen from Figure 7, the tire in this case
- the QRS component is extracted and the left and right parent QRS components are suppressed.
- the positional relationship between the fetal QRS complex and the parent QRS complex is the most common positional relationship in the abdominal wall ECG and the most easily extracted positional relationship.
- the fetal QRS complex is adjacent to the parent QRS complex:
- the fetal QRS complex is adjacent to the parent QRS complex, as indicated by the mark b in Fig. 7. It can be seen from Fig. 7 that the fetal QRS component is extracted in this case, and the adjacent QRS component is suppressed.
- the identification of fetal ECG QRS complexes has an impact.
- the positional relationship between the fetal QRS complex and the parent QRS complex is a small number of positional relationships in the abdominal wall ECG, and it is also difficult to extract the positional relationship.
- the fetal QRS complex is completely coincident with the parent QRS complex, as indicated by the reference c in Fig. 7, and it can be seen from Fig. 7 that the fetal QRS component is extracted in this case, and the coincident parent QRS component is suppressed.
- the positional relationship between the fetal QRS complex and the parent QRS complex is a small number of positional relationships in the abdominal wall ECG and the most difficult to extract positional relationship.
- the fetal QRS complex is far from the parent QRS complex, the fetal QRS complex is adjacent to the parent QRS complex, the fetal QRS complex, and the parent QRS complex are completely Coincident three cases, have a very good extraction effect.
- the algorithm only needs one way of the abdominal wall ECG signal to complete the adaptive fetal signal separation.
- two signals are needed, namely: one abdominal wall ECG signal (main input signal), one chest lead ECG signal (reference input signal), through the main input signal
- the maternal QRS group smoothing process is performed to eliminate the fetal component of the signal segment between the parent QRS groups, and the amplitude correction of the parent QRS group is performed to estimate the abdominal wall signal waveform, and the signal is input into the adaptive filter while being unprocessed.
- the synchronized abdominal wall ECG signal serves as the desired signal for adaptive filtering. Therefore, the present invention requires only one abdominal wall ECG signal, which greatly simplifies the system lead system.
- the separation effect is good: through the smooth and maternal QRS wave group amplitude correction processing, estimate the maternal abdominal wall ECG signal, this method avoids maternal electrocardiogram inhibition caused by the uncorrelated maternal component in the chest lead and abdominal wall lead Incomplete effects have a better separation effect.
- the algorithm design idea is simple and intuitive, for abnormal clinical data, easy to find the cause and add correction module.
- the embodiment of the present invention further provides a computer readable program, wherein when the program is executed, the program causes a computer to execute the fetal electrocardiographic separation method according to the embodiment of the present invention.
- the embodiment of the present invention further provides a storage medium storing a computer readable program, wherein the computer readable program causes the computer to execute the fetal electrocardiographic separation method according to the embodiment of the present invention.
- the above apparatus and method of the present invention may be implemented by hardware or by hardware in combination with software.
- the present invention relates to a computer readable program that, when executed by a logic component, enables the logic component to implement the apparatus or components described above, or to cause the logic component to implement the various methods described above Or steps.
- the present invention also relates to a storage medium for storing the above program, such as a hard disk, a magnetic disk, an optical disk, a DVD, a flash memory, or the like.
- embodiments of the present invention can be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or a combination of software and hardware. Moreover, the invention can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) including computer usable program code.
- computer-usable storage media including but not limited to disk storage, CD-ROM, optical storage, etc.
- the computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing device to operate in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture comprising the instruction device.
- the apparatus implements the functions specified in one or more blocks of a flow or a flow and/or block diagram of the flowchart.
- These computer program instructions can also be loaded onto a computer or other programmable data processing device such that a series of operational steps are performed on a computer or other programmable device to produce computer-implemented processing for execution on a computer or other programmable device.
- the instructions provide steps for implementing the functions specified in one or more of the flow or in a block or blocks of a flow diagram.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Provided are a fetal electrocardiosignal separation method and a device. The fetal electrocardiosignal separation method comprises: detecting QRS complex on a received maternal abdominal wall electrocardiosignal to determine a starting position of the QRS complex and a starting position between QRS complexes (S101); smooth processing a signal segment between the QRS complexes using a smooth window (S102); linear averaging the QRS complex to generate an average template (S103); according to correlation between the maternal QRS complex in the average template and maternal QRS complex corresponding to the current electrocardiocycle, correcting the dynamic amplitude of the maternal QRS complex in the average template to generate a maternal abdominal wall electrocardio estimation signal (S104); taking the maternal abdominal wall electrocardio estimation signal as a reference signal, and the maternal abdominal wall electrocardiosignal as a desired signal, separating the fetal electrocardiosignal from the maternal abdominal wall electrocardiosignal by means of adaptive algorithm (S105). The fetal electrocardiosignal separation method and the device involve a simple fetal electrocardiosignal separation algorithm, in which only one stream of maternal abdominal wall electrocardiosignal is needed, such that single lead fetal electrocardiosignal separation can be achieved.
Description
本发明涉及生物医学工程技术领域,尤其涉及一种胎儿心电分离方法及装置,具体的讲是一种基于平滑窗的胎儿心电分离方法及装置。The invention relates to the field of biomedical engineering technology, in particular to a method and a device for separating fetal ECG, in particular to a method and a device for separating fetal ECG based on a smooth window.
本部分旨在为权利要求书中陈述的本发明的实施方式提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。This section is intended to provide a background or context for the embodiments of the invention set forth in the claims. The description herein is not admitted to be prior art as it is included in this section.
胎儿心电信号(Fetal electrocardiogram,FECG)是胎儿宫内生理活动最主要的客观指标之一,通过分析胎儿心电信号,能够早期诊断胎儿宫内缺氧、宫内窘迫以及先天性心脏病等,有效地降低新生儿各类疾病的发病率和死亡率。因此,围产期监护仪的临床应用将是保证围产期各个阶段的胎儿安全,保障胎儿优生的有力措施之一。Fetal electrocardiogram (FECG) is one of the most important objective indicators of intrauterine physiological activity. By analyzing fetal ECG signals, it can diagnose fetal intrauterine hypoxia, intrauterine distress and congenital heart disease. Effectively reduce the morbidity and mortality of various types of diseases in newborns. Therefore, the clinical application of the perinatal monitor will be one of the powerful measures to ensure the safety of the fetus at all stages of the perinatal period and to ensure the eugenics of the fetus.
在围产期,胎儿在母体腹部由几个不同的电导率的包裹层包围,如胎盘、羊水、母体子宫和母体腹部皮下组织等,其中羊水的导电能力最高,而胎脂的导电能力最低,此外孕妇的皮肤和皮下脂肪也有微弱的导电性,但是导电能力比肌肉组织约小十倍。在采集胎儿心电时,由于胎儿心电信号需要经过上述导电层,本来就很微弱的胎儿心电信号更容易受外界的干扰。这些不同的组织,形成了所谓的“人体容积导体”,这个容积导体不是一个稳定的导体,其导电性和几何形状随着整个妊娠期不断改变。因此,这将导致母体腹部心电混合信号具有如下特点:1)信噪比低,胎儿心电分量太弱,仅为10~50uV左右,而母亲心电信号则可达750uV,两者可差十倍以上。肌电噪声与其它干扰也很强,信号完全淹没于噪声之中。2)时域中胎儿QRS波与母体QRS波有10~30%的重叠。3)频域中胎儿频谱与母体的频谱大部分重叠。4)信号随机性很强,是非平稳过程,特别是由于胎儿的自由移动造成胎儿心电传输通道的时变性。During the perinatal period, the fetus is surrounded by several layers of different conductivity in the mother's abdomen, such as the placenta, amniotic fluid, maternal uterus and the subcutaneous tissue of the mother's abdomen. Among them, amniotic fluid has the highest electrical conductivity, while fetal fat has the lowest electrical conductivity. In addition, the skin and subcutaneous fat of pregnant women also have weak electrical conductivity, but the electrical conductivity is about ten times smaller than that of muscle tissue. When collecting fetal ECG, since the fetal ECG signal needs to pass through the above conductive layer, the originally weak fetal ECG signal is more susceptible to external interference. These different tissues form the so-called "human body volume conductor", which is not a stable conductor whose conductivity and geometry change throughout the gestation period. Therefore, this will lead to the following characteristics of the maternal abdominal ECG mixed signal: 1) low signal-to-noise ratio, fetal ECG component is too weak, only about 10 ~ 50uV, and mother ECG signal can reach 750uV, the difference between the two Ten times more. Myoelectric noise and other disturbances are also strong, and the signal is completely submerged in the noise. 2) The fetal QRS wave in the time domain overlaps with the parent QRS wave by 10 to 30%. 3) The fetal spectrum in the frequency domain mostly overlaps with the spectrum of the mother. 4) The signal is very random and is a non-stationary process, especially due to the time-varying of the fetal ECG transmission channel due to the free movement of the fetus.
母体腹部心电混合信号的上述特点,为获得高精度的胎儿心电信号制造了比较大的障碍,因此从母体腹壁混合信号中提取出胎儿心电信号,以实现围产期胎儿生命信息的监护,是目前广泛研究和关注的问题。The above characteristics of the maternal abdominal ECG mixed signal create a relatively large obstacle for obtaining high-precision fetal ECG signals. Therefore, the fetal ECG signal is extracted from the maternal abdominal wall mixed signal to realize the monitoring of perinatal fetal life information. Is currently a widely studied and concerned issue.
目前,从母体腹壁混合信号中提取出胎儿心电信号的方法主要有以下几种:At present, there are mainly the following methods for extracting fetal ECG signals from the mixed signals of the maternal abdominal wall:
1)匹配滤波:匹配滤波主要利用的是信号的自相关和互相关等统计特性,但是在实际环境中所要面临的情况远比假设情况要复杂的多,采用匹配滤波获得的胎儿心电信号
不是很理想。1) Matched filtering: Matching filtering mainly uses statistical characteristics such as autocorrelation and cross-correlation of signals, but the situation to be faced in the actual environment is much more complicated than the assumption, and the fetal ECG signal obtained by matched filtering is used.
Not very ideal.
2)独立成分分析法(ICA):此方法需要多路腹壁信号,无法实现单导联胎儿心电信号的提取。2) Independent Component Analysis (ICA): This method requires multiple abdominal wall signals and cannot achieve single-lead fetal ECG signal extraction.
3)小波阈值去噪的胎儿心电信号提取方法滤波:利用小波变换结合先验信号知识能够获得效果不错的胎儿心电信号,但是提取过程中容易受到幅度和胎儿心电信号类似信号的干扰,往往导致胎儿心电信号中R波漏检和误检等情况产生,在计算胎儿心率时会引起误差,不利于于临床诊断。3) Wavelet threshold denoising fetal ECG signal extraction method filtering: Using wavelet transform combined with prior signal knowledge can obtain a good fetal ECG signal, but the extraction process is susceptible to amplitude and fetal ECG signal-like interference. It often leads to the occurrence of R wave miss detection and misdetection in the fetal ECG signal. It will cause errors when calculating the fetal heart rate, which is not conducive to clinical diagnosis.
4)神经网络法:此方法可以根据神经网络的预测误差检测心电信号的QRS复合波。PT段信号平缓,可以被较好的预测,而对于QRS复合波,相邻点差异较大,会产生预测误差,利用这个预测误差训练网络,对QRS复合波进行预测,预测时的人工神经网络具有多层结构,利用神经网络学习训练过程需要花费较长时间,而且要求训练样本具有一定的代表性,这在实际腹壁心电信号降噪应用中很难实现。4) Neural network method: This method can detect the QRS complex of ECG signals according to the prediction error of the neural network. The PT segment signal is flat and can be predicted better. For the QRS complex, the difference between adjacent points is large, and the prediction error will be generated. Using this prediction error training network, the QRS complex is predicted, and the artificial neural network is predicted. With a multi-layer structure, it takes a long time to learn the training process by using the neural network, and the training sample is required to have a certain representativeness, which is difficult to realize in the actual abdominal wall ECG signal denoising application.
5)自适应滤波法:自适应滤波方法包括使用一个或几个参考通道训练自适应滤波消除母体心电图。现有的自适应滤波方法去除母体ECG成分,需要形态上类似于被淹没FECG信号的母体ECG参考通道,但是现有的自适应对消方法要求额外提供参考信号,也就是说,除了腹壁电极外,还需要配备额外的母体胸导连电极,自适应方法对消掉胸导联和腹壁导联中相关的母体分量,这意味着母体腹壁心电信号和母体胸导联心电信号需要相关,才能够被较好的对消掉,并且,母体心电信号的传导路径非常复杂,某些情况下,腹壁中的母体心电信号和胸导联中的母体心电信号相关性差。因此,利用现有的自适应滤波算法对消提取胎儿心电信号会有大量的母体心电泄露,对正确识别胎儿心电分量造成了很大干扰。并且,现有的自适应滤波对消无法实现单导联的胎儿心电信号提取。5) Adaptive filtering method: The adaptive filtering method includes training the adaptive filtering to eliminate the mother electrocardiogram using one or several reference channels. The existing adaptive filtering method removes the parent ECG component and requires a parent ECG reference channel that is morphologically similar to the submerged FECG signal, but the existing adaptive cancellation method requires an additional reference signal, that is, in addition to the abdominal wall electrode An additional maternal chest lead electrode is also required, and an adaptive method is used to eliminate the maternal component associated with the chest lead and the abdominal wall lead, which means that the maternal abdominal wall ECG signal and the maternal chest lead ECG signal need to be correlated, Only the better pair can be eliminated, and the conduction path of the maternal ECG signal is very complicated. In some cases, the maternal ECG signal in the abdominal wall is poorly correlated with the maternal ECG signal in the chest lead. Therefore, using the existing adaptive filtering algorithm to cancel the extraction of fetal ECG signals will have a large number of maternal ECG leakage, which will cause great interference to the correct identification of fetal ECG components. Moreover, the existing adaptive filtering cancellation cannot achieve single-lead fetal ECG signal extraction.
发明内容Summary of the invention
本发明提出一种胎儿心电分离方法及装置,用以解决现有的胎儿心电分离算法计算复杂且无法实现单导联的问题。The invention provides a method and a device for separating fetal electrocardiogram, which are used to solve the problem that the existing fetal electrocardiographic separation algorithm is complicated in calculation and cannot realize single lead.
为了达到上述目的,本发明实施例提供一种胎儿心电分离方法,包括:对接收到的母体腹壁心电信号进行QRS波群检测,确定QRS波群的起始位置和QRS波群间的起始位置;利用平滑窗对所述QRS波群间的信号段进行平滑处理;对所述QRS波群进行线性平均,生成平均模板;根据所述平均模板中的母体QRS波群与所述当前心电周期对应
的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号;将所述母体腹壁心电估计信号作为参考信号,将所述母体腹壁心电信号作为期望信号,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号。In order to achieve the above object, an embodiment of the present invention provides a method for separating a fetal electrocardiogram, comprising: performing QRS complex detection on a received abdominal wall electrocardiographic signal, determining a starting position of a QRS complex and a QRS complex; a start position; smoothing a signal segment between the QRS complexes by using a smoothing window; performing linear averaging on the QRS complex to generate an average template; and according to the parent QRS complex and the current center in the average template Electrical cycle correspondence
Correlation of the parent QRS complex, correcting the dynamic amplitude of the parent QRS complex in the average template, generating a maternal abdominal wall ECG estimation signal; using the maternal abdominal wall ECG estimation signal as a reference signal, The maternal abdominal wall ECG signal is used as a desired signal, and the fetal ECG signal is separated from the maternal abdominal wall ECG signal by an adaptive algorithm.
为了达到上述目的,本发明实施例还提供一种胎儿心电分离装置,包括:波群检测模块,用于对接收到的母体腹壁心电信号进行QRS波群检测,确定QRS波群的起始位置和QRS波群间的起始位置;平滑处理模块,用于利用平滑窗对所述QRS波群间的信号段进行平滑处理;平均处理模块,用于对所述QRS波群进行线性平均,生成平均模板;波群修正模块,用于根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号;分离模块,用于将所述母体腹壁心电估计信号作为参考信号,将所述母体腹壁心电信号作为期望信号,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号。In order to achieve the above object, an embodiment of the present invention further provides a fetal electrocardiographic separation device, comprising: a wave group detecting module, configured to perform QRS complex detection on a received parent abdominal wall electrocardiographic signal, and determine a start of a QRS complex; a starting position between the position and the QRS complex; a smoothing module for smoothing the signal segments between the QRS complexes using a smoothing window; and an averaging processing module for linearly averaging the QRS complexes, Generating an average template; a group correction module, configured to: according to the correlation between the parent QRS complex in the average template and the parent QRS complex corresponding to the current ECG cycle, the parent QRS complex in the average template The dynamic amplitude correction is performed to generate a maternal abdominal wall ECG estimation signal; the separation module is configured to use the maternal abdominal wall ECG estimation signal as a reference signal, and use the maternal abdominal wall ECG signal as a desired signal through an adaptive algorithm. The fetal ECG signal is separated from the maternal abdominal wall ECG signal.
本发明提出的胎儿心电分离方法及装置,通过平滑母体QRS波群间信号,消除母体QRS波群间信号段的胎儿分量,以及对母体QRS波群进行幅值修正,从而估计母体腹壁心电信号,此方法避免了胸导联和腹壁导联中母体分量非相关所带来的母体心电抑制不彻底的影响,具有更好的分离效果;本发明的实现只需要一路母体腹壁心电信号,故可以完成单导联胎儿心电信号分离;并且,本发明的胎儿分离算法设计思路简单直观,针对于异常的临床数据,易于查找原因和添加修正模块。The fetal electrocardiographic separation method and device proposed by the invention can eliminate the fetal component of the signal segment between the parent QRS groups by smoothing the signal between the parent QRS groups, and correct the amplitude of the parent QRS group to estimate the maternal abdominal wall electrocardiogram. Signal, this method avoids the incomplete influence of the maternal electrocardiogram caused by the uncorrelated maternal component in the chest lead and the abdominal wall lead, and has better separation effect; the realization of the present invention only needs one parent abdominal wall electrocardiogram signal Therefore, the single-lead fetal ECG signal separation can be completed; and the fetal separation algorithm design method of the present invention is simple and intuitive, and it is easy to find the cause and add the correction module for the abnormal clinical data.
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings used in the description of the embodiments will be briefly described below. It is obvious that the drawings in the following description are only some embodiments of the present invention. Other drawings may also be obtained from those of ordinary skill in the art in light of the inventive work. In the drawing:
图1为本发明实施例的胎儿心电分离方法的系统结构示意图;1 is a schematic structural diagram of a system for a fetal electrocardiographic separation method according to an embodiment of the present invention;
图2为本发明实施例的基于平滑窗的腹壁母体心电估计机制的示意框图;2 is a schematic block diagram of a mother window electrocardiographic estimation mechanism based on a smooth window according to an embodiment of the present invention;
图3为本发明实施例的胎儿心电分离方法的处理流程图;3 is a flowchart of processing of a fetal electrocardiogram separation method according to an embodiment of the present invention;
图4为本发明实施例的胎儿心电分离装置的结构示意图;4 is a schematic structural view of a fetal electrocardiographic separation device according to an embodiment of the present invention;
图5为图4所示实施例的波群修正模块104的结构示意图;
FIG. 5 is a schematic structural diagram of a group correction module 104 of the embodiment shown in FIG. 4;
图6为本发明的具体实施例采用的5路母体腹壁混合心电信号的波形图;6 is a waveform diagram of a 5-channel mother abdominal wall mixed ECG signal used in a specific embodiment of the present invention;
图7为本发明的具体实施例的实验结果分析图。Figure 7 is an analysis of experimental results of a specific embodiment of the present invention.
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。The embodiments of the present invention will be further described in detail below with reference to the accompanying drawings. The illustrative embodiments of the present invention and the description thereof are intended to explain the present invention, but are not intended to limit the invention.
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention are clearly and completely described in the following with reference to the accompanying drawings in the embodiments of the present invention. It is obvious that the described embodiments are only a part of the embodiments of the present invention, but not all embodiments. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative efforts are within the scope of the present invention.
本领域技术技术人员知道,本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此,本公开可以具体实现为以下形式,即:完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。Those skilled in the art will appreciate that embodiments of the present invention can be implemented as a system, apparatus, device, method, or computer program product. Accordingly, the present disclosure may be embodied in the form of full hardware, complete software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.
下面参考本发明的若干代表性实施方式,详细阐释本发明的原理和精神。The principles and spirit of the present invention are explained in detail below with reference to a few representative embodiments of the invention.
本发明提出了一种基于平滑窗的自适应单导联胎儿心电分离方法及装置,对母体腹壁信号中的QRS波群间信号采用平滑窗进行平滑处理,以平滑掉QRS波群间信号段的胎儿心电节拍点,对腹壁心电信号中的母体QRS波群分割提取,对提取后的N个连续母体QRS波群以母体R波波峰位置为定位点,进行线性平均,生成平均QRS波群模板,截取当前输入的腹壁信号中的母体QRS波群和平均模板中的母体QRS波群,利用二者的相关系数,对平均模板中的母体QRS波群进行复制修正,以得到动态的母体腹壁心电估计信号。并且,通过自适应自主学习,对消掉腹壁信号中的母体分量,该方案避免了胸导联中母体分量和腹壁导联中母体分量非相关性带来的母体分量泄露影响,同时可以实现单导联胎儿心电分离。The invention provides an adaptive single-lead fetal electrocardiographic separation method and device based on smooth window, which uses a smoothing window to smooth the QRS inter-group signals in the maternal abdominal wall signal to smooth out the QRS inter-group signal segment. The fetal heartbeat beat point is extracted from the parent QRS complex in the abdominal wall ECG signal. The extracted N consecutive parent QRS complexes are linearly averaged by the parent R wave peak position, and the average QRS wave is generated. The group template intercepts the parent QRS complex in the input abdominal wall signal and the parent QRS complex in the average template, and uses the correlation coefficients of the two to copy and correct the parent QRS complex in the average template to obtain a dynamic matrix. Abdominal wall ECG estimation signal. Moreover, by adaptive autonomous learning, the maternal component in the abdominal wall signal is eliminated, and the scheme avoids the influence of the maternal component leakage caused by the maternal component of the chest lead and the non-correlation of the maternal component in the abdominal wall lead, and can realize the single Lead fetal ECG separation.
图1为本发明实施例的胎儿心电分离方法的系统结构示意图。如图1所示,y(n)表示采集到的母体腹壁心电信号,f(n)为腹壁信号中的胎儿心电信号,x(n)为腹壁信号中的母体心电信号,为自适应滤波器的期望信号,x'(n)为通过腹壁母体心电估计机制估计的腹壁母体心电,自适应滤波器对输入信号不断的进行“学习”和“跟踪”,使输出信号最大限度的逼近腹壁心电信号,从而使误差信号最小,即为分离出的胎儿心电信号。
FIG. 1 is a schematic structural diagram of a system of a fetal electrocardiographic separation method according to an embodiment of the present invention. As shown in Figure 1, y(n) represents the acquired abdominal wall ECG signal, f(n) is the fetal ECG signal in the abdominal wall signal, and x(n) is the maternal ECG signal in the abdominal wall signal. The desired signal of the adaptive filter, x'(n) is the abdominal wall of the abdominal wall estimated by the abdominal wall ECG estimation mechanism, and the adaptive filter continuously "learns" and "tracks" the input signal to make the output signal Maximum approximation of the abdominal wall ECG signal, thus making the error signal Minimal, This is the separated fetal ECG signal.
也就是说,基于自适应的胎儿心电分离算法,需要两路输入信号,即:一路母体腹壁心电信号(主输入信号),一路估计出的腹壁母体心电信号(参考输入信号),这相对于现有技术中的需要接入一路胸导联心电信号作为参考输入信号,实现了单导联。That is to say, based on the adaptive fetal ECG separation algorithm, two input signals are needed, namely: a parent abdominal wall ECG signal (main input signal), and an estimated abdominal wall ECG signal (reference input signal). Compared with the prior art, it is required to access a chest lead ECG signal as a reference input signal, and a single lead is realized.
本发明的核心是如何计算出腹壁母体心电估计信号,以作为自适应算法的参考输入信号。图2为本发明的腹壁母体心电估计机制的示意框图,如图2所示,将QRS波群信号和QRS波群间信号分开处理,采用基于平滑窗的方法进行母体QRS波群间信号进行平滑处理,消除母体QRS波群间信号段的胎儿心电节拍点,利用平均相关法对母体QRS波群进行幅值修正,从而估计出腹壁母体信号波形。The core of the invention is how to calculate the abdominal wall mother ECG estimation signal as a reference input signal to the adaptive algorithm. 2 is a schematic block diagram of the estimation mechanism of the abdominal wall mother body of the present invention. As shown in FIG. 2, the QRS complex signal and the QRS complex signal are separately processed, and the smoothing window based method is used to perform the signal between the parent QRS complexes. The smoothing process eliminates the fetal ECG beat point of the signal segment between the parent QRS complexes, and uses the average correlation method to correct the amplitude of the parent QRS complex to estimate the abdominal wall signal waveform.
图3为本发明实施例的胎儿心电分离方法的处理流程图。如图3所示,包括:FIG. 3 is a flowchart of processing of a fetal electrocardiographic separation method according to an embodiment of the present invention. As shown in Figure 3, it includes:
步骤S101,对接收到的母体腹壁心电信号进行QRS波群检测,确定QRS波群的起始位置和QRS波群间的起始位置;Step S101: performing QRS complex detection on the received abdominal wall ECG signal to determine a starting position between the QRS complex and a starting position between the QRS complexes;
步骤S102,利用平滑窗对所述QRS波群间的信号段进行平滑处理;Step S102: Smoothing a signal segment between the QRS groups by using a smoothing window;
步骤S103,对所述QRS波群进行线性平均,生成平均模板;Step S103, performing linear averaging on the QRS complex to generate an average template;
步骤S104,根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号;Step S104, according to the correlation between the parent QRS complex in the average template and the parent QRS complex corresponding to the current ECG cycle, correct the dynamic amplitude of the parent QRS complex in the average template to generate a matrix. Abdominal wall ECG estimation signal;
步骤S105,将所述母体腹壁心电估计信号作为参考信号,将所述母体腹壁心电信号作为期望信号,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号。Step S105, using the mother abdominal wall ECG estimation signal as a reference signal, using the mother abdominal wall ECG signal as a desired signal, and separating the fetal ECG signal from the mother abdominal wall ECG signal by an adaptive algorithm.
在步骤S102中,采用固定长度窗,对所述QRS波群间的信号段进行平滑处理,所述固定长度窗的取值根据信号采样率和胎儿心电波群长度来确定。In step S102, a signal segment between the QRS complexes is smoothed by using a fixed length window, and the value of the fixed length window is determined according to a signal sampling rate and a fetal ECG group length.
如图2所示,用固定长度窗M,对QRS波群间信号段进行平滑处理,固定窗长度约为胎儿心电QRS波长度,胎儿心电波群长度一般为47ms至85ms,故M取值可为(0.047*d,0.085*d),其中参数d为信号采样率。As shown in Fig. 2, the signal segment between QRS groups is smoothed by a fixed length window M. The length of the fixed window is about the length of the fetal ECG QRS wave, and the length of the fetal ECG group is generally 47ms to 85ms, so the value of M is It can be (0.047*d, 0.085*d), where the parameter d is the signal sampling rate.
QRS波群间信号可表示为:母体心电幅值一般是胎儿心电幅值的十几倍,经过平滑后,母体QRS波群间几乎不再含有胎儿QRS波群。The QRS inter-group signal can be expressed as: The maternal ECG amplitude is generally more than ten times the fetal ECG amplitude. After smoothing, the maternal QRS complex almost no longer contains the fetal QRS complex.
在步骤S103中,对所述QRS波群进行线性平均,生成平均片段,包括:从所述母体腹壁心电信号中,提取出母体心电R波波峰位置;以及以所述R波波峰位置为中心定位点,对所述QRS波群进行线性平均,生成所述平均模板。
In step S103, linearly averaging the QRS complex to generate an average segment, comprising: extracting a parent ECG R wave peak position from the parent abdominal wall ECG signal; and using the R wave peak position as A central anchor point that linearly averages the QRS complex to generate the average template.
具体实施时,需要确定母体腹壁心电信号中的R波波峰位置,该方法为:对接收到的母体腹壁心电信号进行差分运算,以获取腹壁心电信号的斜率曲线,对斜率曲线做绝对滑动积分运算,记录每秒钟内的最大值,对前8秒内的最大值做平均,则R波判定阈值=系数比值*前8秒内的最大值平均,对差分值大于R波判定阈值的段,进行R波波峰位置最大值查找,以确定最终的R波波峰位置点。In the specific implementation, it is necessary to determine the position of the R wave peak in the mother abdominal wall ECG signal. The method is: performing a differential operation on the received abdominal wall ECG signal to obtain the slope curve of the abdominal wall ECG signal, and making the slope curve absolutely Sliding integral operation, recording the maximum value in each second, averaging the maximum value in the first 8 seconds, then the R wave determination threshold = coefficient ratio * the maximum value in the first 8 seconds, the difference value is greater than the R wave determination threshold The segment is subjected to a maximum search of the R wave peak position to determine the final R wave peak position point.
然后以腹壁心电信号中R波波峰位置为中心定位点,截取母体QRS波群心电波形y1,y1=[y1(1),y1(2),...,y1(P)],其中P为母体QRS波群长度。Then take the position of the R wave peak in the abdominal wall ECG signal as the central positioning point, and intercept the parent QRS wave group ECG waveform y 1 , y 1 =[y 1 (1), y 1 (2),...,y 1 ( P)], where P is the parent QRS complex length.
对当前心电周期对应的心电信号片段以及前面M-1个连续的心电信号片段进行线性平均,得到平均模板为:
The linear average of the ECG signal segment corresponding to the current ECG cycle and the previous M-1 consecutive ECG signal segments is obtained, and the average template is:
在本发明实施例中,能否有效的抑制掉母体QRS波群分量,对正确识别分离结果中的胎儿心电分量起着决定性作用。平均模板中的母体心电波形,是对较早的心电波形和最新波形的等权累加和线性平均,估计出的母体心电片段不具有动态性。因为,本发明的步骤S104根据最新母体QRS波群和平均片段中母体QRS波群的相关性,对平均模板中的母体QRS波群进行动态幅值修正。In the embodiment of the present invention, whether the parent QRS complex component can be effectively suppressed can play a decisive role in correctly identifying the fetal ECG component in the separation result. The maternal ECG waveform in the average template is an equal weighted and linear average of the earlier ECG waveform and the latest waveform, and the estimated maternal ECG fragment is not dynamic. Because step S104 of the present invention performs dynamic amplitude correction on the parent QRS complex in the average template based on the correlation between the latest parent QRS complex and the parent QRS complex in the average segment.
具体实施时,方法如下:For specific implementation, the method is as follows:
首先,计算所述当前心电周期对应的QRS波群片段与所述平均模板中的QRS波群片段的相关性,生成相关系数,具体步骤为:First, calculating a correlation between the QRS complex segment corresponding to the current ECG cycle and the QRS complex segment in the average template, and generating a correlation coefficient, the specific steps are:
(1)提取当前心电周期中的母体QRS波群片段:y3=[y3(1),y3(2),...,y3(P)],其中P为母体QRS波群片段的长度;(1) Extract the parent QRS complex segment in the current ECG cycle: y 3 = [y 3 (1), y 3 (2), ..., y 3 (P)], where P is the parent QRS complex The length of the segment;
(2)提取出平均模板中母体QRS波群片段:y2=[y2(1),y2(2),...,y2(P)];(2) Extract the parent QRS complex segment in the average template: y 2 = [y 2 (1), y 2 (2), ..., y 2 (P)];
(3)计算量片段的相关系数r23:(3) Correlation coefficient r 23 of the calculated quantity segment:
其次,根据所述相关系数,用所述当前心电周期对应的QRS波群片段对所述平均模板中的母体QRS波群进行动态幅值修正,所采用的表达式如下:Secondly, according to the correlation coefficient, the QRS complex corresponding to the current ECG cycle is used to perform dynamic amplitude correction on the parent QRS complex in the average template, and the expression used is as follows:
由于平均段生成的平均模板具有当前腹壁心电周期内所对应母体心电波形的信号,并且基本消除掉了其中的胎儿心电波形,通过相关系数,对平均后的母体心电QRS波群段进行动态幅值修正,能最大限度的估计母体腹壁心电信号。Since the average template generated by the average segment has the signal of the corresponding parent ECG waveform in the current abdominal wall electrocardiogram, and the fetal electrocardiogram waveform is substantially eliminated, the averaged parent ECG QRS cluster segment is obtained by the correlation coefficient. Dynamic amplitude correction can be used to estimate the maternal abdominal wall ECG signal to the maximum extent.
在步骤S105中,腹部电极取得主输入信号,将估计出的母体腹壁心电信号作为输入信号,此信号中基本不含有胎儿心电信号,只包含欲消除的噪声和母体心电信号;腹壁心电信号为期望信号,通过自适应算法,消除了输入信号与期望信号中相关的部分,保留不相关的部分,因此腹壁心电信号中的胎儿心电信号分离可以被最大化的提取出来。In step S105, the abdominal electrode obtains the main input signal, and the estimated mother abdominal wall ECG signal is used as an input signal, and the signal substantially does not contain the fetal ECG signal, and only contains the noise to be eliminated and the mother ECG signal; The electrical signal is the desired signal. The adaptive algorithm eliminates the correlation between the input signal and the desired signal, and preserves the uncorrelated portion. Therefore, the fetal ECG signal separation in the abdominal wall ECG signal can be maximized.
本发明实施例中,自适应算法可以为自适应的最小均方(LMS)算法、自适应的递归最小二乘方(RLS)算法、最小二乘格形(LSL)自适应算法等。In the embodiment of the present invention, the adaptive algorithm may be an adaptive least mean square (LMS) algorithm, an adaptive recursive least squares (RLS) algorithm, a least squares lattice (LSL) adaptive algorithm, or the like.
应当注意,尽管在附图中以特定顺序描述了本发明方法的操作,但是,这并非要求或者暗示必须按照该特定顺序来执行这些操作,或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地,可以省略某些步骤,将多个步骤合并为一个步骤执行,和/或将一个步骤分解为多个步骤执行。It should be noted that although the operations of the method of the present invention are described in a particular order in the drawings, this is not required or implied that such operations must be performed in that particular order, or that all illustrated operations must be performed to achieve the desired results. . Additionally or alternatively, certain steps may be omitted, multiple steps being combined into one step, and/or one step being broken down into multiple steps.
在介绍了本发明示例性实施方式的方法之后,接下来,参考图6对本发明示例性实施方式的胎儿心电分离装置进行介绍。该装置的实施可以参见上述方法的实施,重复之处不再赘述。以下所使用的术语“模块”和“单元”,可以是实现预定功能的软件和/或硬件。尽管以下实施例所描述的模块较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。Having described the method of the exemplary embodiment of the present invention, next, a fetal electrocardiographic separation device according to an exemplary embodiment of the present invention will be described with reference to FIG. For the implementation of the device, reference may be made to the implementation of the above method, and the repeated description will not be repeated. The terms "module" and "unit" as used hereinafter may be software and/or hardware that implement a predetermined function. Although the modules described in the following embodiments are preferably implemented in software, hardware, or a combination of software and hardware, is also possible and contemplated.
图4为本发明实施例的胎儿心电分离装置的结构示意图。如图所示,包括:4 is a schematic structural view of a fetal electrocardiographic separation device according to an embodiment of the present invention. As shown, it includes:
波群检测模块101,用于对接收到的母体腹壁心电信号进行QRS波群检测,确定QRS波群的起始位置和QRS波群间的起始位置;平滑处理模块102,用于利用平滑窗对所述QRS波群间的信号段进行平滑处理;平均处理模块103,用于对所述QRS波群进行线性平均,生成平均模板;波群修正模块104,用于根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号;分离模块105,用于将所述母体腹壁心电估计信号作为参考信号,将所述母体腹壁心电信号作为期望信号,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号。The wave group detecting module 101 is configured to perform QRS complex detection on the received parent abdominal wall electrocardiographic signal to determine a starting position between the QRS complex and the QRS complex; and a smoothing module 102 for utilizing smoothing The window smoothes the signal segments between the QRS groups; the averaging processing module 103 is configured to linearly average the QRS groups to generate an average template; the group correction module 104 is configured to be used according to the average template Correlation of the parent QRS complex with the parent QRS complex corresponding to the current ECG cycle, correcting the dynamic amplitude of the parent QRS complex in the average template, and generating a maternal abdominal wall ECG estimation signal; the separation module 105 The parental abdominal wall ECG estimation signal is used as a reference signal, and the mother abdominal wall ECG signal is used as a desired signal, and the fetal ECG signal is separated from the mother abdominal wall ECG signal by an adaptive algorithm.
进一步地,在本发明实施例中,所述平滑处理模块102利用平滑窗对所述QRS波群间的信号段进行平滑处理,具体包括:采用固定长度窗,对所述QRS波群间的信号段进行平滑处理,所述固定长度窗的取值根据信号采样率和胎儿心电波群长度来确定。
Further, in the embodiment of the present invention, the smoothing processing module 102 performs smoothing on the signal segments between the QRS groups by using a smoothing window, and specifically includes: using a fixed length window, the signal between the QRS groups The segment is smoothed, and the value of the fixed length window is determined according to the signal sampling rate and the length of the fetal ECG group.
进一步地,在本发明实施例中,所述平均处理模块103对所述QRS波群进行线性平均,生成平均模板,包括:从所述母体腹壁心电信号中,提取出母体心电R波波峰位置;以所述R波波峰位置为中心定位点,对所述QRS波群进行线性平均,生成所述平均模板。Further, in the embodiment of the present invention, the averaging processing module 103 linearly averages the QRS complex to generate an average template, including: extracting a parent ECG R wave peak from the mother abdominal wall ECG signal Position: positioning the point of the R wave peak as a center, and linearly averaging the QRS group to generate the average template.
进一步地,在本发明实施例中,所述波群修正模块104根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号,如图5所示,具体包括:Further, in the embodiment of the present invention, the group correction module 104 compares the correlation between the parent QRS group in the average template and the parent QRS group corresponding to the current ECG cycle. In the parent QRS wave group, the dynamic amplitude correction is performed to generate the mother abdominal wall ECG estimation signal, as shown in FIG. 5, which specifically includes:
相关系数计算单元1041,用于计算所述当前心电周期对应的QRS波群片段与所述平均模板中的QRS波群片段的相关性,生成相关系数;The correlation coefficient calculation unit 1041 is configured to calculate a correlation between the QRS group segment corresponding to the current ECG cycle and the QRS group segment in the average template, and generate a correlation coefficient;
动态修正单元1042,用于根据所述相关系数,用所述当前心电周期对应的QRS波群片段对所述平均模板中的母体QRS波群进行动态幅值修正,生成母体腹壁心电估计信号。The dynamic correction unit 1042 is configured to perform dynamic amplitude correction on the parent QRS group in the average template by using the QRS complex corresponding to the current ECG cycle according to the correlation coefficient to generate a maternal abdominal wall ECG estimation signal. .
进一步地,在本发明实施例中,所述分离模块105通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号,所述自适应算法包括自适应的最小均方算法、自适应的递归最小二乘方算法、最小二乘格形自适应算法。Further, in the embodiment of the present invention, the separation module 105 separates a fetal ECG signal from the mother abdominal wall ECG signal by an adaptive algorithm, where the adaptive algorithm includes an adaptive least mean square algorithm, Adaptive recursive least squares algorithm, least squares lattice adaptive algorithm.
此外,尽管在上文详细描述中提及了胎儿心电分离装置的若干单元,但是这种划分仅仅并非强制性的。实际上,根据本发明的实施方式,上文描述的两个或更多单元的特征和功能可以在一个单元中具体化。同样,上文描述的一个单元的特征和功能也可以进一步划分为由多个单元来具体化。Moreover, although several elements of the fetal electrocardiographic separation device are mentioned in the detailed description above, such division is merely not mandatory. Indeed, in accordance with embodiments of the present invention, the features and functions of two or more units described above may be embodied in one unit. Also, the features and functions of one unit described above may be further divided into a plurality of units.
以下通过一个具体实施例来说明采用本发明的基于平滑窗的胎儿心电分离方法实现胎儿心电信号分离的效果。The effect of using the smooth window-based fetal electrocardiographic separation method of the present invention to achieve fetal ECG signal separation will be described below by way of a specific embodiment.
参看图6所示,为该具体实施例采用的5路母体腹壁混合心电信号的波形图。如图6所示,每路信号的数据长度为2500个样本点,采样率250Hz。由图6可知,第1路腹壁信号为胎母心电信噪比最大的一路信号,即此路信号效果最好,第4路和第5路腹壁心电信号用肉眼已经无法识别出胎儿心电R波峰位置,此两路信号提取效果已不能为本发明的算法所用,故不在本发明的算法的考虑范围之内。Referring to Figure 6, a waveform diagram of a mixed 5-channel maternal abdominal wall ECG signal used in this embodiment is shown. As shown in Figure 6, the data length of each signal is 2500 sample points, and the sampling rate is 250 Hz. It can be seen from Fig. 6 that the first road abdominal wall signal is the one signal with the largest telecommunication noise ratio of the fetal heart, that is, the signal of this road is the best, and the fourth and fifth abdominal wall ECG signals can not recognize the fetal heart with the naked eye. The position of the electric R peak, which is not used by the algorithm of the present invention, is not considered by the algorithm of the present invention.
本实施例选择第2路信号作为算法验证数据,自适应滤波采用自适应递归最小二乘方算法,其实验结果分析参看图7。In this embodiment, the second path signal is selected as the algorithm verification data, and the adaptive filtering adopts the adaptive recursive least squares algorithm, and the experimental result analysis is shown in FIG. 7.
1、胎儿QRS波群远离母体QRS波群:1. The fetal QRS complex is far from the parent QRS complex:
胎儿QRS波群远离母体QRS波群,如图7中标注a,由图7可知,这种情况下的胎
儿QRS波分量被提取出,左右的母体QRS波分量被抑制掉。此种胎儿QRS波群和母体QRS波群的位置关系,是腹壁心电中最普遍的位置关系,也是最容易提取出的位置关系。The fetal QRS complex is far from the parent QRS complex, as labeled a in Figure 7, as can be seen from Figure 7, the tire in this case
The QRS component is extracted and the left and right parent QRS components are suppressed. The positional relationship between the fetal QRS complex and the parent QRS complex is the most common positional relationship in the abdominal wall ECG and the most easily extracted positional relationship.
2、胎儿QRS波群紧邻母体QRS波群:2. The fetal QRS complex is adjacent to the parent QRS complex:
胎儿QRS波群紧邻母体QRS波群,如图7中标注b所示,由图7可知,这种情况下的胎儿QRS波分量被提取出,紧邻的母体QRS波分量被抑制掉,并未对胎儿心电QRS波群的识别造成影响。此种胎儿QRS波群和母体QRS波群的位置关系,是腹壁心电中少数比例存在的位置关系,也是较难提取出的位置关系。The fetal QRS complex is adjacent to the parent QRS complex, as indicated by the mark b in Fig. 7. It can be seen from Fig. 7 that the fetal QRS component is extracted in this case, and the adjacent QRS component is suppressed. The identification of fetal ECG QRS complexes has an impact. The positional relationship between the fetal QRS complex and the parent QRS complex is a small number of positional relationships in the abdominal wall ECG, and it is also difficult to extract the positional relationship.
3、胎儿QRS波群和母体QRS波群完全重合:3. The fetal QRS complex and the parent QRS complex are completely coincident:
胎儿QRS波群完全和母体QRS波群重合,如图7中标注c所示,由图7可知,这种情况下的胎儿QRS波分量被提取出,重合的母体QRS波分量被抑制掉。此种胎儿QRS波群和母体QRS波群的位置关系,是腹壁心电中少数比例存在的位置关系,也是最难提取出的位置关系。The fetal QRS complex is completely coincident with the parent QRS complex, as indicated by the reference c in Fig. 7, and it can be seen from Fig. 7 that the fetal QRS component is extracted in this case, and the coincident parent QRS component is suppressed. The positional relationship between the fetal QRS complex and the parent QRS complex is a small number of positional relationships in the abdominal wall ECG and the most difficult to extract positional relationship.
由以上结果分析可知,通过本发明的基于平滑窗的胎儿心电分离方法,胎儿QRS波群远离母体QRS波群、胎儿QRS波群紧邻母体QRS波群、胎儿QRS波群和母体QRS波群完全重合三种情况,均有非常理想的提取效果。It can be seen from the above result analysis that, by the smooth window-based fetal electrocardiographic separation method of the present invention, the fetal QRS complex is far from the parent QRS complex, the fetal QRS complex is adjacent to the parent QRS complex, the fetal QRS complex, and the parent QRS complex are completely Coincident three cases, have a very good extraction effect.
本发明提出的胎儿心电分离方法及装置,具有以下有益效果:The fetal electrocardiographic separation method and device proposed by the invention have the following beneficial effects:
1、单导联:算法仅需要一路母体腹壁心电信号,即可完成基于自适应的胎儿信号分离。传统的采用基于自适应的胎儿心电信号分离算法,需要两路信号,即:一路腹壁心电信号(主输入信号),一路胸导联心电信号(参考输入信号),通过对主输入信号进行母体QRS波群间平滑处理,消除母体QRS波群间信号段的胎儿分量,对母体QRS波群进行幅值修正,从而估计出腹壁母体信号波形,信号被输入自适应滤波器,同时未处理的同步腹壁心电信号作为自适应滤波的期望信号,因此,本发明只需要一路腹壁心电信号,这将大大简化了系统的导联系统。1. Single lead: The algorithm only needs one way of the abdominal wall ECG signal to complete the adaptive fetal signal separation. Traditionally, based on adaptive fetal ECG signal separation algorithm, two signals are needed, namely: one abdominal wall ECG signal (main input signal), one chest lead ECG signal (reference input signal), through the main input signal The maternal QRS group smoothing process is performed to eliminate the fetal component of the signal segment between the parent QRS groups, and the amplitude correction of the parent QRS group is performed to estimate the abdominal wall signal waveform, and the signal is input into the adaptive filter while being unprocessed. The synchronized abdominal wall ECG signal serves as the desired signal for adaptive filtering. Therefore, the present invention requires only one abdominal wall ECG signal, which greatly simplifies the system lead system.
2、分离效果好:通过平滑和母体QRS波波群幅值修正处理,估计母体腹壁心电信号,此方法避免了胸导联和腹壁导联中母体分量非相关所带来的母体心电抑制不彻底的影响,具有更好的分离效果。2, the separation effect is good: through the smooth and maternal QRS wave group amplitude correction processing, estimate the maternal abdominal wall ECG signal, this method avoids maternal electrocardiogram inhibition caused by the uncorrelated maternal component in the chest lead and abdominal wall lead Incomplete effects have a better separation effect.
3、算法设计思路简单直观,针对于异常的临床数据,易于查找原因和添加修正模块。
3, the algorithm design idea is simple and intuitive, for abnormal clinical data, easy to find the cause and add correction module.
本发明实施例还提供一种计算机可读程序,其中当执行所述程序时,所述程序使得计算机执行本发明实施例所述的胎儿心电分离方法。The embodiment of the present invention further provides a computer readable program, wherein when the program is executed, the program causes a computer to execute the fetal electrocardiographic separation method according to the embodiment of the present invention.
本发明实施例还提供一种存储有计算机可读程序的存储介质,其中所述计算机可读程序使得计算机执行本发明实施例所述的胎儿心电分离方法。The embodiment of the present invention further provides a storage medium storing a computer readable program, wherein the computer readable program causes the computer to execute the fetal electrocardiographic separation method according to the embodiment of the present invention.
本发明以上的装置和方法可以由硬件实现,也可以由硬件结合软件实现。本发明涉及这样的计算机可读程序,当该程序被逻辑部件所执行时,能够使该逻辑部件实现上文所述的装置或构成部件,或使该逻辑部件实现上文所述的各种方法或步骤。本发明还涉及用于存储以上程序的存储介质,如硬盘、磁盘、光盘、DVD、flash存储器等。The above apparatus and method of the present invention may be implemented by hardware or by hardware in combination with software. The present invention relates to a computer readable program that, when executed by a logic component, enables the logic component to implement the apparatus or components described above, or to cause the logic component to implement the various methods described above Or steps. The present invention also relates to a storage medium for storing the above program, such as a hard disk, a magnetic disk, an optical disk, a DVD, a flash memory, or the like.
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art will appreciate that embodiments of the present invention can be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or a combination of software and hardware. Moreover, the invention can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) including computer usable program code.
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present invention has been described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (system), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or FIG. These computer program instructions can be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing device to produce a machine for the execution of instructions for execution by a processor of a computer or other programmable data processing device. Means for implementing the functions specified in one or more of the flow or in a block or blocks of the flow chart.
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。The computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing device to operate in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture comprising the instruction device. The apparatus implements the functions specified in one or more blocks of a flow or a flow and/or block diagram of the flowchart.
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions can also be loaded onto a computer or other programmable data processing device such that a series of operational steps are performed on a computer or other programmable device to produce computer-implemented processing for execution on a computer or other programmable device. The instructions provide steps for implementing the functions specified in one or more of the flow or in a block or blocks of a flow diagram.
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人
员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
The principles and embodiments of the present invention have been described in connection with the specific embodiments of the present invention. The description of the above embodiments is only for the understanding of the method of the present invention and its core ideas;
In the above, the description of the present invention is not limited to the scope of the present invention.
Claims (10)
- 一种胎儿心电分离方法,其特征在于,包括:A method for separating fetal ECG, characterized in that it comprises:对接收到的母体腹壁心电信号进行QRS波群检测,确定QRS波群的起始位置和QRS波群间的起始位置;Performing QRS complex detection on the received abdominal wall ECG signal to determine the starting position of the QRS complex and the starting position between the QRS complexes;利用平滑窗对所述QRS波群间的信号段进行平滑处理;Smoothing the signal segments between the QRS complexes using a smoothing window;对所述QRS波群进行线性平均,生成平均模板;Performing a linear average on the QRS complex to generate an average template;根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号;And according to the correlation between the parent QRS complex in the average template and the parent QRS complex corresponding to the current ECG cycle, correcting the dynamic amplitude of the parent QRS complex in the average template to generate the maternal abdominal wall electrocardiogram Estimated signal;将所述母体腹壁心电估计信号作为参考信号,将所述母体腹壁心电信号作为期望信号,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号。The maternal abdominal wall ECG estimation signal is used as a reference signal, and the maternal abdominal wall electrocardiographic signal is used as a desired signal, and the fetal ECG signal is separated from the maternal abdominal wall electrocardiographic signal by an adaptive algorithm.
- 根据权利要求1所述的胎儿心电分离方法,其特征在于,利用平滑窗对所述QRS波群间的信号段进行平滑处理,包括:The fetal electrocardiographic separation method according to claim 1, wherein the smoothing of the signal segments between the QRS complexes by using a smoothing window comprises:采用固定长度窗,对所述QRS波群间的信号段进行平滑处理,所述固定长度窗的取值根据信号采样率和胎儿心电波群长度来确定。The signal segment between the QRS complexes is smoothed by a fixed length window, and the value of the fixed length window is determined according to a signal sampling rate and a fetal ECG group length.
- 根据权利要求1所述的胎儿心电分离方法,其特征在于,对所述QRS波群进行线性平均,生成平均模板,包括:The fetal electrocardiographic separation method according to claim 1, wherein linearly averaging the QRS complex to generate an average template comprises:从所述母体腹壁心电信号中,提取出母体心电R波波峰位置;Extracting the peak position of the parent ECG R wave from the maternal abdominal wall electrocardiographic signal;以所述R波波峰位置为中心定位点,对所述QRS波群进行线性平均,生成所述平均模板。The QRS complex is linearly averaged with the R wave peak position as a central positioning point to generate the average template.
- 根据权利要求3所述的胎儿心电分离方法,其特征在于,根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号,包括:The fetal electrocardiographic separation method according to claim 3, wherein said average template is based on a correlation between a parent QRS complex in said average template and a parent QRS complex corresponding to said current ECG cycle In the parent QRS complex, the dynamic amplitude correction is performed to generate the maternal abdominal wall ECG estimation signal, including:计算所述当前心电周期对应的QRS波群片段与所述平均模板中的QRS波群片段的相关性,生成相关系数;Calculating a correlation between the QRS complex segment corresponding to the current ECG cycle and the QRS complex segment in the average template to generate a correlation coefficient;根据所述相关系数,用所述当前心电周期对应的QRS波群片段对所述平均模板中的母体QRS波群进行动态幅值修正,生成母体腹壁心电估计信号。And according to the correlation coefficient, the QRS complex corresponding to the current ECG cycle is used to perform dynamic amplitude correction on the parent QRS complex in the average template to generate a maternal abdominal wall ECG estimation signal.
- 根据权利要求1所述的胎儿心电分离方法,其特征在于,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号,包括: The fetal electrocardiographic separation method according to claim 1, wherein the fetal ECG signal is separated from the maternal abdominal wall electrocardiographic signal by an adaptive algorithm, including:所述自适应算法包括自适应的最小均方算法、自适应的递归最小二乘方算法、最小二乘格形自适应算法。The adaptive algorithm includes an adaptive least mean square algorithm, an adaptive recursive least squares algorithm, and a least squares lattice adaptive algorithm.
- 一种胎儿心电分离装置,其特征在于,包括:A fetal electrocardiographic separation device, comprising:波群检测模块,用于对接收到的母体腹壁心电信号进行QRS波群检测,确定QRS波群的起始位置和QRS波群间的起始位置;The wave group detecting module is configured to perform QRS complex detection on the received abdominal wall ECG signal to determine the starting position of the QRS group and the starting position between the QRS groups;平滑处理模块,用于利用平滑窗对所述QRS波群间的信号段进行平滑处理;a smoothing processing module, configured to smooth a signal segment between the QRS groups by using a smoothing window;平均处理模块,用于对所述QRS波群进行线性平均,生成平均模板;An averaging processing module, configured to perform linear averaging on the QRS complex to generate an average template;波群修正模块,用于根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号;a group correction module, configured to: according to a correlation between a parent QRS group in the average template and a parent QRS group corresponding to the current ECG period, a dynamic amplitude of a parent QRS group in the average template Correction, generating a maternal abdominal wall ECG estimation signal;分离模块,用于将所述母体腹壁心电估计信号作为参考信号,将所述母体腹壁心电信号作为期望信号,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号。a separation module, configured to use the mother abdominal wall ECG estimation signal as a reference signal, and use the mother abdominal wall ECG signal as a desired signal to separate a fetal ECG signal from the mother abdominal wall ECG signal through an adaptive algorithm .
- 根据权利要求6所述的胎儿心电分离装置,其特征在于,所述平滑处理模块利用平滑窗对所述QRS波群间的信号段进行平滑处理,具体包括:The fetal electrocardiographic apparatus according to claim 6, wherein the smoothing processing module smoothes the signal segments between the QRS complexes by using a smoothing window, and specifically includes:采用固定长度窗,对所述QRS波群间的信号段进行平滑处理,所述固定长度窗的取值根据信号采样率和胎儿心电波群长度来确定。The signal segment between the QRS complexes is smoothed by a fixed length window, and the value of the fixed length window is determined according to a signal sampling rate and a fetal ECG group length.
- 根据权利要求6所述的胎儿心电分离装置,其特征在于,所述平均处理模块对所述QRS波群进行线性平均,生成平均模板,包括:The fetal electrocardiographic apparatus according to claim 6, wherein the averaging processing module linearly averages the QRS complex to generate an average template, including:从所述母体腹壁心电信号中,提取出母体心电R波波峰位置;Extracting the peak position of the parent ECG R wave from the maternal abdominal wall electrocardiographic signal;以所述R波波峰位置为中心定位点,对所述QRS波群进行线性平均,生成所述平均模板。The QRS complex is linearly averaged with the R wave peak position as a central positioning point to generate the average template.
- 根据权利要求8所述的胎儿心电分离装置,其特征在于,所述波群修正模块根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号,包括:The fetal electrocardiographic apparatus according to claim 8, wherein said group correction module is based on a correlation between a parent QRS complex in said average template and a parent QRS complex corresponding to said current ECG cycle And correcting the dynamic amplitude of the parent QRS wave group in the average template to generate a maternal abdominal wall ECG estimation signal, including:相关系数计算单元,用于计算所述当前心电周期对应的QRS波群片段与所述平均模板中的QRS波群片段的相关性,生成相关系数;a correlation coefficient calculation unit, configured to calculate a correlation between the QRS complex segment corresponding to the current ECG cycle and a QRS complex segment in the average template, to generate a correlation coefficient;动态修正单元,用于根据所述相关系数,用所述当前心电周期对应的QRS波群片段对所述平均模板中的母体QRS波群进行动态幅值修正,生成母体腹壁心电估计信号。 The dynamic correction unit is configured to perform dynamic amplitude correction on the parent QRS group in the average template by using the QRS complex corresponding to the current ECG cycle according to the correlation coefficient to generate a maternal abdominal wall ECG estimation signal.
- 根据权利要求6所述的胎儿心电分离装置,其特征在于,所述分离模块通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号,包括:The fetal electrocardiographic apparatus according to claim 6, wherein the separation module separates the fetal electrocardiographic signal from the maternal abdominal wall electrocardiographic signal by an adaptive algorithm, including:所述自适应算法包括自适应的最小均方算法、自适应的递归最小二乘方算法、最小二乘格形自适应算法。 The adaptive algorithm includes an adaptive least mean square algorithm, an adaptive recursive least squares algorithm, and a least squares lattice adaptive algorithm.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2016/093525 WO2018023697A1 (en) | 2016-08-05 | 2016-08-05 | Fetal electrocardiosignal separation method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2016/093525 WO2018023697A1 (en) | 2016-08-05 | 2016-08-05 | Fetal electrocardiosignal separation method and device |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2018023697A1 true WO2018023697A1 (en) | 2018-02-08 |
Family
ID=61073325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2016/093525 WO2018023697A1 (en) | 2016-08-05 | 2016-08-05 | Fetal electrocardiosignal separation method and device |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2018023697A1 (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110604564A (en) * | 2019-09-16 | 2019-12-24 | 深圳市理邦精密仪器股份有限公司 | Electrode system, fetal electrocardiosignal acquisition and measurement method and device |
CN110613444A (en) * | 2019-09-16 | 2019-12-27 | 深圳市理邦精密仪器股份有限公司 | Method and device for displaying mother-fetus mixed electrocardiosignals |
CN111772627A (en) * | 2019-04-04 | 2020-10-16 | 中山大学 | An online fetal ECG signal extraction device and method |
CN112120688A (en) * | 2019-06-25 | 2020-12-25 | 深圳市理邦精密仪器股份有限公司 | Electrocardiosignal processing method, electrocardiosignal processing equipment and computer-readable storage medium |
CN114767126A (en) * | 2022-04-29 | 2022-07-22 | 广东粤港澳大湾区国家纳米科技创新研究院 | Fetal electrocardiogram monitoring method and system based on marked interference signals |
CN116269430A (en) * | 2023-03-01 | 2023-06-23 | 华南师范大学 | Electrocardiosignal R peak position extraction method, module, monitoring device and storage medium |
CN117617984A (en) * | 2022-08-09 | 2024-03-01 | 荣耀终端有限公司 | Electrocardiosignal processing method and electronic equipment |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5908393A (en) * | 1998-05-04 | 1999-06-01 | Cambridge Heart, Inc. | Reducing noise in a biological signal |
US20050267376A1 (en) * | 2004-05-28 | 2005-12-01 | Dorothee Marossero | Maternal-fetal monitoring system |
EP2016894A1 (en) * | 2007-07-20 | 2009-01-21 | Stichting Voor De Technische Wetenschappen | ECG signal processing |
CN105411577A (en) * | 2015-12-30 | 2016-03-23 | 深圳先进技术研究院 | Method and system for separating fetal ECG (electrocardiogram) |
CN105640545A (en) * | 2015-12-31 | 2016-06-08 | 深圳先进技术研究院 | Fetal electrocardiosignal extraction method and device |
-
2016
- 2016-08-05 WO PCT/CN2016/093525 patent/WO2018023697A1/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5908393A (en) * | 1998-05-04 | 1999-06-01 | Cambridge Heart, Inc. | Reducing noise in a biological signal |
US20050267376A1 (en) * | 2004-05-28 | 2005-12-01 | Dorothee Marossero | Maternal-fetal monitoring system |
EP2016894A1 (en) * | 2007-07-20 | 2009-01-21 | Stichting Voor De Technische Wetenschappen | ECG signal processing |
CN105411577A (en) * | 2015-12-30 | 2016-03-23 | 深圳先进技术研究院 | Method and system for separating fetal ECG (electrocardiogram) |
CN105640545A (en) * | 2015-12-31 | 2016-06-08 | 深圳先进技术研究院 | Fetal electrocardiosignal extraction method and device |
Non-Patent Citations (1)
Title |
---|
ZHOU, LIGAO ET AL.: "The Development of Fetal ECG Monitor Using Match Filtering", CHINESE JOURNAL OF MEDICAL INSTRUMENTATION, vol. 13, no. 5, 30 September 1989 (1989-09-30) * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111772627A (en) * | 2019-04-04 | 2020-10-16 | 中山大学 | An online fetal ECG signal extraction device and method |
CN111772627B (en) * | 2019-04-04 | 2023-08-15 | 中山大学 | Online fetal electrocardiosignal extraction device and method |
CN112120688A (en) * | 2019-06-25 | 2020-12-25 | 深圳市理邦精密仪器股份有限公司 | Electrocardiosignal processing method, electrocardiosignal processing equipment and computer-readable storage medium |
CN110604564A (en) * | 2019-09-16 | 2019-12-24 | 深圳市理邦精密仪器股份有限公司 | Electrode system, fetal electrocardiosignal acquisition and measurement method and device |
CN110613444A (en) * | 2019-09-16 | 2019-12-27 | 深圳市理邦精密仪器股份有限公司 | Method and device for displaying mother-fetus mixed electrocardiosignals |
CN110613444B (en) * | 2019-09-16 | 2022-06-07 | 深圳市理邦精密仪器股份有限公司 | Method and device for displaying mother-fetus mixed electrocardiosignals |
CN114767126A (en) * | 2022-04-29 | 2022-07-22 | 广东粤港澳大湾区国家纳米科技创新研究院 | Fetal electrocardiogram monitoring method and system based on marked interference signals |
CN117617984A (en) * | 2022-08-09 | 2024-03-01 | 荣耀终端有限公司 | Electrocardiosignal processing method and electronic equipment |
CN116269430A (en) * | 2023-03-01 | 2023-06-23 | 华南师范大学 | Electrocardiosignal R peak position extraction method, module, monitoring device and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2018023697A1 (en) | Fetal electrocardiosignal separation method and device | |
He et al. | Automatic detection of QRS complexes using dual channels based on U-Net and bidirectional long short-term memory | |
Liu et al. | A multi-step method with signal quality assessment and fine-tuning procedure to locate maternal and fetal QRS complexes from abdominal ECG recordings | |
CN112826513B (en) | A fetal heart rate detection system on FECG based on deep learning and specific correction | |
US20200337582A1 (en) | Method for Accurately Extracting Abnormal Potential within QRS | |
WO2006081447A2 (en) | Blind adaptive filter extraction of fetal electrocardiogram signal estimate | |
Abel et al. | A comprehensive survey on signal processing and machine learning techniques for non-invasive fetal ECG extraction | |
CN104305992B (en) | A kind of interactive fast automatic extracting method of Fetal ECG | |
CN101257843A (en) | Nonlinear Signal Separation Methods Using Nonlinear State-Space Projection Methods | |
Sedighian et al. | Pediatric heart sound segmentation using Hidden Markov Model | |
CN108272451A (en) | A kind of QRS wave recognition methods based on improvement wavelet transformation | |
CN103054572B (en) | Blind source separation method on basis of late potential signals of single-lead cardiac ventricle | |
CN104644160A (en) | Electrocardiogram pseudo-difference signal identification method and device | |
CN107622259A (en) | A T wave detection method, electrocardiographic data analysis method and device | |
Fotiadou et al. | Deep convolutional long short-term memory network for fetal heart rate extraction | |
CN107684423A (en) | A kind of Fetal ECG separation method and device | |
Wang et al. | Fetal ECG signal extraction from long-term abdominal recordings based on adaptive QRS removal and joint blind source separation | |
WO2018023698A1 (en) | Fetal-electrocardiogram separation method and device | |
CN109009088B (en) | TWA multichannel fusion estimation method based on tensor decomposition | |
US11622710B2 (en) | Efficient fetal-maternal ECG signal separation from two maternal abdominal leads via diffusion-based channel selection | |
Dessì et al. | Identification of fetal QRS complexes in low density non-invasive biopotential recordings | |
Katebi et al. | Unsupervised hidden semi-Markov model for automatic beat onset detection in 1D Doppler ultrasound | |
CN107684422A (en) | A kind of Fetal ECG separation method and device | |
Souriau et al. | Fetal heart rate monitoring by fusion of estimations from two modalities: A modified Viterbi’s algorithm | |
Golpaygani et al. | Detection and identification of S1 and S2 heart sounds using wavelet decomposition method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 16911292 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
32PN | Ep: public notification in the ep bulletin as address of the adressee cannot be established |
Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205 DATED 04.06.2019) |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 16911292 Country of ref document: EP Kind code of ref document: A1 |