Theoretical and experimental study on the detection limit of the micro-ring resonator based ultrasound point detectors

Combining the diffusive laser excitation and the photoacoustic signals detection, photoacoustic computed tomography (PACT) is uniquely suited for deep tissue imaging. A diffraction-limited ultrasound point detector is highly desirable for maximizing the spatial resolution and the field-of-view of the reconstructed volumetric images. Among all the available ultrasound detectors, micro-ring resonator (MRR) based ultrasound detectors offer the lowest area-normalized limit of detection (nLOD) in a miniature form-factor, making it an ideal candidate as an ultrasound point detector. However, despite their wide adoption for photoacoustic imaging, the underlying signal transduction process has not been systematically studied yet. Here we report a comprehensive theoretical model capturing the transduction of incident acoustic signals into digital data, and the associated noise propagation process, using experimentally calibrated key process parameters. The theoretical model quantifies the signal-to-noise ratio (SNR) and the nLOD under the influence of the key process variables, including the quality factor (Q-factor) of the MRR and the driving wavelength. While asserting the need for higher Q-factors, the theoretical model further quantifies the optimal driving wavelength for optimizing the nLOD. Given the MRR with a Q-factor of 1 × 105, the theoretical model predicts an optimal SNR of 30.1 dB and a corresponding nLOD of 3.75 × 10-2 mPa mm2/Hz1/2, which are in good agreement with the experimental measurements of 31.0 dB and 3.39 × 10-2 mPa mm2/Hz1/2, respectively. The reported theoretical model can be used in guiding the optimization of MRR-based ultrasonic detectors and PA experimental conditions, in attaining higher imaging resolution and contrast. The optimized operating condition has been further validated by performing PACT imaging of a human hair phantom.

Adaptive polarization photoacoustic computed tomography for biological anisotropic tissue imaging

Most photoacoustic computed tomography (PACT) systems usually ignore the anisotropy of the tissue absorption coefficient, which will lead to the lack of information in reconstructed images. In this work, the effect is addressed of the possible optical absorption anisotropy of tissue on PACT images. The functional relationship is derived between the photoacoustic response and the polarization angle of the excitation light. An adaptive polarized light photoacoustic imaging (AP-PACT) approach is proposed and shown to make up for the lack of imaging information and achieve optimal image contrast when imaging samples with anisotropic optical absorption, by utilizing the standard deviation of photoacoustic response as the feedback signal in an adaptive data acquisition process. The method is implemented both on phantom and in vitro experiments, which show that AP-PACT can recover anisotropic absorption-related information from reconstructed images and thus significantly improve their quality.

The influences of finite aperture size in photoacoustic computed tomography

In photoacoustic computed tomography (PACT), the "finite aperture effect" is often characterized as a tangential resolution that increases proportionally with the distance from the rotation center. However, this conclusion is based on the inaccurate point-detector assumption used in image reconstruction. In this study, we appropriately modeled the finite size of the acoustic detector in the back-projection (BP) based image reconstruction to improve the accuracy of the time delay calculation and systematically investigated its effects. Our results showed that the main effect of the finite aperture size is the creation of a limited high-quality imaging region (HQIR) around the scanning center, due to the directional sensitivity of the detector. We also demonstrated that the "finite aperture effect" can reduce the optimal number of detectors required for spatial anti-aliasing. These new findings provide novel perspectives for optimizing PACT systems and corresponding reconstruction methods.

Adaptive dual-speed ultrasound and photoacoustic computed tomography

Full-ring dual-modal ultrasound and photoacoustic computed tomography has unique advantages of nearly isotropic spatial resolution, complementary contrast, deep penetration, and full-view detection. However, the imaging quality may be deteriorated by the inaccurate sound speed estimation. Automatic determining and compensation for sound speed has been a long-standing problem in image reconstruction. Here, we present new adaptive dual-speed ultrasound and photoacoustic computed tomography (ADS-USPACT) to address this challenge. The system features full-view coverage (360°), high-speed dual-modal imaging (10-Hz), automated dual sound speed correction, and synergistic high imaging quality. To correct the sound speed, we develop a two-compartment method that can automatically segment the sample boundary and search for the optimal sound speed based on the rich ultrasonic pulse-echo signals. The method does not require the operator's intervention. We validate this technique in numerical simulation, phantom study, and in vivo experiments. The ADS-USPACT represents significant progress in dual-modal imaging.

Spectral crosstalk in photoacoustic computed tomography

