ECG-based gating in ultra high field cardiovascular magnetic resonance using an independent component analysis approach
Journal of Cardiovascular Magnetic Resonance volume 15, Article number: 104 (2013)
In Cardiovascular Magnetic Resonance (CMR), the synchronization of image acquisition with heart motion is performed in clinical practice by processing the electrocardiogram (ECG). The ECG-based synchronization is well established for MR scanners with magnetic fields up to 3 T. However, this technique is prone to errors in ultra high field environments, e.g. in 7 T MR scanners as used in research applications. The high magnetic fields cause severe magnetohydrodynamic (MHD) effects which disturb the ECG signal. Image synchronization is thus less reliable and yields artefacts in CMR images.
A strategy based on Independent Component Analysis (ICA) was pursued in this work to enhance the ECG contribution and attenuate the MHD effect. ICA was applied to 12-lead ECG signals recorded inside a 7 T MR scanner. An automatic source identification procedure was proposed to identify an independent component (IC) dominated by the ECG signal. The identified IC was then used for detecting the R-peaks. The presented ICA-based method was compared to other R-peak detection methods using 1) the raw ECG signal, 2) the raw vectorcardiogram (VCG), 3) the state-of-the-art gating technique based on the VCG, 4) an updated version of the VCG-based approach and 5) the ICA of the VCG.
ECG signals from eight volunteers were recorded inside the MR scanner. Recordings with an overall length of 87 min accounting for 5457 QRS complexes were available for the analysis. The records were divided into a training and a test dataset. In terms of R-peak detection within the test dataset, the proposed ICA-based algorithm achieved a detection performance with an average sensitivity (Se) of 99.2%, a positive predictive value (+P) of 99.1%, with an average trigger delay and jitter of 5.8 ms and 5.0 ms, respectively. Long term stability of the demixing matrix was shown based on two measurements of the same subject, each being separated by one year, whereas an averaged detection performance of Se = 99.4% and +P = 99.7% was achieved.
Compared to the state-of-the-art VCG-based gating technique at 7 T, the proposed method increased the sensitivity and positive predictive value within the test dataset by 27.1% and 42.7%, respectively.
The presented ICA-based method allows the estimation and identification of an IC dominated by the ECG signal. R-peak detection based on this IC outperforms the state-of-the-art VCG-based technique in a 7 T MR scanner environment.
Cardiac applications have a steadily gaining popularity in Magnetic Resonance Imaging (MRI). Nowadays, Cardiovascular Magnetic Resonance (CMR) is applied in clinical routine in MR scanners with field strengths up to 3 T. Ultra high field (UHF) (i.e. MRI scanners with 7 T or 9.4 T) CMR applications are becoming more and more popular in the research community[2, 3].
Challenges in CMR are the synchronization of the cardiac activity with the acquisition of the MR images and the compensation of motion artefacts during the cardiac cycle. The synchronization is called gating and can be performed prospectively and retrospectively. The state-of-the-art methods are based on the electrocardiogram (ECG) signal[5, 6]. These methods rely on the detection of the R-wave, which is a feature more or less easily extracted and occurring while the ventricle cells are depolarized. Other techniques for synchronisation are based on Doppler ultrasound, finger plethysmography, self gating techniques or a recently developed phonocardiogram technique.
Inside the MR scanner, ECG signals are severely affected by the magnetohydrodynamic (MHD) effect. The MHD effect results from the interaction between the pulsatile blood flow – which is caused by the rhythmic action of the heart – and the static magnetic field of the MR scanner, B0. The MHD effect leads to an induction of a voltage across blood vessels which are perpendicular to the MR scanner’s static magnetic field. The resulting body surface potentials of the MHD effect superimpose the ECG signal which exacerbates the detection of the R-peaks. It is assumed that the majority of the MHD-related distortions are caused by the flow of blood in the aortic arch[11, 12]. The MHD effect is highest during the ejection phase, i.e. during the ventricular systole, where the blood from the ventricles flows into the aorta and the pulmonary arteries (Figure1)a.
ECG-based gating methods rely on a spatial representation of the heart’s electrical activity: the vectorcardiogram (VCG). These methods assume that the R-wave can be spatially discriminated from the MRI induced artefacts and more specifically from the magnetohydrodynamic (MHD) effect[5, 6, 13].
However, the VCG-based gating approaches are prone to errors in MR scanners with magnetic field strengths of 7 T and above due to the increased MHD effect at these fields strengths[9, 14, 15]. In such cases the image acquisition is based on false detections, i.e. missed R-waves and the detection of false positives can occur. Consequently, the resulting MR image is contaminated with motion artefacts and blurring effects.
This work investigates the use of Independent Component Analysis (ICA) for the enhancement of the ECG contribution over the MHD effect in order to improve the R-peak detection and hence the gating ability of the ECG in UHF MRI scanners resulting in an increase of the intrinsic MR image quality. ICA is a statistical technique applied to multidimensional data in order to blindly separate the underlying signal components based on a cost function which maximizes the statistical independence of the resulting components[16, 17]. In ECG signal processing, ICA was used for many applications such as the extraction of the foetal ECG, for artefact removal or for the extraction of atrial activity[18–20]. In CMR, ICA was used to reduce the gradient artefacts from 3-lead ECG signals. In this work, ICA is applied to 12-lead ECG signals recorded inside a 7 T MR scanner. The independent component (IC) which is dominated by the ECG signal is automatically identified by a simple classification rule based on a template matching method. The identified IC is used for QRS detection.
ECG database acquisition
A standard 12-lead Holter ECG (CardioMem CM3000-12, GETEMED, Germany) with a sampling rate of 1024 Hz, a resolution of 12 bits, an input voltage range of ±6 mV and an analogue bandwidth ranging from 0.05 Hz to 100 Hz was used for recording the ECGs. The recorded ECGs consisted of the limb leads (I, II, III), of the augmented limb leads (aVR, aVL, aVF) and of the six precordial leads (V1-V6). Figure2(a) shows the electrode positions of a 12-lead ECG as it was used during the experiments. In contrast, Figure2(b) depicts an orthogonal 2-lead ECG configuration as it is typically used in CMR.
The ECG signals referred to in this work were recorded in a 7.0 T MR scanner (Magnetom, Siemens, Germany). MR imaging was switched off during all ECG measurements and hence gradient and radiofrequency magnetic fields did not exist. The ECG device was placed outside the MR scanner room using lead cable extensions. Standard MR-safe ECG Electrodes were used (H34SG, Covidien, Ireland). The study was approved by the local ethics committee, and written informed consent was obtained from all subjects prior to the ECG measurements. The ECG signals were recorded from eight healthy volunteers (seven male, one female) aged between 24 and 31 years and were recorded outside and inside the MR scanner during normal breathing in supine position. Inside the MR scanner, the ECG measurements were made in a Feet first (Ff) and in a Head first (Hf) position.
A change between the Ff and Hf position inverts the orientation of the MR scanner’s static magnetic field with respect to the volunteer. Hence, the polarity of the MHD effect changes[10, 22]. The thorax was positioned in the centre of the MR scanner’s bore during the measurements. The ECG signals were acquired for 2-3 min outside and 4-5 min for each position inside the scanner. Nine different ECG datasets D1- D9 were recorded inside the MR scanner. Each of the nine datasets consisted of two subsets: one subset containing the Ff (D1(Ff) - D9(Ff)) and one subset containing the Hf (D1(Hf) - D9(Hf)) measurements. For the development and evaluation of the ICA-based gating method, the ECG records were separated into a training and a test dataset. D1- D5 were used as training datasets and D6- D9 were used as test datasets. Hence, ten ECG records were available from the training dataset (D1(Ff)- D5(Ff) and D1(Hf)- D5(Hf)) and eight records were available from the test dataset (D6(Ff)- D9(Ff) and D6(Hf)- D9(Hf)). The ten subsets of D1- D5 (training dataset) had a total length of 47 min corresponding to 2853 R-peaks. The eight subsets of D6- D9 (test dataset) had a total length of 40 min corresponding to 2604 R-peaks. Test dataset D9 was recorded from the same volunteer as the training dataset D5, but D5 was recorded one year prior to D9.
Additional ECG datasets recorded outside the MR scanner from each volunteer were used for the generation of QRS templates, for comparing different signal properties and for the application of the VCG-based algorithms. Figure3(a)-3(f) show ECG signals of leads II and V3 measured outside and inside the 7 T MR scanner. Figure4 depicts the variations of the MHD effect of the precordial lead V3 in four different datasets recorded in the Hf position.
ICA-based suppression of the MHD effect
The ECG signals recorded inside the MRI scanner were contaminated by the MHD and other noise components. This linear mixture of the different signal components can be described by
where A is the so-called (linear) mixing matrix, s k is a vector containing the source signals and x k is the measurement signal vector containing the mixture of the source signals at time instant k. In order to separate the ECG from the MHD and other noise components, an estimate of the real source signals,, was required. Therefore, ICA was applied to the measured data vector x k to find a demixing matrix W so that:
where are the estimated source signals or independent components (ICs). Several algorithms have been proposed to solve this problem which differ mainly from the cost function used for measuring the Gaussianity (which is a measure of independence) of the signals. In this work the FastICA algorithm was used. The FastICA algorithm was not adapted to the specific problem discussed in this work. Where required, the dimensionality of x k was reduced by limiting the number of eigenvalues or principal components (PCs) prior to the application of ICA.
From the estimated ICs, one IC which was dominated by the R-peak or QRS complex needed to be identified. This pivotal step is described in the following section.
Identification of an independent component for gating
One crucial step with the application of ICA to the given problem was the identification of the IC. This IC was required to represent the cardiac activity, i.e. it should contain the R-peak so that standard R-peak/QRS detection algorithms could be applied. For the practical application of the proposed method, an automated identification of was required.
A template matching algorithm depicted in Figure5 was used for an automated identification of. The demixing matrix W which was estimated from the disturbed ECG signal x k was applied to the ECG xk,out acquired outside the MR scanner resulting in a set of linear combinations. For each of these linear combinations, a template uk,qrs of the QRS complex was generated using the following procedure. First, a QRS detection algorithm was applied to the clean ECG signal in order to estimate the positions of the QRS complexes. The QRS template for each linear lead combination was then estimated from segments of 80 ms length centred around the detected QRS positions and averaged over ten consecutive QRS complexes. For the identification of, the cross-correlation functions between each QRS template in uk,qrs and the corresponding IC in were calculated. Based on the maximum value of the cross-correlation functions, one IC in was identified as. The IC corresponding to in is referred to as.
Description of the proposed method
The extraction of from the initial ECG signals and the subsequent R-peak detection were applied in a two-stage process (Figure6).
In the first stage, ICA was applied to a 30 s reference segment xRef,k of the 12-lead ECG signal which was recorded inside the MR scanner. The resulting demixing matrix W was used to transform the reference vector xRef,k into the vector containing the estimated sources, i.e.. This is described by Eq. 2. The estimated signals contained in were dominated by different (physiological) signal sources, e.g. by the ECG signal itself, by the MHD effect, by respiratory motion artefacts or by measurement noise. After the identification of, the demixing matrix W was reduced to a demixing vector w since only was required for the R-peak detection. The estimation of can be summarized as the scalar product of w and x k :
where w is a row vector and x k is a column vector with each row in x k corresponding to one of the twelve ECG leads at time instant k. The estimated demixing vector w remained unchanged, i.e. it was only estimated once for each dataset.
In the second stage, the demixing vector w was applied to the most recent ECG sample x k resulting in a real-time estimation of. The demixing process was followed by the application of a R-peak/QRS detection algorithm to.
This R-peak detection algorithm uses the first derivative of the ECG signal to enhance the QRS complex. To ensure a reliable R-peak detection inside the MR scanner, the slope of the QRS complex has to be preserved in. This was investigated quantitatively by comparing the upslope differences of the QRS complex in the identified IC and in the corresponding IC. After testing both distributions (upslopes in and in) for normality using a Lilliefors test, a two-sample t-test (α = 0.05) was conducted to test if the mean slopes of and showed significant differences.
ECG lead configurations used for estimating
For the estimation of, ICA was applied to different combinations of the recorded ECG leads. The aim of this experiment was to investigate how the additional precordial ECG leads affected the quality of the R-peak detection results.
At the first step, ICA was applied to 8 of the 12 leads contained in each ECG dataset. These were the limb leads I and II and the precordial leads V1-V6. The remaining 4 leads (III, aVR, aVL and aVF) were excluded because they can be derived from linear combinations of limb leads I and II. Hence, these 4 leads contained no additional information. The usage of 8 leads for the estimation of is referred to as lead configuration LC1.
In the second step, ICA was applied to the same datasets but fewer ECG leads were used. These datasets consisted of combinations of the limb leads I and II and of different precordial leads. Furthermore, the influence of decreasing the dimensionality of the datasets was investigated. This was achieved by a reduction of the principal components (PCs) before ICA was applied to the datasets. These combinations of varying ECG leads and PCs are referred to as lead configurations LC2a-LC2h. Results for methods LC2 using leads I, II and at least two precordial lead are given for the best lead/PC combinations. Results for different combinations of lead I and II and of one precordial lead are computed to highlight the influence of a single precordial lead on the R-peak detection (LC2e-LC2g). Finally, ICA was applied to lead I and II without any additional precordial lead (LC2h).
Other R-peak detection methods (M1-M5)
Besides the ICA-based method for QRS detection, other methods for detecting the R-peaks in the contaminated datasets were investigated. The five other R-peak detection techniques, called methods M1-M5, were applied to the 12-lead ECG signals which were acquired using the electrode configuration shown in Figure2(a). The 3D VCG v(k) = (v x (k),v y (k),v z (k)) required for methods M2-M5 was obtained from the 12-lead ECG signals using the inverse Dower matrix. An overview of the five different methods M1-M5 which were implemented for comparison with the proposed ICA-based method is given below. Methods M1, M2 and M5 used the R-peak detector from. R-peak detection results using methods M1-M5 are given for the training and test dataset.
R-peak detection in a single ECG lead (M1)
The QRS detection algorithm was successively applied to each of the 12 leads in all datasets. Finally, one lead was chosen which gave the best results averaged over all datasets. Results are given for this single lead.
R-peak detection in a single VCG lead (M2)
QRS detection was performed in all three VCG leads v(k) = (v x (k),v y (k),v z (k)). Results are presented for the same lead in all datasets which gave the best average detection results.
3D VCG-based R-peak detection (M3)
The VCG-based gating method was implemented according to[5, 13]. A reference vector rout = (r x ,r y ,r z ) which points at the spatial position of the R-peak was obtained for each dataset based on the clean ECG signal acquired outside the MR scanner. The distance and angle between rout and the VCG vector v(k) was calculated. Those positions of the vector v(k) where the distance and angle between rout and v(k) fell below a threshold were defined as R-peaks.
Modified 3D VCG-based R-peak detection (M4)
As explained in the introduction, the R-peak detection based on the 3D VCG (method M3) is prone to errors in 7 T MR scanners due to the altered amplitude and angle of the R-peak vector in the actual VCG v(k). To cope with this problem, a modified version of the VCG-based method was implemented. This method calculated a reference vector rin from the ECG signal which was recorded inside the MR scanner instead of obtaining it from the ECG recorded outside the scanner as described by method M3. The remaining procedure of the VCG-based algorithm remained unchanged.
ICA of the VCG for R-peak detection (M5)
To further analyse the usability of the VCG signal v(k), ICA was applied to the three VCG leads v x (k),v y (k) and v z (k) of each dataset. As for the ICA-based method, one IC was identified using the template matching process.
Evaluation of the different gating methods
R-peak detection performance
The QRS detection performance using the different methods was evaluated by the sensitivity (Se) and the positive predictive value (+P):
where TP is the number of true positives, FP the number false positives and FN the number of false negatives. These parameters were estimated according to the ANSI/AAMI EC57 standard recommendations. A manual annotation of the R-peaks was used as gold standard for performance estimation of the different R-peak detection approaches. The manual annotations were made by two ECG experts.
Mean trigger propagation delay μpd and jitter σpd
The stability of the trigger timing was evaluated based on the propagation delay and the jitter. The mean trigger propagation delay μpd was defined as the difference between the position of the manually annotated R-peak and the R-peak detected by the QRS detection algorithm and was averaged over each dataset.
Long term stability and robustness of the proposed ICA-based gating method
In order to calculate the demixing matrix W only once per volunteer at the beginning of the measurement, the long term stability of W had to be investigated. The training set D5 and the test set D9 – which were acquired from the same volunteer but were separated by one year – were used for this experiment. Two demixing matrices and obtained from the datasets D5(Ff) and D5(Hf) were applied to the two corresponding datasets D9(Ff) and D9(Hf) resulting in two ICs () which were used for R-peak detection. R-peak detection results are given in terms of Se and +P for both datasets.
Estimation of, identification of and slope of the QRS complex in and
Four exemplary ICs from two different datasets obtained by FastICA are shown in Figure7: two ICs mainly dominated by the ECG, one dominated by the ECG and the MHD effect and one IC mainly dominated by the MHD effect. The IC dominated by the MHD effect is characterized by an attenuation of the R-peak’s amplitude. None of the ICs could be clearly assigned to a certain signal source since each IC contained a mixture of ECG, MHD and noise components.
The IC required for R-peak detection was identified in all datasets based on the template matching process described in the methods section.
Figure8(a)-(b) depict the IC used for gating and its first derivative. In Figure8(c), the IC obtained from the ECG acquired outside the MR scanner is shown. The first derivative or upslope of is depicted in Figure8(d). Using the Lillifors test it was shown that the upslopes of the QRS complex measured in and exhibit a normal distribution. Figure9 compares the mean and standard deviations of the QRS-upslopes in to those of the corresponding ICs. The average absolute relative deviations of the slopes in compared to were 14.1% in the training set and 16.7% in the test set. The two-sample t-test revealed a significant change of the slope’s mean value.
R-peak detection in using lead configurations LC1 and LC2
Detailed results for the R-peak detection in using lead configuration LC1 are given in Table1 for the training dataset and in Table2 for the test dataset. For the training dataset, an average sensitivity Se of 99.5% and an average positive predictive value +P of 99.4% were achieved. The average values of the mean propagation delay μpd and jitter σpd were 3.5 ms and 6.5 ms, respectively. Dataset D4(Ff) shows a higher jitter when compared to the other datasets. Comparable R-peak detection results were achieved for the test dataset with Se = 99.2%, +P = 99.1%, μpd = 5.8 ms and σpd = 5.0 ms. Exemplary sections of R-peaks in are shown in Figures7(a),7(c) and8(a).
Table3 summarizes the R-peak detection results in using lead configurations LC2a-LC2h. Se and +P slightly decreased when less precordial leads were used. Compared to lead configuration LC1, the mean propagation delay μpd and jitter σpd increased. The results for the 3-lead configurations LC2e-LC2g revealed the impact of the precordial leads V2-V4 on the overall R-peak detection quality. Out of these 3-lead configurations, LC2g using the precordial lead V4 gave the best R-peak detection results (Se = 98.2%, +P = 97.3% for the training dataset and Se = 98.5%, +P = 97.2% for the test dataset). Without using any of the precordial leads but only limb leads I and II as in configuration LC2h, Se and +P decreased to 77.1% and 66.3% for the training dataset and to 78.9% and 65.8% for the test dataset.
Long term stability of the demixing matrix
Two different ICs were obtained after applying the demixing matrices and to the corresponding ECG datasets D9(Ff) and D9(Hf) acquired from the same volunteer. Using these two ICs, the following R-peak detection results were achieved:
D9(Ff): Se = 99.7%, +P = 99.6%
D9(Hf): Se = 99.1%, +P = 99.7%
Other R-peak detection methods
Table4 compares the results obtained from the gating methods M1-M5 to the results of the ICA-based approach using lead configuration LC1.
The results given for the ECG-based method (M1) were estimated from the precordial lead V4 in all datasets because this ECG lead yielded the best R-peak detection in terms of the average Se and +P. Whereas the MHD induced contaminations in the limb leads resulted in a lower R-peak detection quality, the precordial leads V3-V5 showed less contaminations and hence better R-peak detection results. Figure3 depicts this difference between the limb lead II and the precordial lead V3.
The usage of the VCG lead x for R-peak detection in method M2 lead to an increased Se and +P when compared the detection results in method M1.
The state-of-the-art VCG-based gating method (M3) resulted in a Se and +P within the test dataset of 72.1% and 56.4%, respectively. The usage of a modified reference vector rin in method M4 lead to an increase of +P of approximately 21%. A comparison of both methods for one dataset is shown in Figure10. Figure10(a) depicts the results from the unmodified VCG-based approach (M3). This approach resulted in multiple peaks per R-peak which lead to the low +P given in Table4. In contrast, the results from the VCG-based approach depicted in Figure10(b) which used a modified reference vector rin (M4) showed less peaks per R-peak resulting in a higher +P.
The application of ICA to the three VCG leads in method M5 resulted in a decreased Se and +P when compared to method M2 where QRS detection was directly performed in one of the VCG leads.
ICA-based gating method
Although it was shown previously that ICA is not suitable for recovering the diagnostic information contained in the MHD-contaminated ECG signal, ICA was successfully applied to 12-lead ECG signals for the real-time estimation of an IC which was dominated by the R-peak. This IC was then used for R-peak detection. The influence of different ECG leads on the results achieved by ICA was investigated by applying ICA to different lead combinations. The best R-peak detection results were achieved with lead configuration LC1 where eight ECG leads (I,II,V1-V6) were used for the estimation of. For the combinations of only three ECG leads used with ICA, it was shown that the best average R-peak detection results were achieved by a combination of the two limb leads I and II and of the precordial lead V4. These experiments emphasized the importance of the precordial leads – especially of lead V4 – for the proposed method. Not using the precordial leads as in configuration LC2h caused a substantial decrease in R-peak detection performance (Table3).
For 7 T CMR application based on a prospective gating scheme, the mean trigger propagation delay μpd should be less than 20 ms because the systole’s mechanical contraction begins 30 ms-70 ms after the R-peak. A large jitter σpd could lead to blurring of the CMR image, to ghost artefacts in the phase encoding direction or to a false estimation of the blood flow rate in phase contrast MRI. For the estimation of blood flow rates, the jitter should be less than 15 ms. Using the proposed ICA-based method with lead configuration LC1, a mean detection delay μpd of 5.8 ms and a jitter σpd of 5.0 ms were achieved in the test dataset using eight ECG leads in lead configuration LC1. Hence, the usage of this method would enable the acquisition of proper CMR images.
The demixing matrix W was only estimated once for each dataset using a 30 s segment of the recorded ECG signal and was reduced to a demixing vector w after was identified. This allowed a demixing of the different signal components and the estimation of in real-time. W was stable over a long period of time (1 year) which was shown using demixing matrices from previously recorded 7 T datasets of the same volunteer – even though the ECG electrode positions slightly varied between both measurements. Besides, the demixing matrix W was obtained from and applied to datasets which were acquired during normal breathing. Large baseline shifts caused by the physical movement of the ECG electrodes during breathing can affect the ECG signals. These variations can complicate the R-peak detection. The results obtained in this study show that the proposed method for the estimation of can be applied over a long period of time and is robust against variations caused by normal breathing.
The demixing matrices W varied between the different datasets D1-D9. It was however noticed that the precordial leads were contributing more for the generation of the IC. Nevertheless, the combination of the precordial leads varied over the subject population, and it was also shown in this study that the limb leads add information and allow for a better separation of the MHD and the ECG. Two reasons can be given for the variation of the demixing matrices. Firstly, the anatomy and physiology of the different volunteers leads to high variations of the MHD signals as it is shown in Figure4 and was previously discussed in. Hence, the ECG leads were not affected in the same way in the datasets D1-D9 and different combinations of the measured ECG leads x k were required to obtain an IC which was suitable for R-peak detection. Secondly, this phenomenon can be explained by the positions of the ECG electrodes with respect to the position and orientation of the heart’s electrical axis. These anatomical parameters vary between the different volunteers or datasets. As a result, the ECG traces – either recorded outside or inside the MR scanner – vary between the different volunteers. This requires one initial computation of the demixing matrix W for each new subject. In this work it was shown that ICA provides an elegant way of estimating W for different subjects without further assumptions.
A similar ICA-based method was utilized previously for the suppression of gradient artefacts from 3-lead ECG signals acquired during MRI. These gradient induced artefacts exhibit temporal variations due to different MRI sequences, slice locations and orientations. Hence, a continuous update of the demixing matrix was required during runtime. As described in this work, an update of the demixing matrix for the suppression of the MHD effect using the 12-lead ECG signals was not necessary.
Given the results discussed above, the presented ICA-based method is suitable to be used in 7 T CMR imaging for the real-time estimation of an IC which can then be used for R-peak detection and gating.
Comparison to other gating methods
From the other R-peak detection techniques M1-M5, method M2 using VCG lead x gave the best results. For method M1 which utilized the original ECG signal, the best R-peak detection results were achieved using the precordial ECG lead V4. Together with the results obtained by the ICA-based method, this again emphasizes the importance of the precordial leads for R-peak detection in 7 T CMR.
The results obtained for the state-of-the-art VCG-based gating method (M3) proposed in[5, 13] showed low Se and +P values which are insufficient for gating at 7 T CMR. Based on the functional principle of the VCG method, a modified reference vector was suggested which increased +P within the test dataset by approximately 21%. Although this modification improved the VCG-based technique, results are still not sufficient for a proper gating in CMR. Furthermore, the required manual R-peak annotation makes this modification impractical for a clinical application.
ECG recording hardware
The results achieved in this work are based on 12-lead ECG signals. However, the hardware which was used to record the 12-lead ECGs is not MR safe and hence it cannot be used during MR imaging. Up to now, there is no MR safe 12-lead ECG system available on the market but a commercial 12-lead ECG system was made MR conditional previously[32, 33]. With the results obtained in this work showing the superior signal quality of the precordial leads, such hardware could be promising for future applications – whether for gating in UHF CMR or for ECG-based diagnostics during MRI.
The demixing matrix W was only estimated once for each dataset. Where this approach was successfully applied to nine datasets D1-D9 of healthy volunteers, problems could arise for patients with cardiac pathologies. For such cases, a recalculation of W could be required during runtime. This needs to be evaluated with additional datasets of volunteers or patients suffering from myocardial pathologies.
Due to the lack of MR-conditional 12-lead ECG hardware devices, contaminations caused by the switched gradient magnetic fields of the MR scanner have not been investigated. In this case, however, the demixing matrix W would be computed before playing the gradients. The resulting IC is a linear combination of the original ECG leads. Based on these linear combinations, the gradient artefacts in will still respect the linear time invariant assumption needed for the application of previously published methods for suppressing the gradient artefacts[21, 34, 35].
Due to the MHD effect, the state-of-the-art VCG-based gating technique is prone to errors in 7 T CMR. To cope with this problem, an approach based on ICA was developed and applied to 12-lead ECG signals recorded inside a 7 T MR scanner. An IC which was dominated by the R-peak was estimated by ICA. This IC was automatically identified and used for R-peak detection. Given the high R-peak detection quality within the test dataset (Se = 99.2%, +P = 99.1%), the proposed ICA-based technique outperforms the state-of-the-art VCG-based gating method (Se = 72.1%, +P = 56.4%).
Previous investigations revealed that the precordial leads of the ECG are less affected by the distortions caused by the MHD effect[22, 36]. The importance of the precordial leads for the proposed method was demonstrated by applying ICA to several lead configurations (LC2) with different combinations of the limb and precordial leads and by using a single ECG lead for R-peak detection (method M1). This finding could be helpful for the development of improved MR-conditional ECG hardware with additional precordial leads or for further investigations regarding the placement of the ECG electrodes.
a Image source (Figure1):https://upload.wikimedia.org/wikipedia/commons/f/f4/Wiggers_Diagram.svg
Cardiovascular magnetic resonance
Independent component analysis
Magnetic resonance imaging
Ultra high field
Attili AK, Schuster A, Nagel E, Reiber JH, van der Geest RJtitle=Quantification in cardiac MRI: advances in image acquisition and processing: Int J Cardiovasc Imaging. 2010, 26: 27-40.
Niendorf T, Sodickson DK, Krombach GA, Schulz-Menger J: Toward cardiovascular MRI at 7 T: clinical needs, technical solutions and research promises. Eur Radiol. 2010, 20 (12): 2806-16. 10.1007/s00330-010-1902-8.
van Ooij P, Kleinloog R, Zwanenburg J, Visser F, Luijten P, Barker A, Markl M, Nederveen A, Majoie C, Regli L, et al: Improved depiction of hemodynamics in intracranial aneurysms by 4D flow MRI at 7T compared to 3T. J Cardiovasc Magn Reson. 2013, 15 (Suppl 1): W12-10.1186/1532-429X-15-S1-W12.
Scott A, Keegan J, Firmin D: Motion in cardiovascular MR imaging. Radiology. 2009, 250 (2): 331-51. 10.1148/radiol.2502071998.
Fischer S, Wickline S, Lorenz C: Novel real-time R-wave detection algorithm based on the vectorcardiogram for accurate gated magnetic resonance acquisitions. Magnet Reson Med. 1999, 42 (2): 361-70. 10.1002/(SICI)1522-2594(199908)42:2<361::AID-MRM18>3.0.CO;2-9.
Nacif M, Zavodni A, Kawel N, Choi E, Lima J, Bluemke D: Cardiac magnetic resonance imaging and its electrocardiographs (ECG): tips and tricks. Int J Cardiovasc Imaging. 2011, 28 (6): 1-11.
Rubin J, Brian Fowlkes J, Prince M, Rhee R, Chenevert T: Doppler US gating of cardiac MR imaging. Acad Radiol. 2000, 7 (12): 1116-22. 10.1016/S1076-6332(00)80065-3.
Nijm G, Sahakian A, Swiryn S, Carr J, Sheehan J, Larson A: Comparison of self-gated cine MRI retrospective cardiac synchronization algorithms. J Magn Reson Imaging. 2008, 28 (3): 767-72. 10.1002/jmri.21514.
Frauenrath T, Hezel F, Renz W, de Geyer dÓrth T, Dieringer M, von Knobelsdorff-Brenkenhoff F, Prothmann M, Schulz-Menger J, Niendorf T: Acoustic cardiac triggering: a practical solution for synchronization and gating of cardiovascular magnetic resonance at 7 Tesla. J Cardiovasc Magn Reson. 2010, 12 (67): 1-14.
Keltner JR, Roos MS, Brakeman PR, Budinger TF: Magnetohydrodynamics of Blood Flow. Magnet Reson Med. 1990, 16: 139-49. 10.1002/mrm.1910160113.
Kinouchi Y, Yamaguchi H, Tenforde TS: Theoretical analysis of magnetic field interactions with aortic blood flow. Bioelectromagnetics. 1996, 17: 21-32. 10.1002/(SICI)1521-186X(1996)17:1<21::AID-BEM3>3.0.CO;2-8. [http://dx.doi.org/3.0.CO;2-8],
Martin V, Drochon A, Fokapu O, Gerbeau J: MagnetoHemoDynamics Effect on Electrocardiograms. FIMH, Lect Notes Comput Sc. 2011, 6666: 325-32.
Chia J, Fischer S, Wickline S, Lorenz C: Performance of QRS detection for cardiac magnetic resonance imaging with a novel vectorcardiographic triggering method. J Magn Reson Imaging. 2000, 12 (5): 678-88. 10.1002/1522-2586(200011)12:5<678::AID-JMRI4>3.0.CO;2-5.
Brandts A, Westenberg J, Versluis M, Kroft L, Smith N, Webb A, de Roos A: Quantitative assessment of left ventricular function in humans at 7 T. Magnet Reson Med. 2010, 64 (5): 1471-7. 10.1002/mrm.22529.
Krug J, Rose G, Stucht D, Clifford G, Oster J: Limitations of VCG based gating methods in ultra high field cardiac MRI. J Cardiovasc Magn Reson. 2013, 15 (Suppl 1): W19-10.1186/1532-429X-15-S1-W19.
Hyvarinen A: Survey on independent component analysis. Neural Comput Surv. 1999, 2 (4): 94-128.
Comon P, Jutten C: Handbook of Blind Source Separation: Independent component analysis and applications. 2010, New York: Academic Press
Sameni R, Jutten C, Shamsollahi MB: What ICA provides for ECG processing: Application to noninvasive fetal ECG extraction. Proc Int Symp Signal Process Inf Technol (ISSPIT). 2006, Canada: Vancouver, 656-661.
He T, Clifford G, Tarassenko L: Application of ICA in Removing Artefacts from the ECG. Neural Comput Appl. 2006, 15 (2): 105-16. 10.1007/s00521-005-0013-y.
Phlypo R, Zarzoso V, Comon P, Dasseler Y, Lemahieu I: Extraction of atrial activity from the ECG by spectrally constrained ICA based on kurtosis sign. Proc 7th Int Conf ICA, LNCS. 2007, 4666: 641-8.
Oster J, Pietquin O, Abächerli R, Kraemer M, Felblinger J: Independent component analysis-based artefact reduction: application to the electrocardiogram for improved magnetic resonance imaging triggering. Physiol Meas. 2009, 30: 1381-97. 10.1088/0967-3334/30/12/007.
Krug J, Rose G: Comput Cardiol. 2012, 38: 769-72.
Hyvarinen A: Fast and robust fixed-point algorithms for independent component analysis. IEEE Trans Neural Netw. 1999, 10 (3): 626-34. 10.1109/72.761722.
Clifford GD: Signal processing methods for heart rate variability [Phd thesis]. England: University of Oxford; 2002
Lilliefors HW: On the Kolmogorov-Smirnov test for normality with mean and variance unknown. J Am Stat Assoc. 1967, 62 (318): 399-402. 10.1080/01621459.1967.10482916.
Edenbrandt L, Pahlm O: Vectorcardiogram synthesized from a 12-lead ECG: superiority of the inverse Dower matrix. J Electrocardiol. 1988, 21 (4): 361-367. 10.1016/0022-0736(88)90113-6.
American National Standard: Testing and reporting performance results of cardiac rhythm and ST segment measurement algorithms. ANSI/AAMI/ISO EC57:1998/(R)2008. Arlington (VA), 1998
Krug J, Rose G, Stucht D, Clifford G, Oster J: Filtering the magnetohydrodynamic effect from 12-lead ECG signals using independent component analysis. Comput Cardiol. 2012, 39: 590-93.
Lentner C: Geigy scientific tables: vol 5. Heart and circulation. 1990, West Caldwell (NJ): CIBA-Geigy AG
Thompson RB, McVeigh ER: Flow-gated phase-contrast MRI using radial acquisitions. Magnet Reson Med. 2004, 52 (3): 598-604. 10.1002/mrm.20187.
Kyriakou A, Neufeld E, Szczerba D, Kainz W, Luechinger R, Kozerke S, McGregor R, Kuster N: Patient-specific simulations and measurements of the magneto-hemodynamic effect in human primary vessels. Physiol Meas. 2012, 33: 117-30. 10.1088/0967-3334/33/2/117.
Tse Z, Dumoulin C, Clifford G, Jerosch-Herald M, Kacher D, Kwong R, Stevenson W, Schmidt E: MRI-compatible 12-lead ECGs with MHD separation: application to cardiac MRI gating, physiological monitoring and noninvasive cardiac-output estimation. 2010, Stockholm: Sweden
Tse Z, Dumoulin C, Clifford G, Schweitzer J, Qin L, Oster J, Jerosch-Herold M, Kwong R, Michaud G, Stevenson W, Schmidt E: A 1.5T MRI-conditional 12-lead electrocardiogram for MRI and intra-MR intervention. Magnet Reson Med. 2013, Epub ahead of print.http://dx.doi.org/10.1002/mrm.24744,
Felblinger J, Slotboom J, Kreis R, Jung B, Boesch C, et al: Restoration of electrophysiological signals distorted by inductive effects of magnetic field gradients during MR sequences. Magnet Reson Med. 1999, 41 (4): 715-21. 10.1002/(SICI)1522-2594(199904)41:4<715::AID-MRM9>3.0.CO;2-7.
Oster J, Pietquin O, Kraemer M, Felblinger J: Nonlinear Bayesian filtering for denoising of electrocardiograms acquired in a magnetic resonance environment. IEEE Trans Biomed Eng. 2010, 57 (7): 1628-38.
Schmidt EJ, Clifford G, Jerosch-Herald M, Kwong RY, Epstein L, Kacher D, Dumoulin CL, Jolesz F: MRI-compatible 12-lead ECG: improved MHD suppression, ischemia monitoring, and non-invasive cardiac output. 2009, Hawaii: Honolulu
JWK is funded by the Federal Ministry of Education and Research (BMBF, Germany) in context of the 'INKA’ project (03IP710). JO is a Newton international fellow (round 2010 from the Royal Academy of Engineering). The ECG recording hardware was kindly provided by GETEMED AG, Germany. The authors would like to thank Daniel Stucht for supporting the ECG data acquisition.
The authors declare that they have no competing interests.
The initial manuscript was read critically and edited by all authors. JWK, GDC and JO developed and evaluated the algorithms. JWK designed and conducted the measurements and implemented the algorithms. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Krug, J.W., Rose, G., Clifford, G.D. et al. ECG-based gating in ultra high field cardiovascular magnetic resonance using an independent component analysis approach. J Cardiovasc Magn Reson 15, 104 (2013). https://doi.org/10.1186/1532-429X-15-104