A Novel Technique for Fetal ECG Extraction Using Single-Channel Abdominal Recording

Non-invasive fetal electrocardiograms (FECGs) are an alternative method to standard means of fetal monitoring which permit long-term continual monitoring. However, in abdominal recording, the FECG amplitude is weak in the temporal domain and overlaps with the maternal electrocardiogram (MECG) in the spectral domain. Research in the area of non-invasive separations of FECG from abdominal electrocardiograms (AECGs) is in its infancy and several studies are currently focusing on this area. An adaptive noise canceller (ANC) is commonly used for cancelling interference in cases where the reference signal only correlates with an interference signal, and not with a signal of interest. However, results from some existing studies suggest that propagation of electrocardiogram (ECG) signals from the maternal heart to the abdomen is nonlinear, hence the adaptive filter approach may fail if the thoracic and abdominal MECG lack strict waveform similarity. In this study, singular value decomposition (SVD) and smooth window (SW) techniques are combined to build a reference signal in an ANC. This is to avoid the limitation that thoracic MECGs recorded separately must be similar to abdominal MECGs in waveform. Validation of the proposed method with r01 and r07 signals from a public dataset, and a self-recorded private dataset showed that the proposed method achieved F1 scores of 99.61%, 99.28% and 98.58%, respectively for the detection of fetal QRS. Compared with four other single-channel methods, the proposed method also achieved higher accuracy values of 99.22%, 98.57% and 97.21%, respectively. The findings from this study suggest that the proposed method could potentially aid accurate extraction of FECG from MECG recordings in both clinical and commercial applications.


Introduction
Currently the most widespread techniques to monitor fetal health are Doppler ultrasound, invasive fetal electrocardiogram (FECG) monitoring, and non-invasive FECG monitoring. The monitoring method that adopts Doppler ultrasound is routinely applied during pregnancy and delivery [1]. However, the ultrasound technique has not been proven to be fully safe for the fetus [2]. In addition, the method does not allow for long-term monitoring since it is sensitive to maternal and fetal motion [3], and only provides one monitoring parameter that is the fetal heart rate (FHR), which is not sufficient for accurate thoracic MECG recording and it was used as the reference signal in ANC. The limitation associated with the lack of strict waveform similarity between the interference (MECG in AECG) and the reference signal (MECG recorded in maternal thorax) could be avoided completely in the proposed method. The method only require one abdominal recording. Therefore, the electrode configuration allows for building light and comfortable recording systems for continuous long-term fetus monitoring.
This paper is organized as follows: in Section 2, the theory and implemented detail of our proposed method are described. In Section 3, the two databases used in the current study and statistical assessment are presented. In Section 4, the theory and implementation detail of the four typical single-channel methods are presented in detail. Section 5 presents the Results while Section 6 discusses the results. Moreover, the paper is concluded in Section 7.

Basic Structure
The block diagram of the proposed method is depicted in Figure 1. The abdominal recording is considered to be AECG in this paper (the noise in the abdominal recording is ignored). The SWSVD scheme estimates the MECG of the abdominal signal that serves as input in form of the reference signal to the adaptive filter. The adaptive filter removes the MECG in the abdominal recording according to the estimated reference. For this purpose, an adaptive algorithm is built to model learning to minimize the interference signal (MECG) in abdominal recording with inputx(n) and output x (n), where x (n) matches the target signal y(n) as closely as possible, andf (n) is the target signal and the extracted FECG represents the output of the adaptive filter. the reference signal (MECG recorded in maternal thorax) could be avoided completely in the proposed method. The method only require one abdominal recording. Therefore, the electrode configuration allows for building light and comfortable recording systems for continuous long-term fetus monitoring. This paper is organized as follows: in Section 2, the theory and implemented detail of our proposed method are described. In Section 3, the two databases used in the current study and statistical assessment are presented. In Section 4, the theory and implementation detail of the four typical single-channel methods are presented in detail. Section 5 presents the Results while Section 6 discusses the results. Moreover, the paper is concluded in Section 7.

Basic Structure
The block diagram of the proposed method is depicted in Figure 1. The abdominal recording is considered to be AECG in this paper (the noise in the abdominal recording is ignored). The SWSVD scheme estimates the MECG of the abdominal signal that serves as input in form of the reference signal to the adaptive filter. The adaptive filter removes the MECG in the abdominal recording according to the estimated reference. For this purpose, an adaptive algorithm is built to model learning to minimize the interference signal (MECG) in abdominal recording with input ( ) and output ′( ), where ′( ) matches the target signal ( ) as closely as possible, and ( ) is the target signal and the extracted FECG represents the output of the adaptive filter.

