Radiological Society of North America/Quantitative Imaging Biomarker Alliance Shear Wave Speed Bias Quantification in Elastic and Viscoelastic Phantoms

Mitigating high frame rate demands in shear wave elastography using radial basis function-based reconstruction: An experimental phantom study

BACKGROUND

Shear wave elastography (SWE) is a technique that quantifies tissue stiffness by assessing the speed of shear waves propagating after being excited by acoustic radiation force. SWE allows the quantification of elastic tissue properties and serves as an adjunct to conventional ultrasound techniques, aiding in tissue characterization. To capture this transient propagation of the shear wave, the ultrasound device must be able to reach very high frame rates.

METHODOLOGY

In this paper, our aim is to relax the high frame rate requirement for SWE imaging. To this end, we propose lower frame rate SWE imaging followed by employing a 2-dimensional (2D) radial basis functions (RBF)-based interpolation. More specifically, the process involves obtaining low frame rate data and then temporal upsampling to reach a synthetic high frame rate data by inserting the 'UpS-1' image frames with missing values between two successive image frames (UpS: Upsampling rate). Finally, we apply the proposed interpolation technique to reconstruct the missing values within the incomplete high frame rate data.

RESULTS AND CONCLUSION

The results obtained from employing the proposed model on two experimental datasets indicate that we can relax the frame rate requirement of SWE imaging by a factor of 4 while maintaining shear wave speed (SWS), group velocity, and phase velocity estimates closely align with the high frame rate SWE model so that the error is less than 3%. Furthermore, analysis of the structural similarity index (SSIM) and root mean squared error (RMSE) on the 2D-SWS maps highlights the efficacy of the suggested technique in enhancing local SWS estimates, even at a downsampling (DS) factor of 4. For DS≤4, the SSIM values between the 2D-SWS maps produced by the proposed technique and those generated by the original high frame rate data consistently remain above 0.94. Additionally, the RMSE values is below 0.37 m/s, indicating promising performance of the proposed technique in reconstruction of SWS values.

Assessing measurement consistency of a diffusion tensor imaging (DTI) quality control (QC) anisotropy phantom

OBJECTIVES

We evaluated a quality control (QC) phantom designed to mimic diffusion characteristics and white matter fiber tracts in the brain. We hypothesized that acquisition of diffusion tensor imaging (DTI) data on different vendors and over multiple repeated measures would not contribute to significant variability in calculated diffusion tensor scalar metrics such as fractional anisotropy (FA) and mean diffusivity (MD).

MATERIALS AND METHODS

The DTI QC phantom was scanned using a 32-direction DTI sequence on General Electric (GE), Siemens, and Philips 3 Tesla scanners. Motion probing gradients (MPGs) were investigated as a source of variance in our statistical design, and data were acquired on GE and Siemens scanners using GE, Siemens, and Philips vendor MPGs for 32 directions. In total, 8 repeated scans were made for each GE/Siemens combination of vendor and MPGs with 8 repeated scans on a Philips machine using its stock DTI sequence. Data were analyzed using 2-way ANOVAs to investigate repeat scan and vendor variances and 3-way ANOVAs with repeat, MPG, and vendor as factors.

RESULTS

No statistical differences (i.e., P > 0.05) were found in any DTI scalar metrics (FA, MD) or for any factor, suggesting system constancy across imaging platforms and the specified phantom's reliability and reproducibility across vendors and conditions.

DISCUSSION

A DTI QC phantom demonstrates that DTI measurements maintain their consistency across different MRI systems and can contribute to a standard that is more reliable for quantitative MRI analyses.

Rapid wideband characterization of viscoelastic material properties by Bessel function-based time harmonic ultrasound elastography (B-THE)

Elastography is an emerging diagnostic technique that uses conventional imaging modalities such as sonography or magnetic resonance imaging to quantify tissue stiffness. However, different elastography methods provide different stiffness values, which require calibration using well-characterized phantoms or tissue samples. A comprehensive, fast, and cost-effective elastography technique for phantoms or tissue samples is still lacking. Therefore, we propose ultrasound Bessel-fit-based time harmonic elastography (B-THE) as a novel tool to provide rapid feedback on stiffness-related shear wave speed (SWS) and viscosity-related wave penetration rate (PR) over a wide range of harmonic vibration frequencies. The method relies on external induction and B-mode capture of cylindrical shear waves that satisfy the Bessel wave equation for efficient fit-based parameter recovery. B-THE was demonstrated in polyacrylamide phantoms in the frequency range of 20-200 Hz and was cross-validated by magnetic resonance elastography (MRE) using clinical 3-T MRI and compact 0.5-T tabletop MRI scanners. Frequency-independent material parameters were derived from rheological models and validated by numerical simulations. B-THE quantified frequency-resolved SWS and PR 13 to 176 times faster than more expensive clinical MRE and tabletop MRE and have a good accuracy (relative deviation to reference: 6 %, 10 % and 4 % respectively). Simulations of liver-mimicking material phantoms showed that a simultaneous fit of SWS and PR based on the fractional Maxwell rheological model outperformed a fit on PR solely. B-THE provides a comprehensive and fast elastography technique for the quantitative characterization of the viscoelastic behavior of soft tissue mimicking materials.

Advancements of non‐invasive imaging technologies for the diagnosis and staging of liver fibrosis: Present and future

