 Research
 Open Access
 Published:
Reduction of CPR artifacts in the ventricular fibrillation ECG by coherent line removal
BioMedical Engineering OnLine volume 9, Article number: 2 (2010)
Abstract
Background
Interruption of cardiopulmonary resuscitation (CPR) impairs the perfusion of the fibrillating heart, worsening the chance for successful defibrillation. Therefore ECGanalysis during ongoing chest compression could provide a considerable progress in comparison with standard analysis techniques working only during "handsoff" intervals.
Methods
For the reduction of CPRrelated artifacts in ventricular fibrillation ECG we use a localized version of the coherent line removal algorithm developed by Sintes and Schutz. This method can be used for removal of periodic signals with sufficiently coupled harmonics, and can be adapted to specific situations by optimal choice of its parameters (e.g., the number of harmonics considered for analysis and reconstruction). Our testing was done with 14 different human ventricular fibrillation (VF) ECGs, whose fibrillation band lies in a frequency range of [1 Hz, 5 Hz]. The VFECGs were mixed with 12 different ECGCPRartifacts recorded in an animal experiment during asystole. The length of each of the ECGdata was chosen to be 20 sec, and testing was done for all 168 = 14 × 12 pairs of data. VFtoCPR ratio was chosen as 20 dB, 15 dB, 10 dB, 5 dB, 0 dB, 5 dB and 10 dB. Here 20 dB corresponds to the highest level of CPRartifacts.
Results
For nonoptimized coherent line removal based on signals with a VFtoCPR ratio of 20 dB, 15 dB, 10 dB, 5 dB and 0 dB, the signaltonoise gains (SNRgains) were 9.3 ± 2.4 dB, 9.4 ± 2.4 dB, 9.5 ± 2.5 dB, 9.3 ± 2.5 dB and 8.0 ± 2.7 (mean ± std, n = 168), respectively. Characteristically, an original VFtoCPR ratio of 10 dB, corresponds to a variance ratio var(VF):var(CPR) = 1:10. An improvement by 9.5 dB results in a restored VFtoCPR ratio of 0.5 dB, corresponding to a variance ratio var(VF):var(CPR) = 1:1.1, the variance of the CPR in the signal being reduced by a factor of 8.9.
Discussion
The localized coherent line removal algorithm uses the information of a single ECG channel. In contrast to multichannel algorithms, no additional information such as thorax impedance, blood pressure, or pressure exerted on the sternum during CPR is required. Predictors of defibrillation success such as mean and median frequency of VFECGs containing CPRartifacts are prone to being governed by the harmonics of the artifacts. Reduction of CPRartifacts is therefore necessary for determining reliable values for estimators of defibrillation success.
Conclusions
The localized coherent line removal algorithm reduces CPRartifacts in VFECG, but does not eliminate them. Our SNRimprovements are in the same range as offered by multichannel methods of Rheinberger et al., Husoy et al. and Aase et al. The latter two authors dealt with different ventricular rhythms (VF and VT), whereas here we dealt with VF, only. Additional developments are necessary before the algorithm can be tested in real CPR situations.
Background
Frequent interruptions of chest compressions (CC) as part of cardiopulmonary resuscitation (CPR) during ventricular fibrillation (VF) and pulseless ventricular tachycardia (VT) impair myocardial perfusion and worsen the chance for successful defibrillation with stable return of spontaneous circulation [1, 2]. Eilevstjonn et al. reported that "noflow times" (NFT) comprise about 50% of time during resuscitation [3], and gave suggestions on how to reduce NFT. On the other hand, analysis of the ECG for fibrillation detection necessarily requires interruption of CC, at least with the ECGanalysis algorithms currently implemented in defibrillators available on the market.
Therefore ECGanalysis during ongoing chest compression can provide a considerable progress in comparison with standard analysis techniques working only during "handsoff" intervals [4]. These intervals have recently been found to be unexpectedly long and harmful [1, 5, 6]. In addition to avoiding "handsoff" times, new analysis techniques could become the prerequisite for prediction of defibrillation success probability [7–18]. These analysis techniques would allow to avoid unpromising and therefore ultimately damaging defibrillator shocks.
Aase, Husoy, Eilevstjonn et al. [19–22] have developed adaptive filtering approaches to realtime separation of VF/VT and CPR for a multichannelcontext. Berger et al. suggested an adaptive noise cancellation technique [23] and recently Kalman filtering techniques were used [4, 24]. We note also the contributions by Aramendi et al. [25] and Irusta et al. [26]. In the present work, we concentrate again on timefrequency methods [27, 28] and on the situation where only one ECG channel is available, without any additional information concerning blood pressure or concerning the pressure on the sternum applied during resuscitation. This is particularly adapted to the current use of automated external defibrillators (AEDs), where  so far  no such additional information is available.
In this work we present a method based on a timefrequency analysis. The method makes use of a windowed Fourier transform that captures characteristic features of VF signals and CPR artifacts. A commonly used criterion to assess an algorithm is the improvement of the signaltonoise ratio (SNR). The SNR is expressed as the variance of the "proper" VFECG without CPRartifacts (signal) divided by the variance of the CPRartifacts (noise) in the ECG. Coherent line removal can be used for reduction of CPRrelated artifacts in VF signals, because the fibrillation ECG does not contain a line spectrum (at some particular frequency) which would unintentionally be removed by the algorithm, but many different frequencies which are continuously distributed in a "fibrillation (frequency) band".
Mere improvement of signaltonoise ratio is not a guarantee for a better estimate of the "proper" VFECG. It is of additional importance, that typical ECGbased parameters like the median frequency (of the ECG) show similar results in estimate_{VF}(t) as compared to ECG_{VF}(t). Furthermore, it is important that an artifactfree ECGsignal (containing VF only) shows approximately the same median frequency before and after application of the CPRfiltering algorithm. Median frequency is considered to be an interesting parameter for prediction of defibrillation success [8, 10, 29], even though not unequivocally [30, 31].
In the present work, we illustrate by examples how the proposed method affects the median frequency and present numerical results for the SNRimprovement.
Methods
(A) Data
We used one exemplary dataset from a VFexperiment in a pig model for illustration of the effect of our CPRreduction algorithm on the power spectrum. Data in this animal experiment were recorded with 12 bit and 1000 Hz sampling frequency. For the present purpose, we used a downsampled version of the ECG recorded with 200 Hz sampling frequency.
The actual testing of our CPRreduction algorithm was done with 14 different human ventricular fibrillation (VF) ECGs, which were mutually mixed with 12 different ECGCPRartifacts recorded in an animal experiment during asystole with an applied CPRfrequency between 80/min and 120/min. The length of each of the ECGdata was chosen to be 20 sec, and testing was done for all 168 = 14 × 12 pairs of data.
The 14 different human ECGs have been collected using a Welch Allyn PIC 50 defibrillator, recorded with 12 bit and 375 Hz sampling frequency during real outofhospital CPR situations. For the present purpose we chose human ECGs with frequencies lying (roughly) in the range [1 Hz, 5 Hz].
The pig experiments were conducted according to Utsteinstyle guidelines [32] and approved by the Federal Austrian Animal Experiment Committee. The recording of the human data was approved by the local Ethics Committee of Innsbruck Medical University.
(B) Data processing and quality assessment for CPRreduction algorithms
Data were processed using MATLAB (The Mathworks, Natick (MA), version R2007b). For computation of mean frequency, median frequency and dominant frequency the upper cutoff frequency was 30 Hz. The lower cutoff frequency was 4.33 Hz (for ECGdata including CPRartifacts) and 2.2 Hz (for ECGdata purged from CPRartifacts).
For spectrograms shown in Figures the frequency range has been restricted to [0 Hz, 15 Hz] to improve visibility (the spectrogram above 15 Hz does not contain much interesting information). The "typical frequencies" corresponding to ventricular fibrillation ECG are called the "fibrillation band". This "fibrillation band" changes in the course of time. For human VFECGs, the fibrillation band typically is within the frequency window [1 Hz, 6 Hz].
In an ECGsignal
composed of a CPRartifactfree VFECG and an ECG containing only CPRrelated artifacts, the SNR is defined as
where var(ECG_{VF}) is the variance of the proper VFECG and var(ECG_{CPR}) is the variance of the CPRrelated artifacts in the ECG. The acronym SNR refers to "signaltonoise ratio", where the signal is ECG_{VF} and the noise is ECG_{CPR}. Usually it is expressed in decibel units (dB), i.e., as 10 times the logarithm (with basis 10) of this SNR. For very strong artifacts, the signaltonoise ratio is around 20 dB, whereas 0 dB corresponds to rather small CPRrelated artifacts.
An algorithm for extraction of CPRartifacts in an ECG gives rise to a decomposition
The quality of such an algorithm can be assessed by looking at the variance of error
or by considering the restored signaltonoise ratio
In particular, we use the SNRgain, i.e., the difference (rSNR  SNR).
(C) Coherent line removal
For the reduction of CPRrelated artifacts we used the coherent line removal algorithm developed by Sintes and Schutz [28, 33, 34]. This method can be used for removal of periodic signals with sufficiently strong harmonics, and is presented in more detail in the Appendix. Before applying coherent line removal to a signal s = s(t), the approximate CPRfrequency is estimated from the sum of power at (each) frequency f and its harmonics, taking the frequency f _{0} at which the following function maximizes:
The function ŝ = ŝ(f) is the Fourier transform of the ECGsignal s = s(t). Furthermore f is some chosen frequency, and (k f), k = 2,3, ..., M, are its harmonics. An example, plotting this function for an ECG containing CPRartifacts is shown in Fig 1.
Coherent line removal has a few parameters which can be adapted to the signal in question:

