Diffusive optical tomography in the Bayesian framework
This work provides a theoretical foundation for using diffusion approximations in Bayesian inverse problems for optical tomography, relevant to medical imaging.
The paper develops a Bayesian framework for diffusive optical tomography, proving that the inverse problem for the radiative transfer equation converges to that for the diffusion equation in the strong scattering regime, with convergence measured via Hellinger distance and Kullback-Leibler divergence.
Many naturally-occuring models in the sciences are well-approximated by simplified models, using multiscale techniques. In such settings it is natural to ask about the relationship between inverse problems defined by the original problem and by the multiscale approximation. We develop an approach to this problem and exemplify it in the context of optical tomographic imaging. Optical tomographic imaging is a technique for infering the properties of biological tissue via measurements of the incoming and outgoing light intensity; it may be used as a medical imaging methodology. Mathematically, light propagation is modeled by the radiative transfer equation (RTE), and optical tomography amounts to reconstructing the scattering and the absorption coefficients in the RTE from boundary measurements. We study this problem in the Bayesian framework, focussing on the strong scattering regime. In this regime the forward RTE is close to the diffusion equation (DE). We study the RTE in the asymptotic regime where the forward problem approaches the DE, and prove convergence of the inverse RTE to the inverse DE in both nonlinear and linear settings. Convergence is proved by studying the distance between the two posterior distributions using the Hellinger metric, and using Kullback-Leibler divergence.