Multispectral photoacoustic (PA) imaging faces two major challenges: the spectral coloring effect, which has been studied extensively as an optical inversion problem, and the spectral crosstalk, which is basically a result of non-ideal acoustic inversion. So far, there is no systematic work to analyze the spectral crosstalk because acoustic inversion and spectroscopic measurement are always treated as decoupled. In this work, we theorize and demonstrate through a series of simulations and experiments how imperfect acoustic inversion induces inaccurate PA spectrum measurement. We provide detailed analysis to elucidate how different factors, including limited bandwidth, limited view, light attenuation, out-of-plane signal, and image reconstruction schemes, conspire to render the measured PA spectrum inaccurate. We found that the model-based reconstruction outperforms universal back-projection in suppressing the spectral crosstalk in some cases.

Photoacoustic imaging of squirrel monkey cortical and subcortical brain regions during peripheral electrical stimulation

The investigation of neuronal activity in non-human primate models is of critical importance due to their genetic similarity to human brains. In this study, we tested the feasibility of using photoacoustic imaging for the detection of cortical and subcortical responses due to peripheral electrical stimulation in a squirrel monkey model. Photoacoustic computed tomography and photoacoustic microscopy were applied on squirrel monkeys for real-time deep subcortical imaging and optical-resolution cortical imaging, respectively. The electrically evoked hemodynamic changes in primary somatosensory cortex, premotor cortices, primary motor cortex, and underlying subcortical areas were measured. Hemodynamic responses were observed in both cortical and subcortical brain areas at the cortices during external stimulation, demonstrating the feasibility of photoacoustic technique for functional imaging of non-human primate brain.

High-Frequency 3D Photoacoustic Computed Tomography Using an Optical Microring Resonator

3D photoacoustic computed tomography (3D-PACT) has made great advances in volumetric imaging of biological tissues, with high spatial-temporal resolutions and large penetration depth. The development of 3D-PACT requires high-performance acoustic sensors with a small size, large detection bandwidth, and high sensitivity. In this work, we present a new high-frequency 3D-PACT system that uses a micro-ring resonator (MRR) as the acoustic sensor. The MRR sensor has a size of 80 μm in diameter, and was fabricated using the nanoimprint lithography technology. Using the MRR sensor, we have developed a transmission-mode 3D-PACT system that has achieved a detection bandwidth of ~23 MHz, an imaging depth of ~8 mm, a lateral resolution of 114 μm, and an axial resolution of 57 μm. We have demonstrated the 3D PACT's performance on in vitro phantoms, ex vivo mouse brain, and in vivo mouse ear and tadpole. The MRR-based 3D-PACT system can be a promising tool for structural, functional, and molecular imaging of biological tissues at depths.

Progress of clinical translation of handheld and semi-handheld photoacoustic imaging

Photoacoustic imaging (PAI), featuring rich contrast, high spatial/temporal resolution and deep penetration, is one of the fastest-growing biomedical imaging technology over the last decade. To date, numbers of handheld and semi-handheld photoacoustic imaging devices have been reported with corresponding potential clinical applications. Here, we summarize emerged handheld and semi-handheld systems in terms of photoacoustic computed tomography (PACT), optoacoustic mesoscopy (OAMes), and photoacoustic microscopy (PAM). We will discuss each modality in three aspects: laser delivery, scanning protocol, and acoustic detection. Besides new technical developments, we also review the associated clinical studies, and the advantages/disadvantages of these new techniques. In the end, we propose the challenges and perspectives of miniaturized PAI in the future.

A three-dimensional modeling method for quantitative photoacoustic breast imaging with handheld probe

By providing complementary functional information, photoacoustic (PA) breast imaging based on the handheld ultrasound (US) probe has demonstrated promising potential for breast cancer diagnosis. However, the quantitative PA imaging primarily relies on the knowledge of the optical fluence distribution in the three-dimensional (3D) heterogeneous breast tissue. Previous studies based on the handheld system generally provided two-dimensional (2D) B-scan results, which contains limited anatomical information of the tissue and the lesion. This study proposed a method to perform 3D modeling of the photon transportation for dual-modality PA/US system based on the local 3D breast anatomical information by scanning US probe. Then the calculated optical fluence distribution can be used for PA imaging. Our phantom and clinical pilot study results demonstrated that this method has potential to improve the accuracy of the quantitative PA breast imaging, and it can also be used in other clinical implementations.

Adaptive photoacoustic computed tomography

For many optical imaging modalities, image qualities are inevitably degraded by wavefront distortions caused by varying light speed. In optical microscopy and astronomy, adaptive optics (AO) has long been applied to compensate for such unwanted aberrations. Photoacoustic computed tomography (PACT), despite relying on the ultrasonic wave for image formation, suffers from the acoustic version of the same problem. However, this problem has traditionally been regarded as an inverse problem of jointly reconstructing both the initial pressure and the sound speed distributions. In this work, we proposed a method similar to indirect wavefront sensing in AO. We argued that wavefront distortions can be extracted and corrected by a frequency domain analysis of local images. In addition to an adaptively reconstructed aberration-free image, the speed of sound map can be subsequently estimated. We demonstrated the method by in silico, phantom, and in vivo experiments.

