Imaging Performance for Two Row-Column Arrays.

This study evaluates the volumetric imaging performance of two prototyped 62 + 62 row-column-addressed (RCA) 2-D array transducer probes using three synthetic aperture imaging (SAI) emission sequences and two different beamformers. The probes are fabricated using capacitive micromachined ultrasonic transducer (CMUT) and piezoelectric transducer (PZT) technology. Both have integrated apodization to reduce ghost echoes and are designed with similar acoustical features, i.e., 3-MHz center frequency, λ /2 pitch, and [Formula: see text] active footprint. Raw RF data are obtained using an experimental research ultrasound scanner, SARUS. The SAI sequences are designed for imaging down to 14 cm at a volume rate of 88 Hz. Two beamforming methods, spatial matched filtering and row-column adapted delay-and-sum, are used for beamforming the RF data. The imaging quality is investigated through simulations and phantom measurements. Both probes on average have similar lateral full-width at half-maximum (FWHM) values, but the PZT probe has 20% smaller cystic resolution values and 70% larger contrast-to-noise ratio (CNR) compared to the capacitive micromachined ultrasonic transducer (CMUT) probe. The CMUT probe can penetrate down to 15 cm, and the PZT probe down to 30 cm. The CMUT probe has 17% smaller axial FWHM values. The matched filter focusing shows an improved B-mode image for measurements on a cyst phantom with an improved speckle pattern and better visualization of deeper lying cysts. The results of this study demonstrate the potentials of RCA 2-D arrays against fully addressed 2-D arrays, which are low channel count (e.g., 124 instead of 3844), low acoustic intensity mechanical index (MI ≤ 0.88 and spatial-peak-temporal-average intensity [Formula: see text]), and high penetration depth (down to 30 cm), which makes 3-D imaging at high volume rates possible with equipment in the price range of conventional 2-D imaging.

Use of a Real-Time Stress Map for Assessment of Applied Stress for Strain Elastography: Utility in Training and Computation of Strain Ratios.

OBJECTIVES: There is a significant learning curve in strain elastography. Uniform appropriate levels of stress must be applied for accurate elastograms. If the stress is not applied appropriately, inaccurate results will be obtained, particularly when strain ratios are being estimated. This paper describes a new technique which allows the real-time visualization of the applied stress with a color-coded stress map. The potential use of this map is discussed. METHODS: Ten patients (5 breast, 5 thyroid) and phantoms were scanned using the stress map. The stress applied was varied and the resultant change in the strain image evaluated. RESULTS: The stress map was able to document if appropriate stress was applied when performing strain elastography. When inappropriate stress was applied or physiological process effected the strain image the stress map demonstrated the areas of inaccurate measurements in the stress map. CONCLUSIONS: The display of a stress map that depicts the degree and uniformity of applied stress would be helpful both for training of the appropriate technique, for confirming that the elastogram is appropriate for evaluation, and that strain ratio estimates are accurate.

Curvilinear 3-D Imaging Using Row-Column-Addressed 2-D Arrays With a Diverging Lens: Phantom Study.

A double-curved diverging lens over the flat row-column-addressed (RCA) 2-D array can extend its inherent rectilinear 3-D imaging field of view (FOV) to a curvilinear volume region, which is necessary for applications such as abdominal and cardiac imaging. Two concave lenses with radii of 12.7 and 25.4 mm were manufactured using RTV664 silicone. The diverging properties of the lenses were evaluated based on simulations and measurements on several phantoms. The measured FOV for both lenses in contact with tissue mimicking phantom was less than 15% different from the theoretical predictions, i.e., a curvilinear FOV of and for the 12.7- and 25.4-mm radii lenses. A synthetic aperture imaging sequence with single-element transmissions was designed for imaging down to 140 mm at a volume rate of 88 Hz. The performance was evaluated in terms of signal-to-noise ratio, FOV, and full-width at half-maximum (FWHM) of a focused beam. The penetration depths in a tissue mimicking phantom with 0.5-dB/(cm MHz) attenuation were 100 and 125 mm for the lenses with radii of 12.7 and 25.4 mm. The azimuth, elevation, and radial FWHM at 43-mm depth were (5.8, 5.8, 1) and (6, 6, 1) . The results of this study confirm that the proposed lens approach is an effective method for increasing the FOV, when imaging with RCA 2-D arrays.

Curvilinear 3-D Imaging Using Row-Column-Addressed 2-D Arrays With a Diverging Lens: Feasibility Study.

