Adaptive focusing for transcranial ultrasound imaging using dual arrays

Device for Multifocal Delivery of Ultrasound Into Deep Brain Regions in Humans

Low-intensity focused ultrasound provides the means to noninvasively stimulate or release drugs in specified deep brain targets. However, successful clinical translations require hardware that maximizes acoustic transmission through the skull, enables flexible electronic steering, and provides accurate and reproducible targeting while minimizing the use of MRI. We have developed a device that addresses these practical requirements. The device delivers ultrasound through the temporal and parietal skull windows, which minimize the attenuation and distortions of the ultrasound by the skull. The device consists of 252 independently controlled elements, which provides the ability to modulate multiple deep brain targets at a high spatiotemporal resolution, without the need to move the device or the subject. And finally, the device uses a mechanical registration method that enables accurate deep brain targeting both inside and outside of the MRI. Using this method, a single MRI scan is necessary for accurate targeting; repeated subsequent treatments can be performed reproducibly in an MRI-free manner. We validated these functions by transiently modulating specific deep brain regions in two patients with treatment-resistant depression.

Transcranial Phase Correction Using Pulse-Echo Ultrasound and Deep Learning: A 2-D Numerical Study

Phase aberration caused by human skulls severely degrades the quality of transcranial ultrasound images, posing a major challenge in the practical application of transcranial ultrasound techniques in adults. Aberration can be corrected if the skull profile (i.e., thickness distribution) and speed of sound (SOS) are known. However, accurately estimating the skull profile and SOS using ultrasound with a physics-based approach is challenging due to the complexity of the interaction between ultrasound and the skull. A deep learning approach is proposed herein to estimate the skull profile and SOS using ultrasound radiofrequency (RF) signals backscattered from the skull. A numerical study was performed to test the approach's feasibility. Realistic numerical skull models were constructed from computed tomography (CT) scans of five ex vivo human skulls in this numerical study. Acoustic simulations were performed on 3595 skull segments to generate array-based ultrasound backscattered signals. A deep learning model was developed and trained to estimate skull thickness and SOS from RF channel data. The trained model was shown to be highly accurate. The mean absolute error (MAE) was 0.15 mm (2% error) for thickness estimation and 13 m/s (0.5% error) for SOS estimation. The Pearson correlation coefficient between the estimated and ground-truth values was 0.99 for thickness and 0.95 for SOS. Aberration correction performed using deep-learning-estimated skull thickness and SOS values yielded significantly improved beam focusing (e.g., narrower beams) and transcranial imaging quality (e.g., improved spatial resolution and reduced artifacts) compared with no aberration correction. The results demonstrate the feasibility of the proposed approach for transcranial phase aberration correction.

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.

Dual-Probe Transcranial Full-Waveform Inversion: A Brain Phantom Feasibility Study

OBJECTIVE

Despite being a low-cost, portable and safe medical imaging technique, transcranial ultrasound imaging is not used widely in adults because of the severe degradation and distortion of signals caused by the skull. Full-waveform inversion (FWI) has recently been found to have potential as an effective method for transcranial ultrasound tomography to obtain high-quality, subwavelength-resolution acoustic models of the brain using low-frequency ultrasound data. In this study is the first demonstration of this method in recovering a high-resolution 2-D reconstruction of a brain and skull ultrasound imaging phantom using experimentally acquired data.

METHODS

A 2:5 scale brain phantom encased within a 3-D-printed skull-mimicking layer was created to simulate a clinical transcranial imaging target. To obtain tomographic ultrasound data on the brain and skull phantom, a tomographic ultrasound acquisition system was designed and implemented using commercially available low-frequency cardiac probes. FWI reconstructions of the brain and skull phantom were performed using the acquired tomographic data and were compared with corresponding synthetic reconstructions. This comparison was used to evaluate the feasibility of the proposed imaging system when employing different transducer array configurations.

RESULTS

We demonstrate the successful FWI reconstruction of the brain phantom within the skull mimic from experimentally acquired tomographic ultrasound data. To mitigate the effects of the skull-mimicking material, a reflection-matching algorithm was applied to model the morphology of the skull layer prior to performing the inversion.

CONCLUSION

The findings of this study provide a promising step toward the clinical use of FWI for transcranial ultrasound imaging in adults.

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

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

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

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

Mathematical Model of Ultrasound Attenuation With Skull Thickness for Transcranial-Focused Ultrasound

Transcranial-focused ultrasound (tFUS) has potential for both neuromodulation and neuroimaging. Due to the influence of head tissue, especially the skull, its attenuation is a key issue affecting precise focusing. The objective of the present study was to construct a mathematical model of ultrasound attenuation inclusive of skull thickness. First, combined with real skull phantom experiments and simulation experiments, tFUS attenuation of different head tissues was investigated. Furthermore, based on the system identification method, a mathematical model of ultrasound attenuation was constructed taking skull thickness into account. Finally, the performance of the mathematical model was tested, and its potential applications were investigated. For different head tissues, including scalp, skull, and brain tissue, the skull was found to be the biggest influencing factor for ultrasound attenuation, the attenuation caused by it being 4.70 times and 7.06 times that of attenuation caused by the brain and scalp, respectively. Consistent with the results of both the simulation and phantom experiments, the attenuation of the mathematical model increased as the skull thickness increased. The average error of the mathematical model was 1.87% in the phantom experiment. In addition, the experimental results show that the devised mathematical model is suitable for different initial pressures and different skulls with correlation coefficients higher than 0.99. Both simulation and phantom experiments validated the effectiveness of the proposed mathematical model. It can be concluded from this experiment that the proposed mathematical model can accurately calculate the tFUS attenuation and can significantly contribute to further research and application of tFUS.

Full-Matrix Phase Shift Migration Method for Transcranial Ultrasonic Imaging

A spectrum-domain method, called full-matrix phase shift migration (FM-PSM), is presented for transcranial ultrasound phase correction and imaging with ideal synthetic aperture focusing technology. The simulated data obtained using the pseudospectral time-domain method are used to evaluate the feasibility of the method. The experimental data measured from a 3-D printed skull phantom are used to evaluate the algorithm performance in terms of resolution, contrast-to-noise ratio (CNR), and eccentricity comparing with the classical ray-tracing delay and sum (DAS) method. In wire imaging experiment, FM-PSM has a lateral resolution of 0.22 mm and ray-tracing DAS has a lateral resolution of 0.24 mm measured at -6-dB drop using a transducer with a center frequency of 6.25 MHz. In cylinder imaging experiment, FM-PSM has a CNR of 2.14 and ray-tracing DAS has a CNR of 1.82, which illustrates about 17% improvement. For a J -element array and an output image with pixels M ×N (lateral × axial), the computational cost of the DAS is of O(J ×M²×N²) ; on the contrary, it reduces to O(J ×M ×N²) with the proposed FM-PSM. The results suggest that FM-PSM is an efficiency method for transcranial ultrasonic imaging.

Ultrasound for the Brain: A Review of Physical and Engineering Principles, and Clinical Applications

The emergence of new ultrasound technologies has improved our understanding of the brain functions and offered new opportunities for the treatment of brain diseases. Ultrasound has become a valuable tool in preclinical animal and clinical studies as it not only provides information about the structure and function of brain tissues but can also be used as a therapy alternative for brain diseases. High-resolution cerebral flow images with high sensitivity can be acquired using novel functional ultrasound and super-resolution ultrasound imaging techniques. The noninvasive treatment of essential tremors has been clinically approved and it has been demonstrated that the ultrasound technology can revolutionize the currently existing treatment methods. Microbubble-mediated ultrasound can remotely open the blood-brain barrier enabling targeted drug delivery in the brain. More recently, ultrasound neuromodulation received a great amount of attention due to its noninvasive and deep penetration features and potential therapeutic benefits. This review provides a thorough introduction to the current state-of-the-art research on brain ultrasound and also introduces basic knowledge of brain ultrasound including the acoustic properties of the brain/skull and engineering techniques for ultrasound. Ultrasound is expected to play an increasingly important role in the diagnosis and therapy of brain diseases.

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

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

The use of transcranial ultrasound and clinical assessment to diagnose ischaemic stroke due to large vessel occlusion in remote and rural areas

Rapid endovascular thrombectomy, which can only be delivered in specialist centres, is the most effective treatment for acute ischaemic stroke due to large vessel occlusion (LVO). Pre-hospital selection of these patients is challenging, especially in remote and rural areas due to long transport times and limited access to specialist clinicians and diagnostic facilities. We investigated whether combined transcranial ultrasound and clinical assessment (“TUCA” model) could accurately triage these patients and improve access to thrombectomy. We recruited consecutive patients within 72 hours of suspected stroke, and performed non-contrast transcranial colour-coded ultrasonography within 24 hours of brain computed tomography. We retrospectively collected clinical information, and used hospital discharge diagnosis as the “gold standard”. We used binary regression for diagnosis of haemorrhagic stroke, and an ordinal regression model for acute ischaemic stroke with probable LVO, without LVO, transient ischaemic attacks (TIA) and stroke mimics. We calculated sensitivity, specificity, positive and negative predictive values and performed a sensitivity analysis. We recruited 107 patients with suspected stroke from July 2017 to December 2019 at two study sites: 13/107 (12%) with probable LVO, 50/107 (47%) with acute ischaemic stroke without LVO, 18/107 (17%) with haemorrhagic stroke, and 26/107 (24%) with stroke mimics or TIA. The model identified 55% of cases with probable LVO who would have correctly been selected for thrombectomy and 97% of cases who would not have required this treatment (sensitivity 55%, specificity 97%, positive and negative predictive values 75% and 93%, respectively). Diagnostic accuracy of the proposed model was superior to the clinical assessment alone. These data suggest that our model might be a useful tool to identify pre-hospital patients requiring mechanical thrombectomy, however a larger sample is required with the use of CT angiogram as a reference test.

A novel transcranial ultrasound imaging method with diverging wave transmission and deep learning approach

Real time brain transcranial ultrasound imaging is extremely intriguing because of its numerous applications. However, the skull causes phase distortion and amplitude attenuation of ultrasound signals due to its density: the speed of sound is significantly different in bone tissue than in soft tissue. In this study, we propose an ultrafast transcranial ultrasound imaging technique with diverging wave (DW) transmission and a deep learning approach to achieve large field-of-view with high resolution and real time brain ultrasound imaging. DW transmission provides a frame rate of several kiloHz and a large field of view that is suitable for human brain imaging via a small acoustic window. However, it suffers from poor image quality because the diverging waves are all unfocused. Here, we adopted adaptive beamforming algorithms to improve both the image contrast and the lateral resolution. Both simulated and in situ experiments with a human skull resulted in significant image improvements. However, the skull still introduces a wavefront offset and distortion, which degrades the image quality even when adaptive beamforming methods are used. Thus, we also employed a U-Net neural network to detect the contour and position of the skull directly from the acquired RF signal matrix. This approach avoids the need for beamforming, image reconstruction, and image segmentation, making it more suitable for clinical use.

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.

Catheter Hydrophone Aberration Correction for Transcranial Histotripsy Treatment of Intracerebral Hemorrhage: Proof-of-Concept

Histotripsy is a minimally invasive ultrasound therapy that has shown rapid liquefaction of blood clots through human skullcaps in an in vitro intracerebral hemorrhage model. However, the efficiency of these treatments can be compromised if the skull-induced aberrations are uncorrected. We have developed a catheter hydrophone which can perform aberration correction (AC) and drain the liquefied clot following histotripsy treatment. Histotripsy pulses were delivered through an excised human skullcap using a 256-element, 500-kHz hemisphere array transducer with a 15-cm focal distance. A custom hydrophone was fabricated using a mm PZT-5h crystal interfaced to a coaxial cable and integrated into a drainage catheter. An AC algorithm was developed to correct the aberrations introduced between histotripsy pulses from each array element. An increase in focal pressure of up to 60% was achieved at the geometric focus and 27%-62% across a range of electronic steering locations. The sagittal and axial -6-dB beam widths decreased from 4.6 to 2.2 mm in the sagittal direction and 8 to 4.4 mm in the axial direction, compared to 1.5 and 3 mm in the absence of aberration. After performing AC, lesions with diameters ranging from 0.24 to 1.35 mm were generated using electronic steering over a mm grid in a tissue-mimicking phantom. An average volume of 4.07 ± 0.91 mL was liquefied and drained after using electronic steering to treat a 4.2-mL spherical volume in in vitro bovine clots through the skullcap.

Image-guided ultrasound phased arrays are a disruptive technology for non-invasive therapy

Focused ultrasound offers a non-invasive way of depositing acoustic energy deep into the body, which can be harnessed for a broad spectrum of therapeutic purposes, including tissue ablation, the targeting of therapeutic agents, and stem cell delivery. Phased array transducers enable electronic control over the beam geometry and direction, and can be tailored to provide optimal energy deposition patterns for a given therapeutic application. Their use in combination with modern medical imaging for therapy guidance allows precise targeting, online monitoring, and post-treatment evaluation of the ultrasound-mediated bioeffects. In the past there have been some technical obstacles hindering the construction of large aperture, high-power, densely-populated phased arrays and, as a result, they have not been fully exploited for therapy delivery to date. However, recent research has made the construction of such arrays feasible, and it is expected that their continued development will both greatly improve the safety and efficacy of existing ultrasound therapies as well as enable treatments that are not currently possible with existing technology. This review will summarize the basic principles, current statures, and future potential of image-guided ultrasound phased arrays for therapy.

Computational exploration of wave propagation and heating from transcranial focused ultrasound for neuromodulation

OBJECTIVE

While ultrasound is largely established for use in diagnostic imaging, its application for neuromodulation is relatively new and crudely understood. The objective of the present study was to investigate the effects of tissue properties and geometry on the wave propagation and heating in the context of transcranial neuromodulation.

APPROACH

A computational model of transcranial-focused ultrasound was constructed and validated against empirical data. The models were then incrementally extended to investigate a number of issues related to the use of ultrasound for neuromodulation, including the effect on wave propagation of variations in geometry of skull and gyral anatomy as well as the effect of multiple tissue and media layers, including scalp, skull, CSF, and gray/white matter. In addition, a sensitivity analysis was run to characterize the influence of acoustic properties of intracranial tissues. Finally, the heating associated with ultrasonic stimulation waveforms designed for neuromodulation was modeled.

MAIN RESULTS

The wave propagation of a transcranially focused ultrasound beam is significantly influenced by the cranial domain. The half maximum acoustic beam intensity profiles are insensitive overall to small changes in material properties, though the inclusion of sulci in models results in greater peak intensity values compared to a model without sulci (1%-30% greater). Finally, heating using currently employed stimulation parameters in humans is highest in bone (0.16 °C) and is negligible in brain (4.27 × 10(-3) °C) for a 0.5 s exposure.

SIGNIFICANCE

Ultrasound for noninvasive neuromodulation holds great promise and appeal for its non-invasiveness, high spatial resolution and deep focal lengths. Here we show gross brain anatomy and biological material properties to have limited effect on ultrasound wave propagation and to result in safe heating levels in the skull and brain.

Investigation of the correlation between diffuse infrared and ultrasound for transcranial ultrasound

Over the past two decades the feasibility for using transcranial ultrasound as both a therapeutic and diagnostic tool has been established. Various aberration-correction techniques have been proposed to achieve transcranial focusing, including CT-derived model based corrections, ultrasound-derived model based corrections, magnetic resonance acoustic radiation force (MR-ARFI) techniques, and techniques involving the invasive introduction of an acoustic source or receiver into the brain. Here, we investigate the correlation between transcranial infrared light (IR) and transcranial ultrasound, where we examine whether IR could be an indicator of any of the key acoustic properties that affect transcranial transmission (signal attenuation, speed of sound, and bone density). Nine human skull samples were utilized in the study. The interior of each sample was illuminated over its inner surface using a diffuse light source. Light transmitted to the outer surface was detected by a 3-mm diameter 940-nm infrared sensor. Acoustic measurements were likewise obtained in a water tank using a 12.7-mm diameter 1-MHz source and a needle hydrophone receiver. Results reveal a positive correlation between the acoustic time-of-flight and optical intensity (the correlation coefficient is between 0.5 and 0.9). Subsequent investigation shows this correlation to hold independent of the presence or absence of dura mater on the samples. Poor correlation is observed between acoustic amplitude and optical intensity (the correlation coefficient is between 0.1 and 0.7).

Comparison of analytical and numerical approaches for CT-based aberration correction in transcranial passive acoustic imaging

Computed tomography (CT)-based aberration corrections are employed in transcranial ultrasound both for therapy and imaging. In this study, analytical and numerical approaches for calculating aberration corrections based on CT data were compared, with a particular focus on their application to transcranial passive imaging. Two models were investigated: a three-dimensional full-wave numerical model (Connor and Hynynen 2004 IEEE Trans. Biomed. Eng. 51 1693-706) based on the Westervelt equation, and an analytical method (Clement and Hynynen 2002 Ultrasound Med. Biol. 28 617-24) similar to that currently employed by commercial brain therapy systems. Trans-skull time delay corrections calculated from each model were applied to data acquired by a sparse hemispherical (30 cm diameter) receiver array (128 piezoceramic discs: 2.5 mm diameter, 612 kHz center frequency) passively listening through ex vivo human skullcaps (n  =  4) to emissions from a narrow-band, fixed source emitter (1 mm diameter, 516 kHz center frequency). Measurements were taken at various locations within the cranial cavity by moving the source around the field using a three-axis positioning system. Images generated through passive beamforming using CT-based skull corrections were compared with those obtained through an invasive source-based approach, as well as images formed without skull corrections, using the main lobe volume, positional shift, peak sidelobe ratio, and image signal-to-noise ratio as metrics for image quality. For each CT-based model, corrections achieved by allowing for heterogeneous skull acoustical parameters in simulation outperformed the corresponding case where homogeneous parameters were assumed. Of the CT-based methods investigated, the full-wave model provided the best imaging results at the cost of computational complexity. These results highlight the importance of accurately modeling trans-skull propagation when calculating CT-based aberration corrections. Although presented in an imaging context, our results may also be applicable to the problem of transmit focusing through the skull.

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.

Analysis and measurement of the modulation transfer function of harmonic shear wave induced phase encoding imaging

Shear wave induced phase encoding (SWIPE) imaging generates ultrasound backscatter images of tissue-like elastic materials by using traveling shear waves to encode the lateral position of the scatters in the phase of the received echo. In contrast to conventional ultrasound B-scan imaging, SWIPE offers the potential advantages of image formation without beam focusing or steering from a single transducer element, lateral resolution independent of aperture size, and the potential to achieve relatively high lateral resolution with low frequency ultrasound. Here a Fourier series description of the phase modulated echo signal is developed, demonstrating that echo harmonics at multiples of the shear wave frequency reveal target k-space data at identical multiples of the shear wavenumber. Modulation transfer functions of SWIPE imaging systems are calculated for maximum shear wave acceleration and maximum shear constraints, and compared with a conventionally focused aperture. The relative signal-to-noise ratio of the SWIPE method versus a conventionally focused aperture is found through these calculations. Reconstructions of wire targets in a gelatin phantom using 1 and 3.5 MHz ultrasound and a cylindrical shear wave source are presented, generated from the fundamental and second harmonic of the shear wave modulation frequency, demonstrating weak dependence of lateral resolution with ultrasound frequency.

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.

