Computational and Mathematical Methods in Medicine

Volume 2014 (2014), Article ID 784862, 10 pages

http://dx.doi.org/10.1155/2014/784862

Research Article

New Estimators and Guidelines for Better Use of Fetal Heart Rate Estimators with Doppler Ultrasound Devices

Iulian Voicu,1,2 Sébastien Ménigot,1,2 Denis Kouamé,³ and Jean-Marc Girault1,2

¹Signal & Imaging Group, University François Rabelais of Tours, PRES Loire Valley University, UMR INSERM U930, 7 Avenue Marcel Dassault, 37200 Tours Cedex, France

²Inserm, U930, 10 Boulevard Tonnellé, BP 3223, 37032 Tours Cedex, France

³IRIT UMR 5505, Université Paul Sabatier Toulouse 3, 118 Route de Narbonne, 31062 Toulouse Cedex 9, France

Received 28 May 2013; Revised 25 October 2013; Accepted 4 November 2013; Published 29 January 2014

Academic Editor: Brynjar Karlsson

Copyright © 2014 Iulian Voicu 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.

Characterizing fetal wellbeing with a Doppler ultrasound device requires computation of a score based on fetal parameters. In order to analyze the parameters derived from the fetal heart rate correctly, an accuracy of 0.25 beats per minute is needed. Simultaneously with the lowest false negative rate and the highest sensitivity, we investigated whether various Doppler techniques ensure this accuracy. We found that the accuracy was ensured if directional Doppler signals and autocorrelation estimation were used. Our best estimator provided sensitivity of 95.5%, corresponding to an improvement of 14% compared to the standard estimator.

1. Introduction

Continuous monitoring of fetal parameters has shown their advantages in estimating fetal wellbeing [1]. According to the report of the Society for Maternal-Fetal Medicine [2], continuous fetal heart rate monitoring has reduced infant mortality. The development or the improvement of noninvasive methods dedicated to continuous fetal monitoring is therefore of major interest.

One important parameter in assessing fetal wellbeing is the variability of the fetal heart rate. This parameter, which corresponds to the variation between intervals of two consecutive heart beats, is an indicator of central nervous system development [3–5]. It characterizes fetal behavior states [6–8] and can be an indicator of further neurological evolution [9]. Variability analysis provides a good indication of fetal distress [10] and identifies fetuses with intrauterine growth retardation [11]. According to Dawes criteria [12], variability of 4 ms is a predictor of the lack of acidosis, while a value of 2.6 ms is critical for the fetus. In the normal fetal heart rate range (110–160 bpm), a time variability of 4 ms corresponds to a cardiac frequency variability of 0.81 bpm, while a time variability of 2.6 ms corresponds to a cardiac frequency variability of 0.53 bpm (the values of 0.53 bpm and 0.81 bpm are obtained as follows: bpm and bpm, resp.). Other authors [13] suggest that it is necessary to estimate the heart rate with an accuracy of 0.25 bpm in order to analyze fetal heart rate variability correctly. Reliable estimation of fetal heart rate and hence of heart rate variability is therefore essential.

Several methods are available to assess the fetal heart rate. These methods differ both at the Doppler signal level (directional or nondirectional) and at the level of the algorithm that estimates the heart rate. For example, existing devices on the market such as Sonicaid Oxford (Oxford Sonicaid Instruments, Abington, UK) [14], Hewlett-Packard 8030A (Hewlett-Packard, Palo Alto, CA, USA) [15], and Philips Avalon F40 (Philips, Amsterdam, Netherlands) use the envelope of the nondirectional Doppler signal. Other authors [16, 17] have used the envelope of the directional Doppler signal. Several algorithms based on autocorrelation are commonly used to estimate the heart rate (Oxford Sonicaid, Hewlett-Packard 8030A, and Philips Avalon F40). These algorithms have been applied either directly to the Doppler signal envelope (directional or nondirectional) or to the signals resulting from discrete wavelet decomposition of the envelope. In the latter case, the final estimated heart rate involves the combination of different estimates [16].

