Viscoelastic response (VisR) imaging for assessment of viscoelasticity in Voigt materials

Quantitative Viscoelastic Response (QVisR): Direct Estimation of Viscoelasticity With Neural Networks

We present a machine learning method to directly estimate viscoelastic moduli from displacement time-series profiles generated by viscoelastic response (VisR) ultrasound excitations. VisR uses two colocalized acoustic radiation force (ARF) pushes to approximate tissue viscoelastic creep response and tracks displacements on-axis to measure the material relaxation. A fully connected neural network is trained to learn a nonlinear mapping from VisR displacements, the push focal depth, and the measurement axial depth to the material elastic and viscous moduli. In this work, we assess the validity of quantitative VisR (QVisR) in simulated materials, propose a method of domain adaption to phantom VisR displacements, and show in vivo estimates from a clinically acquired dataset.

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.

Quantitative Evaluation of the Lumbar Multifidus Muscle by Shear Wave Dispersion in Healthy Adults: A Preliminary Study

OBJECTIVES

To investigate the application value of shear wave dispersion (SWD) in healthy adults with the lumbar multifidus muscle (LMM), to determine the range of normal reference values, and to analyze the influences of factors on the parameter.

METHODS

Ninety-five healthy volunteers participated in the study, from whom 2-dimensional, shear wave elastography (SWE), and SWD images of the bilateral LMM were acquired in three positions (prone, standing, and anterior flexion). Subcutaneous fat thickness (SFH), SWE velocity, and SWD slope were measured accordingly for analyses.

RESULTS

The mean SWD slope of the bilateral LMM in the prone position was as follows: left: 14.8 ± 3.1 (m/second)/kHz (female) and 13.0 ± 2.5 (m/second)/kHz (male); right: 14.8 ± 3.7 (m/second)/kHz (female) and 14.2 ± 3.4 (m/second)/kHz (male). In the prone position, there was a weak negative correlation between the bilateral LMM SWD slope of activity level 2 and level 1 (β = -1.5 (2 versus 1, left), -1.9 (2 versus 1, right), all P < .05), and between the left SWD slope of activity level 3 and level 1 (β = -2.3 [3 versus 1, left], P < .05). The correlation between SWE velocity and SWD slope value changed with the position: there was a weak positive correlation in the prone position (r = 0.3 [left], 0.37 [right], both P < .05), and a moderate positive correlation in the standing and anterior flexed positions (r = 0.49-0.74, both P < .001). SFH was moderately negatively correlated with bilateral SWD slope values in the anterior flexion (left: r = -0.4, P = .01; right: r = -0.7, P < .01).

CONCLUSIONS

SWD imaging can be used as an adjunct tool to aid in the assessment of viscosity in LMM. Further, activity level, and position are influencing factors that should be considered in clinical practice.

Real-Time Nondestructive Viscosity Measurement of Soft Tissue Based on Viscoelastic Response Optical Coherence Elastography

Viscoelasticity of the soft tissue is an important mechanical factor for disease diagnosis, biomaterials testing and fabrication. Here, we present a real-time and high-resolution viscoelastic response-optical coherence elastography (VisR-OCE) method based on acoustic radiation force (ARF) excitation and optical coherence tomography (OCT) imaging. The relationship between displacements induced by two sequential ARF loading-unloading and the relaxation time constant of the soft tissue-is established for the Kelvin-Voigt material. Through numerical simulation, the optimal experimental parameters are determined, and the influences of material parameters are evaluated. Virtual experimental results show that there is less than 4% fluctuation in the relaxation time constant values obtained when various Young's modulus and Poisson's ratios were given for simulation. The accuracy of the VisR-OCE method was validated by comparing with the tensile test. The relaxation time constant of phantoms measured by VisR-OCE differs from the tensile test result by about 3%. The proposed VisR-OCE method may provide an effective tool for quick and nondestructive viscosity testing of biological tissues.

Comparing Focused-Tracked and Plane Wave-Tracked ARFI Log(VoA) In Silico and in Application to Human Carotid Atherosclerotic Plaque, Ex Vivo

A significant risk factor for ischemic stroke is carotid atherosclerotic plaque that is susceptible to rupture, with rupture potential conveyed by plaque morphology. Human carotid plaque composition and structure have been delineated noninvasively and in vivo by evaluating log(VoA), a parameter derived as the decadic log of the second time derivative of displacement induced by an acoustic radiation force impulse (ARFI). In prior work, ARFI-induced displacement was measured using conventional focused tracking; however, this requires a long data acquisition period, thereby reducing framerate. We herein evaluate if ARFI log(VoA) framerate can be increased without a reduction in plaque imaging performance using plane wave tracking instead. In silico, both focused- and plane wave-tracked log(VoA) decreased with increasing echobrightness, quantified as signal-to-noise ratio (SNR), but did not vary with material elasticity for SNRs below 40 dB. For SNRs of 40-60 dB, both focused- and plane wave-tracked log(VoA) varied with SNR and material elasticity. Above 60 dB SNR, both focused- and plane wave-tracked log(VoA) varied with material elasticity alone. This suggests that log(VoA) discriminates features according to a combination of their echobrightness and mechanical property. Further, while both focused- and plane-wave tracked log(VoA) values were artifactually inflated by mechanical reflections at inclusion boundaries, plane wave-tracked log(VoA) was more strongly impacted by off-axis scattering. Applied to three excised human cadaveric carotid plaques with spatially aligned histological validation, both log(VoA) methods detected regions of lipid, collagen, and calcium (CAL) deposits. These findings support that plane wave tracking performs comparably to focused tracking for log(VoA) imaging and that plane wave-tracked log(VoA) is a viable approach to discriminating clinically relevant atherosclerotic plaque features at a 30-fold higher framerate than by focused tracking.

Improving in situ acoustic intensity estimates using MR acoustic radiation force imaging in combination with multifrequency MR elastography

PURPOSE

Magnetic resonance acoustic radiation force imaging (MR-ARFI) enables focal spot localization during nonablative transcranial ultrasound therapies. As the acoustic radiation force is proportional to the applied acoustic intensity, measured MR-ARFI displacements could potentially be used to estimate the acoustic intensity at the target. However, variable brain stiffness is an obstacle. The goal of this study was to develop and assess a method to accurately estimate the acoustic intensity at the focus using MR-ARFI displacements in combination with viscoelastic properties obtained with multifrequency MR elastography (MRE).

METHODS

Phantoms with a range of viscoelastic properties were fabricated, and MR-ARFI displacements were acquired within each phantom using multiple acoustic intensities. Voigt model parameters were estimated for each phantom based on storage and loss moduli measured using multifrequency MRE, and these were used to predict the relationship between acoustic intensity and measured displacement.

RESULTS

Using assumed viscoelastic properties, MR-ARFI displacements alone could not accurately estimate acoustic intensity across phantoms. For example, acoustic intensities were underestimated in phantoms stiffer than the assumed stiffness and overestimated in phantoms softer than the assumed stiffness. This error was greatly reduced using individualized viscoelasticity measurements obtained from MRE.

CONCLUSION

We demonstrated that viscoelasticity information from MRE could be used in combination with MR-ARFI displacements to obtain more accurate estimates of acoustic intensity. Additionally, Voigt model viscosity parameters were found to be predictive of the relaxation rate of each phantom's time-varying displacement response, which could be used to optimize patient-specific MR-ARFI pulse sequences.

Acoustic Radiation Force: A Review of Four Mechanisms for Biomedical Applications

Radiation force is a universal phenomenon in any wave motion where the wave energy produces a static or transient force on the propagation medium. The theory of acoustic radiation force (ARF) dates back to the early 19th century. In recent years, there has been an increasing interest in the biomedical applications of ARF. Following a brief history of ARF, this article describes a concise theory of ARF under four physical mechanisms of radiation force generation in tissue-like media. These mechanisms are primarily based on the dissipation of acoustic energy of propagating waves, the reflection of the incident wave, gradients of the compressional wave speeds, and the spatial variations of energy density in standing acoustic waves. Examples describing some of the practical applications of ARF under each mechanism are presented. This article concludes with a discussion on selected ideas for potential future applications of ARF in biomedicine.

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.

Effects of Loading and Boundary Conditions on the Performance of Ultrasound Compressional Viscoelastography: A Computational Simulation Study to Guide Experimental Design

