Next Article in Journal
Towards the Internet of Smart Trains: A Review on Industrial IoT-Connected Railways
Next Article in Special Issue
A Smartphone-Based Automatic Measurement Method for Colorimetric pH Detection Using a Color Adaptation Algorithm
Previous Article in Journal
Analysis of Maneuvering Targets with Complex Motions by Two-Dimensional Product Modified Lv’s Distribution for Quadratic Frequency Modulation Signals
Previous Article in Special Issue
Design and Validation of a Breathing Detection System for Scuba Divers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single-Lead Fetal ECG Extraction Based on a Parallel Marginalized Particle Filter

1
College of Electronics and Information, Hangzhou Dianzi University; Hangdian Smart City Research Center of Zhejiang Province, Hangzhou 310018, China
2
School of Communication Engineering, Hangzhou Dianzi University, Hangzhou 310018, China
3
Women’s Hospital School of Medicine Zhejiang University, Hangzhou 310006, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(6), 1456; https://doi.org/10.3390/s17061456
Submission received: 2 May 2017 / Revised: 15 June 2017 / Accepted: 17 June 2017 / Published: 21 June 2017

Abstract

:
This paper presents a novel method for extracting the fetal ECG (FECG) from a single-lead abdominal signal. A dynamical model for a modified abdominal signal is proposed, in which both the maternal ECG (MECG) and the FECG are modeled, and then a parallel marginalized particle filter (par-MPF) is used for tracking the abdominal signal. Finally, the FECG and MECG are simultaneously separated. Several experiments are conducted using both simulated and clinical signals. The results indicate that the method proposed in this paper effectively extracts the FECG and outperforms other Bayesian filtering algorithms.

1. Introduction

Fetal ECG (FECG) reflects changes in fetal heart activity and is used as a primary method for evaluating fetal heart status. However, there is no uniform, effective, non-invasive FECG signal elicitation and processing technique, partly because the FECG obtained with an abdominal sensor has a low signal-to-noise ratio and is contaminated by strong interference [1], including from the maternal ECG (MECG), fetal brain activity, and movement artifacts. An additional problem is the lack of a standard FECG database as well as a lack of knowledge of fetal cardiac function. Hence, the analysis of FECG is still in its infancy.
Currently, sophisticated and precise biomedical amplifiers can be used to obtain ECG signals. Accordingly, extracting the FECG from these recordings has gained considerable attention [2]. A number of effective algorithms have been proposed for FECG extraction, and most of these techniques require a multi-lead input. However, the number of leads affects the performance of the algorithm. Extracting the FECG with a single lead would both be more convenient and have a lower cost. For this purpose, many methods have been examined.
Adaptive filtering [3,4] is a simple and fast approach for FECG extraction, and conventional adaptive filtering is based on training an adaptive filter to remove the MECG using a maternal reference signal. However, the fetal ECG extracted using this algorithm is still influenced by the maternal ECG and other disturbances. An adaptive comb filter [5] is used for the FECG estimation from the abdominal signal, and the filter can adjust itself to the temporal variations in the fundamental frequency, which makes it qualified for the estimation of a quasi-periodic component from the abdominal signal. An important advantage of the adaptive comb filter over conventional adaptive filters is its ability to separate the MECG and FECG with a temporal overlap; additionally, it requires only one sensor.
Wavelet-based [6,7,8] methods have shown promising results. The properties of the wavelet transformation allow the extraction of the ECG waveforms from noise and artifact cancellation. These methods require the proper selection of the mother wavelet and scale to identify the frequency components in the signal. However, the selection of the mother wavelet is extremely important for the FECG and MECG delineation. An optimal wavelet should consider the edge variation of the ECG waveforms, which may be adapted to the signal morphology because a different mother wavelet can provide different results. In the case of waveform reconstruction, EMD [9] for respiratory signal evaluation shows better performance.
Extended Kalman filtering (EKF)/extended Kalman smoother (EKS) [10,11,12] has also been applied to extract FECG, which is a synthetic dynamic ECG model within an EKF framework that has been extended to jointly model several ECGs to extract the desired ECGs from an abdominal recording of maternal and fetal ECGs and noise. The results show that model can obtain the FECG and has a high signal-to-noise ratio. However, these methods should approximate the posterior density with a Gaussian density. If the assumption does not hold, the prevailing Bayesian filtering methods that are based on the Gaussian model [13] may fail to address the non-Gaussian and nonlinear FECG signals. Thus, another technique, referred to as MPF, is proposed in this paper.
In this paper, a modified abdominal signal dynamical model containing MECG and FECG is proposed. Both of the new ECG dynamical models have 17 state variables that allow the artificial ECG to adapt to normal or abnormal morphologies. A parallel marginalized particle filter is then used for tracking the abdominal signal. Finally, the FECG and MECG are separated at the same time from an abdominal signal. To evaluate the proposed Bayesian method, both simulated signal experiments and clinical signal experiments were conducted in this paper. Additionally, the performance of the filter was compared with several other Bayesian algorithms.
The rest of this paper is organized as follows: in Section 2, the previously developed synthetic ECG dynamical models are described. The proposed methodology is fully described in Section 3. In Section 4, a description of the performance indices that have been used to estimate the performance of the proposed methods is presented. The results and discussion of the different filtering methods are presented in Section 5; finally, in Section 6, some of the main conclusions are noted.

