Photoacoustic Tomography with Time-Dependent Damping: Theoretical and a Convolutional Neural Network-Guided Numerical Inversion Procedure

Abstract

Abstract. In photoacoustic tomography (PAT), a hybrid imaging modality that is based on the acoustic detection of optical absorption from biological tissue exposed to a pulsed laser, a short pulse laser generates an initial pressure proportional to the absorbed optical energy, which then propagates acoustically and is measured on the boundary. To account for the significant signal distortion caused by acoustic attenuation in biological tissue, we model PAT in heterogeneous media using a damped wave equation featuring spatially varying sound speed and a time-dependent damping term. Under natural assumptions, we show that the initial pressure is uniquely determined by the boundary measurements using a harmonic extension of the boundary data with energy decay. For constant damping, an expansion in Dirichlet eigenfunctions of [Formula: see text] leads to an explicit series reconstruction formula for the initial pressure. Finally, we develop a gradient free numerical method based on the Pontryagin’s maximum principle to provide a robust and computationally viable approach to image reconstruction in attenuating PAT.