Jha Abhinav K, Kupinski Matthew A, Barrett Harrison H, Clarkson Eric, Hartman John H
College of Optical Sciences, University of Arizona, Tucson, Arizona 85721, USA.
J Opt Soc Am A Opt Image Sci Vis. 2012 Sep 1;29(9):1885-99. doi: 10.1364/JOSAA.29.001885.
We present the implementation, validation, and performance of a three-dimensional (3D) Neumann-series approach to model photon propagation in nonuniform media using the radiative transport equation (RTE). The RTE is implemented for nonuniform scattering media in a spherical harmonic basis for a diffuse-optical-imaging setup. The method is parallelizable and implemented on a computing system consisting of NVIDIA Tesla C2050 graphics processing units (GPUs). The GPU implementation provides a speedup of up to two orders of magnitude over non-GPU implementation, which leads to good computational efficiency for the Neumann-series method. The results using the method are compared with the results obtained using the Monte Carlo simulations for various small-geometry phantoms, and good agreement is observed. We observe that the Neumann-series approach gives accurate results in many cases where the diffusion approximation is not accurate.
我们展示了一种三维（3D）诺伊曼级数方法的实现、验证及其性能，该方法使用辐射传输方程（RTE）对非均匀介质中的光子传播进行建模。在漫射光学成像设置中，RTE在球谐基下针对非均匀散射介质实现。该方法可并行化，并在由英伟达Tesla C2050图形处理单元（GPU）组成的计算系统上实现。与非GPU实现相比，GPU实现提供了高达两个数量级的加速，这为诺伊曼级数方法带来了良好的计算效率。使用该方法得到的结果与使用蒙特卡罗模拟针对各种小尺寸几何模型得到的结果进行了比较，观察到两者吻合良好。我们发现，在许多扩散近似不准确的情况下，诺伊曼级数方法都能给出准确的结果。
