Assessment of elastic parameters of human skin using dynamic elastography

Sonoelastography and transient elastography are two ultrasound-based techniques that facilitate noninvasive characterization of the viscoelastic properties of soft tissues by investigating their response to shear mechanical excitation. Young's modulus is the principle assessment parameter. Because it defines local tissue stiffness, it is of major interest for the medical imaging and cosmetic industries as it could replace subjective palpation by yielding local, quantitative information. In this paper, we describe a new high-resolution device capable of measuring local Young's modulus in very thin layers (1-5 mm) and devoted to the in vivo evaluation of the elastic properties of human skin. It uses an ultrasonic probe (50 MHz) for tracking the displacements induced by a 300 Hz shear wave generated by a ring surrounding the transducer. The displacements are measured using a conventional cross-correlation technique between successive ultrasonic back-scattered echoes. First, this noninvasive technique has been experimentally proven to be accurate for investigating elasticity in different skin-mimicking phantoms. Second, data were acquired in vivo on human forearms. As expected, Young's modulus was found to be higher in the dermis than in the hypodermis and other soft tissues.     

Quantitative imaging of nonlinear shear modulus by combining static elastography and shear wave elastography

The study of new tissue mechanical properties such as shear nonlinearity could lead to better tissue characterization and clinical diagnosis. This work proposes a method combining static elastography and shear wave elastography to derive the nonlinear shear modulus by applying the acoustoelasticity theory in quasi-incompressible soft solids. Results demonstrate that by applying a moderate static stress at the surface of the investigated medium, and by following the quantitative evolution of its shear modulus, it is possible to accurately and quantitatively recover the local Landau (A) coefficient characterizing the shear nonlinearity of soft tissues.     

Application of 1-D transient elastography for the shear modulus assessment of thin-layered soft tissue: comparison with supersonic shear imaging technique

Elasticity estimation of thin-layered soft tissues has gained increasing interest propelled by medical applications like skin, corneal, or arterial wall shear modulus assessment. In this work, the authors propose one-dimensional transient elastography (1DTE) for the shear modulus assessment of thin-layered soft tissue. Experiments on three phantoms with different elasticities and plate thicknesses were performed. First, using 1DTE, the shear wave speed dispersion curve inside the plate was obtained and validated with finite difference simulation. No dispersive effects were observed and the shear wave speed was directly retrieved from time-of-flight measurements. Second, the supersonic shear imaging (SSI) technique (considered to be a gold standard) was performed. For the SSI technique, the propagating wave inside the plate is guided as a Lamb wave. Experimental SSI dispersion curves were compared with finite difference simulation and fitted using a generalized Lamb model to retrieve the plate bulk shear wave speed. Although they are based on totally different mechanical sources and induce completely different diffraction patterns for the shear wave propagation, the 1DTE and SSI techniques resulted in similar shear wave speed estimations. The main advantage of the 1DTE technique is that bulk shear wave speed can be directly retrieved without requiring a dispersion model.     

Assessment of viscous and elastic properties of sub-wavelength layered soft tissues using shear wave spectroscopy: theoretical framework and in vitro experimental validation

In elastography, quantitative imaging of soft tissue elastic properties is provided by local shear wave speed estimation. Shear wave imaging in a homogeneous medium thicker than the shear wavelength is eased by a simple relationship between shear wave speed and local shear modulus. In thin layered organs, the shear wave is guided and thus undergoes dispersive effects. This case is encountered in medical applications such as elastography of skin layers, corneas, or arterial walls. In this work, we proposed and validated shear wave spectroscopy as a method for elastic modulus quantification in such layered tissues. Shear wave dispersion curves in thin layers were obtained by finite-difference simulations and numerical solving of the boundary conditions. In addition, an analytical approximation of the dispersion equation was derived from the leaky Lamb wave theory. In vitro dispersion curves obtained from phantoms were consistent with numerical studies (deviation <1.4%). The least-mean-squares fitting of the dispersion curves enables a quantitative and accurate (error < 5% of the transverse speed) assessment of the elasticity. Dispersion curves were also found to be poorly influenced by shear viscosity. This phenomenon allows independent recovery of the shear modulus and the viscosity, using, respectively, the dispersion curve and the attenuation estimation along the propagation axis.     

超声编码激励在瞬时弹性成像检测中的应用

瞬时弹性(Transient Elastography,TE)成像广泛应用于肝硬化研究。然而,在临床应用中,对于肥胖病人,该方法很难实现对深度位置的瞬时剪切波进行检测。研究了将超声编码激励应用于瞬时弹性成像剪切波检测的可行性,选用7位巴克码进行编码检测研究。通过剪切波信噪比和检测穿透力两个指标对编码检测与传统短脉冲检测结果进行量化和对比。弹性仿体实验表明,编码检测可以提供比传统短脉冲检测更高的剪切波信噪比和检测深度。离体猪肝实验表明将编码激励应用于组织检测时同样可以实现高信噪比剪切波检测。这些结果表明编码检测应用于瞬时弹性成像检测是可行的,可以增加其检测深度。     

Shear modulus imaging with 2-D transient elastography

In previous work, we have shown that time-resolved 2-D transient elastography is a promising technique for characterizing the elasticity of soft tissues. It involves the measurement of the displacements induced by the propagation of low frequency (LF) pulsed shear waves in biological tissues. In this paper, we present a novel apparatus that contains a LF vibrating device surrounding a linear array of 128 ultrasonic transducers that performs ultrafast ultrasonic imaging (up to 10,000 frames/s) and that is able to follow in real time the propagation of a LF shear wave in the human body. The vibrating device is made of two rods, fixed to electromagnetic vibrators, that produce in the ultrasonic image area a large amplitude shear wave. The geometry has been chosen both to enhance the sensitivity and to create a quasi linear shear wave front in the imaging plane. An inversion algorithm is used to recover the shear modulus map from the spatio-temporal data, and the first experimental results obtained from tissue-equivalent materials are presented.     

Imaging feedback of histotripsy treatments using ultrasound shear wave elastography

Histotripsy is a cavitation-based ultrasound therapy that mechanically fractionates soft solid tissues into fluid-like homogenates. This paper investigates the feasibility of imaging the tissue elasticity change during the histotripsy process as a tool to provide feedback for the treatments. The treatments were performed on agar tissue phantoms and ex vivo kidneys using 3-cycle ultrasound pulses delivered by a 750-kHz therapeutic array at peak negative/positive pressure of 17/108 MPa and a repetition rate of 50 Hz. Lesions with different degrees of damage were created with increasing numbers of therapy pulses from 0 to 2000 pulses per treatment location. The elasticity of the lesions was measured with ultrasound shear wave elastography, in which a quasi-planar shear wave was induced by acoustic radiation force generated by the therapeutic array, and tracked with ultrasound imaging at 3000 frames per second. Based on the shear wave velocity calculated from the sequentially captured frames, the Young's modulus was reconstructed. Results showed that the lesions were more easily identified on the shear wave velocity images than on B-mode images. As the number of therapy pulses increased from 0 to 2000 pulses/location, the Young's modulus decreased exponentially from 22.1 ± 2.7 to 2.1 ± 1.1 kPa in the tissue phantoms (R2 = 0.99, N = 9 each), and from 33.0 ± 7.1 to 4.0 ± 2.5 kPa in the ex vivo kidneys (R2 = 0.99, N = 8 each). Correspondingly, the tissues transformed from completely intact to completely fractionated as examined via histology. A good correlation existed between the lesions' Young's modulus and the degree of tissue fractionation as examined with the percentage of remaining structurally intact cell nuclei (R2 = 0.91, N = 8 each). These results indicate that lesions produced by histotripsy can be detected with high sensitivity using shear wave elastography. Because the decrease in the tissue elasticity corresponded well with the morphological and histological change, this study provides a basis for predicting the local treatment outcomes from tissue elasticity change.     

On the effects of reflected waves in transient shear wave elastography

In recent years, novel quantitative techniques have been developed to provide noninvasive and quantitative stiffness images based on shear wave propagation. Using radiation force and ultrafast ultrasound imaging, the supersonic shear imaging technique allows one to remotely generate and follow a transient plane shear wave propagating in vivo in real time. The tissue shear modulus, i.e., its stiffness, can then be estimated from the shear wave local velocity. However, because the local shear wave velocity is estimated using a time-of- flight approach, reflected shear waves can cause artifacts in the estimated shear velocity because the incident and reflected waves propagate in opposite directions. Such effects have been reported in the literature as a potential drawback of elastography techniques based on shear wave speed, particularly in the case of high stiffness contrasts, such as in atherosclerotic plaque or stiff lesions. In this letter, we present our implementation of a simple directional filter, previously used for magnetic resonance elastography, which separates the forward- and backward-propagating waves to solve this problem. Such a directional filter could be applied to many elastography techniques based on the local estimation of shear wave speed propagation, such as acoustic radiation force imaging (ARFI), shearwave dispersion ultrasound vibrometry (SDUV), needle-based elastography, harmonic motion imaging, or crawling waves when the local propagation direction is known and high-resolution spatial and temporal data are acquired.     