Constructing a double-curved row-column-addressed (RCA) 2-D array or applying a diverging lens over the flat RCA 2-D array can extend the imaging field-of-view (FOV) to a curvilinear volume without increasing the aperture size, which is necessary for applications, such as abdominal and cardiac imaging. Extended FOV and low channel count of double-curved RCA 2-D arrays make 3-D imaging possible with equipment in the price range of conventional 2-D imaging. This paper proposes a delay-and-sum beamformation scheme specific to double-curved RCA 2-D arrays and validates its focusing ability based on simulations. A synthetic aperture imaging sequence with single element transmissions is designed for imaging down to 14 cm at a volume rate of 88 Hz. Using a diverging lens with an f-number of -1 circumscribing the underlying RCA array, the imaging quality of a double-curved λ/2 -pitch 3-MHz 62 + 62 RCA 2-D array is investigated as a function of depth within a curvilinear FOV of 60 °×60° . The simulated double-curved 2-D array exhibits the same full-width-at-half-maximum values for a point scatterer within its curvilinear FOV at a fixed radial distance compared with a flat 2-D array within its rectilinear FOV. The results of this paper demonstrate that the proposed beamforming approach is accurate for achieving correct time-of-flight calculations, and hence avoids geometrical distortions.

Ultrasound Vector Flow Imaging-Part I: Sequential Systems.

This paper gives a review of the most important methods for blood velocity vector flow imaging (VFI) for conventional sequential data acquisition. This includes multibeam methods, speckle tracking, transverse oscillation, color flow mapping derived VFI, directional beamforming, and variants of these. The review covers both 2-D and 3-D velocity estimation and gives a historical perspective on the development along with a summary of various vector flow visualization algorithms. The current state of the art is explained along with an overview of clinical studies conducted and methods for presenting and using VFI. A number of examples of VFI images are presented, and the current limitations and potential solutions are discussed.

Ultrasound Vector Flow Imaging-Part II: Parallel Systems.

This paper gives a review of the current state-of-the-art in ultrasound parallel acquisition systems for flow imaging using spherical and plane waves emissions. The imaging methods are explained along with the advantages of using these very fast and sensitive velocity estimators. These experimental systems are capable of acquiring thousands of images per second for fast moving flow as well as yielding the estimates of low velocity flow. These emerging techniques allow the vector flow systems to assess highly complex flow with transitory vortices and moving tissue, and they can also be used in functional ultrasound imaging for studying brain function in animals. This paper explains the underlying acquisition and estimation methods for fast 2-D and 3-D velocity imaging and gives a number of examples. Future challenges and the potentials of parallel acquisition systems for flow imaging are also discussed.

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.

Implementation of a versatile research data acquisition system using a commercially available medical ultrasound scanner.

This paper describes the design and implementation of a versatile, open-architecture research data acquisition system using a commercially available medical ultrasound scanner. The open architecture will allow researchers and clinicians to rapidly develop applications and move them relatively easy to the clinic. The system consists of a standard PC equipped with a camera link and an ultrasound scanner equipped with a research interface. The ultrasound scanner is an easy-to-use imaging device that is capable of generating high-quality images. In addition to supporting the acquisition of multiple data types, such as B-mode, M-mode, pulsed Doppler, and color flow imaging, the machine provides users with full control over imaging parameters such as transmit level, excitation waveform, beam angle, and focal depth. Beamformed RF data can be acquired from regions of interest throughout the image plane and stored to a file with a simple button press. For clinical trials and investigational purposes, when an identical image plane is desired for both an experimental and a reference data set, interleaved data can be captured. This form of data acquisition allows switching between multiple setups while maintaining identical transducer, scanner, region of interest, and recording time. Data acquisition is controlled through a graphical user interface running on the PC. This program implements an interface for third-party software to interact with the application. A software development toolkit is developed to give researchers and clinicians the ability to utilize third-party software for data analysis and flexible manipulation of control parameters. Because of the advantages of speed of acquisition and clinical benefit, research projects have successfully used the system to test and implement their customized solutions for different applications. Three examples of system use are presented in this paper: evaluation of synthetic aperture sequential beamformation, transverse oscillation for blood velocity estimation, and acquisition of spectral velocity data for evaluating aortic aneurysms.

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.

Three-dimensional synthetic aperture focusing using a rocking convex array transducer.

Volumetric imaging can be performed using 1-D arrays in combination with mechanical motion. Outside the elevation focus of the array, the resolution and contrast quickly degrade compared with the lateral plane, because of the fixed transducer focus. This paper shows the feasibility of using synthetic aperture focusing for enhancing the elevation focus for a convex rocking array. The method uses a virtual source (VS) for defocused multi-element transmit, and another VS in the elevation focus point. This allows a direct time-of-flight to be calculated for a given 3-D point. To avoid artifacts and increase SNR at the elevation VS, a plane-wave VS approach has been implemented. Simulations and measurements using an experimental scanner with a convex rocking array show an average improvement in resolution of 26% and 33%, respectively. This improvement is also seen in in vivo measurements. An evaluation of how a change in transducer design will affect the resolution improvement shows a potential for using a modified transducer for 3-D imaging with improved elevation focusing and contrast.

