- Open Access
Semi-automatic segmentation of myocardium at risk in T2-weighted cardiovascular magnetic resonance
© Sjögren et al; licensee BioMed Central Ltd. 2012
- Received: 23 August 2011
- Accepted: 31 January 2012
- Published: 31 January 2012
T2-weighted cardiovascular magnetic resonance (CMR) has been shown to be a promising technique for determination of ischemic myocardium, referred to as myocardium at risk (MaR), after an acute coronary event. Quantification of MaR in T2-weighted CMR has been proposed to be performed by manual delineation or the threshold methods of two standard deviations from remote (2SD), full width half maximum intensity (FWHM) or Otsu. However, manual delineation is subjective and threshold methods have inherent limitations related to threshold definition and lack of a priori information about cardiac anatomy and physiology. Therefore, the aim of this study was to develop an automatic segmentation algorithm for quantification of MaR using anatomical a priori information.
Forty-seven patients with first-time acute ST-elevation myocardial infarction underwent T2-weighted CMR within 1 week after admission. Endocardial and epicardial borders of the left ventricle, as well as the hyper enhanced MaR regions were manually delineated by experienced observers and used as reference method. A new automatic segmentation algorithm, called Segment MaR, defines the MaR region as the continuous region most probable of being MaR, by estimating the intensities of normal myocardium and MaR with an expectation maximization algorithm and restricting the MaR region by an a priori model of the maximal extent for the user defined culprit artery. The segmentation by Segment MaR was compared against inter observer variability of manual delineation and the threshold methods of 2SD, FWHM and Otsu.
MaR was 32.9 ± 10.9% of left ventricular mass (LVM) when assessed by the reference observer and 31.0 ± 8.8% of LVM assessed by Segment MaR. The bias and correlation was, -1.9 ± 6.4% of LVM, R = 0.81 (p < 0.001) for Segment MaR, -2.3 ± 4.9%, R = 0.91 (p < 0.001) for inter observer variability of manual delineation, -7.7 ± 11.4%, R = 0.38 (p = 0.008) for 2SD, -21.0 ± 9.9%, R = 0.41 (p = 0.004) for FWHM, and 5.3 ± 9.6%, R = 0.47 (p < 0.001) for Otsu.
There is a good agreement between automatic Segment MaR and manually assessed MaR in T2-weighted CMR. Thus, the proposed algorithm seems to be a promising, objective method for standardized MaR quantification in T2-weighted CMR.
- Cardiovascular Magnetic Resonance
- Left Anterior Descend
- Left Ventricular Mass
- Right Coronary Artery
- Inter Observer Variability
Myocardium at risk (MaR) is defined as the ischemic myocardium during coronary artery occlusion and is the region that will be subject to infarction if the blood flow is not restored. Myocardium at risk can be measured using T2-weighted cardiovascular magnetic resonance(CMR)  due to the myocardial edema occurring in the ischemic myocardium [2, 3] up to one week after percutaneous coronary intervention (PCI) . By determining MaR using T2-weighted CMR and myocardial infarction (MI) size using late gadolinium enhancement (LGE), the efficacy of reperfusion therapy can be assessed as myocardial salvage in a single CMR session.
In the event of an acute coronary occlusion, a single artery is usually affected. As a consequence of the occlusion, transmural ischemia occurs within the affected coronary artery's perfusion territory [5, 6]. The myocardium subjected to ischemia becomes edematous and shows an increased signal intensity in T2-weighted CMR compared to non-ischemic myocardium. Several techniques have been proposed for quantitative assessment of MaR in T2-weighted CMR, such as manual delineation , and threshold methods of two standard deviations (2SD) from remote [5, 8], full width half maximum (FWHM) intensity  and Otsu . Human observers delineating MaR take into account both regional intensity differences and a priori knowledge on perfusion territories and transmurality, which may improve accuracy of MaR quantification. However, manual delineation is subjective and time consuming. Semi-automatic and automatic threshold methods such as 2SD, FWHM and Otsu have been proposed as more objective methods. A more advanced automatic algorithm for quantification of edema in T2-weighted CMR has recently been developed by Johnstone et al . Their algorithm shows promising result for an automatic segmentation approach of edema and thereby MaR in T2-weighted CMR by incorporating regional analysis. However, neither the threshold methods, (2SD, FWHM and Otsu), nor the algorithm by Johnstone et al. uses a priori knowledge on the appearance of MaR and the cardiac anatomy, which is considered when performing manual delineations.
Therefore, the aim of this study was to develop an automatic segmentation algorithm for quantification of MaR in T2-weighted CMR images which uses a priori knowledge on the appearance of MaR and cardiac anatomy.
The study was approved by the local ethics committee and all patients gave their written informed consent. Forty seven patients (age 60.3 ± 9.8 years, range 39 - 83, 39 males) with first-time acute ST-elevation myocardial infarction (STEMI) due to a single occluded coronary artery confirmed by angiography were prospectively included in the study. All patients were treated with primary percutaneous coronary intervention (PCI) with coronary stenting, resulting in TIMI grade 3 flow in the culprit artery.
Within a week after admission patients were imaged in the supine position using either a 1.5 T system (Magnetom Vision, Siemens, Erlangen, Germany) with a CP body array coil or a 1.5 T system (Philips Intera CV or Achieva, Philips, Best, the Netherlands) with a cardiac synergy coil. Initial scout images were acquired to locate the heart, and a T2-weighted triple inversion turbo spin echo sequence (STIR) was employed to depict the myocardium at risk. T2-weighted CMR images were acquired in the short-axis view, covering the left ventricle from the base to apex. Imaging parameters were: echo time 43 ms (Siemens) or 100 ms (Philips); repetition time 2 heart beats; number of averages 2; inversion time 180 ms; typical image resolution 1.5 × 1.5 mm (Siemens) or 1.4 × 1.7 mm reconstructed to 0.7 × 0.7 mm (Philips); slice thickness 8 mm with a typical slice gap of 2 mm. When acquiring images with the cardiac synergy coil no parallel imaging was performed (SENSE = 1).
The MaR was manually delineated according to the method previously described by Carlsson et al . In short, endocardial and epicardial borders of the LV were traced in all short-axis slices by an experienced observer and the papillaries were excluded from the myocardium. Regions of hyper enhanced myocardium was manually delineated as myocardium at risk (MaR) by an experienced observer and expressed as percent of left ventricular mass (LVM). Hypo-intense myocardium within the area of increased signal intensity was regarded as microvascular obstruction  and was included in the MaR.
The new segmentation algorithm, called Segment MaR, was implemented in the freely available cardiac image analysis software Segment (http://segment.heiberg.se)  and will be made available at time of publication. Segment was also used for manual delineation and implementation of the threshold methods (2SD, FWHM and Otsu).
Automatic segmentation algorithm, Segment MaR
The automatic segmentation algorithm, Segment MaR, defines the MaR within the manually delineated left ventricular myocardium based on the culprit artery defined by the user. The MaR region is defined as a continuous region which has a higher probability of being MaR compared to normal myocardium, based on the signal intensity, and fulfills a priori criteria for MaR regarding transmurality, shape, size and extent within the perfusion territory of the culprit artery.
User input as culprit artery and orientation
Estimate intensity distribution for normal myocardium and MaR
Define MaR probability from intensity distributions
Calculate MaR probability in 24 sectors for each short axis slice
Find region of connected sectors with high probability of MaR which fulfills regional criteria
Interpolate shape of MaR region over slices by normalized averaging
Define MaR region in the short-axis view from sector based segmentation
In step 1) the user defines the culprit artery as either, LAD, LCx, RCA or LM based on the overall appearance of the hyperenhanced region and indicates the orientation of the heart by indicating the inferior insertion point for the right ventricle (Figure 2:1).
In step 2) the intensity distribution for normal myocardium and MaR is estimated from the intensity histogram by an expectation maximization (EM) algorithm slice by slice . The intensity distributions were analyzed slice by slice since the intensity can vary between slices in T2-weighted CMR images. The EM-algorithm estimates the mean and standard deviation of the intensity distributions for MaR and normal myocardium by refining an initial estimation. The initial estimation of mean and standard deviation for MaR is calculated from the intensities within the normal extent model for the culprit artery and for normal myocardium calculated from the intensities outside the maximal extent model for the culprit artery. The intensity distributions for MaR and normal myocardium are assumed to be Gaussian and are thereby defined by their mean and standard deviation as estimated by the EM-algorithm. Figure 2:2 shows the intensity histogram of the myocardium and the estimated Gaussian intensity distributions for normal myocardium and MaR.
In step 3) a MaR probability is defined from the intensity distributions for normal myocardium and MaR. The MaR probability is defined from the Bayesian probability of MaR given that the intensity is either MaR or normal myocardium by dividing the Gaussian intensity distribution for MaR by the sum of the Gaussian distribution for MaR and normal myocardium. Thus, the probability function is in the range from 0 to 1 and values above 0.5 indicate higher probability of MaR than normal myocardium. A probability function and its corresponding Gaussian distributions are shown in Figure 2:3.
In step 4) the MaR probability is calculated for 24 sectors in each short axis slice. The probability function in step 3) is defined for each short axis slice and mapped to each pixel in the slice. The probability value is then averaged for 24 sectors in each slice resulting in a sector-based MaR probability. The sector-based MaR probability is shown in Figure 2:4 where bright colors indicate high probability of MaR and the red border indicates sectors with a probability value above 0.5.
In step 5) a region of connected sectors with a high probability for MaR which fulfills regional criteria is identified. The criteria to be fulfilled is a) sectors with a probability value above 0.5 should be connected to its nearest neighboring sector within the slice or in an adjacent slice in a 4-neighbourhood to constitute a region, b) sectors should be localized within the maximal extent model for the culprit artery, c) in the slices with outflow tract only sectors on the anterior side of the outflow tract is considered MaR for LAD and LM and only sectors on the inferior side for LCx and RCA. Finally, the MaR probability of each region is calculated by summing the probability value of being MaR for each pixel within the region and summing the probability of being normal i.e. 1 minus the probability of being MaR for each pixel outside the region. If multiple connected regions are found, the region with the highest probability is chosen and the other regions are eliminated from the MaR region. The outer boundary of the new MaR region is indicated with a red border in Figure 2:5.
In step 6), the shape of the MaR region is refined by interpolating the outer boundary of the MaR region over slices by normalized averaging .The normalized averaging interpolates the outer boundary by using certainty values for the outer boundary of each slice and a narrow kernel with width of 3 slices in both apical and basal direction. The certainty value is lowered if a) the region is close to the maximal extent model, b) if the difference in extent deviates from normal difference between slices and c) if the intensity appearance does not match the boundary of the MaR region. The certainty based on closeness to the maximal extent, a), is calculated by a linear function from one to zero between the normal and maximal extent model in each slice. The certainty based on difference in extent between slices, b), is calculated from a Gaussian function with standard deviation of two sectors and a mean of increasing two sectors from base to apex for LAD and LM respectively a mean of decreasing one sector for RCA and LCx. The certainty based on intensity appearance, c), was calculated as the MaR probability, as defined in step 5, for each slice. The new boundary defined by the normalized averaging is interpolated over slices to give a smooth appearance of the MaR region (Figure 2:6).
In step 7) the MaR region is defined in a short axis view from the sector based MaR region. This is done by defining all pixels within a sector as MaR if the sector is within the outer boundary of the sector based MaR region (Figure 2:7).
Comparison to other segmentation methods
The new automatic segmentation method, Segment MaR, was compared to inter observer variability of manual delineation and three threshold methods, 2 standard deviations from remote (2SD), full width half maximum intensity (FWHM) and Otsu. All methods used the same manual delineation of endocardium and epicardium and were applied slice by slice. A slice by slice approach was chosen since the intensity varies between slices in T2-weighted images.
A second experienced and independent observer manually delineated MaR in all subjects for inter observer analysis. The 2SD threshold method estimates an intensity threshold from a remote region as the mean plus two standard deviations of the intensity within the remote region. The remote region was defined as the region outside the maximal extent model of the culprit artery, indicated by the white sectors in the maximal extent model in Figure 1. The FWHM threshold method  estimates an intensity threshold from a remote region as midway between the mean intensity within the remote region and the maximal intensity within the myocardium. The remote region was defined in the same way as for 2SD. The threshold method of Otsu  estimates the intensity threshold from the histogram of all intensities to get minimal variance both above and below the threshold. The intensity threshold was calculated and applied slice by slice for all three threshold methods.
The quantification of MaR by the automatic Segment MaR algorithm, manual second observer delineation, the threshold methods of 2SD, FWHM and Otsu were all compared against the reference observer using Bland-Altman bias (mean ± standard deviation), paired t-test and linear regression analysis (correlation coefficient and p-value). Regional agreement to manual delineation by the reference observer was evaluated by calculating the Dice similarity coefficient (DSC) , which can be derived from the kappa statistics for the classification of pixels . The DSC is calculated as two times the volume of the intersection of the MaR regions divided by the sum of the volumes of the two MaR regions. The DSC is therefore 0 if the regions do not overlap and 1 if the regions overlap perfectly. The DSC was calculated against the reference observer, for Segment MaR, second observer delineation, 2SD, FWHM and Otsu, for each patient and expressed as mean ± standard deviation.
Results for all five segmentation methods compared to reference delineation
[% of LVM]
-1.9 ± 6.4
0.85 ± 0.07
Second observer delineation
-2.3 ± 4.9
0.90 ± 0.08
-7.7 ± 11.4
0.69 ± 0.14
-21.0 ± 9.9
0.46 ± 0.14
5.3 ± 9.6
0.68 ± 0.10
In eight of the forty-seven patients multiple hyperenhanced regions were detected in step 5 of the algorithm. In seven out of those the same region as by manual delineation was identified. In the one case, a dark artifact divided the MaR region into two disconnected regions, resulting in the Segment MaR algorithm only identifying one of these and subsequently underestimating the MaR region.
This study has presented an automatic segmentation algorithm for quantification of MaR from T2-weighted CMR, based on the EM-algorithm and a priori information on normal and maximal perfusion territories for the culprit artery. Compared to manual delineation by a reference observer, the new algorithm, Segment MaR, performed better than previously suggested threshold methods (2SD, FWHM and Otsu). Quantitative bias and regional agreement for Segment MaR were similar to inter observer variability of manual delineation.
The new automatic segmentation algorithm, Segment MaR, estimates an intensity based probability of MaR from all intensity information in each short axis slice by an EM-algorithm. The use of all intensity information may make the estimate more robust to noise, artifacts and variation in signal homogeneity. The use of a priori information in the EM-algorithm showed an added value compared to using the EM-algorithm without a priori information in the bias analysis (Figure 5). The constraint of the extent by a maximal extent model eliminates artifacts located outside the perfusion territory and together with the use of only one region of connected sectors step 5 showed an added value in the bias analysis. Small non-transmural artifacts within the perfusion territory of the culprit artery may also be less likely to be considered as MaR by requiring transmural regions from the connected sectors. The bias analysis showed no added value by use of the interpolation although it is needed for a physiological appearance of the MaR region.
In the edema algorithm by Johnstone et al.  an EM-algorithm was used to estimate the intensity distributions and a threshold was defined as 2 standard deviations above the mean intensity of normal myocardium. Thereby the threshold is similar to that of 2SD threshold with the exception that the intensity of normal myocardium is defined from the EM-algorithm instead of a remote region. Using a threshold of 2SD does not utilize the intensity information on MaR which has been estimated by the EM-algorithm and thereby the variation in signal intensity may not be taken into consideration. Another difference between Segment MaR and the edema algorithm by Johnstone is that the Segment MaR algorithm is based on 24 sectors in each slice instead of pixel wise segmentation. The pixel-wise segmentation does not consider transmurality and may give better precision on the boundary of MaR. Pixel-wise segmentation may, however, be more sensitive to artifacts. Quantitative results for the edema algorithm, reported by Johnstone, were a bias of 1.1 ± 10.1% of LVM and DSC of 0.50 ± 0.27 to their reference of manual delineation and an aim for future work was set to reach DSC > 0.7 which has been stated to indicate excellent regional agreement by Zijdenbos et al. . In the present study, Segment MaR showed a DSC value of 0.85 ± 0.07 and bias of -1.9 ± 6.4% of LVM compared to the manual delineation of the reference observer. Note-able is also the difference in DSC for inter observer variability of manual delineation which was 0.72 ± 0.14 in the study by Johnstone compared to 0.90 ± 0.08 in the present study. This may indicate that the regional agreement may be lower and the quantitative variability of manual delineation may be larger between observers than reported in the present study and thereby the use of an automatic algorithm such as Segment MaR, which has a low bias and high regional agreement to manual delineation, may decrease the inter observer variability.
The weak correlation to manual delineation by the reference observer for the threshold methods of 2SD, FWHM and Otsu may be explained by the fact that the methods are solely based on a threshold and by weaknesses in defining the threshold for the different methods. Using a fixed number of standard deviations as in 2SD does not account for the variability of intensity for MaR and more importantly the contrast in T2-weighted CMR is lower than for other CMR . The threshold by FWHM may be sensitive to artifacts since the brightest pixel intensity is used to find the threshold. Both 2SD and FWHM are sensitive to the definition of the remote region which is currently not standardized. The remote region in this study was automatically defined as the myocardium outside the maximal extent model of the culprit artery. This definition of the remote myocardium was chosen to obtain an objective and standardized representation of normal myocardium. This strategy may, however, result in artifacts being included in the remote myocardium. An overestimation of the threshold may for 2SD be caused by the large standard deviation of the intensity within the remote region. The border zone between normal myocardium and remote myocardium influences the remote myocardium proportionally more for larger MaR regions which may result in a overestimation of the threshold and subsequent underestimation of the MaR region and thus explain the trend seen in bias for 2SD (Figure 3). For FWHM an overestimation of the threshold may be caused by artifacts within the myocardium and the remote region. Bright artifacts may result in a threshold which only identifies the artifact as the MaR region and this can explain the trend in bias seen for FWHM since it will result in a larger underestimation for larger MaR regions. The threshold defined by Otsu does not depend on any remote region but may instead be unstable in the definition of the threshold since it assumes that an optimal threshold should be found in each short axis slice. This implies that both MaR and normal myocardium should be present in each slice, which is not the case in most patients as for example in the basal slices of an LAD occlusion or apical slices of an RCA occlusion. The Otsu threshold may thereby overestimate the MaR region in slices lacking MaR and underestimate MaR in slices lacking normal myocardium. This may explain the large overestimation for small MaR regions and large underestimation for large MaR regions. The EM-algorithm used in Segment MaR also assumes two intensity distributions as in the threshold method of Otsu but is accompanied by a priori information both as initialization to the EM-algorithm and in the post processing of finding connected sectors.
Manual monitoring and possibly manual corrections are as for all automatic segmentation algorithms needed for research and clinical use. The use of Segment MaR may, however, decrease the degree of variability introduced by the subjectivity of manual delineation since the Segment MaR algorithm showed a low bias and high correlation to manual delineation regarding quantitative assessment and an excellent regional agreement according to DSC.
One limitation in this study is the lack of ground truth for in vivo quantification of MaR. Manual delineation according to the methodology used by Carlsson et al.  when validating T2-weighted CMR for MaR to SPECT was chosen as reference method and in this study there was a good inter observer agreement. Due to the limited number of patients in the study it was not possible to use a separate training and test set and the parameters in the automatic Segment MaR algorithm could not be optimized. The Segment MaR algorithm has not been specifically designed for the imaging systems of Philips and Siemens. Further research is suggested to investigate the performance of the algorithm in a larger cohort of patients and possibly optimize and improve the algorithm for specific imaging systems.
In this study, an automatic segmentation algorithm, called Segment MaR, for quantification of myocardium at risk (MaR) in T2-weighted CMR has been presented and showed to have a good agreement to manual delineation. Both the quantitative and regional agreement to manual delineation was better for Segment MaR than for the threshold methods of two standard deviations from remote, full width half maximum intensity and Otsu. The Segment MaR algorithm seems to be a promising, objective method for standardized measurement of MaR in T2-weighted CMR.
The authors would like to acknowledge Ann-Helen Arvidsson and Christel Carlander, both with the Lund Cardiac MR Group, for their skillful assistance with image acquisition. This study has been funded by the Swedish Research Council (2008-2461, 2008-2949), The Swedish Heart and Lung Foundation, The Medical Faculty of Lund University, Sweden, and Region of Scania, Sweden.
- Aletras AH, Tilak GS, Natanzon A, Hsu LY, Gonzalez FM, Hoyt RF, Arai AE: 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-1870.View ArticlePubMedGoogle Scholar
- Fishbein MC, Maclean D, Maroko PR: The histopathologic evolution of myocardial infarction. Chest. 1978, 73: 843-849.View ArticlePubMedGoogle Scholar
- Jennings RB, Schaper J, Hill ML, Steenbergen C, Reimer KA: Effect of reperfusion late in the phase of reversible ischemic injury. Changes in cell volume, electrolytes, metabolites, and ultrastructure. Circ Res. 1985, 56: 262-278.View ArticlePubMedGoogle Scholar
- Carlsson M, Ubachs JF, Hedstrom E, Heiberg E, Jovinge S, Arheden H: Myocardium at risk after acute infarction in humans on cardiac magnetic resonance: quantitative assessment during follow-up and validation with single-photon emission computed tomography. JACC Cardiovasc Imaging. 2009, 2: 569-576.View ArticlePubMedGoogle 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-1587.View ArticlePubMedGoogle Scholar
- Stork A, Muellerleile K, Bansmann PM, Graessner J, Kaul M, Kemper J, Adam G, Lund GK: Value of T2-weighted, first-pass and delayed enhancement, and cine CMR to differentiate between acute and chronic myocardial infarction. Eur Radiol. 2007, 17: 610-617.View ArticlePubMedGoogle Scholar
- Abdel-Aty H, Zagrosek A, Schulz-Menger J, Taylor AJ, Messroghli D, Kumar A, Gross M, Dietz R, Friedrich MG: Delayed enhancement and T2-weighted cardiovascular magnetic resonance imaging differentiate acute from chronic myocardial infarction. Circulation. 2004, 109: 2411-2416.View ArticlePubMedGoogle Scholar
- Wright J, Adriaenssens T, Dymarkowski S, Desmet W, Bogaert J: Quantification of myocardial area at risk with T2-weighted CMR: comparison with contrast-enhanced CMR and coronary angiography. JACC Cardiovasc Imaging. 2009, 2: 825-831.View ArticlePubMedGoogle Scholar
- Tilak GS, Hsu LY, Hoyt RF, Arai AE, Aletras AH: In vivo T2-weighted magnetic resonance imaging can accurately determine the ischemic area at risk for 2-day-old nonreperfused myocardial infarction. Invest Radiol. 2008, 43: 7-15.View ArticlePubMedGoogle Scholar
- Burchell T, Flett A, Petersen S, Davies L, Mohiddin S, Mathur A, Westwood M: Comparing analysis methods for quantification of myocardial oedema in patients following reperfused ST-elevation MI. Journal of Cardiovascular Magnetic Resonance. 2011, 13: M11-PubMed CentralView ArticleGoogle 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. 2011Google Scholar
- Beek AM, Nijveldt R, van Rossum AC: Intramyocardial hemorrhage and microvascular obstruction after primary percutaneous coronary intervention. Int J Cardiovasc Imaging. 2010, 26: 49-55.PubMed CentralView ArticlePubMedGoogle Scholar
- Heiberg E, Sjogren J, Ugander M, Carlsson M, Engblom H, Arheden H: Design and validation of Segment--freely available software for cardiovascular image analysis. BMC Med Imaging. 2010, 10: 1-PubMed CentralView ArticlePubMedGoogle Scholar
- Cerqueira MD, Weissman NJ, Dilsizian V, Jacobs AK, Kaul S, Laskey WK, Pennell DJ, Rumberger JA, Ryan T, Verani MS: Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart: a statement for healthcare professionals from the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association. Circulation. 2002, 105: 539-542.View ArticlePubMedGoogle Scholar
- Dempster AP, Laird NM, Rubin DB: Maximum Likelihood from Incomplete Data Via Em Algorithm. J Roy Stat Soc B Met. 1977, 39: 1-38.Google Scholar
- Granlund G, Knutsson H: Signal processing for computer vision. 1995, Linkoping: Kluwer Academic PublishersView ArticleGoogle 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
- Otsu N: Threshold Selection Method from Gray-Level Histograms. Ieee T Syst Man Cyb. 1979, 9: 62-66.View ArticleGoogle Scholar
- Dice LR: Measures of the Amount of Ecologic Association between Species. Ecology. 1945, 26: 297-302.View ArticleGoogle Scholar
- Zijdenbos AP, Dawant BM, Margolin RA, Palmer AC: Morphometric Analysis of White-Matter Lesions in Mr-Images-Method and Validation. Ieee T Med Imaging. 1994, 13: 716-724.View ArticleGoogle Scholar
- Abdel-Aty H, Simonetti O, Friedrich MG: T2-weighted cardiovascular magnetic resonance imaging. J Magn Reson Imaging. 2007, 26: 452-459.View ArticlePubMedGoogle Scholar
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.