Non-Negative Spherical Deconvolution (NNSD) for Fiber Orientation Distribution Function Estimation

International audienceIn diffusion Magnetic Resonance Imaging (dMRI), Spherical Deconvolution (SD) is a commonly used approach for estimating the fiber Orientation Distribution Function (fODF). As a Probability Density Function (PDF) that characterizes the distribution of fiber orientations, the fODF is expected to be non-negative and to integrate to unity on the continuous unit sphere S2 . However, many existing approaches, despite using continuous representation such as Spherical Harmonics (SH), impose non-negativity only on discretized points of S2. Therefore, non-negativity is not guaranteed on the whole S2 Existing approaches are also known to .exhibit false positive fODF peaks, especially in regions with low anisotropy, causing an over-estimation of the number of fascicles that traverse each voxel. This paper proposes a novel approach, called Non-Negative SD (NNSD), to overcome the above limitations. NNSD offers the following advantages. First, NNSD is the first SH based method that guarantees non-negativity of the fODF throughout the unit sphere. Second, unlike approaches such as Maximum Entropy SD (MESD), Cartesian Tensor Fiber Orientation Distribution (CT-FOD), and discrete representation based SD (DR-SD) techniques, the SH representation allows closed form of spherical integration, efficient computation in a low dimensional space resided by the SH coefficients, and accurate peak detection on the continuous domain defined by the unit sphere. Third, NNSD is significantly less susceptible to producing false positive peaks in regions with low anisotropy. Evaluations of NNSD in comparison with Constrained SD (CSD), MESD, and DR-SD (implemented using L1-regularized least-squares with non-negative constraint), indicate that NNSD yields improved performance for both synthetic and real data. The performance gain is especially prominent for high resolution (1.25 mm)^3 data

Single- and Multiple-Shell Uniform Sampling Schemes for Diffusion MRI Using Spherical Codes

In diffusion MRI (dMRI), a good sampling scheme is important for efficient acquisition and robust reconstruction. Diffusion weighted signal is normally acquired on single or multiple shells in q-space. Signal samples are typically distributed uniformly on different shells to make them invariant to the orientation of structures within tissue, or the laboratory coordinate frame. The Electrostatic Energy Minimization (EEM) method, originally proposed for single shell sampling scheme in dMRI, was recently generalized to multi-shell schemes, called Generalized EEM (GEEM). GEEM has been successfully used in the Human Connectome Project (HCP). However, EEM does not directly address the goal of optimal sampling, i.e., achieving large angular separation between sampling points. In this paper, we propose a more natural formulation, called Spherical Code (SC), to directly maximize the minimal angle between different samples in single or multiple shells. We consider not only continuous problems to design single or multiple shell sampling schemes, but also discrete problems to uniformly extract sub-sampled schemes from an existing single or multiple shell scheme, and to order samples in an existing scheme. We propose five algorithms to solve the above problems, including an incremental SC (ISC), a sophisticated greedy algorithm called Iterative Maximum Overlap Construction (IMOC), an 1-Opt greedy method, a Mixed Integer Linear Programming (MILP) method, and a Constrained Non-Linear Optimization (CNLO) method. To our knowledge, this is the first work to use the SC formulation for single or multiple shell sampling schemes in dMRI. Experimental results indicate that SC methods obtain larger angular separation and better rotational invariance than the state-of-the-art EEM and GEEM. The related codes and a tutorial have been released in DMRITool.Comment: Accepted by IEEE transactions on Medical Imaging. Codes have been released in dmritool https://diffusionmritool.github.io/tutorial_qspacesampling.htm

Increasing the Analytical Accessibility of Multishell and Diffusion Spectrum Imaging Data Using Generalized Q-Sampling Conversion

Many diffusion MRI researchers, including the Human Connectome Project (HCP), acquire data using multishell (e.g., WU-Minn consortium) and diffusion spectrum imaging (DSI) schemes (e.g., USC-Harvard consortium). However, these data sets are not readily accessible to high angular resolution diffusion imaging (HARDI) analysis methods that are popular in connectomics analysis. Here we introduce a scheme conversion approach that transforms multishell and DSI data into their corresponding HARDI representations, thereby empowering HARDI-based analytical methods to make use of data acquired using non-HARDI approaches. This method was evaluated on both phantom and in-vivo human data sets by acquiring multishell, DSI, and HARDI data simultaneously, and comparing the converted HARDI, from non-HARDI methods, with the original HARDI data. Analysis on the phantom shows that the converted HARDI from DSI and multishell data strongly predicts the original HARDI (correlation coefficient > 0.9). Our in-vivo study shows that the converted HARDI can be reconstructed by constrained spherical deconvolution, and the fiber orientation distributions are consistent with those from the original HARDI. We further illustrate that our scheme conversion method can be applied to HCP data, and the converted HARDI do not appear to sacrifice angular resolution. Thus this novel approach can benefit all HARDI-based analysis approaches, allowing greater analytical accessibility to non-HARDI data, including data from the HCP

Non-Local Non-Negative Spherical Deconvolution for Single and Multiple Shell Diffusion MRI

International audienceIn diffusion MRI (dMRI), Spherical Deconvolution (SD) is a category of methods which estimate the fiber Orientation Distribution Function (fODF). Existing SD methods, including the widely used Constrained SD [1], normally have two common limitations: 1) the non-negativity constraint of the fODFs is not satisfied in the continuous sphere; 2) many spurious peaks are detected, especially in the regions with low anisotropy; In [2], we proposed a novel SD method, called Non-Negative SD (NNSD), to avoid these two limitations. NNSD guarantees the non-negativity constraint of fODFs in the continuous sphere S2 , and it is robust to the false positive peaks. In this abstract, we propose Non-Local NNSD (NLNNSD) which considers non-local spatial information and Rician noise in NNSD, and apply it to the testing data in ISBI contest

