Skip to main content

125 Characterization of infarct heterogeneity by fuzzy clustering


Infarcted tissue interspersed with surviving bundles of myocytes can cause ventricular arrhythmias via the reentry phenomena. Infarct heterogeneity as assessed by delayed-enhancement MRI (DE-MRI) has been shown to correlate with inducibility for ventricular tachycardia (VT) [1]. However, the separation of the infarct core from the periphery (or "gray zone") using conventional DE-MRI images relies on an arbitrary selection of a signal intensity (SI) cut-off. DE-MRI images generally have a low SNR; therefore, noise may influence the size of the gray zone. Recently, a new method using an inversion-recovery SSFP (IR-SSFP) sequence has been described to produce multiple infarct-enhanced images at various inversion times [2].


In this work a new method for assessing infarct heterogeneity based on parameterization and fuzzy clustering of IR-SSFP images was developed.

Methods and results

The IR-SSFP sequence was implemented on a 1.5 T scanner (CV/i, GE Healthcare, Milwaukee, WI). Short-axis images were acquired in patients with suspected myocardial infarcts 10–20 minutes after the injection of 0.2 mmol/kg of Gd-DTPA (Magnevist, Berlex Inc., Wayne, NJ). The parameters for the IR-SSFP sequence were: BW = ± 125 kHz, TR/TE = 2.7/1.3 ms, SSFP flip angle = 45° and spatial resolution of 1.6 × 1.6 × 8 mm over a 32 cm FOV. An eight-channel cardiac array was used for signal reception. The IR-SSFP sequence yielded 20 images throughout the heart cycle each at a different inversion time. The first nine images after the inversion pulse occur during diastole and thus have minimal motion; these nine images were used for a pixel-by-pixel analysis of the inversion recovery behaviour. For each pixel, the time course of the SI during its signal recovery was extracted (Fig 1). The SIs were then fitted to the signal recovery equation for an IR-SSFP sequence [3]. The two unknown parameters in the signal recovery equation are: the apparent relaxation time T1* (shorter than the true T1 because of the effects of the SSFP readout pulses) and the steady-state (SS) value. From the regression analysis, pixel-by-pixel T1* and SS parameter maps were calculated. The parameter maps in Figure 1 show a non-transmural antero-septal infarct as an area of reduced T1* (dark band in the T1* map) and increased SS. A scatter plot of the T1* versus SS value for all pixels can then be displayed (Fig 2). This scatter plot was used as the input to a clustering algorithm designed to automatically detect and classify each pixel. The clustering algorithm used was the Gustafson-Kessel (GK) modification of the fuzzy C-means algorithm [4], which uses only the raw data and the number of clusters as inputs. For each pixel, the algorithm calculated the probabilities of that pixel being associated with infarct, healthy myocardium, and blood. Pixels having a probability of greater than 75% for belonging to infarct or myocardium were classified as infarct core or healthy myocardium, respectively; pixels with probabilities between 25–75% for both infarct and myocardium were classified as heterogeneous infarct and constitute the gray zone. Clustering results from one patient are shown in Figure 2 and demonstrate the infarct core (in green) and gray zone (in yellow).

Figure 1
figure 1

Signal intensities (left) for one pixel each in healthy myocardium, infarct and blood from nine IR-SSFP images, each at a different inversion time. The fitted recovery curves on a pixel-by-pixel basis yield parameter maps for T1* (middle) and steady-state values (right).

Figure 2
figure 2

A T1* versus SS scatter plot (left) and corresponding image (right) showing the fuzzy clustering results with infarct core (green), gray zone (yellow), myocardium (blue) and blood (red). Inversion-recovery SSFP yields multiple delayed enhancement images at different inversion times that have been used to create pixel-by-pixel parameter maps. Fuzzy clustering methods were developed to use the paramter maps to detect areas of heterogeneous infarct.


Fuzzy clustering applied to IR-SSFP images can be used to classify pixels as infarct core, infarct grey zone, healthy myocardium and blood. Future work is needed to determine if heterogeneity measured by this method is predictive of VT, and to compare this method to the SI threshold method using conventional DE-MRI images.


  1. Schmidt A: Circ. 2007, 115: 2006-14. 10.1161/CIRCULATIONAHA.106.653568.

    Article  Google Scholar 

  2. Detsky JS: Magn Reson Med. 2007, 58: 365-72. 10.1002/mrm.21291.

    Article  CAS  PubMed  Google Scholar 

  3. Schmitt P: Magn Reson Med. 2004, 51: 661-7. 10.1002/mrm.20058.

    Article  PubMed  Google Scholar 

  4. Gustafson DE: Proc of IEEE CDC. 1979, 761-6.

    Google Scholar 

Download references

Author information

Authors and Affiliations


Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution 2.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Detsky, J.S., Stainsby, J.A., Dick, A.J. et al. 125 Characterization of infarct heterogeneity by fuzzy clustering. J Cardiovasc Magn Reson 10 (Suppl 1), A26 (2008).

Download citation

  • Published:

  • DOI: