A three-dimensional finite elements approach for the coupled radiative transfer equation and diffusion approximation modeling in fluorescence imaging

During the last few years a quite large number of fluorescence molecular imaging applications have been reported in the literature, as one of the most challenging aspects in medical imaging is to “see” a tumor embedded into tissue, which is a turbid medium, by using fluorescent probes for tumor labeling. However, the forward solvers, required for the successful convergence of the inverse problem, are still lacking accuracy and time feasibility. Moreover, initialization of these solvers may be proven even more difficult than solving the inverse problem itself. This paper describes in depth a coupled radiative transfer equation and diffusion approximation model for solving the forward problem in fluorescence imaging. The theoretical confrontation of these solvers comprises the model deployment, its Galerkin finite elements approximation and the domain discretization scheme. Finally, a new optical properties mapping algorithm, based on super-ellipsoid models, is implemented, providing a fully automated simulation target construction within feasible time.