Adaptive multi-element synthetic aperture imaging with motion and phase aberration correction

Improved Spatiotemporal Resolution in Echocardiography Using Mixed Geometry Imaging Sequences

Cardiac ultrasound seeks to image the most dynamic environment in the body-the moving heart. Many modern ultrasound imaging techniques address the tradeoff between spatial and temporal resolution using either narrow focused beams or with broad beam, synthetic aperture (SA) sequences that have been shown to suffer from motion artifacts. Retrospective encoding for conventional ultrasound sequences (REFoCUS) unifies the processing of these various geometric sequences, but the motion sensitivity of this approach has yet to be investigated. We hypothesize that a "mixed sequence" enabled by the REFoCUS method incorporating several beam geometries may better resolve cardiac motion over a wide field of view (FOV) and at a high frame rate. First, the motion sensitivity of REFoCUS was evaluated in simulation for several focused and broad transmit profiles. Focused transmissions resolve both lateral and axial motion much more effectively than broad transmissions, with performance similar to conventional beamforming techniques. Second, a mixed sequence was designed that insonifies the full field-of-view with plane wave (PW) transmissions and key moving targets with focused transmissions. This mixed sequence was tested in simulation and in vivo and was used to image the heart as well as the liver, a low-motion control. By combining a sparse PW sequence ( n=60 ) with a small group of targeted focused transmissions ( n = 10), the anterior mitral valve leaflet (AML) at its peak observed velocity was better resolved. We believe that mixed sequences have strong potential to resolve cardiac motion at clinically relevant frame rates.

Implementation of constrained swept synthetic aperture using a mechanical fixture

Resolution and target detectability in ultrasound imaging are directly tied to the size of the imaging array. This is especially important for imaging at depth, such as in the detection and diagnosis of hepatocellular carcinoma and other lesions in the liver. Swept synthetic aperture (SSA) imaging has shown promise for building large effective apertures from small physical arrays using motion, but has required bulky fixtures and external motion tracking for precise positioning. In this study we present an approach that constrains the transducer motion with a simple linear sliding fixture and estimates motion from the ultrasound data itself using either speckle tracking or channel correlation. We demonstrate in simulation and phantom experiments the ability of both techniques to accurately estimate lateral transducer motion and form SSA images with improved resolution and target detectability. We observed errors under 83 μm across a 50 mm sweep in simulation and found improvements of up to 61% in resolution and up to 33% in lesion detectability experimentally even imaging through ex vivo tissue layers. This approach will increase the accessibility of SSA imaging and allow us to test its use in clinical settings.

High-Frequency Endoscopic Ultrasound Imaging With Phase-Corrected-and-Sum and Coherence Factor Weighting

High-frequency endoscopic ultrasound (HFEUS) imaging is an important tool commonly used in clinical practice for imaging hollow organs. The virtual source synthetic aperture (VSSA) method is effective in improving the imaging quality of HFEUS. However, interference from the motor control unit severely affects the accuracy of the conventional delay and sum (DAS) method, thus compromising the effectiveness of VSSA. In this article, a new computational method based on phase correction was proposed to overcome these shortcomings, which is named phase-corrected-and-sum (PCAS). Meanwhile, the parameters of coherence factor weighting (CFW) can be obtained from the correlation coefficient of the superimposed signals to further increase the imaging quality. Three kinds of imaging experiments were designed to evaluate the proposed method. Compared with the conventional method, the results show that the PCAS-CFW method improves the lateral resolution by about 10% and the contrast-to-noise ratio (CNR) by about 44%. Therefore, this proposed method is capable of significantly improving HFEUS image quality, and this method can be easily integrated into current HFEUS imaging systems, showing great potential for clinical applications.

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.

Extending Retrospective Encoding for Robust Recovery of the Multistatic Data Set

Robust recovery of multistatic synthetic aperture data from conventional ultrasound sequences can enable complete transmit-and-receive focusing at all points in the field of view without the drawbacks of virtual-source synthetic aperture and further enables more advanced imaging applications, such as backscatter coherence, sound speed estimation, and phase aberration correction. Recovery of the multistatic data set has previously been demonstrated on a steered transmit sequence for phased arrays using an adjoint-based method. We introduce two methods to improve the accuracy of the multistatic data set. We first modify the original technique used for steered transmit sequences by applying a ramp filter to compensate for the nonuniform frequency scaling introduced by the adjoint-based method. Then, we present a regularized inversion technique that allows additional aperture specification and is intended to work for both steered transmit and walking aperture sequences. The ramp-filtered adjoint and regularized inversion techniques, respectively, improve the correlation of the recovered signal with the ground truth from 0.9404 to 0.9774 and 0.9894 in steered transmit sequences and 0.4610 to 0.4733 and 0.9936 in walking aperture sequences.

