Fetal ECG Extraction from Abdominal Signals: A Review on Suppression of Fundamental Power Line Interference Component and Its Harmonics

Interference of power line (PLI) (fundamental frequency and its harmonics) is usually present in biopotential measurements. Despite all countermeasures, the PLI still corrupts physiological signals, for example, electromyograms (EMG), electroencephalograms (EEG), and electrocardiograms (ECG). When analyzing the fetal ECG (fECG) recorded on the maternal abdomen, the PLI represents a particular strong noise component, being sometimes 10 times greater than the fECG signal, and thus impairing the extraction of any useful information regarding the fetal health state. Many signal processing methods for cancelling the PLI from biopotentials are available in the literature. In this review study, six different principles are analyzed and discussed, and their performance is evaluated on simulated data (three different scenarios), based on five quantitative performance indices.


Introduction
The fetal heart rate (fHR) and the morphological analysis of the fetal electrocardiogram (fECG) are two of the most important tools used nowadays in clinical investigations to examine the health state of the fetus during pregnancy. The fHR is the mostly used parameter in fetal monitoring, since 1818 [1]. While the fHR track shows a predictive value of almost 99% for the fetal well being investigation, an abnormal fHR has a predictive value of only 50%. Hence, it provides relatively poor specificity in detecting the fetal distress [2]. Additional information about the fetal well being can be obtained by analyzing the morphology of the fECG signal, which was recently introduced in clinical practice for fetal monitoring. Its clinical relevance was demonstrated by a series of clinical studies [3], randomized controlled trials [4][5][6][7][8] and prospective observational studies [9][10][11][12][13][14][15][16][17][18], which prove that clinical fetal monitoring based on both fHR and fECG morphology analysis, especially the ST waveform analysis, leads to the reduction in the number of operative vaginal deliveries, smaller rate of metabolic acidosis at birth, less blood samples performed during labor, and fetal morbidity reduction.
The standard procedure to record the fHR is the cardiotocography (CTG), sometimes known as electronic fetal monitoring [19]. When necessary to investigate both the instantaneous fHR and the fECG morphology, an invasive fetal monitoring method that uses a wire electrode attached to the fetal scalp [20], after the membrane rupture, is preferred. However, both methods have important drawbacks: (i) the fHR obtained via CTG has the potential problems of reliability and accuracy [21,22]; in addition, the beat-tobeat variability of fHR is not present in the CTG traces [23,24]; hence, rapid variations of the fHR cannot be detected; (ii) the second recording technique is invasive [20]; thus, it can put the life of both the mother and the fetus in danger (e.g., possible infections can lead to different complications). An alternative method to obtain the instantaneous fHR and the fECG morphology is the abdominal recording of the fECG which considers an array of electrodes placed on the maternal abdomen. This recording procedure overcomes the main drawbacks of the methods used in clinical routine for fetal monitoring. However, the limitation of this technique is the very low signal-to-noise-ratio (SNR) of the available recorded fECG. This is mainly due to the fact that the fECG signal is generated by a small source (fetus heart). In addition, it has to propagate through different attenuating media to reach the maternal belly surface. Hence, the fECG signals contained in abdominal signals (ADSs) provide an amplitude of about 10 V [24] which becomes still smaller around 28th until 32nd weeks of gestational age due to the appearance of the insulating layer called vernix caseosa.
Furthermore, the signal of interest, that is, the fECG, is only one (weak) component of the ADS mixture; other (disturbing) signals with higher power that also exist are the electromyogram (EMG) of the abdominal muscles, the electrohysterogram (EHG), the maternal ECG (mECG), the baseline wander basically due to the maternal respiration, and the power line interference (PLI). Among them, the PLI, with the fundamental PLI component of 50 Hz/60 Hz, and its harmonics is one of the most disturbing noise sources, because it can reach amplitudes much greater than the abdominal fECG signal, making its analysis almost impossible (see Figure 1).
The PLI is determined by the power supply network, and its appearance in the abdominal recordings is explained by (i) the electrostatic induction and parasitic capacitance coupling between the body and the ground; and (ii) the electromagnetic induction through loops of the recording cables, where a time-varying magnetic field generates a voltage proportional to the loop area (depending on its orientation) and to the strength of the magnetic field [25]. If the cables are twisted, the induced voltage is reduced [25], but still significant for the fECG analysis.
The fundamental PLI is definitely a problem in fECG analysis, and its harmonics, usually present, make the PLI cancelling problem even more complex. The harmonics are usually generated by connected nonlinear loads: neon lamps, TVs, microwaves ovens, fridges, air conditioning devices, computers, and basically almost any power electronics device connected to a single-phase distribution system. The disturbing sources are in fact the rectifiers and semiconductor switches present in almost all of these nonlinear loads which introduce distortions in the power supply waveforms [26,27]. Surprisingly, the 3rd harmonic, that is, 150 Hz/180 Hz, is the most powerful PLI harmonic [28].
Although there are many practical solutions to reduce the PLI, for example, the cable twisting and shielding, the use of differential recording involving an instrumentation amplifier with high common mode rejection ratio (CMRR) at power line frequency, and the proper skin preparation to reduce the electrode imbalance, the PLI still affects the biopotential measurements. For example, the advantage of high CMRR instrumentation amplifiers is limited in real life applications, since a slight imbalance in the electrode-skin impedance leads to the divider effect [25]; thus, the PLI is partly transformed into a differential signal [25,29,30] which passes through the amplifier together with the signal of interest, when no notch filter is available in the amplifier circuit to suppress the PLI components. However, the spectrum of the ECG signal for neonates ranges between 0.01 Hz and 250 Hz [31], while the abdominal fECG frequencies are supposed to be up to 500 Hz; thus, such a notch filter affects also the signal of interest. Finally, any PLI disturbance, even if much attenuated, impairs the morphological analysis of the abdominal fECG due to its very low amplitude.
From this large variety of PLI suppression methods, six representative approaches are selected in this review study: digital notch filters (DNF), adaptive filters (AF), Hilbert Huang transform (HHT), wavelet transform (WT), blind source separation (BSS), and neural networks (NN); they are briefly described in the following sections which discuss their advantages and disadvantages. A recent and representative algorithm is implemented for each PLI cancelling approach and is evaluated using simulated data. The signal processing method with the best performance in PLI reduction, considering the minimal distortion of the original signal as the evaluation criterion, is identified.

PLI Cancellation Using Digital Fixed Notch Filters.
Notch filters are used to cancel narrow band interferences, one common application being the PLI suppression in biopotential measurements. Digital fixed notch filters can be designed to remove multiple frequencies, having the advantage of being designed to remove the fundamental frequency and also its harmonics (multiple-notch filters or comb filters [34,59]).
However, the main problem of multiple-notch filters, when used for cancelling the PLI signal from fECG signals, is the fECG and PLI spectral overlapping. Thus, the notch filter should have a very narrow bandwidth in order to suppress mainly the 50 Hz and its harmonic components, and not the useful information contained in the fECG spectrum. But this requirement comes into conflict with the fact that actually the real PLI signal does not have a fixed fundamental frequency, but rather a frequency that varies around the value of 50 Hz which requires a wide bandwidth of the multiplenotch filter. Moreover, the transient time introduced by the notch filter can be too long, in which case the fECG can be considerably distorted. Hamilton [60] has investigated the effect of the transient time of the notch filters, which increases much when the bandwidth is decreased. They observed the ringing effect appearing near the QRS complex and ST segment when narrow bandwidth notch filters are used. This distortion decreases when the transient time increases [60,61]. However, in practice, a long transient time reduces the capacity of the filter to track the noise level changes [60]. The latest available international standards for ECG acquisition, American National Standard Association for the Advancement of Medical Instrumentation (AAMI) [62] and International Standard IEC 60601 [63], do not specify any requirements for the transition band of the notch filters. The only specification present in the IEC 60601 standard states that "notch filters for line frequency interference suppression shall not introduce on the ECG record more than 25 V peak ringing noise" [63]. Nevertheless, the notch filter should have a short transient time, minimal distortion, and very narrow bandwidth.
Pei and Tseng [64] propose a method to decrease the transient time of multiple IIR notch filters. This technique uses the vector projection in order to find better initial values for the IIR notch filter. A more recent paper [32] reports better results in suppressing the transient time than the ones obtained when applying the method introduced in [64]. Piskorowski proposes a time-variant multiple-notch IIR filter. The transient time is reduced by varying the pole radius with time and thus the filter is able to cancel the fixed frequencies PLI components as fast as possible, with no longterm selectivity impairment [32]. This type of filter should not be confused with adaptive filters, which are varying the notch central frequency, as is explained in the next section.
The time-varying multiple-notch IIR filter (TVMNF) proposed in [32] is chosen as representative for this category of power line PLI suppression methods. Thus, the general transfer function of the multiple-notch filter is where is the number of notches, Ω is the central frequency of the notch, is the order of the harmonics, ( ) is a symmetrical polynomial, and is the pole radius. The selectivity of the filter increases when is increased, but this results also in a longer transition. Thus, the transition time is influenced by the radius and in order to improve the time domain filter response, the is varying in time. The difference equation of the IIR multiple-notch filter with a time-varying parameter is where the variation of the pole radius varies is described by [65] with variation range = (0)/ and = lim → ∞ ( ); V includes the exponential variation of ( ) in (3), and is the sampling frequency [32]. The value chosen for is critical; a very high value generates narrow notches, but their transition time is increased, while a small value leads to a decreased transition time, but less selective notches are obtained, which results in filtering out important fECG frequency components. Thus, ( ) has an exponential variation, from an initial value, (0), to the desired one, [32].
The filter is implemented with one notch frequency and three notch frequencies, respectively, choosing = 0.9 and V = 2, as suggested in [32]. In Figure 2 a block diagram is used to describe the TVMNF algorithm.

Adaptive Filtering in PLI Cancellation.
As previously described, fixed notch filters have the main drawback that the central frequency of the notch cannot be modified. This makes the PLI cancellation difficult when the PLI fundamental frequency has slight variations, which is often the case in real applications. To overcome this problem adaptive filters are introduced, which have the ability to adapt their notch frequency, tracking the changes in the PLI fundamental frequency. The first adaptive filter was introduced by Widrow et al., [41] and according to the review [60], this type of filters introduces less distortion than the fixed notch filters, having also a shorter transition time. The adaptive noise cancellers assume simultaneous recording of the noise source by an additional channel, that is, the reference signal. The reference is supposed to be uncorrelated with the signal of interest, but correlated with the disturbing signal. The filter adapts its parameter in order to make the reference signal as close as possible to the noise that disturbs the signal of interest, by minimizing the output error, considering the least mean square (LMS) criterion.
Many adaptive filters that suppress the PLI in ECG recordings are available in the literature. Wan et al. [67] propose a LMS adaptive algorithm with variable step size and suggest that faster convergence rate and smaller mean square error are obtained, as compared to the traditional approach. Costa and Tavares [68] come with an improvement of the basic adaptive canceller by providing also harmonics' suppression, with a minimal increase in computational complexity. Thus, the algorithm is suitable for low cost acquisition systems. Liangling et al. [42] exploit the capability of the adaptive filters to cancel the PLI component and baseline wander from ECG tracings [67].
However, these adaptive filters have a practical drawback. They all need an additional recording of the disturbing signal, that is, the reference signal. That is not always possible (e.g., in abdominal fECG recording using portable devices). Thus, adaptive filters with no reference signal are more suitable for practical applications. Ziarani and Konrad [69] propose a filter which is able to estimate the amplitude, the phase, and the frequency of the PLI components. An improved version of these adaptive filters is presented by Martens et al. [39] (Improved Adaptive Canceller-IAC), who developed an algorithm able to suppress both the PLI fundamental frequency and its harmonics. The algorithm considers the fact that large QRS amplitudes can distort the estimation of the PLI components and produce large transient segments and thus the adaptive process is blocked in such situations. The algorithm is compared with the classical adaptive filter and two notch filters with large and narrow bandwidth, respectively. The algorithm proposed by Martens et al. outperforms the other techniques, showing a stable behavior even in the worst conditions. This algorithm is chosen as representative for the adaptive filtering PLI cancellation approach and is implemented in the current study.
In Figure 3 a general scheme of the adaptive notch filter is presented.

Blind Source Separation Applied in PLI Suppression.
This PLI cancelling method is based on a completely different concept than digital filters, considering the statistical properties of a mixture of signals. Each signal source is extracted from the mixture, as long as they satisfy some conditions. Many approaches to estimate the ICA parameters exist: maximization of nongaussianity [65,70], maximum likelihood estimation [71,72], tensorial methods [73], and so forth. Different research groups use ICA algorithms to extract the fECG from abdominal recorded signals: Zarzoso and Nandi 2001 [74], Vrins et al. [75], Sameni et al. [76,77], Lathauwer et al. [73], Camargo-Olivares et al. [78], Cardoso [79], and so forth. All ICA studies report that the used ICA methods perform reasonably well in extracting the fECG signals from ADS, that is, separating the abdominal fECG from other types of noise signals, including the PLI component, present in the ADSs. It should be noticed that the interpretation of the abdominal fECG obtained via BSS methods does not have a clear physical explanation, since ICA does not take into account the position of the electrodes and other physical parameters.
The algorithm proposed in [65], FastICA, is chosen as representative for this PLI cancelling approach and is considered in the current study. In Figure 4 the general block diagram of the ICA concept is depicted.

Hilbert Huang Transform Applied in PLI Reduction.
The Hilbert Huang Transform (HHT) is a powerful method for analyzing nonlinear and nonstationary time series and it was introduced by Huang et al. [80]. The method overcomes the shortcomings of the Fourier transform which is valid just for stationary time series. It is proven that Fourier transform offers a wrong energy-frequency distribution with no physical meaning when applied on nonstationary time series [80]. Taking into account that most real signals, and especially the biopotentials, are nonstationary, the HHT is suitable for their analysis.
The method has two steps: (i) generation of Intrinsic Mode Functions (IMFs) through Empirical Mode Decomposition (EMD); (ii) Hilbert analysis [80]. IMFs are fully generated from the data set and must satisfy two conditions: (a) the number of zero crossing and the number of extrema have to be equal or to differ at maximum by one; (b) the mean value of the envelope including the local maxima and the envelope defined by the local minima is zero at any point.
The generation of IMFs is fully data driven and is obtained by decomposing the time series (the process is called "sifting") using the EMD; high frequency components are decomposed into the first IMFs, while the low frequency components are found in the higher order IMFs.
This method is reported to cancel the PLI signal from the ECG [81][82][83]. The basic idea is to discard the IMFs which contain noise and to reconstruct the signal from the remaining IMFs. The main problem is to identify which IMFs contain just the PLI signal. In conventional EMD the first IMF is considered to be related to the PLI signal since it contains the higher frequency and is therefore discarded, which is not necessarily true, since the sifting process is not perfect, allowing high frequency ECG components in the first IMF. Moreover, if the PLI contains harmonics, then the number of IMFs containing PLI components is higher, since the harmonics are decomposed into different IMFs. Pal and Mitra [66] propose an algorithm that identifies the IMFs containing the PLI components. It computes the IMFs' cumulative mean and their powers, and using these parameters identifies which IMFs contain PLI components, on a threshold basis. However, the algorithm is not robust when the power of the PLI signal is high as compared to the signal of interest, the fECG. Therefore, the current study improves the original HHT based PLI cancelling method, by a more precise identification of the noise IMFs being obtained. The Hilbert transform is computed for each IMF and the instantaneous frequency is derived thereafter. The IMFs containing the PLI components are then discarded; as already mentioned, these IMFs can still contain information about the high frequency ECG components, that is, the QRS complex. In order to recover the QRS complex, the IMFs with QRS complex are firstly detected, based on the instantaneous frequency, and then the algorithm described by Kabir and Shahnaz [47] is applied as follows: (a) the QRS complex boundaries are identified: (a.1) the peak locations are detected; (a.2) two nearest local minima, located on both sides of peak are found; (a.3) one zero-crossing point on the left-hand side of its left minimum and the other one on the righthand side of the right minimum are detected. The boundaries of the QRS complex are assumed to be between these two points; (b) a Tukey window centered on the peaks, which spans to cover the QRS complex, is applied, that is, multiplies the selected IMF. This window offers a flat gain at the wave and decreases gradually to zero ensuring a smooth transition with minimal distortion. Thus, the information of the QRS complex is preserved, allowing the estimation of the fECG. In Figure 5 the block diagram of the algorithm is depicted.

PLI Cancellation by Applying Neural Networks.
The classical application of NN in cardiac signals processing is the classification of ECG signals, pattern recognition [84,85], and fECG extraction from ADS [86]. Methods for fECG SNR improvement are described in [87] where a Functional Link Artificial Neural Network (FLANN) is proposed to remove the Gaussian and baseline wander noise. Zhang and Benveniste [88] and Poungponsri and Yu [89] use NN combined with Wavelet transform for better results. However, in a recent article Poungponsri and Yu [51] come with an improvement of the method in [89] and the algorithm is tested also on PLI cancellation (Wavelet Neural Network-WNN). The NN based adaptive filtering approach proposed in [51] for ECG signal noise reduction removes the PLI signal by applying firstly the wavelet decomposition. The wavelet coefficients are further applied to a neural network trained to reconstruct the denoised ECG (see Figure 6). The algorithm was initially developed by the authors to cancel all the noise sources overlapping the ECG signal recorded at 360 Hz (the PLI fundamental frequency is 60 Hz).
The algorithm presented by Suranai et al. is chosen as representative for this category of methods and its steps are as follows: (a) A real signal, not affected by the PLI, is applied to a feed-forward NN with 64 inputs, 2 hidden layers, and 1 output corresponding to the denoised ECG signal. The inputs of the NN are obtained by applying the Wavelet Packet Decomposition (WPD) using the Debauchies 4 wavelet; the thresholding (soft threshold) is then applied to reduce the high frequency noise. The resulting 64 wavelet coefficients are the inputs of the NN having 56 hidden neurons on the first hidden layer and 12 neurons on the second one. The hyperbolic tangent activation function is used for all the neurons: The first 3350 samples of the "clean" selected signal resampled at 300 Hz (record 220 from the Massachusetts Institute of Technology (MIT-BIH) database [90]) are used for this training phase involving the back propagation algorithm (4000, for the raw signal, recorded at 360 Hz). In order to allow the NN to remove the PLI from signals recorded with different amplifier gains, the signals are initially normalized as follows: where max and min are the maximum, and, respectively, the minimum values of the signal. The training is then continued for the same segment, affected by the PLI, so that the NN learn to cancel the PLI noise.
(b) The noisy signals are then applied to the network (i.e., their wavelet coefficients, after thresholding); the denoised signal is considered as the output of the network.
It should be noticed that the algorithm presented above is adapted in this study to remove the PLI signal with the fundamental frequency of 50 Hz and for a sampling frequency of 1000 Hz.

