Correlation-based full-waveform shear wave elastography.

Objective. With the ultimate goal of reconstructing 3D elasticity maps from ultrasound particle velocity measurements in a plane, we present in this paper a methodology of inverting for 2D elasticity maps from measurements on a single line.Approach. The inversion approach is based on gradient optimization where the elasticity map is iteratively modified until a good match is obtained between simulated and measured responses. Full-wave simulation is used as the underlying forward model to accurately capture the physics of shear wave propagation and scattering in heterogeneous soft tissue. A key aspect of the proposed inversion approach is a cost functional based on correlation between measured and simulated responses.Main results. We illustrate that the correlation-based functional has better convexity and convergence properties compared to the traditional least-squares functional, and is less sensitive to initial guess, robust against noisy measurements and other errors that are common in ultrasound elastography. Inversion with synthetic data illustrates the effectiveness of the method to characterize homogeneous inclusions as well as elasticity map of the entire region of interest.Significance. The proposed ideas lead to a new framework for shear wave elastography that shows promise in obtaining accurate maps of shear modulus using shear wave elastography data obtained from standard clinical scanners.

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.

Measurement of wave propagation through a tube using dual transducers for elastography in arteries.

Objective.Measuring waves induced with acoustic radiation force (ARF) in arteries has been studied over the last decade. To date, it remains a challenge to quantitatively assess the local arterial biomechanical properties. The cylindrical shape and waveguide behavior of waves propagating in the arterial wall pose complexities to determining the mechanical properties of the artery.Approach. In this paper, an artery-mimicking tube in water is examined utilizing three-dimensional measurements. The cross-section of the tube is measured while a transducer is translated over 41 different positions along the length of the tube. Motion in the radial direction is calculated using two components of motion which are measured from the two orthogonal views of the cross-section. This enables more accurate estimation of motion along the circumference of tube.Main results. The results provide more information to categorize the motion in tube wall into two types of responses: a transient response and a steady state response. The transient response is caused by ARF application and the waves travel along the length of the tube for a relatively short period of time. This corresponds to the axial and circumferential propagating waves. The two circumferential waves travel along the circumference of tube in CW (clockwise) and CCW (counter-clockwise) direction and result in a standing wave. By using a directional filter, the two waves were successfully separated, and their propagation was more clearly visualized. As a steady state response, a circumferential mode is generated showing a symmetric motion (i.e. the proximal and distal walls move in the opposite direction) following the transient response.Significance.This study provides a more comprehensive understanding of the waves produced in an artery-mimicking tube with ARF application, which will provide opportunities for improving measurement of arterial mechanical properties.

Full wave simulation of arterial response under acoustic radiation force.

With the ultimate goal of estimating arterial viscoelasticity using shear wave elastography, this paper presents a practical methodology to simulate the response of a human carotid artery under acoustic radiation force (ARF). The artery is idealized as a nearly incompressible viscoelastic hollow cylinder submerged in incompressible, inviscid fluid. For this idealization, we develop a multi-step methodology for efficient computation of three-dimensional response under complex ARF excitation, while capturing the fluid-structure interaction between the arterial wall and the surrounding fluid. The specific steps include (a) performing dimensional reduction through semi-analytical finite element formulation, (b) efficient finite element discretization using traditional and recent techniques. The computational efficiency is further enhanced by utilizing (c) modal superposition, followed by, where appropriate, (d) impulse response function. In addition to developing the methodology, convergence analysis is performed for a typical arterial geometry, leading to recommendations on various discretization parameters. At the end, the computational effort is shown to be several orders of magnitude less than the traditional, fully three-dimensional analysis using finite element methods, leading to a practical yet accurate simulation of arterial response under ARF excitations.

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 anin vivovalidation 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.

Toward improved accuracy in shear wave elastography of arteries through controlling the arterial response to ultrasound perturbation in-silico and in phantoms.

Dispersion-based inversion has been proposed as a viable direction for materials characterization of arteries, allowing clinicians to better study cardiovascular conditions using shear wave elastography. However, these methods rely ona prioriknowledge of the vibrational modes dominating the propagating waves induced by acoustic radiation force excitation: differences between anticipated and real modal content are known to yield errors in the inversion. We seek to improve the accuracy of this process by modeling the artery as a fluid-immersed cylindrical waveguide and building an analytical framework to prescribe radiation force excitations that will selectively excite certain waveguide modes using ultrasound acoustic radiation force. We show that all even-numbered waveguide modes can be eliminated from the arterial response to perturbation, and confirm the efficacy of this approach within silicotests that show that odd modes are preferentially excited. Finally, by analyzing data from phantom tests, we find a set of ultrasound focal parameters that demonstrate the viability of inducing the desired odd-mode response in experiments.

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.

