Model-Based Reconstruction of Large Three-Dimensional Optoacoustic Datasets

Compressed Sensing for Biomedical Photoacoustic Imaging: A Review

Photoacoustic imaging (PAI) is a rapidly developing emerging non-invasive biomedical imaging technique that combines the strong contrast from optical absorption imaging and the high resolution from acoustic imaging. Abnormal biological tissues (such as tumors and inflammation) generate different levels of thermal expansion after absorbing optical energy, producing distinct acoustic signals from normal tissues. This technique can detect small tissue lesions in biological tissues and has demonstrated significant potential for applications in tumor research, melanoma detection, and cardiovascular disease diagnosis. During the process of collecting photoacoustic signals in a PAI system, various factors can influence the signals, such as absorption, scattering, and attenuation in biological tissues. A single ultrasound transducer cannot provide sufficient information to reconstruct high-precision photoacoustic images. To obtain more accurate and clear image reconstruction results, PAI systems typically use a large number of ultrasound transducers to collect multi-channel signals from different angles and positions, thereby acquiring more information about the photoacoustic signals. Therefore, to reconstruct high-quality photoacoustic images, PAI systems require a significant number of measurement signals, which can result in substantial hardware and time costs. Compressed sensing is an algorithm that breaks through the Nyquist sampling theorem and can reconstruct the original signal with a small number of measurement signals. PAI based on compressed sensing has made breakthroughs over the past decade, enabling the reconstruction of low artifacts and high-quality images with a small number of photoacoustic measurement signals, improving time efficiency, and reducing hardware costs. This article provides a detailed introduction to PAI based on compressed sensing, such as the physical transmission model-based compressed sensing method, two-stage reconstruction-based compressed sensing method, and single-pixel camera-based compressed sensing method. Challenges and future perspectives of compressed sensing-based PAI are also discussed.

Multi-Scale Volumetric Dynamic Optoacoustic and Laser Ultrasound (OPLUS) Imaging Enabled by Semi-Transparent Optical Guidance

Major biological discoveries are made by interrogating living organisms with light. However, the limited penetration of un-scattered photons within biological tissues limits the depth range covered by optical methods. Deep-tissue imaging is achieved by combining light and ultrasound. Optoacoustic imaging exploits the optical generation of ultrasound to render high-resolution images at depths unattainable with optical microscopy. Recently, laser ultrasound has been suggested as a means of generating broadband acoustic waves for high-resolution pulse-echo ultrasound imaging. Herein, an approach is proposed to simultaneously interrogate biological tissues with light and ultrasound based on layer-by-layer coating of silica optical fibers with a controlled degree of transparency. The time separation between optoacoustic and ultrasound signals collected with a custom-made spherical array transducer is exploited for simultaneous 3D optoacoustic and laser ultrasound (OPLUS) imaging with a single laser pulse. OPLUS is shown to enable large-scale anatomical characterization of tissues along with functional multi-spectral imaging of chromophores and assessment of cardiac dynamics at ultrafast rates only limited by the pulse repetition frequency of the laser. The suggested approach provides a flexible and scalable means for developing a new generation of systems synergistically combining the powerful capabilities of optoacoustics and ultrasound imaging in biology and medicine.

Deep learning framework for three-dimensional surface reconstruction of object of interest in photoacoustic tomography

Photoacoustic tomography (PAT) is a non-ionizing hybrid imaging technology of clinical importance that combines the high contrast of optical imaging with the high penetration of ultrasonic imaging. Two-dimensional (2D) tomographic images can only provide the cross-sectional structure of the imaging target rather than its overall spatial morphology. This work proposes a deep learning framework for reconstructing three-dimensional (3D) surface of an object of interest from a series of 2D images. It achieves end-to-end mapping from a series of 2D images to a 3D image, visually displaying the overall morphology of the object. The framework consists of four modules: segmentation module, point cloud generation module, point cloud completion module, and mesh conversion module, which respectively implement the tasks of segmenting a region of interest, generating a sparse point cloud, completing sparse point cloud and reconstructing 3D surface. The network model is trained on simulation data sets and verified on simulation, phantom, and in vivo data sets. The results showed superior 3D reconstruction performance both visually and on the basis of quantitative evaluation metrics compared to the state-of-the-art non-learning and learning approaches. This method potentially enables high-precision 3D surface reconstruction from the tomographic images output by the preclinical PAT system without changing the imaging system. It provides a general deep learning scheme for 3D reconstruction from tomographic scanning data.

