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

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.

[1]  G. C. Pomraning,et al.  Linear Transport Theory , 1967 .

[2]  M. Schweiger,et al.  Gauss–Newton method for image reconstruction in diffuse optical tomography , 2005, Physics in medicine and biology.

[3]  M. Herranz,et al.  Optical Imaging in Breast Cancer Diagnosis: The Next Evolution , 2012, Journal of oncology.

[4]  Tom Goldstein,et al.  The Split Bregman Method for L1-Regularized Problems , 2009, SIAM J. Imaging Sci..

[5]  B. Pogue,et al.  Optical image reconstruction using frequency-domain data: simulations and experiments , 1996 .

[6]  Guillaume Bal,et al.  Transport- and diffusion-based optical tomography in small domains: a comparative study. , 2007, Applied optics.

[7]  Phaneendra K. Yalavarthy,et al.  Sparse Recovery Methods Hold Promise for Diffuse Optical Tomographic Image Reconstruction , 2014, IEEE Journal of Selected Topics in Quantum Electronics.

[8]  Ge Wang,et al.  A Theoretical Framework of X-Ray Dark-Field Tomography , 2011, SIAM J. Appl. Math..

[9]  Bo Han,et al.  Image Reconstruction for Diffuse Optical Tomography Based on Radiative Transfer Equation , 2015, Comput. Math. Methods Medicine.

[10]  Hua-bei Jiang.  Diffuse Optical Tomography , 2010 .

[11]  E. Lewis,et al.  Computational Methods of Neutron Transport , 1993 .

[12]  Elizabeth M C Hillman,et al.  Optical brain imaging in vivo: techniques and applications from animal to man. , 2007, Journal of biomedical optics.

[13]  R. Kress,et al.  Inverse Acoustic and Electromagnetic Scattering Theory , 1992 .

[14]  Huabei Jiang,et al.  Enhancing finite element-based photoacoustic tomography using total variation minimization , 2011 .

[15]  Tony F. Chan,et al.  High-Order Total Variation-Based Image Restoration , 2000, SIAM J. Sci. Comput..

[16]  Weimin Han,et al.  RTE-based bioluminescence tomography: A theoretical study , 2011 .

[17]  Hongming Zhang,et al.  Speckle noise reduction algorithm with total variation regularization in optical coherence tomography. , 2015, Optics express.

[18]  S Arridge,et al.  A gradient-based optimisation scheme foroptical tomography. , 1998, Optics express.

[19]  C. Clason,et al.  Total variation regularization for nonlinear fluorescence tomography with an augmented Lagrangian splitting approach. , 2010, Applied optics.

[20]  S R Arridge,et al.  Recent advances in diffuse optical imaging , 2005, Physics in medicine and biology.

[21]  Rachid Deriche,et al.  Regularizing Flows for Constrained Matrix-Valued Images , 2004 .

[22]  Shuang Zhang,et al.  Reconstruction Method for Optical Tomography Based on the Linearized Bregman Iteration with Sparse Regularization , 2015, Comput. Math. Methods Medicine.

[23]  M. Nikolova.  An Algorithm for Total Variation Minimization and Applications , 2004 .

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

[25]  Hamid Dehghani,et al.  Near infrared optical tomography using NIRFAST: Algorithm for numerical model and image reconstruction. , 2009, Communications in numerical methods in engineering.

[26]  Yiqiu Dong,et al.  An algorithm for total variation regularized photoacoustic imaging , 2015, Adv. Comput. Math..

[27]  P. Lions,et al.  Image recovery via total variation minimization and related problems , 1997 .

[28]  Ronald F. Gariepy,et al.  Measure Theory and Fine Properties of Functions, Revised Edition , 1865 .

[29]  Hongkai Zhao,et al.  Multilevel bioluminescence tomography based on radiative transfer equation part 2: total variation and l1 data fidelity. , 2010, Optics express.

[30]  Hongkai Zhao,et al.  A Fast-Forward Solver of Radiative Transfer Equation , 2009 .

[31]  Stanley Osher,et al.  Bregman methods in quantitative photoacoustic tomography , 2010 .

[32]  Weimin Han,et al.  Discrete-Ordinate Discontinuous Galerkin Methods for Solving the Radiative Transfer Equation , 2010, SIAM J. Sci. Comput..

[33]  Jie Tian,et al.  Total variation regularization for bioluminescence tomography with the split Bregman method. , 2012, Applied optics.

[34]  Weimin Han,et al.  A theoretical study for RTE-based parameter identification problems , 2013 .

[35]  TWO-WEEK Loan COpy,et al.  University of California , 1886, The American journal of dental science.

[36]  Jiaying Xiao,et al.  Detection of cervical cancer based on photoacoustic imaging-the in-vitro results. , 2015, Biomedical optics express.