Sparse Rectangular and Spiral Array Designs for 3D Medical Ultrasound Imaging

High frequency ultrasound 2D array design and fabrication with 3D printed interposers at 200 μm pitch

High-frequency 2D ultrasound arrays (10-50 MHz) can adjust the focus position in three-dimensional space to produce highly focused beams and have potential use in ultrasound elastography and ophthalmology applications. However, these applications remain unfulfilled due to the difficulty in realizing the array element pitch requirements of high-frequency ultrasound 2D arrays with existing technologies. In recent years, the rapid development of 3D printing technology has made complex and fine structures possible. By utilizing high-resolution 3D printing technology, we can effectively solve the connection problem between circuits and acoustic modules in the ultrasound 2D array fabrication process, making it possible to fabricate high-frequency ultrasound 2D arrays. In this work, we demonstrate two alternative fabrication processes for high-frequency 2D ultrasound arrays with a realized pitch of 200 µm and 256 elements per array. The fabricated arrays achieved a yield of 82 % and 75 %, with an average fractional bandwidth of 43 %. The element-to-element variability in bandwidth was 18 %, which reflects the inherent challenges of fine-pitch high-frequency array fabrication. These results highlight the feasibility of scalable fabrication techniques while identifying areas for future optimization.

Ultrafast 3D synthetic aperture imaging with Hadamard-encoded aperiodic interval codes and aperiodic sparse arrays with separate transmitters and receivers

3D synthetic aperture (SA) imaging of volumes can be obtained using sparse 2D ultrasound arrays. However, even with just 256 elements, the volumetric imaging rate can be relatively slow due to having to transmit on each element in succession. Hadamard Aperiodic Interval (HAPI) codes can be used to image the full SA dataset in one extended transmit to speed up the synthetic aperture imaging, but their long nature produces large deadzones if the same elements are used as both transmitters and receivers. In this simulation study, we use a 2D Costas sparse array with separate transmitters and receivers to remedy the deadzone problem, and use it with the HAPI-coded imaging scheme to obtain fully transmit-receive focused, wide field-of-view 3D volumes with high-resolution and high SNR at ultrafast volumetric imaging rates of more than 500 volumes per second, almost nine times faster than non-coded SA imaging with the same imaging parameters. We show similar PSF performance compared to non-coded SA, and a 26 dB improvement in SNR with order-256 HAPI codes. We also present cyst simulations showing similar contrast for the HAPI-coded SA method compared to non-coded SA in the context of no noise, and improved contrast in the context of noise.

Design and Evaluation of a Weighted Periodic Sparse Array for Low-Complexity 1-D Phased Array Ultrasound Imaging Systems

A sparse array offers a significant reduction in the complexity of ultrasonic imaging systems by decreasing the number of active elements and associated electrical circuits needed to form a focused beam. Consequently, for 1-D arrays, it has been adopted in the development of miniaturized systems such as portable, handheld, or smartphone-based systems. Previously, we developed an analytic method that can design a pair of 1-D periodic sparse arrays (PSAs) satisfying three specific constraints, which are the array size, desired grating lobe level, and sparseness factor (SF). In this study, we further developed our method by incorporating aperture weighting functions, which take the form of tapered rectangular functions to introduce null points on the beam pattern. These null points effectively suppress grating lobes generated by a matching pair of arrays. The design process commences with determining transmit and receive PSA patterns, followed by deriving corresponding aperture weighting functions. First, aperture functions of a base and weighting arrays are convolved, which is then upsampled to the targeted array size. Finally, the upsampled aperture is convolved to an aperture function of a subarray, resulting in weighted PSAs (wPSAs). Pulsed wave (PW) simulation confirmed improved grating lobe suppression with wPSAs compared to PSAs. Phantom imaging experiments using a 1-D phased array validated the enhanced contrast due to suppressed grating lobes but at the cost of small degradation in lateral resolution. The signal-to-noise ratio (SNR) also gradually declined with the greater SFs, but no significant difference in SNR was observed between wPSAs and PSAs. Finally, in vivo echocardiography imaging highlighted the clinical potential of wPSAs, particularly with high SFs. Overall, these results suggest that wPSAs can effectively enhance contrast compared to PSAs under the given SF or, alternatively, wPSA with greater SFs can achieve comparable image quality to PSAs with lower SFs. In conclusion, the wPSA approach holds promise for further reducing the complexity of ultrasound imaging systems.

An Adaptive Harmonic Separation Technique for Ultrasound Harmonic Imaging

OBJECTIVE

An adaptive harmonic separation (HS) technique is proposed to overcome the limitations in conventional filtering techniques for ultrasound (US) tissue harmonic imaging (THI).

METHODS

Based on expectation-maximization source separation, the proposed HS technique adaptively models the depth-varying fundamental and harmonic components in the frequency domain and separates the two by applying their calculated posterior probabilities. Phantom experiments with a Tx center frequency of 2 MHz are conducted to evaluate the proposed HS-based US THI schemes.

RESULTS

The phantom images show that the proposed single-pulse THI scheme utilizing the HS technique provides not only an average improvement of 19.2% in axial resolution compared to the conventional bandpass filtering scheme but also similar image quality to that of the conventional pulse-inversion (PI) scheme which requires two Tx/Rx sequences for each scan line. Furthermore, when combined with the PI technique, the HS technique provides a uniform axial resolution over the entire 170 mm imaging depth with an average improvement of 17.1% compared to the conventional PI scheme.

CONCLUSION

These results show that the proposed adaptive HS technique is capable of improving both the frame rate and the image quality of US THI.

Design and Implementation of Analog-Digital Hybrid Beamformers for Low-Complexity Ultrasound Systems: A Feasibility Study

Low-complexity ultrasound systems are increasingly desired for both wearable, point-of-care ultrasound and high-end massive-channel ultrasound for 3-D matrix imaging. However, the imaging capabilities, including spatial resolution and contrast, could suffer as low complexity systems are pursued, which remains as an unresolved tradeoff. To mitigate this limitation, this study revisits the general structures of analog and digital beamformers and introduces a hybrid approach, referred to as analog-digital hybrid beamforming, to implement efficient ultrasound systems. The suggested hybrid beamforming takes two stages sequentially, where the first analog stage partially beamforms M-channel RF signals to N sum-out data (i.e., M-to-N beamforming), and the second digital stage beamforms N partial sums to single final beamformed data (i.e., N-to-1 beamforming). Our approach was systematically designed and implemented with only four major integrated circuits, which was capable of driving full 64-channel transmission and reception. The developed system was demonstrated with a customized 64-channel 1-D phased array using a commercial tissue mimicking phantom. From the phantom imaging results, signal-to-noise ratio, contrast-to-noise ratio, and full beam width at half maximum values were quantitatively evaluated. The demonstrated results indicate that the analog-digital hybrid beamforming can be applied to any type of array for sophisticated 3-D imaging and tiny wearable ultrasound applications.

2-D Array Design and Fabrication With Pitch-Shifting Interposer at Frequencies From 4 MHz up to 10 MHz

High element density and strict constraints of the element's size have significantly limited the design and fabrication of 2-D ultrasonic arrays, especially fully sampled 2-D arrays. Recently, 3-D printing technology has been one of the most rapidly developing fields. Along with the great progress of 3-D printing technology, complex and detailed 3-D structures have become readily available with a short iteration cycle, which allows us to reduce the complexity of routing and helps to ameliorate assembly problems in 2-D ultrasound array fabrication. In this work, we designed and fabricated 2-D ultrasound arrays for an array of applications with a pitch-shifting interposer, which allowed us to fit different array designs with the same circuit design and significantly reduce the requirements in routing and connection for 2-D array fabrication at frequencies from 4 to 10 MHz. Results demonstrated that this design would make 2-D arrays more available and affordable.

