Elastography is emerging as an imaging modality that can distinguish normal versus diseased tissues via their biomechanical properties. This paper reviews current approaches to elastography in three areas—quasi-static, harmonic and transient—and describes inversion schemes for each elastographic imaging approach. Approaches include first-order approximation methods; direct and iterative inversion schemes for linear elastic; isotropic materials and advanced reconstruction methods for recovering parameters that characterize complex mechanical behavior. The paper's objective is to document efforts to develop elastography within the framework of solving an inverse problem, so that elastography may provide reliable estimates of shear modulus and other mechanical parameters. We discuss issues that must be addressed if model-based elastography is to become the prevailing approach to quasi-static, harmonic and transient elastography: (1) developing practical techniques to transform the ill-posed problem with a well-posed one; (2) devising better forward models to capture the complex mechanical behavior of soft tissues and (3) developing better test procedures to evaluate the performance of modulus elastograms.

/pdf/model-based-elastography-a-survey-of-approaches-to-the-1vm8fucbux.pdf

Model-based elastography: a survey of approaches to the inverse elasticity problem

The technique of MR elastography (MRE) has emerged as a useful modality for quantitatively imaging the mechanical properties of soft tissues in vivo. Recently, MRE has been introduced as a clinical tool for evaluating chronic liver disease, but many other potential applications are being explored. These applications include measuring tissue changes associated with diseases of the liver, breast, brain, heart, and skeletal muscle including both focal lesions (e.g., hepatic, breast, and brain tumors) and diffuse diseases (e.g., fibrosis and multiple sclerosis). The purpose of this review article is to summarize some of the recent developments of MRE and to highlight some emerging applications.

https://onlinelibrary.wiley.com/doi/pdf/10.1002/jmri.23597

Review of MR elastography applications and recent developments.

A signal-to-noise ratio (SNR) measure based on the octahedral shear strain (the maximum shear strain in any plane for a 3D state of strain) is presented for magnetic resonance elastography (MRE), where motion-based SNR measures are commonly used. The shear strain, γ, is directly related to the shear modulus, μ, through the definition of shear stress, τ = μγ. Therefore, noise in the strain is the important factor in determining the quality of motion data, rather than the noise in the motion. Motion and strain SNR measures were found to be correlated for MRE of gelatin phantoms and the human breast. Analysis of the stiffness distributions of phantoms reconstructed from the measured motion data revealed a threshold for both strain and motion SNR where MRE stiffness estimates match independent mechanical testing. MRE of the feline brain showed significantly less correlation between the two SNR measures. The strain SNR measure had a threshold above which the reconstructed stiffness values were consistent between cases, whereas the motion SNR measure did not provide a useful threshold, primarily due to rigid body motion effects.

/pdf/an-octahedral-shear-strain-based-measure-of-snr-for-3d-mr-3hryid7xet.pdf

An octahedral shear strain-based measure of SNR for 3D MR elastography

Purpose: Nonlinear inversion (NLI) in MR elastography requires discretization of the displacement field for a finite element(FE) solution of the “forward problem”, and discretization of the unknown mechanical property field for the iterative solution of the “inverse problem”. The resolution requirements for these two discretizations are different: the forward problem requires sufficient resolution of the displacement FE mesh to ensure convergence, whereas lowering the mechanical property resolution in the inverse problem stabilizes the mechanical property estimates in the presence of measurement noise. Previous NLI implementations use the same FE mesh to support the displacement and property fields, requiring a trade-off between the competing resolution requirements. Methods: This work implements and evaluates multiresolution FE meshes for NLI elastography, allowing independent discretizations of the displacements and each mechanical property parameter to be estimated. The displacement resolution can then be selected to ensure mesh convergence, and the resolution of the property meshes can be independently manipulated to control the stability of the inversion. Results: Phantom experiments indicate that eight nodes per wavelength (NPW) are sufficient for accurate mechanical property recovery, whereas mechanical property estimation from 50 Hzin vivobrain data stabilizes once the displacement resolution reaches 1.7 mm (approximately 19 NPW). Viscoelasticmechanical property estimates of in vivobrain tissue show that subsampling the loss modulus while holding the storage modulus resolution constant does not substantially alter the storage modulusimages. Controlling the ratio of the number of measurements to unknown mechanical properties by subsampling the mechanical property distributions (relative to the data resolution) improves the repeatability of the property estimates, at a cost of modestly decreased spatial resolution. Conclusions: Multiresolution NLI elastography provides a more flexible framework for mechanical property estimation compared to previous single mesh implementations.

