A Comparison of Finite Element-Based Inversion Algorithms, Local Frequency Estimation, and Direct Inversion Approach Used in MRE

Comparison of ultrasound elastography, magnetic resonance elastography and finite element model to quantify nonlinear shear modulus

Objective. The aim of this study is to validate the estimation of the nonlinear shear modulus (A) from the acoustoelasticity theory with two experimental methods, ultrasound (US) elastography and magnetic resonance elastography (MRE), and a finite element method.Approach. Experiments were performed on agar (2%)-gelatin (8%) phantom considered as homogeneous, elastic and isotropic. Two specific setups were built to ensure a uniaxial stress step by step on the phantom, one for US and a nonmagnetic version for MRE. The stress was controlled identically in both imaging techniques, with a water tank placed on the top of the phantom and filled with increasing masses of water during the experiment. In US, the supersonic shear wave elastography was implemented on an ultrafast US device, driving a 6 MHz linear array to measure shear wave speed. In MRE, a gradient-echo sequence was used in which the three spatial directions of a 40 Hz continuous wave displacement generated with an external driver were encoded successively. Numerically, a finite element method was developed to simulate the propagation of the shear wave in a uniaxially stressed soft medium.Main results. Similar shear moduli were estimated at zero stress using experimental methods,μ0US= 12.3 ± 0.3 kPa andμ0MRE= 11.5 ± 0.7 kPa. Numerical simulations were set with a shear modulus of 12 kPa and the resulting nonlinear shear modulus was found to be -58.1 ± 0.7 kPa. A very good agreement between the finite element model and the experimental models (AUS= -58.9 ± 9.9 kPa andAMRE= -52.8 ± 6.5 kPa) was obtained.Significance. These results show the validity of such nonlinear shear modulus measurement quantification in shear wave elastography. This work paves the way to develop nonlinear elastography technique to get a new biomarker for medical diagnosis.

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.

Phase-regularized and displacement-regularized compressed sensing for fast magnetic resonance elastography

Liver magnetic resonance elastography (MRE) is a noninvasive stiffness measurement technique that captures the tissue displacement in the phase of the signal. To limit the scanning time to a single breath-hold, liver MRE usually involves advanced readout techniques such as simultaneous multislice (SMS) or multishot methods. Furthermore, all these readout techniques require additional in-plane acceleration using either parallel imaging capabilities, such as sensitivity encoding (SENSE), or k -space undersampling, such as compressed sensing (CS). However, these methods apply a single regularization function on the complex image. This study aims to design and evaluate methods that use separate regularization on the magnitude and phase of MRE to exploit their distinct spatiotemporal characteristics. Specifically, we introduce two compressed sensing methods. The first method, termed phase-regularized compressed sensing (PRCS), applies a two-dimensional total variation (TV) prior to the magnitude and two-dimensional wavelet regularization to the phase. The second method, termed displacement-regularized compressed sensing (DRCS), exploits the spatiotemporal redundancy using 3D total variation on the magnitude. Additionally, DRCS includes a displacement fitting function to apply wavelet regularization to the displacement phasor. Both DRCS and PRCS were evaluated with different levels of compression factors in three datasets: an in silico abdomen dataset, an in vitro tissue-mimicking phantom, and an in vivo liver dataset. The reconstructed images were compared with the full sampled reconstruction, zero-filling reconstruction, wavelet-regularized compressed sensing, and a low rank plus sparse reconstruction. The metrics used for quantitative evaluation were the structural similarity index (SSIM) of magnitude (M-SSIM), displacement (D-SSIM), and shear modulus (S-SSIM), and mean shear modulus. Results from highly undersampled in silico and in vitro datasets demonstrate that the DRCS method provides higher reconstruction quality than the conventional compressed sensing method for a wide range of stiffness values. Notably, DRCS provides 24% and 22% increase in D-SSIM compared with CS for the in silico and in vitro datasets, respectively. Comparison with liver stiffness measured from full sampled data and highly undersampled data (CR=4) demonstrates that the DRCS method provided the strongest correlation ( R 2 =0.95), second-lowest mean bias (-0.18 kPa, lowest for CS with -0.16 kPa), and lowest coefficient of variation (CV=3.6%). Our results demonstrate the potential of using DRCS to improve the reconstruction quality of accelerated MRE.