2. ECG Dynamical Model

The ECG dynamical equations, which describe the states used in the KF, were proposed by McSharry in 2003 [14]. The dynamical model is used for generating a synthetic ECG with Cartesian coordinates and consists of three nonlinear dynamical state equations, as follows:
{ x ˙ = α x ω y y ˙ = α y + ω x z ˙ = i { Q , P , R , S , T } a i Δ θ i exp ( Δ θ i 2 b i ) ( z z 0 )
where α = 1 x 2 + y 2 , Δ θ i = ( θ ξ i ) mod ( 2 π ) , θ = arctan 2 ( y , x ) , π arctan 2 ( y , x ) π , ω is the angular velocity of the trajectory as it moves around the limit cycle, and z 0 is used for modeling the baseline wander of the ECG signal. Each of the P, Q, R, S, and T waves of the ECG waveform is modeled with a Gaussian function, and all are located at specific angular positions ξ i , and a i , b i in Equation (1), and they correspond to the amplitude and width parameters of the Gaussian function.
Then, a simplified dynamic model in its discrete form was used, as proposed by Sameni [15], with the assumption of a small sampling period as follows:
{ θ k + 1 = ( θ k + ω δ ) mod ( 2 π ) z k + 1 = i { Q , P , R , S , T } δ α i ω b i 2 Δ θ i , k exp ( Δ θ i , k 2 2 b i 2 ) + z k + η k
where θ is the phase, z is the amplitude of the ECG, and η is a random additive noise.

3. Proposed Methodology for FECG Extraction

3.1. Modified Abdominal Signal Dynamical Model

It is well understood that the abdominal electrocardiogram (AECG) is a mixture of the MECG, FECG and noise. Thus, we consider the amplitude z and phase θ in Equation (3) of the FECG and the amplitude z and phase θ in Equation (4) of the MECG as state variables so that the new dynamic model can be used to estimate the MECG and FECG simultaneously:
{ θ k + 1 f = ( θ k f + ω f δ ) mod ( 2 π ) z k + 1 f = i = 1 5 δ α i f ω f ( b i f ) 2 Δ θ i , k f exp ( ( Δ θ i , k f ) 2 2 ( b i f ) 2 ) + z k f + η k f
{ θ k + 1 m = ( θ k m + ω m δ ) mod ( 2 π ) z k + 1 f = i = 1 5 δ α i f ω f ( b i f ) 2 Δ θ i , k f exp ( ( Δ θ i , k f ) 2 2 ( b i f ) 2 ) + z k f + η k f
In addition, most previous Bayesian ECG analyses directly use Equation (2) as the state space representation instead of the basic parameters of the Gaussian wave. As a result, occasional morphologic disturbances from random muscle artifacts have considerable effects on the performance of the filter. In these cases, the filter does not have sufficient time to adapt the Gaussian parameters, and a slight phase error of the model can lead to large errors in the Gaussian parameters [16]. The morphology of the ECG is determined by the Gauss parameters α i , b i , and ξ i ; therefore, α i , b i , and ξ i are also considered state variables following the random noise ε as follows:
{ α k + 1 f = α k f + η α , k b k + 1 f = b k f + η b , k ξ k + 1 f = ξ k f + η ξ , k
{ α k + 1 m = α k m + η α , k b k + 1 m = b k m + η b , k ξ k + 1 m = ξ k m + η ξ , k
where the superscripts f and m in Equations (5) and (6) denote the fetus and mother, respectively. The modified abdominal ECG dynamical model is shown in Equation (7), x k represents state variables, and η k represents process noise. The initial values of the dynamical model parameters can be obtained using the automatic off-line parameter selection procedure proposed in [15]:
{ x k = [ θ k f , z k f , ( α k f ) T , ( b k f ) T , ( ξ k f ) T , θ k m , z k m , ( α k m ) T , ( b k m ) T , ( ξ k m ) T ] T η k = [ η θ , k f , η z , k f , ( η α , k f ) T , ( η b , k f ) T , ( η ξ , k f ) T , η θ , k m , η z , k m , ( η α , k m ) T , ( η b , k m ) T , ( η ξ , k m ) T ] T x k + 1 = x k + η k

3.2. Observation Equation of the Dynamical Model

The phase observation and the noise ECG measurements are related to the measurement vectors described by Sameni [15]. The MECGs can be between five and 10 times greater in intensity than the FECGs, and thus, the extracted FECG signals will be distorted when only the abdominal signal is considered an observation. Thus, in this paper, two auxiliary observation equations are added as observations. These equations are derived from Equation (2), and both are deformation formulas of Equation (2):
s k z k 1 f + i = 1 5 δ α i f ω f ( b i f ) 2 Δ θ i , k 1 f exp ( ( Δ θ i , k 1 f ) 2 2 ( b i f ) 2 ) = z k m + v k m
s k z k 1 m + i = 1 5 δ α i m ω m ( b i m ) 2 Δ θ i , k 1 m exp ( ( Δ θ i , k 1 m ) 2 2 ( b i m ) 2 ) = z k f + v k f
where Equation (8) is the auxiliary FECG observation equation, and Equation (9) is the auxiliary MECG observation equation. Both equations can be used to estimate the MECG and FECG more accurately. The corresponding observation equations are as follows:
{ φ k f = θ k f + u k f s k z k 1 m + i = 1 5 δ α i m ω m ( b i m ) 2 Δ θ i , k 1 m exp ( ( Δ θ i , k 1 m ) 2 2 ( b i m ) 2 ) = z k f + v k f φ k m = θ k m + u k m s k z k 1 f + i = 1 5 δ α i f ω f ( b i f ) 2 Δ θ i , k 1 f exp ( ( Δ θ i , k 1 f ) 2 2 ( b i f ) 2 ) = z k m + v k m
where u k f , v k f , u k m , and v k m are observation noise. The advantage of MPF is that it can approximate non-Gaussian filtering problems; however, to compare it with the Kalman filter, we assumed that the observation noise is Gaussian white noise.

3.3. Parallel-Marginalized Particle Filter

The basic idea of a marginalized particle filter (Rao-Blackwellized particle filter) divides the state vector into linear components and nonlinear components [13]. Using the Rao-Blackwell theorem, we marginalized out the linear state variables and estimated them using the Kalman filter, while the particle filter was used to estimate the nonlinear state variables [17]. Thus, the optimal estimation of the state variables was obtained in the linear condition by the Kalman filter, while the dimension of the state variables was reduced in the nonlinear condition by the particle filter. Therefore, the marginalized particle filter algorithm greatly improves the accuracy when using the same number of particles.
In this paper, the par-MPF algorithm is proposed based on the improved AECG dynamic model, which was used to estimate the abdominal signal. Finally, the FECG and MECG can be obtained in parallel. We consider Gaussian amplitude parameters α k f , α k m as linear variables and the other parameters as the nonlinear variables [17]. Therefore, the state variable x k can be split into x k L and x k N L , which are referred to as linear state and nonlinear variables, respectively:
{ x k L = [ ( α k f ) T , ( α k m ) T ] T x k N L = [ θ k f , z k f , ( b k f ) T , ( ξ k f ) T , θ k m , z k m , ( b k m ) T , ( ξ k m ) T ] T
Using the Bayesian theorem, we can calculate the following:
p ( x k N L , x k L | Y k ) = p ( x k L | x k N L , Y k ) p ( x k N L , Y k )
where p ( x k L | x k N L , Y k ) is estimated by the Kalman filter, and p ( x k N L , Y k ) is estimated by the particle filter.
In summary, the par-MPF algorithm presented in this paper includes the following steps:
  • Step 1: Initialize the M particles:
    { { x 0 | 1 N L , ( i ) } i = 1 M ~ p ( x 0 N L ) { x 0 | 1 L , ( i ) } i = 1 M ~ p ( x 0 L ) ,
    and set k = 0 .
  • Step 2: Calculate the importance:
    weight:
    w k ( i ) = p ( Y k | x k | k 1 N L , ( i ) )   i = 1 , 2 , , M ;
    normalization:
    w k ( i ) = w k ( i ) / i = 1 M w k ( i ) .
  • Step 3: Resample the M particles:
    p ( x k | k N L , ( i ) = x k | k 1 N L , ( j ) ) = w k ( j )
  • Step 4: Update the Kalman filter measurement:
    { α ˜ k | k f / m , ( i ) = α ˜ k | k 1 f / m , ( i ) P k | k f / m = P k | k 1 f / m ,
    where the superscript f / m means that there are two cases, i.e., the fetal ECG and the maternal ECG estimates.
  • Step 5: Update the particle filter time:
    p ( x k + 1 | k f / m N L , ( i ) | x k | k f / m N L , ( i ) ) = N ( μ k f / m , Σ k f / m )
    where:
    { μ k f / m = [ θ k f / m , ( i ) + w f / m δ , g T ( x k f / m N L , ( i ) ) α k f / m , ( i ) + z k f / m , ( i ) , b k f / m , ( i ) , ξ k f / m , ( i ) ] T Σ k f / m = d i a g ( σ θ f / m 2 , σ z f / m 2 , σ b f / m 2 , σ ξ f / m 2 ) g ( x k N L ) = [ g ( x k , 1 N L ) , g ( x k , 2 N L ) , g ( x k , 3 N L ) , g ( x k , 4 N L ) , g ( x k , 5 N L ) ] g ( x k , j N L ) = δ α j , k ω b j , k 2 Δ θ j , k exp ( Δ θ j , k 2 2 b j , k 2 )
  • Step 6: Update the Kalman filter time:
    α ˜ k + 1 | k f / m , ( i ) = α ˜ k | k f / m , ( i ) + L k f / m ( z ˜ k + 1 | k f / m , ( i ) z k f / m , ( i ) g T ( x k + 1 | k f / m N L , ( i ) ) α ˜ k | k f / m , ( i ) )
    where:
    { F k f / m = g T ( x ˜ k + 1 | k f / m N L , ( i ) ) P k | k f / m , ( i ) g ( x ˜ k + 1 | k f / m N L , ( i ) ) + Q k f / m N L L k f / m = P k | k f / m , ( i ) g ( x ˜ k + 1 | k f / m N L , ( i ) ) ( F k f / m ) 1 ,
    P is the process noise covariance matrix, and Q is the observation noise covariance matrix.
  • Step 7:
    k = k + 1 ,
    and iterate from step 2.

4. Performance Index of the FECG Extraction

4.1. Performance Index of the Simulated Data

For the simulated data, the percentage root mean-square difference ( P R D ) and the signal-to-noise ratio enhancement ( S N R e n h ) were used to estimate the performance of filtering algorithms. The percentage root mean-square difference is defined as follows [18]:
P R D = n ( f ( n ) f ^ ( n ) ) 2 n ( f ( n ) ) 2
where f ( n ) is the original FECG signal, f ^ ( n ) denotes the extracted FECG, and PRD is the most frequently used interference measurement tool that enumerates the error between the original signal and the reconstructed signal. The nearer the PRD is to zero, the more it indicates the similarity between the two signals and, hence, indicates the signal quality.
The signal-to-noise ratio enhancement is defined as follows [19]:
S N R e n h = 10 log n ( a ( n ) f ( n ) ) 2 n ( f ( n ) f ^ ( n ) ) 2
where a ( n ) is the abdominal signal, f ( n ) is the original FECG signal, and f ^ ( n ) denotes the extracted FECG. S N R e n h denotes the difference between the SNR before and after filtering. A larger S N R e n h indicates a better performance of the algorithm.

4.2. Performance Index of the Clinical Data

For the clinical data, the signal-to-noise ratio is based on an eigenvalue analysis and a cross correlation can be used to estimate the performance of filtering algorithms. Jagannath D.J. et al. [18], used SNR, including S N R svd , the quotient between the first singular value (fetus signal energy) and the sum of the rest of the singular values (noise energy), and S N R c o r , as an indicator of the goodness of the extraction method. For a square matrix, similarly, the eigenvalue was used to represent signal energy instead. Both of the SNR values can be calculated by the following steps [6]:
  • Step 1: The extracted FECG is divided into N pieces crossing to the R peak, and each piece includes M samples and one QRS complex.
  • Step 2: All of the pieces are stored in columns of an M by N matrix U , and the vectors u ( k ) are the zero mean and are normalized to the unit length; that is,
    u T ( k ) u ( k ) = 1 .
  • Step 3: A signal-to-noise ratio based on eigenvalues can be calculated as:
    S N R e i g = λ max N λ max
    where λ contains the eigenvalues of U .
  • Step 4: A signal-to-noise ratio based on the cross-correlation coefficients can be calculated as
    S N R c o r = η 1 η
    where
    η = 2 M ( M 1 ) i = 0 M 2 k = i + 1 M 1 f ( i ) T f ( k ) ,
    and f is the FECG.

5. Results and Analysis

To verify the feasibility and reliability of the proposed algorithm, both simulated and clinical data were used to study the performance of the proposed method, and other methods were tested with the same signals, including EKS, EKF and the Bayesian adaptive neuro fuzzy inference system [18]. The main ideas of these algorithms are shown in Table 1.

5.1. FECG Extraction on the Simulated Data

5.1.1. Simulated Data without Noise

An artificial FECG-MECG mixture simulator, as proposed by Behar [20], was used for creating the simulated data, where the MECG was simulated with a beat rate of 70 bpm, and the FECG was simulated with a beat rate of 120 bpm. The sampling frequency was set to 250 Hz, and the acquisition time was 10 s. The parameter S N R f m = n f e c g ( n ) 2 n m e c g ( n ) 2 = 10 dB , and a synthetic AECG signal is shown in Figure 1. The simulated FECG signal is shown in Figure 2.
The FECG signals extracted are shown in Figure 3. It is clear that all three algorithms correctly separate the mixed signal, and from the extracted FECG, we can see that unlike par-MPF, it does not fail when the MECG and FECG overlap. This performance was particularly noted for the samples between 700 and 800 and between 2200 and 2300 in Figure 3, in which some parts of the FECG signal deteriorated during the MECG extraction by the EKS and EKF methods. In contrast, the proposed par-MPF jointly models the FECG and MECG, resulting in a better estimate of the FECG. The reason for this is that some parts of the FECG signal deteriorated during the MECG extraction by the EKS and EKF methods. In contrast, the proposed method obtains the MECG and FECG at the same time and results in a better estimate of the FECG.
Another 10 FECG datasets, in which S N R f m = 2 , 4 , 18 , 20 dB , were used to further verify the proposed algorithm. Three filtering algorithms were imposed on these signals to split the MECG and FECG, and the P R D and S N R e n h of the extracted FECG are shown in Figure 4a,b, respectively.
As shown in Figure 4, when the S N R f m increases, the P R D value of the par-MPF is smaller than the EKS and EKF, while the S N R e n h is higher than the EKS and EKF. Overall, the proposed method has a stable performance; thus, the algorithm can be applied to the simulated data with noise.

5.1.2. Simulated Data with Noise

The artificial FECG-MECG mixture simulator was used again to generate the AECG with a motion artifact, as shown in Figure 5, in which S N R f m = 10   dB , and S N R m n = 6   dB , keeping the other conditions the same:
S N R m n = n [ f e c g ( n ) 2 + m e c g ( n ) 2 ] n n o i s e ( n ) 2
To verify the feasibility of the proposed algorithm, the AECG in Figure 5 was processed by the three filtering methods. The extracted FECG is shown in Figure 6. The EKS and EKF also fail when the MECG and FECG overlap, while the par-MPF does not fail. This experiment shows that the algorithm is feasible for the FECG separation of the clinical signal.
We changed the values of S N R m n , after which the P R D and S N R e n h were calculated to evaluate the performance of these filtering methods. The results of the P R D and S N R e n h calculations are shown in Figure 7a,b, respectively, when S N R m n = 2 , 4 , , 18 , 20   dB .
Figure 7 clearly demonstrates that when S N R m n increases, the P R D of the par-MPF is smaller than the EKS and EKF, while the S N R e n h is higher than the EKS and EKF. In summary, when the AECG is mixed with different noise, the proposed algorithm shows a perfect filtering performance, and the extracted FECG is high quality. Therefore, the application of the filtering algorithm to the clinical AECG should be considered.

5.2. FECG Extraction on the Different Database.

5.2.1. Database for the Identification of Systems

The clinical data are from the Database for the Identification of Systems (DaISy), which is a well-known ECG measured from a pregnant woman. A total of eight channels are available, which are sampled at 250 Hz with a duration of 10 s [21]. Five abdominal signals are shown in Figure 8. The FECG components in lead_1 were the most obvious, while the lead_4 signal had so much noise that the FECG components could hardly be seen.
We ran the proposed algorithm and the other three algorithms. Figure 9 presents the FECG extracted using the lead_1. The proposed method resulted in a superior estimation of the FECG compared to the other methods. When the MECG and FECG completely overlap in samples between 1600 and 1700, there is a small effect on the ANFIS-EKS and ANFIS-EKF outputs, while the FECGs suffer a severe loss when extracted by the EKS and EKF. As described above, there are two steps in extracting the FECG signal with the EKS and EKF algorithms. The FECG is obtained from the residual signal, which will reduce the QRS wave of the FECG component, while the par-MPF algorithm is based on the improved AECG model and extracts the MCG and FECG in a parallel manner. Thus, the FECG is estimated more accurately, and the performance of the algorithm is greatly improved.
The FECG signals obtained using the par-MPF from three other abdominal signals are shown in Figure 10. To further illustrate the advantages of our algorithm, another experiment was conducted using another AECG without lead_4. The S N R e i g and S N R c o r are calculated in Table 2.
From the experimental results, we can see that the signal-to-noise ratio of our algorithm is higher than that of the other four algorithms; furthermore, the difference between the signal-to-noise ratio of the par-MPF is small. It can be concluded that the performance of the algorithm is better than that of the other Bayesian algorithms.

5.2.2. Abdominal and Direct Fetal ECG Database

The clinical data used in this section are from the Abdominal and Direct Fetal ECG Database. The five-minute multichannel fetal ECG recordings, with cardiologist-verified annotations of all fetal heartbeats, were from five women in labor at the Medical University of Silesia, Poland. Each record includes four signals from the maternal abdomen and a simultaneously recorded reference direct fetal ECG from the fetal scalp; all of the signals were sampled at 1 KHz [22,23]. The first 5000 points were used to plot the ECG waveform as shown in Figure 11, where (a) is the fetal scalp ECG signal, and (b) to (e) are the abdominal signals.
The result of the filtering algorithms using the signal Abdomen_1 are shown in Figure 12. These fives algorithms adapted to the different shapes of the AECG, which included the interference of noise. However, the performance of the EKF was the worst among these algorithms, and the FECG was mixed with more noise. The results of the FECG extraction using EKS, ANFIS-EKS, and ANFIS-EKF had relatively low noise. In contrast, the FECG extraction using the par-MPF had the best quality.
The FECG signals obtained using the par-MPF by the three other abdominal signals of r01 are shown in Figure 13; the algorithm performance was evaluated by calculating the S N R e i g and S N R c o r using these three signals of r01, and the results are shown in Table 3. All of the algorithms extracted the FECG with different signal qualities; however, the proposed algorithm not only extracted the FECG without noise, but the quality of the signal was finer, and thus, it is helpful for the follow-up extraction of the R wave and FHRV analysis.
Finally, the signals r04, r07, r08 and r10 in the Abdominal and Direct Fetal ECG Database were processed using the above filtering algorithm [24,25,26]. Because the signals of Abdomen_2 of r07, Abdomen_1 of r08 and Abdomen_2 of r10 were poor, they were excluded from the experiments. The comparison of the S N R e i g and S N R c o r among the five algorithms is shown in Figure 14 and Figure 15, respectively.
In Figure 14 and Figure 15, the X-axis and Y-axis indicate the number of the signal used in these experiments and the SNR, respectively. From Figure 14 and Figure 15, we can clearly see that for all signals tested, the SNR of the proposed method is higher than that for KS, EKF, ANFIS-EKS and ANFIS-EKF. The higher SNR here means a better extraction performance. Figure 14 and Figure 15 show that the performance indices of the proposed method are obviously higher than those of the other algorithms, and no matter how the abdominal ECG changes, the algorithm maintains a good performance. The EKF achieved the lowest SNR and showed the poorest performance. The results indicate that the algorithm in this paper effectively extracted a high-quality FECG from different pregnant women. Therefore, this algorithm has great prospects for clinical application.

6. Conclusions

A Bayesian filtering framework for extracting the FECG from a single-channel mixture of the abdominal signal MECG and FECG, including noise, is proposed in this paper. A modified abdominal signal dynamical model is proposed to improve the filtering accuracy. The FECG and MECG are obtained at the same time, substantially reducing the running time of the algorithm. In addition, simulated and clinical data were used to verify the performance of the filter. All the results show that the proposed method successfully extracted the FECG in many scenarios, and the extracted FECG had a high quality.
Finally, the proposed method was compared with some other Bayesian filter algorithms, such as the EKS, EKF and the Bayesian adaptive neuro fuzzy inference system. The proposed method provided a superior estimation of the FECGs compared with the other methods, and when the MECG and FECG completely overlapped, there was little effect on the ANFIS-EKS and ANFIS-EKF. However, the FECGs suffered a severe loss when extracted by the EKS and EKF methods. Moreover, the performance indices of the proposed method were obviously higher than those of the other algorithms. Taken together, these experiments reveal that the strength of the proposed method is the best with respect to signal extraction. The proposed par-MPF algorithm has been verified to obtain a high signal-to-noise ratio FECG signal from a single-channel abdominal signal. However, whether the algorithm is suited for the extraction of twin FECGs is an issue for further study.

Acknowledgments

This research was supported by National Natural Science Foundation of China (Grant No. 61571173); the research was supported by Zhejiang Provincial Natural Science Foundation of China under Grant No. LY14F010021. The research was also supported by the Welfare Project of the Science Technology Department of Zhejiang Province (Grant No. 2015C31085, 2016C33079) and Smart City Collaborative Innovation Center of Zhejiang Province.

Author Contributions

Zhidong Zhao and Wen Xu conceived and performed the theoretical derivation; Huiling Tong and Yanjun Deng designed and performed the experiments; Wen Xu, Yefei Zhang and Haihui Ye designed and performed the analogue simulations; Wen Xu wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lukosevicius, M.; Marozas, V. Noninvasive fetal QRS detection using an echo state network and dynamic programming. Physiol. Meas. 2014, 35, 1685–1697. [Google Scholar] [CrossRef] [PubMed]
  2. Ayat, M.; Assaleh, K.; Al-Nashash, H. Extracting fetal ECG from a single maternal abdominal record. In Proceedings of the IEEE 8th GCC Conference and Exhibition (GCCCE), Muscat, Oman, 1–4 February 2015; pp. 1–4. [Google Scholar]
  3. Puthusserypady, S. Extraction of fetal electrocardiogram using H∞ adaptive algorithms. Med. Biol. Eng. Comput. 2007, 45, 927–937. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, W.; Liu, H.X.; Cheng, J.C. Adaptive filtering in phase space for foetal electrocardiogram estimation from an abdominal electrocardiogram signal and a thoracic electrocardiogram signal. IET Signal Process. 2012, 6, 171–177. [Google Scholar]
  5. Zheng, W.; Wei, X.Y.; Zhong, J.J.; Liu, H.X. Noninvasive fetal ECG estiomation using adaptive comb filter. Comput. Methods Progr. Biomed. 2013, 112, 125–134. [Google Scholar]
  6. Wu, S.; Shen, Y.; Zhou, Z.; Lin, L.; Zeng, Y.; Gao, X. Research of fetal ECG extraction using wavelet analysis and adaptive filtering. Comput. Biol. Med. 2013, 43, 1622–1627. [Google Scholar] [CrossRef] [PubMed]
  7. Elmansouri, K.; Latif, R.; Maoulainine, F. Efficient fetal heart rate extraction using undecimated wavelet transform. In Proceedings of the 2nd World Conference on Complex Systems (WCCS), Agadir, Morocco, 10–12 November 2014; pp. 696–701. [Google Scholar]
  8. Desai, K.D.; Sankhe, M.S. A real-time fetal ECG feature extraction using multiscale discrete wavelet transform. In Proceedings of the 5th International Conference on Biomedical Engineering and Informatics (BMEI), Chongqing, China, 16–18 October 2012; pp. 407–412. [Google Scholar]
  9. Campolo, M.; Labate, D.; Foresta, F.L.; Morabito, F.C. ECG-derived respiratory signal using Empirical Mode Decomposition. In Proceedings of the IEEE International Workshop on Medical Measurements and Applications, Bari, Italy, 30–31 May 2011; pp. 399–403. [Google Scholar]
  10. Niknazar, M.; Rivet, B.; Jutten, C. Fetal ECG extraction by extended state kalman filtering based on single-channel recordings. IEEE Trans. Biomed. Eng. 2013, 60, 1345–1352. [Google Scholar] [CrossRef] [PubMed]
  11. Zeng, X.; Zhou, X.; Li, G.; Liu, Q. Robust adaptive fetal heart rate estimation for single-channel abdominal ECG recording. In Proceedings of the 5th International Conference on Biomedical Engineering and Informatics (BMEI), Chongqing, China, 16–18 October 2012; pp. 587–591. [Google Scholar]
  12. Panigrahy, D.; Sahu, P.K. Extraction of fetal electrocardiogram (ECG) by extended state Kalman filtering and adaptive neuro-fuzzy inference system (ANFIS) based on single channel abdominal recording. Sadhana 2015, 40, 1091–1104. [Google Scholar] [CrossRef]
  13. Schon, T.; Gustafsson, F.; Nordlund, P. Marginalized particle filters for mixed linear/nonlinear state-space models. IEEE Trans. Signal Process. 2005, 53, 2279–2289. [Google Scholar] [CrossRef]
  14. McSharry, P.E.; Clifford, G.D.; Tarassenko, L.; Smith, L.A. A dynamical model for generating synthetic electrocardiogram signals. IEEE Trans. Biomed. Eng. 2003, 50, 289–292. [Google Scholar] [CrossRef] [PubMed]
  15. Sameni, R.; Shamsollahi, M.B.; Jutten, C.; Clifford, G.D. A nonlinear Bayesian filtering framework for ECG denoising. IEEE Trans. Biomed. Eng. 2007, 54, 2172–2185. [Google Scholar] [CrossRef] [PubMed]
  16. Li, G.; Zeng, X.; Lin, J.; Zhou, X. Genetic particle filtering for denoising of ECG corrupted by muscle artifacts. In Proceedings of the 8th International Conference on Natural Computation, Chongqing, China, 29–31 May 2012; pp. 562–565. [Google Scholar]
  17. Lin, C.; Bugallo, M.; Mailhes, C.; Tourneret, J.Y. ECG denoising using a dynamical model and a marginalized particle filter. In Proceedings of the Conference Record of the 45th Asilomar Conference on Signals, Systems and Computers (ASILOMAR), Pacific Grove, CA, USA, 6–9 November 2011; pp. 1679–1683. [Google Scholar]
  18. Jagannath, D.J.; Selvakumar, A.I. Superior foetal electrocardiogram signal elicitation using a novel artificial intelligent Bayesian methodology. Appl. Soft Comput. 2015, 37, 1002–1017. [Google Scholar] [CrossRef]
  19. Martinez, M.; Calpe, J.; Soria, E.; Guerrero, J.F. Methods to evaluate the performance of fetal electrocardiogram extraction algorithms. In Proceedings of the Computers in Cardiology 2001, Rotterdam, The Netherlands, 23–26 September 2001; pp. 253–256. [Google Scholar]
  20. Behar, J.; Andreotti, F.; Zaunseder, S.; Li, Q.; Oster, J.; Clifford, G.D. An ECG simulator for generating maternal-foetal activity mixtures on abdominal ECG recordings. Physiol. Meas. 2014, 35, 1537–1550. [Google Scholar] [CrossRef] [PubMed]
  21. L.LathauwerDaISy: Database for the Identification of Systems. Available online: http://www.esat.kuleuven.ac.be./sista/daisy (accessed on 18 June 2017).
  22. Jezewski, J.; Matonia, A.; Kupka, T.; Roj, D.; Czabanski, R. Determination of 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. [Google Scholar] [CrossRef] [PubMed]
  23. Abdominal and Direct Foetal ECG Database. Available online: https://www.physionet.org/physiobank/database/adfecgdb (accessed on 18 June 2017).
  24. Karvounis, E.C.; Tsipouras, M.G.; Fotiadis, D.I. Detection of Fetal Heart Rate through 3-D Phase Space Analysis from Multivariate Abdominal Recordings. IEEE Trans. Biomed. Eng. 2009, 56, 1394–1406. [Google Scholar] [CrossRef] [PubMed]
  25. Castillo, E.; Morales, D.P.; Botella, G.; Garcia, A.; Parrilla, L.; Palma, A.J. Efficient wavelet-based ECG processing for single-lead FHR extraction. Digit. Signal Process. 2013, 23, 1897–1909. [Google Scholar] [CrossRef]
  26. Castillo, E.; Morales, D.P.; Garcia, A.; Martinez-Marti, F.; Parrilla, L.; Palma, A.J. Noise Suppression in ECG Signals through Efficient One-Step Wavelet Processing Techniques. J. Appl. Math. 2013, 2013, 4. [Google Scholar] [CrossRef]
Figure 1. Simulated data.
Figure 1. Simulated data.
Sensors 17 01456 g001
Figure 2. Simulated FECG.
Figure 2. Simulated FECG.
Sensors 17 01456 g002
Figure 3. FECG estimation: (a) FECG extracted by par-MPF; (b) FECG extracted by EKS; (c) FECG extracted by EKF.
Figure 3. FECG estimation: (a) FECG extracted by par-MPF; (b) FECG extracted by EKS; (c) FECG extracted by EKF.
Sensors 17 01456 g003
Figure 4. Comparison of performance indices: (a) PRD; (b) S N R e n h .
Figure 4. Comparison of performance indices: (a) PRD; (b) S N R e n h .
Sensors 17 01456 g004
Figure 5. Simulated data with noise.
Figure 5. Simulated data with noise.
Sensors 17 01456 g005
Figure 6. FECG estimation: (a) FECG extracted by par-MPF; (b) FECG extracted by EKS; (c) FECG extracted by EKF.
Figure 6. FECG estimation: (a) FECG extracted by par-MPF; (b) FECG extracted by EKS; (c) FECG extracted by EKF.
Sensors 17 01456 g006
Figure 7. Comparison of performance indices: (a) PRD; (b) S N R e n h .
Figure 7. Comparison of performance indices: (a) PRD; (b) S N R e n h .
Sensors 17 01456 g007
Figure 8. AECG: (a) Lead_1; (b) Lead_2; (c) Lead_3; (d) Lead_4; (e) Lead_5.
Figure 8. AECG: (a) Lead_1; (b) Lead_2; (c) Lead_3; (d) Lead_4; (e) Lead_5.
Sensors 17 01456 g008aSensors 17 01456 g008b
Figure 9. FECG estimation: (a) FECG extracted by par-MPF; (b) FECG extracted by EKS; (c) FECG extracted by EKF; (d) FECG extracted by ANFIS-EKS; (e) FECG extracted by ANFIS-EKF.
Figure 9. FECG estimation: (a) FECG extracted by par-MPF; (b) FECG extracted by EKS; (c) FECG extracted by EKF; (d) FECG extracted by ANFIS-EKS; (e) FECG extracted by ANFIS-EKF.
Sensors 17 01456 g009aSensors 17 01456 g009b
Figure 10. FECG estimation: (a) FECG extracted from lead_2; (b) FECG extracted from lead_3; (c) FECG extracted from lead_5.
Figure 10. FECG estimation: (a) FECG extracted from lead_2; (b) FECG extracted from lead_3; (c) FECG extracted from lead_5.
Sensors 17 01456 g010
Figure 11. Signal r01: (a) Direct_1; (b) Abdomen_1; (c) Abdomen_2; (d) Abdomen_3; (e) Abdomen_4.
Figure 11. Signal r01: (a) Direct_1; (b) Abdomen_1; (c) Abdomen_2; (d) Abdomen_3; (e) Abdomen_4.
Sensors 17 01456 g011
Figure 12. FECG estimation: (a) FECG extracted by par-MPF; (b) FECG extracted by EKS; (c) FECG extracted by EKF; (d) FECG extracted by ANFIS-EKS; (e) FECG extracted by ANFIS-EKF.
Figure 12. FECG estimation: (a) FECG extracted by par-MPF; (b) FECG extracted by EKS; (c) FECG extracted by EKF; (d) FECG extracted by ANFIS-EKS; (e) FECG extracted by ANFIS-EKF.
Sensors 17 01456 g012
Figure 13. FECG estimation: (a) FECG extracted from Abdomen_2; (b) FECG extracted from Abdomen_3; (c) FECG extracted from Abdomen_4.
Figure 13. FECG estimation: (a) FECG extracted from Abdomen_2; (b) FECG extracted from Abdomen_3; (c) FECG extracted from Abdomen_4.
Sensors 17 01456 g013
Figure 14. Comparison of the S N R e i g .
Figure 14. Comparison of the S N R e i g .
Sensors 17 01456 g014
Figure 15. Comparison of the S N R c o r .
Figure 15. Comparison of the S N R c o r .
Sensors 17 01456 g015
Table 1. Filtering algorithms used in the experiments.
Table 1. Filtering algorithms used in the experiments.
Filtering AlgorithmMain Idea
par-MPFMarginalized particle filter
EKS [10]Extended Kalman filter + Smooth
EKF [11]Extended Kalman filter
ANFIS-EKS [12]EKS + Adaptive neuro fuzzy inference system
ANFIS-EKF [12]EKF + Adaptive neuro fuzzy inference system
Table 2. The SNR based on the eigenvalues and the cross-correlation coefficients.
Table 2. The SNR based on the eigenvalues and the cross-correlation coefficients.
MethodsLead 1Lead 2Lead 3Lead 4
S N R e i g S N R c o r S N R e i g S N R c o r S N R e i g S N R c o r S N R e i g S N R c o r
par-MPF2.08641.68202.08491.71202.41192.19262.10331.7037
EKS [10]1.74971.56771.82501.37141.39030.77971.14140.9533
EKF [11]1.22451.15061.23670.83921.35900.74650.62580.4130
ANFIS-EKS [12]1.83341.61521.87091.48891.44521.49171.18370.9856
ANFIS-EKF [12]1.50971.31601.54911.25571.37670.92470.84650.6551
Table 3. SNR based on the eigenvalues and the cross-correlation coefficients of r01.
Table 3. SNR based on the eigenvalues and the cross-correlation coefficients of r01.
MethodAbdomen-1Abdomen-2Abdomen-3Abdomen-4
S N R e i g S N R c o r S N R e i g S N R c o r S N R e i g S N R c o r S N R e i g S N R c o r
par-MPF12.889814.733212.146413.962411.469513.214413.35215.3952
EKS [10]3.50253.61333.86094.11445.03895.12866.74886.9799
EKF [11]2.42222.39022.2362.25963.58683.25834.26513.8168
ANFIS-EKS [12]4.22854.58794.14384.62565.54275.55497.20397.2104
ANFIS-EKF [12]3.00713.08732.55732.74173.90213.84445.57715.4508

Share and Cite

MDPI and ACS Style

Zhao, Z.; Tong, H.; Deng, Y.; Xu, W.; Zhang, Y.; Ye, H. Single-Lead Fetal ECG Extraction Based on a Parallel Marginalized Particle Filter. Sensors 2017, 17, 1456. https://doi.org/10.3390/s17061456

AMA Style

Zhao Z, Tong H, Deng Y, Xu W, Zhang Y, Ye H. Single-Lead Fetal ECG Extraction Based on a Parallel Marginalized Particle Filter. Sensors. 2017; 17(6):1456. https://doi.org/10.3390/s17061456

Chicago/Turabian Style

Zhao, Zhidong, Huiling Tong, Yanjun Deng, Wen Xu, Yefei Zhang, and Haihui Ye. 2017. "Single-Lead Fetal ECG Extraction Based on a Parallel Marginalized Particle Filter" Sensors 17, no. 6: 1456. https://doi.org/10.3390/s17061456

APA Style

Zhao, Z., Tong, H., Deng, Y., Xu, W., Zhang, Y., & Ye, H. (2017). Single-Lead Fetal ECG Extraction Based on a Parallel Marginalized Particle Filter. Sensors, 17(6), 1456. https://doi.org/10.3390/s17061456

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop