Estimation of white matter fiber parameters from compressed multiresolution diffusion MRI using sparse Bayesian learning.

We present a sparse Bayesian unmixing algorithm BusineX: Bayesian Unmixing for Sparse Inference-based Estimation of Fiber Crossings (X), for estimation of white matter fiber parameters from compressed (under-sampled) diffusion MRI (dMRI) data. BusineX combines compressive sensing with linear unmixing and introduces sparsity to the previously proposed multiresolution data fusion algorithm RubiX, resulting in a method for improved reconstruction, especially from data with lower number of diffusion gradients. We formulate the estimation of fiber parameters as a sparse signal recovery problem and propose a linear unmixing framework with sparse Bayesian learning for the recovery of sparse signals, the fiber orientations and volume fractions. The data is modeled using a parametric spherical deconvolution approach and represented using a dictionary created with the exponential decay components along different possible diffusion directions. Volume fractions of fibers along these directions define the dictionary weights. The proposed sparse inference, which is based on the dictionary representation, considers the sparsity of fiber populations and exploits the spatial redundancy in data representation, thereby facilitating inference from under-sampled q-space. The algorithm improves parameter estimation from dMRI through data-dependent local learning of hyperparameters, at each voxel and for each possible fiber orientation, that moderate the strength of priors governing the parameter variances. Experimental results on synthetic and in-vivo data show improved accuracy with a lower uncertainty in fiber parameter estimates. BusineX resolves a higher number of second and third fiber crossings. For under-sampled data, the algorithm is also shown to produce more reliable estimates.

Sparse Bayesian Inference of White Matter Fiber Orientations from Compressed Multi-resolution Diffusion MRI.

The RubiX [1] algorithm combines high SNR characteristics of low resolution data with high spacial specificity of high resolution data, to extract microstructural tissue parameters from diffusion MRI. In this paper we focus on estimating crossing fiber orientations and introduce sparsity to the RubiX algorithm, making it suitable for reconstruction from compressed (under-sampled) data. We propose a sparse Bayesian algorithm for estimation of fiber orientations and volume fractions from compressed diffusion MRI. The data at high resolution is modeled using a parametric spherical deconvolution approach and represented using a dictionary created with the exponential decay components along different possible directions. Volume fractions of fibers along these orientations define the dictionary weights. The data at low resolution is modeled using a spatial partial volume representation. The proposed dictionary representation and sparsity priors consider the dependence between fiber orientations and the spatial redundancy in data representation. Our method exploits the sparsity of fiber orientations, therefore facilitating inference from under-sampled data. Experimental results show improved accuracy and decreased uncertainty in fiber orientation estimates. For under-sampled data, the proposed method is also shown to produce more robust estimates of fiber orientations.

Estimation of the CSA-ODF using Bayesian compressed sensing of multi-shell HARDI.

PURPOSE: Diffusion MRI provides important information about the brain white matter structures and has opened new avenues for neuroscience and translational research. However, acquisition time needed for advanced applications can still be a challenge in clinical settings. There is consequently a need to accelerate diffusion MRI acquisitions. METHODS: A multi-task Bayesian compressive sensing (MT-BCS) framework is proposed to directly estimate the constant solid angle orientation distribution function (CSA-ODF) from under-sampled (i.e., accelerated image acquisition) multi-shell high angular resolution diffusion imaging (HARDI) datasets, and accurately recover HARDI data at higher resolution in q-space. The proposed MT-BCS approach exploits the spatial redundancy of the data by modeling the statistical relationships within groups (clusters) of diffusion signal. This framework also provides uncertainty estimates of the computed CSA-ODF and diffusion signal, directly computed from the compressive measurements. Experiments validating the proposed framework are performed using realistic multi-shell synthetic images and in vivo multi-shell high angular resolution HARDI datasets. RESULTS: Results indicate a practical reduction in the number of required diffusion volumes (q-space samples) by at least a factor of four to estimate the CSA-ODF from multi-shell data. CONCLUSION: This work presents, for the first time, a multi-task Bayesian compressive sensing approach to simultaneously estimate the full posterior of the CSA-ODF and diffusion-weighted volumes from multi-shell HARDI acquisitions. It demonstrates improvement of the quality of acquired datasets by means of CS de-noising, and accurate estimation of the CSA-ODF, as well as enables a reduction in the acquisition time by a factor of two to four, especially when "staggered" q-space sampling schemes are used. The proposed MT-BCS framework can naturally be combined with parallel MR imaging to further accelerate HARDI acquisitions.

Pushing spatial and temporal resolution for functional and diffusion MRI in the Human Connectome Project.

