Computational and Mathematical Methods in Medicine

Volume 2018, Article ID 7061456, 7 pages

https://doi.org/10.1155/2018/7061456

## An Improved FastICA Method for Fetal ECG Extraction

College of Life Science and Bioengineering, Beijing University of Technology, Beijing, China

Correspondence should be addressed to Shuicai Wu; nc.ude.tujb@iaciuhsuw

Received 8 January 2018; Revised 3 April 2018; Accepted 22 April 2018; Published 17 May 2018

Academic Editor: Chuangyin Dang

Copyright © 2018 Li Yuan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

*Objective*. The fast fixed-point algorithm for independent component analysis (FastICA) has been widely used in fetal electrocardiogram (ECG) extraction. However, the FastICA algorithm is sensitive to the initial weight vector, which affects the convergence of the algorithm. In order to solve this problem, an improved FastICA method was proposed to extract fetal ECG.* Methods*. First, the maternal abdominal mixed signal was centralized and whitened, and the overrelaxation factor was incorporated into Newton’s iterative algorithm to process the initial weight vector randomly generated. The improved FastICA algorithm was used to separate the source components, selected the best maternal ECG from the separated source components, and detected the R-wave location of the maternal ECG. Finally, the maternal ECG component in each channel was removed by the singular value decomposition (SVD) method to obtain a clean fetal ECG signal.* Results*. An annotated clinical fetal ECG database was used to evaluate the improved algorithm and the conventional FastICA algorithm. The average number of iterations of the algorithm was reduced from 35 before the improvement to 13. Correspondingly, the average running time was reduced from 1.25 s to 1.04 s when using the improved algorithm. The signal-to-noise ratio (SNR) based on eigenvalues of the improved algorithm was 1.55, as compared to 0.99 of the conventional FastICA algorithm. The SNR based on cross-correlation coefficients of the conventional algorithm was also improved from 0.59 to 2.02. The sensitivity, positive predictive accuracy, and harmonic mean () of the improved method were 99.37%, 99.00%, and 99.19%, respectively, while these metrics of the conventional FastICA method were 99.03%, 98.53%, and 98.78%, respectively.* Conclusions*. The proposed improved FastICA algorithm based on the overrelaxation factor, while maintaining the rate of convergence, relaxes the requirement of initial weight vector, avoids the unbalanced convergence, reduces the number of iterations, and improves the convergence performance.

#### 1. Introduction

Electrocardiogram (ECG) is an important tool used by physicians for identifying abnormalities in the human heart activity [1]. Similarly, fetal ECG signal can also reflect electrophysiological activity of the fetal heart. Physicians can detect in time fetal abnormalities during fetal development through the fetal ECG waveform analysis, such as fetal distress and intrauterine hypoxia. In addition, a small number of abnormal fetal ECG waveforms are also a manifestation of congenital heart disease, for which early measures can be taken to reduce neonatal morbidity and mortality.

Currently, there are two ways for obtaining fetal ECG. One is the invasive scalp electrode method, which can directly measure the pure fetal ECG signal. However, it can only detect fetal ECG signal during the time of birth, and it is invasive so it may cause harm to the mother and the fetus. The other method is noninvasive, abdominal electrode method. The signals from the abdominal body surface are collected by placing an electrode patch in the abdomen of the mother, which allows for long-term monitoring during pregnancy without harming the mother or the fetus. However, the signals from maternal abdomen surface are very complex, which not only contain weak fetal ECG and maternal ECG but also include the mother’s respiratory noise, frequency interference, and other signals [2]. In particular, the magnitude of the maternal ECG detected in the abdomen is about 2–10 times that of the fetal ECG [3], which makes the extraction of fetal ECG difficult. Therefore, it is necessary to develop a noninvasive method that can extract fetal ECG effectively.

At present, fetal ECG extraction algorithms mainly include adaptive filtering [4, 5], wavelet analysis [6], matched filtering [7], blind source separation [8], independent component analysis (ICA) [9], neural network [10, 11], and singular value decomposition (SVD) [12]. Among these methods, ICA can separate the source signals from the mixed signals under the assumption that the source signals are statistically independent of each other, without needing any information regarding the source signals or the mixed matrix. Therefore, ICA is considered as a promising method for extracting fetal ECG. In recent years, researchers have proposed many improved ICA algorithms, which can separate non-Gaussian signals. Among them, because of its fast convergence, the fast fixed-point algorithm for independent component analysis (FastICA) [13] has been widely used in the extraction of fetal ECG. However, the FastICA algorithm is sensitive to the initial weight vector, and different initial weight vectors may lead to different convergence performances of the algorithm.