PLI Reduction by Applying the Wavelet Transform.
In the recent years, discrete wavelet transforms and thresholding techniques have been used for ECG denoising [91]. Wavelet based noise cancelling techniques became very popular because they are able to decompose the signal into timefrequency domain which is appropriate for the analysis of nonstationary signals. It is reported in the literature that discrete wavelet transform does not introduce any artificial information to the original signals; the threshold is generated based on the attributes extracted from the signal [92,93].
The main problem is the identification of the mother wavelet, the level of decomposition, and the optimal threshold. Garg et al. [94] compare different mother wavelet functions for ECG denoising and conclude that the recovery of the ECG with minimal artifacts is obtained when using Sym10 decomposition at level 5 and hard shrinkage function with either rigorous SURE or heuristic SURE threshold [94]. More recently, Galiana-Merino et al. [43] use the discrete stationary wavelet packet transform (DSWPT) to suppress the PLI signal and its harmonics from electromyographic (EMG) signals (the DSWPT algorithm). Basically it is a shift invariant transform to isolate the 50 Hz and its harmonics, with the disturbing sine signals being reconstructed thereafter using the DSWPT coefficients. In Figure 7 the steps of the DSWPT algorithm are described as follows. (e) The refined sine disturbances are subtracted from the signal that has to be denoised.
The main results regarding the performance of the selected algorithms are summarized in Table 1.