Dynamic modeling of photoacoustic sensor data to classify human blood samples

The photoacoustic effect is an attractive tool for diagnosis in several biomedical applications. Analyzing photoacoustic signals, however, is challenging to provide qualitative results in an automated way. In this work, we introduce a dynamic modeling scheme of photoacoustic sensor data to classify blood samples according to their physiological status. Thirty-five whole human blood samples were studied with a state-space model estimated by a subspace method. Furthermore, the samples are classified using the model parameters and the linear discriminant analysis algorithm. The classification performance is compared with time- and frequency-domain features and an autoregressive-moving-average model. As a result, the proposed analysis can predict five blood classes: healthy women and men, microcytic and macrocytic anemia, and leukemia. Our findings indicate that the proposed method outperforms conventional signal processing techniques to analyze photoacoustic data for medical diagnosis. Hence, the method is a promising tool in point-of-care devices to detect hematological diseases in clinical scenarios.

Signal domain adaptation network for limited-view optoacoustic tomography

Optoacoustic (OA) imaging is based on optical excitation of biological tissues with nanosecond-duration laser pulses and detection of ultrasound (US) waves generated by thermoelastic expansion following light absorption. The image quality and fidelity of OA images critically depend on the extent of tomographic coverage provided by the US detector arrays. However, full tomographic coverage is not always possible due to experimental constraints. One major challenge concerns an efficient integration between OA and pulse-echo US measurements using the same transducer array. A common approach toward the hybridization consists in using standard linear transducer arrays, which readily results in arc-type artifacts and distorted shapes in OA images due to the limited angular coverage. Deep learning methods have been proposed to mitigate limited-view artifacts in OA reconstructions by mapping artifactual to artifact-free (ground truth) images. However, acquisition of ground truth data with full angular coverage is not always possible, particularly when using handheld probes in a clinical setting. Deep learning methods operating in the image domain are then commonly based on networks trained on simulated data. This approach is yet incapable of transferring the learned features between two domains, which results in poor performance on experimental data. Here, we propose a signal domain adaptation network (SDAN) consisting of i) a domain adaptation network to reduce the domain gap between simulated and experimental signals and ii) a sides prediction network to complement the missing signals in limited-view OA datasets acquired from a human forearm by means of a handheld linear transducer array. The proposed method showed improved performance in reducing limited-view artifacts without the need for ground truth signals from full tomographic acquisitions.

Deep learning on photoacoustic tomography to remove image distortion due to inaccurate measurement of the scanning radius

Photoacoustic tomography (PAT) is a non-invasive, non-ionizing hybrid imaging modality that holds great potential for various biomedical applications and the incorporation with deep learning (DL) methods has experienced notable advancements in recent times. In a typical 2D PAT setup, a single-element ultrasound detector (USD) is used to collect the PA signals by making a 360° full scan of the imaging region. The traditional backprojection (BP) algorithm has been widely used to reconstruct the PAT images from the acquired signals. Accurate determination of the scanning radius (SR) is required for proper image reconstruction. Even a slight deviation from its nominal value can lead to image distortion compromising the quality of the reconstruction. To address this challenge, two approaches have been developed and examined herein. The first framework includes a modified version of dense U-Net (DUNet) architecture. The second procedure involves a DL-based convolutional neural network (CNN) for image classification followed by a DUNet. The first protocol was trained with heterogeneous simulated images generated from three different phantoms to learn the relationship between the reconstructed and the corresponding ground truth (GT) images. In the case of the second scheme, the first stage was trained with the same heterogeneous dataset to classify the image type and the second stage was trained individually with the appropriate images. The performance of these architectures has been tested on both simulated and experimental images. The first method can sustain SR deviation up to approximately 6% for simulated images and 5% for experimental images and can accurately reproduce the GTs. The proposed DL-approach extends the limits further (approximately 7% and 8% for simulated and experimental images, respectively). Our results suggest that classification-based DL method does not need a precise assessment of SR for accurate PAT image formation.

Towards quantitative in vivo dosimetry using x-ray acoustic computed tomography

BACKGROUND

Radiation dosimetry is essential for radiation therapy (RT) to ensure that radiation dose is accurately delivered to the tumor. Despite its wide use in clinical intervention, the delivered radiation dose can only be planned and verified via simulation. This makes precision radiotherapy challenging while in-line verification of the delivered dose is still absent in the clinic. X-ray-induced acoustic computed tomography (XACT) has recently been proposed as an imaging tool for in vivo dosimetry.

PURPOSE