Refraction correction in 3D transcranial ultrasound imaging

We present the first correction of refraction in three-dimensional (3D) ultrasound imaging using an iterative approach that traces propagation paths through a two-layer planar tissue model, applying Snell's law in 3D. This approach is applied to real-time 3D transcranial ultrasound imaging by precomputing delays offline for several skull thicknesses, allowing the user to switch between three sets of delays for phased array imaging at the push of a button. Simulations indicate that refraction correction may be expected to increase sensitivity, reduce beam steering errors, and partially restore lost spatial resolution, with the greatest improvements occurring at the largest steering angles. Distorted images of cylindrical lesions were created by imaging through an acrylic plate in a tissue-mimicking phantom. As a result of correcting for refraction, lesions were restored to 93.6% of their original diameter in the lateral direction and 98.1% of their original shape along the long axis of the cylinders. In imaging two healthy volunteers, the mean brightness increased by 8.3% and showed no spatial dependency.

A review of numerical and experimental compensation techniques for skull-induced phase aberrations in transcranial focused ultrasound

The development of phased array transducers and their integration with magnetic resonance (MR) guidance and thermal monitoring has established transcranial MR-guided focused ultrasound (tcMRgFUS) as an attractive non-invasive modality for neurosurgical interventions. The presence of the skull, however, compromises the efficiency of transcranial FUS (tcFUS) therapy, as its heterogeneous nature and acoustic characteristics induce significant phase aberrations and energy attenuation, especially at the higher acoustic frequencies employed in tcFUS thermal therapy. These aberrations may distort and shift the acoustic focus as well as induce heating at the patient's scalp and skull bone. Phased array transducers feature hundreds of elements that can be driven individually, each with its own phase and amplitude. This feature allows for compensation of skull-induced aberrations by calculation and application of appropriate phase and amplitude corrections. In this paper, we illustrate the importance of precise refocusing and provide a comprehensive review of the wide variety of numerical and experimental techniques that have been used to estimate these corrections.

Transcranial ultrasound imaging with speed of sound-based phase correction: a numerical study

This paper presents a numerical study for ultrasound transcranial imaging. To correct for the phase aberration from the skull, two critical steps are needed prior to brain imaging. In the first step, the skull shape and speed of sound are acquired by either CT scans or ultrasound scans. In the ultrasound scan approach, phased array and double focusing technique are utilized, which are able to estimate the thickness of the skull with a maximum error of around 10% and the average speed of sound in the skull is underestimated by less than 2%. In the second step, the fast marching method is used to compute the phase delay based on the known skull shape and sound speed from the first step, and the computation can be completed in seconds for 2D problems. The computed phase delays are then used in combination with the conventional delay-and-sum algorithm for generating B-mode images. Images of wire phantoms with CT or ultrasound scan-based phase correction are shown to have much less artifact than the ones without correction. Errors of deducing speed of sound from CT scans are also discussed regarding its effect on the transcranial ultrasound images. Assuming the speed of sound grows linearly with the density, this study shows that, the CT-based phase correction approach can provide clear images of wire phantoms even if the speed of sound is overestimated by 400 m s(-1), or the linear coefficient is overestimated by 40%. While in this study, ultrasound scan-based phase correction performs almost equally well with the CT-based approach, potential problems are identified and discussed.

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.

Contrast-enhanced transcranial two-dimensional ultrasound imaging using shear-mode conversion at low frequency

The distortion and attenuation of transcranial ultrasound (US) signals are significant problems in US imaging of the brain. Of the variety of proposed solutions, shear-mode transmission through the skull is one of the more recent options and has been shown to reduce distortion of the US beam. This study examined the effects of transcranial shear-mode transmission on the images of a contrast-agent-filled polytetrafluoroethylene tube produced by a 32-element 750 kHz linear phased array transducer through an ex vivo human skull section. Although the tube was successfully imaged using shear-mode transmission with subharmonic imaging in 6 of 9 cases, the tube was visible in only 1 of 9 cases for both the fundamental and the second harmonic frequencies. Some improvement in the location of the axial image was seen at the fundamental frequency using shear mode. No improvement was seen at the other two frequencies, but this may be due to low transducer sensitivity. As well, neither the presence of the skull nor the incident angle changed the distance at which signals from the two tubes could be resolved. With this transducer, these distances were found to be 5 mm laterally and 3 mm axially for the fundamental and second harmonic images, and 10 mm and 5 mm for the subharmonic images. The results show that the subharmonic signal was the most successful of the three examined in penetrating a thick skull but that the success comes at the cost of image resolution.

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.

The ultrasound brain helmet: new transducers and volume registration for in vivo simultaneous multi-transducer 3-D transcranial imaging

Because stroke remains an important and time-sensitive health concern in developed nations, we present a system capable of fusing 3-D transcranial ultrasound volumes acquired from two sides of the head. This system uses custom sparse array transducers built on flexible multilayer circuits that can be positioned for simultaneous imaging through both temporal acoustic windows, allowing for potential registration of multiple real-time 3-D scans of cerebral vasculature. We examine hardware considerations for new matrix arrays-transducer design and interconnects-in this application. Specifically, it is proposed that SNR may be increased by reducing the length of probe cables. This claim is evaluated as part of the presented system through simulation, experimental data, and in vivo imaging. Ultimately, gains in SNR of 7 dB are realized by replacing a standard probe cable with a much shorter flex interconnect; higher gains may be possible using ribbon-based probe cables. In vivo images are presented, showing cerebral arteries with and without the use of microbubble contrast agent; they have been registered and fused using a simple algorithm which maximizes normalized cross-correlation.

In vivo mapping of brain elasticity in small animals using shear wave imaging

A combination of radiation force and ultrafast ultrasound imaging is used to both generate and track the propagation of a shear wave in the brain whose local speed is directly related to stiffness, characterized by the dynamic shear modulus G*. When performed on trepanated rats, this approach called shear wave imaging (SWI) provides 3-D brain elasticity maps reaching a spatial resolution of 0.7 mm×1 mm×0.4 mm with a good reproducibility (<13%). The dynamic shear modulus of brain tissues exhibits values in the 2-25 kPa range with a mean value of 12 kPa and is quantified for different anatomical regions. The anisotropy of the shear wave propagation is studied and the first in vivo anisotropy map of brain elasticity is provided. The propagation is found to be isotropic in three gray matter regions but highly anisotropic in two white matter regions. The good temporal resolution (~10 ms per acquisition) of SWI also allows a dynamic estimation of brain elasticity to within a single cardiac cycle, showing that brain pulsatility does not transiently modify local elasticity. SWI proves its potential for the study of pathological modifications of brain elasticity both in small animal models and in clinical intra-operative imaging.

Locating an acoustic point source scattered by a skull phantom via time reversal matched filtering

