External vibration multi-directional ultrasound shearwave elastography (EVMUSE): application in liver fibrosis staging

Continuous Shear Wave Elastography for Liver Using Frame-to-Frame Equalization of Complex Amplitude

This study addresses a crucial necessity in the field of noninvasive liver fibrosis diagnosis by introducing the concept of continuous shear wave elastography (C-SWE), utilizing an external vibration source and color Doppler imaging. However, an application of C-SWE to assess liver elasticity, a deep region within the human body, arises an issue of signal instability in the obtained data. To tackle this challenge, this work proposes a method involving the acquisition of multiple frames of datasets, which are subsequently compressed. Furthermore, the proposed frame-to-frame equalization method compensates discrepancies in the initial phase that might exist among multiple-frame datasets, thereby significantly enhancing signal stability. The experimental validation of this approach encompasses both phantom tests and in vivo experiments. In the phantom tests, the proposed technique is validated through a comparison with the established shear wave elastography (SWE) technique. The results demonstrate a remarkable agreement, with an error in shear wave velocity of less than 4.2%. Additionally, the efficacy of the proposed method is confirmed through in vivo tests. As a result, the stabilization of observed shear waves using the frame-to-frame equalization technique exhibits promising potential for accurately assessing human liver elasticity. These findings collectively underscore the viability of C-SWE as a potential diagnostic instrument for liver fibrosis.

Examination of rapid adjustment system based on screen score obtained using continuous shear wave elastography

PURPOSE

Continuous shear wave elastography (C-SWE) can be expected to be applied to portable muscle elasticity diagnosis. To establish diagnostic technology, it will be necessary to improve measurement techniques and quantitative measurement accuracy.

METHODS

In this study, we investigated two screen scores: the quality index (Q-index), which determines whether the intensity of a power Doppler image is appropriate, and the shear wave propagation direction index (SWDI), which determines the uniformity of shear wave propagation.

RESULTS

First, we performed numerical simulations with white noise and found that the coefficient of variation of shear wave velocity estimation was less than 5% when the normalized Q-index was greater than 0.27. Furthermore, regarding the SWDI, we clarified the relationship between the standard deviation in shear wave propagation direction and the SWDI. Next, the relationship between the Q-index and coefficient of variation of estimated shear wave velocity was evaluated through experiments using a tissue-mimicking phantom. The results showed that there was a negative correlation between the Q-index and the coefficient of variation, and the fluctuation of the propagation velocity could be inferred from the Q-index. Finally, we showed the results of applying the screen scores to muscle relaxation monitoring and confirmed its usefulness in clinical applications.

CONCLUSION

By applying the screen scores, we showed improved stability in speed estimation in C-SWE, and demonstrated the possibility of clinical applicability.

Feasibility of Phase Velocity Imaging Using Multi Frequency Oscillation-Shear Wave Elastography

OBJECTIVE

To assess viscoelasticity, a pathologically relevant biomarker, shear wave elastography (SWE) generally uses phase velocity (PV) dispersion relationship generated via pulsed acoustic radiation force (ARF) excitation pulse. In this study, a multi-frequency oscillation (MFO)- excitation pulse with higher weight to higher frequencies is proposed to generate PV images via the generation of motion with energy concentrated at the target frequencies in contrast to the broadband frequency motion generated in pulsed SWE (PSWE).

METHODS

The feasibility of MFO-SWE to generate PV images at 100 to 1000 Hz in steps of 100 Hz was investigated by imaging 6 and 70 kPa inclusions with 6.5 and 10.4 mm diameter and ex vivo bovine liver with and without the presence of an aberration layer and chicken muscle ex vivo, and 4T1 mouse breast tumor, in vivo with comparisons to PSWE.

RESULTS

MFO-SWE-derived CNR was statistically higher than PSWE for 6 kPa (both with and without aberration) and 70 kPa (with aberration) inclusions and derived SNR of the liver was statistically higher than PSWE at higher frequency (600-1000 Hz). Quantitatively, at 600-1000 Hz, MFO-SWE improved CNR of inclusions (without and with) aberration on an average by (8.2 and 156)% and of the tumor by 122%, respectively, and improved SNR of the liver (without and with) aberration by (20.2 and 51.5)% and of chicken muscle by 72%, respectively compared to the PSWE.

CONCLUSIONS AND SIGNIFICANCE

These results indicate the advantages of MFO-SWE to improve PV estimation at higher frequencies which could improve viscoelasticity quantification and feature delineation.

Shear wave generation from non-spherical bubble collapse in a tissue phantom

Elastography is a non-invasive technique to detect tissue anomalies via the local elastic modulus using shear waves. Commonly shear waves are produced via acoustic focusing or the use of mechanical external sources, shear waves may result also naturally from cavitation bubbles during medical intervention, for example from thermal ablation. Here, we measure the shear wave emitted from a well-controlled single laser-induced cavitation bubble oscillating near a rigid boundary. The bubbles are generated in a transparent tissue-mimicking hydrogel embedded with tracer particles. High-speed imaging of the tracer particles and the bubble shape allow quantifying the shear wave and relate it to the bubble dynamics. It is found that different stages of the bubble dynamics contribute to the shear wave generation and the mechanism of shear wave emission, its direction and the efficiency of energy converted into the shear wave depend crucially on the bubble to wall stand-off distance.