The Human Connectome Project (HCP) relies primarily on three complementary magnetic resonance (MR) methods. These are: 1) resting state functional MR imaging (rfMRI) which uses correlations in the temporal fluctuations in an fMRI time series to deduce 'functional connectivity'; 2) diffusion imaging (dMRI), which provides the input for tractography algorithms used for the reconstruction of the complex axonal fiber architecture; and 3) task based fMRI (tfMRI), which is employed to identify functional parcellation in the human brain in order to assist analyses of data obtained with the first two methods. We describe technical improvements and optimization of these methods as well as instrumental choices that impact speed of acquisition of fMRI and dMRI images at 3T, leading to whole brain coverage with 2 mm isotropic resolution in 0.7 s for fMRI, and 1.25 mm isotropic resolution dMRI data for tractography analysis with three-fold reduction in total dMRI data acquisition time. Ongoing technical developments and optimization for acquisition of similar data at 7 T magnetic field are also presented, targeting higher spatial resolution, enhanced specificity of functional imaging signals, mitigation of the inhomogeneous radio frequency (RF) fields, and reduced power deposition. Results demonstrate that overall, these approaches represent a significant advance in MR imaging of the human brain to investigate brain function and structure.

A Framework for Linear and Non-Linear Registration of Diffusion-Weighted MRIs Using Angular Interpolation.

Registration of diffusion-weighted magnetic resonance images (DW-MRIs) is a key step for population studies, or construction of brain atlases, among other important tasks. Given the high dimensionality of the data, registration is usually performed by relying on scalar representative images, such as the fractional anisotropy (FA) and non-diffusion-weighted (b0) images, thereby ignoring much of the directional information conveyed by DW-MR datasets itself. Alternatively, model-based registration algorithms have been proposed to exploit information on the preferred fiber orientation(s) at each voxel. Models such as the diffusion tensor or orientation distribution function (ODF) have been used for this purpose. Tensor-based registration methods rely on a model that does not completely capture the information contained in DW-MRIs, and largely depends on the accurate estimation of tensors. ODF-based approaches are more recent and computationally challenging, but also better describe complex fiber configurations thereby potentially improving the accuracy of DW-MRI registration. A new algorithm based on angular interpolation of the diffusion-weighted volumes was proposed for affine registration, and does not rely on any specific local diffusion model. In this work, we first extensively compare the performance of registration algorithms based on (i) angular interpolation, (ii) non-diffusion-weighted scalar volume (b0), and (iii) diffusion tensor image (DTI). Moreover, we generalize the concept of angular interpolation (AI) to non-linear image registration, and implement it in the FMRIB Software Library (FSL). We demonstrate that AI registration of DW-MRIs is a powerful alternative to volume and tensor-based approaches. In particular, we show that AI improves the registration accuracy in many cases over existing state-of-the-art algorithms, while providing registered raw DW-MRI data, which can be used for any subsequent analysis.

Hierarchical topological network analysis of anatomical human brain connectivity and differences related to sex and kinship.

Modern non-invasive brain imaging technologies, such as diffusion weighted magnetic resonance imaging (DWI), enable the mapping of neural fiber tracts in the white matter, providing a basis to reconstruct a detailed map of brain structural connectivity networks. Brain connectivity networks differ from random networks in their topology, which can be measured using small worldness, modularity, and high-degree nodes (hubs). Still, little is known about how individual differences in structural brain network properties relate to age, sex, or genetic differences. Recently, some groups have reported brain network biomarkers that enable differentiation among individuals, pairs of individuals, and groups of individuals. In addition to studying new topological features, here we provide a unifying general method to investigate topological brain networks and connectivity differences between individuals, pairs of individuals, and groups of individuals at several levels of the data hierarchy, while appropriately controlling false discovery rate (FDR) errors. We apply our new method to a large dataset of high quality brain connectivity networks obtained from High Angular Resolution Diffusion Imaging (HARDI) tractography in 303 young adult twins, siblings, and unrelated people. Our proposed approach can accurately classify brain connectivity networks based on sex (93% accuracy) and kinship (88.5% accuracy). We find statistically significant differences associated with sex and kinship both in the brain connectivity networks and in derived topological metrics, such as the clustering coefficient and the communicability matrix.

Comparative study of semi-implicit schemes for nonlinear diffusion in hyperspectral imagery.

Nonlinear diffusion has been successfully employed over the past two decades to enhance images by reducing undesirable intensity variability within the objects in the image, while enhancing the contrast of the boundaries (edges) in scalar and, more recently, in vector-valued images, such as color, multispectral, and hyperspectral imagery. In this paper, we show that nonlinear diffusion can improve the classification accuracy of hyperspectral imagery by reducing the spatial and spectral variability of the image, while preserving the boundaries of the objects. We also show that semi-implicit schemes can speedup significantly the evolution of the nonlinear diffusion equation with respect to traditional explicit schemes.