Efficient inclusion of total variation type priors in quantitative photoacoustic tomography

A. Hannukainen, N. Hyvönen, H. Majander, T. Tarvainen

Research output: Contribution to journalArticleScientificpeer-review

5 Citations (Scopus)
197 Downloads (Pure)


Quantitative photoacoustic tomography is an emerging imaging technique aimed at estimating the distribution of optical parameters inside tissues from photoacoustic images, which are formed by combining optical information and ultrasonic propagation. This optical parameter estimation problem is ill-posed and needs to be approached within the framework of inverse problems. Photoacoustic images are three-dimensional and high-resolution. Furthermore, high-resolution reconstructions of the optical parameters are targeted. Therefore, in order to provide a practical method for quantitative photoacoustic tomography, the inversion algorithm needs to be able to perform successfully with problems of prominent size. In this work, an efficient approach for the inverse problem of quantitative photoacoustic tomography is proposed, assuming an edge-preferring prior for the optical parameters. The method is based on iteratively combining priorconditioned LSQR with a lagged diffusivity step and a linearization of the measurement model, with the needed multiplications by Jacobians performed in a matrix-free manner. The algorithm is tested with three-dimensional numerical simulations. The results show that the approach can be used to produce accurate and high quality estimates of absorption and diffusion in complex three-dimensional geometries with moderate computation time and cost.

Original languageEnglish
Pages (from-to)1132-1153
Number of pages22
Issue number3
Publication statusPublished - 16 Aug 2016
MoE publication typeA1 Journal article-refereed


  • LSQR
  • Matrix-free implementation
  • Perona-Malik
  • Priorconditioning
  • Quantitative photoacoustic tomography
  • Total variation

Fingerprint Dive into the research topics of 'Efficient inclusion of total variation type priors in quantitative photoacoustic tomography'. Together they form a unique fingerprint.

Cite this