Automatic segmentation of the left ventricle cavity and myocardium in MRI data

A novel approach for the automatic segmentation has been developed to extract the epi-cardium and endo-cardium boundaries of the left ventricle (lv) of the heart. The developed segmentation scheme takes multi-slice and multi-phase magnetic resonance (MR) images of the heart, transversing the short-axis length from the base to the apex. Each image is taken at one instance in the heart's phase. The images are segmented using a diffusion-based filter followed by an unsupervised clustering technique and the resulting labels are checked to locate the (lv) cavity. From cardiac anatomy, the closest pool of blood to the lv cavity is the right ventricle cavity. The wall between these two blood-pools (interventricular septum) is measured to give an approximate thickness for the myocardium. This value is used when a radial search is performed on a gradient image to find appropriate robust segments of the epi-cardium boundary. The robust edge segments are then joined using a normal spline curve. Experimental results are presented with very encouraging qualitative and quantitative results and a comparison is made against the state-of-the art level-sets method

Doctor of Philosophy

dissertationComputational simulation has become an indispensable tool in the study of both basic mechanisms and pathophysiology of all forms of cardiac electrical activity. Because the heart is comprised of approximately 4 billion electrically active cells, it is not possible to geometrically model or computationally simulate each individual cell. As a result computational models of the heart are, of necessity, abstractions that approximate electrical behavior at the cell, tissue, and whole body level. The goal of this PhD dissertation was to evaluate several aspects of these abstractions by exploring a set of modeling approaches in the field of cardiac electrophysiology and to develop means to evaluate both the amplitude of these errors from a purely technical perspective as well as the impacts of those errors in terms of physiological parameters. The first project used subject specific models and experiments with acute myocardial ischemia to show that one common simplification used to model myocardial ischemia-the simplest form of the border zone between healthy and ischemic tissue-was not supported by the experimental results. We propose a alternative approximation of the border zone that better simulates the experimental results. The second study examined the impact of simplifications in geometric models on simulations of cardiac electrophysiology. Such models consist of a connected mesh of polygonal elements and must often capture complex external and internal boundaries. A conforming mesh contains elements that follow closely the shapes of boundaries; nonconforming meshes fit the boundaries only approximately and are easier to construct but their impact on simulation accuracy has, to our knowledge, remained unknown. We evaluated the impact of this simplification on a set of three different forms of bioelectric field simulations. The third project evaluated the impact of an additional geometric modeling error; positional uncertainty of the heart in simulations of the ECG. We applied a relatively novel and highly efficient statistical approach, the generalized Polynomial Chaos-Stochastic Collocation method (gPC-SC), to a boundary element formulation of the electrocardiographic forward problem to carry out the necessary comprehensive sensitivity analysis. We found variations large enough to mask or to mimic signs of ischemia in the ECG

An Automatic Level Set Based Liver Segmentation from MRI Data Sets

A fast and accurate liver segmentation method is a challenging work in medical image analysis area. Liver segmentation is an important process for computer-assisted diagnosis, pre-evaluation of liver transplantation and therapy planning of liver tumors. There are several advantages of magnetic resonance imaging such as free form ionizing radiation and good contrast visualization of soft tissue. Also, innovations in recent technology and image acquisition techniques have made magnetic resonance imaging a major tool in modern medicine. However, the use of magnetic resonance images for liver segmentation has been slow when we compare applications with the central nervous systems and musculoskeletal. The reasons are irregular shape, size and position of the liver, contrast agent effects and similarities of the gray values of neighbor organs. Therefore, in this study, we present a fully automatic liver segmentation method by using an approximation of the level set based contour evolution from T2 weighted magnetic resonance data sets. The method avoids solving partial differential equations and applies only integer operations with a two-cycle segmentation algorithm. The efficiency of the proposed approach is achieved by applying the algorithm to all slices with a constant number of iteration and performing the contour evolution without any user defined initial contour. The obtained results are evaluated with four different similarity measures and they show that the automatic segmentation approach gives successful results

Multiframe Temporal Estimation of Cardiac Nonrigid Motion

A robust, flexible system for tracking the point to point nonrigid motion of the left ventricular (LV) endocardial wall in image sequences has been developed. This system is unique in its ability to model motion trajectories across multiple frames. The foundation of this system is an adaptive transversal filter based on the recursive least-squares algorithm. This filter facilitates the integration of models for periodicity and proximal smoothness as appropriate using a contour-based description of the object’s boundaries. A set of correspondences between contours and an associated set of correspondence quality measures comprise the input to the system. Frame-to-frame relationships from two different frames of reference are derived and analyzed using synthetic and actual images. Two multiframe temporal models, both based on a sum of sinusoids, are derived. Illustrative examples of the system’s output are presented for quantitative analysis. Validation of the system is performed by comparing computed trajectory estimates with the trajectories of physical markers implanted in the LV wall. Sample case studies of marker trajectory comparisons are presented. Ensemble statistics from comparisons with 15 marker trajectories are acquired and analyzed. A multiframe temporal model without spatial periodicity constraints was determined to provide excellent performance with the least computational cost. A multiframe spatiotemporal model provided the best performance based on statistical standard deviation, although at significant computational expense.National Heart, Lung, and Blood InstituteAir Force of Scientific ResearchNational Science FoundationOffice of Naval ResearchR01HL44803F49620-99-1-0481F49620-99-1-0067MIP-9615590N00014-98-1-054

