Invariant heart beat span versus variant heart beat intervals and its application to fetal ECG extraction
© Yan et al.; licensee BioMed Central. 2014
Received: 16 September 2014
Accepted: 5 December 2014
Published: 12 December 2014
The fundamental assumptions for various kinds of fetal electrocardiogram (fECG) extraction methods are not consistent with each other, which is a very important problem needed to be ascertained.
Based on two public databases, the regularity on ECG wave durations for normal sinus rhythm is investigated statistically. Taking the ascertained regularity as an assumption, a new fECG extraction algorithm is proposed, called Partial R-R interval Resampling (PRR).
Both synthetic and real abdominal ECG signals are used to test the algorithm. The results indicate that the PRR algorithm has better performance over the whole R-R interval resampling based comb filtering method (RR) and linear template method (LP), which takes advantages of both LP and RR.
The final drawn conclusion is: (1) the proposition should be true that the individual’s heart beat span is invariable for normal sinus rhythm; (2) the proposed PRR fetal ECG extraction algorithm can estimate the maternal ECG (mECG) more accurately and stably even in the condition of large HRV, finally resulting in better fetal ECG extraction.
KeywordsFetal ECG extraction HRV Template Comb filter
Non-invasive fetal electrocardiogram (fECG) extraction is a classic dilemma in biomedical signal processing. Researchers put forward various kinds of methods to solve this problem. These methods can be classified into: (1) multi-channel methods, like principal component analysis and its variants [1–3], independent component analysis [4, 5], periodic component analysis (piCA) [6–9], some wavelet based methods [10, 11], adaptive method [12–14] and so on; (2) single-channel methods, such as singular value decomposition (SVD) , R-R interval resampling based comb filter (RR) , linear template methods (LP) [17–21] and so on; (3) and fusion methods [22–25]. A detailed description about these methods can be found in literature .
Among the above three kinds of methods, the single-channel methods are always used to subtract the maternal electrocardiogram (mECG) component from the mixed abdominal ECG (aECG) recording. Since mECG component usually takes up the highest proportion in aECG, the accurate estimation of mECG component is a key procedure to obtain final high-quality fECG.
When using SVD  and RR  methods to estimate mECG component, R-R interval resampling is required because of the periodicity and the heart rate variability (HRV) of mECG signal. In essence, an underlying hypothesis is taken that the lengths of various waves in each cardiac cycle are proportional to the length of the R-R interval (heart beat interval).
Linear template (LP) [17–21] methods are simple, effective and widely used mECG estimation methods. In LP methods, the wave complexes in current cardiac cycle are estimated by summing the weighted wave complexes from neighbour cardiac cycles, depending on at least the peak detection of each QRS complex for different heart cycles. In essence, it is assumed that when heart beat interval is time-variant, the heart beat span (i.e. the length from the beginning of P wave to the end of T wave) is stable. So no resampling step is taken in current LP methods.
For the above two assumptions, there is no clear conclusion on which one is correct. In [27, 28], it is said that the length from the beginning of P wave to the peak of R wave (Ps-R) is weakly dependent on heart rate, while in literatures [21, 29] the authors think it depends on the instantaneous heart rate. As for the length from the peak of R wave to the end of the T wave (R-Te), literatures [27, 28] have the opinion that the T wave is strongly dependent on heart rate, becoming narrower and closer to the QRS complex at rapid heart rates, while the opposite view is given in . In , it is considered that the length from the end of T wave to the beginning of next P wave (Te-Ps) is strongly dependent on the heartbeat duration. So far, we have not seen literatures to ascertain such disagreement problem for normal sinus rhythm.
Obviously, the above disagreement problem is related to whether SVD, RR and LP methods are used correctly. It is a very crucial thing for designing a high-quality fECG extraction method, especially when the mECG has large HRV.
Based on MIT Normal Sinus Rhythm Database (NSRDB as Set A)  and the MIT-BIH Arrhythmia Database (ADB as Set B) [30, 31], the problem is investigated in this paper and a regularity comes out that heart beat span is invariant for variable heart beat intervals. Then we propose a new corresponding single channel fECG extraction algorithm based on the regularity. Both HRV and periodicity are considered in the method. And its effectiveness is verified by synthetic and real recordings.
Investigation on regularity of ECG waves duration
The databases selected as research resources are Set A  and Set B [30, 31]. Set A is a normal sinus rhythm database, including 18 long-time records, with two channels for each record, and the sampling frequency is 128 Hz. And Set B is a database containing 48 two-channel recordings with different degree of arrhythmia; the sampling frequency is 360 Hz.
The above result on record “19093” presents that the t Ps-R and t R-Te are quite stable compared with t Ps-Ps . That is, there exists so called invariant heart beat span with the variant heart beat intervals. A further validation with the databases is given below.
First, because of the difficulty in detecting Ps and Te points accurately, for convenience we replace Ps-R by P-R interval, R-Te by R-T interval, Te-Ps by T-P interval and Ps-Ps by R-R interval for statistics’ validation. If the above conclusion from data “19093” is correct, the both substitute intervals P-R and R-T also should be stable, and T-P interval and R-R interval should always vary accordingly. Then, twelve records from Set A and seven from Set B with obvious P and T waves are selected. The first 500 seconds (s) data for each selected record are cut out for statistical analysis. Note that pieces that P or T wave is not clear are thrown away.
t P-R , t R-T , t T-P and t R-R for the nineteen records in Set A and Set B
Partial R-R interval resampling based fECG extraction Algorithm (PRR)
The proposed core steps of PRR are further described as follows:
Step 1 in Figure 5: segmentation. Detect the maternal R peaks in aECG. Take the Ps point as the point T Ps-R s before each R peak and Te point as the point T R-Te s after each R peak. Then a complete cardiac cycle can be divided as Ps-Te-Ps. Ps-Te period is the previously mentioned heart beat span and Te-Ps period is the diastolic period.
Step 2 in Figure 5: Te-Ps interval up-resample. Denote the lengths of diastolic period as l 1p , l 2p ,…, l np-1 , l np (points). Each diastolic period is resampled to have the same number of samples l vp . l vp can be the maximum one of l 1p , l 2p ,…, l np-1 , l np . Thus every diastolic period is up resampled. And every heart beat span keeps unchanged. The total length of a cardiac cycle after resampling is L = (T Ps-R + T R-Te )*fs + l vp , here fs is the sampling rate .
Step 3 in Figure 5: comb filtering. The up-resampled periodic signal convolutes with a comb filter. The designing detail for a comb filter can be found in literature . The result is an estimation of up-resampled mECG.Step 4 in Figure 5: Te-Ps interval down-resample. Resample diastolic periods of the estimated up-resampled mECG to recover their original lengths. It turns out to be the estimation to mECG. The residual ECG (rECG) is obtained after subtracting it from the original aECG. The rECG can be seen as the estimation to fECG; of course, the noise can be weakened in a further step to get a clearer fECG which is not discussed in the paper.
In above Step 1, the T Ps-R and T R-Te values have slight difference for different individuals. T Ps-R lies between 0.15 and 0.20 s, while T R-Te lies 0.20 to 0.35 s . In reality, it is difficult to detect the points Ps and Te steadily online, since P wave and T wave may be invisible in aECG. Therefore, we set T Ps-R =0.2 s and T R-Te = 0.4 s fixedly in the proposed PRR. A further illustration about this can be found in the Discussion Section.
The performance of the PRR algorithm is tested on synthetic mECG and aECG data and real antenatal abdominal recordings. The RR and LP are also implemented for comparison. In LP, the template window is chosen as the window between R peak - 5/12 T and R peak + 7/12 T, here T is the average length of several R-R intervals [17, 18].
Tests on synthetic data
For the simulated synthetic data showed in Figure 6(a), its ECG component is estimated by PRR, RR and LP respectively. The estimation errors are calculated and showed in Figure 6(b-d) respectively, in which error = eECG-ECG, with eECG representing the estimation of ECG. The root mean squares (RMS) of the estimation errors are also calculated, which are 0.0159, 0.0234 and 0.0419 respectively.
WPRs using PRR, RR and LP for 20 synthetic aECGs
PRR × 10-2
RR × 10-2
LP × 10-2
PRR × 10-2
RR × 10-2
LP × 10-2
From Figure 6(b-d) and Figure 8(b-d), it’s quite clear that: (1) RR method always has quite serious erroneous estimations on P and T waves; (2) LP method may have some transient errors at the edges of a diastolic period; (3) for PRR and LP, their estimations to the QRS complexes are almost the same. From Figure 7 it’s showed that: when r = 0, i.e. ECG has no HRV, the RMSs of errors for synthetic ECGs are the same for the three method; for others, RMSs for RR and LP increase with r but RMSs for PRR stay the same. Generally speaking, the tests on synthetic ECG show that PRR can estimate the ECG signal more accurately especially when it has large HRV.
Experiments on real data
Tests are taken on data from MIT ADFED [30, 31] which contains five real-life records. Each record has five channels - one direct fECG and four abdominal ECG. The sampling rate is 1 KHz and their sampling durations are five minutes.
WPR s for rECGs of record “r01” (*indicate the signal piece showed in Figure 9 )
On invariant heart beat span
Based on the lengths of ECG waves for recordings in Set A and Set B, we get a conclusion that for normal sinus rhythm when the instant heart beat interval changes the heart beat span keeps almost unchanged, while the length of diastolic period fluctuates with heart beat interval.
HRV is always explained by the autonomic nerve’s regulation of heart and circulation system . The activity of parasympathetic or sympathetic nervous system affects the sinus node’s beating rhythm [26, 27], and further affect the heart beat to have HRV. However, autonomic nerve may affect HRV just by determining the beat beginning point, i.e., Ps point on ECG, don't control the duration of heart being excited, i.e., heart beat span. The heart beat span for a person should only be determined by his or her own heart's physiology and healthy condition itself. Therefore, if there is no conduction block disease, the HRV for normal sinus rhythm just affects the variation of the length of Te-Ps interval. The lengths of P wave, P-R interval, QRS wave and S-T interval are independent of the heart beat rate , which constituting the stable heart beat span.
Individual difference’s effect on PRR performance
In Method Section we set T Ps-R and T R-Te to be fixed. But actually different individuals may have slight difference on them. Now we discuss the effect below.
The real lengths of Ps-R and R-Te intervals and their offsets to the fixed values for the 11 synthetic signals
t Ps-R (s)
t R-Te (s)
(t Ps-R -T Ps-R ) + (t R-Te -T R-Te )
With statistical analysis to the lengths of ECG waves from two databases, the conclusion comes out: heart beat span is invariant while the heart beat interval changes all the time. According to this prior information, a new modified fECG extraction algorithm is proposed called PRR. It copes with the mECG estimation inaccuracy caused by HRV and enhances the signal to noise ratio of the remaining fECG signal. Experiments show that the larger the HRV, the more obvious the advantage of PRR over other methods. So the proposed PRR method is outstanding on robustness, and it is believed it can be used into the practical fECG extraction. The computation complexity of this method can be studied in further work.
This paper is supported by the National Natural Science Foundation of China (No. 61271079).
- Langley P, Bowers EJ, Murray A: Principal component analysis as a tool for analyzing beat-to-beat changes in ECG features: application to ECG-derived respiration. IEEE Trans Biomed Eng 2010, 57: 821–829.View ArticleGoogle Scholar
- Jukka A, Lipponen , Tarvainen MP: Principal component model for maternal ECG extraction in fetal QRS detection. Physiol Meas 2014, 35: 1637–1648. 10.1088/0967-3334/35/8/1637View ArticleGoogle Scholar
- Liu C, Li P, Di Maria C, Zhao L, Zhang H, Chen Z: A multi-step method with signal quality assessment and fine-tuning procedure to locate maternal and fetal QRS complexes from abdominal ECG recordings. Physiol Meas 2014, 35: 1665–1683. 10.1088/0967-3334/35/8/1665View ArticleGoogle Scholar
- Wei L, Rajapakse JC: ICA with Reference. Neurocomputing 2006, 69: 2244–2257. 10.1016/j.neucom.2005.06.021View ArticleGoogle Scholar
- Martin-Clemente R, Jose Luis C-O, Susana H-M, Mar E, Isabel R: Fast Technique for noninvasive fetal ECG extraction. IEEE Trans Biomed Eng 2011, 58: 227–230.View ArticleGoogle Scholar
- Shi Z, Zhang C: Semi-blind source extraction for fetal electrocardiogram extraction by combining non-Gaussianity and time-correlation. Neurocomputing 2007, 70: 1547–1581. 10.1016/j.neucom.2006.11.007View ArticleGoogle Scholar
- Zhang H, Shi Z, Guo C, Feng E: Semi-blind source extraction algorithm for fetal electrocardiogram based on generalized autocorrelations and reference signals. J Comput Appl Math 2009, 223: 409–420. 10.1016/j.cam.2008.01.031MathSciNetView ArticleGoogle Scholar
- Sameni R, Jutten C, Shamsollahi MB: Multichannel Electrocardiogram Decomposition Using Periodic Component Analysis. IEEE Trans Biomed Eng 2008, 55: 1935–1940.View ArticleGoogle Scholar
- Tsalaile T, Sameni R, Sanei S, Jutten C, Chambers J: Sequential blind source extraction for quasi-periodic signals with time-varying period. IEEE Trans Biomed Eng 2009, 56: 646–655.View ArticleGoogle Scholar
- Shuicai W, Yanni S, Zhuhuang Z, Lan L, Yanjun Z, Xiaofeng G: Research of fetal ECG extraction using wavelet analysis and adaptive filtering. Comput Biol Med 2013, 43: 1622–1627. 10.1016/j.compbiomed.2013.07.028View ArticleGoogle Scholar
- Almeida R, Goncalves H, Bernardes J, Ana Paula R: Fetal QRS detection and heart rate estimation: a wavelet-based approach. Physiol Meas 2014, 35: 1723–1735. 10.1088/0967-3334/35/8/1723View ArticleGoogle Scholar
- Puthusserypady S: Extraction of fetal electrocardiogram using H(infinity) adaptive algorithms. Med Bio Eng Comput 2007, 45: 927–937. 10.1007/s11517-007-0235-5View ArticleGoogle Scholar
- Assaleh K: Extraction of fetal electrocardiogram using adaptive neuro-fuzzy inference systems. IEEE Trans Biomed Eng 2007, 54: 59–68.View ArticleGoogle Scholar
- Assaleh K, Al-Nashash H: A novel technique for the extraction of fetal ECG using polynomial networks. IEEE Trans Biomed Eng 2006, 52: 1148–1152.View ArticleGoogle Scholar
- Zeng XP, Li SH, GuoJun Li Y, Zhou DHM: Fetal ECG extraction by combining single-channel SVD and cyclostationarity-based blind source separation. Int J Signal Process, Image Process Pattern Recognit 2013, 6: 367–376. 10.14257/ijsip.2013.6.5.32View ArticleGoogle Scholar
- Wei Z, Hongxing L, Aijun H, Xinbao N, Jianchun C: Single-lead fetal electrocardiogram estimation by means of combining R-peak detection, resampling and comb filter. Med Eng Phys 2010, 32: 708–719. 10.1016/j.medengphy.2010.04.012View ArticleGoogle Scholar
- Comani S, Mantini D, Lagatta A, Esposito F, Di Luzio S, Romani GL: Time course reconstruction of fetal cardiac signals from fMCG: independent component analysis versus adaptive maternal beat subtraction. Physio Meas 2004, 25: 1305–1321. 10.1088/0967-3334/25/5/019View ArticleGoogle Scholar
- Mihaela U, Johannes WM B, Swan Guid O, Rodica S: Fetal ECG extraction during labor using an adaptive maternal beat subtraction technique. Biomed Tech 2007, 52: 56–60. 10.1515/BMT.2007.011View ArticleGoogle Scholar
- Cerutti S, Baselli G, Civardi S, Ferrazzi E, Anna Maria M, Massimo P, Giorgio P: Variability analysis of fetal heart rate signals as obtained from abdominal electrocardiographic recordings. J Perinat Med 1986, 14: 445–452. 10.1515/jpme.19188.8.131.525View ArticleGoogle Scholar
- Martens SMM, Chiara R, Massimo M, Sluijter RJ: A robust fetal ECG detection method for abdominal recordings. Physiol Meas 2007, 28: 373–388. 10.1088/0967-3334/28/4/004View ArticleGoogle Scholar
- Vullings R, Peters CHL, Sluijter RJ, Mischi M, Oei SG, Bergmans JWM: Dynamic segmentation and linear prediction for maternal ECG removal in antenatal abdominal recordings. Physiol Meas 2009, 30: 291–307. 10.1088/0967-3334/30/3/005View ArticleGoogle Scholar
- Andreotti F, Riedl M, Himmelsbach T, Wedekind D, Wessel N, Stepan H, Schmieder C, Jank A, Malberg H, Zaunseder S: Robust fetal ECG extraction and detection from abdominal leads. Physiol Meas 2014, 35: 1551–1567. 10.1088/0967-3334/35/8/1551View ArticleGoogle Scholar
- Christov I, Simova I, Abacherli R: Extraction of the fetal ECG in noninvasive recordings by signal decompositions. Physiol Meas 2014, 35: 1713–1721. 10.1088/0967-3334/35/8/1713View ArticleGoogle Scholar
- Varanini M, Tartarisco G, Billeci L, Macerata A, Pioggia G, Balocchi R: An efficient unsuperivsed fetal QRS complex detection from abdominal maternal ECG. Physiol Meas 2014, 35: 16.7–1619. 10.1088/0967-3334/35/2/167View ArticleGoogle Scholar
- Haghpanahi M, Borkholder DA: Fetal QRS extraction from abdominal recordings via model-based signal processing and intelligent signal merging. Physiol Meas 2014, 35: 1591–1605. 10.1088/0967-3334/35/8/1591View ArticleGoogle Scholar
- Clifford GD, Ikaro S, Joachim B, Moody GB: Non-invasive fetal ECG analysis. Physiol Meas 2014, 35: 1521–1536. 10.1088/0967-3334/35/8/1521View ArticleGoogle Scholar
- Rajendra Acharya U, Paul Joseph K, Kannathal N, Choo Min L, Suri JS: Hear rate variability: a review. Med Bio Eng Comput 2004, 44: 1031–1051.View ArticleGoogle Scholar
- Leif S, Pablo L: Bioelectrical Signal Processing in Cardiac and Neurological Applications. Academic Press; 2005.Google Scholar
- Guyton AC, Hall JE: Textbook of Medical Physiology. 10th edition. Philadelphia: W B Saunders; 2000.Google Scholar
- Goldberger AL, Amaral LAN, Glass L, Hausdorff JM, Ivanov PC, Mark RG, Mietus JE, Moody GB, Peng C-K, Stanley HE: PhysioBank, PhysioToolkit, and PhysioNet: Components of a New Research Resource for Complex Physiologic Signals. Circulation 2000, 101: e215-e220. Circulation Electronic Pages [http://circ.ahajournals.org/content/101/23/e215.full] 10.1161/01.CIR.101.23.e215View ArticleGoogle Scholar
- Jezewski J, Matonia A, Kupka T, Roj D, Czabanski R: Determination of the fetal heart rate from abdominal signals: evaluation of beat-to-beat accuracy in relation to the direct fetal electrocardiogram. Biomed Tech 2012, 57: 383–394.View ArticleGoogle Scholar
- Mao L, Zhang G, Sun J: P wave detection in ECG signal based on location estimation and recognition post-processing. Signal Process 2009, 25: 948–952.Google Scholar
- McSharry PE, Clifford GD, Tarassenko L, Smith L: A dynamical model for generating synthetic electrocardiogram signals. IEEE Trans Biomed Eng 2003, 50: 289–294. 10.1109/TBME.2003.808805View ArticleGoogle Scholar
- Philip Saul J: Beat-to-beat variations of heart rate reflect modulation of cardiac autonomic outflow. News Physiol Sci 1990, 5: 32–37.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.