On diffusion tensor estimation

Visualization of human optic nerve by diffusion tensor mapping and degree of neuropathy

Diffusion-weighted magnetic resonance imaging of the human optic nerve and tract is technically difficult because of its small size, the inherent strong signal generated by the surrounding fat and the cerebrospinal fluid, and due to eddy current-induced distortions and subject movement artifacts. The effects of the bone canal through which the optic nerve passes, and the proximity of blood vessels, muscles and tendons are generally unknown. Also, the limited technical capabilities of the scanners and the minimization of acquisition times result in poor quality diffusion-weighted images. It is challenging for current tractography methods to accurately track optic pathway fibers that correspond to known anatomy. Despite these technical limitations and low image resolution, here we show how to visualize the optic nerve and tract and quantify nerve atrophy. Our visualization method based on the analysis of the diffusion tensor shows marked differences between a healthy male subject and a male subject with progressive optic nerve neuropathy. These differences coincide with diffusion scalar metrics and are not visible on standard morphological images. A quantification of the degree of optic nerve atrophy in a systematic way is provided and it is tested on 9 subjects from the Human Connectome Project.

Brain Waves: Emergence of Localized, Persistent, Weakly Evanescent Cortical Loops

An inhomogeneous anisotropic physical model of the brain cortex is presented that predicts the emergence of nonevanescent (weakly damped) wave-like modes propagating in the thin cortex layers transverse to both the mean neural fiber direction and the cortex spatial gradient. Although the amplitude of these modes stays below the typically observed axon spiking potential, the lifetime of these modes may significantly exceed the spiking potential inverse decay constant. Full-brain numerical simulations based on parameters extracted from diffusion and structural MRI confirm the existence and extended duration of these wave modes. Contrary to the commonly agreed paradigm that the neural fibers determine the pathways for signal propagation in the brain, the signal propagation because of the cortex wave modes in the highly folded areas will exhibit no apparent correlation with the fiber directions. Nonlinear coupling of those linear weakly evanescent wave modes then provides a universal mechanism for the emergence of synchronized brain wave field activity. The resonant and nonresonant terms of nonlinear coupling between multiple modes produce both synchronous spiking-like high-frequency wave activity as well as low-frequency wave rhythms. Numerical simulation of forced multiple-mode dynamics shows that, as forcing increases, there is a transition from damped to oscillatory regime that can then transition quickly to a nonoscillatory state when a critical excitation threshold is reached. The resonant nonlinear coupling results in the emergence of low-frequency rhythms with frequencies that are several orders of magnitude below the linear frequencies of modes taking part in the coupling. The localization and persistence of these weakly evanescent cortical wave modes have significant implications in particular for neuroimaging methods that detect electromagnetic physiological activity, such as EEG and magnetoencephalography, and for the understanding of brain activity in general, including mechanisms of memory.

Universal theory of brain waves: from linear loops to nonlinear synchronized spiking and collective brain rhythms

An inhomogeneous anisotropic physical model of the brain cortex is presented that predicts the emergence of non-evanescent (weakly damped) wave-like modes propagating in the thin cortex layers transverse to both the mean neural fiber direction and to the cortex spatial gradient. Although the amplitude of these modes stays below the typically observed axon spiking potential, the lifetime of these modes may significantly exceed the spiking potential inverse decay constant. Full brain numerical simulations based on parameters extracted from diffusion and structural MRI confirm the existence and extended duration of these wave modes. Contrary to the standard paradigm that the neural fibers determine the pathways for signal propagation in the brain, the signal propagation due to the cortex wave modes in highly folded areas will exhibit no apparent correlation with the fiber directions. The results are consistent with numerous recent experimental animal and human brain studies demonstrating the existence of electrostatic field activity in the form of traveling waves (including studies where neuronal connections were severed) and with wave loop induced peaks observed in EEG spectra. In addition, we demonstrate that the resonant and non-resonant terms of the nonlinear coupling between multiple modes produce both synchronous spiking-like high frequency wave activity as well as low frequency wave rhythms as a result of their unique dispersion properties. Numerical simulation of forced multiple mode dynamics shows that as forcing increases there is a transition from damped to oscillatory regime that subsequently decays away as over-excitation is reached. The resonant nonlinear coupling results in the emergence of low frequency rhythms with frequencies that are several orders of magnitude below the linear frequencies of modes taking part in the coupling. The localization and persistence of these cortical wave modes, and this new mechanism for understanding the nature of spiking behavior, have significant implications in particular for neuroimaging methods that detect electromagnetic physiological activity, such as EEG and MEG, and in general for the understanding of brain activity, including mechanisms of memory.

A robust variational approach for simultaneous smoothing and estimation of DTI

Estimating diffusion tensors is an essential step in many applications - such as diffusion tensor image (DTI) registration, segmentation and fiber tractography. Most of the methods proposed in the literature for this task are not simultaneously statistically robust and feature preserving techniques. In this paper, we propose a novel and robust variational framework for simultaneous smoothing and estimation of diffusion tensors from diffusion MRI. Our variational principle makes use of a recently introduced total Kullback-Leibler (tKL) divergence for DTI regularization. tKL is a statistically robust dissimilarity measure for diffusion tensors, and regularization by using tKL ensures the symmetric positive definiteness of tensors automatically. Further, the regularization is weighted by a non-local factor adapted from the conventional non-local means filters. Finally, for the data fidelity, we use the nonlinear least-squares term derived from the Stejskal-Tanner model. We present experimental results depicting the positive performance of our method in comparison to competing methods on synthetic and real data examples.

