Spherical Deconvolution of MultichannelDiffusion MRI Data with Non-Gaussian NoiseModels and Spatial Regularization
View/ Open
Date
2015-10-15Author
Canales Rodríguez, Erick Jorge
Daducci, Alessandro
Sotiropoulos, Stamatios N.
Caruyer, Emmanuel
Aja Fernández, Santiago
Radua, Joaquim
Yurramendi Mendizabal, Yosu
Iturria Medina, Yasser
Melie García, Lester
Alemán Gómez, Yasser
Thiran, Jean-Philippe
Sarró, Salvador
Pomarol-Clotet, Edith
Salvador, Raymond
Metadata
Show full item record
PLOS ONE 10(10) : (2015) // Article ID e0138910
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.
Collections
Except where otherwise noted, this item's license is described as 2015 Canales-Rodríguez et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited