Software beamforming: comparison between a phased array and synthetic transmit aperture

Elevational Synthetic Aperture Focusing for Three-Dimensional Photoacoustic Imaging Using a Clinical One-Dimensional Array Transducer

OBJECTIVE

Two-dimensional (2D) photoacoustic (PA) imaging based on array transducers provide high spatial resolution in the lateral direction by adopting receive dynamic focusing. However, the quality of PA image is often deteriorated by poor elevational resolution which is achieved by an acoustic lens. To overcome this limitation, we present a three-dimensional (3D) image reconstruction method using a commercial one-dimensional (1D) array transducer.

METHODS

In the method, the elevational resolution is improved by applying synthetic aperture focusing (SAF) technique along the elevational direction. For this, a commercially available 1D array transducer with an acoustic lens is modeled and appropriate synthetic focusing delay that can minimize the effect of the acoustic lens is derived by mathematical analysis.

RESULTS

From the simulation and experiment results, it was demonstrated that the proposed method can enhance the image quality of PA imaging, i.e., elevational resolution and signal-to-noise ratio (SNR).

CONCLUSION

3D PA images with improved elevational resolution were achieved using a clinical 1D array transducer.

SIGNIFICANCE

The presented method may be useful for clinical application such as detecting microcalcification, imaging of tumor vasculature and guidance of biopsy in real time.

Real-Time Volumetric Synthetic Aperture Software Beamforming of Row-Column Probe Data

Two delay-and-sum beamformers for 3-D synthetic aperture imaging with row-column addressed arrays are presented. Both beamformers are software implementations for graphics processing unit (GPU) execution with dynamic apodizations and third-order polynomial subsample interpolation. The first beamformer was written in the MATLAB programming language and the second was written in C/C++ with the compute unified device architecture (CUDA) extensions by NVIDIA. Performance was measured as volume rate and sample throughput on three different GPUs: a 1050 Ti, a 1080 Ti, and a TITAN V. The beamformers were evaluated across 112 combinations of output geometry, depth range, transducer array size, number of virtual sources, floating-point precision, and Nyquist rate or in-phase/quadrature beamforming using analytic signals. Real-time imaging defined as more than 30 volumes per second was attained by the CUDA beamformer on the three GPUs for 13, 27, and 43 setups, respectively. The MATLAB beamformer did not attain real-time imaging for any setup. The median, single-precision sample throughput of the CUDA beamformer was 4.9, 20.8, and 33.5 Gsamples/s on the three GPUs, respectively. The throughput of CUDA beamformer was an order of magnitude higher than that of the MATLAB beamformer.

Synthetic Aperture Imaging Using High-Frequency Convex Array for Ophthalmic Ultrasound Applications

High-frequency ultrasound (HFUS) imaging has emerged as an essential tool for pre-clinical studies and clinical applications such as ophthalmic and dermatologic imaging. HFUS imaging systems based on array transducers capable of dynamic receive focusing have considerably improved the image quality in terms of spatial resolution and signal-to-noise ratio (SNR) compared to those by the single-element transducer-based one. However, the array system still suffers from low spatial resolution and SNR in out-of-focus regions, resulting in a blurred image and a limited penetration depth. In this paper, we present synthetic aperture imaging with a virtual source (SA-VS) for an ophthalmic application using a high-frequency convex array transducer. The performances of the SA-VS were evaluated with phantom and ex vivo experiments in comparison with the conventional dynamic receive focusing method. Pre-beamformed radio-frequency (RF) data from phantoms and excised bovine eye were acquired using a custom-built 64-channel imaging system. In the phantom experiments, the SA-VS method showed improved lateral resolution (>10%) and sidelobe level (>4.4 dB) compared to those by the conventional method. The SNR was also improved, resulting in an increased penetration depth: 16 mm and 23 mm for the conventional and SA-VS methods, respectively. Ex vivo images with the SA-VS showed improved image quality at the entire depth and visualized structures that were obscured by noise in conventional imaging.

In Vivo Evaluation of Plane Wave Imaging for Abdominal Ultrasonography

