An Ultrasound Matrix Transducer for High-Frame-Rate 3-D Intra-cardiac Echocardiography

OBJECTIVE

Described here is the development of an ultrasound matrix transducer prototype for high-frame-rate 3-D intra-cardiac echocardiography.

METHODS

The matrix array consists of 16 × 18 lead zirconate titanate elements with a pitch of 160 µm × 160 µm built on top of an application-specific integrated circuit that generates transmission signals and digitizes the received signals. To reduce the number of cables in the catheter to a feasible number, we implement subarray beamforming and digitization in receive and use a combination of time-division multiplexing and pulse amplitude modulation data transmission, achieving an 18-fold reduction. The proposed imaging scheme employs seven fan-shaped diverging transmit beams operating at a pulse repetition frequency of 7.7 kHz to obtain a high frame rate. The performance of the prototype is characterized, and its functionality is fully verified.

RESULTS

The transducer exhibits a transmit efficiency of 28 Pa/V at 5 cm per element and a bandwidth of 60% in transmission. In receive, a dynamic range of 80 dB is measured with a minimum detectable pressure of 10 Pa per element. The element yield of the prototype is 98%, indicating the efficacy of the manufacturing process. The transducer is capable of imaging at a frame rate of up to 1000 volumes/s and is intended to cover a volume of 70° × 70° × 10 cm.

CONCLUSION

These advanced imaging capabilities have the potential to support complex interventional procedures and enable full-volumetric flow, tissue, and electromechanical wave tracking in the heart.

Computation of ultrasound propagation in a population of nonlinearly oscillating microbubbles including multiple scattering

In contrast-enhanced echography, the simulation of nonlinear propagation of ultrasound through a population of oscillating microbubbles imposes a computational challenge. Also, the numerical complexity increases because each scatterer has individual properties. To address these problems, the Iterative Nonlinear Contrast Source (INCS) method has been extended to include a large population of nonlinearly responding microbubbles. The original INCS method solves the Westervelt equation in a four-dimensional spatiotemporal domain by generating increasingly accurate field corrections to iteratively update the acoustic pressure. The field corrections are computed by the convolution of a nonlinear contrast source with the Green's function of the linear background medium. Because the convolution integral allows a coarse discretization, INCS can efficiently deal with large-scale problems. To include a population of microbubbles, these are considered as individual contrast point sources with their own nonlinear response. The field corrections are computed as before, but now, in each iteration, the temporal signature of each contrast point source is computed by solving the bubble's Marmottant equation. Physically, each iteration adds an order of multiple scattering. Here, the performance of the extended INCS method and the significance of multiple scattering is demonstrated through various results from different configurations.

Improving Lateral Resolution in 3-D Imaging With Micro-beamforming Through Adaptive Beamforming by Deep Learning

There is an increased desire for miniature ultrasound probes with small apertures to provide volumetric images at high frame rates for in-body applications. Satisfying these increased requirements makes simultaneous achievement of a good lateral resolution a challenge. As micro-beamforming is often employed to reduce data rate and cable count to acceptable levels, receive processing methods that try to improve spatial resolution will have to compensate the introduced reduction in focusing. Existing beamformers do not realize sufficient improvement and/or have a computational cost that prohibits their use. Here we propose the use of adaptive beamforming by deep learning (ABLE) in combination with training targets generated by a large aperture array, which inherently has better lateral resolution. In addition, we modify ABLE to extend its receptive field across multiple voxels. We illustrate that this method improves lateral resolution both quantitatively and qualitatively, such that image quality is improved compared with that achieved by existing delay-and-sum, coherence factor, filtered-delay-multiplication-and-sum and Eigen-based minimum variance beamformers. We found that only in silica data are required to train the network, making the method easily implementable in practice.

Coupling Two Ultra-high-Speed Cameras to Elucidate Ultrasound Contrast-Mediated Imaging and Therapy

Ultrasound contrast-mediated medical imaging and therapy both rely on the dynamics of micron- and nanometer-sized ultrasound cavitation nuclei, such as phospholipid-coated microbubbles and phase-change droplets. Ultrasound cavitation nuclei respond non-linearly to ultrasound on a nanosecond time scale that necessitates the use of ultra-high-speed imaging to fully visualize these dynamics in detail. In this study, we developed an ultra-high-speed optical imaging system that can record up to 20 million frames per second (Mfps) by coupling two small-sized, commercially available, 10-Mfps cameras. The timing and reliability of the interleaved cameras needed to achieve 20 Mfps was validated using two synchronized light-emitting diode strobe lights. Once verified, ultrasound-activated microbubble responses were recorded and analyzed. A unique characteristic of this coupled system is its ability to be reconfigured to provide orthogonal observations at 10 Mfps. Acoustic droplet vaporization was imaged from two orthogonal views, by which the 3-D dynamics of the phase transition could be visualized. This optical imaging system provides the temporal resolution and experimental flexibility needed to further elucidate the dynamics of ultrasound cavitation nuclei to potentiate the clinical translation of ultrasound-mediated imaging and therapy developments.

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.

Algorithm to Correct Measurement Offsets Introduced by Inactive Elements of Transducer Arrays in Ultrasonic Flow Metering

Ultrasonic flow meters (UFMs) based on transducer arrays offer several advantages. With electronic beam steering, it is possible to tune the steering angle of the beam for optimal signal-tonoise ratio (SNR) upon reception. Moreover, multiple beams can be generated to propagate through different travel paths, covering a wider section of the flow profile. Furthermore, in a clamp-on configuration, UFMs based on transducer arrays can perform self-calibration. In this manner, userinput is minimized and measurement repeatability is increased. In practice, transducer array elements may break down. This could happen due to aging, exposure to rough environments, and/or rough mechanical contact. As a consequence of inactive array elements, the measured transit time difference contains two offsets. One offset originates from non-uniform spatial sampling of the generated wavefield. Another offset originates from the ill-defined beam propagating through a travel path different from the intended one. In this paper, an algorithm is proposed that corrects for both of these offsets. The algorithm also performs a filtering operation in the frequency-wavenumber domain of all spurious (i.e., flow-insensitive) wave modes. The advantage of implementing the proposed algorithm is demonstrated on simulations and measurements, showing improved accuracy and precision of the transit time differences compared to the values obtained when the algorithm is not applied. The proposed algorithm can be implemented in both in-line and clamp-on configuration of UFMs based on transducer arrays.

Imaging Scheme for 3-D High-Frame-Rate Intracardiac Echography: A Simulation Study

Atrial fibrillation (AF) is the most common cardiac arrhythmia and is normally treated by RF ablation. Intracardiac echography (ICE) is widely employed during RF ablation procedures to guide the electrophysiologist in navigating the ablation catheter, although only 2-D probes are currently clinically used. A 3-D ICE catheter would not only improve visualization of the atrium and ablation catheter, but it might also provide the 3-D mapping of the electromechanical wave (EW) propagation pattern, which represents the mechanical response of cardiac tissue to electrical activity. The detection of this EW needs 3-D high-frame-rate imaging, which is generally only realizable in tradeoff with channel count and image quality. In this simulation-based study, we propose a high volume rate imaging scheme for a 3-D ICE probe design that employs 1-D micro-beamforming in the elevation direction. Such a probe can achieve a high frame rate while reducing the channel count sufficiently for realization in a 10-Fr catheter. To suppress the grating-lobe (GL) artifacts associated with micro-beamforming in the elevation direction, a limited number of fan-shaped beams with a wide azimuthal and narrow elevational opening angle are sequentially steered to insonify slices of the region of interest. An angular weighted averaging of reconstructed subvolumes further reduces the GL artifacts. We optimize the transmit beam divergence and central frequency based on the required image quality for EW imaging (EWI). Numerical simulation results show that a set of seven fan-shaped transmission beams can provide a frame rate of 1000 Hz and a sufficient spatial resolution to visualize the EW propagation on a large 3-D surface.

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.

Accelerated 2-D Real-Time Refraction-Corrected Transcranial Ultrasound Imaging

In a recent study, we proposed a technique to correct aberration caused by the skull and reconstruct a transcranial B-mode image with a refraction-corrected synthetic aperture imaging (SAI) scheme. Given a sound speed map, the arrival times were calculated using a fast marching technique (FMT), which solves the Eikonal equation and, therefore, is computationally expensive for real-time imaging. In this article, we introduce a two-point ray tracing method, based on Fermat's principle, for fast calculation of the travel times in the presence of a layered aberrator in front of the ultrasound probe. The ray tracing method along with the reconstruction technique is implemented on a graphical processing unite (GPU). The point spread function (PSF) in a wire phantom image reconstructed with the FMT and the GPU implementation was studied with numerical synthetic data and experiments with a bone-mimicking plate and a sagittally cut human skull. The numerical analysis showed that the error on travel times is less than 10% of the ultrasound temporal period at 2.5 MHz. As a result, the lateral resolution was not significantly degraded compared with images reconstructed with FMT-calculated travel times. The results using the synthetic, bone-mimicking plate, and skull dataset showed that the GPU implementation causes a lateral/axial localization error of 0.10/0.20, 0.15/0.13, and 0.26/0.32 mm compared with a reference measurement (no aberrator in front of the ultrasound probe), respectively. For an imaging depth of 70 mm, the proposed GPU implementation allows reconstructing 19 frames/s with full synthetic aperture (96 transmission events) and 32 frames/s with multiangle plane wave imaging schemes (with 11 steering angles) for a pixel size of [Formula: see text]. Finally, refraction-corrected power Doppler imaging is demonstrated with a string phantom and a bone-mimicking plate placed between the probe and the moving string. The proposed approach achieves a suitable frame rate for clinical scanning while maintaining the image quality.

Design and Proof-of-Concept of a Matrix Transducer Array for Clamp-On Ultrasonic Flow Measurements

Common clamp-on ultrasonic flow meters consist of two single-element transducers placed on the pipe wall. Flow speed is measured noninvasively, i.e., without interrupting the flow and without perforating the pipe wall, which also minimizes safety risks and avoids pressure drops inside the pipe. However, before metering, the transducers have to be carefully positioned along the pipe axis to correctly align the acoustic beams and obtain a well-calibrated flowmeter. This process is done manually, is dependent on the properties of the pipe and the liquid, does not account for pipe imperfections, and becomes troublesome on pipelines with an intricate shape. Matrix transducer arrays are suitable to dynamically steer acoustic beams and realize self-alignment upon reception, without user input. In this work, the design of a broadband 37×17 matrix array (center frequency of 1 MHz) to perform clamp-on ultrasonic flow measurements over a wide range of liquids ( c=1000-2000 m/s, α ≤ 1 dB/MHz · cm) and pipe sizes is presented. Three critical aspects were assessed: efficiency, electronic beam steering, and wave mode conversion in the pipe wall. A prototype of a proof-of-concept flowmeter consisting of two 36-element linear arrays (center frequency of 1.1 MHz) was fabricated and placed on a 1-mm-thick, 40-mm inner diameter stainless steel pipe in a custom-made flow loop filled with water. At resonance, simulated and measured efficiencies in water of the linear arrays compared well: 0.88 and 0.81 kPa/V, respectively. Mean flow measurements were achieved by electronic beam steering of the acoustic beams and using both compressional and shear waves generated in the pipe wall. Correlation coefficients of between measured and reference flow speeds were obtained, thus showing the operational concept of an array-based clamp-on ultrasonic flowmeter.

Measurement of Pipe and Liquid Parameters Using the Beam Steering Capabilities of Array-Based Clamp-On Ultrasonic Flow Meters

Clamp-on ultrasonic flow meters (UFMs) are installed on the outside of the pipe wall. Typically, they consist of two single-element transducers mounted on angled wedges, which are acoustically coupled to the pipe wall. Before flow metering, the transducers are placed at the correct axial position by manually moving one transducer along the pipe wall until the maximum amplitude of the relevant acoustic pulse is obtained. This process is time-consuming and operator-dependent. Next to this, at least five parameters of the pipe and the liquid need to be provided manually to compute the flow speed. In this work, a method is proposed to obtain the five parameters of the pipe and the liquid required to compute the flow speed. The method consists of obtaining the optimal angles for different wave travel paths by varying the steering angle of the emitted acoustic beam systematically. Based on these optimal angles, a system of equations is built and solved to extract the desired parameters. The proposed method was tested experimentally with a custom-made clamp-on UFM consisting of two linear arrays placed on a water-filled stainless steel pipe. The obtained parameters of the pipe and the liquid correspond very well with the expected (nominal) values. Furthermore, the performed experiment also demonstrates that a clamp-on UFM based on transducer arrays can achieve self-alignment without the need to manually move the transducers.

Error analysis and reliability of zero-order Lamb mode inversion for waveguide characterization

In recent years, several fitting techniques have been presented to reconstruct the parameters of a plate from its Lamb wave dispersion curves. Published studies show that these techniques can yield high accuracy results and have the potential of reconstructing several parameters at once. The precision with which parameters can be reconstructed by inverting Lamb wave dispersion curves, however, remains an open question of fundamental importance to many applications. In this work, we introduce a method of analyzing dispersion curves that yields quantitative information on the precision with which the parameters can be extracted. In our method, rather than employing error minimization algorithms, we compare a target dispersion curve to a database of theoretical ones that covers a given parameter space. By calculating a measure of dissimilarity (error) for every point in the parameter space, we reconstruct the distribution of the error in that space, beside the location of its minimum. We then introduce dimensionless quantities that describe the distribution of this error, thus yielding information about the spread of similar curves in the parameter space. We demonstrate our approach by considering both idealized and realistic scenarios, analyzing the dispersion curves obtained numerically for a plate and experimentally for a pipe. Our results show that the precision with which each parameter is reconstructed depends on the mode used, as well as the frequency range in which it is considered.

Time-resolved absolute radius estimation of vibrating contrast microbubbles using an acoustical camera