In this study, we first verified whether the pulsed Doppler techniques currently used in commercial devices ensure such an accuracy of 0.25 bpm, and we propose here some recommendations regarding parameter settings. We also compared the techniques used in commercial devices with other pulsed Doppler techniques that use directional Doppler signals. We evaluated the efficacy of all these techniques empirically (error of estimation of the fetal heart rate, sensitivity, and false negative rate).

The originality of this study lies in the recommendations on the parameter settings of the system and in the description of the individual limitations of each technique. Finally, to improve detection probability, a new method based on the combination of heart rates obtained from directional signals is proposed.

2. Materials and Methods

In this section, we describe the Doppler system we developed, the patients, and the synthetic signals used for the comparison of each estimator. The synthetic signals were inferred from real signals. Finally, we present the various existing techniques for estimation of fetal heart rate and we present a new technique based on a combined procedure.

2.1. The Doppler System

In order to evaluate fetal wellbeing objectively and to classify the fetus, we codeveloped the pulsed, multitransducer, multichannel Doppler Actifoetus unit with Althaïs Technologies (Tours, France).

Our system comprised a personal computer (PC) and our Actifoetus unit. The Actifoetus unit contained three groups of four transducers and a Doppler acquisition board. The detailed operating functions of the acquisition board were presented in [18].

The transducers exploring the fetal heart were non-focused and monoelement. They were circular in shape, with a diameter of 13.5 mm and an acoustic power of 1 mW/cm². Geometrically, the transducers were located at the center of gravity and at the top of an equilateral triangle with sides measuring 40.7 mm.

The transducers were placed on the mother’s abdomen. They transmitted a sinusoidal pulse at 2.25 MHz with a pulse repetition frequency (PRF) of 1 kHz. Note that a theoretical accuracy of bpm can be achieved with this value of 1 kHz and accuracy can be still further improved by performing interpolation of the correlation function. The wave was propagated through the mother’s abdomen towards the fetal heart. The backscattered signal was recorded from five different depths, annotated . Note that only one channel was considered in the present study.

The ultrasound signal received was converted into an electrical signal and amplified to compensate for the attenuation of 1 dB/cm/MHz. The signal was then demodulated in phase () and quadrature () [19]. After demodulation, the signals were digitized. The digital outputs of the converters represented the digital Doppler signal.

2.2. Patients

The Doppler signals were acquired at the CHRU “Bretonneau” Tours, France. The consent of each patient was obtained and the study was approved by the Ethics Committee of the Clinical Investigation Centre for Innovative Technology of Tours (CIC-IT 806 CHRU of Tours). Patients were older than eighteen years and all pregnancies were single. The recordings were made during the twenty-fifth and fortieth gestational weeks. Evolution during pregnancy was normal for all fetuses.

2.3. Simulation

Because it was difficult to quantify the effectiveness of the estimation techniques directly on real signals and because there was no suitable model, we generated synthetic signals. These synthetic signals were used as a ground-truth to evaluate the effectiveness of each estimator. To make these signals as realistic as possible, we proceeded in two stages: an analyzing stage deducing the characteristics of the real Doppler signal envelope and a synthetic phase providing realistic simulated signals. Figures 1(a) and 1(b) show 1000 ms of the envelope of a real nondirectional Doppler signal and the two envelopes of corresponding directional signals obtained from the and signals [19, 20]. The synthetic envelope signal was calculated as follows:

where and are the envelopes of directional Doppler signals produced by the scatters that approach and move away from the transducer, respectively.

(a) Envelope of the real nondirectional Doppler signal

(b) Envelopes of real directional signals

Figure 1: A real Doppler signal (kHz) of 1000 ms, recorded with the second transducer in the fourth-channel: (a) Doppler signal (dashed line) and its envelope (solid line); (b) the envelopes of directional signals corresponding to ultrasound scatters approaching to the transducer (solid line) and moving away from the transducer (dash-dot line), respectively.