Shear elasticity probe for soft tissues with 1-D transient elastography

Important tissue parameters such as elasticity can be deduced from the study of the propagation of low frequency shear waves. A new method for measuring the shear velocity in soft tissues is presented in this paper. Unlike conventional transient elastography in which the ultrasonic transducer and the low frequency vibrator are two separated parts, the new method relies on a probe that associates the vibrator and the transducer, which is built on the axis of the vibrator. This setup is easy to use. The low frequency shear wave is driven by the transducer itself that acts as a piston while it is used in pulse echo mode to acquire ultrasonic lines. The results obtained with the new method are in good agreement with those obtained with the conventional one.     

1-D elasticity assessment in soft solids from shear wave correlation: the time-reversal approach

One-channel time-reversal (TR) experiments allow focalization of waves in reverberant cavities. According to the Rayleigh criterion, the focal spot width is directly related to the wavelength and therefore depends on the mechanical properties of the medium. Thus, the general idea of this work is to extract quantitative estimations of these mechanical properties using a time-reversal approach based on cross-correlations of the wave field. An external source creates mechanical waves in the audible frequency range. One component of the vectorial field is measured along a line as function of time with signal processing developed in the field of 1-D elastography. The shear wavelength information is deduced from these mechanical waves using spatiotemporal correlations and interpreted in the frame of the time-reversal symmetry. The impact of wave attenuation in soft solids is reduced using a spatial average of the correlation field. The result is shown to be suitable for global elasticity estimation. The advantage is that the technique is almost independent of the source kind, shape, and time excitation function. This robustness as regard to shear wave source allows translation of this technique to applications in the medical field, including deep or moving organs.     

The role of viscosity in the impulse diffraction field of elastic waves induced by the acoustic radiation force

Several ultrasound-based techniques for the estimation of soft tissue elasticity are currently being investigated. Most of them study the medium response to dynamic excitations. Such responses are usually modeled in a purely elastic medium using a Green's function solution of the motion equation. However, elasticity by itself is not necessarily a discriminant parameter for malignancy diagnosis. Modeling viscous properties of tissues could also be of great interest for tumor characterization. We report in this paper an explicit derivation of the Green's function in a viscous and elastic medium taking into account shear, bulk, and coupling waves. From this theoretical calculation, 3D simulations of mechanical waves in viscoelastic soft tissues are presented. The relevance of the viscoelastic Green's function is validated by comparing simulations with experimental data. The experiments were conducted using the supersonic shear imaging (SSI) technique which dynamically and remotely excites tissues using acoustic radiation force. We show that transient shear waves generated with SSI are modeled very precisely by the Green's function formalism. The combined influences of out-of-plane diffraction, beam shape, and shear viscosity on the shape of transient waves are carefully studied as they represent a major issue in ultrasound-based viscoelasticity imaging techniques.     

Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography

The emergence of ultrafast frame rates in ultrasonic imaging has been recently made possible by the development of new imaging modalities such as transient elastography. Data acquisition rates reaching more than thousands of images per second enable the real-time visualization of shear mechanical waves propagating in biological tissues, which convey information about local viscoelastic properties of tissues. The first proposed approach for reaching such ultrafast frame rates consists of transmitting plane waves into the medium. However, because the beamforming process is then restricted to the receive mode, the echographic images obtained in the ultrafast mode suffer from a low quality in terms of resolution and contrast and affect the robustness of the transient elastography mode. It is here proposed to improve the beamforming process by using a coherent recombination of compounded plane-wave transmissions to recover high-quality echographic images without degrading the high frame rate capabilities. A theoretical model is derived for the comparison between the proposed method and the conventional B-mode imaging in terms of contrast, signal-to-noise ratio, and resolution. Our model predicts that a significantly smaller number of insonifications, 10 times lower, is sufficient to reach an image quality comparable to conventional B-mode. Theoretical predictions are confirmed by in vitro experiments performed in tissue-mimicking phantoms. Such results raise the appeal of coherent compounds for use with standard imaging modes such as B-mode or color flow. Moreover, in the context of transient elastography, ultrafast frame rates can be preserved while increasing the image quality compared with flat insonifications. Improvements on the transient elastography mode are presented and discussed.     

Estimation of polyvinyl alcohol cryogel mechanical properties with four ultrasound elastography methods and comparison with gold standard testings

