Experimental 3-D Ultrasound Imaging with 2-D Sparse Arrays using Focused and Diverging Waves

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.

Assessing Therapeutic Nanoparticle Accumulation in Tumors Using Nanobubble-Based Contrast-Enhanced Ultrasound Imaging

This study explores the challenges associated with nanoparticle-based drug delivery to the tumor parenchyma, focusing on the widely utilized enhanced permeability and retention effect (EPR). While EPR has been a key strategy, its inconsistent clinical success lacks clear mechanistic understanding and is hindered by limited tools for studying relevant phenomena. This work introduces an approach that employs multiparametric dynamic contrast-enhanced ultrasound (CEUS) with a nanoscale contrast agent for noninvasive, real-time examination of tumor microenvironment characteristics. We demonstrate that CEUS imaging can: (1) evaluate tumor microenvironment features, (2) be used to help predict the distribution of doxorubicin-loaded liposomes in the tumor parenchyma, and (3) be used to predict nanotherapeutic efficacy. CEUS using nanobubbles (NBs) was carried out in two tumor types of high (LS174T) and low (U87) vascular permeability. LS174T tumors consistently showed significantly different time intensity curve (TIC) parameters, including area under the rising curve (AUCR, 2.7×) and time to peak intensity (TTP, 1.9×) compared to U87 tumors. Crucially, a recently developed decorrelation time (DT) parameter specific to NB CEUS dynamics successfully predicted the distribution of doxorubicin-loaded liposomes within the tumor parenchyma (r = 0.86 ± 0.13). AUCR, TTP, and DT were used to correlate imaging findings to nanotherapeutic response with 100% accuracy in SKOV-3 tumors. These findings suggest that NB-CEUS parameters can effectively discern tumor vascular permeability, serving as a biomarker for identifying tumor characteristics and predicting the responsiveness to nanoparticle-based therapies. The observed differences between LS174T and U87 tumors and the accurate prediction of nanotherapeutic efficacy in SKOV-3 tumors indicate the potential utility of this method in predicting treatment efficacy and evaluating EPR in diseases characterized by pathologically permeable vasculature. Ultimately, this research contributes valuable insights into refining drug delivery strategies and assessing the broader applicability of EPR-based approaches.

TOBE-Costas Arrays for Fast High-Resolution 3-D Power Doppler Imaging

Two-dimensional sparse arrays and row-column arrays are both alternatives to 2-D fully addressed arrays with lower channel counts. Row-column arrays have recently demonstrated fast 3-D structural and flow imaging but commonly suffer from high grating lobes or require multiplexing to achieve better quality. Two-dimensional sparse arrays enable full-volume acquisitions for each transmit event, but plane-wave transmissions with them usually lack quality in terms of uniformity of wavefronts. Here, we propose a novel architecture that combines both types of these arrays in one aperture, enabling imaging using row-column or sparse arrays alone or a hybrid imaging scheme where the row-column array is used in transmission and a 2-D sparse array in reception. This hybrid imaging scheme can potentially solve the shortcomings of each of these approaches. The sparse array layout chosen is a Costas array, characterized by having only one element per row and column, facilitating its integration with row-column arrays. We simulate images acquired with TOBE-Costas arrays using the hybrid imaging scheme and compare them to row-column and sparse spiral arrays of equivalent aperture size (128λ × 128λ at 7.5 MHz) in ultrafast plane-wave imaging of point targets and 3-D power Doppler imaging of synthetic flow phantoms. Our simulation results show that TOBE-Costas arrays exhibit superior resolution and lower sidelobe levels compared with plane-wave compounding with row-column arrays. Compared with density-tapered spiral arrays, they provide a larger field of view and finer resolution.

Row Transmission for High Volume-Rate Ultrasound Imaging With a Matrix Array

The widely used Vermon 1024-element matrix array for 3-D ultrasound imaging has three blank rows in the elevational direction, which breaks the elevation periodicity, thus degrading volumetric image quality. To bypass the blank rows in elevation while maintaining the steering capability in azimuth, we proposed a row-transmission (RT) scheme to improve 3-D spatial resolution. Specifically, we divided the full array into four apertures, each with multiple rows along the elevation. Each multirow aperture (MRA) was further divided into subapertures to transmit diverging waves (DWs) sequentially. Coherent DW compounding (CDWC) was realized in azimuth, while the elevation was multielement synthetic aperture (M-SA) imaging by regarding each row as an array of dashed line elements. An in-house spatiotemporal coding strategy, cascaded synthetic aperture (CaSA), was incorporated into the RT scheme as RT-CaSA to increase the signal-to-noise ratio (SNR). We compared the proposed RT with conventional bank-by-bank transmission-reception (Bank) and sparse-random-aperture compounding (SRAC) in a wire phantom and the in vivo human abdominal aorta (AA) to assess the performance of anatomical imaging and aortic wall motion estimation. Phantom results demonstrated superior lateral resolution achieved by our RT scheme (+19.52% and +16.88% versus Bank, +15.32% and +19.72% versus SRAC, in the azimuth-depth and elevation-depth planes, respectively). Our RT-CaSA showed excellent contrast ratios (CRs) (+8.19 and +8.08 dB versus Bank, +6.81 and +5.85 dB versus SRAC, +0.99 and +0.90 dB versus RT) and the highest in vivo aortic wall motion estimation accuracy. The RT scheme was demonstrated to have potential for various matrix array-based 3-D imaging research.

Branched Convolutional Neural Networks for Receiver Channel Recovery in High-Frame-Rate Sparse-Array Ultrasound Imaging

For high-frame-rate ultrasound imaging, it remains challenging to implement on compact systems as a sparse imaging configuration with limited array channels. One key issue is that the resulting image quality is known to be mediocre not only because unfocused plane-wave excitations are used but also because grating lobes would emerge in sparse-array configurations. In this article, we present the design and use of a new channel recovery framework to infer full-array plane-wave channel datasets for periodically sparse arrays that operate with as few as one-quarter of the full-array aperture. This framework is based on a branched encoder-decoder convolutional neural network (CNN) architecture, which was trained using full-array plane-wave channel data collected from human carotid arteries (59 864 training acquisitions; 5-MHz imaging frequency; 20-MHz sampling rate; plane-wave steering angles between -15° and 15° in 1° increments). Three branched encoder-decoder CNNs were separately trained to recover missing channels after differing degrees of channelwise downsampling (2, 3, and 4 times). The framework's performance was tested on full-array and downsampled plane-wave channel data acquired from an in vitro point target, human carotid arteries, and human brachioradialis muscle. Results show that when inferred full-array plane-wave channel data were used for beamforming, spatial aliasing artifacts in the B-mode images were suppressed for all degrees of channel downsampling. In addition, the image contrast was enhanced compared with B-mode images obtained from beamforming with downsampled channel data. When the recovery framework was implemented on an RTX-2080 GPU, the three investigated degrees of downsampling all achieved the same inference time of 4 ms. Overall, the proposed framework shows promise in enhancing the quality of high-frame-rate ultrasound images generated using a sparse-array imaging setup.

Increasing abdominal aortic aneurysm curvature visibility using 3D dual probe bistatic ultrasound imaging combined with probe translation

