What If Each Voxel Were Measured With a Different Diffusion Protocol?
Authors
Affiliations (3)
Affiliations (3)
- Center for Biomedical Imaging and Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, New York University School of Medicine, New York, NY, USA.
- GE HealthCare Technology and Innovation Center, Niskayuna, New York, NY, USA.
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, South Carolina, USA.
Abstract
Expansion of diffusion MRI (dMRI) both into the realm of strong gradients and into accessible imaging with portable low-field devices brings about the challenge of gradient nonlinearities. Spatial variations of the diffusion gradients make diffusion weightings and directions non-uniform across the field of view, and deform perfect shells in the <math xmlns="http://www.w3.org/1998/Math/MathML"> <semantics><mrow><mi>q</mi></mrow> <annotation>$$ q $$</annotation></semantics> </math> -space designed for isotropic directional coverage. Such imperfections hinder parameter estimation: Anisotropic shells hamper the deconvolution of the fiber orientation distribution function (fODF), while brute-force retraining of a nonlinear regressor for each unique set of directions and diffusion weightings is computationally inefficient. Here, we propose a protocol-independent parameter estimation (PIPE) method that enables fast parameter estimation for the most general case where each voxel is measured with a different protocol in <math xmlns="http://www.w3.org/1998/Math/MathML"> <semantics><mrow><mi>q</mi></mrow> <annotation>$$ q $$</annotation></semantics> </math> -space. PIPE applies to any spherical convolution-based dMRI model, irrespective of its complexity, which makes it suitable both for white and gray matter in the brain or spinal cord, and for other tissues where fiber bundles have the same properties (fiber response) within a voxel, but are distributed with an arbitrary fODF. We also derive a parsimonious representation that isolates isotropic and anisotropic effects of gradient nonlinearities on multidimensional diffusion encodings. Applied to in vivo human MRI with linear tensor encoding on a high-performance gradient system, PIPE evaluates fiber response and fODF parameters for the whole brain in the presence of significant gradient nonlinearities in under 3 min. PIPE enables fast parameter estimation in the presence of arbitrary gradient nonlinearities, eliminating the need to arrange dMRI in shells or to retrain the estimator for different protocols in each voxel. PIPE applies to any model based on a convolution of a voxel-wise fiber response and fODF, and data from varying <math xmlns="http://www.w3.org/1998/Math/MathML"> <semantics><mrow><mi>b</mi></mrow> <annotation>$$ b $$</annotation></semantics> </math> -values, diffusion/echo times, and other scan parameters.