Nonparametric tests of structure for high angular resolution diffusion imaging in Q-space

High angular resolution diffusion imaging data is the observed characteristic function for the local diffusion of water molecules in tissue. This data is used to infer structural information in brain imaging. Nonparametric scalar measures are proposed to summarize such data, and to locally characterize spatial features of the diffusion probability density function (PDF), relying on the geometry of the characteristic function. Summary statistics are defined so that their distributions are, to first-order, both independent of nuisance parameters and also analytically tractable. The dominant direction of the diffusion at a spatial location (voxel) is determined, and a new set of axes are introduced in Fourier space. Variation quantified in these axes determines the local spatial properties of the diffusion density. Nonparametric hypothesis tests for determining whether the diffusion is unimodal, isotropic or multi-modal are proposed. More subtle characteristics of white-matter microstructure, such as the degree of anisotropy of the PDF and symmetry compared with a variety of asymmetric PDF alternatives, may be ascertained directly in the Fourier domain without parametric assumptions on the form of the diffusion PDF. We simulate a set of diffusion processes and characterize their local properties using the newly introduced summaries. We show how complex white-matter structures across multiple voxels exhibit clear ellipsoidal and asymmetric structure in simulation, and assess the performance of the statistics in clinically-acquired magnetic resonance imaging data.Comment: Published in at http://dx.doi.org/10.1214/10-AOAS441 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Optimizing Magnetic Resonance Imaging for Image-Guided Radiotherapy

Magnetic resonance imaging (MRI) is playing an increasingly important role in image-guided radiotherapy. MRI provides excellent soft tissue contrast, and is flexible in characterizing various tissue properties including relaxation, diffusion and perfusion. This thesis aims at developing new image analysis and reconstruction algorithms to optimize MRI in support of treatment planning, target delineation and treatment response assessment for radiotherapy. First, unlike Computed Tomography (CT) images, MRI cannot provide electron density information necessary for radiation dose calculation. To address this, we developed a synthetic CT generation algorithm that generates pseudo CT images from MRI, based on tissue classification results on MRI for female pelvic patients. To improve tissue classification accuracy, we learnt a pelvic bone shape model from a training dataset, and integrated the shape model into an intensity-based fuzzy c-menas classification scheme. The shape-regularized tissue classification algorithm is capable of differentiating tissues that have significant overlap in MRI intensity distributions. Treatment planning dose calculations using synthetic CT image volumes generated from the tissue classification results show acceptably small variations as compared to CT volumes. As MRI artifacts, such as B1 filed inhomogeneity (bias field) may negatively impact the tissue classification accuracy, we also developed an algorithm that integrates the correction of bias field into the tissue classification scheme. We modified the fuzzy c-means classification by modeling the image intensity as the true intensity corrupted by the multiplicative bias field. A regularization term further ensures the smoothness of the bias field. We solved the optimization problem using a linearized alternating direction method of multipliers (ADMM) method, which is more computational efficient over existing methods. The second part of this thesis looks at a special MR imaging technique, diffusion-weighted MRI (DWI). By acquiring a series of DWI images with a wide range of b-values, high order diffusion analysis can be performed using the DWI image series and new biomarkers for tumor grading, delineation and treatment response evaluation may be extracted. However, DWI suffers from low signal-to-noise ratio at high b-values, and the multi-b-value acquisition makes the total scan time impractical for clinical use. In this thesis, we proposed an accelerated DWI scheme, that sparsely samples k-space and reconstructs images using a model-based algorithm. Specifically, we built a 3D block-Hankel tensor from k-space samples, and modeled both local and global correlations of the high dimensional k-space data as a low-rank property of the tensor. We also added a phase constraint to account for large phase variations across different b-values, and to allow reconstruction from partial Fourier acquisition, which further accelerates the image acquisition. We proposed an ADMM algorithm to solve the constrained image reconstruction problem. Image reconstructions using both simulated and patient data show improved signal-to-noise ratio. As compared to clinically used parallel imaging scheme which achieves a 4-fold acceleration, our method achieves an 8-fold acceleration. Reconstructed images show reduced reconstruction errors as proved on simulated data and similar diffusion parameter mapping results on patient data.PHDElectrical Engineering: SystemsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/143919/1/llliu_1.pd

