It looks like you're using Internet Explorer 6: Features on this site are not supported by that browser version. Please upgrade to the latest version of Internet Explorer.

You are viewing a version of the Ocean Optics Web Book site specially-formatted for printing.

# level 2

## Beam and Point Spread Functions

 Page updated: Aug 3, 2018 Principal author: Curtis Mobley
This page defines beam and point spread functions and introduces their use in image analysis.

#### The Beam Spread Function

Consider a well collimated source emitting spectral power in direction , as illustrated in Fig. . As the beam passes through the medium, scattering will spread out the beam as illustrated by the green arrows in the figure. Absorption will reduce the beam power. The combined effects of scattering and absorption give some spectral irradiance on the surface of a sphere of radius at an angle of relative to the direction of the emitted beam. The beam spread function (BSF) is then defined as the detected irradiance normalized by the emitted power:

 (1)

The irradiance sensor used to measure is a plane irradiance sensor, which has a cosine response for angles relative to the normal to the detector surface.

#### The Point Spread Function

Now consider a source that emits spectral intensity in a cosine pattern as illustrated by the short green arrows in Fig. . The angular distribution of the emitted intensity is

Integrating this intensity over all solid angles shows that the total emitted spectral power is :

Here is the azimuthal angle of the emitted intensity, which is independent of .

As illustrated by the long green arrows in Fig. , the emitted intensity will give rise to a radiance , where is the distance from the source and is direction measured from the axis of the emitted intensity.

With these definitions, the point spread function (PSF) is defined as the detected radiance normalized by the maximum of the emitted intensity:

 (2)

Figure shows a photograph of a street light seen through a foggy atmosphere. Although a street light is not a cosine-emitting point source, the angular pattern of the "glow" seen around the light gives a qualitative feeling for the PSF, which is just the pattern of the radiance distribution centered on the source.

#### Equivalence of the PSF and BSF

Figures and are drawn to highlight the symmetry between the BSF and the PSF. The principle of electromagnetic reciprocity says that if light follows a particular path from one point A to another point B, then light leaving B in the opposite direction will retrace the original path back to A. This is illustrated by the green arrows connecting the sources and receivers in the two figures. This reciprocity suggests that the BSF and PSF contain equivalent information. Figure shows PSF and BSF measurements from two very similar water bodies obtained by Mertens and Replogle (1977). The instruments were mounted on a 20 m long underwater frame, which allowed a maximum range between source and detector of approximately 19 m. Although the water IOPs were slightly different for these two measurements, the closeness of the PSF and BSF over several orders of magnitude suggest that the BSF and PSF are numerically equal. Mertens and Replogle stated this equality without proof. Gordon (1994b) started with the radiative transfer equation and used reciprocity to show that the BSF and PSF are indeed numerically equal.

It is important to understand the difference between the volume scattering function (VSF) and the beam or point spread function. The VSF describes the redirection of light by a single scattering event. The VSF depends only on the scattering angle. The PSF describes the cumulative effects of scattering and absorption as the light travels from the source to the detector. The PSF therefore depends on both the scattering and absorption properties of the medium, and on the distance from the source to the detector.

#### The PSF and Image Prediction

The PSF plays a fundamental role in the prediction of how an object appears when seen through an absorbing and scattering medium such as water. Consider the image of a shark seen in Fig. . Let represent the "bright-dark" pattern of the shark image when seen at distance , or as seen through a vacuum. Thus for a gray-scale digital image and 8 bit resolution, for a black pixel and for white pixel. When viewing a particular point in the image through a vacuum at some distance from the shark, the detected radiance comes only from point and has the value . However, when viewing point in the through a scattering medium, every point in the scene contributes to the radiance seen at . The concentric white circles in the figure represent contours of the . Point thus contributes a radiance of , where angle . Likewise point contributes radiance , where angle . The total radiance at point is given by summing up the contributions from all points in the image:

 (3)

where the PSF argument is shorthand for distance and angle

when the shark is seen from distance . Thus, given the actual appearance of an object and the PSF of the medium, the appearance of the object as seen through the medium can be computed.

Integrals of the form are known as convolution integrals. Equation () is a 2D convolution. These integrals have the property that the Fourier transform of a convolution is the product of the Fourier transforms of the two functions. This leads to a very important procedure for image analysis. Rather than compute the double integral of Eq. () for every point of the scene, the Fourier transforms of and can be computed, multiplied, and then the inverse transform of the product gives the final image. This topic will be developed on other pages; meanwhile see Goodman (1996). (Incidentally, the term "convolution" is frequently misused. An integral of the form is not a convolution; it is just the integral of the product of two functions and or, if you wish, of one function weighed by the other. Such integrals do not occur in image analysis and do not have a convenient Fourier transform.)

Finally, note that in a vacuum and with a perfect optical system, the PSF reduces to a 2D Dirac delta function: , in which case

Thus the image at distance is the same as at distance 0.

Because of the importance of the PSF in image analysis, considerable effort has been expended to develop models of the PSF as a function of water IOPs. A number of such models are reviewed in Hou et al. (2008) The page on The Lidar Equation shows the fundamental role of the BSF in lidar remote sensing.