Effects of motion on correlations of pulse-echo ultrasound signals: Applications in delay estimation and aperture coherence

The correlation between two pulse-echo ultrasound signals is used to achieve a wide range of ultrasound techniques, such as Doppler imaging and elastography. Prior theoretical descriptions of pulse-echo correlations were restricted to stationary scatterers. Here, a theory for the correlation of moving scatterers is presented. An expression is derived for the correlation of two pulse-echo signals with arbitrary transmit and receive apertures acquired from a medium undergoing bulk motion using the Fresnel approximation. The derivation is shown to coincide with prior derivations in the absence of scatterer motion. The theory was compared against simulations in applications of phase-shift estimation and aperture coherence measurements. The phase-shift estimate and jitter were accurately predicted under axial and transverse motion for focused transmit apertures and for sequential and interleaved synthetic transmit apertures. The theory also accurately predicted how motion affects the correlation coefficient between receive aperture elements for a synthetic transmit aperture. The presented theory provides a framework for analyzing the correlations of arbitrary pulse-echo configurations for applications in which scatterer motion is expected.

Doppler-Based Motion Compensation Strategies for 3-D Diverging Wave Compounding and Multiplane-Transmit Beamforming: A Simulation Study

Fast imaging of the heart has shown promise toward bringing new diagnostic information. Although most studies to date have been based on 2-D imaging technology, the ultimate diagnostic tool would enable fast 3-D echocardiography. Hereto, 3-D diverging wave compounding (DWC) and 3-D multiline-transmit (MLT) beamforming have recently been proposed. Moreover, in our recent study, a hybrid technique was proposed in which multiple planar diverging waves were transmitted [i.e., multiplane-transmit (MPT)]. The proposed 3MPT sequence was demonstrated to outperform $9 \times 9$ DWC and 16MLT-4MLA (i.e., multiline acquisition) while imaging moving targets. However, none of the investigated beamforming techniques made use of motion compensation (MoCo) strategies. In this paper, we therefore propose Doppler-based MoCo strategies for 3-D DWC and MPT and test them via computer simulations. It is demonstrated that the MoCo strategies proposed for both DWC and MPT are effective and significantly restore image quality. Moreover, the MPT beamforming with MoCo outperforms $9 \times 9$ DWC with MoCo in terms of contrast ratio and contrast-to-noise ratio. The proposed MPT beamforming with MoCo thus provides volumetric images with relatively high temporal resolution (~66 Hz) and high image quality that is minimally affected by motion artifacts.

Compressed Sensing Based Synthetic Transmit Aperture for Phased Array Using Hadamard Encoded Diverging Wave Transmissions

Previously, we proposed compressed sensing based synthetic transmit aperture (CS-STA) to improve the contrast and frame rate of STA while maintaining its spatial resolution in linear array by choosing uniform random matrix as the measurement matrix and transmitting the plane waves (PWs). In this paper, to extend CS-STA for phased array imaging and further improve its performance, we design four types of CS-STA implementations with different combinations of measurement matrices (i.e., uniform random and Hadamard matrices) and transmitted waves [i.e., PW and diverging wave (DW)]. Through simulations and phantom experiments with a 3 MHz, 64-element phased array, we find that type-IV CS-STA with the combination of a Hadamard matrix and DW outperforms the other three implementations including the previously proposed type-I CS-STA in terms of image quality and reconstruction time. Specifically, PW transmission produces visible discontinuity and the reconstruction time with uniform random matrix is about 100-fold longer than that with the Hadamard matrix. Compared with STA, with eightfold higher frame rate, type-IV CS-STA achieves 8.2 and 12.3 dB higher contrast-to-noise ratio and signal-to-noise ratio in the simulations, respectively. These improvements are slightly lower in the phantom experiments, which are 6.2 and 6.6 dB, respectively. In addition, CS-STA does not deteriorate the spatial resolution of STA, with the maximum deterioration being smaller than 1/8 wavelength. These results demonstrate that type-IV CS-STA can achieve phased array imaging with high image quality at high frame rate and may be beneficial to cardiac imaging.

Recovery of the Complete Data Set From Focused Transmit Beams

