TY - GEN A1 - Buchmann, Jens A1 - Kaplan, Bernhard A1 - Powell, Samuel A1 - Prohaska, Steffen A1 - Laufer, Jan T1 - 3D quantitative photoacoustic tomography using an adjoint radiance Monte Carlo model and gradient descent N2 - Quantitative photoacoustic tomography aims to recover maps of the local concentrations of tissue chromophores from multispectral images. While model-based inversion schemes are promising approaches, major challenges to their practical implementation include the unknown fluence distribution and the scale of the inverse problem. This paper describes an inversion scheme based on a radiance Monte Carlo model and an adjoint-assisted gradient optimization that incorporates fluence-dependent step sizes and adaptive moment estimation. The inversion is shown to recover absolute chromophore concentrations, blood oxygen saturation and the Grüneisen parameter from in silico 3D phantom images for different radiance approximations. The scattering coefficient was assumed to be homogeneous and known a priori. T3 - ZIB-Report - 19-16 KW - quantitative photoacoustic imaging KW - blood oxygen saturation KW - inverse problem KW - model-based inversion KW - Monte Carlo KW - spectral unmixing Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-72995 SN - 1438-0064 ER - TY - GEN A1 - Buchmann, Jens A1 - Kaplan, Bernhard A1 - Powell, Samuel A1 - Prohaska, Steffen A1 - Laufer, Jan T1 - Quantitative PA tomography of high resolution 3-D images: experimental validation in tissue phantoms N2 - Quantitative photoacoustic tomography aims recover the spatial distribution of absolute chromophore concentrations and their ratios from deep tissue, high-resolution images. In this study, a model-based inversion scheme based on a Monte-Carlo light transport model is experimentally validated on 3-D multispectral images of a tissue phantom acquired using an all-optical scanner with a planar detection geometry. A calibrated absorber allowed scaling of the measured data during the inversion, while an acoustic correction method was employed to compensate the effects of limited view detection. Chromophore- and fluence-dependent step sizes and Adam optimization were implemented to achieve rapid convergence. High resolution 3-D maps of absolute concentrations and their ratios were recovered with high accuracy. Potential applications of this method include quantitative functional and molecular photoacoustic tomography of deep tissue in preclinical and clinical studies. T3 - ZIB-Report - 19-60 KW - quantitative photoacoustic imaging KW - blood oxygen saturation KW - inverse problem KW - model-based inversion KW - Monte Carlo KW - spectral unmixing Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-75416 SN - 1438-0064 ER - TY - GEN A1 - Kaplan, Bernhard A1 - Laufer, Jan A1 - Prohaska, Steffen A1 - Buchmann, Jens T1 - Monte-Carlo-based inversion scheme for 3D quantitative photoacoustic tomography N2 - 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. T3 - ZIB-Report - 17-04 KW - quantitative photoacoustic tomography KW - model-based inversion KW - oxygen saturation KW - chromophore concentration KW - photoacoustic imaging KW - Monte Carlo methods for light transport KW - boundary conditions KW - coordinate search Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-62318 SN - 1438-0064 ER -