Quantitative Evaluation of Bioluminescence Imaging Based on Radiative Transfer Equation
KAN Xing-xing1, CHEN Chun-xiao1, GONG Rong-fang2
1. Department of Biomedical Engineering，Nanjing University of Aeronautics & Astronautics, Jiangsu Province Nanjing 211106, China; 2. Department of Mathematics，Nanjing University of Aeronautics & Astronautics, Jiangsu Province Nanjing 211106, China
Abstract Bioluminescence imaging is a kind of emerging detection technology at cellular, molecular and genetic level. The most popular bioluminescence imaging model is diffusion approximation (DA). However, because of the ill-posedness of the DA-based inverse problem and the instability of reconstruction algorithms, the location accuracy of the reconstructed sources is low. Radiative transfer equation (RTE), which considers the direction of the photon migration and the effect of absorption and scattering in tissues, can accurately express the transmission of bioluminescent photons through the tissues. In this paper, we studied the bioluminescence imaging based on the RTE. 2D simulations were performed, and quantitative evaluation was given by the absolute source position error, the relative source area error and the minimum bounding box. The results of the experiment showed that the imaging quality based on RTE was better than that one based on DA.
KAN Xing-xing,CHEN Chun-xiao,GONG Rong-fang. Quantitative Evaluation of Bioluminescence Imaging Based on Radiative Transfer Equation[J]. Chinese Journal of Biomedical Engineering, 2016, 25(1): 39-46.
[1] Ahamed F, Ahmed F, Guo YH. A high order accurate nodal discontinuous Galerkin method (DGM) for numerical solution of hyperbolic equation[J]. International Journal of Applied Physics and Mathematics, 2012, 2(5):362-364 . [2] Chen L, Gao J. Progress of in vivo optical molecular imaging technology[J]. Medical Recapitulate, 2010, 16(24):3800-3802. [3] Mezzanotte L, Fazzina R, Michelini, et al. In vivo bioluminescence imaging of murine xenograft cancer models with a red-shifted thermostable luciferase[J]. Molecular Imaging and Biology, 2010,12(4):406-414. [4] Kim JB, Urban K, Cochran E. Non-invasive detection of a small number of bioluminescent cancer cells in vivo[J]. Plos One, 2010, 5(2): 9364 -9364. [5] Lim E, Modi KD, Kim J. In vivo bioluminescent imaging of mammary tumors using IVIS spectrum[J]. Journal of Visualized Experiments, 2009, 2009(26): 89-96. [6] Zhang H, Sun L, Xiao E. Intraperitoneal chemotherapy with mitomycin C bound to activated carbon nanoparticles using bioluminescence imaging technology[J]. Military Medical Sciences, 2011, 35(4): 299-302. [7] Li X, Gao K. Comparison between bioluminescent imaging and small-animal PET-CT in Xenotransplantation tumor model with Luc-expressing cell[J]. Chinese Journal of Comparative Medicine, 2011, 21(3):10-13. [8] Kuo C, Troy CTL, Xu L, et al. Three-dimensional reconstruction of in vivo bioluminescent sources based on multispectral imaging[J]. Journal of Biomedical Optics, 2007, 12(2):024007. [9] Kai L, Jie T, Wei Y, et al. Applications of Monte Carlo method in simulating diffuse optical imaging[J]. Journal of Software, 2009, 20(5):1216-1225. [10] Arridge SR, Schweiger M, Hiraoka M. A finite element approach for modeling photon transport in tissue[J]. Medical Physics, 1993, 20(2):299-309. [11] Schweiger M, Arridge SR, Hiraoka M. The finite element method for the propagation of light in scattering media: boundary and source conditions[J]. Medical Physics, 1995, 22(11):1779-1792. [12] Schweiger M, Arridge SR. The finite element method for the propagation of light in scattering media: frequency domain case[J]. Medical Physics, 1997, 24(6):895-902. [13] Arridge SR, Dehghani H, Schweiger M. The finite element model for the propagation of light in scattering media: a direct method for domains with nonscattering regions[J]. Medical Physics, 2000, 27(1):252-264. [14] Cong W, Wang G, Kumar D. Practical reconstruction method for bioluminescence tomography[J]. Optics Express, 2005, 13(18):6756-6771. [15] Yu J, Liu F, He X. Adaptive parameter selection for Tikhonov regularization in bioluminescence tomography[C]. IEEE 2010 3rd International Conference on Biomedical Engineering and Informatics, 2010:86-90. [16] Dehghani H, Davis SC, Pogue BW. Spectrally resolved bioluminescence tomography using the reciprocity approach[J]. Medical Physics, 2008, 35(11):4863-4871. [17] Feng J, Jia K. An optimal permissible source region strategy for multispectral bioluminescence tomography[J]. Optics Express, 2008, 16(20):15641-15654. [18] Yu J, He X, Geng G, et al. Hybrid Multilevel sparse reconstruction for a whole domain bioluminescence tomography using adaptive finite element[J]. Computational and Mathematical Methods in Medicine, 2013, 2013(1):548491. [19] Bal G, Tamasan A. Inverse source problems in transport equations[J]. SIAM Journal on Mathematical Analysis, 2007, 39(1):57-76. [20] Charette A, Boulanger J, Kim HK. An overview on recent radiation transport algorithm development for optical tomography imaging[J]. Journal of Quantitative Spectroscopy & Radiative Transfer, 2008, 109(s17-18): 2743-2766. [21] Hao G, Hongkai Z. Multilevel bioluminescence tomography based on radiative transfer equation Part 1:l1 regularization[J]. Optical Express, 2010, 18(3):1854-1871. [22] Hao G, Hongkai Z. Multilevel bioluminescence tomography based on radiative transfer equation Part 2: total variation and l1 data fidelity[J]. Optical Express, 2010, 18(3):2894-2912. [23] Klose AD, Ntziachristos V, Hielscher AH. The inverse source problem based on the radiative transfer equation in optical molecular imaging[J]. Journal of Computational Physics, 2005, 202(1):323-345 . [24] Han W, Huang J, Eichholz JA. Discrete-ordinate discontinuous Galerkin methods for solving the radiative transport equation[J]. SIAM Journal on Scientific Computing, 2010, 32(2):477-497 . [25] Lewis EE, Miller WF. Computational methods of neutron transport[M]. New York: John Wiley & Sons Inc, 1983. [26] Hao G, Hongkai Z. A fast forward solver of radiative transport equation[J]. Transport Theory and Statistical Physics, 2009, 38(3):149-192 . [27] Carlson BG, Lathrop KD. Computing methods in reactor physics[M] (Eds. H.Greenspan, C.N. Kelber, and D. Okrent).New York: Gordon & Breach, 1968: 167. [28] He X, Zhang Z, Yu J. Trust region method for solving the bioluminescence tomography inverse problem[C]. IEEE 2010 Symposium on Photonics and Optoelectronic (SOPO), 2010:1-4. [29] Zhang B, Yang X, Qin C. A trust region method in adaptive finite element framework for bioluminescence tomography[J]. Optical Society of America, 2010, 18(7):6477-6491. [30] Tikhonov AN, Aresenin VY. Solutions of ill-posed problems[M]. Washington DC: V. H. Winston & Sons, 1977. [31] Hansen PC. Analysis of discrete ill-posed problems by means of the L-curve[J]. SIAM Review, 1992, 34(4):561-580.