Most of the XACT studies focus on localizing the radiation beam. However, it has not been studied for its potential for quantitative dosimetry. The aim of this study was to investigate the feasibility of using XACT for quantitative in vivo dose reconstruction during radiotherapy.

METHODS

Varian Eclipse system was used to generate simulated uniform and wedged 3D radiation field with a size of 4 cm× $ \times \ $ 4 cm. In order to use XACT for quantitative dosimetry measurements, we have deconvoluted the effects of both the x-ray pulse shape and the finite frequency response of the ultrasound detector. We developed a model-based image reconstruction algorithm to quantify radiation dose in vivo using XACT imaging, and universal back-projection (UBP) reconstruction is used as comparison. The reconstructed dose was calibrated before comparing it to the percent depth dose (PDD) profile. Structural similarity index matrix (SSIM) and root mean squared error (RMSE) are used for numeric evaluation. Experimental signals were acquired from 4 cm× $ \times \ $ 4 cm radiation field created by Linear Accelerator (LINAC) at depths of 6, 8, and 10 cm beneath the water surface. The acquired signals were processed before reconstruction to achieve accurate results.

RESULTS

Applying model-based reconstruction algorithm with non-negative constraints successfully reconstructed accurate radiation dose in 3D simulation study. The reconstructed dose matches well with the PDD profile after calibration in experiments. The SSIMs between the model-based reconstructions and initial doses are over 85%, and the RMSEs of model-based reconstructions are eight times lower than the UBP reconstructions. We have also shown that XACT images can be displayed as pseudo-color maps of acoustic intensity, which correspond to different radiation doses in the clinic.

CONCLUSION

Our results show that the XACT imaging by model-based reconstruction algorithm is considerably more accurate than the dose reconstructed by UBP algorithm. With proper calibration, XACT is potentially applicable to the clinic for quantitative in vivo dosimetry across a wide range of radiation modalities. In addition, XACT's capability of real-time, volumetric dose imaging seems well-suited for the emerging field of ultrahigh dose rate "FLASH" radiotherapy.

Full-view LED-based optoacoustic tomography

Optoacoustic tomography is commonly performed with bulky and expensive short-pulsed solid-state lasers providing high per-pulse energies in the millijoule range. Light emitting diodes (LEDs) represent a cost-effective and portable alternative for optoacoustic signal excitation that can additionally provide excellent pulse-to-pulse stability. Herein, we introduce a full-view LED-based optoacoustic tomography (FLOAT) system for deep tissue in vivo imaging. It is based on a custom-made electronic unit driving a stacked array of LEDs, which attains 100 ns pulse width and highly stable (0.62 % standard deviation) total per-pulse energy of 0.48 mJ. The illumination source is integrated into a circular array of cylindrically-focused ultrasound detection elements to result in a full-view tomographic configuration, which plays a critical role in circumventing limited-view effects, enhancing the effective field-of-view and image quality for cross-sectional (2D) imaging. We characterized the FLOAT performance in terms of pulse width, power stability, excitation light distribution, signal-to-noise and penetration depth. FLOAT of the human finger revealed a comparable imaging performance to that achieved with the standard pulsed Nd:YAG laser. It is anticipated that this compact, affordable and versatile illumination technology will facilitate optoacoustic imaging developments in resource-limited settings for biological and clinical applications.

Model-Based 3-D X-Ray Induced Acoustic Computerized Tomography

X-ray-induced acoustic (XA) computerized tomography (XACT) is an evolving imaging technique that aims to reconstruct the X-ray energy deposition from XA measurements. Main challenges in XACT are the poor signal-to-noise ratio and limited field-of-view, which cause artifacts in the images. We demonstrate the efficacy of model-based (MB) algorithms for three-dimensional XACT and compare with the traditional algorithms. The MB algorithm is based on iterative, matrix-free approach for regularized-least-squares minimization corresponding to XACT. The matrix-free-LSQR (MF-LSQR) and the non-iterative model-backprojection (MBP) reconstructions were evaluated and compared with universal backprojection (UBP), time-reversal (TR) and fast-Fourier transform (FFT)-based reconstructions for numerical and experimental XACT datasets. The results demonstrate the capability of MF-LSQR algorithm to reduce noisy artifacts thus yielding better reconstructions. MBP and MF-LSQR algorithms perform particularly well with the experimental XACT dataset, where noise in signals significantly affects the reconstruction of the target in UBP and FFT-based reconstructions. The TR reconstruction for experimental XACT are comparable to MF-LSQR, but takes thrice as much time and filters the frequency components greater than maximum frequency supported by the grid, resulting loss of resolution. The MB algorithms are able to overcome the challenges in XACT and hence are vital for the clinical translation of XACT.