Exploring Local White Matter Geometric Structure in diffusion MRI Using Director Field Analysis

International audienceIn this abstract, inspired by microscopic theoretical treatment of phases in liquid crystals, we introduce a novel mathematical framework, called Director Field Analysis (DFA), to study local geometric structural information of white matter. DFA extracts some meaningful scalar indices related with the degree of orientational alignment, dispersion, and orientational distortion, from the Orientation Distribution Function (ODF) field reconstructed by Diffusion Tensor Imaging (DTI) or High Angular Resolution Diffusion Imaging (HARDI)

Director Field Analysis to Explore Local White Matter Geometric Structure in diffusion MRI

International audienceIn diffusion MRI, a tensor field or a spherical function field, e.g., an Orientation Distribution Function (ODF) field, are estimated from measured diffusion weighted images. In this paper, inspired by microscopic theoretical treatment of phases in liquid crystals, we introduce a novel mathematical framework, called Director Field Analysis (DFA), to study local geometric structural information of white matter from the estimated tensor field or spherical function field. 1) We propose Orientational Order (OO) and Orientational Dispersion (OD) indices to describe the degree of alignment and dispersion of a spherical function in each voxel; 2) We estimate a local orthogonal coordinate frame in each voxel with anisotropic diffusion; 3) Finally, we define three indices to describe three types of orientational distortion (splay, bend, and twist) in a local spatial neighborhood, and a total distortion index to describe distortions of all three types. To our knowledge, this is the first work to quantitatively describe orientational distortion (splay, bend, and twist) in diffusion MRI. The proposed scalar indices are useful to detect local geometric changes of white matter for voxel-based or tract-based analysis in both DTI and HARDI acquisitions

Estimation of Fiber Orientations Using Neighborhood Information

Data from diffusion magnetic resonance imaging (dMRI) can be used to reconstruct fiber tracts, for example, in muscle and white matter. Estimation of fiber orientations (FOs) is a crucial step in the reconstruction process and these estimates can be corrupted by noise. In this paper, a new method called Fiber Orientation Reconstruction using Neighborhood Information (FORNI) is described and shown to reduce the effects of noise and improve FO estimation performance by incorporating spatial consistency. FORNI uses a fixed tensor basis to model the diffusion weighted signals, which has the advantage of providing an explicit relationship between the basis vectors and the FOs. FO spatial coherence is encouraged using weighted l1-norm regularization terms, which contain the interaction of directional information between neighbor voxels. Data fidelity is encouraged using a squared error between the observed and reconstructed diffusion weighted signals. After appropriate weighting of these competing objectives, the resulting objective function is minimized using a block coordinate descent algorithm, and a straightforward parallelization strategy is used to speed up processing. Experiments were performed on a digital crossing phantom, ex vivo tongue dMRI data, and in vivo brain dMRI data for both qualitative and quantitative evaluation. The results demonstrate that FORNI improves the quality of FO estimation over other state of the art algorithms.Comment: Journal paper accepted in Medical Image Analysis. 35 pages and 16 figure

Higher-Order Tensors and Differential Topology in Diffusion MRI Modeling and Visualization

Diffusion Weighted Magnetic Resonance Imaging (DW-MRI) is a noninvasive method for creating three-dimensional scans of the human brain. It originated mostly in the 1970s and started its use in clinical applications in the 1980s. Due to its low risk and relatively high image quality it proved to be an indispensable tool for studying medical conditions as well as for general scientific research. For example, it allows to map fiber bundles, the major neuronal pathways through the brain. But all evaluation of scanned data depends on mathematical signal models that describe the raw signal output and map it to biologically more meaningful values. And here we find the most potential for improvement. In this thesis we first present a new multi-tensor kurtosis signal model for DW-MRI. That means it can detect multiple overlapping fiber bundles and map them to a set of tensors. Compared to other already widely used multi-tensor models, we also add higher order kurtosis terms to each fiber. This gives a more detailed quantification of fibers. These additional values can also be estimated by the Diffusion Kurtosis Imaging (DKI) method, but we show that these values are drastically affected by fiber crossings in DKI, whereas our model handles them as intrinsic properties of fiber bundles. This reduces the effects of fiber crossings and allows a more direct examination of fibers. Next, we take a closer look at spherical deconvolution. It can be seen as a generalization of multi-fiber signal models to a continuous distribution of fiber directions. To this approach we introduce a novel mathematical constraint. We show, that state-of-the-art methods for estimating the fiber distribution become more robust and gain accuracy when enforcing our constraint. Additionally, in the context of our own deconvolution scheme, it is algebraically equivalent to enforcing that the signal can be decomposed into fibers. This means, tractography and other methods that depend on identifying a discrete set of fiber directions greatly benefit from our constraint. Our third major contribution to DW-MRI deals with macroscopic structures of fiber bundle geometry. In recent years the question emerged, whether or not, crossing bundles form two-dimensional surfaces inside the brain. Although not completely obvious, there is a mathematical obstacle coming from differential topology, that prevents general tangential planes spanned by fiber directions at each point to be connected into consistent surfaces. Research into how well this constraint is fulfilled in our brain is hindered by the high precision and complexity needed by previous evaluation methods. This is why we present a drastically simpler method that negates the need for precisely finding fiber directions and instead only depends on the simple diffusion tensor method (DTI). We then use our new method to explore and improve streamsurface visualization.<br /