In this paper, an improved FastICA method was proposed to solve the above problem. By incorporating an overrelaxation factor into the iterative algorithm, the initial weight vector generated randomly can be relaxed. By choosing the appropriate overrelaxation factor, the iterative algorithm with slower convergence rate can converge, and the divergent iterative algorithm may become convergent.

#### 2. Methods

##### 2.1. FastICA Algorithm

FastICA is a fixed-point iterative algorithm, minimizing mutual information between estimated components [14]. Separation of independent components is accomplished when the maximum of non-Gaussianity is attained. There are different kinds of FastICA, including those based on kurtosis, based on the maximum likelihood, and based on the maximum negentropy (MNE), and so forth. In this paper, the MNE-based FastICA algorithm was used, which took the maximization of negentropy as a search direction and extracted each independent source signal in turn.

Before using the FastICA algorithm, the observed signal was centralized and whitened. The mean removal process was conducted to subtract the mean vector of the signal from the observed signals so that the observed signal became zero mean, simplifying the FastICA algorithm. The observation signal was whitened using the principal component analysis (PCA) whitening algorithm so that the components after the whitening were uncorrelated. The purpose of the FastICA algorithm based on fixed-point iterative structure was to make have the maximum non-Gaussianity, where was a row of the separation matrix . The objective function was set as where was the expectation operator and was a Gaussian random variable with zero mean and unit variance. It was assumed that also had zero mean and unit variance. was a nonquadratic function.

According to the Kuhn-Tucker condition, the optimization of could be obtained by (2) under the constraint of :where was a constant and could be obtained by , where was the initial value of and was a nonlinear function, which was the derivative of . We chose the nonlinear function . The Newton iterative method was employed to solve (2). The left part of (2) was denoted as , and the Jacobian matrix was

In order to simplify the computation of the inverse of the matrix, (3) was approximated. Because the data were whitened, (3) could be simplified as

The Jacobin matrix was a diagonal matrix, and its inverse matrix could be simply calculated. Similarly, replace the value of with the current value of for the constant . Therefore, we could obtain the approximated Newton iterative formula as follows:where and represented the updated value of . In order to improve the stability of the algorithm, was normalized by after an iteration. To simplify (5), we obtained an iterative formula of simplified FastICA algorithm:

##### 2.2. FastICA Algorithm Improvement

###### 2.2.1. Selection of Overrelaxation Factor

In order to address the problem that FastICA is sensitive to the initial weight vector, we introduced an overrelaxation factor into the iterative algorithm. If is guaranteed to have a decreasing property for a given norm, (7) was satisfied:that is,

By introducing the overrelaxation factor , it was guaranteed that can enter the convergence area of the Newton iteration algorithm from a certain value of , so as to ensure that the algorithm can achieve the convergence effect under any circumstance.

There are many methods for choosing the overrelaxation factors, such as the golden section method and step-by-step experimental method. In this paper, we used step-by-step experiment to obtain the optimal overrelaxation factor. The value of is , where . All the values of that satisfied (9) were recorded. Among these values of , the one that satisfied (8) was found and denoted as . Then, would be the optimal overrelaxation factor .where .

The selection of the overrelaxation factor included the following steps.

*Step 1. *Set the initial value , , and create an empty vector .

*Step 2. *.

*Step 3. *If the value of satisfies (9), then is used to calculate the target function TF, , and the result of TF is added to .

*Step 4. *.

*Step 5. *If , repeat Steps 2–4. Otherwise, the value of corresponding to the smallest value of TF in the vector is found and denoted as . Then, is the optimal overrelaxation factor .

###### 2.2.2. Improved FastICA

The convergence performance of the FastICA algorithm was affected by the initial weight vector. Since the initial weight vector of the FastICA algorithm was selected randomly, the efficiency of each iteration was different due to different initial weight vector. The obtained independent components would also be slightly different. Although it was also possible to select the principal component obtained during the whitening process as the initial weight vector, the algorithm easily converged to the initial value of the whitening. The FastICA algorithm should relax the initial weight vector requirements, that is, to achieve convergence in a wide range.

In order to improve the FastICA algorithm with respect to the initial weight vector , Xu et al. [15] proposed the incorporation of relaxation factor (low relaxation factor ) into the FastICA iteration. Because the negative gradient direction was the fastest declining direction, the gradient value was chosen as the optimal relaxation factor in [15]. However, the computational of the relaxation factor was of high complexity [15]. In this paper, the overrelaxation factor () was introduced into the FastICA iteration to deal with the initial weight vector generated randomly, and the requirement of the initial weight vector of the FastICA algorithm was relaxed. Compared with the low relaxation iteration [15], the proposed overrelaxation iteration converged much faster and was of low computational complexity.