Most biomaterials and tissues are viscoelastic; thus, evaluating viscoelastic properties is important for numerous biomedical applications. Compressional viscoelastography is an ultrasound imaging technique used for measuring the viscoelastic properties of biomaterials and tissues. It analyzes the creep behavior of a material under an external mechanical compression. The aim of this study is to use finite element analysis to investigate how loading conditions (the distribution of the applied compressional pressure on the surface of the sample) and boundary conditions (the fixation method used to stabilize the sample) can affect the measurement accuracy of compressional viscoelastography. The results show that loading and boundary conditions in computational simulations of compressional viscoelastography can severely affect the measurement accuracy of the viscoelastic properties of materials. The measurement can only be accurate if the compressional pressure is exerted on the entire top surface of the sample, as well as if the bottom of the sample is fixed only along the vertical direction. These findings imply that, in an experimental validation study, the phantom design should take into account that the surface area of the pressure plate must be equal to or larger than that of the top surface of the sample, and the sample should be placed directly on the testing platform without any fixation (such as a sample container). The findings indicate that when applying compressional viscoelastography to real tissues in vivo, consideration should be given to the representative loading and boundary conditions. The findings of the present simulation study will provide a reference for experimental phantom designs regarding loading and boundary conditions, as well as guidance towards validating the experimental results of compressional viscoelastography.

Hydrophone Spatial Averaging Correction for Acoustic Exposure Measurements From Arrays-Part I: Theory and Impact on Diagnostic Safety Indexes

This article reports underestimation of mechanical index (MI) and nonscanned thermal index for bone near focus (TIB) due to hydrophone spatial averaging effects that occur during acoustic output measurements for clinical linear and phased arrays. TIB is the appropriate version of thermal index (TI) for fetal imaging after ten weeks from the last menstrual period according to the American Institute of Ultrasound in Medicine (AIUM). Spatial averaging is particularly troublesome for highly focused beams and nonlinear, nonscanned modes such as acoustic radiation force impulse (ARFI) and pulsed Doppler. MI and variants of TI (e.g., TIB), which are displayed in real-time during imaging, are often not corrected for hydrophone spatial averaging because a standardized method for doing so does not exist for linear and phased arrays. A novel analytic inverse-filter method to correct for spatial averaging for pressure waves from linear and phased arrays is derived in this article (Part I) and experimentally validated in a companion article (Part II). A simulation was developed to estimate potential spatial-averaging errors for typical clinical ultrasound imaging systems based on the theoretical inverse filter and specifications for 124 scanner/transducer combinations from the U.S. Food and Drug Administration (FDA) 510(k) database from 2015 to 2019. Specifications included center frequency, aperture size, acoustic output parameters, hydrophone geometrical sensitive element diameter, etc. Correction for hydrophone spatial averaging using the inverse filter suggests that maximally achievable values for MI, TIB, thermal dose ( t ₄₃ ), and spatial-peak-temporal-average intensity ( [Formula: see text]) for typical clinical systems are potentially higher than uncorrected values by (means ± standard deviations) 9% ± 4% (ARFI MI), 19% ± 15% (ARFI TIB), 50% ± 41% (ARFI t ₄₃ ), 43% ± 39% (ARFI [Formula: see text]), 7% ± 5% (pulsed Doppler MI), 15% ± 11% (pulsed Doppler TIB), 42% ± 31% (pulsed Doppler t ₄₃ ), and 33% ± 27% (pulsed Doppler [Formula: see text]). These values correspond to frequencies of 3.2 ± 1.3 (ARFI) and 4.1 ± 1.4 MHz (pulsed Doppler), and the model predicts that they would increase with frequency. Inverse filtering for hydrophone spatial averaging significantly improves the accuracy of estimates of MI, TIB, t ₄₃ , and [Formula: see text] for ARFI and pulsed Doppler signals.

Resonant acoustic rheometry for non-contact characterization of viscoelastic biomaterials

Resonant Acoustic Rheometry (RAR) is a new, non-contact technique to characterize the mechanical properties of soft and viscoelastic biomaterials, such as hydrogels, that are used to mimic the extracellular matrix in tissue engineering. RAR uses a focused ultrasound pulse to generate a microscale perturbation at the sample surface and tracks the ensuing surface wave using pulse-echo ultrasound. The frequency spectrum of the resonant surface waves is analyzed to extract viscoelastic material properties. In this study, RAR was used to characterize fibrin, gelatin, and agarose hydrogels. Single time point measurements of gelled samples with static mechanical properties showed that RAR provided consistent quantitative data and measured intrinsic material characteristics independent of ultrasound parameters. RAR was also used to longitudinally track dynamic changes in viscoelastic properties over the course of fibrin gelation, revealing distinct phase and material property transitions. Application of RAR was verified using finite element modeling and the results were validated against rotational shear rheometry. Importantly, RAR circumvents some limitations of conventional rheology methods and can be performed in a high-throughput manner using conventional labware. Overall, these studies demonstrate that RAR can be a valuable tool to noninvasively quantify the viscoelastic mechanical properties of soft hydrogel biomaterials.

Ultrasound neuromodulation depends on pulse repetition frequency and can modulate inhibitory effects of TTX

Ultrasound is gaining traction as a neuromodulation method due to its ability to remotely and non-invasively modulate neuronal activity with millimeter precision. However, there is little consensus about optimal ultrasound parameters required to elicit neuromodulation and how specific parameters drive mechanisms that underlie ultrasound neuromodulation. We address these questions in this work by performing a study to determine effective ultrasound parameters in a transgenic mouse brain slice model that enables calcium imaging as a quantitative readout of neuronal activity for ultrasound neuromodulation. We report that (1) calcium signaling increases with the application of ultrasound; (2) the neuronal response rate to ultrasound is dependent on pulse repetition frequency (PRF); and (3) ultrasound can reversibly alter the inhibitory effects of tetrodotoxin (TTX) in pharmacological studies. This study offers mechanistic insight into the PRF dependence of ultrasound neuromodulation and the nature of ultrasound/ion channel interaction.

Viscoelastic Response Ultrasound Derived Relative Elasticity and Relative Viscosity Reflect True Elasticity and Viscosity: In Silico and Experimental Demonstration

Viscoelastic response (VisR) ultrasound characterizes the viscoelastic properties of tissue by fitting acoustic radiation force (ARF)-induced displacements in the region of ARF excitation to a 1-D mass-spring-damper (MSD) model. Elasticity and viscosity are calculated separately but relative to the applied ARF amplitude. We refer to these parameters as "relative elasticity (RE)" and "relative viscosity (RV)." We herein test the hypothesis that RE and RV linearly correlate to true elasticity and viscosity in tissue. VisR imaging was simulated in 144 homogeneous viscoelastic materials with varying elasticities and viscosities. Derived RE linearly correlated with material elasticity and varied by an average of 2.52% when the material viscosity changed from 0.1 to 1.3 Pa · s. Derived RV linearly correlated with material viscosity but varied by an average of 102.5% when material elasticity changed from 3.33 to 20 kPa. The effect of elasticity on RV measurement was compensated using the slope of the linear relationship between RV and natural frequency ( ωtextn ). After compensation, RV [Formula: see text] (elasticity compensated RV) linearly correlated with material viscosity and varied by less than 1.00% on average when the modeled shear elastic modulus changed from 3.3 to 20 kPa. In addition to elasticity compensation, variation in ARF amplitude over depth was compensated, yielding REDC and [Formula: see text]. REDC and [Formula: see text] successfully contrasted elastic and viscous inclusions, respectively, in three simulated phantoms. Experimentally, in the homogeneous oil-in-gelatin phantoms and excised livers, REDC linearly correlated with shear wave dispersion ultrasound vibrometry (SDUV) derived shear elastic modulus, and [Formula: see text] linearly correlated with SDUV-derived shear viscosity. In excised livers containing viscoelastic oil-in-gelatin inclusions, the inclusions were successfully contrasted from the liver background by both REDC and [Formula: see text]. These results suggest that RE and RV are relevant for qualitatively assessing the elastic and viscous properties of tissue.

Spectral Quantification of Nonlinear Elasticity Using Acoustoelasticity and Shear-Wave Dispersion

Tissue biomechanical properties are known to be sensitive to pathological changes. Accordingly, various techniques have been developed to estimate tissue mechanical properties. Shear-wave elastography (SWE) measures shear-wave speed (SWS) in tissues, which can be related to shear modulus. Although viscosity or stress-strain nonlinearity may act as confounder of SWE, their explicit characterization may also provide additional information about tissue composition as a contrast modality. Viscosity can be related to frequency dispersion of SWS, which can be characterized using multi-frequency measurements, herein called spectral SWE (SSWE). Additionally, nonlinear shear modulus can be quantified and parameterized based on SWS changes with respect to applied stress, a phenomenon called acoustoelasticity (AE). In this work, we characterize the nonlinear parameters of tissue as a function of excitation frequency by utilizing both AE and SSWE together. For this, we apply incremental amounts of quasi-static stress on a medium, while imaging and quantifying SWS dispersion via SSWE. Results from phantom and ex vivo porcine liver experiments demonstrate the feasibility of measuring frequency-dependent nonlinear parameters using the proposed method. SWS propagation in porcine liver tissue was observed to change from 1.8 m/s at 100 Hz to 3.3 m/s at 700 Hz, while increasing by approximately 25% from a strain of 0% to 12% across these frequencies.