The focused transmit beam is a standard tool for clinical ultrasound scanning, concentrating energy from a number of array elements toward an imaging target. However, above and below the transmit focus, much of the energy in the beam is spread in a broadened main lobe and long off-axis tails that are ignored by conventional beamforming methods. This paper proposes a method to decompose a set of focused transmit beams into their constituent components-diverging waves from individual array elements. The recovery of this complete data set enables synthetic transmit focusing at all points in the field of view without beam shape or focal depth artifacts commonly associated with virtual source synthetic aperture beamforming. An efficient frequency-domain linear decoding implementation is introduced. The principles of the method are demonstrated both in transmit field simulations and experimental imaging. At depth, up to a 9.6-dB improvement in electronic signal-to-noise ratio and 8.9-dB improvement in contrast were observed in comparison with conventional dynamic receive beamforming. The proposed method is broadly applicable to existing scan sequences and requires only channel data for processing.

Feasibility of Multiplane-Transmit Beamforming for Real-Time Volumetric Cardiac Imaging: A Simulation Study

Today's 3-D cardiac ultrasound imaging systems suffer from relatively low spatial and temporal resolution, limiting their applicability in daily clinical practice. To address this problem, 3-D diverging wave imaging with spatial coherent compounding (DWC) as well as 3-D multiline-transmit (MLT) imaging have recently been proposed. Currently, the former improves the temporal resolution significantly at the expense of image quality and the risk of introducing motion artifacts, whereas the latter only provides a moderate gain in volume rate but mostly preserves quality. In this paper, a new technique for real-time volumetric cardiac imaging is proposed by combining the strengths of both approaches. Hereto, multiple planar (i.e., 2-D) diverging waves are simultaneously transmitted in order to scan the 3-D volume, i.e., multiplane transmit (MPT) beamforming. The performance of a 3MPT imaging system was contrasted to that of a 3-D DWC system and that of a 3-D MLT system by computer simulations during both static and moving conditions of the target structures while operating at similar volume rate. It was demonstrated that for stationary targets, the 3MPT imaging system was competitive with both the 3-D DWC and 3-D MLT systems in terms of spatial resolution and sidelobe levels (i.e., image quality). However, for moving targets, the image quality quickly deteriorated for the 3-D DWC systems while it remained stable for the 3MPT system while operating at twice the volume rate of the 3-D-MLT system. The proposed MPT beamforming approach was thus demonstrated to be feasible and competitive to state-of-the-art methodologies.

Effects of line fiducial parameters and beamforming on ultrasound calibration

Ultrasound (US)-guided interventions are often enhanced via integration with an augmented reality environment, a necessary component of which is US calibration. Calibration requires the segmentation of fiducials, i.e., a phantom, in US images. Fiducial localization error (FLE) can decrease US calibration accuracy, which fundamentally affects the total accuracy of the interventional guidance system. Here, we investigate the effects of US image reconstruction techniques as well as phantom material and geometry on US calibration. It was shown that the FLE was reduced by 29% with synthetic transmit aperture imaging compared with conventional B-mode imaging in a Z-bar calibration, resulting in a 10% reduction of calibration error. In addition, an evaluation of a variety of calibration phantoms with different geometrical and material properties was performed. The phantoms included braided wire, plastic straws, and polyvinyl alcohol cryogel tubes with different diameters. It was shown that these properties have a significant effect on calibration error, which is a variable based on US beamforming techniques. These results would have important implications for calibration procedures and their feasibility in the context of image-guided procedures.

Error analysis of subaperture processing in 1-D ultrasound arrays

To simplify the medical ultrasound system and reduce the cost, several techniques have been proposed to reduce the interconnections between the ultrasound probe and the back-end console. Among them, subaperture processing (SAP) is the most straightforward approach and is widely used in commercial products. This paper reviews the most important error sources of SAP, such as static focusing, delay quantization, linear delay profile, and coarse apodization, and the impacts introduced by these errors are shown. We propose to use main lobe coherence loss as a simple classification of the quality of the beam profile for a given design. This figure-ofmerit (FoM) is evaluated by simulations with a 1-D ultrasound subaperture array setup. The analytical expressions and the coherence loss can work as a quick guideline in subaperture design by equalizing the merit degradations from different error sources, as well as minimizing the average or maximum loss over ranges. For the evaluated 1-D array example, a good balance between errors and cost was achieved using a subaperture size of 5 elements, focus at 40 mm range, and a delay quantization step corresponding to a phase of π/4.

Comparison of 3-D synthetic aperture phased-array ultrasound imaging and parallel beamforming

