OSA's Digital Library

Journal of the Optical Society of America A

Journal of the Optical Society of America A


  • Editor: Stephen A. Burns
  • Vol. 23, Iss. 9 — Sep. 1, 2006
  • pp: 2105–2118

Perturbation theory for the diffusion equation by use of the moments of the generalized temporal point-spread function. I. Theory

Angelo Sassaroli, Fabrizio Martelli, and Sergio Fantini  »View Author Affiliations

JOSA A, Vol. 23, Issue 9, pp. 2105-2118 (2006)

View Full Text Article

Enhanced HTML    Acrobat PDF (172 KB)

Browse Journals / Lookup Meetings

Browse by Journal and Year


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools



We approach the perturbative solution to the diffusion equation for the case of absorbing inclusions embedded in a heterogeneous scattering medium by using general properties of the radiative transfer equation and the solution of the Fredholm equation of the second kind given by the Neumann series. The terms of the Neumann series are used to obtain the expression of the moments of the generalized temporal point-spread function derived in transport theory. The moments are calculated independently by using Monte Carlo simulations for validation of the theory. While the mixed moments are correctly derived from the theory by using the solution of the diffusion equation in the geometry of interest, in order to obtain the self moments we should reframe the problem in transport theory and use a suitable solution of the radiative transfer equation for the calculation of the multiple integrals of the corresponding Neumann series. Since the rigorous theory leads to impractical formulas, in order to simplify and speed up the calculation of the self moments, we propose a heuristic method based on the calculation of only a single integral and some scaling parameters. We also propose simple quadrature rules for the calculation of the mixed moments for speeding up the computation of perturbations due to multiple defects. The theory can be developed in the continuous-wave domain, the time domain, and the frequency domain. In a companion paper [ J. Opt. Soc. Am. A 23, 2119 (2006) ] we discuss the conditions of applicability of the theory in practical cases found in diffuse optical imaging of biological tissues.

© 2006 Optical Society of America

OCIS Codes
(170.3660) Medical optics and biotechnology : Light propagation in tissues
(170.5280) Medical optics and biotechnology : Photon migration
(290.7050) Scattering : Turbid media
(300.1030) Spectroscopy : Absorption

ToC Category:
Medical Optics and Biotechnology

Original Manuscript: October 25, 2005
Revised Manuscript: March 8, 2006
Manuscript Accepted: March 25, 2006

Virtual Issues
Vol. 1, Iss. 10 Virtual Journal for Biomedical Optics

Angelo Sassaroli, Fabrizio Martelli, and Sergio Fantini, "Perturbation theory for the diffusion equation by use of the moments of the generalized temporal point-spread function. I. Theory," J. Opt. Soc. Am. A 23, 2105-2118 (2006)