Mechanical Anisotropy Assessment in Kidney Cortex Using ARFI Peak Displacement: Preclinical Validation and Pilot In Vivo Clinical Results in Kidney Allografts

The kidney is an anisotropic organ, with higher elasticity along versus across nephrons. The degree of mechanical anisotropy in the kidney may be diagnostically relevant if properly exploited; however, if improperly controlled, anisotropy may confound stiffness measurements. The purpose of this study is to demonstrate the clinical feasibility of acoustic radiation force (ARF)-induced peak displacement (PD) measures for both exploiting and obviating mechanical anisotropy in the cortex of human kidney allografts, in vivo. Validation of the imaging methods is provided by preclinical studies in pig kidneys, in which ARF-induced PD values were significantly higher ( , Wilcoxon) when the transducer executing asymmetric ARF was oriented across versus along the nephrons. The ratio of these PD values obtained with the transducer oriented across versus along the nephrons strongly linearly correlated ( R2 = 0.95 ) to the ratio of shear moduli measured by shear wave elasticity imaging. On the contrary, when a symmetric ARF was implemented, no significant difference in PD was observed ( p > 0.01 ). Similar results were demonstrated in vivo in the kidney allografts of 14 patients. The symmetric ARF produced PD measures with no significant difference ( p > 0.01 ) between along versus across alignments, but the asymmetric ARF yielded PD ratios that remained constant over a six-month observation period post-transplantation, consistent with stable serum creatinine level and urine protein-to-creatinine ratio in the same patient population ( p > 0.01 ). The results of this pilot in vivo clinical study suggest the feasibility of 1) implementing symmetrical ARF to obviate mechanical anisotropy in the kidney cortex when anisotropy is a confounding factor and 2) implementing asymmetric ARF to exploit mechanical anisotropy when mechanical anisotropy is a potentially relevant biomarker.

Production of acoustic radiation force using ultrasound: methods and applications

INTRODUCTION

Acoustic radiation force (ARF) is used in many biomedical applications. The transfer of momentum in acoustic waves can be used in a multitude of ways to perturb tissue and manipulate cells.

AREAS COVERED

This review will briefly cover the acoustic theory related to ARF, particularly that related to application in tissues. The use of ARF in measurement of mechanical properties will be treated in detail with emphasis on the spatial and temporal modulation of the ARF. Additional topics covered will be the manipulation of particles with ARF, correction of phase aberration with ARF, modulation of cellular behavior with ARF, and bioeffects related to ARF use.

EXPERT COMMENTARY

The use of ARF can be tailored to specific applications for measurements of mechanical properties or correction of focusing for ultrasound beams. Additionally, noncontact manipulation of particles and cells with ARF enables a wide array of applications for tissue engineering and biosensing.

Viscoelastic parameters as discriminators of breast masses: Initial human study results

Shear wave elastography is emerging as a clinically valuable diagnostic tool to differentiate between benign and malignant breast masses. Elastography techniques assume that soft tissue can be modelled as a purely elastic medium. However, this assumption is often violated as soft tissue exhibits viscoelastic properties. In order to explore the role of viscoelastic parameters in suspicious breast masses, a study was conducted on a group of patients using shear wave dispersion ultrasound vibrometry in the frequency range of 50–400 Hz. A total of 43 female patients with suspicious breast masses were recruited before their scheduled biopsy. Of those, 15 patients did not meet the data selection criteria. Voigt model based shear elasticity showed a significantly (p = 7.88x10^-6) higher median value for the 13 malignant masses (16.76±13.10 kPa) compared to 15 benign masses (1.40±1.12 kPa). Voigt model based shear viscosity was significantly different (p = 4.13x10^-5) between malignant (8.22±3.36 Pa-s) and benign masses (2.83±1.47 Pa-s). Moreover, the estimated time constant from the Voigt model, which is dependent on both shear elasticity and viscosity, differed significantly (p = 6.13x10^-5) between malignant (0.68±0.33 ms) and benign masses (3.05±1.95 ms). Results suggest that besides elasticity, viscosity based parameters like shear viscosity and time constant can also be used to differentiate between malignant and benign breast masses.

