Local Sound Speed Estimation for Pulse-Echo Ultrasound in Layered Media

Spatial Ambiguity Correction in Coherence-Based Average Sound Speed Estimation

Sound speed estimation can potentially correct the focusing errors in medical ultrasound. Maximizing the echo spatial coherence as a function of beamforming sound speed is a known technique to estimate the average sound speed. However, beamformation with changing sound speed causes a spatial shift of the echo signals resulting in noise and registration errors in the average sound speed estimates. We show that the spatial shift can be predicted and corrected, leading to superior sound speed estimates. Methods are presented for axial and 2-D location correction. Methods were evaluated using simulations and experimental phantom data. The location correction strategies improved the variance of sound speed estimates and reduced artifacts in the presence of strong backscatter variations. Limitations of the proposed methods and potential improvement strategies were evaluated.

2-D Slicewise Waveform Inversion of Sound Speed and Acoustic Attenuation for Ring Array Ultrasound Tomography Based on a Block LU Solver

Ultrasound tomography is an emerging imaging modality that uses the transmission of ultrasound through tissue to reconstruct images of its mechanical properties. Initially, ray-based methods were used to reconstruct these images, but their inability to account for diffraction often resulted in poor resolution. Waveform inversion overcame this limitation, providing high-resolution images of the tissue. Most clinical implementations, often directed at breast cancer imaging, currently rely on a frequency-domain waveform inversion to reduce computation time. For ring arrays, ray tomography was long considered a necessary step prior to waveform inversion in order to avoid cycle skipping. However, in this paper, we demonstrate that frequency-domain waveform inversion can reliably reconstruct high-resolution images of sound speed and attenuation without relying on ray tomography to provide an initial model. We provide a detailed description of our frequency-domain waveform inversion algorithm with open-source code and data that we make publicly available.

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.

Windowed Radon Transform for Robust Speed-of-Sound Imaging With Pulse-Echo Ultrasound

In recent years, methods estimating the spatial distribution of tissue speed of sound with pulse-echo ultrasound are gaining considerable traction. They can address limitations of B-mode imaging, for instance in diagnosing fatty liver diseases. Current state-of-the-art methods relate the tissue speed of sound to local echo shifts computed between images that are beamformed using restricted transmit and receive apertures. However, the aperture limitation affects the robustness of phase-shift estimations and, consequently, the accuracy of reconstructed speed-of-sound maps. Here, we propose a method based on the Radon transform of image patches able to estimate local phase shifts from full-aperture images. We validate our technique on simulated, phantom and in-vivo data acquired on a liver and compare it with a state-of-the-art method. We show that the proposed method enhances the stability to changes of beamforming speed of sound and to a reduction of the number of insonifications. In particular, the deployment of pulse-echo speed-of-sound estimation methods onto portable ultrasound devices can be eased by the reduction of the number of insonifications allowed by the proposed method.

Investigating pulse-echo sound speed estimation in breast ultrasound with deep learning

Ultrasound is an adjunct tool to mammography that can quickly and safely aid physicians in diagnosing breast abnormalities. Clinical ultrasound often assumes a constant sound speed to form diagnostic B-mode images. However, the components of breast tissue, such as glandular tissue, fat, and lesions, differ in sound speed. Given a constant sound speed assumption, these differences can degrade the quality of reconstructed images via phase aberration. Sound speed images can be a powerful tool for improving image quality and identifying diseases if properly estimated. To this end, we propose a supervised deep-learning approach for sound speed estimation from analytic ultrasound signals. We develop a large-scale simulated ultrasound dataset that generates representative breast tissue samples by modeling breast gland, skin, and lesions with varying echogenicity and sound speed. We adopt a fully convolutional neural network architecture trained on a simulated dataset to produce an estimated sound speed map. The simulated tissue is interrogated with a plane wave transmit sequence, and the complex-value reconstructed images are used as input for the convolutional network. The network is trained on the sound speed distribution map of the simulated data, and the trained model can estimate sound speed given reconstructed pulse-echo signals. We further incorporate thermal noise augmentation during training to enhance model robustness to artifacts found in real ultrasound data. To highlight the ability of our model to provide accurate sound speed estimations, we evaluate it on simulated, phantom, and in-vivo breast ultrasound data.

Estimation error in speed of sound caused by rotation of measured cross-section from short-axis plane of blood vessels: a preliminary study

PURPOSE

Estimating the speed of sound (SoS) in ultrasound propagation media is important for improving the quality of B-mode images and for quantitative tissue characterization. We have been studying a method for estimating the SoS by measuring the reception time distribution of waves scattered from a scatterer at the elements in a probe. Previously, the measurement cross section was assumed to be perpendicular to the long axis of the blood vessel. In this study, we experimentally investigated the relationship between rotation angle [Formula: see text] of the probe relative to the short-axis plane of the blood vessel and the estimated SoS, [Formula: see text].

METHODS

Water tank and phantom experiments were conducted to investigate the characteristics of [Formula: see text] and element signals when the probe was rotated.

RESULTS

The received signal powers at the elements around both edges greatly decreased as [Formula: see text] increased. We introduced a parameter representing the decrease in power, [Formula: see text], in the received signal at the elements at both edges relative to the center element. [Formula: see text] was estimated to be larger as [Formula: see text] increased, especially for [Formula: see text]. [Formula: see text] also increased as [Formula: see text] increased. An approximately proportional relationship existed between the errors in [Formula: see text] and [Formula: see text].

CONCLUSION

Based on these results, we can distinguish between the presence and the absence of SoS misestimations using the difference in power among the elements in the received signal. In the absence of misestimation, we can obtain the true SoS, even if the target has a non-negligible size, by applying our previously proposed methods.

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.

Real-time assessment of high-intensity focused ultrasound heating and cavitation with hybrid optoacoustic ultrasound imaging

High-intensity focused ultrasound (HIFU) enables localized ablation of biological tissues by capitalizing on the synergistic effects of heating and cavitation. Monitoring of those effects is essential for improving the efficacy and safety of HIFU interventions. Herein, we suggest a hybrid optoacoustic-ultrasound (OPUS) approach for real-time assessment of heating and cavitation processes while providing an essential anatomical reference for accurate localization of the HIFU-induced lesion. Both effects could clearly be observed by exploiting the temperature dependence of optoacoustic (OA) signals and the strong contrast of gas bubbles in pulse-echo ultrasound (US) images. The differences in temperature increase and its rate, as recorded with a thermal camera for different HIFU pressures, evinced the onset of cavitation at the expected pressure threshold. The estimated temperatures based on OA signal variations were also within 10-20 % agreement with the camera readings for temperatures below the coagulation threshold (∼50 °C). Experiments performed in excised tissues as well as in a post-mortem mouse demonstrate that both heating and cavitation effects can be effectively visualized and tracked using the OPUS approach. The good sensitivity of the suggested method for HIFU monitoring purposes was manifested by a significant increase in contrast-to-noise ratio within the ablated region by > 10 dB and > 5 dB for the OA and US images, respectively. The hybrid OPUS-based monitoring approach offers the ease of handheld operation thus can readily be implemented in a bedside setting to benefit several types of HIFU treatments used in the clinics.

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.

Speed-of-sound estimation in ultrasound propagation medium by considering size of target scatterer

PURPOSE

Accurate speed-of-sound (SoS) estimation in an ultrasound propagation medium improves imaging quality and contributes to better diagnosis of diseases. In conventional time-delay-based SoS estimation approaches studied by several groups, a received wave is assumed to be scattered from an ideal point scatterer. In these approaches, the SoS is overestimated when the target scatterer has a non-negligible size. In this paper, we propose the SoS estimation method that considers target size.

METHODS

In the proposed method, the error ratio of the estimated SoS using the conventional time-delay-based approach is determined from measurable parameters using the geometric relationship between the received elements and target. Subsequently, the SoS erroneously estimated using conventional estimation, assuming the ideal point scatterer as a target, is corrected by the determined estimation error ratio. To validate the proposed method, the SoS in water was estimated for several wire sizes.

RESULTS

The SoS in the water was overestimated using the conventional SoS estimation method, with a maximum positive error of 38 m/s. The proposed method corrected the SoS estimates, and the errors were suppressed to within 6 m/s, irrespective of the wire diameter.

CONCLUSION

The present results demonstrate that the proposed method can estimate the SoS by considering the target size without using information on the true SoS, true target depth, and true target size, which is applicable to in vivo measurements.

Noninvasive estimation of local speed of sound by pulse-echo ultrasound in a rat model of nonalcoholic fatty liver

Objective. Speed of sound has previously been demonstrated to correlate with fat concentration in the liver. However, estimating speed of sound in the liver noninvasively can be biased by the speed of sound of the tissue layers overlying the liver. Here, we demonstrate a noninvasive local speed of sound estimator, which is based on a layered media assumption, that can accurately capture the speed of sound in the liver. We validate the estimator using an obese Zucker rat model of non-alcoholic fatty liver disease and correlate the local speed of sound with liver steatosis.Approach.We estimated the local and global average speed of sound noninvasively in 4 lean Zucker rats fed a normal diet and 16 obese Zucker rats fed a high fat diet for up to 8 weeks. The ground truth speed of sound and fat concentration were measured from the excised liver using established techniques.Main Results. The noninvasive, local speed of sound estimates of the livers were similar in value to their corresponding 'ground truth' measurements, having a slope ± standard error of the regression of 0.82 ± 0.15 (R²= 0.74 andp< 0.001). Measurement of the noninvasive global average speed of sound did not reliably capture the 'ground truth' speed of sound in the liver, having a slope of 0.35 ± 0.07 (R²= 0.74 andp< 0.001). Decreasing local speed of sound was observed with increasing hepatic fat accumulation (approximately -1.7 m s^-1per 1% increase in hepatic fat) and histopathology steatosis grading (approximately -10 to -13 m s^-1per unit increase in steatosis grade). Local speed of sound estimates were highly correlated with steatosis grade, having Pearson and Spearman correlation coefficients both ranging from -0.87 to -0.78. In addition, a lobe-dependent speed of sound in the liver was observed by theex vivomeasurements, with speed of sound differences of up to 25 m s^-1(p< 0.003) observed between lobes in the liver of the same animal.Significance.The findings of this study suggest that local speed of sound estimation has the potential to be used to predict or assist in the measurement of hepatic fat concentration and that the global average speed of sound should be avoided in hepatic fat estimation due to significant bias in the speed of sound estimate.