Design, Implementation, and Medical Applications of 2-D Ultrasound Sparse Arrays

An ultrasound sparse array consists of a sparse distribution of elements over a 2-D aperture. Such an array is typically characterized by a limited number of elements, which in most cases is compatible with the channel number of the available scanners. Sparse arrays represent an attractive alternative to full 2-D arrays that may require the control of thousands of elements through expensive application-specific integrated circuits (ASICs). However, their massive use is hindered by two main drawbacks: the possible beam profile deterioration, which may worsen the image contrast, and the limited signal-to-noise ratio (SNR), which may result too low for some applications. This article reviews the work done for three decades on 2-D ultrasound sparse arrays for medical applications. First, random, optimized, and deterministic design methods are reviewed together with their main influencing factors. Then, experimental 2-D sparse array implementations based on piezoelectric and capacitive micromachined ultrasonic transducer (CMUT) technologies are presented. Sample applications to 3-D (Doppler) imaging, super-resolution imaging, photo-acoustic imaging, and therapy are reported. The final sections discuss the main shortcomings associated with the use of sparse arrays, the related countermeasures, and the next steps envisaged in the development of innovative arrays.

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.

Design of 2D Planar Sparse Binned Arrays Based on the Coarray Analysis

The analysis of the beampattern is the base of sparse arrays design process. However, in the case of bidimensional arrays, this analysis has a high computational cost, turning the design process into a long and complex task. If the imaging system development is considered a holistic process, the aperture is a sampling grid that must be considered in the spatial domain through the coarray structure. Here, we propose to guide the aperture design process using statistical parameters of the distribution of the weights in the coarray. We have studied three designs of sparse matrix binned arrays with different sparseness degrees. Our results prove that there is a relationship between these parameters and the beampattern, which is valuable and improves the array design process. The proposed methodology reduces the computational cost up to 58 times with respect to the conventional fitness function based on the beampattern analysis.

Design of Ultrasonic Synthetic Aperture Imaging Systems Based on a Non-Grid 2D Sparse Array

This work provides a guide to design ultrasonic synthetic aperture systems for non-grid two-dimensional sparse arrays such as spirals or annular segmented arrays. It presents an algorithm that identifies which elements have a more significant impact on the beampattern characteristics and uses this information to reduce the number of signals, the number of emitters and the number of parallel receiver channels involved in the beamforming process. Consequently, we can optimise the 3D synthetic aperture ultrasonic imaging system for a specific sparse array, reducing the computational cost, the hardware requirements and the system complexity. Simulations using a Fermat spiral array and experimental data based on an annular segmented array with 64 elements are used to assess this algorithm.

A two-steps implementation of 3D ultrasound imaging in frequency domain with 1D array transducer

Compared with B-mode imaging, three-dimensional (3D) ultrasound imaging is more helpful in research and clinical application. At present, the 3D ultrasound images can be acquired directly with two-dimensional (2D) array transducer or reconstructed from a series of B-mode images obtained with one-dimensional (1D) array transducer. Imaging with 2D array transducer can achieve a high frame rate, but suffering from the complexity of the imaging system, such as the large amount of channels, and high computational complexity. Reconstructing 3D images from a series of B-mode images can be implemented by recording the position and orientation of the slice images. This is a low-cost and flexible imaging method, but usually suffering from the low imaging quality and low frame rate. In our previous work, a novel 3D ultrasound imaging method in frequency domain with a moved 1D array transducer is presented. This method can reduce the computational complexity with FFT, and get improved imaging quality and frame rate to some extent. Besides, this method can be adopted to construct images with a row-column-addressed 2D array, which can reduce the amount of channels effectively. In this paper, a two-steps implementation of this imaging method is proposed, in which the combined implementation of the 3D imaging is decomposed to two steps of 2D imaging processes in Frequency domain. In the first step, the received echoes of the 1D array transducer at each position are processed with a 2D imaging processes in the lateral- axial planes. In the second step, a 2D imaging processes is preformed in the planes of orthogonal to the transducer. Simulation results show that the two-steps implementation can achieve almost the same imaging quality to the previous work. Compared with the implementation of 3D imaging in our previous work, the proposed two-steps implementation can be carried out with parallel process to improve the computational efficiency, or carried out with loop to reduce the hardware cost. Besides, the first step can be performed with a conventional DAS imaging method when a cylindrical wave is adopted for imaging. The influence of the spread angle of the field is also discussed.