the time window,

the number har of harmonics (default = 4),

the frequency width delta for line removal (default = 4),

and the number of harmonics NumHar considered for reconstruction (default = 10),
We choose a time window of 2048 data points, corresponding to 10.24 sec (at a sampling frequency of 200 Hz), and determined the optimal choice of har, delta and NumHar in a grid search, varying har and NumHar from 2 to 10, and delta from 2 to 5 (or took default values). The "grid search" simply computes the variance of the error var(ECG_{CPR}  estimate_{CPR}) and looks for the smallest error (among all different sets of parameters) for our 168 = 14 × 12 datasets. Minimization of var(ECG_{CPR}  estimate_{CPR}) is equivalent to maximum SNR improvement.
Results
Fig 2 shows the windowed Fourier transform of an ECG for a VFexperiment in a pig model. Fibrillation in this example starts at ~80 sec, and cardiopulmonary resuscitation (CPR) after ~310 sec. CPR is performed with approximately ~110/min = ~1.8 Hz. The respective CPRrelated artifacts at ~1.8 Hz are clearly visible in the spectrogram, together with the harmonics at ~3.7 Hz, ~5.5 Hz etc. The mean and median frequency computed for the frequency window [4.33 Hz, 30 Hz] are good parameters for the "fibrillation band" as long as no CPR is performed [8, 10]. As soon as CPR starts, the mean and median frequency do not reflect the course of the "fibrillation band" any more. The "fibrillation band" is increasing after start of CPR, whereas the mean and median frequency are decreasing after start of CPR. Starting from ~1000 sec in this example, the median frequency is completely dominated by the second harmonic of the CPRartifacts at ~5.5 Hz.
Fig 3 shows the windowed Fourier transform of the ECG in the same pig experiment, but now after purging the ECG from CPRrelated artifacts using localized coherent line removal. The mean and median frequency can now be computed with respect to the extended frequency window [2 Hz, 30 Hz] and follow the "fibrillation band" apart from the last stage of the experiment (> 1100 sec), where no "fibrillation band" is visible any more (probably due to asystole). Fig 4 also shows estimated VFtoCPR ratios which typically are in a range of 15 dB to 5 dB.
In Fig 4, the median VFECG frequency is shown for the original ECG as compared with the median frequency of the VFECG purged from CPRrelated artifacts by coherent line removal. Before start of CPR, these two median frequencies coincide almost perfectly which indicates that coherent line removal does not distort VFECG when the latter does not contain CPRrelated artifacts. After start of CPR, the two median frequencies are rather different and visual inspection shows that the median frequency of the original VFECG (containing CPRrelated artifacts) is governed by the harmonics of the CPRartifacts. In the present example, the second harmonic at 5.5 Hz is particularly dominating. The median frequency of the VFECG purged from CPRartifacts by coherent line removal, is not governed by the harmonics of CPRrelated artifacts, but by the "fibrillation band" which reflects the ventricular fibrillation of the heart of the animal investigated.
For Fig 5 and Table 1, results from human outofhospital VFECGs are shown, which were mixed with CPRrelated ECG artifacts. In Fig 5, the respective VFtoCPR ratio was chosen to be 10 dB. The reconstructed CPRECG is shown for optimal SNRimprovement, i.e., use of optimized parameters for the coherent line removal.
In this particular example, the human ECG mixed with CPRECG had a "fibrillation band" in the frequency range [1 Hz, 5 Hz], which is expected to be difficult to separate from CPRrelated artifacts with ~1.8 Hz and harmonics at ~3.7 Hz, ~5.5 Hz etc.
Table 1 refers to 14 different human outofhospital VFECG datasets, mutually mixed with 12 different CPRECGs, but now using different choices of VFtoCPR ratios of 20 dB, 15 dB, 10 dB, 5 dB, 0 dB, 5 dB, and 10 dB. The optimal SNRgains (for optimized parameters of the coherent line removal algorithm) were 10.3 ± 2.0 dB, 10.5 ± 2.2 dB, 10.5 ± 2.6 dB, 10.2 ± 2.8 dB, 8.9 ± 2.9 dB, 6.3 ± 3.7 dB, and 1.7 ± 4.0 dB (mean ± std, n = 168) in the example used. Optimization by a grid search was done with respect to the parameters which can be chosen for implementation of the coherent line removal algorithm (such as the number of harmonics used for the estimation of the CPRrelated artifacts in the ECG). If the algorithm is not optimized, but used with the fixed default values of the parameters, the SNRgains are 9.3 ± 2.4 dB, 9.4 ± 2.4 dB, 9.5 ± 2.5 dB, 9.3 ± 2.5 dB, 8.0 ± 2.7 dB, 4.9 ± 3.7 dB, and 1.0 ± 4.0 dB.
Discussion
An important point in the development of CPRartifact reduction algorithms is the use of human VFECG data with typical frequencies lying in the same range as the CPRartifacts themselves. This is not the case in animal experimental data: there often the frequencies of ventricular fibrillation are much higher than the frequencies of CPRartifacts, simplifying the separation of VFECG and CPRartifacts. Here, we deliberately took human outofhospital VFECG data with a frequency range of about [1 Hz, 5 Hz], i.e., much lower than in animal data. This frequency range matches well with the frequency range of cardiopulmonary resuscitation of [1.3 Hz, 2 Hz] (= [80/min, 120/min]) and its harmonics. These are compression frequencies near those given in the present resuscitation guidelines (100/min). The animal VFECG as used in Figs 1, 2, 3 and 4 served merely the purpose to clearly illustrate the strategy of the coherent line removal method.
The VFtoCPR ratio (expressed in dB) is negative for very strong CPRrelated artifacts and positive, if the CPRrelated artifacts are very weak. In case of typical VFtoCPR ratios which range between 15 dB to 5 dB, the coherent line removal algorithm achieves an SNRimprovement of ~9.5.
A typical example would be an original VFtoCPR ratio of 15 dB, corresponding to a variance ratio var(VF):var(CPR) = 1:31.6. An improvement by 9.5 dB would result in a restored VFtoCPR ratio of 5.5 dB, corresponding to a variance ratio var(VF):var(CPR) = 1:3.5, the variance of the CPR in the signal being reduced by a factor of 8.9.
Another typical example would be an original VFtoCPR ratio of 10 dB, corresponding to a variance ratio var(VF):var(CPR) = 1:10. An improvement by 9.5 dB would result in a restored VFtoCPR ratio of 0.5 dB, corresponding to a variance ratio var(VF):var(CPR) = 1:1.1, the variance of the CPR in the signal being reduced by a factor of 8.9.
Many studies on CPR artifacts removal used SNR improvement to measure the quality of CPR artifact suppression [19, 21, 24, 35–39]. All these studies were based on the additive data model, i.e., adding a pure artifact signal to the artifactfree ECG, that we adopted. The SNR improvement that we achieved on our testing data is close to that obtained in studies that were based on a twochannel setting, i.e., including additional information such as blood pressure [24, 39] and/or compression depth [19, 21]. For instance, at the SNR level of 10 dB, we obtained the average SNR improvement of 10.5 dB which is smaller than 12.4 dB and 11.84 dB obtained in refs [19] and [39], respectively, but superior to 10.2 dB obtained in [24]. None of the previous studies has analyzed the SNR improvement for SNR level smaller than 10 dB, i.e., for very large artifacts. It is interesting that in our study the SNR improvement for the SNR levels 15 dB and 20 dB did not increase further and stayed in the same range as the SNR improvement obtained for 5 dB and 10 dB, see Table 1. In contrast, SNR improvement decreases with increasing SNR levels that exceed 5 dB as it was also reported in the previous studies stated above.
Here we used a time window of about 10 sec for application of the coherent line removal algorithm (i.e., typically 2 timewindows for a dataset of 20 sec length), which is the reason that single strong CPRartifacts (with a duration of ~0.6 sec, as shown in Fig 4 in the time window [21 sec, 22 sec]) were not eliminated.
During performance of CPR, automatic analysis of the ECG is difficult, because the predominant part of the signal consists of CPRartifacts. This is illustrated in Fig 2 by computing the parameters mean frequency, the median frequency and the dominant frequency in the frequency window [4.33 Hz, 30 Hz] [8, 10, 29, 30]. Since this frequency window excludes the CPRartifacts at ~1.8 Hz and its first harmonic ~3.7 Hz, one might be tempted to consider these parameters (such as the median frequency) to be a good indicator of the "fibrillation band" [7, 8, 11, 16, 17, 40–45]. This is not the case. As illustrated in Fig 2, only before start of CPR the mean frequency (blue), median frequency (magenta) and dominant frequency (black) in Fig 2 follow very well the "fibrillation band". After start of CPR, these parameters are much more governed by the CPRartifacts than by the "fibrillation band" itself. The dominant frequency, in particular, coincides with the second harmonic of CPRartifacts at ~5.5 Hz during the time period [~310 sec, 1410 sec]. During this time period the "fibrillation band" changes dramatically. Nevertheless the dominant frequency is entirely independent of these changes in the "fibrillation band". For the mean frequency and the median frequency, the situation is not so pronounced, but still very unsatisfactory. The median frequency, for example, coincides with the second harmonic of CPRartifacts during the time period [~908 sec, ~1410 sec]. We therefore conclude that parameters such as mean, median and dominant frequency are blurred by CPRartifacts and are therefore not particularly adapted to serve as quantitative indicators for the location of the "fibrillation band".
Fig 3 shows the same pig experiment as Fig 2 again, but now the ECG has been purged from CPRartifacts by local coherent line removal. Again the mean frequency (blue), the median frequency (magenta) and the dominant frequency (black) are shown. The dominant frequency still does not seem to be helpful. The mean and median frequency, on the other hand, now follow closely the "fibrillation band" and can therefore be considered as parameters describing the location of the fibrillation band for the time period [~80 sec, ~1100 sec]. Later than 1100 sec the "fibrillation band" in this pig experiment does not seem to exist any more due to asystole of the heart of the animal. The lower panel of Fig 3 shows that the proportion of the power in the frequency band [0.33 Hz, 2 Hz] as compared to the frequency band [0.33 Hz, 30 Hz] is much lower after filtering of CPRartifacts than without it. For a large time period [~80 sec, ~1000 sec] this proportion is only about 10%  20% (instead of 80%  90%).
Fig 4 shows a comparison (again for the same pig experiment) of the median frequencies before and after CPRfiltering. Observe, in particular, that during the time period [~80 sec, ~310 sec], during which CPR is not performed, there is almost a perfect match between the median frequencies computed with/without CPRfiltering. This is to say that the filtering algorithm used (namely localized coherent line removal) does not remove anything from the "fibrillation band". During the time period [~80 sec ~1400 sec] where CPR is performed, the median frequency after CPRfiltering (magenta line in Fig 4) follows much better the "fibrillation band" than the median frequency before CPRfiltering (black line in Fig 4).
Optimization of the parameters of the coherent line removal algorithm do not give much better results than just using default values (see Table 1). We consider this to be an important result, showing that optimization of the coherent line removal algorithm is not important. Also the default values like har and NumHar could be changed as long as they are not obviously too small. All values for har and NumHar larger than our default values would give results rather equal in quality.
When dealing with CPR artefacts of real outofhospital cardiac arrests a wide variety of compressions rates are found. Coherent line removal as presented here depends on constant compression rate within the respective chosen timewindow.
It is common standard that, when introducing a new method, simple data models (such as additive simulations of CPR artefacts) are used for preliminary testing. Methods need to perform adequately on such simplified data before extending testing on a heterogenous dataset. Most of the authors hitherto have used additive models of VF. Also, to the best of our knowledge, the SNR is the basic criterion for classifying the quality of a CPRremoval algorithm. In our recently published paper [46], we investigated the efficiency of various twochannel methods for CPRartefact removal in nonshockable rhythms: for such ECGs, twochannel methods could not reduce CPR artefacts without affecting the rhythm analysis for shock recommendation. For the coherent line removal algorithm the same is true, i.e., the efficiency and quality of CPRartefact removal for nonshockable rhythms is not satisfactory. Possibly a better understanding of the spectral distribution of rhythms other than VF is needed for future adaption of the coherent line algorithm.
From a computational point of view, the coherent line removal algorithm is very fast. It could probably be implemented on a Coldfire processor containing a floating point unit. Such Coldfire processors are commercially available, but only when ordering large amounts (> 10000 units). The coherent line removal algorithm is, in particular, much faster than the Kalman filter algorithm presented in ref [24].
Conclusions
The localized coherent line removal algorithm reduces CPRartifacts in VFECG, but does not eliminate them. It uses the ECGchannel only, without any additional information (like blood pressure). Our SNRimprovements are in the same range as offered by multichannel methods of Rheinberger et al. [24], Husoy et al. [20] and Aase et al. [21]. In refs. [20, 21] the authors dealt with different rhythms (VF and VT) whereas we dealt with VF, exclusively. Additional developments are necessary before the algorithm can be tested in real CPR situations.
Appendix: Short description of the coherent line removal algorithm
In Ref [33], the coherent part y(t) of a signal s(t) = y(t) + n(t) is described as
where the overbar denotes complex conjugation and where α _{ k }are appropriate coefficients. Here m(t) is a nearly monochromatic signal,
with slowly varying amplitude r(t) and frequency f _{0}(t), and with . The number M of considered harmonics is one of the parameters of the coherent line removal algorithm.
In the present application, the coherent part y(t) of the ECGsignal is identified with its CPRartifacts, corresponding to the function estimate_{CPR}(t) in the Methods Section. The frequency f _{0}(t) is the frequency of resuscitation (e.g., 110/min = 1.83 Hz).
We use a localized version of the approach from Ref [33], using timewindows of 10.24 sec length. Localization is necessary, because CPRrelated artifacts change in amplitude and frequency, and may disappear at all when no resuscitation is performed. Localization means that the coefficients α _{ k }depend on the specific timewindow, e.g., [0 sec, 10.24 sec], [10.24 sec, 20.48 sec], [20.48 sec, 30.72 sec] etc. If no resuscitation artifacts are present in a timewindow, then we have α _{ k }= 0 for all harmonics k = 1, 2, ..., M.
The stepinfrequency is given by
The width of the frequency window used for the estimation of the parameters α _{ k }is related to the parameter δ by
For our default value δ = 4, the width of the frequency window is therefore given as
In the frequency space, the decomposition s(t) = y(t) + n(t) is given as
where the hat ŝ = ŝ(v) indicates the Fourier transform (depending on frequency v) of a function s = s(t). Here n = n(t) corresponds to the proper VFECG without the CPRrelated artifacts. Restricting oneself to the frequency window around the k^{th}harmonic, one gets
Applying the inverse Fourier transform, this leads to
Using the abbreviations
and
one gets
Considering the VFECG n(t) and the respective restrictions n _{ k }(t) to a frequency window around the k^{th}harmonic as a stochastic processes with ensemble mean value zero, ⟨n _{ k }(t)⟩ = 0, leading to
Hence multiplication of the stochastic processes B _{ k }(t) by appropriate scalars Γ_{ k }leads to stochastic processes
which all have the same ensemble mean value
The values of the scalars Γ_{ k }can be obtained from a least square method, comparing the first harmonic with the other harmonics considered, and taking Γ_{ k }as the scalar which leads to the minimal value of the ensemble expectation . This is achieved by setting
From the values of Γ_{ k }one may compute the values for the α _{ k }, namely α _{ k }= (α _{1}Γ_{ k })^{k}.
Finally the interference m(t) is reconstructed as a linear combination
such that it has the same mean and minimum variance. Here . Note that the number M of harmonics for analysis may be different from the number N of harmonics for reconstruction. The parameters M and N may be used for optimization of the coherent line removal algorithm. This leads to
which allows to estimate the coherent part of the original signal.
References
 1.
