Photoacoustic tomography with a model-based approach involving realistic detector properties

Abstract

A computational and experimental study is conducted to examine how directivity associated with a finite aperture sensor affects photoacoustic tomography (PAT) image reconstruction. Acoustic signals for the simulation work were computed using a discrete particle approach from three numerical phantoms including a vasculature. The theoretical framework and a Monte Carlo approach for construction of a tissue configuration are discussed in detail. While simulating forward data, the directivity of the sensor was taken into account. The image reconstruction was accomplished using system matrix based methods like l2 norm Tikhonov regularization, l1 norm regularization and total variation (TV) minimization. Accordingly, two different system matrices were constructed- (i) assuming transducer as a point detector (PD) and (ii) retaining properties of a finite detector with directivity (FDWD). Image reconstruction was also performed utilizing experimentally measured PA signals. Both the computational and experimental results demonstrate that blur-free PAT imaging can be achieved with the FDWD method. Additionally, TV minimization provides marginally better image reconstruction compared to the other schemes.