Ultrasound (US) contrast agents consist of microbubbles ranging from 1 to 10 μm in size. The acoustical response of individual microbubbles can be studied with high-frame-rate optics or an "acoustical camera" (AC). The AC measures the relative microbubble oscillation while the optical camera measures the absolute oscillation. In this article, the capabilities of the AC are extended to measure the absolute oscillations. In the AC setup, microbubbles are insonified with a high- (25 MHz) and low-frequency US wave (1-2.5 MHz). Other than the amplitude modulation (AM) from the relative size change of the microbubble (employed in Renaud, Bosch, van der Steen, and de Jong (2012a). "An 'acoustical camera' for in vitro characterization of contrast agent microbubble vibrations," Appl. Phys. Lett. 100(10), 101911, the high-frequency response from individual vibrating microbubbles contains a phase modulation (PM) from the microbubble wall displacement, which is the extension described here. The ratio of PM and AM is used to determine the absolute radius, R₀. To test this sizing, the size distributions of two monodisperse microbubble populations ( _R = 2.1 and 3.5 μm) acquired with the AC were matched to the distribution acquired with a Coulter counter. As a result of measuring the absolute size of the microbubbles, this "extended AC" can capture the full radial dynamics of single freely floating microbubbles with a throughput of hundreds of microbubbles per hour.

Internalization of targeted microbubbles by endothelial cells and drug delivery by pores and tunnels

Ultrasound insonification of microbubbles can locally enhance drug delivery by increasing the cell membrane permeability. To aid development of a safe and effective therapeutic microbubble, more insight into the microbubble-cell interaction is needed. In this in vitro study we aimed to investigate the initial 3D morphology of the endothelial cell membrane adjacent to individual microbubbles (n = 301), determine whether this morphology was affected upon binding and by the type of ligand on the microbubble, and study its influence on microbubble oscillation and the drug delivery outcome. High-resolution 3D confocal microscopy revealed that targeted microbubbles were internalized by endothelial cells, while this was not the case for non-targeted or IgG1-κ control microbubbles. The extent of internalization was ligand-dependent, since α_vβ₃-targeted microbubbles were significantly more internalized than CD31-targeted microbubbles. Ultra-high-speed imaging (~17 Mfps) in combination with high-resolution confocal microscopy (n = 246) showed that microbubble internalization resulted in a damped microbubble oscillation upon ultrasound insonification (2 MHz, 200 kPa peak negative pressure, 10 cycles). Despite damped oscillation, the cell's susceptibility to sonoporation (as indicated by PI uptake) was increased for internalized microbubbles. Monitoring cell membrane integrity (n = 230) showed the formation of either a pore, for intracellular delivery, or a tunnel (i.e. transcellular perforation), for transcellular delivery. Internalized microbubbles caused fewer transcellular perforations and smaller pore areas than non-internalized microbubbles. In conclusion, studying microbubble-mediated drug delivery using a state-of-the-art imaging system revealed receptor-mediated microbubble internalization and its effect on microbubble oscillation and resulting membrane perforation by pores and tunnels.

An iterative method to evaluate one-dimensional pulsed nonlinear elastic wavefields and mixing of elastic waves in solids

Over the last 15 years, literature on nondestructive testing has shown that the generation of higher harmonics and nonlinear mixing of waves could be used to obtain the nonlinearity parameters of an elastic medium and thereby gather information about its state, e.g., aging and fatigue. To design ultrasound measurement setups based on these phenomena, efficient numerical modeling tools are needed. In this paper, the iterative nonlinear contrast source method for numerical modeling of nonlinear acoustic waves is extended to the one-dimensional elastic case. In particular, nonlinear mixing of two collinear bulk waves (one compressional, one shear) in a homogeneous, isotropic medium is considered, taking into account its third-order elastic constants ( A,  B,  and  C). The obtained results for nonlinear propagation are in good agreement with a benchmark solution based on the modified Burgers equation. The results for the resonant waves that are caused by the one-way and two-way mixing of primary waves are in quantitative agreement with the results in the literature [Chen, Tang, Zhao, Jacobs, and Qu, J. Acoust. Soc. Am. 136(5), 2389-2404 (2014)]. The contrast source approach allows the identification of the propagating and evanescent components of the scattered wavefield in the wavenumber-frequency domain, which provides physical insight into the mixing process and explains the propagation direction of the resonant wave.

Refraction-Corrected Transcranial Ultrasound Imaging Through the Human Temporal Window Using a Single Probe

Transcranial ultrasound imaging (TUI) is a diagnostic modality with numerous applications, but unfortunately, it is hindered by phase aberration caused by the skull. In this article, we propose to reconstruct a transcranial B-mode image with a refraction-corrected synthetic aperture imaging (SAI) scheme. First, the compressional sound velocity of the aberrator (i.e., the skull) is estimated using the bidirectional headwave technique. The medium is described with four layers (i.e., lens, water, skull, and water), and a fast marching method calculates the travel times between individual array elements and image pixels. Finally, a delay-and-sum algorithm is used for image reconstruction with coherent compounding. The point spread function (PSF) in a wire phantom image and reconstructed with the conventional technique (using a constant sound speed throughout the medium), and the proposed method was quantified with numerical synthetic data and experiments with a bone-mimicking plate and a human skull, compared with the PSF achieved in a ground truth image of the medium without the aberrator (i.e., the bone plate or skull). A phased-array transducer (P4-1, ATL/Philips, 2.5 MHz, 96 elements, pitch = 0.295 mm) was used for the experiments. The results with the synthetic signals, the bone-mimicking plate, and the skull indicated that the proposed method reconstructs the scatterers with an average lateral/axial localization error of 0.06/0.14 mm, 0.11/0.13 mm, and 1.0/0.32 mm, respectively. With the human skull, an average contrast ratio (CR) and full-width-half-maximum (FWHM) of 37.1 dB and 1.75 mm were obtained with the proposed approach, respectively. This corresponds to an improvement of CR and FWHM by 7.1 dB and 36% compared with the conventional method, respectively. These numbers were 12.7 dB and 41% with the bone-mimicking plate.

A spatial and temporal characterisation of single proton acoustic waves in proton beam cancer therapy

An in vivo range verification technology for proton beam cancer therapy, preferably in real-time and with submillimeter resolution, is desired to reduce the present uncertainty in dose localization. Acoustical imaging technologies exploiting possible local interactions between protons and microbubbles or nanodroplets might be an interesting option. Unfortunately, a theoretical model capable of characterising the acoustical field generated by an individual proton on nanometer and micrometer scales is still missing. In this work, such a model is presented. The proton acoustic field is generated by the adiabatic expansion of a region that is locally heated by a passing proton. To model the proton heat deposition, secondary electron production due to protons has been quantified using a semi-empirical model based on Rutherford's scattering theory, which reproduces experimentally obtained electronic stopping power values for protons in water within 10% over the full energy range. The electrons transfer energy into heat via electron-phonon coupling to atoms along the proton track. The resulting temperature increase is calculated using an inelastic thermal spike model. Heat deposition can be regarded as instantaneous, thus, stress confinement is ensured and acoustical initial conditions are set. The resulting thermoacoustic field in the nanometer and micrometer range from the single proton track is computed by solving the thermoacoustic wave equation using k-space Green's functions, yielding the characteristic amplitudes and frequencies present in the acoustic signal generated by a single proton in an aqueous medium. Wavefield expansion and asymptotic approximations are used to extend the spatial and temporal ranges of the proton acoustic field.

Measurement of Pipe and Fluid Properties With a Matrix Array-Based Ultrasonic Clamp-On Flow Meter

Current ultrasonic clamp-on flow meters consist of a pair of single-element transducers that are carefully positioned before use. This positioning process consists of manually finding the distance between the transducer elements, along the pipe axis, for which maximum signal-to-noise ratio (SNR) is achieved. This distance depends on the sound speed, thickness, and diameter of the pipe and on the sound speed of the liquid. However, these parameters are either known with low accuracy or completely unknown during positioning, making it a manual and troublesome process. Furthermore, even when sensor positioning is done properly, uncertainty about the mentioned parameters, and therefore on the path of the acoustic beams, limits the final accuracy of flow measurements. In this research, we address these issues using an ultrasonic clamp-on flow meter consisting of two matrix arrays, which enables the measurement of pipe and liquid parameters by the flow meter itself. Automatic parameter extraction, combined with the beam-steering capabilities of transducer arrays, yields a sensor capable of compensating for pipe imperfections. Three parameter extraction procedures are presented. In contrast to similar literature, the procedures proposed here do not require that the medium be submerged nor do they require a priori information about it. First, axial Lamb waves are excited along the pipe wall and recorded with one of the arrays. A dispersion curve-fitting algorithm is used to extract bulk sound speeds and wall thickness of the pipe from the measured dispersion curves. Second, circumferential Lamb waves are excited, measured, and corrected for dispersion to extract the pipe diameter. Third, pulse-echo measurements provide the sound speed of the liquid. The effectiveness of the first two procedures has been evaluated using simulated and measured data of stainless steel and aluminum pipes, and the feasibility of the third procedure has been evaluated using simulated data.

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.

Exploiting nonlinear wave propagation to improve the precision of ultrasonic flow meters

Acoustic wave propagation in ultrasonic flow measurements is typically assumed to be linear and reciprocal. However, if the transmitting transducer generates a sufficiently high pressure, nonlinear wave propagation effects become significant. In flow measurements, this would translate into more information to estimate the flow and therefore a higher precision relative to the linear case. In this work, we investigate how the generated harmonics can be used to measure flow. Measurements in a custom-made flow loop and simulations using the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation will show that the second harmonic component provides similar transit time differences to those obtained from the fundamental component, their linear combination results in more precise flow measurements compared to the estimations with the fundamental component alone.

Tapering of the interventricular septum can affect ultrasound shear wave elastography: An in vitro and in silico study

Shear wave elastography (SWE) has the potential to determine cardiac tissue stiffness from non-invasive shear wave speed measurements, important, e.g., for predicting heart failure. Previous studies showed that waves traveling in the interventricular septum (IVS) may display Lamb-like dispersive behaviour, introducing a thickness-frequency dependency in the wave speed. However, the IVS tapers across its length, which complicates wave speed estimation by introducing an additional variable to account for. The goal of this work is to assess the impact of tapering thickness on SWE. The investigation is performed by combining in vitro experiments with acoustic radiation force (ARF) and 2D finite element simulations, to isolate the effect of the tapering curve on ARF-induced and natural waves in the heart. The experiments show a 11% deceleration during propagation from the thick to the thin end of an IVS-mimicking tapered phantom plate. The numerical analysis shows that neglecting the thickness variation in the wavenumber-frequency domain can introduce errors of more than 30% in the estimation of the shear modulus, and that the exact tapering curve, rather than the overall thickness reduction, determines the dispersive behaviour of the wave. These results suggest that septal geometry should be accounted for when deriving cardiac stiffness with SWE.

Motion-compensated noninvasive periodontal health monitoring using handheld and motor-based photoacoustic-ultrasound imaging systems

Simultaneous visualization of the teeth and periodontium is of significant clinical interest for image-based monitoring of periodontal health. We recently reported the application of a dual-modality photoacoustic-ultrasound (PA-US) imaging system for resolving periodontal anatomy and periodontal pocket depths in humans. This work utilized a linear array transducer attached to a stepper motor to generate 3D images via maximum intensity projection. This prior work also used a medical head immobilizer to reduce artifacts during volume rendering caused by motion from the subject (e.g., breathing, minor head movements). However, this solution does not completely eliminate motion artifacts while also complicating the imaging procedure and causing patient discomfort. To address this issue, we report the implementation of an image registration technique to correctly align B-mode PA-US images and generate artifact-free 2D cross-sections. Application of the deshaking technique to PA phantoms revealed 80% similarity to the ground truth when shaking was intentionally applied during stepper motor scans. Images from handheld sweeps could also be deshaken using an LED PA-US scanner. In ex vivo porcine mandibles, pigmentation of the enamel was well-estimated within 0.1 mm error. The pocket depth measured in a healthy human subject was also in good agreement with our prior study. This report demonstrates that a modality-independent registration technique can be applied to clinically relevant PA-US scans of the periodontium to reduce operator burden of skill and subject discomfort while showing potential for handheld clinical periodontal imaging.

Lamb Waves and Adaptive Beamforming for Aberration Correction in Medical Ultrasound Imaging

Phase aberration in transcranial ultrasound imaging (TUI) caused by the human skull leads to an inaccurate image reconstruction. In this article, we present a novel method for estimating the speed of sound and an adaptive beamforming technique for phase aberration correction in a flat polyvinylchloride (PVC) slab as a model for the human skull. First, the speed of sound of the PVC slab is found by extracting the overlapping quasi-longitudinal wave velocities of symmetrical Lamb waves in the frequency-wavenumber domain. Then, the thickness of the plate is determined by the echoes from its front and back side. Next, an adaptive beamforming method is developed, utilizing the measured sound speed map of the imaging medium. Finally, to minimize reverberation artifacts caused by strong scatterers (i.e., needles), a dual probe setup is proposed. In this setup, we image the medium from two opposite directions, and the final image can be the minimum intensity projection of the inherently co-registered images of the opposed probes. Our results confirm that the Lamb wave method estimates the longitudinal speed of the slab with an error of 3.5% and is independent of its shear wave speed. Benefiting from the acquired sound speed map, our adaptive beamformer reduces (in real time) a mislocation error of 3.1, caused by an 8 mm slab, to 0.1 mm. Finally, the dual probe configuration shows 7 dB improvement in removing reverberation artifacts of the needle, at the cost of only 2.4-dB contrast loss. The proposed image formation method can be used, e.g., to monitor deep brain stimulation procedures and localization of the electrode(s) deep inside the brain from two temporal bones on the sides of the human skull.

Erratum: Lamb Waves and Adaptive Beamforming for Aberration Correction in Medical Ultrasound Imaging

In our paper titled "Lamb Waves and Adaptive Beamforming for Aberration Correction in Medical Ultrasound Imaging" [1], we mentioned that the superposition of the different symmetric (S) modes in the frequency-wavenumber (f-k) domain results in a high intensity region where its slope corresponds to the longitudinal wave speed in the slab. However, we have recently understood that this high intensity region belongs to the propagation of a wave called lateral wave or head wave [2-5]. It is generated if the longitudinal sound speed of the aberrator (i.e. the PVC slab) is larger than that of water and if the incident wavefront is curved. When the incidence angle at the interface between water and PVC is near the critical angle, the refracted wave in PVC re-radiates a small part of its energy into the fluid (i.e. the head wave). As discussed in [4], if the thickness of the waveguide is larger than the wavelength, the first arriving signal is the head wave. This is also the case in our study [1] where the ultrasound wavelength of a compressional wave in PVC was close to 1 mm, and a PVC slab with a thickness of 8 mm was used.

Parasternal Versus Apical View in Cardiac Natural Mechanical Wave Speed Measurements

Shear wave speed measurements can potentially be used to noninvasively measure myocardial stiffness to assess the myocardial function. Several studies showed the feasibility of tracking natural mechanical waves induced by aortic valve closure in the interventricular septum, but different echocardiographic views have been used. This article systematically studied the wave propagation speeds measured in a parasternal long-axis and in an apical four-chamber view in ten healthy volunteers. The apical and parasternal views are predominantly sensitive to longitudinal or transversal tissue motion, respectively, and could, therefore, theoretically measure the speed of different wave modes. We found higher propagation speeds in apical than in the parasternal view (median of 5.1 m/s versus 3.8 m/s, , n = 9 ). The results in the different views were not correlated ( r = 0.26 , p = 0.49 ) and an unexpectedly large variability among healthy volunteers was found in apical view compared with the parasternal view (3.5-8.7 versus 3.2-4.3 m/s, respectively). Complementary finite element simulations of Lamb waves in an elastic plate showed that different propagation speeds can be measured for different particle motion components when different wave modes are induced simultaneously. The in vivo results cannot be fully explained with the theory of Lamb wave modes. Nonetheless, the results suggest that the parasternal long-axis view is a more suitable candidate for clinical diagnosis due to the lower variability in wave speeds.

High-Resolution Imaging of Intracellular Calcium Fluctuations Caused by Oscillating Microbubbles

Ultrasound insonification of microbubbles can locally enhance drug delivery, but the microbubble-cell interaction remains poorly understood. Because intracellular calcium (Ca_i²⁺) is a key cellular regulator, unraveling the Ca_i²⁺ fluctuations caused by an oscillating microbubble provides crucial insight into the underlying bio-effects. Therefore, we developed an optical imaging system at nanometer and nanosecond resolution that can resolve Ca_i²⁺ fluctuations and microbubble oscillations. Using this system, we clearly distinguished three Ca_i²⁺ uptake profiles upon sonoporation of endothelial cells, which strongly correlated with the microbubble oscillation amplitude, severity of sonoporation and opening of cell-cell contacts. We found a narrow operating range for viable drug delivery without lethal cell damage. Moreover, adjacent cells were affected by a calcium wave propagating at 15 μm/s. With the unique optical system, we unraveled the microbubble oscillation behavior required for drug delivery and Ca_i²⁺ fluctuations, providing new insight into the microbubble-cell interaction to aid clinical translation.

Ultrasound-Responsive Cavitation Nuclei for Therapy and Drug Delivery

Therapeutic ultrasound strategies that harness the mechanical activity of cavitation nuclei for beneficial tissue bio-effects are actively under development. The mechanical oscillations of circulating microbubbles, the most widely investigated cavitation nuclei, which may also encapsulate or shield a therapeutic agent in the bloodstream, trigger and promote localized uptake. Oscillating microbubbles can create stresses either on nearby tissue or in surrounding fluid to enhance drug penetration and efficacy in the brain, spinal cord, vasculature, immune system, biofilm or tumors. This review summarizes recent investigations that have elucidated interactions of ultrasound and cavitation nuclei with cells, the treatment of tumors, immunotherapy, the blood-brain and blood-spinal cord barriers, sonothrombolysis, cardiovascular drug delivery and sonobactericide. In particular, an overview of salient ultrasound features, drug delivery vehicles, therapeutic transport routes and pre-clinical and clinical studies is provided. Successful implementation of ultrasound and cavitation nuclei-mediated drug delivery has the potential to change the way drugs are administered systemically, resulting in more effective therapeutics and less-invasive treatments.

Fundamental modeling of wave propagation in temporally relaxing media with applications to cardiac shear wave elastography

Shear wave elastography (SWE) might allow non-invasive assessment of cardiac stiffness by relating shear wave propagation speed to material properties. However, after aortic valve closure, when natural shear waves occur in the septal wall, the stiffness of the muscle decreases significantly, and the effects of such temporal variation of medium properties on shear wave propagation have not been investigated yet. The goal of this work is to fundamentally investigate these effects. To this aim, qualitative results were first obtained experimentally using a mechanical setup, and were then combined with quantitative results from finite difference simulations. The results show that the amplitude and period of the waves increase during propagation, proportional to the relaxation of the medium, and that reflected waves can originate from the temporal stiffness variation. These general results, applied to literature data on cardiac stiffness throughout the heart cycle, predict as a major effect a period increase of 20% in waves propagating during a healthy diastolic phase, whereas only a 10% increase would result from the impaired relaxation of an infarcted heart. Therefore, cardiac relaxation can affect the propagation of waves used for SWE measurements and might even provide direct information on the correct relaxation of a heart.

Suppression of Lamb wave excitation via aperture control of a transducer array for ultrasonic clamp-on flow metering

During ultrasonic clamp-on flow metering, Lamb waves propagating in the pipe wall may limit the measurement accuracy by introducing absolute errors in the flow estimates. Upon reception, these waves can interfere with the up and downstream waves refracting from the liquid, and disturb the measurement of the transit time difference that is used to obtain the flow speed. Thus, suppression of the generation of Lamb waves might directly increase the accuracy of a clamp-on flow meter. Existing techniques apply to flow meters with single element transducers. This paper considers the application of transducer arrays and presents a method to achieve a predefined amount of suppression of these spurious Lamb waves based on appropriate amplitude weightings of the transducer elements. Finite element simulations of an ultrasonic clamp-on flow measurement setting will be presented to show the effect of array aperture control on the suppression of the Lamb waves in a 1-mm-thick stainless steel pipe wall. Furthermore, a proof-of-principle experiment will be shown that demonstrates a good agreement with the simulations.

Impact of Bit Errors in Digitized RF Data on Ultrasound Image Quality

This article quantitatively analyzes the impact of bit errors in digitized RF data on ultrasound image quality. The quality of B-mode images in both linear array and phased array imaging is evaluated by means of three objective image quality metrics: peak signal-to-noise ratio, structural similarity index, and contrast-to-noise ratio, when bit errors are introduced to the RF data with different bit-error rates (BERs). The effectiveness of coding schemes for forward error detection and correction to improve the image quality is also studied. The results show that ultrasound imaging is inherently resilient to high BER. The image quality suffers unnoticeable degradation for BER lower than 1E-6. Simple 1-bit parity coding with 9% added redundancy helps to retain similar image quality for BER up to 1E-4, and Hamming coding with 33.3% added redundancy allows the BER to increase to 1E-3. These results can serve as a guideline in the datalink design for ultrasound probes with in-probe receive digitization. With much more relaxed BER requirements than in typical datalinks, the design can be optimized by allowing fewer cables with higher data rate per cable or lower power consumption with the same cable count.

Enhanced contrast acoustic-resolution photoacoustic microscopy using double-stage delay-multiply-and-sum beamformer for vasculature imaging

In acoustic-resolution photoacoustic microscopy (AR-PAM) systems, the lateral resolution in the focal zone of the ultrasound (US) transducer is determined by the numerical aperture (NA) of the transducer. To have a high lateral resolution, a large NA is used. However, the larger the NA, the smaller the depth of focus [DOF]. As a result, the lateral resolution is deteriorated at depths out of the focal region. The synthetic aperture focusing technique (SAFT) along with a beamformer can be used to improve the resolution outside the focal region. In this work, for image formation in AR-PAM, we propose the double-stage delay-multiply-and-sum (DS_DMAS) algorithm to be combined with SAFT. The proposed method is evaluated experimentally using hair targets and in vivo vasculature imaging. It is shown that DS_DMAS provides a higher resolution and contrast compared to other methods. For the B-mode images obtained using the hair phantom, the proposed method reduces the average noise level for all the depths by about 134%, 57% and 23%, compared to the original low- resolution, SAFT+DAS and SAFT+DMAS methods, respectively. All the results indicate that the proposed method can be an appropriate algorithm for image formation in AR-PAM systems.

System identification of ankle joint dynamics based on plane-wave ultrasound muscle imaging

Effective treatment of movement disorders requires thorough understanding of human limb control. Joint dynamics can be assessed using robotic manipulators and system identification. Due to tendon compliance, joint angle and muscle length are not proportional. This study uses plane-wave ultrasound imaging to investigate the dynamic relation between ankle joint angle and muscle fiber stretch. The first goal is to determine the feasibility of using ultrasound imaging with system identification; the second goal is to assess the relation between ankle angle, muscle stretch, and reflex size. Soleus and gastrocnemius muscle stretches were assessed with ultrasound imaging and image tracking. For small (1° SD) continuous motions, muscle stretch was proportional to ankle angle during a relax task, but images were too noisy to make that assessment during an active position task. For transient perturbations with high velocity (> 90°/s) the muscle length showed oscillations that were not present in the ankle angle, demonstrating a non-proportional relationship and muscle-tendon interaction. The gastrocnemius velocity predicted the size of the short-latency reflex better than the ankle angle velocity. Concluding, plane-wave ultrasound muscle imaging is feasible for system identification experiments and shows that muscle length and ankle angle are proportional during a relax task with small continuous perturbations.

Two-Stage Beamforming for Phased Array Imaging Using the Fast Hankel Transform

An ultrasound scan generates a huge amount of data. To form an image, this data has to be transferred to the imaging system. This is an issue for applications where the data transfer capacity is limited such as hand-held systems, wireless probes, and miniaturized array probes. Two-stage beamforming methods can be used to significantly reduce the data transfer requirements. In the first stage, which is applied in-probe, the amount of data is reduced from channel to scanline data. In the imaging system, the data are then beamformed to obtain images, which are synthetically focused over the entire image. Currently, two approaches exist for the second stage. The first approach is a time-of-flight (TOF) approach called synthetic aperture sequential beamforming (SASB), which has been developed for both linear and phased arrays. SASB does, however, introduce artifacts in the image that can be reduced by tapering the first-stage scanlines at the cost of lateral resolution. The second approach is based on the wave equation, but a computationally efficient method for phased arrays that is producing sector scan data is lacking. Here, we propose an algorithm that uses the fast Hankel transform to obtain a fast algorithm. The imaging performance of this method is evaluated with simulations and experiments. Compared with PSASB, which is an adaption of SASB for phased arrays, our method requires a similar amount of operations to construct the entire image and there is no tradeoff between resolution and artifacts. These results show the advantage of using the wave equation instead of a TOF approach.

Acoustic Characterization of the CLINIcell for Ultrasound Contrast Agent Studies

Ultrasound contrast agents consist of gas-filled coated microbubbles that oscillate upon ultrasound insonification. Their characteristic oscillatory response provides contrast enhancement for imaging and has the potential to locally enhance drug delivery. Since microbubble response depends on the local acoustic pressure, an ultrasound compatible chamber is needed to study their behavior and the underlying drug delivery pathways. In this study, we determined the amplitude of the acoustic pressure in the CLINIcell, an optically transparent chamber suitable for cell culture. The pressure field was characterized based on microbubble response recorded using the Brandaris 128 ultrahigh-speed camera and an iterative processing method. The results were compared to a control experiment performed in an OptiCell, which is conventionally used in microbubble studies. Microbubbles in the CLINIcell responded in a controlled manner, comparable to those in the OptiCell. For frequencies from 1 to 4 MHz, the mean pressure amplitude was -5.4 dB with respect to the externally applied field. The predictable ultrasound pressure demonstrates the potential of the CLINIcell as an optical, ultrasound, and cell culture compatible device to study microbubble oscillation behavior and ultrasound-mediated drug delivery.

A 2-D Ultrasound Transducer With Front-End ASIC and Low Cable Count for 3-D Forward-Looking Intravascular Imaging: Performance and Characterization

Intravascular ultrasound (IVUS) is an imaging modality used to visualize atherosclerosis from within the inner lumen of human arteries. Complex lesions like chronic total occlusions require forward-looking IVUS (FL-IVUS), instead of the conventional side-looking geometry. Volumetric imaging can be achieved with 2-D array transducers, which present major challenges in reducing cable count and device integration. In this work, we present an 80-element lead zirconium titanate matrix ultrasound transducer for FL-IVUS imaging with a front-end application-specific integrated circuit (ASIC) requiring only four cables. After investigating optimal transducer designs, we fabricated the matrix transducer consisting of 16 transmit (TX) and 64 receive (RX) elements arranged on top of an ASIC having an outer diameter of 1.5 mm and a central hole of 0.5 mm for a guidewire. We modeled the transducer using finite-element analysis and compared the simulation results to the values obtained through acoustic measurements. The TX elements showed uniform behavior with a center frequency of 14 MHz, a -3-dB bandwidth of 44%, and a transmit sensitivity of 0.4 kPa/V at 6 mm. The RX elements showed center frequency and bandwidth similar to the TX elements, with an estimated receive sensitivity of /Pa. We successfully acquired a 3-D FL image of three spherical reflectors in water using delay-and-sum beamforming and the coherence factor method. Full synthetic-aperture acquisition can be achieved with frame rates on the order of 100 Hz. The acoustic characterization and the initial imaging results show the potential of the proposed transducer to achieve 3-D FL-IVUS imaging.

Fast Volumetric Imaging Using a Matrix Transesophageal Echocardiography Probe with Partitioned Transmit-Receive Array

We describe a 3-D multiline parallel beamforming scheme for real-time volumetric ultrasound imaging using a prototype matrix transesophageal echocardiography probe with diagonally diced elements and separated transmit and receive arrays. The elements in the smaller rectangular transmit array are directly wired to the ultrasound system. The elements of the larger square receive aperture are grouped in 4 × 4-element sub-arrays by micro-beamforming in an application-specific integrated circuit. We propose a beamforming sequence with 85 transmit-receive events that exhibits good performance for a volume sector of 60° × 60°. The beamforming is validated using Field II simulations, phantom measurements and in vivo imaging. The proposed parallel beamforming achieves volume rates up to 59 Hz and produces good-quality images by angle-weighted combination of overlapping sub-volumes. Point spread function, contrast ratio and contrast-to-noise ratio in the phantom experiment closely match those of the simulation. In vivo 3-D imaging at 22-Hz volume rate in a healthy adult pig clearly visualized the cardiac structures, including valve motion.

Acoustic Characterization of a Vessel-on-a-Chip Microfluidic System for Ultrasound-Mediated Drug Delivery

Ultrasound in the presence of gas-filled microbubbles can be used to enhance local uptake of drugs and genes. To study the drug delivery potential and its underlying physical and biological mechanisms, an in vitro vessel model should ideally include 3-D cell culture, perfusion flow, and membrane-free soft boundaries. Here, we propose an organ-on-a-chip microfluidic platform to study ultrasound-mediated drug delivery: the OrganoPlate. The acoustic propagation into the OrganoPlate was determined to assess the feasibility of controlled microbubble actuation, which is required to study the microbubble-cell interaction for drug delivery. The pressure field in the OrganoPlate was characterized non-invasively by studying experimentally the well-known response of microbubbles and by simulating the acoustic wave propagation in the system. Microbubble dynamics in the OrganoPlate were recorded with the Brandaris 128 ultrahigh-speed camera (17 million frames/s) and a control experiment was performed in an OptiCell, an in vitro monolayer cell culture chamber that is conventionally used to study ultrasound-mediated drug delivery. When insonified at frequencies between 1 and 2 MHz, microbubbles in the OrganoPlate experienced larger oscillation amplitudes resulting from higher local pressures. Microbubbles responded similarly in both systems when insonified at frequencies between 2 and 4 MHz. Numerical simulations performed with a 3-D finite-element model of ultrasound propagation into the OrganoPlate and the OptiCell showed the same frequency-dependent behavior. The predictable and homogeneous pressure field in the OrganoPlate demonstrates its potential to develop an in vitro 3-D cell culture model, well suited to study ultrasound-mediated drug delivery.

Dual stage beamforming in the absence of front-end receive focusing

Ultrasound front-end receive designs for miniature, wireless, and/or matrix transducers can be simplified considerably by direct-element summation in receive. In this paper we develop a dual-stage beamforming technique that is able to produce a high-quality image from scanlines that are produced with focused transmit, and simple summation in receive (no delays). We call this non-delayed sequential beamforming (NDSB). In the first stage, low-resolution RF scanlines are formed by simple summation of element signals from a running sub-aperture. In the second stage, delay-and-sum beamforming is performed in which the delays are calculated considering the transmit focal points as virtual sources emitting spherical waves, and the sub-apertures as large unfocused receive elements. The NDSB method is validated with simulations in Field II. For experimental validation, RF channel data were acquired with a commercial research scanner using a 5 MHz linear array, and were subsequently processed offline. For NDSB, good average lateral resolution (0.99 mm) and low grating lobe levels (<-40 dB) were achieved by choosing the transmit [Formula: see text] as 0.75 and the transmit focus at 15 mm. NDSB was compared with conventional dynamic receive focusing (DRF) and synthetic aperture sequential beamforming (SASB) with their own respective optimal settings. The full width at half maximum of the NDSB point spread function was on average 20% smaller than that of DRF except for at depths  <30 mm and 10% larger than SASB considering all the depths. NDSB showed only a minor degradation in contrast-to-noise ratio and contrast ratio compared to DRF and SASB when measured on an anechoic cyst embedded in a tissue-mimicking phantom. In conclusion, using simple receive electronics front-end, NDSB can attain an image quality better than DRF and slightly inferior to SASB.

Improving the Performance of a 1-D Ultrasound Transducer Array by Subdicing

In medical ultrasound transducer design, the geometry of the individual elements is crucial since it affects the vibration mode of each element and its radiation impedance. For a fixed frequency, optimal vibration (i.e., uniform surface motion) can be achieved by designing elements with very small width-to-thickness ratios. However, for optimal radiation impedance (i.e., highest radiated power), the width should be as large as possible. This leads to a contradiction that can be solved by subdicing wide elements. To systematically examine the effect of subdicing on the performance of a 1-D ultrasound transducer array, we applied finite-element simulations. We investigated the influence of subdicing on the radiation impedance, on the time and frequency response, and on the directivity of linear arrays with variable element widths. We also studied the effect of varying the depth of the subdicing cut. The results show that, for elements having a width greater than 0.6 times the wavelength, subdicing improves the performance compared with that of nonsubdiced elements: the emitted pressure may be increased up to a factor of three, the ringing time may be reduced by up to 50%, the bandwidth increased by up to 77%, and the sidelobes reduced by up to 13 dB. Moreover, this simulation study shows that all these improvements can already be achieved by subdicing the elements to a depth of 70% of the total element thickness. Thus, subdicing can improve important transducer parameters and, therefore, help in achieving images with improved signal-to-noise ratio and improved resolution.

A Broadband Polyvinylidene Difluoride-Based Hydrophone with Integrated Readout Circuit for Intravascular Photoacoustic Imaging

Intravascular photoacoustic (IVPA) imaging can visualize the coronary atherosclerotic plaque composition on the basis of the optical absorption contrast. Most of the photoacoustic (PA) energy of human coronary plaque lipids was found to lie in the frequency band between 2 and 15 MHz requiring a very broadband transducer, especially if a combination with intravascular ultrasound is desired. We have developed a broadband polyvinylidene difluoride (PVDF) transducer (0.6 × 0.6 mm, 52 μm thick) with integrated electronics to match the low capacitance of such a small polyvinylidene difluoride element (<5 pF/mm(2)) with the high capacitive load of the long cable (∼100 pF/m). The new readout circuit provides an output voltage with a sensitivity of about 3.8 μV/Pa at 2.25 MHz. Its response is flat within 10 dB in the range 2 to 15 MHz. The root mean square (rms) output noise level is 259 μV over the entire bandwidth (1-20 MHz), resulting in a minimum detectable pressure of 30 Pa at 2.25 MHz.

A Prototype PZT Matrix Transducer With Low-Power Integrated Receive ASIC for 3-D Transesophageal Echocardiography

This paper presents the design, fabrication, and experimental evaluation of a prototype lead zirconium titanate (PZT) matrix transducer with an integrated receive ASIC, as a proof of concept for a miniature three-dimensional (3-D) transesophageal echocardiography (TEE) probe. It consists of an array of 9 ×12 piezoelectric elements mounted on the ASIC via an integration scheme that involves direct electrical connections between a bond-pad array on the ASIC and the transducer elements. The ASIC addresses the critical challenge of reducing cable count, and includes front-end amplifiers with adjustable gains and micro-beamformer circuits that locally process and combine echo signals received by the elements of each 3 ×3 subarray. Thus, an order-of-magnitude reduction in the number of receive channels is achieved. Dedicated circuit techniques are employed to meet the strict space and power constraints of TEE probes. The ASIC has been fabricated in a standard 0.18-μm CMOS process and consumes only 0.44 mW/channel. The prototype has been acoustically characterized in a water tank. The ASIC allows the array to be presteered across ±37° while achieving an overall dynamic range of 77 dB. Both the measured characteristics of the individual transducer elements and the performance of the ASIC are in good agreement with expectations, demonstrating the effectiveness of the proposed techniques.

F-k Domain Imaging for Synthetic Aperture Sequential Beamforming

Spatial resolution in medical ultrasound images is a key component in image quality and an important factor for clinical diagnosis. In early systems, the lateral resolution was optimal in the focus but rapidly decreased outside the focal region. Improvements have been found in, e.g., dynamic-receive beamforming, in which the entire image is focused in receive, but this requires complex processing of element data and is not applicable for mechanical scanning of single-element images. This paper exploits the concept of two-stage beamforming based on virtual source-receivers, which reduces the front-end computational load while maintaining a similar data rate and frame rate compared to dynamic-receive beamforming. We introduce frequency-wavenumber domain data processing to obtain fast second-stage data processing while having similarly high lateral resolution as dynamic-receive beamforming and processing in time-space domain. The technique is very suitable in combination with emerging technologies such as application-specific integrated circuits (ASICs), hand-held devices, and wireless data transfer. The suggested method consists of three steps. In the first step, single-focused RF line data are shifted in time to relocate the focal point to a new origin t' = 0, z' = 0. This new origin is considered as an array of virtual source/receiver pairs, as has been suggested previously in literature. In the second step, the dataset is efficiently processed in the wavenumber-frequency domain to form an image that is in focus throughout its entire depth. In the third step, the data shift is undone to obtain a correct depth axis in the image. The method has been tested first with a single-element scanning system and second in a tissue-mimicking phantom using a linear array. In both setups, the method resulted in a −6-dB lateral point spread function (PSF) which was constant over the entire depth range, and similar to dynamic-receive beamforming and synthetic aperture sequential beamforming. The signal-to-noise ratio increased by 6 dB in both the near field and far field. These results show that the second-stage processing algorithm effectively produces a focused image over the entire depth range from a single-focused ultrasound field.

Calibrating Doppler imaging of preterm intracerebral circulation using a microvessel flow phantom

Introduction: Preterm infants are born during critical stages of brain development, in which the adaptive capacity of the fetus to extra-uterine environment is limited. Inadequate brain perfusion has been directly linked to preterm brain damage. Advanced high-frequency ultrasound probes and processing algorithms allow visualization of microvessels and depiction of regional variation. To assess whether visualization and flow velocity estimates of preterm cerebral perfusion using Doppler techniques are accurate, we conducted an in vitro experiment using a microvessel flow phantom.

Materials and Methods: An in-house developed flow phantom containing two microvessels (inner diameter 200 and 700 μm) with attached syringe pumps, filled with blood-mimicking fluid, was used to generate non-pulsatile perfusion of variable flow. Measurements were performed using an Esaote MyLab70 scanner.

Results: Microvessel mimicking catheters with velocities as low as 1 cm/s were adequately visualized with a linear ultrasound probe. With a convex probe, velocities <2 cm/s could not be depicted. Within settings, velocity and diameter measurements were highly reproducible [intra-class correlation 0.997 (95% CI 0.996–0.998) and 0.914 (0.864–0.946)]. Overall, mean velocity was overestimated up to threefold, especially in high velocity ranges. Significant differences were seen in velocity measurements when using steer angle correction and in vessel diameter estimation (p < 0.05).

Conclusion: Visualization of microvessel-size catheters mimicking small brain vessels is feasible. Reproducible velocity and diameter results can be obtained, although important overestimation of the values is observed. Before velocity estimates of microcirculation can find its use in clinical practice, calibration of the ultrasound machine for any specific Doppler purpose is essential. The ultimate goal is to develop a sonographic tool that can be used for objective study of regional perfusion in routine practice.

Computation of nonlinear ultrasound fields using a linearized contrast source method

Nonlinear ultrasound is important in medical diagnostics because imaging of the higher harmonics improves resolution and reduces scattering artifacts. Second harmonic imaging is currently standard, and higher harmonic imaging is under investigation. The efficient development of novel imaging modalities and equipment requires accurate simulations of nonlinear wave fields in large volumes of realistic (lossy, inhomogeneous) media. The Iterative Nonlinear Contrast Source (INCS) method has been developed to deal with spatiotemporal domains measuring hundreds of wavelengths and periods. This full wave method considers the nonlinear term of the Westervelt equation as a nonlinear contrast source, and solves the equivalent integral equation via the Neumann iterative solution. Recently, the method has been extended with a contrast source that accounts for spatially varying attenuation. The current paper addresses the problem that the Neumann iterative solution converges badly for strong contrast sources. The remedy is linearization of the nonlinear contrast source, combined with application of more advanced methods for solving the resulting integral equation. Numerical results show that linearization in combination with a Bi-Conjugate Gradient Stabilized method allows the INCS method to deal with fairly strong, inhomogeneous attenuation, while the error due to the linearization can be eliminated by restarting the iterative scheme.

Single pulse frequency compounding protocol for superharmonic imaging

Second harmonic imaging is currently accepted as the standard in commercial echographic systems. A new imaging technique, coined as superharmonic imaging (SHI), combines the third till the fifth harmonics, arising during nonlinear sound propagation. It could further enhance the resolution and quality of echographic images. To meet the bandwidth requirement for SHI a dedicated phased array has been developed: a low frequency subarray, intended for transmission, interleaved with a high frequency subarray, used in reception. As the bandwidth of the elements is limited, the spectral gaps in between the harmonics cause multiple (ghost) reflection artifacts. A dual-pulse frequency compounding method aims at suppressing those artifacts at a price of a reduced frame rate. In this study we explore a possibility of performing frequency compounding within a single transmission. The traditional frequency compounding method suppresses the ripples by consecutively emitting two short Gaussian bursts with a slightly different center frequency. In the newly proposed method, the transmit aperture is divided into two parts: the first half is used to send a pulse at the lower center frequency, while the other half simultaneously transmits at a slightly higher center frequency. The suitability of the protocol for medical imaging applications in terms of the steering capabilities was performed in a simulation study with INCS and the hydrophone measurements. Moreover, an experimental study was carried out to find the optimal parameters for the clinical imaging protocol. The latter was subsequently used to obtain the images of a tissue mimicking phantom containing strongly reflecting wires. Additionally, the images of a human heart in the parasternal projection were acquired. The scanning aperture with the developed protocol amounts to approximately 90°, which is sufficient to capture the cardiac structures in the standard anatomical projections. The theoretically estimated and experimentally measured grating lobe levels are equal to -28.3 dB and -35.9 dB, respectively. A considerable improvement in the axial resolution of the SHI component (0.73 mm) at -6 dB in comparison with the third harmonic (2.23 mm) was observed. A similar comparison in terms of the lateral resolution slightly favored the superharmonic component by 0.2 mm. Additionally, the images of the tissue mimicking phantom exhibited the absence of the multiple reflection artifacts. The in-vivo acquisition allows one to clearly observe the dynamic of the mitral valve leaflets. The new method is equally effective in eliminating the ripple artifacts associated with SHI as the dual-pulse technique, while the full frame rate is maintained.

Comparison of fundamental, second harmonic, and superharmonic imaging: a simulation study

In medical ultrasound, fundamental imaging (FI) uses the reflected echoes from the same spectral band as that of the emitted pulse. The transmission frequency determines the trade-off between penetration depth and spatial resolution. Tissue harmonic imaging (THI) employs the second harmonic of the emitted frequency band to construct images. Recently, superharmonic imaging (SHI) has been introduced, which uses the third to the fifth (super) harmonics. The harmonic level is determined by two competing phenomena: nonlinear propagation and frequency dependent attenuation. Thus, the transmission frequency yielding the optimal trade-off between the spatial resolution and the penetration depth differs for THI and SHI. This paper quantitatively compares the concepts of fundamental, second harmonic, and superharmonic echocardiography at their optimal transmission frequencies. Forward propagation is modeled using a 3D-KZK implementation and the iterative nonlinear contrast source (INCS) method. Backpropagation is assumed to be linear. Results show that the fundamental lateral beamwidth is the narrowest at focus, while the superharmonic one is narrower outside the focus. The lateral superharmonic roll-off exceeds the fundamental and second harmonic roll-off. Also, the axial resolution of SHI exceeds that of FI and THI. The far-field pulse-echo superharmonic pressure is lower than that of the fundamental and second harmonic. SHI appears suited for echocardiography and is expected to improve its image quality at the cost of a slight reduction in depth-of-field.

A contrast source method for nonlinear acoustic wave fields in media with spatially inhomogeneous attenuation

Experimental data reveals that attenuation is an important phenomenon in medical ultrasound. Attenuation is particularly important for medical applications based on nonlinear acoustics, since higher harmonics experience higher attenuation than the fundamental. Here, a method is presented to accurately solve the wave equation for nonlinear acoustic media with spatially inhomogeneous attenuation. Losses are modeled by a spatially dependent compliance relaxation function, which is included in the Westervelt equation. Introduction of absorption in the form of a causal relaxation function automatically results in the appearance of dispersion. The appearance of inhomogeneities implies the presence of a spatially inhomogeneous contrast source in the presented full-wave method leading to inclusion of forward and backward scattering. The contrast source problem is solved iteratively using a Neumann scheme, similar to the iterative nonlinear contrast source (INCS) method. The presented method is directionally independent and capable of dealing with weakly to moderately nonlinear, large scale, three-dimensional wave fields occurring in diagnostic ultrasound. Convergence of the method has been investigated and results for homogeneous, lossy, linear media show full agreement with the exact results. Moreover, the performance of the method is demonstrated through simulations involving steered and unsteered beams in nonlinear media with spatially homogeneous and inhomogeneous attenuation.

Optimization of a phased-array transducer for multiple harmonic imaging in medical applications: frequency and topology

Second-harmonic imaging is currently one of the standards in commercial echographic systems for diagnosis, because of its high spatial resolution and low sensitivity to clutter and near-field artifacts. The use of nonlinear phenomena mirrors is a great set of solutions to improve echographic image resolution. To further enhance the resolution and image quality, the combination of the 3rd to 5th harmonics--dubbed the superharmonics--could be used. However, this requires a bandwidth exceeding that of conventional transducers. A promising solution features a phased-array design with interleaved low- and high-frequency elements for transmission and reception, respectively. Because the amplitude of the backscattered higher harmonics at the transducer surface is relatively low, it is highly desirable to increase the sensitivity in reception. Therefore, we investigated the optimization of the number of elements in the receiving aperture as well as their arrangement (topology). A variety of configurations was considered, including one transmit element for each receive element (1/2) up to one transmit for 7 receive elements (1/8). The topologies are assessed based on the ratio of the harmonic peak pressures in the main and grating lobes. Further, the higher harmonic level is maximized by optimization of the center frequency of the transmitted pulse. The achievable SNR for a specific application is a compromise between the frequency-dependent attenuation and nonlinearity at a required penetration depth. To calculate the SNR of the complete imaging chain, we use an approach analogous to the sonar equation used in underwater acoustics. The generated harmonic pressure fields caused by nonlinear wave propagation were modeled with the iterative nonlinear contrast source (INCS) method, the KZK, or the Burger's equation. The optimal topology for superharmonic imaging was an interleaved design with 1 transmit element per 6 receive elements. It improves the SNR by ~5 dB compared with the interleaved (1/2) design reported in literature. The optimal transmit frequency for superharmonic echocardiography was found to be 1.0 to 1.2 MHz. For superharmonic abdominal imaging this frequency was found to be 1.7 to 1.9 MHz. For 2nd-harmonic echocardiography, the optimal transmit frequency of 1.8 MHz reported in the literature was corroborated with our simulation results.