Sort:  Author  |  Year  |  Journal  |  Reset  


  1. S. R. Arridge, M. Cope, and D. T. Delpy, "The theoretical basis for the determination of optical pathlength in tissue: temporal and frequency analysis," Phys. Med. Biol. 37, 1531-1560 (1992). [CrossRef] [PubMed]
  2. A. Kienle, "Light diffusion through a turbid parallelepiped," J. Opt. Soc. Am. A 22, 1883-1888 (2005). [CrossRef]
  3. F. Martelli, A. Sassaroli, S. Del Bianco, Y. Yamada, and G. Zaccanti, "Solution of the time-dependent diffusion equation for layered random media by the eigenfunction method," Phys. Rev. E 67, 056623 (2003). [CrossRef]
  4. D. A. Boas, M. A. O'Leary, B. Chance, and A. G. Yodh, "Scattering of diffuse photon density waves by spherical inhomogeneities within turbid media: Analytic solution and applications," Proc. Natl. Acad. Sci. U.S.A. 91, 4887-4891 (1994). [CrossRef] [PubMed]
  5. D. A. Boas, M. A. O'Leary, B. Chance, and A. G. Yodh, "Detection and characterization of optical inhomogeneities with diffuse photon density waves: a signal to noise analysis," Appl. Opt. 36, 75-92 (1997). [CrossRef] [PubMed]
  6. P. N. den Outer, Th. M. Nieuwenhuizen, and A. Lagendijk, "Location of objects in multiple-scattering media," J. Opt. Soc. Am. A 10, 1209-1218 (1993). [CrossRef]
  7. K. M. Case, F. de Hoffman, and G. Placzek, Introduction to the Theory of Neutron Diffusion (U.S. Government Printing Office, 1953).
  8. R. Graaf and K. Rinzema, "Practical improvements on photon diffusion theory: application to isotropic scattering," Phys. Med. Biol. 46, 3043-3050 (2001). [CrossRef]
  9. J. C. J. Paasschens, "Solution of the time-dependent Boltzmann equation," Phys. Rev. E 56, 1135-1141 (1997). [CrossRef]
  10. R. Graaf, J. G. Aarnoudse, F. F. de Mul, and H. W. Jentink, "Light propagation parameters for anisotropically scattering media based on a rigorous solution of the transport equation," Appl. Opt. 26, 2273-2279 (1989). [CrossRef]
  11. K. Rinzema, B. J. Hoenders, H. A. Ferwerda, and J. J. ten Bosch, "Analytic calculation of the radiance in an anisotropically scattering turbid medium close to a source," Pure Appl. Opt. 4, 629-642 (1995). [CrossRef]
  12. K. Rinzema, L. H. P. Murrer, and W. M. Star, "Direct experimental verification of light transport theory in an optical phantom," J. Opt. Soc. Am. A 15, 2078-2088 (1998). [CrossRef]
  13. B. W. Pogue, M. S. Patterson, H. Jiang, and K. D. Paulsen, "Initial assessment of a simple system for frequency-domain diffuse optical tomography," Phys. Med. Biol. 40, 1709-1729 (1995). [CrossRef] [PubMed]
  14. Y. Yao, Y. Wang, Y. Pei, W. Zhu, and R. L. Barbour, "Frequency domain optical imaging of absorption and scattering distributions by a Born iterative method," J. Opt. Soc. Am. A 14, 325-342 (1997). [CrossRef]
  15. H. Jiang, K. D. Paulsen, U. L. Osterberg, B. W. Pogue, and M. S. Patterson, "Optical image reconstruction using frequency-domain data: simulations and experiments," J. Opt. Soc. Am. A 13, 253-266 (1996). [CrossRef]
  16. S. R. Arridge, M. Schweiger, M. Hiraoka, and D. T. Delpy, "A finite element approach to modelling photon transport in tissue," Med. Phys. 20, 299-309 (1993). [CrossRef] [PubMed]
  17. S. R. Arridge and M. Schweiger, "Photon-measurement density functions. Part 2. Finite-element-method calculations," Appl. Opt. 34, 8026-8037 (1995). [CrossRef] [PubMed]
  18. H. Zhao, F. Gao, Y. Tanikawa, K. Homma, and Y. Yamada, "Time-resolved diffuse optical tomography imaging for the provision of both anatomical and functional information about biological tissue," Appl. Opt. 44, 1905-1916 (2005). [CrossRef] [PubMed]
  19. R. L. Barbour, H. L. Graber, Y. Pei, S. Zhong, and C. H. Schmitz, "Optical tomographic imaging of dynamic features of dense-scattering media," J. Opt. Soc. Am. A 18, 3018-3036 (2001). [CrossRef]
  20. K. Ren, G. S. Abdoulaev, G. Bal, and A. H. Hielscher, "Algorithm for solving the equation of radiative transfer in the frequency domain," Opt. Lett. 29, 578-580 (2004). [CrossRef] [PubMed]
  21. G. S. Abdoulaev and A. H. Hielscher, "Three-dimensional optical tomography with the equation of radiative transfer," J. Electron. Imaging 12, 594-600 (2003). [CrossRef]
  22. R. Graaf, M. H. Koelink, F. F. M. de Mul, W. G. Zijlstra, A. C. M. Dassel, and J. G. Aarnoudse, "Condensed Monte Carlo simulations for the description of light transport," Appl. Opt. 32, 426-434 (1993). [CrossRef]
  23. B. C. Wilson and G. Adam, "A Monte Carlo model for the absorption and flux distributions of light in tissue," Med. Phys. 10, 824-830 (1983). [CrossRef] [PubMed]
  24. L. Wang, S. L. Jacques, and L. Zheng, "MCML-Monte Carlo modeling of light transport in multi-layered tissues," Comput. Methods Programs Biomed. 47, 131-146 (1995). [CrossRef] [PubMed]
  25. A. Sassaroli, C. Blumetti, F. Martelli, L. Alianelli, D. Contini, A. Ismaelli, and G. Zaccanti, "Monte Carlo procedure for investigating light propagation and imaging of highly scattering media," Appl. Opt. 37, 7392-7400 (1998). [CrossRef]
  26. S. R. Arridge and M. Schweiger, "Photon-measurement density functions. Part 1. Analytical forms," Appl. Opt. 34, 7395-7409 (1995). [CrossRef] [PubMed]
  27. J. C. Hebden and S. R. Arridge, "Imaging through scattering media by the use of an analytical model of perturbation amplitudes in the time domain," Appl. Opt. 35, 6788-6796 (1996). [CrossRef] [PubMed]
  28. M. A. O'Leary, D. A. Boas, B. Chance, and A. G. Yodh, "Experimental images of inhomogeneous turbid media by frequency-domain diffusing- photon tomography," Opt. Lett. 20, 426-428 (1995). [CrossRef] [PubMed]
  29. D. A. Boas, "A fundamental limitation of linearized algorithms for diffuse optical tomography," Opt. Express 1, 404-413 (1997). [CrossRef] [PubMed]
  30. M. R. Ostermeyer and S. L. Jacques, "Perturbation theory for diffuse light transport in complex biological tissues," J. Opt. Soc. Am. A 14, 255-261 (1997). [CrossRef]
  31. A. Sassaroli, F. Martelli, and S. Fantini, "Perturbation theory for the diffusion equation by use of the moments of the generalized temporal point-spread function. II. Continuous-wave results," J. Opt. Soc. Am. A 23, 2119-2131 (2006). [CrossRef]
  32. M. Hiraoka, M. Firbank, M. Essenpreis, M. Cope, S. R. Arridge, P. van der Zee, and D. T. Delpy, "A Monte Carlo investigation of optical path length in inhomogeneous tissue and its application to near-infrared spectroscopy," Phys. Med. Biol. 38, 1859-1876 (1993). [CrossRef] [PubMed]
  33. Y. Tsuchiya, "Photon path distribution in inhomogeneous turbid media: theoretical analysis and a method of calculation," J. Opt. Soc. Am. A 19, 1383-1389 (2002). [CrossRef]
  34. J. Steinbrink, H. Wabnitz, H. Obrig, A. Villringer, and H. Rinneberg, "Determining changes in NIR absorption using a layered model of the human head," Phys. Med. Biol. 46, 879-896 (2001). [CrossRef] [PubMed]
  35. G. Strangman, M. A. Franceschini, and D. A. Boas, "Factors affecting the accuracy of near-infrared spectroscopy concentration calculations for focal changes in oxygenation parameters," Neuroimage 18, 865-879 (2003). [CrossRef] [PubMed]
  36. A. Ishimaru, Wave Propagation and Scattering in Random Media (Academic, 1978).
  37. K. M. Case and P. F. Zweifel, Linear Transport Theory (Addison-Wesley, 1967).
  38. J. J. Duderstadt and W. R. Martin, Transport Theory (Wiley, 1979).
  39. S. Carraresi, T. S. Mohamed Shatir, F. Martelli, and G. Zaccanti, "Accuracy of a perturbation model to predict the effect of scattering and absorbing inhomogeneities on photon migration," Appl. Opt. 40, 4622-4632 (2001). [CrossRef]
  40. K. Furutsu and Y. Yamada, "Diffusion approximation for a dissipative random medium and the applications," Phys. Rev. E 50, 3634-3640 (1994). [CrossRef]
  41. G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists, 4th ed. (Academic, 1995), p. 934.
  42. G. Zaccanti, L. Alianelli, C. Blumetti, and C. Carraresi, "Method for measuring the mean time of flight spent by photons inside a volume element of a highly diffusing medium," Opt. Lett. 24, 1290-1292 (1999). [CrossRef]
  43. R. Aronson, "Radiative transfer implies a modified reciprocity relation," J. Opt. Soc. Am. A 14, 486-490 (1997). [CrossRef]
  44. S. Feng, F. Zeng, and B. Chance, "Photon migration in the presence of a single defect: a perturbation analysis," Appl. Opt. 34, 3826-3837 (1995). [CrossRef] [PubMed]
  45. M. Bassani, F. Martelli, and G. Zaccanti, "Independence of the diffusion coefficient from absorption: experimental and numerical evidence," Opt. Lett. 22, 853-855 (1997). [CrossRef] [PubMed]
  46. D. A. Boas, J. P. Culver, J. J. Stott, and A. K. Dunn, "Three-dimensional Monte Carlo code for photon migration through complex heterogeneous media including the adult human head," Opt. Express 10, 159-170 (2002). [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


Fig. 1

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited