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

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

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

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.     

SNR Weighting for Shear Wave Speed Reconstruction in Tomoelastography

In tomoelastography, to achieve a final wave speed map by combining reconstructions obtained from all spatial directions and excitation frequencies, the use of weights is inevitable. Here, a new weighting scheme, which maximizes the signal‐to‐noise ratio (SNR) of the final wave speed map, has been proposed. To maximize the SNR of the final wave speed map, the use of squares of estimated SNR values of reconstructed individual maps has been proposed. Therefore, derivations of the SNR of the reconstructed wave speed maps have become necessary. Considering the noise on the complex MRI signal, the SNR of the reconstructed wave speed map was formulated by an analytical approach assuming a high SNR, and the results were verified using Monte Carlo simulations (MCSs). It has been assumed that the noise remains approximately Gaussian when the image SNR is high enough, despite the nonlinear operations in tomoelastography inversion. Hence, the SNR threshold was determined by comparing the SNR computed by MCSs and analytical approximations. The weighting scheme was evaluated for accuracy, spatial resolution and SNR performances on simulated phantoms. MR elastography (MRE) experiments on two different phantoms were conducted. Wave speed maps were generated for simulated 3D human abdomen MRE data and experimental human abdomen MRE data. The simulation results demonstrated that the SNR‐weighted inversion improved the SNR performance of the wave speed map by a factor of two compared to the performance of the original (i.e., amplitude‐weighted) reconstruction. In the case of a low SNR, no bias occurred in the wave speed map when SNR weighting was used, whereas 10% bias occurred when the original weighting (i.e., amplitude weighting) was used. Thus, while not altering the accuracy or spatial resolution of the wave speed map with the proposed weighting method, the SNR of the wave speed map has been significantly improved.     

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.     

Unipolar MR elastography: Theory,numerical analysis and implementation

In MR elastography (MRE), zeroth moment balanced motion‐encoding gradients (MEGs) are incorporated into MRI sequences to induce a phase shift proportional to the local displacement caused by external actuation. To maximize the signal‐to‐noise ratio (SNR), fractional encoding is employed, i.e., the MEG duration is reduced below the wave period. Here, gradients encode primarily the velocity of the motion‐reducing encoding efficiency. Thus, in GRE‐MRE, T₂* decay and motion sensitivity have to be balanced, imposing a lower limit on repetition times (TRs). We propose to use a single trapezoidal gradient, a “unipolar gradient”, to directly encode spin displacement. Such gradients cannot be used in conventional sequences as they exhibit a large zeroth moment and dephase magnetization. By time‐reversing a spoiled SSFP sequence, the spoiling gradient becomes an efficient unipolar MEG. The proposed “unipolar MRE” technique benefits from this approach in three ways: first, displacement encoding is split over multiple TRs increasing motion sensitivity; second, spoiler and MEG coincide, allowing a reduction in TR; third, motion sensitivity of a typical unipolar lobe is of an order of magnitude higher than a bipolar MEG of equal duration. In this work, motion encoding using unipolar MRE is analyzed using the extended phase graph (EPG) formalism with a periodic motion propagator. As an approximation, the two‐transverse TR approximation for diffusion‐weighted SSFP is extended to incorporate cyclic motion. A complex encoding efficiency metric is introduced to compare the displacement fields of unipolar and conventional GRE‐MRE sequences in both magnitude and phase. The derived theoretical encoding equations are used to characterize the proposed sequence using an extensive parameter study. Unipolar MRE is validated against conventional GRE‐MRE in a phantom study showing excellent agreement between measured displacement fields. In addition, unipolar MRE yields significantly increased octahedral shear strain‐SNR relative to conventional GRE‐MRE and allows for the recovery of high stiffness inclusions, where conventional GRE‐MRE fails.     

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.     

Magnetic resonance elastography of the lungs: A repeatability and reproducibility study

