Abstract:Bioluminescence tomography(BLT) is an instrumental molecular imaging modality designed for the 3D location and quantification of bioluminescent sources distribution in vivo. In our context, the diffusion approximation(DA) to radiative transfer equation(RTE) is utilized to model the forward process of light propagation. Mathematically, the solution uniqueness does not hold for DA-based BLT which is an inverse source problem of partial differential equations and hence is highly ill-posed. In the current work, we concentrate on a general regularization framework for BLT with Bregman distance as data fidelity and total variation(TV) as regularization. Two specializations of the Bregman distance, the least squares(LS) distance and Kullback-Leibler(KL) divergence, which correspond to the Gaussian and Poisson environments respectively, are demonstrated and the resulting regularization problems are denoted as LS+TV and KL+TV. Based on the constrained Landweber(CL) scheme and expectation maximization(EM) algorithm for BLT, iterative algorithms for the LS+ TV and KL+TV problems in the context of BLT are developed, which are denoted as CL-TV and EM-TV respectively. They are both essentially gradient-based algorithms alternatingly performing the standard CL or EM iteration step and the TV correction step which requires the solution of a weighted ROF model. Chambolle's duality-based approach is adapted and extended to solving the weighted ROF subproblem. Numerical experiments for a 3D heterogeneous mouse phantom are carried out and preliminary results are reported to verify and evaluate the proposed algorithms. It is found that for piecewise-constant sources both CL-TV and EM-TV outperform the conventional CL and EM algorithms for BLT.

Local Convergence of an Em-Like Image Reconstruction Method for Diffuse Optical Tomography

An Em-Like Optimization Scheme For Diffuse Optical Tomography

Meshfree strategy for diffuse optical tomography image reconstruction

Iterative Total Least-Squares Image Reconstruction Algorithm for Optical Tomography by the Conjugate Gradient Method

Applying Weierstrass approximation theorem for trigonometric polynomials to limited-view tomographic image reconstruction

Experimental Demonstration of an Analytic Method for Image Reconstruction in Optical Tomography with Large Data Sets

Analytic Image Reconstruction in Local Phase-Contrast Tomography

Spatial Location Weighted Optimization Scheme for DC Optical Tomography.

Investigation of algorithms for image reconstruction in local tomography

Image Reconstruction Based on Sequential Monte Carlo Principle for EMT

Image reconstruction based on nonlinear diffusion model for limited-angle computed tomography

A Wavelet-Based Multiresolution Regularized Least Squares Reconstruction Approach for Optical Tomography

Deep Learning Image Reconstruction Simulation for Electromagnetic Tomography

An Improved Finite-Element Based Reconstruction Algorithm for Fluorescence Tomography

A Reconstruction Algorithm for Fluorescence Molecular Tomography Based on the Total Least-Squares Rule and the Normalized Conjugate Gradient Method

Iterative Reconstruction for Bioluminescence Tomography with Total Variation Regularization

A higher order diffusion model for three-dimensional photon migration and image reconstruction in optical tomography.

Image Reconstruction Based on Sparsity and Sensitive Field Optimization for EMT

Image Reconstruction Algorithms of Computed Tomography

A finite-element-based reconstruction method for 3D fluorescence tomography

Image Reconstruction for Bioluminescence Tomography from Partial Measurement