Reference Formation Scheme SWSVD
A detailed description of the SWSVD scheme is presented in Figure 2. The aim of SWSVD is to create an artificial reference signal (estimated abdominal MECG) using a single channel abdominal recording. An ECG cycle basically consists of a P-wave, T-wave, and QRS complex (which can be subdivided into a separate Q-wave, R-wave, and S-wave) [29][30][31]. The T-wave, P-wave, and the segment between them are referred to as inter-QRS waveform in the current study. Basically, the estimation of the MECG from abdominal recording involves two steps which are dynamic estimation of MECG QRS segments using SVD method, and smoothening of FECG QRS in inter-QRS segments by average filter. The estimated and smoothed segments were concatenated as MECG wave which was used to build an artificial reference signal for the adaptive filter.

Reference Formation Scheme SWSVD
A detailed description of the SWSVD scheme is presented in Figure 2. The aim of SWSVD is to create an artificial reference signal (estimated abdominal MECG) using a single channel abdominal recording. An ECG cycle basically consists of a P-wave, T-wave, and QRS complex (which can be subdivided into a separate Q-wave, R-wave, and S-wave) [29][30][31]. The T-wave, P-wave, and the segment between them are referred to as inter-QRS waveform in the current study. Basically, the estimation of the MECG from abdominal recording involves two steps which are dynamic estimation of MECG QRS segments using SVD method, and smoothening of FECG QRS in inter-QRS segments by average filter. The estimated and smoothed segments were concatenated as MECG wave which was used to build an artificial reference signal for the adaptive filter. segment between them are referred to as inter-QRS waveform in the current study. Basically, the estimation of the MECG from abdominal recording involves two steps which are dynamic estimation of MECG QRS segments using SVD method, and smoothening of FECG QRS in inter-QRS segments by average filter. The estimated and smoothed segments were concatenated as MECG wave which was used to build an artificial reference signal for the adaptive filter.  The first step in estimation of MECG signals is to extract the AECG QRS complex from abdominal recordings. For this purpose, SVD method used in [32] was adopted. N preceding extracted MECG QRS complexes waves are arranged in the form of a designed matrix Y. The consecutive MECG QRS complex waveforms occupy consecutive rows of the later matrix, while the R peak lies in the same column. Hence, matrix Y is constructed as in Equation (1). Maternal R-peak was detected according the method proposed in the works of Pan and Tompkins [33]. The MECG QRS complex was extracted with fixed length before the maternal R-peak, W 1 , and after the maternal R-peak, W 2 , as shown in Figure 3: where, W is the duration of AECG QRS cycle and N is the number of AECG QRS complex cycles: The MECG QRS complex component (principle components) of design matrix Y can be obtained using SVD. Parameter N was set to 20 in this study. The MECG QRS complex in the abdominal recoding can be estimated as follows: where q is the principle component number of the abdominal recording and W j is the direction vector of the j-th principle component. An appropriate reverse transformation is performed on the data reconstruction of the estimated MECG matrix Y'. Each data segment {y(·)} is transformed to {x'(·)}. As used in some earlier studies [30,31], the duration of the QRS complex segment is set and fixed as 0.1 s. In this current study, QRS wave segment is set to be the samples between 0.04 s and 0.06 s after the R peak was detected. In the later test section, W 1 was set to be 40 sample points while W 2 was set to be 60 sample points using a test data recorded at sample frequency of 1000 Hz.
In the second step, a smoothening window method was applied to smooth the fetal component associated with the inter-QRS wave. In order to approximate the maternal inter-QRS waveform, an average filter was used to average a number of points from the inter QRS signal to produce each point in the output signal. In equation form, this can be expressed as: where x(i) is AECG,x(i) is the estimated MECG, and M(i) is the number of points used in the moving average filter. in the output signal. In equation form, this can be expressed as: (4) where x(i) is AECG, ( ) is the estimated MECG, and M(i) is the number of points used in the moving average filter. The group of points from the input signal is chosen symmetrically around the output point. Symmetrical averaging requires that M(i) be an odd number. If the sampling sequence x(i) is the sequence before point A and after point D of a previous cycle shown in Figure 3, then M(i) is set to be a constant L which is approximately equal to the length of FECG QRS complex, that is the number of sample points. The length of FECG QRS complex is about 26-61 ms when the pregnancy is between 17 to 41 weeks [34]. In this study, L is set to 30 samples in a test data with frequency of 1000 Hz.
Jitter in the reference signal introduced by mismatch in the estimated QRS and Inter-QRS will cause a significant output error. In order to avoid such error in the estimated reference signal, the length of the average filter was reduced linearly when the input signal x(i) was spanning point A to B or point D to C.
That is to say, if AECG of x(i) is in the segment between points A and B or the segment between points C and D, the M(i) is determined by the following rules: where i is the location of AECG x(i); k and k + 1 represent the two adjacent R peak locations.

Adaptive Filter
For the adaptive filter, we adopted the recursive least square (RLS) method. In contrast to the least mean squares filter that adapts its coefficients by considering only the current error value, the RLS filter considers the total error from the initially supplied signal values to the current data point. In the context of FECG component extraction, the RLS filter tends to converge faster and is usually more accurate. The operational procedures of the RLS method are described as follows: , be the input signal vector and ω(n) = [ω 1 (n), ω 2 (n), . . . , ω H (n)], be the coefficient vector of the adaptive filter. H is the number of filter coefficients and the output of the filter is given by: The observable error is determined by using Equation (7): In addition, the cost function was minimized by using the formulae presented in Equation (8): where λ ∈ [0, 1] and it denotes the forgetting factor which determines the weight of previous data when updating filter coefficients.
To minimize the cost function ε(n), partial derivatives for all the entries of the coefficient vector were taken and results were initialized to zero as shown in Equation (9): At each iteration, the filter coefficients of the RLS algorithm were updated based on the following equations: The formula in Equation (12) was derived based on Equations (10) and (11): The parameters λ and H were set to 0.99 and 20, respectively, in the study.

Database
Two different datasets were used in the current study to assess the performance of the proposed method. The first dataset was collected from an online clinical database (Abdominal and Direct Fetal ECG Database [35]) provided by the PhysioNet. The Abdominal and Direct Fetal ECG Database consist of five multichannel datasets, that is: r01, r04, r07, r08 and r10, recorded during delivery processes. Each multichannel dataset consists of four AECG recordings acquired from the maternal abdomen and one direct FECG recording acquired from the fetal head simultaneously. The configuration of the abdominal electrodes comprised of four electrodes placed around the navel, a reference electrode placed above the pubic symphysis, and a common mode reference electrode (with active-ground signal) placed on the left leg. In order to reduce the skin impedance, the areas under the Ag-AgCl electrodes (Red Dot 2271, 3M, St. Paul, MN, USA) were abraded. The first AECG recording of r01 and the fourth AECG recording of r07 are termed as public dataset and denoted as DB1 in this study.
Meanwhile, the second set of data was obtained with BIOPAC Acquisition Systems, Inc. (Goleta, CA, USA), provided with a ECG100C module at the PLA Navy Anqing Hospital, Anhui, China. The acquisition procedures includes placing the two active electrodes below the mother's navel, thus to be as closer as possible to the heart of the fetus as in Figure 4. The electrodes were of type "offset-electrode" and this was chosen to ensure high-quality recording [36]. This dataset contains a single abdominal recording that was captured from the abdomen of a pregnant woman at 40 weeks of gestation. The signal were sampled at 1000 Hz and the private database is denoted DB2.
Both the public and private datasets were used to evaluate the performance of the proposed method and the four typical single channel FECG extraction methods considered in this study as shown in Section 5.
CA, USA), provided with a ECG100C module at the PLA Navy Anqing Hospital, Anhui, China. The acquisition procedures includes placing the two active electrodes below the mother's navel, thus to be as closer as possible to the heart of the fetus as in Figure 4. The electrodes were of type "offsetelectrode" and this was chosen to ensure high-quality recording [36]. This dataset contains a single abdominal recording that was captured from the abdomen of a pregnant woman at 40 weeks of gestation. The signal were sampled at 1000 Hz and the private database is denoted DB2. Both the public and private datasets were used to evaluate the performance of the proposed method and the four typical single channel FECG extraction methods considered in this study as shown in Section 5.

Statistical Assessment
Statistical analysis was conducted to compare the performance of the proposed method with four existing typical single channel methods in terms of FECG extraction accuracy. The adapted version [37] of the Pan and Tompkins algorithm [33] was used to detect the fetal QRS in this study. The performance of the FECG QRS complex detection was assessed by beat-to-beat comparisons between the detected FECG QRS complex and the annotated. In accordance with the ANSI/AAMI guideline [38], sensitivity (Se), positive predictive value (PPV), and accuracy (ACC) metrics [39] were used for the assessment and there formulae are presented below: where TP, FP and FN are the number of true positives (correctly detected FECG QRS complex), false positives (falsely detected non-existent R peaks), and false negative (missed FECG QRS complex detections) respectively. T = TP + FN is the number of annotated FECG QRS complex. The F1 is the overall probability that the FECG QRS complex are detected correctly and could alternatively be used as a measure of the method's accuracy: The adult ECG acceptance interval between the detected and annotated locations of QRS complex was defined as 150 ms [40], and a matching interval of 50 ms was used in the FECG QRS complex detection due to higher FHR [41]. PhysioNet provides an implementation of FECG QRS detection which is available online [42]. In the implementation, fetal heartbeats were detected with an energy threshold of 0.42 and the refractory period between two R-peaks is 0.15 s.

Compared Methods
Four typical single-channel methods were implemented in order to evaluate performance of the proposed method. Procedures of the four methods are discussed briefly in this section.

Cerutti Method
The Cerutti method proposed in [43] builds a maternal template to each cycle by averaging and scaling procedures. Afterwards, the maternal templates are subtracted from each abdominal cycle, leaving residual FECG and noise. As presented in other studies [30,31], the P wave, QRS complex and T wave durations are normally about 0.2 s, 0.1 s and 0.4 s, respectively. The maternal templates are defined as signals from 0.25 s before the detected maternal R-peak and 0.45 s after the R peak. An averaging template → x is then made as mean of N preceding maternal templates synchronized on mother R peak and this is updated on every maternal cycle. The parameter N is set to 20 in our implementation of this method. Subsequently, the averaging template is scaled with a constant factor γ. Hence, the estimated maternal template is expressed as Since maternal and fetal QRS are not synchronized, and fetal QRS is larger than maternal QRS in one order amplitude, no FECG is present in → x . Therefore, the FECG is left when estimated template → x is subtracted from abdominal cycle → y .

Kanjial Method
The Kanjial method, as described in [32], use SVD to estimate principle components related to MECG in an abdominal recording. The estimated MECG is cancelled from the abdominal recording to produce FECG. The abdominal sample data is arranged as the in a design matrix X such that consecutive MECG cycles occupy consecutive rows of the matrix, and peaks of the MECG cycles lie in the same column of the matrix: where m and n are the number of MECG cycles and samples sizes in MECG cycle, respectively. SVD is performed on design matrix X and the most dominant component X mother = U ∑ V T is the estimated maternal component, forming the fetal component X f etus = X − X mother , where U ∈ pxp , V ∈ txt , U T U = I, V T V = I, ∑ = diag σ 1 , . . . , σ g . . . 0 , g = min(p,t).
The design matrix is updated every 20 MECG cycles and the duration of MECG cycles is 0.7 s (0.25 s before and 0.45 s after the detected maternal R peaks).

Suzanna Method
The Suzanna method [30] is an improvement of Cerutti method. In the estimating maternal ECG procedure, Cerutti method scale the average maternal ECG template as a whole while Suzanna method performs a separate scaling of the P wave, the QRS complex and the T wave. Refer to the paper [30], the maternal template → x is windowed with a total length of 0.7 s. In this method, QRS wave segment is the samples within 0.05 s before and after the detected R peak. P wave segment and the T wave segment is the samples within 0.2 s before and 0.4 s after the QRS wave segment. The segments vectors are assembled in a matrix M as follows: The scaling parameters γ P , γ QRS , γ T correspond to the P wave segment, T wave segment and QRS wave segment. The scaling vector: The estimated maternal template is given by The scaling vector is solved based on the search of min least squares error minS → r − → y : In the averaging procedure, the averaging maternal template is built using 20 maternal cycles and updated on every maternal cycle.

Vullings Method
The Vullings method, proposed in [3], only estimate the MECG by linear combination of preceding MECG cycles with different weight coefficients. X is the prediction matrix and synchronized on maternal R peaks.
where m and n are the number of MECG cycles and sample sizes in cycle. The MECG cycle vector → x could be estimated by the prediction matrix X and weighted coefficient vector . . λ m ] as in Equation (24): The weight coefficient vector is calculated by minimizing the mean squared error (MSE) between the estimated MECG cycle vector → x and actual MECG cycle vector → x as in Equation (25): The parameter m is set to seven maternal cycles and the duration of MECG cycles is 0.7 s, that is, 0.25 s before and 0.45 s after the detected maternal R peaks.

