Methods for Grating Lobe Suppression in Ultrasound Plane Wave Imaging

An efficient ultrasonic wavenumber-domain plane wave imaging method towards the inspection of curved structures

Ultrasonic phased array testing is commonly employed for inspecting curved structures. Conventional plane wave imaging techniques, based on delay-and-sum in the time-domain, offer high image quality and inspection accuracy but suffer from low frame rates due to their high computational complexity. In this work, an efficient wavenumber-domain imaging method that combines non-stationary wavefield extrapolation and f-k migration is proposed for curved structure inspection. Special emission focal laws are designed to generate a sequence of steered plane waves through the curved interface. The raw data is then extrapolated to the top boundary of the region of interest, followed by f-k migration to reconstruct images with high time efficiency. Simulation and experimental evaluations demonstrate a time reduction by a factor of up to 32.24 compared to conventional time-domain plane wave image reconstruction with equivalent image quality, highlighting its potential for monitoring flaws in real-time.

Nonlinear beamforming for intracardiac echocardiography: a comparative study

Intracardiac echocardiography (ICE) enables cardiac imaging with a wide field of view, deep imaging depth, and high frame rate during surgery. However, strong sidelobe and grating lobe artifacts created by the ultra-compact transducer degrade its image quality, making diagnosis and monitoring of treatment difficult. Conventionally, aperture apodization algorithms are often used to suppress sidelobe and grating lobe artifacts at the expense of lateral resolution, which is undesirable in ICE. In this study, we present comparative results of the beamforming methods specifically in ICE application. We demonstrate and compare five nonlinear beamforming algorithms in ICE: nonlinear pth root delay and sum (NL-p-DAS), nonlinear pth root spectral magnitude scaling (NL-p-SMS), delay-and-sum with coherence factors (DAS + SCF), delay and sum with apodization (DAS + apodization) and delay and sum (DAS). Phantom and ex-vivo experiment compare the performance of each algorithm in static and dynamic conditions. DAS + SCF shows the best lateral resolution, and all four algorithms improve the image contrast and sidelobe suppression over conventional DAS. NL-p-SMS stands out for the best axial resolution and suppression of grating lobe artifacts. For motion tracking, NL-p-SMS shows better temporal resolution than other methods. Overall, all the beamforming algorithms other than DAS showed improved image quality. Among them, NL-p-SMS, which has a high temporal resolution, showed the potential for providing more accurate information regards movement tracking.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13534-024-00352-9.

Frame rate improvement in ultrafast coherent plane wave compounding

Coherent plane wave compounding (CPWC), as an ultrafast ultrasound imaging technique, makes a significant breakthrough in frame rate enhancement. However, there exists a compromise between the quality of the final image and the frame rate in CPWC. In this paper, we propose an efficient method to minimize the number of required emissions, and consequently, improve the frame rate, while maintaining the image quality. To this end, we down-sample the angle interval using two specific sampling factors. More precisely, we construct two different subsets, each of which consists of a few numbers of emissions. The optimal values of the angle intervals are achieved based on the beampattern that corresponds to the reference case (that is, the case where all plane waves are used). Finally, in order to keep the image quality comparable with the reference case, we apply some modifications to the image reconstruction procedure. In the proposed algorithm, the Delay-and-Sum beamformed images of two considered subsets are convolved to achieve the final reconstructed image. The obtained results confirm the efficiency of the proposed method in terms of frame rate improvement compared to the reference case. In particular, by using the proposed method, the required emissions in PICMUS data reduce to 16, which is 4.6 times smaller compared to the reference case. Also, the gCNR values of the proposed method and the reference case are obtained as 0.98 and 0.97, respectively, for in-vivo dataset. This demonstrates that the proposed method successfully preserves the quality of the reconstructed image by using much fewer emissions.

Frequency-Dependent F-Number Suppresses Grating Lobes and Improves the Lateral Resolution in Coherent Plane-Wave Compounding

Ultrafast imaging modes, such as coherent plane-wave compounding (CPWC), increase image uniformity and reduce grating lobe artifacts by dynamic receive apertures. The focal length and the desired aperture width maintain a given ratio, which is called the F -number. Fixed F -numbers, however, exclude useful low-frequency components from the focusing and reduce the lateral resolution. Herein, this reduction is avoided by a frequency-dependent F -number. This F -number derives from the far-field directivity pattern of a focused aperture and can be expressed in closed form. The F -number, at low frequencies, widens the aperture to improve the lateral resolution. The F -number, at high frequencies, narrows the aperture to avoid lobe overlaps and suppress grating lobes. Phantom and in vivo experiments with a Fourier-domain beamforming algorithm validated the proposed F -number in CPWC. The lateral resolution, which was measured by the median lateral full-widths at half-maximum of wires, improved by up to 46.8% and 14.9% in a wire and a tissue phantom, respectively, in comparison to fixed F -numbers. Grating lobe artifacts, which were measured by the median peak signal-to-noise ratios of wires, reduced by up to 9.9 dB in comparison to the full aperture. The proposed F -number thus outperformed F -numbers that were recently derived from the directivity of the array elements.

On the physics of ultrasound transmission for in-plane needle tracking in guided interventions

Objective.In ultrasound (US) guided interventions, the accurate visualization and tracking of needles is a critical challenge, particularly during in-plane insertions. An inaccurate identification and localization of needles lead to severe inadvertent complications and increased procedure times. This is due to the inherent specular reflections from the needle with directivity depending on the angle of incidence of the US beam, and the needle inclination.Approach.Though several methods have been proposed for improved needle visualization, a detailed study emphasizing the physics of specular reflections resulting from the interaction of transmitted US beam with the needle remains to be explored. In this work, we discuss the properties of specular reflections from planar and spherical wave US transmissions respectively through multi-angle plane wave (PW) and synthetic transmit aperture (STA) techniques for in-plane needle insertion angles between 15°-50°.Main Results.The qualitative and quantitative results from simulations and experiments reveal that the spherical waves enable better visualization and characterization of needles than planar wavefronts. The needle visibility in PW transmissions is severely degraded by the receive aperture weighting during image reconstruction than STA due to greater deviation in reflection directivity. It is also observed that the spherical wave characteristics starts to alter to planar characteristics due to wave divergence at large needle insertion depths.Significance.The study highlights that synergistic transmit-receive imaging schemes addressing the physical properties of reflections from the transmit wavefronts are imperative for the precise imaging of needle interfaces and hence have strong potential in elevating the quality of outcomes from US guided interventional practices.

Grating Lobe Reduction in Plane-Wave Imaging With Angular Compounding Using Subtraction of Coherent Signals

Plane-wave imaging (PWI) with angular compounding has gained in popularity over recent years, because it provides high frame rates and good image properties. However, most linear arrays used in clinical practice have a pitch that is equal to than the wavelength of ultrasound. Hence, the presence of grating lobes is a concern for PWI using multiple transmit angles. The presence of grating lobes produces clutter in images and reduces the ability to observe tissue contrast. Techniques to reduce or eliminate the presence of grating lobes for PWI using multiple angles will result in improved image quality. Null subtraction imaging (NSI) is a nonlinear beamforming technique that has been explored for improving the lateral resolution of ultrasonic imaging. However, the apodization scheme used in NSI also eliminates or greatly reduces the presence of grating lobes. Imaging tasks using NSI were evaluated in simulations and physical experiments involving tissue-mimicking phantoms and rat tumors in vivo. Images created with NSI were compared with images created using traditional delay and sum (DAS) with Hann apodization and images created using a generalized coherence factor (GCF). NSI was observed to greatly reduce the presence of grating lobes in ultrasonic images, compared to DAS with Hann and GCF, while maintaining spatial resolution and contrast in the images. Therefore, NSI can provide a novel means of creating images using PWI with multiple steering angles on clinically available linear arrays while reducing the adverse effects associated with grating lobes.

Optimization of a random linear ultrasonic therapeutic array based on a genetic algorithm

Given their advantage of suppressing grating lobes, randomly arranged linear arrays have potential for use in ultrasonic treatment. The current work proposes a method based on genetic algorithm to optimize the random arrangement of array elements, so that the suppression effect of grating lobes can be significantly improved with reduced calculating time. The maximum and average kerfs of array elements are used as genes, and the ratio of the maximum to the secondary maximum sound pressure at the focal plane is used as the optimized target. Typically, the calculation requirements of the current method can be reduced to ∼ 25% of the traversing method. We further discuss how the kerf width, the effective ratio of element areas and the ratio of focal distance to array aperture affect the suppression of grating lobes. For a typical linear array with 32 elements (1-MHz operating frequency, 1.5-mm element width and 150-mm focal distance), the results suggest that the grating lobes are suppressed well when (1) the ratio of maximum width to average width of the element is between 5 and 8, (2) the ratio of the effective element area to the area of the whole array is between 0.5 and 0.9, and (3) the ratio of the effective emission aperture to the actual emission aperture of the array is as large as possible. Based on optimized parameters, an experimental array was fabricated and the measured results of corresponding sound field were entirely consistent with the simulated results (Given her role as an Associate Editor of this journal, Juan Tu had no involvement in the peer-review of articles for which she was an author and had no access to information regarding the peer-review. Full responsibility for the peer-review process for this article was delegated to another Editor of this journal.).

Experimental Study of Aperiodic Plane Wave Imaging for Ultrafast 3-D Ultrasound Imaging

OBJECTIVE

Although plane wave imaging (PWI) with multiple plane waves (PWs) steered at different angles enables ultrafast three-dimensional (3-D) ultrasonic imaging, there is still a challenging tradeoff between image quality and frame rate. To address this challenge, we recently proposed the aperiodic PWI (APWI) with mathematical analysis and simulation study. In this paper, we demonstrate the feasibility of APWI and evaluate the performance with phantom and in vivo experiments.

METHODS

APWI with a concentric ring angle pattern (APWI-C) and APWI with a sunflower pattern (APWI-S) are evaluated. For experimental verification of the methods, the experimental results are compared with simulation results in terms of the mainlobe-to-sidelobe ratio. In addition, the performance of APWI is compared with that of conventional PWI by using a commercial phantom. To examine the potential for clinical use of APWI, a gallstone-mimicking phantom study and an in vivo carotid artery experiment are also conducted.

RESULTS

In the phantom study, the APWI methods provide a contrast ratio approximately 23 dB higher than that of PWI. In a gallstone mimicking experiment, the proposed methods yield 3-D rendered stone images more similar to the real stones than PWI. In the in vivo carotid artery images, APWI reduces the clutter artifacts inside the artery.

CONCLUSION

Phantom and in vivo studies show that the APWI enhances the contrast without compromising the spatial resolution and frame rate.

SIGNIFICANCE

This study experimentally demonstrates the feasibility and advantage of APWI for ultrafast 3-D ultrasonic imaging.

Pixel-Oriented Adaptive Apodization for Plane-Wave Imaging Based on Recovery of the Complete Dataset

In theory, coherent plane-wave compounding (CPWC) enables ultrafast ultrasound imaging while maintaining a high imaging quality that is comparable to conventional B-mode imaging based on focused beam transmissions. However, in practice, due to the imperfect synthetization of transmit focusing (e.g., heterogeneous speed of sound in tissue and limited range of steering angle), CPWC suffers from a variety of imaging artifacts resulting from side lobes, grating lobes, and axial lobes. This study focuses on addressing the issues of axial lobes for CPWC, which constitutes an important source of clutter that leads to the degradation of contrast ratio and contrast-to-noise ratio (CR and CNR) of CPWC. We first investigated the source of the axial lobes based on plane-wave propagation and the delay-and-sum (DAS) beamforming. We then proposed a new method that is based on pixel-oriented adaptive apodization (POAA) to eliminate the axial lobes throughout the entire field of view (FOV). POAA was first validated in a simulation study, followed by in vitro phantom experiments and an in vivo case study on a carotid artery from a healthy volunteer. In the simulation study, suppression of axial lobes by 120 dB was observed from wire targets, and an improvement of CNR by up to 60% was found in a cyst-mimicking digital phantom. In the phantom experiment, POAA showed an improvement in CNR by around 20% over conventional methods. The effectiveness of axial lobe suppression was finally demonstrated in vivo, where POAA showed a substantial suppression of clutters throughout the entire FOV.