We verified (Figures 1(a) and 1(b)) that the envelope of the real nondirectional Doppler signal had the signatures of both envelopes of the directional signals. For example, at around 500 ms, the envelope of the nondirectional signal was mainly influenced by scatters that approached the transducer, while at around 400 ms, we observed the influence of movements away from the transducer. The alternating influence of these two movements determined the envelope of the nondirectional signal.

2.3.1. Analysis of Real Directional Signals

Figure 2 shows 2000 ms of the envelope of a real directional Doppler signal. In order to find the important parameters required for the synthesis of this signal, we extracted its intrinsic features (the number and amplitudes of the peaks, the lags between the peaks, and the differences in amplitude between the peaks). The values of these parameters were evaluated by considering a quasiconstant fetal heart rate.

Figure 2: Envelope of a real directional Doppler signal of 2000 ms. The parameters defining the synthetic signal are the amplitudes and durations of peaks, the lag, and the differences in amplitude of two consecutive peaks over time.

Figure 2 represents a sequence of several patterns. These patterns were made up of peaks that corresponded to cardiac wall and valve movements of the fetal heart. As suggested by Shakespeare et al. [21], although six peaks (atrial contraction, ventricular contraction, opening and closing of the mitral valves, and opening and closing of the aortic valves) could be detected theoretically, only a few peaks were in practice detected in the nondirectional Doppler signal. From our analysis, it appeared that the most likely pattern was that with four peaks. Note that this signature composed of four peaks could vary considerably from one beat to another, and it was similar to that identified by Jezewski et al. [13]. Among all these patterns, the most likely was the pattern with peaks in the order 2143; that is, the highest peak was in second position, the second highest peak was in first position, and so forth. For the 2143-pattern, we evaluated the amplitudes of each peak (, , , and ), the peak durations (, , , and ), the lags between two consecutive peaks (, , and ), and the differences in amplitude between two consecutive peaks (, , and ). The results of this statistical analysis are reported in Table 1.

Table 1: Statistics evaluated using () patterns: represents the statistics of the maxima, where ,  , , represent the differences between the positions of two adjacent maxima over time; , , represent statistical differences between two adjacent maxima over time, and represents the statistics of the peak durations. We found that the statistics of the four peaks of , , were identical.

Law

Gaussian
Uniform

(ms)

As the patterns observed in Figure 2 were noisy, we decided to assess the noise level in order to simulate noisy synthetic Doppler signals. We assessed the signal to noise ratio (SNR) as follows:

where and are the powers of the active and passive regions, respectively. We considered the active region as the area containing the pattern peaks, whereas there were none in the passive region. Using (2), we found that the SNR calculated on our real signals corresponded to a Gaussian law: .

2.3.2. Synthesis of Synthetic Directional Signals

Analysis of the envelope of directional signals showed the presence of four peaks, which appeared in order 2143 inside the periodic patterns. The synthesis of such a signal must account for these characteristics. Equation (3) shows the two possible components of such a signal:

where is the noise, is the peak amplitude, is the peak frequency, is the unit rectangular function centered on with width and , is the pattern duration, and is the synthetic signal period. We set a constant interval between the highest peaks of two consecutive patterns of the synthetic signal, as illustrated in Figure 2. We also chose the pattern period as 50% of the synthetic cardiac cycle period , since this period can vary between 40 and 60% [22].

Using (3), we generated two synthetic noisy envelopes corresponding to the envelopes of the directional signals. The envelope of the nondirectional synthetic signal was modeled using (1), being the sum of the envelopes of the both directional synthetic signals. In order to simplify our study, only was calculated, being a delayed and amplified version of . To simulate realistic signals, we introduced a lag between the directional components:

where and are the envelopes of directional signals and is the lag between the two envelopes. is a factor that represents the amplitude ratio between the two types of envelope. From Figure 1,   ms and .

