Synthetic Aperture Real-time Ultrasound System Research Articles

Development of Super-Resolution Sharpness-Based Axial Localization for Ultrasound Imaging

Super-resolution ultrasound mostly uses image-based methods for the localization of single scatterers. These methods are largely based on the center of mass (COM) calculation. Sharpness-based localization is an alternative to COM for scatterer localization in the axial direction. Simulated ultrasound point scatterer data (center frequency $f_{0} = 7$ MHz and wavelength $\lambda = 220\,\,\mu \text{m}$ ) showed that the normalized sharpness method can provide scatterer axial localization with an accuracy down to $2~\mu \text{m}$ ( $ ), which is a two-order of magnitude improvement compared to that achievable by the conventional imaging ( $\approx \lambda $ ), and a five-fold improvement compared to the COM estimate ( $\approx 10~\mu \text{m}$ or $0.05\lambda $ ). Similar results were obtained experimentally using wire-target data acquired by the Synthetic Aperture Real-time Ultrasound System. The performance of the proposed method was also found to be consistent across different types of ultrasound transmission. The localization precision deteriorates in the presence of noise, but even in very low signal-to-noise ratio (SNR = 0 dB), the uncertainty was not higher than $6~\mu \text{m}$ , which outperforms the COM estimate. The method can be implemented in image data as well as by using the raw signals. It is proposed that the signal-derived localization should replace the image-based equivalent, as it provides at least 10 times improved accuracy.

Resolving Ultrasound Contrast Microbubbles Using Minimum Variance Beamforming.

Minimum Variance (MV) beamforming is known to improve the lateral resolution of ultrasound images and enhance the separation of isolated point scatterers. This paper aims to evaluate the adaptive beamformer's performance with flowing microbubbles (MBs) which are relevant to super-resolution ultrasound imaging. Simulations using point scatterer data from single emissions were complemented by an experimental investigation performed using a capillary tube phantom and the Synthetic Aperture Real-time Ultrasound System (SARUS). The MV performance was assessed by the minimum distance that allows the display of two scatterers positioned side-by-side, the lateral Full-Width-at-Half-Maximum (FWHM), and the Peak-Sidelobe-Level (PSL). In the tube, scatterer responses separated by down to [Formula: see text] (or 1.05λ ) were distinguished by the MV method, while the standard Delay-And-Sum (DAS) beamformers were unable to achieve such separation. Up to ninefold FWHM decrease was also measured in favor of the MV beamformer for individual echoes from MBs. The lateral distance between two scatterers impacted on their FWHM value, and additional differences in the scatterers' axial or out-of-plane position also impacted on their size and appearance. The simulation and experimental results were in agreement in terms of lateral resolution. The point scatterer study showed that the proposed MV imaging scheme provided clear resolution benefits compared to DAS. Current super-resolution methods mainly depend on DAS beamformers. Instead, the use of the MV method may provide a larger number of detected, and potentially better localized, MB scatterers.

Directional Transverse Oscillation Vector Flow Estimation.

A method for estimating vector velocities using transverse oscillation (TO) combined with directional beamforming is presented. In directional TO (DTO), a normal focused field is emitted and the received signals are beamformed in the lateral direction transverse to the ultrasound beam to increase the amount of data for vector velocity estimation. The approach is self-calibrating as the lateral oscillation period is estimated from the directional signal through a Fourier transform to yield quantitative velocity results over a large range of depths. The approach was extensively simulated using Field IIpro and implemented on the experimental Synthetic Aperture Real-time Ultrasound System (SARUS) scanner in connection with a BK Medical 8820e convex array transducer. Velocity estimates for DTO are found for beam-to-flow angles of 60°, 75°, and 90°, and vessel depths from 24 to 156 mm. Using 16 emissions, the standard deviation (SD) for angle estimation at depths ranging from 24 to 104 mm is between 6.01° and 0.93° with a mean SD of 2.8°. The mean relative SD for the lateral velocity component is 9.2% and the mean relative bias -3.4% or four times lower than for traditional TO. The approach also works for deeper lying vessels with a slight increase in SD to 15.7%, but a maintained bias of -3.5% from 126 to 156 mm. Data for a pulsating flow have also been acquired for 15 cardiac cycles using a CompuFlow 1000 pump. The relative SD was here 7.4% for a femoral artery waveform.

Synthetic Aperture Ultrasound Fourier Beamformation Using Virtual Sources.

An efficient Fourier beamformation algorithm is presented for multistatic synthetic aperture ultrasound imaging using virtual sources. The concept is based on the frequency domain wavenumber algorithm from radar and sonar and is extended to a multielement transmit/receive configuration using virtual sources. Window functions are used to extract the azimuth processing bandwidths and weight the data to reduce side lobes in the final image. Field II simulated data and SARUS (Synthetic Aperture Real-time Ultrasound System) measured data are used to evaluate the results in terms of point spread function, resolution, contrast, signal-to-noise ratio, and processing time. Lateral resolutions of 0.53 and 0.66 mm are obtained for Fourier Beamformation Using Virtual Sources (FBV) and delay and sum (DAS) on point target simulated data. Corresponding axial resolutions are 0.21 mm for FBV and 0.20 mm for DAS. The results are also consistent over different depths evaluated using a simulated phantom containing several point targets at different depths. FBV shows a better lateral resolution at all depths, and the axial and cystic resolutions of -6, -12, and -20 dB are almost the same for FBV and DAS. To evaluate the cyst phantom metrics, three different criteria of power ratio, contrast ratio, and contrast-to-noise ratio have been used. Results show that the algorithms have a different performance in the cyst center and near the boundary. FBV has a better performance near the boundary; however, DAS is better in the more central area of the cyst. Measured data from phantoms are also used for evaluation. The results confirm applicability of FBV in ultrasound, and 20 times less processing time is attained in comparison with DAS. Evaluating the results over a wide variety of parameters and having almost the same results for simulated and measured data demonstrates the ability of FBV in preserving the quality of image as DAS, while providing a more efficient algorithm with 20 times less computations.

Safety Assessment of Advanced Imaging Sequences I: Measurements.

A method for rapid measurement of intensities (I(spta)), mechanical index (MI), and probe surface temperature for any ultrasound scanning sequence is presented. It uses the scanner's sampling capability to give an accurate measurement of the whole imaging sequence for all emissions to yield the true distributions. The method is several orders of magnitude faster than approaches using an oscilloscope, and it also facilitates validating the emitted pressure field and the scanner's emission sequence software. It has been implemented using the experimental synthetic aperture real-time ultrasound system (SARUS) scanner and the Onda AIMS III intensity measurement system (Onda Corporation, Sunnyvale, CA, USA). Four different sequences have been measured: a fixed focus emission, a duplex sequence containing B-mode and flow emissions, a vector flow sequence with B-mode and flow emissions in 17 directions, and finally a SA duplex flow sequence. A BK8820e (BK Medical, Herlev, Denmark) convex array probe is used for the first three sequences and a BK8670 linear array probe for the SA sequence. The method is shown to give the same intensity values within 0.24% of the AIMS III Soniq 5.0 (Onda Corporation, Sunnyvale, CA, USA) commercial intensity measurement program. The approach can measure and store data for a full imaging sequence in 3.8-8.2 s per spatial position. Based on I(spta), MI, and probe surface temperature, the method gives the ability to determine whether a sequence is within U.S. FDA limits, or alternatively indicate how to scale it to be within limits.

SARUS: A synthetic aperture real-time ultrasound system

The Synthetic Aperture Real-time Ultrasound System (SARUS) for acquiring and processing synthetic aperture (SA) data for research purposes is described. The specifications and design of the system are detailed, along with its performance for SA, nonlinear, and 3-D flow estimation imaging. SARUS acquires individual channel data simultaneously for up to 1024 transducer elements for a couple of heart beats, and is capable of transmitting any kind of excitation. The 64 boards in the system house 16 transmit and 16 receive channels each, where sampled channel data can be stored in 2 GB of RAM and processed using five field-programmable gate arrays (FPGAs). The fully parametric focusing unit calculates delays and apodization values in real time in 3-D space and can produce 350 million complex samples per channel per second for full non-recursive synthetic aperture B-mode imaging at roughly 30 high-resolution images/s. Both RF element data and beamformed data can be stored in the system for later storage and processing. The stored data can be transferred in parallel using the system's sixty-four 1-Gbit Ethernet interfaces at a theoretical rate of 3.2 GB/s to a 144-core Linux cluster.