Monte-Carlo-based inversion scheme for 3D quantitative photoacoustic tomography

Please always quote using this URN: urn:nbn:de:0297-zib-62318
  • The goal of quantitative photoacoustic tomography (qPAT) is to recover maps of the chromophore distributions from multiwavelength images of the initial pressure. Model-based inversions that incorporate the physical processes underlying the photoacoustic (PA) signal generation represent a promising approach. Monte-Carlo models of the light transport are computationally expensive, but provide accurate fluence distributions predictions, especially in the ballistic and quasi-ballistic regimes. Here, we focus on the inverse problem of 3D qPAT of blood oxygenation and investigate the application of the Monte-Carlo method in a model-based inversion scheme. A forward model of the light transport based on the MCX simulator and acoustic propagation modeled by the k-Wave toolbox was used to generate a PA image data set acquired in a tissue phantom over a planar detection geometry. The combination of the optical and acoustic models is shown to account for limited-view artifacts. In addition, the errors in the fluence due to, for example, partial volume artifacts and absorbers immediately adjacent to the region of interest are investigated. To accomplish large-scale inversions in 3D, the number of degrees of freedom is reduced by applying image segmentation to the initial pressure distribution to extract a limited number of regions with homogeneous optical parameters. The absorber concentration in the tissue phantom was estimated using a coordinate descent parameter search based on the comparison between measured and modeled PA spectra. The estimated relative concentrations using this approach lie within 5 % compared to the known concentrations. Finally, we discuss the feasibility of this approach to recover the blood oxygenation from experimental data.

Download full text files

Export metadata

  • Export Bibtex
Author:Bernhard Kaplan, Jan Laufer, Steffen Prohaska, Jens Buchmann
Document Type:ZIB-Report
Tag:Monte Carlo methods for light transport; boundary conditions; chromophore concentration; coordinate search; model-based inversion; oxygen saturation; photoacoustic imaging; quantitative photoacoustic tomography
MSC-Classification:68-XX COMPUTER SCIENCE (For papers involving machine computations and programs in a specific mathematical area, see Section -04 in that area)
CCS-Classification:J. Computer Applications
Date of first Publication:2017/03/24
Series (Serial Number):ZIB-Report (17-04)
Published in:Appeared in: Proc. SPIE 10064, Photons Plus Ultrasound: Imaging and Sensing 2017, 100645J (March 23, 2017)

$Rev: 13581 $