Three-dimensional Neumann-series approach to model light transport in nonuniform media

Abhinav K. Jha, Matthew A. Kupinski, Harrison H. Barrett, Eric Clarkson, John H. Hartman

Research output: Contribution to journalArticlepeer-review

18 Scopus citations


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.

Original languageEnglish (US)
Pages (from-to)1885-1899
Number of pages15
JournalJournal of the Optical Society of America A: Optics and Image Science, and Vision
Issue number9
StatePublished - Sep 1 2012
Externally publishedYes

ASJC Scopus subject areas

  • Electronic, Optical and Magnetic Materials
  • Atomic and Molecular Physics, and Optics
  • Computer Vision and Pattern Recognition


Dive into the research topics of 'Three-dimensional Neumann-series approach to model light transport in nonuniform media'. Together they form a unique fingerprint.

Cite this