Showing papers in "IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control in 2017"

Acoustic Wave Filter Technology–A Review

Abstract: Today, acoustic filters are the filter technology to meet the requirements with respect to performance dictated by the cellular phone standards and their form factor. Around two billion cellular phones are sold every year, and smart phones are of a very high percentage of approximately two-thirds. Smart phones require a very high number of filter functions ranging from the low double-digit range up to almost triple digit numbers in the near future. In the frequency range up to 1 GHz, surface acoustic wave (SAW) filters are almost exclusively employed, while in the higher frequency range, bulk acoustic wave (BAW) and SAW filters are competing for their shares. Prerequisites for the success of acoustic filters were the availability of high-quality substrates, advanced and highly reproducible fabrication technologies, optimum filter techniques, precise simulation software, and advanced design tools that allow the fast and efficient design according to customer specifications. This paper will try to focus on innovations leading to high volume applications of intermediate frequency (IF) and radio frequency (RF) acoustic filters, e.g., TV IF filters, IF filters for cellular phones, and SAW/BAW RF filters for the RF front-end of cellular phones.

High-Performance SAW Resonator on New Multilayered Substrate Using LiTaO 3 Crystal

Abstract: To develop the high-performance filters and duplexers required for recent long-term evolution frequency bands in mobile handsets, a surface acoustic wave (SAW) resonator is needed that has a higher quality (Q) and a lower temperature coefficient of frequency (TCF). To achieve this, the authors focused on acoustic energy confinement in the depth direction for a rotated Y-X LiTaO3 (LT) substrate. Characteristics of multilayered substrates with low-impedance and high-impedance layers under LT layer were studied numerically in terms of acoustic energy distribution, phase velocity, coupling coefficient, and temperature characteristics employing a finite-element method simulation. After several calculations, a novel multilayered structure was developed that uses SiO2 for a low-impedance layer and AlN for a high-impedance layer under the thin LT layer. A one-port resonator using the new substrate was fabricated, and its experimental results showed that the developed resonator had a Bode-Q over 4000 and TCF of −8 ppm/°C, which are four times higher than and one-fifth as small as those of a conventional 4° YX-LT SAW resonator, respectively. By applying this technology, a band 25 duplexer with very narrow duplex gap was successfully developed, which shows extremely low insertion loss, steep cutoff characteristics, and stable temperature characteristics.

Quantitative Frequency-Domain Passive Cavitation Imaging

Abstract: Passive cavitation detection has been an instrumental technique for measuring cavitation dynamics, elucidating concomitant bioeffects, and guiding ultrasound therapies. Recently, techniques have been developed to create images of cavitation activity to provide investigators with a more complete set of information. These techniques use arrays to record and subsequently beamform received cavitation emissions, rather than processing emissions received on a single-element transducer. In this paper, the methods for performing frequency-domain delay, sum, and integrate passive imaging are outlined. The method can be applied to any passively acquired acoustic scattering or emissions, including cavitation emissions. To compare data across different systems, techniques for normalizing Fourier transformed data and converting the data to the acoustic energy received by the array are described. A discussion of hardware requirements and alternative imaging approaches is additionally outlined. Examples are provided in MATLAB.

High-Quality Plane Wave Compounding Using Convolutional Neural Networks

Abstract: Single plane wave (PW) imaging produces ultrasound images of poor quality at high frame rates (ultrafast). High-quality PW imaging usually relies on the coherent compounding of several successive steered emissions (typically more than ten), which in turn results in a decreased frame rate. We propose a new strategy to reduce the number of emitted PWs by learning a compounding operation from data, i.e., by training a convolutional neural network to reconstruct high-quality images using a small number of transmissions. We present experimental evidence that this approach is promising, as we were able to produce high-quality images from only three PWs, competing in terms of contrast ratio and lateral resolution with the standard compounding of 31 PWs ( $10\times $ speedup factor).

Global Time-Delay Estimation in Ultrasound Elastography

Abstract: A critical step in quasi-static ultrasound elastography is the estimation of time delay between two frames of radio-frequency (RF) data that are obtained while the tissue is undergoing deformation. This paper presents a novel technique for time-delay estimation (TDE) of all samples of RF data simultaneously , thereby exploiting all the information in RF data for TDE. A nonlinear cost function that incorporates similarity of RF data intensity and prior information of displacement continuity is formulated. Optimization of this function involves searching for TDE of all samples of the RF data, rendering the optimization intractable with conventional techniques given that the number of variables can be approximately one million. Therefore, the optimization problem is converted to a sparse linear system of equations, and is solved in real time using a computationally efficient optimization technique. We call our method GLobal Ultrasound Elastography (GLUE), and compare it to dynamic programming analytic minimization (DPAM) and normalized cross correlation (NCC) techniques. Our simulation results show that the contrast-to-noise ratio (CNR) values of the axial strain maps are 4.94 for NCC, 14.62 for DPAM, and 26.31 for GLUE. Our results on experimental data from tissue mimicking phantoms show that the CNR values of the axial strain maps are 1.07 for NCC, 16.01 for DPAM, and 18.21 for GLUE. Finally, our results on in vivo data show that the CNR values of the axial strain maps are 3.56 for DPAM and 13.20 for GLUE.

Ultrasonic Fingerprint Sensor With Transmit Beamforming Based on a PMUT Array Bonded to CMOS Circuitry

Abstract: In this paper, we present a single-chip $65 \times 42$ element ultrasonic pulse-echo fingerprint sensor with transmit (TX) beamforming based on piezoelectric micromachined ultrasonic transducers directly bonded to a CMOS readout application-specific integrated circuit (ASIC) The readout ASIC was realized in a standard 180-nm CMOS process with a 24-V high-voltage transistor option Pulse-echo measurements are performed column-by-column in sequence using either one column or five columns to TX the ultrasonic pulse at 20 MHz TX beamforming is used to focus the ultrasonic beam at the imaging plane where the finger is located, increasing the ultrasonic pressure and narrowing the 3-dB beamwidth to $50~\mu \text{m}$ , a factor of 64 narrower than nonbeamformed measurements The surface of the sensor is coated with a poly-dimethylsiloxane (PDMS) layer to provide good acoustic impedance matching to skin Scanning laser Doppler vibrometry of the PDMS surface was used to map the ultrasonic pressure field at the imaging surface, demonstrating the expected increase in pressure, and reduction in beamwidth Imaging experiments were conducted using both PDMS phantoms and real fingerprints The average image contrast is increased by a factor of 15 when beamforming is used

Ultrasonic Shear Wave Elasticity Imaging Sequencing and Data Processing Using a Verasonics Research Scanner

Abstract: Ultrasound elasticity imaging has been developed over the last decade to estimate tissue stiffness. Shear wave elasticity imaging (SWEI) quantifies tissue stiffness by measuring the speed of propagating shear waves following acoustic radiation force excitation. This paper presents the sequencing and data processing protocols of SWEI using a Verasonics system. The selection of the sequence parameters in a Verasonics programming script is discussed in detail. The data processing pipeline to calculate group shear wave speed (SWS), including tissue motion estimation, data filtering, and SWS estimation, is demonstrated. In addition, the procedures for calibration of beam position, scanner timing, and transducer face heating are provided to avoid SWS measurement bias and transducer damage.

Architecture of an Ultrasound System for Continuous Real-Time High Frame Rate Imaging

Abstract: High frame rate (HFR) imaging methods based on the transmission of defocused or plane waves rather than focused beams are increasingly popular. However, the production of HFR images poses severe requirements both in the transmission and the reception sections of ultrasound scanners. In particular, major technical difficulties arise if the images must be continuously produced in real-time, i.e., without any acquisition interruption nor loss of data. This paper presents the implementation of the real-time HFR-compounded imaging application in the ULA-OP 256 research platform. The beamformer sustains an average output sample rate of 470 MSPS. This allows continuously producing coherently compounded images, each of 64 lines by 1280 depths (here corresponding to 15.7 mm width and 45 mm depth, respectively), at frame rates up to 5.3 kHz. Imaging tests addressed to evaluate the achievable speed and quality performance were conducted on phantom. Results obtained by real-time compounding frames obtained with different numbers of steering angles between +7.5° and −7.5° are presented.

Microbubble Axial Localization Errors in Ultrasound Super-Resolution Imaging

Abstract: Acoustic super-resolution imaging has allowed the visualization of microvascular structure and flow beyond the diffraction limit using standard clinical ultrasound systems through the localization of many spatially isolated microbubble signals. The determination of each microbubble position is typically performed by calculating the centroid, finding a local maximum, or finding the peak of a 2-D Gaussian function fit to the signal. However, the backscattered signal from a microbubble depends not only on diffraction characteristics of the waveform, but also on the microbubble behavior in the acoustic field. Here, we propose a new axial localization method by identifying the onset of the backscattered signal. We compare the accuracy of localization methods using in vitro experiments performed at 7-cm depth and 2.3-MHz center frequency. We corroborate these findings with simulation results based on the Marmottant model. We show experimentally and in simulations that detecting the onset of the returning signal provides considerably increased accuracy for super-resolution. Resulting experimental cross-sectional profiles in super-resolution images demonstrate at least 5.8 times improvement in contrast ratio and more than 1.8 times reduction in spatial spread (provided by 90% of the localizations) for the onset method over centroiding, peak detection, and 2-D Gaussian fitting methods. Simulations estimate that these latter methods could create errors in relative bubble positions as high as $900~\mu \text{m}$ at these experimental settings, while the onset method reduced the interquartile range of these errors by a factor of over 2.2. Detecting the signal onset is, therefore, expected to considerably improve the accuracy of super-resolution.

Accelerated Singular Value-Based Ultrasound Blood Flow Clutter Filtering With Randomized Singular Value Decomposition and Randomized Spatial Downsampling

Abstract: Singular value decomposition (SVD)-based ultrasound blood flow clutter filters have recently demonstrated substantial improvement in clutter rejection for ultrafast plane wave microvessel imaging, and have become the commonly used clutter filtering method for many novel ultrafast imaging applications such as functional ultrasound and super-resolution imaging. At present, however, the computational burden of SVD remains as a major hurdle for practical implementation and clinical translation of this method. To address this challenge, in the study we present two blood flow clutter filtering methods based on randomized SVD (rSVD) and randomized spatial downsampling to accelerate SVD clutter filtering with minimal compromise to the clutter filter performance. rSVD accelerates SVD computation by approximating the $k$ largest singular values, while random downsampling accelerates both full SVD and rSVD by decomposing the original large data matrix into small matrices that can be processed in parallel. An in vitro blood flow phantom study with the presence of heavy tissue clutter showed significantly improved computational performance using the proposed methods with minimal deterioration to the clutter filter performance (less than 3-dB reduction in blood to clutter ratio, less than 0.2-cm2/s2 increase in flow mean squared error, less than 0.1-cm/s increase in the standard deviation of the vessel blood flow signal, and less than 0.3-cm/s increase in tissue clutter velocity for both full SVD and rSVD when the downsampling factor was less than $20\times$ ). The maximum acceleration was about threefold from randomized spatial downsampling, and approximately another threefold from rSVD. An in vivo rabbit kidney perfusion study showed that rSVD provided comparable performance to full SVD in clutter rejection in vivo (maximum difference of blood to clutter ratio was less than 0.6 dB), and random downsampling provided artifact-free perfusion imaging results when combined with both full SVD and rSVD. The blood to clutter ratio was still above 10 dB with a downsampling factor of $60\times$ . We also demonstrated real-time microvessel imaging feasibility (~40-ms processing time) by combining rSVD with random downsampling. The findings and methods presented in this paper may greatly facilitate the new area of ultrafast microvessel imaging research.

3-D Nonlinear Acoustic Inverse Scattering: Algorithm and Quantitative Results

Abstract: We describe a novel 3-D ultrasound technology, the quantitative transmission ultrasound system and algorithm to image a pendent breast in a water bath. Quantitative accuracy is verified using phantoms. Morphological accuracy is verified using cadaveric breast and in vivo images, and spatial resolution is estimated. This paper generalizes an earlier 2-D algorithm to a full 3-D inversion algorithm and shows the importance of such a 3-D algorithm for artifact suppression as compared with the 2-D algorithm. The resultant high-resolution ultrasound images, along with quantitative information regarding tissue speed of sound/stiffness, provide a more accurate depiction of the breast anatomy and lesions, contributing to improved breast care.

Methods of Generating Submicrometer Phase-Shift Perfluorocarbon Droplets for Applications in Medical Ultrasonography

Abstract: Continued advances in the field of ultrasound and ultrasound contrast agents have created new approaches to imaging and medical intervention. Phase-shift perfluorocarbon droplets, which can be vaporized by ultrasound energy to transition from the liquid to the vapor state, are one of the most highly researched alternatives to clinical ultrasound contrast agents (i.e., microbubbles). In this paper, part of a special issue on methods in biomedical ultrasonics, we survey current techniques to prepare ultrasound-activated nanoscale phase-shift perfluorocarbon droplets, including sonication, extrusion, homogenization, microfluidics, and microbubble condensation. We provide example protocols and discuss advantages and limitations of each approach. Finally, we discuss best practice in characterization of this class of contrast agents with respect to size distribution and ultrasound activation.

Lung Ultrasound Surface Wave Elastography: A Pilot Clinical Study

Abstract: A lung ultrasound surface wave elastography (LUSWE) technique is developed to measure superficial lung tissue elastic properties. The purpose of this paper was to translate LUSWE into clinical studies for assessing patients with interstitial lung disease (ILD) and present the pilot data from lung measurements on 10 healthy subjects and 10 patients with ILD. ILD includes multiple lung disorders in which the lung tissue is distorted and stiffened by tissue fibrosis. Chest radiography and computed tomography are the most commonly used techniques for assessing lung disease, but they are associated with radiation and cannot directly measure lung elastic properties. LUSWE provides a noninvasive and nonionizing technique to measure the elastic properties of superficial lung tissue. LUSWE was used to measure regions of both lungs through six intercostal spaces for patients and healthy subjects. The data are presented as wave speed at 100, 150, and 200 Hz at the six intercostal spaces. As an example, the surface wave speeds are, respectively, 1.88 ± 0.11 m/s at 100 Hz, 2.74 ± 0.26 m/s at 150 Hz, and 3.62 ± 0.13 m/s at 200 Hz for a healthy subject in the upper right lung; this is in comparison to measurements from an ILD patient of 3.3 ± 0.37 m/s at 100 Hz, 4.38 ± 0.33 m/s at 150 Hz, and 5.24 ± 0.44 m/s at 200 Hz in the same lung space. Significant differences in wave speed between healthy subjects and ILD patients were found. LUSWE is a safe and noninvasive technique which may be useful for assessing ILD.

Design of HIFU Transducers for Generating Specified Nonlinear Ultrasound Fields

Abstract: Various clinical applications of high-intensity focused ultrasound have different requirements for the pressure levels and degree of nonlinear waveform distortion at the focus. The goal of this paper is to determine transducer design parameters that produce either a specified shock amplitude in the focal waveform or specified peak pressures while still maintaining quasi-linear conditions at the focus. Multiparametric nonlinear modeling based on the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation with an equivalent source boundary condition was employed. Peak pressures, shock amplitudes at the focus, and corresponding source outputs were determined for different transducer geometries and levels of nonlinear distortion. The results are presented in terms of the parameters of an equivalent single-element spherically shaped transducer. The accuracy of the method and its applicability to cases of strongly focused transducers were validated by comparing the KZK modeling data with measurements and nonlinear full diffraction simulations for a single-element source and arrays with 7 and 256 elements. The results provide look-up data for evaluating nonlinear distortions at the focus of existing therapeutic systems as well as for guiding the design of new transducers that generate specified nonlinear fields.

Erratum to “A 200–1380 kHz Quadrifrequency Focused Ultrasound Transducer for Neurostimulation in Rodents and Primates: Transcranial In Vitro Calibration and Numerical Study of the Influence of Skull Cavity”

Abstract: Low intensity transcranial focused ultrasound has been demonstrated to produce neuromodulation in both animals and humans. Primarily for technical reasons, frequency is one of the most poorly investigated critical wave parameters. We propose the use of a quadri-band transducer capable of operating at 200, 320, 850, and 1380 kHz for further investigation of the frequency dependence of neuromodulation efficacy while keeping the position of the transducer fixed with respect to the subject’s head. This paper presents the results of the transducer calibration in water, in vitro transmission measurements through a monkey skull flap, 3-D simulations based on both a $\mu $ -computed tomography ( $\mu $ CT)-scan of a rat and on CT-scans of two macaques. A maximum peak pressure greater than 0.52 MPa is expected at each frequency in rat and macaque heads. According to the literature, our transducer can achieve neuromodulation in rodents and primates at each four frequencies. The impact of standing waves is shown to be most prominent at the lowest frequencies.

Guidelines for Finite-Element Modeling of Acoustic Radiation Force-Induced Shear Wave Propagation in Tissue-Mimicking Media

Abstract: Ultrasound shear wave elastography is emerging as an important imaging modality for evaluating tissue material properties. In its practice, some systematic biases have been associated with ultrasound frequencies, focal depths and configuration, and transducer types (linear versus curvilinear), along with displacement estimation and shear wave speed estimation algorithms. Added to that, soft tissues are not purely elastic, so shear waves will travel at different speeds depending on their spectral content, which can be modulated by the acoustic radiation force (ARF) excitation focusing, duration, and the frequency-dependent stiffness of the tissue. To understand how these different acquisition and material property parameters may affect the measurements of shear wave velocity, the simulations of the propagation of shear waves generated by ARF excitations in viscoelastic media are a very important tool. This paper serves to provide an in-depth description of how these simulations are performed. The general scheme is broken into three components: 1) simulation of the 3-D ARF push beam; 2) applying that force distribution to a finite-element model; and 3) extraction of the motion data for post-processing. All three components will be described in detail and combined to create a simulation platform that is powerful for developing and testing algorithms for academic and industrial researchers involved in making quantitative shear-wave-based measurements of tissue material properties.

The Practicalities of Obtaining and Using Hydrophone Calibration Data to Derive Pressure Waveforms

Abstract: This paper considers the means by which calibration data are used to convert hydrophone output voltage into pressure. Hydrophone frequency responses are complex-valued quantities, and only by correcting for the magnitude and phase variations, is it possible to accurately recover the original pressure waveform. The limitations of current hydrophone calibration techniques are discussed, and a new method of obtaining hydrophone phase data is presented. Magnitude and phase information is measured via both coarse increment (1 MHz) and fine increment (50 kHz) calibration techniques for three exemplar hydrophones (0.5 mm needle, 0.2 mm needle, and 0.4 mm membrane). Frequently hydrophone calibration data are available at frequency increments that do not match that required by the deconvolution process. Therefore, a variety of methods to interpolate the calibrated system response to obtain correctly spaced data are considered, and two spline interpolation methods are found to offer best performance. Data preconditioning and filtering to address artifacts above and below the 1 to 40 MHz bandwidth of the coarse frequency increment calibration are also investigated, and a simple procedure for selecting an appropriate low-pass filter is presented. The revised calibration data are used to deconvolve the hydrophone frequency response for experimentally derived waveforms. Standard ultrasonic output parameters (such as peak compressional and peak rarefactional pressures, pulse intensity integral, and temporal peak and pulse average acoustic intensities) are calculated from these waveforms. Although the three hydrophones used in this paper are of different types and have a range of active element sizes, all output parameters derived from the deconvolved waveforms have <5% variation from their respective population means (with the majority being within <2%).

Efficient Strategies for Estimating the Spatial Coherence of Backscatter

Abstract: The spatial coherence of ultrasound backscatter has been proposed to reduce clutter in medical imaging, to measure the anisotropy of the scattering source, and to improve the detection of blood flow. These techniques rely on correlation estimates that are obtained using computationally expensive strategies. In this paper, we assess the existing spatial coherence estimation methods and propose three computationally efficient modifications: a reduced kernel, a downsampled receive aperture, and the use of an ensemble correlation coefficient. The proposed methods are implemented in simulation and in vivo studies. Reducing the kernel to a single sample improved computational throughput and improved axial resolution. Downsampling the receive aperture was found to have negligible effect on estimator variance, and improved computational throughput by an order of magnitude for a downsample factor of 4. The ensemble correlation estimator demonstrated lower variance than the currently used average correlation. Combining the three methods, the throughput was improved 105-fold in simulation with a downsample factor of 4- and 20-fold in vivo with a downsample factor of 2.

In Vivo Intracardiac Vector Flow Imaging Using Phased Array Transducers for Pediatric Cardiology

Abstract: Two-dimensional blood speckle tracking (ST) has shown promise for measuring complex flow patterns in neonatal hearts using linear arrays and high-frame-rate plane wave imaging. For general pediatric applications, however, the need for phased array probes emerges due to the limited intercostal acoustic window available. In this paper, a clinically approved real-time duplex imaging setup with phased array probes was used to investigate the potential of blood ST for the 2-D vector flow imaging of children with congenital heart disease. To investigate transmit beam pattern and tracking accuracy, straight tubes with parabolic flow were simulated at three depths (4.5, 7, and 9.5 cm). Due to the small aperture available, diffraction effects could be observed when approaching 10 cm, which limited the number of parallel receive beams that could be utilized. Moving to (slightly) diverging beams was shown to solve this issue at the expense of a loss in signal-to-noise ratio. To achieve consistent estimates, a forward–backward tracking scheme was introduced to avoid measurement bias occurring due to tracking kernel averaging artifacts at flow domain boundaries. Promising results were observed for depths <10 cm in two pediatric patients, where complex cardiac flow patterns could be estimated and visualized. As a loss in penetration compared with color flow imaging is expected, a larger clinical study is needed to establish the clinical feasibility of this approach.

Characterization of Contrast Agent Microbubbles for Ultrasound Imaging and Therapy Research

Abstract: The high efficiency with which gas microbubbles can scatter ultrasound compared with the surrounding blood pool or tissues has led to their widespread employment as contrast agents in ultrasound imaging In recent years, their applications have been extended to include super-resolution imaging and the stimulation of localized bio-effects for therapy The growing exploitation of contrast agents in ultrasound and in particular these recent developments have amplified the need to characterize and fully understand microbubble behavior The aim in doing so is to more fully exploit their utility for both diagnostic imaging and potential future therapeutic applications This paper presents the key characteristics of microbubbles that determine their efficacy in diagnostic and therapeutic applications and the corresponding techniques for their measurement In each case, we have presented information regarding the methods available and their respective strengths and limitations, with the aim of presenting information relevant to the selection of appropriate characterization methods First, we examine methods for determining the physical properties of microbubble suspensions and then techniques for acoustic characterization of both suspensions and single microbubbles The next section covers characterization of microbubbles as therapeutic agents, including as drug carriers for which detailed understanding of their surface characteristics and drug loading capacity is required Finally, we discuss the attempts that have been made to allow comparison across the methods employed by various groups to characterize and describe their microbubble suspensions and promote wider discussion and comparison of microbubble behavior

