An efficient and robust reconstruction method for optical tomography with the time-domain radiative transfer equation

Diffuse optical tomography is a novel molecular imaging technology for small animal studies. Most known reconstruction methods use the diffusion equation (DA) as forward model, although the validation of DA breaks down in certain situations. In this work, we use the radiative transfer equation as forward model which provides an accurate description of the light propagation within biological media and investigate the potential of sparsity constraints in solving the diffuse optical tomography inverse problem. The feasibility of the sparsity reconstruction approach is evaluated by boundary angular-averaged measurement data and internal angular-averaged measurement data. Simulation results demonstrate that in most of the test cases the reconstructions with sparsity regularization are both qualitatively and quantitatively more reliable than those with standardL2regularization. Results also show the competitive performance of the split Bregman algorithm for the DOT image reconstruction with sparsity regularization compared with other existingL1algorithms.

Download Full-text

Diffuse Optical Tomography Based on Radiative Transfer Equation with L∞ Data Fidelity and Total Variation Penalty Regularization

Journal of Medical Imaging and Health Informatics ◽

10.1166/jmihi.2017.2154 ◽

2017 ◽

Vol 7 (6) ◽

pp. 1212-1218 ◽

Cited By ~ 1

Author(s):

Bo Bi ◽

Jinping Tang

Keyword(s):

Radiative Transfer ◽

Total Variation ◽

Transfer Equation ◽

Radiative Transfer Equation ◽

Optical Tomography ◽

Diffuse Optical Tomography ◽

Data Fidelity

Download Full-text

Mixed Total Variation and L1 Regularization Method for Optical Tomography Based on Radiative Transfer Equation

Computational and Mathematical Methods in Medicine ◽

10.1155/2017/2953560 ◽

2017 ◽

Vol 2017 ◽

pp. 1-15 ◽

Cited By ~ 3

Author(s):

Jinping Tang ◽

Bo Han ◽

Weimin Han ◽

Bo Bi ◽

Li Li

Keyword(s):

Radiative Transfer ◽

Total Variation ◽

Transfer Equation ◽

Radiative Transfer Equation ◽

Imaging Modality ◽

Optical Tomography ◽

Identification Problem ◽

Split Bregman ◽

Computation Efficiency ◽

Reconstruction Methods

Optical tomography is an emerging and important molecular imaging modality. The aim of optical tomography is to reconstruct optical properties of human tissues. In this paper, we focus on reconstructing the absorption coefficient based on the radiative transfer equation (RTE). It is an ill-posed parameter identification problem. Regularization methods have been broadly applied to reconstruct the optical coefficients, such as the total variation (TV) regularization and the L1 regularization. In order to better reconstruct the piecewise constant and sparse coefficient distributions, TV and L1 norms are combined as the regularization. The forward problem is discretized with the discontinuous Galerkin method on the spatial space and the finite element method on the angular space. The minimization problem is solved by a Jacobian-based Levenberg-Marquardt type method which is equipped with a split Bregman algorithms for the L1 regularization. We use the adjoint method to compute the Jacobian matrix which dramatically improves the computation efficiency. By comparing with the other imaging reconstruction methods based on TV and L1 regularizations, the simulation results show the validity and efficiency of the proposed method.

Download Full-text

Propagation of Ultra-Short-Pulse Radiation in Participating Media: A Laguerre-Galerkin Solution

Volume 11: Micro and Nano Systems, Parts A and B ◽

10.1115/imece2007-41911 ◽

2007 ◽

Author(s):

Tuba Okutucu ◽

Yaman Yener

Keyword(s):

Radiative Transfer ◽

Transfer Equation ◽

Radiative Transfer Equation ◽

Transfer Problem ◽

Short Pulse ◽

Optical Tomography ◽

Time Derivative ◽

Product Analysis ◽

Participating Media ◽

Radiative Transfer Problem

Transient analysis of the radiative transfer problem in participating media has become essential due to the recent applications involving extremely small time scales. In classical radiation problems, the time derivative term in the radiative transfer equation has a negligible order of magnitude compared to the others. Lasers of pico- to femtosecond pulse durations are now being used to investigate the properties of scattering and absorbing media in such applications as, optical tomography, combustion product analysis, and remote sensing. For such applications, the time derivative in the radiative transfer equation can no longer be neglected. Numerous approaches such as, integral formulation, direct numerical approach, discrete ordinates method, Monte Carlo simulations, and Galerkin technique have been introduced for the solution of transient radiative transfer problems in participating media. In the present work, Laguerre-Galerkin solutions for both rectangular and Gaussian incident pulse profiles are presented.

Download Full-text