A Riemannian approach for estimating orientation distribution function (ODF) images from high-angular resolution diffusion imaging (HARDI)

High-angular resolution diffusion imaging (HARDI) is a magnetic resonance technique estimating the direction of self-diffusion of water molecules in biological tissue. HARDI encodes at each pixel (voxel) the orientation distribution function (ODF) of water diffusion molecules, i.e. the probability distribution function of finding a water molecule which moved in a certain direction during the observation time. As a consequence ODF images differ from usual gray scale images with respect to their underlying geometry as well as with respect to their error distribution. We present a Bayesian estimator for ODF images considering these differences. To this end, we derive a likelihood function based on the Rician distribution of the NMR signals and propose prior distributions considering ODFs as Riemanian manifolds. Utilizing properties of spherical harmonics and the square root representation of ODFs allows us to effectively reconstruct and regularize ODF images in one step within this Riemannian framework. Experiments demonstrate the merits of our approach on synthetic as well as on real data.

[1]  Hanno Scharr,et al.  Riemannian Bayesian estimation of diffusion tensor images , 2009, 2009 IEEE 12th International Conference on Computer Vision.

[2]  Alan Connelly,et al.  Robust determination of the fibre orientation distribution in diffusion MRI: Non-negativity constrained super-resolved spherical deconvolution , 2007, NeuroImage.

[3]  Nikos Paragios,et al.  A Convex Semi-definite Positive Framework for DTI Estimation and Regularization , 2007, ISVC.

[4]  R. Vidal,et al.  A nonparametric Riemannian framework for processing high angular resolution diffusion images (HARDI) , 2009, 2009 IEEE Conference on Computer Vision and Pattern Recognition.

[5]  Christophe Lenglet,et al.  ODF reconstruction in q-ball imaging with solid angle consideration , 2009, 2009 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[6]  Rachid Deriche,et al.  A Riemannian Framework for Orientation Distribution Function Computing , 2009, MICCAI.

[7]  R. Deriche,et al.  Regularized, fast, and robust analytical Q‐ball imaging , 2007, Magnetic resonance in medicine.

[8]  Baba C. Vemuri,et al.  A Unified Computational Framework for Deconvolution to Reconstruct Multiple Fibers From Diffusion Weighted MRI , 2007, IEEE Transactions on Medical Imaging.

[9]  N. Makris,et al.  High angular resolution diffusion imaging reveals intravoxel white matter fiber heterogeneity , 2002, Magnetic resonance in medicine.

[10]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[11]  D. Tuch.  Q‐ball imaging , 2004, Magnetic resonance in medicine.

[12]  A. R. Edmonds.  Angular Momentum in Quantum Mechanics , 1957 .

[13]  Anuj Srivastava,et al.  Riemannian Analysis of Probability Density Functions with Applications in Vision , 2007, 2007 IEEE Conference on Computer Vision and Pattern Recognition.

[14]  Luc Brun,et al.  Robust Variational Estimation of PDF functions from Diffusion MR Signal , 2008 .

[15]  Christophe Lenglet,et al.  Estimating Orientation Distribution Functions with Probability Density Constraints and Spatial Regularity , 2009, MICCAI.