Coherent Plane Wave Compounding Combined With Tensor Completion Applied for Ultrafast Imaging

To solve the problem of resolution and contrast in plane wave imaging (PWI), coherent plane wave compounding (CPWC) was introduced, in which scanning was performed at different angles, which can achieve the desired image quality by combining the images obtained from PWI at different angles. However, the application of this idea reduces the frame rate in proportion to the number of plane waves (PWs) or angles, so that in this modality, when dealing with some applications such as shear wave imaging (SWI) and strain imaging, there is always a compromise between the frame rate and the image quality. Tensor completion (TC) is a powerful technique to recover missing information of a low-rank tensor from limited observations based on rank minimization. In this article, we present an idea based on TC to make this compromise lighter; in other words, with a smaller number of angles, we can achieve the desired quality of the output image. To evaluate the proposed idea, plane wave imaging challenge in medical ultrasound (PICMUS) datasets was used, which were recorded at 75 different angles. The results of the resolution evaluation showed that using 20% of the coherent PWs and reconstructing other 80% by TC, compared with the situation of using only 20% of the coherent PWs provided a resolution improvement of 14.97% and 17.4% in the simulated and experimental point targets, respectively. Also, the results of the contrast investigation showed that the contrast ratio (CR) improved by 72.6%, 62.9%, and 111.4% in the simulated cyst target data, experimental cyst targets, and in vivo carotid cross section, respectively. The results confirmed that using 20% of the coherent PWs and reconstructing other 80% by TC, the image quality is very close to that obtained by considering all 75 angles, so that the difference in resolution is less than 2% and the difference in contrast to noise ratio (CNR) is less than 5 dB. Therefore, with this idea, it can be said that less compromise is needed; in other words, despite having a higher frame rate, an acceptable quality can be achieved.

In Vivo Evaluation of Plane Wave Imaging for Abdominal Ultrasonography

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

Analysis of the Time and Phase Delay Resolutions in Ultrasound Baseband I/Q Beamformers

OBJECTIVE

The ultrasound baseband in-phase/quadrature beamformer (IQBF) has been widely employed in medical ultrasound imaging to reduce the amount of channel data or to decrease the data rate of the beamforming process. The aim of this study is to assess the effect of the time and phase delay compensation accuracies on the IQBF and thereby to suggest the criteria for selecting the delay resolutions of the IQBF.

METHODS

Mathematical models of the gain loss (GL) and sidelobe level (SL) in closed form are suggested, and the relationships between the parameters (time and phase delay resolutions of the IQBF and the signal bandwidth) and the errors (GL and SL) are investigated. The performance of the IQBF is compared with that of the traditional radio-frequency beamformer (RFBF). Simulation and phantom and in vivo experimental results are shown to corroborate the theoretical analysis.

RESULTS AND CONCLUSION

Theoretical analysis and simulation and experimental results show that a phase delay resolution with a quantization step of 2π/64 is sufficient for phase compensation and that a time delay resolution with a sampling rate of 4f₀ and 2f₀ in the IQBF is sufficient for data with a -6 dB bandwidth of 50% and 25%, respectively, for similar performance as the RFBF with a sampling rate of 16f₀, where f₀ is the center frequency of the ultrasound signal.

SIGNIFICANCE

The suggested criteria have the potential to be used for designing an efficient IQBF satisfying the desired specifications and beamforming accuracy.

Design of a Volumetric Imaging Sequence Using a Vantage-256 Ultrasound Research Platform Multiplexed With a 1024-Element Fully Sampled Matrix Array

Ultrasound imaging using a matrix array allows real-time multi-planar volumetric imaging. To enhance image quality, the matrix array should provide fast volumetric ultrasound imaging with spatially consistent focusing in the lateral and elevational directions. However, because of the significantly increased data size, dealing with massive and continuous data acquisition is a significant challenge. We have designed an imaging acquisition sequence that handles volumetric data efficiently using a single 256-channel Verasonics ultrasound research platform multiplexed with a 1024-element matrix array. The developed sequence has been applied for building an ultrasonic pupilometer. Our results demonstrate the capability of the developed approach for structural visualization of an ex vivo porcine eye and the temporal response of the modeled eye pupil with moving iris at the volume rate of 30 Hz. Our study provides a fundamental ground for researchers to establish their own volumetric ultrasound imaging platform and could stimulate the development of new volumetric ultrasound approaches and applications.

