Fast and High-Resolution Three-Dimensional Hybrid-Domain Photoacoustic Imaging Incorporating Analytical-Focused Transducer Beam Amplitude

Influence mechanism of the temporal duration of laser irradiation on photoacoustic technique: a review

Significance

In the photoacoustic (PA) technique, the laser irradiation in the time domain (i.e., laser pulse duration) governs the characteristics of PA imaging-it plays a crucial role in the optical-acoustic interaction, the generation of PA signals, and the PA imaging performance.

Aim

We aim to provide a comprehensive analysis of the impact of laser pulse duration on various aspects of PA imaging, encompassing the signal-to-noise ratio, the spatial resolution of PA imaging, the acoustic frequency spectrum of the acoustic wave, the initiation of specific physical phenomena, and the photothermal-PA (PT-PA) interaction/conversion.

Approach

By surveying and reviewing the state-of-the-art investigations, we discuss the effects of laser pulse duration on the generation of PA signals in the context of biomedical PA imaging with respect to the aforementioned aspects.

Results

First, we discuss the impact of laser pulse duration on the PA signal amplitude and its correlation with the lateral resolution of PA imaging. Subsequently, the relationship between the axial resolution of PA imaging and the laser pulse duration is analyzed with consideration of the acoustic frequency spectrum. Furthermore, we examine the manipulation of the pulse duration to trigger physical phenomena and its relevant applications. In addition, we elaborate on the tuning of the pulse duration to manipulate the conversion process and ratio from the PT to PA effect.

Conclusions

We contribute to the understanding of the physical mechanisms governing pulse-width-dependent PA techniques. By gaining insight into the mechanism behind the influence of the laser pulse, we can trigger the pulse-with-dependent physical phenomena for specific PA applications, enhance PA imaging performance in biomedical imaging scenarios, and modulate PT-PA conversion by tuning the pulse duration precisely.

A Review of Emerging Electromagnetic-Acoustic Sensing Techniques for Healthcare Monitoring

Conventional electromagnetic (EM) sensing techniques such as radar and LiDAR are widely used for remote sensing, vehicle applications, weather monitoring, and clinical monitoring. Acoustic techniques such as sonar and ultrasound sensors are also used for consumer applications, such as ranging and in vivo medical/healthcare applications. It has been of long-term interest to doctors and clinical practitioners to realize continuous healthcare monitoring in hospitals and/or homes. Physiological and biopotential signals in real-time serve as important health indicators to predict and prevent serious illness. Emerging electromagnetic-acoustic (EMA) sensing techniques synergistically combine the merits of EM sensing with acoustic imaging to achieve comprehensive detection of physiological and biopotential signals. Further, EMA enables complementary fusion sensing for challenging healthcare settings, such as real-world long-term monitoring of treatment effects at home or in remote environments. This article reviews various examples of EMA sensing instruments, including implementation, performance, and application from the perspectives of circuits to systems. The novel and significant applications to healthcare are discussed. Three types of EMA sensors are presented: (1) Chip-based radar sensors for health status monitoring, (2) Thermo-acoustic sensing instruments for biomedical applications, and (3) Photoacoustic (PA) sensing and imaging systems, including dedicated reconstruction algorithms were reviewed from time-domain, frequency-domain, time-reversal, and model-based solutions. The future of EMA techniques for continuous healthcare with enhanced accuracy supported by artificial intelligence (AI) is also presented.

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.

High-fidelity deconvolution for acoustic-resolution photoacoustic microscopy enabled by convolutional neural networks