/pdf/multiresolution-mr-elastography-using-nonlinear-inversion-5gv1s1atky.pdf

Multiresolution MR elastography using nonlinear inversion

Magnetic resonance elastography (MRE) is used to quantify the viscoelastic shear modulus, G*, of human and animal tissues. Previously, values of G* determined by MRE have been compared to values from mechanical tests performed at lower frequencies. In this study, a novel dynamic shear test (DST) was used to measure G* of a tissue-mimicking material at higher frequencies for direct comparison to MRE. A closed-form solution, including inertial effects, was used to extract G* values from DST data obtained between 20 and 200 Hz. MRE was performed using cylindrical 'phantoms' of the same material in an overlapping frequency range of 100-400 Hz. Axial vibrations of a central rod caused radially propagating shear waves in the phantom. Displacement fields were fit to a viscoelastic form of Navier's equation using a total least-squares approach to obtain local estimates of G*. DST estimates of the storage G' (Re[G*]) and loss modulus G″ (Im[G*]) for the tissue-mimicking material increased with frequency from 0.86 to 0.97 kPa (20-200 Hz, n = 16), while MRE estimates of G' increased from 1.06 to 1.15 kPa (100-400 Hz, n = 6). The loss factor (Im[G*]/Re[G*]) also increased with frequency for both test methods: 0.06-0.14 (20-200 Hz, DST) and 0.11-0.23 (100-400 Hz, MRE). Close agreement between MRE and DST results at overlapping frequencies indicates that G* can be locally estimated with MRE over a wide frequency range. Low signal-to-noise ratio, long shear wavelengths and boundary effects were found to increase residual fitting error, reinforcing the use of an error metric to assess confidence in local parameter estimates obtained by MRE.

/pdf/viscoelastic-properties-of-soft-gels-comparison-of-magnetic-46w7rw6yvi.pdf

Viscoelastic properties of soft gels: comparison of magnetic resonance elastography and dynamic shear testing in the shear wave regime

Purpose: The clinical efficacy of breast elastography may be limited when the authors employ the assumption that soft tissues exhibit linear, frequency-independent isotropic mechanical properties during the recovery of shear modulus. Thus, the purpose of this research was to evaluate the degradation in performance incurred when linear-elastic MR reconstruction methods are applied to phantoms that are fabricated using viscoelastic materials.

Methods: To develop phantoms with frequency-dependent mechanical properties, the authors measured the complex modulus of two groups of cylindrical-shaped gelatin samples over a wide frequency range (up to 1 kHz) with the established principles of time-temperature superposition (TTS). In one group of samples, the authors added varying amounts of agar (1%–4%); in the other group, the authors added varying amounts of sucrose (2.5%–20%). To study how viscosity affected the performance of the linear-elastic reconstruction method, the authors constructed an elastically heterogeneous MR phantom to simulate the case where small viscoelastic lesions were surrounded by relatively nonviscous breast tissue. The breast phantom contained four linear, viscoelastic spherical inclusions (10 mm diameter) that were embedded in normal gelatin. The authors imaged the breast phantom with a clinical prototype of a MRE system and recovered the shear-modulus distribution using the overlapping-subzone-linear-elastic image-reconstruction method. The authors compared the recovered shear modulus to that measured using the TTS method.

Results: The authors demonstrated that viscoelastic phantoms could be fabricated by including sucrose in the gelation process and that small viscoelastic inclusions were visible in MR elastograms recovered using a linear-elastic MR reconstruction process; however, artifacts that degraded contrast and spatial resolution were more prominent in highly viscoelastic inclusions. The authors also established that the accuracy of the MR elastograms depended on the degree of viscosity that the inclusion exhibited.

Conclusions: The results demonstrated that reconstructing shear modulus from other constitutive laws, such as viscosity, should improve both the accuracy and quality of MR elastograms of the breast.

Keith M. Paulsen

Papers

The performance of steady-state harmonic magnetic resonance elastography when applied to viscoelastic materials