Results
The proposed method was implemented in Matlab 8.6.1 (R2015b) on a computer equipped with a 2 GHz Pentium ® processor. This system configuration is capable of removing MECG from the abdominal signals in real time.

Preprocess
In this study, five-order low-pass and five-order high-pass Butterworth bidirectional filters were used to remove the high frequency and low frequency signals, respectively. Figure 5 shows the different filtered results for one segment of r01 recording in DB1 using three choices of the low cut-off frequency, namely 1, 5 and 8 Hz. In all cases, the same high cut-off frequency of 100 Hz was used.

Statistical Assessment
Statistical analysis was conducted to compare the performance of the proposed method with four existing typical single channel methods in terms of FECG extraction accuracy. The adapted version [37] of the Pan and Tompkins algorithm [33] was used to detect the fetal QRS in this study. The performance of the FECG QRS complex detection was assessed by beat-to-beat comparisons between the detected FECG QRS complex and the annotated. In accordance with the ANSI/AAMI guideline [38], sensitivity (Se), positive predictive value (PPV), and accuracy (ACC) metrics [39] were used for the assessment and there formulae are presented below: The adult ECG acceptance interval between the detected and annotated locations of QRS complex was defined as 150 ms [40], and a matching interval of 50 ms was used in the FECG QRS complex detection due to higher FHR [41]. PhysioNet provides an implementation of FECG QRS

Performance Based on DB1
DB1 contains two signals: signal r01 and signal r07. The r01 signal has a duration of 5 min with a total of 642 fetal heartbeats while signal r07 also has a duration of 5 min but a total of 625 fetal heartbeats. Both signals were preprocessed with the same condition, that is, low cutoff frequency (f b ) of 8 Hz and high cutoff frequency (f h ) of 100 Hz. This was done to evaluate the performance of the proposed method and the four typical single-channel FECG extraction methods. Tables 1 and 2 show the FECG QRS detector statistical assessment for signals r01 and r07, respectively. Direct FECG channels in DB1 provide the fetal heart location accurately. Figure 6 shows a single segment performance using signal r01 of our proposed method and four typical single-channel methods with same preprocessing procedure, that is, low cutoff frequency (f b ) of 8 Hz and high cutoff frequency (f h ) of 100 Hz. proposed method and the four typical single-channel FECG extraction methods. Tables 1 and 2 show the FECG QRS detector statistical assessment for signals r01 and r07, respectively. Direct FECG channels in DB1 provide the fetal heart location accurately. Figure 6 shows a single segment performance using signal r01 of our proposed method and four typical single-channel methods with same preprocessing procedure, that is, low cutoff frequency (fb) of 8 Hz and high cutoff frequency (fh) of 100 Hz.

Performance Based on DB2
Figure 7 presents one segment performance using the private database (DB2). The duration of AECG recordings in the database is 1 min 20 s, and it contains a total of 177 fetal heartbeats. Locations of the fetal heartbeats in the dataset from DB2 were labeled by two doctors in the PLA Navy Anqing Hospital, Anhui, China. The database was used to test the performances of the proposed method and other four typical single-channel methods explained in previous section. Table 3 is the FECG QRS detected from statistical assessment based on DB2. The results presented in Table 3 and Figure 7 are based on the same preprocessing condition (low cutoff frequency f b = 7 Hz and high cutoff frequency f h = 100 Hz).  Figure 7. A segment performance using DB2 of our proposed method and other four typical single-channel methods with same preprocessed (low cutoff frequency f b = 7 Hz and high cutoff frequency f h = 100 Hz).

