Variational segmentation framework in prolate spheroidal coordinates for 3D real-time echocardiography

This paper presents a new formulation of a deformable model segmentation in prolate spheroidal coordinates for segmentation of 3D cardiac echocardiography data. The prolate spheroidal coordinate system enables a representation of the segmented surface with descriptors specifically adapted to the "ellipsoidal" shape of the ventricle. A simple data energy term, based on gray-level information, guides the segmentation. The segmentation framework provides a very fast and simple algorithm to evolve an initial ellipsoidal object towards the endocardial surface of the myocardium with near real-time deformations. With near real-time performance, additional constraints on landmark points, can be used interactively to prevent leakage of the surface

Augmenting CT cardiac roadmaps with segmented streaming ultrasound

Static X-ray computed tomography (CT) volumes are often used as anatomic roadmaps during catheter-based cardiac interventions performed under X-ray fluoroscopy guidance. These CT volumes provide a high-resolution depiction of soft-tissue structures, but at only a single point within the cardiac and respiratory cycles. Augmenting these static CT roadmaps with segmented myocardial borders extracted from live ultrasound (US) provides intra-operative access to real-time dynamic information about the cardiac anatomy. In this work, using a customized segmentation method based on a 3D active mesh, endocardial borders of the left ventricle were extracted from US image streams (4D data sets) at a frame rate of approximately 5 frames per second. The coordinate systems for CT and US modalities were registered using rigid body registration based on manually selected landmarks, and the segmented endocardial surfaces were overlaid onto the CT volume. The root-mean squared fiducial registration error was 3.80 mm. The accuracy of the segmentation was quantitatively evaluated in phantom and human volunteer studies via comparison with manual tracings on 9 randomly selected frames using a finite-element model (the US image resolutions of the phantom and volunteer data were 1.3 x 1.1 x 1.3 mm and 0.70 x 0.82 x 0.77 mm, respectively). This comparison yielded 3.70±2.5 mm (approximately 3 pixels) root-mean squared error (RMSE) in a phantom study and 2.58±1.58 mm (approximately 3 pixels) RMSE in a clinical study. The combination of static anatomical roadmap volumes and dynamic intra-operative anatomic information will enable better guidance and feedback for image-guided minimally invasive cardiac interventions

Feature Tracking Cardiac Magnetic Resonance via Deep Learning and Spline Optimization

Feature tracking Cardiac Magnetic Resonance (CMR) has recently emerged as an area of interest for quantification of regional cardiac function from balanced, steady state free precession (SSFP) cine sequences. However, currently available techniques lack full automation, limiting reproducibility. We propose a fully automated technique whereby a CMR image sequence is first segmented with a deep, fully convolutional neural network (CNN) architecture, and quadratic basis splines are fitted simultaneously across all cardiac frames using least squares optimization. Experiments are performed using data from 42 patients with hypertrophic cardiomyopathy (HCM) and 21 healthy control subjects. In terms of segmentation, we compared state-of-the-art CNN frameworks, U-Net and dilated convolution architectures, with and without temporal context, using cross validation with three folds. Performance relative to expert manual segmentation was similar across all networks: pixel accuracy was ~97%, intersection-over-union (IoU) across all classes was ~87%, and IoU across foreground classes only was ~85%. Endocardial left ventricular circumferential strain calculated from the proposed pipeline was significantly different in control and disease subjects (-25.3% vs -29.1%, p = 0.006), in agreement with the current clinical literature.Comment: Accepted to Functional Imaging and Modeling of the Heart (FIMH) 201

State of the Art of Level Set Methods in Segmentation and Registration of Medical Imaging Modalities

Segmentation of medical images is an important step in various applications such as visualization, quantitative analysis and image-guided surgery. Numerous segmentation methods have been developed in the past two decades for extraction of organ contours on medical images. Low-level segmentation methods, such as pixel-based clustering, region growing, and filter-based edge detection, require additional pre-processing and post-processing as well as considerable amounts of expert intervention or information of the objects of interest. Furthermore the subsequent analysis of segmented objects is hampered by the primitive, pixel or voxel level representations from those region-based segmentation. Deformable models, on the other hand, provide an explicit representation of the boundary and the shape of the object. They combine several desirable features such as inherent connectivity and smoothness, which counteract noise and boundary irregularities, as well as the ability to incorporate knowledge about the object of interest. However, parametric deformable models have two main limitations. First, in situations where the initial model and desired object boundary differ greatly in size and shape, the model must be re-parameterized dynamically to faithfully recover the object boundary. The second limitation is that it has difficulty dealing with topological adaptation such as splitting or merging model parts, a useful property for recovering either multiple objects or objects with unknown topology. This difficulty is caused by the fact that a new parameterization must be constructed whenever topology change occurs, which requires sophisticated schemes. Level set deformable models, also referred to as geometric deformable models, provide an elegant solution to address the primary limitations of parametric deformable models. These methods have drawn a great deal of attention since their introduction in 1988. Advantages of the contour implicit formulation of the deformable model over parametric formulation include: (1) no parameterization of the contour, (2) topological flexibility, (3) good numerical stability, (4) straightforward extension of the 2D formulation to n-D. Recent reviews on the subject include papers from Suri. In this chapter we give a general overview of the level set segmentation methods with emphasize on new frameworks recently introduced in the context of medical imaging problems. We then introduce novel approaches that aim at combining segmentation and registration in a level set formulation. Finally we review a selective set of clinical works with detailed validation of the level set methods for several clinical applications