In this paper, we introduced the overrelaxation factor into the iteration of FastICA algorithm. The value of was calculated according to Section 2.2.1. The iteration of became as follows:

##### 2.3. Fetal ECG Extraction by Improved FastICA

To achieve long-term monitoring of pregnant women during the perinatal period, noninvasive abdominal methods were preferred to collect mixed signals from the maternal abdomen. The collected mixed signals contained not only maternal ECG and fetal ECG, but also some noises. Before using the improved FastICA algorithm to extract fetal ECG, the observed signals were preprocessed. The third-order low-pass Butterworth filter was used to estimate the baseline signal of each channel, and the low-pass filter cut-off frequency was 5 Hz. The baseline drift removal signal was obtained by subtracting the baseline signal from the observed signal. Then, the processed signal was centralized and whitened to make the signal uncorrelated, thus reducing the complexity of the algorithm. The improved FastICA was used to separate each component of the source signal. As the amplitude of the maternal ECG signal is 2–10 times that of the fetal ECG signal, the extracted fetal ECG signal usually contained some maternal ECG interference. In order to remove the maternal ECG interference, the maternal ECG signal channel was selected from the separated signal components, and R-wave detection was conducted on the maternal ECG signal. Then, the SVD method was used to remove the maternal ECG components in each channel to obtain clean fetal ECG without maternal ECG interference [16].

Fetal ECG extraction steps were as follows.

*Step 1. *Remove baseline drift and centralize the observed signal , .

*Step 2. *.

*Step 3. *Initialize random vector , and set the error of convergence .

*Step 4. *According to Section 2.2.1, calculate overrelaxation factor .

*Step 5. *According to (10), adjust .

*Step 6. *Calculate , and normalize it.

*Step 7. *If , the algorithm does not converge; return to Step 5 until the algorithm converges.

*Step 8. *Using the resulting separation matrix , all the estimates of the source signal are obtained. The maternal ECG component is selected, and R-wave detection is performed.

*Step 9. *Use SVD to remove the maternal ECG components in each channel to get a clean fetal ECG signal.

#### 3. Fetal ECG Extraction Algorithm Evaluation

The signal-to-noise ratio (SNR) and statistical evaluation were used to compare the performances of the proposed algorithm and the conventional FastICA algorithm in terms of fetal ECG extraction performance.

##### 3.1. Signal-to-Noise Ratio Evaluation

The signal-to-noise ratio (SNR) based on eigenvalues and the SNR based on cross-correlation coefficients proposed by Outram [17] were used to evaluate the performance of the proposed algorithm and the conventional FastICA algorithm.

First, R-wave detection of the extracted fetal ECG was performed. Each R-wave was taken as a reference to cut off signal segments whose length is . Each signal segment contained a complete QRS wave. Then, the signal segments were used to form a matrix , and the normalization of each column of the matrix with zero mean and unit variance was performed. The SNR based on eigenvalues was defined as follows:where was the eigenvalues of the matrix and was the maximum of the eigenvalues of the matrix .

The SNR based on cross-correlation coefficients is defined in where , where is the fetal ECG signal that contains intact QRS waves.

##### 3.2. Statistical Evaluation

The performance of the methods was evaluated on the total length of fetal QRS signals, using sensitivity (Sens), positive predictive accuracy (PPA) [18], and harmonic mean () [19]:where TP, FP, and FN are the number of true positive (correct detection of fetal ECG QRS complexes), false positive (falsely detected absence of R-peak), and false-negative (false fetal ECG QRS complex detection), respectively. is the overall probability that the fetal ECG QRS complex is correctly detected and can be used as a measure of the accuracy of the proposed method. The Pan and Tompkins algorithm [20] was used to detect the fetal QRS in this work.

#### 4. Results

The improved FastICA algorithm was applied to fetal ECG signal extraction. The algorithm was validated using the clinical database Abdominal and Direct Fetal Electrocardiogram Database (ADFECGDB) [21]. The data were collected from five pregnant women of childbirth at 38–41 weeks of gestation. Each record contains four signals obtained from the maternal abdomen and one acquired directly from the fetal head at a sampling rate of 1 kHz and a collection time of 5 minutes.

Figure 1 shows the original record of the maternal abdomen signal database named r01. Figures 2 and 3 represent the residual signals after the cancelling of the maternal ECG component. Figure 2 shows the four ECG components of r01 extracted by the improved FastICA algorithm, where S1 is the extracted fetal ECG signal. Figure 3 shows the four ECG components of r01 extracted by the conventional FastICA algorithm, where S2 is the extracted fetal ECG signal. The amplitudes of R-waves (denoted by black circles in S2 of Figure 3) are small, which may render the R-waves undetected.