Signal restoration algorithm for photoacoustic imaging systems

In a photoacoustic (PA) imaging system, the detectors are bandwidth-limited. Therefore, they capture PA signals with some unwanted ripples. This limitation degrades the resolution/contrast and induces sidelobes and artifacts in the reconstructed images along the axial direction. To compensate for the limited bandwidth effect, we present a PA signal restoration algorithm, where a mask is designed to extract the signals at the absorber positions and remove the unwanted ripples. This restoration improves the axial resolution and contrast in the reconstructed image. The restored PA signals can be considered as the input of the conventional reconstruction algorithms (e.g., Delay-and-sum (DAS) and Delay-multiply-and-sum (DMAS)). To compare the performance of the proposed method, DAS and DMAS reconstruction algorithms were performed with both the initial and restored PA signals on numerical and experimental studies (numerical targets, tungsten wires, and human forearm). The results show that, compared with the initial PA signals, the restored PA signals can improve the axial resolution and contrast by 45% and 16.1 dB, respectively, and suppress background artifacts by 80%.

Enhancing optoacoustic mesoscopy through calibration-based iterative reconstruction

Optoacoustic mesoscopy combines rich optical absorption contrast with high spatial resolution at tissue depths beyond reach for microscopic techniques employing focused light excitation. The mesoscopic imaging performance is commonly hindered by the use of inaccurate delay-and-sum reconstruction approaches and idealized modeling assumptions. In principle, image reconstruction performance could be enhanced by simulating the optoacoustic signal generation, propagation, and detection path. However, for most realistic experimental scenarios, the underlying total impulse response (TIR) cannot be accurately modelled. Here we propose to capture the TIR by scanning of a sub-resolution sized absorber. Significant improvement of spatial resolution and depth uniformity is demonstrated over 3 mm range, outperforming delay-and-sum and model-based reconstruction implementations. Reconstruction performance is validated by imaging subcutaneous murine vasculature and human skin in vivo. The proposed experimental calibration and reconstruction paradigm facilitates quantitative inversions while averting complex physics-based simulations. It can readily be applied to other imaging modalities employing TIR-based reconstructions.

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.

Learned regularization for image reconstruction in sparse-view photoacoustic tomography

Constrained data acquisitions, such as sparse view measurements, are sometimes used in photoacoustic computed tomography (PACT) to accelerate data acquisition. However, it is challenging to reconstruct high-quality images under such scenarios. Iterative image reconstruction with regularization is a typical choice to solve this problem but it suffers from image artifacts. In this paper, we present a learned regularization method to suppress image artifacts in model-based iterative reconstruction in sparse view PACT. A lightweight dual-path network is designed to learn regularization features from both the data and the image domains. The network is trained and tested on both simulation and in vivo datasets and compared with other methods such as Tikhonov regularization, total variation regularization, and a U-Net based post-processing approach. Results show that although the learned regularization network possesses a size of only 0.15% of a U-Net, it outperforms other methods and converges after as few as five iterations, which takes less than one-third of the time of conventional methods. Moreover, the proposed reconstruction method incorporates the physical model of photoacoustic imaging and explores structural information from training datasets. The integration of deep learning with a physical model can potentially achieve improved imaging performance in practice.

Expediting Image Acquisition in Reflection Ultrasound Computed Tomography

Reflection ultrasound computed tomography (RUCT) attains optimal image quality from objects that can be fully accessed from multiple directions, such as the human breast or small animals. Owing to the full-view tomography approach based on the compounding of images taken from multiple angles, RUCT effectively mitigates several deficiencies afflicting conventional pulse-echo ultrasound (US) systems, such as speckle patterns and interuser variability. On the other hand, the small interelement pitch required to fulfill the spatial sampling criterion in the circular transducer configuration used in RUCT typically implies the use of an excessive number of independent array elements. This increases the system's complexity and costs, and limits the achievable imaging speed. Here, we explore acquisition schemes that enable RUCT imaging with the reduced number of transmit/receive elements. We investigated the influence of the element size in transmission and reception in a ring array geometry. The performance of a sparse acquisition approach based on partial acquisition from a subset of the elements has been further assessed. A larger element size is shown to preserve contrast and resolution at the center of the field of view (FOV), while a reduced number of elements is shown to cause uniform loss of contrast and resolution across the entire FOV. The tradeoffs of achievable FOV, contrast-to-noise ratio, and temporal and spatial resolutions are assessed in phantoms and in vivo mouse experiments. The experimental analysis is expected to aid the development of optimized hardware and image acquisition strategies for RUCT and, thus, result in more affordable imaging systems facilitating wider adoption.

