Influence of electrode placement on signal quality for ambulatory pregnancy monitoring.
Monitoring of the fetal heart rate (fHR), fetal electrocardiogram (fECG) morphology, and uterine activity is important for fetal health assessment during pregnancy and delivery. Current monitoring methods are, however, not suited for long-term observation throughout pregnancy. The most commonly used method for fetal monitoring during pregnancy uses Doppler ultrasound and a strain gauge to determine the fHR and uterine activity, respectively. While allowing for noninvasive measurements, Doppler ultrasound measurements need continuous attention of a trained physician, making them unsuitable for long-term observation . Furthermore, in many situations these measurements do not provide conclusive information for accurate assessment of fetal health and, therefore, additional information is needed for clinical decision making . During delivery, a fetal scalp electrode and an intrauterine pressure catheter (IUPC) can be used as part of an internal cardiotocogram (CTG) .
Although this method provides accurate fECG and IUP readings, thus improving decision making, it requires rupturing of the membranes; it can, therefore, be applied during delivery only.
Recently published methods enable noninvasive extraction of the fECG from abdominal recordings, the results of which have shown to be more reliable and accurate than ultrasound [4-11]. Additionally, based on the electrohysterogram (EHG), which is a recording of the electrical activity causing uterine muscle contractions, an intrauterine pressure (IUP) estimate can also be obtained [12, 13]. As both methods use abdominal contact electrodes, a single set of electrodes can enable prolonged noninvasive monitoring of both fECG and IUP. This allows for determination of the fetal health and onset of labor throughout pregnancy with a relatively high patient comfort.
Numerous electrode configurations have been presented for both fECG and EHG measurements, using a wide variety of leads ranging in numbers from only a few to one hundred [6, 7, 12, 14-17]. Although abdominal electrode placement determines the measured signal quality to a great extent, no consensus has been reached on standard electrode placement. Because of the prolonged measurement periods in ambulatory monitoring, comfort in placing and wearing the sensors is important, imposing additional requirements on the electrode spacing and number of sensors used. Also, the number of sensors greatly influences the total power consumption of the monitoring system, which is a major concern in ambulatory monitoring [15,18].
The goal of this paper is to determine an electrode grid that allows for a dynamic selection of the optimal electrode positions for an assessment of the fetal health and onset of labor in an ambulatory setting. More specifically, we focus on determination of the fetal heart-rate and uterine contraction frequency near the end of the pregnancy, when monitoring these parameters is more relevant. To this end, the influence of abdominal electrode placement on the signal quality of the fetal and maternal electrocardiogram (ECG) as well as the EHG is explored. Various bipolar electrode placements are evaluated based on both simulated and measured signals to determine the influence of scale and direction on measurement quality for all signals of interest. For simulation, a dipole-based model is used to calculate the potential on the thoracic surface, with a maternal and fetal vector cardiogram (VCG) as signal sources. These VCGs can be extracted from any standard thoracic or abdominal ECG recording, respectively, to simulate a multitude of electrode placements. For the measurements, an electrode placement protocol is employed that enables measurement of the maternal and fetal ECG as well as the EHG from a plurality of bipolar leads in multiple directions and interelectrode scales. A bipolar setup is used because it reduces the common mode interference, while enhancing the relative contribution of the signals of interest, depending on the lead orientation [6, 19]. Both simulations and measurements are used to determine the influence of electrode placement on the recorded signal quality, defined by its signal-to-noise ratio (SNR). Results show that both electrode orientation and distance significantlyinfluence the measurement quality of the signals of interest.
First, an introduction of all physiological parameters influencing abdominal measurements is presented in Section 2. Next, Section 3 outlines a strategy to find an optimal electrode grid for ambulatory fetal health monitoring, where a model for simulation of the maternal and fetal ECG is proposed in Section 3.1, followed by the measurement setup description in Section 3.2. Finally, Section 3.3 introduces all criteria for choosing electrode positions for an optimal electrode grid. The simulation and measurement results are shown in Section 4. Section 5 contains a discussion on the obtained results and provides guidelines for choosing an optimal electrode grid.
2. Physiological Parameters
Various physiological signals play an important role in pregnancy monitoring. The fECG can be used to determine the fHR with beat-to-beat accuracy, but it also allows for analysis of the fECG morphology . The EHG on the other hand can be used to estimate the IUP and recognize the first signs of approaching delivery . During labor, the fHR can then be used in combination with the EHG. The pressure exerted on the fetus by the contracting uterus acts as a stress test, and the fetal heart rate variability gives an indication of the level of fetal distress . Multiple noise sources are also present in abdominally measured signals, the most noticable of which are the maternal electrocardiogram (mECG) and electromyographic signals from the abdominal muscles . A list of the main signals that are present in abdominal recordings is presented hereafter.
2.1. Fetal ECG. The fECG is a representation of the potential map in the fetal heart resulting from de- and repolarization of myocytes over time. For a fetus with a gestational age between 32 to 41 weeks, the fetal HR is in the range of 60240 BPM with an expected HR around 140 BPM, while a higher instantaneous HR is possible due to arrhythmias . Most of the QRS spectral energy is contained in the 20-60 Hz range, with a peak-to-peak amplitude in the range of 3-25 [micro]V.
2.2. Electrohysterogram. The EHG is a bioelectric signal associated with the propagation of action potentials through uterine muscle cells. It is the primary cause of uterine muscle contractions. Therefore, the EHG can be used as an estimator of uterine activity and permits estimation of the IUP associated with each contraction . Contrary to the myocardium, the action potential waves propagating through theuterusdonot have afixedpacemaker location, resulting in an erratic direction of propagation. Research has shown that the main EHG frequency components are in the 0.3-0.8 Hz band with a peak-to-peak amplitude >100 [micro]V 14, 20, 24].
2.3. Maternal ECG. The mECG is generally the most prominent signal present in abdominal ECG recordings. Therefore, the mECG needs to be removed to enable reliable fECG detection. This procedure typically requires accurate detection of the maternal R-peaks and delineation of the the P-wave, QRS complex, and T-wave . For pregnant women, the HR can range from 50 to 210 BPM with an expected HR around 80 BPM, while most of the QRS spectral energy is contained in the 10-30 Hz frequency range [25,26]. Signal amplitude and morphology both depend on electrode placement. The amplitude is in the order of 100 [micro]V when measured on the abdomen, which is 4-20 times that of the fECG [6, 27].
2.4. Additional Noise Sources. Abdominal ECG and EHG recordings are affected by various noise sources. Apart from the mECG, abdominal recordings can contain various high-amplitude noise sources, which are typically highly nonstationary in time and colored in spectrum 28, 29].
Skeletal muscle artifacts introduce high frequency components between 10 Hz and 500 Hz during skeletal muscle activity, when recorded abdominally. Inparticular during a contraction, the skeletal muscle electromyogram can be regarded as the major noise contribution [24, 30]. Electrode motion artifacts are transient baseline changes which originate from changes in the electrode-skin impedance and changes of skin potential due to deformation of the skin . On average, electrode motion artifacts are 100-500 ms in length with an amplitude of up to five times that of the mECG .
Fetal movements can cause skeletal muscle artifacts, and more importantly, can result in a different orientation of the fetal heart vector with respect to the electrode grid, changing the amplitude and morphology of the measured signals .
The vernix caseosa is a thin fatty layer which starts developing on the fetal skin from the 28th week of gestation, changing the conductivity in the maternal volume conductor. This severely reduces the measured fECG amplitude until the vernix caseosa partly dissolves, which typically starts in the 32nd week of gestation .
In an effort to determine the optimal electrode placement for ambulatory fetal health monitoring, the influence of electrode placement on the signal quality of the mECG, fECG, and EHG is determined. To this end, an electrophysiological model for simulation of the abdominal skin potential and a set of abdominal measurements are introduced.
The model described in Section 3.1 can simulate the electrical potential on the body surface as a result of electrical field vectors of the maternal and fetal hearts. This model can be used to simulate any electrode configuration on the maternal thorax based on standard thoracic and abdominal measurements and used to determine and verify the optimal measurement positions for both mECG and fECG. A similar simulation of the EHG is not possible due to the unknown origin and propagation direction of the signal .
Section 3.2 describes a set of abdominal measurements and a method of extracting bipolar leads at various interelectrode scales and in various orientations. These bipolar leads allow for comparison of the relative signal quality of all ECG and EHG signals at different positions and orientations.
Section 3.3 defines the criteria for selection of optimal electrode sets for maternal and fetal ECG as well as EHG measurement based on their SNR. The leads with the highest SNR for each patient are selected based on the results for both simulated and measured signals. Here, the SNR calculation of the EHG signal is based solely on abdominal measurements. The statistical significance of the SNR improvement over the other bipolar leads is determined based on the P value to determine which changes in direction and interelectrode distance result in significant SNR reduction.
3.1. Electrophysiological Model
3.1.1. The Maternal Heart. Cardiac contractions are the result of the generation and conduction of electrical pulses, referred to as action potentials, through the myocardium. Propagation of these action potentials acts as the movement of numerous electrical dipoles, which can be modeled as a single dipole.
This results in an electrical cardiac field vector [??] [35, 36], which can be represented as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)
where [q.sub.i] and [d.sub.i] are the size of the ith charge and its distance from the origin, respectively, and [[??].sub.x], [[??].sub.y], and [[??].sub.z] are three orthonormal vectors related to the three body axes. This cardiac vector originates in the center of the heart and varies over time in both amplitude and orientation defining a specific trajectory, which is referred to as the VCG. Using the Moore-Penrose inverse of the Dower matrix [37, 38], also called the inverse Dower transform, it is possible to extract p from a standard 12-lead ECG [39, 40]. The inverse Dower transform takes into account the standard locations of the recording leads and the attenuation of the body volume conductor for each electrode; it is defined as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)
where [V.sub.1] - [V.sub.6] are the six precordial leads and [D.sub.I] and [D.sub.II] are Einthoven leads I and II, respectively. The resulting cardiac vector [??] can now be obtained from a standard ECG recording and used to simulate the potential map on the skin surface. Assuming the body volume conductor to be a homogeneous conductor, a linear projection of [??] can be used to estimate the ECG signal. In case of a homogeneous infinite-volume conductor, the potential v on the body surface generated by [??] is defined by the simplified model
v = [1/4[pi][sigm]][[??] x [??]/[[??].sup.3] (3)
where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is the vector connecting the origin of the cardiac vector p with the electrode location and <r is the conductivity constant of the conductive medium. Figure 1 shows the relative locations of both the maternal and fetal cardiac vector, [[??].sub.m] and [[??].sub.j], the approximate electrode placement on the maternal abdomen, and [??] connecting both.
The model can now be used to calculate the potential of a simulated ECG anywhere on the thorax using the VCG obtained from a standard 12-lead ECG. Here, recordings I12 and I32 from the St. Petersburg INCART 12-lead Arrhythmia Database were selected for maternal VCG extraction based on gender and age of the subjects. Nonarrhythmic episodes with a total length of 30 min were selected for the female patients, 39 and 41 years of age. Figure 2 shows the coronal, sagittal, and transverse projections of a VCG extracted from recording I12, as well as a plot of Einthoven lead II, both measured and simulated using (3).
3.1.2. The Fetal Heart. The fetal heart model is similar to that of the maternal heart with the exception of its position and orientation relative to the electrodes. However, the used fetal VCG (fVCG) differs from the VCG used for simulation of the maternal cardiac activity. A matrix operation similar to the inverse Dower matrix in (2) can be used to extract the fVCG from the abdominal electrodes. To this end, the matrix coefficients are adjusted based on the spatial locations of the electrodes with respect to the fetal heart, as discussed in . For the measurement setup as described in , the adjusted matrix to obtain the fVCG [[??].sup.f] is defined as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)
where [L.sub.1] - [L.sub.8] are leads 1-8 as described in . Two abdominal measurements with a high SNR measured at a gestational age of 22 weeks and 24 weeks, respectively, were selected from the dataset described in  for extraction of a fVCG. Figure 3 shows the three orthogonal projections of a fVCG obtained from these recordings, as well as a plot of the 2nd electrode ([L.sub.2]) of the abdominal recording, both measured and simulated. It can be noted that the fVCG is rotated with respect to the visualization planes compared to the maternal VCG, due to the orientation of the fetal body. Additionally, the fVCG is much more erratic due to the lower SNR.
3.1.3. Noise Sources. Reduction of the recorded ECG signals to three dimensions using the inverse Dower transform results in reduction of the noise, as can be observed in Figure 3(d), where the simulated signal is smoother than the original measurement. In these simulations the remaining noise from the original ECG, in fact, originates from the same spatial location as the ECG signal. Therefore, the noise amplitude in the simulation scales at the same rate as the ECG with changing electrode distance. Hence, the SNR will remain equal irrespective of the interelectrode distance. Electrode motion artifacts and skeletal muscle artifacts in real measurements are localized electrical signals and, therefore, the amplitude of these artifacts is not influenced by interelectrode distance. As a consequence, the ECG signal amplitude typically scales up faster with increasing interelectrode distance than the noise, resulting in an increased SNR. As the frequency spectrum of skeletal muscle artifacts shows the largest overlap with that of the fetal and maternal ECG, modeled skeletal muscle noise is added to each of the signals at the simulated electrode positions.
The skeletal muscle artifacts are modeled by coloring white noise based on the typical electromyogram power spectrum. The filter used for noise coloring is a 2nd order band-pass filter combined with a low-pass filter tuned to match the recorded power spectrum in . As the amplitude of artifacts due to fetal movements scales with changing interelectrode distance, this type of noise is not added to the model.
The modeled skin potential [v.sub.0] can now be described for each electrode by
[v.sub.0] (t) = [v.sub.m] (t) + [v.sub.f] (t) + [alpha] x (h x w) (t), (5)
where [v.sub.m] and [v.sub.j] are the maternal and fetal simulation results, respectively, h(t) is the impulse response of the skeletalmuscle-artifacts-coloring filter, w(t) is white noise, a is a scaling factor for the skeletal muscle artifacts, and * is the convolution operator. The scaling factor a is chosen such that the SNR of the best bipolar lead is comparable to that of recordings as described in Section 3.2.
3.2. Abdominal Measurements. A set of test measurements was performed to determine optimal electrode placement for ambulatory fetal monitoring. The dataset contains five abdominal measurements of 20 min each on pregnant women at full gestation (39 weeks and 4 days [+ or -] 12 days) with a total of 7,148 maternal and 13,612 fetal heart beats, as well as 35 uterine contractions. Ultrasonography by medical personnel showed all fetuses have a cephalic presentation and no fetuses making any mayor thoracic movements. Simultaneously with each measurement, an external cardiotocogram was recorded to monitor both fHR and uterine activity. All measurements were performed with an ethical approval at the Maxima Medical Center in Veldhoven (The Netherlands), after patients signed an informed consent.
All measurements were performed using a Refa system (TMS International, Enschede, the Netherlands), comprising a multichannel amplifier for electrophysiological signals, with a noise level <1 [micro]V RMS and a discretization level of 18.4nV/bit. An electrode grid with a triangular layout, as shown in Figure 4(a), was used for all measurements. This layout was chosen because a minimum of three electrodes are needed to register a two-dimensional signal, and because it is easily scalable to a different grid size while reusing a large part of the electrodes. In addition to the electrodes shown in Figure 4, ground and reference electrodes were placed on the hip [6, 12]. The triangular grid shown in Figure 4(a) can be used to analyze recorded signals at three different interelectrode scales, each consisting of six electrodes as shown in Figure 4(b). At each of the three indicated scales, bipolar derivations were calculated in 6 different directions, resulting in a total of 18 bipolar signals for each measurement. Because a triangular grid is used to determine bipolar signals at 30[degrees] increments, the interelectrode distance within a single scale varies. Bipolar leads II, IV, and VI are a fraction [[delta].sub.[alpha]] = 1 - [square root 3/4] shorter than those of leads I, III, and V.
3.3. Quality Measures. To define the ability of the electrode grid to capture all electrophysiological characteristics of interest, a signal quality criterion for both the fetal and maternal ECG as well as the EHG is defined.
3.3.1. Maternal and Fetal ECG Quality Measure. An estimate of the maternal and fetal SNR is used as the quality measure for the maternal and fetal ECG, respectively. For the mECG, the abdominal ECG is first preprocessed to only include noise in the same frequency band as the ECG, which can influence ECG detection and extraction results. Band-pass filtering with cut-off frequencies of 2 Hz and 70 Hz is used to remove baseline wander and high-frequency noise. Powerline interference is removed by means of a notch filter. Next, the SNR is calculated for a single maternal QRS complex based on the method described in . This SNR is defined as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)
where [V.sub.mQRS,PP] is the peak-to-peak voltage of the maternal QRS-complex in a 100 ms wide window around the annotated peak location, [S.sub.noise] is the noise signal in the 2 segments between the surrounding QRS-complexes, and N is the number of samples in [S.sub.noise]. Here, the signal power can be derived from the square of [V.sub.mQRS,PP] by dividing it by 8, when approximating the [Q.sub.RS] complex as a sinusoidal signal . Because [V.sub.mQRS] PP is determined from the measured signal, which includes noise, the calculated SNR defined above is actually a ratio of the signal-plus-noise to the noise and as a consequence is typically over-estimated.
SNR calculations for the fECG are similar to those for the mECG, with the exception of the removal of the mECG prior to fetal SNR (SNRy) calculation and changes in the band-pass cut-off frequencies. Various methods for removal of the mECG have been presented in the literature [4-9]. As the method proposed by Vullings et al. outperforms the other techniques in both mECG removal and fHR detection, the maternal ECG was removed using dynamic segmentation and linear prediction . Based on the assumption that all maternal ECG components are removed prior to fetal R-peak detection using , the amplitude of the maternal R-peaks becomes irrelevant when calculating the fetal SNR.
Additionally, the cut-off frequencies of the filter are set at 5 Hz and 100 Hz, respectively, to account for the different spectral content of the fECG. For the calculation of both maternal SNR and fetal SNR, the R-peak locations must be known. To this end, all R-peaks in the dataset have been manually annotated.
3.3.2. EHG Quality Measure. Similar to the ECG, the EHG quality measure is defined by its SNR, where signal and noise segments were defined as contraction and noncontraction segments in the abdominal EHG signal, respectively. The contraction timing was defined as periods with an IUP increase of over 10% compared to the baseline pressure, as measured with the IUPC reference signal. Signal segments with a 10% rise over the baseline shorter than 20 s were rejected. A SNR was obtained for each contraction, as defined by
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)
where [S.sub.EHG] and [S.sub.noise] are the signal and noise power in the 0.3-0.8 Hz frequency band, respectively, in line with , and M is the number of samples in [S.sub.EHG].
Figure 5 shows a 4 s segment of both simulation results of the model described in Section 3.1 and a recording as described in Section 3.2, respectively. Both use the electrode grid as introduced in Figure 4(b) with an interelectrode distance of 16 cm. In both plots the mECG is an order of magnitude larger than the fECG and that the orientation of the bipolar measurement influences the amplitude and morphology of both maternal and ECG signals. In both cases, leads I and III are the leads with the largest QRS energy for the maternal and fetal ECG, respectively.
Tables 1 and 2 show the SNR in simulated signals for the maternal and fetal ECG, respectively. Next to the three scales at which measurements were performed, two additional scales with 12 cm and 20 cm interelectrode distance were simulated. The bold faced SNR in each row indicates the bipolar channel with the highest SNR at the current scale, for each of the simulations. The SNR for both mECG and fECG shows an increase with increasing interelectrode distance.
Tables 3 and 4 show the SNR, for the maternal and fetal ECG, respectively, for the abdominal measurements as described in Section 3.2. The bold faced SNR in each row indicates the bipolar channel with the highest SNR at the current scale, for each of the patients. Similar to the simulations, the SNR for both maternal and fetal ECG shows a clear increase with increasing interelectrode distance. For the maternal ECG, a clear preferential measurement direction can also be observed.
Figure 6 shows part of the fECG SNR in six bipolar derivations at the 16 cm scale. For display purposes, the SNR has been smoothened taking the running average of the SNR in a 20 s window. Even though the SNR in all leads fluctuates over time, Lead IV clearly has the highest SNR throughout the recording. Leads I and VI, which are at a large angle with lead IV, have the lowest SNR.
Figure 7 shows the mean and standard deviation of the EHG SNR for each of the six possible bipolar lead orientations and each of the three interelectrode scales. From this figure it can be observed that, on average, a horizontal measurement direction gives the best EHG measurement results, and the SNR with respect to the interelectrode distance contains an optimum between 8 cm and 16 cm.
The statistical significance of the obtained signal quality at various measured electrode locations is expressed by its P value, obtained from a number of t-tests. A value below 0.05 is considered to be statistically significant. Figure 8 shows the statistical significance of a relative change of interelectrode distance and rotation, with respect to the optimal electrode placement. Additionally, the 0.05-threshold line is shown in both figures indicating which changes in interelectrode distance or orientation are significant.
5. Discussion and Conclusion
Long-term ambulatory fetal monitoring based on noninvasive abdominal measurement of the mECG, fECG, and EHG is an important step closer to becoming feasible when using integrated on-board signal processing to reduce power consumption of wireless transmission to a negligible level. Efficiency of the signal processing tasks is proportional to the number of leads and signal quality, making the choice of optimal electrode number and placement increasingly important. In this paper, the influence of abdominal electrode placement on the measurement quality of mECG, fECG, and EHG signals is explored. The aim is to obtain an electrode grid for ambulatory measurements that guarantees, at all time, optimal electrode pairs for extraction of the fetal heartrate and uterine contraction frequency in the final stage of pregnancy. Both simulated and measured signals are used to evaluate various electrode placements in an effort to determine the influence of scale and direction on measurement quality for all signals of interest. As the orientation of the fetus is a priori unknown and, as a rule, the optimal measurement direction does not equal that of the mECG, an electrode grid is needed for optimal detection of both signals. Based on these results, an electrode grid is proposed, which allows for optimal registration of all signals of interest.
For evaluation of the maternal and fetal ECG, a dipole-based model of the heart using a VCG as input was introduced to simulate the potential due to cardiac activity on the skin surface. The model allows for realistic simulation of the potential at any position on the maternal abdomen and simulation results can be used to determine the influence of electrode placement on the SNR of maternal and fetal ECG signals. Simulation results based on 2 standard 12-lead ECG recordings and 2 in-house abdominal recordings show a strong influence of both the interelectrode distance and orientation of bipolar electrode pairs on the SNR. Depending on the morphology of the ECG, a 90[degrees] rotation of the bipolar lead orientation or a reduction of the interelectrode distance by half can lower the SNR by as much as 5 dB.
A set of test measurements was performed to validate and verify the obtained simulation results for the mECG and fECG and obtain quantitative information on the EHG. To this end, a triangular measurement grid was used, which enables measurement of all signals of interest from a plurality of bipolar leads in multiple directions and interelectrode scales. For each of the bipolar derivations, as described in Section 3.2, the SNR was determined for the mECG, fECG, and EHG signals.
For the ECG signals, the influence of interelectrode distance and lead orientation observed in the measurement results clearly show the same trend as the simulation results. The optimal measurement direction on the abdomen for the mECG is horizontal, where an increase in the SNR can be observed with increasing interelectrode distance up to 16 cm. The SNR of the fECG also increases with increasing interelectrode distance, but, contrary to the mECG, it does not show a single electrode pair clearly having the highest SNR for all recordings. This can be related to the changing fetal orientation, resulting in morphological changes in the recorded signal due to a varying angle between the fetal cardiac vector and the bipolar lead. For simulation 1 and patients 2 to 5 in the measurement set, the approximate direction of lead IV can be observed to contain the highest SNR. Simulation 2 and patient 1 show a different optimal direction, as those fetuses displayed a sacrum anterior and occipitoposterior presentation, respectively.
To determine the significance of a change in interelectrode distance or orientation, a f-test was performed on every combination of interelectrode distance and orientation with respect to the electrode pair with the highest SNR. A change in interelectrode distance by only 20% produces a statistically significant SNR reduction for the mECG, while for the fECG no significant result was obtained. This fECG result can most probably be ascribed to the low overall SNR, which results in a high relative variation in SNR. A change in measurement orientation by an angle >55[degrees] compared to the optimum results in a significant drop in SNR for both ECG signals. Using more than 3 electrodes in a circular layout is therefore not expected to significantly increase the SNR of either the maternal or fetal ECG.
Contrary to the ECG, the SNR of the EHG signal does not consistently improve with increasing interelectrode distance but shows a clear optimum between 8 cm and 14 cm. A f-test shows that interelectrode distances of both 4 cm and 16 cm yield a significantly lower SNR. The optimum is introduced because bipolar measurements of propagating signals near the electrodes introduce dips in the signal power spectrum, the frequencies of which are dependent on interelectrode distance and propagation velocity . With an interelectrode distance of 10 cm and an assumed maximal propagation velocity of 10cm/s , the bandwidth of the main lobe of this spatial filter is 1 Hz. Increasing the interelectrode distance beyond 10 cm results in the bandwidth reducing below the EHG frequency band. On average, a horizontal measurement direction of the EHG shows the highest signal quality, although no statistically significant results are obtained. Given the complex underlying mechanisms, different factors, for example, respiration, may contribute to these results . Additionally, the a priori unknown origin and direction of propagation of the uterine electrical activity reduce the statistical significance . As the optimal measurement direction during each contraction is different, measuring multiple bipolar lead directions simultaneously is expected to improve reliable contraction detection compared to a single bipolar pair.
Based on our results for the mECG, fECG, and EHG, an electrode grid for optimal measurement of all signals of interest can be defined. In the final stage of pregnancy, when the fetus typically has a cephalic presentation, an electrode grid as shown in Figure 9(a) offers optimal measurement results, considering the reduced freedom of movement of the fetus. Earlier throughout the pregnancy the fetus has much more room to move around, and additional fECG measurement directions may be required to guarantee optimal signal quality. Tracking the current fetal position using an algorithm as presented in  will allow for optimal electrode selection. An electrode grid with six electrodes at a scale of 20 cm, as shown in Figure 9(b), offers 360[degrees] fECG measurements and various EHG signals. While giving optimal results for the EHG measurements, simulation results show that also the maternal and fetal ECG are obtained with a best possible quality given the posed size constraints. Dynamic electrode selection might be used to only record leads with the highest SNR for mECG, fECG, and EHG, reducing power consumption in the analog front-end and digital signal processing.
Future work will focus on methods for efficient and dynamic electrode selection to reduce power consumption in a device aimed at ambulatory pregnancy monitoring. Additional insights could be provided by an extended dataset including measurements at lower gestational ages.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was supported by the Dutch technology foundation STW
 H. Kieler, S. Cnattingius, B. Haglund, J. Palmgren, and O. Axelsson, "Sinistrality-a side-effect of prenatal sonography: a comparative study of young men," Epidemiology, vol. 12, no. 6, pp. 618-623, 2001.
 I. Amer-Wahlin, C. Hellsten, H. Noren et al., "Cardiotocography only versus cardiotocography plus ST analysis of fetal electrocardiogram for intrapartum fetal monitoring: a Swedish randomised controlled trial," The Lancet, vol. 358, no. 9281, pp. 534-538, 2001.
 E. H. Hon, R. H. Paul, and R. W. Hon, "Electronic evaluation of fetal heart rate. XI. Description of a spiral electrode," Obstetrics and Gynecology, vol. 40, no. 3, pp. 362-365, 1972.
 L. De Lathauwer, B. De Moor, and J. Vandewalle, "Fetal electrocardiogram extraction by blind source subspace separation," IEEE Transactions on Biomedical Engineering, vol. 47, no. 5, pp. 567-572, 2000.
 V. Zarzoso and A. K. Nandi, "Noninvasive fetal electrocardiogram extraction: blind separation versus adaptive noise cancellation," IEEE Transactions on Biomedical Engineering, vol. 48, no. 1, pp. 12-18, 2001.
 S. M. M. Martens, C. Rabotti, M. Mischi, and R. J. Sluijter, "A robust fetal ECG detection method for abdominal recordings," Physiological Measurement, vol. 28, no. 4, pp. 373-388, 2007.
 R. Vullings, C. H. L. Peters, R. J. Sluijter, M. Mischi, S. G. Oei, and J. W. M. Bergmans, "Dynamic segmentation and linear prediction for maternal ECG removal in antenatal abdominal recordings," Physiological Measurement, vol. 30, no. 3, pp. 291307, 2009.
 E. C. Karvounis, M. G. Tsipouras, D. I. Fotiadis, and K. K. Naka, "An automated methodology for fetal heart rate extraction from the abdominal electrocardiogram," IEEE Transactions on Information Technology in Biomedicine, vol. 11, no. 6, pp. 628638, 2007.
 M. Ungureanu, J. W. M. Bergmans, S. G. Oei, and R. Strungaru, "Fetal ECG extraction during labor using an adaptive maternal beat subtraction technique," Biomedizinische Technik, vol. 52, no. 1, pp. 56-60, 2007.
 G. Clifford, R. Sameni, J. Ward, J. Robinson, and A. J. Wolfberg, "Clinically accurate fetal ECG parameters acquired from maternal abdominal sensors," American Journal of Obstetrics and Gynecology, vol. 205, no. 1, pp. 47.e1-47.e5, 2011.
 W. R. Cohen, S. Ommani, S. Hassan et al., "Accuracy and reliability of fetal heart rate monitoring using maternal abdominal surface electrodes," Acta Obstetricia et Gynecologica Scandinavica, vol. 91, no. 11, pp. 1306-1313, 2012.
 C. Rabotti, M. Mischi, J. O. E. H. Van Laar, G. S. Oei, and J. W. M. Bergmans, "Estimation of internal uterine pressure by joint amplitude and frequency analysis of electrohysterographic signals," Physiological Measurement, vol. 29, no. 7, pp. 829-841, 2008.
 G. Haran, M. Elbaz, M. D. Fejgin, and T. Biron-Shental, "A comparison of surface acquired uterine electromyography and intrauterine pressure catheter to assess uterine activity," American Journal of Obstetrics and Gynecology, vol. 206, no. 5, pp. 412.e1-417.e5, 2012.
 C. Rabotti, M. Mischi, S. G. Oei, and J. W. M. Bergmans, "Noninvasive estimation of the electrohysterographic action-potential conduction velocity," IEEE Transactions on Biomedical Engineering, vol. 57, no. 9, pp. 2178-2187, 2010.
 V Vigneron, A. Azancot, C. Herail, M. Schmidt, O. Sibony, and C. Jutten, "Sensor array for fetal ECG. part 2 : sensor selections," in Proceedings of 2nd International Conference on Computational Intelligence in Medicine and Healthcare (CIMED '05), pp. 434-442, Lisbon, Portugal, June 2005.
 F. Vrins, C. Jutten, and M. Verleysen, "Sensor array and electrode selection for non-invasive fetal electrocardiogram extraction by independent component analysis," in Independent Component Analysis and Blind Signal Separation, vol. 3195, pp. 1017-1024, Springer, Berlin, Germany, 2004.
 B. R. Hayes-gill, C. W. Barratt, and J.-F. Pieri, "Electrode and electrode positioning arrangement for abdominal fetal electrocardiogram detection," U.S. Patent 20 110 306 862, 2011.
 B. Murmann, "Limits on ADC power dissipation," in Analog Circuit Design, M. Steyaert, J. H. Huijsing, and A. H. van Roermund, Eds., pp. 351-367, Springer, Amsterdam, The Netherlands, 2006.
 S. Graczyk, J. Jezewski, J. Wrobel, and A. Gacek, "Abdominal electrohysterogram data acquisition problems and their source
of origin," in Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '95), pp. PS13-PS14, February 1995.
 R. E. Garfield and W. L. Maner, "Physiology and electrical activity of uterine contractions," Seminars in Cell and Developmental Biology, vol. 18, no. 3, pp. 289-295, 2007.
 G. Buonocore, R. Bracci, and M. Weindling, Neonatology: A Practical Approach to Neonatal Diseases, Springer, New York, NY, USA, 2011.
 E. M. Symonds, D. Sahota, and A. Chang, Fetal Electrocardiography, World Scientific, River Edge, NJ, USA, 2001.
 G. S. Dawes, C. R. S. Houghton, C. W. G. Redman, and G. H. A. Visser, "Pattern of the normal human fetal heart rate," British Journal of Obstetrics and Gynaecology, vol. 89, no. 4, pp. 276284, 1982.
 C. Buhimschi, M. B. Boyle, and R. E. Garfield, "Electrical activity of the human uterus during pregnancy as recorded from the abdominal surface," Obstetrics and Gynecology, vol. 90, no. 1, pp. 102-111, 1997.
 R. A. Robergs and R. Landwehr, "The surprising history of the "HRmax=220-age" equation," Journal of Exercise Physiology Online, vol. 5, no. 2, pp. 1-10, 2002.
 B.-U. Kohler, C. Hennig, and R. Orglmeister, "The principles of software QRS detection," IEEE Engineering in Medicine and Biology Magazine, vol. 21, no. 1, pp. 42-57, 2002.
 S. Abboud and D. Sadeh, "Spectral analysis of the fetal electrocardiogram," Computers in Biology and Medicine, vol. 19, no. 6, pp. 409-415, 1989.
 R. Sameni, M. B. Shamsollahi, C. Jutten, and G. D. Clifford, "A nonlinear Bayesian filtering framework for ECG denoising," IEEE Transactions on Biomedical Engineering, vol. 54, no. 12, pp. 2172-2185, 2007.
 G. Friesen, T. C. Jannett, M. Afify Jadallah, S. L. Yates, S. R. Qu int, and H. Troy Nagle, "A comparison of the noise sensitivity of nine QRS detection algorithms," IEEE Transactions on Biomedical Engineering, vol. 37, no. 1, pp. 85-98, 1990.
 B. Bigland-Ritchie, E. F. Donovan, and C. S. Roussos, "Conduction velocity and EMG power spectrum changes in fatigue of sustained maximal efforts," Journal of Applied Physiology Respiratory Environmental and Exercise Physiology, vol. 51, no. 5, pp. 1300-1305, 1981.
 H. de Talhouet and J. G. Webster, "The origin of skin-stretchcaused motion artifacts under electrodes," Physiological Measurement, vol. 17, no. 2, pp. 81-93, 1996.
 R. Vullings, C. Peters, M. Mischi, S. Oei, and J. Bergmans, "Fetal movement quantification by fetal vectorcardiography: a preliminary study," in Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '08), pp. 1056-1059, August 2008.
 T. Oostendorp, Modelling the fetal ECG: een Wetenschappelijke Proeve op het Gebied van de Geneeskunde en Tandheelkunde, Katholieke Universiteit Nijmegen, Nijmegen, The Netherlands, 1989.
 S. Wray, "Uterine contraction and physiological mechanisms of modulation," American Journal of Physiology, vol. 264, no. 1, pp. C1-C18, 1993.
 D. B. Geselowitz, "Dipole theory in electrocardiography," The American Journal of Cardiology, vol. 14, no. 3, pp. 301-306, 1964.
 H. C. Burger and J. B. V Milaan, "Heart-vector and leads," British Heart Journal, vol. 8, no. 3, pp. 157-161, 1946.
 E. H. Moore, "On the reciprocal of the general algebraic matrix," Bulletin of the American Mathematical Society, vol. 26, pp. 394395, 1920.
 R. Penrose, "A generalized inverse for matrices," Mathematical Proceedings of the Cambridge Philosophical Society, vol. 51, pp. 406-413, 1955.
 G. E. Dower, H. B. Machado, and J. A. Osborne, "On deriving the electrocardiogram from vectorcardiographic leads," Clinical Cardiology, vol. 3, no. 2, pp. 87-95, 1980.
 L. Edenbrandt and O. Pahlm, "Vectorcardiogram synthesized from a 12-lead ECG: superiority of the inverse Dower matrix," Journal of Electrocardiology, vol. 21, no. 4, pp. 361-367, 1988.
 G. Moody, "WFDB applications guide (10.5.9)," 2010, http://physionet.org/physiotools/wag/.
 P Zipp, "Effect of electrode parameters on the bandwidth of the surface e.m.g. power-density spectrum," Medical and Biological Engineering and Computing, vol. 16, no. 5, pp. 537-541, 1978.
 C. Rabotti, R. Bijloo, G. Oei, and M. Mischi, "Vectorial analysis of the electrohysterogram for prediction of preterm delivery: a preliminary study," in Proceedings of the 33rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '11), pp. 3880-3883, September 2011.
 C. Rabotti, H. de Lau, N. Haazen, G. Oei, and M. Mischi, "Ultrasound analysis of the uterine wall movement for improved electrohysterographic measurement and modeling," in Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '13), pp. 7436-7439, 2013.
 R. Vullings, M. Mischi, S. Oei, and J. Bergmans, "Novel Bayesian vectorcardiographic loop alignment for improved monitoring of ECG and fetal movement," IEEE Transactions on Biomedical Engineering, vol. 60, no. 6, pp. 1580-1588, 2013.
M. J. Rooijakkers, (1) S. Song, (1) C. Rabotti, (1) S. G. Oei, (2) J. W. M. Bergmans, (1) E. Cantatore, (1) and M. Mischi (1)
(1) Faculty of Electrical Engineering, University of Technology Eindhoven, 5612 AZ Eindhoven, The Netherlands
(2) Maxima Medical Center, 5504 DB Veldhoven, The Netherlands
Correspondence should be addressed to M. J. Rooijakkers; firstname.lastname@example.org
Received 31 May 2013; Revised 13 November 2013; Accepted 15 November 2013; Published 3 February 2014
Academic Editor: Maria Gabriellla Signorini
TABLE 1: Mean [SNR.sub.m] in dB for both simulations. I II III IV V VI Sim. 1 4 cm 5.8 4.8 3.1 3.5 4.6 4.7 8 cm 9.4 4.1 3.6 7.1 9.7 10.8 12 cm 12.4 8.1 4.9 6.3 12.1 12.4 16 cm 14.6 10.9 5.7 7.6 14.3 14.2 20 cm 16.1 12.0 6.2 8.5 15.5 15.8 Sim. 2 4 cm 5.1 4.3 4.5 4.9 5.6 4.8 8 cm 8.0 3.8 5.9 7.8 9.3 9.3 12 cm 10.9 7.8 6.8 9.0 10.8 10.7 16 cm 12.6 9.8 7.9 9.6 12.1 12.0 20 cm 14.3 11.2 0.0 10.4 13.6 13.6 TABLE 2: Mean [SNR.sub.f] in dB for both simulations. I II III Simulation 1 4 cm 1.9 2.2 3.2 8 cm 1.8 2.3 4.2 12 cm 1.5 2.0 4.0 16 cm 1.1 1.8 3.5 20 cm 1.3 1.7 3.4 Simulation 2 4 cm 2.2 2.1 2.0 8 cm 2.7 2.9 2.1 12 cm 3.0 3.3 2.6 16 cm 3.2 3.9 3.0 20 cm 3.4 4.4 3.3 IV V VI Simulation 1 4 cm 3.3 3.2 2.5 8 cm 4.7 4.2 2.0 12 cm 5.5 3.9 1.9 16 cm 6.3 3.2 1.2 20 cm 6.9 3.1 1.4 Simulation 2 4 cm 1.9 2.0 2.1 8 cm 1.7 1.7 1.9 12 cm 1.5 1.8 2.3 16 cm 1.2 1.0 2.0 20 cm 1.4 1.3 2.6 TABLE 3: Mean [SNR.sub.m] in dB for all measurements. I II III Patient 1 4 cm -3.9 3.1 1.4 8 cm 7.2 -0.6 3.5 16 cm 11.9 78 6.3 Patient 2 4 cm 73 8.5 7.5 8 cm 12.5 10.8 11.2 16 cm 11.5 12.7 12.5 Patient 3 4 cm 7.9 7.4 3.3 8 cm 10.5 8.9 6.3 16 cm 12.3 10.8 9.3 Patient 4 4 cm 6.4 4.9 3.3 8 cm 8.7 6.6 8.5 16 cm 14.3 13.2 10.8 Patient 5 4 cm 1.0 -0.6 -1.1 8 cm 6.9 3.4 3.0 16 cm 12.3 11.3 10.1 IV V VI Patient 1 4 cm 0.5 -3.4 -4.2 8 cm 3.0 5.0 -0.2 16 cm 7.4 10.3 11.2 Patient 2 4 cm 3.9 1.9 5.5 8 cm 7.1 7.4 6.4 16 cm 9.9 9.5 11.1 Patient 3 4 cm 1.4 2.3 0.2 8 cm 4.3 6.0 6.2 16 cm 8.7 11.4 12.0 Patient 4 4 cm 2.2 4.6 3.5 8 cm 2.6 8.5 8.8 16 cm 5.7 12.4 11.6 Patient 5 4 cm -0.3 0.8 0.2 8 cm 2.9 6.3 6.6 16 cm 8.5 10.9 11.2 TABLE 4: Mean [SNR.sub.f] in dB for all measurements. I II III Patient 1 4 cm -2.7 1.2 0.9 8 cm 1.1 -1.1 0.8 16 cm 3.7 5.0 4.8 Patient 2 4 cm -0.1 2.0 3.1 8 cm 1.2 3.0 5.9 16 cm -0.7 2.0 5.2 Patient 3 4 cm 0.9 2.7 1.4 8 cm 1.1 4.1 5.2 16 cm -0.7 3.3 5.7 Patient 4 4 cm -0.3 -0.7 -0.2 8 cm -1.7 -2.1 -1.1 16 cm -6.2 -3.5 2.5 Patient 5 4 cm -2.8 -2.1 1.6 8 cm -1.1 -2.3 6.6 16 cm -0.6 -4.5 3.0 IV V VI Patient 1 4 cm 0.4 -1.5 -3.0 8 cm 0.6 0.9 -0.9 16 cm 3.9 0.1 0.1 Patient 2 4 cm 2.2 1.2 -0.8 8 cm 4.2 4.3 5.0 16 cm 6.0 1.0 -2.5 Patient 3 4 cm 0.7 0.3 -2.2 8 cm 4.0 1.9 -0.7 16 cm 7.4 0.6 -3.4 Patient 4 4 cm 0.6 0.3 -0.6 8 cm 1.3 1.6 0.6 16 cm 4.9 0.7 -1.7 Patient 5 4 cm 1.2 1.4 -0.9 8 cm 5.0 6.6 3.3 16 cm 8.5 8.3 3.4
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Rooijakkers, M.J.; Song, S.; Rabotti, C.; Oei, S.G.; Bergmans, J.W.M.; Cantatore, E.; Mischi, M.|
|Publication:||Computational and Mathematical Methods in Medicine|
|Date:||Jan 1, 2014|
|Previous Article:||ADLD: a novel graphical representation of protein sequences and its application.|
|Next Article:||The influence of tone inventory on ERP without focal attention: a cross-language study.|