Unsupervised 4D myocardium segmentation with a Markov Random Field based deformable model

Automatic detection of myocardial ischaemia using generalisable spatio-temporal hierarchical Bayesian modelling of DCE-MRI

Dynamic Contrast Enhanced Magnetic Resonance Imaging (DCE-MRI) can be used as a non-invasive method for the assessment of myocardial perfusion. The acquired images can be utilised to analyse the spatial extent and severity of myocardial ischaemia (regions with impaired microvascular blood flow). In the present paper, we propose a novel generalisable spatio-temporal hierarchical Bayesian model (GST-HBM) to automate the detection of ischaemic lesions and improve the in silico prediction accuracy by systematically integrating spatio-temporal context information. We present a computational inference procedure with an adequate trade-off between accuracy and computational efficiency, whereby model parameters are sampled from the posterior distribution with Gibbs sampling, while lower-level hyperparameters are selected using model selection strategies based on the Watanabe Akaike information criterion (WAIC). We have assessed our method on both synthetic (in silico) data with known gold-standard and 12 sets of clinical first-pass myocardial perfusion DCE-MRI datasets. We have also carried out a comparative performance evaluation with four established alternative methods: Gaussian mixture model (GMM), opening and closing operations based on Gaussian mixture model (GMMC&Omax), Markov random field constrained Gaussian mixture model (GMM-MRF) and model-based hierarchical Bayesian model (M-HBM). Our results show that the proposed GST-HBM method achieves much higher in silico prediction accuracy than the established alternative methods. Furthermore, this method appears to provide a more robust delineation of ischaemic lesions in datasets affected by spatially variant noise.

Topology-sensitive weighting model for myocardial segmentation

Accurate myocardial segmentation is crucial for the diagnosis of various heart diseases. However, segmentation results often suffer from topology structural errors, such as broken connections and holes, especially in cases of poor image quality. These errors are unacceptable in clinical diagnosis. We proposed a Topology-Sensitive Weight (TSW) model to keep both pixel-wise accuracy and topological correctness. Specifically, the Position Weighting Update (PWU) strategy with the Boundary-Sensitive Topology (BST) module can guide the model to focus on positions where topological features are sensitive to pixel values. The Myocardial Integrity Topology (MIT) module can serve as a guide for maintaining myocardial integrity. We evaluate the TSW model on the CAMUS dataset and a private echocardiography myocardial segmentation dataset. The qualitative and quantitative experimental results show that the TSW model significantly enhances topological accuracy while maintaining pixel-wise precision.

An automatic segmentation method with self-attention mechanism on left ventricle in gated PET/CT myocardial perfusion imaging

OBJECTIVES

We aimed to propose an automatic segmentation method for left ventricular (LV) from 16 electrocardiogram (ECG) -gated ¹³N-NH₃ PET/CT myocardial perfusion imaging (MPI) to improve the performance of LV function assessment.

METHODS

Ninety-six cases with confirmed or suspected obstructive coronary artery disease (CAD) were enrolled in this research. The LV myocardial contours were delineated by physicians as ground truth. We developed an automatic segmentation method, which introduces the self-attention mechanism into 3D U-Net to capture global information of images so as to achieve fine segmentation of LV. Three cross-validation tests were performed on each gate (64 vs. 32 for training vs. validation). The effectiveness was validated by quantitative metrics (modified hausdorff distance, MHD; dice ratio, DR; 3D MHD) as well as cardiac functional parameters (end-systolic volume, ESV; end-diastolic volume, EDV; ejection fraction, EF). Furthermore, the feasibility of the proposed method was also evaluated by intra- and inter-observers with DR and 3D-MHD.

RESULTS

Compared with backbone network, the proposed approach improved the average DR from 0.905 ± 0.0193 to 0.9202 ± 0.0164, and decreased the average 3D MHD from 0.4611 ± 0.0349 to 0.4304 ± 0.0339. The average relative error of LV volume between proposed method and ground truth is 1.09±3.66%, and the correlation coefficient is 0.992 ± 0.007 (P < 0.001). The EDV, ESV, EF deduced from the proposed approach were highly correlated with ground truth (r ≥ 0.864, P < 0.001), and the correlation with commercial software is fair (r ≥ 0.871, P < 0.001). DR and 3D MHD of contours and myocardium from two observers are higher than 0.899 and less than 0.5194.

CONCLUSION

The proposed approach is highly feasible for automatic segmentation of the LV cavity and myocardium, with potential to benefit the precision of LV function assessment.

Ultrafast Computation of Left Ventricular Ejection Fraction by Using Temporal Intensity Variation in Cine Cardiac Magnetic Resonance