High frame rate ultrasound (US) imaging techniques in 3D are promising tools for capturing abdominal aortic aneurysms (AAAs) over time, however, with the limited number of channel-to-element connections current footprints are small, which limits the field of view. Moreover, the maximal steering angle of the ultrasound beams in transmit and the maximal receptance angle in receive are insufficient for capturing the curvy shape of the AAA. Therefore, an approach is needed towards large arrays. In this study, high frame rate bistatic 3D US data (17 Hz) were acquired with two synchronized matrix arrays positioned at different locations (multi-aperture imaging) using a translation stage to simulate what a larger array with limited channel-to-element connections can potentially achieve. Acquisitions were performed along an AAA shaped phantom with different probe tilting angles (0 up to ± 30°). The performance of different multi-aperture configurations was quantified using the generalized contrast-to-noise ratio of the wall and lumen (gCNR). Furthermore, a parametric model of the multi-aperture system was used to estimate in which AAA wall regions the contrast is expected to be high. This was evaluated for AAAs with increasing diameters and curvature. With an eight-aperture 0° probe angle configuration a 69 % increase in field of view was measured in the longitudinal direction compared to the field of view of a single aperture configuration. When increasing the number of apertures from two to eight, the gCNR improved for the upper wall and lower wall by 35 % and 13 % (monostatic) and by 36 % and 13 % (bistatic). Contrast improvements up to 22 % (upper wall) and 12 % (lower wall) are achieved with tilted probe configurations compared to non-tilted configurations. Moreover, with bistatic imaging with tilted probe configurations gCNR improvements up to 4 % (upper wall) and 7 % (lower wall) are achieved compared to monostatic imaging. Furthermore, imaging with a larger inter-probe distance improved the gCNR for a ± 15° probe angle configuration. The gCNR has an expected pattern over time, where the contrast is lower when there is more wall motion (systole) and higher when motion is reduced (diastole). Furthermore, a higher frame rate (45 Hz) yields a lower gCNR, because fewer compound angles are used. The results of the parametric model suggest that a flat array is suitable for imaging AAA shapes with limited curvature, but that it is not suitable for imaging larger AAA shapes with more curvature. According to the model, tilted multi-aperture configurations combined with bistatic imaging can achieve a larger region with high contrast compared to non-tilted configurations. The findings of the model are in agreement with experimental findings. To conclude, this study demonstrates the vast improvements in field of view and AAA wall visibility that a large, sparsely populated 3D array can potentially achieve when imaging AAAs compared to single or dual aperture imaging. In the future, larger arrays, less thermal noise, more steering, and more channel-to-element connections combined with carefully chosen orientations of (sub-) apertures will likely advance 3D imaging of AAAs.

Split-based elevational localization of photoacoustic guidewire tip by 1D array probe using spatial impulse response

Objective. Photoacoustic emitters on the tip of a therapeutic device have been intensively studied for echo-guided intervention purposes. In this study, a novel method for localizing the guidewire tip emitter in the elevation direction using a 1D array probe is proposed to resolve the issue of the tip potentially deviating from the ultrasound-imaged plane.Approach. Our method uses the 'interference split' that appears when the emitter is off-plane. Here, a point source from the emitter splits into two points in images. Based on the split, 'split-based elevation localization (SEL)' is introduced to estimate the absolute elevation position of the emitter. Additionally, 'Signed SEL' incorporates an asymmetric feature into the 1D probe to obtain the sign of the elevation localization. An attenuative coupler is attached to the half side of the probe to control the interference split. In SEL and Signed SEL, we propose a modeled split matching (MSM) algorithm to localize the tip position. MSM performs pattern matching of a measured split waveform with modeled split waveforms corresponding to all emitter positions in a region of interest. The modeled waveforms are precalculated using the spatial impulse response. The proposed method is numerically and experimentally validated.Main results. Numerical simulations for time-domain wave propagation clearly demonstrated the interference split phenomena. In the experimental validation with a vessel-mimicking phantom, the proposed methods successfully estimated the elevation positions,yb.SEL exhibited a root-mean-squared error (RMSE) of 2.0 mm for the range of 0 mm ≤yb≤ 30 mm, while Signed SEL estimated the absolute position with an RMSE of 2.4 mm and the sign with an accuracy of 80.8% for the range of -30 mm ≤yb≤ 30 mm.Significance.These results suggest that the proposed method could provide approximate tip positions and help sonographers track it by fanning the probe.

A Multiplexed 32 × 32 2D Matrix Array Transducer for Flexible Sub-Aperture Volumetric Ultrasound Imaging

A fully-sampled two-dimensional (2D) matrix array ultrasonic transducer is essential for fast and accurate three-dimensional (3D) volumetric ultrasound imaging. However, these arrays, usually consisting of thousands of elements, not only face challenges of poor performance and complex wiring due to high-density elements and small element sizes but also put high requirements for electronic systems. Current commercially available fully-sampled matrix arrays, dividing the aperture into four fixed sub-apertures to reduce system channels through multiplexing are widely used. However, the fixed sub-aperture configuration limits imaging flexibility and the gaps between sub-apertures lead to reduced imaging quality. In this study, we propose a high-performance multiplexed matrix array by the design of 1-3 piezocomposite and gapless sub-aperture configuration, as well as optimized matching layer materials. Furthermore, we introduce a sub-aperture volumetric imaging method based on the designed matrix array, enabling high-quality and flexible 3D ultrasound imaging with a low-cost 256-channel system. The influence of imaging parameters, including the number of sub-apertures and steering angle on imaging quality was investigated by simulation, in vitro and in vivo imaging experiments. The fabricated matrix array has a center frequency of 3.4 MHz and a -6 dB bandwidth of above 70%. The proposed sub-aperture volumetric imaging method demonstrated a 10% improvement in spatial resolution, a 19% increase in signal-to-noise ratio, and a 57.7% increase in contrast-to-noise ratio compared with the fixed sub-aperture array imaging method. This study provides a new strategy for high-quality volumetric ultrasound imaging with a low-cost system.

Sub-aperture ultrafast volumetric ultrasound imaging for fully sampled dual-mode matrix array

Fully sampled dual-mode matrix array ultrasound transducer is capable of performing imaging and therapeutic ultrasound in three dimensions (3D). It is a promising tool for many clinical applications because of its precise multi-focus therapy with imaging guidance by itself. Our team previously designed a 256-element fully sampled dual-mode matrix array transducer, while its imaging quality needs to be further improved. In this work, we propose a high-contrast sub-aperture volumetric imaging strategy to improve the imaging quality of the dual-mode matrix array. We first analyzed the effect of various parameters of sub-aperture imaging on the imaging quality by Field II. Based on the optimized parameters, we compared the resolution and signal to noise ratio (SNR) of sub-aperture imaging with those of full aperture imaging on phantoms and rabbit brain. The experimental results showed the proposed sub-aperture imaging method could obtain a comparable resolution to full aperture imaging. Moreover, the average intensity of noise signal near the wire phantom decreased by about 5 dB and the SNR of tissue phantom image increased by 8 %. The proposed sub-aperture imaging method also enabled clearer and more accurate imaging of the rabbit brain. The obtained results indicate the proposed sub-aperture imaging is a promising method for practical use of a fully sampled dual-mode matrix array for volumetric ultrasound imaging.

High-Frame-Rate Volumetric Porcine Renal Vasculature Imaging

OBJECTIVE

The aim of this study was to assess the feasibility and imaging options of contrast-enhanced volumetric ultrasound kidney vasculature imaging in a porcine model using a prototype sparse spiral array.

METHODS

Transcutaneous freehand in vivo imaging of two healthy porcine kidneys was performed according to three protocols with different microbubble concentrations and transmission sequences. Combining high-frame-rate transmission sequences with our previously described spatial coherence beamformer, we determined the ability to produce detailed volumetric images of the vasculature. We also determined power, color and spectral Doppler, as well as super-resolved microvasculature in a volume. The results were compared against a clinical 2-D ultrasound machine.

RESULTS

Three-dimensional visualization of the kidney vasculature structure and blood flow was possible with our method. Good structural agreement was found between the visualized vasculature structure and the 2-D reference. Microvasculature patterns in the kidney cortex were visible with super-resolution processing. Blood flow velocity estimations were within a physiological range and pattern, also in agreement with the 2-D reference results.

CONCLUSION

Volumetric imaging of the kidney vasculature was possible using a prototype sparse spiral array. Reliable structural and temporal information could be extracted from these imaging results.

Assessing Tumor Microenvironment Characteristics and Stratifying EPR with a Nanobubble Companion Nanoparticle via Contrast-Enhanced Ultrasound Imaging

The tumor microenvironment is characterized by dysfunctional endothelial cells, resulting in heightened vascular permeability. Many nanoparticle-based drug delivery systems attempt to use this enhanced permeability combined with impaired lymphatic drainage (a concept known as the 'enhanced permeability and retention effect' or EPR effect) as the primary strategy for drug delivery, but this has not proven to be as clinically effective as anticipated. The specific mechanisms behind the inconsistent clinical outcomes of nanotherapeutics have not been clearly articulated, and the field has been hampered by a lack of accessible tools to study EPR-associated phenomena in clinically relevant scenarios. While medical imaging has tremendous potential to contribute to this area, it has not been broadly explored. This work examines, for the first time, the use of multiparametric dynamic contrast-enhanced ultrasound (CEUS) with a novel nanoscale contrast agent to examine tumor microenvironment characteristics noninvasively and in real-time. We demonstrate that CEUS imaging can: (1) evaluate tumor microenvironment features and (2) be used to help predict the distribution of doxorubicin-loaded liposomes in the tumor parenchyma. CEUS using nanobubbles (NBs) was carried out in two tumor types of high (LS174T) and low (U87) vascular permeability, and time-intensity curve (TIC) parameters were evaluated in both models prior to injection of doxorubicin liposomes. Consistently, LS174T tumors showed significantly different TIC parameters, including area under the rising curve (2.7x), time to peak intensity (1.9x) and decorrelation time (DT, 1.9x) compared to U87 tumors. Importantly, the DT parameter successfully predicted tumoral nanoparticle distribution (r = 0.86 ± 0.13). Ultimately, substantial differences in NB-CEUS generated parameters between LS174T and U87 tumors suggest that this method may be useful in determining tumor vascular permeability and could be used as a biomarker for identifying tumor characteristics and predicting sensitivity to nanoparticle-based therapies. These findings could ultimately be applied to predicting treatment efficacy and to evaluating EPR in other diseases with pathologically permeable vasculature.

Optimized Virtual Sources Distributions for 3-D Ultrafast Diverging Wave Compounding Imaging: A Simulation Study

Ultrafast ultrasound imaging allows observing rapid phenomena; combined with 3-D imaging it has the potential to provide a more accurate analysis of organs which leads, in the end, to better diagnosis. Coherent compounding using diverging waves is commonly used to reconstruct high-quality images on large volumes while keeping the frame rate high enough to allow dynamic analysis. In practice, the virtual sources (VSs) that drive the diverging waves are often distributed in a deterministic way: following a regular grid, concentric rings, and spirals. Even though those deterministic distributions can offer various tradeoffs in terms of imaging performance, other distributions can be considered to improve imaging performance. It is herein suggested to look at alternative VSs distributions for optimizing the lateral resolution and the secondary lobes level (SLL) on several point spread functions (PSFs) by means of a multiobjective genetic algorithm. The optimization framework has led to seven pseudo-irregular distributions of VSs distributions that have not yet been found in the literature. An analysis of the imaging performance with a simulated phantom shows that these new distributions offer different tradeoffs between lateral resolution and contrast, respectively, measured on point-like reflectors and anechoic cysts. As an example, one of these optimized distributions improves the lateral resolution by 16% and gives equivalent contrast values on cysts and PSF isotropy properties, when compared to a concentric-rings-based distribution.

High-Volume-Rate 3-D Ultrasound Imaging Using Fast-Tilting and Redirecting Reflectors

Three-dimensional ultrasound imaging has many advantages over 2-D imaging such as more comprehensive tissue evaluation and less operator dependence. However, developing a low-cost and accessible 3-D ultrasound solution with high volume rate and imaging quality remains a challenging task. Recently, we proposed a 3-D ultrasound imaging technique: fast acoustic steering via tilting electromechanical reflectors (FASTER), which uses a fast-tilting acoustic reflector to steer ultrafast plane waves elevationally to achieve high-volume-rate 3-D imaging with conventional 1-D transducers. However, the initial FASTER implementation requires a water tank for acoustic wave conduction and cannot be conveniently used for regular handheld scanning. To address these limitations, here, we developed a novel ultrasound probe clip-on device that encloses a fast-tilting reflector, a redirecting reflector, and an acoustic wave conduction medium. The new FASTER 3-D imaging device can be easily attached to or removed from clinical ultrasound transducers, allowing rapid transformation from 2-D to 3-D imaging. In vitro B-mode studies demonstrated that the proposed method provided comparable imaging quality to conventional, mechanical-translation-based 3-D imaging while offering a much faster volume rate (e.g., 300 versus  ∼  10 Hz). We also demonstrated 3-D power Doppler (PD) and 3-D super-resolution ultrasound localization microscopy (ULM) with the FASTER device. An in vivo imaging study showed that the FASTER device could clearly visualize the 3-D anatomy of the basilic vein. These results suggest that the newly developed redirecting reflector and the clip-on device could overcome key hurdles for future clinical translation of the FASTER 3-D imaging technology.

Costas Sparse 2-D Arrays for High-Resolution Ultrasound Imaging

Two-dimensional arrays enable volumetric ultrasound imaging but have been limited to small aperture size and hence low resolution due to the high cost and complexity of fabrication, addressing, and processing associated with large fully addressed arrays. Here, we propose Costas arrays as a gridded sparse 2-D array architecture for volumetric ultrasound imaging. Costas arrays have exactly one element for every row and column, such that the vector displacement between any pair of elements is unique. These properties ensure aperiodicity, which helps eliminate grating lobes. Compared with previously reported works, we studied the distribution of active elements based on an order-256 Costas layout on a wider aperture ( 96 λ×96 λ at 7.5 MHz center frequency) for high-resolution imaging. Our investigations with focused scanline imaging of point targets and cyst phantoms showed that Costas arrays exhibit lower peak sidelobe levels compared with random sparse arrays of the same size and offer comparable performance in terms of contrast compared with Fermat spiral arrays. In addition, Costas arrays are gridded, which could ease the manufacturing and has one element for each row/column, which enables simple interconnection strategies. Compared with state-of-the-art matrix probes, which are commonly 32×32 , the proposed sparse arrays achieve higher lateral resolution and a wider field of view.

Coherent Bistatic 3-D Ultrasound Imaging Using Two Sparse Matrix Arrays

In the last decade, many advances have been made in high frame rate 3-D ultrasound imaging, including more flexible acquisition systems, transmit (TX) sequences, and transducer arrays. Compounding multiangle transmits of diverging waves has shown to be fast and effective for 2-D matrix arrays, where heterogeneity between transmits is key in optimizing the image quality. However, the anisotropy in contrast and resolution remains a drawback that cannot be overcome with a single transducer. In this study, a bistatic imaging aperture is demonstrated that consists of two synchronized matrix ( 32×32 ) arrays, allowing for fast interleaved transmits with a simultaneous receive (RX). First, for a single array, the aperture efficiency for high volume rate imaging was evaluated between sparse random arrays and fully multiplexed arrays. Second, the performance of the bistatic acquisition scheme was analyzed for various positions on a wire phantom and was showcased in a dynamic setup mimicking the human abdomen and aorta. Sparse array volume images were equal in resolution and lower in contrast compared to fully multiplexed arrays but can efficiently minimize decorrelation during motion for multiaperture imaging. The dual-array imaging aperture improved the spatial resolution in the direction of the second transducer, reducing the average volumetric speckle size with 72% and the axial-lateral eccentricity with 8%. In the aorta phantom, the angular coverage increased by a factor of 3 in the axial-lateral plane, raising the wall-lumen contrast with 16% compared to single-array images, despite accumulation of thermal noise in the lumen.

3-D Coherent Multi-Transducer Ultrasound Imaging with Sparse Spiral Arrays

Coherent multi-transducer ultrasound (CoMTUS) creates an extended effective aperture through the coherent combination of multiple arrays, which results in images with enhanced resolution, extended field-of-view, and higher sensitivity. The subwavelength localization accuracy of the multiple transducers required to coherently beamform the data is achieved by using the echoes backscattered from targeted points. In this study, CoMTUS is implemented and demonstrated for the first time in 3-D imaging using a pair of 256-element 2-D sparse spiral arrays, which keep the channel count low and limit the amount of data to be processed. The imaging performance of the method was investigated using both simulations and phantom tests. The feasibility of free-hand operation is also experimentally demonstrated. Results show that, in comparison to a single dense array system using the same total number of active elements, the proposed CoMTUS system improves spatial resolution (up to 10 times) in the direction where both arrays are aligned, contrast-to-noise-ratio (CNR, up to 46%), and generalized CNR (up to 15%). Overall, CoMTUS shows a narrower main lobe and higher contrast-to-noise ratio, which results in an increased dynamic range and better target detectability.