Fast acquisition of propagating waves in humans with low-field MRI: Toward accessible MR elastography

Most commonly used at clinical magnetic fields (1.5 to 3 T), magnetic resonance elastography (MRE) captures mechanical wave propagation to reconstruct the mechanical properties of soft tissue with MRI. However, in terms of noninvasively assessing disease progression in a broad range of organs (e.g., liver, breast, skeletal muscle, and brain), its accessibility is limited and its robustness is challenged when magnetic susceptibility differences are encountered. Low-field MRE offers an opportunity to overcome these issues, and yet it has never been demonstrated in vivo in humans with magnetic fields <1.5 T mainly because of the long acquisition times required to achieve a sufficient signal-to-noise ratio. Here, we describe a method to accelerate 3D motion-sensitized MR scans at 0.1 T using only 10% k-space sampling combined with a high-performance detector and an efficient encoding acquisition strategy. Its application is demonstrated in vivo in the human forearm for a single motion-encoding direction in less than 1 min.

Numerical simulation of wave propagation through interfaces using the extended finite element method for magnetic resonance elastography

Magnetic resonance elastography (MRE) is an elasticity imaging technique for quantitatively assessing the stiffness of human tissues. In MRE, finite element method (FEM) is widely used for modeling wave propagation and stiffness reconstruction. However, in front of inclusions with complex interfaces, FEM can become burdensome in terms of the model partition and computationally expensive. In this work, we implement a formulation of FEM, known as the eXtended finite element method (XFEM), which is a method used for modeling discontinuity like crack and heterogeneity. Using a level-set method, it makes the interface independent of the mesh, thus relieving the meshing efforts. We investigate this method in two studies: wave propagation across an oblique linear interface and stiffness reconstruction of a random-shape inclusion. In the first study, numerical results by XFEM and FEM models revealing the wave conversion rules at linear interface are presented and successfully compared to the theoretical predictions. The second study, investigated in a pseudo-practical application, demonstrates further the applicability of XFEM in MRE and the convenience, accuracy, and speed of XFEM with respect to FEM. XFEM can be regarded as a promising alternative to FEM for inclusion modeling in MRE.

Design, Construction, and Implementation of a Magnetic Resonance Elastography Actuator for Research Purposes

Magnetic resonance elastography (MRE) is a technique for determining the mechanical response of soft materials using applied harmonic deformation of the material and a motion-sensitive magnetic resonance imaging sequence. This technique can elucidate significant information about the health and development of human tissue such as liver and brain, and has been used on phantom models (e.g., agar, silicone) to determine their suitability for use as a mechanical surrogate for human tissues in experimental models. The applied harmonic deformation used in MRE is generated by an actuator, transmitted in bursts of a specified duration, and synchronized with the magnetic resonance signal excitation. These actuators are most often a pneumatic design (common for human tissues or phantoms) or a piezoelectric design (common for small animal tissues or phantoms). Here, we describe how to design and assemble both a pneumatic and a piezoelectric MRE actuator for research purposes. For each of these actuator types, we discuss displacement requirements, end-effector options and challenges, electronics and electronic-driving requirements and considerations, and full MRE implementation. We also discuss how to choose the actuator type, size, and power based on the intended material and use. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Design, construction, and implementation of a convertible pneumatic MRE actuator for use with tissues and phantom models Basic Protocol 2: Design, construction, and implementation of a piezoelectric MRE actuator for localized excitation in phantom models.

In silicoevaluation and optimisation of magnetic resonance elastography of the liver

Objective.Magnetic resonance elastography (MRE) is widely adopted as a biomarker of liver fibrosis. However,in vivoMRE accuracy is difficult to assess.Approach.Finite element model (FEM) simulation was employed to evaluate liver MRE accuracy and inform methodological optimisation. MRE data was simulated in a 3D FEM of the human torso including the liver, and compared with spin-echo echo-planar imaging MRE acquisitions. The simulated MRE results were compared with the ground truth magnitude of the complex shear modulus (∣G*∣) for varying: (1) ground truth liver ∣G*∣; (2) simulated imaging resolution; (3) added noise; (4) data smoothing. Motion and strain-based signal-to-noise (SNR) metrics were evaluated on the simulated data as a means to select higher-quality voxels for preparation of acquired MRE summary statistics of ∣G*∣.Main results.The simulated MRE accuracy for a given ground truth ∣G*∣ was found to be a function of imaging resolution, motion-SNR and smoothing. At typical imaging resolutions, it was found that due to under-sampling of the MRE wave-field, combined with motion-related noise, the reconstructed simulated ∣G*∣ could contain errors on the scale of the difference between liver fibrosis stages, e.g. 54% error for ground truth ∣G*∣ = 1 kPa. Optimum imaging resolutions were identified for given ground truth ∣G*∣ and motion-SNR levels.Significance.This study provides important knowledge on the accuracy and optimisation of liver MRE. For example, for motion-SNR ≤ 5, to distinguish between liver ∣G*∣ of 2 and 3 kPa (i.e. early-stage liver fibrosis) it was predicted that the optimum isotropic voxel size is 4-6 mm.