Steen S, Liao Q, Pierre L, Paskevicius A, Sjoberg T: Evaluation of LUCAS, a new device for automatic mechanical compression and active decompression resuscitation. Resuscitation 2002, 55: 285–299. 10.1016/S03009572(02)00271X
 2.
Wik L, KramerJohansen J, Myklebust H, Sorebo H, Svensson L, Fellows B, Steen PA: Quality of cardiopulmonary resuscitation during outofhospital cardiac arrest. Jama 2005, 293: 299–304. 10.1001/jama.293.3.299
 3.
Eilevstjonn J, KramerJohansen J, Eftestol T, Stavland M, Myklebust H, Steen PA: Reducing no flow times during automated external defibrillation. Resuscitation 2005, 67: 95–101. 10.1016/j.resuscitation.2005.04.009
 4.
Ruiz de Gauna S, Ruiz J, Irusta U, Aramendi E, Eftestol T, KramerJohansen J: A method to remove CPR artefacts from human ECG using only the recorded ECG. Resuscitation 2008, 76: 271–278. 10.1016/j.resuscitation.2007.08.002
 5.
Sato Y, Weil MH, Sun S, Tang W, Xie J, Noc M, Bisera J: Adverse effects of interrupting precordial compression during cardiopulmonary resuscitation. Crit Care Med 1997, 25: 733–736. 10.1097/0000324619970500000005
 6.
