Respiratory organ motion in interventional MRI : tracking, guiding and modeling

Respiratory organ motion is one of the major challenges in interventional MRI, particularly in interventions with therapeutic ultrasound in the abdominal region. High-intensity focused ultrasound found an application in interventional MRI for noninvasive treatments of different abnormalities. In order to guide surgical and treatment interventions, organ motion imaging and modeling is commonly required before a treatment start. Accurate tracking of organ motion during various interventional MRI procedures is prerequisite for a successful outcome and safe therapy. In this thesis, an attempt has been made to develop approaches using focused ultrasound which could be used in future clinically for the treatment of abdominal organs, such as the liver and the kidney. Two distinct methods have been presented with its ex vivo and in vivo treatment results. In the first method, an MR-based pencil-beam navigator has been used to track organ motion and provide the motion information for acoustic focal point steering, while in the second approach a hybrid imaging using both ultrasound and magnetic resonance imaging was combined for advanced guiding capabilities. Organ motion modeling and four-dimensional imaging of organ motion is increasingly required before the surgical interventions. However, due to the current safety limitations and hardware restrictions, the MR acquisition of a time-resolved sequence of volumetric images is not possible with high temporal and spatial resolution. A novel multislice acquisition scheme that is based on a two-dimensional navigator, instead of a commonly used pencil-beam navigator, was devised to acquire the data slices and the corresponding navigator simultaneously using a CAIPIRINHA parallel imaging method. The acquisition duration for four-dimensional dataset sampling is reduced compared to the existing approaches, while the image contrast and quality are improved as well. Tracking respiratory organ motion is required in interventional procedures and during MR imaging of moving organs. An MR-based navigator is commonly used, however, it is usually associated with image artifacts, such as signal voids. Spectrally selective navigators can come in handy in cases where the imaging organ is surrounding with an adipose tissue, because it can provide an indirect measure of organ motion. A novel spectrally selective navigator based on a crossed-pair navigator has been developed. Experiments show the advantages of the application of this novel navigator for the volumetric imaging of the liver in vivo, where this navigator was used to gate the gradient-recalled echo sequence

Deep learning-based fully automatic segmentation of wrist cartilage in MR images

The study objective was to investigate the performance of a dedicated convolutional neural network (CNN) optimized for wrist cartilage segmentation from 2D MR images. CNN utilized a planar architecture and patch-based (PB) training approach that ensured optimal performance in the presence of a limited amount of training data. The CNN was trained and validated in twenty multi-slice MRI datasets acquired with two different coils in eleven subjects (healthy volunteers and patients). The validation included a comparison with the alternative state-of-the-art CNN methods for the segmentation of joints from MR images and the ground-truth manual segmentation. When trained on the limited training data, the CNN outperformed significantly image-based and patch-based U-Net networks. Our PB-CNN also demonstrated a good agreement with manual segmentation (Sorensen-Dice similarity coefficient (DSC) = 0.81) in the representative (central coronal) slices with large amount of cartilage tissue. Reduced performance of the network for slices with a very limited amount of cartilage tissue suggests the need for fully 3D convolutional networks to provide uniform performance across the joint. The study also assessed inter- and intra-observer variability of the manual wrist cartilage segmentation (DSC=0.78-0.88 and 0.9, respectively). The proposed deep-learning-based segmentation of the wrist cartilage from MRI could facilitate research of novel imaging markers of wrist osteoarthritis to characterize its progression and response to therapy

Customizable tubular model for n-furcating blood vessels and its application to 3D reconstruction of the cerebrovascular system

Understanding the 3D cerebral vascular network is one of the pressing issues impacting the diagnostics of various systemic disorders and is helpful in clinical therapeutic strategies. Unfortunately, the existing software in the radiological workstation does not meet the expectations of radiologists who require a computerized system for detailed, quantitative analysis of the human cerebrovascular system in 3D and a standardized geometric description of its components. In this study, we show a method that uses 3D image data from magnetic resonance imaging with contrast to create a geometrical reconstruction of the vessels and a parametric description of the reconstructed segments of the vessels. First, the method isolates the vascular system using controlled morphological growing and performs skeleton extraction and optimization. Then, around the optimized skeleton branches, it creates tubular objects optimized for quality and accuracy of matching with the originally isolated vascular data. Finally, it optimizes the joints on n-furcating vessel segments. As a result, the algorithm gives a complete description of shape, position in space, position relative to other segments, and other anatomical structures of each cerebrovascular system segment. Our method is highly customizable and in principle allows reconstructing vascular structures from any 2D or 3D data. The algorithm solves shortcomings of currently available methods including failures to reconstruct the vessel mesh in the proximity of junctions and is free of mesh collisions in high curvature vessels. It also introduces a number of optimizations in the vessel skeletonization leading to a more smooth and more accurate model of the vessel network. We have tested the method on 20 datasets from the public magnetic resonance angiography image database and show that the method allows for repeatable and robust segmentation of the vessel network and allows to compute vascular lateralization indices. Graphical abstract: [Figure not available: see fulltext.]</p