Evaluating Renal Transplant Status Using Viscoelastic Response (VisR) Ultrasound

Chronic kidney disease is most desirably and cost-effectively treated by renal transplantation, but graft survival is a major challenge. Although irreversible graft damage can be averted by timely treatment, intervention is delayed when early graft dysfunction goes undetected by standard clinical metrics. A more sensitive and specific parameter for delineating graft health could be the viscoelastic properties of the renal parenchyma, which are interrogated non-invasively by Viscoelastic Response (VisR) ultrasound, a new acoustic radiation force (ARF)-based imaging method. Assessing the performance of VisR imaging in delineating histologically confirmed renal transplant pathologies in vivo is the purpose of the study described here. VisR imaging was performed in patients with (n = 19) and without (n = 25) clinical indication for renal allograft biopsy. The median values of VisR outcome metrics (τ, relative elasticity [RE] and relative viscosity [RV]) were calculated in five regions of interest that were manually delineated in the parenchyma (outer, center and inner) and in the pelvis (outer and inner). The ratios of a given VisR metric for all possible region-of-interest combinations were calculated, and the corresponding ratios were statistically compared between biopsied patients subdivided by diagnostic categories versus non-biopsied, control allografts using the two-sample Wilcoxon test (p <0.05). Although τ ratios non-specifically differentiated allografts with vascular disease, tubular/interstitial scarring, chronic allograft nephropathy and glomerulonephritis from non-biopsied control allografts, RE distinguished only allografts with vascular disease and tubular/interstitial scarring, and RV distinguished only vascular disease. These results suggest that allografts with scarring and vascular disease can be identified using non-invasive VisR RE and RV metrics.

Acoustic Radiation Force-Induced Creep-Recovery (ARFICR): A Noninvasive Method to Characterize Tissue Viscoelasticity

Ultrasound shear wave elastography is a promising noninvasive, low cost, and clinically viable tool for liver fibrosis staging. Current shear wave imaging technologies on clinical ultrasound scanners ignore shear wave dispersion and use a single group velocity measured over the shear wave bandwidth to estimate tissue elasticity. The center frequency and bandwidth of shear waves induced by acoustic radiation force depend on the ultrasound push beam (push duration, -number, etc.) and the viscoelasticity of the medium, and therefore are different across scanners from different vendors. As a result, scanners from different vendors may give different tissue elasticity measurements within the same patient. Various methods have been proposed to evaluate shear wave dispersion to better estimate tissue viscoelasticity. A rheological model such as the Kelvin-Voigt model is typically fitted to the shear wave dispersion to solve for the elasticity and viscosity of tissue. However, these rheological models impose strong assumptions about frequency dependence of elasticity and viscosity. Here, we propose a new method called Acoustic Radiation Force Induced Creep-Recovery (ARFICR) capable of quantifying rheological model-independent measurements of elasticity and viscosity for more robust tissue health assessment. In ARFICR, the creep-recovery time signal at the focus of the push beam is used to calculate the relative elasticity and viscosity (scaled by an unknown constant) over a wide frequency range. Shear waves generated during the ARFICR measurement are also detected and used to calculate the shear wave velocity at its center frequency, which is then used to calibrate the relative elasticity and viscosity to absolute elasticity and viscosity. In this paper, finite-element method simulations and experiments in tissue mimicking phantoms are used to validate and characterize the extent of viscoelastic quantification of ARFICR. The results suggest that ARFICR can measure tissue viscoelasticity reliably. Moreover, the results showed the strong frequency dependence of viscoelastic parameters in tissue mimicking phantoms and healthy liver.

Viscoelasticity Mapping by Identification of Local Shear Wave Dynamics

Estimation of soft tissue elasticity is of interest in several clinical applications. For instance, tumors and fibrotic lesions are notoriously stiff compared with benign tissue. A fully quantitative measure of lesion stiffness can be obtained by shear wave (SW) elastography. This method uses an acoustic radiation force to produce laterally propagating SWs that can be tracked to obtain the velocity, which in turn is related to Young's modulus. However, not only elasticity, but also viscosity plays an important role in the propagation process of SWs. In fact, viscosity itself is a parameter of diagnostic value for the detection and characterization of malignant lesions. In this paper, we describe a new method that enables imaging viscosity from SW elastography by local model-based system identification. By testing the method on simulated data sets and performing in vitro experiments, we show that the ability of the proposed technique to generate parametric maps of the viscoelastic material properties from SW measurements, opening up new possibilities for noninvasive tissue characterization.

