Diffeomorphic Metric Mapping of High Angular Resolution Diffusion Imaging based on Riemannian Structure of Orientation Distribution Functions

In this paper, we propose a novel large deformation diffeomorphic registration algorithm to align high angular resolution diffusion images (HARDI) characterized by orientation distribution functions (ODFs). Our proposed algorithm seeks an optimal diffeomorphism of large deformation between two ODF fields in a spatial volume domain and at the same time, locally reorients an ODF in a manner such that it remains consistent with the surrounding anatomical structure. To this end, we first review the Riemannian manifold of ODFs. We then define the reorientation of an ODF when an affine transformation is applied and subsequently, define the diffeomorphic group action to be applied on the ODF based on this reorientation. We incorporate the Riemannian metric of ODFs for quantifying the similarity of two HARDI images into a variational problem defined under the large deformation diffeomorphic metric mapping (LDDMM) framework. We finally derive the gradient of the cost function in both Riemannian spaces of diffeomorphisms and the ODFs, and present its numerical implementation. Both synthetic and real brain HARDI data are used to illustrate the performance of our registration algorithm

A linear and regularized ODF estimation algorithm to recover multiple fibers in Q-Ball imaging

Due the well-known limitations of diffusion tensor imaging (DTI), high angular resolution diffusion imaging is currently of great interest to characterize voxels containing multiple fiber crossings. In particular, Q-ball imaging (QBI) is now a popular reconstruction method to obtain the orientation distribution function (ODF) of these multiple fiber distributions. The latter captures all important angular contrast by expressing the probability that a water molecule will diffuse into any given solid angle. However, QBI and other high order spin displacement estimation methods involve non-trivial numerical computations and lack a straightforward regularization process. In this paper, we propose a simple linear and regularized analytic solution for the Q-ball reconstruction of the ODF. First, the signal is modeled with a physically meaningful high order spherical harmonic series by incorporating the Laplace-Beltrami operator in the solution. This leads to an elegant mathematical simplification of the Funk-Radon transform using the Funk-Hecke formula. In doing so, we obtain a fast and robust model-free ODF approximation. We validate the accuracy of the ODF estimation quantitatively using the multi-tensor synthetic model where the exact ODF can be computed. We also demonstrate that the estimated ODF can recover known multiple fiber regions in a biological phantom and in the human brain. Another important contribution of the paper is the development of ODF sharpening methods. We show that sharpening the measured ODF enhances each underlying fiber compartment and considerably improves the extraction of fibers. The proposed techniques are simple linear transformations of the ODF and can easily be computed using our spherical harmonics machinery

Ball and rackets: inferring fiber fanning from diffusion-weighted MRI

A number of methods have been proposed for resolving crossing fibers from diffusion-weighted (DW) MRI. However, other complex fiber geometries have drawn minimal attention. In this study, we focus on fiber orientation dispersion induced by within-voxel fanning. We use a multi-compartment, model-based approach to estimate fiber dispersion. Bingham distributions are employed to represent continuous distributions of fiber orientations, centered around a main orientation, and capturing anisotropic dispersion. We evaluate the accuracy of the model for different simulated fanning geometries, under different acquisition protocols and we illustrate the high SNR and angular resolution needs. We also perform a qualitative comparison between our parametric approach and five popular non-parametric techniques that are based on orientation distribution functions (ODFs). This comparison illustrates how the same underlying geometry can be depicted by different methods. We apply the proposed model on high-quality, post-mortem macaque data and present whole-brain maps of fiber dispersion, as well as exquisite details on the local anatomy of fiber distributions in various white matter regions

Improving Estimation of Fiber Orientations in Diffusion MRI Using Inter-Subject Information Sharing

Diffusion magnetic resonance imaging is widely used to investigate diffusion patterns of water molecules in the human brain. It provides information that is useful for tracing axonal bundles and inferring brain connectivity. Diffusion axonal tracing, namely tractography, relies on local directional information provided by the orientation distribution functions (ODFs) estimated at each voxel. To accurately estimate ODFs, data of good signal-to-noise ratio and sufficient angular samples are desired. This is however not always available in practice. In this paper, we propose to improve ODF estimation by using inter-subject image correlation. Specifically, we demonstrate that diffusion-weighted images acquired from different subjects can be transformed to the space of a target subject to drastically increase the number of angular samples to improve ODF estimation. This is largely due to the incoherence of the angular samples generated when the diffusion signals are reoriented and warped to the target space. To reorient the diffusion signals, we propose a new spatial normalization method that directly acts on diffusion signals using local affine transforms. Experiments on both synthetic data and real data show that our method can reduce noise-induced artifacts, such as spurious ODF peaks, and yield more coherent orientations

Efficient Sampling for Accelerated Diffusion Magnetic Resonance Imaging

Diffusion magnetic resonance imaging (dMRI) is a non-invasive method that allows connectivity mapping of the brain. However, despite major advances in this field, accurate inference of these patterns and its applicability within a clinical context is still in its early stages. This thesis describes a conceptually novel method for reconstructing neuronal pathways inside the brain from diffusion-weighted imaging (DWI) measurements with high angular resolution and short data acquisition time. The proposed method combines recent theoretical advances on spherical sampling and noise reduction techniques from the field of compressed sensing. Numerical simulations were performed to study the best sampling strategy under a novel sampling theorem on the sphere in order to reduce the acquisition time during dMRI scans. Furthermore, these results were combined with the recently proposed spherical deconvolution technique to reconstruct the distribution of neuronal tracts (or fibers) within one voxel with high angular resolution between multiple crossing fibers. The spherical deconvolution problem was hereby formulated as an inverse problem and solved using techniques adopted from the field of compressed sensing. Since the result of the spherical deconvolution step is sparse in nature, the basis pursuit denoising formulation of the inverse problem is optimal within this context. Finally, the resulting fiber orientation reconstruction was compared with diffusion spectrum imaging (DSI) – a classic model-free acquisition method. Simulations revealed that the proposed approach is superior to DSI in terms of both, acquisition time and angular resolution of crossing fibers (>=40° with at least 90% sensitivity). Our investigations showed that the application of spherical deconvolution stated as a basis pursuit denoising problem holds great promise for high angular resolution dMRI