Practical and clinical applications of pancreatic magnetic resonance elastography: a systematic review

Magnetic resonance elastography (MRE) is a non-invasive technique suitable for assessing mechanical properties of tissues, i.e., stiffness. MRE of the pancreas is relatively new, but recently an increasing number of studies have successfully assessed pancreas diseases with MRE aiming to differentiate healthy from pathological pancreatic tissue with or without fibrosis. This review will systematically describe the practical and clinical applications of pancreatic MRE. We conducted a systematic literature search with a pre-specified search strategy using PubMed and Embase according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. English peer-reviewed articles applying MRE of the pancreas were included. Two independent reviewers assessed the studies. The literature search yielded 14 studies. The pancreatic stiffness for healthy volunteers ranged from 1.11. to 1.21 kPa at a driver frequency of 40 Hz. In benign tumors, the stiffness values were slightly higher or sometimes even lower (range 0.78 to 2.00 kPa), compared to the healthy pancreas parenchyma whereas, in malignant tumors, the stiffness values tended to be higher (1.42 to 6.06 kPa). The pancreatic stiffness was increased in both acute (median: 1.99 kPa) and chronic pancreatitis (> 1.50 kPa). MRE is a promising technique for detecting and quantifying pancreatic stiffness. It is related to fibrosis and seems to be useful in assessing treatment response and clinical follow-up of pancreatic diseases. However, most of the described practical settings were characterized by a lack of uniformity and inconsistency in reporting standards across studies. Harmonization between centers is necessary to achieve more consensus and optimization of pancreatic MRE protocols.

Plane wave elastography: a frequency-domain ultrasound shear wave elastography approach

In this paper, we propose plane wave elastography (PWE), a novel ultrasound shear wave elastography (SWE) approach. Currently, commercial methods for SWE rely on directional filtering based on the prior knowledge of the wave propagation direction, to remove complicated wave patterns formed due to reflection and refraction. The result is a set of decomposed directional waves that are separately analyzed to construct shear modulus fields that are then combined through compounding. Instead, PWE relies on a rigorous representation of the wave propagation using the frequency-domain scalar wave equation to automatically select appropriate propagation directions and simultaneously reconstruct shear modulus fields. Specifically, assuming a homogeneous, isotropic, incompressible, linear-elastic medium, we represent the solution of the wave equation using a linear combination of plane waves propagating in arbitrary directions. Given this closed-form solution, we formulate the SWE problem as a nonlinear least-squares optimization problem which can be solved very efficiently. Through numerous phantom studies, we show that PWE can handle complicated waveforms without prior filtering and is competitive with state-of-the-art that requires prior filtering based on the knowledge of propagation directions.

Three-Dimensional Multi-Frequency Shear Wave Absolute Vibro-Elastography (3D S-WAVE) With a Matrix Array Transducer: Implementation and Preliminary In Vivo Study of the Liver

