Gradient-Based Optimization for Poroelastic and Viscoelastic MR Elastography

We describe an efficient gradient computation for solving inverse problems arising in magnetic resonance elastography (MRE). The algorithm can be considered as a generalized 'adjoint method' based on a Lagrangian formulation. One requirement for the classic adjoint method is assurance of the self-adjoint property of the stiffness matrix in the elasticity problem. In this paper, we show this property is no longer a necessary condition in our algorithm, but the computational performance can be as efficient as the classic method, which involves only two forward solutions and is independent of the number of parameters to be estimated. The algorithm is developed and implemented in material property reconstructions using poroelastic and viscoelastic modeling. Various gradient- and Hessian-based optimization techniques have been tested on simulation, phantom and in vivo brain data. The numerical results show the feasibility and the efficiency of the proposed scheme for gradient calculation.

[1]  Jorge Nocedal,et al.  Global Convergence Properties of Conjugate Gradient Methods for Optimization , 1992, SIAM J. Optim..

[2]  K D Paulsen,et al.  An overlapping subzone technique for MR‐based elastic property reconstruction , 1999, Magnetic resonance in medicine.

[3]  J. Shirron,et al.  On the noninvasive determination of material parameters from a knowledge of elastic displacements theory and numerical simulation , 1998, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[4]  Peter M. Pinsky,et al.  Recovery of shear modulus in elastography using an adjoint method with B-spline representation , 2005 .

[5]  Jeffrey C Bamber,et al.  Evaluation of the adjoint equation based algorithm for elasticity imaging , 2004, Physics in medicine and biology.

[6]  C. Vogel.  Computational Methods for Inverse Problems , 1987 .

[7]  K D Paulsen,et al.  Elasticity reconstruction from experimental MR displacement data: initial experience with an overlapping subzone finite element inversion process. , 2000, Medical physics.

[8]  C. Jack,et al.  Decreased brain stiffness in Alzheimer's disease determined by magnetic resonance elastography , 2011, Journal of magnetic resonance imaging : JMRI.

[9]  Caroline Maleke,et al.  Quantitative viscoelastic parameters measured by harmonic motion imaging , 2009, Physics in medicine and biology.

[10]  P. Wolfe.  Convergence Conditions for Ascent Methods. II , 1969 .

[11]  R L Ehman,et al.  Complex‐valued stiffness reconstruction for magnetic resonance elastography by algebraic inversion of the differential equation , 2001, Magnetic resonance in medicine.

[12]  Assad A. Oberai,et al.  INVERSE PROBLEMS PII: S0266-5611(03)54272-1 Solution of inverse problems in elasticity imaging using the adjoint method , 2003 .

[13]  K D Paulsen,et al.  Effects of frequency- and direction-dependent elastic materials on linearly elastic MRE image reconstructions , 2010, Physics in medicine and biology.

[14]  A. Manduca,et al.  Magnetic resonance elastography by direct visualization of propagating acoustic strain waves. , 1995, Science.

[15]  Keith D. Paulsen,et al.  Spatially-Resolved Hydraulic Conductivity Estimation Via Poroelastic Magnetic Resonance Elastography , 2014, IEEE Transactions on Medical Imaging.

[16]  Mickael Tanter,et al.  MR elastography of breast lesions: Understanding the solid/liquid duality can improve the specificity of contrast‐enhanced MR mammography , 2007, Magnetic resonance in medicine.

[17]  M. McGarry,et al.  Rayleigh Damped Magnetic Resonance Elastograpy , 2008 .

[18]  Mickael Tanter,et al.  Viscoelastic shear properties of in vivo breast lesions measured by MR elastography. , 2005, Magnetic resonance imaging.

[19]  K D Paulsen,et al.  Magnetic resonance elastography using 3D gradient echo measurements of steady-state motion. , 2001, Medical physics.

[20]  James G Chase,et al.  Silicone breast phantoms for elastographic imaging evaluation. , 2013, Medical physics.

[21]  Keith D. Paulsen,et al.  MR elastography of hydrocephalus , 2009, Medical Imaging.

[22]  J. Nocedal,et al.  Global Convergence of a Class of Quasi-newton Methods on Convex Problems, Siam Some Global Convergence Properties of a Variable Metric Algorithm for Minimization without Exact Line Searches, Nonlinear Programming, Edited , 1996 .

[23]  J. Shirron,et al.  Evaluation of a material parameter extraction algorithm using MRI-based displacement measurements , 2000, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[24]  Dieter Klatt,et al.  Assessment of liver viscoelasticity using multifrequency MR elastography , 2008, Magnetic resonance in medicine.

[25]  C. Pichot,et al.  Inverse scattering: an iterative numerical method for electromagnetic imaging , 1991 .

[26]  D. Plewes,et al.  Elastic moduli of normal and pathological human breast tissues: an inversion-technique-based investigation of 169 samples , 2007, Physics in medicine and biology.

[27]  J. Ophir,et al.  Elastography: A Quantitative Method for Imaging the Elasticity of Biological Tissues , 1991, Ultrasonic imaging.

[28]  Ralph Sinkus,et al.  Magnetic resonance elastography for the noninvasive staging of liver fibrosis. , 2008, Gastroenterology.

[29]  Keith D. Paulsen,et al.  Estimating Hydraulic Conductivity in Vivo Using Magnetic Resonance Elastography , 2011 .

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

[31]  Keith D. Paulsen,et al.  Modeling of Soft Poroelastic Tissue in Time-Harmonic MR Elastography , 2009, IEEE Transactions on Biomedical Engineering.

[32]  Huan Wang,et al.  Local mechanical properties of white matter structures in the human brain , 2013, NeuroImage.

