Automated segmentation of myocardial scar in late enhancement MRI using combined intensity and spatial information

Multisize Electrode Field-of-View: Validation by High Resolution Gadolinium-Enhanced Cardiac Magnetic Resonance

BACKGROUND

Voltage mapping to detect ventricular scar is important for guiding catheter ablation, but the field-of-view of unipolar, bipolar, conventional, and microelectrodes as it relates to the extent of viable myocardium (VM) is not well defined.

OBJECTIVES

The purpose of this study was to evaluate electroanatomic voltage-mapping (EAVM) with different-size electrodes for identifying VM, validated against high-resolution ex-vivo cardiac magnetic resonance (HR-LGE-CMR).

METHODS

A total of 9 swine with early-reperfusion myocardial infarction were mapped with the QDOT microcatheter. HR-LGE-CMR (0.3-mm slices) were merged with EAVM. At each EAVM point, the underlying VM in multisize transmural cylinders and spheres was quantified from ex vivo CMR and related to unipolar and bipolar voltages recorded from conventional and microelectrodes.

RESULTS

In each swine, 220 mapping points (Q1, Q3: 216, 260 mapping points) were collected. Infarcts were heterogeneous and nontransmural. Unipolar and bipolar voltage increased with VM volumes from >175 mm3 up to >525 mm3 (equivalent to a 5-mm radius cylinder with height >6.69 mm). VM volumes in subendocardial cylinders with 1- or 3-mm depth correlated poorly with all voltages. Unipolar voltages recorded with conventional and microelectrodes were similar (difference 0.17 ± 2.66 mV) and correlated best to VM within a sphere of radius 10 and 8 mm, respectively. Distance-weighting did not improve the correlation.

CONCLUSIONS

Voltage increases with transmural volume of VM but correlates poorly with small amounts of VM, which limits EAVM in defining heterogeneous scar. Microelectrodes cannot distinguish thin from thick areas of subendocardial VM. The field-of-view for unipolar recordings for microelectrodes and conventional electrodes appears to be 8 to 10 mm, respectively, and unexpectedly similar.

Deep learning for classification of late gadolinium enhancement lesions based on the 16-segment left ventricular model

PURPOSE

This study aimed to develop and validate a deep learning-based method that allows for segmental analysis of myocardial late gadolinium enhancement (LGE) lesions.

METHODS

Cardiac LGE data from 170 patients with coronary artery disease and non-ischemic heart disease were used for training, validation, and testing. Short-axis images were transformed to polar space after identification of the left ventricular (LV) center point and anterior right ventricular (RV) insertion point. Images were obtained after dividing the polar transformed images into segments based on the 16-segment LV model. Five different deep convolutional neural network (CNN) models were developed and validated using the labeled data, where the image after the division corresponded to a segment, and the lesion labeling was based on the 16-segment LV model. Unseen testing data were used to evaluate the performance of the lesion classification.

RESULTS

Without manual lesion segmentation and annotation, the proposed method showed an area under the curve (AUC) of 0.875, and a precision, recall, and F1-score of 0.723, 0.783, and 0.752, respectively for the lesion class when the pretrained ResNet50 model was tested for all slice images. The two pretrained models of ResNet50 and EfficientNet-B0 outperformed the three non-pretrained CNN models in terms of AUCs (0.873-0.875 vs. 0.834-0.841).

CONCLUSION

The proposed method is based on learning a deep CNN model from polar transformed images to predict LGE lesions with good accuracy and does not require time-consuming annotation procedures such as lesion segmentation.

MyoPS: A benchmark of myocardial pathology segmentation combining three-sequence cardiac magnetic resonance images

Assessment of myocardial viability is essential in diagnosis and treatment management of patients suffering from myocardial infarction, and classification of pathology on the myocardium is the key to this assessment. This work defines a new task of medical image analysis, i.e., to perform myocardial pathology segmentation (MyoPS) combining three-sequence cardiac magnetic resonance (CMR) images, which was first proposed in the MyoPS challenge, in conjunction with MICCAI 2020. Note that MyoPS refers to both myocardial pathology segmentation and the challenge in this paper. The challenge provided 45 paired and pre-aligned CMR images, allowing algorithms to combine the complementary information from the three CMR sequences for pathology segmentation. In this article, we provide details of the challenge, survey the works from fifteen participants and interpret their methods according to five aspects, i.e., preprocessing, data augmentation, learning strategy, model architecture and post-processing. In addition, we analyze the results with respect to different factors, in order to examine the key obstacles and explore the potential of solutions, as well as to provide a benchmark for future research. The average Dice scores of submitted algorithms were 0.614±0.231 and 0.644±0.153 for myocardial scars and edema, respectively. We conclude that while promising results have been reported, the research is still in the early stage, and more in-depth exploration is needed before a successful application to the clinics. MyoPS data and evaluation tool continue to be publicly available upon registration via its homepage (www.sdspeople.fudan.edu.cn/zhuangxiahai/0/myops20/).

Artificial Intelligence as a Diagnostic Tool in Non-Invasive Imaging in the Assessment of Coronary Artery Disease

Coronary artery disease (CAD) remains a leading cause of mortality and morbidity worldwide, and it is associated with considerable economic burden. In an ageing, multimorbid population, it has become increasingly important to develop reliable, consistent, low-risk, non-invasive means of diagnosing CAD. The evolution of multiple cardiac modalities in this field has addressed this dilemma to a large extent, not only in providing information regarding anatomical disease, as is the case with coronary computed tomography angiography (CCTA), but also in contributing critical details about functional assessment, for instance, using stress cardiac magnetic resonance (S-CMR). The field of artificial intelligence (AI) is developing at an astounding pace, especially in healthcare. In healthcare, key milestones have been achieved using AI and machine learning (ML) in various clinical settings, from smartwatches detecting arrhythmias to retinal image analysis and skin cancer prediction. In recent times, we have seen an emerging interest in developing AI-based technology in the field of cardiovascular imaging, as it is felt that ML methods have potential to overcome some limitations of current risk models by applying computer algorithms to large databases with multidimensional variables, thus enabling the inclusion of complex relationships to predict outcomes. In this paper, we review the current literature on the various applications of AI in the assessment of CAD, with a focus on multimodality imaging, followed by a discussion on future perspectives and critical challenges that this field is likely to encounter as it continues to evolve in cardiology.

Automatic development of 3D anatomical models of border zone and core scar regions in the left ventricle

Patients with myocardial infarction are at elevated risk of sudden cardiac death, and scar tissue arising from infarction is known to play a role. The accurate identification of scars therefore is crucial for risk assessment, quantification and guiding interventions. Typically, core scars and grey peripheral zones are identified by radiologists and clinicians based on cardiac late gadolinium enhancement magnetic resonance images (LGE-MRI). Scar regions from LGE-MRI vary in size, shape, heterogeneity, artifacts, and image resolution. Thus, manual segmentation is time consuming, and influenced by the observer's experience (bias effect). We propose a fully automatic framework that develops 3D anatomical models of the left ventricle with border zone and core scar regions that are free from bias effect. Our myocardium (SOCRATIS), border scar and core scar (BZ-SOCRATIS) segmentation pipelines were evaluated using internal and external validation datasets. The automatic myocardium segmentation framework performed a Dice score of 81.9% and 70.0% in the internal and external validation dataset. The automatic scar segmentation pipeline achieved a Dice score of 60.9% for the core scar segmentation and 43.7% for the border zone scar segmentation in the internal dataset and in the external dataset a Dice score of 44.2% for the core scar segmentation and 54.8% for the border scar segmentation respectively. To the best of our knowledge, this is the first study outlining a fully automatic framework to develop 3D anatomical models of the left ventricle with border zone and core scar regions. Our method exhibits high performance without the need for training or tuning in an unseen cohort (unsupervised).

