Diffusive optical tomography in the Bayesian framework

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.

[1]  Andrew M. Stuart,et al.  Inverse problems: A Bayesian perspective , 2010, Acta Numerica.

[2]  Georg Stadler,et al.  Extreme-scale UQ for Bayesian inverse problems governed by PDEs , 2012, 2012 International Conference for High Performance Computing, Networking, Storage and Analysis.

[3]  G. G. Stokes "J." , 1890, The New Yale Book of Quotations.

[4]  J. Sylvester,et al.  A global uniqueness theorem for an inverse boundary value problem , 1987 .

[5]  Erkki Somersalo,et al.  Linear inverse problems for generalised random variables , 1989 .

[6]  Alexander D. Klose,et al.  Reprint of: Optical tomography using the time-independent equation of radiative transfer --- Part 1: Forward model , 2010 .

[7]  E. Somersalo,et al.  Statistical and computational inverse problems , 2004 .

[8]  Mourad Choulli,et al.  INVERSE SCATTERING AND INVERSE BOUNDARY VALUE PROBLEMS FOR THE LINEAR BOLTZMANN EQUATION , 1996 .

[9]  A. Stuart,et al.  The Bayesian Approach to Inverse Problems , 2013, 1302.6989.

[10]  R. Alcouffe,et al.  Comparison of finite-difference transport and diffusion calculations for photon migration in homogeneous and heterogeneous tissues. , 1998, Physics in medicine and biology.

[11]  G. Roberts,et al.  MCMC Methods for Functions: ModifyingOld Algorithms to Make Them Faster , 2012, 1202.0709.

[12]  N. Holmer,et al.  Electrical Impedance Tomography , 1991 .

[13]  J. Zabczyk,et al.  Stochastic Equations in Infinite Dimensions , 2008 .

[14]  Guillaume Bal,et al.  Inverse transport with isotropic sources and angularly averaged measurements , 2008 .

[15]  Guillaume Bal,et al.  Time-dependent angularly averaged inverse transport , 2009, 0902.3432.

[16]  C. Bardos,et al.  DIFFUSION APPROXIMATION AND COMPUTATION OF THE CRITICAL SIZE , 1984 .

[17]  Y. Marzouk,et al.  A stochastic collocation approach to Bayesian inference in inverse problems , 2009 .

[18]  P. Bassanini,et al.  Elliptic Partial Differential Equations of Second Order , 1997 .

[19]  Qin Li,et al.  Validity and Regularization of Classical Half-Space Equations , 2016, 1606.01311.

[20]  Kui Ren,et al.  Recent Developments in Numerical Techniques for Transport-Based Medical Imaging Methods , 2010 .

[21]  Lei Wu,et al.  Geometric Correction for Diffusive Expansion of Steady Neutron Transport Equation , 2014, 1404.2583.

[22]  Gideon Simpson,et al.  Algorithms for Kullback-Leibler Approximation of Probability Measures in Infinite Dimensions , 2014, SIAM J. Sci. Comput..

[23]  S. Arridge Optical tomography in medical imaging , 1999 .

[24]  Angela A. Eick,et al.  Mechanisms of light scattering from biological cells relevant to noninvasive optical-tissue diagnostics. , 1998, Applied optics.

[25]  Gunther Uhlmann,et al.  Inverse Problems for the Stationary Transport Equation in the Diffusion Scaling , 2018, SIAM J. Appl. Math..

[26]  Giovanni Alessandrini,et al.  Stable determination of conductivity by boundary measurements , 1988 .

[27]  Li Wang,et al.  Stability of stationary inverse transport equation in diffusion scaling , 2017, 1703.00097.

[28]  Jenn-Nan Wang Stability estimates of an inverse problem for the stationary transport equation , 1999 .

[29]  Andrew M. Stuart,et al.  Approximation of Bayesian Inverse Problems for PDEs , 2009, SIAM J. Numer. Anal..

[30]  Guillaume Bal,et al.  Inverse transport theory and applications , 2009 .