A pseudo-dynamic sub-optimal filter for elastography under static loading and measurements

We propose a pseudo-dynamic form of a sub-optimal Kalman filter for elastography of plane-strain models of soft tissues under strictly static deformations and partial measurements. Since the tissue material is nearly incompressible and is thus prone to volumetric locking via standard displacement-based finite element formulations, we use a Cosserat point approach for deriving the static equilibrium equations. A pseudo-dynamical form of the equilibrium equations, with added noise and appropriate augmentation by the discretized shear modulus as additional states, is then adopted as the process equation such that its steady-state solution approaches the static response of the plane-strain model. A fictitious noise of small intensity is also added to the measurement equation and, following linearization of the process equation, a Kalman filter is applied to reconstruct the shear modulus profile. We present several numerical experiments, some of which also bring forth the relative advantages of the proposed approach over a deterministic reconstruction based on a quasi-Newton search.

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

[2]  Bernard Hanzon,et al.  A differential geometric approach to nonlinear filtering: the projection filter , 1998, IEEE Trans. Autom. Control..

[3]  A.R. Skovoroda,et al.  Tissue elasticity reconstruction based on ultrasonic displacement and strain images , 1995, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[4]  J. C. Simo,et al.  Variational and projection methods for the volume constraint in finite deformation elasto-plasticity , 1985 .

[5]  T. Hughes Generalization of selective integration procedures to anisotropic and nonlinear media , 1980 .

[6]  Paul E. Barbone,et al.  Elastic modulus imaging: On the uniqueness and nonuniqueness of the elastography inverse problem in , 2004 .

[7]  Peter Wriggers,et al.  Response of a nonlinear elastic general Cosserat brick element in simulations typically exhibiting locking and hourglassing , 2005 .

[8]  Marc Bonnet,et al.  Inverse problems in elasticity , 2005 .

[9]  A. Corigliano,et al.  Parameter identification in explicit structural dynamics: performance of the extended Kalman filter , 2004 .

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

[11]  Y. Saad,et al.  Krylov Subspace Methods on Supercomputers , 1989 .

[12]  Roger B. Sidje,et al.  Expokit: a software package for computing matrix exponentials , 1998, TOMS.

[13]  A. Maniatty,et al.  Shear modulus reconstruction in dynamic elastography: time harmonic case , 2006, Physics in medicine and biology.

[14]  Faouzi Kallel,et al.  Tissue elasticity reconstruction using linear perturbation method , 1996, IEEE Trans. Medical Imaging.

[15]  M. Rubin,et al.  A new 3-D finite element for nonlinear elasticity using the theory of a Cosserat point , 2003 .

[16]  C. Sumi,et al.  Estimation of shear modulus distribution in soft tissue from strain distribution , 1995, IEEE Transactions on Biomedical Engineering.

[17]  Simon J. Godsill,et al.  On sequential Monte Carlo sampling methods for Bayesian filtering , 2000, Stat. Comput..

[18]  J. Bamber,et al.  Quantitative elasticity imaging: what can and cannot be inferred from strain images. , 2002, Physics in medicine and biology.

[19]  Yi Liu,et al.  Tomography-based 3-D anisotropic elastography using boundary measurements , 2005, IEEE Transactions on Medical Imaging.

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

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

[22]  C. S. Manohar,et al.  A sequential importance sampling filter with a new proposal distribution for state and parameter estimation of nonlinear dynamical systems , 2008, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[23]  J. C. Simo,et al.  A CLASS OF MIXED ASSUMED STRAIN METHODS AND THE METHOD OF INCOMPATIBLE MODES , 1990 .

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

[25]  M. Doyley,et al.  Evaluation of an iterative reconstruction method for quantitative elastography , 2000 .

[26]  G. Milstein Numerical Integration of Stochastic Differential Equations , 1994 .

[27]  Jonathan Ophir,et al.  Enhancing the performance of model-based elastography by incorporating additional a priori information in the modulus image reconstruction process , 2006, Physics in medicine and biology.

[28]  K. R. Raghavan,et al.  Forward and inverse problems in elasticity imaging of soft tissues , 1994 .

[29]  S. Arridge,et al.  State-estimation approach to the nonstationary optical tomography problem. , 2003, Journal of the Optical Society of America. A, Optics, image science, and vision.

[30]  Peter Wriggers,et al.  A new finite element based on the theory of a Cosserat point—extension to initially distorted elements for 2D plane strain , 2007 .

[31]  A Galerkin/least-square finite element formulation for nearly incompressible elasticity/stokes flow , 2007 .

[32]  G. Kallianpur Stochastic Filtering Theory , 1980 .