Ultrashort echo-time MRI versus CT for skull aberration correction in MR-guided transcranial focused ultrasound: In vitro comparison on human calvaria

Model-based navigation of transcranial focused ultrasound neuromodulation in humans: Application to targeting of the amygdala and thalamus

BACKGROUND

Transcranial focused ultrasound (tFUS) neuromodulation has shown promise in animals but is challenging to translate to humans because of the thicker skull that heavily scatters ultrasound waves.

OBJECTIVE

We develop and disseminate a model-based navigation (MBN) tool for acoustic dose delivery in the presence of skull aberrations that is easy to use by non-specialists.

METHODS

We pre-compute acoustic beams for thousands of virtual transducer locations on the scalp of the subject under study. We use the hybrid angular spectrum solver mSOUND, which runs in ∼4 seconds per solve per CPU yielding pre-computation times under one hour for scalp meshes with up to 4,000 faces and a parallelization factor of 5. We combine this pre-computed set of beam solutions with optical tracking, thus allowing real-time display of the tFUS beam as the operator freely navigates the transducer around the subject' scalp. We assess the impact of MBN versus line-of-sight targeting (LOST) positioning in simulations of 13 subjects.

RESULTS

Our navigation tool has a display refresh rate of ∼10 Hz. In our simulations, MBN increased the acoustic dose in the thalamus and amygdala by 8-67% compared to LOST and avoided complete target misses that affected 10-20% of LOST cases. MBN yields a lower variability of the deposited dose across subjects than LOST.

CONCLUSIONS

MBN may yield greater and more consistent (less variable) ultrasound dose deposition than transducer placement with line-of-sight targeting, and thus may become a helpful tool to improve the efficacy of tFUS neuromodulation.

Robust deep learning estimation of cortical bone porosity from MR T1-weighted images for individualized transcranial focused ultrasound planning

Objective

Transcranial focused ultrasound (tFUS) is an emerging neuromodulation approach that has been demonstrated in animals but is difficult to translate to humans because of acoustic attenuation and scattering in the skull. Optimal dose delivery requires subject-specific skull porosity estimates which has traditionally been done using CT. We propose a deep learning (DL) estimation of skull porosity from T1-weighted MRI images which removes the need for radiation-inducing CT scans.

Approach

We evaluate the impact of different DL approaches, including network architecture, input size and dimensionality, multichannel inputs, data augmentation, and loss functions. We also propose back-propagation in the mask (BIM), a method whereby only voxels inside the skull mask contribute to training. We evaluate the robustness of the best model to input image noise and MRI acquisition parameters and propagate porosity estimation errors in thousands of beam propagation scenarios.

Main results

Our best performing model is a cGAN with a ResNet-9 generator with 3D 64×64×64 inputs trained with L1 and L2 losses. The model achieved a mean absolute error of 6.9% in the test set, compared to 9.5% with the pseudo-CT of Izquierdo et al. (38% improvement) and 9.4% with the generic pixel-to-pixel image translation cGAN pix2pix (36% improvement). Acoustic dose distributions in the thalamus were more accurate with our approach than with the pseudo-CT approach of both Burgos et al. and Izquierdo et al, resulting in near-optimal treatment planning and dose estimation at all frequencies compared to CT (reference).

Significance

Our DL approach porosity estimates with ~7% error, is robust to input image noise and MRI acquisition parameters (sequence, coils, field strength) and yields near-optimal treatment planning and dose estimates for both central (thalamus) and lateral brain targets (amygdala) in the 200-1000 kHz frequency range.

Aberration correction in abdominal histotripsy

In transabdominal histotripsy, ultrasound pulses are focused on the body to noninvasively destroy soft tissues via cavitation. However, the ability to focus is limited by phase aberration, or decorrelation of the ultrasound pulses due to spatial variation in the speed of sound throughout heterogeneous tissue. Phase aberration shifts, broadens, and weakens the focus, thereby reducing the safety and efficacy of histotripsy therapy. This paper reviews and discusses aberration effects in histotripsy and in related therapeutic ultrasound techniques (e.g., high intensity focused ultrasound), with an emphasis on aberration by soft tissues. Methods for aberration correction are reviewed and can be classified into two groups: model-based methods, which use segmented images of the tissue as input to an acoustic propagation model to predict and compensate phase differences, and signal-based methods, which use a receive-capable therapy array to detect phase differences by sensing acoustic signals backpropagating from the focus. The relative advantages and disadvantages of both groups of methods are discussed. Importantly, model-based methods can correct focal shift, while signal-based methods can restore substantial focal pressure, suggesting that both methods should be combined in a 2-step approach. Aberration correction will be critical to improving histotripsy treatments and expanding the histotripsy treatment envelope to enable non-invasive, non-thermal histotripsy therapy for more patients.

Generating Patient-Specific Acoustic Simulations for Transcranial Focused Ultrasound Procedures Based on Optical Tracking Information

Optical tracking is a real-time transducer positioning method for transcranial focused ultrasound (tFUS) procedures, but the predicted focus from optical tracking typically does not incorporate subject-specific skull information. Acoustic simulations can estimate the pressure field when propagating through the cranium but rely on accurately replicating the positioning of the transducer and skull in a simulated space. Here, we develop and characterize the accuracy of a workflow that creates simulation grids based on optical tracking information in a neuronavigated phantom with and without transmission through an ex vivo skull cap. The software pipeline could replicate the geometry of the tFUS procedure within the limits of the optical tracking system (transcranial target registration error (TRE): 3.9 ± 0.7 mm). The simulated focus and the free-field focus predicted by optical tracking had low Euclidean distance errors of 0.5±0.1 and 1.2±0.4 mm for phantom and skull cap, respectively, and some skull-specific effects were captured by the simulation. However, the TRE of simulation informed by optical tracking was 4.6±0.2, which is as large or greater than the focal spot size used by many tFUS systems. By updating the position of the transducer using the original TRE offset, we reduced the simulated TRE to 1.1 ± 0.4 mm. Our study describes a software pipeline for treatment planning, evaluates its accuracy, and demonstrates an approach using MR-acoustic radiation force imaging as a method to improve dosimetry. Overall, our software pipeline helps estimate acoustic exposure, and our study highlights the need for image feedback to increase the accuracy of tFUS dosimetry.

Evaluation of synthetically generated computed tomography for use in transcranial focused ultrasound procedures

Purpose

Transcranial focused ultrasound (tFUS) is a therapeutic ultrasound method that focuses sound through the skull to a small region noninvasively and often under magnetic resonance imaging (MRI) guidance. CT imaging is used to estimate the acoustic properties that vary between individual skulls to enable effective focusing during tFUS procedures, exposing patients to potentially harmful radiation. A method to estimate acoustic parameters in the skull without the need for CT is desirable.

Approach

We synthesized CT images from routinely acquired T1-weighted MRI using a 3D patch-based conditional generative adversarial network and evaluated the performance of synthesized CT (sCT) images for treatment planning with tFUS. We compared the performance of sCT with real CT (rCT) images for tFUS planning using Kranion and simulations using the acoustic toolbox, k-Wave. Simulations were performed for 3 tFUS scenarios: (1) no aberration correction, (2) correction with phases calculated from Kranion, and (3) phase shifts calculated from time reversal.

Results

From Kranion, the skull density ratio, skull thickness, and number of active elements between rCT and sCT had Pearson's correlation coefficients of 0.94, 0.92, and 0.98, respectively. Among 20 targets, differences in simulated peak pressure between rCT and sCT were largest without phase correction (12.4 % ± 8.1 % ) and smallest with Kranion phases (7.3 % ± 6.0 % ). The distance between peak focal locations between rCT and sCT was < 1.3    mm for all simulation cases.

Conclusions

Real and synthetically generated skulls had comparable image similarity, skull measurements, and acoustic simulation metrics. Our work demonstrated similar results for 10 testing cases comparing MR-sCTs and rCTs for tFUS planning. Source code and a docker image with the trained model are available at https://github.com/han-liu/SynCT_TcMRgFUS.

A simulation study on the sensitivity of transcranial ray-tracing ultrasound modeling to skull properties

In transcranial focused ultrasound therapies, such as treating essential tremor via thermal ablation in the thalamus, acoustic energy is focused through the skull using a phased-array transducer. Ray tracing is a computationally efficient method that can correct skull-induced phase aberrations via per-element phase delay calculations using patient-specific computed tomography (CT) data. However, recent studies show that variations in CT-derived Hounsfield unit may account for only 50% of the speed of sound variability in human skull specimens, potentially limiting clinical transcranial ultrasound applications. Therefore, understanding the sensitivity of treatment planning methods to material parameter variations is essential. The present work uses a ray-tracing simulation model to explore how imprecision in model inputs, arising from clinically significant uncertainties in skull properties or considerations of acoustic phenomena, affects acoustic focusing quality through the skull. We propose and validate new methods to optimize ray-tracing skull simulations for clinical treatment planning, relevant for predicting intracranial target's thermal rise, using experimental data from ex-vivo human skulls.

Acoustoelectric Time-Reversal for Ultrasound Phase-Aberration Correction

Acoustoelectric imaging (AEI) is a technique that combines ultrasound (US) with radio frequency recording to detect and map local current source densities. This study demonstrates a new method called acoustoelectric time reversal (AETR), which uses AEI of a small current source to correct for phase aberrations through a skull or other US-aberrating layers with applications to brain imaging and therapy. Simulations conducted at three different US frequencies (0.5, 1.5, and 2.5 MHz) were performed through media layered with different sound speeds and geometries to induce aberrations of the US beam. Time delays of the acoustoelectric (AE) signal from a monopole within the medium were calculated for each element to enable corrections using AETR. Uncorrected aberrated beam profiles were compared with those after applying AETR corrections, which demonstrated a strong recovery (29%-100%) of lateral resolution and increases in focal pressure up to 283%. To further demonstrate the practical feasibility of AETR, we further conducted bench-top experiments using a 2.5 MHz linear US array to perform AETR through 3-D-printed aberrating objects. These experiments restored lost lateral restoration up to 100% for the different aberrators and increased focal pressure up to 230% after applying AETR corrections. Cumulatively, these results highlight AETR as a powerful tool for correcting focal aberrations in the presence of a local current source with applications to AEI, US imaging, neuromodulation, and therapy.

Current and Emerging Systems for Focused Ultrasound-Mediated Blood-Brain Barrier Opening

With an ever-growing list of neurological applications of focused ultrasound (FUS), there has been a consequent increase in the variety of systems for delivering ultrasound energy to the brain. Specifically, recent successful pilot clinical trials of blood-brain barrier (BBB) opening with FUS have generated substantial interest in the future applications of this relatively novel therapy, with divergent, purpose-built technologies emerging. With many of these technologies at various stages of pre-clinical and clinical investigation, this article seeks to provide an overview and analysis of the numerous medical devices in active use and under development for FUS-mediated BBB opening.

Transcranial ultrasound simulations: A review

Transcranial ultrasound is more and more used for therapy and imaging of the brain. However, the skull is a highly attenuating and aberrating medium, with different structures and acoustic properties among samples and even within a sample. Thus, case-specific simulations are needed to perform transcranial focused ultrasound interventions safely. In this article, we provide a review of the different methods used to model the skull and to simulate ultrasound propagation through it.

Classical and Learned MR to Pseudo-CT Mappings for Accurate Transcranial Ultrasound Simulation

Model-based treatment planning for transcranial ultrasound therapy typically involves mapping the acoustic properties of the skull from an X-ray computed tomography (CT) image of the head. Here, three methods for generating pseudo-CT (pCT) images from magnetic resonance (MR) images were compared as an alternative to CT. A convolutional neural network (U-Net) was trained on paired MR-CT images to generate pCT T images from either T1-weighted or zero-echo time (ZTE) MR images (denoted tCT and zCT, respectively). A direct mapping from ZTE to pCT was also implemented (denoted cCT). When comparing the pCT and ground-truth CT images for the test set, the mean absolute error was 133, 83, and 145 Hounsfield units (HU) across the whole head, and 398, 222, and 336 HU within the skull for the tCT, zCT, and cCT images, respectively. Ultrasound simulations were also performed using the generated pCT images and compared to simulations based on CT. An annular array transducer was used targeting the visual or motor cortex. The mean differences in the simulated focal pressure, focal position, and focal volume were 9.9%, 1.5 mm, and 15.1% for simulations based on the tCT images; 5.7%, 0.6 mm, and 5.7% for the zCT; and 6.7%, 0.9 mm, and 12.1% for the cCT. The improved results for images mapped from ZTE highlight the advantage of using imaging sequences, which improves the contrast of the skull bone. Overall, these results demonstrate that acoustic simulations based on MR images can give comparable accuracy to those based on CT.

Comparison between MR and CT imaging used to correct for skull-induced phase aberrations during transcranial focused ultrasound

Transcranial focused ultrasound with the InSightec Exablate system uses thermal ablation for the treatment of movement and mood disorders and blood brain barrier disruption for tumor therapy. The system uses computed tomography (CT) images to calculate phase corrections that account for aberrations caused by the human skull. This work investigates whether magnetic resonance (MR) images can be used as an alternative to CT images to calculate phase corrections. Phase corrections were calculated using the gold standard hydrophone method and the standard of care InSightec ray tracing method. MR binary image mask, MR-simulated-CT (MRsimCT), and CT images of three ex vivo human skulls were supplied as inputs to the InSightec ray tracing method. The degassed ex vivo human skulls were sonicated with a 670 kHz hemispherical phased array transducer (InSightec Exablate 4000). 3D raster scans of the beam profiles were acquired using a hydrophone mounted on a 3-axis positioner system. Focal spots were evaluated using six metrics: pressure at the target, peak pressure, intensity at the target, peak intensity, positioning error, and focal spot volume. Targets at the geometric focus and 5 mm lateral to the geometric focus were investigated. There was no statistical difference between any of the metrics at either target using either MRsimCT or CT for phase aberration correction. As opposed to the MRsimCT, the use of CT images for aberration correction requires registration to the treatment day MR images; CT misregistration within a range of ± 2 degrees of rotation error along three dimensions was shown to reduce focal spot intensity by up to 9.4%. MRsimCT images used for phase aberration correction for the skull produce similar results as CT-based correction, while avoiding both CT to MR registration errors and unnecessary patient exposure to ionizing radiation.

Numerical investigation of the energy distribution of Low-intensity transcranial focused ultrasound neuromodulation for hippocampus

OBJECTIVE

Ultrasonic neuromodulation as a safe and non-invasive brain stimulation method that delivers a low-intensity, focused ultrasound to nervous system tissue in a targeted area of the brain. The objective of this study is to numerically investigate the ultrasound wave propagation and the energy distribution within the brain tissues using customized single element focused ultrasound transducers (SEFT), targeting the hippocampus.

METHODS

A high resolution detailed human head model with seven tissue types was constructed from magnetic resonance imaging (MRI). A full-wave finite-difference time-domain simulation platform, Sim4life, was then used to simulate a 3D non-linear ultrasound wave equation to the specific region of interest, the hippocampus. Three customized SEFT were used to test the effect of transducer positions, and another customized transducer was used to compare the sensitivity effect on heterogeneous and homogeneous brain models. Finally, the sensitivity and performance of low intensity focusing ultrasound stimulation were evaluated.

RESULTS

An optimized application of SEFT was customized to deliver 100 W/m² intensity of energy deposition at the hippocampus region. About 85.65% of the generated volume beam was delivered to the targeted hippocampus region and the beam overlap parameter was affected by different transducer positions. Deflection angle changes of SEFT at the range of ± 5% did not have a significant effect on energy delivery and position displacement. Only 0.5% of peak pressure change was observed between heterogeneous and homogeneous brain models. The sensitivity analysis also showed that the sound speed is the most influential acoustic parameter.

SIGNIFICANCE

This study demonstrated that ultrasound neuromodulation targeting the depth brain tissue of the hippocampus could be a potential and promising alternative method to some non-acoustic brain stimulation modalities. In the numerical study of ultrasound brain stimulations, ultrasound parameters and the brain model need to be properly determined to simulate the ultrasonic neuromodulations.

Zero-TE MRI: principles and applications in the head and neck

Zero echo-time (ZTE) MRI is a novel imaging technique that utilizes ultrafast readouts to capture signal from short-T2 tissues. Additional sequence advantages include rapid imaging times, silent scanning, and artifact resistance. A robust application of this technology is imaging of cortical bone without the use of ionizing radiation, thus representing a viable alternative to CT for both rapid screening and "one-stop-shop" MRI. Although ZTE is increasingly used in musculoskeletal and body imaging, neuroimaging applications have historically been limited by complex anatomy and pathology. In this article, we review the imaging physics of ZTE including pulse sequence options, practical limitations, and image reconstruction. We then discuss optimization of settings for ZTE bone neuroimaging including acquisition, processing, segmentation, synthetic CT generation, and artifacts. Finally, we examine clinical utility of ZTE in the head and neck with imaging examples including malformations, trauma, tumors, and interventional procedures.

Ultrasound stimulation of the motor cortex during tonic muscle contraction

Transcranial ultrasound stimulation (tUS) shows potential as a noninvasive brain stimulation (NIBS) technique, offering increased spatial precision compared to other NIBS techniques. However, its reported effects on primary motor cortex (M1) are limited. We aimed to better understand tUS effects in human M1 by performing tUS of the hand area of M1 (M1hand) during tonic muscle contraction of the index finger. Stimulation during muscle contraction was chosen because of the transcranial magnetic stimulation-induced phenomenon known as cortical silent period (cSP), in which transcranial magnetic stimulation (TMS) of M1hand involuntarily suppresses voluntary motor activity. Since cSP is widely considered an inhibitory phenomenon, it presents an ideal parallel for tUS, which has often been proposed to preferentially influence inhibitory interneurons. Recording electromyography (EMG) of the first dorsal interosseous (FDI) muscle, we investigated effects on muscle activity both during and after tUS. We found no change in FDI EMG activity concurrent with tUS stimulation. Using single-pulse TMS, we found no difference in M1 excitability before versus after sparsely repetitive tUS exposure. Using acoustic simulations in models made from structural MRI of the participants that matched the experimental setups, we estimated in-brain pressures and generated an estimate of cumulative tUS exposure experienced by M1hand for each subject. We were unable to find any correlation between cumulative M1hand exposure and M1 excitability change. We also present data that suggest a TMS-induced MEP always preceded a near-threshold cSP.

Acoustic properties across the human skull

Transcranial ultrasound is emerging as a noninvasive tool for targeted treatments of brain disorders. Transcranial ultrasound has been used for remotely mediated surgeries, transient opening of the blood-brain barrier, local drug delivery, and neuromodulation. However, all applications have been limited by the severe attenuation and phase distortion of ultrasound by the skull. Here, we characterized the dependence of the aberrations on specific anatomical segments of the skull. In particular, we measured ultrasound propagation properties throughout the perimeter of intact human skulls at 500 kHz. We found that the parietal bone provides substantially higher transmission (average pressure transmission 31 ± 7%) and smaller phase distortion (242 ± 44 degrees) than frontal (13 ± 2%, 425 ± 47 degrees) and occipital bone regions (16 ± 4%, 416 ± 35 degrees). In addition, we found that across skull regions, transmission strongly anti-correlated (R=-0.79) and phase distortion correlated (R=0.85) with skull thickness. This information guides the design, positioning, and skull correction functionality of next-generation devices for effective, safe, and reproducible transcranial focused ultrasound therapies.

Transcranial MR-Guided Histotripsy System

Histotripsy has been previously shown to treat a wide range of locations through excised human skulls in vitro. In this article, a transcranial magnetic resonance (MR)-guided histotripsy (tcMRgHt) system was developed, characterized, and tested in the in vivo pig brain through an excised human skull. A 700-kHz, 128-element MR-compatible phased-array ultrasound transducer with a focal depth of 15 cm was designed and fabricated in-house. Support structures were also constructed to facilitate transcranial treatment. The tcMRgHt array was acoustically characterized with a peak negative pressure up to 137 MPa in free field, 72 MPa through an excised human skull with aberration correction, and 48.4 MPa without aberration correction. The electronic focal steering range through the skull was 33.5 mm laterally and 50 mm axially, where a peak negative pressure above the 26-MPa cavitation intrinsic threshold can be achieved. The MR compatibility of the tcMRgHt system was assessed quantitatively using SNR, B0 field map, and B1 field map in a clinical 3T magnetic resonance imaging (MRI) scanner. Transcranial treatment using electronic focal steering was validated in red blood cell phantoms and in vivo pig brain through an excised human skull. In two pigs, targeted cerebral tissue was successfully treated through the human skull as confirmed by MRI. Excessive bleeding or edema was not observed in the peri-target zones by the time of pig euthanasia. These results demonstrated the feasibility of using this preclinical tcMRgHt system for in vivo transcranial treatment in a swine model.

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.

Implementation of a Skull-Conformal Phased Array for Transcranial Focused Ultrasound Therapy

OBJECTIVE

To implement a skull-conformal phased array for ultrasound-guided transcranial focused ultrasound therapy with improved patient comfort.

METHODS

Using patient-specific computed tomography and MRI neuroimaging data, tightly-conforming helmet scaffolds were designed computationally. The helmet scaffolds were designed to hold reusable transducer modules at near-normal incidence in an optimal configuration for the treatment location(s) of interest. Numerical simulations of trans-skull ultrasound propagation were performed to evaluate different conformal array designs and to compare with hemispherical arrays similar to those employed clinically. A 4096-element phased array was constructed by 3D printing a helmet scaffold optimised for an ex vivo human skullcap, and its performance was evaluated via benchtop and in vivo experiments.

RESULTS

Acoustic field measurements confirmed the system's ability to focus through human skull bone using simulation-based transcranial aberration corrections. Preliminary in vivo testing demonstrated safe trans-human skull blood-brain barrier (BBB) opening in rodents.

CONCLUSION

Patient-specific conformal ultrasound phased arrays appear to be a feasible and safe approach for conducting transcranial BBB opening procedures.

SIGNIFICANCE

Skull-conformal phased arrays stand to improve patient comfort and have the potential to accelerate the adoption of transcranial FUS therapy by improving access to the technology.

Acoustic Attenuation: Multifrequency Measurement and Relationship to CT and MR Imaging

Transcranial magnetic resonance-guided focused ultrasound (tcMRgFUS) is gaining significant acceptance as a noninvasive treatment for motion disorders and shows promise for novel applications such as blood-brain barrier opening for tumor treatment. A typical procedure relies on CT-derived acoustic property maps to simulate the transfer of ultrasound through the skull. Accurate estimates of the acoustic attenuation in the skull are essential to accurate simulations, but there is no consensus about how attenuation should be estimated from CT images and there is interest in exploring MR as a predictor of attenuation in the skull. In this study, we measure the acoustic attenuation at 0.5, 1, and 2.25 MHz in 89 samples taken from two ex vivo human skulls. CT scans acquired with a variety of X-ray energies, reconstruction kernels, and reconstruction algorithms, and MR images acquired with ultrashort and zero echo time sequences are used to estimate the average Hounsfield unit value, MR magnitude, and T2^* value in each sample. The measurements are used to develop a model of attenuation as a function of frequency and each individual imaging parameter.

Experimental and Computational Investigation of Guided Waves in a Human Skull

We investigate guided (Lamb) waves in a human cadaver skull through experiments and computational simulations. Ultrasonic wedge transducers and scanning laser Doppler vibrometry are used respectively to excite and measure Lamb waves propagating in the cranial bone of a degassed skull. Measurements are performed over a section of the parietal bone and temporal bone spanning the squamous suture. The experimental data are analyzed for the identification of wave modes and the characterization of dispersion properties. In the parietal bone, for instance, the A₀ wave mode is excited between 200 and 600 kHz, and higher-order Lamb waves are excited from 1 to 1.8 MHz. From the experimental dispersion curves and average thickness extracted from the skull computed tomography scan, we estimate average isotropic material properties that capture the essential dispersion characteristics using a semi-analytical finite-element model. We also explore the leaky and non-leaky wave behavior of the degassed skull with water loading in the cranial cavity. Successful excitation of leaky Lamb waves is confirmed (for higher-order wave modes with phase velocity faster than the speed of sound in water) from 500 kHz to 1.5 MHz, which may find applications in imaging and therapeutics at the brain periphery or skull-brain interface (e.g., for metastases). The non-leaky A₀ Lamb wave mode propagates between 200 and 600 kHz, with or without fluid loading, for potential use in skull-related diagnostics and imaging (e.g., for sutures).

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.

Achilles Subtendon Structure and Behavior as Evidenced From Tendon Imaging and Computational Modeling

The Achilles tendon is the largest and strongest tendon in the human body and is essential for storing elastic energy and positioning the foot for walking and running. Recent research into Achilles tendon anatomy and mechanics has revealed the importance of the Achilles subtendons, which are unique and semi-independent structures arising from each of the three muscular heads of the triceps surae. Of particular importance is the ability for the subtendons to slide, the role that this has in healthy tendons, and the alteration of this property in aging and disease. In this work, we discuss technical approaches that have led to the current understanding of Achilles subtendons, particularly imaging and computational modeling. We introduce a 3D geometrical model of the Achilles subtendons, built from dual-echo UTE MRI. We revisit and discuss computational models of Achilles subtendon twisting suggesting that optimal twist reduces both rupture loads and stress concentrations by distributing stresses. Second harmonic generation imaging shows collagenous subtendons within a rabbit Achilles tendon; a clear absence of signal between the subtendons indicates an inter-subtendon region on the order of 30 μm in our rabbit animal model. Entry of wheat germ agglutinin in both the inter-fascicular and the inter-subtendon regions suggests a glycoprotein-containing inter-subtendon matrix which may facilitate low friction sliding of the subtendons in healthy mammals. Lastly, we present a new computational model coupled with human exercise trials to demonstrate the magnitude of Achilles subtendon sliding which occurs during rehabilitation exercises for Achilles tendinopathy, and shows that specific exercise can maximize subtendon sliding and interface strains, without maximizing subtendon strains. This work demonstrates the value of imaging and computational modeling for probing tendon structure-function relationships and may serve to inform and develop treatments for Achilles tendinopathy.

Computationally Efficient Transcranial Ultrasonic Focusing: Taking Advantage of the High Correlation Length of the Human Skull

The phase correction necessary for transcranial ultrasound therapy requires numerical simulation to noninvasively assess the phase shift induced by the skull bone. Ideally, the numerical simulations need to be fast enough for clinical implementation in a brain therapy protocol and to provide accurate estimation of the phase shift to optimize the refocusing through the skull. In this article, we experimentally performed transcranial ultrasound focusing at 900 kHz on N = 5 human skulls. To reduce the computation time, we propose here to perform the numerical simulation at 450 kHz and use the corresponding phase shifts experimentally at 900 kHz. We demonstrate that a 450-kHz simulation restores 94.2% of the pressure when compared with a simulation performed at 900 kHz and 85.0% of the gold standard pressure obtained by an invasive time reversal procedure based on the signal recorded by a hydrophone placed at the target. From a 900- to 450-kHz simulation, the grid size is divided by 8, and the computation time is divided by 10.

Neuronavigated Repetitive Transcranial Ultrasound Stimulation Induces Long-Lasting and Reversible Effects on Oculomotor Performance in Non-human Primates

Since the late 2010s, Transcranial Ultrasound Stimulation (TUS) has been used experimentally to carryout safe, non-invasive stimulation of the brain with better spatial resolution than Transcranial Magnetic Stimulation (TMS). This innovative stimulation method has emerged as a novel and valuable device for studying brain function in humans and animals. In particular, single pulses of TUS directed to oculomotor regions have been shown to modulate visuomotor behavior of non-human primates during 100 ms ultrasound pulses. In the present study, a sustained effect was induced by applying 20-s trains of neuronavigated repetitive Transcranial Ultrasound Stimulation (rTUS) to oculomotor regions of the frontal cortex in three non-human primates performing an antisaccade task. With the help of MRI imaging and a frame-less stereotactic neuronavigation system (SNS), we were able to demonstrate that neuronavigated TUS (outside of the MRI scanner) is an efficient tool to carry out neuromodulation procedures in non-human primates. We found that, following neuronavigated rTUS, saccades were significantly modified, resulting in shorter latencies compared to no-rTUS trials. This behavioral modulation was maintained for up to 20 min. Oculomotor behavior returned to baseline after 18-31 min and could not be significantly distinguished from the no-rTUS condition. This study is the first to show that neuronavigated rTUS can have a persistent effect on monkey behavior with a quantified return-time to baseline. The specificity of the effects could not be explained by auditory confounds.

Enhanced Numerical Method for the Design of 3-D-Printed Holographic Acoustic Lenses for Aberration Correction of Single-Element Transcranial Focused Ultrasound

