Efficient optical parameter mapping based on time-domain radiative transfer equation combined with parallel programming.

A two-dimensional optical parameter mapping based on the time-domain radiative transfer equation (TD-RTE) is studied in this work. The finite element method with structured and unstructured grids is employed to solve TD-RTE and OpenMP parallel technology is employed to improve the computing efficiency. The sequential quadratic programming algorithm is used as a powerful optimization method to reconstruct absorption and scattering parameter fields and the maximum a posteriori estimation is employed by introducing the regularization term into the objective function to improve the ill-posed inverse problem. In addition, the effects of measurement errors on reconstruction accuracy are investigated thoroughly. All the simulation results demonstrate that the reconstructed scheme we developed is accurate and efficient in optical parameter mapping based on TD-RTE.