A Tiled Ultrasound Matrix Transducer for Volumetric Imaging of the Carotid Artery

High frame rate three-dimensional (3D) ultrasound imaging would offer excellent possibilities for the accurate assessment of carotid artery diseases. This calls for a matrix transducer with a large aperture and a vast number of elements. Such a matrix transducer should be interfaced with an application-specific integrated circuit (ASIC) for channel reduction. However, the fabrication of such a transducer integrated with one very large ASIC is very challenging and expensive. In this study, we develop a prototype matrix transducer mounted on top of multiple identical ASICs in a tiled configuration. The matrix was designed to have 7680 piezoelectric elements with a pitch of 300 μm × 150 μm integrated with an array of 8 × 1 tiled ASICs. The performance of the prototype is characterized by a series of measurements. The transducer exhibits a uniform behavior with the majority of the elements working within the -6 dB sensitivity range. In transmit, the individual elements show a center frequency of 7.5 MHz, a -6 dB bandwidth of 45%, and a transmit efficiency of 30 Pa/V at 200 mm. In receive, the dynamic range is 81 dB, and the minimum detectable pressure is 60 Pa per element. To demonstrate the imaging capabilities, we acquired 3D images using a commercial wire phantom.

Three-Dimensional Ultrasound Localization Microscopy with Bipartite Graph-Based Microbubble Pairing and Kalman-Filtering-Based Tracking on a 256-Channel Verasonics Ultrasound System with a 32 × 32 Matrix Array

Three-dimensional (3D) ultrasound localization microscopy (ULM) using a 2-D matrix probe and microbubbles (MBs) has been recently proposed to visualize microvasculature beyond the ultrasound diffraction limit in three spatial dimensions. However, 3D ULM suffers from several limitations: (1) high system complexity due to numerous channel counts, (2) complex MB flow dynamics in 3D, and (3) extremely long acquisition time. To reduce the system complexity while maintaining high image quality, we used a sub-aperture process to reduce received channel counts. To address the second issue, a 3D bipartite graph-based method with Kalman filtering-based tracking was used in this study for MB tracking. An MB separation approach was incorporated to separate high concentration MB data into multiple, sparser MB datasets, allowing better MB localization and tracking for a limited acquisition time. The proposed method was first validated in a flow channel phantom, showing improved spatial resolutions compared with the contrasted enhanced power Doppler image. Then the proposed method was evaluated with an in vivo chicken embryo brain dataset. Results showed that the reconstructed 3D super-resolution image achieved a spatial resolution of around 52 μm (smaller than the wavelength of around 200 μm). Microvessels that cannot be resolved clearly using localization only, can be well identified with the tailored 3D pairing and tracking algorithms. To sum up, the feasibility of the 3D ULM is shown, indicating the great possibility in clinical applications.

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.

Unfocused Field Analysis of a Density-Tapered Spiral Array for High-Volume-Rate 3-D Ultrasound Imaging

Spiral array transducers with a sparse 2-D aperture have demonstrated their potential in realizing 3-D ultrasound imaging with reduced data rates. Nevertheless, their feasibility in high-volume-rate imaging based on unfocused transmissions has yet to be established. From a metrology standpoint, it is essential to characterize the acoustic field of unfocused transmissions from spiral arrays not only to assess their safety but also to identify the root cause of imaging irregularities due to the array's sparse aperture. Here, we present a field profile analysis of unfocused transmissions from a density-tapered spiral array transducer (256 hexagonal elements, 220- [Formula: see text] element diameter, and 1-cm aperture diameter) through both simulations and hydrophone measurements. We investigated plane- and diverging-wave transmissions (five-cycle, 7.5-MHz pulses) from 0° to 10° steering for their beam intensity characteristics and wavefront arrival time profiles. Unfocused firings were also tested for B-mode imaging performance (ten compounded angles, -5° to 5° span). The array was found to produce unfocused transmissions with a peak negative pressure of 93.9 kPa at 2 cm depth. All transmissions steered up to 5° were free of secondary lobes within 12 dB of the main beam peak intensity. All wavefront arrival time profiles were found to closely match the expected profiles with maximum root-mean-squared errors of [Formula: see text] for plane wave (PW) and [Formula: see text] for diverging wave. The B-mode images showed good spatial resolution with a penetration depth of 22 mm in PW imaging. Overall, these results demonstrate that the density-tapered spiral array can facilitate unfocused transmissions below regulatory limits (mechanical index: 0.034; spatial-peak, pulse-average intensity: 0.298 W/cm2) and with suppressed secondary lobes while maintaining smooth wavefronts.

Sparse 2-D PZT-on-PCB Arrays With Density Tapering

Two-dimensional (2-D) arrays offer volumetric imaging capabilities without the need for probe translation or rotation. A sparse array with elements seeded in a tapering spiral pattern enables one-to-one connection to an ultrasound machine, thus allowing flexible transmission and reception strategies. To test the concept of sparse spiral array imaging, we have designed, realized, and characterized two prototype probes designed at 2.5-MHz low-frequency (LF) and 5-MHz high-frequency (HF) center frequencies. Both probes share the same electronic design, based on piezoelectric ceramics and rapid prototyping with printed circuit board substrates to wire the elements to external connectors. Different center frequencies were achieved by adjusting the piezoelectric layer thickness. The LF and HF prototype probes had 88% and 95% of working elements, producing peak pressures of 21 and 96 kPa/V when focused at 5 and 3 cm, respectively. The one-way -3-dB bandwidths were 26% and 32%. These results, together with experimental tests on tissue-mimicking phantoms, show that the probes are viable for volumetric imaging.

A theranostic 3D ultrasound imaging system for high resolution image-guided therapy

Microbubble contrast agents are a diagnostic tool with broad clinical impact and an increasing number of indications. Many therapeutic applications have also been identified. Yet, technologies for ultrasound guidance of microbubble-mediated therapy are limited. In particular, arrays that are capable of implementing and imaging microbubble-based therapy in three dimensions in real-time are lacking. We propose a system to perform and monitor microbubble-based therapy, capable of volumetric imaging over a large field-of-view. To propel the promise of the theranostic treatment strategies forward, we have designed and tested a unique array and system for 3D ultrasound guidance of microbubble-based therapeutic protocols based on the frequency, temporal and spatial requirements. Methods: Four 256-channel plane wave scanners (Verasonics, Inc, WA, USA) were combined to control a 1024-element planar array with 1.3 and 2.5 MHz therapeutic and imaging transmissions, respectively. A transducer aperture of ~40×15 mm was selected and Field II was applied to evaluate the point spread function. In vitro experiments were performed on commercial and custom phantoms to assess the spatial resolution, image contrast and microbubble-enhanced imaging capabilities. Results: We found that a 2D array configuration with 64 elements separated by λ-pitch in azimuth and 16 elements separated by 1.5λ-pitch in elevation ensured the required flexibility. This design, of 41.6 mm × 16 mm, thus provided both an extended field-of-view, up to 11 cm x 6 cm at 10 cm depth and steering of ±18° in azimuth and ±12° in elevation. At a depth of 16 cm, we achieved a volume imaging rate of 60 Hz, with a contrast ratio and resolution, respectively, of 19 dB, 0.8 mm at 3 cm and 20 dB and 2.1 mm at 12.5 cm. Conclusion: A single 2D array for both imaging and therapeutics, integrated with a 1024 channel scanner can guide microbubble-based therapy in volumetric regions of interest.

Exploring 3D elastic-wave scattering at interfaces using high-resolution phased-array system