Magnetic resonance elastography (MRE) is commonly regarded as the imaging-based gold-standard for liver fibrosis staging, comparable to biopsy. While ultrasound-based elastography methods for liver fibrosis staging have been developed, they are confined to a 1D or a 2D region of interest and to a limited depth. 3D Shear Wave Absolute Vibro-Elastography (S-WAVE) is a steady-state, external excitation, volumetric elastography technique that is similar to MRE, but has the additional advantage of multi-frequency excitation. We present a novel ultrasound matrix array implementation of S-WAVE that takes advantage of 3D imaging. We use a matrix array transducer to sample axial multi-frequency steady-state tissue motion over a volume, using a Color Power Angiography sequence. Tissue motion with the frequency components {40,50,60} and {45,55,65} Hz are acquired over a (90° lateral) × (40° elevational) × (16 cm depth) sector with an acquisition time of 12 seconds. We compute the elasticity map in 3D using local spatial frequency estimation. We characterize this new approach in tissue phantoms against measurements obtained with transient elastography and MRE. Six healthy volunteers and eight patients with chronic liver disease were imaged. Their MRE and S-WAVE volumes were aligned using T1 to B-mode registration for direct comparison in common regions of interest. S-WAVE and MRE results are correlated with R² = 0.92, while MRE and TE results are correlated with R² = 0.71. Our findings show that S-WAVE with matrix array has the potential to deliver a similar assessment of liver fibrosis as MRE in a more accessible, inexpensive way, to a broader set of patients.

Enhanced complex local frequency elastography method for tumor viscoelastic shear modulus reconstruction

BACKGROUND AND OBJECTIVES

The Mayo Clinic provides a magnetic resonance (MR) elastography software named MRE Wave, which uses the conventional local frequency elastography (LFE) method. However, MRE Wave is unable to supply complex viscoelasticity maps for elastography. We sought to improve the local frequency estimation algorithm used in LFE, which we refer to as the Enhanced Complex Local Frequency Elastography (EC-LFE) algorithm.

METHODS

The proposed algorithm uses wave equations under the hypotheses of being linear, isotropic, and locally homogeneous. Two 2D simulation models were used to investigate the accuracy and sensitivity of the EC-LFE algorithm for detecting small tumors. The corresponding statistical parameters were the relative root mean square (RMS) error and contrast-to-noise ratio (CNR). EC-LFE was investigated with two different parameter sets, one with an optimally chosen parameter ξ (EC-LFE Adj, for short) and the other with ξ = 0 (EC-LFE0). We compared the MRE Wave and the EC-LFE using series signal-to-noise (SNR) wave data.

RESULTS

The elasticity RMS error of the MRE Wave software was about 1%, and that of the EC-LFE0 and EC-LFE Adj were about 0.2%. The elasticity standard deviation of the MRE Wave software was about 3% of the mean value, and those of the EC-LFE0 and EC-LFE Adj were about 1% of the mean value. The elasticity CNR value of EC-LFE0 reached 1.93 times that of the MRE Wave in the region of small tumors (less than 10-point sampling). The viscosity RMS errors of the EC-LFE0 could be less than 5%.

CONCLUSION

Compared to conventional methods, the EC-LFE was more accurate and sensitive for small tumor detection and exhibited higher noise immunity. The improved algorithm output more parameters and outperformed than the MRE Wave, thereby rendering them more suitable for clinical applications.

Requirements for accurate estimation of shear modulus by magnetic resonance elastography: A computational comparative study

BACKGROUND

Magnetic resonance (MR) elastography is a non-destructive method of measuring biological tissue and is conducive to the early detection of tumors. Researchers usually set different assumptions according to different research objects, then establish and solve wave equations to estimate the shear modulus. Establishing a more reasonable model for a measured object estimates a more accurate shear modulus. Different assumptions of the mathematical model, and the method used to solve the wave equation causes deviation of the estimation.

OBJECTIVE

This study focused on shear modulus deviations caused by differences in calculation methods. The author demonstrated a method to ensure that the measuring range of the selected reconstruction algorithm with selected drive frequency covers the elasticity range of the target tissue. It is hoped to arouse the interest of researchers to introduce new transform domain methods to the field of MR elastography.

METHOD

In linear, isotropic and local homogeneity assumptions, the typical representative of two different calculation methods are algebraic inversion of the differential equation (AIDE) algorithm and local frequency elastography (LFE) algorithm. To compare the accuracy of these calculation methods, the author adopted a digital phantom that can set the parameter values accurately. It is assumed that the phantom tissue was linear and isotropic, and that the driving wave was sinusoidal. The displacement distribution of waves in the tissue was calculated by the finite element simulation method in two different resolutions with the signal-to-noise ratio (SNR) set to 40 dB and the threshold of relative mean error (RME) no more than 10%. The wavelength-to-pixel-size ratios of the two methods under the setting threshold of RME were compared.

RESULTS