QUANTIFICATION OF MYOCARDIAL MECHANICS IN LEFT VENTRICLES UNDER INOTROPIC STIMULATION AND IN HEALTHY RIGHT VENTRICLES USING 3D DENSE CMR

Statistical data from clinical studies indicate that the death rate caused by heart disease has decreased due to an increased use of evidence-based medical therapies. This includes the use of magnetic resonance imaging (MRI), which is one of the most common non-invasive approaches in evidence-based health care research. In the current work, I present 3D Lagrangian strains and torsion in the left ventricle of healthy and isoproterenol-stimulated rats, which were investigated using Displacement ENcoding with Stimulated Echoes (DENSE) cardiac magnetic resonance (CMR) imaging. With the implementation of the 12-segment model, a detailed profile of regional cardiac mechanics was reconstructed for each subject. Statistical analysis revealed that isoproterenol induced a significant change in the strains and torsion in certain regions at the mid-ventricle level. In addition, I investigated right ventricular cardiac mechanics with the methodologies developed for the left ventricle. This included a comparison of different regions within the basal and mid-ventricular regions. Despite no regional variation found in the peak circumferential strain, the peak longitudinal strain exhibited regional variation at the anterior side of the RV due to the differences in biventricular torsion, mechanism of RV free wall contraction, and fiber architecture at RV insertions. Future applications of the experimental work presented here include the construction and validation of biventricular finite element models. Specifically, the strains predicted by the models will be statistically compared with experimental strains. In addition, the results of the present study provide an essential reference of RV baseline evaluated with DENSE MRI, a highly objective technique

Foetal echocardiographic segmentation

Congenital heart disease affects just under one percentage of all live births [1]. Those defects that manifest themselves as changes to the cardiac chamber volumes are the motivation for the research presented in this thesis. Blood volume measurements in vivo require delineation of the cardiac chambers and manual tracing of foetal cardiac chambers is very time consuming and operator dependent. This thesis presents a multi region based level set snake deformable model applied in both 2D and 3D which can automatically adapt to some extent towards ultrasound noise such as attenuation, speckle and partial occlusion artefacts. The algorithm presented is named Mumford Shah Sarti Collision Detection (MSSCD). The level set methods presented in this thesis have an optional shape prior term for constraining the segmentation by a template registered to the image in the presence of shadowing and heavy noise. When applied to real data in the absence of the template the MSSCD algorithm is initialised from seed primitives placed at the centre of each cardiac chamber. The voxel statistics inside the chamber is determined before evolution. The MSSCD stops at open boundaries between two chambers as the two approaching level set fronts meet. This has significance when determining volumes for all cardiac compartments since cardiac indices assume that each chamber is treated in isolation. Comparison of the segmentation results from the implemented snakes including a previous level set method in the foetal cardiac literature show that in both 2D and 3D on both real and synthetic data, the MSSCD formulation is better suited to these types of data. All the algorithms tested in this thesis are within 2mm error to manually traced segmentation of the foetal cardiac datasets. This corresponds to less than 10% of the length of a foetal heart. In addition to comparison with manual tracings all the amorphous deformable model segmentations in this thesis are validated using a physical phantom. The volume estimation of the phantom by the MSSCD segmentation is to within 13% of the physically determined volume

Shape Analysis Based Strategies for Evaluation of Adaptations in In Vivo Right Ventricular Geometry and Mechanics as Effected by Pulmonary Hypertension

Pulmonary hypertension (PH) is a deadly disease, which as it progresses over time alters many aspects of the afflicted heart, and particularly the right ventricle (RV), such as its size, shape, and mechanical material properties. However, due to the limitations of what can be measured noninvasively in a standard clinical setting and the difficulty caused by the intrinsic complexity of the human RV, there has been little success to-date to identify clinically obtainable metrics of RV shape, deformation, or material properties that are quantitatively linked to the onset and progression of PH. Towards addressing this challenge, this work proposes the use of the shape and shape change of the RV, which is measurable from standard clinical imaging, along with statistical analysis and inverse material characterization strategies to identify new metrics of RV mechanical function that will be uniquely predictive of the state of the heart subject to PH. Thus, this thesis can be broken into two components: the first is statistical shape analysis of the RV, and the second is inverse characterization of heart wall mechanical material properties from RV shape change and measurable hemodynamics. For the statistical shape analysis investigation, a custom approach using harmonic mapping and proper orthogonal decomposition is applied to determine the fundamental components of shape (i.e., modes) from a dataset of 50 patients with varying states of PH, including some without PH at all. For the inverse characterization work, a novel method was developed to estimate the heterogeneous properties of a structure, given only the target shape of that structure, after a known excitation is applied to deform the structure. Lastly, the inverse characterization algorithm was extended to be applicable to actual in vivo cardiac data, particularly through the inclusion of a registration step to account for the organ-scale rotation and translation of the heart. Future work remains to expand on the computational efficiency of this inverse solution estimation procedure, and to further evaluate and improve upon the consistency and clinical interpretability of the material property estimates