Liver fibrosis is a reparative response triggered by liver injury. Non‐invasive assessment and staging of liver fibrosis in patients with chronic liver disease are of paramount importance, as treatment strategies and prognoses depend significantly on the degree of fibrosis. Although liver fibrosis has traditionally been staged through invasive liver biopsy, this method is prone to sampling errors, particularly when biopsy sizes are inadequate. Consequently, there is an urgent clinical need for an alternative to biopsy, one that ensures precise, sensitive, and non‐invasive diagnosis and staging of liver fibrosis. Non‐invasive imaging assessments have assumed a pivotal role in clinical practice, enjoying growing popularity and acceptance due to their potential for diagnosing, staging, and monitoring liver fibrosis. In this comprehensive review, we first delved into the current landscape of non‐invasive imaging technologies, assessing their accuracy and the transformative impact they have had on the diagnosis and management of liver fibrosis in both clinical practice and animal models. Additionally, we provided an in‐depth exploration of recent advancements in ultrasound imaging, computed tomography imaging, magnetic resonance imaging, nuclear medicine imaging, radiomics, and artificial intelligence within the field of liver fibrosis research. We summarized the key concepts, advantages, limitations, and diagnostic performance of each technique. Finally, we discussed the challenges associated with clinical implementation and offer our perspective on advancing the field, hoping to provide alternative directions for the future research.

Quality Assurance of Point and 2D Shear Wave Elastography through the Establishment of Baseline Data Using Phantoms

Ultrasound elastography has been available on most modern systems; however, the implementation of quality processes tends to be ad hoc. It is essential for a medical physicist to benchmark elastography measurements on each system and track them over time, especially after major software upgrades or repairs. This study aims to establish baseline data using phantoms and monitor them for quality assurance in elastography. In this paper, we utilized two phantoms: a set of cylinders, each with a composite material with varying Young's moduli, and an anthropomorphic abdominal phantom containing a liver modeled to represent early-stage fibrosis. These phantoms were imaged using three ultrasound manufacturers' elastography functions with either point or 2D elastography. The abdominal phantom was also imaged using magnetic resonance elastography (MRE) as it is recognized as the non-invasive gold standard for staging liver fibrosis. The scaling factor was determined based on the data acquired using MR and US elastography from the same vendor. The ultrasound elastography measurements showed inconsistency between different manufacturers, but within the same manufacturer, the measurements showed high repeatability. In conclusion, we have established baseline data for quality assurance procedures and specified the criteria for the acceptable range in liver fibrosis phantoms during routine testing.

Comparison of 2D Shear Wave Elastography and Transient Elastography in Non-Invasive Evaluation of Liver Fibrosis in Hepatitis C Virus-Related Chronic Liver Disease

Background: Transient Elastography (TE) is widely regarded as the most reliable non-invasive method for evaluating liver fibrosis. Recently, new techniques such as 2D Shear Wave Elastography (2D-SWE) have been developed. This study aimed to evaluate the correlation between TE and 2D-SWE in patients with HCV-related chronic liver disease and to redefine the cut-off values of 2D-SWE for predicting different stages of fibrosis based on our results. Methods: Both TE (Fibroscan, Echosens, Paris, France) and 2D-SWE (SuperSonic Imagine) were performed simultaneously in 170 patients, including those with active and eradicated HCV infection. Spearman's rank correlation coefficient was used to assess the correlation between the two measurements, and the concordance between the assigned METAVIR classes was calculated using Cohen's kappa coefficient. ROC curves were constructed to determine the optimal cut-off values for 2D-SWE. Results: Ten patients were excluded for invalid measurements. In the remaining 160 patients, TE and 2D-SWE demonstrated a high correlation (ρ = 0.83, p < 0.0001) and good agreement in METAVIR classification (k = 0.74). The optimal cut-off values identified for 2D-SWE were as follows: ≥ 7 kPa for F ≥ 2, ≥ 8.3 kPa for F ≥ 3, and ≥ 9.4 kPa for F4. Conclusions: 2D-SWE is a viable alternative to TE for patients with HCV-related chronic liver disease. Our data suggest that the currently accepted 2D-SWE cut-off values for cirrhosis (F4) should be reconsidered and potentially lowered.

Clinical, Safety, and Engineering Perspectives on Wearable Ultrasound Technology: A Review

Wearable ultrasound has the potential to become a disruptive technology enabling new applications not only in traditional clinical settings, but also in settings where ultrasound is not currently used. Understanding the basic engineering principles and limitations of wearable ultrasound is critical for clinicians, scientists, and engineers to advance potential applications and translate the technology from bench to bedside. Wearable ultrasound devices, especially monitoring devices, have the potential to apply acoustic energy to the body for far longer durations than conventional diagnostic ultrasound systems. Thus, bioeffects associated with prolonged acoustic exposure as well as skin health need to be carefully considered for wearable ultrasound devices. This article reviews emerging clinical applications, safety considerations, and future engineering and clinical research directions for wearable ultrasound technology.

Performance evaluation of commercial and non-commercial shear wave elastography implementations for vascular applications

BACKGROUND

Shear wave elastography (SWE) is mainly used for stiffness estimation of large, homogeneous tissues, such as the liver and breasts. However, little is known about its accuracy and applicability in thin (∼0.5-2 mm) vessel walls. To identify possible performance differences among vendors, we quantified differences in measured wave velocities obtained by commercial SWE implementations of various vendors over different imaging depths in a vessel-mimicking phantom. For reference, we measured SWE values in the cylindrical inclusions and homogeneous background of a commercial SWE phantom. Additionally, we compared the accuracy between a research implementation and the commercially available clinical SWE on an Aixplorer ultrasound system in phantoms and in vivo in patients.

METHODS

SWE measurements were performed over varying depths (0-35 mm) using three ultrasound machines with four ultrasound probes in the homogeneous 20 kPa background and cylindrical targets of 10, 40, and 60 kPa of a multi-purpose phantom (CIRS-040GSE) and in the anterior and posterior wall of a homogeneous polyvinyl alcohol vessel-mimicking phantom. These phantom data, along with in vivo SWE data of carotid arteries in 23 patients with a (prior) head and neck neoplasm, were also acquired in the research and clinical mode of the Aixplorer ultrasound machine. Machine-specific estimated phantom stiffness values (CIRS phantom) or wave velocities (vessel phantom) over all depths were visualized, and the relative error to the reference values and inter-frame variability (interquartile range/median) were calculated. Correlations between SWE values and target/vessel wall depth were explored in phantoms and in vivo using Spearman's correlations. Differences in wave velocities between the anterior and posterior arterial wall were assessed with Wilcoxon signed-rank tests. Intra-class correlation coefficients were calculated for a sample of ten patients as a measure of intra- and interobserver reproducibility of SWE analyses in research and clinical mode.