Multi-angle data acquisition to compensate transducer finite size in photoacoustic tomography

In photoacoustic tomography (PAT) systems, the tangential resolution decreases due to the finite size of the transducer as the off-center distance increases. To address this problem, we propose a multi-angle detection approach in which the transducer used for data acquisition rotates around its center (with specific angles) as well as around the scanning center. The angles are calculated based on the central frequency and diameter of the transducer and the radius of the region-of-interest (ROI). Simulations with point-like absorbers (for point-spread-function evaluation) and a vasculature phantom (for quality assessment), and experiments with ten 0.5 mm-diameter pencil leads and a leaf skeleton phantom are used for evaluation of the proposed approach. The results show that a location-independent tangential resolution is achieved with 150 spatial sampling and central rotations with angles of ±8°/±16°. With further developments, the proposed detection strategy can replace the conventional detection (rotating a transducer around ROI) in PAT.

Model-Based X-Ray-Induced Acoustic Computed Tomography

X-ray-induced acoustic computed tomography (XACT) provides X-ray absorption-based contrast with acoustic detection. For its clinical translation, XACT imaging often has a limited field of view. This can result in image artifacts and overall loss of quantification accuracy. In this article, we aim to demonstrate model-based XACT image reconstruction to address these problems. An efficient matrix-free implementation of the regularized LSQR (MF-LSQR)-based minimization scheme and a noniterative model back-projection (MBP) scheme for computing XACT reconstructions have been demonstrated in this article. The proposed algorithms have been numerically validated and then used to perform reconstructions from experimental measurements obtained from an XACT setup. While the commonly used back-projection (BP) algorithm produces limited-view and noisy artifacts in the region of interest (ROI), model-based LSQR minimization overcomes these issues. The model-based algorithms also reduce the ring artifacts caused due to the nonuniformity response of the multichannel data acquisition. Using the model-based reconstruction algorithms, we are able to obtain reasonable XACT reconstructions for acoustic measurements of up to 120° view. Although the MBP is more efficient than the model-based LSQR algorithm, it provides only the structural information of the ROI. Overall, it has been demonstrated that the model-based image reconstruction yields better image quality for XACT than the standard BP. Moreover, the combination of model-based image reconstruction with different regularization methods can solve the limited-view problem for XACT imaging (in many realistic cases where the full-view dataset is unavailable), and hence pave the way for future clinical translation.

Photoacoustic computed tomography of mechanical HIFU-induced vascular injury

Mechanical high-intensity focused ultrasound (HIFU) has been used for cancer treatment and drug delivery. Existing monitoring methods for mechanical HIFU therapies such as MRI and ultrasound imaging often suffer from high cost, poor spatial-temporal resolution, and/or low sensitivity to tissue's hemodynamic changes. Evaluating vascular injury during mechanical HIFU treatment, therefore, remains challenging. Photoacoustic computed tomography (PACT) is a promising tool to meet this need. Intrinsically sensitive to optical absorption, PACT provides high-resolution imaging of blood vessels using hemoglobin as the endogenous contrast. In this study, we have developed an integrated HIFU-PACT system for detecting vascular rupture in mechanical HIFU treatment. We have demonstrated singular value decomposition for enhancing hemorrhage detection. We have validated the HIFU-PACT performance on phantoms and in vivo animal tumor models. We expect that PACT-HIFU will find practical applications in oncology research using small animal models.

Optoacoustic imaging of the skin

Optoacoustic (OA, photoacoustic) imaging capitalizes on the synergistic combination of light excitation and ultrasound detection to empower biological and clinical investigations with rich optical contrast while effectively bridging the gap between micro and macroscopic imaging realms. State-of-the-art OA embodiments consistently provide images at micron-scale resolution through superficial tissue layers by means of focused illumination that can be smoothly exchanged for acoustic-resolution images at diffuse light depths of several millimetres to centimetres via ultrasound beamforming or tomographic reconstruction. Taken together, this unique multi-scale imaging capacity opens unprecedented capabilities for high-resolution in vivo interrogations of the skin at scalable depths. Moreover, diverse anatomical and functional information is retrieved via dynamic mapping of endogenous chromophores such as haemoglobin, melanin, lipids, collagen, water and others. This, along with the use of non-ionizing radiation, facilitates a clinical translation of the OA modalities. We review recent progress in OA imaging of the skin in preclinical and clinical studies exploiting the rich contrast provided by endogenous substances in tissues. The imaging capabilities of existing approaches are discussed in the context of initial translational studies on skin cancer, inflammatory skin diseases, wounds and other conditions.