Computerized Analysis of Magnetic Resonance Images to Study Cerebral Anatomy in Developing Neonates

The study of cerebral anatomy in developing neonates is of great importance for the understanding of brain development during the early period of life. This dissertation therefore focuses on three challenges in the modelling of cerebral anatomy in neonates during brain development. The methods that have been developed all use Magnetic Resonance Images (MRI) as source data. To facilitate study of vascular development in the neonatal period, a set of image analysis algorithms are developed to automatically extract and model cerebral vessel trees. The whole process consists of cerebral vessel tracking from automatically placed seed points, vessel tree generation, and vasculature registration and matching. These algorithms have been tested on clinical Time-of- Flight (TOF) MR angiographic datasets. To facilitate study of the neonatal cortex a complete cerebral cortex segmentation and reconstruction pipeline has been developed. Segmentation of the neonatal cortex is not effectively done by existing algorithms designed for the adult brain because the contrast between grey and white matter is reversed. This causes pixels containing tissue mixtures to be incorrectly labelled by conventional methods. The neonatal cortical segmentation method that has been developed is based on a novel expectation-maximization (EM) method with explicit correction for mislabelled partial volume voxels. Based on the resulting cortical segmentation, an implicit surface evolution technique is adopted for the reconstruction of the cortex in neonates. The performance of the method is investigated by performing a detailed landmark study. To facilitate study of cortical development, a cortical surface registration algorithm for aligning the cortical surface is developed. The method first inflates extracted cortical surfaces and then performs a non-rigid surface registration using free-form deformations (FFDs) to remove residual alignment. Validation experiments using data labelled by an expert observer demonstrate that the method can capture local changes and follow the growth of specific sulcus

Investigation of Flow Disturbances and Multi-Directional Wall Shear Stress in the Stenosed Carotid Artery Bifurcation Using Particle Image Velocimetry

Hemodynamics and shear forces are associated with pathological changes in the vascular wall and its function, resulting in the focal development of atherosclerosis. Flow complexities that develop in the presence of established plaques create environments favourable to thrombosis formation and potentially plaque rupture leading to stroke. The carotid artery bifurcation is a common site of atherosclerosis development. Recently, the multi-directional nature of shear stress acting on the endothelial layer has been highlighted as a risk factor for atherogenesis, emphasizing the need for accurate measurements of shear stress magnitude as well direction. In the absence of comprehensive patient specific datasets numerical simulations of hemodynamics are limited by modeling assumptions. The objective of this thesis was to investigate the relative contributions of various factors - including geometry, rheology, pulsatility, and compliance – towards the development of disturbed flow and multi-directional wall shear stress (WSS) parameters related to the development of atherosclerosis An experimental stereoscopic particle image velocimetry (PIV) system was used to measure instantaneous full-field velocity in idealized asymmetrically stenosed carotid artery bifurcation models, enabling the extraction of bulk flow features and turbulence intensity (TI). The velocity data was combined with wall location information segmented from micro computed tomography (CT) to obtain phase-averaged maps of WSS magnitude and direction. A comparison between Newtonian and non-Newtonian blood-analogue fluids demonstrated that the conventional Newtonian viscosity assumption underestimates WSS magnitude while overestimating TI. Studies incorporating varying waveform pulsatility demonstrated that the levels of TI and oscillatory shear index (OSI) depend on the waveform amplitude in addition to the degree of vessel constriction. Local compliance resulted in a dampening of disturbed flow due to volumetric capacity of the upstream vessel, however wall tracking had a negligible effect on WSS prediction. While the degree of stenosis severity was found to have a dominant effect on local hemodynamics, comparable relative differences in metrics of flow and WSS disturbances were found due to viscosity model, waveform pulsatility and local vessel compliance

Imaging Biomarkers of Pulmonary Structure and Function

