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.
Perturbation theory for the diffusion equation by use of the moments of the generalized temporal point-spread function. I. Theory / A. Sassaroli; F. Martelli; S. Fantini. - In: JOURNAL OF THE OPTICAL SOCIETY OF AMERICA. A, OPTICS, IMAGE SCIENCE, AND VISION. - ISSN 1084-7529. - STAMPA. - 23:(2006), pp. 2105-2118. [10.1364/JOSAA.23.002105]
Perturbation theory for the diffusion equation by use of the moments of the generalized temporal point-spread function. I. Theory
MARTELLI, FABRIZIOMembro del Collaboration Group
;
2006
Abstract
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.File | Dimensione | Formato | |
---|---|---|---|
josa_pert1.pdf
Accesso chiuso
Descrizione: Articolo principale
Tipologia:
Pdf editoriale (Version of record)
Licenza:
Tutti i diritti riservati
Dimensione
172.21 kB
Formato
Adobe PDF
|
172.21 kB | Adobe PDF | Richiedi una copia |
I documenti in FLORE sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.