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 is chemic 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

Synthetic generation of myocardial blood-oxygen-level-dependent MRI time series via structural sparse decomposition modeling

This paper aims to identify approaches that generate appropriate synthetic data (computer generated) for cardiac phase-resolved blood-oxygen-level-dependent (CP-BOLD) MRI. CP-BOLD MRI is a new contrast agent- and stress-free approach for examining changes in myocardial oxygenation in response to coronary artery disease. However, since signal intensity changes are subtle, rapid visualization is not possible with the naked eye. Quantifying and visualizing the extent of disease relies on myocardial segmentation and registration to isolate the myocardium and establish temporal correspondences and ischemia detection algorithms to identify temporal differences in BOLD signal intensity patterns. If transmurality of the defect is of interest pixel-level analysis is necessary and thus a higher precision in registration is required. Such precision is currently not available affecting the design and performance of the ischemia detection algorithms. In this work, to enable algorithmic developments of ischemia detection irrespective to registration accuracy, we propose an approach that generates synthetic pixel-level myocardial time series. We do this by 1) modeling the temporal changes in BOLD signal intensity based on sparse multi-component dictionary learning, whereby segmentally derived myocardial time series are extracted from canine experimental data to learn the model; and 2) demonstrating the resemblance between real and synthetic time series for validation purposes. We envision that the proposed approach has the capacity to accelerate development of tools for ischemia detection while markedly reducing experimental costs so that cardiac BOLD MRI can be rapidly translated into the clinical arena for the noninvasive assessment of ischemic heart disease

Dictionary Learning Based Image Descriptor for Myocardial Registration of CP-BOLD MR:Medical Image Computing and Computer-Assisted Intervention -- MICCAI 2015

Cardiac Phase-resolved Blood Oxygen-Level-Dependent (CP-BOLD) MRI is a new contrast agent- and stress-free imaging technique for the assessment of myocardial ischemia at rest. The precise registration among the cardiac phases in this cine type acquisition is essential for automating the analysis of images of this technique, since it can potentially lead to better specificity of ischemia detection. However, inconsistency in myocardial intensity patterns and the changes in myocardial shape due to the heart’s motion lead to low registration performance for state-of-the-art methods. This low accuracy can be explained by the lack of distinguishable features in CP-BOLD and inappropriate metric definitions in current intensity-based registration frameworks. In this paper, the sparse representations, which are defined by a discriminative dictionary learning approach for source and target images, are used to improve myocardial registration. This method combines appearance with Gabor and HOG features in a dictionary learning framework to sparsely represent features in a low dimensional space. The sum of absolute differences of these distinctive sparse representations are used to define a similarity term in the registration framework. The proposed approach is validated on a dataset of CP-BOLD MR and standard CINE MR acquired in baseline and ischemic condition across 10 canines

Fast dynamic magnetic resonance imaging using sparse recovery methods and novel signal encoding formulations

Magnetic resonance imaging (MRI) is a non-invasive imaging modality that provides excellent soft tissue contrast without using ionizing radiations. These qualities/properties make MRI the preferred imaging modality for critical organs like heart and brain. Over the past decade, the advancement in hardware and image reconstruction algorithms has led to substantial improvements in MRI in terms of imaging speeds, quality and reliability. However, MRI speeds need to be further improved while retaining/maintaining the image quality given that the emerging medical diagnostic procedures are increasingly relying on detailed characterization of physiological functions that evolve on time scales too fast to be captured using conventional MRI methods. This dissertation starts with presenting a sparse signal recovery based fast MRI method. This method synergistically combines a data redundancy scheme for high frequency details with a novel and physically realizable MR signal encoding formulation. The new signal encoding formulation uses clinically deployed tagging radio frequency pulses to mix information in the spatial frequency domain prior to acquisition. Thus, the new formulation leads to a more uniform coverage of spatial frequency information even at high accelerations. The synergistic combination of image-detail redundancy encoding with tagging based signal encoding allows recovery of edges and fine structures with unprecedented quality. Next, this dissertation evaluates the use of fast spiral trajectories for high spatial resolution functional imaging of human superior colliculus. Gradient efficient and motion-robust spiral trajectories are used to keep fMRI scan durations short. . However, high resolution imaging of human subcortical structures using these trajectories is limited due to the weak functional responses of SC structures and also low signal-to-noise ratio associated with small voxels. To improve the functional sensitivity of spiral trajectories, dual echo variants are used. Combination of two echoes of the dual-echo variants reduces noise and thereby improves the functional sensitivity of high resolution fMRI. Lastly, this dissertation presents a novel formulation for fast dynamic MRI which combines the generic linear dynamical system model with sparse recovery techniques. Specifically, the formulation uses a known prior spatio-temporal model to predict the underlying image and uses sparse recovery techniques to recover the residual image. The spatio-temporal evolution model inherently encodes for coupled data redundancies in the spatial- and temporal-dimensions. Also, the generalizability of the formulation in choosing the evolution model allows it to be applicable to various physiological functions.Electrical and Computer Engineerin

