Statistical Inversion Using Sparsity and Total Variation Prior And Monte Carlo Sampling Method For Diffuse Optical Tomography

Abstract: In this paper, we formulate the reconstruction problem in diffuse optical tomography (DOT) in a statistical setting for determining the optical parameters, scattering and absorption, from boundary photon density measurements. A special kind of adaptive Metropolis algorithm for the reconstruction procedure using sparsity and total variation prior is presented. Finally, a simulation study of this technique with different regularization functions and its comparison to the deterministic Iteratively Regularized Gauss Newton method shows the effectiveness and stability of the method.