This paper demonstrates that synthetic aperture imaging (SAI) can be used to achieve real-time 3-D ultrasound phased-array imaging. It investigates whether SAI increases the image quality compared with the parallel beamforming (PB) technique for real-time 3-D imaging. Data are obtained using both simulations and measurements with an ultrasound research scanner and a commercially available 3.5- MHz 1024-element 2-D transducer array. To limit the probe cable thickness, 256 active elements are used in transmit and receive for both techniques. The two imaging techniques were designed for cardiac imaging, which requires sequences designed for imaging down to 15 cm of depth and a frame rate of at least 20 Hz. The imaging quality of the two techniques is investigated through simulations as a function of depth and angle. SAI improved the full-width at half-maximum (FWHM) at low steering angles by 35%, and the 20-dB cystic resolution by up to 62%. The FWHM of the measured line spread function (LSF) at 80 mm depth showed a difference of 20% in favor of SAI. SAI reduced the cyst radius at 60 mm depth by 39% in measurements. SAI improved the contrast-to-noise ratio measured on anechoic cysts embedded in a tissue-mimicking material by 29% at 70 mm depth. The estimated penetration depth on the same tissue-mimicking phantom shows that SAI increased the penetration by 24% compared with PB. Neither SAI nor PB achieved the design goal of 15 cm penetration depth. This is likely due to the limited transducer surface area and a low SNR of the experimental scanner used.

Coherent plane wave compounding for very high frame rate ultrasonography of rapidly moving targets

Coherent plane wave compounding is a promising technique for achieving very high frame rate imaging without compromising image quality or penetration. However, this approach relies on the hypothesis that the imaged object is not moving during the compounded scan sequence, which is not the case in cardiovascular imaging. This work investigates the effect of tissue motion on retrospective transmit focusing in coherent compounded plane wave imaging (PWI). Two compound scan sequences were studied based on a linear and alternating sequence of tilted plane waves, with different timing characteristics. Simulation studies revealed potentially severe degradations in the retrospective focusing process, where both radial and lateral resolution was reduced, lateral shifts of the imaged medium were introduced, and losses in signal-to-noise ratio (SNR) were inferred. For myocardial imaging, physiological tissue displacements were on the order of half a wavelength, leading to SNR losses up to 35 dB, and reductions of contrast by 40 dB. No significant difference was observed between the different tilt sequences. A motion compensation technique based on cross-correlation was introduced, which significantly recovered the losses in SNR and contrast for physiological tissue velocities. Worst case losses in SNR and contrast were recovered by 35 dB and 27-35 dB, respectively. The effects of motion were demonstrated in vivo when imaging a rat heart. Using PWI, very high frame rates up to 463 fps were achieved at high image quality, but a motion correction scheme was then required.

Sequential beamforming for synthetic aperture imaging

Synthetic aperture sequential beamforming (SASB) is a novel technique which allows to implement synthetic aperture beamforming on a system with a restricted complexity, and without storing RF-data. The objective is to improve lateral resolution and obtain a more depth independent resolution compared to conventional ultrasound imaging. SASB is a two-stage procedure using two separate beamformers. The initial step is to construct and store a set of B-mode image lines using a single focal point in both transmit and receive. The focal points are considered virtual sources and virtual receivers making up a virtual array. The second stage applies the focused image lines from the first stage as input data, and take advantage of the virtual array in the delay and sum beamforming. The size of the virtual array is dynamically expanded and the image is dynamically focused in both transmit and receive and a range independent lateral resolution is obtained. The SASB method has been investigated using simulations in Field II and by off-line processing of data acquired with a commercial scanner. The lateral resolution increases with a decreasing F#. Grating lobes appear if F#≤2 for a linear array with λ-pitch. The performance of SASB with the virtual source at 20mm and F#=1.5 is compared with conventional dynamic receive focusing (DRF). The axial resolution is the same for the two methods. For the lateral resolution there is improvement in FWHM of at least a factor of 2 and the improvement at -40dB is at least a factor of 3. With SASB the resolution is almost constant throughout the range. For DRF the FWHM increases almost linearly with range and the resolution at -40dB is fluctuating with range. The theoretical potential improvement in SNR of SASB over DRF has been estimated. An improvement is attained at the entire range, and at a depth of 80mm the improvement is 8dB.

In vivo evaluation of synthetic aperture sequential beamforming

Ultrasound in vivo imaging using synthetic aperture sequential beamformation (SASB) is compared with conventional imaging in a double blinded study using side-by-side comparisons. The objective is to evaluate if the image quality in terms of penetration depth, spatial resolution, contrast and unwanted artifacts is comparable to conventional imaging. In vivo data was acquired using a ProFocus ultrasound scanner (BK Medical, Herlev, Denmark) and a 192-element 3.5 MHz convex array transducer (Sound Technology Inc., PA, USA). Data were acquired interleaved, ensuring that the exact same anatomical locations were scanned. Eighteen volunteers were scanned abdominally resulting in 85 image sequence pairs. Evaluation of image quality was performed by five medical doctors. Results show that image quality using SASB was significantly better than conventional imaging (p value: <0.01). There was not a significant difference in penetration depth (p value: 0.55). The study supports that in vivo ultrasound imaging using SASB is feasible for abdominal imaging.

Coded tissue harmonic imaging with nonlinear chirp signals

Coded tissue harmonic imaging with pulse inversion (CTHI-PI) based on a linear chirp signal can improve the signal-to-noise ratio with minimizing the peak range sidelobe level (PRSL), which is the main advantage over CTHI with bandpass filtering (CTHI-BF). However, the CTHI-PI technique could suffer from motion artifacts due to decreasing frame rate caused by two firings of opposite phase signals for each scanline. In this paper, a new CTHI method based on a nonlinear chirp signal (CTHI-NC) is presented, which can improve the separation of fundamental and harmonic components without sacrificing frame rate. The nonlinear chirp signal is designed to minimize the PRSL value by optimizing its frequency sweep rate and time duration. The performance of the CTHI-NC method was evaluated by measuring the PRSL and mainlobe width after compression. From the in vitro experiments, the CTHI-NC provided the PRSL of -40.6 dB and the mainlobe width of 2.1 μs for the transmit quadratic nonlinear chirp signal with the center frequency of 2.1 MHz, the fractional bandwidth at -6 dB of 0.6 and the time duration of 15 μs. These results indicate that the proposed method could be used for improving frame rates in CTHI while providing comparable image quality to CTHI-PI.

Least-squares estimation of imaging parameters for an ultrasonic array using known geometric image features

Ultrasonic array images are adversely affected by errors in the assumed or measured imaging parameters. For non-destructive testing and evaluation, this can result in reduced defect detection and characterization performance. In this paper, an autofocus algorithm is presented for estimating and correcting imaging parameter errors using the collected echo data and a priori knowledge of the image geometry. Focusing is achieved by isolating a known geometric feature in the collected data and then performing a weighted leastsquares minimization of the errors between the data and a feature model, with respect to the unknown parameters. The autofocus algorithm is described for the estimation of element positions in a flexible array coupled to a specimen with an unknown surface profile. Experimental results are shown using a prototype flexible array and it is demonstrated that (for an isolated feature and a well-prescribed feature model) the algorithm is capable of generating autofocused images that are comparable in quality to benchmark images generated using accurately known imaging parameters.

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.

Coarray synthesis based on polynomial decomposition

Synthetic aperture (SA) techniques have been frequently used to reduce the volume and complexity of the imaging systems. A useful tool for designing synthetic aperture configurations is the coarray. This is the virtual aperture that produces in one way the same beam pattern as the SA system in emission and reception. In this correspondence, we propose a new algorithm, based on the polynomial decomposition, that allows to obtain any wanted coarray on a linear array using whatever synthetic aperture configuration. With this fast and simple algorithm, the desired coarray is decomposed into a set of sub-apertures, whose length is determined by the requirements and resources of the system. The result is the set of weights that have to be applied on the sub-apertures to get the desired coarray, and consequently, the wanted beam pattern.

Minimally redundant 2-D array designs for 3-D medical ultrasound imaging

In real-time ultrasonic 3-D imaging, in addition to difficulties in fabricating and interconnecting 2-D transducer arrays with hundreds of elements, there are also challenges in acquiring and processing data from a large number of ultrasound channels. The coarray (spatial convolution of the transmit and receive arrays) can be used to find efficient array designs that capture all of the spatial frequency content (a transmit-receive element combination corresponds to a spatial frequency) with a reduced number of active channels and firing events. Eliminating the redundancies in the transmit-receive element combinations and firing events reduces the overall system complexity and improves the frame rate. Here we explore four reduced redundancy 2-D array configurations for miniature 3-D ultrasonic imaging systems. Our approach is based on 1) coarray design with reduced redundancy using different subsets of linear arrays constituting the 2-D transducer array, and 2) 3-D scanning using fan-beams (narrow in one dimension and broad in the other dimension) generated by the transmit linear arrays. We form the overall array response through coherent summation of the individual responses of each transmit-receive array pairs. We present theoretical and simulated point spread functions of the array configurations along with quantitative comparison in terms of the front-end complexity and image quality.

Timing-error-difference calibration using reciprocal signals