The correction of transcranial focused ultrasound aberrations is a relevant issue for enhancing various non-invasive medical treatments. The emission through multi-element phased arrays has been the most widely accepted method to improve focusing in recent years; however, the number and size of transducers represent a bottleneck that limits the focusing accuracy of the technique. To overcome this limitation, a new disruptive technology, based on 3-D-printed acoustic lenses, has recently been proposed. As the submillimeter precision of the latest generation of 3-D printers has been proven to overcome the spatial limitations of phased arrays, a new challenge is to improve the accuracy of the numerical simulations required to design this type of ultrasound lens. In the study described here, we evaluated two improvements in the numerical model applied in previous works for the design of 3-D-printed lenses: (i) allowing the propagation of shear waves in the skull by means of its simulation as an isotropic solid and (ii) introduction of absorption into the set of equations that describes the dynamics of the wave in both fluid and solid media. The results obtained in the numerical simulations are evidence that the inclusion of both s-waves and absorption significantly improves focusing.

Rapid dual-RF, dual-echo, 3D ultrashort echo time craniofacial imaging: A feasibility study

PURPOSE

To develop a dual-radiofrequency (RF), dual-echo, 3D ultrashort echo-time (UTE) pulse sequence and bone-selective image reconstruction for rapid high-resolution craniofacial MRI.

METHODS

The proposed pulse sequence builds on recently introduced dual-RF UTE imaging. While yielding enhanced bone specificity by exploiting high sensitivity of short T₂ signals to variable RF pulse widths, the parent technique exacts a 2-fold scan time penalty relative to standard dual-echo UTE. In the proposed method, the parent sequence's dual-RF scheme was incorporated into dual-echo acquisitions while radial view angles are varied every pulse-to-pulse repetition period. The resulting 4 echoes (2 for each RF) were combined by view-sharing to construct 2 sets of k-space data sets, corresponding to short and long TEs, respectively, leading to a 2-fold increase in imaging efficiency. Furthermore, by exploiting the sparsity of bone signals in echo-difference images, acceleration was achieved by solving a bone-sparsity constrained image reconstruction problem. In vivo studies were performed to evaluate the effectiveness of the proposed acceleration approaches in comparison to the parent method.

RESULTS

The proposed technique achieves 1.1-mm isotropic skull imaging in 3 minutes without visual loss of image quality, compared to the parent technique (scan time = 12 minutes). Bone-specific images and corresponding 3D renderings of the skull were found to depict the expected craniofacial anatomy over the entire head.

CONCLUSION

The proposed method is able to achieve high-resolution volumetric craniofacial images in a clinically practical imaging time, and thus may prove useful as a potential alternative to computed tomography.

Design and Implementation of a Transmit/Receive Ultrasound Phased Array for Brain Applications

Focused ultrasound phased array systems have attracted increased attention for brain therapy applications. However, such systems currently lack a direct and real-time method to intraoperatively monitor ultrasound pressure distribution for securing treatment. This study proposes a dual-mode ultrasound phased array system design to support transmit/receive operations for concurrent ultrasound exposure and backscattered focal beam reconstruction through a spherically focused ultrasound array. A 256-channel ultrasound transmission system was used to transmit focused ultrasonic energy (full 256 channels), with an extended implementation of multiple-channel receiving function (up to 64 channels) using the same 256-channel ultrasound array. A coherent backscatter-received beam formation algorithm was implemented to map the point spread function (PSF) and focal beam distribution under a free-field/transcranial environment setup, with the backscattering generated from a strong scatterer (a point reflector or a microbubble-perfused tube) or a weakly scattered tissue-mimicking graphite phantom. Our results showed that PSF and focal beam can be successfully reconstructed and visualized in free-field conditions and can also be transcranially reconstructed following skull-induced aberration correction. In vivo experiments were conducted to demonstrate its capability to preoperatively and semiquantitatively map a focal beam to guide blood-brain barrier opening. The proposed system may have potential for real-time guidance of ultrasound brain intervention, and may facilitate the design of a dual-mode ultrasound phased array for brain therapeutic applications.

Simulating transvertebral ultrasound propagation with a multi-layered ray acoustics model

The simulation accuracy of transvertebral ultrasound propagation using a multi-layered ray acoustics model based on CT-derived vertebral geometry was investigated through comparison with experimental measurements of pressure fields in ex vivo human vertebral foramen. A spherically focused transducer (5 cm diameter, f-number 1.2, 514 kHz) was geometrically focused to the centre of individual thoracic vertebral foramen, through the posterior bony elements. Transducer propagation paths through the laminae and the spinous processes were tested. Simulation transducer-vertebra configurations were registered to experiment transducer-vertebra configurations, and simulation accuracy of the simulation model was evaluated for predicting maximum transmitted pressure to the canal, voxel pressure in the canal, and focal distortion. Accuracy in predicting maximum transmitted pressure was calculated by vertebra, and it is shown that simulation predicts maximum pressure with a greater degree of accuracy than a vertebra-specific insertion loss. Simulation error in voxel pressure was evaluated using root-mean-square error and cross-correlation, and found to be similar to the water-only case. Simulation accuracy in predicting focal distortion was evaluated by comparing experiment and simulation maximum pressure location and weighted  >50% focal volume location. Average simulation error across all measurements and simulations in maximum pressure location and weighted  >50% focal volume location were 2.3 mm and 1.5 mm, respectively. These errors are small relative to the dimensions of the transducer focus (4.9 mm full width half maximum), the spinal cord (10 mm diameter), and vertebral canal diameter (15-20 mm diameter). These results suggest that ray acoustics can be applied to simulating transvertebral ultrasound propagation.

Inside/outside the brain binary cavitation localization based on the lowpass filter effect of the skull on the harmonic content: a proof of concept study