Abella BS, Alvarado JP, Myklebust H, Edelson DP, Barry A, O'Hearn N, Hoek TL, Becker LB: Quality of cardiopulmonary resuscitation during inhospital cardiac arrest. Jama 2005, 293: 305–310. 10.1001/jama.293.3.305
 7.
Achleitner U, Wenzel V, Strohmenger HU, Krismer AC, Lurie KG, Lindner KH, Amann A: The effects of repeated doses of vasopressin or epinephrine on ventricular fibrillation in a porcine model of prolonged cardiopulmonary resuscitation. Anesth Analg 2000, 90: 1067–1075. 10.1097/0000053920000500000012
 8.
Amann A, Rheinberger K, Achleitner U: Algorithms to analyze ventricular fibrillation signals. Curr Opin Crit Care 2001, 7: 152–156. 10.1097/0007519820010600000003
 9.
Eftestol T, Sunde K, Aase SO, Husoy JH, Steen PA: "Probability of successful defibrillation" as a monitor during CPR in outofhospital cardiac arrested patients. Resuscitation 2001, 48: 245–254. 10.1016/S03009572(00)002665
 10.
Amann A, Rheinberger K, Achleitner U, Krismer AC, Lingnau W, Lindner KH, Wenzel V: The prediction of defibrillation outcome using a new combination of mean frequency and amplitude in porcine models of cardiac arrest. Anesth Analg 2002, 95: 716–722. 10.1097/0000053920020900000034
 11.