The elastic-wave scattering at interfaces, such as cracks, is essential for nondestructive inspections, and hence, understanding the phenomenon is crucial. However, the elastic-wave scattering at cracks is very complex in three dimensions since microscopic asperities of crack faces can be multiple scattering sources. We propose a method for exploring 3D elastic-wave scattering based on our previously developed high-resolution 3D phased-array system, the piezoelectric and laser ultrasonic system (PLUS). We describe the principle of PLUS, which combines a piezoelectric transmitter and a 2D mechanical scan of a laser Doppler vibrometer, enabling us to resolve a crack into a collection of scattring sources. Subsequently, we show how the 3D elastic-wave scattering in the vicinity of each response can be extracted. Here, we experimentally applied PLUS to a fatigue-crack specimen. We found that diverse 3D elastic-wave scattering occurred in a manner depending on the responses within the fatigue crack. This is significant because access to such information will be useful for optimizing inspection conditions, designing ultrasonic measurement systems, and characterizing cracks. More importantly, the described methodology is very general and can be applied to not only metals but also other materials such as composites, concrete, and rocks, leading to progress in many fields.

Boosting transducer matrix sensitivity for 3D large field ultrasound localization microscopy using a multi-lens diffracting layer: a simulation study

Mapping blood microflows of the whole brain is crucial for early diagnosis of cerebral diseases. Ultrasound localization microscopy (ULM) was recently applied to map and quantify blood microflows in 2D in the brain of adult patients down to the micron scale. Whole brain 3D clinical ULM remains challenging due to the transcranial energy loss which significantly reduces the imaging sensitivity. Large aperture probes with a large surface can increase both resolution and sensitivity. However, a large active surface implies thousands of acoustic elements, with limited clinical translation. In this study, we investigate via simulations a new high-sensitive 3D imaging approach based on large diverging elements, combined with an adapted beamforming with corrected delay laws, to increase sensitivity. First, pressure fields from single elements with different sizes and shapes were simulated. High directivity was measured for curved element while maintaining high transmit pressure. Matrix arrays of 256 elements with a dimension of 10 × 10 cm with small (λ/2), large (4λ), and curved elements (4λ) were compared through point spread functions analysis. A large synthetic microvessel phantom filled with 100 microbubbles per frame was imaged using the matrix arrays in a transcranial configuration. 93% of the bubbles were detected with the proposed approach demonstrating that the multi-lens diffracting layer has a strong potential to enable 3D ULM over a large field of view through the bones.

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.

Signal-to-noise ratio of diverging waves in multiscattering media: Effects of signal duration and divergence angle

In this paper, SNR maximization in coded diverging waves is studied, and experimental verification of the results is presented. Complementary Golay sequences and binary phase shift keying modulation are used to code the transmitted signal. The SNR in speckle and pin targets is maximized with respect to chip signal length. The maximum SNR is obtained in diverging wave transmission when the chip signal is as short a duration as the array permits. We determined the optimum diverging wave profile to confine the transmitted ultrasound energy in the imaging sector. The optimized profile also contributes to the SNR maximization. The SNR performances of the optimized coded diverging wave and conventional single-focused phased array imaging are compared on a single frame basis. The SNR of the optimized coded diverging wave is higher than that of the conventional single-focused phased array imaging at all depths and regions.

3D Transcranial Ultrasound Localization Microscopy in the Rat Brain with a Multiplexed Matrix Probe

OBJECTIVE

Ultrasound Localization Microscopy (ULM) provides images of the microcirculation in-depth in living tissue. However, its implementation in two-dimension is limited by the elevation projection and tedious plane-by-plane acquisition. Volumetric ULM alleviates these issues and can map the vasculature of entire organs in one acquisition with isotropic resolution. However, its optimal implementation requires many independent acquisition channels, leading to complex custom hardware.

METHODS

In this article, we implemented volumetric ultrasound imaging with a multiplexed 32 x 32 probe driven by a single commercial ultrasound scanner. We propose and compare three different sub-aperture multiplexing combinations for localization microscopy in silico and in vitro with a flow of microbubbles in a canal. Finally, we evaluate the approach for micro-angiography of the rat brain.The "light" combination allows a higher maximal volume rate than the "full" combination while maintaining the field of view and resolution.

RESULTS

In the rat brain, 100,000 volumes were acquired within 7 min with a dedicated ultrasound sequence and revealed vessels down to 31 m in diameter with flows from 4.3 mm/s to 28.4 mm/s.

CONCLUSION

This work demonstrates the ability to perform a complete angiography with unprecedented resolution in the living rats brain with a simple and light setup through the intact skull.

SIGNIFICANCE

We foresee that it might contribute to democratize 3D ULM for both preclinical and clinical studies.

High Frame Rate Volumetric Imaging of Microbubbles Using a Sparse Array and Spatial Coherence Beamforming

Volumetric ultrasound imaging of blood flow with microbubbles enables a more complete visualization of the microvasculature. Sparse arrays are ideal candidates to perform volumetric imaging at reduced manufacturing complexity and cable count. However, due to the small number of transducer elements, sparse arrays often come with high clutter levels, especially when wide beams are transmitted to increase the frame rate. In this study, we demonstrate with a prototype sparse array probe and a diverging wave transmission strategy, that a uniform transmission field can be achieved. With the implementation of a spatial coherence beamformer, the background clutter signal can be effectively suppressed, leading to a signal to background ratio improvement of 25 dB. With this approach, we demonstrate the volumetric visualization of single microbubbles in a tissue-mimicking phantom as well as vasculature mapping in a live chicken embryo chorioallantoic membrane.

Real-Time Volumetric Synthetic Aperture Software Beamforming of Row-Column Probe Data

Two delay-and-sum beamformers for 3-D synthetic aperture imaging with row-column addressed arrays are presented. Both beamformers are software implementations for graphics processing unit (GPU) execution with dynamic apodizations and third-order polynomial subsample interpolation. The first beamformer was written in the MATLAB programming language and the second was written in C/C++ with the compute unified device architecture (CUDA) extensions by NVIDIA. Performance was measured as volume rate and sample throughput on three different GPUs: a 1050 Ti, a 1080 Ti, and a TITAN V. The beamformers were evaluated across 112 combinations of output geometry, depth range, transducer array size, number of virtual sources, floating-point precision, and Nyquist rate or in-phase/quadrature beamforming using analytic signals. Real-time imaging defined as more than 30 volumes per second was attained by the CUDA beamformer on the three GPUs for 13, 27, and 43 setups, respectively. The MATLAB beamformer did not attain real-time imaging for any setup. The median, single-precision sample throughput of the CUDA beamformer was 4.9, 20.8, and 33.5 Gsamples/s on the three GPUs, respectively. The throughput of CUDA beamformer was an order of magnitude higher than that of the MATLAB beamformer.

Ultrasound Planar Array Imaging Metric Analysis

Planar array design makes the tradeoff between the 3-D ultrasound image quality and the system complexity based on the imaging metrics. The -6 dB mainlobe width (MW), mainlobe-to-sidelobe energy ratio (MSR), peak sidelobe level (PSL), and average sidelobe level (ASL) are the common imaging metrics for linear array design. MW is used for lateral resolution evaluation, while MSR, PSL, and ASL are adopted for contrast resolution evaluation. However, simulation results show that these metrics cannot fully evaluate the planar array performance. This article proposes several new imaging metrics for planar array: averaged mainlobe acoustic energy level and mainlobe energy density curve are the lateral resolution metrics, while mainlobe-to-sidelobe energy density ratio is the contrast resolution metric. The new metrics take into account the influence of the mainlobe area on the planar array performance evaluation. PSF analysis and simulated images show that the proposed metrics can evaluate planar array performance more accurately than the existing metrics. Moreover, uniform planar arrays with different scales and random sparse arrays are tested to show how to use the proposed metrics in planar array design.

Sparse Convolutional Beamforming for 3-D Ultrafast Ultrasound Imaging

