Three dimensional photoacoustic tomography in Bayesian framework
Self archived versionpublished version
MetadataShow full item record
CitationTick, J. Pulkkinen, A. Lucka, F. Ellwood, R. Cox, BT. Kaipio, JP. Arridge, SR. Tarvainen, T. (2018). Three dimensional photoacoustic tomography in Bayesian framework. JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA, 144 (4) , 2061. 10.1121/1.5057109.
The image reconstruction problem (or inverse problem) in photoacoustic tomography is to resolve the initial pressure distribution from detected ultrasound waves generated within an object due to an illumination by a short light pulse. Recently, a Bayesian approach to photoacoustic image reconstruction with uncertainty quantification was proposed and studied with two dimensional numerical simulations. In this paper, the approach is extended to three spatial dimensions and, in addition to numerical simulations, experimental data are considered. The solution of the inverse problem is obtained by computing point estimates, i.e., maximum a posteriori estimate and posterior covariance. These are computed iteratively in a matrix-free form using a biconjugate gradient stabilized method utilizing the adjoint of the acoustic forward operator. The results show that the Bayesian approach can produce accurate estimates of the initial pressure distribution in realistic measurement geometries and that the reliability of these estimates can be assessed.