2.4. Estimators

In this section, we describe the different estimators used in our study. Each estimator that was based on the autocorrelation function was denoted by , , as illustrated in Figure 3. Each estimator operated on different signals: , , and .

Figure 3: General diagram of the Doppler data processing acquired using one transducer and one channel. and represent the envelopes of the directional signals, and represents the envelope of the nondirectional signal. , are the two autocorrelation estimators.

Devices existing on the market currently use the envelope and autocorrelation. The estimators that used these configurations were and . The mathematical expression of the two autocorrelation estimators is given in [23] and is represented thereafter by

where is the size of the analyzing window, is the time for which the estimator is computed, and is the lag. represents one of the signals analyzed (, , or ).

We tested other estimators () which used directional signals and , together with and .

2.4.1. Algorithm

The algorithm to estimate the fetal heart rate was the same for all three signals (, , and ). The steps of the algorithm were as follows.

(1)Extract from each signal under consideration (, , or ) a limited number of samples, being the window size.(2)Compute and .(3)Using an empirical threshold, detect the position of   peaks in and .(4)From the position of peaks of and , determine the durations between consecutive peaks with .(5)Calculate the cardiac frequencies with , . This conditional test limits the number of cardiac frequencies estimated from or in the average computation. This conditional test also permits removal of cardiac frequency estimates that are half the expected value, as are sometimes observed (Shakespeare et al. [21]).(6)Calculate the average cardiac frequency (FHR) from not exceeding 35 bpm [13]:

As an illustration, consider a window of  s. Whenever the cardiac frequency was bpm, 16 peaks were observed in the autocorrelation function. Using an empirically set threshold, the duration between each peak was measured ( s s) and cardiac frequencies of bpmbpm were estimated. The average cardiac frequency was obtained by bpm. Note that 4 peaks were observed with 60 bpm and  s were estimated. The average cardiac frequency was obtained by bpm.

Thus the proposed algorithm correctly worked in the range of 60–240 bpm. However, the standard deviation of the FHR estimation was not the same for its extreme values since in one case the average was obtained with three values whereas the average was obtained with fifteen values in the other.

Note that such an algorithm is not perfect since it is hypothesised that the FHR is constant during the process. Sometimes the second peak of the autocorrelation can be lower than the third and the FHR estimate is incorrect. A process must be performed to remove outliers.

2.4.2. Elimination of Outliers

In order to eliminate outlier estimates associated with estimator dysfunction, we introduced a postprocessing step. This postprocessing step was applied only in the case of real signals. An estimate was considered to be an outlier if it laid outside the statistic computed from 40 previous estimates, or if it differed between two consecutive analysis windows by 35 bpm.

2.4.3. Combination

In order to improve the effectiveness of the FHR estimation, we combined estimations. For and , the two values of the fetal heart rate estimated on signals and were combined. The estimate on the two signals was achieved using or . The combination rule we used was as follows:

(i)if the heart rate was detected on a single signal, the combined value took this value;(ii)if the heart rate was detected on both signals, the combined value was the average of the two values.

Note that, in contrast to Kret’s study [16], we combined the fetal heart rates estimated on both directional Doppler signals, while Kret’s technique was based on combination of fetal heart rate estimations computed after discrete wavelet decomposition of the envelope of the directional Doppler signal. Since Kret’s technique was applied only for continuous Doppler signals, it was not taken into consideration in our study.

3. Results

To find the best estimators and the conditions in which they could be used, we performed a series of simulations and experiments. Using simulations, we sought configurations that ensured an error of estimation, that is, the expected accuracy below 0.25 bpm, the highest sensitivity, and the lowest false negative rate. Experimentally, we sought the best configuration for optimal use of the estimator.

3.1. Simulated Signals

We present the results from two types of simulation. In the first series of simulations, we sought parameter settings that ensured the desired accuracy of 0.25 bpm. In the second series of simulations, we evaluated the effectiveness of each estimator in terms of true positive rate and false negative rate.