Optimized automated cardiac MR scar quantification with GAN-based data augmentation

BACKGROUND

The clinical utility of late gadolinium enhancement (LGE) cardiac MRI is limited by the lack of standardization, and time-consuming postprocessing. In this work, we tested the hypothesis that a cascaded deep learning pipeline trained with augmentation by synthetically generated data would improve model accuracy and robustness for automated scar quantification.

METHODS

A cascaded pipeline consisting of three consecutive neural networks is proposed, starting with a bounding box regression network to identify a region of interest around the left ventricular (LV) myocardium. Two further nnU-Net models are then used to segment the myocardium and, if present, scar. The models were trained on the data from the EMIDEC challenge, supplemented with an extensive synthetic dataset generated with a conditional GAN.

RESULTS

The cascaded pipeline significantly outperformed a single nnU-Net directly segmenting both the myocardium (mean Dice similarity coefficient (DSC) (standard deviation (SD)): 0.84 (0.09) vs 0.63 (0.20), p < 0.01) and scar (DSC: 0.72 (0.34) vs 0.46 (0.39), p < 0.01) on a per-slice level. The inclusion of the synthetic data as data augmentation during training improved the scar segmentation DSC by 0.06 (p < 0.01). The mean DSC per-subject on the challenge test set, for the cascaded pipeline augmented by synthetic generated data, was 0.86 (0.03) and 0.67 (0.29) for myocardium and scar, respectively.

CONCLUSION

A cascaded deep learning-based pipeline trained with augmentation by synthetically generated data leads to myocardium and scar segmentations that are similar to the manual operator, and outperforms direct segmentation without the synthetic images.

MA-SOCRATIS: An automatic pipeline for robust segmentation of the left ventricle and scar