Acoustic-resolution photoacoustic microscopy (AR-PAM) image resolution is determined by the point spread function (PSF) of the imaging system. Previous algorithms, including Richardson-Lucy (R-L) deconvolution and model-based (MB) deconvolution, improve spatial resolution by taking advantage of the PSF as prior knowledge. However, these methods encounter the problems of inaccurate deconvolution, meaning the deconvolved feature size and the original one are not consistent (e.g., the former can be smaller than the latter). We present a novel deep convolution neural network (CNN)-based algorithm featuring high-fidelity recovery of multiscale feature size to improve lateral resolution of AR-PAM. The CNN is trained with simulated image pairs of line patterns, which is to mimic blood vessels. To investigate the suitable CNN model structure and elaborate on the effectiveness of CNN methods compared with non-learning methods, we select five different CNN models, while R-L and directional MB methods are also applied for comparison. Besides simulated data, experimental data including tungsten wires, leaf veins, and in vivo blood vessels are also evaluated. A custom-defined metric of relative size error (RSE) is used to quantify the multiscale feature recovery ability of different methods. Compared to other methods, enhanced deep super resolution (EDSR) network and residual in residual dense block network (RRDBNet) model show better recovery in terms of RSE for tungsten wires with diameters ranging from 30 μ m to 120 μ m . Moreover, AR-PAM images of leaf veins are tested to demonstrate the effectiveness of the optimized CNN methods (by EDSR and RRDBNet) for complex patterns. Finally, in vivo images of mouse ear blood vessels and rat ear blood vessels are acquired and then deconvolved, and the results show that the proposed CNN method (notably RRDBNet) enables accurate deconvolution of multiscale feature size and thus good fidelity.

Achieving depth-independent lateral resolution in AR-PAM using the synthetic-aperture focusing technique

Acoustic-resolution photoacoustic microscopy (AR-PAM) is a promising imaging modality that renders images with ultrasound resolution and extends the imaging depth beyond the optical ballistic regime. To achieve a high lateral resolution, a large numerical aperture (NA) of a focused transducer is usually applied for AR-PAM. However, AR-PAM fails to hold its performance in the out-of-focus region. The lateral resolution and signal-to-noise ratio (SNR) degrade substantially, thereby leading to a significantly deteriorated image quality outside the focal area. Based on the concept of the synthetic-aperture focusing technique (SAFT), various strategies have been developed to address this challenge. These include 1D-SAFT, 2D-SAFT, adaptive-SAFT, spatial impulse response (SIR)-based schemes, and delay-multiply-and-sum (DMAS) strategies. These techniques have shown progress in achieving depth-independent lateral resolution, while several challenges remain. This review aims to introduce these developments in SAFT-based approaches, highlight their fundamental mechanisms, underline the advantages and limitations of each approach, and discuss the outlook of the remaining challenges for future advances.

Image reconstruction for endoscopic photoacoustic tomography including effects of detector responses

In photoacoustic tomography (PAT), conventional image reconstruction methods are generally based on the assumption of an ideal point-like ultrasonic detector. This assumption is appropriate when the receiving surface of the detector is sufficiently small and/or the distance between the imaged object and the detector is large enough. However, it does not hold in endoscopic applications of PAT. In this study, we propose a model-based image reconstruction method for endoscopic photoacoustic tomography (EPAT), considering the effect of detector responses on image quality. We construct a forward model to physically describe the imaging process of EPAT, including the generation of the initial pressure due to optical absorption and thermoelastic expansion, the propagation of photoacoustic waves in tissues, and the acoustic measurement. The model outputs the theoretical sampling voltage signal, which is the response of the ultrasonic detector to the acoustic pressure reaching its receiving surface. The images representing the distribution map of the optical absorption energy density on cross-sections of the imaged luminal structures are reconstructed from the sampling voltage signals output by the detector through iterative inversion of the forward model. Compared with the conventional approaches based on back-projection and other imaging models, our method improved the quality and spatial resolution of the resulting images.

3-D ultrasonic image reconstruction in frequency domain using a virtual transducer model

In ultrasonic non-destructive testing, image reconstruction is essential to restore the diffracted ultrasound signals to improve the lateral resolution of images. Some reconstruction methods, like DAS-based synthetic aperture imaging, are inefficient, especially for reconstructing three-dimensional (3-D) images. Other methods do not provide high-resolution results, because they neglect the distortion effect introduced by transducer geometry. To overcome these disadvantages, we propose a 3-D ultrasonic image reconstruction method based on synthetic aperture wavenumber algorithm. It considers wave diffraction and transducer geometry effects, and can refocus the reflectors even in non-focal zone, which suits for large depth range imaging. This method builds a virtual transducer model in frequency domain by treating the focused transducer as a virtual planar transducer on its focal plane. In addition, the method uses non-uniform fast Fourier transform and deconvolution operation to achieve the 3-D image reconstruction, which has remarkably improved the efficiency and accuracy. According to the experimental results, the lateral resolution of an image reconstructed by the proposed method can reach 290.2 μm, exceeding the lateral resolution limitation of the 15 MHz focused transducer (523.24 μm). Furthermore, the proposed method only takes 0.744 s to reconstruct a 3-D image with 1000×100×100 pixels, while the time domain SAFT takes about 1163.8 s. It shows the potential for real-time 3-D imaging under advanced hardware.

An Excitation-Reception Collinear Probe for Ultrasonic, Photoacoustic, and Thermoacoustic Tri-Modal Volumetric Imaging

Imaging systems that integrate multiple modalities can reveal complementary anatomic and functional information as they exploit different contrast mechanisms, which have shown great application potential and advantages in preclinical studies. A portable and easy-to-use imaging probe will be more conducive to transfer to clinical practice. Here, we present a tri-modal ultrasonic (US), photoacoustic (PA), and thermoacoustic (TA) imaging system with an excitation-reception collinear probe. The acoustic field, light field, and electric field of the probe were designed to be coaxial, realizing homogeneous illumination and high-sensitivity detection at the same detection position. US images can provide detailed information about structures, PA images can delineate the morphology of blood vessels in tissues, and TA images can reveal dielectric properties of the tissues. Moreover, phantoms and in vivo human finger experiments were performed by the tri-modal imaging system to demonstrate its performance. The results show that the tri-modal imaging system with the proposed probe has the ability to detect small breast tumors with a radius of only 2.5 mm and visualize the anatomical structure of the finger in three dimensions. Our work confirms that the tri-modal imaging system equipped with a collinear probe can be applied to a variety of different scenarios, which lays a solid foundation for the application of the tri-modality system in clinical trials.

Super-Resolution Photoacoustic Microscopy Using Structured-Illumination

A novel super-resolution volumetric photoacoustic microscopy, based on the theory of structured-illumination, is proposed in this paper. The structured-illumination will be introduced in order to surpass the diffraction limit in a photoacoustic microscopy (PAM) structure. Through optical excitation of the targeted object with a sinusoidal spatial fringe pattern, the object's frequency spectrum is forced to shift in the spatial frequency domain. The shifting in the desired direction leads to the passage of the high-frequency contents of the object through the passband of the acoustic diffraction frequency response. Finally, combining the low-frequency image with the high-frequency parts in four regular orientations in the spatial frequency domain is equivalent to imaging the targeted object with an imaging system of two-fold bandwidth and thus half lateral resolution. In order to obtain the image of out-of-focus regions and improve the lateral resolution outside the focal region of a PAM imaging system, Fourier-domain reconstruction algorithm based on the synthetic aperture focusing technique (SAFT) using the virtual detector concept is utilized for reduction in the required computational load and time. The performance of the proposed imaging system is validated with in vivo and ex vivo targets. The experimental results obtained from several tungsten filaments in the depth range of 1.2 mm, show an improvement of -6 dB lateral resolution from 55- [Formula: see text] to 25- [Formula: see text] and also an improvement of signal-to-noise ratio (SNR) from 16-22 dB to 27-33 dB in the proposed system.

Photoacoustic imaging in monitoring of compartmental syndrome in rat extremities

Muscle ischemia injury is the essence of compartment syndrome (CS). Photoacoustic (PA) imaging can monitor hemoglobin concentration changes in ischemic tissue by determining the state of light-absorbing molecules. This study investigated whether PA imaging can provide accurate CS monitoring. Rats received compression on the lower hind limb for 3 h to induce ischemia injury, followed by PA imaging of desired muscles for 24 h. PA intensities of the injured group were significantly lower than that in the control group. Histology findings correlated well with the PA findings. The results demonstrated that PA imaging could be a noninvasive and timely tool for clinically monitoring CS.

Evaluation of Reconstruction Methodology for Helical Scan Guided Photoacoustic Endoscopy

Photoacoustic endoscopy (PAE), combining both advantages of optical contrast and acoustic resolution, can visualize the chemical-specific optical information of tissues inside human-body. Recently, its corresponding reconstruction methods have been extensively researched. However, most of them are limited on cylindrical scan trajectories, rather than a helical scan which is more clinically practical. On this note, this article proposes a methodology of imaging reconstruction and evaluation for helical scan guided PAE. Different from traditional reconstruction method, synthetic aperture focusing technique (SAFT), our method reconstructs image using wavefield extrapolation which significantly improves computational efficiency and even takes only 0.25 seconds for 3-D reconstructions. In addition, the proposed evaluation methodology can estimate the resolutions and deviations of reconstructed images in advance, and then can be used to optimize the PAE scan parameters. Groups of simulations as well as ex-vivo experiments with different scan parameters are provided to fully demonstrate the performance of the proposed techniques. The quantitatively measured angular resolutions and deviations agree well with our theoretical derivation results D√{r_s² +h²} / [1.25(r_s r_d +h²)] (rad) and -h l / (r_s r_d +h²) (rad), respectively D,r_d, r_s,h and l represent transducer diameter, radius of scan trajectory, radius of source position, unit helical pitch and the distance from targets to helical scan plane, respectively). This theoretical result also suits for circular and cylindrical scan in case of h = 0 .

Three-dimensional synthetic aperture focusing photoacoustic microscopy based on the acoustic simulation generated delay time and weighted factor

Acoustic resolution photoacoustic microscopy (ARPAM) is a promising imaging tool in biomedical applications for its advantage of penetration over other optical imaging techniques. However, the lateral resolution of ARPAM deteriorates significantly in the out-of-focus region. The synthetic aperture focusing technique (SAFT) is required to restore this kind of focus-related imaging distortion. The conventional SAFT method is based on the virtual detector (VD) conception, in which the phase of the received photoacoustic (PA) signal is calculated by assuming the focus of the transducer as a VD. Nevertheless, the phase of the received PA signal is not only determined by the geometrical parameters of the transducer, but also by the transducer's electromechanic response and the original PA signal. Ignoring these two factors will reduce the quality of the imaging results. In this work, a new SAFT method, which is based on acoustic simulation, is proposed for ARPAM. The measured PA signal from a point target at the focus is employed to evaluate the convolution of the transducer's electromechanic response and the original PA signal. This measured signal is used as the excitation in an acoustic simulation. The simulation, which is based on the geometrical and acoustic parameters of the transducer, is employed to calculate the delay time and weighted coefficient for the SAFT calculation. The phantom experiments with point and line targets indicate that the proposed method obtains imaging results with better lateral resolution and improved signal-noise ratio compared with the widely used VD-based SAFT method.

Low-Cost Multi-Wavelength Photoacoustic Imaging Based on Portable Continuous-Wave Laser Diode Module

Photoacoustic imaging (PAI), an emerging imaging technique, exploits the merits of both optical and ultrasound imaging, equipped with optical contrast and deep penetration. Typical linear PAI relies on a nanosecond laser pulse to induce photoacoustic signals. To construct a multi-wavelength PAI system, a multi-wavelength nano-second laser source is required, which greatly increases the cost of the PAI system. However, according to the nonlinear photoacoustic effect, the amplitude of the photoacoustic signals will vary with different base temperatures of the tissue. Therefore, using continuous-wave lasers with different wavelengths to induce different temperature variations at the same point of the tissue, and then using a single-wavelength pulsed laser to induce photoacoustic signals has been an alternative method to achieve multi-wavelength PAI. In this paper, based on the nonlinear photoacoustic effect, we developed a continuous-wave multi-wavelength laser source to cut down the cost of the conventional multi-wavelength PAI system. The principle will be introduced firstly, followed by qualitative and quantitative experiments.

Rapid Three-Dimensional Photoacoustic Imaging Reconstruction for Irregularly Layered Heterogeneous Media