Real-time 3-D ultrasound (US) provides a complete visualization of inner body organs and blood vasculature, crucial for diagnosis and treatment of diverse diseases. However, 3-D systems require massive hardware due to the huge number of transducer elements and consequent data size. This increases cost significantly and limit both frame rate and image quality, thus preventing the 3-D US from being common practice in clinics worldwide. A recent study presented a technique called sparse convolutional beamforming algorithm (SCOBA), which obtains improved image quality while allowing notable element reduction in the context of 2-D focused imaging. In this article, we build upon previous work and introduce a nonlinear beamformer for 3-D imaging, called COBA-3D, consisting of 2-D spatial convolution of the in-phase and quadrature received signals. The proposed technique considers diverging-wave transmission and achieves improved image resolution and contrast compared with standard delay-and-sum beamforming while enabling a high frame rate. Incorporating 2-D sparse arrays into our method creates SCOBA-3D: a sparse beamformer that offers significant element reduction and, thus, allows performing 3-D imaging with the resources typically available for 2-D setups. To create 2-D thinned arrays, we present a scalable and systematic way to design 2-D fractal sparse arrays. The proposed framework paves the way for affordable ultrafast US devices that perform high-quality 3-D imaging, as demonstrated using phantom and ex-vivo data.

A Study on Capacitive Micromachined Ultrasonic Transducer Periodic Sparse Array

Capacitive micromachined ultrasonic transducer (CMUT) is an ultrasonic transducer based on the microelectromechanical system (MEMS). CMUT elements are easily made into a high-density array, which will increase the hardware complexity. In order to reduce the number of active channels, this paper studies the grating lobes generated by CMUT periodic sparse array (PSA) pairs. Through the design of active element positions in the transmitting and receiving processes, the simulation results of effective aperture and beam patterns show that the common grating lobes (CGLs) generated by the transmit and receive array are eliminated. On the basis of point targets imaging, a CMUT linear array with 256 elements is used to carry out the PSA pairs experiment. Under the same sparse factor (SF), the optimal sparse array configuration can be selected to reduce the imaging artifacts. This conclusion is of great significance for the application of CMUT in three-dimensional ultrasound imaging.

Architecture for an Ultrasound Advanced Open Platform With an Arbitrary Number of Independent Channels

Ultrasound open platforms are programmable and flexible tools for the development and test of novel methods. In most cases, they embed the electronics for the independent control of (maximum) 256 probe elements. However, a higher number of channels is needed for the control of 2-D array probes. This paper presents a system architecture that, through the hardware and software synchronization of multiple ULA-OP 256 scanners, may implement advanced open platforms with an arbitrary number of channels. The proposed solution needs a single personal computer, maintains real-time features, and preserves portability. A prototype demonstrator, composed of two ULA-OP 256 scanners connected to 512 elements of a matrix array, was implemented and tested according to different channel configurations. Experiments performed under MATLAB control confirmed that by doubling the number of elements (from 256 to 512) the signal-to-noise and contrast ratios improve by 9 dB and 3 dB, respectively. Furthermore, as a full 512-channel scanner, the demonstrator can produce real-time B-mode images at 18 Hz, high enough for probe positioning during acquisitions. Also, the demonstrator permitted the implementation of a new high frame rate, bi-plane, triplex modality. All probe elements are excited to simultaneously produce two planar, perpendicular diverging waves. Each scanner independently processes the echoes received by the 256 connected elements to beamform 1300 frames per second. For each insonified plane, good quality morphological (B-mode), qualitative (color flow-), and quantitative (spectral-) Doppler images are finally shown in real-time by a dedicated interface.

Sparse channel sampling for ultrasound localization microscopy (SPARSE-ULM)

Ultrasound localization microscopy (ULM) has recently enabled the mapping of the cerebral vasculaturein vivowith a resolution ten times smaller than the wavelength used, down to ten microns. However, with frame rates up to 20000 frames per second, this method requires large amount of data to be acquired, transmitted, stored, and processed. The transfer rate is, as of today, one of the main limiting factors of this technology. Herein, we introduce a novel reconstruction framework to decrease this quantity of data to be acquired and the complexity of the required hardware by randomly subsampling the channels of a linear probe. Method performance evaluation as well as parameters optimization were conductedin silicousing the SIMUS simulation software in an anatomically realistic phantom and then compared toin vivoacquisitions in a rat brain after craniotomy. Results show that reducing the number of active elements deteriorates the signal-to-noise ratio and could lead to false microbubbles detections but has limited effect on localization accuracy. In simulation, the false positive rate on microbubble detection deteriorates from 3.7% for 128 channels in receive and 7 steered angles to 11% for 16 channels and 7 angles. The average localization accuracy ranges from 10.6μm and 9.93μm for 16 channels/3 angles and 128 channels/13 angles respectively. These results suggest that a compromise can be found between the number of channels and the quality of the reconstructed vascular network and demonstrate feasibility of performing ULM with a reduced number of channels in receive, paving the way for low-cost devices enabling high-resolution vascular mapping.

Fast 3-D Velocity Estimation in 4-D Using a 62 + 62 Row-Column Addressed Array

This article presents an imaging scheme capable of estimating the full 3-D velocity vector field in a volume using row-column addressed arrays (RCAs) at a high volume rate. A 62 + 62 RCA array is employed with an interleaved synthetic aperture sequence. It contains repeated emissions with rows and columns interleaved with B-mode emissions. The sequence contains 80 emissions in total and can provide continuous volumetric data at a volume rate above 125 Hz. A transverse oscillation cross correlation estimator determines all three velocity components. The approach is investigated using Field II simulations and measurements using a specially built 3-MHz 62 + 62 RCA array connected to the SARUS experimental scanner. Both the B-mode and flow sequences have a penetration depth of 14 cm when measured on a tissue-mimicking phantom (0.5-dB/[ [Formula: see text]] attenuation). Simulations of a parabolic flow in a 12-mm-diameter vessel at a depth of 30 mm, beam-to-flow angle of 90°, and xy-rotation of 45° gave a standard deviation (SD) of (3.3, 3.4, 0.4)% and bias of (-3.3, -3.9, -0.1)%, for ( v_x , v_y , and v_z ). Decreasing the beam-to-flow angle to 60° gave an SD of (8.9, 9.1, 0.8)% and bias of (-7.6, -9.5, -7.2)%, showing a slight increase. Measurements were carried out using a similar setup, and pulsing at 2 kHz yielded comparable results at 90° with an SD of (5.8, 5.5, 1.1)% and bias of (1.4, -6.4, 2.4)%. At 60°, the SD was (5.2, 4.7 1.2)% and bias (-4.6, 6.9, -7.4)%. Results from measurements across all tested settings showed a maximum SD of 6.8% and a maximum bias of 15.8% for a peak velocity of 10 cm/s. A tissue-mimicking phantom with a straight vessel was used to introduce clutter, tissue motion, and pulsating flow. The pulsating velocity magnitude was estimated across ten pulse periods and yielded an SD of 10.9%. The method was capable of estimating transverse flow components precisely but underestimated the flow with small beam-to-flow angles. The sequence provided continuous data in both time and space throughout the volume, allowing for retrospective analysis of the flow. Moreover, B-mode planes can be selected retrospectively anywhere in the volume. This shows that tensor velocity imaging (full 3-D volumetric vector flow imaging) can be estimated in 4-D ( x, y, z, and t ) using only 62 channels in receive, making 4-D volumetric imaging implementable on current scanner hardware.

Helical toroid phantom for 3D flow imaging investigations

The medical physics community has hitherto lacked an effective calibration phantom to holistically evaluate the performance of three-dimensional (3D) flow imaging techniques. Here, we present the design of a new omnidirectional, three-component (3-C) flow phantom whose lumen is consisted of a helical toroid structure (4 mm lumen diameter; helically winded for 5 revolutions over a torus with 10 mm radius; 5 mm helix radius). This phantom's intraluminal flow trajectory embraces all combinations of x, y, and z directional components, as confirmed using computational fluid dynamics (CFD) simulations. The phantom was physically fabricated via lost-core casting with polyvinyl alcohol cryogel (PVA) as the tissue mimic. 3D ultrasound confirmed that the phantom lumen expectedly resembled a helical toroid geometry. Pulsed Doppler measurements showed that the phantom, when operating under steady flow conditions (3 ml s^-1 flow rate), yielded flow velocity magnitudes that agreed well with those derived from CFD at both the inner torus (-47.6 ± 5.7 versus -52.0 ± 2.2 cm s^-1; mean ± 1 S.D.) and the outer torus (49.5 ± 4.2 versus 48.0 ± 1.7 cm s^-1). Additionally, 3-C velocity vectors acquired from multi-angle pulsed Doppler showed good agreement with CFD-derived velocity vectors (<7% and 10° difference in magnitude and flow angle, respectively). Ultrasound color flow imaging further revealed that the phantom's axial flow pattern was aligned with the CFD-derived flow profile. Overall, the helical toroid phantom has strong potential as an investigative tool in 3D flow imaging innovation endeavors, such as the development of flow vector estimators and visualization algorithms.