[33]  Kim Thomson,et al.  Quantitative shear wave ultrasound elastography: initial experience in solid breast masses , 2010, Breast Cancer Research.

[34]  Keith D. Paulsen,et al.  Including Spatial Information in Nonlinear Inversion MR Elastography Using Soft Prior Regularization , 2013, IEEE Transactions on Medical Imaging.

[35]  M. Nash,et al.  Electromechanical model of excitable tissue to study reentrant cardiac arrhythmias. , 2004, Progress in biophysics and molecular biology.

[36]  J. Georgiadis,et al.  Multiresolution MR elastography using nonlinear inversion. , 2012, Medical physics.

[37]  Assad A. Oberai,et al.  Solution of the nonlinear elasticity imaging inverse problem: The incompressible case. , 2011 .

[38]  Brett E. Bouma,et al.  Tissue Elasticity Estimation with Optical Coherence Elastography: Toward Mechanical Characterization of In Vivo Soft Tissue , 2005, Annals of Biomedical Engineering.

[39]  John B Weaver,et al.  Contrast detection in fluid-saturated media with magnetic resonance poroelastography. , 2010, Medical physics.

[40]  K D Paulsen,et al.  Three‐dimensional subzone‐based reconstruction algorithm for MR elastography , 2001, Magnetic resonance in medicine.

[41]  Dieter Klatt,et al.  In vivo viscoelastic properties of the brain in normal pressure hydrocephalus , 2010, NMR in biomedicine.

[42]  M. Henriksen,et al.  Nonlinear viscoelastic stress analysis—a finite element approach , 1984 .

[43]  K D Paulsen,et al.  Subzone based magnetic resonance elastography using a Rayleigh damped material model. , 2011, Medical Physics (Lancaster).

[44]  Keith D. Paulsen,et al.  Magnetic Resonance Poroelastography: An Algorithm for Estimating the Mechanical Properties of Fluid-Saturated Soft Tissues , 2010, IEEE Transactions on Medical Imaging.

[45]  Carl-Fredrik Westin,et al.  Local multiscale frequency and bandwidth estimation , 1994, Proceedings of 1st International Conference on Image Processing.

[46]  John B Weaver,et al.  Brain mechanical property measurement using MRE with intrinsic activation , 2012, Physics in medicine and biology.

[47]  Frauke Zipp,et al.  MR-elastography reveals degradation of tissue integrity in multiple sclerosis , 2010, NeuroImage.

[48]  Patrick Amestoy,et al.  Hybrid scheduling for the parallel solution of linear systems , 2006, Parallel Comput..

[49]  R. Willinger,et al.  Shear Properties of Brain Tissue over a Frequency Range Relevant for Automotive Impact Situations: New Experimental Results. , 2004, Stapp car crash journal.

[50]  S E Salcudean,et al.  Sparsity regularization in dynamic elastography , 2012, Physics in medicine and biology.

[51]  Jorge Nocedal,et al.  Representations of quasi-Newton matrices and their use in limited memory methods , 1994, Math. Program..

[52]  Sevan Goenezen,et al.  Inverse Problems , 2008 .

[53]  R. Sinkus,et al.  High-resolution tensor MR elastography for breast tumour detection. , 2000, Physics in medicine and biology.

[54]  J. F. Greenleaf,et al.  Magnetic resonance elastography: Non-invasive mapping of tissue elasticity , 2001, Medical Image Anal..

[55]  R. Ehman,et al.  Early detection of nonalcoholic steatohepatitis in patients with nonalcoholic fatty liver disease by using MR elastography. , 2011, Radiology.

[56]  P. Janmey,et al.  Tissue Cells Feel and Respond to the Stiffness of Their Substrate , 2005, Science.

[57]  Steven E. Benzley,et al.  A Comparison of All Hexagonal and All Tetrahedral Finite Element Meshes for Elastic and Elasto-plastic Analysis , 2011 .

[58]  Assad A. Oberai,et al.  Adjoint‐weighted variational formulation for the direct solution of inverse problems of general linear elasticity with full interior data , 2010 .

[59]  Jorge Nocedal,et al.  Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization , 1997, TOMS.

[60]  Eli E. W. Van Houten,et al.  Use of a Rayleigh damping model in elastography , 2008, Medical & Biological Engineering & Computing.

[61]  Patrick Amestoy,et al.  A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling , 2001, SIAM J. Matrix Anal. Appl..

[62]  Curtis L. Johnson,et al.  Magnetic resonance elastography of the brain using multishot spiral readouts with self‐navigated motion correction , 2013, Magnetic resonance in medicine.

[63]  Keith D. Paulsen,et al.  Development of a poroelastic dynamic mechanical analysis technique for biphasic media , 2013, Medical Imaging.

[64]  E. Houten,et al.  Parameter identification in a generalized time-harmonic Rayleigh damping model for elastography. , 2014 .

[65]  John B Weaver,et al.  Initial in vivo experience with steady‐state subzone‐based MR elastography of the human breast , 2003, Journal of magnetic resonance imaging : JMRI.

[66]  Jochen Guck,et al.  Critical review: cellular mechanobiology and amoeboid migration. , 2010, Integrative biology : quantitative biosciences from nano to macro.

[67]  J. Shewchuk.  An Introduction to the Conjugate Gradient Method Without the Agonizing Pain , 1994 .

[68]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[69]  J. Georgiadis,et al.  Suitability of poroelastic and viscoelastic mechanical models for high and low frequency MR elastography. , 2015, Medical physics.

[70]  Matthew D.j. McGarry.  Improvement and evaluation of nonlinear inversion MR elastography , 2013 .