This paper examines the utilization of the time reversal matched filtering method to resolve the location of an acoustic point source beneath a skull phantom (variable thickness layer), without the removal of this layer. This acoustical process is examined experimentally in a water tank immersion system containing an acoustic source, a custom-made skull phantom, and a receiving transducer in a pitch-catch arrangement. The phantom is designed to approximately model the acoustic properties of an average human skull bone (minus the diploe layer), while the variable thickness of the phantom introduces a variable time delay to the acoustic wave, relative to its entry point on the phantom. This variable delay is measured and corrected for, and a matched filtering time reversed process is used to determine the location of the point source. The results of the experiment are examined for various positions of the acoustic source behind the phantom and compared to the reference cases with no phantom present. The average distance between these two cases is found to be 4.36 mm, and within the expected deviation in results due to not accounting for the effects of refraction.

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.

The ultrasound brain helmet: feasibility study of multiple simultaneous 3D scans of cerebral vasculature

We describe early stage experiments to test the feasibility of an ultrasound brain helmet to produce multiple simultaneous real-time three-dimensional (3D) scans of the cerebral vasculature from temporal and suboccipital acoustic windows of the skull. The transducer hardware and software of the Volumetrics Medical Imaging (Durham, NC, USA) real-time 3D scanner were modified to support dual 2.5 MHz matrix arrays of 256 transmit elements and 128 receive elements which produce two simultaneous 64 degrees pyramidal scans. The real-time display format consists of two coronal B-mode images merged into a 128 degrees sector, two simultaneous parasagittal images merged into a 128 degrees x 64 degrees C-mode plane and a simultaneous 64 degrees axial image. Real-time 3D color Doppler scans from a skull phantom with latex blood vessel were obtained after contrast agent injection as a proof of concept. The long-term goal is to produce real-time 3D ultrasound images of the cerebral vasculature from a portable unit capable of internet transmission thus enabling interactive 3D imaging, remote diagnosis and earlier therapeutic intervention. We are motivated by the urgency for rapid diagnosis of stroke due to the short time window of effective therapeutic intervention.

Towards aberration correction of transcranial ultrasound using acoustic droplet vaporization

We report on the first experiments demonstrating the transcranial acoustic formation of stable gas bubbles that can be used for transcranial ultrasound aberration correction. It is demonstrated that the gas bubbles can be formed transcranially by phase-transitioning single, superheated, micron-size, liquid dodecafluoropentane droplets with ultrasound, a process known as acoustic droplet vaporization (ADV). ADV was performed at 550 kHz, where the skull is less attenuating and aberrating, allowing for higher-amplitudes to be reached at the focus. Additionally, it is demonstrated that time-reversal focusing at 1 MHz can be used to correct for transcranial aberrations with a single gas bubble acting as a point beacon. Aberration correction was performed using a synthetic aperture approach and verified by the realignment of the scattered waveforms. Under the conditions described below, time-reversal aberration correction using gas bubbles resulted in a gain of 1.9 +/- 0.3 in an introduced focusing factor. This is a small fraction of the gain anticipated from complete transmit-receive of a fully-populated two-dimensional array with sub-wavelength elements. (E-mail: khaworth@umich.edu).

Compensating for bone interfaces and respiratory motion in high-intensity focused ultrasound

Bursts of focused ultrasound energy a thousand times more intense than diagnostic ultrasound have become a non-invasive option for treating cancer, from breast to prostate or uterine fibroid, during the last decade. Despite this progress, many issues still need to be addressed. First, the distortions caused by defocusing obstacles, such as the skull or ribs, on the ultrasonic therapeutic beam are still being investigated. Multi-element transducer technology must be used in order to achieve such transcranial or transcostal adaptive focusing. Second, the problem of motion artifacts, a key component in the treatment of abdominal lesions, has been shown significantly to influence the efficacy and treatment time. Though many methods have been proposed for the detection of organ motion, little work has been done to develop a comprehensive solution including motion tracking and feedback correction in real time. This paper is a review of the work achieved by authors in transcranial high-intensity focused ultrasound (HIFU), transcostal HIFU and motion compensated HIFU. For these three issues, the optimal solution can be reached using the same technology of multi-element transducers devices able to work both in transmit and receive modes.