Determining Sources in the Bioluminescence Tomography Problem
Ming-Hui Ding, Rongfang Gong, Hongyu Liu, Catharine W. K. Lo
TL;DR
This paper addresses the bioluminescence tomography (BLT) inverse source problem under the diffusion approximation, aiming to reconstruct an internal light source from boundary measurements. It advances the theory by proving uniqueness for sources of the form $q=\varphi\chi_\omega$ when the domain $\omega$ is either a $C^2$-smooth region or a polyhedral/corona-shaped region, using a single boundary datum, and it supports these results with comprehensive numerical experiments. The authors develop a microlocal CGO-based analysis to handle corners and non-smooth geometries, and they implement a Levenberg-Marquardt inversion with Tikhonov regularization to numerically recover both the support and intensity of the sources. The work broadens the class of geometries for which BLT uniqueness can be guaranteed and demonstrates practical recoverability for both smooth and non-smooth source domains, with potential impact on in vivo imaging and related biomedical applications.
Abstract
In this paper, we revisit the bioluminescence tomography (BLT) problem, where one seeks to reconstruct bioluminescence signals (an internal light source) from external measurements of the Cauchy data. As one kind of optical imaging, the BLT has many merits such as high signal-to-noise ratio, non-destructivity and cost-effectiveness etc., and has potential applications such as cancer diagnosis, drug discovery and development as well as gene therapies and so on. In the literature, BLT is extensively studied based on diffusion approximation (DA) equation, where the distribution of peak sources is to be reconstructed and no solution uniqueness is guaranteed without adequate a priori information. Motivated by the solution uniqueness issue, several theoretical results are explored. The major contributions in this work that are new to the literature are two-fold: first, we show the theoretical uniqueness of the BLT problem where the light sources are in the shape of $C^2$ domains or polyhedral- or corona-shaped; second, we support our results with plenty of problem-orientated numerical experiments.
