Closed-form solution for the Wigner phase-space distribution function for diffuse reflection and small-angle scattering in a random medium
JOSA A, Vol. 17, Issue 12, pp. 2464-2474 (2000)
http://dx.doi.org/10.1364/JOSAA.17.002464
Acrobat PDF (232 KB)
Abstract
Within the paraxial approximation, a closed-form solution for the Wigner phase-space distribution function is derived for diffuse reflection and small-angle scattering in a random medium. This solution is based on the extended Huygens–Fresnel principle for the optical field, which is widely used in studies of wave propagation through random media. The results are general in that they apply to both an arbitrary small-angle volume scattering function, and arbitrary (real) ABCD optical systems. Furthermore, they are valid in both the single- and multiple-scattering regimes. Some general features of the Wigner phase-space distribution function are discussed, and analytic results are obtained for various types of scattering functions in the asymptotic limit s≫1, where s is the optical depth. In particular, explicit results are presented for optical coherence tomography (OCT) systems. On this basis, a novel way of creating OCT images based on measurements of the momentum width of the Wigner phase-space distribution is suggested, and the advantage over conventional OCT images is discussed. Because all previous published studies regarding the Wigner function are carried out in the transmission geometry, it is important to note that the extended Huygens–Fresnel principle and the ABCD matrix formalism may be used successfully to describe this geometry (within the paraxial approximation). Therefore for completeness we present in an appendix the general closed-form solution for the Wigner phase-space distribution function in ABCD paraxial optical systems for direct propagation through random media, and in a second appendix absorption effects are included.
© 2000 Optical Society of America
OCIS Codes
(170.1650) Medical optics and biotechnology : Coherence imaging
(170.4500) Medical optics and biotechnology : Optical coherence tomography
(170.7050) Medical optics and biotechnology : Turbid media
(290.4210) Scattering : Multiple scattering
Citation
Harold T. Yura, Lars Thrane, and Peter E. Andersen, "Closed-form solution for the Wigner phase-space distribution function for diffuse reflection and small-angle scattering in a random medium," J. Opt. Soc. Am. A 17, 2464-2474 (2000)
http://www.opticsinfobase.org/josaa/abstract.cfm?URI=josaa-17-12-2464
Sort: Year | Journal | Reset
References
- E. P. Wigner, “On the quantum correction for thermodynamic equilibrium,” Phys. Rev. 40, 749–759 (1932).
- M. G. Raymer, C. Cheng, D. M. Toloudis, M. Anderson, and M. Beck, “Propagation of Wigner coherence functions in multiple scattering media,” in Advances in Optical Imaging and Photon Migration, R. R. Alfano and J. G. Fujimoto, eds., Vol. 2 of OSA Trends in Optics and Photonics Series (Optical Society of America, Washington, D.C., 1996), pp. 236–238; C.-C. Cheng and M. G. Raymer, “Long-range saturation of spatial decoherence in wave-field transport in random multiple-scattering media,” Phys. Rev. Lett. 82, 4807–4810 (1999); M. G. Raymer and C.-C. Cheng, “Propagation of the optical Wigner function in random multiple-scattering media,” in Laser–Tissue Interaction XI: Photochemical, Photothermal, and Photomechanical, D. D. Duncan, J. O. Hollinger, and S. L. Jacques, eds., Proc. SPIE 3914, 376–380 (2000).
- S. John, G. Pang, and Y. Yang, “Optical coherence propagation and imaging in a multiple scattering medium,” J. Biomed. Opt. 1, 180–191 (1996).
- A. Wax and J. E. Thomas, “Measurement of smoothed Wigner phase-space distributions for small-angle scattering in a turbid medium,” J. Opt. Soc. Am. A 15, 1896–1908 (1998).
- R. F. Lutomirski and H. T. Yura, “Propagation of a finite optical beam in an inhomogeneous medium,” Appl. Opt. 10, 1652–1658 (1971).
- M. Hillery, R. F. O’Connell, M. O. Scully, and E. P. Wigner, “Distribution functions in physics: fundamentals,” Phys. Rep. 106, 121–167 (1984).
- It is straightforward to show that within the paraxial approximation, the specific radiance distribution N(P, θ)= k^{2} W(P, kθ) for those cases where the Wigner phase-space distribution is positive definite.
- In random media where the scattering particles are large compared with the wavelength and the index of refraction ratio is near unity, the bulk backscattering efficiency is much smaller than the scattering efficiency. Moreover, the scattering is primarily in the forward direction, which is the basis of using the paraxial approximation. Therefore the bulk backscattering may be neglected when one is considering the light propagation problem, since its contribution is mall. An example of this is skin tissue (cell sizes of 5–10-μm diameter and index of refraction ratio of 1.45/1.4=1.04).
- V. A. Banakh and V. L. Mironov, LIDAR in a Turbulent Atmosphere (Artech House, Boston, Mass., 1987).
- L. Thrane, H. T. Yura, and P. E. Andersen, “Analysis of optical coherence tomography systems based on the extended Huygens–Fresnel principle,” J. Opt. Soc. Am. A 17, 484–490 (2000).
- H. T. Yura and S. G. Hanson, “Optical beam wave propagation through complex optical systems,” J. Opt. Soc. Am. A 4, 1931–1948 (1987).
- Because we are dealing with “real” ABCD optical systems, we tacitly assume that B≠0.
- V. I. Tatarskii, The Effects of the Turbulent Atmosphere on Wave Propagation (National Technical Information Service, Springfield, Va., 1971).
- H. T. Yura and S. G. Hanson, “Effects of receiver optics contamination on the performance of laser velocimeter systems,” J. Opt. Soc. Am. A 13, 1891–1902 (1996).
- One can obtain corresponding results that are valid for propagation through the turbulent atmosphere by formally replacing the exponent in Eq. (4) by one-half the corresponding point-source wave structure function.13
- For completeness, the corresponding mutual coherence function is given by 〈U(P +p /2)U* (P +p /2)〉 =(η/πB^{2})|Γ_{pt} (p)|^{2} K(−p)exp (−ikDp ⋅P /B).
- J. W. Goodman, Statistical Optics (Wiley, New York, 1985).
- D. Arnush, “Underwater light-beam propagation in the small-angle-scattering approximation,” J. Opt. Soc. Am. 62, 1109–1111 (1972).
- I. Dror, A. Sandrov, and N. S. Kopeika, “Experimental investigation of the influence of the relative position of the scattering layer on image quality: the shower curtain effect,” Appl. Opt. 37, 6495–6499 (1998).
- L. G. Henyey and J. L. Greenstein, “Diffuse radiation in the galaxy,” Astrophys. J. 93, 70–83 (1941).
- See, e.g., S. L. Jacques and L. Wang, “Monte Carlo modeling of light transport in tissues,” Optical-Thermal Response of Laser-Irradiated Tissue, A. J. Welch and M. J. C. van Gemert, eds. (Plenum, New York, 1995), Chap. 4.
- C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983).
- It is straightforward to show that the quantity θ_{rms} appearing in Eq. (15) can be written as θ_{rms} =2(1−g).
- J. M. Schmitt and A. Knüttel, “Model of optical coherence tomography of heterogeneous tissue,” J. Opt. Soc. Am. A 14, 1231–1242 (1997).
- M. J. C. van Gemert, S. L. Jacques, H. J. C. M. Sterenborg, and W. M. Star, “Skin optics,” IEEE Trans. Biomed. Eng. 36, 1146–1154 (1989).
- H. T. Yura, “A multiple scattering analysis of the propagation of radiance through the atmosphere,” in Proceedings of the Union Radio-Scientifique Internationale Open Symposium (Union Radio-Scientifique Internationale, Ghent, Belgium, 1977), pp. 65–69.
- A. E. Siegman, Lasers (Oxford U. Press, Oxford, UK, 1986).
- For axially symmetric scattering functions, the integral over scattering angles reduces to 2π0 ∞ d θθσ(θ)J_{0} kz^{′} Z r +1−z^{′} Z p, where J_{0} (⋅) is the Bessel function of the first kind, of order zero.
- For axially symmetric scattering functions, the integral over scattering angles reduces to2π0 ∞ d θθσ(θ)J_{0} kB(z^{′}) B r +1−B(z^{′}) B p.
- For a spatially uniform medium of index of refraction n, we have A=D=1, B=Z/n, B(z^{′})=z^{′} /n, and Eq. (A7) reduces to Eq. (A1).
- Note that the limits on the z^{′} integration of the jth term in the summation are now from z_{j} to z_{j} +Δz_{j}.
- The corresponding Wigner function is obtained by replacing θ by q/k_{m}.
- As expected physically, the corresponding irradiance, obtained by integrating the radiance pattern over all solid angle, is given by exp (−μ_{A} Z)I_{0}.
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.
« Previous Article | Next Article »
OSA is a member of CrossRef.