Magnetic Resonance Imaging of the Brain in Moving Subjects. Application of Fetal, Neonatal and Adult Brain Studies

Imaging in the presence of subject motion has been an ongoing challenge for magnetic resonance imaging (MRI). Motion makes MRI data inconsistent, causing artifacts in conventional anatomical imaging as well as invalidating diffusion tensor imaging (DTI) reconstruction. In this thesis some of the important issues regarding the acquisition and reconstruction of anatomical and DTI imaging of moving subjects are addressed; methods to achieve high resolution and high signalto- noise ratio (SNR) volume data are proposed. An approach has been developed that uses multiple overlapped dynamic single shot slice by slice imaging combined with retrospective alignment and data fusion to produce self consistent 3D volume images under subject motion. We term this method as snapshot MRI with volume reconstruction or SVR. The SVR method has been performed successfully for brain studies on subjects that cannot stay still, and in some cases were moving substantially during scanning. For example, awake neonates, deliberately moved adults and, especially, on fetuses, for which no conventional high resolution 3D method is currently available. Fine structure of the in-utero fetal brain is clearly revealed for the first time with substantially improved SNR. The SVR method has been extended to correct motion artifacts from conventional multi-slice sequences when the subject drifts in position during data acquisition. Besides anatomical imaging, the SVR method has also been further extended to DTI reconstruction when there is subject motion. This has been validated successfully from an adult who was deliberately moving and then applied to inutero fetal brain imaging, which no conventional high resolution 3D method is currently available. Excellent fetal brain 3D apparent diffusion coefficient (ADC) maps in high resolution have been achieved for the first time as well as promising fractional Anisotropy (FA) maps. Pilot clinical studies using SVR reconstructed data to study fetal brain development in-utero have been performed. Growth curves for the normally developing fetal brain have been devised by the quantification of cerebral and cerebellar volumes as well as some one dimensional measurements. A Verhulst model is proposed to describe these growth curves, and this approach has achieved a correlation over 0.99 between the fitted model and actual data

Realistic simulation of artefacts in diffusion MRI for validating post-processing correction techniques

AbstractIn this paper we demonstrate a simulation framework that enables the direct and quantitative comparison of post-processing methods for diffusion weighted magnetic resonance (DW-MR) images. DW-MR datasets are employed in a range of techniques that enable estimates of local microstructure and global connectivity in the brain. These techniques require full alignment of images across the dataset, but this is rarely the case. Artefacts such as eddy-current (EC) distortion and motion lead to misalignment between images, which compromise the quality of the microstructural measures obtained from them. Numerous methods and software packages exist to correct these artefacts, some of which have become de-facto standards, but none have been subject to rigorous validation. In the literature, improved alignment is assessed using either qualitative visual measures or quantitative surrogate metrics. Here we introduce a simulation framework that allows for the direct, quantitative assessment of techniques, enabling objective comparisons of existing and future methods. DW-MR datasets are generated using a process that is based on the physics of MRI acquisition, which allows for the salient features of the images and their artefacts to be reproduced. We apply this framework in three ways. Firstly we assess the most commonly used method for artefact correction, FSL's eddy_correct, and compare it to a recently proposed alternative, eddy. We demonstrate quantitatively that using eddy_correct leads to significant errors in the corrected data, whilst eddy is able to provide much improved correction. Secondly we investigate the datasets required to achieve good correction with eddy, by looking at the minimum number of directions required and comparing the recommended full-sphere acquisitions to equivalent half-sphere protocols. Finally, we investigate the impact of correction quality by examining the fits from microstructure models to real and simulated data