Improving Fiber Alignment in HARDI by Combining Contextual PDE Flow with Constrained Spherical Deconvolution

We propose two strategies to improve the quality of tractography results computed from diffusion weighted magnetic resonance imaging (DW-MRI) data. Both methods are based on the same PDE framework, defined in the coupled space of positions and orientations, associated with a stochastic process describing the enhancement of elongated structures while preserving crossing structures. In the first method we use the enhancement PDE for contextual regularization of a fiber orientation distribution (FOD) that is obtained on individual voxels from high angular resolution diffusion imaging (HARDI) data via constrained spherical deconvolution (CSD). Thereby we improve the FOD as input for subsequent tractography. Secondly, we introduce the fiber to bundle coherence (FBC), a measure for quantification of fiber alignment. The FBC is computed from a tractography result using the same PDE framework and provides a criterion for removing the spurious fibers. We validate the proposed combination of CSD and enhancement on phantom data and on human data, acquired with different scanning protocols. On the phantom data we find that PDE enhancements improve both local metrics and global metrics of tractography results, compared to CSD without enhancements. On the human data we show that the enhancements allow for a better reconstruction of crossing fiber bundles and they reduce the variability of the tractography output with respect to the acquisition parameters. Finally, we show that both the enhancement of the FODs and the use of the FBC measure on the tractography improve the stability with respect to different stochastic realizations of probabilistic tractography. This is shown in a clinical application: the reconstruction of the optic radiation for epilepsy surgery planning

On the Reliability of Diffusion Neuroimaging

Over the last years, diffusion imaging techniques like DTI, DSI or Q-Ball received increasin

Quantitative evaluation of 10 tractography algorithms on a realistic diffusion MR phantom.

International audienceAs it provides the only method for mapping white matter fibers in vivo, diffusion MRI tractography is gaining importance in clinical and neuroscience research. However, despite the increasing availability of different diffusion models and tractography algorithms, it remains unclear how to select the optimal fiber reconstruction method, given certain imaging parameters. Consequently, it is of utmost importance to have a quantitative comparison of these models and algorithms and a deeper understanding of the corresponding strengths and weaknesses. In this work, we use a common dataset with known ground truth and a reproducible methodology to quantitatively evaluate the performance of various diffusion models and tractography algorithms. To examine a wide range of methods, the dataset, but not the ground truth, was released to the public for evaluation in a contest, the "Fiber Cup". 10 fiber reconstruction methods were evaluated. The results provide evidence that: 1. For high SNR datasets, diffusion models such as (fiber) orientation distribution functions correctly model the underlying fiber distribution and can be used in conjunction with streamline tractography, and 2. For medium or low SNR datasets, a prior on the spatial smoothness of either the diffusion model or the fibers is recommended for correct modelling of the fiber distribution and proper tractography results. The phantom dataset, the ground truth fibers, the evaluation methodology and the results obtained so far will remain publicly available on: http://www.lnao.fr/spip.php?rubrique79 to serve as a comparison basis for existing or new tractography methods. New results can be submitted to [email protected] and updates will be published on the webpage

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

Defining Meyer's loop-temporal lobe resections, visual field deficits and diffusion tensor tractography

Anterior temporal lobe resection is often complicated by superior quadrantic visual field deficits (VFDs). In some cases this can be severe enough to prohibit driving, even if a patient is free of seizures. These deficits are caused by damage to Meyer's loop of the optic radiation, which shows considerable heterogeneity in its anterior extent. This structure cannot be distinguished using clinical magnetic resonance imaging sequences. Diffusion tensor tractography is an advanced magnetic resonance imaging technique that enables the parcellation of white matter. Using seed voxels antero-lateral to the lateral geniculate nucleus, we applied this technique to 20 control subjects, and 21 postoperative patients. All patients had visual fields assessed with Goldmann perimetry at least three months after surgery. We measured the distance from the tip of Meyer's loop to the temporal pole and horn in all subjects. In addition, we measured the size of temporal lobe resection using postoperative T1-weighted images, and quantified VFDs. Nine patients suffered VFDs ranging from 22% to 87% of the contralateral superior quadrant. In patients, the range of distance from the tip of Meyer's loop to the temporal pole was 24–43 mm (mean 34 mm), and the range of distance from the tip of Meyer's loop to the temporal horn was –15 to +9 mm (mean 0 mm). In controls the range of distance from the tip of Meyer's loop to the temporal pole was 24–47 mm (mean 35 mm), and the range of distance from the tip of Meyer's loop to the temporal horn was –11 to +9 mm (mean 0 mm). Both quantitative and qualitative results were in accord with recent dissections of cadaveric brains, and analysis of postoperative VFDs and resection volumes. By applying a linear regression analysis we showed that both distance from the tip of Meyer's loop to the temporal pole and the size of resection were significant predictors of the postoperative VFDs. We conclude that there is considerable variation in the anterior extent of Meyer's loop. In view of this, diffusion tensor tractography of the optic radiation is a potentially useful method to assess an individual patient's risk of postoperative VFDs following anterior temporal lobe resection

DeepTract: A Probabilistic Deep Learning Framework for White Matter Fiber Tractography

We present DeepTract, a deep-learning framework for estimating white matter fibers orientation and streamline tractography. We adopt a data-driven approach for fiber reconstruction from diffusion weighted images (DWI), which does not assume a specific diffusion model. We use a recurrent neural network for mapping sequences of DWI values into probabilistic fiber orientation distributions. Based on these estimations, our model facilitates both deterministic and probabilistic streamline tractography. We quantitatively evaluate our method using the Tractometer tool, demonstrating competitive performance with state-of-the art classical and machine learning based tractography algorithms. We further present qualitative results of bundle-specific probabilistic tractography obtained using our method. The code is publicly available at: https://github.com/itaybenou/DeepTract.git