Efficient inclusion of total variation type priors in quantitative photoacoustic tomography

Research output: Scientific - peer-reviewArticle

Standard

Efficient inclusion of total variation type priors in quantitative photoacoustic tomography. / Hannukainen, A.; Hyvönen, N.; Majander, H.; Tarvainen, T.

In: SIAM JOURNAL ON IMAGING SCIENCES, Vol. 9, No. 3, 16.08.2016, p. 1132-1153.

Research output: Scientific - peer-reviewArticle

Harvard

APA

Vancouver

Author

Hannukainen, A.; Hyvönen, N.; Majander, H.; Tarvainen, T. / Efficient inclusion of total variation type priors in quantitative photoacoustic tomography.

In: SIAM JOURNAL ON IMAGING SCIENCES, Vol. 9, No. 3, 16.08.2016, p. 1132-1153.

Research output: Scientific - peer-reviewArticle

Bibtex - Download

@article{376fd3ccfd7a41309e6b7b40a706e469,
title = "Efficient inclusion of total variation type priors in quantitative photoacoustic tomography",
keywords = "LSQR, Matrix-free implementation, Perona-Malik, Priorconditioning, Quantitative photoacoustic tomography, Total variation",
author = "A. Hannukainen and N. Hyvönen and H. Majander and T. Tarvainen",
year = "2016",
month = "8",
doi = "10.1137/15M1051737",
volume = "9",
pages = "1132--1153",
journal = "SIAM JOURNAL ON IMAGING SCIENCES",
issn = "1936-4954",
publisher = "Society for Industrial and Applied Mathematics Publications",
number = "3",

}

RIS - Download

TY - JOUR

T1 - Efficient inclusion of total variation type priors in quantitative photoacoustic tomography

AU - Hannukainen,A.

AU - Hyvönen,N.

AU - Majander,H.

AU - Tarvainen,T.

PY - 2016/8/16

Y1 - 2016/8/16

N2 - 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.

AB - 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.

KW - LSQR

KW - Matrix-free implementation

KW - Perona-Malik

KW - Priorconditioning

KW - Quantitative photoacoustic tomography

KW - Total variation

UR - http://www.scopus.com/inward/record.url?scp=84989338130&partnerID=8YFLogxK

U2 - 10.1137/15M1051737

DO - 10.1137/15M1051737

M3 - Article

VL - 9

SP - 1132

EP - 1153

JO - SIAM JOURNAL ON IMAGING SCIENCES

T2 - SIAM JOURNAL ON IMAGING SCIENCES

JF - SIAM JOURNAL ON IMAGING SCIENCES

SN - 1936-4954

IS - 3

ER -

ID: 9103271