Estimation of material parameters from slow and fast shear waves in an incompressible, transversely isotropic material

• Published in: Journal of biomechanics Vol. 48; no. 15; pp. 4002 - 4009
• Main Authors: Tweten, Dennis J., Okamoto, Ruth J., Schmidt, John L., Garbow, Joel R., Bayly, Philip V.
• Format: Journal Article
• Language: English
• Published: United States Elsevier Ltd 26.11.2015; Elsevier Limited
• Subjects: Accuracy; Anisotropy; Elasticity Imaging Techniques; Estimates; Filtering; Finite Element Analysis; Finite element method; Humans; Inversion algorithms; Mathematical analysis; Mathematical models; Models, Theoretical; MR elastography; Physical Medicine and Rehabilitation; Reproducibility of Results; Shear; Shear waves; Sound waves; Studies; Transversely isotropic material; Viscoelasticity; Wave propagation; Shear waves; MR elastography; Inversion algorithms; Anisotropy; Transversely isotropic material
• ISSN: 0021-9290; 1873-2380; 1873-2380
• DOI: 10.1016/j.jbiomech.2015.09.009

This paper describes a method to estimate mechanical properties of soft, anisotropic materials from measurements of shear waves with specific polarization and propagation directions. This method is applicable to data from magnetic resonance elastography (MRE), which is a method for measuring shear waves in live subjects or in vitro samples. Here, we simulate MRE data using finite element analysis. A nearly incompressible, transversely isotropic (ITI) material model with three parameters (shear modulus, shear anisotropy, and tensile anisotropy) is used, which is appropriate for many fibrous, biological tissues. Both slow and fast shear waves travel concurrently through such a material with speeds that depend on the propagation direction relative to fiber orientation. A three-parameter estimation approach based on directional filtering and isolation of slow and fast shear wave components (directional filter inversion, or DFI) is introduced. Wave speeds of each isolated shear wave component are estimated using local frequency estimation (LFE), and material properties are calculated using weighted least squares. Data from multiple finite element simulations are used to assess the accuracy and reliability of DFI for estimation of anisotropic material parameters.