Perturbation Monte Carlo Method for Quantitative Photoacoustic Tomography

Extractor-attention-predictor network for quantitative photoacoustic tomography

Quantitative photoacoustic tomography (qPAT) holds great potential in estimating chromophore concentrations, whereas the involved optical inverse problem, aiming to recover absorption coefficient distributions from photoacoustic images, remains challenging. To address this problem, we propose an extractor-attention-predictor network architecture (EAPNet), which employs a contracting-expanding structure to capture contextual information alongside a multilayer perceptron to enhance nonlinear modeling capability. A spatial attention module is introduced to facilitate the utilization of important information. We also use a balanced loss function to prevent network parameter updates from being biased towards specific regions. Our method obtains satisfactory quantitative metrics in simulated and real-world validations. Moreover, it demonstrates superior robustness to target properties and yields reliable results for targets with small size, deep location, or relatively low absorption intensity, indicating its broader applicability. The EAPNet, compared to the conventional UNet, exhibits improved efficiency, which significantly enhances performance while maintaining similar network size and computational complexity.

Application limits of the scaling relations for Monte Carlo simulations in diffuse optics. Part 1: theory

Monte Carlo (MC) is a powerful tool to study photon migration in scattering media, yet quite time-consuming to solve inverse problems. To speed up MC-simulations, scaling relations can be applied to an existing initial MC-simulation to generate a new data-set with different optical properties. We named this approach trajectory-based since it uses the knowledge of the detected photon trajectories of the initial MC-simulation, in opposition to the slower photon-based approach, where a novel MC-simulation is rerun with new optical properties. We investigated the convergence and applicability limits of the scaling relations, both related to the likelihood that the sample of trajectories considered is representative also for the new optical properties. For absorption, the scaling relation contains smoothly converging Lambert-Beer factors, whereas for scattering it is the product of two quickly diverging factors, whose ratio, for NIRS cases, can easily reach ten orders of magnitude. We investigated such instability by studying the probability-distribution for the number of scattering events in trajectories of given length. We propose a convergence test of the scattering scaling relation based on the minimum-maximum number of scattering events in recorded trajectories. We also studied the dependence of MC-simulations on optical properties, most critical in inverse problems, finding that scattering derivatives are ascribed to small deviations in the distribution of scattering events from a Poisson distribution. This paper, which can also serve as a tutorial, helps to understand the physics of the scaling relations with the causes of their limitations and devise new strategies to deal with them.

Carbon nanomaterial-involved EMT and CSC in cancer

Carbon nanomaterials (CNMs) are ubiquitous in our daily lives because of the outstanding physicochemical properties. CNMs play curial parts in industrial and medical fields, however, the risks of CNMs exposure to human health are still not fully understood. In view of, it is becoming extremely difficult to ignore the existence of the toxicity of CNMs. With the increasing exploitation of CNMs, it's necessary to evaluate the potential impact of these materials on human health. In recent years, more and more researches have shown that CNMs are contributed to the cancer formation and metastasis after long-term exposure through epithelial-mesenchymal transition (EMT) and cancer stem cells (CSCs) which is associated with cancer progression and invasion. This review discusses CNMs properties and applications in industrial and medical fields, adverse effects on human health, especially the induction of tumor initiation and metastasis through EMT and CSCs procedure.

Quantitative photoacoustic tomography with light fluence compensation based on radiance Monte Carlo model

Objective. Photoacoustic tomography (PAT) is a rapidly evolving imaging modality that provides images with high contrast and spatial resolution showing the optical properties of biological tissues. The photoacoustic pressure is proportional to the product of the optical absorption coefficient and the local light fluence. The essential challenge in reconstructing quantitative images representing spatially varying absorption coefficients is the unknown light fluence. In addition, optical attenuation induces spatial variations in the light fluence, and the heterogeneity of the fluence determines the limits of reconstruction quality and depth.Approach.In this work, a reconstruction enhancement scheme is proposed to compensate for the variation in the light fluence in the absorption coefficient recovery. The inverse problem of the radiance Monte Carlo model describing light transport through the tissue is solved by using an alternating optimization strategy. In the iteration, the absorption coefficients and photon weights are alternately updated.Main results.The method provides highly accurate quantitative images of absorption coefficients in simulations, phantoms, andin vivostudies. The results show that the method has great potential for improving the accuracy of absorption coefficient recovery compared to conventional reconstruction methods that ignore light fluence variations. Comparison with state-of-the-art fluence compensation methods shows significant improvements in root mean square error, normalized mean square absolute distance, and structural similarity metrics.Significance.This method achieves high precision quantitative imaging by compensating for nonuniform light fluence without increasing the complexity and operation of the imaging system.

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.