Internal-Illumination Photoacoustic Tomography Enhanced by a Graded-Scattering Fiber Diffuser

The penetration depth of photoacoustic imaging in biological tissues has been fundamentally limited by the strong optical attenuation when light is delivered externally through the tissue surface. To address this issue, we previously reported internal-illumination photoacoustic imaging using a customized radial-emission optical fiber diffuser, which, however, has complex fabrication, high cost, and non-uniform light emission. To overcome these shortcomings, we have developed a new type of low-cost fiber diffusers based on a graded-scattering method in which the optical scattering of the fiber diffuser is gradually increased as the light travels. The graded scattering can compensate for the optical attenuation and provide relatively uniform light emission along the diffuser. We performed Monte Carlo numerical simulations to optimize several key design parameters, including the number of scattering segments, scattering anisotropy factor, divergence angle of the optical fiber, and reflective index of the surrounding medium. These optimized parameters collectively result in uniform light emission along the fiber diffuser and can be flexibly adjusted to accommodate different applications. We fabricated and characterized the prototype fiber diffuser made of agarose gel and intralipid. Equipped with the new fiber diffuser, we performed thorough proof-of-concept studies on ex vivo tissue phantoms and an in vivo swine model to demonstrate the deep-imaging capability (~10 cm achieved ex vivo) of photoacoustic tomography. We believe that the internal light delivery via the optimized fiber diffuser is an effective strategy to image deep targets (e.g., kidney) in large animals or humans.

Sound Out the Deep Colors: Photoacoustic Molecular Imaging at New Depths

Photoacoustic tomography (PAT) has become increasingly popular for molecular imaging due to its unique optical absorption contrast, high spatial resolution, deep imaging depth, and high imaging speed. Yet, the strong optical attenuation of biological tissues has traditionally prevented PAT from penetrating more than a few centimeters and limited its application for studying deeply seated targets. A variety of PAT technologies have been developed to extend the imaging depth, including employing deep-penetrating microwaves and X-ray photons as excitation sources, delivering the light to the inside of the organ, reshaping the light wavefront to better focus into scattering medium, as well as improving the sensitivity of ultrasonic transducers. At the same time, novel optical fluence mapping algorithms and image reconstruction methods have been developed to improve the quantitative accuracy of PAT, which is crucial to recover weak molecular signals at larger depths. The development of highly-absorbing near-infrared PA molecular probes has also flourished to provide high sensitivity and specificity in studying cellular processes. This review aims to introduce the recent developments in deep PA molecular imaging, including novel imaging systems, image processing methods and molecular probes, as well as their representative biomedical applications. Existing challenges and future directions are also discussed.

Cross-sectional photoacoustic tomography image reconstruction with a multi-curve integration model

BACKGROUND AND OBJECTIVE

In acoustic inversion of photoacoustic tomography (PAT), an imaging model that precisely describes both the ultrasonic wave propagation and the detector properties is of crucial importance. Inspired by the multi-stripe integration model in clinical X-ray computed tomography systems, in this work, we introduce the Multi-Curve-Integration-based acoustic inversion for cross-sectional Photoacoustic Tomography (MCI-PAT).

METHODS

We assumed that in cross-sectional PAT system, the three-dimensional (3-D) wave propagation problem could be reduced to a two-dimensional (2-D) problem in a limited, yet sufficient field of view. Under such condition, the MCI-PAT imaging model is generated by integrating several circular acoustic curves, the centers of which are points evenly distributed on the finite-size ultrasonic transducer surface. In this way, the spatial detector response is taken into account, while its computational burden does not largely increase because the integration process is performed only on a 2-D plane.

RESULTS

As proven by simulation, phantom and in vivo small animal experiments, the MCI-PAT method leads to promising improvement in PAT image quality. Comparing to traditional imaging models that considered only a single acoustic curve, our proposed method successfully improved the visibility of small structures and achieved evident enhancement on signal-to-noise ratio.

CONCLUSIONS

The performance of the MCI-PAT method demonstrates that for cross-sectional PAT, a 2-D simplification of the propagation of multiple photoacoustic waves is feasible. Due to its simplicity, our method can be used as an add-on to current system models considering only a single acoustic curve.