Discussion
The pre-processing is a crucial step in FECG extraction from abdominal recordings. The low frequency cut-off recommended is 0.05 Hz for the diagnostic electrocardiographs [44]. In the context of non-invasive FECG extraction, it is common to use a low cut-off frequency value of more than 0.05 Hz. Figure 5 shows the different abdominal ECG waves with different baseline cut-off frequency f b . It is obvious that some FECG QRS is hardly visible at very low amplitude or overlaid totally by noise in the original AECG. The FECG QRS have been enhanced by a low frequency cut-off which is greater than 5 Hz. It is reasonable that such low cut-off frequency affects the morphology of the ECG signals but this is acceptable since the objective is to analyze FHR [45]. For all approaches, higher f b (f b = 10 Hz) have been shown to perform better than a lower f b (f b = 2 Hz) with up to 3% increase with respect to F1 measure [3].
In order to maximally suppress the artifacts and minimize distorted FECG signal, preprocessing conditions with low (f b = 8 Hz) and high (f h = 100 Hz) cut-off frequencies were used to the dataset obtained from DB1. Meanwhile the dataset from DB2 was preprocessed with low (f b = 7 Hz) and high (f h = 100 Hz) cut-off frequencies. Figure 6 shows an example of a single segment with the proposed method and other typical single-channel methods by using signal r01 obtained from DB1. The performance of our FECG extraction method and the existing methods was assessed from two aspects of the segment. The first aspect shows the ability of the methods to extract the fetal QRS component. It was observed that the fetal heartbeat was successfully extracted even though the fetal QRS overlapped with the maternal QRS. In the second aspect, the ability of the methods with respect to removing the MECG QRS complex was shown. From both aspects, the proposed method was observed to have better performance than the existing methods by removing the maternal QRS from the AECG recordings. In addition, the maternal QRS component could not be cancelled completely, and this would make it difficult to accurately discriminate the fetal QRS. As shown in Figure 6, all the five single-channel FECG extraction methods performed well. However, our proposed method have the best performance in removing MECG component. From Table 1, it could be observed that our proposed method only recorded two false positive fetal heartbeats detection and three false negative fetal heartbeats detection based on signal r01 which contains a total of 642 fetal heartbeats. Suzanna method achieved an F1 score of 99.53% at low and high cut-off frequencies of f b = 8 Hz and f h = 100 Hz respectively, while our proposed method recorded a slightly higher F1 score of 99.61%. Also from Table 2, the proposed method only recorded four false positive fetal heartbeats detection and five false negative fetal heartbeats detection based on signal r07 which contains a total of 625 fetal heartbeats. Hence, our proposed method has the best performance in fetal heartbeats detection using the signals in DB1.
The DB2 just contain one channel AECG and the signal length is 1 min and 20 s. A total of 177 fetal heartbeats is contained in the DB2. Compared to DB1, the FECG QRS is interfered by a lot of noise which is very close with FECG QRS in frequency domain. In addition, the morphology of the MECG QRS complex displays time-variation and amplitude-variation largely. That is to say, signal quality of DB2 is much worse than the signals in DB1. From Figure 7, it is obvious that the MECG QRS component cannot be removed completely by the four existing single channel methods considered in this study. The residual MECG components have severely affected the identification of beat locations in the fetal ECG. The FECG QRS statistical assessment in Table 3 also indicates that the residual MECG QRS component resulted to more false positive and false negative fetal heartbeats detection in DB2. However, the F1 of our proposed method is 0.9858 while the F1 of Cerutti, Kanjial, Suzanna, and Vullings methods are 94.65%, 94.89%, 92.35% and 82.62%, respectively. It is obvious that our proposed method outperformed other single-channel FECG extraction methods with respect to DB2. One possible reason why the proposed method performed better than the existing single channel methods is that it integrates an extensive ANC structure that can effectively cancel the interferences in terms of MECG in the AECG. Therefore, when dealing with signals of poor quality (like DB2), our proposed extensive ANC structure could be effective during FECG extraction.
The four methods (Cerutti, Kanjial, Suzanna, and Vullings) that were considered for comparison in this study rely on the construction of an adaptive MECG template using the preceding MECG cycle. The estimated MECG template is subtracted coherently from AECG and FECG is left. These temporal methods have been widely used in literatures. Our work focuses on such methods because they only require a relatively low amount of dataset (a single abdominal channel) and are of low computational complexity that would allow long-term monitoring. However, the statistical analyses based on DB2 showed that the existing methods could not remove the MECG QRS component from AECG completely in situations where there are large variation in the morphology of the continuous MECG QRS complexes. The proposed extensive ANC structure performed better than the four methods in estimating MECG using continual preceding MECG QRS waveforms. Therefore, the proposed extensive ANC structure using average window technique and SVD method could effectively estimate the MECG in FECG extraction for FHR.

Conclusions
In this study, a novel framework based on ANC technique that uses a single channel abdominal recording for FECG extraction is proposed. The proposed method integrates averaging window and SVD techniques to build reference signal for the ANC based extraction framework. The proposed method could avoid the limitation of recording at least one thoracic MECG which should be linear with an abdominal MECG recording, and this is common with the classical ANC methods. Datasets consisting of abdominal and direct fetal ECG from two different sources were used to assess the performance of the proposed system and four typical single channel methods. Compared with the four single-channel methods, the proposed method performed better in terms of ACC and F1 metric. Also, the statistical results shows that the proposed method is able to produce reliable fetal heartbeat extraction using a single abdominal recording.