Sparse and Optimal Acquisition Design for Diffusion MRI and Beyond

PURPOSE Diffusion magnetic resonance imaging (MRI) in combination with functional MRI promises a whole new vista for scientists to investigate noninvasively the structural and functional connectivity of the human brain-the human connectome, which had heretofore been out of reach. As with other imaging modalities, diffusion MRI data are inherently noisy and its acquisition time-consuming. Further, a faithful representation of the human connectome that can serve as a predictive model requires a robust and accurate data-analytic pipeline. The focus of this paper is on one of the key segments of this pipeline-in particular, the development of a sparse and optimal acquisition (SOA) design for diffusion MRI multiple-shell acquisition and beyond. METHODS The authors propose a novel optimality criterion for sparse multiple-shell acquisition and quasimultiple-shell designs in diffusion MRI and a novel and effective semistochastic and moderately greedy combinatorial search strategy with simulated annealing to locate the optimum design or configuration. The goal of the optimality criteria is threefold: first, to maximize uniformity of the diffusion measurements in each shell, which is equivalent to maximal incoherence in angular measurements; second, to maximize coverage of the diffusion measurements around each radial line to achieve maximal incoherence in radial measurements for multiple-shell acquisition; and finally, to ensure maximum uniformity of diffusion measurement directions in the limiting case when all the shells are coincidental as in the case of a single-shell acquisition. The approach taken in evaluating the stability of various acquisition designs is based on the condition number and the A-optimal measure of the design matrix. RESULTS Even though the number of distinct configurations for a given set of diffusion gradient directions is very large in general-e.g., in the order of 10(232) for a set of 144 diffusion gradient directions, the proposed search strategy was found to be effective in finding the optimum configuration. It was found that the square design is the most robust (i.e., with stable condition numbers and A-optimal measures under varying experimental conditions) among many other possible designs of the same sample size. Under the same performance evaluation, the square design was found to be more robust than the widely used sampling schemes similar to that of 3D radial MRI and of diffusion spectrum imaging (DSI). CONCLUSIONS A novel optimality criterion for sparse multiple-shell acquisition and quasimultiple-shell designs in diffusion MRI and an effective search strategy for finding the best configuration have been developed. The results are very promising, interesting, and practical for diffusion MRI acquisitions.

[1]  J Hennig,et al.  RARE imaging: A fast imaging method for clinical MR , 1986, Magnetic resonance in medicine.

[2]  B. Hermann,et al.  Children with new-onset epilepsy exhibit diffusion abnormalities in cerebral white matter in the absence of volumetric differences , 2010, Epilepsy Research.

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

[4]  Rachid Deriche,et al.  Multiple q-shell diffusion propagator imaging , 2011, Medical Image Anal..

[5]  Charles H Cunningham,et al.  Temporal stability of adaptive 3D radial MRI using multidimensional golden means , 2009, Magnetic resonance in medicine.

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

[7]  Carl-Fredrik Westin,et al.  Sparse Multi-Shell Diffusion Imaging , 2011, MICCAI.

[8]  D L Parker,et al.  Comparison of gradient encoding schemes for diffusion‐tensor MRI , 2001, Journal of magnetic resonance imaging : JMRI.

[9]  James E. pLebensohn Geometry and the Imagination , 1952 .

[10]  P. Basser,et al.  MR diffusion tensor spectroscopy and imaging. , 1994, Biophysical journal.

[11]  M. Horsfield,et al.  Optimal strategies for measuring diffusion in anisotropic systems by magnetic resonance imaging , 1999, Magnetic resonance in medicine.

[12]  Peter Börnert,et al.  Three‐dimensional radial ultrashort echo‐time imaging with T2 adapted sampling , 2006, Magnetic resonance in medicine.

[13]  P. Basser,et al.  A unifying theoretical and algorithmic framework for least squares methods of estimation in diffusion tensor imaging. , 2006, Journal of magnetic resonance.

[14]  S Skare,et al.  Condition number as a measure of noise performance of diffusion tensor data acquisition schemes with MRI. , 2000, Journal of magnetic resonance.

[15]  Cheng Guan Koay,et al.  Analytically exact spiral scheme for generating uniformly distributed points on the unit sphere , 2011, J. Comput. Sci..

[16]  T. Hendler,et al.  High b‐value q‐space analyzed diffusion‐weighted MRI: Application to multiple sclerosis , 2002, Magnetic resonance in medicine.

[17]  Wolfhard Semmler,et al.  Sodium MRI using a density‐adapted 3D radial acquisition technique , 2009, Magnetic resonance in medicine.

[18]  Daniel C Alexander,et al.  Optimal acquisition orders of diffusion‐weighted MRI measurements , 2007, Journal of magnetic resonance imaging : JMRI.

[19]  T D Nguyen,et al.  Optimization of view ordering for motion artifact suppression. , 2001, Magnetic resonance imaging.

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

[21]  C-M Lait,et al.  True three-dimensional image reconstruction by nuclear magnetic resonance zeugmatography , 1981 .

[22]  D. Tank,et al.  Brain magnetic resonance imaging with contrast dependent on blood oxygenation. , 1990, Proceedings of the National Academy of Sciences of the United States of America.

[23]  James G Pipe,et al.  Generation and visualization of four-dimensional MR angiography data using an undersampled 3-D projection trajectory , 2006, IEEE Transactions on Medical Imaging.

