Optimal focusing by spatio-temporal inverse filter. I. Basic principles

Innovative aberration correction in ultrasound diagnostics with direct phase estimation for enhanced image quality

The paper addresses a crucial challenge in medical radiology and introduces a novel general approach, which utilises applied mathematics and information technology techniques, for aberration correction in ultrasound diagnostics. Ultrasound imaging of inhomogeneous media inherently suffers from variations in ultrasonic speed between tissue. The characteristics of aberrations are unique to each patient due to tissue morphology. This study proposes a new phase aberration correction method based on the Fourier transform and leveraging of the synthetic aperture mode. The proposed method enables correction after the emission and reception of ultrasonic wave, allowing for the estimation of aberration profiles for different parts of the sonogram. To demonstrate the method's performance, this study included the conducting of experiments using a commercially available quality control phantom, an ex-vivo temporal human bone, and specially designed distortion layers. At a frequency of 2 MHz, the experiments demonstrated an increase of two-and-three-quarters in echo signal intensity and a decrease of nearly two-fold in the width of the angular distribution compared to the pre-correction state. However, it is important to note that the implementation of the method has a limitation, as it requires an aperture synthesis mode and access to raw RF data, which restricts use in common scanners. To ensure the reproducibility of the results, this paper provides public access to an in-house C +  + code for aberration correction following the proposed method, as well as the dataset used in this study.

Field and Temperature Shaping for Microwave Hyperthermia: Recent Treatment Planning Tools to Enhance SAR-Based Procedures

The aim of the article is to provide a summary of the work carried out in the framework of a research project funded by the Italian Ministry of Research. The main goal of the activity was to introduce multiple tools for reliable, affordable, and high-performance microwave hyperthermia for cancer therapy. The proposed methodologies and approaches target microwave diagnostics, accurate in vivo electromagnetic parameters estimation, and improvement in treatment planning using a single device. This article provides an overview of the proposed and tested techniques and shows their complementarity and interconnection. To highlight the approach, we also present a novel combination of specific absorption rate optimization via convex programming with a temperature-based refinement method implemented to mitigate the effect of thermal boundary conditions on the final temperature map. To this purpose, numerical tests were carried out for both simple and anatomically detailed 3D scenarios for the head and neck region. These preliminary results show the potential of the combined technique and improvements in the temperature coverage of the tumor target with respect to the case wherein no refinement is adopted.

Ultrasonic array imaging of highly attenuative materials with spatio-temporal singular value decomposition

With the increasing use of high density polyethylene (HDPE) pipes in nuclear industry which puts safety at the forefront, ultrasonic array imaging methods play a vital role in the structural integrity of HDPE pipe materials. However, the viscoelastic attenuation of HDPE pipe materials significantly decreases the level of signals, leading to a low signal-to-noise ratio caused by electronic noise. In this work, a domain-adapted spatio-temporal singular value decomposition (STSVD) processing algorithm combined with the total focusing method is proposed to improve the ultrasonic array image quality. First, the real-valued radio frequency (RF) data or A-scan signals are demodulated into the complex analytic signals containing in-phase and quadrature (I/Q) components. Then, the STSVD processing algorithm is used to filter the I/Q data, and the filtered I/Q data is converted into RF signals. Finally, the total focusing method is applied to the processed RF signals to produce the image of the region under detection as a stage of post-processing. Experiments are carried out with an ultrasonic linear phased array in contact with the HDPE pipe materials containing multiple side-drilled holes and through-wall notches. Results show that the proposed method can produce images with high quality to provide good inspection and characterization of defects in highly attenuative materials, especially the deeper defects.

The physics of knocking over LEGO minifigures with time reversal focused vibrations for use in a museum exhibit

Time reversal (TR) is a method of focusing wave energy at a point in space. The optimization of a TR demonstration is described, which knocks over one selected LEGO minifigure among other minifigures by focusing the vibrations within an aluminum plate at the target minifigure. The aim is to achieve a high repeatability of the demonstration along with reduced costs to create a museum exhibit. By comparing the minifigure's motion to the plate's motion directly beneath its feet, it is determined that a major factor inhibiting the repeatability is that the smaller vibrations before the focal event cause the minifigure to bounce repeatedly and it ends up being in the air during the main vibrational focal event, which was intended to launch the minifigure. The deconvolution TR technique is determined to be optimal in providing the demonstration repeatability. The amplitude, frequency, and plate thickness are optimized in a laboratory setting. An eddy current sensor is then used to reduce the costs, and the impact on the repeatability is determined. A description is given of the implementation of the demonstration for a museum exhibit. This demonstration illustrates the power of the focusing acoustic waves, and the principles learned by optimizing this demonstration can be applied to other real-world applications.

Imaging through a scattering medium: the Fisher information and the generalized Abbe limit

Enhanced-resolution imaging in complex scattering media is revisited from a parameter estimation perspective. A suitably defined Fisher information is shown to offer useful insights into the limiting precision of parameter estimation in a scattering environment and, hence, into the limiting spatial resolution that can be achieved in imaging-through-scattering settings. The Fisher information that defines this resolution limit via the Cramér-Rao lower bound is shown to scale with the number of adaptively controlled space-time modes of the probe field, suggesting a physically intuitive generalization of the Abbe limit to the spatial resolution attainable for complex scattering systems. In a conventional, direct-imaging microscopy setting, this bound is shown to converge to the canonical Abbe limit.

The SVD beamformer with diverging waves: a proof-of-concept for fast aberration correction

The rise of ultrafast ultrasound imaging-with plane or diverging waves - paved the way to new applications of ultrasound in biomedical applications. However, propagation through complex layers (typically fat, muscle, and bone) hinder considerably the image quality, especially because of sound speed heterogeneities. In difficult-to-image patients, in the case of the hepatic steatosis for instance, a good image and a reliable sound speed quantification are crucial to provide a powerful non-invasive diagnosis tool. In this work, we proposed to adapt the singular value decomposition (SVD) beamformer method for diverging waves and thus present a novel aberration correction approach for widely used curved arrays. We probed its efficiency experimentally bothin vitroandin vivo. Besides the proposed matrix formalism, we explored the physical meaning of the SVD of ultrafast data. Finally, we demonstrated the ability of the technique to improve the image quality and offer new perspectives particularly in quantitative liver imaging.

Comparison Between Ray-Tracing and Full-Wave Simulation for Transcranial Ultrasound Focusing on a Clinical System Using the Transfer Matrix Formalism

Only one high-intensity focused ultrasound device has been clinically approved for transcranial brain surgery at the time of writing. The device operates within 650 and 720 kHz and corrects the phase distortions induced by the skull of each patient using a multielement phased array. Phase correction is estimated adaptively using a proprietary algorithm based on computed-tomography (CT) images of the patient's skull. In this article, we assess the performance of the phase correction computed by the clinical device and compare it to: 1) the correction obtained with a previously validated full-wave simulation algorithm using an open-source pseudo-spectral toolbox and 2) a hydrophone-based correction performed invasively to measure the aberrations induced by the skull at 650 kHz. For the full-wave simulation, three different mappings between CT Hounsfield units and the longitudinal speed of sound inside the skull were tested. All methods are compared with the exact same setup due to transfer matrices acquired with the clinical system for N = 5 skulls and T = 2 different targets for each skull. We show that the clinical ray-tracing software and the full-wave simulation restore, respectively, 84% ± 5% and 86% ± 5% of the pressure obtained with hydrophone-based correction for targets located in central brain regions. On the second target (off-center), we also report that the performance of both algorithms degrades when the average incident angles of the acoustic beam at the skull surface increase. When incident angles are higher than 20°, the restored pressure drops below 75% of the pressure restored with hydrophone-based correction.

The SVD Beamformer: Physical Principles and Application to Ultrafast Adaptive Ultrasound

A shift of paradigm is currently underway in biomedical ultrasound thanks to plane or diverging waves coherent compounding for faster imaging. One remaining challenge consists in handling phase and amplitude aberrations induced during the ultrasonic propagation through complex layers. Unlike conventional line-per-line imaging, ultrafast ultrasound provides backscattering information from the whole imaged area for each transmission. Here, we take benefit from this feature and propose an efficient approach to perform fast aberration correction. Our method is based on the Singular Value Decomposition (SVD) of an ultrafast compound matrix containing backscattered data for several plane wave transmissions. First, we explain the physical signification of SVD and associated singular vectors within the ultrafast matrix formalism. We theoretically demonstrate that the separation of spatial and angular variables, rendered by SVD on ultrafast data, provides an elegant and straightforward way to optimize the angular coherence of backscattered data. In heterogeneous media, we demonstrate that the first spatial and angular singular vectors retrieve respectively the non-aberrated image of a region of interest, and the phase and amplitude of its aberration law. Numerical, in vitro and in vivo results prove the efficiency of the image correction, but also the accuracy of the aberrator determination. Based on spatial and angular coherence, we introduce a complete methodology for adaptive beamforming of ultrafast data, performed on successive isoplanatism patches undergoing SVD beamforming. The simplicity of this method paves the way to real-time adaptive ultrafast ultrasound imaging and provides a theoretical framework for future quantitative ultrasound applications.

Multitouch Vibrotactile Feedback on a Tactile Screen by the Inverse Filter Technique: Vibration Amplitude and Spatial Resolution

Nowadays, tactile surfaces, such as smartphones, provide haptic feedback to signify that a task has been performed correctly or more generally to enrich the interaction. However, this haptic feedback induces vibrations in the surface that propagate to the whole surface, reverberate and attenuate, thus making multi-finger interaction, with different feedbacks, difficult. Recently, the Inverse Filter Method has been proposed to control the propagation of these vibrations, and thus enable to product localized multitouch on a glass surface. This way, a user can put several fingers on a tactile surface and yet feel stimuli independently on his/her different fingers. This article continues this work and demonstrates that a localized multitouch haptic feedback can be delivered in real time using a capacitive screen. To achieve this, this article presents the two necessary steps: a calibration step and an interpolation calculation in order to save calculation and learning time. Furthermore, the paper describes the performance of the device through a study on the behaviour of the screen subjected to the Inverse Filter Method, indicating the movement of the whole screen and the voltage requirement for any haptic feedback.

Confinement of Vibrotactile Stimuli in Narrow Plates: Principle and Effect of Finger Loading

On a touch surface, providing a local vibrotactile feedback enables multiuser and multitouch interactions. While the vibration propagation usually impedes this localization, we show in this article that narrow strip-shaped plates constitute waveguides in which bending waves below a cut-off frequency do not propagate. We provide a theoretical explanation of the phenomenon and experimental validations. We thus show that vibrations up to a few kHz are well confined on top of the actuated area with vibration amplitude over 1 micrometer. The principle was validated with piezoelectric actuators of various shapes and a linear resonant actuator (LRA). Investigation of the effect of a fingertip load on the system through theory and experimentation was conducted and revealed that almost no attenuation was brought by the fingertip when using low frequency evanescent waves. Finally, a perceptual validation was conducted and showed dynamic stimuli with a large frequency spectrum could be felt and distinguished.

A Noninvasive Ultrasound Resonance Method for Detecting Skull Induced Phase Shifts May Provide a Signal for Adaptive Focusing

OBJECTIVE

There may be a need to perform dynamic skull aberration corrections during the non-invasive high-intensity transcranial treatment with magnetic resonance imaging (MRI) -guided focused ultrasound in order to accurately and rapidly restore the focus in the brain.

METHODS

This could possibly be accomplished by using an ultrasound-based correction method based on the skulls' thickness resonance frequencies. The focus of a 500 kHz transducer was centered in the ex vivo human skull caps at different temperatures. The pulse-echoed signals reflected from the skulls were analyzed in the frequency domain to reveal the resonance frequencies for the phase shift calculation. The accuracy was compared to both hydrophone and computed tomography (CT) based analytical methods.

RESULTS

Around 73% of the measurements (n = 784) were in the optimal constructive interference region, with a 15° decrease in the average phase error compared to the previous study. In the best implementation, it performed approximately the same or better than the CT based analytical method currently in clinical use. Linear correlation was found between the resonance frequencies or skull induced phase shifts and the skull temperature with an average rate of -0.4 kHz/°C and 2.6 deg/°C, respectively.

CONCLUSION

The ultrasound based resonance method has shown the feasibility of detecting heating-induced changes of skull phase shift non-invasively and accurately.

SIGNIFICANCE

Since the technique can be made MRI compatible and integrated in the therapy arrays, it may enable temperature tracking and adaptive focusing during high-intensity transcranial ultrasound treatments, to prevent skull overheating and preserve the transcranial focusing integrity.

The potential of constrained SAR focusing for hyperthermia treatment planning: analysis for the head & neck region

Clinical trials have shown that hyperthermia is a potent adjuvant to conventional cancer treatments, but the temperatures currently achieved in the clinic are still suboptimal. Hyperthermia treatment planning simulations have potential to improve the heating profile of phased-array applicators. An important open challenge is the development of an effective optimization procedure that enables uniform heating of the target region while keeping temperature below a threshold in healthy tissues. In this work, we analyzed the effectiveness and efficiency of a recently proposed optimization approach, i.e. focusing via constrained power optimization (FOCO), using 3D simulations of twelve clinical patient specific models. FOCO performance was compared against a clinically used particle swarm based optimization approach. Evaluation metrics were target coverage at the 25% iso-SAR level, target hotspot quotient, median target temperature (T50) and computational requirements. Our results show that, on average, constrained power focusing performs slightly better than the clinical benchmark ([Formula: see text]T50 [Formula: see text] °C), but outperforms this clinical benchmark for large target volumes ([Formula: see text]40 cm[Formula: see text], [Formula: see text]T50 [Formula: see text] °C). In addition, the results are achieved in a shorter time ([Formula: see text]%) and are repeatable because the approach is formulated as a convex optimization problem.

Theoretical characterization of annular array as a volumetric optoacoustic ultrasound handheld probe

Optoacoustic ultrasound (OPUS) is a promising hybridized technique for simultaneous acquisition of functional and morphological data. The optical specificity of optoacoustic leverages the diagnostic aptitude of ultrasonography beyond anatomy. However, this integration has been rarely practiced for volumetric imaging. The challenge lies in the effective imaging probes that preserve the functionality of both modalities. The potentials of a sparse annular array for volumetric OPUS imaging are theoretically investigated. In order to evaluate and optimize the performance characteristics of the probe, series of analysis in the framework of system model matrix was carried out. The two criteria of voxel crosstalk and eigenanalysis have been employed to unveil information about the spatial sensitivity, aliasing, and number of definable spatial frequency components. Based on these benchmarks, the optimal parameters for volumetric handheld probe are determined. In particular, the number, size, and the arrangement of the elements and overall aperture dimension were investigated. The result of the numerical simulation suggests that the segmented-annular array of 128 negatively focused elements with 1λ  ×  20λ size, operating at 5-MHz central frequency showcases a good agreement with the physical requirement of both imaging systems. We hypothesize that these features enable a high-throughput volumetric passive/active ultrasonic imaging system with great potential for clinical applications.

Time reversal focusing of high amplitude sound in a reverberation chamber

Time reversal (TR) is a signal processing technique that can be used for intentional sound focusing. While it has been studied in room acoustics, the application of TR to produce a high amplitude focus of sound in a room has not yet been explored. The purpose of this study is to create a virtual source of spherical waves with TR that are of sufficient intensity to study nonlinear acoustic propagation. A parameterization study of deconvolution, one-bit, clipping, and decay compensation TR methods is performed to optimize high amplitude focusing and temporal signal focus quality. Of all TR methods studied, clipping is shown to produce the highest amplitude focal signal. An experiment utilizing eight horn loudspeakers in a reverberation chamber is done with the clipping TR method. A peak focal amplitude of 9.05 kPa (173.1 dB peak re 20 μPa) is achieved. Results from this experiment indicate that this high amplitude focusing is a nonlinear process.

3D-printed adaptive acoustic lens as a disruptive technology for transcranial ultrasound therapy using single-element transducers

The development of multi-element arrays for better control of the shape of ultrasonic beams has opened the way for focusing through highly aberrating media, such as the human skull. As a result, the use of brain therapy with transcranial-focused ultrasound has rapidly grown. Although effective, such technology is expensive. We propose a disruptive, low-cost approach that consists of focusing a 1 MHz ultrasound beam through a human skull with a single-element transducer coupled with a tailored silicone acoustic lens cast in a 3D-printed mold and designed using computed tomography-based numerical acoustic simulation. We demonstrate on N  =  3 human skulls that adding lens-based aberration correction to a single-element transducer increases the deposited energy on the target 10 fold.

An integrated and highly sensitive ultrafast acoustoelectric imaging system for biomedical applications

Direct imaging of the electrical activation of the heart is crucial to better understand and diagnose diseases linked to arrhythmias. This work presents an ultrafast acoustoelectric imaging (UAI) system for direct and non-invasive ultrafast mapping of propagating current densities using the acoustoelectric effect. Acoustoelectric imaging is based on the acoustoelectric effect, the modulation of the medium's electrical impedance by a propagating ultrasonic wave. UAI triggers this effect with plane wave emissions to image current densities. An ultrasound research platform was fitted with electrodes connected to high common-mode rejection ratio amplifiers and sampled by up to 128 independent channels. The sequences developed allow for both real-time display of acoustoelectric maps and long ultrafast acquisition with fast off-line processing. The system was evaluated by injecting controlled currents into a saline pool via copper wire electrodes. Sensitivity to low current and low acoustic pressure were measured independently. Contrast and spatial resolution were measured for varying numbers of plane waves and compared to line per line acoustoelectric imaging with focused beams at equivalent peak pressure. Temporal resolution was assessed by measuring time-varying current densities associated with sinusoidal currents. Complex intensity distributions were also imaged in 3D. Electrical current densities were detected for injected currents as low as 0.56 mA. UAI outperformed conventional focused acoustoelectric imaging in terms of contrast and spatial resolution when using 3 and 13 plane waves or more, respectively. Neighboring sinusoidal currents with opposed phases were accurately imaged and separated. Time-varying currents were mapped and their frequency accurately measured for imaging frame rates up to 500 Hz. Finally, a 3D image of a complex intensity distribution was obtained. The results demonstrated the high sensitivity of the UAI system proposed. The plane wave based approach provides a highly flexible trade-off between frame rate, resolution and contrast. In conclusion, the UAI system shows promise for non-invasive, direct and accurate real-time imaging of electrical activation in vivo.

Ex vivo optimisation of a heterogeneous speed of sound model of the human skull for non-invasive transcranial focused ultrasound at 1 MHz

Transcranial brain therapy has recently emerged as a non-invasive strategy for the treatment of various neurological diseases, such as essential tremor or neurogenic pain. However, treatments require millimetre-scale accuracy. The use of high frequencies (typically ≥1 MHz) decreases the ultrasonic wavelength to the millimetre scale, thereby increasing the clinical accuracy and lowering the probability of cavitation, which improves the safety of the technique compared with the use of low-frequency devices that operate at 220 kHz. Nevertheless, the skull produces greater distortions of high-frequency waves relative to low-frequency waves. High-frequency waves require high-performance adaptive focusing techniques, based on modelling the wave propagation through the skull. This study sought to optimise the acoustical modelling of the skull based on computed tomography (CT) for a 1 MHz clinical brain therapy system. The best model tested in this article corresponded to a maximum speed of sound of 4000 m.s^-1 in the skull bone, and it restored 86% of the optimal pressure amplitude on average in a collection of six human skulls. Compared with uncorrected focusing, the optimised non-invasive correction led to an average increase of 99% in the maximum pressure amplitude around the target and an average decrease of 48% in the distance between the peak pressure and the selected target. The attenuation through the skulls was also assessed within the bandwidth of the transducers, and it was found to vary in the range of 10 ± 3 dB at 800 kHz and 16 ± 3 dB at 1.3 MHz.

Time reversal focusing of elastic waves in plates for an educational demonstration

The purpose of this research is to develop a visual demonstration of time reversal focusing of vibrations in a thin plate. Various plate materials are tested to provide optimal conditions for time reversal focusing. Specifically, the reverberation time in each plate and the vibration coupling efficiency from a shaker to the plate are quantified to illustrate why a given plate provides the best spatially confined focus as well as the highest focal amplitude possible. A single vibration speaker and a scanning laser Doppler vibrometer (SLDV) are used to provide the time reversal focusing. Table salt is sprinkled onto the plate surface to allow visualization of the high amplitude, spatially localized time reversal focus; the salt is thrown upward only at the focal position. Spatial mapping of the vibration focusing on the plate using the SLDV is correlated to the visual salt jumping demonstration. The time reversal focusing is also used to knock over an object when the object is placed at the focal position; some discussion of optimal objects to use for this demonstration are given.

SAW Synthesis With IDTs Array and the Inverse Filter: Toward a Versatile SAW Toolbox for Microfluidics and Biological Applications

Surface acoustic waves (SAWs) are versatile tools to manipulate fluids at small scales for microfluidics and biological applications. A nonexhaustive list of operations that can be performed with SAW includes sessile droplet displacement, atomization, division, and merging but also the actuation of fluids embedded in microchannels or the manipulation of suspended particles. However, each of these operations requires a specific design of the wave generation system, the so-called interdigitated transducers (IDTs). Depending on the application, it might indeed be necessary to generate focused or plane, propagating or standing, and aligned or shifted waves. Furthermore, the possibilities offered by more complex wave fields such as acoustical vortices for particle tweezing and liquid twisting cannot be explored with classical IDTs. In this paper, we show that the inverse filter technique coupled with an IDTs array enables us to synthesize all classical wave fields used in microfluidics and biological applications with a single multifunctional platform. It also enables us to generate swirling SAWs, whose potential for the on-chip synthesis of tailored acoustical vortices has been demonstrated lately. The possibilities offered by this platform are illustrated by performing many operations successively on sessile droplets with the same system.

Observation of a Single-Beam Gradient Force Acoustical Trap for Elastic Particles: Acoustical Tweezers

We demonstrate the trapping of elastic particles by the large gradient force of a single acoustical beam in three dimensions. Acoustical tweezers can push, pull and accurately control both the position and the forces exerted on a unique particle. Forces in excess of 1 micronewton were exerted on polystyrene beads in the submillimeter range. A beam intensity less than 50  W/cm^{2} was required, ensuring damage-free trapping conditions. The large spectrum of frequencies covered by coherent ultrasonic sources provides a wide variety of manipulation possibilities from macroscopic to microscopic length scales. Our observations could open the way to important applications, in particular, in biology and biophysics at the cellular scale and for the design of acoustical machines in microfluidic environments.

HIFU Tissue Ablation: Concept and Devices

High intensity focused ultrasound (HIFU) is rapidly gaining clinical acceptance as a technique capable of providing non-invasive heating and ablation for a wide range of applications. Usually requiring only a single session, treatments are often conducted as day case procedures, with the patient either fully conscious, lightly sedated or under light general anesthesia. HIFU scores over other thermal ablation techniques because of the lack of necessity for the transcutaneous insertion of probes into the target tissue. Sources placed either outside the body (for treatment of tumors or abnormalities of the liver, kidney, breast, uterus, pancreas brain and bone), or in the rectum (for treatment of the prostate), provide rapid heating of a target tissue volume, the highly focused nature of the field leaving tissue in the ultrasound propagation path relatively unaffected. Numerous extra-corporeal, transrectal and interstitial devices have been designed to optimize application-specific treatment delivery for the wide-ranging areas of application that are now being explored with HIFU. Their principle of operation is described here, and an overview of their design principles is given.

MR-Guided Transcranial Focused Ultrasound

Previous chapters introduced the ability of using focused ultrasound to ablate tissues. It has led to various clinical applications in the treatment of uterine fibroid, prostate or liver cancers. Nevertheless, treating the brain non-invasively with focused ultrasound has been considered beyond reach for almost a century: The skull bone protects the brain from mechanical injuries, but it also reflects and refracts ultrasound, making it difficult to target the brain with focused ultrasound. Fortunately, aberration correction techniques have been developed recently and thermal lesioning in the thalamus has been achieved clinically. This chapter introduces the aberration effect of the skull bone and how it can be corrected non-invasively. It also presents the latest clinical results obtained with thermal ablation and introduces novel non-thermal approaches that could revolutionize brain therapy in the future.

Facilitation of Drug Transport across the Blood-Brain Barrier with Ultrasound and Microbubbles

Medical treatment options for central nervous system (CNS) diseases are limited due to the inability of most therapeutic agents to penetrate the blood–brain barrier (BBB). Although a variety of approaches have been investigated to open the BBB for facilitation of drug delivery, none has achieved clinical applicability. Mounting evidence suggests that ultrasound in combination with microbubbles might be useful for delivery of drugs to the brain through transient opening of the BBB. This technique offers a unique non-invasive avenue to deliver a wide range of drugs to the brain and promises to provide treatments for CNS disorders with the advantage of being able to target specific brain regions without unnecessary drug exposure. If this method could be applied for a range of different drugs, new CNS therapeutic strategies could emerge at an accelerated pace that is not currently possible in the field of drug discovery and development. This article reviews both the merits and potential risks of this new approach. It assesses methods used to verify disruption of the BBB with MRI and examines the results of studies aimed at elucidating the mechanisms of opening the BBB with ultrasound and microbubbles. Possible interactions of this novel delivery method with brain disease, as well as safety aspects of BBB disruption with ultrasound and microbubbles are addressed. Initial translational research for treatment of brain tumors and Alzheimer’s disease is presented.

Random calibration for accelerating MR-ARFI guided ultrasonic focusing in transcranial therapy

Transcranial focused ultrasound is a promising therapeutic modality. It consists of placing transducers around the skull and emitting shaped ultrasound waves that propagate through the skull and then concentrate on one particular location within the brain. However, the skull bone is known to distort the ultrasound beam. In order to compensate for such distortions, a number of techniques have been proposed recently, for instance using Magnetic Resonance Imaging feedback. In order to fully determine the focusing distortion due to the skull, such methods usually require as many calibration signals as transducers, resulting in a lengthy calibration process. In this paper, we investigate how the number of calibration sequences can be significantly reduced, based on random measurements and optimization techniques. Experimental data with six human skulls demonstrate that the number of measurements can be up to three times lower than with the standard methods, while restoring 90% of the focusing efficiency.

Improving the air coupling of bulk piezoelectric transducers with wedges of power-law profiles: a numerical study

An air-coupled ultrasonic transducer is created by bonding a bulk piezoelectric element onto the surface of a thick plate with a wedge of power-law profile. The wedge is used to improve the ultrasonic radiation efficiency. The power-law profile provides a smooth, impedance-matching transition for the mechanical energy to be transferred from the thick plate to the air, through the large-amplitude flexural waves observed in the thinnest region of the wedge. The performance of the proposed transducer is examined numerically and compared to that of a design where the piezoelectric element is isolated and where it is affixed to a thin plate of uniform thickness. The numerical analysis is first focused on the free-field radiation of the transducers. Then, time-reversal experiments are simulated by placing the transducers inside a cavity of arbitrary shape with some perfectly reflecting boundaries. In addition to time-reversal mirrors, the proposed concept could be integrated in the design of phased arrays and parametric arrays.

A comparison of methods for focusing the field of a HIFU array transducer through human ribs

A forward model, which predicts the scattering by human ribs of a multi-element high-intensity focused ultrasound transducer, was used to investigate the efficacy of a range of focusing approaches described in the literature. This forward model is based on the boundary element method and was described by Gélat et al (2011 Phys. Med. Biol. 56 5553-81; 2012 Phys. Med. Biol. 57 8471-97). The model has since been improved and features a complex surface impedance condition at the surface of the ribs. The inverse problem of focusing through the ribs was implemented on six transducer array-rib topologies and five methods of focusing were investigated, including spherical focusing, binarized apodization based on geometric ray tracing, phase conjugation and the decomposition of the time-reversal operator method. The excitation frequency was 1 MHz and the array was of spherical-section type. Both human and idealized rib topologies were considered. The merit of each method of focusing was examined. It was concluded that the constrained optimization approach offers greater potential than the other focusing methods in terms of maximizing the ratio of acoustic pressure magnitudes at the focus to those on the surface of the ribs whilst taking full advantage of the dynamic range of the phased array.

3-D transcranial ultrasound imaging with bilateral phase aberration correction of multiple isoplanatic patches: a pilot human study with microbubble contrast enhancement

With stroke currently the second-leading cause of death globally, and 87% of all strokes classified as ischemic, the development of a fast, accessible, cost-effective approach for imaging occlusive stroke could have a significant impact on health care outcomes and costs. Although clinical examination and standard computed tomography alone do not provide adequate information for understanding the complex temporal events that occur during an ischemic stroke, ultrasound imaging is well suited to the task of examining blood flow dynamics in real time and may allow for localization of a clot. A prototype bilateral 3-D ultrasound imaging system using two matrix array probes on either side of the head allows for correction of skull-induced aberration throughout two entire phased array imaging volumes. We investigated the feasibility of applying this custom correction technique in five healthy volunteers with Definity microbubble contrast enhancement. Subjects were scanned simultaneously via both temporal acoustic windows in 3-D color flow mode. The number of color flow voxels above a common threshold increased as a result of aberration correction in five of five subjects, with a mean increase of 33.9%. The percentage of large arteries visualized by 3-D color Doppler imaging increased from 46% without aberration correction to 60% with aberration correction.

Comparison and visualization of focusing wave fields from various time reversal techniques in elastic media

