Dataset Information

Spherical Deconvolution of Multichannel Diffusion MRI Data with Non-Gaussian Noise Models and Spatial Regularization.

ABSTRACT: Spherical deconvolution (SD) methods are widely used to estimate the intra-voxel white-matter fiber orientations from diffusion MRI data. However, while some of these methods assume a zero-mean Gaussian distribution for the underlying noise, its real distribution is known to be non-Gaussian and to depend on many factors such as the number of coils and the methodology used to combine multichannel MRI signals. Indeed, the two prevailing methods for multichannel signal combination lead to noise patterns better described by Rician and noncentral Chi distributions. Here we develop a Robust and Unbiased Model-BAsed Spherical Deconvolution (RUMBA-SD) technique, intended to deal with realistic MRI noise, based on a Richardson-Lucy (RL) algorithm adapted to Rician and noncentral Chi likelihood models. To quantify the benefits of using proper noise models, RUMBA-SD was compared with dRL-SD, a well-established method based on the RL algorithm for Gaussian noise. Another aim of the study was to quantify the impact of including a total variation (TV) spatial regularization term in the estimation framework. To do this, we developed TV spatially-regularized versions of both RUMBA-SD and dRL-SD algorithms. The evaluation was performed by comparing various quality metrics on 132 three-dimensional synthetic phantoms involving different inter-fiber angles and volume fractions, which were contaminated with noise mimicking patterns generated by data processing in multichannel scanners. The results demonstrate that the inclusion of proper likelihood models leads to an increased ability to resolve fiber crossings with smaller inter-fiber angles and to better detect non-dominant fibers. The inclusion of TV regularization dramatically improved the resolution power of both techniques. The above findings were also verified in human brain data.

SUBMITTER: Canales-Rodriguez EJ

PROVIDER: S-EPMC4607500 | biostudies-literature | 2015

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Spherical Deconvolution of Multichannel Diffusion MRI Data with Non-Gaussian Noise Models and Spatial Regularization.

Canales-Rodríguez Erick J EJ Daducci Alessandro A Sotiropoulos Stamatios N SN Caruyer Emmanuel E Aja-Fernández Santiago S Radua Joaquim J Yurramendi Mendizabal Jesús M JM Iturria-Medina Yasser Y Melie-García Lester L Alemán-Gómez Yasser Y Thiran Jean-Philippe JP Sarró Salvador S Pomarol-Clotet Edith E Salvador Raymond R

PloS one 20151015 10

Spherical deconvolution (SD) methods are widely used to estimate the intra-voxel white-matter fiber orientations from diffusion MRI data. However, while some of these methods assume a zero-mean Gaussian distribution for the underlying noise, its real distribution is known to be non-Gaussian and to depend on many factors such as the number of coils and the methodology used to combine multichannel MRI signals. Indeed, the two prevailing methods for multichannel signal combination lead to noise pat ...[more]

PMID: 26470024

Similar Datasets

Project description:Virtual dissection of white matter (WM) using diffusion MRI tractography is confounded by its poor reproducibility. Despite the increased adoption of advanced reconstruction models, early region-of-interest driven protocols based on diffusion tensor imaging (DTI) remain the dominant reference for virtual dissection protocols. Here we bridge this gap by providing a comprehensive description of typical WM anatomy reconstructed using a reproducible automated subject-specific parcellation-based approach based on probabilistic constrained-spherical deconvolution (CSD) tractography. We complement this with a WM template in MNI space comprising 68 bundles, including all associated anatomical tract selection labels and associated automated workflows. Additionally, we demonstrate bundle inter- and intra-subject variability using 40 (20 test-retest) datasets from the human connectome project (HCP) and 5 sessions with varying b-values and number of b-shells from the single-subject Multiple Acquisitions for Standardization of Structural Imaging Validation and Evaluation (MASSIVE) dataset. The most reliably reconstructed bundles were the whole pyramidal tracts, primary corticospinal tracts, whole superior longitudinal fasciculi, frontal, parietal and occipital segments of the corpus callosum and middle cerebellar peduncles. More variability was found in less dense bundles, e.g., the fornix, dentato-rubro-thalamic tract (DRTT), and premotor pyramidal tract. Using the DRTT as an example, we show that this variability can be reduced by using a higher number of seeding attempts. Overall inter-session similarity was high for HCP test-retest data (median weighted-dice = 0.963, stdev = 0.201 and IQR = 0.099). Compared to the HCP-template bundles there was a high level of agreement for the HCP test-retest data (median weighted-dice = 0.747, stdev = 0.220 and IQR = 0.277) and for the MASSIVE data (median weighted-dice = 0.767, stdev = 0.255 and IQR = 0.338). In summary, this WM atlas provides an overview of the capabilities and limitations of automated subject-specific probabilistic CSD tractography for mapping white matter fasciculi in healthy adults. It will be most useful in applications requiring a reproducible parcellation-based dissection protocol, and as an educational resource for applied neuroimaging and clinical professionals.

Project description:Parkinson's disease (PD) is a progressive neurodegenerative disorder that affects extensive regions of the central nervous system. In this work, we evaluated the structural connectome of patients with PD, as mapped by diffusion-weighted MRI tractography and a multi-shell, multi-tissue (MSMT) constrained spherical deconvolution (CSD) method to increase the precision of tractography at tissue interfaces. The connectome was mapped with probabilistic MSMT-CSD in 21 patients with PD and in 21 age- and gender-matched controls. Mapping was also performed by deterministic single-shell, single tissue (SSST)-CSD tracking and probabilistic SSST-CSD tracking for comparison. A support vector machine was trained to predict diagnosis based on a linear combination of graph metrics. We showed that probabilistic MSMT-CSD could detect significantly reduced global strength, efficiency, clustering, and small-worldness, and increased global path length in patients with PD relative to healthy controls; by contrast, probabilistic SSST-CSD only detected the difference in global strength and small-worldness. In patients with PD, probabilistic MSMT-CSD also detected a significant reduction in local efficiency and detected clustering in the motor, frontal temporoparietal associative, limbic, basal ganglia, and thalamic areas. The network-based statistic identified a subnetwork of reduced connectivity by MSMT-CSD and probabilistic SSST-CSD in patients with PD, involving key components of the cortico-basal ganglia-thalamocortical network. Finally, probabilistic MSMT-CSD had superior diagnostic accuracy compared with conventional probabilistic SSST-CSD and deterministic SSST-CSD tracking. In conclusion, probabilistic MSMT-CSD detected a greater extent of connectome pathology in patients with PD, including those with cortico-basal ganglia-thalamocortical network disruptions. Connectome analysis based on probabilistic MSMT-CSD may be useful when evaluating the extent of white matter connectivity disruptions in PD.

Dataset Information

Spherical Deconvolution of Multichannel Diffusion MRI Data with Non-Gaussian Noise Models and Spatial Regularization.

Publications

Spherical Deconvolution of Multichannel Diffusion MRI Data with Non-Gaussian Noise Models and Spatial Regularization.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets