Light transport in turbid media with non-scattering, low-scattering and high absorption heterogeneities based on hybrid simplified spherical harmonics with radiosity model.

Modeling light propagation in the whole body is essential and necessary for optical imaging. However, non-scattering, low-scattering and high absorption regions commonly exist in biological tissues, which lead to inaccuracy of the existing light transport models. In this paper, a novel hybrid light transport model that couples the simplified spherical harmonics approximation (SPN) with the radiosity theory (HSRM) was presented, to accurately describe light transport in turbid media with non-scattering, low-scattering and high absorption heterogeneities. In the model, the radiosity theory was used to characterize the light transport in non-scattering regions and the SPN was employed to handle the scattering problems, including subsets of low-scattering and high absorption. A Neumann source constructed by the light transport in the non-scattering region and formed at the interface between the non-scattering and scattering regions was superposed into the original light source, to couple the SPN with the radiosity theory. The accuracy and effectiveness of the HSRM was first verified with both regular and digital mouse model based simulations and a physical phantom based experiment. The feasibility and applicability of the HSRM was then investigated by a broad range of optical properties. Lastly, the influence of depth of the light source on the model was also discussed. Primary results showed that the proposed model provided high performance for light transport in turbid media with non-scattering, low-scattering and high absorption heterogeneities.

[1]  S Arridge,et al.  3D optical tomography in the presence of void regions. , 2000, Optics express.

[2]  Simon R. Arridge,et al.  Variable order spherical harmonic expansion scheme for the radiative transport equation using finite elements , 2011, J. Comput. Phys..

[3]  Ge Wang,et al.  A mouse optical simulation environment (MOSE) to investigate bioluminescent phenomena in the living mouse with the Monte Carlo method. , 2004, Academic radiology.

[4]  Stephen J. Wright,et al.  Reconstruction in optical tomography using the PN approximations , 2006 .

[5]  Jie Tian,et al.  Multimodality Molecular Imaging , 2008, IEEE Engineering in Medicine and Biology Magazine.

[6]  Ville Kolehmainen,et al.  Hybrid radiative-transfer-diffusion model for optical tomography. , 2005, Applied optics.

[7]  A. Chatziioannou,et al.  Tomographic bioluminescence imaging by use of a combined optical-PET (OPET) system: a computer simulation feasibility study , 2005, Physics in medicine and biology.

[8]  Jorge Ripoll,et al.  Tomographic imaging with polarized light. , 2012, Journal of the Optical Society of America. A, Optics, image science, and vision.

[9]  R. Leahy,et al.  Digimouse: a 3D whole body mouse atlas from CT and cryosection data , 2007, Physics in medicine and biology.

[10]  S R Arridge,et al.  An investigation of light transport through scattering bodies with non-scattering regions. , 1996, Physics in medicine and biology.

[11]  S. Arridge,et al.  Photon migration in non-scattering tissue and the effects on image reconstruction. , 1999, Physics in medicine and biology.

[12]  J. Ripoll,et al.  Free-space propagation of diffuse light: theory and experiments. , 2003, Physical review letters.

[13]  Jie Tian,et al.  Spectrally resolved bioluminescence tomography with the third-order simplified spherical harmonics approximation. , 2009, Physics in medicine and biology.

[14]  Jie Tian,et al.  Evaluation of the simplified spherical harmonics approximation in bioluminescence tomography through heterogeneous mouse models. , 2010, Optics express.

[15]  S. Arridge,et al.  Coupled radiative transfer equation and diffusion approximation model for photon migration in turbid medium with low-scattering and non-scattering regions , 2005, Physics in medicine and biology.

[16]  S. Arridge,et al.  Optical tomography: forward and inverse problems , 2009, 0907.2586.

[17]  Jie Tian,et al.  GPU-based Monte Carlo simulation for light propagation in complex heterogeneous tissues. , 2010, Optics express.

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

[19]  Jie Tian,et al.  A multilevel adaptive finite element algorithm for bioluminescence tomography. , 2006, Optics express.

[20]  Huabei Jiang,et al.  A higher order diffusion model for three-dimensional photon migration and image reconstruction in optical tomography , 2009, Physics in medicine and biology.

[21]  V A Markel,et al.  Inverse scattering for the diffusion equation with general boundary conditions. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[22]  Y Liu,et al.  Flux vector formulation for photon propagation in the biological tissue. , 2007, Optics letters.

[23]  Vasilis Ntziachristos,et al.  Looking and listening to light: the evolution of whole-body photonic imaging , 2005, Nature Biotechnology.

[24]  Yujie Lu,et al.  Experimental bioluminescence tomography with fully parallel radiative-transfer-based reconstruction framework. , 2009, Optics express.

[25]  M. Schweiger,et al.  The finite element model for the propagation of light in scattering media: a direct method for domains with nonscattering regions. , 2000, Medical physics.

[26]  Jie Tian,et al.  Comparisons of hybrid radiosity-diffusion model and diffusion equation for bioluminescence tomography in cavity cancer detection. , 2012, Journal of biomedical optics.

[27]  Alexander D Klose,et al.  The forward and inverse problem in tissue optics based on the radiative transfer equation: a brief review. , 2010, Journal of quantitative spectroscopy & radiative transfer.

[28]  Stefan Andersson-Engels,et al.  Evaluation of a radiative transfer equation and diffusion approximation hybrid forward solver for fluorescence molecular imaging , 2012, Journal of biomedical optics.

[29]  Arridge,et al.  Optical tomography in the presence of void regions , 2000, Journal of the Optical Society of America. A, Optics, image science, and vision.

[30]  Jie Tian,et al.  Influence investigation of a void region on modeling light propagation in a heterogeneous medium. , 2013, Applied optics.

[31]  E. Okada,et al.  Hybrid Monte Carlo-diffusion method for light propagation in tissue with a low-scattering region. , 2003, Applied optics.

[32]  J. Willmann,et al.  Molecular imaging in drug development , 2008, Nature Reviews Drug Discovery.

[33]  Geoffrey McLennan,et al.  Practical reconstruction method for bioluminescence tomography. , 2005, Optics express.

[34]  Edward W. Larsen,et al.  Light transport in biological tissue based on the simplified spherical harmonics equations , 2006, J. Comput. Phys..

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

[36]  Jae Hoon Lee,et al.  Modeling of diffuse-diffuse photon coupling via a nonscattering region: a comparative study. , 2004, Applied optics.

[37]  Jie Tian,et al.  Graphics processing unit parallel accelerated solution of the discrete ordinates for photon transport in biological tissues. , 2011, Applied optics.