Cavitation activity induced by ultrasound may occur during high intensity focused ultrasound (HIFU) treatment, due to bubble nucleation under high peak negative pressure, and during blood-brain-barrier (BBB) disruption, due to injected ultrasound contrast agents (UCAs). Such microbubble activity has to be monitored to assess the safety and efficiency of ultrasonic brain treatments. In this study, we aim at assessing whether cavitation occurs within cerebral tissue by binary discriminating cavitation activity originating from the inside or the outside of the skull. The results were obtained from both in vitro experiments mimicking BBB opening, by using UCA flow, and in vitro thermal necrosis in calf brain samples. The sonication was applied using a 1 MHz focused transducer and the acoustic response of the microbubbles was recorded with a wideband passive cavitation detector. The spectral content of the recorded signal was used to localize microbubble activity. Since the skull acts as a low pass filter, the ratio of high harmonics to low harmonics is lower for cavitation events located inside the skull compared to events outside the skull. Experiments showed that the ratio of the 5/2 ultraharmonic to the 1/2 subharmonic for binary localization cavitation activity achieves 100% sensitivity and specificity for both monkey and human skulls. The harmonic ratio of the fourth to the second harmonic provided 100% sensitivity and 96% and 46% specificity on a non-human primate for thermal necrosis and BBB opening, respectively. Nonetheless, the harmonic ratio remains promising for human applications, as the experiments showed 100% sensitivity and 100% specificity for both thermal necrosis and BBB opening through the human skull. The study requires further validation on a larger number of skull samples.

Measurements of the Relationship Between CT Hounsfield Units and Acoustic Velocity and How It Changes With Photon Energy and Reconstruction Method

Transcranial magnetic resonance-guided focused ultrasound continues to gain traction as a noninvasive treatment option for a variety of pathologies. Focusing ultrasound through the skull can be accomplished by adding a phase correction to each element of a hemispherical transducer array. The phase corrections are determined with acoustic simulations that rely on speed of sound estimates derived from CT scans. While several studies have investigated the relationship between acoustic velocity and CT Hounsfield units (HUs), these studies have largely ignored the impact of X-ray energy, reconstruction method, and reconstruction kernel on the measured HU, and therefore the estimated velocity, and none have measured the relationship directly. In this paper, 91 ex vivo human skull fragments from two skulls are imaged by 80 CT scans with a variety of energies and reconstruction methods. The average HU from each fragment is found for each scan and correlated with the speed of sound measured using a through transmission technique in that fragment. As measured by the -squared value, the results show that CT is able to account for 23%-53% of the variation in velocity in the human skull. Both the X-ray energy and the reconstruction technique significantly alter the -squared value and the linear relationship between HU and speed of sound in bone. Accounting for these variations will lead to more accurate phase corrections and more efficient transmission of acoustic energy through the skull.

Enhanced Energy Localization in Hyperthermia Treatment Based on Hybrid Electromagnetic and Ultrasonic System: Proof of Concept with Numerical Simulations

This paper proposes a hybrid hyperthermia treatment system, utilizing two noninvasive modalities for treating brain tumors. The proposed system depends on focusing electromagnetic (EM) and ultrasound (US) energies. The EM hyperthermia subsystem enhances energy localization by incorporating a multichannel wideband setting and coherent-phased-array technique. A genetic algorithm based optimization tool is developed to enhance the specific absorption rate (SAR) distribution by reducing hotspots and maximizing energy deposition at tumor regions. The treatment performance is also enhanced by augmenting an ultrasonic subsystem to allow focused energy deposition into deep tumors. The therapeutic faculty of ultrasonic energy is assessed by examining the control of mechanical alignment of transducer array elements. A time reversal (TR) approach is then investigated to address challenges in energy focus in both subsystems. Simulation results of the synergetic effect of both modalities assuming a simplified model of human head phantom demonstrate the feasibility of the proposed hybrid technique as a noninvasive tool for thermal treatment of brain tumors.

Accurate simulation of transcranial ultrasound propagation for ultrasonic neuromodulation and stimulation

Non-invasive, focal neurostimulation with ultrasound is a potentially powerful neuroscientific tool that requires effective transcranial focusing of ultrasound to develop. Time-reversal (TR) focusing using numerical simulations of transcranial ultrasound propagation can correct for the effect of the skull, but relies on accurate simulations. Here, focusing requirements for ultrasonic neurostimulation are established through a review of previously employed ultrasonic parameters, and consideration of deep brain targets. The specific limitations of finite-difference time domain (FDTD) and k-space corrected pseudospectral time domain (PSTD) schemes are tested numerically to establish the spatial points per wavelength and temporal points per period needed to achieve the desired accuracy while minimizing the computational burden. These criteria are confirmed through convergence testing of a fully simulated TR protocol using a virtual skull. The k-space PSTD scheme performed as well as, or better than, the widely used FDTD scheme across all individual error tests and in the convergence of large scale models, recommending it for use in simulated TR. Staircasing was shown to be the most serious source of error. Convergence testing indicated that higher sampling is required to achieve fine control of the pressure amplitude at the target than is needed for accurate spatial targeting.

Sensitivity of simulated transcranial ultrasound fields to acoustic medium property maps

High intensity transcranial focused ultrasound is an FDA approved treatment for essential tremor, while low-intensity applications such as neurostimulation and opening the blood brain barrier are under active research. Simulations of transcranial ultrasound propagation are used both for focusing through the skull, and predicting intracranial fields. Maps of the skull acoustic properties are necessary for accurate simulations, and can be derived from medical images using a variety of methods. The skull maps range from segmented, homogeneous models, to fully heterogeneous models derived from medical image intensity. In the present work, the impact of uncertainties in the skull properties is examined using a model of transcranial propagation from a single element focused transducer. The impact of changes in bone layer geometry and the sound speed, density, and acoustic absorption values is quantified through a numerical sensitivity analysis. Sound speed is shown to be the most influential acoustic property, and must be defined with less than 4% error to obtain acceptable accuracy in simulated focus pressure, position, and volume. Changes in the skull thickness of as little as 0.1 mm can cause an error in peak intracranial pressure of greater than 5%, while smoothing with a 1 [Formula: see text] kernel to imitate the effect of obtaining skull maps from low resolution images causes an increase of over 50% in peak pressure. The numerical results are confirmed experimentally through comparison with sonications made through 3D printed and resin cast skull bone phantoms.

Rapid full-wave phase aberration correction method for transcranial high-intensity focused ultrasound therapies

Background

Non-invasive high-intensity focused ultrasound (HIFU) can be used to treat a variety of disorders, including those in the brain. However, the differences in acoustic properties between the skull and the surrounding soft tissue cause aberrations in the path of the ultrasonic beam, hindering or preventing treatment.

Methods

We present a method for correcting these aberrations that is fast, full-wave, and allows for corrections at multiple treatment locations. The method is simulation-based: an acoustic model is built based on high-resolution CT scans, and simulations are performed using the hybrid angular spectrum (HAS) method to determine the phases needed for correction.

Results

Computation of corrections for clinically applicable resolutions can be achieved in approximately 15 min. Experimental results with a plastic model designed to mimic the aberrations caused by the skull show that the method can recover 95 % of the peak pressure obtained using hydrophone-based time-reversal methods. Testing using an ex vivo human skull flap resulted in recovering up to 70 % of the peak pressure at the focus and 61 % when steering (representing, respectively, a 1.52- and 1.19-fold increase in the peak pressure over the uncorrected case). Additionally, combining the phase correction method with rapid HAS simulations allows evaluation of such treatment metrics as the effect of misregistration on resulting pressure levels.

Conclusions

The method presented here is able to rapidly compute phases required to improve ultrasound focusing through the skull at multiple treatment locations. Combining phase correction with rapid simulation techniques allows for evaluation of various treatment metrics such as the effect of steering on pressure levels. Since the method computes 3D pressure patterns, it may also be suitable for predicting off-focus hot spots during treatments—a primary concern for transcranial HIFU. Additionally, the plastic-skull method presented here may be a useful tool in evaluating the effectiveness of phase correction methods.

A 3D model of the Achilles tendon to determine the mechanisms underlying nonuniform tendon displacements

The Achilles is the thickest tendon in the body and is the primary elastic energy-storing component during running. The form and function of the human Achilles is complex: twisted structure, intratendinous interactions, and differential motor control from the triceps surae muscles make Achilles behavior difficult to intuit. Recent in vivo imaging of the Achilles has revealed nonuniform displacement patterns that are not fully understood and may result from complex architecture and musculotendon interactions. In order to understand which features of the Achilles tendon give rise to the nonuniform deformations observed in vivo, we used computational modeling to predict the mechanical contributions from different features of the tendon. The aims of this study are to: (i) build a novel computational model of the Achilles tendon based on ultrashort echo time MRI, (ii) compare simulated displacements with published in vivo ultrasound measures of displacement, and (iii) use the model to elucidate the effects of tendon twisting, intratendon sliding, retrocalcaneal insertion, and differential muscle forces on tendon deformation. Intratendon sliding and differential muscle forces were found to be the largest factors contributing to displacement nonuniformity between tendon regions. Elimination of intratendon sliding or muscle forces reduced displacement nonuniformity by 96% and 85%, respectively, while elimination of tendon twist and the retrocalcaneal insertion reduced displacement nonuniformity by only 35% and 3%. These results suggest that changes in the complex internal structure of the tendon alter the interaction between muscle forces and tendon behavior and therefore may have important implications on muscle function during movement.

MR thermometry for focused ultrasound monitoring utilizing model predictive filtering and ultrasound beam modeling

Background

A major challenge in using magnetic resonance temperature imaging (MRTI) to monitor focused ultrasound (FUS) applications is achieving high spatio-temporal resolution over a large field of view (FOV). This is important to accurately monitor all ultrasound (US) power depositions. Magnetic resonance (MR) subsampling in conjunction with thermal model-based reconstruction of the MRTI utilizing Pennes bioheat transfer equation (PBTE) is one promising approach. The thermal properties used in the thermal model are often estimated from a pre-treatment, low-power sonication.

Methods

In this proof-of-concept study we investigate the use of US simulations computed using the hybrid angular spectrum (HAS) method to estimate the US power deposition density Q, thereby avoiding the pre-treatment sonication and any potential tissue damage. MRTI reconstructions are performed using a thermal model-based reconstruction method called model predictive filtering (MPF). Experiments are performed in a homogeneous gelatin phantom and in a gelatin phantom with embedded plastic skull. MPF reconstructions are compared to separate sonications imaged with fully sampled data over a smaller FOV. Temperature root-mean-square errors (RMSE) and focal spot positions and shapes are evaluated.

Results

HAS simulations accurately predict the location of the focal spot (to within 1 mm) in both phantoms. Accurate temperature maps (RMSE below 1 °C), where the location of the focal spot agrees well with fully sampled “truth” (to within 1 mm), are also achieved in both phantoms.

Conclusions

HAS simulations can be used to accurately predict the focal spot location in homogeneous media and when focusing through an aberrating plastic skull. The HAS simulated power deposition (Q) patterns can be used in the MPF thermal model-based reconstruction to obtain accurate temperature maps with high spatio-temporal resolution over large FOVs.

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.

Ultrasound treatment of neurological diseases--current and emerging applications

Like cardiovascular disease and cancer, neurological disorders present an increasing challenge for an ageing population. Whereas nonpharmacological procedures are routine for eliminating cancer tissue or opening a blocked artery, the focus in neurological disease remains on pharmacological interventions. Setbacks in clinical trials and the obstacle of access to the brain for drug delivery and surgery have highlighted the potential for therapeutic use of ultrasound in neurological diseases, and the technology has proved useful for inducing focused lesions, clearing protein aggregates, facilitating drug uptake, and modulating neuronal function. In this Review, we discuss milestones in the development of therapeutic ultrasound, from the first steps in the 1950s to recent improvements in technology. We provide an overview of the principles of diagnostic and therapeutic ultrasound, for surgery and transient opening of the blood-brain barrier, and its application in clinical trials of stroke, Parkinson disease and chronic pain. We discuss the promising outcomes of safety and feasibility studies in preclinical models, including rodents, pigs and macaques, and efficacy studies in models of Alzheimer disease. We also consider the challenges faced on the road to clinical translation.

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.