Amann A, Achleitner U, Antretter H, Bonatti JO, Krismer AC, Lindner KH, Rieder J, Wenzel V, Voelckel WG, Strohmenger HU: Analysing ventricular fibrillation ECGsignals and predicting defibrillation success during cardiopulmonary resuscitation employing N(alpha)histograms. Resuscitation 2001, 50: 77–85. 10.1016/S03009572(01)003227
 12.
Watson JN, Uchaipichat N, Addison PS, Clegg GR, Robertson CE, Eftestol T, Steen PA: Improved prediction of defibrillation success for outofhospital VF cardiac arrest using wavelet transform methods. Resuscitation 2004, 63: 269–275. 10.1016/j.resuscitation.2004.06.012
 13.
Yang Z, Yang Z, Lu W, Harrison RG, Eftestol T, Steen PA: A probabilistic neural network as the predictive classifier of outofhospital defibrillation outcomes. Resuscitation 2005, 64: 31–36. 10.1016/j.resuscitation.2004.07.002
 14.
Eftestol T, Losert H, KramerJohansen J, Wik L, Sterz F, Steen PA: Independent evaluation of a defibrillation outcome predictor for outofhospital cardiac arrested patients. Resuscitation 2005, 67: 55–61. 10.1016/j.resuscitation.2005.05.006
 15.