Regularized Dual Averaging Image Reconstruction for Full-Wave Ultrasound Computed Tomography

Abstract: Ultrasound computed tomography (USCT) holds great promise for breast cancer screening. Waveform inversion-based image reconstruction methods account for higher order diffraction effects and can produce high-resolution USCT images, but are computationally demanding. Recently, a source encoding technique has been combined with stochastic gradient descent (SGD) to greatly reduce image reconstruction times. However, this method bundles the stochastic data fidelity term with the deterministic regularization term. This limitation can be overcome by replacing SGD with a structured optimization method, such as the regularized dual averaging method, that exploits knowledge of the composition of the cost function. In this paper, the dual averaging method is combined with source encoding techniques to improve the effectiveness of regularization while maintaining the reduced reconstruction times afforded by source encoding. It is demonstrated that each iteration can be decomposed into a gradient descent step based on the data fidelity term and a proximal update step corresponding to the regularization term. Furthermore, the regularization term is never explicitly differentiated, allowing nonsmooth regularization penalties to be naturally incorporated. The wave equation is solved by the use of a time-domain method. The effectiveness of this approach is demonstrated through computer simulation and experimental studies. The results suggest that the dual averaging method can produce images with less noise and comparable resolution to those obtained by the use of SGD.

Noise Equalization for Ultrafast Plane Wave Microvessel Imaging

Abstract: Ultrafast plane wave microvessel imaging significantly improves ultrasound Doppler sensitivity by increasing the number of Doppler ensembles that can be collected within a short period of time. The rich spatiotemporal plane wave data also enable more robust clutter filtering based on singular value decomposition. However, due to the lack of transmit focusing, plane wave microvessel imaging is very susceptible to noise. This paper was designed to: 1) study the relationship between ultrasound system noise (primarily time gain compensation induced) and microvessel blood flow signal and 2) propose an adaptive and computationally cost-effective noise equalization method that is independent of hardware or software imaging settings to improve microvessel image quality.

A Prototype Therapy System for Transcutaneous Application of Boiling Histotripsy

Abstract: Boiling histotripsy (BH) is a method of focused ultrasound surgery that noninvasively applies millisecond-length pulses with high-amplitude shock fronts to generate liquefied lesions in tissue. Such a technique requires unique outputs compared to a focused ultrasound thermal therapy apparatus, particularly to achieve high in situ pressure levels through intervening tissue. This paper describes the design and characterization of a system capable of producing the necessary pressure to transcutaneously administer BH therapy through clinically relevant overlying tissue paths using pulses with duration up to 10 ms. A high-voltage electronic pulser was constructed to drive a 1-MHz focused ultrasound transducer to produce shock waves with amplitude capable of generating boiling within the pulse duration in tissue. The system output was characterized by numerical modeling with the 3-D Westervelt equation using boundary conditions established by acoustic holography measurements of the source field. Such simulations were found to be in agreement with directly measured focal waveforms. An existing derating method for nonlinear therapeutic fields was used to estimate in situ pressure levels at different tissue depths. The system was tested in ex vivo bovine liver samples to create BH lesions at depths up to 7 cm. Lesions were also created through excised porcine body wall (skin, adipose, and muscle) with 3–5 cm thickness. These results indicate that the system is capable of producing the necessary output for transcutaneous ablation with BH.

A Frequency-Shift Method to Measure Shear-Wave Attenuation in Soft Tissues

Abstract: In vivo quantification of shear-wave attenuation in soft tissues may help to better understand human tissue rheology and lead to new diagnostic strategies. Attenuation is difficult to measure in acoustic radiation force elastography because the shear-wave amplitude decreases due to a combination of diffraction and viscous attenuation. Diffraction correction requires assuming a cylindrical wavefront and an isotropic propagation medium, which may not be the case in some applications. In this paper, the frequency-shift method, used in ultrasound imaging and seismology, was adapted for shear-wave attenuation measurement in elastography. This method is not sensitive to diffraction effects. For a linear frequency dependence of the attenuation, a closed-form relation was obtained between the decrease in the peak frequency of the gamma-distributed wave amplitude spectrum and the attenuation coefficient of the propagation medium. The proposed method was tested against a plane-wave reference method in homogeneous agar–gelatin phantoms with 0%, 10%, and 20% oil concentrations, and hence different attenuations of 0.117, 0.202, and 0.292 $\text {Np}\cdot \text {m}^{-1}$ /Hz, respectively. Applicability to biological tissues was demonstrated with two ex vivo porcine liver samples (0.79 and 1.35 $\text {Np} \,\cdot \, \text {m}^{-1}$ /Hz) and an in vivo human muscle, measured along (0.43 $\text {Np}\,\cdot \, \text {m}^{-1}$ /Hz) and across (1.77 $\text {Np}\cdot \text {m}^{-1}$ /Hz) the tissue fibers. In all cases, the data supported the assumptions of a gamma-distributed spectrum for the source and linear frequency attenuation for the tissue. This method provides tissue attenuation, which is relevant diagnostic information to model viscosity, in addition to shear-wave velocity used to assess elasticity. Data processing is simple and could be performed automatically in real time for clinical applications.

