A novel approach for the time-domain fluorescence imaging of a semi-infinite turbid medium: Monte Carlo evaluation

We investigated the feasibility of a two-step scheme for reconstruction of a fluorophore target embedded in a semi-infinite medium. In this scheme, we neglected the presence of the fluorophore target for the excitation light and used an analytical solution of the time-dependent radiative transfer equation (RTE) for the excitation light in a homogeneous semi-infinite media instead of solving the RTE numerically in the forward calculation. In the first step of this reconstruction scheme, we implemented a pixel-based reconstruction using the Landweber method with adjoint fields. The second step uses this result as an initial guess for solving the shape and contrast value reconstruction problem using the level set method. Numerical experiments using Monte Carlo data measurements, show that the proposed scheme provides reconstructions of shape, location and contrast value of the target with rather good accuracy. The computation times of the solution of the forward problem and the whole reconstruction process were reduced by about forty and fifteen percent, respectively.

[1]  Rolf Rannacher,et al.  Numerical methods in multidimensional radiative transfer , 2009 .

[2]  Sung-Ho Han,et al.  Estimating the depth and lifetime of a fluorescent inclusion in a turbid medium using a simple time-domain optical method. , 2008, Optics letters.

[3]  Guillaume Bal,et al.  Inverse transport theory and applications , 2009 .

[4]  Alessandro Torricelli,et al.  Heuristic Green's function of the time dependent radiative transfer equation for a semi-infinite medium. , 2007, Optics express.

[5]  P. Yalavarthy,et al.  Analytical solutions for diffuse fluorescence spectroscopy/imaging in biological tissues. Part I: zero and extrapolated boundary conditions. , 2013, Journal of the Optical Society of America. A, Optics, image science, and vision.

[6]  Bertrand Zavidovique,et al.  From Data to Images: a Shape Based Approach for Fluorescence Tomography , 2011 .

[7]  Sung-Ho Han,et al.  Analysis of the fluorescence temporal point-spread function in a turbid medium and its application to optical imaging. , 2008, Journal of biomedical optics.

[8]  Otmar Scherzer,et al.  KACZMARZ METHODS FOR REGULARIZING NONLINEAR ILL-POSED EQUATIONS I: CONVERGENCE ANALYSIS , 2007 .

[9]  Alwin Kienle,et al.  Hybrid Green’s function of the time-dependent radiative transfer equation for anisotropically scattering semi-infinite media , 2013, Journal of biomedical optics.

[10]  H. Engl,et al.  Regularization of Inverse Problems , 1996 .

[11]  O. Dorn,et al.  Sparsity and level set regularization for diffuse optical tomography using a transport model in 2D , 2016 .

[12]  S R Arridge,et al.  Reconstructing absorption and diffusion shape profiles in optical tomography by a level set technique. , 2006, Optics letters.

[13]  Simon R. Arridge,et al.  Parameter and structure reconstruction in optical tomography , 2008 .

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

[15]  Yukio Yamada,et al.  Diffuse optical tomography: Present status and its future , 2014 .

[16]  Yong Wang,et al.  Simple time-domain optical method for estimating the depth and concentration of a fluorescent inclusion in a turbid medium. , 2004, Optics letters.

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

[18]  Oliver Dorn,et al.  A transport-backtransport method for optical tomography , 1998 .

[19]  Dominique Lesselier,et al.  Level Set Methods for Structural Inversion and Image Reconstruction , 2015, Handbook of Mathematical Methods in Imaging.

[20]  Simon R. Arridge,et al.  Fluorescence lifetime imaging in turbid media , 2014 .

[21]  Oliver Dorn,et al.  Shape reconstruction in scattering media with voids using a transport model and level sets , 2002 .

[22]  H Dehghani,et al.  A three-dimensional finite element model and image reconstruction algorithm for time-domain fluorescence imaging in highly scattering media , 2011, Physics in medicine and biology.

[23]  Barbara Kaltenbacher,et al.  Iterative Regularization Methods for Nonlinear Ill-Posed Problems , 2008, Radon Series on Computational and Applied Mathematics.

[24]  Goro Nishimura,et al.  A new scheme of the time-domain fluorescence tomography for a semi-infinite turbid medium , 2017 .

[25]  Fabrizio Martelli,et al.  Light Propagation Through Biological Tissue and Other Diffusive Media: Theory, Solutions, and Software , 2009 .

[26]  M. Patterson,et al.  Improved solutions of the steady-state and the time-resolved diffusion equations for reflectance from a semi-infinite turbid medium. , 1997, Journal of the Optical Society of America. A, Optics, image science, and vision.

[27]  Zhixiong Guo,et al.  FLUORESCENCE IMAGE RECONSTRUCTION FOR OPTICAL TOMOGRAPHY BASED ON TRANSIENT RADIATION TRANSFER EQUATION , 2003 .

[28]  D. Boas,et al.  Experimental images of heterogeneous turbid media by frequency-domain diffusing-photon tomography. , 1995, Optics letters.