3D synthetic aperture imaging with a therapeutic spherical random phased array for transcostal applications

Experimental validation of a synthetic aperture imaging technique using a therapeutic random phased array is described, demonstrating the dual nature of imaging and therapy of such an array. The transducer is capable of generating both continuous wave high intensity beams for ablating the tumor and low intensity ultrasound pulses to image the target area. Pulse-echo data is collected from the elements of the phased array to obtain B-mode images of the targets. Since therapeutic arrays are optimized for therapy only with concave apertures having low f-number and large directive elements often coarsely sampled, imaging can not be performed using conventional beamforming. We show that synthetic aperture imaging is capable of processing the acquired RF data to obtain images of the field of interest. Simulations were performed to compare different synthetic aperture imaging techniques to identify the best algorithm in terms of spatial resolution. Experimental validation was performed using a 1 MHz, 256-elements, spherical random phased array with 130 mm radius of curvature. The array was integrated with a research ultrasound scanner via custom connectors to acquire raw RF data for variety of targets. Imaging was implemented using synthetic aperture beamforming to produce images of a rib phantom and ex vivo ribs. The array was shown to resolve spherical targets within ±15 mm of either side of the axis in the focal plane and obtain 3D images of the rib phantom up to ±40 mm of either side of the central axis and at a depth of 3-9 cm from the array surface. The lateral and axial full width half maximum was 1.15 mm and 2.75 mm, respectively. This study was undertaken to emphasize that both therapy and image guidance with a therapeutic random phased array is possible and such a system has the potential to address some major limitations in the existing high intensity focused ultrasound (HIFU) systems. The 3D images obtained with a therapeutic array can be used to identify and locate strong scattering objects aiding to image guidance and treatment planning of the HIFU procedure.

Real-Time Passive Acoustic Mapping Using Sparse Matrix Multiplication

Passive acoustic mapping enables the spatiotemporal monitoring of cavitation with circulating microbubbles during focused ultrasound (FUS)-mediated blood-brain barrier opening. However, the computational load for processing large data sets of cavitation maps or more complex algorithms limit the visualization in real-time for treatment monitoring and adjustment. In this study, we implemented a graphical processing unit (GPU)-accelerated sparse matrix-based beamforming and time exposure acoustics in a neuronavigation-guided ultrasound system for real-time spatiotemporal monitoring of cavitation. The system performance was tested in silico through benchmarking, in vitro using nonhuman primate (NHP) and human skull specimens, and demonstrated in vivo in NHPs. We demonstrated the stability of the cavitation map for integration times longer than 62.5 [Formula: see text]. A compromise between real-time displaying and cavitation map quality obtained from beamformed RF data sets with a size of 2000 ×128 ×30 (axial [Formula: see text]) was achieved for an integration time of [Formula: see text], which required a computational time of 0.27 s (frame rate of 3.7 Hz) and could be displayed in real-time between pulses at PRF = 2 Hz. Our benchmarking tests show that the GPU sparse-matrix algorithm processed the RF data set at a computational rate of [Formula: see text]/pixel/sample, which enables adjusting the frame rate and the integration time as needed. The neuronavigation system with real-time implementation of cavitation mapping facilitated the localization of the cavitation activity and helped to identify distortions due to FUS phase aberration. The in vivo test of the method demonstrated the feasibility of GPU-accelerated sparse matrix computing in a close to a clinical condition, where focus distortions exemplify problems during treatment. These experimental conditions show the need for spatiotemporal monitoring of cavitation with real-time capability that enables the operator to correct or halt the sonication in case substantial aberrations are observed.

High-frame-rate volume imaging using sparse-random-aperture compounding

High-frame-rate volume imaging (HFR-VI) aims to provide high-quality images with high-temporal information. Despite its potential, HFR-VI translation into clinical applications has been challenging due to the high cost of the equipment required to drive matrix probes with a large number of elements. The goal of this study is to introduce and test sparse-random-aperture compounding (SRAC), a technique that allows use of matrix probes with an ultrasound system that has fewer channels while maintaining high frame rates. Four scanning methods were implemented with a 256-channel system using a 4-to-1 multiplexer and a 3 MHz matrix probe with 1024 elements. These methods used three types of waves, either single-diverging waves (SDW), multiplane-diverging waves (MDW) or wide beams (WB); and were driven using one to four SRAC. All methods were also implemented in a 1024-channel multisystem. The main-lobe-to-side-lobe ratio (MLSLR) and the contrast ratio (CR) were studied using a string phantom and a CIRS phantom, respectively. The results showed an increase in the MLSLR and CR as a function of the number of SRAC. The multisystem provided the best results for the MLSLR. However, four SRAC outperformed the multisystem with respect to CR. The method using SDW provided the highest frame rates (i.e. 1875 and 7500 Hz for four and one SRAC, respectively), however it provided the lowest image quality. The two methods using MDWs showed a good compromise between image quality and frame rate (i.e. 187 to 750 Hz for four and one SRAC). WB provided the best image quality at the expense of frame rate (i.e. 18 to 75 Hz for four and one SRAC). Our results suggest that SRAC in combination with the tested scanning methods can provide a low-channel count alternative for HFR-VI systems and allows a tunable tradeoff between image quality and frame rate guided by the desired application.

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.

Source Density Apodization: Image Artifact Suppression Through Source Pitch Nonuniformity

Conventional ultrasound imaging probes typically comprise finite-sized arrays of periodically spaced transducer elements which, in the case of phased arrays, can result in severe grating and sidelobe artifacts. Whereas side lobes can be effectively suppressed through amplitude apodization ("AmpA"), grating lobes arising from periodicity in transducer placement can only be suppressed by decreasing the element pitch, which is technologically challenging and costly. In this work, we present source density apodization ("SDA") as an alternative apodization scheme, where the spatial source density (and, hence, the element pitch) is varied across the imaging aperture. Using an all-optical ultrasound imaging setup capable of video-rate 2-D imaging as well as dynamic and arbitrary reconfiguration of the source array geometry, we show both numerically and experimentally how SDA and AmpA are equivalent for large numbers of sources. For low numbers of sources, SDA is shown to yield superior image quality as both side and grating lobes are effectively suppressed. In addition, we demonstrate how asymmetric SDA schemes can be used to locally and dynamically improve the image quality. Finally, we demonstrate how a nonsmoothly varying spatial source density (such as that obtained for randomized arrays or in the presence of source positioning uncertainty or inaccuracy) can yield severe image artifacts. The application of SDA can, thus, yield high image quality even for low channel counts, which can ultimately result in higher imaging frame rates using acquisition systems of reduced complexity.

Sparse hand-held probe for optoacoustic ultrasound volumetric imaging: an experimental proof-of-concept study

We present an experimental proof-of-concept study on the performance of a sparse segmented annular array for optoacoustic imaging. A capacitive micromachined ultrasonic transducer was equipped with a negatively focused acoustic lens and scanned in an annular fashion to exploit the performance of the sparse array geometry proposed in our recent numerical studies [Biomed. Opt. Express10, 1545 (2019)BOEICL2156-708510.1364/BOE.10.001545; J. Biomed. Opt.23, 025004 (2018)JBOPFO1083-366810.1117/1.JBO.23.2.025004]. A dedicated water tank was made using a 3D printer for light delivery and mounting the sample. A phantom experiment was carried out to showcase the possibility of full-field optoacoustic ultrasound (OPUS) imaging and confirm the earlier numerical results. This proof of concept opens the door towards a prototype of OPUS imaging for (pre-) clinical studies.

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.

