Diffuse optical tomography by simulated annealing via a spin Hamiltonian

Diffuse optical tomography (DOT) is an imaging modality that uses near-infrared light. Although iterative numerical schemes are commonly used for its inverse problem, correct solutions are not obtained unless good initial guesses are chosen. We propose a numerical scheme of DOT, which works even when good initial guesses of optical parameters are not available. We use simulated annealing (SA), which is a method of the Markov-chain Monte Carlo. To implement SA for DOT, a spin Hamiltonian is introduced in the cost function, and the Metropolis algorithm or single-component Metropolis-Hastings algorithm is used. By numerical experiments, it is shown that an initial random spin configuration is brought to a converged configuration by SA, and targets in the medium are reconstructed. The proposed numerical method solves the inverse problem for DOT by finding the ground state of a spin Hamiltonian with SA. (c) 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

Automated summary

Diffuse optical tomography uses near-infrared light and an iterative numerical scheme, with simulated annealing proposed as a method to find solutions even without good initial guesses. The proposed numerical method successfully reconstructs targets in the medium by finding the ground state of a spin Hamiltonian using simulated annealing.