Simultaneous identification of elastic properties, thickness, and diameter of arteries excited with ultrasound radiation force

Full waveform inversion for arterial viscoelasticity

Objective. Arterial viscosity is emerging as an important biomarker, in addition to the widely used arterial elasticity. This paper presents an approach to estimate arterial viscoelasticity using shear wave elastography (SWE).Approach. While dispersion characteristics are often used to estimate elasticity from SWE data, they are not sufficiently sensitive to viscosity. Driven by this, we develop a full waveform inversion (FWI) methodology, based on directly matching predicted and measured wall velocity in space and time, to simultaneously estimate both elasticity and viscosity. Specifically, we propose to minimize an objective function capturing the correlation between measured and predicted responses of the anterior wall of the artery.Results. The objective function is shown to be well-behaving (generally convex), leading us to effectively use gradient optimization to invert for both elasticity and viscosity. The resulting methodology is verified with synthetic data polluted with noise, leading to the conclusion that the proposed FWI is effective in estimating arterial viscoelasticity.Significance. Accurate estimation of arterial viscoelasticity, not just elasticity, provides a more precise characterization of arterial mechanical properties, potentially leading to a better indicator of arterial health.

The influence of acoustic radiation force beam shape and location on wave spectral content for arterial dispersion ultrasound vibrometry

Objective. Arterial dispersion ultrasound vibrometry (ADUV) relies on the use of guided waves in arterial geometries for shear wave elastography measurements. Both the generation of waves through the use of acoustic radiation force (ARF) and the techniques employed to infer the speed of the resulting wave motion affect the spectral content and accuracy of the measurement. In particular, the effects of the shape and location of the ARF beam in ADUV have not been widely studied. In this work, we investigated how such variations of the ARF beam affect the induced motion and the measurements in the dispersive modes that are excited. Approach. The study includes an experimental evaluation on an arterial phantom and an in vivo validation of the observed trends, observing the two walls of the waveguide, simultaneously, when subjected to variations in the ARF beam extension (F/N) and focus location. Main results. Relying on the theory of guided waves in cylindrical shells, the shape of the beam controls the selection and nature of the induced modes, while the location affects the measured dispersion curves (i.e. variation of phase velocity with frequency or wavenumber, multiple modes) across the waveguide walls. Significance. This investigation is important to understand the spectral content variations in ADUV measurements and to maximize inversion accuracy by tuning the ARF beam settings in clinical applications.

Multimodal guided wave inversion for arterial stiffness: methodology and validation in phantoms

Arterial stiffness is an important biomarker for many cardiovascular diseases. Shear wave elastography is a recent technique aimed at estimating local arterial stiffness using guided wave inversion (GWI), i.e. matching the computed and measured wave dispersion. This paper develops and validates a new GWI approach by synthesizing various recent observations and algorithms: (a) refinements to signal processing to obtain more accurate experimental dispersion curves; (b) an efficient forward model to compute theoretical dispersion curves for immersed, incompressible cylindrical waveguides; (c) an optimization framework based on the recent observation that the measured dispersion curve is multimodal, i.e. it matches for not one but two different wave modes in two different frequency ranges. The resulting inversion approach is validated using extensive experimental data from rubber tube phantoms, not only for modulus estimation but also to simultaneously estimate modulus and wall thickness. The observations indicate that the modulus estimates are best performed with the information on wall thickness. The approach, which takes less than half a minute to run, is shown to be accurate, with the modulus estimated with less than 4% error for 70% of the experiments.

Tissue Characterization by Low-Frequency Acoustic Waves Generated by a Single High-Frequency Focused Ultrasound Beam

The mechanical properties of biological tissues are fingerprints of certain pathologic processes. Ultrasound systems have been used as a non-invasive technique to both induce kilohertz-frequency mechanical vibrations and detect waves resulting from interactions with biological structures. However, existing methodologies to produce kilohertz-frequency mechanical vibrations using ultrasound require the use of variable-frequency, dual-frequency and high-power systems. Here, we propose and demonstrate the use of bursts of megahertz- frequency acoustic radiation to observe kilohertz-frequency mechanical responses in biological tissues. Femoral bones were obtained from 10 healthy mice and 10 mice in which osteoporosis had been induced. The bones' porosity, trabecular number, trabecular spacing, connectivity and connectivity density were determined using micro-computed tomography (μCT). The samples were irradiated with short, focused acoustic radiation pulses (f = 3.1 MHz, t = 15 μs), and the low-frequency acoustic response (1-100 kHz) was acquired using a dedicated hydrophone. A strong correlation between the spectral maps of the acquired signals and the μCT data was found. In a subsequent evaluation, soft tissue stiffness measurements were performed with a gel wax-based tissue-mimicking phantom containing three spherical inclusions of the same type of gel but different densities and Young's moduli, yet with approximately the same echogenicity. Conventional B-mode ultrasound was unable to image the inclusions, while the novel technique proposed here showed good image contrast.

Experimental Investigation of Guided Wave Imaging in Thin Soft Media under Various Coupling Conditions

Guided wave imaging for the artery remains in its infancy in clinical practice mainly because of complex arterial microstructure, hemodynamics and boundary conditions. Despite the theoretically known potential effect of the surrounding medium on guided wave propagation in thin media in non-destructive testing, experimental evidence pertaining to thin soft materials, such as the artery, is relatively scarce in the relevant literature. Therefore, this study first evaluated the propagating guided wave generated by acoustic radiation force in polyvinyl alcohol-based hydrogel plates differing in thickness and stiffness under various material coupling conditions (water and polyvinyl alcohol bulk). A thin-walled polyvinyl alcohol hollow cylindrical phantom coupled by softer gelatin-agar phantoms and an excised porcine aorta surrounded by water and pork belly were further examined. Guided waves in the thin structure and shear waves in the bulk media were captured by ultrafast ultrasound imaging, and guided wave dispersion as a function of the frequency-thickness product was analyzed using the zero-order anti-symmetric Lamb wave model to estimate the shear modulus of each thin medium studied. Results confirmed the deviated shear modulus estimates from the ground truth for thin plates, the thin-walled hollow cylindrical phantom and the porcine aorta bounded by stiffness-unmatched bulk medium. The findings indicated the need for (i) careful interpretation of estimated shear moduli of thin structure bounded by bulk media and (ii) a generalized guided wave model that takes into account the effect of coupling medium.

Dispersive Radon transform

Dispersion results in the spreading and overlapping of the wave-packets, which often limits the capability of signal interpretation; on the other hand, such a phenomenon can also be used for structure or media evaluation. In this study, the authors propose an original dispersive Radon transform (DRT), which is formulated as integration transform along a set of dispersion curves. Multichannel dispersive signals of each individual mode can be concentrated to a well localized region in the DRT domain. The proposed DRT establishes the sparse projection of the dispersive components and provides an efficient solution for mode separation, noise filtering, and missing data reconstruction. Particularly the DRT method allows projecting the temporal signals of dispersive waves on the space of parameters of interest, which can be used to solve the inverse problem for waveguide or media property estimation. The least-square procedure and sparse scheme of the DRT are introduced. A high-resolution DRT is designed based on an iterative reweighting inversion scheme, which resembles the infinite-aperture velocity gather. The proposed method is applied by analyzing ultrasonic guided waves in plate-like structures and in a human radius specimen. The results suggest that the DRT method can significantly enhance the interpretation of dispersive signals.