Observation of Guided Acoustic Waves in a Human Skull

Ray theory-based compounded plane wave ultrasound imaging for aberration corrected transcranial imaging: Phantom experiments and simulations

Compounded plane wave imaging (CPWI) allows high-frame-rate measurement and has been one of the most promising modalities for real-time brain imaging. However, ultrasonic brain imaging using the CPWI modality is usually performed with a worn thin or removal of the skull layer. Otherwise, the skull layer is expected to distort the ultrasonic wavefronts and significantly decrease intracranial imaging quality. The motivation of this study is to investigate a CPWI method for transcranial brain imaging with the skull layer. A coordinate transformation ray-tracing (CTRT) approach was proposed to track the distorted ultrasonic wavefronts and calculate the time delays for the ultrasound plane wave passing through the skull layer. With an accurate correction for the time delays in beamforming, the CTRT-based CPWI could achieve high-quality intracranial images with the presence of skulls. The proposed CTRT-based CPWI method was verified using a simplified three-layer transcranial model. The full-wave simulation demonstrated that CTRT could accurately (i.e., relative percentage error less than0.18%) track the distorted transmitting wavefront through skull. Compared with the CPWI without aberration correction, the CTRT-based CPWI provided high-quality intracranial imaging and could accurately localize intracranial point scatterers; specifically, positioning error decreases from 0.5 mm to 0.1 mm on average in the axial direction and from 0.7 mm to 0.1 mm on average in the lateral direction. As the compounded angles increased in the CTRT-based CPWI, the contrast improved by 16.2 dB on average for the region of interest, and the array performance indicator (representing resolution) decreased by 4.0 on average for the intracranial point scatterers. The CTRT is of low computational cost compared with full wave simulation. This study suggested that the proposed CTRT-based CPWI might have the potential for real-time and non-invasive transcranial aberration-corrected imaging.

Optimizing transcranial ultrasound delivery at large incident angles by leveraging cranial leaky guided wave dispersion

We investigate the role of leaky guided waves in transcranial ultrasound transmission in temporal and parietal bones at large incidence angles. Our numerical and experimental results show that the dispersion characteristics of the fundamental leaky guided wave mode with longitudinal polarization can be leveraged to estimate the critical angle above which efficient shear mode conversion takes place, and below which major transmission drops can be expected. Simulations that employ a numerical propagator matrix and a Semi-Analytical approach establish the transcranial dispersion characteristics and transmission coefficients at different incident angles. Experimental transmission tests conducted at 500 kHz and radiation tests performed in the 200-800 kHz range confirm the numerical findings in terms of transmitted peak pressure and frequency-radiation angle spectra, based on which the connection between critical angles, dispersion and transmission is demonstrated. Our results support the identification of transcranial ultrasound strategies that leverage shear mode conversion, which is less sensitive to phase aberrations compared to normal incidence ultrasound. These findings can also enable higher transmission rates in cranial bones with low porosity by leveraging dispersion information extracted through signal processing, without requiring measurement of geometric and mechanical properties of the cranial bone.

Coregistered transcranial optoacoustic and magnetic resonance angiography of the human brain

Imaging modalities capable of visualizing the human brain have led to major advances in neurology and brain research. Multi-spectral optoacoustic tomography (MSOT) has gained importance for studying cerebral function in rodent models due to its unique capability to map changes in multiple hemodynamic parameters and to directly visualize neural activity within the brain. The technique further provides molecular imaging capabilities that can facilitate early disease diagnosis and treatment monitoring. However, transcranial imaging of the human brain is hampered by acoustic attenuation and other distortions introduced by the skull. Here, we demonstrate non-invasive transcranial MSOT angiography of pial veins through the temporal bone of an adult healthy volunteer. Time-of-flight (TOF) magnetic resonance angiography (MRA) and T₁-weighted structural magnetic resonance imaging (MRI) were further acquired to facilitate anatomical registration and interpretation. The superior middle cerebral vein in the temporal cortex was identified in the MSOT images, matching its location observed in the TOF-MRA images. These initial results pave the way toward the application of MSOT in clinical brain imaging.

Multimodal Noninvasive Functional Neurophotonic Imaging of Murine Brain-Wide Sensory Responses

Modern optical neuroimaging approaches are expanding the ability to elucidate complex brain function. Diverse imaging contrasts enable direct observation of neural activity with functional sensors along with the induced hemodynamic responses. To date, decoupling the complex interplay of neurovascular coupling and dynamical physiological states has remained challenging when employing single-modality functional neuroimaging readings. A hybrid fluorescence optoacoustic tomography platform combined with a custom data processing pipeline based on statistical parametric mapping is devised, attaining the first noninvasive observation of simultaneous calcium and hemodynamic activation patterns using optical contrasts. Correlated changes in the oxy- and deoxygenated hemoglobin, total hemoglobin, oxygen saturation, and rapid GCaMP6f fluorescence signals are observed in response to peripheral sensory stimulation. While the concurrent epifluorescence serves to corroborate and complement the functional optoacoustic observations, the latter further aids in decoupling the rapid calcium responses from the slowly varying background in the fluorescence recordings mediated by hemodynamic changes. The hybrid imaging platform expands the capabilities of conventional neuroimaging methods to provide more comprehensive functional readings for studying neurovascular and neurometabolic coupling mechanisms and related diseases.

Deep Learning-Based Photoacoustic Imaging of Vascular Network Through Thick Porous Media

Photoacoustic imaging is a promising approach used to realize in vivo transcranial cerebral vascular imaging. However, the strong attenuation and distortion of the photoacoustic wave caused by the thick porous skull greatly affect the imaging quality. In this study, we developed a convolutional neural network based on U-Net to extract the effective photoacoustic information hidden in the speckle patterns obtained from vascular network images datasets under porous media. Our simulation and experimental results show that the proposed neural network can learn the mapping relationship between the speckle pattern and the target, and extract the photoacoustic signals of the vessels submerged in noise to reconstruct high-quality images of the vessels with a sharp outline and a clean background. Compared with the traditional photoacoustic reconstruction methods, the proposed deep learning-based reconstruction algorithm has a better performance with a lower mean absolute error, higher structural similarity, and higher peak signal-to-noise ratio of reconstructed images. In conclusion, the proposed neural network can effectively extract valid information from highly blurred speckle patterns for the rapid reconstruction of target images, which offers promising applications in transcranial photoacoustic imaging.

Increasing the transmission efficiency of transcranial ultrasound using a dual-mode conversion technique based on Lamb waves

Transcranial focused ultrasound (FUS) is a noninvasive treatment for brain tumors and neuromodulation. Based on normal incidence, conventional FUS techniques use a focused or an array of ultrasonic transducers to overcome the attenuation and absorption of ultrasound in the skull; however, this remains the main limitation of using FUS. A dual-mode conversion technique based on Lamb waves is proposed to achieve high transmission efficiency. This concept was validated using the finite element analysis (FEA) and experiments based on changes in the incident angle. Aluminum, plexiglass, and a human skull were used as materials with different attenuations. The transmission loss was calculated for each material, and the results were compared with the reflectance function of the Lamb waves. Oblique incidence based on dual-mode conversion exhibited a better transmission efficiency than that of a normal incidence for all of the specimens. The total transmission losses for the materials were 13.7, 15.46, and 3.91 dB less than those associated with the normal incidence. A wedge transducer was designed and fabricated to implement the proposed method. The results demonstrated the potential applicability of the dual-mode conversion technique for the human skull.

Effect of Skull Porous Trabecular Structure on Transcranial Ultrasound Imaging in the Presence of Elastic Wave Mode Conversion at Varying Incidence Angle