RESULTS

There was a high variability in obtained SWE values among ultrasound machines, probes, and, in some cases, with depth. Compared to the homogeneous CIRS-background, this variation was more pronounced for the inclusions and the vessel-mimicking phantom. Furthermore, higher stiffnesses were generally underestimated. In the vessel-mimicking phantom, anterior wave velocities were (incorrectly) higher than posterior wave velocities (3.4-5.6 m/s versus 2.9-5.9 m/s, p ≤ 0.005 for 3/4 probes) and remarkably correlated with measurement depth for most machines (Spearman's ρ = -0.873-0.969, p < 0.001 for 3/4 probes). In the Aixplorer's research mode, this difference was smaller (3.3-3.9 m/s versus 3.2-3.6 m/s, p = 0.005) and values did not correlate with measurement depth (Spearman's ρ = 0.039-0.659, p ≥ 0.002). In vivo, wave velocities were higher in the posterior than the anterior vessel wall in research (left p = 0.001, right p < 0.001) but not in clinical mode (left: p = 0.114, right: p = 0.483). Yet, wave velocities correlated with vessel wall depth in clinical (Spearman's ρ = 0.574-0.698, p < 0.001) but not in research mode (Spearman's ρ = -0.080-0.466, p ≥ 0.003).

CONCLUSIONS

We observed more variation in SWE values among ultrasound machines and probes in tissue with high stiffness and thin-walled geometry than in low stiffness, homogeneous tissue. Together with a depth-correlation in some machines, where carotid arteries have a fixed location, this calls for caution in interpreting SWE results in clinical practice for vascular applications.

Ultrasound Viscoelastography by Acoustic Radiation Force: A State-of-the-Art Review

Ultrasound elastography (USE) is a promising tool for tissue characterization as several diseases result in alterations of tissue structure and composition, which manifest as changes in tissue mechanical properties. By imaging the tissue response to an applied mechanical excitation, USE mimics the manual palpation performed by clinicians to sense the tissue elasticity for diagnostic purposes. Next to elasticity, viscosity has recently been investigated as an additional, relevant, diagnostic biomarker. Moreover, since biological tissues are inherently viscoelastic, accounting for viscosity in the tissue characterization process enhances the accuracy of the elasticity estimation. Recently, methods exploiting different acquisition and processing techniques have been proposed to perform ultrasound viscoelastography. After introducing the physics describing viscoelasticity, a comprehensive overview of the currently available USE acquisition techniques is provided, followed by a structured review of the existing viscoelasticity estimators classified according to the employed processing technique. These estimators are further reviewed from a clinical usage perspective, and current outstanding challenges are discussed.

Review of the Low-Temperature Acoustic Properties of Water, Aqueous Solutions, Lipids, and Soft Biological Tissues

Existing data on the acoustic properties of low-temperature biological materials is limited and widely dispersed across fields. This makes it difficult to employ this information in the development of ultrasound applications in the medical field, such as cryosurgery and rewarming of cryopreserved tissues. In this review, the low-temperature acoustic properties of biological materials, and the measurement methods used to acquire them were collected from a range of scientific fields. The measurements were reviewed from the acoustic setup to thermal methodologies for samples preparation, temperature monitoring, and system insulation. The collected data contain the longitudinal and shear velocity, and attenuation coefficient of biological soft tissues and biologically relevant substances-water, aqueous solutions, and lipids-in the temperature range down to -50 °C and in the frequency range from 108 kHz to 25 MHz. The multiple reflection method (MRM) was found to be the preferred method for low-temperature samples, with a buffer rod inserted between the transducer and sample to avoid direct contact. Longitudinal velocity changes are observed through the phase transition zone, which is sharp in pure water, and occurs more slowly and at lower temperatures with added solutes. Lipids show longer transition zones with smaller sound velocity changes; with the longitudinal velocity changes observed during phase transition in tissues lying between these two extremes. More general conclusions on the shear velocity and attenuation coefficient at low-temperatures are restricted by the limited data. This review enhance knowledge guiding for further development of ultrasound applications in low-temperature biomedical fields, and may help to increase the precision and standardization of low-temperature acoustic property measurements.

Systematic quantification of differences in shear wave elastography estimates between linear-elastic and viscoelastic material assumptionsa)

Quantitative, accurate, and standardized metrics are important for reliable shear wave elastography (SWE)-based biomarkers. For over two decades, the linear-elastic material assumption has been employed in SWE modes. In recent years, viscoelasticity estimation methods have been adopted in a few clinical systems. The current study aims to systematically quantify differences in SWE estimates obtained using linear-elastic and viscoelastic material assumptions. An acousto-mechanical simulation framework of acoustic radiation force impulse-based SWE was created to elucidate the effect of material viscosity and shear modulus on SWE estimates. Shear modulus estimates exhibited errors up to 72% when a numerical viscoelastic phantom was assessed as linearly elastic. Shear modulus estimates of polyvinyl alcohol phantoms between rheometry and SWE following the Kelvin-Voigt viscoelastic model assumptions were not significantly different. However, the percentage difference in shear modulus estimates between rheometry and SWE using the linear-elastic assumption was 50.1%-62.1%. In ex vivo liver, the percentage difference in shear modulus estimates between linear-elastic and viscoelastic methods was 76.1%. These findings provide a direct and systematic quantification of the potential error introduced when viscoelastic tissues are imaged with SWE following the linear-elastic assumption. This work emphasizes the need to utilize viscoelasticity estimation methods for developing robust quantitative imaging biomarkers.

Ultrasound Shear Wave Elastography in Cardiology

The advent of high-frame rate imaging in ultrasound allowed the development of shear wave elastography as a noninvasive alternative for myocardial stiffness assessment. It measures mechanical waves propagating along the cardiac wall with speeds that are related to stiffness. The use of cardiac shear wave elastography in clinical studies is increasing, but a proper understanding of the different factors that affect wave propagation is required to correctly interpret results because of the heart's thin-walled geometry and intricate material properties. The aims of this review are to give an overview of the general concepts in cardiac shear wave elastography and to discuss in depth the effects of age, hemodynamic loading, cardiac morphology, fiber architecture, contractility, viscoelasticity, and system-dependent factors on the measurements, with a focus on clinical application. It also describes how these factors should be considered during acquisition, analysis, and reporting to ensure an accurate, robust, and reproducible measurement of the shear wave.

Viscoelastic Characterization of Phantoms for Ultrasound Elastography Created Using Low- and High-Viscosity Poly(vinyl alcohol) with Ethylene Glycol as the Cryoprotectant

Ultrasound elastography enables noninvasive characterization of the tissue mechanical properties. Phantoms are widely used in ultrasound elastography for developing, testing, and validating imaging techniques. Creating phantoms with a range of viscoelastic properties relevant to human organs and pathological conditions remains an active area of research. Poly(vinyl alcohol) (PVA) cryogel phantoms offer a long shelf life, robustness, and convenient handling and storage. The goal of this study was to develop tunable phantoms using PVA with a clinically relevant range of viscoelastic properties. We combined low- and high-viscosity PVA to tune the viscoelastic properties of the phantom. Further, phantoms were created with an ethylene glycol-based cryoprotectant to determine whether it reduces the variability in the viscoelastic properties. Scanning electron microscopy (SEM) was performed to evaluate the differences in microstructure between phantoms. The density, longitudinal sound speed, and acoustic attenuation spectra (5-20 MHz) of the phantoms were measured. The phantoms were characterized using a shear wave viscoelastography approach assuming the Kelvin-Voigt model. Microstructural differences were revealed by SEM between phantoms with and without a cryoprotectant and with different PVA mixtures. The longitudinal sound speed and attenuation power-law fit exponent of the phantoms were within the clinical range (1510-1571 m/s and 1.23-1.38, respectively). The measured shear modulus (G) ranged from 3.3 to 17.7 kPa, and the viscosity (η) ranged from 2.6 to 7.3 Pa·s. The phantoms with the cryoprotectant were more homogeneous and had lower shear modulus and viscosity (G = 2.17 ± 0.2 kPa; η = 2.0 ± 0.05 Pa·s) than those without a cryoprotectant (G = 3.93 ± 0.7 kPa; η = 2.6 ± 0.14 Pa·s). Notably, phantoms with relatively constant viscosities and varying shear moduli were achieved by this method. These findings advance the development of well-characterized viscoelastic phantoms for use in elastography.

On the Role of Coherent Plane Wave Compounding in Shear Wave Elasticity Imaging: The Convolution Effect and Its Implications

OBJECTIVE

The clinical applicability of shear wave elasticity imaging (SWEI) has been confounded by its appreciable inter-system variability and unsatisfactory sensitivity. While SWEI relies on plane wave imaging (PWI) to achieve real-time rendering, it has been rarely noticed that PWI can affect SWEI's performance. This work is aimed at demonstrating that the use of coherent plane wave compounding (CPWC) can be a factor causing SWEI's underperformance.

METHODS

We presented a model to formally describe the slow-time behavior of CPWC in motion tracking. This model reveals that CPWC introduces temporal convolution on the observed motion, making the motion sampling process a low-pass filter (LPF). For validation, shear waves were produced in a phantom in the same way but sampled via PWI using different compounding numbers (CN) and pulse repetition frequencies (PRF), with the obtained signals compared with the inferences drawn from our model. Similar experiments were performed to reconstruct two small targets in the phantom in order to appraise the impact of CPWC on SWEI's sensitivity.

DISCUSSION

The validation experiment shows that the measurements match well with the model inferences, which verifies the LPF nature of CPWC. The phantom study also shows that either increasing CN or decreasing PRF can cause the loss of high-frequency motion information, leading to blurred target delineation by SWEI.

CONCLUSION

The convolution effect can help understand the variability of SWEI. Researchers should beware this effect when working on SWEI standardization. Clinicians using SWEI should also be cautious because this effect makes it harder to identify small lesions.

Influence of Measurement Depth and Acquisition Parameters on Shear Wave Speed and Shear Wave Dispersion in Certified Phantoms Using a Canon Aplio Clinical Ultrasound Scanner

OBJECTIVE

The aim of the work described here was to investigate the relative contribution of confounding factors on liver shear wave speed (SWS) and shear wave dispersion slope (SWDS) measurements in three certified phantoms using a Canon Aplio clinical ultrasound scanner.

METHODS

A Canon Aplio i800 i-series ultrasound system (Canon Medical Systems Corporation, Otawara, Tochigi, Japan) with i8CX1 convex array (center frequency = 4 MHz) was used to examine dependencies caused by the depth, width and height of the acquisition box (AQB), the depth and size of the region of interest (ROI), the AQB angle and the pressure of the ultrasound probe on the surface of the phantom.

RESULTS

Results revealed that depth is the most significant confounder in both SWS and SWDS measurements. AQB angle, height and width and ROI size exhibited minimal confounding effects on measurements. For SWS, the most consistent measurement depth is when the top of the AQB is placed between 2 and 4 cm, and the ROI is located between 3 and 7 cm deep. For SWDS, results indicate that measurement values significantly decrease with depth from the surface of the phantom until approximately 7 cm deep, and consequently no stable area of AQB placement or ROI depth exists.

CONCLUSION

In contrast to SWS, the same ideal acquisition depth range cannot necessarily be applied to SWDS measurements because of a significant depth dependency.

Evaluating Variability of Commercial Liver Fibrosis Elastography Phantoms

OBJECTIVE

Liver fibrosis has been found to increase the mechanical stiffness of the liver. To mimic different stages of liver fibrosis, commercially available phantoms (Model 039, CIRS, Inc.) have been produced for clinical quality assurance and research purposes. The purpose of this study was to investigate the mechanical property variability of the phantoms in two lots of CIRS Model 039 phantoms.

METHODS

Each lot consisted of phantoms of four stiffness types, and there were 8-10 phantoms of each type. Shear wave elastography measurements were conducted on each phantom at 10 different angles. Group velocity measurements and phase velocity curves were calculated for every SWE acquisition. Multilevel functional principal component analysis (MFPCA) was performed on phase velocity data, which decomposes each phase velocity curve into the sum of eigenfunctions of two levels. The variance of the component scores of levels 1 and 2 were used to represent inter-phantom and intra-phantom variability, respectively. The 95% confidence intervals of phase velocity in a phantom type were calculated to reflect curve variability.

DISCUSSION

The standard deviations of the group velocity for phantoms of any type were less than 0.04 and 0.02 m/s for lots 1 and 2, respectively. For both lots, in every type, the phase velocity curves of most individual phantoms fall within the 95% confidence interval.

CONCLUSION

MFPCA is an effective tool for analyzing the inter- and intra-phantom variability of phase velocity curves. Given the known variability of a fully tested lot, estimation of the variability of a new lot can be performed with a reduced number of phantoms tested.

A New Estimation Scheme for Improving the Performance of Shear Wave Elasticity Imaging

Shear wave velocity (SWV) reconstruction based on time-of-flight (TOF) is widely adopted to realize shear wave elasticity imaging (SWEI). It typically breaks down the reconstruction of a SWV image into many kernels and treats them independently. We hypothesized that information exchange among kernels improves the performance of SWEI. Therefore, we propose the approach of iterative re-weighted least squares based on inter-kernel communication (IKC-IRLS). We also hypothesized that time-to-peak (TTP) is superior to cross-correlation (CC) in visualizing small targets because TTP uses higher shear wave frequencies than CC. To examine the hypotheses, IKC-IRLS was combined with TTP data and compared with four established methods. The five methods were tested by imaging several small-size stiff targets (2.5, 4.0 and 6.4 mm in diameter) using different kernel sizes in the simulation and real experiments. The results indicate that the IKC-IRLS approach can mitigate speckle noise and is robust to TTP outliers. Consequently, the proposed method achieves the highest contrast-to-noise ratio and the lowest mean absolute percentage error of target in almost all tested cases.

Noninvasive ultrasound assessment of tissue internal pressure using dual mode elasticity imaging: a phantom study

Objective.Tissue internal pressure, such as interstitial fluid pressure in solid tumors and intramuscular pressure in compartment syndrome, is closely related to the pathological state of tissues. It is of great diagnostic value to measure and/or monitor the internal pressure of targeted tissues. Because most of the current methods for measuring tissue pressure are invasive, noninvasive methods are highly desired. In this study, we developed a noninvasive method for qualitative assessment of tissue internal pressure based on a combination of two ultrasound elasticity imaging methods: strain imaging and shear wave elasticity imaging.Approach.The method was verified through experimental investigation using two tissue-mimicking phantoms each having an inclusion confined by a membrane, in which hydrostatic pressures can be applied and maintained. To examine the sensitivity of the elasticity imaging methods to pressure variation, strain ratio and shear modulus ratio (SMR) between the inclusion and background of phantom were obtained.Main results.The results first experimentally prove that pressure, in addition to elasticity, is a contrast mechanism of strain imaging, and further demonstrate that a comparative analysis of strain ratio and SMR is an effective method for noninvasive tissue internal pressure detection.Significance.This work provides a new perspective in interpreting the strain ratio data in medical diagnosis, and it also provides a noninvasive alternative for assessing tissue internal pressure, which could be valuable for the diagnosis of pressure-related diseases.

US Backscatter for Liver Fat Quantification: An AIUM-RSNA QIBA Pulse-Echo Quantitative Ultrasound Initiative

Nonalcoholic fatty liver disease (NAFLD) is believed to affect one-third of American adults. Noninvasive methods that enable detection and monitoring of NAFLD have the potential for great public health benefits. Because of its low cost, portability, and noninvasiveness, US is an attractive alternative to both biopsy and MRI in the assessment of liver steatosis. NAFLD is qualitatively associated with enhanced B-mode US echogenicity, but visual measures of B-mode echogenicity are negatively affected by interobserver variability. Alternatively, quantitative backscatter parameters, including the hepatorenal index and backscatter coefficient, are being investigated with the goal of improving US-based characterization of NAFLD. The American Institute of Ultrasound in Medicine and Radiological Society of North America Quantitative Imaging Biomarkers Alliance are working to standardize US acquisition protocols and data analysis methods to improve the diagnostic performance of the backscatter coefficient in liver fat assessment. This review article explains the science and clinical evidence underlying backscatter for liver fat assessment. Recommendations for data collection are discussed, with the aim of minimizing potential confounding effects associated with technical and biologic variables.

Inter-platform Variability of Liver Elastography: Pairwise Comparisons of Four Devices

This study was aimed at determining whether liver stiffness measurements by 2-D shear wave elastography using GE's (2D-SWE-GE) and Canon's (2D-SWE-Canon) newest apparatus and vibration-controlled transient elastography (VCTE) share the same distribution of values compared with Hologic Supersonic Imagine (2D-SWE-SSI). In participants with chronic liver disease recruited in two university centers from August 2020 to February 2021, liver stiffness was measured the same day by the same operator with 2D-SWE-SSI plus one of the following devices: 2D-SWE-GE (n = 314), 2D-SWE-Canon (n = 311), and VCTE-M probe (n = 812). VCTE-M and 2D-SWE-SSI values shared the highest correlation and concordance coefficients (0.933 and 0.920, respectively) and a coefficient of variation below 20%, whatever the range of values. 2D-SWE-GE had the lowest variations, with 2D-SWE-SSI values below 13 kPa. However, both 2D-SWE-GE and 2D-SWE-Canon exhibited a frank underestimation of the high percentiles' 2D-SWE-SSI values with coefficients of variation of -21.7% and -25.8% from 13- to 17-kPa values, and -44.3% and -32.4% from 17-kPa values, respectively. In conclusion, knowledge of the vendor-specific distribution of values is mandatory for interpreting results obtained with different machines. If all four techniques behave closely in low values allowing excluding advanced chronic liver diseases in larger populations, discrepancies are observed in high percentile values.

Liver fibrosis assessment: MR and US elastography

Elastography has emerged as a preferred non-invasive imaging technique for the clinical assessment of liver fibrosis. Elastography methods provide liver stiffness measurement (LSM) as a surrogate quantitative biomarker for fibrosis burden in chronic liver disease (CLD). Elastography can be performed either with ultrasound or MRI. Currently available ultrasound-based methods include strain elastography, two-dimensional shear wave elastography (2D-SWE), point shear wave elastography (pSWE), and vibration-controlled transient elastography (VCTE). MR Elastography (MRE) is widely available as two-dimensional gradient echo MRE (2D-GRE-MRE) technique. US-based methods provide estimated Young's modulus (eYM) and MRE provides magnitude of the complex shear modulus. MRE and ultrasound methods have proven to be accurate methods for detection of advanced liver fibrosis and cirrhosis. Other clinical applications of elastography include liver decompensation prediction, and differentiation of non-alcoholic steatohepatitis (NASH) from simple steatosis (SS). In this review, we briefly describe the different elastography methods, discuss current clinical applications, and provide an overview of advances in the field of liver elastography.

Development of a Robotic Shear Wave Elastography System for Noninvasive Staging of Liver Disease in Murine Models

Shear wave elastography (SWE) is an ultrasound-based stiffness quantification technology that is used for noninvasive liver fibrosis assessment. However, despite widescale clinical adoption, SWE is largely unused by preclinical researchers and drug developers for studies of liver disease progression in small animal models due to significant experimental, technical, and reproducibility challenges. Therefore, the aim of this work was to develop a tool designed specifically for assessing liver stiffness and echogenicity in small animals to better enable longitudinal preclinical studies. A high-frequency linear array transducer (12-24 MHz) was integrated into a robotic small animal ultrasound system (Vega; SonoVol, Inc., Durham, NC) to perform liver stiffness and echogenicity measurements in three dimensions. The instrument was validated with tissue-mimicking phantoms and a mouse model of nonalcoholic steatohepatitis. Female C57BL/6J mice (n = 40) were placed on choline-deficient, L-amino acid-defined, high-fat diet and imaged longitudinally for 15 weeks. A subset was sacrificed after each imaging timepoint (n = 5) for histological validation, and analyses of receiver operating characteristic (ROC) curves were performed. Results demonstrated that robotic measurements of echogenicity and stiffness were most strongly correlated with macrovesicular steatosis (R²  = 0.891) and fibrosis (R²  = 0.839), respectively. For diagnostic classification of fibrosis (Ishak score), areas under ROC (AUROCs) curves were 0.969 for ≥Ishak1, 0.984 for ≥Ishak2, 0.980 for ≥Ishak3, and 0.969 for ≥Ishak4. For classification of macrovesicular steatosis (S-score), AUROCs were 1.00 for ≥S2 and 0.997 for ≥S3. Average scanning and analysis time was <5 minutes/liver. Conclusion: Robotic SWE in small animals is feasible and sensitive to small changes in liver disease state, facilitating in vivo staging of rodent liver disease with minimal sonographic expertise.

Comparison of inversion methods in MR elastography: An open-access pipeline for processing multifrequency shear-wave data and demonstration in a phantom, human kidneys, and brain

PURPOSE

Magnetic resonance elastography (MRE) maps the viscoelastic properties of soft tissues for diagnostic purposes. However, different MRE inversion methods yield different results, which hinder comparison of values, standardization, and establishment of quantitative MRE markers. Here, we introduce an expandable, open-access, webserver-based platform that offers multiple inversion techniques for multifrequency, 3D MRE data.

METHODS

The platform comprises a data repository and standard MRE inversion methods including local frequency estimation (LFE), direct-inversion based multifrequency dual elasto-visco (MDEV) inversion, and wavenumber-based (k-) MDEV. The use of the platform is demonstrated in phantom data and in vivo multifrequency MRE data of the kidneys and brains of healthy volunteers.

RESULTS

Detailed maps of stiffness were generated by all inversion methods showing similar detail of anatomy. Specifically, the inner renal cortex had higher shear wave speed (SWS) than renal medulla and outer cortex without lateral differences. k-MDEV yielded higher SWS values than MDEV or LFE (full kidney/brain k-MDEV: 2.71 ± 0.19/1.45 ± 0.14 m/s, MDEV: 2.14 ± 0.16/0.99 ± 0.11 m/s, LFE: 2.12 ± 0.15/0.89 ± 0.06 m/s).

CONCLUSION

The freely accessible platform supports the comparison of MRE results obtained with different inversion methods, filter thresholds, or excitation frequencies, promoting reproducibility in MRE across community-developed methods.

Spatiotemporal Deconvolution of Hydrophone Response for Linear and Nonlinear Beams-Part I: Theory, Spatial-Averaging Correction Formulas, and Criteria for Sensitive Element Size

This article reports spatiotemporal deconvolution methods and simple empirical formulas to correct pressure and beamwidth measurements for spatial averaging across a hydrophone sensitive element. Readers who are uninterested in hydrophone theory may proceed directly to Appendix A for an easy method to estimate spatial-averaging correction factors. Hydrophones were modeled as angular spectrum filters. Simulations modeled nine circular transducers (1-10 MHz; F/1.4-F/3.2) driven at six power levels and measured with eight hydrophones (432 beam/hydrophone combinations). For example, the model predicts that if a 200- [Formula: see text] membrane hydrophone measures a moderately nonlinear 5-MHz beam from an F/1 transducer, spatial-averaging correction factors are 33% (peak compressional pressure or p_c ), 18% (peak rarefactional pressure or p ), and 18% (full width half maximum or FWHM). Theoretical and experimental estimates of spatial-averaging correction factors to were in good agreement (within 5%) for linear and moderately nonlinear signals. Criteria for maximum appropriate hydrophone sensitive element size as functions of experimental parameters were derived. Unlike the oft-cited International Electrotechnical Commission (IEC) criterion, the new criteria were derived for focusing rather than planar transducers and can accommodate nonlinear signals in addition to linear signals. Responsible reporting of hydrophone-based pressure and beamwidth measurements should always acknowledge spatial-averaging considerations.

Liquid-Liver Phantom: Mimicking the Viscoelastic Dispersion of Human Liver for Ultrasound- and MRI-Based Elastography

OBJECTIVES

Tissue stiffness can guide medical diagnoses and is exploited as an imaging contrast in elastography. However, different elastography devices show different liver stiffness values in the same subject, hindering comparison of values and establishment of system-independent thresholds for disease detection. There is a need for standardized phantoms that specifically address the viscosity-related dispersion of stiffness over frequency. To improve standardization of clinical elastography across devices and platforms including ultrasound and magnetic resonance imaging (MRI), a comprehensively characterized phantom is introduced that mimics the dispersion of stiffness of the human liver and can be generated reproducibly.

MATERIALS AND METHODS

The phantom was made of linear polymerized polyacrylamide (PAAm) calibrated to the viscoelastic properties of healthy human liver in vivo as reported in the literature. Stiffness dispersion was analyzed using the 2-parameter springpot model fitted to the dispersion of shear wave speed of PAAm, which was measured by shear rheometry, ultrasound-based time-harmonic elastography, clinical magnetic resonance elastography (MRE), and tabletop MRE in the frequency range of 5 to 3000 Hz. Imaging parameters for ultrasound and MRI, reproducibility, aging behavior, and temperature dependency were assessed. In addition, the frequency bandwidth of shear wave speed of clinical elastography methods (Aplio i900, Canon; Acuson Sequoia, Siemens; FibroScan, EchoSense) was characterized.

RESULTS

Within the entire frequency range analyzed in this study, the PAAm phantom reproduced well the stiffness dispersion of human liver in vivo despite its fluid properties under static loading (springpot stiffness parameter, 2.14 [95% confidence interval, 2.08-2.19] kPa; springpot powerlaw exponent, 0.367 [95% confidence interval, 0.362-0.373]). Imaging parameters were close to those of liver in vivo with only slight variability in stiffness values of 0.5% (0.4%, 0.6%), 4.1% (3.9%, 4.5%), and -0.63% (-0.67%, -0.58%), respectively, between batches, over a 6-month period, and per °C increase in temperature.

CONCLUSIONS

The liquid-liver phantom has useful properties for standardization and development of liver elastography. First, it can be used across clinical and experimental elastography devices in ultrasound and MRI. Second, being a liquid, it can easily be adapted in size and shape to specific technical requirements, and by adding inclusions and scatterers. Finally, because the phantom is based on noncrosslinked linear PAAm constituents, it is easy to produce, indicating potential widespread use among researchers and vendors to standardize liver stiffness measurements.

Shear wave speed measurement bias in a viscoelastic phantom across six ultrasound elastography systems: a comparative study with transient elastography and magnetic resonance elastography

Purpose

To quantify the bias of shear wave speed (SWS) measurements in a viscoelastic phantom across six different ultrasound (US) systems and to compare the SWS with those from transient elastography (TE) and magnetic resonance elastography (MRE).

Methods

A viscoelastic phantom of stiffness representing fibrotic liver or healthy thyroid was measured with nine (linear probe) and 10 (convex probe) modes of six different US-based shear wave elastography (SWE) systems using linear and convex probes. SWS measurements of three regions of interest were repeated thrice at two focal depths, coupling the probe to the phantom using a jig. An MRE system using three motion-encoding gradient frequencies of 60, 90, and 120 Hz and TE were also used to measure the stiffness of the phantom.

Results

The SWS from different SWE systems had mean coefficients of variation of 9.0–9.2% and 5.4–5.6% with linear and convex probes, respectively, in viscoelastic phantom measurement. The focal depth was a less significant source of SWS variability than the system. The total average SWS obtained with US-SWE systems was 19.9% higher than that obtained with MRE at 60 Hz, which is commonly used in clinical practice, and 31.5% higher than that obtained with TE using the M probe.

Conclusions

Despite the measurement biases associated with the SWE systems, biases were not necessarily consistent, and they changed with the probes used and depth measured. The SWS of the viscoelastic phantom obtained using different modalities increased according to the shear wave frequency used.

Full Characterization of in vivo Muscle as an Elastic, Incompressible, Transversely Isotropic Material Using Ultrasonic Rotational 3D Shear Wave Elasticity Imaging

Using a 3D rotational shear wave elasticity imaging (SWEI) setup, 3D shear wave data were acquired in the vastus lateralis of a healthy volunteer. The innate tilt between the transducer face and the muscle fibers results in the excitation of multiple shear wave modes, allowing for more complete characterization of muscle as an elastic, incompressible, transversely isotropic (ITI) material. The ability to measure both the shear vertical (SV) and shear horizontal (SH) wave speed allows for measurement of three independent parameters needed for full ITI material characterization: the longitudinal shear modulus μL , the transverse shear modulus μT , and the tensile anisotropy χE . Herein we develop and validate methodology to estimate these parameters and measure them in vivo, with μL = 5.77±1.00 kPa, μT = 1.93±0.41 kPa (giving shear anisotropy χμ = 2.11±0.92 ), and χE = 4.67±1.40 in a relaxed vastus lateralis muscle. We also demonstrate that 3D SWEI can be used to more accurately characterize muscle mechanical properties as compared to 2D SWEI.

Assessing cardiac stiffness using ultrasound shear wave elastography

Shear wave elastography offers a new dimension to echocardiography: it measures myocardial stiffness. Therefore, it could provide additional insights into the pathophysiology of cardiac diseases affecting myocardial stiffness and potentially improve diagnosis or guide patient treatment. The technique detects fast mechanical waves on the heart wall with high frame rate echography, and converts their propagation velocity into a stiffness value. A proper interpretation of shear wave data is required as the shear wave interacts with the intrinsic, yet dynamically changing geometrical and material characteristics of the heart under pressure. This dramatically alters the wave physics of the propagating wave, demanding adapted processing methods compared to other shear wave elastography applications as breast tumor and liver stiffness staging. Furthermore, several advanced analysis methods have been proposed to extract supplementary material features such as viscosity and anisotropy, potentially offering additional diagnostic value. This review explains the general mechanical concepts underlying cardiac shear wave elastography and provides an overview of the preclinical and clinical studies within the field. We also identify the mechanical and technical challenges ahead to make shear wave elastography a valuable tool for clinical practice.

Sources of Variability in Shear Wave Speed and Dispersion Quantification with Ultrasound Elastography: A Phantom Study

There is a growing interest in quantifying shear-wave dispersion (SWD) with ultrasound shear-wave elastography (SWE). Recent studies suggest that SWD complements shear-wave speed (SWS) in diffuse liver disease diagnosis. To accurately interpret these metrics in clinical practice, we analyzed the impact of operator-dependent acquisition parameters on SWD and SWS measurements. Considered parameters were the acquisition depth, lateral position and size of the region of interest (ROI), as well as the size of the SWE acquisition box. Measurements were performed using the Canon Aplio i800 system (Canon Medical Systems, Otawara, Tochigi, Japan) and four homogeneous elasticity phantoms with certified stiffness values ranging from 3.7 to 44 kPa. In general, SWD exhibited two to three times greater variability than SWS. The acquisition depth was the main variance-contributing factor for both SWS and SWD, which decayed significantly with depth. The lateral ROI position contributed as much as the acquisition depth to the total variance in SWD. Locations close to the initial shear-wave excitation pulse were more robust to biases because of inaccurate probe-phantom coupling. The size of the ROI and acquisition box did not introduce significant variations. These results suggest that future guidelines on multiparametric elastography should account for the depth- and lateral-dependent variability of measurements.

Evaluation of Robustness of Local Phase Velocity Imaging in Homogenous Tissue-Mimicking Phantoms

Shear wave elastography (SWE) is a method of evaluating mechanical properties of soft tissues. Most current implementations of SWE report the group velocity for shear wave velocity, which assumes an elastic, isotropic, homogenous and incompressible tissue. Local phase velocity imaging (LPVI) is a novel method of phase velocity reconstruction that allows for accurate evaluation of shear wave velocity at specified frequencies. This method's robustness was evaluated in 11 elastic and 8 viscoelastic phantoms using linear and curvilinear arrays. We acquired data with acoustic radiation force push beams with different focal depths and F-numbers and reconstructed phase velocity images over a wide range of frequencies. Regardless of phantom, push beam focal depth and reconstruction frequency, an F-number around 3.0 was found to produce the largest usable area in the phase velocity reconstructions. For elastic phantoms scanned with a linear array, the optimal focal depth, frequency range and maximum region of interest (ROI) were 20-30 mm, 100-400 Hz and 2.70 cm², respectively. For viscoelastic phantoms scanned with a linear array, the optimal focal depth, frequency and maximum ROI were 20-30 mm, 100-300 Hz and 1.54 cm², respectively. For the curvilinear array in the same phantoms, optimal focal depth, frequency range and maximum ROIs were 45-60 mm, 100-400 and 100-300 Hz and 1.54 cm², respectively. In further work, LPVI reconstructions from inclusion phantoms will be evaluated to simulate non-homogeneous tissues. Additionally, LPVI will be evaluated in larger-volume phantoms to account for wave reflection from the containers when using the curvilinear array.

Novel Uses of Ultrasound to Assess Kidney Mechanical Properties

Ultrasound is a key imaging tool for evaluating the kidney. Over the last two decades, methods to measure the mechanical properties of soft tissues have been developed and used in clinical practice, although use in the kidney has not been as widespread as for other applications. The mechanical properties of the kidney are determined by the structure and composition of the renal parenchyma and perfusion characteristics. Because pathologic processes change these factors, the mechanical properties change and can be used for diagnostic purposes and for monitoring treatment or disease progression. Ultrasound-based elastography methods for evaluating the mechanical properties of the kidney use focused ultrasound beams to perturb the kidney and then high frame-rate ultrasound methods are used to measure the resulting motion. The motion is analyzed to estimate the mechanical properties. This review will describe the principles of these methods and discuss several seminal studies related to characterizing the kidney. Additionally, an overview of the clinical use of elastography methods in native and kidney allografts will be provided. Perspectives on future developments and uses of elastography technology along with other complementary ultrasound imaging modalities will be provided.