Multi-atlas segmentation of cardiac regions and total infarct scar (MA-SOCRATIS) is an unsupervised automatic pipeline to segment left ventricular myocardium and scar from late gadolinium enhanced MR images (LGE-MRI) of the heart. We implement two different pipelines for myocardial and scar segmentation from short axis LGE-MRI. Myocardial segmentation has two steps; initial segmentation and re-estimation. The initial segmentation step makes a first estimate of myocardium boundaries by using multi-atlas segmentation techniques. The re-estimation step refines the myocardial segmentation by a combination of k-means clustering and a geometric median shape variation technique. An active contour technique determines the unhealthy and healthy myocardial wall. The scar segmentation pipeline is a combination of a Rician-Gaussian mixture model and full width at half maximum (FWHM) thresholding, to determine the intensity pixels in scar regions. Following this step a watershed method with an automatic seed-points framework segments the final scar region. MA-SOCRATIS was evaluated using two different datasets. In both datasets ground truths were based on manual segmentation of short axis images from LGE-MRI scans. The first dataset included 40 patients from the MS-CMRSeg 2019 challenge dataset (STACOM at MICCAI 2019). The second is a collection of 20 patients with scar regions that are challenging to segment. MA-SOCRATIS achieved robust and accurate performance in automatic segmentation of myocardium and scar regions without the need of training or tuning in both cohorts, compared with state-of-the-art techniques (intra-observer and inter observer myocardium segmentation: 81.9% and 70% average Dice value, and scar (intra-observer and inter observer segmentation: 70.5% and 70.5% average Dice value).

Recent Advances in Fibrosis and Scar Segmentation From Cardiac MRI: A State-of-the-Art Review and Future Perspectives

Segmentation of cardiac fibrosis and scars is essential for clinical diagnosis and can provide invaluable guidance for the treatment of cardiac diseases. Late Gadolinium enhancement (LGE) cardiovascular magnetic resonance (CMR) has been successful in guiding the clinical diagnosis and treatment reliably. For LGE CMR, many methods have demonstrated success in accurately segmenting scarring regions. Co-registration with other non-contrast-agent (non-CA) modalities [e.g., balanced steady-state free precession (bSSFP) cine magnetic resonance imaging (MRI)] can further enhance the efficacy of automated segmentation of cardiac anatomies. Many conventional methods have been proposed to provide automated or semi-automated segmentation of scars. With the development of deep learning in recent years, we can also see more advanced methods that are more efficient in providing more accurate segmentations. This paper conducts a state-of-the-art review of conventional and current state-of-the-art approaches utilizing different modalities for accurate cardiac fibrosis and scar segmentation.

Myocardial Infarction Quantification from Late Gadolinium Enhancement MRI Using Top-Hat Transforms and Neural Networks

Late gadolinium enhancement (LGE) MRI is the gold standard technique for myocardial viability assessment. Although the technique accurately reflects the damaged tissue, there is no clinical standard to quantify myocardial infarction (MI). Moreover, commercial software used in clinical practice are mostly semi-automatic, and hence require direct intervention of experts. In this work, a new automatic method for MI quantification from LGE-MRI is proposed. Our novel segmentation approach is devised for accurately detecting not only hyper-enhanced lesions, but also microvascular obstruction areas. Moreover, it includes a myocardial disease detection step which extends the algorithm for working under healthy scans. The method is based on a cascade approach where firstly, diseased slices are identified by a convolutional neural network (CNN). Secondly, by means of morphological operations a fast coarse scar segmentation is obtained. Thirdly, the segmentation is refined by a boundary-voxel reclassification strategy using an ensemble of very light CNNs. We tested the method on a LGE-MRI database with healthy (n = 20) and diseased (n = 80) cases following a 5-fold cross-validation scheme. Our approach segmented myocardial scars with an average Dice coefficient of 77.22 ± 14.3% and with a volumetric error of 1.0 ± 6.9 cm3. In a comparison against nine reference algorithms, the proposed method achieved the highest agreement in volumetric scar quantification with the expert delineations (p< 0.001 when compared to the other approaches). Moreover, it was able to reproduce the scar segmentation intra- and inter-rater variability. Our approach was shown to be a good first attempt towards automatic and accurate myocardial scar segmentation, although validation over larger LGE-MRI databases is needed.

Reproducibility of clinical late gadolinium enhancement magnetic resonance imaging in detecting left atrial scar after atrial fibrillation ablation

INTRODUCTION

Late gadolinium enhancement (LGE) cardiac magnetic resonance imaging (MRI) can be used to detect postablation atrial scar (PAAS) but its reproducibility and reliability in clinical scans across different magnetic flux densities and scar detection methods are unknown.

METHODS

Patients (n = 45) having undergone two consecutive MRIs (3 months apart) on 3T and 1.5T scanners were studied. We compared PAAS detection reproducibility using four methods of thresholding: simple thresholding, Otsu thresholding, 3.3 standard deviations (SD) above blood pool (BP) mean intensity, and image intensity ratio (IIR). We performed a texture study by dividing the left atrial wall intensity histogram into deciles and evaluated the correlation of the same decile of the two scans as well as to a randomized distribution of intensities, quantified using Dice Similarity Coefficient (DSC).

RESULTS

The choice of scanner did not significantly affect the reproducibility. The scar detection performed by Otsu thresholding (DSC of 71.26 ± 8.34) resulted in a better correlation of the two scans compared with the methods of 3.3 SD above BP mean intensity (DSC of 57.78 ± 21.2, p < .001) and IIR above 1.61 (DSC of 45.76 ± 29.55, p <.001). Texture analysis showed that correlation only for voxels with intensities in deciles above the 70th percentile of wall intensity histogram was better than random distribution (p < .001).

CONCLUSIONS

Our results demonstrate that clinical LGE-MRI can be reliably used for visualizing PAAS across different magnetic flux densities if the threshold is greater than 70th percentile of the wall intensity distribution. Also, atrial wall-based thresholding is better than BP-based thresholding for reproducible PAAS detection.

Three-dimensional Deep Convolutional Neural Networks for Automated Myocardial Scar Quantification in Hypertrophic Cardiomyopathy: A Multicenter Multivendor Study

Background Cardiac MRI late gadolinium enhancement (LGE) scar volume is an important marker for outcome prediction in patients with hypertrophic cardiomyopathy (HCM); however, its clinical application is hindered by a lack of measurement standardization. Purpose To develop and evaluate a three-dimensional (3D) convolutional neural network (CNN)-based method for automated LGE scar quantification in patients with HCM. Materials and Methods We retrospectively identified LGE MRI data in a multicenter (n = 7) and multivendor (n = 3) HCM study obtained between November 2001 and November 2011. A deep 3D CNN based on U-Net architecture was used for LGE scar quantification. Independent CNN training and testing data sets were maintained with a 4:1 ratio. Stacks of short-axis MRI slices were split into overlapping substacks that were segmented and then merged into one volume. The 3D CNN per-site and per-vendor performances were evaluated with respect to manual scar quantification performed in a core laboratory setting using Dice similarity coefficient (DSC), Pearson correlation, and Bland-Altman analyses. Furthermore, the performance of 3D CNN was compared with that of two-dimensional (2D) CNN. Results This study included 1073 patients with HCM (733 men; mean age, 49 years ± 17 [standard deviation]). The 3D CNN-based quantification was fast (0.15 second per image) and demonstrated excellent correlation with manual scar volume quantification (r = 0.88, P < .001) and ratio of scar volume to total left ventricle myocardial volume (%LGE) (r = 0.91, P < .001). The 3D CNN-based quantification strongly correlated with manual quantification of scar volume (r = 0.82-0.99, P < .001) and %LGE (r = 0.90-0.97, P < .001) for all sites and vendors. The 3D CNN identified patients with a large scar burden (>15%) with 98% accuracy (202 of 207) (95% confidence interval [CI]: 95%, 99%). When compared with 3D CNN, 2D CNN underestimated scar volume (r = 0.85, P < .001) and %LGE (r = 0.83, P < .001). The DSC of 3D CNN segmentation was comparable among different vendors (P = .07) and higher than that of 2D CNN (DSC, 0.54 ± 0.26 vs 0.48 ± 0.29; P = .02). Conclusion In the hypertrophic cardiomyopathy population, a three-dimensional convolutional neural network enables fast and accurate quantification of myocardial scar volume, outperforms a two-dimensional convolutional neural network, and demonstrates comparable performance across different vendors. © RSNA, 2019 Online supplemental material is available for this article.

Myocardial Radiomics in Cardiac MRI

OBJECTIVE. The purpose of this article is to review the nascent field of radiomics in cardiac MRI. CONCLUSION. Cardiac MRI produces a large number of images in a fairly inefficient manner with sometimes limited clinical application. In the era of precision medicine, there is increasing need for imaging to account for a broader array of diseases in an efficient and objective manner. Radiomics, the extraction and analysis of quantitative imaging features from medical imaging, may offer potential solutions to this need.

Algorithms for left atrial wall segmentation and thickness - Evaluation on an open-source CT and MRI image database

Structural changes to the wall of the left atrium are known to occur with conditions that predispose to Atrial fibrillation. Imaging studies have demonstrated that these changes may be detected non-invasively. An important indicator of this structural change is the wall's thickness. Present studies have commonly measured the wall thickness at few discrete locations. Dense measurements with computer algorithms may be possible on cardiac scans of Computed Tomography (CT) and Magnetic Resonance Imaging (MRI). The task is challenging as the atrial wall is a thin tissue and the imaging resolution is a limiting factor. It is unclear how accurate algorithms may get and how they compare in this new emerging area. We approached this problem of comparability with the Segmentation of Left Atrial Wall for Thickness (SLAWT) challenge organised in conjunction with MICCAI 2016 conference. This manuscript presents the algorithms that had participated and evaluation strategies for comparing them on the challenge image database that is now open-source. The image database consisted of cardiac CT (n=10) and MRI (n=10) of healthy and diseased subjects. A total of 6 algorithms were evaluated with different metrics, with 3 algorithms in each modality. Segmentation of the wall with algorithms was found to be feasible in both modalities. There was generally a lack of accuracy in the algorithms and inter-rater differences showed that algorithms could do better. Benchmarks were determined and algorithms were ranked to allow future algorithms to be ranked alongside the state-of-the-art techniques presented in this work. A mean atlas was also constructed from both modalities to illustrate the variation in thickness within this small cohort.

Microstructural Infarct Border Zone Remodeling in the Post-infarct Swine Heart Measured by Diffusion Tensor MRI

Introduction: Computational models of the heart increasingly require detailed microstructural information to capture the impact of tissue remodeling on cardiac electromechanics in, for example, hearts with myocardial infarctions. Myocardial infarctions are surrounded by the infarct border zone (BZ), which is a site of electromechanical property transition. Magnetic resonance imaging (MRI) is an emerging method for characterizing microstructural remodeling and focal myocardial infarcts and the BZ can be identified with late gadolinium enhanced (LGE) MRI. Microstructural remodeling within the BZ, however, remains poorly characterized by MRI due, in part, to the fact that LGE and DT-MRI are not always available for the same heart. Diffusion tensor MRI (DT-MRI) can evaluate microstructural remodeling by quantifying the DT apparent diffusion coefficient (ADC, increased with decreased cellularity), fractional anisotropy (FA, decreased with increased fibrosis), and tissue mode (decreased with increased fiber disarray). The purpose of this work was to use LGE MRI in post-infarct porcine hearts (N = 7) to segment remote, BZ, and infarcted myocardium, thereby providing a basis to quantify microstructural remodeling in the BZ and infarcted regions using co-registered DT-MRI. Methods: Chronic porcine infarcts were created by balloon occlusion of the LCx. 6-8 weeks post-infarction, MRI contrast was administered, and the heart was potassium arrested, excised, and imaged with LGE MRI (0.33 × 0.33 × 0.33 mm) and co-registered DT-MRI (1 × 1 × 3 mm). Myocardium was segmented as remote, BZ, or infarct by LGE signal intensity thresholds. DT invariants were used to evaluate microstructural remodeling by quantifying ADC, FA, and tissue mode. Results: The BZ significantly remodeled compared to both infarct and remote myocardium. BZ demonstrated a significant decrease in cellularity (increased ADC), significant decrease in tissue organization (decreased FA), and a significant increase in fiber disarray (decreased tissue mode) relative to remote myocardium (all p < 0.05). Microstructural remodeling in the infarct was similar, but significantly larger in magnitude (all p < 0.05). Conclusion: DT-MRI can identify regions of significant microstructural remodeling in the BZ that are distinct from both remote and infarcted myocardium.

Empirical Mode Decomposition and Monogenic Signal-Based Approach for Quantification of Myocardial Infarction From MR Images

Quantification of myocardial infarction on late Gadolinium enhancement cardiovascular magnetic resonance (LGE-CMR) images into heterogeneous infarct periphery (or gray zone) and infarct core plays an important role in cardiac diagnosis, especially in identifying patients at high risk of cardiovascular mortality. However, quantification task is challenging due to noise corrupted in cardiac MR images, the contrast variation, and limited resolution of images. In this study, we propose a novel approach for automatic myocardial infarction quantification, termed DEMPOT, which consists of three key parts: Decomposition of image into intrinsic modes, monogenic phase performing on combined dominant modes, and multilevel Otsu thresholding on the phase. In particular, inspired by the Hilbert-Huang transform, we perform the multidimensional ensemble empirical mode decomposition and 2-D generalization of the Hilbert transform known as the Riesz transform on the MR image to obtain the monogenic phase that is robust to noise and contrast variation. Then, a two-stage algorithm using multilevel Otsu thresholding is accomplished on the monogenic phase to automatically quantify the myocardium into healthy, gray zone, and infarct core regions. Experiments on LGE-CMR images with myocardial infarction from 82 patients show the superior performance of the proposed approach in terms of reproducibility, robustness, and effectiveness.

Myocardium segmentation from DE MRI with guided random walks and sparse shape representation

PURPOSE

For patients with myocardial infarction (MI), delayed enhancement (DE) cardiovascular magnetic resonance imaging (MRI) is a sensitive and well-validated technique for the detection and visualization of MI. The myocardium viability assessment with DE MRI is important in diagnosis and treatment management, where myocardium segmentation is a prerequisite. However, few academic works have focused on automated myocardium segmentation from DE images. In this study, we aim to develop an automatic myocardium segmentation algorithm that targets DE images.

METHODS

We propose a segmentation framework based on both prior shape knowledge and image intensity. Instead of the strong request of the pre-segmentation of cine MRI in the same session, we use the sparse representation method to model the myocardium shape. Data from the Cardiac MR Left Ventricle Segmentation Challenge (2009) are used to build the shape template repository. The method of guided random walks is used to integrate the shape model and intensity information. An iterative approach is used to gradually improve the results.

RESULTS

The proposed method was tested on the DE MRI data from 30 MI patients. The proposed method achieved Dice similarity coefficients (DSC) of 74.60 ± 7.79% with 201 shape templates and 73.56 ± 6.32% with 56 shape templates, which were close to the inter-observer difference (73.94 ± 5.12%). To test the generalization of the proposed method to routine clinical images, the DE images of 10 successive new patients were collected, which were unseen during the method development and parameter tuning, and a DSC of 76.02 ± 7.43% was achieved.

CONCLUSION

The authors propose a novel approach for the segmentation of myocardium from DE MRI by using the sparse representation-based shape model and guided random walks. The sparse representation method effectively models the prior shape with a small number of shape templates, and the proposed method has the potential to achieve clinically relevant results.

Myocardium Segmentation From DE MRI Using Multicomponent Gaussian Mixture Model and Coupled Level Set

Objective: In this paper, we propose a fully automatic framework for myocardium segmentation of delayed-enhancement (DE) MRI images without relying on prior patient-specific information. Methods: We employ a multicomponent Gaussian mixture model to deal with the intensity heterogeneity of myocardium caused by the infarcts. To differentiate the myocardium from other tissues with similar intensities, while at the same time maintain spatial continuity, we introduce a coupled level set (CLS) to regularize the posterior probability. The CLS, as a spatial regularization, can be adapted to the image characteristics dynamically. We also introduce an image intensity gradient based term into the CLS, adding an extra force to the posterior probability based framework, to improve the accuracy of myocardium boundary delineation. The prebuilt atlases are propagated to the target image to initialize the framework. Results: The proposed method was tested on datasets of 22 clinical cases, and achieved Dice similarity coefficients of 87.43 ± 5.62% (endocardium), 90.53 ± 3.20% (epicardium) and 73.58 ± 5.58% (myocardium), which have outperformed three variants of the classic segmentation methods. Conclusion: The results can provide a benchmark for the myocardial segmentation in the literature. Significance: DE MRI provides an important tool to assess the viability of myocardium. The accurate segmentation of myocardium, which is a prerequisite for further quantitative analysis of myocardial infarction (MI) region, can provide important support for the diagnosis and treatment management for MI patients.Objective: In this paper, we propose a fully automatic framework for myocardium segmentation of delayed-enhancement (DE) MRI images without relying on prior patient-specific information. Methods: We employ a multicomponent Gaussian mixture model to deal with the intensity heterogeneity of myocardium caused by the infarcts. To differentiate the myocardium from other tissues with similar intensities, while at the same time maintain spatial continuity, we introduce a coupled level set (CLS) to regularize the posterior probability. The CLS, as a spatial regularization, can be adapted to the image characteristics dynamically. We also introduce an image intensity gradient based term into the CLS, adding an extra force to the posterior probability based framework, to improve the accuracy of myocardium boundary delineation. The prebuilt atlases are propagated to the target image to initialize the framework. Results: The proposed method was tested on datasets of 22 clinical cases, and achieved Dice similarity coefficients of 87.43 ± 5.62% (endocardium), 90.53 ± 3.20% (epicardium) and 73.58 ± 5.58% (myocardium), which have outperformed three variants of the classic segmentation methods. Conclusion: The results can provide a benchmark for the myocardial segmentation in the literature. Significance: DE MRI provides an important tool to assess the viability of myocardium. The accurate segmentation of myocardium, which is a prerequisite for further quantitative analysis of myocardial infarction (MI) region, can provide important support for the diagnosis and treatment management for MI patients.

Cardiac imaging: working towards fully-automated machine analysis & interpretation

INTRODUCTION

Non-invasive imaging plays a critical role in managing patients with cardiovascular disease. Although subjective visual interpretation remains the clinical mainstay, quantitative analysis facilitates objective, evidence-based management, and advances in clinical research. This has driven developments in computing and software tools aimed at achieving fully automated image processing and quantitative analysis. In parallel, machine learning techniques have been used to rapidly integrate large amounts of clinical and quantitative imaging data to provide highly personalized individual patient-based conclusions. Areas covered: This review summarizes recent advances in automated quantitative imaging in cardiology and describes the latest techniques which incorporate machine learning principles. The review focuses on the cardiac imaging techniques which are in wide clinical use. It also discusses key issues and obstacles for these tools to become utilized in mainstream clinical practice. Expert commentary: Fully-automated processing and high-level computer interpretation of cardiac imaging are becoming a reality. Application of machine learning to the vast amounts of quantitative data generated per scan and integration with clinical data also facilitates a move to more patient-specific interpretation. These developments are unlikely to replace interpreting physicians but will provide them with highly accurate tools to detect disease, risk-stratify, and optimize patient-specific treatment. However, with each technological advance, we move further from human dependence and closer to fully-automated machine interpretation.

Myocardial Infarct Segmentation From Magnetic Resonance Images for Personalized Modeling of Cardiac Electrophysiology

Accurate representation of myocardial infarct geometry is crucial to patient-specific computational modeling of the heart in ischemic cardiomyopathy. We have developed a methodology for segmentation of left ventricular (LV) infarct from clinically acquired, two-dimensional (2D), late-gadolinium enhanced cardiac magnetic resonance (LGE-CMR) images, for personalized modeling of ventricular electrophysiology. The infarct segmentation was expressed as a continuous min-cut optimization problem, which was solved using its dual formulation, the continuous max-flow (CMF). The optimization objective comprised of a smoothness term, and a data term that quantified the similarity between image intensity histograms of segmented regions and those of a set of training images. A manual segmentation of the LV myocardium was used to initialize and constrain the developed method. The three-dimensional geometry of infarct was reconstructed from its segmentation using an implicit, shape-based interpolation method. The proposed methodology was extensively evaluated using metrics based on geometry, and outcomes of individualized electrophysiological simulations of cardiac dys(function). Several existing LV infarct segmentation approaches were implemented, and compared with the proposed method. Our results demonstrated that the CMF method was more accurate than the existing approaches in reproducing expert manual LV infarct segmentations, and in electrophysiological simulations. The infarct segmentation method we have developed and comprehensively evaluated in this study constitutes an important step in advancing clinical applications of personalized simulations of cardiac electrophysiology.

Evaluation of state-of-the-art segmentation algorithms for left ventricle infarct from late Gadolinium enhancement MR images

Studies have demonstrated the feasibility of late Gadolinium enhancement (LGE) cardiovascular magnetic resonance (CMR) imaging for guiding the management of patients with sequelae to myocardial infarction, such as ventricular tachycardia and heart failure. Clinical implementation of these developments necessitates a reproducible and reliable segmentation of the infarcted regions. It is challenging to compare new algorithms for infarct segmentation in the left ventricle (LV) with existing algorithms. Benchmarking datasets with evaluation strategies are much needed to facilitate comparison. This manuscript presents a benchmarking evaluation framework for future algorithms that segment infarct from LGE CMR of the LV. The image database consists of 30 LGE CMR images of both humans and pigs that were acquired from two separate imaging centres. A consensus ground truth was obtained for all data using maximum likelihood estimation. Six widely-used fixed-thresholding methods and five recently developed algorithms are tested on the benchmarking framework. Results demonstrate that the algorithms have better overlap with the consensus ground truth than most of the n-SD fixed-thresholding methods, with the exception of the Full-Width-at-Half-Maximum (FWHM) fixed-thresholding method. Some of the pitfalls of fixed thresholding methods are demonstrated in this work. The benchmarking evaluation framework, which is a contribution of this work, can be used to test and benchmark future algorithms that detect and quantify infarct in LGE CMR images of the LV. The datasets, ground truth and evaluation code have been made publicly available through the website: https://www.cardiacatlas.org/web/guest/challenges.

Prospective, randomized comparison of gadopentetate and gadobutrol to assess chronic myocardial infarction applying cardiovascular magnetic resonance

Background

We hypothesized that the contrast medium gadobutrol is not inferior compared to Gd-DTPA in identifying and quantifying ischemic late gadolinium enhancement (LGE), even by using a lower dose.

Methods

We prospectively enrolled 30 patients with chronic myocardial infarction as visualized by LGE during clinical routine scan at 1.5 T with 0.20 mmol/kg Gd-DTPA. Participants were randomized to either 0.15 mmol/kg gadobutrol (group A) or 0.10 mmol/kg gadobutrol (group B). CMR protocol was identical in both exams.

LGE was quantified using a semiautomatic approach. Signal intensities of scar, remote myocardium, blood and air were measured. Signal to noise (SNR) and contrast to noise ratios (CNR) were calculated.

Results

Signal intensities were not different between Gd-DTPA and gadobutrol in group A, whereas significant differences were detected in group B. SNR of injured myocardium (53.5+/−21.4 vs. 30.1+/−10.4, p = 0.0001) and CNR between injured and remote myocardium (50.3+/−20.3 vs. 27.3+/−9.3, p < 0.0001) were lower in gadobutrol. Infarct size was lower in both gadobutrol groups compared to Gd-DTPA (group A: 16.8+/−10.2 g vs. 12.8+/−6.8 g, p = 0.03; group B: 18.6+/−12.0 g vs. 14.0+/−9.9 g, p = 0.0016).

Conclusions

Taking application of 0.2 mmol/kg Gd-DTPA as the reference, the delineation of infarct scar was similar with 0.15 mmol/kg gadobutrol, whereas the use 0.10 mmol/kg gadobutrol led to reduced tissue contrast.

Trial registration

The study had been registered under EudraCT Number: 2010-020775-22. Registration date: 2010.08.10

Cardiac magnetic resonance image-based classification of the risk of arrhythmias in post-myocardial infarction patients

INTRODUCTION

Patients surviving myocardial infarction (MI) can be divided into high and low arrhythmic risk groups. Distinguishing between these two groups is of crucial importance since the high-risk group has been shown to benefit from implantable cardioverter defibrillator insertion; a costly surgical procedure with potential complications and no proven advantages for the low-risk group. Currently, markers such as left ventricular ejection fraction and myocardial scar size are used to evaluate arrhythmic risk.

METHODS

In this paper, we propose quantitative discriminative features extracted from late gadolinium enhanced cardiac magnetic resonance images of post-MI patients, to distinguish between 20 high-risk and 34 low-risk patients. These features include size, location, and textural information concerning the scarred myocardium. To evaluate the discriminative power of the proposed features, we used several built-in classification schemes from matrix laboratory (MATLAB) and Waikato environment for knowledge analysis (WEKA) software, including k-nearest neighbor (k-NN), support vector machine (SVM), decision tree, and random forest.

RESULTS

In Experiment 1, the leave-one-out cross-validation scheme is implemented in MATLAB to classify high- and low-risk groups with a classification accuracy of 94.44%, and an AUC of 0.965 for a feature combination that captures size, location and heterogeneity of the scar. In Experiment 2 with the help of WEKA, nested cross-validation is performed with k-NN, SVM, adjusting decision tree and random forest classifiers to differentiate high-risk and low-risk patients. SVM classifier provided average accuracy of 92.6%, and AUC of 0.921 for a feature combination capturing location and heterogeneity of the scar. Experiment 1 and Experiment 2 show that textural features from the scar are important for classification and that localization features provide an additional benefit.

CONCLUSION

These promising results suggest that the discriminative features introduced in this paper can be used by medical professionals, or in automatic decision support systems, along with the recognized risk markers, to improve arrhythmic risk stratification in post-MI patients.

Temporal changes in interpapillary muscle dynamics as an active indicator of mitral valve and left ventricular interaction in ischemic mitral regurgitation

BACKGROUND

Regional subpapillary myocardial hypokinesis may impair lateral reduction in the interpapillary muscle distance (IPMD) from diastole to systole, and adversely affect mitral valve geometry and tethering.

OBJECTIVES

The goal of this study was to investigate the impact of impaired lateral shortening in the interpapillary muscle distance on mitral valve geometry and function in ischemic heart disease.

METHODS

To quantify ventricular size/shape, regional myocardial contraction, lateral shortening of the IPMD, mitral valve geometry, and severity of mitral regurgitation, 67 patients with ischemic heart disease underwent cardiac magnetic resonance imaging, and a correlation analysis of measured parameters was performed. The impact of reduced IPMD shortening on mitral valve (dys)function was confirmed in swine and in a physiological computational mitral valve model.

RESULTS

Lateral shortening of the IPMD from diastole to systole was severely reduced in patients with moderate/severe ischemic mitral regurgitation (9.6 ± 2.8 mm), but preserved in mild IMR (11.5 ± 3.4 mm). Left ventricular size and ejection fraction did not differ between the groups. In swine with subpapillary infarction and impaired IPMD, mitral regurgitation was evident within 1 week, compared to those pigs with a nonpapillary infarction and preserved IPMD. In the controlled computational valve model, IPMD had the maximal impact on regurgitation, and was exacerbated with additional annular dilation.

CONCLUSIONS

By using cardiac magnetic resonance imaging in humans, we demonstrated that it is the impairment of lateral shortening between the papillary muscles, and not passive ventricular size, that governs the severity of mitral regurgitation. Loss of lateral shortening of IPMD tethers the leaflet edges and impairs their systolic closure, resulting in mitral regurgitation, even in small ventricles. Understanding the lateral dynamics of ventricular-valve interactions could aid the development of new repair techniques for ischemic mitral regurgitation.

Preprocedural magnetic resonance imaging for image-guided catheter ablation of scar-related ventricular tachycardia

To present and validate a highly automated MRI analysis workflow for image-guided catheter ablation of scar-related ventricular tachycardia (VT) ablation procedures. A cohort of 15 post-infarction patients underwent MRI prior to VT ablation. The MRI study included a black-blood turbo spin echo sequence for visualizing the aortic root and ostium of the left main (LM) coronary artery, and a 3D late gadolinium enhanced sequence for visualizing the LV anatomy and myocardial scar substrate. Semi-automated segmentation of the LV, aortic root and ostium of LM was performed, followed by fully automated segmentation of myocardial scar. All segmented structures were aligned using an automated image registration algorithm to remove inter-scan displacement. MRI was integrated at the beginning of the procedure after mapping a single LM point. The integration performance was compared to that of the traditional iterative closest point (ICP) method. The proposed method required a single LM mapping point only, compared to 255 ± 43 points with the ICP method. The single-point method achieved a mean point-to-surface distance of 4.9 ± 1.5 mm on the LV surface and 5.1 ± 1.7 mm on the aorta surface (ICP: 3.7 ± 0.8 and 9.2 ± 7.2 mm, P < 0.05). The Cohen's kappa coefficient between the MRI-defined and EAM-defined scar was 0.36 ± 0.16 for the presented method, significantly higher than that of ICP method (0.23 ± 0.21, P = 0.03), indicating more accurate scar substrate localization during integration. This study demonstrated the feasibility of preprocedural MRI integration into the VT ablation procedure, with highly automated image analysis workflow and minimal mapping effort.

Super-resolution reconstruction of late gadolinium-enhanced MRI for improved myocardial scar assessment

PURPOSE

To develop and validate a method for improving image resolution of late gadolinium-enhanced (LGE) magnetic resonance imaging (MRI) for accurate assessment of myocardial scar.

MATERIALS AND METHODS

In a cohort of 37 postinfarction patients, LGE was performed prior to ventricular tachycardia catheter ablation therapy at 1.5T. A super-resolution reconstruction (SRR) technique was applied to the three anisotropic views: short-axis (SA), two-chamber, and four-chamber, to reconstruct a single isotropic volume. For compensation of the interscan heart motion, a joint localized gradient-correlation-based scheme was developed. Scar was identified as either core or gray zone in both the SRR and original SA volumes, and evaluated based on the clinically established bipolar voltage range of the in vivo electroanatomical voltage mapping (EAVM).

RESULTS

Compared to the SA volume, the SRR method resulted in significantly (P < 0.05) reduced myocardial scar gray zone sizes (10.5 ± 8.8 g vs. 9.2 ± 8.1 g) and improved agreement of the bipolar voltage range of scar gray zone (0.99 ± 0.65 mV vs. 1.46 ± 1.15 mV).

CONCLUSION

We propose an SRR method to automatically reconstruct a high-quality isotropic LGE volume from three orthogonal views. Analysis of the in vivo EAVM demonstrated improved myocardial scar assessment from the SRR volume compared with the SA LGE alone.

Myocardial scar identification based on analysis of Look-Locker and 3D late gadolinium enhanced MRI

The aim of this study is to introduce and evaluate an approach for objective and reproducible scar identification from late gadolinium enhanced (LGE) MR by analysis of LGE data with post-contrast T(1) mapping from a routinely acquired T(1) scout Look-Locker (LL) sequence. In 90 post-infarction patients, a LL sequence was acquired prior to a three-dimensional LGE sequence covering the entire left ventricle. In 60/90 patients (training set), the T(1) relaxation rates of remote myocardium and dense myocardial scar were linearly regressed to that of blood. The learned linear relationship was applied to 30/90 patients (validation set) to identify the remote myocardium and dense scar, and to normalize the LGE signal intensity to a range from 0 to 100 %. A 50 % threshold was applied to identify myocardial scar. In the validation set, two observers independently performed manual scar identification, annotated reference regions for the full-width-half-maxima (FWHM) and standard deviation (SD) method, and analyzed the LL sequence for the proposed method. Compared with the manual, FWHM, and SD methods, the proposed method demonstrated the highest inter-class correlation coefficient (0.997) and Dice overlap index (98.7 ± 1.3 %) between the two observers. The proposed method also showed excellent agreement with the gold-standard manual scar identification, with a Dice index of 89.8 ± 7.5 and 90.2 ± 6.6 % for the two observers, respectively. Combined analysis of LL and LGE sequences leads to objective and reproducible myocardial scar identification in post-infarction patients.

A Method to Standardize Quantification of Left Atrial Scar From Delayed-Enhancement MR Images

Delayed-enhancement magnetic resonance imaging (DE-MRI) is an effective technique for detecting left atrial (LA) fibrosis both pre and postradiofrequency ablation for the treatment of atrial fibrillation. Fixed thresholding models are frequently utilized clinically to segment and quantify scar in DE-MRI due to their simplicity. These methods fail to provide a standardized quantification due to interobserver variability. Quantification of scar can be used as an endpoint in clinical studies and therefore standardization is important. In this paper, we propose a segmentation algorithm for LA fibrosis quantification and investigate its performance. The algorithm was validated using numerical phantoms and 15 clinical data sets from patients undergoing LA ablation. We demonstrate that the approach produces good concordance with expert manual delineations. The method offers a standardized quantification technique for evaluation and interpretation of DE-MRI scans.

Fast, accurate, and fully automatic segmentation of the right ventricle in short-axis cardiac MRI

This paper presents a fully automatic method to segment the right ventricle (RV) from short-axis cardiac MRI. A combination of a novel window-constrained accumulator thresholding technique, binary difference of Gaussian (DoG) filters, optimal thresholding, and morphology are utilized to drive the segmentation. A priori segmentation window constraints are incorporated to guide and refine the process, as well as to ensure appropriate area confinement of the segmentation. Training and testing were performed using a combined 48 patient datasets supplied by the organizers of the MICCAI 2012 right ventricle segmentation challenge, allowing for unbiased evaluations and benchmark comparisons. Marked improvements in speed and accuracy over the top existing methods are demonstrated.

Interactive Hierarchical-Flow Segmentation of Scar Tissue From Late-Enhancement Cardiac MR Images

We propose a novel multi-region image segmentation approach to extract myocardial scar tissue from 3-D whole-heart cardiac late-enhancement magnetic resonance images in an interactive manner. For this purpose, we developed a graphical user interface to initialize a fast max-flow-based segmentation algorithm and segment scar accurately with progressive interaction. We propose a partially-ordered Potts (POP) model to multi-region segmentation to properly encode the known spatial consistency of cardiac regions. Its generalization introduces a custom label/region order constraint to Potts model to multi-region segmentation. The combinatorial optimization problem associated with the proposed POP model is solved by means of convex relaxation, for which a novel multi-level continuous max-flow formulation, i.e., the hierarchical continuous max-flow (HMF) model, is proposed and studied. We demonstrate that the proposed HMF model is dual or equivalent to the convex relaxed POP model and introduces a new and efficient hierarchical continuous max-flow based algorithm by modern convex optimization theory. In practice, the introduced hierarchical continuous max-flow based algorithm can be implemented on the parallel GPU to achieve significant acceleration in numerics. Experiments are performed in 50 whole heart 3-D LE datasets, 35 with left-ventricular and 15 with right-ventricular scar. The experimental results are compared to full-width-at-half-maximum and Signal-threshold to reference-mean methods using manual expert myocardial segmentations and operator variabilities and the effect of user interaction are assessed. The results indicate a substantial reduction in image processing time with robust accuracy for detection of myocardial scar. This is achieved without the need for additional region constraints and using a single optimization procedure, substantially reducing the potential for error.

Multi-Stencil Streamline Fast Marching: A General 3-D Framework to Determine Myocardial Thickness and Transmurality in Late Enhancement Images

We propose a fully 3-D methodology for the computation of myocardial nonviable tissue transmurality in contrast enhanced magnetic resonance images. The outcome is a continuous map defined within the myocardium where not only current state-of-the-art measures of transmurality can be calculated, but also information on the location of nonviable tissue is preserved. The computation is done by means of a partial differential equation framework we have called multi-stencil streamline fast marching. Using it, the myocardial and scarred tissue thickness is simultaneously computed. Experimental results show that the proposed 3-D method allows for the computation of transmurality in myocardial regions where current 2-D methods are not able to as conceived, and it also provides more robust and accurate results in situations where the assumptions on which current 2-D methods are based-i.e., there is a visible endocardial contour and its corresponding epicardial points lie on the same slice-, are not met.

Pre to Intraoperative Data Fusion Framework for Multimodal Characterization of Myocardial Scar Tissue

Merging multimodal information about myocardial scar tissue can help electrophysiologists to find the most appropriate target during catheter ablation of ventricular arrhythmias. A framework is presented to analyze and combine information from delayed enhancement magnetic resonance imaging (DE-MRI) and electro-anatomical mapping data. Using this information, electrical, mechanical, and image-based characterization of the myocardium are performed. The presented framework allows the left ventricle to be segmented by DE-MRI and the scar to be characterized prior to the intervention based on image information. It allows the electro-anatomical maps obtained during the intervention from a navigation system to be merged together with the anatomy and scar information extracted from DE-MRI. It also allows for the estimation of endocardial motion and deformation to assess cardiac mechanics. Therefore, electrical, mechanical, and image-based characterization of the myocardium can be performed. The feasibility of this approach was demonstrated on three patients with ventricular tachycardia associated to ischemic cardiomyopathy by integrating images from DE-MRI and electro-anatomical maps data in a common framework for intraoperative myocardial tissue characterization. The proposed framework has the potential to guide and monitor delivery of radio frequency ablation of ventricular tachycardia. It is also helpful for research purposes, facilitating the study of the relationship between electrical and mechanical properties of the tissue, as well as with tissue viability from DE-MRI.

Segmentation of Scarred Myocardium in Cardiac Magnetic Resonance Images

The segmentation of scarred and nonscarred myocardium in Cardiac Magnetic Resonance (CMR) is obtained using different features and feature combinations in a Bayes classifier. The used features are found as a local average of intensity values and the underlying texture information in scarred and nonscarred myocardium. The segmentation classifier was trained and tested with different experimental setups and parameter combinations and was cross validated due to limited data. The experimental results show that the intensity variations are indeed an important feature for good segmentation, and the average area under the Receiver Operating Characteristic (ROC) curve, that is, the AUC, is 91.58 ± 3.2%. The segmentation using texture features also gives good segmentation with average AUC values at 85.89 ± 5.8%, that is, lower than the direct current (DC) feature. However, the texture feature gives robust performance compared to a local mean (DC) feature in a test set simulated from the original CMR data. The segmentation of scarred myocardium is comparable to manual segmentation in all the cross validation cases.

Evaluation of current algorithms for segmentation of scar tissue from late gadolinium enhancement cardiovascular magnetic resonance of the left atrium: an open-access grand challenge

BACKGROUND

Late Gadolinium enhancement (LGE) cardiovascular magnetic resonance (CMR) imaging can be used to visualise regions of fibrosis and scarring in the left atrium (LA) myocardium. This can be important for treatment stratification of patients with atrial fibrillation (AF) and for assessment of treatment after radio frequency catheter ablation (RFCA). In this paper we present a standardised evaluation benchmarking framework for algorithms segmenting fibrosis and scar from LGE CMR images. The algorithms reported are the response to an open challenge that was put to the medical imaging community through an ISBI (IEEE International Symposium on Biomedical Imaging) workshop.

METHODS

The image database consisted of 60 multicenter, multivendor LGE CMR image datasets from patients with AF, with 30 images taken before and 30 after RFCA for the treatment of AF. A reference standard for scar and fibrosis was established by merging manual segmentations from three observers. Furthermore, scar was also quantified using 2, 3 and 4 standard deviations (SD) and full-width-at-half-maximum (FWHM) methods. Seven institutions responded to the challenge: Imperial College (IC), Mevis Fraunhofer (MV), Sunnybrook Health Sciences (SY), Harvard/Boston University (HB), Yale School of Medicine (YL), King's College London (KCL) and Utah CARMA (UTA, UTB). There were 8 different algorithms evaluated in this study.

RESULTS

Some algorithms were able to perform significantly better than SD and FWHM methods in both pre- and post-ablation imaging. Segmentation in pre-ablation images was challenging and good correlation with the reference standard was found in post-ablation images. Overlap scores (out of 100) with the reference standard were as follows: Pre: IC = 37, MV = 22, SY = 17, YL = 48, KCL = 30, UTA = 42, UTB = 45; Post: IC = 76, MV = 85, SY = 73, HB = 76, YL = 84, KCL = 78, UTA = 78, UTB = 72.

CONCLUSIONS

The study concludes that currently no algorithm is deemed clearly better than others. There is scope for further algorithmic developments in LA fibrosis and scar quantification from LGE CMR images. Benchmarking of future scar segmentation algorithms is thus important. The proposed benchmarking framework is made available as open-source and new participants can evaluate their algorithms via a web-based interface.

Probability mapping of scarred myocardium using texture and intensity features in CMR images

Background

The myocardium exhibits heterogeneous nature due to scarring after Myocardial Infarction (MI). In Cardiac Magnetic Resonance (CMR) imaging, Late Gadolinium (LG) contrast agent enhances the intensity of scarred area in the myocardium.

Methods

In this paper, we propose a probability mapping technique using Texture and Intensity features to describe heterogeneous nature of the scarred myocardium in Cardiac Magnetic Resonance (CMR) images after Myocardial Infarction (MI). Scarred tissue and non-scarred tissue are represented with high and low probabilities, respectively. Intermediate values possibly indicate areas where the scarred and healthy tissues are interwoven. The probability map of scarred myocardium is calculated by using a probability function based on Bayes rule. Any set of features can be used in the probability function.

Results

In the present study, we demonstrate the use of two different types of features. One is based on the mean intensity of pixel and the other on underlying texture information of the scarred and non-scarred myocardium. Examples of probability maps computed using the mean intensity of pixel and the underlying texture information are presented. We hypothesize that the probability mapping of myocardium offers alternate visualization, possibly showing the details with physiological significance difficult to detect visually in the original CMR image.

Conclusion

The probability mapping obtained from the two features provides a way to define different cardiac segments which offer a way to identify areas in the myocardium of diagnostic importance (like core and border areas in scarred myocardium).

Three-dimensional segmentation of the left ventricle in late gadolinium enhanced MR images of chronic infarction combining long- and short-axis information

Automatic segmentation of the left ventricle (LV) in late gadolinium enhanced (LGE) cardiac MR (CMR) images is difficult due to the intensity heterogeneity arising from accumulation of contrast agent in infarcted myocardium. In this paper, we present a comprehensive framework for automatic 3D segmentation of the LV in LGE CMR images. Given myocardial contours in cine images as a priori knowledge, the framework initially propagates the a priori segmentation from cine to LGE images via 2D translational registration. Two meshes representing respectively endocardial and epicardial surfaces are then constructed with the propagated contours. After construction, the two meshes are deformed towards the myocardial edge points detected in both short-axis and long-axis LGE images in a unified 3D coordinate system. Taking into account the intensity characteristics of the LV in LGE images, we propose a novel parametric model of the LV for consistent myocardial edge points detection regardless of pathological status of the myocardium (infarcted or healthy) and of the type of the LGE images (short-axis or long-axis). We have evaluated the proposed framework with 21 sets of real patient and four sets of simulated phantom data. Both distance- and region-based performance metrics confirm the observation that the framework can generate accurate and reliable results for myocardial segmentation of LGE images. We have also tested the robustness of the framework with respect to varied a priori segmentation in both practical and simulated settings. Experimental results show that the proposed framework can greatly compensate variations in the given a priori knowledge and consistently produce accurate segmentations.

A comprehensive 3-D framework for automatic quantification of late gadolinium enhanced cardiac magnetic resonance images

Late gadolinium enhanced (LGE) cardiac magnetic resonance (CMR) can directly visualize nonviable myocardium with hyperenhanced intensities with respect to normal myocardium. For heart attack patients, it is crucial to facilitate the decision of appropriate therapy by analyzing and quantifying their LGE CMR images. To achieve accurate quantification, LGE CMR images need to be processed in two steps: segmentation of the myocardium followed by classification of infarcts within the segmented myocardium. However, automatic segmentation is difficult usually due to the intensity heterogeneity of the myocardium and intensity similarity between the infarcts and blood pool. Besides, the slices of an LGE CMR dataset often suffer from spatial and intensity distortions, causing further difficulties in segmentation and classification. In this paper, we present a comprehensive 3-D framework for automatic quantification of LGE CMR images. In this framework, myocardium is segmented with a novel method that deforms coupled endocardial and epicardial meshes and combines information in both short- and long-axis slices, while infarcts are classified with a graph-cut algorithm incorporating intensity and spatial information. Moreover, both spatial and intensity distortions are effectively corrected with specially designed countermeasures. Experiments with 20 sets of real patient data show visually good segmentation and classification results that are quantitatively in strong agreement with those manually obtained by experts.

Improved myocardial scar characterization by super-resolution reconstruction in late gadolinium enhanced MRI

Image resolution is an important factor for accurate myocardial scar assessment from late gadolinium enhanced (LGE) MR. It has been shown that the conventionally used short-axis (SA) LGE acquisition with anisotropic resolution may overestimate the scar size due to partial volume effect, undermining the prognostic and diagnostic accuracy of LGE MRI in critical clinical applications. In this work, we present a method for combining three complementary anisotropic orthogonal LGE sequences of the heart region into a single isotropic volume. Our algorithm is based on the super-resolution reconstruction technique and employs joint localized gradient-correlation-based technique for compensation of breathing motion. The proposed method was validated on the gold standard electroanatomical voltage mapping (EAVM) data of 15 post-infarction patients. The reconstructed myocardial scar image demonstrated improved agreement with the EAVM compared to the conventional SA image, especially at the clinically significant gray zone region.

Automated quantification of myocardial infarction using graph cuts on contrast delayed enhanced magnetic resonance images

In this work, we propose a semi-automated myocardial infarction quantification method for cardiac contrast delayed enhancement magnetic resonance images (DE-MRI). Advantages of this method include that it reduces manual contouring of the left ventricle, obviates a remote myocardium region, and automatically distinguishes infarct, healthy and heterogeneous ("gray zone") tissue despite variability in intensity and noise across images. Quantitative evaluation results showed that the automatically determined infarct core and gray zone size have high correlation with that derived from the averaged results of the manual full width at half maximum (FWHM) methods (R(2)=0.99 for infarct core and gray zone size). Compared with the manual method, a much better reproducibility was achieved with the proposed algorithm and it shortens the evaluation time to one second per image, compared with 2-5 min per image for the manual method.

Automatic classification of scar tissue in late gadolinium enhancement cardiac MRI for the assessment of left-atrial wall injury after radiofrequency ablation

Radiofrequency ablation is a promising procedure for treating atrial fibrillation (AF) that relies on accurate lesion delivery in the left atrial (LA) wall for success. Late Gadolinium Enhancement MRI (LGE MRI) at three months post-ablation has proven effective for noninvasive assessment of the location and extent of scar formation, which are important factors for predicting patient outcome and planning of redo ablation procedures. We have developed an algorithm for automatic classification in LGE MRI of scar tissue in the LA wall and have evaluated accuracy and consistency compared to manual scar classifications by expert observers. Our approach clusters voxels based on normalized intensity and was chosen through a systematic comparison of the performance of multivariate clustering on many combinations of image texture. Algorithm performance was determined by overlap with ground truth, using multiple overlap measures, and the accuracy of the estimation of the total amount of scar in the LA. Ground truth was determined using the STAPLE algorithm, which produces a probabilistic estimate of the true scar classification from multiple expert manual segmentations. Evaluation of the ground truth data set was based on both inter- and intra-observer agreement, with variation among expert classifiers indicating the difficulty of scar classification for a given a dataset. Our proposed automatic scar classification algorithm performs well for both scar localization and estimation of scar volume: for ground truth datasets considered easy, variability from the ground truth was low; for those considered difficult, variability from ground truth was on par with the variability across experts.