Precise time-of-flight calculation for 3-D synthetic aperture focusing.

Conventional linear arrays can be used for 3-D ultrasound imaging by moving the array in the elevation direction and stacking the planes in a volume. The point-spread function is larger in the elevation plane, because the aperture is smaller and has a fixed elevation focus. Resolution improvements in elevation can be achieved by applying synthetic aperture focusing to the beamformed-in-plane RF data. The proposed method uses a virtual source placed at the elevation focus for postbeamforming. This has previously been done in 2 steps, in-plane focusing followed by synthetic aperture postfocusing in elevation, due to lack of a simple expression for the exact time of flight. This paper presents a new single step method for calculating the time of flight for a 3-D case using a linear array. The new method is more flexible and is able to beamform a fewer number of points much more efficiently. The method is evaluated using both simulated data and phantom measurements using the RASMUS experimental scanner. Computational cost of the method is higher than the 2-step method for a full volume beamforming, but it allows for a reduction of an order-of-magnitude if 3 planes are used for real-time visualization. In addition, the need for a temporary storage of beamformed data is removed.

Fast parametric beamformer for synthetic aperture imaging.

This paper describes the design and implementation of a real-time delay-and-sum synthetic aperture beamformer. The beamforming delays and apodization coefficients are described parametrically. The image is viewed as a set of independent lines that are defined in 3D by their origin, direction, and inter-sample distance. The delay calculation is recursive and inspired by the coordinate rotation digital computer (CORDIC) algorithm. Only 3 parameters per channel and line are needed for their generation. The calculation of apodization coefficients is based on a piece- wise linear approximation. The implementation of the beamformer is optimized with respect to the architecture of a novel synthetic aperture real-time ultrasound scanner (SARUS), in which 4 channels are processed by the same set of field-programmable gate arrays (FPGA). In synthetic transmit aperture imaging, low-resolution images are formed after every emission. Summing all low-resolution images produces a perfectly focused high-resolution image. The design of the beamformer is modular, and a single beamformation unit can produce 4600 low-resolution images per second, each consisting of 32 lines and 1024 complex samples per line. In its present incarnation, 3 such modules fit in a single device. The summation of low-resolution images is performed internally in the FPGA to reduce the required bandwidth. The delays are calculated with a precision of 1/16th of a sample, and the apodization coefficients with 7-bit precision. The accumulation of low-resolution images is performed with 24-bit precision. The level of the side- and grating lobes, introduced by the use of integer numbers in the calculations and truncation of intermediate results, is below -86 dB from the peak.

Synthetic aperture ultrasound imaging.

The paper describes the use of synthetic aperture (SA) imaging in medical ultrasound. SA imaging is a radical break with today's commercial systems, where the image is acquired sequentially one image line at a time. This puts a strict limit on the frame rate and the possibility of acquiring a sufficient amount of data for high precision flow estimation. These constrictions can be lifted by employing SA imaging. Here data is acquired simultaneously from all directions over a number of emissions, and the full image can be reconstructed from this data. The paper demonstrates the many benefits of SA imaging. Due to the complete data set, it is possible to have both dynamic transmit and receive focusing to improve contrast and resolution. It is also possible to improve penetration depth by employing codes during ultrasound transmission. Data sets for vector flow imaging can be acquired using short imaging sequences, whereby both the correct velocity magnitude and angle can be estimated. A number of examples of both phantom and in vivo SA images will be presented measured by the experimental ultrasound scanner RASMUS to demonstrate the many benefits of SA imaging.

Ultrasound research scanner for real-time synthetic aperture data acquisition.

Conventional ultrasound systems acquire ultrasound data sequentially one image line at a time. The architecture of these systems is therefore also sequential in nature and processes most of the data in a sequential pipeline. This often makes it difficult to implement radically different imaging strategies on the platforms and makes the scanners less accessible for research purposes. A system designed for imaging research flexibility is the prime concern. The possibility of sending out arbitrary signals and the storage of data from all transducer elements for 5 to 10 seconds allows clinical evaluation of synthetic aperture and 3D imaging. This paper describes a real-time system specifically designed for research purposes. The system can acquire multichannel data in real-time from multi-element ultrasound transducers, and can perform some real-time processing on the acquired data. The system is capable of performing real-time beamforming for conventional imaging methods using linear, phased, and convex arrays. Image acquisition modes can be intermixed, and this makes it possible to perform initial trials in a clinical environment with new imaging modalities for synthetic aperture imaging, 2D and 3D B-mode, and velocity imaging using advanced coded emissions. The system can be used with 128-element transducers and can excite 128 transducer elements and receive and sample data from 64 channels simultaneously at 40 MHz with 12-bit precision. Two-to-one multiplexing in receive can be used to cover 128 receive channels. Data can be beamformed in real time using the system's 80 signal processing units, or it can be stored directly in RAM. The system has 16 Gbytes RAM and can, thus, store more than 3.4 seconds of multichannel data. It is fully software programmable and its signal processing units can also be reconfigured under software control. The control of the system is done over a 100-Mbits/s Ethernet using C and Matlab. Programs for doing, e.g., B-mode imaging can be written directly in Matlab and executed on the system over the net from any workstation running Matlab. The overall system concept is presented along with its implementation and examples of B-mode and in vivo synthetic aperture flow imaging.