Shear wave dispersion analysis of incompressible waveguides.

A suite of methodologies is presented to compute shear wave dispersion in incompressible waveguides encountered in biomedical imaging; plate, tube, and general prismatic waveguides, all immersed in an incompressible fluid, are considered in this effort. The developed approaches are based on semi-analytical finite element methods in the frequency domain with a specific focus on the complexities associated with the incompressibility of the solid media as well as the simplification facilitated by the incompressibility of the surrounding fluid. The proposed techniques use the traditional idea of selective reduced integration for the solid medium and the more recent idea of perfectly matched discrete layers for the surrounding fluid. Also, used is the recently developed complex-length finite element method for platelike structures. Several numerical examples are presented to illustrate the practicality and effectiveness of the developed techniques in computing shear wave dispersion in a variety of waveguides.

Characterizing blood clots using acoustic radiation force optical coherence elastography and ultrasound shear wave elastography.

Thromboembolism in a cerebral blood vessel is associated with high morbidity and mortality. Mechanical thrombectomy (MT) is one of the emergenc proceduresperformed to remove emboli. However, the interventional approaches such as aspiration catheters or stent retriever are empirically selected. An inappropriate selection of surgical devices can influence the success rate during embolectomy, which can lead to an increase in brain damage. There has been growing interest in the study of clot composition and using a priori knowledge of clot composition to provide guidance for an appropriate treatment strategy for interventional physicians. Developing imaging tools which can allow interventionalists to understand clot composition could affect management and device strategy. In this study, we investigated how clots of different compositions can be characterized by using acoustic radiation force optical coherence elastography (ARF-OCE) and compared with ultrasound shear wave elastography (SWE). Five different clots compositions using human blood were fabricated into cylindrical forms from fibrin-rich (21% red blood cells, RBCs) to RBC-rich (95% RBCs). Using the ARF-OCE and SWE, we characterized the wave velocities measured in the time-domain. In addition, the semi-analytical finite element model was used to explore the relationship between the phase velocities with various frequency ranges and diameters of the clots. The study demonstrated that the wave group velocities generally decrease as RBC content increases in ARF-OCE and SWE. The correlation of the group velocities from the OCE and SWE methods represented a good agreement as RBC composition is larger than 39%. Using the phase velocity dispersion analysis applied to ARF-OCE data, we estimated the shear wave velocities decoupling the effects of the geometry and material properties of the clots. The study demonstrated that the composition of the clots can be characterized by elastographic methods using ARF-OCE and SWE, and OCE demonstrated better ability to discriminate between clots of different RBC compositions, compared to the ultrasound-based approach, especially in clots with low RBC compositions.

Arterial waveguide model for shear wave elastography: implementation and in vitro validation.

Arterial stiffness is found to be an early indicator of many cardiovascular diseases. Among various techniques, shear wave elastography has emerged as a promising tool for estimating local arterial stiffness through the observed dispersion of guided waves. In this paper, we develop efficient models for the computational simulation of guided wave dispersion in arterial walls. The models are capable of considering fluid-loaded tubes, immersed in fluid or embedded in a solid, which are encountered in in vitro/ex vivo, and in vivo experiments. The proposed methods are based on judiciously combining Fourier transformation and finite element discretization, leading to a significant reduction in computational cost while fully capturing complex 3D wave propagation. The developed methods are implemented in open-source code, and verified by comparing them with significantly more expensive, fully 3D finite element models. We also validate the models using the shear wave elastography of tissue-mimicking phantoms. The computational efficiency of the developed methods indicates the possibility of being able to estimate arterial stiffness in real time, which would be beneficial in clinical settings.

A high resolution computer model for sound propagation in the human thorax based on the Visible Human data set.

A parallel supercomputer model based on realistic tissue data is developed for sound propagation in the human thorax and the sound propagation behavior is analyzed under various conditions using artificial sound sources. The model uses the Visible Human male data set for a realistic representation of the human thorax. The results were analyzed in time and frequency domains. The analysis suggests that lower frequencies of around 100 Hz are more effectively transmitted through the thorax and that the spatial confinement of sound waves within the thorax results in a resonance effect at around 1500 Hz. The results confirm previous studies that show the size of the thorax plays a significant role in the type of sound generated at the chest wall.