Distributed Aberration Correction Techniques Based on Tomographic Sound Speed Estimates

Real-Time Speed-of-Sound Estimation In Vivo via Steered Plane Wave Ultrasound

Speed-of-sound (SoS) is an intrinsic acoustic property of human tissues and has been regarded as a potential biomarker of tissue health. To foster the clinical use of this emerging biomarker in medical diagnostics, it is important for SoS estimates to be derived and displayed in real time. Here, we demonstrate that concurrent global SoS estimation and B-mode imaging can be achieved live on a portable ultrasound scanner. Our innovation is hinged upon the design of a novel pulse-echo SoS estimation framework that is based on steered plane wave imaging. It has accounted for the effects of refraction and imaging depth when the medium SoS differs from the nominal value of 1540 m/s that is conventionally used in medical imaging. The accuracy of our SoS estimation framework was comparatively analyzed with through-transmit time-of-flight measurements in vitro on 15 custom agar phantoms with different SoS values (1508-1682 m/s) and in vivo on human calf muscles ( N = 9 ; SoS range: 1560-1586 m/s). Our SoS estimation framework has a mean signed difference (MSD) of - 0.6 ± 2.3 m/s in vitro and - 2.2 ± 11.2 m/s in vivo relative to the reference measurements. In addition, our real-time system prototype has yielded simultaneous SoS estimates and B-mode imaging at an average frame rate of 18.1 fps. Overall, by realizing real-time tissue SoS estimation with B-mode imaging, our innovation can foster the use of tissue SoS as a biomarker in medical ultrasound diagnostics.

Improved Transcranial Plane-Wave Imaging With Learned Speed-of-Sound Maps

Although transcranial ultrasound plane-wave imaging (PWI) has promising clinical application prospects, studies have shown that variable speed-of-sound (SoS) would seriously damage the quality of ultrasound images. The mismatch between the conventional constant velocity assumption and the actual SoS distribution leads to the general blurring of ultrasound images. The optimization scheme for reconstructing transcranial ultrasound image is often solved using iterative methods like full-waveform inversion. These iterative methods are computationally expensive and based on prior magnetic resonance imaging (MRI) or computed tomography (CT) information. In contrast, the multi-stencils fast marching (MSFM) method can produce accurate time travel maps for the skull with heterogeneous acoustic speed. In this study, we first propose a convolutional neural network (CNN) to predict SoS maps of the skull from PWI channel data. Then, use these maps to correct the travel time to reduce transcranial aberration. To validate the performance of the proposed method, numerical, phantom and intact human skull studies were conducted using a linear array transducer (L11-5v, 128 elements, pitch = 0.3 mm). Numerical simulations demonstrate that for point targets, the lateral resolution of MSFM-restored images increased by 65%, and the center position shift decreased by 89%. For the cyst targets, the eccentricity of the fitting ellipse decreased by 75%, and the center position shift decreased by 58%. In the phantom study, the lateral resolution of MSFM-restored images was increased by 49%, and the position shift was reduced by 1.72 mm. This pipeline, termed AutoSoS, thus shows the potential to correct distortions in real-time transcranial ultrasound imaging, as demonstrated by experiments on the intact human skull.

Investigation of a method to estimate the average speed of sound using phase variances of element signals for ultrasound compound imaging

PURPOSE

In the receive beamforming of an ultrasonography system, a B-mode image is reconstructed by assuming an average speed of sound (SoS) as a constant value. In our previous studies, we proposed a method for estimating the average SoS based on the coherence factor (CF) and the reciprocal of phase variances of element signals in delay-and-sum (DAS) beamforming. In this paper, we investigate the accuracy of estimation of the average SoS for compound imaging.

METHODS

For this purpose, two numerical simulations were performed with k-Wave software. Also, the estimation methods based on the CF and the reciprocal were applied to in vivo data from the common carotid artery, and B-mode images were reconstructed using the estimated average SoS.

RESULTS

In the first numerical simulation using an inhomogeneous phantom, the relationship between the accuracy and the transmission angles for the estimation was investigated, and the root mean squared errors (RMSEs) of estimates obtained based on the CF and the reciprocal of the phase variance were 1.25 ± 0.09, and 0.765 ± 0.17% at the transmission sequence of steering angles of (- 10°, - 5°, 0°, 5°, 10°), respectively. In the second numerical simulation using a cyst phantom, lateral resolutions were improved by reconstructing the image using the estimates obtained using the proposed strategy (reciprocal). By the proposed strategy, improvement of the continuity of the lumen-intima interface in the lateral direction was observed in the in vivo experiment.

CONCLUSION

Consequently, the results indicated that the proposed strategy was beneficial for estimation of the average SoS and image reconstruction.

Aberration correction in diagnostic ultrasound: A review of the prior field and current directions