Lung diseases are one of the leading causes of death worldwide, from which four million people die annually. Lung diseases are associated with changes in the mechanical properties of the lungs. Several studies have shown the feasibility of using magnetic resonance elastography (MRE) to quantify the lungs' shear stiffness. The aim of this study is to investigate the reproducibility and repeatability of lung MRE, and its shear stiffness measurements, obtained using a modified spin echo‐echo planar imaging (SE‐EPI) MRE sequence. In this study, 21 healthy volunteers were scanned twice by repositioning the volunteers to image right lung both at residual volume (RV) and total lung capacity (TLC) to assess the reproducibility of lung shear stiffness measurements. Additionally, 19 out of the 21 volunteers were scanned immediately without moving the volunteers to test the repeatability of the modified SE‐EPI MRE sequence. A paired t‐test was performed to determine the significant difference between stiffness measurements obtained at RV and TLC. Concordance correlation and Bland–Altman's analysis were performed to determine the reproducibility and repeatability of the SE‐EPI MRE‐derived shear stiffness measurements. The SE‐EPI MRE sequence is highly repeatable with a concordance correlation coefficient (CCC) of 0.95 at RV and 0.96 at TLC. Similarly, the stiffness measurements obtained across all volunteers were highly reproducible with a CCC of 0.95 at RV and 0.92 at TLC. The mean shear stiffness of the lung at RV was 0.93 ± 0.22 kPa and at TLC was 1.41 ± 0.41 kPa. TLC showed a significantly higher mean shear stiffness (P = 0.0004) compared with RV. Lung MRE stiffness measurements obtained using the SE‐EPI sequence were reproducible and repeatable, both at RV and TLC. Lung shear stiffness changes across respiratory cycle with significantly higher stiffness at TLC than RV.     

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.     

Three-dimensional analysis of shear wave propagation observed by in vivo magnetic resonance elastography of the brain

Dynamic magnetic resonance elastography (MRE) is a non-invasive method for the quantitative determination of the mechanical properties of soft tissues in vivo. In MRE, shear waves are generated in the tissue and visualized using phase-sensitive MR imaging methods. The resulting two-dimensional (2-D) wave images can reveal in-plane elastic properties when possible geometrical biases of the wave patterns are taken into account. In this study, 3-D MRE experiments of in vivo human brain are analyzed to gain knowledge about the direction of wave propagation and to deduce in-plane elastic properties. The direction of wave propagation was determined using a new algorithm which identifies minimal wave velocities along rays from the surface into the brain. It was possible to quantify biases of the elastic parameters due to projections onto coronal, sagittal and transversal image planes in 2-D MRE. It was found that the in-plane shear modulus is increasingly overestimated when the image slice is displaced from narrow slabs of 2-5cm through the center of the brain. The mean shear modulus of the brain was deduced from 4-D wave data with about 3.5kPa. Using the proposed slice positions in 2-D MRE, this shear modulus can be reproduced with an acceptable error within a fraction of the full 3-D examination time.     

Ristretto MRE: A generalized multi‐shot GRE‐MRE sequence

In order to acquire consistent k‐space data in MR elastography, a fixed temporal relationship between the MRI sequence and the underlying period of the wave needs to be ensured. To this end, conventional GRE‐MRE enforces synchronization through repeated triggering of the transducer and forcing the sequence repetition time to be equal to an integer multiple of the wave period. For wave frequencies below 100 Hz, however, this leads to prolonged acquisition times, as the repetition time scales inversely with frequency. A previously developed multi‐shot approach (eXpresso MRE) to multi‐slice GRE‐MRE tackles this issue by acquiring an integer number of slices per wave period, which allows acquisition to be accelerated in typical scenarios by a factor of two or three. In this work, it is demonstrated that the constraints imposed by the eXpresso scheme are overly restrictive. We propose a generalization of the sequence in three steps by incorporating sequence delays into imaging shots and allowing for interleaved wave‐phase acquisition. The Ristretto scheme is compared in terms of imaging shot and total scan duration relative to eXpresso and conventional GRE‐MRE and is validated in three different phantom studies. First, the agreement of measured displacement fields in different stages of the sequence generalization is shown. Second, performance is compared for 25, 36, 40, and 60 Hz actuation frequencies. Third, the performance is assessed for the acquisition of different numbers of slices (13 to 17). In vivo feasibility is demonstrated in the liver and the breast. Here, Ristretto is compared with an optimized eXpresso sequence, leading to scan accelerations of 15% and 5%, respectively, without compromising displacement field and stiffness estimates in general. The Ristretto concept allows us to choose imaging shot durations on a fine grid independent of the number of slices and the wave frequency, permitting 2‐ to 4.5‐fold acceleration of conventional GRE‐MRE acquisitions.     

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.     

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.     

