Inter-Scanner Harmonization of High Angular Resolution DW-MRI using Null Space Deep Learning

Diffusion-weighted magnetic resonance imaging (DW-MRI) allows for non-invasive imaging of the local fiber architecture of the human brain at a millimetric scale. Multiple classical approaches have been proposed to detect both single (e.g., tensors) and multiple (e.g., constrained spherical deconvolution, CSD) fiber population orientations per voxel. However, existing techniques generally exhibit low reproducibility across MRI scanners. Herein, we propose a data-driven tech-nique using a neural network design which exploits two categories of data. First, training data were acquired on three squirrel monkey brains using ex-vivo DW-MRI and histology of the brain. Second, repeated scans of human subjects were acquired on two different scanners to augment the learning of the network pro-posed. To use these data, we propose a new network architecture, the null space deep network (NSDN), to simultaneously learn on traditional observed/truth pairs (e.g., MRI-histology voxels) along with repeated observations without a known truth (e.g., scan-rescan MRI). The NSDN was tested on twenty percent of the histology voxels that were kept completely blind to the network. NSDN significantly improved absolute performance relative to histology by 3.87% over CSD and 1.42% over a recently proposed deep neural network approach. More-over, it improved reproducibility on the paired data by 21.19% over CSD and 10.09% over a recently proposed deep approach. Finally, NSDN improved gen-eralizability of the model to a third in vivo human scanner (which was not used in training) by 16.08% over CSD and 10.41% over a recently proposed deep learn-ing approach. This work suggests that data-driven approaches for local fiber re-construction are more reproducible, informative and precise and offers a novel, practical method for determining these models.

[1]  G. Ball,et al.  Voxel-wise comparisons of cellular microstructure and diffusion-MRI in mouse hippocampus using 3D Bridging of Optically-clear histology with Neuroimaging Data (3D-BOND) , 2018, Scientific Reports.

[2]  Bennett A. Landman,et al.  Comparison of multi-fiber reproducibility of PAS-MRI and Q-ball with empirical multiple b-value HARDI , 2017, Medical Imaging.

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

[4]  Denis Le Bihan,et al.  Looking into the functional architecture of the brain with diffusion MRI , 2003, Nature Reviews Neuroscience.

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

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

[7]  F. Calamante,et al.  How many diffusion gradient directions are required for HARDI , 2009 .

[8]  Camilo Bermudez,et al.  Deep Learning Captures More Accurate Diffusion Fiber Orientations Distributions than Constrained Spherical Deconvolution , 2019, ArXiv.

[9]  Stefan Skare,et al.  How to correct susceptibility distortions in spin-echo echo-planar images: application to diffusion tensor imaging , 2003, NeuroImage.

[10]  Uwe Klose,et al.  Reproducibility and consistency of evaluation techniques for HARDI data , 2009, Magnetic Resonance Materials in Physics, Biology and Medicine.

[11]  J. Alger,et al.  Between-Scanner and Between-Visit Variation in Normal White Matter Apparent Diffusion Coefficient Values in the Setting of a Multi-Center Clinical Trial , 2016, Clinical Neuroradiology.

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

[13]  Geoffrey E. Hinton,et al.  ImageNet classification with deep convolutional neural networks , 2012, Commun. ACM.

[14]  Stamatios N. Sotiropoulos,et al.  An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging , 2016, NeuroImage.

[15]  Prasanna Parvathaneni,et al.  SHARD: spherical harmonic-based robust outlier detection for HARDI methods , 2018, Medical Imaging.

[16]  Jessica A. Turner,et al.  Multi-site study of diffusion metric variability: effects of site, vendor, field strength, and echo time on regions-of-interest and histogram-bin analyses , 2016, SPIE Medical Imaging.

[17]  Ragini Verma,et al.  Harmonization of multi-site diffusion tensor imaging data , 2017, NeuroImage.

[18]  Peter Savadjiev,et al.  Inter-site and inter-scanner diffusion MRI data harmonization , 2016, NeuroImage.

[19]  R. Deriche,et al.  Apparent diffusion coefficients from high angular resolution diffusion imaging: Estimation and applications , 2006, Magnetic resonance in medicine.

[20]  Bennett A. Landman,et al.  Comparison of 3D orientation distribution functions measured with confocal microscopy and diffusion MRI , 2016, NeuroImage.

[21]  Peter Savadjiev,et al.  Harmonizing Diffusion MRI Data Across Multiple Sites and Scanners , 2015, MICCAI.

[22]  Stephen M Smith,et al.  Fast robust automated brain extraction , 2002, Human brain mapping.

[23]  Wu Liu,et al.  A discriminative null space based deep learning approach for person re-identification , 2016, 2016 4th International Conference on Cloud Computing and Intelligence Systems (CCIS).

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

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