title: Quantitative photoacoustic tomography using forward and adjoint Monte Carlo models of radiance.
abstract: Forward and adjoint Monte Carlo (MC) models of radiance are proposed for use in model-based quantitative photoacoustic tomography. A two-dimensional (2-D) radiance MC model using a harmonic angular basis is introduced and validated against analytic solutions for the radiance in heterogeneous media. A gradient-based optimization scheme is then used to recover 2-D absorption and scattering coefficients distributions from simulated photoacoustic measurements. It is shown that the functional gradients, which are a challenge to compute efficiently using MC models, can be calculated directly from the coefficients of the harmonic angular basis used in the forward and adjoint models. This work establishes a framework for transport-based quantitative photoacoustic tomography that can fully exploit emerging highly parallel computing architectures.
authors.affiliation: University College London, Department of Medical Physics and Biomedical Engineering, Malet Place, WC1E 6BT London, United Kingdom.
authors.forename: Roman
authors.initials: R
authors.lastname: Hochuli
authors.affiliation: University College London, Department of Medical Physics and Biomedical Engineering, Malet Place, WC1E 6BT London, United Kingdom.
authors.forename: Samuel
authors.initials: S
authors.lastname: Powell
authors.affiliation: University College London, Department of Computer Science, Malet Place, WC1E 6BT London, United Kingdom.
authors.forename: Simon
authors.initials: S
authors.lastname: Arridge
authors.affiliation: University College London, Department of Medical Physics and Biomedical Engineering, Malet Place, WC1E 6BT London, United Kingdom.
authors.forename: Ben
authors.initials: B
authors.lastname: Cox
mesh_terms: D000465:Algorithms; D008027:Light; D009010:Monte Carlo Method; D061088:Photoacoustic Techniques / Q000379:methods*; D012542:Scattering, Radiation; D014054:Tomography / Q000379:methods*