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

Abstract

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.

Evaluation of wave delivery methodology for brain MRE: Insights from computational simulations

“…In our previous brain MRE simulation work it was found that inversion artifacts occurred close to interfaces between brain tissue and the FTM and ventricles. The conclusion of that earlier study was that errors at the boundaries were caused by a combination of factors: 1) reflection, refraction, and scattering at tissue boundaries leading to wave interference, which results in inversion artifacts at larger sampling steps (3 mm); 2) material heterogeneity bringing about errors in the direct inversion algorithm, which assumes local homogeneity ; and 3) averaging across the tissue boundaries due to interpolation, derivative calculations, and smoothing of the curl vector field during the inversion . However, the findings of this present study emphasize the importance of wave reflection and the resulting interference patterns, as the different wave delivery methods produce different predominant directions of motion accompanied by different magnitudes of inversion error.…”

Focal application of low-dose-rate brachytherapy for prostate cancer: a pilot study

“…The current protocol includes the collection of MR elastography and mpTRUS data. Magnetic resonance elastography is performed as a part of the mpMRI procedure [51,52] and mpTRUS, which also includes TRUS elastography that is collected immediately prior to TTMB [53,54,55]. A sufficiently accurate image-based cancer detection method (e.g.…”

MR elastography measurement of the effect of passive warmup prior to eccentric exercise on thigh muscle mechanical properties

“…The CIRS 049 phantom had been previously used to verify MRE measurements 55, 56. In both cases different inversion approaches to the ones used here (i.e., finite element modeling) were applied; however, the recovered stiffness values were found to be lower than values reported by the manufacturer for the stiffest target and higher for the softest target.…”