光子学报, 2010, 39 (9): 1594, 网络出版: 2010-11-04
二维稳态辐射传输方程的有限差分求解法
A Finite-Difference-Method Solution to Two-dimensional Steady-State Radiative Transfer Equation
辐射传输方程 有限差分方法 扩散光学层析 正向问题 Radiative transfer equation Finite difference method Diffuse optical tomography Forward model
摘要
针对扩散光学层析在小动物成像中的应用问题并基于混浊介质空间光子三维散射的实际物理效应,提出的二维稳态辐射传输方程的有限差分数值求解新方法.在此基础上,研究了不同的空间剖分网格和角度离散密度对求解准确度的影响,并通过将所提方法与蒙特卡洛模拟进行比对,验证方法的正确性.研究表明:在均匀组织体内,当离散角度达到一定数量时,由辐射传输方程的有限差分解获得的透射面和侧面的光子密度对空间网格大小并不敏感,而在反射面上光子密度计算则需要较密的空间网格才能够达到一定准确度.本研究为发展基于辐射传输方程的扩散光学层析理论奠定了基础.
Abstract
A novel finite difference method (FDM) solution to the 2-dimensional time-independent radiative transfer equation (RTE) is proposed for small-animal imaging applications, which is based on the real physics of 3-dimensional scattering of a photon in turbid medium. The influences of the angular discretization and the spatially meshing grids on the accuracy of FDM-RTE solution are investigated and a comparison is made between the FDM-RTE calculations to the relevant Monte-Carlo simulation ones for a validation of the proposed method. The comparison shows that, as the number of the discretized solid angles is over 80, the accuracy of the FDM-RTE calculated forward photon density is insensitive to the spatially meshing grids on the transmission and side surfaces but on the reflection surface. This study provides a viable methodology to the development of the RTE-based Diffuse Optical Tomography(DOT) technology.
金蒙, 高峰, 李娇, 赵会娟. 二维稳态辐射传输方程的有限差分求解法[J]. 光子学报, 2010, 39(9): 1594. JIN Meng, GAO Feng, LI Jiao, ZHAO Hui-juan. A Finite-Difference-Method Solution to Two-dimensional Steady-State Radiative Transfer Equation[J]. ACTA PHOTONICA SINICA, 2010, 39(9): 1594.