Contrast and Volume Rate Enhancement of 3-D Ultrasound Imaging Using Aperiodic Plane Wave Angles: A Simulation Study

Three-dimensional plane wave imaging (PWI) with a 2-D array has been studied for ultrafast volumetric imaging in medical ultrasound. Compared to 2-D PWI, 3-D PWI requires the transmission of an increased number of plane waves (PWs) to scan a volume of interest and achieve transmit dynamic focusing in both the lateral and elevational directions. To reduce the number of PW angles for a given 2-D angular range by mitigating the grating lobe level, we propose two aperiodic patterns of PW angles: concentric rings with a uniform radial interval and the well-known sunflower pattern. Both patterns are validated to provide uniform angle distributions without regular periodicity, and thereby reduce the grating lobe level compared to a periodic angle distribution with the same number of PW angles. Simulation studies show that the aperiodic patterns enhance the contrast of B-mode images by approximately 3-6 dB over all depths. This enhancement implies that the aperiodic angle sets can increase the volume rate by approximately 2-6 times compared to the periodic angle set at the same contrast and spatial resolution.

Efficient Transmit Delay Calculation in Ultrasound Coherent Plane-Wave Compound Imaging for Curved Array Transducers

The recently introduced plane-wave compounding method based on multiple plane-wave excitation has enabled several new applications due to its high frame rate (>1000 Hz). In this paper, a new efficient transmit delay calculation method in plane-wave compound imaging for a curved array transducer is presented. In the proposed method, the transmit delay is only calculated for a steering angle of 0° and is shifted along the element of the transducer to obtain other transmit delays for different steering angles. To evaluate the performance of the proposed method, the computational complexity was measured for various transmission conditions. For the number of elements and plane-wave excitations of 128 and 65, respectively, the number of operations was substantially decreased in the proposed method compared with the conventional method (256 vs. 8320). The benefits of the proposed method were demonstrated with phantom and in vivo experiments, where coherent plane-wave compounding with 65 excitations provided larger CR and CNR values compared to nine excitations (−22.5 dB and 2.7 vs. −11.3 dB and 1.9, respectively). These results indicate the proposed method can effectively reduce the computational complexity for plane-wave compound imaging in curved array transducers.

Distributing Synthetic Focusing Over Multiple Push-Detect Events Enhances Shear Wave Elasticity Imaging Performance

Plane wave (PW) imaging is a commonly used method for tracking waves during shear wave elasticity imaging (SWEI), but its unfocused transmission beam reduces tracking accuracy and precision. Coherent compounding minimizes this problem, but SWEI's stringent frame rate requirement and the coarse pitch of most clinical transducers limit its effectiveness. Synthetic aperture imaging (SAI) is an alternate ultrasound imaging approach with a much tighter focus than PW imaging, but its lower transmission power has deterred researchers from using SAI in SWEI. Hadamard-encoded multielement SAI can overcome this limitation. However, only a limited number of subapertures (3-5) can be transmitted in a single push-detect event. We have developed methods to distribute more subapertures or more compounding angles over multiple push-detect events. In this paper, we report the results of experiments conducted on phantoms to assess SWEI's performance when using Hadamard-encoded distributed-multielement synthetic aperture (HDMSA) imaging or distributed plane wave compounding (DPWC) to track shear waves. Tracking shear waves with HDMSA improved the elastographic signal-to-noise ratio (SNR_e) by 61.6%-89.5% depending on the phantom employed. Similarly, DPWC tracking improved SNR_e by 56.2%-93.3% for the same group of phantoms. Compared to focused ultrasound tracking (at the focus), SNR_e improved by 28.6% and 32.5% when tracking shear waves with HDMSA and DPWC, respectively. Long acquisitions could introduce decoding errors that decrease the performance when performing HDMSA tracking within the clinical setting. Nevertheless, the results of studies performed on the bicep muscle of three healthy volunteers demonstrate that for stationary organs, tracking shear waves with HDMSA yielded repeatable elastograms that offer better elastographic performance than those produced with current tracking methods.