Transcranial photoacoustic computed tomography based on a layered back-projection method

A major challenge of transcranial human brain photoacoustic computed tomography (PACT) is correcting for the acoustic aberration induced by the skull. Here, we present a modified universal back-projection (UBP) method, termed layered UBP (L-UBP), that can de-aberrate the transcranial PA signals by accommodating the skull heterogeneity into conventional UBP. In L-UBP, the acoustic medium is divided into multiple layers: the acoustic coupling fluid layer between the skull and detectors, the skull layer, and the brain tissue layer, which are assigned different acoustic properties. The transmission coefficients and wave conversion are considered at the fluid-skull and skull-tissue interfaces. Simulations of transcranial PACT using L-UBP were conducted to validate the method. Ex vivo experiments with a newly developed three-dimensional PACT system with 1-MHz center frequency demonstrated that L-UBP can substantially improve the image quality compared to conventional UBP.

Minimally invasive photoacoustic imaging: Current status and future perspectives

Photoacoustic imaging (PAI) is an emerging biomedical imaging modality that is based on optical absorption contrast, capable of revealing distinct spectroscopic signatures of tissue at high spatial resolution and large imaging depths. However, clinical applications of conventional non-invasive PAI systems have been restricted to examinations of tissues at depths less than a few cm due to strong light attenuation. Minimally invasive photoacoustic imaging (miPAI) has greatly extended the landscape of PAI by delivering excitation light within tissue through miniature fibre-optic probes. In the past decade, various miPAI systems have been developed with demonstrated applicability in several clinical fields. In this article, we present an overview of the current status of miPAI and our thoughts on future perspectives.

Review of cost reduction methods in photoacoustic computed tomography

Photoacoustic Computed Tomography (PACT) is a major configuration of photoacoustic imaging, a hybrid noninvasive modality for both functional and molecular imaging. PACT has rapidly gained importance in the field of biomedical imaging due to superior performance as compared to conventional optical imaging counterparts. However, the overall cost of developing a PACT system is one of the challenges towards clinical translation of this novel technique. The cost of a typical commercial PACT system originates from optical source, ultrasound detector, and data acquisition unit. With growing applications of photoacoustic imaging, there is a tremendous demand towards reducing its cost. In this review article, we have discussed various approaches to reduce the overall cost of a PACT system, and provided a cost estimation to build a low-cost PACT system.

Optimizing photoacoustic image reconstruction using cross-platform parallel computation

Three-dimensional (3D) image reconstruction involves the computations of an extensive amount of data that leads to tremendous processing time. Therefore, optimization is crucially needed to improve the performance and efficiency. With the widespread use of graphics processing units (GPU), parallel computing is transforming this arduous reconstruction process for numerous imaging modalities, and photoacoustic computed tomography (PACT) is not an exception. Existing works have investigated GPU-based optimization on photoacoustic microscopy (PAM) and PACT reconstruction using compute unified device architecture (CUDA) on either C++ or MATLAB only. However, our study is the first that uses cross-platform GPU computation. It maintains the simplicity of MATLAB, while improves the speed through CUDA/C++ − based MATLAB converted functions called MEXCUDA. Compared to a purely MATLAB with GPU approach, our cross-platform method improves the speed five times. Because MATLAB is widely used in PAM and PACT, this study will open up new avenues for photoacoustic image reconstruction and relevant real-time imaging applications.

High-resolution deep functional imaging of the whole mouse brain by photoacoustic computed tomography in vivo

Photoacoustic computed tomography (PACT) is a non-invasive imaging technique offering high contrast, high resolution, and deep penetration in biological tissues. We report a PACT system equipped with a high frequency linear transducer array for mapping the microvascular network of a whole mouse brain with the skull intact and studying its hemodynamic activities. The linear array was scanned in the coronal plane to collect data from different angles, and full-view images were synthesized from the limited-view images in which vessels were only partially revealed. We investigated spontaneous neural activities in the deep brain by monitoring the concentration of hemoglobin in the blood vessels and observed strong interhemispherical correlations between several chosen functional regions, both in the cortical layer and in the deep regions. We also studied neural activities during an epileptic seizure and observed the epileptic wave spreading around the injection site and the wave propagating in the opposite hemisphere.

A forward-adjoint operator pair based on the elastic wave equation for use in transcranial photoacoustic computed tomography

