Magnetic resonance elastography (MRE) is a powerful technique for noninvasive determination of the biomechanical properties of tissue, with important applications in disease diagnosis. A typical experimental scenario is to induce waves in the tissue by time-harmonic external mechanical oscillation and then measure the tissue's displacement at fixed spatial positions 8 times during a complete time-period, extracting the dominant frequency signal from the discrete Fourier transform in time. Accurate reconstruction of the tissue's elastic moduli from MRE data is a challenging inverse problem, and we derive and analyze two new methods which address different aspects. The first of these concerns the time signal: using only the dominant frequency component loses information for noisy data and typically gives a complex value for the (real) shear modulus, which is then hard to interpret. Our new reconstruction method is based on the Fourier time-interpolant of the displacement: it uses all the measured information and automatically gives a real value of shear modulus up to rounding error. This derivation is for homogeneous materials, and our second new method (stacked frequency wave inversion, SFWI) concerns the inhomogeneous shear modulus in the time-harmonic case. The underlying problem is ill-conditioned because the coefficient of the shear modulus in the governing equations can be zero or small, and the SFWI approach overcomes this by combining approximations at different frequencies into a single overdetermined matrix--vector equation. Careful numerical tests confirm that both these new algorithms perform well.
- magnetic resonance elastography
- inverse problem