Imaging of Single Transducer-Harmonic Motion Imaging-Derived Displacements at Several Oscillation Frequencies Simultaneously

Mapping of mechanical properties, dependent on the frequency of motion, is relevant in diagnosis, monitoring treatment response, or intra-operative surgical resection planning. While shear wave speeds at different frequencies have been described elsewhere, the effect of frequency on the "on-axis" acoustic radiation force (ARF)-induced displacement has not been previously investigated. Instead of generating single transducer-harmonic motion imaging (ST-HMI)-derived peak-to-peak displacement (P2PD) image at a particular frequency, a novel multi-frequency excitation pulse is proposed to generate P2PD images at 100-1000 Hz simultaneously. The performance of the proposed excitation pulse is compared with the ARFI by imaging 16 different inclusions (Young's moduli of 6, 9, 36, 70 kPa and diameters of 1.6, 2.5, 6.5, and 10.4 mm) embedded in an 18 kPa background. Depending on inclusion size and stiffness, the maximum CNR and contrast were achieved at different frequencies and were always higher than ARFI. The frequency, at which maximum CNR and contrast were achieved, increased with stiffness for fixed inclusion's size and decreased with size for fixed stiffness. In vivo feasibility is tested by imaging a 4T1 breast cancer mouse tumor on Day 6, 12, and 19 post-injection of tumor cells. Similar to phantoms, the CNR of ST-HMI images was higher than ARFI and increased with frequency for the tumor on Day 6. Besides, P2PD at 100-1000 Hz indicated that the tumor became stiffer with respect to the neighboring non-cancerous tissue over time. These results indicate the importance of using a multi-frequency excitation pulse to simultaneously generate displacement at multiple frequencies to better delineate inclusions or tumors.

Three-Dimensional Shear Wave Elastography Using a 2D Row Column Addressing (RCA) Array

Objective. To develop a 3D shear wave elastography (SWE) technique using a 2D row column addressing (RCA) array, with either external vibration or acoustic radiation force (ARF) as the shear wave source. Impact Statement. The proposed method paves the way for clinical translation of 3D SWE based on the 2D RCA, providing a low-cost and high volume rate solution that is compatible with existing clinical systems. Introduction. SWE is an established ultrasound imaging modality that provides a direct and quantitative assessment of tissue stiffness, which is significant for a wide range of clinical applications including cancer and liver fibrosis. SWE requires high frame rate imaging for robust shear wave tracking. Due to the technical challenges associated with high volume rate imaging in 3D, current SWE techniques are typically confined to 2D. Advancing SWE from 2D to 3D is significant because of the heterogeneous nature of tissue, which demands 3D imaging for accurate and comprehensive evaluation. Methods. A 3D SWE method using a RCA array was developed with a volume rate up to 2000 Hz. The performance of the proposed method was systematically evaluated on tissue-mimicking elasticity phantoms and in an in vivo case study. Results. 3D shear wave motion induced by either external vibration or ARF was successfully detected with the proposed method. Robust 3D shear wave speed maps were reconstructed for phantoms and in vivo. Conclusion. The high volume rate 3D imaging provided by the 2D RCA array provides a robust and practical solution for 3D SWE with a clear pathway for future clinical translation.

A Dual-modal Imaging Method Combining Ultrasound and Electromagnetism for Simultaneous Measurement of Tissue Elasticity and Electrical Conductivity

The mechanical and electrical properties of soft tissues are relative to soft tissues' pathological state. Modern medical imaging devices have shown a trend to multi-modal imaging, which will provide complementary functional information to improve the accuracy of disease diagnosis. However, no method or system can simultaneously measure the mechanical and electrical properties of the soft tissue. In this study, we proposed a novel dual-modal imaging method integrated by shear wave elasticity imaging (SWEI) and Magneto-acousto-electrical tomography (MAET) to measure soft tissue's elasticity and conductivity simultaneously. A dual-modal imaging system based on a linear array transducer is built, and the imaging performances of MAET and SWEI were respectively evaluated by phantoms experiment and \textit{in vitro} experiment. Conductivity phantom experiments show that the MAET in this dual-modal system can image conductivity gradient as low as 0.4 S/m. The phantom experiments show that the reconstructed 2-D elasticity maps of the phantoms with inclusions with a diameter larger than 5 mm are relatively accurate. \textit{In vitro} experiments show that the elasticity parameter can significantly distinguish the changes in tissue before and after heating. This study first proposes a method that can simultaneously obtain tissue elasticity and electrical conductivity to the best of our knowledge. Although this paper just carried out the proof of concept experiments of the new method, it demonstrates great potential for disease diagnosis in the future.

Controllable Angle Shear Wavefront Reconstruction Based on Image Fusion Method for Shear Wave Elasticity Imaging

The generation and measurement of shear waves are critical in ultrasonic elasticity imaging. Generally, the resulting wavefront direction is very important for accurately measuring the shear speed and estimating the medium elasticity. In this article, the proposed method can generate a compound shear wavefront with the same direction as speed reconstruction and zero angles between the wavefront and the focus direction, which can improve the estimation accuracy of shear wave velocity. Also, this method, called time-division multipoint excitation image fusion (TDMPEIF), can reconstruct the shear wave propagation images acquired at different depths of a medium according to the frame sequence to produce the shear waves front with a regulable angle. Moreover, the shear wave speed and the elasticity of a medium can be mapped quantitatively with this method. The results demonstrate that the TDMPEIF can improve the quality of the shear wave velocity images, which has wide application value and good promotion prospects for quantitative evaluation of tissue elasticity.

Multiparametric Quantitative US Examination of Liver Fibrosis: A Feature-engineering and Machine-learning Based Analysis

Quantitative ultrasound (QUS), which is commonly used to extract quantitative features from the ultrasound radiofrequency (RF) data or the RF envelope signals for tissue characterization, is becoming a promising technique for noninvasive assessments of liver fibrosis. However, the number of feature variables examined and finally used in the existing QUS methods is typically small, to some extent limiting the diagnostic performance. Therefore, this paper devises a new multiparametric QUS (MP-QUS) method which enables the extraction of a large number of feature variables from US RF signals and allows for the use of feature-engineering and machinelearning based algorithms for liver fibrosis assessment. In the MP-QUS, eighty-four feature variables were extracted from multiple QUS parametric maps derived from the RF signals and the envelope data. Afterwards, feature reduction and selection were performed in turn to remove the feature redundancy and identify the best combination of features in the reduced feature set. Finally, a variety of machine-learning algorithms were tested for classifying liver fibrosis with the selected features, based on the results of which the optimal classifier was established and used for final classification. The performance of the proposed MPQUS method for staging liver fibrosis was evaluated on an animal model, with histologic examination as the reference standard. The mean accuracy, sensitivity, specificity and area under the receiver-operating-characteristic curve achieved by MP-QUS are respectively 83.38%, 86.04%, 80.82% and 0.891 for recognizing significant liver fibrosis, and 85.50%, 88.92%, 85.24% and 0.924 for diagnosing liver cirrhosis. The proposed MP-QUS method paves a way for its future extension to assess liver fibrosis in human subjects.

A Robust and Fast Method for 2-D Shear Wave Speed Calculation

We developed a new method, called the tangent plane method (TPM), for more efficiently and accurately estimating 2-D shear wave speed (SWS) from any direction of wave propagation. In this technique, we estimate SWS by solving the Eikonal equation because this approach is more robust to noise. To further enhance the performance, we computed the tangent plane of the arrival time surface. To evaluate the approach, we performed simulations and also conducted phantom studies. Simulation studies showed that TPM was more robust to noise than the conventional methods such as 2-D cross correlation (CC) and the distance method. The contrast/CNR for an inclusion (69 kPa; manufacturer provided stiffness) of a phantom is 0.54/4.17, 0.54/1.82, and 0.46/1.22. SWS results [mean and standard deviation (SD)] were 4.41 ± 0.49, 4.62 ± 0.85, and 3.66 ± 0.99 m/s, respectively, while the manufacturer's reported value (mean and range) is 4.81 ± 0.49 m/s. This shows that TPM has the higher CNR and lower SD than other methods. To increase the computation speed, an iterative version of TPM (ITPM) was also developed, which calculated the time-of-flight iteratively. ITPM reduced the computation time to 3.6%, i.e., from 748 to 27 s. In vivo case analysis demonstrated the feasibility of using the conventional ultrasound scanner for the proposed 2-D SWS algorithms.

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.

Quantitative Shear Wave Speed Assessment for Muscles with the Diagnosis of Taut Bands and/or Myofascial Trigger Points Using Probe Oscillation Shear Wave Elastography: A Pilot Study

OBJECTIVE

To use probe oscillation shear wave elastography (PROSE) with two vibration sources to generate two shear waves in the imaging plane to quantitatively assess the shear wave speeds (SWSs) of muscles with and without the diagnosis of taut bands (TB) and/or myofascial trigger points (MTrPs).

METHODS

Thirty-three patients were scanned with the PROSE technique. Shear waves were generated through continuous vibration of the ultrasound probe, while the shear wave motions were detected using the same probe. SWSs for the sides with and without TBs and/or MTrPs were computed and compared. The pressure pain thresholds (PPTs) were measured as an indicator of maximum pain tolerance of patients. The statistical differences between the SWSs with and without TBs and/or MTrPs with different PPT values were analyzed using the nonparametric Wilcoxon rank-sum test.

RESULTS

The mean SWSs for the sides with TBs and/or MTrPs are faster than that of the contralateral side without TBs and/or MTrPs. A significant difference was observed between mean SWSs with and without TBs and/or MTrPs without any information of PPT, with rank-sum test P <  .005. Additionally, with the information of PPT, a significant difference was observed between mean SWSs for the sides with and without TBs and/or MTrPs, for PPT values between 0 and 50 N/cm² (P <  .005), but for PPT values between 50 and 90 N/cm² , it was difficult to differentiate mean SWSs with and without TBs and/or MTrPs.

CONCLUSION

Our preliminary results show that SWSs measured from patients had a significant difference between the mean SWSs with and without TBs and/or MTrPs.

Two-dimensional (2D) dynamic vibration optical coherence elastography (DV-OCE) for evaluating mechanical properties: a potential application in tissue engineering

Mechanical properties in tissues are an important indicator because they are associated with disease states. One of the well-known excitation sources in optical coherence elastography (OCE) to determine mechanical properties is acoustic radiation force (ARF); however, a complicated focusing alignment cannot be avoided. Another excitation source is a piezoelectric (PZT) stack to obtain strain images via compression, which can affect the intrinsic mechanical properties of tissues in tissue engineering. In this study, we report a new technique called two-dimensional (2D) dynamic vibration OCE (DV-OCE) to evaluate 2D wave velocities without tedious focusing alignment procedures and is a non-contact method with respect to the samples. The three-dimensional (3D) Fourier transform was utilized to transfer the traveling waves (x, y, t) into 3D k-space (k_x, k_y, f). A spatial 2D wavenumber filter and multi-angle directional filter were employed to decompose the waves with omni-directional components into four individual traveling directions. The 2D local wave velocity algorithm was used to calculate a 2D wave velocity map. Six materials, two homogeneous phantoms with 10 mm thickness, two homogeneous phantoms with 2 mm thickness, one heterogeneous phantom with 2 mm diameter inclusion and an ex vivo porcine kidney, were examined in this study. In addition, the ARF-OCE was used to evaluate wave velocities for comparison. Numerical simulations were performed to validate the proposed 2D dynamic vibration OCE technique. We demonstrate that the experimental results were in a good agreement with the results from ARF-OCE (transient OCE) and numerical simulations. Our proposed 2D dynamic vibration OCE could potentially pave the way for mechanical evaluation in tissue engineering and for laboratory translation with easy-to-setup and contactless advantages.

Diffuse shear wave spectroscopy for soft tissue viscoelastic characterization

In order to limit and slow the development of diseases, they have to be diagnosed early as possible to treat patients in a better and more rapid manner. In this paper, we focus on a noninvasive and quick method based on diffuse fields in elastography to detect diseases that affect the stiffness of organs. To validate our method, a phantom experiment numerically pre-validated is designed. It consists of seven vibrators that generate white noises in a bandwidth of [80-300] Hz and then a complex acoustic field in a phantom. Waves are tracked by a linear ultrasound probe L11-4v linked to a Verasonics Vantage System and are converted into a particle velocity 2D map as a function of time. The phase velocity of the shear waves is calculated using a temporal and 2D spatial Fourier transform and an adapted signal-processing method. Shear wave velocity dispersion measurement in the frequency bandwidth of the vibrators enables one to characterize the dynamic hardness of the material through the viscoelastic parameters μ and η in an acquisition time shorter than a second (T_acq = 300 ms). With the aim of estimating the consistency of the method, the experiment is performed N = 10 times. The measured elastic modulus and viscous parameter that quantify the dynamic properties of the medium correspond to the expected values: μ = 1.23 ± 0.05 kPa and η = 0.51 ± 0.09 Pa∙s.

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.

Application of the novel estimation method by shear wave elastography using vibrator to human skeletal muscle

In recent years, non-invasive measurement of tissue stiffness (hardness) using ultrasound elastography has attracted considerable attention. It has been used to evaluate muscle stiffness in the fields of rehabilitation, sports, and orthopedics. However, ultrasonic diagnostic devices with elastography systems are expensive and clinical use of such devices has been limited. In this study, we proposed a novel estimation method for vibration-based shear wave elastography measurement of human skeletal muscle, then determined its reproducibility and reliability. The coefficient of variation and correlation coefficient were used to determine reproducibility and reliability of the method by measuring the shear wave velocities in konjac phantom gels and agar phantom gels, as well as skeletal muscle. The intra-day, day-to-day, and inter-operator reliabilities were good when measuring the shear wave velocities in phantom gels. The intra-day and day-to-day reliabilities were good when measuring the shear wave velocities in skeletal muscle. The findings confirmed adequate reproducibility and reliability of the novel estimation method for vibration-based shear wave elastography. Therefore, the proposed measurement method may be a useful tool for evaluation of muscle stiffness.

Three-dimensional shear wave elastography on conventional ultrasound scanners with external vibration

Two-dimensional (2D) ultrasound shear wave elastography (SWE) has been widely used for soft tissue properties assessment. Given that shear waves propagate in three dimensions (3D), extending SWE from 2D to 3D is important for comprehensive and accurate stiffness measurement. However, implementation of 3D SWE on a conventional ultrasound scanner is challenging due to the low volume rate (tens of Hertz) associated with limited parallel receive capability of the scanner's hardware beamformer. Therefore, we developed an external mechanical vibration-based 3D SWE technique allowing robust 3D shear wave tracking and speed reconstruction for conventional scanners. The aliased shear wave signal detected with a sub-Nyquist sampling frequency was corrected by leveraging the cyclic nature of the sinusoidal shear wave generated by the external vibrator. Shear wave signals from different sub-volumes were aligned in temporal direction to correct time delays from sequential pulse-echo events, followed by 3D speed reconstruction using a 3D local frequency estimation algorithm. The technique was validated on liver fibrosis phantoms with different stiffness, showing good correlation (r = 0.99, p < 0.001) with values measured from a state-of-the-art SWE system (GE LOGIQ E9). The phantoms with different stiffnesses can be well-differentiated regardless of the external vibrator position, indicating the feasibility of the 3D SWE with regard to different shear wave propagation scenarios. Finally, shear wave speed calculated by the 3D method correlated well with magnetic resonance elastography performed on human liver (r = 0.93, p = 0.02), demonstrating the in vivo feasibility. The proposed technique relies on low volume rate imaging and can be implemented on the widely available clinical ultrasound scanners, facilitating its clinical translation to improve liver fibrosis evaluation.

2-D Ultrasound Shear Wave Elastography With Multi-Sphere-Source External Mechanical Vibration: Preliminary Phantom Results

Ultrasound shear wave elastography (SWE) imaging is emerging as a quantitative and non-invasive tissue characterization modality. Shear wave generation using external mechanical vibration (EMV) has received extensive research interest over acoustic radiation force impulse (ARFI) because of its low cost and potential for portability. In this paper, we propose an EMV concept with multiple spherical sources that can be easily reconfigured in three configurations to induce unique shear wave propagation patterns. We introduce two design embodiments of this concept bench test design for proof of concept and a clinically deployable design. The latter is designed to incorporate size, ergonomics, portability and power consumption considerations and constraints. Experimental validation on elasticity phantoms using both EMV designs demonstrates shear wave generation and elasticity reconstruction comparable in performance to ElastQ, a commercial ARFI-based shear elastography technology from Philips. In addition, the local displacement amplitude induced by EMV is 10 times greater than that induced by ARFI at the same given depth. Finally, the multiple configurations of the presented EMV design would allow exploration of advanced elastography methods such as tissue anisotropic elasticity.

Ramp-Creep Ultrasound Viscoelastography for Measuring Viscoelastic Parameters of Materials

Several ultrasound-based methods have been developed to evaluate the viscoelastic properties of materials. The purpose of this study is to introduce a novel viscoelastography method based on ultrasound acoustic radiation force for measuring the parameters relevant to the viscoelastic properties of materials, named ramp-creep ultrasound viscoelastography (RC viscoelastography). RC viscoelastography uses two different ultrasound excitation modes to cause ramp and creep strain responses in the material. By combining and analyzing the information obtained from these two modes of excitation, the viscoelastic parameters of the material can be quantitatively evaluated. Finite element computer simulation demonstrated that RC viscoelastography can accurately evaluate the viscoelastic parameters of the material, including the relaxation and creep time constants as well as the ratio of viscous fluids to solids in the material, except for the region near the top surface of the material. The novelty of RC viscoelastography is that there is no need to know the magnitude of acoustic radiation force and induced stress in the material in order to evaluate the viscoelastic parameters. In the future, experiments are necessary to test the performance of RC viscoelastography in real biomaterials and biological tissues.

Elastographic Tomosynthesis From X-Ray Strain Imaging of Breast Cancer

Noncancerous breast tissue and cancerous breast tissue have different elastic properties. In particular, cancerous breast tumors are stiff when compared to the noncancerous surrounding tissue. This difference in elasticity can be used as a means for detection through the method of elastographic tomosynthesis by means of physical modulation. This paper deals with a method to visualize elasticity of soft tissues, particularly breast tissues, via x-ray tomosynthesis. X-ray tomosynthesis is now used to visualize breast tissues with better resolution than the conventional single-shot mammography. The advantage of X-ray tomosynthesis over X-ray CT is that fewer projections are needed than CT to perform the reconstruction, thus radiation exposure and cost are both reduced. Two phantoms were used for the testing of this method, a physical phantom and an in silico phantom. The standard root mean square error in the tomosynthesis for the physical phantom was 2.093 and the error in the in silico phantom was negligible. The elastographs were created through the use of displacement and strain graphing. A Gaussian Mixture Model with an expectation–maximization clustering algorithm was applied in three dimensions with an error of 16.667%. The results of this paper have been substantial when using phantom data. There are no equivalent comparisons yet in 3D x-ray elastographic tomosynthesis. Tomosynthesis with and without physical modulation in the 3D elastograph can identify feature groupings used for biopsy. The studies have potential to be applied to human test data used as a guide for biopsy to improve accuracy of diagnosis results. Further research on this topic could prove to yield new techniques for human patient diagnosis purposes.

Image reconstruction utilizing median filtering applied to elastography

BACKGROUND

The resources of ultrafast technology can be used to add another analysis to ultrasound imaging: assessment of tissue viscoelasticity. Ultrafast image formation can be utilized to find transitory shear waves propagating in soft tissue, which permits quantification of the mechanical properties of the tissue via elastography. This technique permits simple and noninvasive diagnosis and monitoring of disease.

METHODS

This article presents a method to estimate the viscoelastic properties and rigidity of structures using the ultrasound technique known as shear wave elasticity imaging (SWEI). The Verasonics Vantage 128 research platform and L11-4v transducer were used to acquire radio frequency signals from a model 049A elastography phantom (CIRS, USA), with subsequent processing and analysis in MATLAB.

RESULTS

The images and indexes obtained reflect the qualitative measurements of the different regions of inclusions in the phantom and the respective alterations in the viscoelastic properties of distinct areas. Comparison of the results obtained with this proposed technique and other commonly used techniques demonstrates the characteristics of median filtering in smoothing variations in velocity to form elastographic images. The results from the technique proposed in this study are within the margins of error indicated by the phantom manufacturer for each type of inclusion; for the phantom base and for type I, II, III, and IV inclusions, respectively, in kPa and percentage errors, these are 25 (24.0%), 8 (37.5%), 14 (28.6%), 45 (17.8%), and 80 (15.0%). The values obtained using the method proposed in this study and mean percentage errors were 29.18 (- 16.7%), 10.26 (- 28.2%), 15.64 (- 11.7%), 45.81 (- 1.8%), and 85.21 (- 6.5%), respectively.

CONCLUSIONS

The new technique to obtain images uses a distinct filtering function which considers the mean velocity in the region around each pixel, in turn allowing adjustments according to the characteristics of the phantom inclusions within the ultrasound and optimizing the resulting elastographic images.

In Vivo Open- and Closed-chest Measurements of Left-Ventricular Myocardial Viscoelasticity using Lamb wave Dispersion Ultrasound Vibrometry (LDUV): A Feasibility Study

Diastolic dysfunction causes close to half of congestive heart failures and is associated with increased stiffness in left-ventricular myocardium. A clinical tool capable of measuring viscoelasticity of the myocardium could be beneficial in clinical settings. We used Lamb wave Dispersion Ultrasound Vibrometry (LDUV) for assessing the feasibility of making in vivo non-invasive measurements of myocardial elasticity and viscosity in pigs. In vivo open-chest measurements of myocardial elasticity and viscosity obtained using a Fourier space based analysis of Lamb wave dispersion are reported. The approach was used to perform ECG-gated transthoracic in vivo measurements of group velocity, elasticity and viscosity throughout a single heart cycle. Group velocity, elasticity and viscosity in the frequency range 50-500 Hz increased from diastole to systole, consistent with contraction and relaxation of the myocardium. Systolic group velocity, elasticity and viscosity were 5.0 m/s, 19.1 kPa, 6.8 Pa·s, respectively. In diastole, the measured group velocity, elasticity and viscosity were 1.5 m/s, 5.1 kPa and 3.2 Pa·s, respectively.

Probe Oscillation Shear Wave Elastography: Initial In Vivo Results in Liver

Shear wave elastography methods are able to accurately measure tissue stiffness, allowing these techniques to monitor the progression of hepatic fibrosis. While many methods rely on acoustic radiation force to generate shear waves for 2-D imaging, probe oscillation shear wave elastography (PROSE) provides an alternative approach by generating shear waves through continuous vibration of the ultrasound probe while simultaneously detecting the resulting motion. The generated shear wave field in in vivo liver is complicated, and the amplitude and quality of these shear waves can be influenced by the placement of the vibrating probe. To address these challenges, a real-time shear wave visualization tool was implemented to provide instantaneous visual feedback to optimize probe placement. Even with the real-time display, it was not possible to fully suppress residual motion with established filtering methods. To solve this problem, the shear wave signal in each frame was decoupled from motion and other sources through the use of a parameter-free empirical mode decomposition before calculating shear wave speeds. This method was evaluated in a phantom as well as in in vivo livers from five volunteers. PROSE results in the phantom as well as in vivo liver correlated well with independent measurements using the commercial General Electric Logiq E9 scanner.

Characterization of Viscoelastic Materials Using Group Shear Wave Speeds

Recent investigations of viscoelastic properties of materials have been performed by observing shear wave propagation following localized, impulsive excitations, and Fourier decomposing the shear wave signal to parameterize the frequency-dependent phase velocity using a material model. This paper describes a new method to characterize viscoelastic materials using group shear wave speeds , , and determined from the shear wave displacement, velocity, and acceleration signals, respectively. Materials are modeled using a two-parameter linear attenuation model with phase velocity and dispersion slope at a reference frequency of 200 Hz. Analytically calculated lookup tables are used to determine the two material parameters from pairs of measured group shear wave speeds. Green's function calculations are used to validate the analytic model. Results are reported for measurements in viscoelastic and approximately elastic phantoms and demonstrate good agreement with phase velocities measured using Fourier analysis of the measured shear wave signals. The calculated lookup tables are relatively insensitive to the excitation configuration. While many commercial shear wave elasticity imaging systems report group shear wave speeds as the measures of material stiffness, this paper demonstrates that differences , , and of group speeds are first-order measures of the viscous properties of materials.

Shear Wave Speed Estimation Using Reverberant Shear Wave Fields: Implementation and Feasibility Studies

Elastography is a modality that estimates tissue stiffness and, thus, provides useful information for clinical diagnosis. Attention has focused on the measurement of shear wave propagation; however, many methods assume shear wave propagation is unidirectional and aligned with the lateral imaging direction. Any deviations from the assumed propagation result in biased estimates of shear wave speed. To address these challenges, directional filters have been applied to isolate shear waves with different propagation directions. Recently, a new method was proposed for tissue stiffness estimation involving creation of a reverberant shear wave field propagating in all directions within the medium. These reverberant conditions lead to simple solutions, facile implementation and rapid viscoelasticity estimation of local tissue. In this work, this new approach based on reverberant shear waves was evaluated and compared with another well-known elastography technique using two calibrated elastic and viscoelastic phantoms. Additionally, the clinical feasibility of this technique was analyzed by assessing shear wave speed in human liver and breast tissues, in vivo. The results indicate that it is possible to estimate the viscoelastic properties in each scanned medium. Moreover, a better approach to estimation of shear wave speed was obtained when only the phase information was taken from the reverberant waves, which is equivalent to setting all magnitudes within the bandpass equal to unity: an idealization of a perfectly isotropic reverberant shear wave field.

Two-Dimensional Time-Harmonic Elastography of the Human Liver and Spleen

Measurement of shear wave speed of the liver and spleen by elastography is an established diagnostic procedure for the detection of hepatic fibrosis, portal hypertension and esophageal varices. However, current elastography systems are limited by the size and penetration depth of elastographic windows. In this study, 2D time-harmonic elastography is proposed for generating full field-of-view shear wave speed maps in great depth. Two-dimensional time-harmonic elastography uses external harmonic stimulation at multiple frequencies to create compound shear wave speed maps. The method is tested in a phantom with soft and stiff inclusions and used for elastography of the liver and spleen in 13 asymptomatic volunteers. Each volunteer was scanned twice to determine the sensitivity of the method to physiologic variations: first, after 2 h of fasting, and a second time, 15 min after drinking 1 L of water. The wave speed maps of the phantom clearly identified the soft and stiff inclusions, yielding values that were consistent with those from magnetic resonance elastography. In vivo wave speed values were 1.49 ± 0.11 m/s for the liver and 2.03 ± 0.15 m/s for the spleen in a lower-frequency band centered at 40 Hz and 3.15 ± 0.30 m/s for the spleen in a higher-frequency band centered at 120 Hz. After water intake, wave speed values increased by 6% in the liver (p = 0.002) and decreased in the spleen by 4% (p = 0.021, low-frequency band) and 6% (p = 0.0002, high-frequency band), suggesting the high sensitivity of the method to altered blood flow and perfusion pressure. Two-dimensional time-harmonic elastography of the liver and spleen is a promising method for measuring tissue stiffness at different states of blood flow and perfusion in a large tissue window and at great penetration depth.

Probe Oscillation Shear Elastography (PROSE): A High Frame-Rate Method for Two-Dimensional Ultrasound Shear Wave Elastography

Ultrasound shear wave elastography (SWE) utilizes the propagation of induced shear waves to characterize the shear modulus of soft tissue. Many methods rely on an acoustic radiation force (ARF) "push beam" to generate shear waves. However, specialized hardware is required to generate the push beams, and the thermal stress that is placed upon the ultrasound system, transducer, and tissue by the push beams currently limits the frame-rate to about 1 Hz. These constraints have limited the implementation of ARF to high-end clinical systems. This paper presents Probe Oscillation Shear Elastography (PROSE) as an alternative method to measure tissue elasticity. PROSE generates shear waves using a harmonic mechanical vibration of an ultrasound transducer, while simultaneously detecting motion with the same transducer under pulse-echo mode. Motion of the transducer during detection produces a "strain-like" compression artifact that is coupled with the observed shear waves. A novel symmetric sampling scheme is proposed such that pulse-echo detection events are acquired when the ultrasound transducer returns to the same physical position, allowing the shear waves to be decoupled from the compression artifact. Full field-of-view (FOV) two-dimensional (2D) shear wave speed images were obtained by applying a local frequency estimation (LFE) technique, capable of generating a 2D map from a single frame of shear wave motion. The shear wave imaging frame rate of PROSE is comparable to the vibration frequency, which can be an order of magnitude higher than ARF based techniques. PROSE was able to produce smooth and accurate shear wave images from three homogeneous phantoms with different moduli, with an effective frame rate of 300 Hz. An inclusion phantom study showed that increased vibration frequencies improved the accuracy of inclusion imaging, and allowed targets as small as 6.5 mm to be resolved with good contrast (contrast-to-noise ratio ≥ 19 dB) between the target and background.

Time Harmonic Elastography Reveals Sensitivity of Liver Stiffness to Water Ingestion

The aim of the study was to test the sensitivity of liver stiffness (LS) measured by time harmonic elastography in large tissue windows to water uptake and post-prandial effects. Each subject gave written informed consent to participate in this institutional review board-approved prospective study. LS was measured by time harmonic elastography in 10 healthy volunteers pre- and post-prandially, as well as before, directly after and 2 h after drinking water. The LS-time function during water intake was measured in 14 scans over 3 h in five volunteers. LS increased by 10% (p = 0.0015) post-prandially and by 11% (p = 0.0024) after pure water ingestion, and decreased to normal values after 2 h. LS was lower after overnight fasting than after 2-h fasting (3%, p = 0.04). Over the time course, LS increased to post-water peak values 15 min after drinking 0.25 L water and remained unaffected by further ingestion of water. In conclusion, our study indicates that LS measured by time harmonic elastography represents an effective-medium property sensitive to physiologic changes in vascular load of the liver.

Performance of 2-Dimensional Ultrasound Shear Wave Elastography in Liver Fibrosis Detection Using Magnetic Resonance Elastography as the Reference Standard: A Pilot Study

OBJECTIVES

To investigate the correlation between 2-dimensional (2D) ultrasound shear wave elastography (SWE) and magnetic resonance elastography (MRE) in liver stiffness measurement and the diagnostic performance of 2D SWE for liver fibrosis when imaging from different intercostal spaces and using MRE as the reference standard.

METHODS

Two-dimensional SWE was performed on 47 patients. One patient was excluded from the study. Each of the remaining 46 patients underwent same-day MRE for clinical purposes. The study was compliant with the Health Insurance Portability and Accountability Act and approved by the Institutional Review Board. Informed consent was obtained from each patient. Two-dimensional SWE measurements were acquired from the ninth, eighth, and seventh intercostal spaces. The correlation with MRE was calculated at each intercostal space and multiple intercostal spaces combined. The performance of 2D SWE in diagnosing liver fibrosis was evaluated by receiver operating characteristic curve analysis using MRE as the standard.

RESULTS

The 47 patients who initially underwent 2D SWE included 22 female and 25 male patients (age range, 19-77 years). The highest correlation between 2D SWE and MRE was from the eighth and seventh intercostal spaces (r = 0.68-0.76). The ranges of the areas under the receiver operating characteristic curves for separating normal or inflamed livers from fibrotic livers using MRE as the clinical reference were 0.84 to 0.92 when using the eighth and seventh intercostal spaces individually and 0.89 to 0.90 when combined.

CONCLUSIONS

The results suggest that 2D SWE and MRE are well correlated when SWE is performed at the eighth and seventh intercostal spaces. The ninth intercostal space is less reliable for diagnosing fibrosis with 2D SWE. Combining measurements from multiple intercostal spaces does not significantly improve the performance of 2D SWE for detection of fibrosis.

A new method for shear wave speed estimation in shear wave elastography

Visualization of mechanical properties of tissue can aid in noninvasive pathology diagnosis. Shear wave elastography (SWE) measures the elastic properties of soft tissues by estimation of local shear wave propagation speed. In this paper, a new robust method for estimation of shear wave speed is introduced which has the potential for simplifying continuous filtering and real-time elasticity processing. Shear waves were generated by external mechanical excitation and imaged at a high frame rate. Three homogeneous phantoms of varying elastic moduli and one inclusion phantom were imaged. Waves propagating in separate directions were filtered and shear wave speed was estimated by inversion of the 1-D first-order wave equation. Final 2-D shear wave speed maps were constructed by weighted averaging of estimates from opposite traveling directions. Shear wave speed results for phantoms with gelatin concentrations of 5%, 7%, and 9% were 1.52 ± 0.10 m/s, 1.86 ± 0.10 m/s, and 2.37 ± 0.15 m/s, respectively, which were consistent with estimates computed from three other conventional methods, as well as compression tests done with a commercial texture analyzer. The method was shown to be able to reconstruct a 2-D speed map of an inclusion phantom with good image quality and variance comparable to conventional methods. Suggestions for further work are given.

Coded excitation plane wave imaging for shear wave motion detection

Plane wave imaging has greatly advanced the field of shear wave elastography thanks to its ultrafast imaging frame rate and the large field-of-view (FOV). However, plane wave imaging also has decreased penetration due to lack of transmit focusing, which makes it challenging to use plane waves for shear wave detection in deep tissues and in obese patients. This study investigated the feasibility of implementing coded excitation in plane wave imaging for shear wave detection, with the hypothesis that coded ultrasound signals can provide superior detection penetration and shear wave SNR compared with conventional ultrasound signals. Both phase encoding (Barker code) and frequency encoding (chirp code) methods were studied. A first phantom experiment showed an approximate penetration gain of 2 to 4 cm for the coded pulses. Two subsequent phantom studies showed that all coded pulses outperformed the conventional short imaging pulse by providing superior sensitivity to small motion and robustness to weak ultrasound signals. Finally, an in vivo liver case study on an obese subject (body mass index = 40) demonstrated the feasibility of using the proposed method for in vivo applications, and showed that all coded pulses could provide higher SNR shear wave signals than the conventional short pulse. These findings indicate that by using coded excitation shear wave detection, one can benefit from the ultrafast imaging frame rate and large FOV provided by plane wave imaging while preserving good penetration and shear wave signal quality, which is essential for obtaining robust shear elasticity measurements of tissue.

Multifrequency time-harmonic elastography for the measurement of liver viscoelasticity in large tissue windows

Elastography of the liver for the non-invasive diagnosis of hepatic fibrosis is an established method. However, investigations of obese patients or patients with ascites are often limited by small and superficial elastographic windows. Therefore, multifrequency time-harmonic elastography (THE) based on time-resolved A-line ultrasound has recently been developed for measuring liver viscoelasticity in wide soft tissue windows and at greater depths. In this study, THE was integrated into a clinical B-mode scanner connected to a dedicated actuator bed driven by superimposed vibrations of 30- to 60-Hz frequencies. The resulting shear waves in the liver were captured along multiple profiles 7 to 14 cm in width and automatically processed for reconstruction of mean efficient shear wave speed and shear wave dispersion slope. This new modality was tested in healthy volunteers and 22 patients with clinically proven cirrhosis. Patients could be separated from controls by higher shear wave speeds (3.11 ± 0.64 m/s, 2.14-4.81 m/s, vs. 1.74 ± 0.10 m/s, 1.60-1.91 m/s) without significant degradation of data by high body mass index or ascites. Furthermore, the wave speed dispersion slope was significantly (p = 0.0025) lower in controls (5.2 ± 1.8 m/s/kHz) than in patients (39.1 ± 32.2 m/s/kHz). In conclusion, THE is useful for the diagnosis of cirrhosis in large tissue windows.