Photoacoustic computed tomography (PACT) is an emerging imaging modality that exploits optical contrast and ultrasonic detection principles to form images of the photoacoustically induced initial pressure distribution within tissue. The PACT reconstruction problem corresponds to an inverse source problem in which the initial pressure distribution is recovered from measurements of the radiated wavefield. A major challenge in transcranial PACT brain imaging is compensation for aberrations in the measured data due to the presence of the skull. Ultrasonic waves undergo absorption, scattering and longitudinal-to-shear wave mode conversion as they propagate through the skull. To properly account for these effects, a wave-equation-based inversion method should be employed that can model the heterogeneous elastic properties of the skull. In this work, a forward model based on a finite-difference time-domain discretization of the three-dimensional elastic wave equation is established and a procedure for computing the corresponding adjoint of the forward operator is presented. Massively parallel implementations of these operators employing multiple graphics processing units (GPUs) are also developed. The developed numerical framework is validated and investigated in computer19 simulation and experimental phantom studies whose designs are motivated by transcranial PACT applications.

A numerical analysis of a semi-dry coupling configuration in photoacoustic computed tomography for infant brain imaging

In the application of photoacoustic human infant brain imaging, debubbled ultrasound gel or water is commonly used as a couplant for ultrasonic transducers due to their acoustic properties. The main challenge in using such a couplant is its discomfort for the patient. In this study, we explore the feasibility of a semi-dry coupling configuration to be used in photoacoustic computed tomography (PACT) systems. The coupling system includes an inflatable container consisting of a thin layer of Aqualene with ultrasound gel or water inside of it. Finite element method (FEM) is used for static and dynamic structural analysis of the proposed configuration to be used in PACT for infant brain imaging. The outcome of the analysis is an optimum thickness of Aqualene in order to meet the weight tolerance requirement with the least attenuation and best impedance match to recommend for an experimental setting.

Joint Reconstruction of Absorbed Optical Energy Density and Sound Speed Distributions in Photoacoustic Computed Tomography: A Numerical Investigation

Photoacoustic computed tomography (PACT) is a rapidly emerging bioimaging modality that seeks to reconstruct an estimate of the absorbed optical energy density within an object. Conventional PACT image reconstruction methods assume a constant speed-of-sound (SOS), which can result in image artifacts when acoustic aberrations are significant. It has been demonstrated that incorporating knowledge of an object's SOS distribution into a PACT image reconstruction method can improve image quality. However, in many cases, the SOS distribution cannot be accurately and/or conveniently estimated prior to the PACT experiment. Because variations in the SOS distribution induce aberrations in the measured photoacoustic wavefields, certain information regarding an object's SOS distribution is encoded in the PACT measurement data. Based on this observation, a joint reconstruction (JR) problem has been proposed in which the SOS distribution is concurrently estimated along with the sought-after absorbed optical energy density from the photoacoustic measurement data. A broad understanding of the extent to which the JR problem can be accurately and reliably solved has not been reported. In this work, a series of numerical experiments is described that elucidate some important properties of the JR problem that pertain to its practical feasibility. To accomplish this, an optimization-based formulation of the JR problem is developed that yields a non-linear iterative algorithm that alternatively updates the two image estimates. Heuristic analytic insights into the reconstruction problem are also provided. These results confirm the ill-conditioned nature of the joint reconstruction problem that will present significant challenges for practical applications.

Investigation of the far-field approximation for modeling a transducer's spatial impulse response in photoacoustic computed tomography

When ultrasonic transducers with large detecting areas and/or compact measurement geometries are employed in photoacoustic computed tomography (PACT), the spatial resolution of reconstructed images can be significantly degraded. Our goal in this work is to clarify the domain of validity of an imaging model that mitigates such effects by use of a far-field approximation. Computer-simulation studies are described that demonstrate the far-field-based imaging model is highly accurate for a practical 3D PACT imaging geometry employed in an existing small animal imaging system. For use in special cases where the far-field approximation is violated, an extension of the far-field-based imaging model is proposed that divides the transducer face into a small number of rectangular patches that are each described accurately by use of the far-field approximation.

Biomedical photoacoustics in China

During the last decade, along with its explosive growth globally, biomedical photoacoustics has become a rapidly growing research field in China as well. In particular, photoacoustic tomography (PAT), capable of imaging intact biological tissue in vivo at great depths, has generated intense interest among Chinese researchers. This review briefly summarizes the current status and recent progress of the research in PAT in China. The focus is on the technology development and biomedical applications of three representative embodiments of PAT: photoacoustic microscopy, photoacoustic computed tomography, and photoacoustic endoscopy. In addition, recent development and studies in other related areas are also reviewed shortly.