Spatial HARDI: Improved visualization of complex white matter architecture with Bayesian spatial regularization

Imaging of water diffusion using magnetic resonance imaging has become an important noninvasive method for probing the white matter connectivity of the human brain for scientific and clinical studies. Current methods, such as diffusion tensor imaging (DTI), high angular resolution diffusion imaging (HARDI) such as q-ball imaging, and diffusion spectrum imaging (DSI), are limited by low spatial resolution, long scan times, and low signal-to-noise ratio (SNR). These methods fundamentally perform reconstruction on a voxel-by-voxel level, effectively discarding the natural coherence of the data at different points in space. This paper attempts to overcome these tradeoffs by using spatial information to constrain the reconstruction from raw diffusion MRI data, and thereby improve angular resolution and noise tolerance. Spatial constraints are specified in terms of a prior probability distribution, which is then incorporated in a Bayesian reconstruction formulation. By taking the log of the resulting posterior distribution, optimal Bayesian reconstruction is reduced to a cost minimization problem. The minimization is solved using a new iterative algorithm based on successive least squares quadratic descent. Simulation studies and in vivo results are presented which indicate significant gains in terms of higher angular resolution of diffusion orientation distribution functions, better separation of crossing fibers, and improved reconstruction SNR over the same HARDI method, spherical harmonic q-ball imaging, without spatial regularization. Preliminary data also indicate that the proposed method might be better at maintaining accurate ODFs for smaller numbers of diffusion-weighted acquisition directions (hence faster scans) compared to conventional methods. Possible impacts of this work include improved evaluation of white matter microstructural integrity in regions of crossing fibers and higher spatial and angular resolution for more accurate tractography.

[1]  Alan Connelly,et al.  Direct estimation of the fiber orientation density function from diffusion-weighted MRI data using spherical deconvolution , 2004, NeuroImage.

[2]  Derek K. Jones,et al.  The effect of gradient sampling schemes on measures derived from diffusion tensor MRI: A Monte Carlo study † , 2004, Magnetic resonance in medicine.

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

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

[5]  Baba C. Vemuri,et al.  Resolution of complex tissue microarchitecture using the diffusion orientation transform (DOT) , 2006, NeuroImage.

[6]  Kalvis M. Jansons,et al.  Persistent angular structure: new insights from diffusion magnetic resonance imaging data , 2003 .

[7]  P. Thiran,et al.  Mapping Human Whole-Brain Structural Networks with Diffusion MRI , 2007, PloS one.

[8]  R. Zabih,et al.  Bayesian parallel imaging with edge‐preserving priors , 2007, Magnetic resonance in medicine.

[9]  T. Mareci,et al.  Generalized diffusion tensor imaging and analytical relationships between diffusion tensor imaging and high angular resolution diffusion imaging , 2003, Magnetic resonance in medicine.

[10]  Mark J. Lowe,et al.  An objective method for regularization of fiber orientation distributions derived from diffusion-weighted MRI , 2007, NeuroImage.

[11]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

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

[13]  M. Moseley,et al.  Magnetic Resonance in Medicine 51:924–937 (2004) Characterizing Non-Gaussian Diffusion by Using Generalized Diffusion Tensors , 2022 .

[14]  William H. Press,et al.  Numerical recipes in C , 2002 .

[15]  C. Beaulieu,et al.  The basis of anisotropic water diffusion in the nervous system – a technical review , 2002, NMR in biomedicine.

[16]  Wim Vanduffel,et al.  Q-ball imaging of macaque white matter architecture , 2005, Philosophical Transactions of the Royal Society B: Biological Sciences.

[17]  P. Boesiger,et al.  SENSE: Sensitivity encoding for fast MRI , 1999, Magnetic resonance in medicine.

[18]  Stephen M. Smith,et al.  Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm , 2001, IEEE Transactions on Medical Imaging.

[19]  H. Groemer.  Geometric Applications of Fourier Series and Spherical Harmonics , 1996 .

[20]  R. Gerchberg.  A practical algorithm for the determination of phase from image and diffraction plane pictures , 1972 .

[21]  Alan Barnett,et al.  Theory of Q‐ball imaging redux: Implications for fiber tracking , 2009, Magnetic resonance in medicine.

[22]  Yihong Yang,et al.  Mapping the orientation of intravoxel crossing fibers based on the phase information of diffusion circular spectrum , 2004, NeuroImage.

[23]  K. Kwong,et al.  Parallel imaging reconstruction using automatic regularization , 2004, Magnetic resonance in medicine.

[24]  S. Arridge,et al.  Detection and modeling of non‐Gaussian apparent diffusion coefficient profiles in human brain data , 2002, Magnetic resonance in medicine.

[25]  P. V. van Zijl,et al.  Three‐dimensional tracking of axonal projections in the brain by magnetic resonance imaging , 1999, Annals of neurology.

[26]  Duan Xu,et al.  Q‐ball reconstruction of multimodal fiber orientations using the spherical harmonic basis , 2006, Magnetic resonance in medicine.

[27]  Abbas F. Sadikot,et al.  Flow-based fiber tracking with diffusion tensor and q-ball data: Validation and comparison to principal diffusion direction techniques , 2005, NeuroImage.

[28]  Erick Jorge Canales-Rodríguez,et al.  Mathematical description of q‐space in spherical coordinates: Exact q‐ball imaging , 2009, Magnetic resonance in medicine.

[29]  Pratik Mukherjee,et al.  Development and initial evaluation of 7-T q-ball imaging of the human brain. , 2008, Magnetic resonance imaging.

[30]  Rachid Deriche,et al.  Impact of Rician Adapted Non-Local Means Filtering on HARDI , 2008, MICCAI.

[31]  Zhi-Pei Liang,et al.  Anatomically constrained reconstruction from noisy data , 2008, Magnetic resonance in medicine.

[32]  Paul M. Thompson,et al.  HARDI Denoising: Variational Regularization of the Spherical Apparent Diffusion Coefficient sADC , 2009, IPMI.

[33]  D. LeBihan,et al.  Validation of q-ball imaging with a diffusion fibre-crossing phantom on a clinical scanner , 2005, Philosophical Transactions of the Royal Society B: Biological Sciences.

[34]  Guy B. Williams,et al.  Inference of multiple fiber orientations in high angular resolution diffusion imaging , 2005, Magnetic resonance in medicine.

[35]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[36]  Baba C. Vemuri,et al.  A novel tensor distribution model for the diffusion-weighted MR signal , 2007, NeuroImage.

[37]  Carl-Fredrik Westin,et al.  Estimation of fiber Orientation Probability Density Functions in High Angular Resolution Diffusion Imaging , 2009, NeuroImage.

[38]  A. Anderson.  Measurement of fiber orientation distributions using high angular resolution diffusion imaging , 2005, Magnetic resonance in medicine.

[39]  L. Frank.  Characterization of anisotropy in high angular resolution diffusion‐weighted MRI , 2002, Magnetic resonance in medicine.

[40]  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.