Abstract:Near-infrared spectroscopy (NIRS) including diffuse optical tomography is an imaging modality which makes use of diffuse light propagation in random media. When optical properties of a random medium are investigated from boundary measurements of reflected or transmitted light, iterative inversion schemes such as the Levenberg–Marquardt algorithm are known to fail when initial guesses are not close enough to the true value of the coefficient to be reconstructed. In this paper, we investigate how this weakness of iterative schemes is overcome using Markov chain Monte Carlo. Using time-resolved measurements performed against a polyurethane-based phantom, we present a case that the Levenberg–Marquardt algorithm fails to work but the proposed hybrid method works well. Then, with a toy model of diffuse optical tomography we illustrate that the Levenberg–Marquardt method fails when it is trapped by a local minimum but the hybrid method can escape from local minima by using the Metropolis–Hastings Markov chain Monte Carlo algorithm until it reaches the valley of the global minimum. The proposed hybrid scheme can be applied to different inverse problems in NIRS which are solved iteratively. We find that for both numerical and phantom experiments, optical properties such as the absorption and reduced scattering coefficients can be retrieved without being trapped by a local minimum when Monte Carlo simulation is run only about 100 steps before switching to an iterative method. The hybrid method is compared with simulated annealing. Although the Metropolis–Hastings MCMC arrives at the steady state at about 10,000 Monte Carlo steps, in the hybrid method the Monte Carlo simulation can be stopped way before the burn-in time.

The inverse Rytov series for diffuse optical tomography

Diffuse optical tomography in time domain with the inverse Rytov series

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

Shape-constrained Reconstruction in Diffuse Optical Tomography by Simulated Annealing

An Optical Tomography-Related Time-Domain Inverse Problern for the Diffusion Equation

Acoustical tomography based on the second-order Rytov transform perturbation approximation

The Forward and Inverse Models in Time-Resolved Optical Tomography Imaging and Their Finite-Element Method Solutions

Acoustical diffraction tomography in a finite form based on the Rytov transform.

Numerical modelling and image reconstruction in diffuse optical tomography

Reconstruction Method for Optical Tomography Based on the Linearized Bregman Iteration with Sparse Regularization.

Radiative Transport and Optical Tomography with Large Datasets

Simultaneous uniqueness and numerical inversion for an inverse problem in the time-domain diffuse optical tomography with fluorescence

Performance evaluation of typical approximation algorithms for nonconvex ℓp-minimization in diffuse optical tomography

Wave equation traveltime tomography using Rytov approximation

Image reconstruction algorithm for steady-state diffuse optical tomography with structural priori information

Optical Image Reconstruction of Non-Scattering and Low Scattering Heterogeneities in Turbid Media Based on the Diffusion Approximation Model

An Efficient and Robust Reconstruction Method for Optical Tomography with the Time-Domain Radiative Transfer Equation

A Hybrid Inversion Scheme Combining Markov Chain Monte Carlo and Iterative Methods for Determining Optical Properties of Random Media

Improved Optical Tomography Based on Hybrid Frequency-Domain and Time-Domain Radiative Transfer Model

Inversion formulae for ray transforms in vector and tensor tomography

Compressive diffuse optical tomography: noniterative exact reconstruction using joint sparsity