磁共振弹性成像的初步实验研究

目的：研究磁共振弹性成像（MRE）技术。方法：研制外部激发装置，设计成像脉冲序列，制作模拟人体软组织的体模。激发装置由序列控制，于体模表面产生低频率剪切波。脉冲序列采用梯度回波序列，在x、y或z轴上施加运动敏感梯度（MSG）。剪切波导致的介质内的周期性移位可使接收信号产生周期性相位位移，从测得的相位位移就能计算出每个体素的移位值，直接显示介质内剪切波的传播。通过调整相位偏置，获得一个完整周期内剪切波的动态传播图像。相位图经局部频率估算法（LFE）处理后计算出量化的弹性模景图。实验采用浓度为1．0％和1．5％不同弹性的琼脂凝胶体模，激发频率分别采用150Hz、200Hz、250Hz和300Hz。结果：MRE的相位图显示了剪切波在体模内的传播，剪切波的波长随激发频率和体模弹性变化。波长与激发频率呈反比，与体模弹性呈正比。剪切波的波长在不同激发频率和不同浓度体模之间呈严格的比例关系。计算出的弹性模量图清楚显示了两种浓度介质的弹性对比。结论：MRE的相位图可显示剪切波在介质内的传播，弹性图可量化和显示介质的弹性模量。     

Frequency-dependent viscoelastic parameters of mouse brain tissue estimated by MR elastography

