Reconstruction-oriented multigrid finite element algorithm on bioluminescence tomography incorporating priori information

摘要

Bioluminescence tomography (BLT) is employed to reconstruct internal bioluminescent source to reveal the molecular and cellular information. However, BLT faces many challenges such as quantitative reconstruction. In this paper, a reconstruction-oriented multigrid finite element algorithm is proposed to fully reconstruct the source density. A source permissible region is utilized to raise the numerical stability. The proposed algorithm transfers diffuse equation into liner relationship between the inner source and boundary information on a mesh. A tolerant algorithm for linearly constrained optimization algorithm can solve the box-constrained quadratic optimization problem. By setting a threshold as the average reconstruction densities, we can refine every possible source element and form next mesh until the reconstructed densities are accepted. Numerical simulation of homogeneous and heterogeneous phantom demonstrates the feasibility and potential of the proposed algorithm.