Efficient inversion strategies for estimating optical properties with Monte Carlo radiative transport models

TL;DR

This paper introduces an adaptive Monte Carlo-based inversion method that significantly reduces computational costs in optical property estimation for biomedical imaging, enabling efficient quantitative tomography.

Contribution

It develops a stochastic gradient approach with adaptive accuracy control, improving efficiency in Monte Carlo radiative transport models for inverse problems.

Findings

Achieves comparable or lower computational cost than a single high-accuracy Monte Carlo run.

Enables practical solutions for complex optical property estimation problems.

Demonstrates effectiveness in quantitative photoacoustic and ultrasound-modulated optical tomography.

Abstract

Indirect imaging problems in biomedical optics generally require repeated evaluation of forward models of radiative transport, for which Monte Carlo is accurate yet computationally costly. We develop a novel approach to reduce this bottleneck which has significant implications for quantitative tomographic imaging in a variety of medical and industrial applications. Using Monte Carlo we compute a fully stochastic gradient of an objective function for a given imaging problem. Leveraging techniques from the machine learning community we then adaptively control the accuracy of this gradient throughout the iterative inversion scheme, in order to substantially reduce computational resources at each step. For example problems of Quantitative Photoacoustic Tomography and Ultrasound Modulated Optical Tomography, we demonstrate that solutions are attainable using a total computational expense that is comparable to (or less than) that which is required for a single high accuracy forward run of the same Monte Carlo model. This approach demonstrates significant computational savings when approaching the full non-linear inverse problem of optical property estimation using stochastic methods.