The lower threshold of wavelength-to-pixel-size ratio for AIDE was close to 10, while that for LFE was nearly 2 (the limitation of Shannon's law) under the setting precision. Thus, the measuring range of the AIDE method was less than that of LFE at the same experimental conditions.

CONCLUSION

The driving frequency selection range of the spatial frequency domain method is wider than that of the spatial domain method. It is worthwhile for researchers to devote more time to introducing new transformation domain method for MR elastography.

Nonlinear Inversion MR Elastography With Low-Frequency Actuation

Magnetic resonance elastography (MRE) has been developed to noninvasively reconstruct mechanical properties for tissue and tissue-like materials over a frequency range of 10 ~200 Hz. In this work, low frequency (1~1.5 Hz) MRE activations were employed to estimate mechanical property distributions of simulated data and experimental phantoms. Nonlinear inversion (NLI) MRE algorithms based on viscoelastic and poroelastic material models were used to solve the inverse problems and recover images of the shear modulus and hydraulic conductivity. Data from a simulated phantom containing an inclusion with property contrast was carried out to study the feasibility of our low frequency actuated approach. To verify the stability of NLI algorithms for low frequency actuation, different levels of synthetic noise were added to the displacement data. Spatial distributions and property values were recovered well for noise level less than 5%. For the presented experimental phantom reconstructions with regularizations, the computed storage moduli from viscoelastic and poroelastic MRE gave similar results. Contrast was detected between inclusions and background in recovered hydraulic conductivity images. Results and findings confirm the feasibility of future in vivo neuroimaging examinations using natural cerebrovascular pulsations at cardiac frequencies, which can eliminate specialized equipment for high frequency actuation.

Analytical Minimization-Based Regularized Subpixel Shear-Wave Tracking for Ultrasound Elastography

Ultrasound elastography is a convenient and affordable method for imaging mechanical properties of tissue, which are often correlated with pathologies. An emerging novel elastography technique applies an external acoustic radiation force to generate a shear wave in the tissue and uses ultrasound imaging to track the shear wave. Accurate tracking of the small tissue motion is a critical step in shear-wave elastography (SWE), but it is challenging due to various sources of noise in the ultrasound data. We formulate tissue displacement estimation as an optimization problem and propose two computationally efficient approaches to estimate the displacement field. The first algorithm is referred to as dynamic programming analytic minimization (DPAM), which utilizes first-order Taylor series expansion of the highly nonlinear cost function to allow for its efficient optimization, and was previously proposed for quasistatic elastography. The second algorithm is a novel technique that utilizes second-order derivatives of the nonlinear cost function. We call the new algorithm second-order analytic minimization elastography (SESAME). We compare DPAM and SESAME to the standard normalized cross correlation (NCC) approach in the context of displacement and speed estimation of wave propagation in SWE. The results of micrometer-order displacement estimation in a uniform simulation phantom illustrate that SESAME outperforms DPAM, which in turn outperforms NCC in terms of signal-to-noise ratio (SNR) and jitter. In addition, the relative difference between true and reconstructed shear modulus (averaged over excitations at different focal depths and several scatterer realizations at each depth) is approximately 3.41%, 1.12%, and 1.01%, respectively, for NCC, DPAM, and SESAME. The performance of the proposed methods is also assessed with real data acquired using a tissue-mimicking phantom, wherein, in comparison to NCC, DPAM and SESAME improve the SNR of displacement estimates by 7.6 and 9.5 dB, respectively. Experimental results on a tissue-mimicking phantom also show that shear modulus reconstruction substantially improved with the proposed DPAM technique over NCC and with some further improvement achieved by utilizing the second-order Taylor series approximation in SESAME instead of the first-order DPAM.

A 2D finite element model for shear wave propagation in biological soft tissues: Application to magnetic resonance elastography

Dynamic elastography is a virtual palpation tool that aims at investigating the mechanical response of biological soft tissues in vivo. The objective of this study is to develop a finite element model (FEM) with low computational cost for reproducing realistically wave propagation for magnetic resonance elastography in heterogeneous soft tissues. Based on the first-order shear deformation theory for moderately thick structures, this model is developed and validated through comparison with analytical formulations of wave propagating in heterogeneous, viscoelastic infinite medium. This 2D-FEM is then compared to experimental data and a 3D-FEM using a commercial software. Our FEM is a powerful promising tool for investigations of magnetic resonance elastography.