Cardiac magnetic resonance enables comprehensive cardiac evaluation; however, intense time and labor requirements for data acquisition and processing have discouraged many clinicians from using it. We have developed an alternative image-processing algorithm that requires minimal user interaction: an ultrafast algorithm that computes left ventricular ejection fraction (LVEF) by using temporal intensity variation in cine balanced steady-state free precession (bSSFP) short-axis images, with or without contrast medium. We evaluated the algorithm's performance against an expert observer's analysis for segmenting the LV cavity in 65 study participants (LVEF range, 12%-70%). In 12 instances, contrast medium was administered before cine imaging. Bland-Altman analysis revealed quantitative effects of LV basal, midcavity, and apical morphologic variation on the algorithm's accuracy. Total computation time for the LV stack was <2.5 seconds. The algorithm accurately delineated endocardial boundaries in 1,132 of 1,216 slices (93%). When contours in the extreme basal and apical slices were not adequate, they were replaced with manually drawn contours. The Bland-Altman mean differences were <1.2 mL (0.8%) for end-diastolic volume, <5 mL (6%) for end-systolic volume, and <3% for LVEF. Standard deviation of the difference was ≤4.1% of LV volume for all sections except the midcavity in end-systole (8.3% of end-systolic volume). We conclude that temporal intensity variation-based ultrafast LVEF computation is clinically accurate across a range of LV shapes and wall motions and is suitable for postcontrast cine SSFP imaging. Our algorithm enables real-time processing of cine bSSFP images on a commercial scanner console within 3 seconds in an unobtrusive automated process.

Discrete light sheet microscopic segmentation of left ventricle using morphological tuning and active contours

Left ventricular segmentation using cardiovascular MR scan is required for the diagnosis and further cure of cardiac diseases. Automatic systems for left ventricle segmentation are being studied for attaining more accurate results in a shorter period of time. A novel algorithm introducing discrete segmentation of left ventricle achieves an independent processing of images swiftly. The workflow consists of four segments; first, automated localization is performed on the MR image. Second, performing preprocessing intimately improves and enhances the quality of image using mean contrast adjustment. Central segmentation of endocardium and epicardium layers includes novel MTAC (Morphological tuning using active contours) segmentation algorithm that provides a perfect combination of active contours and morphological tuning to bring an adequate and desirable segmentation. The prospective snake model is a restrained progression, which takes iterations for an impulse throughout the left ventricle contours. At the end, contrast based refining overcomes minor edge problems for both outer and inner boundaries. Proposed algorithm is evaluated via Sunnybrook cardiac MR images by producing an overall average perpendicular distance 2.45 mm, an average dice matrix (endo: 91.3%; epi: 92.16%) and 91.7% dice matrix of overall endocardium and epicardium contours from ground truth contours.

Feasibility of Automatic Seed Generation Applied to Cardiac MRI Image Analysis

We present a method of using interactive image segmentation algorithms to reduce specific image segmentation problems to the task of finding small sets of pixels identifying the regions of interest. To this end, we empirically show the feasibility of automatically generating seeds for GrowCut, a popular interactive image segmentation algorithm. The principal contribution of our paper is the proposal of a method for automating the seed generation method for the task of whole-heart segmentation of MRI scans, which achieves competitive unsupervised results (0.76 Dice on the MMWHS dataset). Moreover, we show that segmentation performance is robust to seeds with imperfect precision, suggesting that GrowCut-like algorithms can be applied to medical imaging tasks with little modeling effort.

Left ventricle segmentation using a Bayesian approach with distance dependent shape priors