With the advancement of aberration correction techniques, transcranial ultrasound imaging has exhibited great potential in applications such as imaging neurological function and guiding therapeutic ultrasound. However, the feasibility of transcranial imaging varies among individuals because of the differences in skull acoustic properties. To better understand the fundamental mechanisms underlying the variation in imaging performance, the effect of the structure of the porous trabecular bone on transcranial imaging performance (i.e., target localization errors and resolution) was investigated for the first time through the use of elastic wave simulations and experiments. Simulation studies using high-resolution computed tomography data from ex vivo skull samples revealed that imaging at large incidence angles reduced the target localization error for skulls having low porosity; however, as skull porosity increased, large angles of incidence resulted in degradation of resolution and increased target localization errors. Experimental results indicate that imaging at normal incidence introduced a localization error of 1.85 ± 0.10 mm, while imaging at a large incidence angle (40°) resulted in an increased localization error of 6.54 ± 1.33 mm and caused a single point target to no longer appear as a single, coherent target in the resulting image, which is consistent with simulation results. This first investigation of the effects of skull microstructure on transcranial ultrasound imaging indicates that imaging performance is highly dependent on the porosity of the skull, particularly at non-normal angles of incidence.

Vibration-based elastic parameter identification of the diploë and cortical tables in dry cranial bones

Various human skull models feature a layered cranial structure composed of homogeneous cortical tables and the inner diploë. However, there is a lack of fundamental validation work of such three-layer cranial bone models by combining high-fidelity computational modeling and rigorous experiments. Here, non-contact vibration experiments are conducted on an assortment of dry bone segments from the largest cranial bone regions (parietal, frontal, occipital, and temporal) to estimate the first handful of modal frequencies and damping ratios, as well as mode shapes, in the audio frequency regime. Numerical models that consider the cortical tables and the diploë as domains with separate isotropic material properties are constructed for each bone segment using a routine that identifies the cortical table-diploë boundaries from micro-computed tomography scan images, and reconstructs a three-dimensional geometry layer by layer. The material properties for cortical tables and diploë are obtained using a Hounsfield Unit-based mass density calculation combined with a parameter identification scheme for Young's modulus estimation. With the identified parameters, the average error between experimental and numerical modal frequencies is 1.3% and the modal assurance criterion values for most modes are above 0.90, indicating that the layered model is suitable for predicting the vibrational behavior of cranial bone. The proposed layered modeling and identified elastic parameters are also useful to support computational modeling of cranial guided waves and mode conversion in medical ultrasound. Additionally, the diploë elastic properties are rarely reported in the literature, making this work a fundamental characterization effort that can guide in the selection of material properties for human head models that consider layered cranial bone.

Radiation Characteristics of Cranial Leaky Lamb Waves

We numerically and experimentally investigate the dispersion properties of leaky Lamb waves in the cranial bone. Cranial Lamb waves leak energy from the skull into the brain when propagating at speeds higher than the speed of sound in the surrounding fluid. The understanding of their radiation mechanism is significantly complicated by the geometric and mechanical characteristics of the cortical tables and the trabecular bone (diploë). Toward such understanding, we here analyze the sub-1.0 MHz radiation angle dispersion spectrum of porous bone phantoms and parietal bone geometries obtained from μ CT scans. Our numerical results show that, when diploic pores are physically modeled, leakage angles computed from time transient finite-element analyses correspond to those predicted by an equivalent three-layered fluid-loaded waveguide model. For the bone geometries analyzed, two main leaky branches are observed in the near-field dispersion spectrum: a fast wave radiated at small angles, which is related to the fastest fundamental Lamb mode supported by the cranial bone, and a slower wave radiated at larger angles. This observation is also confirmed by experimental tests carried out on an immersed parietal bone.

Experimental identification of high order Lamb waves and estimation of the mechanical properties of a dry human skull

We experimentally investigate and characterize high order Lamb wave modes in a dry human skull. Specifically, we show that the diploë supports distinct wave modes in the sub-1.0 MHz frequency regime, and we employ these modes for the estimation of equivalent mechanical properties of cortical and trabecular bones. These modes are efficiently generated in a parietal region by direct contact excitation with a wedge beam transducer, and are recorded via infrared laser vibrometry. Frequency/wavenumber data are estimated using a matrix pencil method applied to wavefield measurements recorded on the outer cortical surface. The semi-analytical finite element model of an equivalent three-layered plate provides the platform for the identification of wave modes based on their through-the-thickness profiles, and supports the estimation of equivalent mechanical properties in conjunction with an optimization algorithm developed for this purpose. The results presented herein illustrate how high order Lamb waves can be used to gain understanding of the wave properties of a human skull and to estimate the orthotropic and equivalent isotropic mechanical properties of cortical and trabecular bones.

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).