Although plane wave imaging (PWI) has been extensively employed for ultrafast ultrasound imaging, its potential for sectorial B-mode imaging with a convex array transducer has not yet been widely recognized. Recently, we reported an optimized PWI approach for sector scanning that exploits the dynamic transmit focusing capability. In this paper, we first report the clinical applicability of the optimized PWI for abdominal ultrasonography by in vivo image and video evaluations and compare it with conventional focusing (CF) and diverging wave imaging (DWI), which is another dynamic transmit focusing technique generally used for sectorial imaging. In vivo images and videos of the liver, kidney, and gallbladder were obtained from 30 healthy volunteers using PWI, DWI, and CF. Three radiologists assessed the phantom images, 156 in vivo images, and 66 in vivo videos. PWI showed significantly enhanced (p < 0.05) spatial resolution, contrast, and noise and artifact reduction, and a 4-fold higher acquisition rate compared to CF and provided similar performances compared to DWI. Because the computations required for PWI are considerably lower than that for DWI, PWI may represent a promising technique for sectorial imaging in abdominal ultrasonography that provides better image quality and eliminates the need for focal depth adjustment.

Efficient Sample Delay Calculation for 2-D and 3-D Ultrasound Imaging

Ultrasound imaging is a reference medical diagnostic technique, thanks to its blend of versatility, effectiveness, and moderate cost. The core computation of all ultrasound imaging methods is based on simple formulae, except for those required to calculate acoustic propagation delays with high precision and throughput. Unfortunately, advanced three-dimensional (3-D) systems require the calculation or storage of billions of such delay values per frame, which is a challenge. In 2-D systems, this requirement can be four orders of magnitude lower, but efficient computation is still crucial in view of low-power implementations that can be battery-operated, enabling usage in numerous additional scenarios. In this paper, we explore two smart designs of the delay generation function. To quantify their hardware cost, we implement them on FPGA and study their footprint and performance. We evaluate how these architectures scale to different ultrasound applications, from a low-power 2-D system to a next-generation 3-D machine. When using numerical approximations, we demonstrate the ability to generate delay values with sufficient throughput to support 10 000-channel 3-D imaging at up to 30 fps while using 63% of a Virtex 7 FPGA, requiring 24 MB of external memory accessed at about 32 GB/s bandwidth. Alternatively, with similar FPGA occupation, we show an exact calculation method that reaches 24 fps on 1225-channel 3-D imaging and does not require external memory at all. Both designs can be scaled to use a negligible amount of resources for 2-D imaging in low-power applications and for ultrafast 2-D imaging at hundreds of frames per second.

Development of a Low-Complexity, Cost-Effective Digital Beamformer Architecture for High-Frequency Ultrasound Imaging

This paper presents a low-complexity, cost-effective digital beamformer architecture for a high-frequency ultrasound imaging system. The proposed beamformer uses a lookup table and linear interpolation methods for computing the dynamic receive focusing delays and a postfiltering technique to minimize hardware complexity. In the postfiltering technique, channel radio-frequency data having the same fractional delay (i.e., 16f₀ resolution) are aggregated prior to interpolation. Thus, only four polyphase structure filters are required in the developed beamformer. In addition, a quadrature bandpass filter that generates an analytic signal is utilized as an interpolation filter; this allows decimation during beam formation and a reduction in computational complexity. The proposed method was evaluated through a 20- [Formula: see text] wire phantom experiment, and the -6-dB lateral and axial resolutions obtained therein were measured and compared with those obtained using a conventional method. The same lateral (165 [Formula: see text]) and axial (80 [Formula: see text]) resolutions at a depth of 5.6 mm were obtained using both the methods, and the proposed method could reduce the beamforming points (i.e., computational complexity) by a factor of the decimation factor (≥4). Images from an excised bovine eye were captured; they showed that the proposed beamformer identified fine anatomical structures such as cornea or iris without compromising the spatial resolution and reduced the computational complexity.

High intensity focused ultrasound (HIFU) focal spot localization using harmonic motion imaging (HMI)

Several ultrasound-based imaging modalities have been proposed for image guidance and monitoring of high-intensity focused ultrasound (HIFU) treatment. However, accurate localization and characterization of the effective region of treatment (focal spot) remain important obstacles in the clinical implementation of HIFU ablation. Harmonic motion imaging for focused ultrasound (HMIFU) is a HIFU monitoring technique that utilizes radiation-force-induced localized oscillatory displacement. HMIFU has been shown to correctly identify the formation and extent of HIFU thermal ablation lesions. However a significant problem remains in identifying the location of the HIFU focus, which is necessary for treatment planning. In this study, the induced displacement was employed to localize the HIFU focal spot inside the tissue prior to treatment. Feasibility was shown with two separate systems. The 1D HMIFU system consisted of a HIFU transducer emitting an amplitude-modulated HIFU beam for mechanical excitation and a confocal single-element, pulse-echo transducer for simultaneous RF acquisition. The 2D HIFU system consists of a HIFU phased array, and a co-axial imaging phased array for simultaneous imaging. Initial feasibility was first performed on tissue-mimicking gelatin phantoms and the focal zone was defined as the region corresponding to the -3dB full width at half maximum of the HMI displacement. Using the same parameters, in vitro experiments were performed in canine liver specimens to compare the defined focal zone with the lesion. In vitro measurements showed good agreement between the HMI predicted focal zone and the induced HIFU lesion location. HMIFU was experimentally shown to be capable of predicting and tracking the focal region in both phantoms and in vitro tissues. The accuracy of focal spot localization was evaluated by comparing with the lesion location in post-ablative tissues, with a R(2) = 0.821 at p < 0.002 in the 2D HMI system. We demonstrated the feasibility of using this HMI-based technique to localize the HIFU focal spot without inducing thermal changes during the planning phase. The focal spot localization method has also been applied on ex vivo human breast tissue ablation and can be fully integrated into any HMI system for planning purposes.