This Letter presents a series of time reversal experiments conducted on the surface of a fused silica glass block. Four different time reversal techniques are compared using three different imaging conditions. The techniques include two classical time reversal experiments: one with a pulse waveform source and one with an impulse response generated from a chirp signal. The other two techniques utilize the deconvolution, or inverse filter, signal processing methods for obtaining the signals to back propagate using a pulse waveform and an impulse response from a chirp. The max-in-time, symmetry, and energy current imaging conditions are compared.

Transcranial passive acoustic mapping with hemispherical sparse arrays using CT-based skull-specific aberration corrections: a simulation study

The feasibility of transcranial passive acoustic mapping with hemispherical sparse arrays (30 cm diameter, 16 to 1372 elements, 2.48 mm receiver diameter) using CT-based aberration corrections was investigated via numerical simulations. A multi-layered ray acoustic transcranial ultrasound propagation model based on CT-derived skull morphology was developed. By incorporating skull-specific aberration corrections into a conventional passive beamforming algorithm (Norton and Won 2000 IEEE Trans. Geosci. Remote Sens. 38 1337-43), simulated acoustic source fields representing the emissions from acoustically-stimulated microbubbles were spatially mapped through three digitized human skulls, with the transskull reconstructions closely matching the water-path control images. Image quality was quantified based on main lobe beamwidths, peak sidelobe ratio, and image signal-to-noise ratio. The effects on the resulting image quality of the source's emission frequency and location within the skull cavity, the array sparsity and element configuration, the receiver element sensitivity, and the specific skull morphology were all investigated. The system's resolution capabilities were also estimated for various degrees of array sparsity. Passive imaging of acoustic sources through an intact skull was shown possible with sparse hemispherical imaging arrays. This technique may be useful for the monitoring and control of transcranial focused ultrasound (FUS) treatments, particularly non-thermal, cavitation-mediated applications such as FUS-induced blood-brain barrier disruption or sonothrombolysis, for which no real-time monitoring techniques currently exist.

Simultaneous bilateral real-time 3-d transcranial ultrasound imaging at 1 MHz through poor acoustic windows

Ultrasound imaging has been proposed as a rapid, portable alternative imaging modality to examine stroke patients in pre-hospital or emergency room settings. However, in performing transcranial ultrasound examinations, 8%-29% of patients in a general population may present with window failure, in which case it is not possible to acquire clinically useful sonographic information through the temporal bone acoustic window. In this work, we describe the technical considerations, design and fabrication of low-frequency (1.2 MHz), large aperture (25.3 mm) sparse matrix array transducers for 3-D imaging in the event of window failure. These transducers are integrated into a system for real-time 3-D bilateral transcranial imaging-the ultrasound brain helmet-and color flow imaging capabilities at 1.2 MHz are directly compared with arrays operating at 1.8 MHz in a flow phantom with attenuation comparable to the in vivo case. Contrast-enhanced imaging allowed visualization of arteries of the Circle of Willis in 5 of 5 subjects and 8 of 10 sides of the head despite probe placement outside of the acoustic window. Results suggest that this type of transducer may allow acquisition of useful images either in individuals with poor windows or outside of the temporal acoustic window in the field.

Pitch-catch phase aberration correction of multiple isoplanatic patches for 3-D transcranial ultrasound imaging

Having previously presented the ultrasound brain helmet, a system for simultaneous 3-D ultrasound imaging via both temporal bone acoustic windows, the scanning geometry of this system is utilized to allow each matrix array to serve as a correction source for the opposing array. Aberration is estimated using cross-correlation of RF channel signals, followed by least mean squares solution of the resulting overdetermined system. Delay maps are updated and real-time 3-D scanning resumes. A first attempt is made at using multiple arrival time maps to correct multiple unique aberrators within a single transcranial imaging volume, i.e., several isoplanatic patches. This adaptive imaging technique, which uses steered unfocused waves transmitted by the opposing, or beacon, array, updates the transmit and receive delays of 5 isoplanatic patches within a 64° x 64° volume. In phantom experiments, color flow voxels above a common threshold have also increased by an average of 92%, whereas color flow variance decreased by an average of 10%. This approach has been applied to both temporal acoustic windows of two human subjects, yielding increases in echo brightness in 5 isoplanatic patches with a mean value of 24.3 ± 9.1%, suggesting that such a technique may be beneficial in the future for performing noninvasive 3-D color flow imaging of cerebrovascular disease, including stroke.

Timing-error-difference calibration of a two-dimensional array imaging system using the overlapping-subaperture algorithm

Timing errors in the transmitting and receiving electronic channels of an imaging system can generate different transmission and reception phase-aberration profiles. To decide if these two profiles need to be measured separately, an overlapping-subaperture algorithm has been proposed in a previous paper to measure the difference between timing errors in transmitting and receiving channels connected to each element in a two-dimensional array. This algorithm has been used to calibrate a custom built imaging system with a curved linear two-dimensional array, and the results are presented in this paper. The experimental results have demonstrated that the overlapping-subaperture algorithm is capable of calibrating the timing-error-difference profile of this imaging system with a standard deviation of only a few nanoseconds. Experimental results have also shown that the time-error-difference profile of this imaging system is smaller than one tenth of a wavelength and there is no need to measure the transmission and reception phase-aberration profiles separately. The derived average phase-aberration profile using the near-field signal-redundancy algorithm can be used to correct phase aberrations for both transmission and reception.

Correction of tissue-motion effects on common-midpoint signals using reciprocal signals

The near field signal redundancy algorithm for phase-aberration correction is sensitive to tissue motion because several separated transmissions are usually needed to acquire a set of common-midpoint signals. If tissues are moving significantly due to, for example, heart beats, the effects of tissue motion on common-midpoint signals need to be corrected before the phase-aberration profile can be successfully measured. Theoretical analyses in this paper show that the arrival-time difference between a pair of common-midpoint signals due to tissue motion is usually very similar to that between the pair of reciprocal signals acquired using the same two transmissions. Based on this conclusion, an algorithm for correcting tissue-motion effects on the peak position of cross-correlation functions between common-midpoint signals is proposed and initial experimental results are also presented.

Optimal spatiotemporal focusing through complex scattering media

