Dynamic estimation of the modeling error statistics in Diffuse Optical Tomography

Diffuse Optical Tomography (DOT) is a non-invasive imaging technique that suffers from a typical large-scale and ill-posed inverse problem with low spatial resolution. In DOT, the inverse problem is computationally intensive and decreasing the computation complexity and making it well-posed is the one of the most challenging research areas. More precisely, one of the well-known complexity reduction techniques is defined as applying modelling error originated from discretization of forward problem. Applying the discretization error in Bayesian inference has already been discussed; the method in which the likelihood is modified by an off-line prior density estimation. This paper implements a new method to enhance the modelling error approach using an iterative scheme to update statistical parameters of modelling discrepancy in DOT. The algorithm is very similar to Ensemble Kalman Filter. Moreover, the reconstruction process in the applied method is conducted by a small sample size rather than off-line method. Hence, the computation complexity is decreased and the algorithm converges in few iterations. The efficiency of the proposed method is illustrated by simulations.