[24]  Cheng Guan Koay,et al.  A signal transformational framework for breaking the noise floor and its applications in MRI. , 2009, Journal of magnetic resonance.

[25]  M. Kraut,et al.  Comparison of weakness progression in inclusion body myositis during treatment with methotrexate or placebo , 2002, Annals of neurology.

[26]  Carlo Pierpaoli,et al.  The Elliptical Cone of Uncertainty and Its Normalized Measures in Diffusion Tensor Imaging , 2008, IEEE Transactions on Medical Imaging.

[27]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[28]  Lawrence Dougherty,et al.  Dynamic MRI with projection reconstruction and KWIC processing for simultaneous high spatial and temporal resolution , 2004, Magnetic resonance in medicine.

[29]  P. Basser,et al.  Simple harmonic oscillator based reconstruction and estimation for three-dimensional q-space MRI , 2009 .

[30]  W. J. Studden,et al.  Theory Of Optimal Experiments , 1972 .

[31]  J Velikina,et al.  Highly constrained backprojection for time‐resolved MRI , 2006, Magnetic resonance in medicine.

[32]  D. Beste Design of Satellite Constellations for Optimal Continuous Coverage , 1978, IEEE Transactions on Aerospace and Electronic Systems.

[33]  O Dössel,et al.  Fast isotropic volumetric coronary MR angiography using free‐breathing 3D radial balanced FFE acquisition , 2004, Magnetic resonance in medicine.

[34]  R. Busse,et al.  Fast spin echo sequences with very long echo trains: Design of variable refocusing flip angle schedules and generation of clinical T2 contrast , 2006, Magnetic resonance in medicine.

[35]  Maria I Altbach,et al.  View‐ordering in radial fast spin‐echo imaging , 2004, Magnetic resonance in medicine.

[36]  J. Dubois,et al.  Optimized diffusion gradient orientation schemes for corrupted clinical DTI data sets , 2006, Magnetic Resonance Materials in Physics, Biology and Medicine.

[37]  Carlo Pierpaoli,et al.  Error Propagation Framework for Diffusion Tensor Imaging via Diffusion Tensor Representations , 2007, IEEE Transactions on Medical Imaging.

[38]  Carlo Pierpaoli,et al.  Probabilistic Identification and Estimation of Noise (piesno): a Self-consistent Approach and Its Applications in Mri , 2009 .

[39]  D. Donoho,et al.  Sparse MRI: The application of compressed sensing for rapid MR imaging , 2007, Magnetic resonance in medicine.

[40]  Olaf Sporns,et al.  The Human Connectome: A Structural Description of the Human Brain , 2005, PLoS Comput. Biol..

[41]  Mario Bertero,et al.  Introduction to Inverse Problems in Imaging , 1998 .

[42]  Andrew L. Alexander,et al.  Hybrid diffusion imaging , 2007, NeuroImage.

[43]  Stephen Wolfram,et al.  The Mathematica Book , 1996 .

[44]  Walter F Block,et al.  Time‐resolved contrast‐enhanced imaging with isotropic resolution and broad coverage using an undersampled 3D projection trajectory , 2002, Magnetic resonance in medicine.

[45]  Andrew L. Alexander,et al.  Diffusion tensor imaging of the corpus callosum in Autism , 2007, NeuroImage.

[46]  Christian Windischberger,et al.  Toward discovery science of human brain function , 2010, Proceedings of the National Academy of Sciences.

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

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

[49]  P. Grenier,et al.  MR imaging of intravoxel incoherent motions: application to diffusion and perfusion in neurologic disorders. , 1986, Radiology.

[50]  A. D. de Crespigny,et al.  Compromised white matter tract integrity in schizophrenia inferred from diffusion tensor imaging. , 1999, Archives of general psychiatry.

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

[52]  S J Riederer,et al.  Fluoroscopically triggered contrast-enhanced three-dimensional MR angiography with elliptical centric view order: application to the renal arteries. , 1997, Radiology.

[53]  K. Scheffler,et al.  Multiecho sequences with variable refocusing flip angles: Optimization of signal behavior using smooth transitions between pseudo steady states (TRAPS) , 2003, Magnetic resonance in medicine.

[54]  V. Wedeen,et al.  Mapping fiber orientation spectra in cerebral white matter with Fourier-transform diffusion MRI , 2000 .

[55]  Kaleem Siddiqi,et al.  Recent advances in diffusion MRI modeling: Angular and radial reconstruction , 2011, Medical Image Anal..

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

[57]  Richard Kijowski,et al.  Effects of refocusing flip angle modulation and view ordering in 3D fast spin echo , 2008, Magnetic resonance in medicine.

[58]  R F Busse,et al.  A flexible view ordering technique for high‐quality real‐time 2DFT MR fluoroscopy , 1999, Magnetic resonance in medicine.

[59]  M. Meyerand,et al.  Extremely efficient and deterministic approach to generating optimal ordering of diffusion MRI measurements. , 2011, Medical physics.

[60]  P. Basser,et al.  Diffusion tensor MR imaging of the human brain. , 1996, Radiology.

[61]  H. Pfeifer Principles of Nuclear Magnetic Resonance Microscopy , 1992 .

[62]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[63]  Cheng Guan Koay,et al.  A simple scheme for generating nearly uniform distribution of antipodally symmetric points on the unit sphere , 2011, J. Comput. Sci..