Effects of frequency- and direction-dependent elastic materials on linearly elastic MRE image reconstructions

MR Elastography With Optimization-Based Phase Unwrapping and Traveling Wave Expansion-Based Neural Network (TWENN)

Magnetic Resonance Elastography (MRE) can characterize biomechanical properties of soft tissue for disease diagnosis and treatment planning. However, complicated wavefields acquired from MRE coupled with noise pose challenges for accurate displacement extraction and modulus estimation. Using optimization-based displacement extraction and Traveling Wave Expansion-based Neural Network (TWENN) modulus estimation, we propose a new pipeline for processing MRE images. An objective function with Dual Data Consistency (Dual-DC) has been used to ensure accurate phase unwrapping and displacement extraction. For the estimation of complex wavenumbers, a complex-valued neural network with displacement covariance as an input has been developed. A model of traveling wave expansion is used to generate training datasets for the network with varying levels of noise. The complex shear modulus map is obtained through fusion of multifrequency and multidirectional data. Validation using brain and liver simulation images demonstrates the practical value of the proposed pipeline, which can estimate the biomechanical properties with minimal root-mean-square errors when compared to state-of-the-art methods. Applications of the proposed method for processing MRE images of phantom, brain, and liver reveal clear anatomical features, robustness to noise, and good generalizability of the pipeline.

Stiffness reconstruction methods for MR elastography

Assessment of tissue stiffness is desirable for clinicians and researchers, as it is well established that pathophysiological mechanisms often alter the structural properties of tissue. Magnetic resonance elastography (MRE) provides an avenue for measuring tissue stiffness and has a long history of clinical application, including staging liver fibrosis and stratifying breast cancer malignancy. A vital component of MRE consists of the reconstruction algorithms used to derive stiffness from wave-motion images by solving inverse problems. A large range of reconstruction methods have been presented in the literature, with differing computational expense, required user input, underlying physical assumptions, and techniques for numerical evaluation. These differences, in turn, have led to varying accuracy, robustness, and ease of use. While most reconstruction techniques have been validated against in silico or in vitro phantoms, performance with real data is often more challenging, stressing the robustness and assumptions of these algorithms. This article reviews many current MRE reconstruction methods and discusses the aforementioned differences. The material assumptions underlying the methods are developed and various approaches for noise reduction, regularization, and numerical discretization are discussed. Reconstruction methods are categorized by inversion type, underlying assumptions, and their use in human and animal studies. Future directions, such as alternative material assumptions, are also discussed.

Gradient-Based Optimization for Poroelastic and Viscoelastic MR Elastography

We describe an efficient gradient computation for solving inverse problems arising in magnetic resonance elastography (MRE). The algorithm can be considered as a generalized 'adjoint method' based on a Lagrangian formulation. One requirement for the classic adjoint method is assurance of the self-adjoint property of the stiffness matrix in the elasticity problem. In this paper, we show this property is no longer a necessary condition in our algorithm, but the computational performance can be as efficient as the classic method, which involves only two forward solutions and is independent of the number of parameters to be estimated. The algorithm is developed and implemented in material property reconstructions using poroelastic and viscoelastic modeling. Various gradient- and Hessian-based optimization techniques have been tested on simulation, phantom and in vivo brain data. The numerical results show the feasibility and the efficiency of the proposed scheme for gradient calculation.

Observation of direction-dependent mechanical properties in the human brain with multi-excitation MR elastography

Magnetic resonance elastography (MRE) has shown promise in noninvasively capturing changes in mechanical properties of the human brain caused by neurodegenerative conditions. MRE involves vibrating the brain to generate shear waves, imaging those waves with MRI, and solving an inverse problem to determine mechanical properties. Despite the known anisotropic nature of brain tissue, the inverse problem in brain MRE is based on an isotropic mechanical model. In this study, distinct wave patterns are generated in the brain through the use of multiple excitation directions in order to characterize the potential impact of anisotropic tissue mechanics on isotropic inversion methods. Isotropic inversions of two unique displacement fields result in mechanical property maps that vary locally in areas of highly aligned white matter. Investigation of the corpus callosum, corona radiata, and superior longitudinal fasciculus, three highly ordered white matter tracts, revealed differences in estimated properties between excitations of up to 33%. Using diffusion tensor imaging to identify dominant fiber orientation of bundles, relationships between estimated isotropic properties and shear asymmetry are revealed. This study has implications for future isotropic and anisotropic MRE studies of white matter tracts in the human brain.

Quantitative 3D magnetic resonance elastography: Comparison with dynamic mechanical analysis

Purpose

Magnetic resonance elastography (MRE) is a rapidly growing noninvasive imaging technique for measuring tissue mechanical properties in vivo. Previous studies have compared two‐dimensional MRE measurements with material properties from dynamic mechanical analysis (DMA) devices that were limited in frequency range. Advanced DMA technology now allows broad frequency range testing, and three‐dimensional (3D) MRE is increasingly common. The purpose of this study was to compare 3D MRE stiffness measurements with those of DMA over a wide range of frequencies and shear stiffnesses.