Eftestol T, Wik L, Sunde K, Steen PA: Effects of cardiopulmonary resuscitation on predictors of ventricular fibrillation defibrillation success during outofhospital cardiac arrest. Circulation 2004, 110: 10–15. 10.1161/01.CIR.0000133323.15565.75
 16.
Neurauter A, Eftestol T, KramerJohansen J, Abella BS, Sunde K, Wenzel V, Lindner KH, Eilevstjonn J, Myklebust H, Steen PA, Strohmenger HU: Prediction of countershock success using single features from multiple ventricular fibrillation frequency bands and feature combinations using neural networks. Resuscitation 2007, 73: 253–263. 10.1016/j.resuscitation.2006.10.002
 17.
Neurauter A, Strohmenger HU: Prediction of countershock success employing single features from multiple ventricular fibrillation frequency bands and feature combinations using neural networks. Resuscitation 2008, 76: 152. 10.1016/j.resuscitation.2007.06.022
 18.
Gundersen K, Kvaloy JT, KramerJohansen J, Eftestol T: Identifying approaches to improve the accuracy of shock outcome prediction for outofhospital cardiac arrest. Resuscitation 2008, 76: 279–284. 10.1016/j.resuscitation.2007.07.019
 19.
Langhelle A, Eftestol T, Myklebust H, Eriksen M, Holten BT, Steen PA: Reducing CPR artefacts in ventricular fibrillation in vitro. Resuscitation 2001, 48: 279–291. 10.1016/S03009572(00)002598
 20.
Husoy JH, Eilevstjonn J, Eftestol T, Aase SO, Myklebust H, Steen PA: Removal of cardiopulmonary resuscitation artifacts from human ECG using an efficient matching pursuitlike algorithm. IEEE Trans Biomed Eng 2002, 49: 1287–1298. 10.1109/TBME.2002.804591
 21.
Aase SO, Eftestol T, Husoy JH, Sunde K, Steen PA: CPR artifact removal from human ECG using optimal multichannel filtering. IEEE Trans Biomed Eng 2000, 47: 1440–1449. 10.1109/10.880095
 22.
Eilevstjonn J, Eftestol T, Aase SO, Myklebust H, Husoy JH, Steen PA: Feasibility of shock advice analysis during CPR through removal of CPR artefacts from the human ECG. Resuscitation 2004, 61: 131–141. 10.1016/j.resuscitation.2003.12.019
 23.
Berger RD, Palazzolo J, Halperin H: Rhythm discrimination during uninterrupted CPR using motion artifact reduction system. Resuscitation 2007, 75: 145–152. 10.1016/j.resuscitation.2007.03.007
 24.
Rheinberger K, Steinberger T, Unterkofler K, Baubin M, Klotz A, Amann A: Removal of CPR artifacts from the ventricular fibrillation ECG by adaptive regression on lagged reference signal. IEEE Trans Biomed Eng 2008, 55: 130–137. 10.1109/TBME.2007.902235
 25.
Aramendi E, Ruiz de Gauna S, Irusta U, Ruiz J, Arcocha MF, Ormaetxe JM: Detection of ventricular fibrillation in the presence of cardiopulmonary resuscitation artefacts. Resuscitation 2007, 72: 115–123. 10.1016/j.resuscitation.2006.05.017
 26.
Irusta U, Ruiz J, de Gauna SR, Eftestol T, KramerJohansen J: A least meansquare filter for the estimation of the cardiopulmonary resuscitation artifact based on the frequency of the compressions. IEEE Trans Biomed Eng 2009, 56: 1052–1062. 10.1109/TBME.2008.2010329
 27.
Klotz A, Feichtinger HG, Amann A: Elimination of CPRartefacts in VFECGs by time frequency methods. Biomedizinische Technik (Biomedical Engineering) 2003, (Suppl 48):218–219. 10.1515/bmte.2003.48.s1.218
 28.
Klotz A, Feichtinger HG, Amann A: Removal of CPR artifacts in ventricular fibrillation ECG by local coherent line removal. Proc of the EUSIPCO 2004, 2203–2206.
 29.
Achleitner U, Wenzel V, Strohmenger HU, Lindner KH, Baubin MA, Krismer AC, Mayr VD, Amann A: The beneficial effect of basic life support on ventricular fibrillation mean frequency and coronary perfusion pressure. Resuscitation 2001, 51: 151–158. 10.1016/S03009572(01)003884
 30.
Stadlbauer KH, Rheinberger K, Wenzel V, Raedler C, Krismer AC, Strohmenger HU, Augenstein S, WagnerBerger HG, Voelckel WG, Lindner KH, Amann A: The effects of nifedipine on ventricular fibrillation mean frequency in a porcine model of prolonged cardiopulmonary resuscitation. Anesth Analg 2003, 97: 226–230. table of contents 10.1213/01.ANE.0000068801.28430.ED
 31.
Holzer M, Behringer W, Sterz F, Kofler J, Oschatz E, Schuster E, Laggner AN: Ventricular fibrillation median frequency may not be useful for monitoring during cardiac arrest treated with endothelin1 or epinephrine. Anesth Analg 2004, 99: 1787–1793. table of contents 10.1213/01.ANE.0000138421.74434.E3
 32.
Idris AH, Becker LB, Ornato JP, Hedges JR, Bircher NG, Chandra NC, Cummins RO, Dick W, Ebmeyer U, Halperin HR, Hazinski MF, Kerber RE, Kern KB, Safar P, Steen PA, Swindle MM, Tsitlik JE, von Planta I, von Planta M, Wears RL, Weil MH: Utsteinstyle guidelines for uniform reporting of laboratory CPR research. A statement for healthcare professionals from a Task Force of the American Heart Association, the American College of Emergency Physicians, the American College of Cardiology, the European Resuscitation Council, the Heart and Stroke Foundation of Canada, the Institute of Critical Care Medicine, the Safar Center for Resuscitation Research, and the Society for Academic Emergency Medicine. Resuscitation 1996, 33: 69–84. 10.1016/S03009572(96)010556
 33.
Sintes AM, Schutz BF: Coherent line removal: filtering out harmonically related line interference from experimental data, with applications to gravitational wave detectors. Phys Rev D 1998, 58: 122003. 10.1103/PhysRevD.58.122003
 34.
Levkov C, Mihov G, Ivanov R, Daskalov I, Christov I, Dotsinsky I: Removal of powerline interference from the ECG: a review of the subtraction procedure. Biomed Eng Online 2005, 4: 50. 10.1186/1475925X450
 35.
Ruiz J, Aramendi E, Ruiz de Gauna S, Lazkano A, Leturiendo LA, Gutierrez JJ: Ventricular fibrillation detection in ventricular fibrillation signals corrupted by cardiopulmonary resuscitation artefacts. Computers in Cardiology 2003, 31: 221–224. full_text
 36.
Ruiz de Gauna S, Ruiz J, Irusta U, Aramendi E, Lazkano A, Gutierrez JJ: CPR artefact removal from VG signals by means of an adaptive Kalman filter using the chest compression frequency as a reference signal. Computers in Cardiology 2005, 32: 175–178. full_text
 37.
Irusta U, Ruiz J, Ruiz de Gauna S, Aramendi A, Lazkano A, Gutierrez JJ: A variable step size LMS alogrithm for the suppression of CPR artefact from a VF signal. Computers in Cardiology 2005, 32: 179–182. full_text
 38.
Aramendi E, Ruiz J, Ruiz de Gauna S, Irusta U, Lazkano A, Gutierrez JJ: A simple effective filtering method for removing CPR caused artefacts from surface ECG signals. Computers in Cardiology 2005, 32: 547–550. full_text
 39.
Werther T, Klotz A, Kracher G, Baubin M, Feichtinger H, Gilly H, Amann A: CPR artifact removal in ventricular fibrillation ECG signals using Gabor multipliers. IEEE Trans Biomed Eng 2009, 56: 320–327. 10.1109/TBME.2008.2003107
 40.
Dzwonczyk R, Brown CG, Werman HA: The median frequency of the ECG during ventricular fibrillation: its use in an algorithm for estimating the duration of cardiac arrest. IEEE Trans Biomed Eng 1990, 37: 640–646. 10.1109/10.55668
 41.
Strohmenger HU, Lindner KH, Keller A, Lindner IM, Pfenninger E, Bothner U: Effects of graded doses of vasopressin on median fibrillation frequency in a porcine model of cardiopulmonary resuscitation: results of a prospective, randomized, controlled trial. Crit Care Med 1996, 24: 1360–1365. 10.1097/0000324619960800000015
 42.
Strohmenger HU, Lindner KH, Keller A, Lindner IM, Pfenninger EG: Spectral analysis of ventricular fibrillation and closedchest cardiopulmonary resuscitation. Resuscitation 1996, 33: 155–161. 10.1016/S03009572(96)010039
 43.
Strohmenger HU, Lindner KH, Lurie KG, Welz A, Georgieff M: Frequency of ventricular fibrillation as a predictor of defibrillation success during cardiac surgery. Anesth Analg 1994, 79: 434–438. 10.1213/0000053919940900000006
 44.
Strohmenger HU, Lindner KH, Prengel AW, Pfenninger EG, Bothner U, Lurie KG: Effects of epinephrine and vasopressin on median fibrillation frequency and defibrillation success in a porcine model of cardiopulmonary resuscitation. Resuscitation 1996, 31: 65–73. 10.1016/03009572(95)008993
 45.
Neurauter A, KramerJohansen J, Eilevstjonn J, Myklebust H, Wenzel V, Lindner KH, Eftestol T, Steen PA, Strohmenger HU: Estimation of the duration of ventricular fibrillation using ECG single feature analysis. Resuscitation 2007, 73: 246–252. 10.1016/j.resuscitation.2006.08.028
 46.
Werther T, Klotz A, Granegger M, Baubin M, Feichtinger HG, Amann A, Gilly H: Strong corruption of electrocardiograms caused by cardiopulmonary resuscitation reduces efficiency of twochannel methods for removing motion artefacts in nonshockable rhythms. Resuscitation 2009, 80: 1301–1307. 10.1016/j.resuscitation.2009.07.020
Acknowledgements
We thank our colleagues and the Red Cross staff at Innsbruck Emergency Medical Service for the recording of preclinical ECGdata. This work was supported by the Austrian Science Fund (FWF) under grant L288, and by the Oesterreichische Nationalbank (OeNB) under grant Jubilaeumsfondsprojekt No 8665.
Author information
Affiliations
Corresponding authors
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
The human study was planned and performed by MB and AA. The animal experiments were planned and performed by WL. The MATLAB mfile for coherent line removal was written by AK, with later modifications by AA. The compilation of human ECGdata (from different traces recorded on the MRLdefibrillator) was performed by TN, who also prepared a respective database and programmed the MATLAB toolbox for scoring and handling of the data. The toolbox was subsequently modified by AKu. The choice of the data and the respective MATLABexperiments were performed by AA with support by AK, TW and MG. The distributed computing of the coherent line algorithm was implemented by AKu. The manuscript was written by AA, TW, MB and WL. All authors approved the final version of the manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Amann, A., Klotz, A., Niederklapfer, T. et al. Reduction of CPR artifacts in the ventricular fibrillation ECG by coherent line removal. BioMed Eng OnLine 9, 2 (2010). https://doi.org/10.1186/1475925X92
Received:
Accepted:
Published:
Keywords
 Ventricular Fibrillation
 Dominant Frequency
 Median Frequency
 Frequency Window
 Adaptive Noise Cancellation