Two-dimensional shear-wave elastography on conventional ultrasound scanners with time-aligned sequential tracking (TAST) and comb-push ultrasound shear elastography (CUSE)

Two-dimensional shear-wave elastography presents 2-D quantitative shear elasticity maps of tissue, which are clinically useful for both focal lesion detection and diffuse disease diagnosis. Realization of 2-D shear-wave elastography on conventional ultrasound scanners, however, is challenging because of the low tracking pulse-repetition-frequency (PRF) of these systems. Although some clinical and research platforms support software beamforming and plane-wave imaging with high PRF, the majority of current clinical ultrasound systems do not have the software beamforming capability, which presents a critical challenge for translating the 2-D shear-wave elastography technique from laboratory to clinical scanners. To address this challenge, this paper presents a time-aligned sequential tracking (TAST) method for shear-wave tracking on conventional ultrasound scanners. TAST takes advantage of the parallel beamforming capability of conventional systems and realizes high-PRF shear-wave tracking by sequentially firing tracking vectors and aligning shear wave data in the temporal direction. The comb-push ultrasound shear elastography (CUSE) technique was used to simultaneously produce multiple shear wave sources within the field-of-view (FOV) to enhance shear wave SNR and facilitate robust reconstructions of 2-D elasticity maps. TAST and CUSE were realized on a conventional ultrasound scanner. A phantom study showed that the shear-wave speed measurements from the conventional ultrasound scanner were in good agreement with the values measured from other 2-D shear wave imaging technologies. An inclusion phantom study showed that the conventional ultrasound scanner had comparable performance to a state-of-the-art shear-wave imaging system in terms of bias and precision in measuring different sized inclusions. Finally, in vivo case analysis of a breast with a malignant mass, and a liver from a healthy subject demonstrated the feasibility of using the conventional ultrasound scanner for in vivo 2-D shear-wave elastography. These promising results indicate that the proposed technique can enable the implementation of 2-D shear-wave elastography on conventional ultrasound scanners and potentially facilitate wider clinical applications with shear-wave elastography.

Sparse matrix beamforming and image reconstruction for 2-D HIFU monitoring using harmonic motion imaging for focused ultrasound (HMIFU) with in vitro validation

Harmonic motion imaging for focused ultrasound (HMIFU) utilizes an amplitude-modulated HIFU beam to induce a localized focal oscillatory motion simultaneously estimated. The objective of this study is to develop and show the feasibility of a novel fast beamforming algorithm for image reconstruction using GPU-based sparse-matrix operation with real-time feedback. In this study, the algorithm was implemented onto a fully integrated, clinically relevant HMIFU system. A single divergent transmit beam was used while fast beamforming was implemented using a GPU-based delay-and-sum method and a sparse-matrix operation. Axial HMI displacements were then estimated from the RF signals using a 1-D normalized cross-correlation method and streamed to a graphic user interface with frame rates up to 15 Hz, a 100-fold increase compared to conventional CPU-based processing. The real-time feedback rate does not require interrupting the HIFU treatment. Results in phantom experiments showed reproducible HMI images and monitoring of 22 in vitro HIFU treatments using the new 2-D system demonstrated reproducible displacement imaging, and monitoring of 22 in vitro HIFU treatments using the new 2-D system showed a consistent average focal displacement decrease of 46.7 ±14.6% during lesion formation. Complementary focal temperature monitoring also indicated an average rate of displacement increase and decrease with focal temperature at 0.84±1.15%/(°)C, and 2.03±0.93%/(°)C , respectively. These results reinforce the HMIFU capability of estimating and monitoring stiffness related changes in real time. Current ongoing studies include clinical translation of the presented system for monitoring of HIFU treatment for breast and pancreatic tumor applications.