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.
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 |
PACS-Classification: | 80.00.00 INTERDISCIPLINARY PHYSICS AND RELATED AREAS OF SCIENCE AND TECHNOLOGY |
Date of first Publication: | 2017/03/24 |
Series (Serial Number): | ZIB-Report (17-04) |
ISSN: | 1438-0064 |
Published in: | Appeared in: Proc. SPIE 10064, Photons Plus Ultrasound: Imaging and Sensing 2017, 100645J (March 23, 2017) |
DOI: | https://doi.org/10.1117/12.2251945 |