Timing errors in transmission and reception electronic channels of medical ultrasound imaging systems are generally smaller than one-tenth of a wavelength and do not influence the focusing quality of the system. However, these errors influence the performance of the near-field-signal-redundancy algorithm for correcting phase-aberrations generated by speed heterogeneity in the medium due to its high sensitivity to errors. The effect of timing errors is to make the transmission and reception phase-aberration profiles different. When the difference is much smaller than the period of the signal, an algorithm has been proposed in a previous work to measure the average of the transmission and reception phase-aberration profiles, and it can be used as an approximation to correct phase-aberrations on both transmission and reception. However, when the difference is large, the transmission and reception phase-aberration profiles need to be measured separately. In this paper, several algorithms that use reciprocal signals are proposed to measure the difference profile of the transmission and reception phase-aberration profiles. Their performances are theoretically analyzed, simulated, and experimentally tested. From the measured average and difference profiles, the transmission and reception phase-aberration profiles can be derived separately and used to correct phase-aberrations on transmission and reception, respectively.

Real-time 3D imaging methods using 2D phased arrays based on synthetic focusing techniques

A fast 3D ultrasound imaging technique using a 2D phased array transducer based on the synthetic focusing method for nondestructive testing or medical imaging is proposed. In the proposed method, each column of a 2D array is fired successively to produce transverse fan beams focused at a fixed depth along a given longitudinal direction and the resulting pulse echoes are received at all elements of a 2D array used. After firing all column arrays, a frame of high-resolution image along a given longitudinal direction is obtained with dynamic focusing employed in the longitudinal direction on receive and in the transverse direction on both transmit and receive. The volume rate of the proposed method can be increased much higher than that of the conventional 2D array imaging by employing an efficient sparse array technique. A simple modification to the proposed method can further increase the volume scan rate significantly. The proposed methods are verified through computer simulations.

The cross algorithm for phase-aberration correction in medical ultrasound images formed with two-dimensional arrays

Common-midpoint signals in the near-field signal-redundancy (NFSR) algorithm for one-dimensional arrays are acquired using three consecutive transducer elements. An all-row-plus-two-column algorithm has been proposed to implement the one-dimensional NFSR algorithm on two dimensional arrays. The disadvantage of this method is that its ambiguity profile is not linear and a timeconsuming iterative method has to be used to linearize the ambiguity profile. An all-row-plus-two-column-and-a-diagonal algorithm has also been proposed. Its ambiguity profile is linear, but it is very sensitive to noise and cannot be used. In this paper, a novel cross algorithm is proposed to implement the NFSR algorithm on two-dimensional arrays. In this algorithm, common-midpoint signals are acquired using four adjacent transducer elements, which is not available in one-dimensional arrays. Its advantage includes a linear ambiguity profile and a higher measurement signal-to-noise ratio. The performance of the cross algorithm is evaluated theoretically. The region of redundancy is analyzed. The procedure for deriving the phaseaberration profile from peak positions of cross-correlation functions between common-midpoint signals is discussed. This algorithm is tested with a simulated data set acquired with a two-dimensional array, and the result shows that the cross algorithm performs better than the all-row plus-twocolumn NFSR algorithm.

Effects influencing focusing in synthetic aperture vector flow imaging

Previously, a synthetic aperture vector velocity estimation method was proposed. Data are beamformed at different directions through a point, where the velocity is estimated. The flow direction is estimated by a search for the direction where the normalized cross-correlation peaks and the velocity magnitude along this direction are found. In this paper, different effects that influence the focusing in this method are investigated. These include the effect of phase errors in the emitted spherical waves, motion effects, and the effect of various interpolation methods in beam-forming. A model based on amplitude drop and phase error for spherical waves created using the virtual source concept is derived. This model can be used to determine the opening angle of a virtual source. Simulations for different virtual source placements are made, and it is recommended that the virtual sources be placed behind the aperture when shallow structures are imaged, and when deeper-lying structures are imaged the virtual sources be placed in front of the aperture. Synthetic aperture methods involve summation of data from numerous emissions. Motion between these emissions results in incoherence and affects resolution, contrast, and the signal-to-noise ratio. The effects of motion on the synthetic aperture vector velocity estimation method are investigated, and it is shown that for both axial and lateral motion, the contrast and signal-to-noise ratio can be seriously affected. A compensation method using the previous vector velocity estimate, when new data are beamformed, is implemented and tested. It is shown from a number of flow phantom experiments that a significant improvement with respect to bias and standard deviation of the velocity estimates can be obtained by using this compensation. Increased performance is gained at the expense of computation time. Different interpolation methods can be used for beam-forming the data. In this paper, the velocity estimation performance using various more complex interpolation schemes are compared to that using linear interpolation. No significant difference in the performance of the method is seen when other interpolation methods are used.

Motion artifacts of extended high frame rate imaging