Dependent Nonparametric Bayesian Group Dictionary Learning for online reconstruction of Dynamic MR images

In this paper, we introduce a dictionary learning based approach applied to the problem of real-time reconstruction of MR image sequences that are highly undersampled in k-space. Unlike traditional dictionary learning, our method integrates both global and patch-wise (local) sparsity information and incorporates some priori information into the reconstruction process. Moreover, we use a Dependent Hierarchical Beta-process as the prior for the group-based dictionary learning, which adaptively infers the dictionary size and the sparsity of each patch; and also ensures that similar patches are manifested in terms of similar dictionary atoms. An efficient numerical algorithm based on the alternating direction method of multipliers (ADMM) is also presented. Through extensive experimental results we show that our proposed method achieves superior reconstruction quality, compared to the other state-of-the- art DL-based methods

Virtual cardiac monolayers for electrical wave propagation

The complex structure of cardiac tissue is considered to be one of the main determinants of an arrhythmogenic substrate. This study is aimed at developing the first mathematical model to describe the formation of cardiac tissue, using a joint in silico-in vitro approach. First, we performed experiments under various conditions to carefully characterise the morphology of cardiac tissue in a culture of neonatal rat ventricular cells. We considered two cell types, namely, cardiomyocytes and fibroblasts. Next, we proposed a mathematical model, based on the Glazier-Graner-Hogeweg model, which is widely used in tissue growth studies. The resultant tissue morphology was coupled to the detailed electrophysiological Korhonen-Majumder model for neonatal rat ventricular cardiomyocytes, in order to study wave propagation. The simulated waves had the same anisotropy ratio and wavefront complexity as those in the experiment. Thus, we conclude that our approach allows us to reproduce the morphological and physiological properties of cardiac tissue

Non-invasive estimation of QLV from the standard 12-lead ECG in patients with left bundle branch block

Background: Cardiac resynchronization therapy (CRT) is a treatment for patients with heart failure and electrical dyssynchrony, i.e., left bundle branch block (LBBB) ECG pattern. CRT resynchronizes ventricular contraction with a right ventricle (RV) and a left ventricle (LV) pacemaker lead. Positioning the LV lead in the latest electrically activated region (measured from Q wave onset in the ECG to LV sensing by the left pacemaker electrode [QLV]) is associated with favorable outcome. However, optimal LV lead placement is limited by coronary venous anatomy and the inability to measure QLV non-invasively before implantation. We propose a novel non-invasive method for estimating QLV in sinus-rhythm from the standard 12-lead ECG. Methods: We obtained 12-lead ECG, LV electrograms and LV lead position in a standard LV 17-segment model from procedural recordings from 135 standard CRT recipients. QLV duration was measured post-operatively. Using a generic heart geometry and corresponding forward model for ECG computation, the electrical activation pattern of the heart was fitted to best match the 12-lead ECG in an iterative optimization procedure. This procedure initialized six activation sites associated with the His-Purkinje system. The initial timing of each site was based on the directions of the vectorcardiogram (VCG). Timing and position of the sites were then changed iteratively to improve the match between simulated and measured ECG. Noninvasive estimation of QLV was done by calculating the time difference between Q-onset on the computed ECG and the activation time corresponding to centroidal epicardial activation time of the segment where the LV electrode is positioned. The estimated QLV was compared to the measured QLV. Further, the distance between the actual LV position and the estimated LV position was computed from the generic ventricular model. Results: On average there was no difference between QLV measured from procedural recordings and non-invasive estimation of QLV ( [Formula: see text] ). Median distance between actual LV pacing site and the estimated pacing site was 18.6 mm (IQR 17.3 mm). Conclusion: Using the standard 12-lead ECG and a generic heart model it is possible to accurately estimate QLV. This method may potentially be used to support patient selection, optimize implant procedures, and to simulate optimal stimulation parameters prior to pacemaker implantation

Sensitivity and specificity of substrate mapping: An in silico framework for the evaluation of electroanatomical substrate mapping strategies

pre-printBackground - Voltage mapping is an important tool for characterizing proarrhythmic electrophysiological substrate, yet it is subject to geometric factors that influence bipolar amplitudes and thus compromise performance. The aim of this study was to characterize the impact of catheter orientation on the ability of bipolar amplitudes to accurately discriminate between healthy and diseased tissues. Methods and Results - We constructed a three-dimensional, in-silico, bidomain model of cardiac tissue containing transmural lesions of varying diameter. A planar excitation wave was stimulated and electrograms were sampled with a realistic catheter model at multiple positions and orientations. We carried out validation studies in animal experiments of acute ablation lesions mapped with a clinical mapping system. Bipolar electrograms sampled at higher inclination angles of the catheter with respect to the tissue demonstrated improvements in both sensitivity and specificity of lesion detection. Removing low voltage electrograms with concurrent activation of both electrodes, suggesting false attenuation of the bipolar electrogram due to alignment with the excitation wavefront, had little effect on the accuracy of voltage mapping. Conclusions - Our results demonstrate possible mechanisms for the impact of catheter orientation on voltage mapping accuracy. Moreover, results from our simulations suggest that mapping accuracy may be improved by selectively controlling the inclination of the catheter to record at higher angles with respect to the tissue