A multiresolution approach to shear wave image reconstruction

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.

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.

A parametric evaluation of shear wave speeds estimated with time-of-flight calculations in viscoelastic media

Shear wave elasticity imaging (SWEI) uses an acoustic radiation force to generate shear waves, and then soft tissue mechanical properties are obtained by analyzing the shear wave data. In SWEI, the shear wave speed is often estimated with time-of-flight (TOF) calculations. To characterize the errors produced by TOF calculations, three-dimensional (3D) simulated shear waves are described by time-domain Green's functions for a Kelvin-Voigt model evaluated for multiple combinations of the shear elasticity and the shear viscosity. Estimated shear wave speeds are obtained from cross correlations and time-to-peak (TTP) calculations applied to shear wave particle velocities and shear wave particle displacements. The results obtained from these 3D shear wave simulations indicate that TTP calculations applied to shear wave particle displacements yield effective estimates of the shear wave speed if noise is absent, but cross correlations applied to shear wave particle displacements are more robust when the effects of noise and shear viscosity are included. The results also show that shear wave speeds estimated with TTP methods and cross correlations using shear wave particle velocities are more sensitive to increases in shear viscosity and noise, which suggests that superior estimates of the shear wave speed are obtained from noiseless or noisy shear wave particle displacements.

A parametric evaluation of shear wave speeds estimated with time-of-flight calculations in viscoelastic media

Shear wave elasticity imaging (SWEI) uses an acoustic radiation force to generate shear waves, and then soft tissue mechanical properties are obtained by analyzing the shear wave data. In SWEI, the shear wave speed is often estimated with time-of-flight (TOF) calculations. To characterize the errors produced by TOF calculations, three-dimensional (3D) simulated shear waves are described by time-domain Green's functions for a Kelvin-Voigt model evaluated for multiple combinations of the shear elasticity and the shear viscosity. Estimated shear wave speeds are obtained from cross correlations and time-to-peak (TTP) calculations applied to shear wave particle velocities and shear wave particle displacements. The results obtained from these 3D shear wave simulations indicate that TTP calculations applied to shear wave particle displacements yield effective estimates of the shear wave speed if noise is absent, but cross correlations applied to shear wave particle displacements are more robust when the effects of noise and shear viscosity are included. The results also show that shear wave speeds estimated with TTP methods and cross correlations using shear wave particle velocities are more sensitive to increases in shear viscosity and noise, which suggests that superior estimates of the shear wave speed are obtained from noiseless or noisy shear wave particle displacements.

Spatial resolution in dynamic optical coherence elastography

Dynamic optical coherence elastography (OCE) tracks elastic wave propagation speed within tissue, enabling quantitative three-dimensional imaging of the elastic modulus. We show that propagating mechanical waves are mode converted at interfaces, creating a finite region on the order of an acoustic wavelength where there is not a simple one-to-one correspondence between wave speed and elastic modulus. Depending on the details of a boundary’s geometry and elasticity contrast, highly complex propagating fields produced near the boundary can substantially affect both the spatial resolution and contrast of the elasticity image. We demonstrate boundary effects on Rayleigh waves incident on a vertical boundary between media of different shear moduli. Lateral resolution is defined by the width of the transition zone between two media and is the limit at which a physical inclusion can be detected with full contrast. We experimentally demonstrate results using a spectral-domain OCT system on tissue-mimicking phantoms, which are replicated using numerical simulations. It is shown that the spatial resolution in dynamic OCE is determined by the temporal and spatial characteristics (i.e., bandwidth and spatial pulse width) of the propagating mechanical wave. Thus, mechanical resolution in dynamic OCE inherently differs from the optical resolution of the OCT imaging system.

Evaluation of Reconstruction Parameters for 2-D Comb-Push Ultrasound Shear Wave Elastography

Shear wave elastography (SWE) is a noninvasive ultrasound imaging modality used in the assessment of the mechanical properties of tissues such as the liver, kidney, skeletal muscle, thyroid, and the breast. Among the methods used to perform SWE is the comb-push ultrasound shear elastography method. This method uses multiple focused ultrasound beams to generate push beams with acoustic radiation force. Applying these push beams generates propagating shear waves. The propagation motion is measured with ultrafast ultrasound imaging. The shear wave motion data are directionally filtered, and a 2-D shear wave velocity (SWV) algorithm is applied to create group velocity maps. This algorithm uses a moving window and a specified patch for performing cross-correlations of time-domain signals. We performed a parametric study of how the choice of the patch and window size affected the reconstruction of the SWV in homogeneous and inclusion phantoms. We quantified the mean velocity and coefficient of variation in the homogeneous phantoms. We measured the contrast-to-noise ratio and bias in the inclusion phantoms. In each of these cases, we found that particular combinations of the patch and window provided optimal values of these evaluation metrics for the phantoms tested. This study provides a basis to construct algorithms to produce optimal SWV reconstructions for various clinical applications.