Based on the high frame rate (HFR) imaging method developed in our lab, an extended high frame rate imaging method with various transmission schemes was developed recently. In this method, multiple, limited-diffraction array beams or steered plane wave transmissions are used to increase image resolution and field of view as well as to reduce sidelobes. Furthermore, the multiple, limited-diffraction array beam transmissions can be approximated with square-wave aperture weightings, allowing one or two transmitters to be used with a multielement array transducer to simplify imaging systems. By varying the number of transmissions, the extended HFR imaging method allows a continuous trade-off between image quality and frame rate. Because multiple transmissions are needed to obtain one frame of image for the method, motion could cause phase misalignment and thus produce artifacts, reducing image contrast and resolution and leading to an inaccurate clinical interpretation of images. Therefore, it is important to study how motion affects the method and provide a useful guidance of using the method properly in various applications. In this paper, computer simulations, in vitro and in vivo experiments were performed to study the effects of motion on the method in different conditions. Results show that a number of factors may affect the motion effects. However, it was found that the extended HFR imaging method is not sensitive to the motions commonly encountered in the clinical applications, as is demonstrated by an in vivo heart experiment, unless the number of transmissions is large and objects are moving at a high velocity near the surface of a transducer.

In-vivo evaluation of convex array synthetic aperture imaging

This paper presents an in-vivo study of synthetic transmit aperture (STA) imaging in comparison with conventional imaging, evaluating whether STA imaging is feasible in-vivo and whether the image quality obtained is comparable with traditional scanned imaging in terms of penetration depth, spatial resolution, contrast resolution and artifacts. Acquisition was performed using our research scanner RASMUS and a 5.5 MHz convex array transducer. STA imaging was acquired using circular wave emulation by 33-element subapertures and a 20 micros linear FM signal as excitation pulse. For conventional imaging, a 64 element aperture was used in transmit and receive with a 1.5 cycle sinusoid excitation pulse. Conventional and STA images were acquired interleaved, ensuring that the exact same anatomical location was scanned. Image sequences were recorded in real time and processed off-line. Seven male volunteers were scanned abdominally and the resulting images were compared by three medical doctors using randomized blinded presentation. Penetration and image quality were scored and evaluated statistically. Results showed slightly but significantly (0.48 cm, p = 0.008) increased penetration using STA. Image quality was also highly significantly (p < 0.001) increased. Results show that in-vivo ultrasound imaging using STA is feasible for abdominal imaging without severe motion artifacts.

Implementation of the near-field signal redundancy phase-aberration correction algorithm on two-dimensional arrays

Near-field signal-redundancy (NFSR) algorithms for phase-aberration correction have been proposed and experimentally tested for linear and phased one-dimensional arrays. In this paper the performance of an all-row-plus-two-column, two-dimensional algorithm has been analyzed and tested with simulated data sets. This algorithm applies the NFSR algorithm for one-dimensional arrays to all the rows as well as the first and last columns of the array. The results from the two column measurements are used to derive a linear term for each row measurement result. These linear terms then are incorporated into the row results to obtain a two-dimensional phase aberration profile. The ambiguity phase aberration profile, which is the difference between the true and the derived phase aberration profiles, of this algorithm is not linear. Two methods, a trial-and-error method and a diagonal-measurement method, are proposed to linearize the ambiguity profile. The performance of these algorithms is analyzed and tested with simulated data sets.

Position and time-delay calibration of transducer elements in a sparse array for underwater ultrasound imaging

This paper describes a novel method for the calibration of the position and time delay of transducer elements in a large, sparse array used for underwater, high-resolution ultrasound imaging. This method is based on the principles used in the global positioning system (GPS). However, unlike GPS, in which the wave propagation speed is generally assumed known, the sound propagation speed in the water usually is unknown and it is calibrated simultaneously in this method to achieve high calibration accuracy. In this method, a high-precision positioning system is used to scan a single hydrophone (used for transmission) in the imaging field of the array. The hydrophone transmits pulses at selected positions, and the transducer elements in the sparse array receive the transmitted signals. Time of flight (TOF) values between transducer elements and hydrophone positions then are measured. From a series of measured TOF values, the position and time delay values for each transducer element as well as the propagation speed can be calibrated. The performances of the calibration algorithm are theoretically analyzed and evaluated with numerical calculations and simulation studies. It is found that this method is capable of calibrating the positions and time delays of transducer elements with high accuracy.

2-D array for 3-D ultrasound imaging using synthetic aperture techniques