Tensor Velocity Imaging With Motion Correction

This article presents a motion compensation procedure that significantly improves the accuracy of synthetic aperture tensor velocity estimates for row-column arrays. The proposed motion compensation scheme reduces motion effects by moving the image coordinates with the velocity field during summation of low-resolution volumes. The velocity field is estimated using a transverse oscillation cross-correlation estimator, and each image coordinate's local tensor velocity is determined by upsampling the field using spline interpolation. The motion compensation procedure is validated using Field II simulations and flow measurements acquired using a 3-MHz row-column addressed probe and the research scanner SARUS. For a peak velocity of 25 cm/s, a pulse repetition frequency of 2 kHz, and a beam-to-flow angle of 60°, the proposed motion compensation procedure was able to reduce the relative bias from -27.0% to -9.4% and the standard deviation from 8.6% to 8.1%. In simulations performed with a pulse repetition frequency of 10 kHz, the proposed method reduces the bias in all cases with beam-to-flow angles of 60° and 75° and peak velocities between 10 and 150 cm/s.

Current Development and Applications of Super-Resolution Ultrasound Imaging

Abnormal changes of the microvasculature are reported to be key evidence of the development of several critical diseases, including cancer, progressive kidney disease, and atherosclerotic plaque. Super-resolution ultrasound imaging is an emerging technology that can identify the microvasculature noninvasively, with unprecedented spatial resolution beyond the acoustic diffraction limit. Therefore, it is a promising approach for diagnosing and monitoring the development of diseases. In this review, we introduce current super-resolution ultrasound imaging approaches and their preclinical applications on different animals and disease models. Future directions and challenges to overcome for clinical translations are also discussed.

Design of 2D Sparse Array Transducers for Anomaly Detection in Medical Phantoms

Aperiodic sparse 2D ultrasonic array configurations, including random array, log spiral array, and sunflower array, have been considered for their potential as conformable transducers able to image within a focal range of 30-80 mm, at an operating frequency of 2 MHz. Optimisation of the imaging performance of potential array patterns has been undertaken based on their simulated far field directivity functions. Two evaluation criteria, peak sidelobe level (PSL) and integrated sidelobe ratio (ISLR), are used to access the performance of each array configuration. Subsequently, a log spiral array pattern with -19.33 dB PSL and 2.71 dB ISLR has been selected as the overall optimal design. Two prototype transducers with the selected log spiral array pattern have been fabricated and characterised, one using a fibre composite element composite array transducer (CECAT) structure, the other using a conventional 1-3 composite (C1-3) structure. The CECAT device demonstrates improved coupling coefficient (0.64 to 0.59), reduced mechanical cross-talk between neighbouring array elements (by 10 dB) and improved operational bandwidth (by 16.5%), while the C1-3 device performs better in terms of sensitivity (~50%). Image processing algorithms, such as Hough transform and morphological opening, have been implemented to automatically detect and dimension particles located within a fluid-filled tube structure, in a variety of experimental scenarios, including bespoke phantoms using tissue mimicking material. Experiments using the fabricated CECAT log spiral 2D array transducer demonstrated that this algorithmic approach was able to detect the walls of the tube structure and stationary anomalies within the tube with a precision of ~0.1 mm.