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

Bernhard A. Kaplan; Jens Buchmann; Steffen Prohaska; Jan Laufer

doi:10.1117/12.2251945

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

Bernhard A. Kaplan, Jens Buchmann, Steffen Prohaska, Jan Laufer

Proceedings Volume 10064, Photons Plus Ultrasound: Imaging and Sensing 2017; 100645J (2017) https://doi.org/10.1117/12.2251945
Event: SPIE BiOS, 2017, San Francisco, California, United States

Abstract

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.

Citation Download Citation

Bernhard A. Kaplan, Jens Buchmann, Steffen Prohaska, and Jan Laufer "Monte-Carlo-based inversion scheme for 3D quantitative photoacoustic tomography", Proc. SPIE 10064, Photons Plus Ultrasound: Imaging and Sensing 2017, 100645J (23 March 2017); https://doi.org/10.1117/12.2251945

ACCESS THE FULL ARTICLE

INSTITUTIONAL
Select your institution to access the SPIE Digital Library.

SELECT YOUR INSTITUTION

PERSONAL
Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.

PERSONAL SIGN IN

No SPIE Account? Create one

PURCHASE THIS CONTENT

SUBSCRIBE TO DIGITAL LIBRARY

50 downloads per 1-year subscription

Members: $195

Non-members: $335 ADD TO CART

25 downloads per 1 - year subscription

Members: $145

Non-members: $250 ADD TO CART

PURCHASE SINGLE ARTICLE

Includes PDF, HTML & Video, when available