The photoacoustic image reconstruction problem (inverse problem) is to estimate an initial acoustic pressure distribution from measurements of ultrasound waves generated within an object due to optical excitation with a short light pulse. In this work, the recently suggested Bayesian approach to photoacoustic tomography is extended to three dimensions and an iterative matrix-free method for the solution of the problem is described. Image reconstruction is investigated with numerical simulations and experimental data. The use of different prior information and noise models in different sensor geometries, including a limited-view setup, is investigated. The results show that the Bayesian approach can produce accurate estimates of the initial pressure distribution even in a limited-view setup provided that prior information and the noise have been properly modelled.

, ,
Centrum Wiskunde & Informatica, Amsterdam (CWI), The Netherlands

Tick, J, Pulkkinen, A, Lucka, F, Ellwood, R, Cox, B.T, Arridge, S, & Tarvainen, T. (2018). Photoacoustic image reconstruction in Bayesian framework. In Progress in Biomedical Optics and Imaging - Proceedings of SPIE. doi:10.1117/12.2288163