Quantitative Reconstruction of Absorption Coefficients for Photoacoustic Tomography

Photoacoustic (PA) tomography (PAT) is a cutting-edge imaging modality for visualizing the internal structure and light-absorption distribution in tissue. However, reconstruction of the absorption distribution has been limited by nonuniform light fluence. This paper introduces a novel method for quantitative reconstruction of the distribution of optical absorption coefficients in tissue. In this method, we implement an iterative algorithm for recovering absorption coefficients from optical absorbed energy maps based on a 3D Monte Carlo simulation of light transport and integrated with fluence compensation to obtain the initialization parameters. In the iteration algorithm, we calculate the deviation between the detected and the computed absorbed energy distribution at each iteration. By minimizing the deviation in the absorbed energy, the recovered values converge to the true absorption distribution. The results of numerical simulation and phantom experiment theoretically and experimentally demonstrate that the proposed method performs an accurately quantitative estimate of the distribution of optical absorption coefficients. This work expects to provide accurate quantitative information for absorbers within tissues or organs, and thereby broaden the clinical applications of PAT.