This talk considers quantitative photoacoustic tomography under the prior assumption that the examined body is characterized by piecewise constant optical parameters. An iterative reconstruction algorithm that aims at finding the MAP estimate for these parameters is introduced. The method is based on iteratively combining a lagged diffusivity step and linearization of the forward model with priorconditioned LSQR. It is numerically demonstrated the algorithm is capable of producing reconstructions on dense, three-dimensional finite element meshes.

This presentation is part of Minisymposium “MS61 - Imaging with Light and Sound (3 parts)”
organized by: Felix Lucka (CWI & UCL) , Tanja Tarvainen (University of Eastern Finland) .

Authors:

Nuutti Hyvönen (Aalto University)

Keywords: