Demonstration of value of optimizing ECG triggering for cardiovascular magnetic resonance in patients with congenital heart disease
© Knesewitsch et al.; licensee BioMed Central Ltd. 2013
Received: 30 July 2012
Accepted: 6 December 2012
Published: 16 January 2013
Optimal ECG triggering is of paramount importance for correct blood flow quantification during cardiovascular magnetic resonance (CMR). However, optimal ECG triggering and therefore blood flow quantification is impaired in many patients with congenital heart disease (CHD) due to complex QRS patterns. Therefore, a new ECG-trigger algorithm was developed to address triggering problems due to complex QRS patterns.
The aim of this study was to test this new ECG-trigger algorithm in routine patients with CHD and its impact on blood flow quantification.
35 consecutive routine patients with CHD undergoing CMR were included in the study. (40% Fallot’s Tetralogy, 20% aortic arch pathology, 14% transposition of the great arteries, 26% others; age 26+/−11 yrs).
In all patients, blood flow in the ascending aorta was quantified using the old ECG-trigger algorithm and the new ECG-trigger algorithm in random order. Blood flow quantified using the old or new ECG-trigger algorithm was compared by Bland-Altman analysis.
Three blinded investigators evaluated the vector clouds and trigger points of both ECG-trigger methods. Evaluation criteria were false positive and false negative triggered QRS complexes (specificity and sensitivity), and accuracy of detection. Accuracy of detection was defined as time scatter of the trigger around the correct trigger point.
Specificity, sensitivity, and accuracy of detection significantly increased using the new ECG-trigger algorithm compared to the old ECG-trigger algorithm.
Blood flow quantification using the old or new ECG-trigger algorithm differed more than 5% in 31% of the cases.
Our results suggest that optimizing ECG triggering during CMR using our new algorithm can avoid errors of >5% in approximately 1/3 of routine patients with congenital heart disease (CHD). We furthermore suggest that incorrect ECG triggering appears to be problematic for blood flow quantification of many patients with CHD undergoing routine CMR.
KeywordsECG triggering Congenital heart disease Blood flow quantification Vectorcardiography (VCG)
Cardiovascular magnetic resonance (CMR) has evolved from a pure imaging method to a powerful tool in the diagnostic management of congenital heart disease (CHD) . This evolution is in a large part due to the capability of CMR to measure flow by phase-velocity CMR . To measure flow by phase-velocity CMR, synchronisation of the images to the heart phase is needed. This synchronisation is realised by ECG triggering. ECG triggering is based on real-time R-wave detection and ensures that each portion of an image is allocated to a specific phase of the cardiac cycle. Accordingly, an imprecise detection leads to an incorrect allocation of the images to the cardiac cycle and ultimately to false diagnostic results [2, 3]. Therefore, optimal ECG triggering is of paramount importance for correct blood flow quantification during CMR.
Reasons for imprecise triggering are known and efficient algorithms have been developed and used for several years. An important reason beside the general noise in normal ECG is the MR-specific environment, especially the magnetohydrodynamic effect. This effect leads to a deformation in the ST segment and may exhibit a T-wave with a larger amplitude than the QRS complex and finally to a wrong triggering.
In normal ECGs these problems seem to have been solved satisfactorily .
However, optimal ECG triggering and therefore blood flow quantification are impaired in many patients with CHD due to complex QRS patterns.
For example patients with CHD such as Tetralogy of Fallot or Ebstein anomaly partially show complex ECGs [4, 5]. Accordingly, in our routine clinical experience incorrect ECG triggering leading to unsatisfactory flow measurements using the old ECG-trigger algorithm occurs frequently. However, the exact extent of incorrect flow measurements is unknown. In a previous study 2 of 347 assessments even failed totally due to unreliable triggering .
Therefore, a new ECG-trigger algorithm was developed to address triggering problems and achieve better performance in patients with complex ECGs.
The new trigger is mainly based on a matched filter . The matched filter is able to detect the R-waves on the rising edge of the R-wave after an initial learning phase.
The old trigger algorithm works with threshold values, which are derived from the ECG.
The aim of this study was to test this new ECG-trigger algorithm in routine clinical patients with CHD and to measure its impact on blood flow quantification.
35 consecutive routine patients with CHD undergoing CMR were included in the study. Fourteen patients had Fallot’s Tetralogy, seven patients aortic arch pathology, five patients atrial switch operation for transposition of the great arteries, and nine patients had other complex congenital heart disease.
The patients consisted of 11 females and 24 males. The average age was 26 years +/−11 years ranging from 12 to 60 years.
Phase-velocity CMR measurements
For cardiovascular magnetic resonance (CMR) a standard cardiac 1.5 Tesla CMR scanner and a standard cardiac 12-channel coil (MAGNETOM Avanto®, software version VB15, Siemens Healthcare, Erlangen, Germany) was used.
To compare the difference of both triggering methods in all patients, the blood flow in the ascending aorta was quantified using the old ECG-trigger algorithm and the new ECG-trigger algorithm in random order directly after each other. The slice orientation and all other acquisition parameters were not changed between the two phase-velocity CMR measurements: free breathing (typical acquisition times: around 3 minutes), retrospective ECG gating, the velocity encoding was set to 200–550 cm/s according to the assumed jet velocity, the slice thickness was 5 mm, repetition time 36.7 ms, echo time 3.09 ms, flip angle 30°, averages 3, segmentation 3, number of phase-encoding steps 192, receiver bandwidth 31.25 kHz, rectangular field of view 260 to 330 × 330 mm, matrix 256 × 256, phase partial Fourier off. Data was reconstructed to provide 30 magnitude (anatomic) and phase (velocity-mapped) images per cardiac cycle. All measurements were automatically compensated for the concomitant gradient effects. It is important to note, that we always imaged the vessel’s region of interest at isocenter to maximize gradient fidelity. Additionally, we took meticulous care of placing the region of interest into the center of the imaging plane. The region of interest was placed into the center of the imaging plane because placing the region of interest out of the center of the imaging plane can lead to large errors .
Furthermore, the ECG was monitored during acquisition. The running acquisition was always aborted when extra systoles were noted.
The phase contrast images were processed with the post-processing software ARGUS (syngo MultiModality Workplace, version VE23B, Siemens Healthcare, Erlangen, Germany) by using manual contour segmentation in all phases of the cardiac cycle.
old ECG-trigger algorithm
The old algorithm is based on two ECG channels acquired simultaneously. The signals are supplied in a first processing branch to a low-pass filter and a derived value sum generator. The output signal compared with a threshold value generates a first comparison result. Additionally, the signals are fed in a second processing branch to a derived value generator. The output signal compared with an upper and lower threshold generates a second comparison result. Those first and second branches exist for both ECG channels. All first and second comparison results are evaluated in a weighted logic circuit, whose output are triggers for the MRI measurement . The mentioned upper and lower threshold values are learned during the initial learning phase.
new ECG-trigger algorithm
The new algorithm learns the shape of the rising edge of the R-wave during the initial learning phase in both ECG channels. Then, the trigger algorithm continuously compares the incoming ECG signal with the learned shape by realtime processing [6, 9]. The comparison is based on the following means: A matched filter which is widely used in telecommunications  and a filter that sums up the squared differences between the incoming ECG signal and the learned shape. The latter corresponds to a pattern-matching mechanism. Additionally, the angle of the VCG vector which is spanned by the two ECG channels at each time instant is utilized for detection of the R-wave. This angle criterion is used as a necessary condition but not as a sufficient condition for the generation of triggers.
The aim of the analysis was to detect failed triggering of the QRS complex and to evaluate the accuracy of the triggering. Both failures will probably lead first to a higher jitter of the RR interval and second to a false flow calculation.
For the evaluation of the triggering the ECG was recorded from the beginning to the end of each sequence, totalling in 192 QRS complexes, respectively.
Failed triggeringA two-channel ECG with the trigger marking was computed to analyse manually that exactly one trigger point is dedicated to every QRS complex. A heart cycle with two trigger points was counted as a false positive and a not triggered QRS-complex as a false negative (see Figures 1, 2 and 3). The results were used to calculate Sensitivity, Specificity and Performance Index (PI) of the method according to the following equations:
Accuracy of triggeringA good accuracy of the triggering means that every trigger point is located at the peak of the R-wave or at least always on the same one point of the R-wave. Bad triggering means that the trigger points are scattered over the ascending and descending part of the R-wave or even in a different part of the ECG, e.g. trigger points in the T-wave instead of the R-wave due to the magnetohydrodynamic effect (see Figure 4). To estimate this kind of error we decided to evaluate the triggering with the help of vectorcardiography (VCG). VCG is a 2D illustration of the electric conduction of the heart. In analogy to the ECG, the VCG consists of a P-, T- and an R-loop. The advantage of this method is that all 192 heart phases and their triggering can be plotted as a vector cloud in one single image (Figure 4). The clouds were evaluated in a blinded fashion as pairs in random order from three different investigators (CM, HS, SF), according to a predefined semiquantitive subjective grading of the ECG-triggering (1: “very good”, small scatter around the same area of the QRS-loop; 2: “good”, large scatter around the same area of the QRS-loop; 3: “bad”, < 5/192 (3%) trigger points in a different area of the QRS-loop; 4: “bad”, ≥ 5/192 (3%) trigger points in a different area of the QRS-loop).
The results of the measurement of the old and new ECG-trigger algorithm were compared by Bland-Altman analysis . The confidence interval was set to the mean value +/−1.96 of the standard deviation. Additional percentage lines were drawn to illustrate measurements with more than 5% and 10% difference.
Correlations between the changes in standard deviation, flow and accuracy
The standard deviation (SD) of the RR interval is a potential marker for unreliable triggering. Failed triggering will obviously lead to a higher SD. And probably a lower accuracy will also lead to a higher SD. Therefore a difference in SD between both measurements should be correlated with the power of the triggering method. We also expect a correlation between the SD and differences in the flow.
To compare these values, the differences in the SD and net flow were calculated as percentage of the mean value SD [%] and net flow [%].
Results of sensitivity, specificity and performance
Performance index (%)
old algorithm (N = 6811)
new algorithm (N = 6802)
Accuracy of triggering
Triggering was found to be more accurate using the new algorithm: Average grading for the old algorithm was 3.2 +/−1.0, average grading for the new algorithm was 1.5 +/−0.6 (p < 0.001). In 18 (51%) cases the accuracy of the new algorithm was remarkably better than the old (with a difference of grading of more than 2) and in 16 (46%) cases slightly better or equal (less than 2,0 difference). Only in one case, the accuracy of the old algorithm seemed to be better.
Correlations between the changes in SD, net flow and accuracy
The SD [%] differed between both algorithms (p = 0.028). The differences of SD [%] also correlated with differences in net flow [%] (p = 0.027 r = 0.373).
Furthermore, differences in grading as a marker for the accuracy correlated with differences in the SD [%] (p = 0.005 r = 0,462).
Due to the positive correlation between our appraisal of the accuracy and the difference in the flow quantifications we can assume that inefficient triggering has a considerable impact on the quality of flow quantification. However, other concomitant factors, which affect the flow in the ascending aorta, cannot totally be eliminated by our method. It is also difficult to achieve objective parameters for the evaluation of the accuracy. The obviously slightly better performance index, mainly based on missed QRS complexes, cannot statistically be correlated with the flow due to the very good performance index of both methods. Therefore, it is not possible to trace a specific difference in flow quantification totally back to an inefficient triggering. Nevertheless, we can say that the new algorithm works better under the same conditions than the old algorithm and because of that we can expect better results.
The difference in the flow quantification underlines the supposed complications. If we presume that a difference of more than 5% between two measurements probably leads to a clinically relevant misinterpretation and that a physical change in aortic flow under stable patient conditions can be neglected, 1/3 of the assessments probably fail due to incorrect measurements. The study shows that remarkable differences between the new and the old algorithms exist. The new algorithm showed significantly better triggering in all tested parameters. We can assume that these differences in triggering are an important source of error in the blood flow quantification. The new algorithm will probably help to achieve better results in many patients with CHD undergoing CMR.
- Fratz S, Hess J, Schuhbaeck A, Buchner C, Hendrich E, Martinoff S, Stern H: Routine clinical cardiovascular magnetic resonance in paediatric and adult congenital heart disease: patients, protocols, questions asked and contributions made. J Cardiovasc Magn Reson. 2008, 10: 46-10.1186/1532-429X-10-46.PubMed CentralView ArticlePubMed
- Powell AJ, Geva T: Blood flow measurement by magnetic resonance imaging in congenital heart disease. Pediatr Cardiol. 2000, 21: 47-58. 10.1007/s002469910007.View ArticlePubMed
- Friesen GM, Jannett TC, Jadallah MA, Yates SL, Quint SR, Nagle HT: A comparison of the noise sensitivity of nine QRS detection algorithms. IEEE Trans Biomed Eng. 1990, 37: 85-98. 10.1109/10.43620.View ArticlePubMed
- Gatzoulis MA, Till JA, Redington AN: Depolarization-repolarization inhomogeneity after repair of tetralogy of fallot. The substrate for malignant ventricular tachycardia?. Circulation. 1997, 95: 401-4. 10.1161/01.CIR.95.2.401.View ArticlePubMed
- Kastor JA, Goldreyer BN, Josephson ME, Perloff JK, Scharf DL, Manchester JH, Shelburne JC, Hirshfeld JW: Electrophysiologic characteristics of Ebstein’s anomaly of the tricuspid valve. Circulation. 1975, 52: 987-95. 10.1161/01.CIR.52.6.987.View ArticlePubMed
- Rößler J, Frank M: United States Patent application publication. Method for identifying a characteristic profile of an R-wave in an EKG signal and a computer program product as well as an electronically readable data medium for performing the method. 2011, Pub No: US 2011/10319780 A1; http://appft1.uspto.gov/netacgi/nph-Parser?Sect1=PTO1&Sect2=HITOFF&d=PG01&p=1&u=/netahtml/PTO/srchnum.html&r=1&f=G&l=50&s1=20110319780.PGNR
- Fratz S, Hess J, Schwaiger M, Martinoff S, Stern HC: More accurate quantification of pulmonary blood flow by magnetic resonance imaging than by lung perfusion scintigraphy in patients with fontan circulation. Circulation. 2002, 106: 1510-3. 10.1161/01.CIR.0000029103.26029.1E.View ArticlePubMed
- Frank M, Merkel S, Mustafa E, Wrobel H: United States Patent application publication. Method for controlling a measurement of a magnetic resonance device on the basis of an electrocardiogram signal. 2008, PUB no: US 2008/0139926 A1, June, 2008; http://appft1.uspto.gov/netacgi/nph-Parser?Sect1=PTO1&Sect2=HITOFF&d=PG01&p=1&u=/netahtml/PTO/srchnum.html&r=1&f=G&l=50&s1=20110319780.PGNR
- Rößler J, Frank M: United States Patent application publication. Method for identifying an R-wave in an ECG signal, ECG measuring device and magnetic resonance scanner. 2010, Pub No: US 2010/0191134 A1, July, 2010; http://appft1.uspto.gov/netacgi/nph-Parser?Sect1=PTO1&Sect2=HITOFF&d=PG01&p=1&u=/netahtml/PTO/srchnum.html&r=1&f=G&l=50&s1=20110319780.PGNR
- Proakis JG: Digital communications. 1995, New York: McGraw-Hill, 3 ed
- Bland JM, Altman DG: Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986, 1: 307-10.View ArticlePubMed
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.