Methods

3D MRE and DMA were performed on eight different polyvinyl chloride samples over 20–205 Hz with stiffness between 3 and 23 kPa. Driving frequencies were chosen to create 1.1, 2.2, 3.3, 4.4, 5.5, and 6.6 effective wavelengths across the diameter of the cylindrical phantoms. Wave images were analyzed using direct inversion and local frequency estimation algorithm with the curl operator and compared with DMA measurements at each corresponding frequency. Samples with sufficient spatial resolution and with an octahedral shear strain signal‐to‐noise ratio > 3 were compared.

Results

Consistency between the two techniques was measured with the intraclass correlation coefficient (ICC) and was excellent with an overall ICC of 0.99.

Conclusions

3D MRE and DMA showed excellent consistency over a wide range of frequencies and stiffnesses. Magn Reson Med 77:1184–1192, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution‐NonCommercial‐NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non‐commercial and no modifications or adaptations are made.

Magnetic resonance elastography (MRE) in cancer: Technique, analysis, and applications

Tissue mechanical properties are significantly altered with the development of cancer. Magnetic resonance elastography (MRE) is a noninvasive technique capable of quantifying tissue mechanical properties in vivo. This review describes the basic principles of MRE and introduces some of the many promising MRE methods that have been developed for the detection and characterization of cancer, evaluation of response to therapy, and investigation of the underlying mechanical mechanisms associated with malignancy.

Non-identifiability of the Rayleigh damping material model in magnetic resonance elastography

Magnetic Resonance Elastography (MRE) is an emerging imaging modality for quantifying soft tissue elasticity deduced from displacement measurements within the tissue obtained by phase sensitive Magnetic Resonance Imaging (MRI) techniques. MRE has potential to detect a range of pathologies, diseases and cancer formations, especially tumors. The mechanical model commonly used in MRE is linear viscoelasticity (VE). An alternative Rayleigh damping (RD) model for soft tissue attenuation is used with a subspace-based nonlinear inversion (SNLI) algorithm to reconstruct viscoelastic properties, energy attenuation mechanisms and concomitant damping behavior of the tissue-simulating phantoms. This research performs a thorough evaluation of the RD model in MRE focusing on unique identification of RD parameters, μI and ρI. Results show the non-identifiability of the RD model at a single input frequency based on a structural analysis with a series of supporting experimental phantom results. The estimated real shear modulus values (μR) were substantially correct in characterising various material types and correlated well with the expected stiffness contrast of the physical phantoms. However, estimated RD parameters displayed consistent poor reconstruction accuracy leading to unpredictable trends in parameter behaviour. To overcome this issue, two alternative approaches were developed: (1) simultaneous multi-frequency inversion; and (2) parametric-based reconstruction. Overall, the RD model estimates the real shear shear modulus (μR) well, but identifying damping parameters (μI and ρI) is not possible without an alternative approach.

Local mechanical properties of white matter structures in the human brain

The noninvasive measurement of the mechanical properties of brain tissue using magnetic resonance elastography (MRE) has emerged as a promising method for investigating neurological disorders. To date, brain MRE investigations have been limited to reporting global mechanical properties, though quantification of the stiffness of specific structures in the white matter architecture may be valuable in assessing the localized effects of disease. This paper reports the mechanical properties of the corpus callosum and corona radiata measured in healthy volunteers using MRE and atlas-based segmentation. Both structures were found to be significantly stiffer than overall white matter, with the corpus callosum exhibiting greater stiffness and less viscous damping than the corona radiata. Reliability of both local and global measures was assessed through repeated experiments, and the coefficient of variation for each measure was less than 10%. Mechanical properties within the corpus callosum and corona radiata demonstrated correlations with measures from diffusion tensor imaging pertaining to axonal microstructure.

Quantitative imaging methods for the development and validation of brain biomechanics models

Rapid deformation of brain tissue in response to head impact or acceleration can lead to numerous pathological changes, both immediate and delayed. Modeling and simulation hold promise for illuminating the mechanisms of traumatic brain injury (TBI) and for developing preventive devices and strategies. However, mathematical models have predictive value only if they satisfy two conditions. First, they must capture the biomechanics of the brain as both a material and a structure, including the mechanics of brain tissue and its interactions with the skull. Second, they must be validated by direct comparison with experimental data. Emerging imaging technologies and recent imaging studies provide important data for these purposes. This review describes these techniques and data, with an emphasis on magnetic resonance imaging approaches. In combination, these imaging tools promise to extend our understanding of brain biomechanics and improve our ability to study TBI in silico.

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

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.

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

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.

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

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.