Higher order positive semidefinite diffusion tensor imaging

Due to the well-known limitations of diffusion tensor imaging, high angular resolution diffusion imaging (HARDI) is used to characterize non-Gaussian diffusion processes. One approach to analyzing HARDI data is to model the apparent diffusion coefficient (ADC) with higher order diffusion tensors. The diffusivity function is positive semidefinite. In the literature, some methods have been proposed to preserve positive semidefiniteness of second order and fourth order diffusion tensors. None of them can work for arbitrarily high order diffusion tensors. In this paper, we propose a comprehensive model to approximate the ADC profile by a positive semidefinite diffusion tensor of either second or higher order. We call this the positive semidefinite diffusion tensor (PSDT) model. PSDT is a convex optimization problem with a convex quadratic objective function constrained by the nonnegativity requirement on the smallest Z-eigenvalue of the diffusivity function. The smallest Z-eigenvalue is a computable measure of the extent of positive definiteness of the diffusivity function. We also propose some other invariants for the ADC profile analysis. Experiment results show that higher order tensors could improve the estimation of anisotropic diffusion and that the PSDT model can depict the characterization of diffusion anisotropy which is consistent with known neuroanatomy. © 2010 Society for Industrial and Applied Mathematics.published_or_final_versio

On High Order Tensor-based Diffusivity Profile Estimation

Diffusion weighted magnetic resonance imaging (dMRI) is used to measure, in vivo, the self-diffusion of water molecules in biological tissues. High order tensors (HOTs) are used to model the apparent diffusion coefficient (ADC) profile at each voxel from the dMRI data. In this paper we propose: (i) A new method for estimating HOTs from dMRI data based on weighted least squares (WLS) optimization; and (ii) A new expression for computing the fractional anisotropy from a HOT that does not suffer from singularities and spurious zeros. We also present an empirical evaluation of the proposed method relative to the two existing methods based on both synthetic and real human brain dMRI data. The results show that the proposed method yields more accurate estimation than the competing methods

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

Doctor of Philosophy

dissertationDiffusion magnetic resonance imaging (dMRI) has become a popular technique to detect brain white matter structure. However, imaging noise, imaging artifacts, and modeling techniques, etc., create many uncertainties, which may generate misleading information for further analysis or applications, such as surgical planning. Therefore, how to analyze, effectively visualize, and reduce these uncertainties become very important research questions. In this dissertation, we present both rank-k decomposition and direct decomposition approaches based on spherical deconvolution to decompose the fiber directions more accurately for high angular resolution diffusion imaging (HARDI) data, which will reduce the uncertainties of the fiber directions. By applying volume rendering techniques to an ensemble of 3D orientation distribution function (ODF) glyphs, which we call SIP functions of diffusion shapes, one can elucidate the complex heteroscedastic structural variation in these local diffusion shapes. Furthermore, we quantify the extent of this variation by measuring the fraction of the volume of these shapes, which is consistent across all noise levels, the certain volume ratio. To better understand the uncertainties in white matter fiber tracks, we propose three metrics to quantify the differences between the results of diffusion tensor magnetic resonance imaging (DT-MRI) fiber tracking algorithms: the area between corresponding fibers of each bundle, the Earth Mover's Distance (EMD) between two fiber bundle volumes, and the current distance between two fiber bundle volumes. Based on these metrics, we discuss an interactive fiber track comparison visualization toolkit we have developed to visualize these uncertainties more efficiently. Physical phantoms, with high repeatability and reproducibility, are also designed with the hope of validating the dMRI techniques. In summary, this dissertation provides a better understanding about uncertainties in diffusion magnetic resonance imaging: where and how much are the uncertainties? How do we reduce these uncertainties? How can we possibly validate our algorithms

A simplified algorithm for inverting higher order diffusion tensors

In Riemannian geometry, a distance function is determined by an inner product on the tangent space. In Riemann-Finsler geometry, this distance function can be determined by a norm. This gives more freedom on the form of the so-called indicatrix or the set of unit vectors. This has some interesting applications, e.g., in medical image analysis, especially in diffusion weighted imaging (DWI). An important application of DWI is in the inference of the local architecture of the tissue, typically consisting of thin elongated structures, such as axons or muscle fibers, by measuring the constrained diffusion of water within the tissue. From high angular resolution diffusion imaging (HARDI) data, one can estimate the diffusion orientation distribution function (dODF), which indicates the relative diffusivity in all directions and can be represented by a spherical polynomial. We express this dODF as an equivalent spherical monomial (higher order tensor) to directly generalize the (second order) diffusion tensor approach. To enable efficient computation of Riemann-Finslerian quantities on diffusion weighted (DW)-images, such as the metric/norm tensor, we present a simple and efficient algorithm to invert even order spherical monomials, which extends the familiar inversion of diffusion tensors, i.e., symmetric matrices.</p

Data augmentation in Rician noise model and Bayesian Diffusion Tensor Imaging

Mapping white matter tracts is an essential step towards understanding brain function. Diffusion Magnetic Resonance Imaging (dMRI) is the only noninvasive technique which can detect in vivo anisotropies in the 3-dimensional diffusion of water molecules, which correspond to nervous fibers in the living brain. In this process, spectral data from the displacement distribution of water molecules is collected by a magnetic resonance scanner. From the statistical point of view, inverting the Fourier transform from such sparse and noisy spectral measurements leads to a non-linear regression problem. Diffusion tensor imaging (DTI) is the simplest modeling approach postulating a Gaussian displacement distribution at each volume element (voxel). Typically the inference is based on a linearized log-normal regression model that can fit the spectral data at low frequencies. However such approximation fails to fit the high frequency measurements which contain information about the details of the displacement distribution but have a low signal to noise ratio. In this paper, we directly work with the Rice noise model and cover the full range of $b$-values. Using data augmentation to represent the likelihood, we reduce the non-linear regression problem to the framework of generalized linear models. Then we construct a Bayesian hierarchical model in order to perform simultaneously estimation and regularization of the tensor field. Finally the Bayesian paradigm is implemented by using Markov chain Monte Carlo.Comment: 37 pages, 3 figure

Higher order positive semidefinite diffusion tensor imaging

2010-2011 > Academic research: refereed > Publication in refereed journalVersion of RecordPublishe

Probing white-matter microstructure with higher-order diffusion tensors and susceptibility tensor MRI.

Diffusion MRI has become an invaluable tool for studying white matter microstructure and brain connectivity. The emergence of quantitative susceptibility mapping and susceptibility tensor imaging (STI) has provided another unique tool for assessing the structure of white matter. In the highly ordered white matter structure, diffusion MRI measures hindered water mobility induced by various tissue and cell membranes, while susceptibility sensitizes to the molecular composition and axonal arrangement. Integrating these two methods may produce new insights into the complex physiology of white matter. In this study, we investigated the relationship between diffusion and magnetic susceptibility in the white matter. Experiments were conducted on phantoms and human brains in vivo. Diffusion properties were quantified with the diffusion tensor model and also with the higher order tensor model based on the cumulant expansion. Frequency shift and susceptibility tensor were measured with quantitative susceptibility mapping and susceptibility tensor imaging. These diffusion and susceptibility quantities were compared and correlated in regions of single fiber bundles and regions of multiple fiber orientations. Relationships were established with similarities and differences identified. It is believed that diffusion MRI and susceptibility MRI provide complementary information of the microstructure of white matter. Together, they allow a more complete assessment of healthy and diseased brains