We present an alternative approach for spatiotemporal focusing through complex scattering media by wave front shaping. Using a nonlinear feedback signal to shape the incident pulsed wave front, we show that the limit of a spatiotemporal matched filter can be achieved; i.e., the wave amplitude at the intended time and focus position is maximized for a given input energy. It is exactly what is also achieved with time reversal. Demonstrated with ultrasound experiments, our method is generally applicable to all types of waves.

A passive inverse filter for Green's function retrieval

Passive methods for the recovery of Green's functions from ambient noise require strong hypotheses, including isotropic distribution of the noise sources. Very often, this distribution is nonisotropic, which introduces bias in the Green's function reconstruction. To minimize this bias, a spatiotemporal inverse filter is proposed. The method is tested on a directive noise field computed from an experimental active seismic data set. The results indicate that the passive inverse filter allows the manipulation of the spatiotemporal degrees of freedom of a complex wave field, and it can efficiently compensate for the noise wavefield directivity.

Filter-based compounded delay estimation with application to strain imaging

Ultrasonic wave interference produces local fluctuations in both the envelope, known as speckle, and phase of echoes. Furthermore, such fluctuations are correlated in space, and subsequent motion estimation from the envelope and/or phase signal produces patterned, correlated errors. Compounding, or combining information from multiple decorrelated looks, reduces such effects. We propose using a filter bank to create multiple looks to produce a compounded motion estimate. In particular, filtering in the lateral direction is shown to preserve delay estimation accuracy in the filtered sub-bands while creating decorrelation between sub-bands at the expense of some lateral resolution. For Gaussian apodization, we explicitly compute the induced signal decorrelation produced by Gabor filters. Furthermore, it is shown that lateral filtering is approximately equivalent to steering, in which filtered sub-bands correspond to signals extracted from shifted sub-apertures. Field II simulation of a point spread function verifies this claim. We use phase zero and its variants as displacement estimators for our compounded result. A simplified deformation model is used to provide computer simulations of deforming an elastic phantom. Simulations demonstrate root mean square error (RMSE) reduction in both displacement and strain of the compounded result over conventional and its laterally blurred versions. Then we apply the methods to experimental data using a commercial elastic phantom, demonstrating an improvement in strain SNR.

Experimental demonstration of the utility of pressure sensitivity kernels in time-reversal

Pressure sensitivity kernels were recently applied to time-reversal acoustics in an attempt to explain the enhanced stability of the time-reversal focal spot [Raghukumar et al., J. Acoust. Soc. Am. 124, 98-112 (2008)]. The theoretical framework developed was also used to derive optimized source functions, closely related to the inverse filter. The use of these optimized source functions results in an inverse filter-like focal spot which is more robust to medium sound speed fluctuations than both time-reversal and the inverse filter. In this paper the theory is applied to experimental data gathered during the Focused Acoustic Fields experiment, conducted in 2005, north of Elba Island in Italy. Sensitivity kernels are calculated using a range-independent sound-speed profile, for a geometry identical to that used in the experiment, and path sensitivities are identified with observed arrivals. The validity of the kernels in tracking time-evolving Green's functions is studied, along with limitations that result from a linearized analysis. An internal wave model is used to generate an ensemble of sound speed profiles, which are then used along with the calculated sensitivity kernels to derive optimized source functions. Focal spots obtained using the observed Green's functions with these optimized source functions are then compared to those obtained using time-reversal and the inverse-filter. It is shown that these functions are able to provide a focal spot superior to time-reversal while being more robust to sound speed fluctuations than the inverse filter or time-reversal.

A method for MRI guidance of intercostal high intensity focused ultrasound ablation in the liver

PURPOSE

High intensity focused ultrasound (HIFU) is a promising method for the noninvasive treatment of liver tumors. However, the presence of ribs in the HIFU beam path remains problematic since it may lead to adverse effects (skin burns) by absorption and reflection of the incident beam at or near the bone surface. This article presents a method based on magnetic resonance (MR) imaging for identification of the ribs in the HIFU beam, and for selection of the transducer elements to deactivate.

METHODS

The ribs are visualized on anatomical images acquired prior to heating and manually segmented. The resulting regions of interest surrounding the ribs are projected onto the transducer surface by ray tracing from the focal point. The transducer elements in the "shadow" of the ribs are then deactivated. The method was validated ex vivo and in vivo in pig liver during breathing under multislice real-time MR thermometry, using the proton resonance frequency shift method.

RESULTS

Ex vivo and in vivo temperature data showed that the temperature increase near the ribs was substantial when HIFU sonications were performed with all elements active, whereas the temperature was reduced with deactivation of the transducer elements located in front of the ribs. The temperature at the focal point was similar with and without deactivation of the transducer elements, indicative of no loss of heat efficiency with the proposed technique.

CONCLUSIONS

This method is simple, rapid, and reliable, and enables intercostal HIFU ablation while sparing ribs and their surrounding tissues.

MR-guided adaptive focusing of ultrasound

Adaptive focusing of ultrasonic waves under the guidance of a magnetic resonance (MR) system is demonstrated for medical applications. This technique is based on the maximization of the ultrasonic wave intensity at one targeted point in space. The wave intensity is indirectly estimated from the local tissue displacement induced at the chosen focus by the acoustic radiation force of ultrasonic beams. Coded ultrasonic waves are transmitted by an ultrasonic array and an MRI scanner is used to measure the resulting local displacements through a motion-sensitive MR sequence. After the transmission of a set of spatially encoded ultrasonic waves, a non-iterative inversion process is employed to accurately estimate the spatial-temporal aberration induced by the propagation medium and to maximize the acoustical intensity at the target.Both programmable and physical aberrating layers introducing strong distortions (up to 2pi radians) were recovered within acceptable errors (<0.8 rad). This noninvasive technique is shown to accurately correct phase aberrations in a phantom gel with negligible heat deposition and limited acquisition time. These refocusing performances demonstrate a major potential in the field of MR-guided ultrasound therapy in particular for transcranial brain high-intensity focused ultrasound.

Multiple-array passive acoustic source localization in urban environments

