A Closed-Form Solution of Rotation Invariant Spherical Harmonic Features in Diffusion MRI

Rotation invariant features are an indispensable tool for characterizing diffusion Magnetic Resonance Imaging (MRI) and in particular for brain tissue microstructure estimation. In this work, we propose a new mathematical framework for efficiently calculating a complete set of such invariants from any spherical function. Specifically, our method is based on the spherical harmonics series expansion of a given function of any order and can be applied directly to the resulting coefficients by performing a simple integral operation analytically. This enable us to derive a general closed-form equation for the invariants. We test our invariants on the diffusion MRI fiber orientation distribution function obtained from the diffusion signal both in-vivo and in synthetic data. Results show how it is possible to use these invariants for characterizing the white matter using a small but complete set of features.

[1]  Gian-Carlo Rota,et al.  Apolarity and Canonical Forms for Homogeneous Polynomials , 1993, Eur. J. Comb..

[2]  J. Veraart,et al.  Mapping orientational and microstructural metrics of neuronal integrity with in vivo diffusion MRI , 2016, 1609.09144.

[3]  Jacob T. Schwartz,et al.  Fast Probabilistic Algorithms for Verification of Polynomial Identities , 1980, J. ACM.

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

[5]  Giuseppe Scotti,et al.  A modified damped Richardson–Lucy algorithm to reduce isotropic background effects in spherical deconvolution , 2010, NeuroImage.

[6]  Rachid Deriche,et al.  Complete Set of Invariants of a 4 th Order Tensor: The 12 Tasks of HARDI from Ternary Quartics , 2014, MICCAI.

[7]  Risi Kondor,et al.  A novel set of rotationally and translationally invariant features for images based on the non-commutative bispectrum , 2007, cs/0701127.

[8]  Steen Moeller,et al.  Advances in diffusion MRI acquisition and processing in the Human Connectome Project , 2013, NeuroImage.

[9]  Pedro M. Q. Aguiar,et al.  Shape representation via elementary symmetric polynomials: A complete invariant inspired by the bispectrum , 2013, 2013 IEEE International Conference on Image Processing.

[10]  Jelle Veraart,et al.  Rotationally-invariant mapping of scalar and orientational metrics of neuronal microstructure with diffusion MRI , 2018, NeuroImage.

[11]  M. Descoteaux,et al.  A Generalized SMT-Based Framework for Diffusion MRI Microstructural Model Estimation , 2018, MICCAI 2018.

[12]  René Vidal,et al.  Rotation Invariant Features for HARDI , 2013, IPMI.

[13]  Jürgen Hennig,et al.  Disentangling micro from mesostructure by diffusion MRI: A Bayesian approach , 2017, NeuroImage.

[14]  Daniel C. Alexander,et al.  Multi-compartment microscopic diffusion imaging , 2016, NeuroImage.

[15]  Ramakrishna Kakarala,et al.  A theory of phase-sensitive rotation invariance with spherical harmonic and moment-based representations , 2010, 2010 IEEE Computer Society Conference on Computer Vision and Pattern Recognition.

[16]  Luke Bloy,et al.  Demons registration of high angular resolution diffusion images , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[17]  Herbert H. H. Homeier,et al.  Some Properties of the Coupling Coefficients of Real Spherical Harmonics and Their Relation to Gaunt Coefficients , 1996 .

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

[19]  Ragini Verma,et al.  On facilitating the use of HARDI in population studies by creating rotation-invariant markers , 2015, Medical Image Anal..

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

[21]  Richard Zippel,et al.  Probabilistic algorithms for sparse polynomials , 1979, EUROSAM.

[22]  Rachid Deriche,et al.  Generalized Invariants of a 4th order tensor: Building blocks for new biomarkers in dMRI , 2012 .