A two-dimensional (2-D) array of 256 X 256 = 65,536 elements, with total area 4 X 4 = 16 cm2, serves as a flexible platform for developing acquisition schemes for 3-D rectilinear ultrasound imaging at 10 MHz using synthetic aperture techniques. This innovative system combines a simplified interconnect scheme and synthetic aperture techniques with a 2-D array for 3-D imaging. A row-column addressing scheme is used to access different elements for different transmit events. This addressing scheme is achieved through a simple interconnect, consisting of one top, one bottom single-layer, flex circuits that, compared to multilayer flex circuits, are simpler to design, cheaper to manufacture, and thinner so their effect on the acoustic response is minimized. We present three designs that prioritize different design objectives: volume acquisiton time, resolution, and sensitivity, while maintaining acceptable figures for the other design objectives. For example, one design overlooks time-acquisition requirements, assumes good noise conditions, and optimizes for resolution, achieving -6 dB and -20 dB beamwidths of less than 0.2 and 0.5 mm, respectively, for an F/2 aperture. Another design can acquire an entire volume in 256 transmit events, with -6 dB and -20 dB beamwidths in the order of 0.4 and 0.8 mm, respectively.

Coherent-array imaging using phased subarrays. Part I: basic principles

The front-end hardware complexity of a coherent array imaging system scales with the number of active array elements that are simultaneously used for transmission or reception of signals. Different imaging methods use different numbers of active channels and data collection strategies. Conventional full phased array (FPA) imaging produces the best image quality using all elements for both transmission and reception, and it has high front-end hardware complexity. In contrast, classical synthetic aperture (CSA) imaging only transmits on and receives from a single element at a time, minimizing the hardware complexity but achieving poor image quality. We propose a new coherent array imaging method--phased subarray (PSA) imaging--that performs partial transmit and receive beam-forming using a subset of adjacent elements at each firing step. This method reduces the number of active channels to the number of subarray elements; these channels are multiplexed across the full array and a reduced number of beams are acquired from each subarray. The low-resolution subarray images are laterally upsampled, interpolated, weighted, and coherently summed to form the final high-resolution PSA image. The PSA imaging reduces the complexity of the front-end hardware while achieving image quality approaching that of FPA imaging.

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.

An efficient motion estimation and compensation method for ultrasound synthetic aperture imaging

This paper describes a method for overcoming motion artifacts in synthetic aperture imaging. The method is based on a computer simulation study on the influence of target motion on synthetic aperture techniques. A region-based motion compensation approach is used in which only the axial motion is estimated and compensated for a given region of interest under the assumption that the whole ROI moves uniformly. The estimated axial motion is calculated with a crosscorrelation method at the point where the focused signal has the maximum energy within the ROI. We also present a method for estimating axial motion using the autocorrelation method that is widely used to estimate average Doppler frequency. Both computer simulations and in vivo experiments show that the proposed crosscorrelation-based method can greatly improve the spatial resolution and SNR of ultrasound imaging by implementing SA techniques for two-way dynamic focusing without motion artifacts. In addition, the autocorrelation-based motion compensation method provides almost the same results as the crosscorrelation-based method, but with a dramatically reduced computational complexity.

A study of synthetic-aperture imaging with virtual source elements in B-mode ultrasound imaging systems

We propose an all point transmit and receive focusing method based on transmit synthetic focusing combined with receive dynamic focusing in a linear array transducer. In the method, on transmit, a virtual source element is assumed to be located at the transmit focal depth of conventional B-mode imaging systems, and transmit synthetic focusing is used in two half planes, one before and the other after the transmit focal depth, using the RF data of each scanline, together with all other relevant RF scanline data previously stored. The proposed new method uses the same data acquisition scheme as the conventional focusing method while maintaining the same frame rate via high-speed signal processing, but it is not suitable for imaging moving objects. It improves upon the lateral resolution and sidelobe level at all imaging depths. Also, it increases the transmit power and image signal-to-noise ratio (SNR), due to transmit field synthesis, and extends the image penetration depth as well. Evaluations with simulation and experimental data show much improvement in resolution and SNR at all imaging depths.

Correlation processing for correction of phase distortions in subaperture imaging

Ultrasonic subaperture imaging combines synthetic aperture and phased array approaches and permits low-cost systems with improved image quality. In subaperture processing, a large array is synthesized using echo signals collected from a number of receive subapertures by multiple firings of a phased transmit subaperture. Tissue inhomogeneities and displacements in subaperture imaging may cause significant phase distortions on received echo signals. Correlation processing on reference echo signals can be used for correction of the phase distortions, for which the accuracy and robustness are critically limited by the signal correlation. In this study, we explore correlation processing techniques for adaptive subaperture imaging with phase correction for motion and tissue inhomogeneities. The proposed techniques use new subaperture data acquisition schemes to produce reference signal sets with improved signal correlation. The experimental test results were obtained using raw radio frequency (RF) data acquired from two different phantoms with 3.5 MHz, 128-element transducer array. The results show that phase distortions can effectively be compensated by the proposed techniques in real-time adaptive subaperture imaging.