Tissue-mimicking phantoms are very useful in the field of tissue characterization and essential in elastography for the purpose of validating motion estimators. This study is dedicated to the characterization of polyvinyl alcohol cryogel (PVA-C) for these types of applications. A strict fabrication procedure was defined to optimize the reproducibility of phantoms having a similar elasticity. Following mechanical stretching tests, the phantoms were used to compare the accuracy of four different elastography methods. The four methods were based on a one-dimensional (1-D) scaling factor estimation, on two different implementations of a 2-D Lagrangian speckle model estimator (quasistatic elastography methods), and on a 1-D shear wave transient elastography technique (dynamic method). Young's modulus was investigated as a function of the number of freeze-thaw cycles of PVA-C, and of the concentration of acoustic scatterers. Other mechanical and acoustic parameters-such as the speed of sound, shear wave velocity, mass density, and Poisson's ratio-also were assessed. The Poisson's ratio was estimated with good precision at 0.499 for all samples, and the Young's moduli varied in a range of 20 kPa for one freeze-thaw cycle to 600 kPa for 10 cycles. Nevertheless, above six freeze-thaw cycles, the results were less reliable because of sample geometry artifacts. However, for the samples that underwent less than seven freeze-thaw cycles, the Young's moduli estimated with the four elastography methods showed good matching with the mechanical tensile tests with a regression coefficient varying from 0.97 to 1.07, and correlations R2 varying from 0.93 to 0.99, depending on the method.     

Ultrafast compound imaging for 2-D motion vector estimation: application to transient elastography

This paper describes a new technique for two-dimensional (2-D) imaging of the motion vector at a very high frame rate with ultrasound. Its potential is experimentally demonstrated for transient elastography. But, beyond this application, it also could be promising for color flow and reflectivity imaging. To date, only axial displacements induced in human tissues by low-frequency vibrators were measured during transient elastography. The proposed technique allows us to follow both axial and lateral displacements during the shear wave propagation and thus should improve Young's modulus image reconstruction. The process is a combination of several ideas well-known in ultrasonic imaging: ultra-fast imaging, multisynthetic aperture beamforming, 1-D speckle tracking, and compound imaging. Classical beamforming in the transmit mode is replaced here by a single plane wave insonification increasing the frame rate by at least a factor of 128. The beamforming is achieved only in the receive mode on two independent subapertures. Comparison of successive frames by a classical 1-D speckle tracking algorithm allows estimation of displacements along two different directions linked to the subapertures beams. The variance of the estimates is finally improved by tilting the emitting plane wave at each insonification, thus allowing reception of successive decorrelated speckle patterns.     

Evolution of ultrasonic velocity and dynamic elastic moduli with shear strain in granular layers

Ultrasonic wave transmission has been used to investigate processes that influence frictional strength, strain localization, fabric development, porosity evolution, and friction constitutive properties in granular materials under a wide range of conditions. We present results from a novel technique using ultrasonic wave propagation to observe the evolution of elastic properties during shear in laboratory experiments conducted at stresses applicable to tectonic faults in Earth’s crust. Elastic properties were measured continuously during loading, compaction, and subsequent shear using piezoelectric transducers fixed within shear forcing blocks in the double-direct-shear configuration. We report high-fidelity measurements of elastic wave properties for normal stresses up to 20 MPa and shear strains up to 500 % in layers of granular quartz, smectite clay, and a quartz-clay mixture. Layers were 0.1–1 cm thick and had nominal contact area of $5 \mathrm{cm} \!\times \! 5 \mathrm{cm}$ . We investigate relationships among frictional strength, granular layer thickness, and ultrasonic wave velocity and amplitude as a function of shear strain and normal stress. For layers of granular quartz, P-wave velocity and amplitude decrease by 20–70 % after a shear strain of 0.5. We find that P-wave velocity increases upon application of shear load for layers of pure clay and for the quartz-clay mixture. The P-wave amplitude of pure clay and quart-clay mixtures first decreases by $\sim $ 50 and 30 %, respectively, and then increases with additional shear strain. Changes in P-wave speed and wave amplitude result from changes in grain contact stiffness, crack density and disruption of granular force chains. Our data indicate that sample dilation and shear localization influence acoustic velocity and amplitude during granular shear.     

平面P1波入射下饱和度对深埋圆形衬砌动应力集中因子的影响