Spherical Array System for High-Precision Transcranial Ultrasound Stimulation and Optoacoustic Imaging in Rodents

Ultrasound can be delivered transcranially to ablate brain tissue, open the blood-brain barrier, or affect neural activity. Transcranial focused ultrasound in small rodents is typically done with low-frequency single-element transducers, which results in unspecific targeting and impedes the concurrent use of fast neuroimaging methods. In this article, we devised a wide-angle spherical array bidirectional interface for high-resolution parallelized optoacoustic imaging and transcranial ultrasound (POTUS) delivery in the same target regions. The system operates between 3 and 9 MHz, allowing to generate and steer focal spots with widths down to [Formula: see text] across a field of view covering the entire mouse brain, while the same array is used to capture high-resolution 3-D optoacoustic data in real time. We showcase the system's versatile beam-forming capacities as well as volumetric optoacoustic imaging capabilities and discuss its potential to noninvasively monitor brain activity and various effects of ultrasound emission.

Anisotropic Longitudinal Wave Propagation in Swine Skull

To understand the in-plane elastic character of ultrasonic waves in the skull, longitudinal wave velocities were studied in the MHz range using a conventional pulse technique. Taking advantage of the thickness of swine skulls, anisotropic in-plane wave velocity changes in the outer and diploe layers were experimentally investigated using structural information measured by X-ray computer tomography (CT). The velocities in the thin inner layer were difficult to measure. The main trabecular alignment (MTA) in the thick swine diploe layer was almost perpendicular to the thickness direction and changed with position inside the skull. The degree of anisotropy of in-plane longitudinal wave velocity ranged 1.07-1.33 in both outer and diploe layers, depending on position and swine sample. The angle of the fastest velocity in the outer layer was different from that in most parts of the diploe layer. Anisotropic character in the diploe layer gradually changed with position in the thickness direction.

Multifocal structured illumination optoacoustic microscopy

Optoacoustic (OA) imaging has the capacity to effectively bridge the gap between macroscopic and microscopic realms in biological imaging. High-resolution OA microscopy has so far been performed via point-by-point scanning with a focused laser beam, thus greatly restricting the achievable imaging speed and/or field of view. Herein we introduce multifocal structured illumination OA microscopy (MSIOAM) that attains real-time 3D imaging speeds. For this purpose, the excitation laser beam is shaped to a grid of focused spots at the tissue surface by means of a beamsplitting diffraction grating and a condenser and is then scanned with an acousto-optic deflector operating at kHz rates. In both phantom and in vivo mouse experiments, a 10 mm wide volumetric field of view was imaged with 15 Hz frame rate at 28 μm spatial resolution. The proposed method is expected to greatly aid in biological investigations of dynamic functional, kinetic, and metabolic processes across multiple scales.

Listening to tissues with new light: recent technological advances in photoacoustic imaging

Photoacoustic tomography (PAT), or optoacoustic tomography, has achieved remarkable progress in the past decade, benefiting from the joint developments in optics, acoustics, chemistry, computing and mathematics. Unlike pure optical or ultrasound imaging, PAT can provide unique optical absorption contrast as well as widely scalable spatial resolution, penetration depth and imaging speed. Moreover, PAT has inherent sensitivity to tissue's functional, molecular, and metabolic state. With these merits, PAT has been applied in a wide range of life science disciplines, and has enabled biomedical research unattainable by other imaging methods. This Review article aims at introducing state-of-the-art PAT technologies and their representative applications. The focus is on recent technological breakthroughs in structural, functional, molecular PAT, including super-resolution imaging, real-time small-animal whole-body imaging, and high-sensitivity functional/molecular imaging. We also discuss the remaining challenges in PAT and envisioned opportunities.