In many situations of interest, obstacles to acoustic wave propagation such as terrain or buildings exist that provide unique challenges to localization. These obstacles introduce multiple propagation paths, reflections, and diffraction into the propagation. In this paper, matched field processing is proposed as an effective method of acoustic localization in a two dimensional scattering environment. Numerical techniques can be used to model complex propagation in a space where analytical solutions are not feasible. Realistically, there is always some uncertainty in model parameters that in turn can adversely affect localization ability. In particular, uncertainty in array location, sound speed, and various parameters affecting inter-array coherence only are investigated. A spatially distributed, multiarray network is shown to mitigate the effects of uncertainty. Multiarray inverse filter processing techniques are evaluated through perturbation of uncertain model parameters. These techniques are more accurate and flexible to implement than other matched field processing methods such as time reversal.

Non-destructive imaging using the time domain topological energy method

The time domain topological gradient has recently been presented as a promising technique for imaging of complex structures [N. Dominguez, et al., Wave Motion 42 (1) (2005) 31-52], showing results obtained on simulated data. The method relies on the adequate combination of two computed ultrasonic fields, one forward and one adjoint. The adjoint field, that carries the information about the defects, is analogous to a time reversal operation. The forward field plays the role of a "photographic developer". The great interest of the method is that the time reversal operation is not done experimentally but numerically, allowing implementation of low cost systems. The expected low sensitivity of the method to the presence of measurement or structure noise had been discussed in a former work. This paper shows results obtained on experimental data in the context of non-destructive testing, showing the efficiency of the method even in complex configuration such that composite material testing.

Energy-based adaptive focusing of waves: application to noninvasive aberration correction of ultrasonic wavefields

An aberration correction method based on the maximization of the wave intensity at the focus of an emitting array is presented. The potential of this new adaptive focusing technique is investigated for ultrasonic focusing in biological tissues. The acoustic intensity is maximized noninvasively through direct measurement or indirect estimation of the beam energy at the focus for a series of spatially coded emissions. For ultrasonic waves, the acoustic energy at the desired focus can be indirectly estimated from the local displacements induced in tissues by the ultrasonic radiation force of the beam. Based on the measurement of these displacements, this method allows determination of the precise estimation of the phase and amplitude aberrations, and consequently the correction of aberrations along the beam travel path. The proof of concept is first performed experimentally using a large therapeutic array with strong electronic phase aberrations (up to 2pi). Displacements induced by the ultrasonic radiation force at the desired focus are indirectly estimated using the time shift of backscattered echoes recorded on the array. The phase estimation is deduced accurately using a direct inversion algorithm which reduces the standard deviation of the phase distribution from sigma = 1.89 radian before correction to sigma = 0.53 radian following correction. The corrected beam focusing quality is verified using a needle hydrophone. The peak intensity obtained through the aberrator is found to be -7.69 dB below the reference intensity obtained without any aberration. Using the phase correction, a sharp focus is restored through the aberrator with a relative peak intensity of -0.89 dB. The technique is tested experimentally using a linear transmit/receive array through a real aberrating layer. The array is used to automatically correct its beam quality, as it both generates the radiation force with coded excitations and indirectly estimates the acoustic intensity at the focus with speckle tracking. This technique could have important implications in the field of high-intensity focused ultrasound even in complex configurations such as transcranial, transcostal, or deep seated organs.

Manipulating spatiotemporal degrees of freedom of waves in random media

We show that all the spatiotemporal degrees of freedom available in a complex medium can be harnessed and converted into spatial ones. This is demonstrated experimentally through an instantaneous spatial inversion, using broadband ultrasonic waves in a multiple scattering sample. We show theoretically that the inversion convergence is governed by the total number of degrees of freedom available in the medium for a fixed bandwidth and demonstrate experimentally its use for complex media investigation. We believe our approach has potential in sensing, imagery, focusing, and telecommunication.

Non-invasive transcranial ultrasound therapy based on a 3D CT scan: protocol validation and in vitro results

A non-invasive protocol for transcranial brain tissue ablation with ultrasound is studied and validated in vitro. The skull induces strong aberrations both in phase and in amplitude, resulting in a severe degradation of the beam shape. Adaptive corrections of the distortions induced by the skull bone are performed using a previous 3D computational tomography scan acquisition (CT) of the skull bone structure. These CT scan data are used as entry parameters in a FDTD (finite differences time domain) simulation of the full wave propagation equation. A numerical computation is used to deduce the impulse response relating the targeted location and the ultrasound therapeutic array, thus providing a virtual time-reversal mirror. This impulse response is then time-reversed and transmitted experimentally by a therapeutic array positioned exactly in the same referential frame as the one used during CT scan acquisitions. In vitro experiments are conducted on monkey and human skull specimens using an array of 300 transmit elements working at a central frequency of 1 MHz. These experiments show a precise refocusing of the ultrasonic beam at the targeted location with a positioning error lower than 0.7 mm. The complete validation of this transcranial adaptive focusing procedure paves the way to in vivo animal and human transcranial HIFU investigations.

A phase and space coherent direct imaging method

A direct imaging algorithm for point and extended targets is presented. The algorithm is based on a physical factorization of the response matrix of a transducer array. The factorization is used to transform a passive target problem to an active source problem and to extract principal components (tones) in a phase consistent way. The multitone imaging function can superpose multiple tones (spatial diversity/aperture of the array) and frequencies (bandwidth of the probing signal) based on phase coherence. The method is a direct imaging algorithm that is simple and efficient since no forward solver or iteration is needed. Robustness of the algorithm with respect to noise is demonstrated via numerical examples.

Shear elasticity estimation from surface wave: the time reversal approach

In this work the shear elasticity of soft solids is measured from the surface wave speed estimation. An external source creates mechanical waves which are detected using acoustic sensors. The surface wave speed estimation is extracted from the complex reverberated elastic field through a time-reversal analysis. Measurements in a hard and a soft gelatin-based phantom are validated by independent transient elastography estimations. In contrast with other elasticity assessment methods, one advantage of the present approach is its low sound technology cost. Experiments performed in cheese and soft phantoms allows one to envision applications in the food industry and medicine.