Medical ultrasound images are reconstructed with simplifying assumptions on wave propagation, with one of the most prominent assumptions being that the imaging medium is composed of a constant sound speed. When the assumption of a constant sound speed are violated, which is true in most in vivoor clinical imaging scenarios, distortion of the transmitted and received ultrasound wavefronts appear and degrade the image quality. This distortion is known as aberration, and the techniques used to correct for the distortion are known as aberration correction techniques. Several models have been proposed to understand and correct for aberration. In this review paper, aberration and aberration correction are explored from the early models and correction techniques, including the near-field phase screen model and its associated correction techniques such as nearest-neighbor cross-correlation, to more recent models and correction techniques that incorporate spatially varying aberration and diffractive effects, such as models and techniques that rely on the estimation of the sound speed distribution in the imaging medium. In addition to historical models, future directions of ultrasound aberration correction are proposed.

A joint method of coherence factor and nonlinear beamforming for synthetic aperture imaging with a ring array

Ultrasound computed tomography (USCT) with a ring array is an emerging diagnostic method for breast cancer. In the literature, synthetic aperture (SA) imaging has employed the delay-and-sum (DAS) beamforming technique for ring-array USCT to obtain isotropic resolution reflection images. However, the images obtained by the conventional DAS beamformer suffer from off-axis clutter and low resolution due to inhomogeneity of the medium and phase distortion. To address these issues, researchers have developed adaptive beamforming methods, such as coherence factor (CF) and convolutional beamforming algorithm (COBA), that improve image quality. In this study, we propose a joint method that combines CF with short-lag COBA (SLCOBA). First, we estimate the average sound speed using CF to address tissue inhomogeneity. Based on the corrected sound speed map, SLCOBA effectively suppresses side lobes and enhances image quality. Numerical results show that the proposed method reduces clutter and noise, improving resolution performance. These findings suggest that the proposed method may be a practical option for breast imaging in inhomogeneous mediums in the future.

Fast Marching based Tissue Adaptive Delay Estimation for Aberration Corrected Delay and Sum Beamforming in Ultrasound Imaging

Conventional ultrasound (US) imaging employs the delay and sum (DAS) receive beamforming with dynamic receive focus for image reconstruction due to its simplicity and robustness. However, the DAS beamforming follows a geometrical method of delay estimation with a spatially constant speed-of-sound (SoS) of 1540 m/s throughout the medium irrespective of the tissue in-homogeneity. This approximation leads to errors in delay estimations that accumulate with depth and degrades the resolution, contrast and overall accuracy of the US image. In this work, we propose a fast marching based DAS for focused transmissions which leverages the approximate SoS map to estimate the refraction corrected propagation delays for each pixel in the medium. The proposed approach is validated qualitatively and quantitatively for imaging depths of upto ∼ 11 cm through simulations, where fat layer-induced aberration is employed to alter the SoS in the medium. To the best of the authors' knowledge, this is the first work considering the effect of SoS on image quality for deeper imaging.Clinical relevance- The proposed approach when employed with an approximate SoS estimation technique can aid in overcoming the fat-induced signal aberrations and thereby in the accurate imaging of various pathologies of liver and abdomen.

Excluding Echo Shift Noise in Real-Time Pulse-Echo Speed-of-Sound Imaging

Computed ultrasound tomography in echo mode (CUTE) allows real-time imaging of the tissue speed of sound (SoS) using handheld ultrasound. The SoS is retrieved by inverting a forward model that relates the spatial distribution of the tissue SoS to echo shift maps detected between varying transmit and receive angles. Despite promising results, in vivo SoS maps often show artifacts due to elevated noise in echo shift maps. To minimize artifacts, we propose a technique where an individual SoS map is reconstructed for each echo shift map separately, as opposed to a single SoS map from all echo shift maps simultaneously. The final SoS map is then obtained as a weighted average over all SoS maps. Due to the partial redundancy between different angle combinations, artifacts that appear only in a subset of the individual maps can be excluded via the averaging weights. We investigate this real-time capable technique in simulations using two numerical phantoms, one with a circular inclusion and one with two layers. Our results demonstrate that the SoS maps reconstructed using the proposed technique are equivalent to the ones using simultaneous reconstruction when considering uncorrupted data but show significantly reduced artifact level for data that are corrupted by noise.

An aberration correction approach for single and dual aperture ultrasound imaging of the abdomen

Abdominal ultrasound image quality is hampered by phase aberration, that is mainly caused by the large speed-of-sound (SoS) differences between fat and muscle tissue in the abdominal wall. The mismatch between the assumed and actual SoS distribution introduces general blurring of the ultrasound images, and acoustic refraction can lead to geometric distortion of the imaged features. Large aperture imaging or dual-transducer imaging can improve abdominal imaging at deep locations by providing increased contrast and resolution. However, aberration effects for large aperture imaging can be even more severe, which limits its full potential. In this study, a model-based aberration correction method for arbitrary acquisition schemes is introduced for delay-and-sum (DAS) beamforming and its performance was analyzed for both single- and dual-transducer ultrasound imaging. The method employs aberration corrected wavefront arrival times, using manually assigned local SoS values. Two wavefront models were compared. The first model is based on a straight ray approximation, and the second model on the Eikonal equation, which is solved by a multi-stencils fast marching method. Their accuracy for abdominal imaging was evaluated in acoustic simulations and phantom experiments involving tissue-mimicking and porcine material with large SoS contrast (∼100 m/s). The lateral resolution was improved by up to 90% in simulations and up to 65% in experiments compared to standard DAS, in which the use of Eikonal beamforming generally outperformed straight ray beamforming. Moreover, geometric distortions were mitigated in multi-aperture imaging, leading to a reduction in position error of around 80%. A study on the sensitivity of the aberration correction to shape and SoS of aberrating layers was performed, showing that even with imperfect segmentations or SoS values, aberration correction still outperforms standard DAS.

Large-Array Deep Abdominal Imaging in Fundamental and Harmonic Mode

Deep abdominal images suffer from poor diffraction-limited lateral resolution. Extending the aperture size can improve resolution. However, phase distortion and clutter can limit the benefits of larger arrays. Previous studies have explored these effects using numerical simulations, multiple transducers, and mechanically swept arrays. In this work, we used an 8.8-cm linear array transducer to investigate the effects of aperture size when imaging through the abdominal wall. We acquired channel data in fundamental and harmonic modes using five aperture sizes. To avoid motion and increase the parameter sampling, we decoded the full-synthetic aperture data and retrospectively synthesized nine apertures (2.9-8.8 cm). We imaged a wire target and a phantom through ex vivo porcine abdominal samples and scanned the livers of 13 healthy subjects. We applied bulk sound speed correction to the wire target data. Although point resolution improved from 2.12 to 0.74 mm at 10.5 cm depth, contrast resolution often degraded with aperture size. In subjects, larger apertures resulted in an average maximum contrast degradation of 5.5 dB at 9-11 cm depth. However, larger apertures often led to visual detection of vascular targets unseen with conventional apertures. An average 3.7-dB contrast improvement over fundamental mode in subjects showed that the known benefits of tissue-harmonic imaging extend to larger arrays.

Sound Speed Estimation for Distributed Aberration Correction in Laterally Varying Media

Spatial variation in sound speed causes aberration in medical ultrasound imaging. Although our previous work has examined aberration correction in the presence of a spatially varying sound speed, practical implementations were limited to layered media due to the sound speed estimation process involved. Unfortunately, most models of layered media do not capture the lateral variations in sound speed that have the greatest aberrative effect on the image. Building upon a Fourier split-step migration technique from geophysics, this work introduces an iterative sound speed estimation and distributed aberration correction technique that can model and correct for aberrations resulting from laterally varying media. We first characterize our approach in simulations where the scattering in the media is known a-priori. Phantom and in-vivo experiments further demonstrate the capabilities of the iterative correction technique. As a result of the iterative correction scheme, point target resolution improves by up to a factor of 4 and lesion contrast improves by up to 10.0 dB in the phantom experiments presented.

Ultrasound Matrix Imaging-Part I: The Focused Reflection Matrix, the F-Factor and the Role of Multiple Scattering

This is the first article in a series of two dealing with a matrix approach for aberration quantification and correction in ultrasound imaging. Advanced synthetic beamforming relies on a double focusing operation at transmission and reception on each point of the medium. Ultrasound matrix imaging (UMI) consists in decoupling the location of these transmitted and received focal spots. The response between those virtual transducers form the so-called focused reflection matrix that actually contains much more information than a confocal ultrasound image. In this paper, a time-frequency analysis of this matrix is performed, which highlights the single and multiple scattering contributions as well as the impact of aberrations in the monochromatic and broadband regimes. Interestingly, this analysis enables the measurement of the incoherent input-output point spread function at any pixel of this image. A fitting process enables the quantification of the single scattering, multiple scattering and noise components in the image. From the single scattering contribution, a focusing criterion is defined, and its evolution used to quantify the amount of aberration throughout the ultrasound image. In contrast to the state-of-the-art coherence factor, this new indicator is robust to multiple scattering and electronic noise, thereby providing a contrasted map of the focusing quality at a much better transverse resolution. After a validation of the proof-of-concept based on time-domain simulations, UMI is applied to the in-vivo study of a human calf. Beyond this specific example, UMI opens a new route for speed-of-sound and scattering quantification in ultrasound imaging.