Data Simulation and Performance Measurements.
The simulated data, used to quantitatively estimate the performance of the proposed algorithms are generated in two steps. Firstly, the fECG is simulated using the dynamic model introduced in [95,96]:̇= − , where is the angular velocity of the time vector as it moves around the limit circle (representing the period ), = 1 − √ 2 + 2 , = tan( / ), Δ = − , contains the amplitudes of the peaks, contains the width of each peak, and are the angles which specify the -, -, -, -,waves/peaks. Secondly, the PLI components, simulated as sinusoids, are added. Usually the PLI fundamental component is supposed to be constant. However, there are some deviations from the fundamental frequency in real applications, mainly due to unstable power sources. Thus, the PLI components can exhibit significant frequency deviation, up to 3% [30,[97][98][99][100] (the deviation differs from country to country, depending on the available power supply technologies). Three data sets are therefore constructed, considering three PLI scenarios (7): (i) the PLI contains just the fundamental power line interference component, of 50 Hz, (ii) the PLI includes both the power line fundamental frequency and its 3rd harmonic (150 Hz), and (iii) the PLI is a sinusoid whose frequency slightly varies in time around 50 Hz. Consider where = 20%. It is clear that the worst scenario is the third one, assuming that the power line fundamental frequency is time-varying.
For each scenario, the SNR defined by (8) is varied; five noise levels are considered: −2 dB, 0 dB, 2 dB, 4 dB, and 5 dB: For the qualitative evaluation of the implemented algorithms, the following performance indices are considered as follows.
(a) Normalized root mean square error, expressed in percentage: (b) Noise retention, expressed in percentage: where is the power of the signal computed with orig fECG = 10 * log 10 ∑ =1 orig fECG( ) 2 .
(c) SNR improvement [47]: where signal is the input signal containing both the fECG and the PLI.
(d) Cross-correlation coefficient, considering the original and the denoised fECG signal:

Results
The results obtained when applying the five selected algorithms are organized as follows. (i) For each performance index a table is constructed. The columns correspond to the evaluated algorithms and the rows to the scenarios, assuming that the SNR is −2, 0, 2, 4, and 5 dB. The best result obtained for the performance index is emphasized in bold for each case and for each scenario. The scenarios for which some algorithms are not working by principle are represented by empty gray cells. (ii) The performance indices for scenario 2, when the SNR is −2 dB, are illustrated for each algorithm (see Figures 8,9,10,and 11

Discussion
In Tables 2, 3, 4, and 5, it can be observed that the worst overall performance is obtained when the WNN is applied. The algorithm is able to reduce the noise if its level is very low, but the QRS complex, containing high frequencies, is disturbed, impairing the fECG morphology analysis. When the signal is hidden by the noise, that is, low level of SNR, the method fails to extract the denoised signal, which can be explained by the thresholding step. The main advantage of the algorithm is the computation time of the testing (denoising) phase, assuming that the neural network is already trained. Because of the overall bad performance the algorithm is excluded from the following discussion. The TVMNF method proposed by [32] shows good results for scenarios 1 and 2 having a RMSD of 2.11%, Figure 8, and a SNR improvement, SRN imp = 76 dB for the scenario 2 (see Figure 11), SNR = −2 dB (worst scenario). No ringing effect is observed near the QRS complexes. However, despite         the good results, the main drawback is that the filter cannot be applied when the fundamental frequency is varying, as expected in real applications. Thus, this type of filter is of limited usage.
As expected, the adaptive methods are able to follow the changes in the frequencies of the PLI components; hence, they can be successfully applied in scenario 3. According to the computed performance indices, the IAC is able to obtain good estimates of the fECG signal when the fundamental frequency is fixed, even when the harmonics are present. However, for the worst SNR, that is, −2 dB, the obtained performance is slightly worse as compared to the performance of the other algorithms (see . For scenario 3 the algorithm has the worst performance (RMSD = 655% and a noise retention factor of NR = 95%). This is due to the fact that the algorithm needs very long adaptation time (approximately 10 s) when step variations in the PLI fundamental frequency appear. Thus, IAC is of limited use in real applications.
The HHT method implemented in this study shows also good results. However, at small SNR values, small oscillations appear near the QRS complex, suggesting that some very low PLI components can be found in the low order IMFs. Moreover, when the method is applied in scenario 3, outliers appear when the PLI fundamental frequency is changing. The main advantage of this method is that it is suitable for nonstationary signals like the biopotentials and that it is fully data driven; that is, no a priori knowledge is necessary. The main disadvantage is that the decomposition does not fully separate the oscillations; thus, some useful information can be found in the IMFs containing the PLI components; in addition, it does not have, by now, a complete mathematical evaluation. However, recent papers present some improvements to the basic method claiming a better decomposition (e.g., Ensemble Empirical Mode Decomposition (EEMD) [101] and Complete Ensemble Empirical Mode Decomposition [102]).
The BSS algorithm is able to separate the fECG from the PLI showing the same performance no matter how the SNR is varying, because by principle it exploits the independence between any two signal sources. However, it should be noted that the comparison with the other method is not quite fair in the current study; in order to have the same simulations, the same signals used for the other algorithms were fed to the Fas-tICA. In real application, the input of ICA algorithm is signals obtained from different channels, meaning different fECG waveforms, which can lower the performance of ICA in fECG extraction. Moreover, the physical relevance of the fECG independent components obtained when applying ICA is still a subject of discussion among researchers. Thus, the method has limited usage in real application, but it can be successfully combined with the adaptive filtering techniques, improving the estimation provided by the PLI reference block.
As theoretically expected, the best performance is obtained when using the DSWPT, if the PLI is stationary and includes exactly the 50 Hz and its harmonics. However, the method does not work in the worst scenario. Since the algorithm estimates the PLI interference assuming that the power line frequency is 50 Hz (in step 2, templates of the sine disturbances are constructed, averaging the segmented signal, using a window of 20 samples), it is expected that the algorithms fail in cancelling the PLI when the frequency is more or less different from 50 Hz, or even worse, when the power line frequency is varying, which is supported by the obtained results.

Conclusions
In this paper a review of PLI cancelling methods applied in fECG signal processing is proposed, revealing the main concepts provided in the literature for suppressing the 50 Hz/60 Hz component and its harmonics from biosignals. The selected algorithms are quantitatively analyzed, using different performance criteria and practical considerations are provided when discussing the PLI cancelling from abdominal fECGs. Three sets of simulated data are constructed and used in the quantitative evaluation of the algorithms, considering the 50 Hz PLI fundamental frequency, the 50 Hz combined with the 150 Hz PLI component, and a varying PLI fundamental frequency. The quantitative performance is monitored using five different indices, corresponding to different SNRs.
While some methods, like WNN, show very bad performances, most of the algorithms have good results, especially in scenarios 1 and 2. The DSWPT has the best performance in scenarios 1 and 2, as depicted in the Figures 8-11, but its main drawback is that it is not suitable for the most realistic scenario, scenario 3. The HHT based algorithm shows the best performance overall, considering the implemented scenarios. Thus, further studies should concentrate on exploiting the capabilities of the HHT method.