3.1.1. Optimal Parameter Settings

The results presented in Figures 4 and 5 were obtained for synthetic signals of 30 s. The parameters that varied in our analysis were the periodicity of the signal , the SNR, the window size , and the lag . We varied the signal periodicity between 1000 and 250 ms, as these values corresponded to the standard range of exploration (60–240 bpm) of different fetal monitors. The SNR range varied between 0 and 14 dB, in order to include our measured SNR values on the real signals and in order to take into account the worst cases. The range of size varied between 512 and 4096 ms. The highest fetal heart rate could be obtained with a window size of 512 ms, although we limited the maximum window to 4096 ms to reduce computation time.

(a)

(b)

Figure 4: Duration of the analyzing window required to reach an error of at least 0.25 bpm with the autocorrelation estimator () and a  dB.

(a)

(b)

Figure 5: Duration of the analyzing window required to reach an error of at least 0.25 bpm with the autocorrelation estimator () and with a  dB.

Figures 4 and 5 show the smallest window size analyzed (ms) of all estimators tested that ensured the expected accuracy of 0.25 bpm in the range of 60–240 bpm and that ensured a SNR at least greater than 0.6 dB. Note that for estimator reported in Figure 5, we showed that there was no size which ensured the desired accuracy, whatever the SNR or the frequency. To test estimator robustness in relation to the increasing complexity of the simulated signals, the lag varied between 0 and 40 ms, this value of 40 ms being taken from Figure 1.

The results derived from Figure 4 showed that accuracy for the envelope signal (estimator ) was no longer achieved for certain frequencies, but it still was for directional signals. Finally, Figure 5, shows the best estimators (, , ) and their respective best parameter settings () that ensured an accuracy of 0.25 bpm with a  dB in the 60–240 bpm range.

To summarize, these first results showed the superiority of compared to and the superiority of the envelope of directional signals compared to that of nondirectional signals. We therefore recommend the use of and the estimators (, , and ) based on the envelope of directional signals.

3.1.2. Performance Levels of Estimators

In this study, the performance levels of estimators we wanted to compute were sensitivity and the false negative rate. Fetal heart rates were evaluated every 250 ms from noisy signals with ms. Sensitivity was computed with the equation: , where was the true positive rate and was the false negative rate. Estimates of simulated heart rate were considered to be false negative if they did not ensure the expected accuracy; otherwise, they were true positive. Sensitivity and the false positive rate were evaluated as the average of 30 values. Each value was determined after analysis of a noisy signal of 30 s where sensitivity and the false positive rate had converged to the highest value and to the lowest value, respectively. Convergence was reached for a minimum SNR of 6 dB.

The results of this second series of simulations are presented in Figures 6 and 7 for ms and dB. Note that the estimation of error of 25 bpm was obtained only for whatever the frequency, whereas accuracy for was obtained only for 100, 150, 200, 220, and 240 bpm.

(a)

(b)

Figure 6: True positive rates for and with  dB, , and an error of estimation of 0.25 bpm. (a) True positive rate for . (b) True positive rate for .

(a)

(b)

Figure 7: False negative rates for and with  dB, . Accuracy of 25 bpm was obtained only for whatever the frequency, whereas for accuracy was obtained for only 100, 150, 200, 220, and 240 bpm. (a) False negative rate for . (b) False negative rate for .

The results set out in Figure 7 show that estimators based on (, , , and ) had an average (average obtained from the cardiac frequency) false negative rate of 1.5%, while those based on (, , , and ) presented a higher average false negative rate of approximately 14.8%. The 97.5% average true positive rate of was slightly lower than that of estimators based on , which was 100%. Finally, when the accuracy of 0.25 bpm was reached, we observed that the estimators based on were generally more accurate than those based on , although the average false negative rate was not zero.

