Compensation of array lens effects for improved co-registration of passive acoustic mapping and B-mode images for cavitation monitoring

Passive Acoustic Mapping for Convex Arrays With the Helical Wave Spectrum Method

Passive acoustic mapping (PAM) has emerged as a valuable imaging modality for monitoring the cavitation activity in focused ultrasound therapies. When it comes to imaging in the human abdomen, convex arrays are preferred due to their large acoustic window. However, existing PAM methods for convex arrays rely on the computationally expensive delay-and-sum (DAS) operation limiting the image reconstruction speed when the field-of-view (FOV) is large. In this work, we propose an efficient and frequency-selective PAM method for convex arrays. This method is based on projecting the helical wave spectrum (HWS) between cylindrical surfaces in the imaging field. Both the in silico and in vitro experiments showed that the HWS method has comparable image quality and similar acoustic cavitation source localization accuracy as the DAS-based methods. Compared to the frequency-domain and time-domain DAS methods, the time-complexity of the HWS method is reduced by one order and two orders of magnitude, respectively. A parallel implementation of the HWS method realized millisecond-level image reconstruction speed. We also show that the HWS method is inherently capable of mapping microbubble (MB) cavitation activity of different status, i.e., no cavitation, stable cavitation, or inertial cavitation. After compensating for the lens effects of the convex array, we further combined PAM formed by the HWS method and B-mode imaging as a real-time dual-mode imaging approach to map the anatomical location where MBs cavitate in a liver phantom experiment. This method may find use in applications where convex arrays are required for cavitation activity monitoring in real time.

Passive Cavitation Imaging Artifact Reduction Using Data-Adaptive Spatial Filtering

Passive cavitation imaging (PCI) with a clinical diagnostic array results in poor axial localization of bubble activity due to the size of the point spread function (PSF). The objective of this study was to determine if data-adaptive spatial filtering improved PCI beamforming performance relative to standard frequency-domain delay, sum, and integrate (DSI) or robust Capon beamforming (RCB). The overall goal was to improve source localization and image quality without sacrificing computation time. Spatial filtering was achieved by applying a pixel-based mask to DSI- or RCB-beamformed images. The masks were derived from DSI, RCB, or phase or amplitude coherence factors (ACFs) using both receiver operating characteristic (ROC) and precision-recall (PR) curve analyses. Spatially filtered passive cavitation images were formed from cavitation emissions based on two simulated sources densities and four source distribution patterns mimicking cavitation emissions induced by an EkoSonic catheter. Beamforming performance was assessed via binary classifier metrics. The difference in sensitivity, specificity, and area under the ROC curve (AUROC) differed by no more than 11% across all algorithms for both source densities and all source patterns. The computational time required for each of the three spatially filtered DSIs was two orders of magnitude less than that required for time-domain RCB and thus this data-adaptive spatial filtering strategy for PCI beamforming is preferable given the similar binary classification performance.

Passive Cavitation Mapping by Cavitation Source Localization From Aperture-Domain Signals-Part I: Theory and Validation Through Simulations

Passive cavitation mapping (PCM) algorithms for diagnostic ultrasound arrays based on time exposure acoustics (TEA) exhibit poor axial resolution, which is in part due to the diffraction-limited point spread function of the imaging system and poor rejection by the delay-and-sum beamformer. In this article, we adapt a method for speed of sound estimation to be utilized as a cavitation source localization (CSL) approach. This method utilizes a hyperbolic fit to the arrival times of the cavitation signals in the aperture domain, and the coefficients of the fit are related to the position of the cavitation source. Wavefronts exhibiting poor fit to the hyperbolic function are corrected to yield improved source localization. We demonstrate through simulations that this method is capable of accurate estimation of the origin of coherent spherical waves radiating from cavitation/point sources. The average localization error from simulated microbubble sources was 0.12 ± 0.12mm ( 0.15 ± 0.14λ₀ for a 1.78-MHz transmit frequency). In simulations of two simultaneous cavitation sources, the proposed technique had an average localization error of 0.2mm ( 0.23λ₀ ), whereas conventional TEA had an average localization error of 0.81mm ( 0.97λ₀ ). The reconstructed PCM-CSL image showed a significant improvement in resolution compared with the PCM-TEA approach.

Dual-Array Passive Acoustic Mapping for Cavitation Imaging With Enhanced 2-D Resolution

Passive acoustic mapping (PAM) techniques have been developed for the purposes of detecting, localizing, and quantifying cavitation activity during therapeutic ultrasound procedures. Implementation with conventional diagnostic ultrasound arrays has allowed planar mapping of bubble acoustic emissions to be overlaid with B-mode anatomical images, with a variety of beamforming approaches providing enhanced resolution at the cost of extended computation times. However, no passive signal processing techniques implemented to date have overcome the fundamental physical limitation of the conventional diagnostic array aperture that results in point spread functions with axial/lateral beamwidth ratios of nearly an order of magnitude. To mitigate this problem, the use of a pair of orthogonally oriented diagnostic arrays was recently proposed, with potential benefits arising from the substantially expanded range of observation angles. This article presents experiments and simulations intended to demonstrate the performance and limitations of the dual-array system concept. The key finding of this study is that source pair resolution of better than 1 mm is now possible in both dimensions of the imaging plane using a pair of 7.5-MHz center frequency conventional arrays at a distance of 7.6cm. With an eye toward accelerating computations for real-time applications, channel count reductions of up to a factor of eight induce negligible performance losses. Modest sensitivities to sound speed and relative array position uncertainties were identified, but if these can be kept on the order of 1% and 1 mm, respectively, then the proposed methods offer the potential for a step improvement in cavitation monitoring capability.

Spatiotemporal Assessment of the Cellular Safety of Cavitation-Based Therapies by Passive Acoustic Mapping

Many useful therapeutic bio-effects can be generated using ultrasound-induced cavitation. However, cavitation is also capable of causing unwanted cellular and vascular damage, which should be monitored to ensure treatment safety. In this work, the unique opportunity provided by passive acoustic mapping (PAM) to quantify cavitation dose across an entire volume of interest during therapy is utilised to provide setup-independent measures of spatially localised cavitation dose. This spatiotemporally quantifiable cavitation dose is then related to the level of cellular damage generated. The cavitation-mediated destruction of equine red blood cells mixed with one of two types of cavitation nuclei at a variety of concentrations is investigated. The blood is placed within a 0.5-MHz ultrasound field and exposed to a range of peak rarefactional pressures up to 2 MPa, with 50 to 50,000 cycle pulses maintaining a 5% duty cycle. Two co-planar linear arrays at 90° to each other are used to generate 400-µm-resolution frequency domain robust capon beamforming PAM maps, which are then used to generate estimates of cavitation dose. A relationship between this cavitation dose and the levels of haemolysis generated was found which was comparable regardless of the applied acoustic pressure, pulse length, cavitation agent type or concentration used. PAM was then used to monitor cellular damage in multiple locations within a tissue phantom simultaneously, with the damage-cavitation dose relationship being similar for the two experimental models tested. These results lay the groundwork for this method to be applied to other measures of safety, allowing for improved ultrasound monitoring of cavitation-based therapies.