Cardiac motion estimation in ultrasound images using a sparse representation and dictionary learning

Les maladies cardiovasculaires sont de nos jours un problème de santé majeur. L'amélioration des méthodes liées au diagnostic de ces maladies représente donc un réel enjeu en cardiologie. Le coeur étant un organe en perpétuel mouvement, l'analyse du mouvement cardiaque est un élément clé pour le diagnostic. Par conséquent, les méthodes dédiées à l'estimation du mouvement cardiaque à partir d'images médicales, plus particulièrement en échocardiographie, font l'objet de nombreux travaux de recherches. Cependant, plusieurs difficultés liées à la complexité du mouvement du coeur ainsi qu'à la qualité des images échographiques restent à surmonter afin d'améliorer la qualité et la précision des estimations. Dans le domaine du traitement d'images, les méthodes basées sur l'apprentissage suscitent de plus en plus d'intérêt. Plus particulièrement, les représentations parcimonieuses et l'apprentissage de dictionnaires ont démontré leur efficacité pour la régularisation de divers problèmes inverses. Cette thèse a ainsi pour but d'explorer l'apport de ces méthodes, qui allient parcimonie et apprentissage, pour l'estimation du mouvement cardiaque. Trois principales contributions sont présentées, chacune traitant différents aspects et problématiques rencontrées dans le cadre de l'estimation du mouvement en échocardiographie. Dans un premier temps, une méthode d'estimation du mouvement cardiaque se basant sur une régularisation parcimonieuse est proposée. Le problème d'estimation du mouvement est formulé dans le cadre d'une minimisation d'énergie, dont le terme d'attache aux données est construit avec l'hypothèse d'un bruit de Rayleigh multiplicatif. Une étape d'apprentissage de dictionnaire permet une régularisation exploitant les propriétés parcimonieuses du mouvement cardiaque, combinée à un terme classique de lissage spatial. Dans un second temps, une méthode robuste de flux optique est présentée. L'objectif de cette approche est de robustifier la méthode d'estimation développée au premier chapitre de manière à la rendre moins sensible aux éléments aberrants. Deux régularisations sont mises en oeuvre, imposant d'une part un lissage spatial et de l'autre la parcimonie des champs de mouvements dans un dictionnaire approprié. Afin d'assurer la robustesse de la méthode vis-à-vis des anomalies, une stratégie de minimisation récursivement pondérée est proposée. Plus précisément, les fonctions employées pour cette pondération sont basées sur la théorie des M-estimateurs. Le dernier travail présenté dans cette thèse, explore une méthode d'estimation du mouvement cardiaque exploitant une régularisation parcimonieuse combinée à un lissage à la fois dans les domaines spatial et temporel. Le problème est formulé dans un cadre général d'estimation de flux optique. La régularisation temporelle proposée impose des trajectoires de mouvement lisses entre images consécutives. De plus, une méthode itérative d'estimation permet d'incorporer les trois termes de régularisations, tout en rendant possible le traitement simultané d'un ensemble d'images. Dans cette thèse, les contributions proposées sont validées en employant des images synthétiques et des simulations réalistes d'images ultrasonores. Ces données avec vérité terrain permettent d'évaluer la précision des approches considérées, et de souligner leur compétitivité par rapport à des méthodes de l'état-del'art. Pour démontrer la faisabilité clinique, des images in vivo de patients sains ou atteints de pathologies sont également considérées pour les deux premières méthodes. Pour la dernière contribution de cette thèse, i.e., exploitant un lissage temporel, une étude préliminaire est menée en utilisant des données de simulation.Cardiovascular diseases have become a major healthcare issue. Improving the diagnosis and analysis of these diseases have thus become a primary concern in cardiology. The heart is a moving organ that undergoes complex deformations. Therefore, the quantification of cardiac motion from medical images, particularly ultrasound, is a key part of the techniques used for diagnosis in clinical practice. Thus, significant research efforts have been directed toward developing new cardiac motion estimation methods. These methods aim at improving the quality and accuracy of the estimated motions. However, they are still facing many challenges due to the complexity of cardiac motion and the quality of ultrasound images. Recently, learning-based techniques have received a growing interest in the field of image processing. More specifically, sparse representations and dictionary learning strategies have shown their efficiency in regularizing different ill-posed inverse problems. This thesis investigates the benefits that such sparsity and learning-based techniques can bring to cardiac motion estimation. Three main contributions are presented, investigating different aspects and challenges that arise in echocardiography. Firstly, a method for cardiac motion estimation using a sparsity-based regularization is introduced. The motion estimation problem is formulated as an energy minimization, whose data fidelity term is built using the assumption that the images are corrupted by multiplicative Rayleigh noise. In addition to a classical spatial smoothness constraint, the proposed method exploits the sparse properties of the cardiac motion to regularize the solution via an appropriate dictionary learning step. Secondly, a fully robust optical flow method is proposed. The aim of this work is to take into account the limitations of ultrasound imaging and the violations of the regularization constraints. In this work, two regularization terms imposing spatial smoothness and sparsity of the motion field in an appropriate cardiac motion dictionary are also exploited. In order to ensure robustness to outliers, an iteratively re-weighted minimization strategy is proposed using weighting functions based on M-estimators. As a last contribution, we investigate a cardiac motion estimation method using a combination of sparse, spatial and temporal regularizations. The problem is formulated within a general optical flow framework. The proposed temporal regularization enforces smoothness of the motion trajectories between consecutive images. Furthermore, an iterative groupewise motion estimation allows us to incorporate the three regularization terms, while enabling the processing of the image sequence as a whole. Throughout this thesis, the proposed contributions are validated using synthetic and realistic simulated cardiac ultrasound images. These datasets with available groundtruth are used to evaluate the accuracy of the proposed approaches and show their competitiveness with state-of-the-art algorithms. In order to demonstrate clinical feasibility, in vivo sequences of healthy and pathological subjects are considered for the first two methods. A preliminary investigation is conducted for the last contribution, i.e., exploiting temporal smoothness, using simulated data

