Application of SQP algorithm for fluorescence tomography with the time-domain equation of radiative transfer

Abstract A reconstruction scheme for the fluorescence tomography is investigated based on the time-domain radiative transfer equation (TD-RTE). Two coupled TD-RTEs, which can provide considerable measurement data, are used as the forward model and solved by the discrete ordinate method. The sequential quadratic programming (SQP) is employed to build the reconstruction scheme for solving the inverse problem. The gradient of objective function is calculated efficiently by the adjoint equation technique. Considering the ill-posed nature of the inverse problem, the regularization term based on the generalized Gaussian Markov random field (GGMRF) model is adopted to enhance the reconstructed image. Influence of the initial guess, contrast, noisy data, and shape of the fluorescent target are analyzed. Simulated results show that the proposed algorithm performs efficiently and accurately on reconstructing the distribution of the fluorescence yield.

[1]  Zhixiong Guo,et al.  Fast 3-d optical imaging with transient fluorescence signals. , 2004, Optics express.

[2]  Alexander D Klose,et al.  Fluorescence tomography with simulated data based on the equation of radiative transfer. , 2003, Optics letters.

[3]  Joan Boulanger,et al.  Numerical developments for short-pulsed Near Infra-Red laser spectroscopy. Part II: inverse treatment , 2005 .

[4]  A. Klose,et al.  Optical tomography using the time-independent equation of radiative transfer-Part 1: Forward model , 2002 .

[5]  Vasilis Ntziachristos,et al.  The inverse source problem based on the radiative transfer equation in optical molecular imaging , 2005 .

[6]  Charles A. Bouman,et al.  Optical diffusion tomography by iterative- coordinate-descent optimization in a Bayesian framework , 1999 .

[7]  Feng Gao,et al.  A linear, featured-data scheme for image reconstruction in time-domain fluorescence molecular tomography. , 2006, Optics express.

[8]  Marco Ferrari,et al.  A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application , 2012, NeuroImage.

[9]  Bin Liu,et al.  Light-field-camera imaging simulation of participatory media using Monte Carlo method , 2016 .

[10]  Huijuan Zhao,et al.  Three-dimensional scheme for time-domain fluorescence molecular tomography based on Laplace transforms with noise-robust factors. , 2008, Optics express.

[11]  Gregory Boverman,et al.  Time resolved fluorescence tomography of turbid media based on lifetime contrast. , 2006, Optics express.

[12]  Joan Boulanger,et al.  An overview on recent radiation transport algorithm development for optical tomography imaging , 2008 .

[13]  David A Boas,et al.  Comparison of frequency-domain and time-domain fluorescence lifetime tomography. , 2008, Optics letters.

[14]  S. Kumar,et al.  Discrete-ordinates solution of short-pulsed laser transport in two-dimensional turbid media. , 2001, Applied optics.

[15]  F Lesage,et al.  Time Domain Fluorescent Diffuse Optical Tomography: analytical expressions. , 2005, Optics express.

[16]  Davide Contini,et al.  Time domain functional NIRS imaging for human brain mapping , 2014, NeuroImage.

[17]  H. Tan,et al.  Multi-start iterative reconstruction of the radiative parameter distributions in participating media based on the transient radiative transfer equation , 2015 .

[18]  L. C. Henyey,et al.  Diffuse radiation in the Galaxy , 1940 .

[19]  John McGhee,et al.  Radiative transport-based frequency-domain fluorescence tomography , 2008, Physics in medicine and biology.

[20]  Masao Fukushima,et al.  A successive quadratic programming algorithm with global and superlinear convergence properties , 1986, Math. Program..

[21]  Andrew K. Dunn,et al.  A Time Domain Fluorescence Tomography System for Small Animal Imaging , 2008, IEEE Transactions on Medical Imaging.

[22]  R. Cubeddu,et al.  Time-resolved fluorescence imaging in biology and medicine , 2002 .

[23]  M. J. D. Powell,et al.  THE CONVERGENCE OF VARIABLE METRIC METHODS FOR NONLINEARLY CONSTRAINED OPTIMIZATION CALCULATIONS , 1978 .

[24]  Hong Qi,et al.  Application of the sequential quadratic programming algorithm for reconstructing the distribution of optical parameters based on the time-domain radiative transfer equation. , 2016, Optics express.

[25]  K. Schittkowski.  The nonlinear programming method of Wilson, Han, and Powell with an augmented Lagrangian type line search function , 1982 .

[26]  Fatmir Asllanaj,et al.  Sensitivity analysis to optical properties of biological tissues subjected to a short-pulsed laser using the time-dependent radiative transfer equation , 2014 .

[27]  J. Laible,et al.  Coupled complex adjoint sensitivities for frequency-domain fluorescence tomography: theory and vectorized implementation , 2003 .

[28]  Fei Liu,et al.  Reconstruction for limited-projection fluorescence molecular tomography based on projected restarted conjugate gradient normal residual. , 2011, Optics letters.

[29]  Hyun Keol Kim,et al.  PDE-Constrained Fluorescence Tomography With the Frequency-Domain Equation of Radiative Transfer , 2010 .

[30]  Brian W. Pogue,et al.  Mathematical model for time-resolved and frequency-domain fluorescence spectroscopy in biological tissues. , 1994, Applied optics.

[31]  Chuanlong Xu,et al.  An inverse method for flue gas shielded metal surface temperature measurement based on infrared radiation , 2016 .