Asthma and chronic obstructive pulmonary disease (COPD) are characterized by airflow limitations resulting from airway obstruction and/or tissue destruction. The diagnosis and monitoring of these pulmonary diseases is primarily performed using spirometry, specifically the forced expiratory volume in one second (FEV1), which measures global airflow obstruction and provides no regional information of the different underlying disease pathologies. The limitations of spirometry and current therapies for lung disease patients have motivated the development of pulmonary imaging approaches, such as computed tomography (CT) and magnetic resonance imaging (MRI). Inhaled hyperpolarized noble gas MRI, specifically using helium-3 (3He) and xenon-129 (129Xe) gases, provides a way to quantify pulmonary ventilation by visualizing lung regions accessed by gas during a breath-hold, and alternatively, regions that are not accessed - coined “ventilation defects.” Despite the strong foundation and many advantages hyperpolarized 3He MRI has to offer research and patient care, clinical translation has been inhibited in part due to the cost and need for specialized equipment, including multinuclear-MR hardware and polarizers, and personnel. Accordingly, our objective was to develop and evaluate imaging biomarkers of pulmonary structure and function using MRI and CT without the use of exogenous contrast agents or specialized equipment. First, we developed and compared CT parametric response maps (PRM) with 3He MR ventilation images in measuring gas-trapping and emphysema in ex-smokers with and without COPD. We observed that in mild-moderate COPD, 3He MR ventilation abnormalities were related to PRM gas-trapping whereas in severe COPD, ventilation abnormalities correlated with both PRM gas-trapping and PRM emphysema. We then developed and compared pulmonary ventilation abnormalities derived from Fourier decomposition of free-breathing proton (1H) MRI (FDMRI) with 3He MRI in subjects with COPD and bronchiectasis. This work demonstrated that FDMRI and 3He MRI ventilation defects were strongly related in COPD, but not in bronchiectasis subjects. In COPD only, FDMRI ventilation defects were spatially related with 3He MRI ventilation defects and emphysema. Based on the FDMRI biomarkers developed in patients with COPD and bronchiectasis, we then evaluated ventilation heterogeneity in patients with severe asthma, both pre- and post-salbutamol as well as post-methacholine challenge, using FDMRI and 3He MRI. FDMRI free-breathing ventilation abnormalities were correlated with but under-estimated 3He MRI static ventilation defects. Finally, based on the previously developed free-breathing MRI approach, we developed a whole-lung free-breathing pulmonary 1H MRI technique to measure regional specific-ventilation and evaluated both asthmatics and healthy volunteers. These measurements not only provided similar information as specific-ventilation measured using plethysmography, but also information about regional ventilation defects that were correlated with 3He MRI ventilation abnormalities. These results demonstrated that whole-lung free-breathing 1H MRI biomarker of specific-ventilation may reflect ventilation heterogeneity and/or gas-trapping in asthma. These important findings indicate that imaging biomarkers of pulmonary structure and function using MRI and CT have the potential to regionally reveal the different pathologies in COPD and asthma without the use of exogenous contrast agents. The development and validation of these clinically meaningful imaging biomarkers are critically required to accelerate pulmonary imaging translation from the research workbench to being a part of the clinical workflow, with the overall goal to improve patient outcomes

Cerebral atrophy in mild cognitive impairment and Alzheimer disease: rates and acceleration.

OBJECTIVE: To quantify the regional and global cerebral atrophy rates and assess acceleration rates in healthy controls, subjects with mild cognitive impairment (MCI), and subjects with mild Alzheimer disease (AD). METHODS: Using 0-, 6-, 12-, 18-, 24-, and 36-month MRI scans of controls and subjects with MCI and AD from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database, we calculated volume change of whole brain, hippocampus, and ventricles between all pairs of scans using the boundary shift integral. RESULTS: We found no evidence of acceleration in whole-brain atrophy rates in any group. There was evidence that hippocampal atrophy rates in MCI subjects accelerate by 0.22%/year2 on average (p = 0.037). There was evidence of acceleration in rates of ventricular enlargement in subjects with MCI (p = 0.001) and AD (p < 0.001), with rates estimated to increase by 0.27 mL/year2 (95% confidence interval 0.12, 0.43) and 0.88 mL/year2 (95% confidence interval 0.47, 1.29), respectively. A post hoc analysis suggested that the acceleration of hippocampal loss in MCI subjects was mainly driven by the MCI subjects that were observed to progress to clinical AD within 3 years of baseline, with this group showing hippocampal atrophy rate acceleration of 0.50%/year2 (p = 0.003). CONCLUSIONS: The small acceleration rates suggest a long period of transition to the pathologic losses seen in clinical AD. The acceleration in hippocampal atrophy rates in MCI subjects in the ADNI seems to be driven by those MCI subjects who concurrently progressed to a clinical diagnosis of AD

Use of radiobiological modeling in treatment plan evaluation and optimization of prostate cancer radiotherapy

There are many tools available that are used to evaluate a radiotherapy treatment plan, such as isodose distribution charts, dose volume histograms (DVH), maximum, minimum and mean doses of the dose distributions as well as DVH point dose constraints. All the already mentioned evaluation tools are dosimetric only without taking into account the radiobiological characteristics of tumors or OARs. It has been demonstrated that although competing treatment plans might have similar mean, maximum or minimum doses they may have significantly different clinical outcomes (Mavroidis et al. 2001). For performing a more complete treatment plan evaluation and comparison the complication-free tumor control probability (P+) and the biologically effective uniform dose (D ) have been proposed (Källman et al. 1992a, Mavroidis et al. 2000). The D concept denotes that any two dose distributions within a target or OAR are equivalent if they produce the same probability for tumor control or normal tissue complication, respectively (Mavroidis et al. 2001)..

Rapid Segmentation Techniques for Cardiac and Neuroimage Analysis

Recent technological advances in medical imaging have allowed for the quick acquisition of highly resolved data to aid in diagnosis and characterization of diseases or to guide interventions. In order to to be integrated into a clinical work flow, accurate and robust methods of analysis must be developed which manage this increase in data. Recent improvements in in- expensive commercially available graphics hardware and General-Purpose Programming on Graphics Processing Units (GPGPU) have allowed for many large scale data analysis problems to be addressed in meaningful time and will continue to as parallel computing technology improves. In this thesis we propose methods to tackle two clinically relevant image segmentation problems: a user-guided segmentation of myocardial scar from Late-Enhancement Magnetic Resonance Images (LE-MRI) and a multi-atlas segmentation pipeline to automatically segment and partition brain tissue from multi-channel MRI. Both methods are based on recent advances in computer vision, in particular max-flow optimization that aims at solving the segmentation problem in continuous space. This allows for (approximately) globally optimal solvers to be employed in multi-region segmentation problems, without the particular drawbacks of their discrete counterparts, graph cuts, which typically present with metrication artefacts. Max-flow solvers are generally able to produce robust results, but are known for being computationally expensive, especially with large datasets, such as volume images. Additionally, we propose two new deformable registration methods based on Gauss-Newton optimization and smooth the resulting deformation fields via total-variation regularization to guarantee the problem is mathematically well-posed. We compare the performance of these two methods against four highly ranked and well-known deformable registration methods on four publicly available databases and are able to demonstrate a highly accurate performance with low run times. The best performing variant is subsequently used in a multi-atlas segmentation pipeline for the segmentation of brain tissue and facilitates fast run times for this computationally expensive approach. All proposed methods are implemented using GPGPU for a substantial increase in computational performance and so facilitate deployment into clinical work flows. We evaluate all proposed algorithms in terms of run times, accuracy, repeatability and errors arising from user interactions and we demonstrate that these methods are able to outperform established methods. The presented approaches demonstrate high performance in comparison with established methods in terms of accuracy and repeatability while largely reducing run times due to the employment of GPU hardware

The Human Connectome Project's neuroimaging approach

Noninvasive human neuroimaging has yielded many discoveries about the brain. Numerous methodological advances have also occurred, though inertia has slowed their adoption. This paper presents an integrated approach to data acquisition, analysis and sharing that builds upon recent advances, particularly from the Human Connectome Project (HCP). The 'HCP-style' paradigm has seven core tenets: (i) collect multimodal imaging data from many subjects; (ii) acquire data at high spatial and temporal resolution; (iii) preprocess data to minimize distortions, blurring and temporal artifacts; (iv) represent data using the natural geometry of cortical and subcortical structures; (v) accurately align corresponding brain areas across subjects and studies; (vi) analyze data using neurobiologically accurate brain parcellations; and (vii) share published data via user-friendly databases. We illustrate the HCP-style paradigm using existing HCP data sets and provide guidance for future research. Widespread adoption of this paradigm should accelerate progress in understanding the brain in health and disease