Adaptive stochastic Gauss-Newton method with optical Monte Carlo for quantitative photoacoustic tomography

SIGNIFICANCE

The image reconstruction problem in quantitative photoacoustic tomography (QPAT) is an ill-posed inverse problem. Monte Carlo method for light transport can be utilized in solving this image reconstruction problem.

AIM

The aim was to develop an adaptive image reconstruction method where the number of photon packets in Monte Carlo simulation is varied to achieve a sufficient accuracy with reduced computational burden.

APPROACH

The image reconstruction problem was formulated as a minimization problem. An adaptive stochastic Gauss-Newton (A-SGN) method combined with Monte Carlo method for light transport was developed. In the algorithm, the number of photon packets used on Gauss-Newton (GN) iteration was varied utilizing a so-called norm test.

RESULTS

The approach was evaluated with numerical simulations. With the proposed approach, the number of photon packets needed for solving the inverse problem was significantly smaller than in a conventional approach where the number of photon packets was fixed for each GN iteration.

CONCLUSIONS

The A-SGN method with a norm test can be utilized in QPAT to provide accurate and computationally efficient solutions.

High-resolution three-dimensional blood flow tomography in the subdiffuse regime using laser speckle contrast imaging

SIGNIFICANCE

Visualizing high-resolution hemodynamics in cerebral tissue over a large field of view (FOV), provides important information in studying disease states affecting the brain. Current state-of-the-art optical blood flow imaging techniques either lack spatial resolution or are too slow to provide high temporal resolution reconstruction of flow map over a large FOV.

AIM

We present a high spatial resolution computational optical imaging technique based on principles of laser speckle contrast imaging (LSCI) for reconstructing the blood flow maps in complex tissue over a large FOV provided that the three-dimensional (3D) vascular structure is known or assumed.

APPROACH

Our proposed method uses a perturbation Monte Carlo simulation of the high-resolution 3D geometry for both accurately deriving the speckle contrast forward model and calculating the Jacobian matrix used in our reconstruction algorithm to achieve high resolution. Given the convex nature of our highly nonlinear problem, we implemented a mini-batch gradient descent with an adaptive learning rate optimization method to iteratively reconstruct the blood flow map. Specifically, we implemented advanced optimization techniques combined with efficient parallelization and vectorization of the forward and derivative calculations to make reconstruction of the blood flow map feasible with reconstruction times on the order of tens of minutes.

RESULTS

We tested our reconstruction algorithm through simulation of both a flow phantom model as well as an anatomically correct murine cerebral tissue and vasculature captured via two-photon microscopy. Additionally, we performed a noise study, examining the robustness of our inverse model in presence of 0.1% and 1% additive noise. In all cases, the blood flow reconstruction error was <2  %   for most of the vasculature, except for the peripheral vasculature which suffered from insufficient photon sampling. Descending vasculature and deeper structures showed slightly higher sensitivity to noise compared with vasculature with a horizontal orientation at the more superficial layers. Our results show high-resolution reconstruction of the blood flow map in tissue down to 500  μm and beyond.

CONCLUSIONS

We have demonstrated a high-resolution computational imaging technique for visualizing blood flow map in complex tissue over a large FOV. Once a high-resolution structural image is captured, our reconstruction algorithm only requires a few LSCI images captured through a camera to reconstruct the blood flow map computationally at a high resolution. We note that the combination of high temporal and spatial resolution of our reconstruction algorithm makes the solution well-suited for applications involving fast monitoring of flow dynamics over a large FOV, such as in functional neural imaging.

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.

Efficient inversion strategies for estimating optical properties with Monte Carlo radiative transport models

SIGNIFICANCE

Indirect imaging problems in biomedical optics generally require repeated evaluation of forward models of radiative transport, for which Monte Carlo is accurate yet computationally costly. We develop an approach to reduce this bottleneck, which has significant implications for quantitative tomographic imaging in a variety of medical and industrial applications.

AIM

Our aim is to enable computationally efficient image reconstruction in (hybrid) diffuse optical modalities using stochastic forward models.

APPROACH

Using Monte Carlo, we compute a fully stochastic gradient of an objective function for a given imaging problem. Leveraging techniques from the machine learning community, we then adaptively control the accuracy of this gradient throughout the iterative inversion scheme to substantially reduce computational resources at each step.

RESULTS

For example problems of quantitative photoacoustic tomography and ultrasound-modulated optical tomography, we demonstrate that solutions are attainable using a total computational expense that is comparable to (or less than) that which is required for a single high-accuracy forward run of the same Monte Carlo model.

CONCLUSIONS

This approach demonstrates significant computational savings when approaching the full nonlinear inverse problem of optical property estimation using stochastic methods.