Photoacoustic imaging (PAI) is susceptible to speed of sound (SOS) differences in heterogeneous media which greatly reduce the resolutions and qualities of the imaging results. Several reconstruction methods have been reported to adapt for heterogenous media, but they are limited by specific deficiencies such as efficiency, accuracy, and model limitation problems. Among them, the plane wave model based on wavefield reconstruction is the most efficient and promising one for high-efficiency three-dimensional PAI. However, the classic plane wave model only suits for planar layered media, severely limiting its applications in practice. To this end, we modify the plane wave model to apply for irregularly layered heterogeneous media and propose a corresponding wavefield extrapolation to reconstruct photoacoustic image. This method employs split-step Fourier to compensate the SOS differences, extrapolates wavefields and reconstructs the image depth by depth. Furthermore, a floating discretization strategy is introduced to control and balance the efficiency and accuracy with a hyperparameter. The simulation and experiment results demonstrate that the proposed method can reconstruct the image with an equivalent resolution to time reversal's and even have higher efficiency and robustness. To reconstruct a three-dimensional image with 50×50×600 pixels, the proposed method takes only 5.5 seconds using a laptop loaded with Intel(R) Core (TM) i7-8550U CPU @1.8GHz.

Reducing Acoustic Inhomogeneity Based on Speed of Sound Autofocus in Microwave Induced Thermoacoustic Tomography

Microwave induced thermoacoustic tomography is a newly developing non-invasive and non-ionizing modality. In practical applications, such as breast tumor detection and brain imaging, the acoustic properties in the tissue to be detected are usually unknown and spatially non-uniform, which results in distortion and blurring of the buried targets. In this paper, a reconstruction method based on speed of sound (SoS) autofocus is proposed to reduce the effect of acoustic inhomogeneity in different soft tissues. According to this method, the number of tissue types, which are referred to as clusters in this work, can be automatically determined by a decision graph. To distinguish the boundaries of different tissues, a Gaussian Mixture Model (GMM) is fitted to the obtained image data for soft clustering instead of traditional hard clustering. Through fixing the tissue centers which are characterized by corresponding data density peaks as the means of Gaussian parameters rather than choosing them randomly, adaptive and robust reconstruction performance can be guaranteed. After performing an iterative GMM optimization, the SoS autofocus is achieved. Image reconstructed by using the updated SoS distribution is with higher accuracy than that with homogeneous assumption. Compared with the existing similar methods, the proposed method strategy obviates the need of extra experiment costs, and possesses good robustness with respect to hard assignment model errors when the medium is relatively complex. Realistic breast model and brain model simulations combined with experiments of agar phantom and pig's brain are provided to demonstrate the effectiveness of the proposed method.

Toward Wearable Healthcare: A Miniaturized 3D Imager With Coherent Frequency-Domain Photoacoustics

Medical monitoring is undergoing a translation from the hospital-based system to the personalized home-based system. With the aim of wearable application of photoacoustic technique, we propose a miniaturized photoacoustic 3D imager for superficial medical imaging. By employing the compact continuous-wave laser diode based optical irradiation and an ultrathin 2D matrix array based photoacoustic detection in the coherent frequency domain, a wearable imaging probe with a size of about 80 × 25 × 24 mm³ and a weight of 21 g is developed. At the backend, an FPGA controlled Howland current source drives the laser diodes to excite linear frequency modulated optical irradiation. Recorded by a portable multichannel data acquisition system, the generated photoacoustic responses are firstly compressed with the coherent frequency domain photoacoustic method and then extrapolated in the wavenumber-frequency domain for fast image reconstruction. With three-wavelength (450 nm, 638 nm, and 808 nm) laser irradiation, photoacoustic imaging can be operated multispectrally, endowing the developed imager with functional imaging capability in 3D space. With the imager worn on the human forearm, hemoglobin oxygen saturation level in superficial arm vasculature can be long-term monitored with high stability. When the imager is applied for imaging in a relatively large area (e.g., early melanoma detection in the human breast), flexible scanning in a handheld manner can be performed. This work opens the application potential of photoacoustic technique in a broad range of areas, including personalized healthcare, home health monitoring, and long-term physiologic monitoring.