The link between tissue elasticity and thermal dose in vivo

The objective of this study was to investigate in vivo the relationship between stiffness and thermal dose. For this purpose, shear wave elastography (SWE)--a novel ultrasound-based technique for real-time mapping of the stiffness of biological soft tissues--is performed in temperature-controlled experiments. Experiments were conducted on nine anesthetized rats. Their right leg was put in a thermo-regulated waterbath. The right leg of each animal was heated at one particular temperature between 38 °C and 48.5 °C for 15 min to 3 h. Shear waves were generated in the muscle using the acoustic radiation force induced by a linear ultrasonic probe. The shear wave propagation was imaged in real time by the probe using an ultrafast scanner prototype (10,000 frames s?1). The local tissue stiffness was derived from the shear wave speed. Two optical fiber sensors were inserted into the muscle to measure in situ the temperature. Stiffness was found to increase strongly during the experiments. When expressed as a function of the thermal dose, the stiffness curves were found to be the same for all experiments. A thermal dose threshold was found at 202 min for an eightfold stiffness increase. Finally, the time-temperature relationship was established for different stiffness ratios. The slope of the time-temperature relationship based on stiffness measurements was found identical to the one obtained for cell death in the seminal paper on the thermal dose by Sapareto and Dewey in 1984 (Int. J. Radiat. Oncol. Biol. Phys. 10 787-800). The present results highlight the stiffness increase as a good indicator of thermal necrosis. SWE imaging can be used in vivo for necrosis threshold determination in thermal therapy.     

Preliminary assessment of one-dimensional MR elastography for use in monitoring focused ultrasound therapy

The purpose of this work is to assess a fast technique that measures tissue stiffness and temperature during focused ultrasound thermal therapy (FUS). A one-dimensional (1D) MR elastography (MRE) pulse sequence was evaluated for the purpose of obtaining rapid measurements of thermally induced changes in tissue stiffness and temperature for monitoring FUS treatments. The accuracy of the 1D measurement was studied by comparing tissue displacements measured by 1D MRE with those measured by the well-established 2D MRE pulse sequence. The reproducibility of the 1D MRE measurement was assessed, in gel phantoms and ex vivo porcine tissue, for varied FUS intensity levels (31.5-199.9 W cm(-2)) and over a range of displacements at the focus (0.1-1 microm). Temperature elevations in agarose gel phantoms were measured using 1D MRE and calibrated using fiberoptic-thermometer-based measurements. The 1D MRE displacement measurements are highly correlated with those obtained with the 2D technique (R(2) = 0.88-0.93), indicating that 1D MRE can successfully measure tissue displacement. Ten repeated trials at each FUS power level yielded a minimum detectable displacement change of 0.2 microm in phantoms and 0.4 microm in tissue (at 95% confidence level). The 1D MRE temperature measurements correlated well with temperature changes measured simultaneously with fiberoptic thermometers (R(2) = 0.97). The 1D MRE technique is capable of detecting tissue displacements as low as 0.4 microm, which is an order of magnitude smaller than 5 microm displacements expected during FUS therapy (Le et al 2005 AIP Conf. Proc.: Ther. Ultrasound 829 186-90). Additionally, 1D MRE was shown to provide adequate measurements of temperature elevations in tissue. These findings indicate that 1D MRE may be an effective tool for monitoring FUS treatments.     

