Comparison of semi-automated methods to quantify infarct size and area at risk by cardiovascular magnetic resonance imaging at 1.5T and 3.0T field strengths
© Khan et al.; licensee BioMed Central. 2015
Received: 30 June 2014
Accepted: 9 February 2015
Published: 25 February 2015
There is currently no gold standard technique for quantifying infarct size (IS) and ischaemic area-at-risk (AAR [oedema]) on late gadolinium enhancement imaging (LGE) and T2-weighted short tau inversion recovery imaging (T2w-STIR) respectively. This study aimed to compare the accuracy and reproducibility of IS and AAR quantification on LGE and T2w-STIR imaging using Otsu’s Automated Technique (OAT) with currently used methods at 1.5T and 3.0T post acute ST-segment elevation myocardial infarction (STEMI).
Ten patients were assessed at 1.5T and 10 at 3.0T. IS was assessed on LGE using 5–8 standard-deviation thresholding (5-8SD), full-width half-maximum (FWHM) quantification and OAT. AAR was assessed on T2w-STIR using 2SD and OAT. Accuracy was assessed by comparison with manual quantification. Interobserver and intraobserver variabilities were assessed using Intraclass Correlation Coefficients and Bland-Altman analysis. IS using each technique was correlated with left ventricular ejection fraction (LVEF).
FWHM and 8SD-derived IS closely correlated with manual assessment at both field strengths (1.5T: 18.3 ± 10.7% LV Mass [LVM] with FWHM, 17.7 ± 14.4% LVM with 8SD, 16.5 ± 10.3% LVM with manual quantification; 3.0T: 10.8 ± 8.2% LVM with FWHM, 11.4 ± 9.0% LVM with 8SD, 11.5 ± 9.0% LVM with manual quantification). 5SD and OAT overestimated IS at both field strengths. OAT, 2SD and manually quantified AAR closely correlated at 1.5T, but OAT overestimated AAR compared with manual assessment at 3.0T. IS and AAR derived by FWHM and OAT respectively had better reproducibility compared with manual and SD-based quantification. FWHM IS correlated strongest with LVEF.
FWHM quantification of IS is accurate, reproducible and correlates strongly with LVEF, whereas 5SD and OAT overestimate IS. OAT accurately assesses AAR at 1.5T and with excellent reproducibility. OAT overestimated AAR at 3.0T and thus cannot be recommended as the preferred method for AAR quantification at 3.0T.
KeywordsMyocardial infarction Late gadolinium enhancement Infarct size T2-weighted STIR Myocardial salvage Myocardial oedema Ischaemic area at risk
Cardiovascular magnetic resonance (CMR)-measured infarct size (IS) [1,2] and myocardial salvage index (MSI) [3,4] are important measures of reperfusion success and predictors of remodelling and prognosis post acute ST-segment elevation myocardial infarction (STEMI). MSI is the proportion of reversibly injured ischaemic area-at-risk (AAR) visualised as myocardium with high signal intensity on T2-weighted images [3-5].
There is currently no gold standard technique for the quantification of IS and AAR on late gadolinium imaging (LGE) and T2-weighted short tau inversion recovery imaging (T2w-STIR) respectively. Semi-automated standard deviation (SD)-based thresholding techniques [4,7], manual (visual) contouring of enhancement [1,2], the full-width half-maximum (FWHM) method [6,8], and recently, automated techniques have been used [9,10]. The heterogeneity of techniques and resulting IS and AAR values hinders comparisons between studies.
Otsu’s Automated Thresholding (OAT) automatically identifies hyperenhanced areas by selecting the grayscale signal intensity threshold giving minimal intraclass variance within enhanced and normal myocardium and is largely user-independent . There are very scarce published data using OAT quantification, of IS  and AAR [13,14].
There are no published studies assessing IS or AAR quantification at 3.0T, or using 7SD and 8SD infarct quantification thresholding at any field strength.
This study aimed to compare the accuracy and reproducibility of IS and AAR quantification on LGE and T2w-STIR using OAT with the currently used quantification methods at 1.5T and 3.0T.
Ten patients were retrospectively, randomly selected using a random number generator  from the cohort of a UK multicentre, prospective CMR study investigating acute STEMI management at 1.5T (Complete Versus culprit-Lesion only PRimary PCI Trial) . Ten further patients were identically selected from a separate multicentre study at 3.0T (Randomized Controlled Trial Comparing Intracoronary Administration of Adenosine or Sodium Nitroprusside to Control for Attenuation of Microvascular Obstruction During Primary Percutaneous Coronary Intervention) . STEMI was diagnosed according to ESC definitions  and patients underwent primary PCI within 12 h of symptom onset. The studies were approved by Trent Research Ethics Committee, conducted according to the Declaration of Helsinki and all participants provided written informed consent.
CMR image acquisition
IS and AAR quantification
Image quality was graded according to a 4-point scale before analysis: 3 = excellent, 2 = good, 1 = moderate and 0 = unanalysable. To remove the confounding variable of image quality on AAR quantification, 26% of studies from the total study population, where T2w-STIR images were deemed non-analysable were excluded from the random number study selection pool. Analysis was performed offline in a central core lab, blinded to patient details using cmr42 (Circle Cardiovascular Imaging, Calgary, Canada). LGE, T2w-STIR and cine images were studied together and co-registered to allow accurate quantification based on all available data. For the assessment of LV volumes and function, IS and AAR, endocardial and epicardial borders were manually contoured on contiguous short-axis LV slices, excluding papillary muscles, trabeculae, epicardial surfaces and blood-pool artefact, and the quantification method applied. The most apical LGE and T2w-STIR slice was excluded to minimize partial volume effect. Total IS and AAR were expressed as percentage of LV mass (LVM).
AAR was quantified on T2w-STIR as hyperenhancement using 2SD thresholding and OAT by 2 blinded readers (JNK, SAN). The ROI size for 2SD was set at 2 cm2. Mean AAR was compared across the techniques and with manual quantification as described above for IS quantification.
Two manual corrections were applied to IS and AAR measurements: [a] inclusion of hypointense regions within enhancement corresponding to microvascular obstruction and intramyocardial haemorrhage in total IS and AAR respectively [4,6]; [b] exclusion of small isolated regions of enhancement without interslice continuity, in non-infarct related artery territories deemed to be noise artefact.
To assess intraobserver variability of the techniques, all images were re-quantified by a single observer after a 2-month interval. We also recorded the time taken to quantify IS and AAR using each of the methods once the endocardial and epicardial contours had been drawn (time taken for [a] quantification of AAR or IS using chosen technique + [b] inclusion of IMH or MVO where appropriate + [c] exclusion of noise artefact).
Normality was assessed using the Shapiro-Wilk test, histograms and Q-Q plots. Normally distributed data were expressed as mean ± standard deviation. IS and AAR by each quantification method were normally distributed and thus compared using paired t-tests, and the accuracy of each method was assessed by comparison with manual assessment using paired t-testing, two-way mixed-effect intraclass correlation coefficient for absolute agreement ICC (three, one)  and Bland-Altman analysis . On ICC, agreement was defined as excellent (ICC ≥ 0.75), good (ICC 0.6-0.74), fair (ICC 0.4-0.59), or poor (ICC < 0.40) . Interobserver and intraobserver variabilities were expressed using ICC (three, one) and Bland-Altman analysis. The significance of differences in reproducibility was assessed using Wilcoxon rank comparison of the squared differences . Statistical tests were performed using SPSS v20 (IBM, USA). p < 0.05 was considered significant.
Baseline demographics by CMR field strength cohort
56.6 ± 14.0
52.6 ± 10.6
LAD IRA (n, %)
LCX IRA (n, %)
RCA IRA (n, %)
IRA-only PCI: n = 3 (30%) Complete revascularisation: n = 7 (70%)
Vasodilator treatment group: n = 7 (70%) Control group: n = 3 (30%)
CMR time post STEMI (d)
3.7 ± 1.3
3.4 ± 2.1
111.6 ± 21.9
107.1 ± 23.1
179.9 ± 33.8
169.3 ± 35.2
94.7 ± 20.9
94.4 ± 32.3
47.2 ± 7.8
45.4 ± 8.1
LGE image quality
2.5 ± 0.6
2.2 ± 0.6
T2w-STIR image quality
2.6 ± 0.5
2.1 ± 0.3
Infarct size (IS) results at 1.5T by quantification method and corresponding reproducibilities
Mean IS (%LVM)
18.3 ± 10.7
25.9 ± 16.1b
22.0 ± 15.8a
19.8 ± 15.3
17.7 ± 14.4
28.2 ± 11.8b
16.5 ± 10.3
ICC v Manual
Mean bias v Manual (±1.96SD LoA)
+1.84 (+10.30, −6.62)
+9.39 (+25.58, −6.81)
+5.57 (+21.65, −10.52)
+3.28 (+18.86, −12.30)
+1.21 (+15.92, −13.50)
+11.71 (+17.39, +6.03)
Interobserver mean bias (±1.96SD LoA)
+0.37 (+9.17, −8.43)
+2.54 (+11.62, −6.54)
+4.43 (+16.01, −7.16)
+4.00 (+15.48, −7.47)
+4.01 (+16.02, −8.01)
+0.55 (+5.82, −4.73)
+5.34 (+14.96, −4.28)
Intraobserver mean bias (±1.96SD LoA)
+0.36 (+3.34, −2.61)
−0.81 (+9.37, −10.99)
+0.01 (+10.28, −10.27)
+0.07 (+11.76, −11.62)
+0.42 (+13.00, −12.17)
+0.81 (+4.21, −2.50)
−1.92 (+0.64, −4.49)
Infarct size (IS) results at 3.0T by quantification method and corresponding reproducibilities
Mean IS (%LVM)
10.8 ± 8.2
17.0 ± 11.2b
14.77 ± 10.4b
13.0 ± 9.7a
11.4 ± 9.0
21.6 ± 9.8b
11.5 ± 9.0
ICC v Manual
Mean bias v Manual (±1.96SD LoA)
+0.22 (+5.09, −4.65)
+6.42 (+14.93, −2.09)
+4.17 (+11.05, −2.71)
+2.38 (+7.92, −3.16)
+0.81 (+5.62, −4.01)
+11.03 (+22.20, −0.15)
Interobserver mean bias (±1.96SD LoA)
−0.49 (+1.74, −2.72)
+0.44 (+7.23, −6.35)
+1.14 (+8.51, −6.23)
+1.40 (+9.25, −6.44)
+1.50 (+9.99, −6.98)
−0.05 (+4.35, −4.44)
+1.97 (+9.48, −5.54)
Intraobserver mean bias (±1.96SD LoA)
+0.20 (+1.49, −1.10)
+0.43 (+2.90, −2.03)
−0.04 (+2.42, −2.50)
+0.10 (+2.19, −1.98)
+0.32 (+2.08, −1.45)
+0.15 (+3.50, −3.21)
+1.14 (+4.35, −2.07)
Area at risk (AAR) by field strength and quantification method and corresponding reproducibilities
Mean value (AAR [% LVM])
34.8 ± 9.8
38.1 ± 13.0
35.4 ± 11.2
35.2 ± 14.4
38.9 ± 9.9a
30.0 ± 8.2
ICC v Manual
Mean bias v Manual (±1.96SD LoA)
+0.31 (+12.20, −11.57)
+3.62 (+11.24, −4.00)
+5.13 (+22.76, −12.50)
+8.92 (+23.15, −5.31)
Interobserver mean bias (±1.96SD LoA)
+3.38 (+9.12, −2.37)
+1.26 (+7.20, −4.68)
+1.31 (+17.11, −14.49)
−0.34 (+15.36, −16.03)
−1.35 (+1.64, −4.33)
−5.20 (+9.54, −19.95)
Intraobserver mean bias (±1.96SD LoA)
+2.80 (+7.28, −1.68)
+0.58 (+3.31, −2.16)
−0.01 (+5.62, −5.64)
+1.00 (+5.50, −3.51)
+0.04 (+2.82, −2.74)
+1.46 (+14.84, −11.93)
Time taken per patient for Infarct Size (IS) and Area at risk (AAR) quantification by field strength and quantification method
17.1 ± 2.4 vs. 20.9 ± 5.5
1.5T: FWHM v 6SD (IS)
17.1 ± 2.4 vs. 19.4 ± 3.1
1.5T: FWHM v 7SD (IS)
17.1 ± 2.4 vs. 19.1 ± 3.7
1.5T: FWHM v 8SD (IS)
17.1 ± 2.4 vs. 19.6 ± 3.2
1.5T: FWHM v OAT (IS)
17.1 ± 2.4 vs. 18.0 ± 2.6
1.5T: FWHM v MANUAL (IS)
17.1 ± 2.4 vs. 21.1 ± 4.7
1.5T: 5SD v OAT (IS)
20.9 ± 5.5 vs. 18.0 ± 2.6
1.5T: 2SD v OAT (AAR)
17.1 ± 2.4 vs. 16.7 ± 2.6
1.5T: 2SD v MANUAL (AAR)
17.1 ± 2.4 vs. 18.3 ± 2.6
1.5T: OAT v MANUAL (AAR)
16.7 ± 2.6 vs. 18.3 ± 2.6
3T: FWHM v 5SD (IS)
18.9 ± 2.7 vs. 24.7 ± 9.1
3T: FWHM v 6SD (IS)
18.9 ± 2.7 vs. 22.2 ± 5.2
3T: FWHM v 7SD (IS)
18.9 ± 2.7 vs. 22.5 ± 5.3
3T: FWHM v 8SD (IS)
18.9 ± 2.7 vs. 21.2 ± 3.0
3T: FWHM v OAT (IS)
18.9 ± 2.7 vs. 20.7 ± 3.2
3T: FWHM v MANUAL (IS)
18.9 ± 2.7 vs. 24.0 ± 3.7
3T: 5SD v OAT (IS)
24.7 ± 9.1 vs. 20.7 ± 3.2
3T: 2SD v OAT (AAR)
19.6 ± 2.7 vs. 18.5 ± 2.6
3T: 2SD v MANUAL (AAR)
19.6 ± 2.7 vs. 18.3 ± 2.6
3T: OAT v MANUAL (AAR)
18.5 ± 2.6 vs. 18.3 ± 2.6
Interobserver and intraobserver variability of IS quantification
Results are displayed in Tables 2 and 3. FWHM and OAT demonstrated extremely high interobserver and intraobserver agreement at both field strengths, with all ICC values >0.922 and mean bias < +1.84%. SD-based techniques demonstrated good interobserver and intraobserver agreement at both field strengths, however lower than for FWHM and OAT, with ICC values >0.888 and mean bias < +4.43%. Interobserver and intraobserver agreement for manual quantification were very high at both field strengths apart from interobserver agreement at 1.5T, which was good (ICC 0.793). Bland-Altman charts for IS are shown in Additional files 2 and 3.
Interobserver agreement for IS at 3.0T was significantly better with FWHM vs. manual quantification (p = 0.037). Intraobserver agreement for IS was significantly better at 1.5T with FWHM vs. 6SD (p = 0.013), 7SD (p = 0.022) and 8SD (p = 0.037), and at 3.0T for FWHM vs. manual (p = 0.047). There was a strong trend towards higher intraobserver agreement for IS at 1.5T with FWHM vs. manual (p = 0.093).
Correlation of myocardial injury with LV ejection fraction
At 1.5T, FWHM and manual quantification demonstrated a strong inverse correlation between IS and LVEF (FWHM: r = −0.745, p = 0.013; manual r = −0.709, p = 0.022). All other methods demonstrated moderate inverse correlation and did not reach statistical significance. At 3.0T, FWHM IS showed a significant, moderate correlation with LVEF (r = −0.673, p = 0.033). The correlation using all other techniques was weaker and not statistically significant.
AAR varied with the quantification method used (Figures 4 and 5). There was no significant difference between 2SD, OAT and manually quantified AAR at 1.5T. At 3.0T, AAR quantified with OAT was larger than that manually contoured (p = 0.004) and similar to that on 2SD. Agreement with manual quantification at 1.5T tended to be higher for OAT than 2SD, with ICC 0.920 and narrower limits of agreement on Bland-Altman analysis. There was no difference in AAR quantification time using OAT, 2SD or manual quantification at 1.5T or 3.0T (Table 5), Additional files 4 and 5.
Interobserver and intraobserver variability of AAR and MSI quantification
OAT had extremely high interobserver and intraobserver agreement for AAR quantification at both field strengths, with all ICC values >0.976. Good interobserver agreement was seen for 2SD quantification of AAR at both field strengths. Manual quantification demonstrated excellent interobserver agreement at 3.0T. Interobserver agreement at 1.5T and intraobserver agreement at both field strengths was good with manual quantification (ICC >0.716).
Interobserver agreement at 3.0T was significantly better for OAT vs. manual quantification (p = 0.017), and at 1.5T was borderline significantly higher for OAT vs. manual (p = 0.059). Intraobserver agreement at 3.0T was significantly better for OAT vs. manual quantification (p = 0.007). The raw datasets for IS and AAR quantification at 1.5T and 3.0T field strengths are available in Additional files 6 and 7.
In this study we assessed IS and AAR quantification in acute STEMI patients with currently used semi-automatic techniques at 1.5T and 3.0T. FWHM and 8SD closely agreed with the reference standard of manual IS quantification at both field strengths, whereas 5SD and OAT led to higher IS values compared to manual quantification. AAR measured by OAT and 2SD were similar to manual quantification only at 1.5T. Interobserver and intraobserver agreement for IS and AAR quantification were better with FWHM and OAT compared with manual quantification respectively, and tended to be better than with SD-based methods. There was an inverse correlation between IS and LVEF for all quantification methods and this was strongest and most significant for FWHM. Our study is the first to assess IS quantification methods using 7SD and 8SD thresholding and to assess IS and AAR quantification at 3.0T.
Mean IS using the quantification techniques
LGE IS quantification in acute MI has been validated in a small number of animal studies. FWHM  and manual quantification  of in-vivo images closely correlated with IS on tetrazolium chloride stained canine hearts. Kim et al.  demonstrated good agreement of 2SD thresholding with tetrazolium chloride stained canine myocardium. However this was on ex-vivo slices with high spatial resolution and in the absence of rhythm and motion artefacts, and may not be generalizable to humans . Indeed, 2SD has been shown to overestimate IS in humans based on functional improvement and IS reduction in enhanced areas.[26,27] There is no histological validation in humans and hence no ‘gold standard’ quantification. We thus used manual assessment as has been used previously [6,12], however derived from the mean of repeated analyses by three experienced CMR cardiologists to increase the robustness of our reference standard.
FWHM and 8SD were the only methods in our study showing good agreement with manual quantification at both field strengths. This may be because they are less prone to IS overestimation resulting from oedema and partial volume effects giving rise to intermediate signal intensities [26,28]. This resulted in negligible requirements for manual exclusion of noise artefact with FWHM and 8SD. This in conjunction with the relative ease in identification of the brightest infarct core compared with deciding on a representative remote ROI is likely to explain the shorter time required for IS quantification using FWHM compared with SD-based techniques.
The greater IS using 5SD compared with manual quantification in our study is in agreement with previous results at 1.5T . These findings indicate that the good agreement between 5SD and manual quantification in chronic ischaemic heart disease , where infarct tends to have a higher and more homogenous SI , cannot be extrapolated to acute STEMI patients. The close correlation of 5SD and in particular OAT with manual assessment shown by Vermes et al.  is in contrast to our findings. IS quantification was only performed on slices with infarct seen visually in that study, thus potentially underestimating IS. In addition, the small remote ROI used for 5SD thresholding (0.5-1 cm2) by Vermes et al. may not adequately represent remote myocardium signal intensity, thus leading to underestimation or overestimation of IS if an excessively bright or dark, isolated region of myocardium is taken as the remote ROI respectively. By setting the ROI size at 2 cm2 for all SD-based methods in our study, we aimed to ensure that the ROI was large enough to represent remote myocardium accurately. Using the same remote ROI for all SD-based methods in our study ensured consistency and removed the effect of ROI size and location when comparing IS between 5–8 SD thresholds. Hence, 6 and 5-SD and 7, 6 and 5-SD quantification overestimated IS at 1.5T and 3.0T respectively due to their intrinsically progressively lower signal intensity thresholds and not due to differences in remote ROI.
We studied IS and AAR quantification early after STEMI. IS decreases with time post PPCI with a reduction of ~30% demonstrated within the first week in some studies [26,30]. The extent of necrosis is overestimated by LGE early post STEMI due to cellular disruption and oedema. As scar resorbs and remodels, IS reduces and scar may become more homogenous in signal intensity and brighter. The relative overestimation of IS by lower standard deviation thresholds and OAT compared with FWHM, 8SD and manual quantification may thus be more significant in acute compared with in chronic infarcts. We chose an early time point to minimise drop-out in the study and most importantly, all the data relating infarct size to subsequent prognosis following STEMI has been based on early measurement of infarct size (usually within 1 week) [4,31]. Whether AAR varies in the first week after STEMI has shown conflicting results [26,27]. As we have only scanned the patients in this study on a single occasion we cannot comment on how the results would have varied if performed at later dates following presentation.
Interobserver and intraobserver variability of IS quantification
The excellent interobserver and intraobserver agreement for FWHM, 5SD and OAT quantification of IS in our study at both field strengths is in agreement with previous studies at 1.5T: FWHM, 5SD [6,27] and FWHM, OAT . Consistent with the work of Flett et al. , we found that the FWHM technique had greater interobserver and intraobserver reproducibility compared with SD-based and manual quantification. This is expected when considering that for each patient there is a single brightest core of infarct. This is in contrast to the remote ROI, which could be drawn on any slice without complete LGE in SD-based quantification, or manual contouring of enhancement, which is completely user-dependent and in the acute phase post STEMI could potentially be affected by partial volume in infarct boundaries and oedematous myocardium.
Mean AAR and MSI using the quantification techniques
T2w-STIR AAR is typically quantified using 2SD thresholding. Validation studies are limited. 2SD-derived AAR on T2w images in vivo correlated with microsphere-assessed AAR in canine myocardium (r = 0.84). There is no gold standard AAR quantification method on T2w-STIR, hence we used manual assessment.
The close agreement between OAT and manually contoured AAR at 1.5T is consistent with the work of McAlindon et al. OAT however demonstrated greater AAR compared with manual quantification at 3.0T. This is in keeping with Sjogren et al.  who showed overestimation of AAR using OAT with a mean bias of +5.3 ± 9.6% compared with manual quantification . The determination of an optimal threshold and quantification of enhancement on every slice with OAT, regardless of oedema is likely to contribute to this. The risk of overestimation of AAR will be greatest in slices with minimal oedema since OAT will deem a proportion of pixels enhanced. This may potentially have contributed to the overestimation of AAR at 3.0T in our study with OAT, since our 3.0T cohort had a smaller AAR than the 1.5T patients. IS was also smaller in our 3.0T cohort and may have contributed to the greater overestimation of IS using OAT at 3.0T compared with 1.5T. Conversely, underestimation of AAR is more likely in slices with complete enhancement since OAT will deem a proportion of pixels unenhanced . T2w-STIR images were of diagnostic quality in all patients in our study, however mean quality control grading was slightly lower at 3.0T (2.1 ± 0.3 [3 T] vs. 2.6 ± 0.5 [1.5T], p = 0.05) and may have potentially contributed to the overestimation of OAT-derived AAR if there was more noise artefact in the AAR or signal intensity drop out in remote regions by reducing the threshold. More work into automated quantification methods is required, in particular at 3.0T. Newer automated techniques, taking into account a priori information about the culprit artery  and including noise and false positive artefact exclusion [10,21] algorithms may improve the accuracy of automated IS and AAR quantification.
The relative degree of AAR overestimation in our study was, however, considerably less than for IS. The predominantly transmural pattern of OAT enhancement for IS and AAR may cause less overestimation of AAR compared with IS, since oedema has been shown to be predominantly transmural in 70-100% of oedematous segments [33,34].
Interobserver and intraobserver variability of AAR quantification
The relatively low interobserver and intraobserver agreement using 2SD compared with OAT at both field strengths is likely to result from varying manual definition of the remote ROI. The extremely high ICC’s obtained with OAT are remarkable considering that these figures still take account differences in manual correction and contouring of endocardial and epicardial borders. Given these results, quantification of AAR with OAT could minimise variability in measurement in multi-centre trials.
The main limitation of our and previous similar studies is the lack of a gold standard for IS and AAR quantification. Different quantification techniques were studied for IS and AAR. FWHM quantification of AAR was not undertaken due to the lower CNR of T2w-STIR imaging, since the vast majority of myocardium would have signal intensity >50% of the maximum at the AAR core, leading to potentially extreme overestimation of AAR and MSI. Indeed, McAlindon et al. demonstrated that FWHM significantly overestimated AAR compared to all other quantification methods tested at 1.5T (2,3,5 SD, OAT, manual quantification) . 5SD thresholding was not assessed for AAR as it has never been validated or correlated with clinical outcomes and the only study to feature it demonstrated that it significantly overestimated AAR compared to all other quantification methods tested at 1.5T (2,3 SD, FWHM, OAT, manual quantification) . 2SD thresholding was not assessed for IS as it has been shown to overestimate IS [6,12] and had the lowest correlation with histological IS on tetrazolium chloride staining using Bland-Altman analysis  compared with all other quantification methods used in studies of IS in acute STEMI (5SD, FWHM, manual quantification). Test-retest reproducibility was not assessed and should be considered in future studies. Infarct heterogeneity and identification of peri-infarct zone (greyzone) was not assessed in this study and may be of interest to assess in future studies using OAT. We deliberately studied patients imaged at different field strengths and with different scanner vendors to represent the situation in multi-centre clinical trials and this should make the results more generalizable. Our sample size (total n = 20) is limited, however is comparable to similar studies in myocardial infarction [8-10,12,14,20] and our findings are supported by their consistency at both field strengths.. Finally our results may not be generalizable to if patients are scanned at different time points following STEMI.
Inter- and intraobserver variability for the quantification of IS with FWHM is excellent at 1.5 and 3.0T and better than when using manual quantification. Only FWHM and 8SD closely agreed with manual delineation of IS at both field strengths. FWHM had better reproducibility, shorter quantification time and closer correlation with LVEF and may be the preferred method for IS quantification in future studies. AAR is similar when assessed with OAT, 2SD and manual quantification at 1.5T, however OAT has excellent intra and interobserver variability and thus has potential in quantification of AAR at 1.5T, especially in multi-centre studies. OAT overestimated AAR at 3.0T compared with manual quantification and thus cannot currently be recommended as the preferred method for AAR quantification at 3.0T.
GPM is supported by a NIHR Intermediate Fellowship. This study was supported by The National Institute for Health Research (funded the CMR scans and software used in image analysis, funding number EME 10/27/01) and British Heart Foundation (funded recruitment of subjects, funding number SP/10/001/28194).
- Klem I, Shah DJ, White RD, Pennell DJ, van Rossum AC, Regenfus M, et al. Prognostic value of routine cardiac magnetic resonance assessment of left ventricular ejection fraction and myocardial damage: an international, multicenter study. Circ Cardiovasc Imaging. 2011;4:610–9.View ArticlePubMedGoogle Scholar
- Hombach V, Grebe O, Merkle N, Waldenmaier S, Höher M, Kochs M, et al. Sequelae of acute myocardial infarction regarding cardiac structure and function and their prognostic significance as assessed by magnetic resonance imaging. Eur Heart J. 2005;26:549–57.View ArticlePubMedGoogle Scholar
- Masci PG, Ganame J, Strata E, Desmet W, Aquaro GD, Dymarkowski S, et al. Myocardial Salvage by CMR Correlates With LV Remodeling and Early ST-Segment Resolution in Acute Myocardial Infarction. JCMG. 2010;3:45–51.Google Scholar
- Eitel I, Desch S, de Waha S, Fuernau G, Gutberlet M, Schuler G, et al. Long-term prognostic value of myocardial salvage assessed by cardiovascular magnetic resonance in acute reperfused myocardial infarction. Heart (British Cardiac Soc). 2011;97:2038–45.View ArticleGoogle Scholar
- Pennell D. Myocardial salvage: retrospection, resolution, and radio waves. Circulation. 2006;113:1821–3.View ArticlePubMedGoogle Scholar
- Flett AS, Hasleton J, Cook C, Hausenloy D, Quarta G, Ariti C, et al. Evaluation of techniques for the quantification of myocardial scar of differing etiology using cardiac magnetic resonance. JACC Cardiovasc Imaging. 2011;4:150–6.View ArticlePubMedGoogle Scholar
- Kim RJ, Fieno DS, Parrish TB, Harris K, Chen E-l, Simonetti O, et al. Relationship of MRI Delayed Contrast Enhancement to Irreversible Injury, Infarct Age, and Contractile Function. Circulation. 1999;100:1992–2002.View ArticlePubMedGoogle Scholar
- Amado LC, Gerber BL, Gupta SN, Rettmann DW, Szarf G, Schock R, et al. Accurate and objective infarct sizing by contrast-enhanced magnetic resonance imaging in a canine myocardial infarction model. J Am Coll Cardiol. 2004;44:2383–9.View ArticlePubMedGoogle Scholar
- Hsu LY, Natanzon A, Kellman P, Hirsch GA, Aletras AH, Arai AE. Quantitative myocardial infarction on delayed enhancement MRI. Part I: animal validation of an automated feature analysis and combined thresholding infarct sizing algorithm. J Magn Reson Imaging. 2006;23:298–308.View ArticlePubMedGoogle Scholar
- Johnstone RI, Greenwood JP, Biglands JD, Plein S, Ridgway JP, Radjenovic A. Assessment of tissue edema in patients with acute myocardial infarction by computer-assisted quantification of triple inversion recovery prepared MRI of the myocardium. Magn Reson Med. 2011;66:564–73.View ArticlePubMedGoogle Scholar
- Otsu N. A Threshold Selection Method from Gray-Level Histograms. IEEE Trans Sys Man Cybern. 1979;5MC-9:62–6.Google Scholar
- Vermes E, Childs H, Carbone I, Barckow P, Friedrich M. Auto-Threshold Quantification of Late Gadolinium Enhancement in Patients With Acute Heart Disease. JMRI. 2013;37:382–90.View ArticlePubMedGoogle Scholar
- Sjogren J, Ubachs JF, Engblom H, Carlsson M, Arheden H, Heiberg E. Semi-automatic segmentation of myocardium at risk in T2-weighted cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2012;14:10.View ArticlePubMed CentralPubMedGoogle Scholar
- McAlindon E, Lawton C, Flett AS, Manghat N, Hamilton MC, Bucciarelli Ducci C. Evaluation of 7 techniques for the quantification of myocardial oedema in STEMI. JCMR. 2013;15:P188.Google Scholar
- Haahr M. Random.org (True Random Number Generator). http://www.random.org 1998, Website (Accessed September 2013).
- Kelly DJ, McCann GP, Blackman D, Curzen NP, Dalby M, Greenwood JP, et al. Complete Versus culprit-Lesion only PRimary PCI Trial (CVLPRIT): a multicentre trial testing management strategies when multivessel disease is detected at the time of primary PCI: rationale and design. EuroIntervention. 2013;8:1190–8.View ArticlePubMedGoogle Scholar
- Nazir S, Khan JN, Mahmoud IZ, Greenwood JP, Blackman D, Kunadian V, et al. The REFLO-STEMI trial comparing intracoronary adenosine, sodium nitroprusside and standard therapy for the attenuation of infarct size and microvascular obstruction during primary percutaneous coronary intervention: study protocol for a randomised controlled trial. Trials. 2014;15:371.View ArticlePubMed CentralPubMedGoogle Scholar
- Steg PG, James SK, Atar D, Badano LP, Lundqvist CB, Borger MA, et al. ESC Guidelines for the management of acute myocardial infarction in patients presenting with STsegment elevation: The Task Force on the management of ST-segment elevation acute myocardial infarction of the European Society of Cardiology (ESC). Eur Heart J. 2012;33(20):2569–619.View ArticlePubMedGoogle Scholar
- Heiberg E, Engblom H, Engvall J, Hedstrom E, Ugander M, Arheden H. Semi-automatic quantification of myocardial infarction from delayed contrast enhanced magnetic resonance imaging. Scand Cardiovasc J. 2005;39:267–75.View ArticlePubMedGoogle Scholar
- Hsu LY, Ingkanisorn WP, Kellman P, Aletras AH, Arai AE. Quantitative myocardial infarction on delayed enhancement MRI. Part II: Clinical application of an automated feature analysis and combined thresholding infarct sizing algorithm. J Magn Reson Imaging. 2006;23:309–14.View ArticlePubMedGoogle Scholar
- Gao H, Kadir K, Payne AR, Soraghan J, Berry C. Highly automatic quantification of myocardial oedema in patients with acute myocardial infarction using bright blood T2-weighted CMR. J Cardiovasc Magn Reson. 2013;15:28.View ArticlePubMed CentralPubMedGoogle Scholar
- Kadir K, Gao H, Payne A, Soraghan J, Berry C. Automatic quantification and 3D visualisation of edema in cardiac MRI. Conf Proc IEEE Eng Med Biol Soc. 2011;2011:8021–4.PubMedGoogle Scholar
- McGraw KO, Wong SP. Forming Inferences About Some Intraclass Correlation Coefficients. Psychol Methods. 1996;1:30–46.View ArticleGoogle Scholar
- Bland JM, Altman DG. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986;1:307–10.View ArticlePubMedGoogle Scholar
- Castillo E, Osman NF, Rosen BD, El-Shehaby I, Pan L, Jerosch-Herold M, et al. Quantitative assessment of regional myocardial function with MR-tagging in a multi-center study: interobserver and intraobserver agreement of fast strain analysis with Harmonic Phase (HARP) MRI. J Cardiovasc Magn Reson. 2005;7:783–91.View ArticlePubMedGoogle Scholar
- Mather AN, Fairbairn TA, Artis NJ, Greenwood JP. Timing of Cardiovascular MR Imaging after Acute Myocardial Infarction : Effect on Estimates of Infarct Characteristics and Prediction of Late Ventricular Remodeling. Radiology. 2011;261:116–26.View ArticlePubMed CentralPubMedGoogle Scholar
- Dall’Armellina E, Karia N, Lindsay AC, Karamitsos TD, Ferreira V, Robson MD, et al. Dynamic changes of edema and late gadolinium enhancement after acute myocardial infarction and their relationship to functional recovery and salvage index. Circ Cardiovasc Imaging. 2011;4:228–36.View ArticlePubMed CentralPubMedGoogle Scholar
- Kim HW, Farzaneh-Far A, Kim RJ. Cardiovascular magnetic resonance in patients with myocardial infarction: current and emerging applications. J Am Coll Cardiol. 2009;55:1–16.View ArticlePubMedGoogle Scholar
- Bondarenko O, Beek AM, Hofman MB, Kuhl HP, Twisk JW, van Dockum WG, et al. Standardizing the definition of hyperenhancement in the quantitative assessment of infarct size and myocardial viability using delayed contrast-enhanced CMR. J Cardiovasc Magn Reson. 2005;7:481–5.View ArticlePubMedGoogle Scholar
- Ibrahim T, Hackl T, Nekolla SG, Breuer M, Feldmair M, Schömig A. Acute Myocardial Infarction : Serial Cardiac MR Imaging Shows a Decrease in Delayed Enhancement of the Myocardium during the 1st Week after Reperfusion. Radiology. 2010;254:88–97.View ArticlePubMedGoogle Scholar
- Wu E, Ortiz JT, Tejedor P, Lee DC, Kansal P, Carr JC, et al. Infarct size by contrast enhanced cardiac magnetic resonance is a stronger predictor of outcomes than left ventricular ejection fraction or end-systolic volume index : prospective cohort study Infarct size by contrast enhanced cardiac magnetic resonance i. Heart. 2008;94:730–6.View ArticlePubMedGoogle Scholar
- Aletras AH, Tilak GS, Natanzon A, Hsu L-Y, Gonzalez FM, Hoyt RF, et al. Retrospective determination of the area at risk for reperfused acute myocardial infarction with T2-weighted cardiac magnetic resonance imaging: histopathological and displacement encoding with stimulated echoes (DENSE) functional validations. Circulation. 2006;113:1865–70.View ArticlePubMedGoogle Scholar
- OH-Ici D, Ridgway JP, Kuehne T, Berger F, Plein S, Sivananthan M, et al. Cardiovascular magnetic resonance of myocardial edema using a short inversion time inversion recovery (STIR) black-blood technique: Diagnostic accuracy of visual and semi-quantitative assessment. JCMR. 2012;14:22.Google Scholar
- Friedrich MG, Abdel-Aty H, Taylor A, Schulz-Menger J, Messroghli D, Dietz R. The salvaged area at risk in reperfused acute myocardial infarction as visualized by cardiovascular magnetic resonance. J Am Coll Cardiol. 2008;51:1581–7.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.