Abstract
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.
© 2012 Optical Society of America
Full Article | PDF ArticleMore Like This
Abhinav K. Jha, Matthew A. Kupinski, Takahiro Masumura, Eric Clarkson, Alexey V. Maslov, and Harrison H. Barrett
J. Opt. Soc. Am. A 29(8) 1741-1757 (2012)
Abhinav K. Jha, Yansong Zhu, Simon Arridge, Dean F. Wong, and Arman Rahmim
Biomed. Opt. Express 9(4) 1389-1407 (2018)
Aki Pulkkinen and Tanja Tarvainen
J. Opt. Soc. Am. A 30(3) 470-478 (2013)