To model bidirectional measurements, a numerical method for computing the incoherent light scattered by a diffusing medium is presented. The results are expressed as a function of the incident and the observer angles (θi, ϕi) and (θf, ϕf), in contrast to the N-flux method, which gives no information about the azimuthal distribution. To solve the multiple-scattering equations, an auxiliary function, expanded on the spherical harmonics, is introduced in the diffusion equation. A set of integral equations on the coefficients are obtained that are well suited for the numerical resolution. The boundary conditions are included in the linear operator of the integral equation, so that each boundary condition is associated with a specific equation. As an illustration, the method is applied to numerical simulations of maps of the light scattered by a thick refractive diffusing layer of refractive index n=1.5, for two directive phase functions and for several incident collimated-beam angles.
© 2002 Optical Society of AmericaPDF Article