  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.
