Fast parametric beamformer for synthetic aperture imaging

A Noniterative Radix-8 CORDIC Algorithm with Low Latency and High Efficiency

An efficient, noniterative Radix-8 (NR-8) coordinate rotation digital computer (CORDIC) algorithm is proposed for low-latency and high-efficiency computation of the functions of sine, cosine, or the phase shift, with which the values of the functions are precisely computed by only using the angle in a narrow range of [0, π/12] rather than in a wide angle range of [0, π/2]. This algorithm is expressed by a formula that simplifies the traditional iterative processes by using a complex multiplier. The results obtained from the simulation and the experiment on an FPGA show that the NR-8 CORDIC algorithm operates well, with which the 16-bit precision output is extremely precise, with only 0.012% of the absolute error for computing the sine or cosine function with a step of 0.001°. Compared with the best conventional CORDIC algorithm, the clock latency of this algorithm significantly decreases down to less than 50%, only needs half of the logic resources and consumes half of the power. This algorithm also takes advantages over other newly improved CORDIC algorithms and requires less than half of the clock latency, even for a 23-bit precision output. Therefore, this algorithm could provide a potential application in real-time systems such as radar digital beamforming.

Acoustic Field Characterization of Medical Array Transducers Based on Unfocused Transmits and Single-Plane Hydrophone Measurements

Medical ultrasonic arrays are typically characterized in controlled water baths using measurements by a hydrophone, which can be translated with a positioning stage. Characterization of 3D acoustic fields conventionally requires measurements at each spatial location, which is tedious and time-consuming, and may be prohibitive given limitations of experimental setup (e.g., the bath and stage) and measurement equipment (i.e., the hydrophone). Moreover, with the development of new ultrasound sequences and modalities, multiple measurements are often required to characterize each imaging mode to ensure performance and clinical safety. Acoustic holography allows efficient characterization of source transducer fields based on single plane measurements. In this work, we explore the applicability of a re-radiation method based on the Rayleigh–Sommerfeld integral to medical imaging array characterization. We show that source fields can be reconstructed at single crystal level at wavelength resolution, based on far-field measurements. This is herein presented for three practical application scenarios: for identifying faulty transducer elements; for characterizing acoustic safety parameters in focused ultrasound sequences from 2D planar measurements; and for estimating arbitrary focused fields based on calibration from an unfocused sound field and software beamforming. The results experimentally show that the acquired pressure fields closely match those estimated using our technique.

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.

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.

Experimental performance assessment of the sub-band minimum variance beamformer for ultrasound imaging

Recent progress in adaptive beamforming techniques for medical ultrasound has shown that current resolution limits can be surpassed. One method of obtaining improved lateral resolution is the Minimum Variance (MV) beamformer. The frequency domain implementation of this method effectively divides the broadband ultrasound signals into sub-bands (MVS) to conform with the narrow-band assumption of the original MV theory. This approach is investigated here using experimental Synthetic Aperture (SA) data from wire and cyst phantoms. A 7MHz linear array transducer is used with the SARUS experimental ultrasound scanner for the data acquisition. The lateral resolution and the contrast obtained, are evaluated and compared with those from the conventional Delay-and-Sum (DAS) beamformer and the MV temporal implementation (MVT). From the wire phantom the Full-Width-at-Half-Maximum (FWHM) measured at a depth of 52mm, is 16.7μm (0.08λ) for both MV methods, while the corresponding values for the DAS case are at least 24 times higher. The measured Peak-Side-lobe-Level (PSL) may reach -41dB using the MVS approach, while the values from the DAS and MVT beamforming are above -24dB and -33dB, respectively. From the cyst phantom, the power ratio (PR), the contrast-to-noise ratio (CNR), and the speckle signal-to-noise ratio (sSNR) measured at a depth of 30mm are at best similar for MVS and DAS, with values ranging between -29dB and -30dB, 1.94 and 2.05, and 2.16 and 2.27 respectively. In conclusion the MVS beamformer is not suitable for imaging continuous targets, and significant resolution gains were obtained only for isolated targets.

Software-based high-level synthesis design of FPGA beamformers for synthetic aperture imaging

Field-programmable gate arrays (FPGAs) can potentially be configured as beamforming platforms for ultrasound imaging, but a long design time and skilled expertise in hardware programming are typically required. In this article, we present a novel approach to the efficient design of FPGA beamformers for synthetic aperture (SA) imaging via the use of software-based high-level synthesis techniques. Software kernels (coded in OpenCL) were first developed to stage-wise handle SA beamforming operations, and their corresponding FPGA logic circuitry was emulated through a high-level synthesis framework. After design space analysis, the fine-tuned OpenCL kernels were compiled into register transfer level descriptions to configure an FPGA as a beamformer module. The processing performance of this beamformer was assessed through a series of offline emulation experiments that sought to derive beamformed images from SA channel-domain raw data (40-MHz sampling rate, 12 bit resolution). With 128 channels, our FPGA-based SA beamformer can achieve 41 frames per second (fps) processing throughput (3.44 × 10(8) pixels per second for frame size of 256 × 256 pixels) at 31.5 W power consumption (1.30 fps/W power efficiency). It utilized 86.9% of the FPGA fabric and operated at a 196.5 MHz clock frequency (after optimization). Based on these findings, we anticipate that FPGA and high-level synthesis can together foster rapid prototyping of real-time ultrasound processor modules at low power consumption budgets.

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.

Modified multi-element synthetic transmit aperture method for ultrasound imaging: A tissue phantom study

The paper presents the modified multi-element synthetic transmit aperture (MSTA) method for ultrasound imaging. It is based on coherent summation of RF echo signals with apodization weights taking into account the finite size of the transmit subaperture and of the receive element. The work presents extension of the previous study where the modified synthetic transmit aperture (STA) method was considered and verified [1]. In the case of MSTA algorithm the apodization weights were calculated for each imaging point and all combinations of the transmit subaperture and receive element using their angular directivity functions (ADFs). The ADFs were obtained from the exact solution of the corresponding mixed boundary-value problem for periodic baffle system modeling the transducer array. Performance of the developed method was tested using Field II simulated synthetic aperture data of point reflectors for 4MHz 128-element transducer array with 0.3mm pitch and 0.02mm kerf to estimate the visualization depth and lateral resolution. Also experimentally determined data of the tissue-mimicking phantom (Dansk Fantom Service, model 571) obtained using 128 elements, 4MHz, linear transducer array (model L14-5/38) and Ultrasonix SonixTOUCH Research platform were used for qualitative assessment of imaging contrast improvement. Comparison of the results obtained by the modified and conventional MSTA algorithms indicated 15dB improvement of the noise reduction in the vicinity of transducer's surface (1mm depth), and concurrent increase in the visualization depth (86% augment of the scattered amplitude at the depth of 90mm). However, this increase was achieved at the expense of minor degradation of the lateral resolution of approximately 8% at the depth of 50mm and 5% at the depth of 90mm.

Synthetic aperture focusing for a single-element transducer undergoing helical motion

This paper describes the application of 3-D synthetic aperture focusing (SAF) to a single-element trans-rectal ultrasound transducer. The transducer samples a 3-D volume by simultaneous rotation and translation, giving a helical motion. Two different 3-D SAF methods are investigated, a direct and a two-step approach. Both methods perform almost identically for simulated scatterers and give a significant improvement in azimuth resolution and a constant resolution in elevation. Side lobes below -60 dB are achievable for both methods. Validation of the method is achieved by scanning a simple wire phantom and a complex phantom containing wires in azimuth and elevation. The simple wire phantom shows the same results as that found through simulation. The complex phantom shows simultaneous focusing in azimuth and elevation for the wire scatterers. Consideration of the processing requirements for both 3-D SAF methods shows that the two-step approach can give equivalent performance using an order of magnitude fewer calculations. This reduction requires a temporary storage of 9.1 GB of data for the investigated setup.

A new synthetic aperture focusing method to suppress the diffraction of ultrasound

Spatial resolution of an ultrasound image is limited by diffraction of ultrasound as it propagates along the axial direction. This paper proposes a method for reducing the diffraction spreading effect of ultrasound by using a synthetic aperture focusing (SAF) method that uses plane waves instead of spherical waves. The new method performs data acquisition and beamforming in the same manner as conventional SAF methods. The main difference is that all array elements are used on each firing to generate a plane wave, the traveling angle of which varies with the position of a receive subaperture. On reception, each scan line is formed by synthesizing RF samples acquired by relevant receive subapertures with delays to force the plane waves to meet at each imaging point. Theoretical analysis and computer simulation with infinite transmit aperture show that the proposed method is capable of suppressing the diffraction of ultrasound and especially causing the lateral beam width to remain unchanged beyond a certain depth determined by the size of a receive subaperture and the maximum traveling angle of plane waves. It is demonstrated that the proposed method is realizable using a linear array transducer. It is also shown that the lateral radiation pattern produced by the proposed method has smaller beam width than that using conventional SAF methods in the region of interest because it suppresses the diffraction of ultrasound.

Novel low-power ultrasound digital preprocessing architecture for wireless display

A complete hardware-based ultrasound preprocessing unit (PPU) is presented as an alternative to available power-hungry devices. Intended to expand the ultrasonic applications, the proposed unit allows replacement of the cable of the ultrasonic probe by a wireless link to transfer data from the probe to a remote monitor. The digital back-end architecture of this PPU is fully pipelined, which permits sampling of ultrasonic signals at a frequency equal to the field-programmable gate array-based system clock, up to 100 MHz. Experimental results show that the proposed processing unit has an excellent performance, an equivalent 53.15 Dhrystone 2.1 MIPS/ MHz (DMIPS/MHz), compared with other software-based architectures that allow a maximum of 1.6 DMIPS/MHz. In addition, an adaptive subsampling method is proposed to operate the pixel compressor, which allows real-time image zooming and, by removing high-frequency noise, the lateral and axial resolutions are enhanced by 25% and 33%, respectively. Realtime images, acquired from a reference phantom, validated the feasibility of the proposed architecture. For a display rate of 15 frames per second, and a 5-MHz single-element piezoelectric transducer, the proposed digital PPU requires a dynamic power of only 242 mW, which represents around 20% of the best-available software-based system. Furthermore, composed by the ultrasound processor and the image interpolation unit, the digital processing core of the PPU presents good power-performance ratios of 26 DMIPS/mW and 43.9 DMIPS/mW at a 20-MHz and 100-MHz sample frequency, respectively.