Figure 8 shows the error of estimation corresponding to different values of SNR when a zero false negative rate was imposed. This zero false negative rate was obtained by modifying the detection threshold, and a direct consequence was an increase in the estimation. The results derived from Figure 8 showed that the zero false negative rate for was ensured for a dB (below the SNR measured on real signals) and for an error of estimation of 0.8 bpm. In the case of and directional signals, we obtained an error of estimation of 4 bpm, whereas for a nondirectional signal it was 6 bpm.

(a)

(b)

(c)

Figure 8: FHR error of estimation (bpm) when the false negative rate was zero, when the , and when the signals tested were , , and .

Table 2 summarizes the effectiveness of each estimator in terms of FHR error of estimation, SNR, and average false negative rate. To reach an error of estimation, that is, the expected accuracy of 0.25 bpm, we recommend , that is, autocorrelation-based estimators (, , , and ), the price to be paid being an average false negative rate of 1.5%. To reach an average false negative rate of 0%, we recommend autocorrelation-based estimators (, , , and ), where the price to be paid is an error of estimation of 0.8 bpm far from the expected accuracy of 0.25 bpm.

Table 2: FHR error of estimation (bpm) obtained by all estimators tested for different configurations of SNR, false negative rate (FNR), and = 4096 ms.

Estimators
FHR error of estimation (bpm)

, , ,
, ,

FNR = 0%, 6 dB > SNR > 2 dB
0.8
4
6

FNR = 1.5%, SNR > 6 dB
0.25
—
—

3.2. Results Obtained on Real Signals

We recorded 580 minutes for the analysis of real Doppler signals. We selected areas where signals had the cardiac activity signature. The performance levels on these signals were evaluated on the envelopes of both the nondirectional and the directional signals. The FHR estimation obtained with a commercial device (Oxford Sonicaid) was used as a reference to evaluate sensitivity which was evaluated for each estimator. All estimators were evaluated using a size of ms.

The results obtained for all the signals are presented in Table 3. The estimators based on directional signals (, , , and ) provided a higher level of sensitivity compared to those which used nondirectional signals (, ). In the case of directional signals, the results of and were close but slightly better for . This result confirmed the results obtained by simulations. We therefore recommend the use of estimators based on calculated on the directional signals (, ).

Table 3: Sensitivity (%) of estimators for = 4096 ms. and are the two autocorrelation estimators, respectively. , are the envelopes of directional signals, and is the envelope of the nondirectional signal. indicates the combined estimator.

Sensitivity ()
Estimators

88.50% ()
88.43% ()

86.63% ()
84.79% ()

81.79% ()
75.05% ()

()
95.48% ()
94.93% ()

Sensitivity was improved using the combination method. The sensitivity of estimator was 95.5% (see Table 3). Using the combination method, the sensitivity increased to about (95.5%–88.5%) 7% compared to directional signals and to about (95.5%–81.8%) 14% compared to nondirectional signals .

4. Discussion and Conclusion

In this study, we focused on different settings of the estimators (window size, lag) that ensured a fetal heart rate estimation with a maximum authorized error of 0.25 bpm. We found in simulation that only estimators based on and directional signals could ensure such an accuracy of 0.25 bpm. The size necessary in this case was ms.

Note that, although we proposed synthetic signals that were as realistic as possible, we are aware that the plotted performance levels are representative only of our simulations and not of all cases encountered in practice. It is likely that the performance levels of the algorithms tested can be reduced in the presence of artefacts. However, the 95% sensitivity obtained from real signals suggests that our proposed estimators may be trusted.

In the case of real signals, the sensitivity was quantified. Since in our study the estimated SNR on the real signals was greater than the threshold of 6 dB (deduced on simulated signals) required to reach the desired accuracy of 0.25 bpm, a denoising filter was not necessary. However, in cases of a SNR lower than 6 dB, a denoising process (Wiener, wavelet) could be introduced.

Sensitivity was quantified using a size of 4096 ms. We found that the estimators , , and based on had slightly greater sensitivity than those based on . We therefore recommend the use of .