Directional synthetic aperture flow imaging.

A method for flow estimation using synthetic aperture imaging and focusing along the flow direction is presented. The method can find the correct velocity magnitude for any flow angle, and full color flow images can be measured using only 32 to 128 pulse emissions. The approach uses spherical wave emissions with a number of defocused elements and a linear frequency-modulated pulse (chirp) to improve the signal-to-noise ratio. The received signals are dynamically focused along the flow direction and these signals are used in a cross-correlation estimator for finding the velocity magnitude. The flow angle is manually determined from the B-mode image. The approach can be used for both tissue and blood velocity determination. The approach was investigated using both simulations and a flow system with a laminar flow. The flow profile was measured with a commercial 7.5 MHz linear array transducer. A plastic tube with an internal diameter of 17 mm was used with an EcoWatt 1 pump generating a laminar, stationary flow. The velocity profile was measured for flow angles of 90 and 60 degrees. The RASMUS research scanner was used for acquiring radio frequency (RF) data from 128 elements of the array, using 8 emissions with 11 elements in each emission. A 20-micros chirp was used during emission. The RF data were subsequently beamformed off-line and stationary echo canceling was performed. The 60-degree flow with a peak velocity of 0.15 m/s was determined using 16 groups of 8 emissions, and the relative standard deviation was 0.36% (0.65 mm/s). Using the same setup for purely transverse flow gave a standard deviation of 1.2% (2.1 mm/s). Variation of the different parameters revealed the sensitivity to number of lines, angle deviations, length of correlation interval, and sampling interval. An in vivo image of the carotid artery and jugular vein of a healthy 29-year-old volunteer was acquired. A full color flow image using only 128 emissions could be made with a high-velocity precision.

In-vivo synthetic aperture flow imaging in medical ultrasound.

A new method for acquiring flow images using synthetic aperture techniques in medical ultrasound is presented. The new approach makes it possible to have a continuous acquisition of flow data throughout the whole image simultaneously, and this can significantly improve blood velocity estimation. Any type of filter can be used for discrimination between tissue and blood flow without initialization, and the number of lines used for velocity estimation is limited only by the nonstationarity of the flow. The new approach is investigated through both simulations and measurements. A flow rig is used for generating a parabolic laminar flow, and a research scanner is used for acquiring RF data from individual transducer elements. A reference profile is calculated from a mass flow meter. The parabolic velocity profile is estimated using the new approach with a relative standard deviation of 2.2% and a mean relative bias of 3.4% using 24 pulse emissions at a flow angle of 45 degrees. The 24 emissions can be used for making a full-color flow map image. An in-vivo image of flow in the carotid artery for a 29-year-old male also is presented. The full image is acquired using 24 emissions.

Real time 3D visualization of ultrasonic data using a standard PC.

This paper describes a flexible, software-based scan converter capable of rendering 3D volumetric data in real time on a standard PC. The display system is used in the remotely accessible and software-configurable multichannel ultrasound sampling system (RASMUS system) developed at the Center for Fast Ultrasound Imaging. The display system is split into two modules: data transfer and display. These two modules are independent and communicate using shared memory and a predefined set of functions. It is, thus, possible to use the display program with a different data-transfer module which is tailored to another source of data (scanner, database, etc.). The data-transfer module of the RASMUS system is based on a digital signal processor from Analog Devices--ADSP 21060. The beamformer is connected to a PC via the link channels of the ADSP. A direct memory access channel transfers the data from the ADSP to a memory buffer. The display module, which is based on OpenGL, uses this memory buffer as a texture map that is passed to the graphics board. The scan conversion, image interpolation, and logarithmic compression are performed by the graphics board, thus reducing the load on the main processor to a minimum. The scan conversion is done by mapping the ultrasonic data to polygons. The format of the image is determined only by the coordinates of the polygons allowing for any kind of geometry to be displayed on the screen. Data from color flow mapping is added by alpha-blending. The 3D data are displayed either as cross-sectional planes, or as a fully rendered 3D volume displayed as a pyramid. All sides of the pyramid can be changed to reveal B-mode or C-mode scans, and the pyramid can be rotated in all directions in real time.