Wideband 2-D Array Design Optimization With Fabrication Constraints for 3-D US Imaging

Abstract: Ultrasound (US) 2-D arrays are of increasing interest due to their electronic steering capability to investigate 3-D regions without requiring any probe movement. These arrays are typically populated by thousands of elements that, ideally, should be individually driven by the companion scanner. Since this is not convenient, the so-called microbeamforming methods, yielding a prebeamforming stage performed in the probe handle by suitable custom integrated circuits, have so far been implemented in a few commercial high-end scanners. A possible approach to implement relatively cheap and efficient 3-D US imaging systems is using 2-D sparse arrays in which a limited number of elements can be coupled to an equal number of independent transmit/receive channels. In order to obtain US beams with adequate characteristics all over the investigated volume, the layout of such arrays must be carefully designed. This paper provides guidelines to design, by using simulated annealing optimization, 2-D sparse arrays capable of fitting specific applications or fabrication/implementation constraints. In particular, an original energy function based on multidepth 3-D analysis of the beam pattern is also exploited. A tutorial example is given, addressed to find the $N_{\mathrm{ e}}$ elements that should be activated in a 2-D fully populated array to yield efficient acoustic radiating performance over the entire volume. The proposed method is applied to a $32 \times 32$ array centered at 3 MHz to select the 128, 192, and 256 elements that provide the best acoustic performance. It is shown that the 256-element optimized array yields sidelobe levels even lower (by 5.7 dB) than that of the reference 716-element circular and (by 10.3 dB) than that of the reference 1024-element array.

Ultrasonic 3-D Vector Flow Method for Quantitative In Vivo Peak Velocity and Flow Rate Estimation

Abstract: Current clinical ultrasound (US) systems are limited to show blood flow movement in either 1-D or 2-D. In this paper, a method for estimating 3-D vector velocities in a plane using the transverse oscillation method, a $32\times32$ element matrix array, and the experimental US scanner SARUS is presented. The aim of this paper is to estimate precise flow rates and peak velocities derived from 3-D vector flow estimates. The emission sequence provides 3-D vector flow estimates at up to 1.145 frames/s in a plane, and was used to estimate 3-D vector flow in a cross-sectional image plane. The method is validated in two phantom studies, where flow rates are measured in a flow-rig, providing a constant parabolic flow, and in a straight-vessel phantom ( ${\oslash }=8$ mm) connected to a flow pump capable of generating time varying waveforms. Flow rates are estimated to be 82.1 ± 2.8 L/min in the flow-rig compared with the expected 79.8 L/min, and to 2.68 ± 0.04 mL/stroke in the pulsating environment compared with the expected 2.57 ± 0.08 mL/stroke. Flow rates estimated in the common carotid artery of a healthy volunteer are compared with magnetic resonance imaging (MRI) measured flow rates using a 1-D through-plane velocity sequence. Mean flow rates were 333 ± 31 mL/min for the presented method and 346 ± 2 mL/min for the MRI measurements.

Cardiac Strain Imaging With Coherent Compounding of Diverging Waves

Abstract: Current methods of cardiac strain imaging at high frame rate suffer from motion matching artifacts or poor lateral resolution. Coherent compounding has been shown to improve echocardiographic image quality while maintaining a high frame rate, but has never been used to image cardiac strain. However, myocardial velocity can have an impact on coherent compounding due to displacements between frames. The objective of this paper was to investigate the feasibility and performance of coherent compounding for cardiac strain imaging at a low and a high myocardial velocity. Left-ventricular contraction in short-axis view was modeled as an annulus with radial thickening and circumferential rotation. Simulated radio-frequency channel data with a cardiac phased array were obtained using three different beamforming methods: single diverging wave, coherent compounding of diverging waves, and conventional focusing. Axial and lateral displacements and strains as well as radial strains were estimated and compared to their true value. In vivo feasibility of cardiac strain imaging with coherent compounding was performed and compared to single diverging wave imaging. At low myocardial velocities, the axial, lateral, and radial strain relative error for nine compounded waves (16.3%, 40.4%, and 18.9%) were significantly lower than those obtained with single diverging wave imaging (19.9%, 80.3%, and 30.6%) and closer to that obtained with conventional focusing (16.7%, 43.7%, and 16%). In vivo left-ventricular radial strains exhibited higher quality with nine compounded waves than with single diverging wave imaging. These results indicate that cardiac strain can be imaged using coherent compounding of diverging waves with a better performance than with single diverging wave imaging while maintaining a high frame rate, and therefore, has the potential to improve diagnosis of myocardial strain-based cardiac diseases.

