Curl-Based Finite Element Reconstruction of the Shear Modulus Without Assuming Local Homogeneity: Time Harmonic Case

In elasticity imaging, the shear modulus is obtained from measured tissue displacement data by solving an inverse problem based on the wave equation describing the tissue motion. In most inversion approaches, the wave equation is simplified using local homogeneity and incompressibility assumptions. This causes a loss of accuracy and therefore imaging artifacts in the resulting elasticity images. In this paper we present a new curl-based finite element method inversion technique that does not rely upon these simplifying assumptions. As done in previous research, we use the curl operator to eliminate the dilatational term in the wave equation, but we do not make the assumption of local homogeneity. We evaluate our approach using simulation data from a virtual tissue phantom assuming time harmonic motion and linear, isotropic, elastic behavior of the tissue. We show that our reconstruction results are superior to those obtained using previous curl-based methods with homogeneity assumption. We also show that with our approach, in the 2-D case, multi-frequency measurements provide better results than single-frequency measurements. Experimental results from magnetic resonance elastography of a CIRS elastography phantom confirm our simulation results and further demonstrate, in a quantitative and repeatable manner, that our method is accurate and robust.