Viscoelastic properties of mouse brain tissue were estimated non-invasively, in vivo, using magnetic resonance elastography (MRE) at 4.7 T to measure the dispersive properties of induced shear waves. Key features of this study include (i) the development and application of a novel MR-compatible actuation system which transmits vibratory motion into the brain through an incisor bar, and (ii) the investigation of the mechanical properties of brain tissue over a 1200 Hz bandwidth from 600-1800 Hz. Displacement fields due to propagating shear waves were measured during continuous, harmonic excitation of the skull. This protocol enabled characterization of the true steady-state patterns of shear wave propagation. Analysis of displacement fields obtained at different frequencies indicates that the viscoelastic properties of mouse brain tissue depend strongly on frequency. The average storage modulus (G') increased from approximately 1.6 to 8 kPa over this range; average loss modulus (G″) increased from approximately 1 to 3 kPa. Both moduli were well approximated by a power-law relationship over this frequency range. MRE may be a valuable addition to studies of disease in murine models, and to pre-clinical evaluations of therapies. Quantitative measurements of the viscoelastic parameters of brain tissue at high frequencies are also valuable for modeling and simulation of traumatic brain injury.     

In vivo brain viscoelastic properties measured by magnetic resonance elastography

Magnetic resonance elastography (MRE) is a non-invasive imaging technique used to visualise and quantify mechanical properties of tissue, providing information beyond what can be currently achieved with standard MR sequences and could, for instance, provide new insight into pathological processes in the brain. This study uses the MRE technique at 3 T to extract the complex shear modulus for in vivo brain tissue utilizing a full three-dimensional approach to reconstruction, removing contributions of the dilatational wave by application of the curl operator. A calibrated phantom is used to benchmark the MRE measurements, and in vivo results are presented for healthy volunteers. The results provide data for in vivo brain storage modulus (G'), finding grey matter (3.1 kPa) to be significantly stiffer than white matter (2.7 kPa). The first in vivo loss modulus (G') measurements show no significant difference between grey matter (2.5 kPa) and white matter (2.5 kPa).     

Tissue characterization using magnetic resonance elastography: preliminary results

The well-documented effectiveness of palpation as a diagnostic technique for detecting cancer and other diseases has provided motivation for developing imaging techniques for noninvasively evaluating the mechanical properties of tissue. A recently described approach for elasticity imaging, using propagating acoustic shear waves and phase-contrast MRI, has been called magnetic resonance elastography (MRE). The purpose of this work was to conduct preliminary studies to define methods for using MRE as a tool for addressing the paucity of quantitative tissue mechanical property data in the literature. Fresh animal liver and kidney tissue specimens were evaluated with MRE at multiple shear wave frequencies. The influence of specimen temperature and orientation on measurements of stiffness was studied in skeletal muscle. The results demonstrated that all of the materials tested (liver, kidney, muscle and tissue-simulating gel) exhibit systematic dependence of shear stiffness on shear rate. These data are consistent with a viscoelastic model of tissue mechanical properties, allowing calculation of two independent tissue properties from multiple-frequency MRE data: shear modulus and shear viscosity. The shear stiffness of tissue can be substantially affected by specimen temperature. The results also demonstrated evidence of shear anisotropy in skeletal muscle but not liver tissue. The measured shear stiffness in skeletal muscle was found to depend on both the direction of propagation and polarization of the shear waves.     

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

The mechanical model commonly used in magnetic resonance elastography (MRE) is linear elasticity. However, soft tissue may exhibit frequency- and direction-dependent (FDD) shear moduli in response to an induced excitation causing a purely linear elastic model to provide an inaccurate image reconstruction of its mechanical properties. The goal of this study was to characterize the effects of reconstructing FDD data using a linear elastic inversion (LEI) algorithm. Linear and FDD phantoms were manufactured and LEI images were obtained from time-harmonic MRE acquisitions with variations in frequency and driving signal amplitude. LEI responses to artificially imposed uniform phase shifts in the displacement data from both purely linear elastic and FDD phantoms were also evaluated. Of the variety of FDD phantoms considered, LEI appeared to tolerate viscoelastic data-model mismatch better than deviations caused by poroelastic and anisotropic mechanical properties in terms of visual image contrast. However, the estimated shear modulus values were substantially incorrect relative to independent mechanical measurements even in the successful viscoelastic cases and the variations in mean values with changes in experimental conditions associated with uniform phase shifts, driving signal frequency and amplitude were unpredictable. Overall, use of LEI to reconstruct data acquired in phantoms with FDD material properties provided biased results under the best conditions and significant artifacts in the worst cases. These findings suggest that the success with which LEI is applied to MRE data in tissue will depend on the underlying mechanical characteristics of the tissues and/or organs systems of clinical interest.     

Non-invasive measurement of brain viscoelasticity using magnetic resonance elastography

The purpose of this work was to develop magnetic resonance elastography (MRE) for the fast and reproducible measurement of spatially averaged viscoelastic constants of living human brain. The technique was based on a phase-sensitive echo planar imaging acquisition. Motion encoding was orthogonal to the image plane and synchronized to intracranial shear vibrations at driving frequencies of 25 and 50 Hz induced by a head-rocker actuator. Ten time-resolved phase-difference wave images were recorded within 60 s and analyzed for shear stiffness and shear viscosity. Six healthy volunteers (six men; mean age 34.5 years; age range 25-44 years) underwent 23-39 follow-up MRE studies over a period of 6 months. Interindividual mean +/- SD shear moduli and shear viscosities were found to be 1.17 +/- 0.03 kPa and 3.1 +/- 0.4 Pas for 25 Hz and 1.56 +/- 0.07 kPa and 3.4 +/- 0.2 Pas for 50 Hz, respectively (P < or = 0.01). The intraindividual range of shear modulus data was 1.01-1.31 kPa (25 Hz) and 1.33-1.77 kPa (50 Hz). The observed modulus dispersion indicates a limited applicability of Voigt's model to explain viscoelastic behavior of brain parenchyma within the applied frequency range. The narrow distribution of data within small confidence intervals demonstrates excellent reproducibility of the experimental protocol. The results are necessary as reference data for future comparisons between healthy and pathological human brain viscoelastic data.     

Multifrequency inversion in magnetic resonance elastography

Time-harmonic shear wave elastography is capable of measuring viscoelastic parameters in living tissue. However, finite tissue boundaries and waveguide effects give rise to wave interferences which are not accounted for by standard elasticity reconstruction methods. Furthermore, the viscoelasticity of tissue causes dispersion of the complex shear modulus, rendering the recovered moduli frequency dependent. Therefore, we here propose the use of multifrequency wave data from magnetic resonance elastography (MRE) for solving the inverse problem of viscoelasticity reconstruction by an algebraic least-squares solution based on the springpot model. Advantages of the method are twofold: (i) amplitude nulls appearing in single-frequency standing wave patterns are mitigated and (ii) the dispersion of storage and loss modulus with drive frequency is taken into account by the inversion procedure, thereby avoiding subsequent model fitting. As a result, multifrequency inversion produces fewer artifacts in the viscoelastic parameter map than standard single-frequency parameter recovery and may thus support image-based viscoelasticity measurement. The feasibility of the method is demonstrated by simulated wave data and MRE experiments on a phantom and in vivo human brain. Implemented as a clinical method, multifrequency inversion may improve the diagnostic value of time-harmonic MRE in a large variety of applications.     

Algebraic Helmholtz inversion in planar magnetic resonance elastography

Magnetic resonance elastography (MRE) is an increasingly used noninvasive modality for diagnosing diseases using the response of soft tissue to harmonic shear waves. We present a study on the algebraic Helmholtz inversion (AHI) applied to planar MRE, demonstrating that the deduced phase speed of shear waves depends strongly on the relative orientations of actuator polarization, motion encoding direction and image plane as well as on the actuator plate size, signal-to-noise ratio and discretization of the wave image. Results from the numerical calculation of harmonic elastic waves due to different excitation directions and simulated plate sizes are compared to experiments on a gel phantom. The results suggest that correct phase speed can be obtained despite these largely uncontrollable influences, if AHI is based on out-of-plane displacements and the actuator is driven at an optimal frequency yielding an optimal pixel per wavelength resolution in the wave image. Assuming plane waves, the required number of pixels per wavelength depends only on the degree of noise.     

Noninvasive assessment of the rheological behavior of human organs using multifrequency MR elastography: a study of brain and liver viscoelasticity

MR elastography (MRE) enables the noninvasive determination of the viscoelastic behavior of human internal organs based on their response to oscillatory shear stress. An experiment was developed that combines multifrequency shear wave actuation with broad-band motion sensitization to extend the dynamic range of a single MRE examination. With this strategy, multiple wave images corresponding to different driving frequencies are simultaneously received and can be analyzed by evaluating the dispersion of the complex modulus over frequency. The technique was applied on the brain and liver of five healthy volunteers. Its repeatability was tested by four follow-up studies in each volunteer. Five standard rheological models (Maxwell, Voigt, Zener, Jeffreys and fractional Zener model) were assessed for their ability to reproduce the observed dispersion curves. The three-parameter Zener model was found to yield the most consistent results with two shear moduli mu(1) = 0.84 +/- 0.22 (1.36 +/- 0.31) kPa, mu(2) = 2.03 +/- 0.19 (1.86 +/- 0.34) kPa and one shear viscosity of eta = 6.7 +/- 1.3 (5.5 +/- 1.6) Pa s (interindividual mean +/- SD) in brain (liver) experiments. Significant differences between the rheological parameters of brain and liver were found for mu(1) and eta (P < 0.05), indicating that human brain is softer and possesses a higher viscosity than liver.     

Viscoelasticity-based MR elastography of skeletal muscle

An in vivo multifrequency magnetic resonance elastography (MRE) protocol was developed for studying the viscoelastic properties of human skeletal muscle in different states of contraction. Low-frequency shear vibrations in the range of 25-62.5 Hz were synchronously induced into the femoral muscles of seven volunteers and measured in a cross-sectional view by encoding the fast-transverse shear wave component parallel to the muscle fibers. The so-called springpot model was used for deriving two viscoelastic constants, μ and α, from the dispersion functions of the complex shear modulus in relaxed and in loaded muscle. Representing the shear elasticity parallel to the muscle fibers, μ increased in all volunteers upon contraction from 2.68 ± 0.23 kPa to 3.87 ± 0.50 kPa. Also α varied with load, indicating a change in the geometry of the mechanical network of muscle from relaxation (α = 0.253 ± 0.009) to contraction (α = 0.270 ± 0.009). These results provide a reference for a future assessment of muscular dysfunction using rheological parameters.     

MR elastography in a murine stroke model reveals correlation of macroscopic viscoelastic properties of the brain with neuronal density

The aim of this study was to investigate the influence of neuronal density on viscoelastic parameters of living brain tissue after ischemic infarction in the mouse using MR elastography (MRE). Transient middle cerebral artery occlusion (MCAO) in the left hemisphere was induced in 20 mice. In vivo 7‐T MRE at a vibration frequency of 900 Hz was performed on days 3, 7, 14 and 28 (n = 5 per group) after MCAO, followed by the analysis of histological markers, such as neuron counts (NeuN). MCAO led to a significant reduction in the storage modulus in the left hemisphere relative to contralateral values (p = 0.03) without changes over time. A correlation between storage modulus and NeuN in both hemispheres was observed, with correlation coefficients of R = 0.648 (p = 0.002, left) and R = 0.622 (p = 0.003, right). The loss modulus was less sensitive to MCAO, but correlated with NeuN in the left hemisphere (R = 0.764, p = 0.0001). In agreement with the literature, these results suggest that the shear modulus in the brain is reduced after transient ischemic insult. Furthermore, our study provides evidence that the in vivo shear modulus of brain tissue correlates with neuronal density. In diagnostic applications, MRE may thus have diagnostic potential as a tool for image‐based quantification of neurodegenerative processes. Copyright © 2013 John Wiley & Sons, Ltd.     

Magnetic resonance elastography: A review

Magnetic resonance elastography (MRE) is a rapidly developing technology for quantitatively assessing the mechanical properties of tissue. The technology can be considered to be an imaging‐based counterpart to palpation, commonly used by physicians to diagnose and characterize diseases. The success of palpation as a diagnostic method is based on the fact that the mechanical properties of tissues are often dramatically affected by the presence of disease processes, such as cancer, inflammation, and fibrosis. MRE obtains information about the stiffness of tissue by assessing the propagation of mechanical waves through the tissue with a special magnetic resonance imaging technique. The technique essentially involves three steps: (1) generating shear waves in the tissue, (2) acquiring MR images depicting the propagation of the induced shear waves, and (3) processing the images of the shear waves to generate quantitative maps of tissue stiffness, called elastograms. MRE is already being used clinically for the assessment of patients with chronic liver diseases and is emerging as a safe, reliable, and noninvasive alternative to liver biopsy for staging hepatic fibrosis. MRE is also being investigated for application to pathologies of other organs including the brain, breast, blood vessels, heart, kidneys, lungs, and skeletal muscle. The purpose of this review article is to introduce this technology to clinical anatomists and to summarize some of the current clinical applications that are being pursued. Clin. Anat. 23:497–511, 2010. © 2010 Wiley‐Liss, Inc.     

Microvasculature alters the dispersion properties of shear waves – a multi‐frequency MR elastography study

Magnetic Resonance Elastography (MRE) uses macroscopic shear wave propagation to quantify mechanical properties of soft tissues. Micro‐obstacles are capable of affecting the macroscopic dispersion properties of shear waves. Since disease or therapy can change the mechanical integrity and organization of vascular structures, MRE should be able to sense these changes if blood vessels represent a source for wave scattering. To verify this, MRE was performed to quantify alteration of the shear wave speed c_s due to the presence of vascular outgrowths using an aortic ring model. Eighteen fragments of rat aorta included in a Matrigel matrix (n=6 without outgrowths, n=6 with a radial outgrowth extent of ~600µm and n=6 with ~850µm) were imaged using a 7 Tesla MR scanner (Bruker, PharmaScan). High resolution anatomical images were acquired in addition to multi‐frequency MRE (ν = 100, 115, 125, 135 and 150 Hz). Average c_s was measured within a ring of ~900µm thickness encompassing the aorta and were normalized to c_s0 of the corresponding Matrigel. The frequency dependence was fit to the power law model c_s ~ν^y. After scanning, optical microscopy was performed to visualize outgrowths. Results demonstrated that in presence of vascular outgrowths (1) normalized c_s significantly increased for the three highest frequencies (Kruskal‐Wallis test, P = 0.0002 at 125 Hz and P = 0.002 at 135 Hz and P = 0.003 at 150 Hz) but not for the two lowest (Kruskal‐Wallis test, P = 0.63 at 100 Hz and P = 0.87 at 115 Hz), and (2) normalized c_s followed a power law behavior not seen in absence of vascular outgrowths (ANOVA test, P < 0.0001). These results showed that vascular outgrowths acted as micro‐obstacles altering the dispersion relationships of propagating shear waves and that MRE could provide valuable information about microvascular changes. Copyright © 2015 John Wiley & Sons, Ltd.     

Calculation of shear stiffness in noise dominated magnetic resonance elastography data based on principal frequency estimation

Magnetic resonance elastography (MRE) is a non-invasive phase-contrast-based method for quantifying the shear stiffness of biological tissues. Synchronous application of a shear wave source and motion encoding gradient waveforms within the MRE pulse sequence enable visualization of the propagating shear wave throughout the medium under investigation. Encoded shear wave-induced displacements are then processed to calculate the local shear stiffness of each voxel. An important consideration in local shear stiffness estimates is that the algorithms employed typically calculate shear stiffness using relatively high signal-to-noise ratio (SNR) MRE images and have difficulties at an extremely low SNR. A new method of estimating shear stiffness based on the principal spatial frequency of the shear wave displacement map is presented. Finite element simulations were performed to assess the relative insensitivity of this approach to decreases in SNR. Additionally, ex vivo experiments were conducted on normal rat lungs to assess the robustness of this approach in low SNR biological tissue. Simulation and experimental results indicate that calculation of shear stiffness by the principal frequency method is less sensitive to extremely low SNR than previously reported MRE inversion methods but at the expense of loss of spatial information within the region of interest from which the principal frequency estimate is derived.     

时域有限差分法研究长骨中的超声导波

目的基于时域有限差分(finite-difference time-domain,FDTD)法研究长骨弹性模量对超声导波传播特性的影响,为早期骨质疲劳的超声导波检测方法提供理论依据。方法对长骨进行理论建模,并对不同弹性模量下的长骨模型进行FDTD仿真;通过对仿真数据的分析,计算出不同模式导波的相速度、群速度、中心频率及能量等特性参数。结果长骨的弹性模量与超声导波的传播特性参数密切相关,各个导波模式的相速度、群速度、中心频率和能量均随弹性模量的减小而减小,其中L(0,5)模式的变化趋势最为显著。结论超声导波的传播特性参量能够反映出长骨弹性模量的变化,进而为长骨早期疲劳诊断提供了一种可能的方法。     

Quantitative sonoelastography for the in vivo assessment of skeletal muscle viscoelasticity

A novel quantitative sonoelastography technique for assessing the viscoelastic properties of skeletal muscle tissue was developed. Slowly propagating shear wave interference patterns (termed crawling waves) were generated using a two-source configuration vibrating normal to the surface. Theoretical models predict crawling wave displacement fields, which were validated through phantom studies. In experiments, a viscoelastic model was fit to dispersive shear wave speed sonoelastographic data using nonlinear least-squares techniques to determine frequency-independent shear modulus and viscosity estimates. Shear modulus estimates derived using the viscoelastic model were in agreement with that obtained by mechanical testing on phantom samples. Preliminary sonoelastographic data acquired in healthy human skeletal muscles confirm that high-quality quantitative elasticity data can be acquired in vivo. Studies on relaxed muscle indicate discernible differences in both shear modulus and viscosity estimates between different skeletal muscle groups. Investigations into the dynamic viscoelastic properties of (healthy) human skeletal muscles revealed that voluntarily contracted muscles exhibit considerable increases in both shear modulus and viscosity estimates as compared to the relaxed state. Overall, preliminary results are encouraging and quantitative sonoelastography may prove clinically feasible for in vivo characterization of the dynamic viscoelastic properties of human skeletal muscle.     

Magnetic resonance elastography of skeletal muscle deep tissue injury

The current state‐of‐the‐art diagnosis method for deep tissue injury in muscle, a subcategory of pressure ulcers, is palpation. It is recognized that deep tissue injury is frequently preceded by altered biomechanical properties. A quantitative understanding of the changes in biomechanical properties preceding and during deep tissue injury development is therefore highly desired. In this paper we quantified the spatial–temporal changes in mechanical properties upon damage development and recovery in a rat model of deep tissue injury. Deep tissue injury was induced in nine rats by two hours of sustained deformation of the tibialis anterior muscle. Magnetic resonance elastography (MRE), T₂‐weighted, and T₂‐mapping measurements were performed before, directly after indentation, and at several timepoints during a 14‐day follow‐up. The results revealed a local hotspot of elevated shear modulus (from 3.30 ± 0.14 kPa before to 4.22 ± 0.90 kPa after) near the center of deformation at Day 0, whereas the T₂ was elevated in a larger area. During recovery there was a clear difference in the time course of the shear modulus and T₂. Whereas T₂ showed a gradual normalization towards baseline, the shear modulus dropped below baseline from Day 3 up to Day 10 (from 3.29 ± 0.07 kPa before to 2.68 ± 0.23 kPa at Day 10, P < 0.001), followed by a normalization at Day 14. In conclusion, we found an initial increase in shear modulus directly after two hours of damage‐inducing deformation, which was followed by decreased shear modulus from Day 3 up to Day 10, and subsequent normalization. The lower shear modulus originates from the moderate to severe degeneration of the muscle. MRE stiffness values were affected in a smaller area as compared with T₂. Since T₂ elevation is related to edema, distributing along the muscle fibers proximally and distally from the injury, we suggest that MRE is more specific than T₂ for localization of the actual damaged area.     

Characterization of the linearly viscoelastic behavior of human tympanic membrane by nanoindentation

Human tympanic membrane (or eardrum) is composed of three membrane layers with collagen fibers oriented in the radial and circumferential directions, and exhibits viscoelastic behavior with membrane (or in-plane) properties different from through-thickness (or out-of-plane) properties. Due to the interaction of bundled fibers and ground substance, which is inhomogeneous, these properties could change with locations. In this paper, we use nanoindentation techniques to measure the viscoelastic functions of four quadrants of tympanic membrane (TM). For measurement of in-plane Young’s relaxation modulus we fixed a sectioned quadrant of the TM on a circular hole and used a spherical nanoindenter tip to apply force at the center of the suspended circular portion of the specimen. An inverse problem solving methodology was employed using finite element method to determine the average in-plane Young’s relaxation modulus of the TM quadrant. Results indicate that the in-plane steady-state Young’s relaxation modulus for four quadrants of the TM does not vary significantly. However, a variation of the modulus from 25.73 MPa to 37.8 MPa was observed with measurements from different individuals. For measurement of Young’s relaxation modulus in the through-thickness direction a spherical indenter tip was used to indent into different locations on the surface of the TM specimen supported by a substrate. Viscoelastic contact mechanics analysis of the load–displacement curve, representative primarily of the through-thickness stiffness of the TM, was conducted to extract the Young’s relaxation modulus in the out-of-plane direction. Results indicate a wide variation in steady-state Young’s relaxation modulus, from 2 MPa to 15 MPa, in the through-thickness direction over the TM.     

In vivo elasticity measurements of extremity skeletal muscle with MR elastography

MR elastography (MRE) has been shown to be capable of non-invasively measuring tissue elasticity even in deep-lying regions. Although limited studies have already been published examining in vivo muscle elasticity, it is still not clear over what range the in vivo elasticity values vary. The present study intends to produce further information by examining four different skeletal muscles in a group of 12 healthy volunteers in the age range of 27-38 years. The examinations were performed in the biceps brachii, the flexor digitorum profundus, the soleus and the gastrocnemius. The average shear modulus was determined to be 17.9 (+/- 5.5), 8.7 (+/- 2.8), 12.5 (+/- 7.3) and 9.9 (+/- 6.8) kPa for each muscle, respectively. To ascertain the reproducibility of the examination, the stiffness measurements in two volunteers were repeated seven times for the biceps brachii. These examinations yielded a mean shear modulus of 11.3 +/-.7 and 13.3 +/- 4.7 kPa for the two subjects. For elasticity reconstruction, an automated reconstruction algorithm is introduced which eliminates variation due to subjective manual image analysis. This study yields new information regarding the expected variation in muscle elasticity in a healthy population, and also reveals the expected variability of the MRE technique in skeletal muscle.