A hybrid inversion scheme of Markov-chain Monte Carlo and iterative methods for diffuse optical tomography

Yu Jiang,Manabu Machida,Gen Nakamura
2018-01-01
Abstract:Diffuse optical tomography is formulated as inverse coefficient problems for the diffusion equation. Iterative inversion schemes such as the Levenberg-Marquardt algorithm are known to fail when initial guesses are not close to the true value of the coefficient to be reconstructed. In this paper, we investigate how this weakness of iterative schemes is overcome by the use of Monte Carlo. We present a toy model of diffuse optical tomography for which the Levenberg-Marquardt algorithm fails to work but the Metropolis-Hastings Markov chain Monte Carlo works. We show that our proposed hybrid scheme solves the inverse problem efficiently by preparing a good initial guess by Monte Carlo and then computing the reconstructed value with the Levenberg-Marquardt algorithm starting from the found initial guess.
What problem does this paper attempt to address?