Magnetic resonance elastography of the lumbar back muscles: A preliminary study

Back pain is associated with increased lumbar paraspinal muscle (LPM) stiffness identified by manual palpation and strain elastography. Recently, magnetic resonance elastography (MRE) has allowed the stiffness of muscle to be characterized noninvasively in vivo, providing quantitative 3D stiffness maps (elastograms). The aim of this study was to characterize the stiffness (shear modulus, SM) of the LPM (multifidus and erector spinae) using MRE. MRE of the lumbar region was performed on seven adults in supine position. MRE was acquired in three muscular states: relaxed with outstretched legs, stretched with passive pelvis flexion, and contracted with outstretched legs and tightened trunk muscles. The mean SM was measured within a region of interest manually defined in the multifidus, erector spinae, and the entire paraspinal compartment. The intermuscular difference and the effects of stretching and contraction were assessed by ANOVA and t‐tests. At rest, the mean SM of the paraspinal compartment was 1.6 ± 0.2 kPa. It increased significantly with stretching to 1.65 ± 0.3 kPa, and with contraction to 2.0 ± 0.7 kPa. Irrespective of muscular state, the erector spinae was significantly stiffer than the multifidus. The multifidus underwent proportionally higher stiffness changes from rest to contraction and stretching. MRE can be used to measure the stiffness of the LPM in different muscular states. We hypothesize that, irrespective of posture, the erector spinae behaves as semi‐rigid beam, and ensures permanent stiffness of the spine. The multifidus behaves as an adaptable muscle that provides segmental flexibility to the spine and tunes the spine stiffness. Clin. Anat. 31:514–520, 2018. © 2018 Wiley Periodicals, Inc.     

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.     

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.     

Sonoelastographic imaging of interference patterns for estimation of the shear velocity of homogeneous biomaterials

The shear wave velocity is one of a few important parameters that characterize the mechanical properties of bio-materials. In this paper, two noninvasive methods are proposed to measure the shear velocity by inspecting the shear wave interference patterns. In one method, two shear wave sources are placed on the opposite two sides of a sample, driven by the identical sinusoidal signals. The shear waves from the two sources interact to create interference patterns, which are visualized by the vibration sonoelastography technique. The spacing between the pattern bands equals half of the shear wavelength. The shear velocity can be obtained by taking the product of the wavelength and the frequency. An alternative method is to drive the two vibration sources at slightly different frequencies. In this case, the interference patterns no longer remain stationary. It is proved that the apparent velocity of the moving patterns is proportional to the shear velocity in the medium. Since the apparent velocity of the patterns can be measured by analysing the video sequence, the shear velocity can be obtained thereafter. These approaches are validated by a conventional shear wave time-of-flight approach, and they are accurate within 4% on various homogeneous tissue-mimicking phantoms.     

激励声束聚焦形状对剪切波速度估计的影响

为了探究运用现有的激励及检测技术实现区域弹性成像的可行性,以及在声辐射力脉冲激励下,聚焦形状、检测点的空间位置等因素对局部剪切波速度估算的影响,运用Filed II以及有限元方法对声辐射力及组织响应情况进行了研究,并选用时域互相关方法对局部剪切波速度进行估计。结果表明,随着检测灵敏度的提高,可以实现发射2~3次激励脉冲即可对检测组织进行区域成像;使用现有的剪切波速度估计方法,浅表处与焦点深度处的剪切波速度估计值相差可达12.7%;f number（聚焦深度/孔径）增大,可以增大成像面积,但是剪切波速度估计的准确性降低。鉴于此,提出波速修正算法,修正后波速估计值更接近理论值。以f number为1.3为例,修正后的平均误差由3.2%降至0.6%。结果证明,现有的声辐射力激励及修正后的剪切波速度估计方法可用于区域弹性成像。