Semi-automated analysis of infarct heterogeneity on DE-MRI using graph cuts
Journal of Cardiovascular Magnetic Resonance volume 14, Article number: T6 (2012)
Two popular methods for determining the threshold values for the infarct core and gray zone on delayed enhancement MR images (DE-MRI) have been proposed previously: full width and half maximum (FWHM) and standard deviation (SD) methods . Major limitations of these methods are:1) three manually drawn contours are needed for endocardial, epicardial and remote myocardium boundaries, which is time consuming and suffers from inter-observer and intra-observer variability; 2) the difficulty in reproducible manual delineation of remote myocardium, is an important contributor to variability in results; and 3) the dependence on the remote region statistics is problematic due to the low SNR of this region . The purpose of this research was to develop a novel algorithm for segmentation of infarct core and gray zone from conventional IR-GRE short-axis MR images with highly robust and reproducible results comparable to the FWHM analysis while eliminating the requirement for drawing a remote myocardial region.
Eleven male patients (age: 63.5±11.8 yr) with known CAD and evidence of LGE and chronic MI had cardiac IR-GRE MR scans with full left ventricle (LV) coverage (7-13 slices). MR imaging was performed on a 1.5 T scanner (CV/i, GE Healthcare) using an 8-channel cardiac coil. DE-MRI was started 10 min after the injection of 0.2 mmol/kg of Gd-DTPA (Magnevist, Berlex). First, the endocardial and epicardial contours were generated from the corresponding SSFP images automatically  with papillary muscles included. (42%(54/136) of the contours need manual adjustment. Next, the graph cuts algorithm  was used to segment the infarct: 1) a two-classes Gaussian Mixture Model (GMM) was created for the myocardial ROI, determined by the endocardial and epicardial contours; 2) each pixel was assigned to the most likely Gaussian component; and 3) a graph was built and the graph cut algorithm finds the optimum classification of healthy and infarcted myocardium. Finally, the segmented infarct is separated into infarct core and gray zone using a threshold of half the maximal signal within the segmented infarct (Fig. 1). Linear regression analysis and Bland-Altman plots were used to compare the FWHM method (requiring manual ROIs for LV and remote myocardium) and our method (Fig. 2).
There were excellent correlations of the infarct size (infarct core 1: R^2 = 0.99; gray zone: R^2 = 0.95) derived from our graph cuts method and the manual FWHM method. The Bland-Altman analysis indicated that there was a small overestimation bias (infarct core: 0.17 g; gray zone: 0.68 g) with limits of agreement of ± 1.40 g (infarct core) and ± 4.27 g (gray zone). This variability is small relative to the reported range of gray zone masses (20 +/- 13 g, N=91 ).
The results for the proposed semi-automated segmentation technique indicate that it will streamline accurate quantification of myocardial infarct on IR-GRE MR infarct images in clinical practice.
The authors thank the Canadian Foundation for Innovation (CFI), the Canadian Institutes of Health Research (CIHR),and MaRS Innovation Proof of Principle Program for their support.
Flett AS, et al: JACC: cardiovascular imaging. 2011, 4 (2): 150-156.
Roes SD, et al: Circulation: cardiovascular imaging. 2009, 2: 183-190.
Lu Y, et al: Segmentation of Left Ventricle in Cardiac Cine MRI: An Automatic Image-Driven Method. Lecture Notes in Computer Science. 2009, 5528/2009: 339-347.
Boykov Y, et al: Fast approximate energy minimization via graph cuts. IEEE transactions PAMI. 2001, 23 (11): 1222-1239. 10.1109/34.969114.
About this article
Cite this article
Lu, Y., Connelly, K.A., Yang, Y. et al. Semi-automated analysis of infarct heterogeneity on DE-MRI using graph cuts. J Cardiovasc Magn Reson 14 (Suppl 1), T6 (2012). https://doi.org/10.1186/1532-429X-14-S1-T6