The Impact of Model-Based Clutter Suppression on Cluttered, Aberrated Wavefronts

Abstract: Recent studies reveal that both phase aberration and reverberation play a major role in degrading ultrasound image quality. We previously developed an algorithm for suppressing clutter, but we have not yet tested it in the context of aberrated wavefronts. In this paper, we evaluate our previously reported algorithm, called aperture domain model image reconstruction (ADMIRE), in the presence of phase aberration and in the presence of multipath scattering and phase aberration. We use simulations to investigate phase aberration corruption and correction in the presence of reverberation. As part of this paper, we observed that ADMIRE leads to suppressed levels of aberration. In order to accurately characterize aberrated signals of interest, we introduced an adaptive component to ADMIRE to account for aberration, referred to as adaptive ADMIRE. We then use ADMIRE, adaptive ADMIRE, and conventional filtering methods to characterize aberration profiles on in vivo liver data. These in vivo results suggest that adaptive ADMIRE could be used to better characterize a wider range of aberrated wavefronts. The aberration profiles’ full-width at half-maximum of ADMIRE, adaptive ADMIRE, and postfiltered data with 0.4- ${\rm mm}^{-1}$ spatial cutoff frequency are 4.0 ± 0.28 mm, 2.8 ± 1.3 mm, and 2.8 ± 0.57 mm, respectively, while the average root-mean square values in the same order are 16 ± 5.4 ns, 20 ± 6.3 ns, and 19 ± 3.9 ns, respectively. Finally, because ADMIRE suppresses aberration, we perform a limited evaluation of image quality using simulations and in vivo data to determine how ADMIRE and adaptive ADMIRE perform with and without aberration correction.

Wireless Power Transfer to Millimeter-Sized Nodes Using Airborne Ultrasound

Abstract: We propose the use of airborne ultrasound for wireless power transfer to mm-sized nodes, with intended application in the next generation of the Internet of Things (IoT). We show through simulation that ultrasonic power transfer can deliver 50 $\mu \text{W}$ to a mm-sized node 0.88 m away from a ~ 50-kHz, 25-cm2 transmitter array, with the peak pressure remaining below recommended limits in air, and with load power increasing with transmitter area. We report wireless power recovery measurements with a precharged capacitive micromachined ultrasonic transducer, demonstrating a load power of 5 $\mu \text{W}$ at a simulated distance of 1.05 m. We present aperture efficiency, dynamic range, and bias-free operation as key metrics for the comparison of transducers meant for wireless power recovery. We also argue that long-range wireless charging at the watt level is extremely challenging with existing technology and regulations. Finally, we compare our acoustic powering system with cutting edge electromagnetically powered nodes and show that ultrasound has many advantages over RF as a vehicle for power delivery. Our work sets the foundation for further research into ultrasonic wireless power transfer for the IoT.

Intracardiac Vortex Dynamics by High-Frame-Rate Doppler Vortography-In Vivo Comparison With Vector Flow Mapping and 4-D Flow MRI.

Abstract: Recent studies have suggested that intracardiac vortex flow imaging could be of clinical interest to early diagnose the diastolic heart function. Doppler vortography has been introduced as a simple color Doppler method to detect and quantify intraventricular vortices. This method is able to locate a vortex core based on the recognition of an antisymmetric pattern in the Doppler velocity field. Because the heart is a fast-moving organ, high frame rates are needed to decipher the whole blood vortex dynamics during diastole. In this paper, we adapted the vortography method to high-frame-rate echocardiography using circular waves. Time-resolved Doppler vortography was first validated in vitro in an ideal forced vortex. We observed a strong correlation between the core vorticity determined by high-frame-rate vortography and the ground-truth vorticity. Vortography was also tested in vivo in ten healthy volunteers using high-frame-rate duplex ultrasonography. The main vortex that forms during left ventricular filling was tracked during two-three successive cardiac cycles, and its core vorticity was determined at a sampling rate up to 80 duplex images per heartbeat. Three echocardiographic apical views were evaluated. Vortography-derived vorticities were compared with those returned by the 2-D vector flow mapping approach. Comparison with 4-D flow magnetic resonance imaging was also performed in four of the ten volunteers. Strong intermethod agreements were observed when determining the peak vorticity during early filling. It is concluded that high-frame-rate Doppler vortography can accurately investigate the diastolic vortex dynamics.