Various cases were considered on the basis of this study, that is, those that do not require a precise estimate of the fetal heart rate and those for which accuracy is critical. The accuracy of fetal heart rate estimation in the first case is not important but the false negative rate should be as low as possible. For example, if the goal of a monitoring system is simply to verify that the fetal heart rate is in the normal range (110–160 bpm), very high accuracy is not needed. In this case, an error of estimation of 0.8 bpm is sufficient. Our computations showed that for an error of estimation of 0.8 bpm, a zero false negative rate in the zones when the rhythm is quasi-constant could be ensured. In the second case, an error of estimation of 0.25 bpm is required for a system in which the goal is not only to estimate the heart rate, but also to evaluate fetal wellbeing. Our study showed that for this type of system, the false negative rate may be slightly higher than zero. It is important to note that this error of estimation was guaranteed for a quasi-constant heart rate. This is not a constraint for such a system, since the variability of fetal heart rate must be evaluated in these ranges to predict fetal distress.

Applied to real signals, the estimators based on provided sensitivity close to those of , and the most efficient of these estimators were those that used directional signals (, , and ).

A 7% increase in sensitivity compared to estimators based on individual directional signals was possible when we combined the two heart rates calculated on the directional signals. A 14% increase in sensitivity compared to estimators based on individual nondirectional signals was possible when we combined the two heart rates calculated on the directional signals. When a combination was used, both signals were processed in parallel, thus doubling the number of operations.

The good levels of performance of our estimator based on this combination suggest first that it can be adapted to multitransducer, multichannel configurations and second that such an estimator will improve fetal diagnosis.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This study was supported by the Agence Nationale de la Research (Project ANR-07-TECSAN-023, Surfoetus). The authors would like to thank the Clinical Investigation Centre for Innovative Technology of Tours (CIC-IT 806 CHRU of Tours) and Professor F. Perrotin’s Obstetric Department Team for their support in recording the signals. The authors would like to thank the anonymous referees for their valuable comments.

References

F. Manning, I. Morrison, I. Lange, C. Harman, and P. Chamberlain, “Fetal assessment based on fetal biophysical profile scoring: experience in 12, 620 referred high-risk pregnancies. I. Perinatal mortality by frequency and etiology,” The American Journal of Obstetrics and Gynecology, vol. 151, no. 3, pp. 343–350, 1985.

H.-Y. Chen, S. P. Chauhan, C. V. Ananth, A. M. Vintzileos, and A. Z. Abuhamad, “Electronic fetal heart rate monitoring and infant mortality: a population-based study in the United States,” The American Journal of Obstetrics and Gynecology, vol. 204, no. 6, pp. S43–S44, 2011.

G. G. Berntson, J. Thomas Bigger Jr., D. L. Eckberg et al., “Heart rate variability: origins methods, and interpretive caveats,” Psychophysiology, vol. 34, no. 6, pp. 623–648, 1997.

M. Malik and A. J. Camm, “Heart rate variability,” Clinical Cardiology, vol. 13, no. 8, pp. 570–576, 1990.

M. David, M. Hirsch, J. Karin, E. Toledo, and S. Akselrod, “An estimate of fetal autonomic state by time-frequency analysis of fetal heart rate variability,” Journal of Applied Physiology, vol. 102, no. 3, pp. 1057–1064, 2007.

B. Arabin and S. Riedewald, “An attempt to quantify characteristics of behavioral states,” The American Journal of Perinatology, vol. 9, no. 2, pp. 115–119, 1992.

B. Frank, B. Pompe, U. Schneider, and D. Hoyer, “Permutation entropy improves fetal behavioural state classification based on heart rate analysis from biomagnetic recordings in near term fetuses,” Medical and Biological Engineering and Computing, vol. 44, no. 3, pp. 179–187, 2006.