We propose a method for segmentation of the left ventricle in magnetic resonance cardiac images. The framework consists of an initial Bayesian segmentation of the central slice of the volume. This segmentation is used to locate a shape prior for the LV myocardial tissue. This shape prior is determined using the fact that the myocardium is approximately annular as seen in the short-axis. Then a second Bayesian segmentation is performed to obtain the final result. This procedure is repeated for the rest of the slices. An extrapolation of the area of the LV is used to determine a stopping criterion. The method was evaluated on the databases of the Cardiac Atlas project. Our results demonstrate a suitable accuracy for myocardial segmentation (≈0.8 Dice's coefficient). For the endocardium and the epicardium the Dice's coefficients are 0.94 and 0.9 respectively. The accuracy was also evaluated in terms of the Hausdorff distance and the average distance. For the myocardium we obtain 8 mm and 2 mm respectively. Our results demonstrate the capability and merits of the proposed method to estimate the structure of the LV. The method requires minimal user input and generates results with quality comparable to more complex approaches. This paper suggests a new efficient approach for automatic LV quantification based on a Bayesian technique with shape priors with errors comparable to state-of-the-art techniques.

Automatic quantification of the LV function and mass: A deep learning approach for cardiovascular MRI

OBJECTIVE

This paper proposes a novel approach for automatic left ventricle (LV) quantification using convolutional neural networks (CNN).

METHODS

The general framework consists of one CNN for detecting the LV, and another for tissue classification. Also, three new deep learning architectures were proposed for LV quantification. These new CNNs introduce the ideas of sparsity and depthwise separable convolution into the U-net architecture, as well as, a residual learning strategy level-to-level. To this end, we extend the classical U-net architecture and use the generalized Jaccard distance as optimization objective function.

RESULTS

The CNNs were trained and evaluated with 140 patients from two public cardiovascular magnetic resonance datasets (Sunnybrook and Cardiac Atlas Project) by using a 5-fold cross-validation strategy. Our results demonstrate a suitable accuracy for myocardial segmentation ( ∼ 0.9 Dice's coefficient), and a strong correlation with the most relevant physiological measures: 0.99 for end-diastolic and end-systolic volume, 0.97 for the left myocardial mass, 0.95 for the ejection fraction and 0.93 for the stroke volume and cardiac output.

CONCLUSION

Our simulation and clinical evaluation results demonstrate the capability and merits of the proposed CNN to estimate different structural and functional features such as LV mass and EF which are commonly used for both diagnosis and treatment of different pathologies.

SIGNIFICANCE

This paper suggests a new approach for automatic LV quantification based on deep learning where errors are comparable to the inter- and intra-operator ranges for manual contouring.

Fast left ventricle tracking using localized anatomical affine optical flow

In daily clinical cardiology practice, left ventricle (LV) global and regional function assessment is crucial for disease diagnosis, therapy selection, and patient follow-up. Currently, this is still a time-consuming task, spending valuable human resources. In this work, a novel fast methodology for automatic LV tracking is proposed based on localized anatomically constrained affine optical flow. This novel method can be combined to previously proposed segmentation frameworks or manually delineated surfaces at an initial frame to obtain fully delineated datasets and, thus, assess both global and regional myocardial function. Its feasibility and accuracy were investigated in 3 distinct public databases, namely in realistically simulated 3D ultrasound, clinical 3D echocardiography, and clinical cine cardiac magnetic resonance images. The method showed accurate tracking results in all databases, proving its applicability and accuracy for myocardial function assessment. Moreover, when combined to previous state-of-the-art segmentation frameworks, it outperformed previous tracking strategies in both 3D ultrasound and cardiac magnetic resonance data, automatically computing relevant cardiac indices with smaller biases and narrower limits of agreement compared to reference indices. Simultaneously, the proposed localized tracking method showed to be suitable for online processing, even for 3D motion assessment. Importantly, although here evaluated for LV tracking only, this novel methodology is applicable for tracking of other target structures with minimal adaptations.

Unsupervised Myocardial Segmentation for Cardiac BOLD

A fully automated 2-D+time myocardial segmentation framework is proposed for cardiac magnetic resonance (CMR) blood-oxygen-level-dependent (BOLD) data sets. Ischemia detection with CINE BOLD CMR relies on spatio-temporal patterns in myocardial intensity, but these patterns also trouble supervised segmentation methods, the de facto standard for myocardial segmentation in cine MRI. Segmentation errors severely undermine the accurate extraction of these patterns. In this paper, we build a joint motion and appearance method that relies on dictionary learning to find a suitable subspace. Our method is based on variational pre-processing and spatial regularization using Markov random fields, to further improve performance. The superiority of the proposed segmentation technique is demonstrated on a data set containing cardiac phase-resolved BOLD MR and standard CINE MR image sequences acquired in baseline and ischemic condition across ten canine subjects. Our unsupervised approach outperforms even supervised state-of-the-art segmentation techniques by at least 10% when using Dice to measure accuracy on BOLD data and performs at par for standard CINE MR. Furthermore, a novel segmental analysis method attuned for BOLD time series is utilized to demonstrate the effectiveness of the proposed method in preserving key BOLD patterns.

A 3D Hermite-based multiscale local active contour method with elliptical shape constraints for segmentation of cardiac MR and CT volumes

Analysis of cardiac images is a fundamental task to diagnose heart problems. Left ventricle (LV) is one of the most important heart structures used for cardiac evaluation. In this work, we propose a novel 3D hierarchical multiscale segmentation method based on a local active contour (AC) model and the Hermite transform (HT) for LV analysis in cardiac magnetic resonance (MR) and computed tomography (CT) volumes in short axis view. Features such as directional edges, texture, and intensities are analyzed using the multiscale HT space. A local AC model is configured using the HT coefficients and geometrical constraints. The endocardial and epicardial boundaries are used for evaluation. Segmentation of the endocardium is controlled using elliptical shape constraints. The final endocardial shape is used to define the geometrical constraints for segmentation of the epicardium. We follow the assumption that epicardial and endocardial shapes are similar in volumes with short axis view. An initialization scheme based on a fuzzy C-means algorithm and mathematical morphology was designed. The algorithm performance was evaluated using cardiac MR and CT volumes in short axis view demonstrating the feasibility of the proposed method.

Automatic initialization and quality control of large-scale cardiac MRI segmentations

Continuous advances in imaging technologies enable ever more comprehensive phenotyping of human anatomy and physiology. Concomitant reduction of imaging costs has resulted in widespread use of imaging in large clinical trials and population imaging studies. Magnetic Resonance Imaging (MRI), in particular, offers one-stop-shop multidimensional biomarkers of cardiovascular physiology and pathology. A wide range of analysis methods offer sophisticated cardiac image assessment and quantification for clinical and research studies. However, most methods have only been evaluated on relatively small databases often not accessible for open and fair benchmarking. Consequently, published performance indices are not directly comparable across studies and their translation and scalability to large clinical trials or population imaging cohorts is uncertain. Most existing techniques still rely on considerable manual intervention for the initialization and quality control of the segmentation process, becoming prohibitive when dealing with thousands of images. The contributions of this paper are three-fold. First, we propose a fully automatic method for initializing cardiac MRI segmentation, by using image features and random forests regression to predict an initial position of the heart and key anatomical landmarks in an MRI volume. In processing a full imaging database, the technique predicts the optimal corrective displacements and positions in relation to the initial rough intersections of the long and short axis images. Second, we introduce for the first time a quality control measure capable of identifying incorrect cardiac segmentations with no visual assessment. The method uses statistical, pattern and fractal descriptors in a random forest classifier to detect failures to be corrected or removed from subsequent statistical analysis. Finally, we validate these new techniques within a full pipeline for cardiac segmentation applicable to large-scale cardiac MRI databases. The results obtained based on over 1200 cases from the Cardiac Atlas Project show the promise of fully automatic initialization and quality control for population studies.

Segmentation of medical images using mean value guided contour

Partial differential equation-based (PDE-based) methods are extensively used in image segmentation, especially in contour models. Difficulties associated with the boundaries, namely troubles with developing initialization, inadequate convergence to boundary concavities, and difficulties connected to saddle points and stationary points of active contours make the contour models suffer from a feeble performance of referring to complex geometries. The present paper is designed to take advantage of mean value theorem rather than minimizing energy function for contours. It is efficiently capable of resolving above-mentioned problems by applying this theorem to the edge map gradient vectors, which is calculated from the image. Since the contour is computed in a straightforward manner from an edge map instead of force balance equation, it varies from other contour-based image segmentation methods. To illustrate the ability of the proposed model in complex geometries and ruptures, several experiments were also provided to validate the model. The experiments' results demonstrated that the proposed method, which is called mean value guided contour (MVGC), is capable of repositioning contours into boundary concavities and has suitable forcefulness in complex geometries.

A novel myocardium segmentation approach based on neutrosophic active contour model

BACKGROUND AND OBJECTIVES

Automatic delineation of the myocardium in echocardiography can assist radiologists to diagnosis heart problems. However, it is still challenging to distinguish myocardium from other tissue due to a low signal-to-noise ratio, low contrast, vague boundary, and speckle noise. The purpose of this study is to automatically detect myocardium region in left ventricle myocardial contrast echocardiography (LVMCE) images to help radiologists' diagnosis and further measurement on infarction size.

METHODS

The LVMCE image is firstly mapped into neutrosophic similarity (NS) domain using the intensity and homogeneity features. Then, a neutrosophic active contour model (NACM) is proposed and the energy function is defined by the NS values. Finally, the ventricle is detected using the curve evolving results. The ventricle's boundary is identified as the endocardium. To speed up the evolution procedure and increase the detection accuracy, a clustering algorithm is employed to obtain the initial ventricle region. The curve evolution procedure in NACM is utilized again to obtain the epicardium, where the initial contour uses the detected endocardium and the anatomy knowledge on the thickness of the myocardium.

RESULTS

Echocardiographic studies are performed on 10 male Sprague-Dawley rats using a Vivid 7 system including 5 normal cases and 5 rats with myocardial infarction. The myocardium boundaries manually outlined by an experienced radiologist are used as the reference standard for the performance evaluation. Two metrics, Hdist and AvgDist, are employed to evaluate the detection results. The NACM method was compared with those from the eliminated particle swarm optimization (EPSO) and active contour model without edges (ACMWE) methods. The mean and standard deviation of the Hdist and AvgDist on endocardium are 6.83 ± 1.12mm and 0.79 ± 0.28mm using EPSO method, 7.12 ± 0.98mm and 0.82 ± 0.32mm using ACMWE method, and 4.55 ± 0.9mm and 0.58 ± 0.18mm using NACM method, respectively. The improvement on epicardium is much more significant, and two metrics are decreased from 7.45 ± 1.24mm, and 1.47 ± 0.34mm using EPSO method, and 8.21±0.43mm, and 1.73±0.47mm using ACMWE method, to 4.94 ± 0.82mm, and 0.84 ± 0.22mm using NACM method, respectively.

CONCLUSIONS

The proposed method can automatically detect myocardium accurately, and is helpful for clinical therapeutics to measure myocardial perfusion and infarct size.

Segmentation of myocardium from cardiac MR images using a novel dynamic programming based segmentation method

PURPOSE

Myocardium segmentation in cardiac magnetic resonance (MR) images plays a vital role in clinical diagnosis of the cardiovascular diseases. Because of the low contrast and large variation in intensity and shapes, myocardium segmentation has been a challenging task. A dynamic programming (DP) based segmentation method, incorporating the likelihood and shape information of the myocardium, is developed for segmenting myocardium in cardiac MR images.

METHODS

The endocardium, i.e., the left ventricle blood cavity, is segmented for initialization, and then the optimal epicardium contour is determined using the polar-transformed image and DP scheme. In the DP segmentation scheme, three techniques are proposed to improve the segmentation performance: (1) the likelihood image of the myocardium is constructed to define the external cost in the DP, thus the cost function incorporates prior probability estimation; (2) the adaptive search range is introduced to determine the polar-transformed image, thereby excluding irrelevant tissues; (3) the connectivity constrained DP algorithm is developed to obtain an optimal closed contour. Four metrics, including the Dice metric (Dice), root mean squared error (RMSE), reliability, and correlation coefficient, are used to assess the segmentation accuracy. The authors evaluated the performance of the proposed method on a private dataset and the MICCAI 2009 challenge dataset. The authors also explored the effects of the three new techniques of the DP scheme in the proposed method.

RESULTS

For the qualitative evaluation, the segmentation results of the proposed method were clinically acceptable. For the quantitative evaluation, the mean (Dice) for the endocardium and epicardium was 0.892 and 0.927, respectively; the mean RMSE was 2.30 mm for the endocardium and 2.39 mm for the epicardium. In addition, the three new techniques in the proposed DP scheme, i.e., the likelihood image of the myocardium, the adaptive search range, and the connectivity constrained DP algorithm, improved the segmentation performance for the epicardium with 0.029, 0.047, and 0.007 in terms of the Dice and 0.98, 1.31, and 0.21 mm in terms of the RMSE, respectively.

CONCLUSIONS

The three techniques (the likelihood image of the myocardium, the adaptive search range, and the connectivity constrained DP algorithm) can improve the segmentation ability of the DP method, and the proposed method with these techniques has the ability to achieve the acceptable segmentation result of the myocardium in cardiac MR images. Therefore, the proposed method would be useful in clinical diagnosis of the cardiovascular diseases.

A review of heart chamber segmentation for structural and functional analysis using cardiac magnetic resonance imaging

Cardiovascular magnetic resonance (CMR) has become a key imaging modality in clinical cardiology practice due to its unique capabilities for non-invasive imaging of the cardiac chambers and great vessels. A wide range of CMR sequences have been developed to assess various aspects of cardiac structure and function, and significant advances have also been made in terms of imaging quality and acquisition times. A lot of research has been dedicated to the development of global and regional quantitative CMR indices that help the distinction between health and pathology. The goal of this review paper is to discuss the structural and functional CMR indices that have been proposed thus far for clinical assessment of the cardiac chambers. We include indices definitions, the requirements for the calculations, exemplar applications in cardiovascular diseases, and the corresponding normal ranges. Furthermore, we review the most recent state-of-the art techniques for the automatic segmentation of the cardiac boundaries, which are necessary for the calculation of the CMR indices. Finally, we provide a detailed discussion of the existing literature and of the future challenges that need to be addressed to enable a more robust and comprehensive assessment of the cardiac chambers in clinical practice.

A combined deep-learning and deformable-model approach to fully automatic segmentation of the left ventricle in cardiac MRI

Segmentation of the left ventricle (LV) from cardiac magnetic resonance imaging (MRI) datasets is an essential step for calculation of clinical indices such as ventricular volume and ejection fraction. In this work, we employ deep learning algorithms combined with deformable models to develop and evaluate a fully automatic LV segmentation tool from short-axis cardiac MRI datasets. The method employs deep learning algorithms to learn the segmentation task from the ground true data. Convolutional networks are employed to automatically detect the LV chamber in MRI dataset. Stacked autoencoders are used to infer the LV shape. The inferred shape is incorporated into deformable models to improve the accuracy and robustness of the segmentation. We validated our method using 45 cardiac MR datasets from the MICCAI 2009 LV segmentation challenge and showed that it outperforms the state-of-the art methods. Excellent agreement with the ground truth was achieved. Validation metrics, percentage of good contours, Dice metric, average perpendicular distance and conformity, were computed as 96.69%, 0.94, 1.81 mm and 0.86, versus those of 79.2-95.62%, 0.87-0.9, 1.76-2.97 mm and 0.67-0.78, obtained by other methods, respectively.

Clinical feasibility of a myocardial signal intensity threshold-based semi-automated cardiac magnetic resonance segmentation method

OBJECTIVES

To assess the accuracy and efficiency of a threshold-based, semi-automated cardiac MRI segmentation algorithm in comparison with conventional contour-based segmentation and aortic flow measurements.

METHODS

Short-axis cine images of 148 patients (55 ± 18 years, 81 men) were used to evaluate left ventricular (LV) volumes and mass (LVM) using conventional and threshold-based segmentations. Phase-contrast images were used to independently measure stroke volume (SV). LV parameters were evaluated by two independent readers.

RESULTS

Evaluation times using the conventional and threshold-based methods were 8.4 ± 1.9 and 4.2 ± 1.3 min, respectively (P < 0.0001). LV parameters measured by the conventional and threshold-based methods, respectively, were end-diastolic volume (EDV) 146 ± 59 and 134 ± 53 ml; end-systolic volume (ESV) 64 ± 47 and 59 ± 46 ml; SV 82 ± 29 and 74 ± 28 ml (flow-based 74 ± 30 ml); ejection fraction (EF) 59 ± 16 and 58 ± 17%; and LVM 141 ± 55 and 159 ± 58 g. Significant differences between the conventional and threshold-based methods were observed in EDV, ESV, and LVM mesurements; SV from threshold-based and flow-based measurements were in agreement (P > 0.05) but were significantly different from conventional analysis (P < 0.05). Excellent inter-observer agreement was observed.

CONCLUSIONS

Threshold-based LV segmentation provides improved accuracy and faster assessment compared to conventional contour-based methods.

KEY POINTS

• Threshold-based left ventricular segmentation provides time-efficient assessment of left ventricular parameters • The threshold-based method can discriminate between blood and papillary muscles • This method provides improved accuracy compared to aortic flow measurements as a reference.

Myocardial Iron Loading Assessment by Automatic Left Ventricle Segmentation with Morphological Operations and Geodesic Active Contour on T2* images

Myocardial iron loading thalassemia patients could be identified using T2* magnetic resonance images (MRI). To quantitatively assess cardiac iron loading, we proposed an effective algorithm to segment aligned free induction decay sequential myocardium images based on morphological operations and geodesic active contour (GAC). Nine patients with thalassemia major were recruited (10 male and 16 female) to undergo a thoracic MRI scan in the short axis view. Free induction decay images were registered for T2* mapping. The GAC were utilized to segment aligned MR images with a robust initialization. Segmented myocardium regions were divided into sectors for a region-based quantification of cardiac iron loading. Our proposed automatic segmentation approach achieve a true positive rate at 84.6% and false positive rate at 53.8%. The area difference between manual and automatic segmentation was 25.5% after 1000 iterations. Results from T2* analysis indicated that regions with intensity lower than 20 ms were suffered from heavy iron loading in thalassemia major patients. The proposed method benefited from abundant edge information of the free induction decay sequential MRI. Experiment results demonstrated that the proposed method is feasible in myocardium segmentation and was clinically applicable to measure myocardium iron loading.

Simplified post processing of cine DENSE cardiovascular magnetic resonance for quantification of cardiac mechanics

BACKGROUND

Cardiovascular magnetic resonance using displacement encoding with stimulated echoes (DENSE) is capable of assessing advanced measures of cardiac mechanics such as strain and torsion. A potential hurdle to widespread clinical adoption of DENSE is the time required to manually segment the myocardium during post-processing of the images. To overcome this hurdle, we proposed a radical approach in which only three contours per image slice are required for post-processing (instead of the typical 30-40 contours per image slice). We hypothesized that peak left ventricular circumferential, longitudinal and radial strains and torsion could be accurately quantified using this simplified analysis.

METHODS AND RESULTS

We tested our hypothesis on a large multi-institutional dataset consisting of 541 DENSE image slices from 135 mice and 234 DENSE image slices from 62 humans. We compared measures of cardiac mechanics derived from the simplified post-processing to those derived from original post-processing utilizing the full set of 30-40 manually-defined contours per image slice. Accuracy was assessed with Bland-Altman limits of agreement and summarized with a modified coefficient of variation. The simplified technique showed high accuracy with all coefficients of variation less than 10% in humans and 6% in mice. The accuracy of the simplified technique was also superior to two previously published semi-automated analysis techniques for DENSE post-processing.

CONCLUSIONS

Accurate measures of cardiac mechanics can be derived from DENSE cardiac magnetic resonance in both humans and mice using a simplified technique to reduce post-processing time by approximately 94%. These findings demonstrate that quantifying cardiac mechanics from DENSE data is simple enough to be integrated into the clinical workflow.

Tracking using motion estimation with physically motivated inter-region constraints

We propose a method for tracking structures (e.g., ventricles and myocardium) in cardiac images (e.g., magnetic resonance) by propagating forward in time a previous estimate of the structures using a new physically motivated motion estimation scheme. Our method estimates motion by regularizing only within structures so that differing motions among different structures are not mixed. It simultaneously satisfies the physical constraints at the interface between a fluid and a medium that the normal component of the fluid's motion must match the normal component of the medium's motion and the No-Slip condition, which states that the tangential velocity approaches zero near the interface. We show that these conditions lead to partial differential equations with Robin boundary conditions at the interface, which couple the motion between structures. We show that propagating a segmentation across frames using our motion estimation scheme leads to more accurate segmentation than traditional motion estimation that does not use physical constraints. Our method is suited to interactive segmentation, prominently used in commercial applications for cardiac analysis, where segmentation propagation is used to predict a segmentation in the next frame. We show that our method leads to more accurate predictions than a popular and recent interactive method used in cardiac segmentation.

Integrated segmentation and interpolation of sparse data

We address the two inherently related problems of segmentation and interpolation of 3D and 4D sparse data and propose a new method to integrate these stages in a level set framework. The interpolation process uses segmentation information rather than pixel intensities for increased robustness and accuracy. The method supports any spatial configurations of sets of 2D slices having arbitrary positions and orientations. We achieve this by introducing a new level set scheme based on the interpolation of the level set function by radial basis functions. The proposed method is validated quantitatively and/or subjectively on artificial data and MRI and CT scans and is compared against the traditional sequential approach, which interpolates the images first, using a state-of-the-art image interpolation method, and then segments the interpolated volume in 3D or 4D. In our experiments, the proposed framework yielded similar segmentation results to the sequential approach but provided a more robust and accurate interpolation. In particular, the interpolation was more satisfactory in cases of large gaps, due to the method taking into account the global shape of the object, and it recovered better topologies at the extremities of the shapes where the objects disappear from the image slices. As a result, the complete integrated framework provided more satisfactory shape reconstructions than the sequential approach.

Groupwise elastic registration by a new sparsity-promoting metric: application to the alignment of cardiac magnetic resonance perfusion images

This paper proposes a methodology for the joint alignment of a sequence of images based on a groupwise registration procedure by using a new family of metrics that exploit the expected sparseness of the temporal intensity curves corresponding to the aligned points. Therefore, this methodology is able to tackle the alignment of temporal sequences of images in which the represented phenomenon varies in time. Specifically, we have applied it to the correction of motion in contrast-enhanced first-pass perfusion cardiac magnetic resonance images. The time sequence is elastically registered as a whole by using the aforementioned family of multi-image metrics and jointly optimizing the parameters of the transformations involved. The proposed metrics are able to cope with dynamic changes in the intensity content of corresponding points in the sequence guided by the assumption that these changes allow for a sparse representation in a properly selected frame. Results have shown the statistically significant improvement in the performance of the proposed metric with respect to previous groupwise registration metrics for the problem at hand, which is especially relevant to correct for elastic deformations.

Automatic functional analysis of left ventricle in cardiac cine MRI

RATIONALE AND OBJECTIVES

A fully automated left ventricle segmentation method for the functional analysis of cine short axis (SAX) magnetic resonance (MR) images was developed, and its performance evaluated with 133 studies of subjects with diverse pathology: ischemic heart failure (n=34), non-ischemic heart failure (n=30), hypertrophy (n=32), and healthy (n=37).

MATERIALS AND METHODS

The proposed automatic method locates the left ventricle (LV), then for each image detects the contours of the endocardium, epicardium, papillary muscles and trabeculations. Manually and automatically determined contours and functional parameters were compared quantitatively.

RESULTS

There was no significant difference between automatically and manually determined end systolic volume (ESV), end diastolic volume (EDV), ejection fraction (EF) and left ventricular mass (LVM) for each of the four groups (paired sample t-test, α=0.05). The automatically determined functional parameters showed high correlations with those derived from manual contours, and the Bland-Altman analysis biases were small (1.51 mL, 1.69 mL, -0.02%, -0.66 g for ESV, EDV, EF and LVM, respectively).

CONCLUSIONS

The proposed technique automatically and rapidly detects endocardial, epicardial, papillary muscles' and trabeculations' contours providing accurate and reproducible quantitative MRI parameters, including LV mass and EF.

Left ventricular segmentation from MRI datasets with edge modelling conditional random fields

BACKGROUND

This paper considers automatic segmentation of the left cardiac ventricle in short axis magnetic resonance images. Various aspects, such as the presence of papillary muscles near the endocardium border, makes simple threshold based segmentation difficult.

METHODS

The endo- and epicardium are modelled as two series of radii which are inter-related using features describing shape and motion. Image features are derived from edge information from human annotated images. The features are combined within a discriminatively trained Conditional Random Field (CRF). Loopy belief propagation is used to infer segmentations when an unsegmented video sequence is given. Powell's method is applied to find CRF parameters by minimizing the difference between ground truth annotations and the inferred contours. We also describe how the endocardium centre points are calculated from a single human-provided centre point in the first frame, through minimization of frame alignment error.

RESULTS

We present and analyse the results of segmentation. The algorithm exhibits robustness against inclusion of the papillary muscles by integrating shape and motion information. Possible future improvements are identified.

CONCLUSIONS

The presented model integrates shape and motion information to segment the inner and outer contours in the presence of papillary muscles. On the Sunnybrook dataset we find an average Dice metric of 0.91 ± 0.02 and 0.93 ± 0.02 for the inner and outer segmentations, respectively. Particularly problematic are patients with hypertrophy where the blood pool disappears from view at end-systole.

A magnetic resonance software simulator for the evaluation of myocardial deformation estimation

This paper proposes a methodology to design a physiologically realistic computer simulator of images of the left ventricle myocardium based on a patient-specific biomechanical model. The simulator takes a magnetic resonance image of a given patient at end diastole, uses a manual segmentation of that image to model the geometry of the myocardium and sets the parameters of the constitutive model used for biomechanical simulation according to a regional labeling of the contractility of the myocardium for that patient. The simulated deformations are used to warp the magnetic resonance dataset throughout the cardiac cycle to generate different image modalities. The simulator is validated by quantifying its ability to model actual deformations in a set of patients affected by an acute myocardial infarction. Specifically a high correlation has been encountered between the ejection fraction derived from the simulated end systolic deformation of the myocardium and the myocardium segmented from actual data. Additionally, most of the parameters that describe the simulated deformation compare well with reported values. Overall, the simulator is intended as a testbed for extensive comparisons of myocardial motion tracking methods due to its ability to relate the impaired myocardial function with the associated ventricular remodeling, a novel contribution in the literature of cardiac image simulators.

Robust estimation of MRI myocardial perfusion parameters

Myocardial perfusion imaging by first-pass contrast-enhanced magnetic resonance allows to asses the viability of a tissue by the study of the contrast agent transit through the cardiac chambers and myocardium. Since visual inspection is difficult and may left aside critical temporal information, the need of automatic quantitative analysis arises. We propose two robust estimators of the parameters that quantify the perfusion according to a prior pharmacokinetic model. The estimators are based on the concentration of the contrast agent inside the tissue and the blood.

4D statistical shape modeling of the left ventricle in cardiac MR images

PURPOSE

Statistical shape models have shown improved reliability and consistency in cardiac image segmentation. They incorporate a sufficient amount of a priori knowledge from the training datasets and solve some major problems such as noise and image artifacts or partial volume effect. In this paper, we construct a 4D statistical model of the left ventricle using human cardiac short-axis MR images.

METHODS

Kernel PCA is utilized to explore the nonlinear variation of a population. The distribution of the landmarks is divided into the inter- and intra-subject subspaces. We compare the result of Kernel PCA with linear PCA and ICA for each of these subspaces. The initial atlas in natural coordinate system is built for the end-diastolic frame. The landmarks extracted from it are propagated to all frames of all datasets. We apply the 4D KPCA-based ASM for segmentation of all phases of a cardiac cycle and compare it with the conventional ASM.

RESULTS

The proposed statistical model is evaluated by calculating the compactness capacity, specificity and generalization ability measures. We investigate the behavior of the nonlinear model for different values of the kernel parameter. The results show that the model built by KPCA is less compact than PCA but more compact than ICA. Although for a constant number of modes the reconstruction error is a little higher for the KPCA-based statistical model, it produces a statistical model with substantially better specificity than PCA- and ICA-based models.

CONCLUSION

Quantitative analysis of the results demonstrates that our method improves the segmentation accuracy.