A nonparametric Riemannian framework for processing high angular resolution diffusion images and its applications to ODF-based morphometry

High angular resolution diffusion imaging (HARDI) has become an important technique for imaging complex oriented structures in the brain and other anatomical tissues. This has motivated the recent development of several methods for computing the orientation probability density function (PDF) at each voxel. However, much less work has been done on developing techniques for filtering, interpolation, averaging and principal geodesic analysis of orientation PDF fields. In this paper, we present a Riemannian framework for performing such operations. The proposed framework does not require that the orientation PDFs be represented by any fixed parameterization, such as a mixture of von Mises-Fisher distributions or a spherical harmonic expansion. Instead, we use a nonparametric representation of the orientation PDF. We exploit the fact that under the square-root re-parameterization, the space of orientation PDFs forms a Riemannian manifold: the positive orthant of the unit Hilbert sphere. We show that various orientation PDF processing operations, such as filtering, interpolation, averaging and principal geodesic analysis, may be posed as optimization problems on the Hilbert sphere, and can be solved using Riemannian gradient descent. We illustrate these concepts with numerous experiments on synthetic, phantom and real datasets. We show their application to studying left/right brain asymmetries.

Diffusion tensor imaging at low SNR: nonmonotonic behaviors of tensor contrasts

Diffusion tensor imaging (DTI) provides measurements of directional diffusivities and has been widely used to characterize changes in the tissue microarchitecture of the brain. DTI is gaining prominence in applications outside of the brain, where resolution, motion and short T2 values often limit the achievable signal-to-noise ratio (SNR). Consequently, it is important to revisit the topic of tensor estimation in low-SNR regimes. A theoretical framework is developed to model noise in DTI, and by using simulations based on this theory, the degree to which the noise, tensor estimation method and acquisition protocol affect tensor-derived quantities, such as fractional anisotropy and apparent diffusion coefficient, is clarified. These results are then validated against clinical data. It is shown that reliability of tensor contrasts depends on the noise level, estimation method, diffusion-weighting scheme and underlying anatomy. The propensity for bias and errors does not monotonically increase with noise. Comparative results are shown in both graphical and tabular forms, so that decisions about suitable acquisition protocols and processing methods can be made on a case-by-case basis without exhaustive experimentation.

Shape-adaptive DCT for denoising of 3D scalar and tensor valued images

During the last ten years or so, diffusion tensor imaging has been used in both research and clinical medical applications. To construct the diffusion tensor images, a large set of direction sensitive magnetic resonance image (MRI) acquisitions are required. These acquisitions in general have a lower signal-to-noise ratio than conventional MRI acquisitions. In this paper, we discuss computationally effective algorithms for noise removal for diffusion tensor magnetic resonance imaging (DTI) using the framework of 3-dimensional shape-adaptive discrete cosine transform. We use local polynomial approximations for the selection of homogeneous regions in the DTI data. These regions are transformed to the frequency domain by a modified discrete cosine transform. In the frequency domain, the noise is removed by thresholding. We perform numerical experiments on 3D synthetical MRI and DTI data and real 3D DTI brain data from a healthy volunteer. The experiments indicate good performance compared to current state-of-the-art methods. The proposed method is well suited for parallelization and could thus dramatically improve the computation speed of denoising schemes for large scale 3D MRI and DTI.

Diffusion tensor analysis with invariant gradients and rotation tangents

Guided by empirically established connections between clinically important tissue properties and diffusion tensor parameters, we introduce a framework for decomposing variations in diffusion tensors into changes in shape and orientation. Tensor shape and orientation both have three degrees-of-freedom, spanned by invariant gradients and rotation tangents, respectively. As an initial demonstration of the framework, we create a tunable measure of tensor difference that can selectively respond to shape and orientation. Second, to analyze the spatial gradient in a tensor volume (a third-order tensor), our framework generates edge strength measures that can discriminate between different neuroanatomical boundaries, as well as creating a novel detector of white matter tracts that are adjacent yet distinctly oriented. Finally, we apply the framework to decompose the fourth-order diffusion covariance tensor into individual and aggregate measures of shape and orientation covariance, including a direct approximation for the variance of tensor invariants such as fractional anisotropy.

Diffusion Tensor Estimation by Maximizing Rician Likelihood

Diffusion tensor imaging (DTI) is widely used to characterize white matter in health and disease. Previous approaches to the estimation of diffusion tensors have either been statistically suboptimal or have used Gaussian approximations of the underlying noise structure, which is Rician in reality. This can cause quantities derived from these tensors - e.g., fractional anisotropy and apparent diffusion coefficient - to diverge from their true values, potentially leading to artifactual changes that confound clinically significant ones. This paper presents a novel maximum likelihood approach to tensor estimation, denoted Diffusion Tensor Estimation by Maximizing Rician Likelihood (DTEMRL). In contrast to previous approaches, DTEMRL considers the joint distribution of all observed data in the context of an augmented tensor model to account for variable levels of Rician noise. To improve numeric stability and prevent non-physical solutions, DTEMRL incorporates a robust characterization of positive definite tensors and a new estimator of underlying noise variance. In simulated and clinical data, mean squared error metrics show consistent and significant improvements from low clinical SNR to high SNR. DTEMRL may be readily supplemented with spatial regularization or a priori tensor distributions for Bayesian tensor estimation.