Estimation of a remote sensing system point-spread function from measured imagery

A general approach to identifying the point spread function (PSF) of a remotely sensed scene is demonstrated in terms of a step function for an abrupt change in the gray level along the row or column of the image data. The estimate of the PSF is made in terms of a finite sum of basis functions, employing a sequence of rectangular pulses covering the spatial extent of the PSF. The approximation, if narrow impulses are employed, provides accurate fidelity to the PSF. The method becomes practical when the geometrical structure of the scene elements producing the measured response is known. The field boundary is obtained through consideration of the differing intensities on each side of the boundary, which is a step discontinuity. The mathematical procedure for the technique is provided, together with a sample problem from Landsat-4 Thematic Mapper data. Atmospheric blurring and electronic effects on the overall PSF and the cubic convolution resampling effects are noted.