Learning-based Algorithms for Inverse Problems in MR Image Reconstruction and Quantitative Perfusion Imaging

Medical imaging has become an integral part of the clinical pipeline through its widespread use in the diagnosis, prognosis and treatment planning of several diseases. Magnetic Resonance Imaging (MRI) is particularly useful because it is free from ionizing radiation and is able to provide excellent soft tissue contrast. However, MRI suffers from drawbacks like long scanning durations that increase the cost of imaging and render the acquired images vulnerable to artifacts like motion. In modalities like Arterial Spin Labeling (ASL), which is used for non-invasive and quantitative perfusion imaging, low signal-to-noise ratio and lack of precision in parameter estimates also present significant problems. In this thesis, we develop and present algorithms whose focus can be divided into two broad categories. First, we investigate the reconstruction of MR images from fewer measurements, using data-driven machine learning to fill in the gaps in acquisition, thereby reducing the scan duration. Specifically, we first combine a supervised and an unsupervised (blind) learned dictionary in a residual fashion as a spatial prior in MR image reconstruction, and then extend this framework to include deep supervised learning. The latter, called blind primed supervised (BLIPS) learning, proposes that there exists synergy between features learned using shallower dictionary-based methods or traditional prior-based image reconstruction and those learned using newer deep supervised learning-based approaches. We show that this synergy can be exploited to yield reconstructions that are approx. 0.5-1 dB better in PSNR (in avg. across undersampling patterns). We also observe that the BLIPS algorithm is more robust to a scarcity of available training data, yielding reconstructions that are approx. 0.8 dB better (in terms of avg. PSNR) compared to strict supervised learning reconstruction when training data is very limited. Secondly, we aim to provide more precise estimates for multiple physiological parameters and tissue properties from ASL scans by estimation-theory-based optimization of ASL scan design, and combination with MR Fingerprinting. For this purpose, we use the Cramer-Rao Lower Bound (CRLB) for optimizing the scan design, and deep learning for regression-based estimation. We also show that regardless of the estimator used, optimization improves the precision in parameter estimates, and enables us to increase the available ‘useful’ information obtained in a fixed scanning duration. Specifically, we successfully improve the theoretical precision of perfusion estimates by 4.6% compared to a scan design where the repetition times are randomly chosen (a popular choice in literature) thereby yielding a 35.2% improvement in the corresponding RMSE in our in-silico experiments. This improvement is also visually evident in our in-vivo studies on healthy human subjects.PHDElectrical and Computer EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/169819/1/anishl_1.pd

Reconstruction algorithms for Magnetic Resonance Imaging

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science, 2013.Cataloged from PDF version of thesis.Includes bibliographical references (p. 135-142).This dissertation presents image reconstruction algorithms for Magnetic Resonance Imaging (MRI) that aims to increase the imaging efficiency. Algorithms that reduce imaging time without sacrificing the image quality and mitigate image artifacts are proposed. The goal of increasing the MR efficiency is investigated across multiple imaging techniques: structural imaging with multiple contrasts preparations, Diffusion Spectrum Imaging (DSI), Chemical Shift Imaging (CSI), and Quantitative Susceptibility Mapping (QSM). The main theme connecting the proposed methods is the utilization of prior knowledge on the reconstructed signal. This prior often presents itself in the form of sparsity with respect to either a prespecified or learned signal transformation.by Berkin Bilgic.Ph.D

Deep Learning in Medical Image Analysis

The accelerating power of deep learning in diagnosing diseases will empower physicians and speed up decision making in clinical environments. Applications of modern medical instruments and digitalization of medical care have generated enormous amounts of medical images in recent years. In this big data arena, new deep learning methods and computational models for efficient data processing, analysis, and modeling of the generated data are crucially important for clinical applications and understanding the underlying biological process. This book presents and highlights novel algorithms, architectures, techniques, and applications of deep learning for medical image analysis