Fitting IVIM with Variable Projection and Simplicial Optimization

Fitting multi-exponential models to Diffusion MRI (dMRI) data has always been challenging due to various underlying complexities. In this work, we introduce a novel and robust fitting framework for the standard two-compartment IVIM microstructural model. This framework provides a significant improvement over the existing methods and helps estimate the associated diffusion and perfusion parameters of IVIM in an automatic manner. As a part of this work we provide capabilities to switch between more advanced global optimization methods such as simplicial homology (SH) and differential evolution (DE). Our experiments show that the results obtained from this simultaneous fitting procedure disentangle the model parameters in a reduced subspace. The proposed framework extends the seminal work originated in the MIX framework, with improved procedures for multi-stage fitting. This framework has been made available as an open-source Python implementation and disseminated to the community through the DIPY project.

[1]  Eric Jones,et al.  SciPy: Open Source Scientific Tools for Python , 2001 .

[2]  Carl Sandrock,et al.  A simplicial homology algorithm for Lipschitz optimisation , 2018, Journal of Global Optimization.

[3]  Joe D. Warren,et al.  The program dependence graph and its use in optimization , 1984, TOPL.

[4]  Bingsheng Huang,et al.  Intravoxel incoherent motion MR imaging in nasopharyngeal carcinoma: comparison and correlation with dynamic contrast enhanced MR imaging , 2017, Oncotarget.

[5]  Stephen P. Boyd,et al.  CVXPY: A Python-Embedded Modeling Language for Convex Optimization , 2016, J. Mach. Learn. Res..

[6]  Maryam Afzali,et al.  MicroLearn : Framework for machine learning , reconstruction , optimization and microstructure modeling , 2018 .

[7]  Richard H. Byrd,et al.  Approximate solution of the trust region problem by minimization over two-dimensional subspaces , 1988, Math. Program..

[8]  Rachid Deriche,et al.  The Dmipy Toolbox: Diffusion MRI Multi-Compartment Modeling and Microstructure Recovery Made Easy , 2019, Front. Neuroinform..

[9]  R. Storn,et al.  Differential Evolution: A Practical Approach to Global Optimization (Natural Computing Series) , 2005 .

[10]  T. Georgiou,et al.  Microstructure Imaging of Crossing (MIX) White Matter Fibers from diffusion MRI , 2016, Scientific Reports.

[11]  Nikolaos V. Sahinidis,et al.  Derivative-free optimization: a review of algorithms and comparison of software implementations , 2013, J. Glob. Optim..

[12]  Denis Le Bihan,et al.  What can we see with IVIM MRI? , 2017, NeuroImage.

[13]  J. Nocedal,et al.  A Limited Memory Algorithm for Bound Constrained Optimization , 1995, SIAM J. Sci. Comput..

[14]  A. Istratov,et al.  Exponential analysis in physical phenomena , 1999 .

[15]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[16]  G. Golub,et al.  Separable nonlinear least squares: the variable projection method and its applications , 2003 .

[17]  Maxime Descoteaux,et al.  Dipy, a library for the analysis of diffusion MRI data , 2014, Front. Neuroinform..

[18]  P. T. While,et al.  A comparative simulation study of bayesian fitting approaches to intravoxel incoherent motion modeling in diffusion‐weighted MRI , 2017, Magnetic resonance in medicine.

[19]  Mats Andersson,et al.  Comparison of methods for estimation of the intravoxel incoherent motion (IVIM) diffusion coefficient (D) and perfusion fraction (f) , 2018, Magnetic Resonance Materials in Physics, Biology and Medicine.

[20]  Thomas F. Coleman,et al.  A Subspace, Interior, and Conjugate Gradient Method for Large-Scale Bound-Constrained Minimization Problems , 1999, SIAM J. Sci. Comput..

[21]  Namkug Kim,et al.  Histogram Analysis of Intravoxel Incoherent Motion for Differentiating Recurrent Tumor from Treatment Effect in Patients with Glioblastoma: Initial Clinical Experience , 2014, American Journal of Neuroradiology.

[22]  R. B. Bapat,et al.  A constructive proof of a permutation-based generalization of Sperner's lemma , 1989, Math. Program..

[23]  Masako Kataoka,et al.  IVIM MRI of the Breast , 2018, Intravoxel Incoherent Motion (IVIM) MRI.

[24]  D. Le Bihan,et al.  Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. , 1988, Radiology.

[25]  G. Golub,et al.  Linear least squares and quadratic programming , 1969 .

[26]  Herbert R. Halvorson.  Pade-Laplace algorithm for sums of exponentials Selecting appropriate exponential model and initial estimates for exponential fitting , 1992 .

[27]  Martijn Froeling,et al.  Comparison of six fit algorithms for the intra-voxel incoherent motion model of diffusion-weighted magnetic resonance imaging data of pancreatic cancer patients , 2018, PloS one.

[28]  Remigijus Paulavičius,et al.  Simplicial Global Optimization , 2014 .

[29]  Jing Yuan,et al.  Liver intravoxel incoherent motion (IVIM) magnetic resonance imaging: a comprehensive review of published data on normal values and applications for fibrosis and tumor evaluation. , 2017, Quantitative imaging in medicine and surgery.

[30]  Utaroh Motosugi,et al.  Diffusion and Intravoxel Incoherent Motion MR Imaging-based Virtual Elastography: A Hypothesis-generating Study in the Liver. , 2017, Radiology.

[31]  Berthold Kiefer,et al.  Intravoxel incoherent motion diffusion-weighted MR imaging for characterization of focal pancreatic lesions. , 2013, Radiology.

[32]  Vimal Bhatia,et al.  The Diffusion-KLMS Algorithm , 2014, 2014 International Conference on Information Technology.

[33]  I. Sobol.  On the distribution of points in a cube and the approximate evaluation of integrals , 1967 .

[34]  Idris Ismail,et al.  Comparative Study of Meta-heuristics Optimization Algorithm using Benchmark Function , 2017 .

[35]  Earl L. Bell.  FITTING MULTI-COMPONENT EXPONENTIAL DECAY CURVES BY DIGITAL COMPUTER , 1965 .

[36]  R. Meuli,et al.  Perfusion Measurement in Brain Gliomas with Intravoxel Incoherent Motion MRI , 2014, American Journal of Neuroradiology.

[37]  Gordon K. Smyth,et al.  A Modified Prony Algorithm for Exponential Function Fitting , 1995, SIAM J. Sci. Comput..