准饱和土是天然土体中比较常见的一种形式,与饱和土相同,其内可传播三种体波:剪切波、快压缩波(P1波)和慢压缩波(P2波)。忽略土颗粒的压缩性,考虑孔隙流体的压缩性以及孔隙流体与土骨架之间的粘性耦合作用,采用修正的B iot模型,运用波函数展开法研究了平面P1波入射时,准饱和土体内深埋圆形衬砌的散射问题,通过数值计算模拟并分析了衬砌的动应力集中问题,结果表明:不论入射频率为低频(100Hz)还是高频(500Hz),对于相同的衬砌,当土体为准饱和时,动应力集中因子沿周向分布曲线的形状和数值基本保持不变;但土体为完全饱和时,虽然饱和度发生了微小的变化(比99%只提高了1%),但动应力集中因子的值却大幅度地降低,其沿周向分布曲线的形状也发生了显著的变化,在工程设计中应重视饱和度对圆形衬砌动应力集中因子的影响。     

Diffraction field of a low frequency vibrator in soft tissues using transient elastography

For the last 10 years, interest has grown in low frequency shear waves that propagate in the human body. However, the generation of shear waves by acoustic vibrators is a relatively complex problem, and the directivity patterns of shear waves produced by the usual vibrators are more complicated than those obtained for longitudinal ultrasonic transducers. To extract shear modulus parameters from the shear wave propagation in soft tissues, it is important to understand and to optimize the directivity pattern of shear wave vibrators. This paper is devoted to a careful study of the theoretical and the experimental directivity pattern produced by a point source in soft tissues. Both theoretical and experimental measurements show that the directivity pattern of a point source vibrator presents two very strong lobes for an angle around 35 degrees . This paper also points out the impact of the near field in the problem of shear wave generation.     

Time-frequency analyses of fluid–solid interaction under sinusoidal translational shear deformation of the viscoelastic rat cerebrum

During normal extracellular fluid (ECF) flow in the brain glymphatic system or during pathological flow induced by trauma resulting from impacts and blast waves, ECF–solid matter interactions result from sinusoidal shear waves in the brain and cranial arterial tissue, both heterogeneous biological tissues with high fluid content. The flow in the glymphatic system is known to be forced by pulsations of the cranial arteries at about 1 Hz. The experimental shear stress response to sinusoidal translational shear deformation at 1 Hz and 25% strain amplitude and either 0% or 33% compression is compared for rat cerebrum and bovine aortic tissue. Time-frequency analyses aim to correlate the shear stress signal frequency components over time with the behavior of brain tissue constituents to identify the physical source of the shear nonlinear viscoelastic response. Discrete fast Fourier transformation analysis and the novel application to the shear stress signal of harmonic wavelet decomposition both show significant 1 Hz and 3 Hz components. The 3 Hz component in brain tissue, whose magnitude is much larger than in aortic tissue, may result from interstitial fluid induced drag forces. The harmonic wavelet decomposition locates 3 Hz harmonics whose magnitudes decrease on subsequent cycles perhaps because of bond breaking that results in easier fluid movement. Both tissues exhibit transient shear stress softening similar to the Mullins effect in rubber. The form of a new mathematical model for the drag force produced by ECF–solid matter interactions captures the third harmonic seen experimentally.     

Measurement of viscoelastic properties of tissue-mimicking material using longitudinal wave excitation

This paper presents an experimental framework for the measurement of the viscoelastic properties of tissuemimicking material. The novelty of the presented framework is in the use of longitudinal wave excitation and the study of the longitudinal wave patterns in finite media for the measurement of the viscoelastic properties. Ultrasound is used to track the longitudinal motions inside a test block. The viscoelastic parameters of the block are then estimated by 2 methods: a wavelength measurement method and a model fitting method. Connections are also made with shear elastography. The viscoelastic parameters are estimated for several homogeneous phantom blocks. The results from the new methods are compared with the conventional rheometry results.     

Rubber toughening of plastics

Mechanisms of fatigue damage in ABS and HIPS polymers were studied by monitoring changes in mechanical properties under sinusoidal and square wave loading at a frequency of 0.033 Hz. Both materials exhibited a large increase in hysteresis and a decrease in modulus. In ABS, these changes occurred to approximately the same extent in tension and compression, whereas in HIPS the changes in tensile properties were much more pronounced. From the shape of the hysteresis loops, and from volumetric measurements, it was concluded that crazing is the dominent fatigue damage mechanism in HIPS, whilst shear yielding is responsible for most of the damage observed in ABS. Large increases in hysteresis caused substantial rises in temperature despite the low frequency, and thus accelerated fatigue damage accumulation.