Bioluminescence tomography with Gaussian prior

Parameterizing the bioluminescent source globally in Gaussians provides several advantages over voxel representation in bioluminescence tomography. It is mathematically unique to recover Gaussians [Med. Phys. 31(8), 2289 (2004)] and practically sufficient to approximate various shapes by Gaussians in diffusive medium. The computational burden is significantly reduced since much fewer unknowns are required. Besides, there are physiological evidences that the source can be modeled by Gaussians. The simulations show that the proposed model and algorithm significantly improves accuracy and stability in the presence of Gaussian or non- Gaussian sources, noisy data or the optical background mismatch. It is also validated through in vivo experimental data.