S. Lange, P. Van Leeuwen, U. Schneider et al., “Heart rate features in fetal behavioural states,” Early Human Development, vol. 85, no. 2, pp. 131–135, 2009.

J. A. DiPietro, D. M. Hodgson, K. A. Costigan, S. C. Hilton, and T. R. B. Johnson, “Fetal neurobehavioral development,” Child Development, vol. 67, no. 5, pp. 2553–2567, 1996.

M. Ferrario, M. G. Signorini, G. Magenes, and S. Cerutti, “Comparison of entropy-based regularity estimators: application to the fetal heart rate signal for the identification of fetal distress,” IEEE Transactions on Biomedical Engineering, vol. 53, no. 1, pp. 119–125, 2006.

A. Kikuchi, T. Shimizu, A. Hayashi et al., “Nonlinear analyses of heart rate variability in normal and growth-restricted fetuses,” Early Human Development, vol. 82, no. 4, pp. 217–226, 2006.

G. S. Dawes, M. O. Lobbs, G. Mandruzzato, M. Moulden, C. W. G. Redman, and T. Wheeler, “Large fetal heart rate decelerations at term associated with changes in fetal heart rate variation,” The American Journal of Obstetrics and Gynecology, vol. 188, no. 1, pp. 105–111, 1993.

J. Jezewski, J. Wrobel, and K. Horoba, “Comparison of Doppler ultrasound and direct electrocardiography acquisition techniques for quantification of fetal heart rate variability,” IEEE Transactions on Biomedical Engineering, vol. 53, no. 5, pp. 855–864, 2006.

J. Pardey, M. Moulden, and C. W. G. Redman, “A computer system for the numerical analysis of nonstress tests,” The American Journal of Obstetrics and Gynecology, vol. 186, no. 5, pp. 1095–1103, 2002.

E. Courtin, W. Ruchay, P. Salfeld, H. Sommer, and A. Versatile, “Semiautomatic fetal monitor for non-technical users,” Hewlett-Packard Journal, vol. 28, no. 5, pp. 16–23, 1977.

T. Kret and K. Kaluzynski, “The fetal heart rate estimation based on continuous ultrasonic Doppler data,” Biocybernetics and Biomedical Engineering, vol. 26, no. 3, pp. 49–56, 2006.

I. Voicu, D. Kouame, M. Fournier-Massignan, and J.-M. Girault, “Estimating fetal heart rate from multiple ultrasound signals,” in Proceedings of the International Conference on Advancements of Medicine and Health Care through Technology, pp. 185–190, Berlin, Germany, September 2009.

A. Kribèche, F. Tranquart, D. Kouame, and L. Pourcelot, “The Actifetus system: a multidoppler sensor system for monitoring fetal movements,” Ultrasound in Medicine and Biology, vol. 33, no. 3, pp. 430–438, 2007.

J. A. Jensen, Estimation of Blood Velocities Using Ultrasound, Cambridge University Press, Cambridge, UK, 1996.

N. Aydin, L. Fan, and D. H. Evans, “Quadrature-to-directional format conversion of Doppler signals using digital methods,” Physiological Measurement, vol. 15, no. 2, article 007, pp. 181–199, 1994.

S. A. Shakespeare, J. A. Crowe, B. R. Hayes-Gill, K. Bhogal, and D. K. James, “The information content of Doppler ultrasound signals from the fetal heart,” Medical and Biological Engineering and Computing, vol. 39, no. 6, pp. 619–626, 2001.

E. Hernandez-Andrade, H. Figueroa-Diesel, C. Kottman et al., “Gestational-age-adjusted reference values for the modified myocardial performance index for evaluation of fetal left cardiac function,” Ultrasound in Obstetrics and Gynecology, vol. 29, no. 3, pp. 321–325, 2007.

A. V. Oppenheim and R. W. Schafer, Discrete-Time Signal Processing, Prentice Hall Signal Processing, Prentice-Hall International, Englewood Cliffs, NJ, USA, 1989.