3-D Super-Resolution Ultrasound Imaging With a 2-D Sparse Array

High-frame-rate 3-D ultrasound imaging technology combined with super-resolution processing method can visualize 3-D microvascular structures by overcoming the diffraction-limited resolution in every spatial direction. However, 3-D super-resolution ultrasound imaging using a full 2-D array requires a system with a large number of independent channels, the design of which might be impractical due to the high cost, complexity, and volume of data produced. In this study, a 2-D sparse array was designed and fabricated with 512 elements chosen from a density-tapered 2-D spiral layout. High-frame-rate volumetric imaging was performed using two synchronized ULA-OP 256 research scanners. Volumetric images were constructed by coherently compounding nine-angle plane waves acquired at a pulse repetition frequency of 4500 Hz. Localization-based 3-D super-resolution images of two touching subwavelength tubes were generated from 6000 volumes acquired in 12 s. Finally, this work demonstrates the feasibility of 3-D super-resolution imaging and super-resolved velocity mapping using a customized 2-D sparse array transducer.

Spectral Doppler Measurements With 2-D Sparse Arrays

The 2-D sparse arrays, in which a few hundreds of elements are distributed on the probe surface according to an optimization procedure, represent an alternative to full 2-D arrays, including thousands of elements usually organized in a grid. Sparse arrays have already been used in B-mode imaging tests, but their application to Doppler investigations has not been reported yet. Since the sparsity of the elements influences the acoustic field, a corresponding influence on the mean frequency (Fm), bandwidth (BW), and signal-to-noise ratio (SNR) of the Doppler spectra is expected. This article aims to assess, by simulations and experiments, to what extent the use of a sparse rather than a full gridded 2-D array has an impact on spectral Doppler measurements. Parabolic flows were investigated by a 3 MHz, 1024-element gridded array and by a sparse array; the latter was obtained by properly selecting a subgroup of 256 elements from the full array. Simulations show that the mean Doppler frequency does not change between the sparse and the full array while there are significant differences on the BW (average reduction of 17.2% for the sparse array, due to different apertures of the two probes) and on the signal power (Ps) (22 dB, due to the different number of active elements). These results are confirmed by flow phantom experiments, which also highlight that the most critical difference between sparse and full gridded array in Doppler measurements is in terms of SNR (-16.8 dB).

High-Frame-Rate Tri-Plane Echocardiography With Spiral Arrays: From Simulation to Real-Time Implementation

Major cardiovascular diseases (CVDs) are associated with (regional) dysfunction of the left ventricle. Despite the 3-D nature of the heart and its dynamics, the assessment of myocardial function is still largely based on 2-D ultrasound imaging, thereby making diagnosis heavily susceptible to the operator's expertise. Unfortunately, to date, 3-D echocardiography cannot provide adequate spatiotemporal resolution in real-time. Hence, tri-plane imaging has been introduced as a compromise between 2-D and true volumetric ultrasound imaging. However, tri-plane imaging typically requires high-end ultrasound systems equipped with fully populated matrix array probes embedded with expensive and little flexible electronics for two-stage beamforming. This article presents an advanced ultrasound system for real-time, high frame rate (HFR), and tri-plane echocardiography based on low element count sparse arrays, i.e., the so-called spiral arrays. The system was simulated, experimentally validated, and implemented for real-time operation on the ULA-OP 256 system. Five different array configurations were tested together with four different scan sequences, including multi-line and planar diverging wave transmission. In particular, the former can be exploited to achieve, in tri-plane imaging, the same temporal resolution currently used in clinical 2-D echocardiography, at the expenses of contrast (-3.5 dB) and signal-to-noise ratio (SNR) (-8.7 dB). On the other hand, the transmission of planar diverging waves boosts the frame rate up to 250 Hz, but further compromises contrast (-10.5 dB), SNR (-9.7 dB), and lateral resolution (+46%). In conclusion, despite an unavoidable loss in image quality and sensitivity due to the limited number of elements, HFR tri-plane imaging with spiral arrays is shown to be feasible in real-time and may enable real-time functional analysis of all left ventricular segments of the heart.

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

In three-dimensional (3D) medical ultrasound imaging with two-dimensional (2D) arrays, sparse 2D arrays have been studied to reduce the number of active channels. Among them, sparse 2D arrays with regular or uniform arrangements of elements have advantages of low side lobe energy and uniform field responses over the entire field of view. This paper presents two uniform sparse array models: sparse rectangular arrays (SRAs) on a rectangular grid and sparse spiral arrays (SSAs) on a sunflower grid. Both arrays can be easily implemented on the commercially available or the custom-made arrays. To suppress the overall grating lobe levels, the transmit (Tx) and receive (Rx) array pairs of both the array models are designed not to have grating lobes at the same locations in the Tx/Rx beam patterns, for which the theoretical design rules are also proposed. Computer simulation results indicate that the proposed array pairs for both the SRAs and the SSAs achieve peak grating lobe levels below -40 dB using about a quarter of the number of elements in the dense rectangular array while maintaining similar beam widths to that of the dense array pair.

Enhancing Image Quality of Photoacoustic Tomography Using Sub-Pitch Array Translation Approach: Simulation and Experimental Validation

OBJECTIVE

The purpose of this paper is to present a more convenient and practical alternate way of increasing the lateral discrete array sampling while using a typical λ pitch linear array transducer at receive for photoacoustic tomography (PAT) application.

METHODS

We have employed a linear translation-based approach, in which the array transducer is translated by sub-pitch amount to create an augmented RF frame data having denser lateral spatial sampling. The denser λ/2 and λ/4 pitch data were reconstructed and compared against conventional λ pitch reconstructed PAT image using simulation and tissue mimicking phantom experiments in terms of improvements in resolution and contrast.

RESULTS

The results from experiments demonstrate a 34.48% improvement in lateral resolution (LR), measured in terms of full-width at half-maximum of the lateral profile of point spread function, and a maximum of 7-dB improvement in contrast is achieved while using a λ/2-pitch configuration when compared to the conventional λ-pitch configuration.

CONCLUSION

It was demonstrated that λ/2- and λ/4-pitch configurations result in better LR and contrast than λ-pitch configuration.

SIGNIFICANCE

Based on the results obtained, the proposed method has the potential to serve as an easy-to-integrate and simple way of achieving better image quality without requiring to increase the system complexity with existing transducer array probe technology in regular clinical scanners.

High-Frame-Rate Contrast-Enhanced Echocardiography Using Diverging Waves: 2-D Motion Estimation and Compensation

Combining diverging ultrasound waves and microbubbles could improve contrast-enhanced echocardiography (CEE), by providing enhanced temporal resolution for cardiac function assessment over a large imaging field of view. However, current image formation techniques using coherent summation of echoes from multiple steered diverging waves (DWs) are susceptible to tissue and microbubble motion artifacts, resulting in poor image quality. In this study, we used correlation-based 2-D motion estimation to perform motion compensation for CEE using DWs. The accuracy of this motion estimation method was evaluated with Field II simulations. The root-mean-square velocity errors were 5.9% ± 0.2% and 19.5% ± 0.4% in the axial and lateral directions, when normalized to the maximum value of 62.8 cm/s which is comparable to the highest speed of blood flow in the left ventricle (LV). The effects of this method on image contrast ratio (CR) and contrast-to-noise ratio (CNR) were tested in vitro using a tissue mimicking rotating disk with a diameter of 10 cm. Compared against the control without motion compensation, a mean increase of 12 dB in CR and 7 dB in CNR were demonstrated when using this motion compensation method. The motion correction algorithm was tested in vivo on a CEE data set acquired with the Ultrasound Array Research Platform II performing coherent DW imaging. Improvement of the B-mode and contrast-mode image quality with cardiac motion and blood flow-induced microbubble motion was achieved. The results of motion estimation were further processed to interpret blood flow in the LV. This allowed for a triplex cardiac imaging technique, consisting of B mode, contrast mode, and 2-D vector flow imaging with a high frame rate of 250 Hz.