Noninvasive evaluation of corneal viscoelasticity based on displacement in response to acoustic radiation force

Noninvasive assessment of corneal mechanical properties in vivo will help to the understanding of the pathogenesis and early diagnosis of ectatic corneal disorders. This study presented a noninvasive method that assesses the corneal biomechanical properties by exciting the cornea with acoustic radiation force and monitoring its displacement using a dual-frequency confocal transducer. A 3.85-MHz pushing element was used to induce localized tissue displacement, and the displacement was detected by the 12.8-MHz detecting elements using pulse-echo methods. Under constant acoustic radiation force, 0the tissue displacement are directly correlated with tissue biomechanical properties, a set of parameters were extracted from local displacement waveform, including relaxation time constant (τ), relative elasticity (RE) and relative viscosity (RV). In order to obtain corneal samples with different mechanical properties, the fresh bovine eyes were performed by collagen cross-linking (CXL). The result indicated that the estimated τ in untreated corneas were statistically significantly different (p<;0.05) from those of treated corneas and the estimation of τ varied significantly with different degrees of CXL. It is possible to develop a non-invasive, effective and efficient method with high spatial resolution.

Efficient 3-D Reconstruction in Ultrasound Elastography via a Sparse Iteration Based on Markov Random Fields

Percutaneous needle-based liver ablation procedures are becoming increasingly common for the treatment of small isolated tumors in hepatocellular carcinoma patients who are not candidates for surgery. Rapid 3-D visualization of liver ablations has potential clinical value, because it can enable interventional radiologists to plan and execute needle-based ablation procedures with real time feedback. Ensuring the right volume of tissue is ablated is desirable to avoid recurrence of tumors from residual untreated cancerous cells. Shear wave velocity (SWV) measurements can be used as a surrogate for tissue stiffness to distinguish stiffer ablated regions from softer untreated tissue. This paper extends the previously reported sheaf reconstruction method to generate complete 3-D visualizations of SWVs without resorting to an approximate intermediate step of reconstructing transverse C planes. The noisy data are modeled using a Markov random field, and a computationally tractable reconstruction algorithm that can handle grids with millions of points is developed. Results from simulated ellipsoidal inclusion data show that this algorithm outperforms standard nearest neighbor interpolation by an order of magnitude in mean squared reconstruction error. Results from the phantom experiments show that it also provides a higher contrast-to-noise ratio by almost 2 dB and better signal-to-noise ratio in the stiff inclusion by over 2 dB compared with nearest neighbor interpolation and has lower computational complexity than linear and spline interpolation.

On the Quantitative Potential of Viscoelastic Response (VisR) Ultrasound Using the One-Dimensional Mass-Spring-Damper Model

Viscoelastic response (VisR) ultrasound is an acoustic radiation force (ARF)-based imaging method that fits induced displacements to a one-dimensional (1-D) mass-spring-damper (MSD) model to estimate the ratio of viscous to elastic moduli, τ, in viscoelastic materials. Error in VisR τ estimation arises from inertia and acoustic displacement underestimation. These error sources are herein evaluated using finite-element method (FEM) simulations, error correction methods are developed, and corrected VisR τ estimates are compared with true simulated τ values to assess VisR's relevance to quantifying viscoelasticity. With regard to inertia, adding a mass term in series with the Voigt model, to achieve the MSD model, accounts for inertia due to tissue mass when ideal point force excitations are used. However, when volumetric ARF excitations are applied, the induced complex system inertia is not described by the single-degree-of-freedom MSD model, causing VisR to overestimate τ. Regarding acoustic displacement underestimation, associated deformation of ARF-induced displacement profiles further distorts VisR τ estimates. However, median error in VisR τ is reduced to approximately -10% using empirically derived error correction functions applied to simulated viscoelastic materials with viscous and elastic properties representative of tissue. The feasibility of corrected VisR imaging is then demonstrated in vivo in the rectus femoris muscle of an adult with no known neuromuscular disorders. These results suggest VisR's potential relevance to quantifying viscoelastic properties clinically.