Surface Roughness Measuring System

This invention relates to novel apparatus for measuring heights of peaks above an arbitrary level by statistical means not limited to roughness or minimum area, and more particularly to a synthetic-aperture radar system for measuring the average height of terrain peaks or ocean waves.

Much work has been done on the remote measurement of ocean wave heights using nonimaging radar, and a review of such techniques is presented by D. E. Barrick, Remote Sensing of Sea State by Radar, Remote Sensing of the Troposphere, edited by V. E. Derr, Ch. 12, U.S. Government Printing Office, Washington, D.C., 1972. G. T. Ruck, D. E. Barrick, and T. Kaliszewski, Bistatic Radar Sea State Monitoring, Battelle Columbus Laboratories, Columbus, Ohio, Rep. 1388-19, 1971, have analyzed the possibility of using two-frequency radar interferometry for measuring wave slopes, and Weissman, et al. have demonstrated the operation of a radar system utilizing this technique both in laboratory wave tank measurements and flights over the Caribbean, [D. E. Weissman, Two Frequency Radar Interferometry Applied to the Measurement of Ocean Wave Height, IEEE Trans. on Ant. and Prop., AP-21, 649=656, 1973; D. E. Weissman, C. T. Swift, W. L. Jones, Jr., J. W. Johnson, W. L. Grantham, J. A. Howell, J. C. Fedors, and J. J. Davis, A Dual Frequency Radar for Ocean Roughness Sampling, URSI Fall Meeting, Boulder, Colorado, 1974.] The possibility of measuring ocean wave heights by studying the return from a short-pulse radar has also been demonstrated by B. D. Zamarayev and A. I. Kalmykov, On the Possibility of Determining the Spatial Structure of an Agitated Ocean Wave by Means of Radar, Izvestia, USSR Academy of Sciences -- Atmospheric and Ocean Physics, 5, 64-66, 1969, and by D. L. Hammone, R. A. Manella and E. J. Walsh, Short Pulse Radar Used to Measure Significant Wave Height and Surface Wind Speed, URSI Fall Meeting, Boulder, Colorado, 1975.

The principal concept involved in the wave height measurement, which is common to the interferometric techniques utilizing radar sensors, is that the resultant intensity for the electric field scattered from the ocean wave at some arbitrary point in space changes as the wavelength or angle of illumination for the electromagnetic radiation varies. This variation is dependent on the height difference between the crest and trough of the wave, and a measurement of the rate of this variation provides an indication of the ocean wave height. This principle has been used before in the design and operation of the tellurometer for the precise determination of range differences, [T. L. Wadley, Electronic Principles of the Tellurometer, Trans. South African Inst. Electr. Engrs., 49, 143-161, 1958.] The various radar techniques for measuring ocean wave heights differ in the handling of the radar data, which is governed by the constraints set by the particular radar instrument being used, the accuracy and spatial resolution of the measurement being related to the bandwidth and the resolution of the radar system involved.

An extensive theoretical and experimental treatment of the utilization of this concept to measure the height of surface features within a resolution cell in images obtained by a tunable dye laser has been performed by Jain, the present inventor [A. Jain, A Wavelength Diversity Technique for Smoothing of Speckle, Ph.D. Thesis, California Institute of Technology, Pasadena, California, 1975], [N. George and A. Jain, "Space and Wavelength Dependence of Speckle Intensity," Appl. Phys., 4, 201-212, 1974], and [George, et al., "Experiments on the Space and Wavelength Dependence of Speckle," Appl. Phys., 1, 157, 1975]. The surface heights obtained by varying the wavelength can be determined quite accurately to the order of a wavelength, the accuracy improving with greater heights.

While the wavelengths of the ocean waves can be determined by inspection of the radar imagery or by obtaining its Fourier transform [W. E. Brown, Jr., C. Elachi, and T. W. Thompson, "Radar Imaging of Ocean Surface Patterns," J. Geophys. Res., Vol. 81, No. 15, pp. 2657 (May 1976); C. Elachi, "Two Dimensional Wavelength Spectrum of Swells Across the North Atlantic," J. Geophys. Res., (in press)], it is also of considerable interest to measure the heights of the ocean waves. A method of processing the images obtained by synthetic aperture radar to determine the heights of ocean waves at localized areas of interest in the images, and experimental results demonstrating this technique are presented.

In accordance with the invention, significant height information is obtained (of peaks in a rough terrain or heights of ocean waves) from a synthetic-aperture radar signal, S(t), by means for compressing the radar signal for intensity images I.sub.1 (x), I.sub.2 (x) . . . I.sub.n (x) over different frequency spacings, .DELTA.f.sub.r0, of the available chirp or Doppler bandwidths, means for cross-correlating these intensity images over a selected area and means for determining the value of a given (e.g., zero) component of this cross-correlation as a function of the spacing .DELTA.f.sub.r0 of the center frequencies, .DELTA.f.sub.r0, for the cross-correlated images as a direct indication of height. In a digital processor, the data compression means is comprised of a plurality of match filters of predetermined width, .DELTA.f.sub.r, each followed by a squarer to provide the intensity images I.sub.1 (x), I.sub.2 (x) . . . I.sub.n (x) which are cross-correlated over selected areas in the x dimension of the image to obtain the values .intg.I.sub.1 (x)I.sub.k (x)dx for k selected from the sequence 1, 2, . . . n. The cross-correlation value for k = 1 yields a normalizing factor that may be used in subsequent cross-correlations for k = 2, 3 . . . n. In an optical processor the data compression means is comprised of a data film on which the radar signal has been recorded, and a tilted plane correlator to produce on a map film the intensity images I.sub.1 (x), I.sub.2 (x) . . . I.sub.n (x) for different frequency spacings, .DELTA.f.sub.r0. The cross-correlating means is then comprised of optical means for superimposing a selected area of an intensity image on a map film at one plane onto a corresponding area selected in the intensity images to be cross-correlated on a duplicate map film in another plane, and illuminating the superimposed areas with coherent light, and means for focusing the light transmitted through the superimposed and selected areas onto an intensity detector on a third plane. In either case, the apparatus thus includes means for compressing synthetic-aperture radar data over different sections of the available chirp or Doppler bandwidths to obtain intensity images, and means for cross-correlating these images to obtain height data. This height data for a predetermined range of frequency spacing, .DELTA.f.sub.r0, is plotted to calibrate the system. Once calibrated, only two intensity images need be cross-correlated at a predetermined spacing to determine height directly.

FIGS. 1(a) and 1(b) show a system diagram for a height processor in which a synthetic-aperture radar signal S(t) is compressed by matched filters f.sub.1 (t,x), f.sub.2 (t,x) . . . f.sub.n (t,x), corresponding to a processing bandwidth of .DELTA.f.sub.r and at chirp or Doppler processing frequencies f.sub.r01, f.sub.r02, . . . f.sub.r0n, respectively, the square of these outputs cross-correlated by the operation x, and the slope of the resulting curve measured to provide the height data.

FIGS. 2(a) and 2(b) show radar imaging geometry useful in understanding the theory of the present invention.

FIG. 3 illustrates the configuration of a tilted plane optical processor, for compressing into a map film the radar signal from a data film, the bandwidth of the chirp or Doppler spectrum being controlled by the size of aperture A of a frequency plane filter F, and the center frequency being controlled by the position of aperture A.

FIG. 4 illustrates the configuration of an optical cross-correlator in which the pupil P limits the area of a map film I.sub.1 at plane I illuminated by a laser L and collimator C, and lenses L.sub.11 and L.sub.22 of focal lengths f(30.5 cm each) image I.sub.1, onto plane II. The lens L.sub.3 of focal length f, Fourier transforms the product I.sub.1 L.sub.2 onto plane III, where the total power is detected by detector D.

FIG. 5 illustrates cross-correlation data curves, R.sup.E (.DELTA.f.sub.r0), for the ocean wave heights of 0.9, 1.1, 1.8, and 4 m as a function of frequency separation for the chirp center frequencies used in the processing using a chirp bandwidth equal to 2 MHz.

Although reference will hereafter be made to the ocean surface in determining its roughness (wave height), it will be apparent that the same techniques may be employed to determine the roughness (peak heights) of land surface. As noted hereinbefore, the use of synthetic aperture imaging radars in oceanography has recently attracted considerable interest (W. E. Brown, Jr., C. Elachi and T. W. Thompson, Oceanographic Observation with Imaging Radar, URSI Fall Meeting, Boulder, Colorado, 1973; L. I. Moskowitz, "The Feasibility of Ocean Current Mapping vis. Synthetic Aperture Radar Methods," Proc. Fall Conv. Amer. Soc. Photogrammetry, 760-771, Oct. 1973; T. R. Larson and J. W. Wright, "Imaging Ocean Current Gradients with Synthetic Aperture Radar," URSI Fall Meeting, Boulder, Colorado, 1974). The radar generates a two-dimensional high resolution image of the ocean surface. The relation between the image brightness distribution and the properties of the surface is still not well understood; however three models (or any combination of them) seem to be the most popular: 1. The image represents the modulation of the small scale surface roughness. This modulation could result from the presence of a swell, internal wave, local wind, slicks, . . . 2. The image represents the modulation of the local surface slope due to the presence of a large wave or ship wake. 3. The image is a result of the modulation of the water velocity by a wave or a current. The last model, if correct, is inherent to synthetic aperture imaging radars, in contrast to real aperture imaging radars, because of the importance of the Doppler information in the generation of the high resolution image (A. W. Rihaczek, Principles of High Resolution Radar, McGraw Hill, N.Y., 1969; J. W. Goodman, Introduction to Fourier Optics, McGraw Hill, San Francisco, 1968; R. O. Harger, Synthetic Aperture Radar Systems, Academic Press, N.Y., 1970). Even though the formation of the image is not yet well understood, radar imagery does show a large variety of ocean surface phenomena which, in most cases, have been verified with ground truth measurements: large gravity waves, internal waves, currents, ship wakes, island wind shadows, surface slicks, weather fronts and others. It should be mentioned that the uniqueness of the synthetic aperture imaging radar is that it provides high resolution imagery (tenths of meters) regardless of weather conditions, cloud coverage, sun illumination and platform altitude (even from orbiting spacecrafts).

The radar image of long ocean waves provides two-dimensional patterns and therefore the wavelength of these waves. But an important parameter of these waves is their height. The system of the present invention shown schematically in the block diagram of FIG. 1 can be used to derive the wave height from the radar image by measuring the rate of change of the speckle pattern in the synthetic aperture radar image as a function of the radar frequency.

A somewhat similar technique has also been used by D. Weisman, "Two Frequency Radar Interferometry Applied to the Measurement of Ocean Wave Height," IEEE Trans. on Ant. and Prop., AP-21, 649-656 (1973) to derive ocean wave height along a flight line using a two-frequency correlation radar technique.

In the art of radar imagery, a speckle pattern is defined as the pattern which results from the interferometric addition of the echoes generated when a coherent electromagnetic wave is scattered from a rough surface. This pattern is familiar to people who observe surfaces using coherent laser light. (This frequency diversity technique has been extensively used in the field of coherent optics to determine surface roughness.) If the surface is smooth there is no speckling. If the surface is slightly rough speckling appears. Then, a change in the electromagnetic wave frequency results in a change of the speckle pattern due to the change in the phase relationships for the electromagnetic field scattered across the rough surface, the magnitude of this change increasing with the height variation on the surface. A detailed statistical analysis gives the resulting value for the cross correlation of the radar images of such a surface at the two different frequencies separated by .DELTA..nu., as ##EQU1## where { } denotes the average value, and I(.nu.) is the intensity distribution in the image formed at frequency .nu., c the velocity of light, h.sub.o the rms height variation on the surface and f a factor depending upon the process used to obtain the imagery and the cross correlation, i.e., film nonlinearities, changes from signal film to image film, image negative to image positive, etc.

To implement the principle of above equation with side-looking synthetic aperture imaging radar 10 shown in FIGS. 2a and 2b, the radar may have a chirped signal with frequency ranging from 1200 MHz to 1210 MHz. The echo is recorded on a data film S which may then be optically processed to generate map films as will be described with reference to FIG. 3 for use in an optical correlator shown in FIG. 4. Alternatively, the echo signal S(t) may be processed digitally to obtain wave height data as shown in FIG. 1.

The principle involved in the ocean wave height measurement using side-looking (or forward-looking) synthetic-aperture imaging radar is that the resultant intensity for the electric field scattered from a rough surface at some arbitrary point in space changes as the wavelength or angle of illumination for the electromagnetic radiation varies. This variation is dependent on the height difference between the crest and trough of the wave, and a measurement of the rate of this variation provides an indication of the ocean wave height (or height of peaks in rough terrain).

The essential elements in the processing of the synthetic-aperture radar data for obtaining significant wave height information are shown in FIG. 1(a). They involve compressing the radar data signal S(t) for ocean images over different sections of the available chirp or Doppler bandwidths f.sub.r01, f.sub.r02, . . . f.sub.r0n, cross-correlating these images and plotting the zero component of this cross-correlation value as a function of the spacings of the center frequencies for the bandwidths. A measure of the value of this curve for a given frequency spacing, .DELTA.f.sub.r, would provide a direct indication of the ocean wave heights. This cross-correlation function, which is the output of the radar wave height processor shown in FIG. 1(a) is equal to the square of the product of the characteristic function of the ocean wave heights and the pupil function used in the chirp spectrum aperture for processing the radar signal. The intermediate calculation of the image produced by the radar system, will take into account the effect of the orbital velocity of the waves when the optical system of FIG. 3 is used with the optical correlator of FIG. 4. A detailed treatment of this effect will be discussed hereinafter. Thus, a knowledge of the chirp spectrum pupil function of the radar, and an experimental determination of the cross-correlation curve by the wave height processor, leads directly to a measurement of the ocean wave height distribution. It has been found that the measured cross-correlation curve is equal to .GAMMA..sub.r (.DELTA.f.sub.r0) exp (-0.12.DELTA.f.sub.r0 .sigma..sub.H), where .sigma..sub.H is the significant wave height, .DELTA.f.sub.r0 the frequency separation for the cross-correlation value, and .GAMMA..sub.r (.DELTA.f.sub.r0) a number depending upon the chirp spectrum aperture function .GAMMA..sub.r for the radar processor.

Referring now to the processor shown schematically in a functional block diagram in FIG. 1(a), a number of functional blocks 1, 2 . . . N receive the chirp return signal S(t) and split it into n frequency channels, each of a predetermined width, .DELTA.f.sub.r, of about 2 MHz. In practice, digital bandwidth filters may be used to provide the n frequency channels for electronic data processing, but in the preferred embodiment to be described with reference to FIGS. 3 and 4, optical filtering techniques are used. Outputs from the functional blocks 1 through N are individually squared in circuits 11, 12 . . . NN to change the electric field signals to intensity signals I.sub.1 (x), I.sub.2 (x) . . . I.sub.n (x). These signals are cross correlated in pairs by a circuit 20 and integrated. A slope detector 21, calibrated in heights, detects the drop-off (width) of the resulting integrated signal to provide as a system output a direct measurement in height. This technique for height determination is based on the fact that the synthetic aperture radar is an imaging instrument, where the resolution is determined by the chirp and doppler bandwidths used for the radar system. Thus, any point on the image contains contributions from the field scattered by an area of the terrain determined by the resolution cell size of the radar. By narrowing the effective doppler or chirp bandwidth utilized by the instrument, the area of the ocean surface contributing to the intensity at a given image point may be widened to any desired magnitude. The total path length that the radar wave has to travel to the image is different for the trough of the wave than for the crest, and the resultant phase difference for the contribution from these two regions of the wave is dependent on the frequency or angle of illumination of the wave. As this frequency or angle of illumination is varied, the fluctuation of the intensity at the image point, which receives contributions from both wave regions, is related to the difference in heights between the crest and trough of the wave, modified by the effect due to the resolution cell size of the radar. By measuring the effective rate of this fluctuation, one can infer the heights (of the ocean waves or of peaks in a rough terrain).

Before processing with a more detailed description of wave (or peak) height determination, the optical system shown in FIGS. 3 and 4 will be described. It has been used with an L-band synthetic aperture imaging radar operating at 25 cm wavelength, with a phased array antenna 75 cm long, 25 cm wide, and 5 cm thick to obtain a range beamwidth of 90.degree., centered 45.degree. off vertical and an azimuth beamwidth of 18.degree., centered on zero doppler. The transmitted peak power is 7 kW, the pulse repetition frequency 1 kHz, the radar platform velocity 250 m/sec, the pulsewidth 1.25 .mu.sec, and the total chirp spectrum bandwidth 10 mHz. The data are recorded optically on a film, shown as a data film 22 in FIG. 3 using a dual-channel recorder, such as the Goodyear 102, dual-channel recorder, and processed optically to convert the signal film into imagery with a resolution of 20 m range (45.degree.) and 10 m along track. Once the data film is developed, it is employed in a tilted plane correlator (TPC) shown schematically in FIG. 3. The data film 22 is illuminated by the output of a HeNe laser (not shown) after it has passed through a spatial filter-collimating lens system. A Fourier transform lens L.sub.1, a frequency plane aperture filter F, and an imaging lens L.sub.2 convert the range data into the range image at a plane on a map film 23. A cylindrical azimuth telescope L.sub.3 has the function of bringing the azimuth data at plane I to an alongtrack image focus on the map film 23. By controlling the position of a data plane aperture 24 and the spatial position of the frequency plane aperture filter F, the chirp frequency and the center bandwidth can be controlled. By setting the value of the frequency plane aperture at a given range bandwidth, the same strip of signal film can be processed for different center frequencies to provide the required images by positioning the location of the frequency plane aperture at the corresponding spatial coordinate in the range dimension. In that manner, the radar return signal recorded on the data film from a synthetic aperture radar in an aircraft 30 (FIG. 2a) is processed by the TCP for different chirp bandwidths f.sub.r01, f.sub.r02 . . . f.sub.rOn of predetermined width (e.g., 2 MHz) to provide on the map film 23 a record of the split radar signal equivalent to that of the different channels 1, 2 . . . n in FIG. 1(a). The map film itself functions as a squarer.

The optical setup used to cross correlate pairs of images of the same area, but obtained for different chirp spectrum center frequencies, is shown in FIG. 4. An argon-ion laser 31, combined with a spatial filter collimator system, C, is used to illuminate the film image I.sub.1 (x, y, f.sub.r0) processed for the particular chirp spectrum frequency f.sub.r0 and the bandwidth f.sub.r. The two Fourier transform lenses L.sub.11 and L.sub.22, each of focal length f, are used to provide unity magnification of an image at plane I onto plane II, where the transparency of transmissivity I.sub.2 (x + .DELTA.x, y + .DELTA.y, f.sub.r0 + .DELTA.f.sub.r0) is situated, which is the image of the same area as I.sub.1 (x, y, f.sub.r0) but processed at the chirp spectrum frequency f.sub.r0 + .DELTA.f.sub.r0 and shifted by (.DELTA.x, .DELTA.y) in the X- and y-directions for the optical system. The distance between L.sub.11 and L.sub.22 is 2f, and between plane I and L.sub.11 and plane II and L.sub.22 is f, so that no phase effects are introduced by the imaging process. The spatial distribution of the electric field at the exit of plane II is proportional to I.sub.1 (x, y, f.sub.r0)I.sub.2 (x + .DELTA.x, y + .DELTA.y, f.sub.r0 + .DELTA.f.sub.r0), and this is Fourier transformed by a lens L.sub.33 of focal length f onto plane III. The power incident on plane III is measured by a photodetector 32 whose aperture is large enough to collect all the Fourier components involved. Using Parseval's theorem, and assuming a unit amplitude illumination at I.sub.1, the power measured by the detector is ##EQU2## where .DELTA.f.sub.r is the bandwidth at which I.sub.1 and I.sub.2 were processed, (.DELTA.x, .DELTA.y) are relative shifts of I.sub.1 with respect to I.sub.2, and f.sub.r0, f.sub.r0 + .DELTA.f.sub.r0 are the two processing frequencies involved. For the purposes of wave height measurement, .DELTA.x and .DELTA.y would be equated to zero and the function W.sub..DELTA.f.sbsb.r (.DELTA.f.sub.r0, 0, 0) plotted as a function of .DELTA.f.sub.r0. In considering I.sub.1 and I.sub.2 to be the two images of the ocean wave, we have assumed that the film nonlinearities of the image film can be ignored and have factored out the gamma of the film. The result is displayed using a CRT 33 or recorder 34, and a slope detector 35 provides a direct indication of ocean wave height after the system has been calibrated.

As noted hereinbefore, as the frequency or angle of illumination is varied, the fluctuation of the intensity of the image point, which receives contributions from both wave regions is related to the difference in heights between the crest and trough of the wave, modified by the effect due to the resolution. A detailed discussion will now follow on the theory by which heights, of ocean waves or peaks of rough terrain, can be inferred from measurements of the effective rate of this fluctuation.

In the linear system idealization for the synthetic aperture radar performance [Radar Handbook, ed. Merrill I. Skolnik, Chap. 23 (by L. J. Cutrona), New York, 1970] the electric field e(x.sub.0, y.sub.0) in the image of the radar output is ##EQU3## wherein .psi.(x.sub.0, y.sub.0 ; x, y) is the generalized ambiguity function as defined by Cutrona [ibid.] for the radar and g(x, y) is the terrain reflectivity.

In the calculation for the ambiguity function, however, it is assumed that the terrain reflectivity does not depend upon the wavelength of illumination or the look angle of the radar. In the case of an uneven surface, this assumption is not applicable, and it is necessary to consider the effects of the wavelength and angular dependence of the radiation scattered by the surface. In order to take this phenomenon into account, calculate the image field for each wavelength component of the chirp spectrum and position of the radar platform and integrate this electric field over the total synthetic antenna length and the chirp spectrum bandwidth available.

The impulse response function for the radar imaging system for a given value of the antenna position, x.sub.A, and the chirp spectrum frequency, f.sub.r, may be expressed as a convolution of the generalized ambiguity function with the Fourier transform of the pupil function for this range and Doppler aperture. This aperture is the Diarc delta function .delta.(f.sub.dd - f.sub.d).delta.(f.sub.rr - f.sub.r), where the phase factor due to he time delay, f.sub.d, is equal to x.sub.A f.sub.r0 /Z, f.sub.r0 being the frequency of operation of the radar and Z the height of the radar platform from the terrain. While the Doppler shift of the return signal is f.sub.d .nu.r/c, for .nu..sub.r the aircraft velocity and c = 3 .times. 10.sup.8 m/sec, it is the quantity f.sub.d that is significant in the matched filter processing of the radar correlator. The Fourier transform of this delta function is defined by ##EQU4## the symbol y.sub.0 being used to specify the particular nature of the chirp processing operation.

The resultant impulse response of the radar system for these particular frequencies (f.sub.d, f.sub.r) is given by the convolution of the impulse response U(x, y.sub.0) with the generalized ambiguity function .psi.(x, y). Thus, for the complex terrain reflectivity g(x, y), the electric field recorded in the radar image is the convolution of g(x, y) with the convolution of .psi.(x, y) with U(x, y.sub.0). Further, for the Doppler spectrum aperture function given by .GAMMA..sub.d (f.sub.d - f.sub.d0), which has an effective width of .DELTA.f.sub.d and is centered at f.sub.d0, and a chirp spectrum aperture function .GAMMA..sub.r (f.sub.r - f.sub.r0) of width .DELTA.f.sub.r centered at the radar frequency f.sub.r0, the electric field in the radar image is the superposition of fields obtained for each Doppler and chirp frequency component weighted by the functions .GAMMA..sub.d (f.sub.d - f.sub.d0), .GAMMA..sub.r (f.sub.r - f.sub.r0). This yields the equation ##EQU5## where * denotes the convolution operation. Eq. (3) describes the field in the image plane of the terrain reflectivity g(x, y) obtained by a synthetic aperture radar defined by the generalized ambiguity function .psi.(x, y) and the processing Doppler and chirp spectrum apertures .GAMMA..sub.d (f.sub.d - f.sub.d0), .GAMMA..sub.r (f.sub.r - f.sub.r0). Since the azimuth aperture for the synthetic aperture radar processor is generally symmetric, it is safe to consider the value of f.sub.d0 to be zero in our further evaluation of Eq. (3).

From the geometry of FIG. 2a the terrain reflectivity function gp(x, y) for a particular area surrounding the point P located at (0, y.sub.p) is: ##EQU6## R(x, y) is the amplitude reflectivity of the terrain due to absorption and dielectric constant variations and h(x, y) is the profile of the terrain resting on the local slopes of angles (.theta..sub.Ptx, .theta..sub.Pty) shown in FIG. 2b for the radar platform and terrain of FIG. 2a which vary slowly compared to the variation of h(x, y). The slope angles .theta..sub.Ptx and .theta..sub.Pty are a function of both the slope angles (.theta..sub.Px, .theta..sub.Py) of the side-looking radar and the tilt (t) of the terrain which, for an ocean wave, is the angle between a tangent to the ocean wave and a horizontal X-Y plane. FIG. 2b illustrates the geometry of the angle .theta..sub.Pty measured in the Y-Z plane and FIG. 2c illustrates the geometry of the angle .theta..sub.Ptx measured in the X-Z plane. Note that in FIG. 2a there is a slight tilt (.theta..sub.tx, .theta..sub.ty) of the terrain with respect to the X and Y axes of the horizontal X-Y plane, to illustrate a typical case in which the roughness of the terrain is the major factor in the tilt angles .theta..sub.Pty, .theta..sub.Ptx. The distance of point P from the radar, .sqroot.Z.sup.2 +y.sub.p.sup.2, is considered to be large, so that the radar wave incident on the area surrounding P is essentially plane and is at angles (.theta..sub.Px, .theta..sub.Py) with respect to the normal to the horizontal plane. The function exp[(i4.pi./c)f.sub.r cos.theta..sub.Pxy h(x, y)] is the phase modulation in the reflected wave due to the small-scale roughness h(x, y). The factor ##EQU7## is the phase modulation in the reflected wave from the tilt variations of the ocean wave and the angles of incidence .theta..sub.Px and .theta..sub.Py of the radar wave and is determined from equating the angle of incidence to the angle of reflection from the surface. The angle .theta..sub.Px is generally small, and for a typical spacecraft radar it is less than 0.1 radian. (.theta..sub.Ptx, .theta..sub.Pty) are the angles of tilt of an ocean wave when g.sub.p (x, y) is an ocean surface and are generally limited to be less than 0.143. Thus, in further computations involving g.sub.P (x, y) it is reasonable to consider sin(.theta..sub.Px - 2.theta..sub.Ptx) as approximately equal to (sin.theta..sub.Px - sin2.theta..sub.Ptx). Since sin.theta..sub.Px is equal to x.sub.A /Z, f.sub.r0 sin.theta..sub.Px is the Doppler frequency f.sub.d, and .theta..sub.Py is generally much larger than .theta..sub.Pty for side-looking radars.

To evaluate the radar performance for a given reflectivity, substitute Eq. (4) into Eq. (3). For h(x, y) = .theta..sub.Ptx = .theta..sub.Pty = 0, which corresponds to an even terrain of reflectivity R(x, y), the image electric field is found to be ##EQU8## Eq. (5) neglects the variation of x sin .theta..sub.Px in comparison to y sin .theta..sub.Py, which is a much larger term, in the function .gamma..sub.r (y). For infinite processing bandwidths, given by the widths of .GAMMA..sub.r (f.sub.r), .GAMMA..sub.d (f.sub.d), the functions .GAMMA..sub.r (x), .gamma..sub.d (y) reduce to delta functions, and Eq. (5) reduces to Eq. (1).

On the other hand, for a radar with an infinitely high resolution, when .psi.(x, y) would be a delta function, the image electric field would be the convolution of the reflectivity with the impulse response functions .GAMMA..sub.d and .GAMMA..sub.r, and the image resolution would be determined primarily by the bandwidths utilized by the radar processor. In general, the resolution of the radar is indeed limited by the data handling capabilities of the processor, and so it is safe to assume, for further analysis, that the generalized ambiguity function is the delta function.

Consider the case of an uneven surface, such as ocean waves, and assume that the function R(x, y) is unity in Eq. (4). We describe h(x, y) as a series of steps whose perimeters are defined by the functions S.sub.ij (x - x.sub.i, y - y.sub.j) centered at the various points (x.sub.i, y.sub.j), and the profile of the surface within this (ij).sup.th step is defined by the function h.sub.ij. The effective length of the S.sub.ij function is given by l.sub.i in the x-direction and l.sub.j in the y-direction. Thus, the surface profile h(x, y) can be written as ##EQU9## For nonoverlapping S.sub.ij, the exponential exp[(-i4.pi./c) cos .theta..sub.Pxy h(x, y)] can be rewritten and the resulting surface reflectivity can be written in the following form: ##EQU10## The first term in Eq. (6) represents the field reflected if the terrain is perfectly flat with only tilt variations. The term identified by the summation consists of two terms, the first being the component of the reflected field due to the surface roughness, and the second being the component that must be subtracted from the field for the flat surface case to provide the reflected field for the smooth component of the terrain.

To obtain the electric field in the image of the field given by Eq. (6), substitute Eq. (6) in Eq. (3), neglect the effect of the variation of f.sub.d.sup.2 compared to f.sub.d in the term cos.theta..sub.Pxy) and assume that the scattering cross-section of the terrain, determined by the Fourier transforms of S.sub.ij (x - x.sub.i, y - y.sub.j) does not vary significantly with the Doppler chirp spectrum frequency variations over the bands used in the radar processor, as compared to the variations of the exponential phase terms containing this (f.sub.d, f.sub.r) dependence. The resultant image electric field is ##EQU11## where .beta..sub.ij is the scattering cross-section factor and is equal to ##EQU12## for s.sub.ij defined by ##EQU13## In the approximation, where S.sub.ij (x, y) are given by the step functions [rect(x/l.sub.j)], the cross-section factor s.sub.ij is equal to ##EQU14## The functions .gamma. .sub.d (2x/c), .gamma. .sub.r [y sin(.theta..sub.Py - 2.theta..sub.Pty)] are the same as those defined in Eq. (5) and, for the chirp and Doppler spectrum pupil functions, approximated by ##EQU15## these functions become ##EQU16## The product .gamma..sub.d .gamma..sub.r on the left-hand side of the convolution sign in Eq. (7) corresponds to the impulse response of the radar, and the resolution cell size is determined to be 2c/[.DELTA.f.sub.r sin(.theta..sub.Py - 2.theta..sub.Pty)] in the range dimension and c/.DELTA.f.sub.d in the azimuth dimension, this being equal to the size of the radar antenna if the full Doppler aperture is used. The function F.sub.ij (x, y; x.sub.i, y.sub.j ; h.sub.ij) describes the effect on the image field due to the phase variations in the scattered field arising from the unevenness of the terrain and is given by ##EQU17## where cos .theta..sub.Py0 is equal to the value of cos .theta..sub.Pyx at the Doppler frequency f.sub.d equal to zero.

The scattering factor .beta. is dependent on the tilts of the ocean waves (.theta..sub.Ptx, .theta..sub.Pty) and also the distance away from nadir as determined by .theta..sub.Py and the size of the small scale roughness (l.sub.i, l.sub.j). A modulation in the wave tilt or roughness size would contribute to a modulation of the intensity in the image. A detailed analysis and experimental evaluation of the effect of tilt and roughness of the terrain on the image modulations for the radar system have been reported by Schaber, et al. [G. G. Schaber, G. L. Berlin, and W. E. Brown, Jr., Geol. Soc. Amer. Bull., March 1975] for the case of desert images and by Elachi and Brown [C. Elachi and W. E. Brown, Jr., IEEE Trans. on Ant. and Prop., in press] for the case of ocean wave imaging.

While in principle, the synthetic antenna length is determined by the area illuminated by the radar, the actual length is modified by the finite angle over which the radiation is backscattered by the small-scale roughness. Since this angle is determined by the spatial lengths of these scatterers, a modulation in the size of the scatterers would also result in a modulation in the azimuth resolution of the radar if the scatterer sizes were larger than the size of the radar antenna.

Note that Eq. (7) does not include the effect of motions of the surface being imaged. Since f.sub.d is the time variable in Eq. (3), the effect of surface motions can be included by making g(x, y) a function of f.sub.d as well. A uniform motion of g(x, y) in the y-direction with a velocity v.sub.g can be included by substituting the function g(x, y - .xi.f.sub.d) for g(x, y) in Eq. (3). .xi. is equal to Zv.sub.g /f.sub.r0 v.sub.r.sup.2, where v.sub.r is the velocity of the radar platform. A resulting evaluation provides no change in the image except for a uniform shift proportional to .xi. in the azimuth direction. A uniform motion in the x-direction is equivalent to a change in the equivalent velocity of the aircraft and thus a change in the focusing parameter for radar imagery.

At any one point on the ocean surface, the height varies according to the orbital velocity, and the cross-section of this point according to the ocean wave length. A patch of constant cross-section moves in the direction of the ocean wave propagation at the ocean wave velocity. Since the radar maps the time delay of the backscattered signal from a patch of constant cross-section, the ocean wave appears to it to be moving in a plane parallel to the radar platform motion and the orbital velocity effects are not important. However, an object of constant cross-section, such as a ship, resting at a particular point on the ocean, and moving vertically at the orbital velocity, will have a modified image and its image will generally be degraded in the azimuth direction. It should be pointed out that the synthetic aperture radar is sensitive to the variation of time delays of the signal scattered by a target, and the actual effect of the instantaneous shift in the frequency of the backscattered signal due to the orbital motion is included in this time delay variation which depends on the variation of position of the scattering patch. In such a consideration, the variation of the velocity and radar cross-section of this patch must also be taken into account, and approximating these quantities by a constant number during the integration time of the radar processor has in the past led to considerable confusion in the interpretation of the synthetic aperture radar ocean wave imagery. A dependence of the ocean wave tilt angle .theta..sub.Ptx on time would lead to a change in the resolution cell size of the radar, the result being a smaller resolution cell size if the tilt angle is increasing and a larger cell size if it is decreasing. This effect, however, would be very small for the normally occurring rates of change of the tilt. The random motion of the small-scale roughness would have negligible overall effect on the radar image since its main contribution is in the backscatter modulation function .beta..sub.ij, a treatment of this negligible effect has been considered by Elachi and Evans [C. Elachi and D. E. Evans, IEEE Trans. on Ant. and Prop., in press].

The ability to determine height of ocean waves, or peaks of a rough terrain, according to the invention derives from Eq. (7). For a particular point P.sub.1 in the area of interest, the image electric field may be rewritten as ##EQU18## where the finite size of .gamma..sub.r (y), .gamma..sub.d (x) limits the number of scatterers that contribute to the field at an image point. In this case, this number is given by the summation over the N scatterers; .alpha. is the uniform background level term and is the first term on the right-hand side of Eq. (7); .beta. is the scattering factor; and (.theta..sub.Pty.sbsb.j, .theta..sub.Ptx.sbsb.i) are the tilt angles at points (x.sub.i, y.sub.j) which have been assumed to be slowly varying compared to the small-scale roughness h(x, y). Assuming that the scatterers are densely packed [L. I. Goldfischer, J. Opt. Soc. Am., 55, 247-253, 1965], N is equal to the resolution cell area divided by the area occupied by a scatterer, and this quantity can be approximated by 2c.sup.2 /(l.sub.i l.sub.j .DELTA.f.sub.r .DELTA.f.sub.d sin .theta..sub.Py) (neglecting the effect on the resolution cell size due to the small tilt angles, the scatterer size modulations, and the heights h.sub.ij).

The small-scale roughness is usually less than the wavelength of the radar and the variation of exp[(i4.pi./c)f.sub.r0 cos .theta..sub.Py0 h.sub.ij ] with respect to f.sub.r0 is negligibly small over the chirp spectrum bandwidth available. Also, since .theta..sub.Pty.sbsb.j is a small number, it is permissible to approximate sin(.theta..sub.Py - 2.theta..sub.Pty.sbsb.j) by (sin .theta..sub.Py - cos .theta..sub.Py sin 2.theta..sub.Pty.sbsb.j). The term exp[(i2.pi./c)f.sub.r0 y.sub.j sin .theta..sub.Py ] represents the variation of the phase of the radar wave incident at angle .theta..sub.Py across the resolution cell size of the radar, and with respect to change in radar frequency, occurs in an interval equal to the chirp bandwidth of the radar processor. The factor ##EQU19## corresponds to the phase variation due to the variations in the tilt angles of the terrain. If the resolution cell is small enough so that the tilt of the waves does not change appreciably across its length, the decorrelation length due to the tilt in the range and azimuth directions is determined by .DELTA.f.sub.r /(2 sin.theta..sub.Pty) .times. (sin.theta..sub.Py /cos.theta..sub.Py) and .DELTA.f.sub.d /(4 sin.theta..sub.Ptx), respectively. The composite decorrelation length would depend upon the sum of the slopes of the individual decorrelation curves for the resolution cell size and the tilts approximated by ##STR1##

The resolution cell size of the radar may be large enough to cover more than half an ocean wave, in which case the tilt angles (.theta..sub.Ptx, .theta..sub.Pty) cannot be assumed to be constant within a resolution cell of the radar, but it is possible to rewrite (y.sub.j cos .theta..sub.Py sin.theta..sub.Pty.sbsb.j + x.sub.i sin2.theta..sub.Ptx.sbsb.i) as 2H.sub.ij, corresponding to the height of the ocean wave. Due to the cos .theta..sub.Py term, the dependence of the resolution cell size on sin .theta..sub.Py and the tilt angle variations, the H.sub.ij as defined is not strictly the ocean wave height parameter, but for the purposes of the statistical analysis, it may be considered to be equivalent. Then calculate the cross-correlation for the image intensities processed for the two frequencies f.sub.r0 and (f.sub.r0 + .DELTA.f.sub.r0) and obtain the following result ##EQU20## The brackets < > denote the mean value operation. .beta. is the average value of .beta..sub.ij within the resolution cell and is the backscattering cross-section of the area of interest. N, the number of scatterers in a resolution cell, is approximately equal to 2c.sup.2 /(.DELTA.f.sub.r .DELTA.f.sub.d sin.theta..sub.Py l.sub.i l.sub.j), where l.sub.i l.sub.j are the sizes of the small-scale scatterers. .sigma..sub.H is the standard deviation of the ocean wave height function H.sub.ij, and .sigma..sub.h is the standard deviation of the heights of the scattering waves. .eta..sub.H and .eta..sub.h are the characteristic functions for the probability density of the wave height H.sub.ij and the roughness h.sub.ij and, for a Gaussian distribution, are equal to exp{-1/2[4.pi.(.DELTA.f.sub.r0 /c) .sigma..sub.H ].sup.2 } and exp{-1/2[4.pi.(.DELTA.f.sub.r0 /c) .sigma..sub.h ].sup.2 }, respectively. For the special case when the resolution cell size is smaller than an ocean wavelength, the characteristic function for the wave height .theta..sub.H in Eq. (9a) may be replaced by the characteristic function for the wave slopes equal to

for .DELTA.x, .DELTA.y, the resolution cell sizes. For a Gaussian density function for the slopes, this function is equal to ##EQU21## where .sigma..sub..theta.ty, .sigma..sub..theta.tx are the standard deviations of the wave slopes in the y- and x-directions, respectively.

.alpha., which corresponds to the specular reflection term, is negligibly small compared to .beta. for the non-nadir angles and may be neglected in Eq. (9a). Further, N is a very large number, and the terms of order N can be neglected in comparison to the terms containing N.sup.2. .sigma..sub.h is a small number compared to .sigma..sub.H, and .eta..sub.h [2.DELTA.f.sub.r0 /c).sigma..sub.h ] may be approximated by unity for the values of .DELTA.f.sub.r0 involved with the radar chirp. The quantity ##EQU22## is defined by R(.DELTA.f.sub.r0) as the normalized cross-correlation function for the radar image. This quantity is simply a product of the square of the radar processor pupil function, .GAMMA..sub.r.sup.2 (.DELTA.f.sub.r0), and the square of the characteristic function of the wave heights or the characteristic function for the wave slopes. This may be written as ##EQU23## for the radar resolution cell much larger than an ocean wavelength, and ##EQU24## for the radar resolution cell smaller than an ocean wavelength. Thus, an experimental determination of the cross-correlation function along with a prior knowledge of the pupil function in range for the radar processor would lead to a determination of the characteristic function for the ocean wave heights from which the wave height spectra of the ocean waves may be determined. When the radar resolution cell is larger than the ocean wavelength, Eq. (9b) can be used to interpret the radar data to obtain information on the distribution of wave heights. When the radar resolution is smaller than an ocean wavelength, Eq. (9c) provides a distribution of wave slopes, and knowledge of the radar resolution cell size can be used to convert this information into wave height information. Also, note by examining the rate of change of R(.DELTA.f.sub.r0) as a function of .sigma..sub.H, that the sensitivity of the measurement of .sigma..sub.H from R(.DELTA.f.sub.r0) increases with .sigma..sub.H. The pupil fraction .GAMMA..sub.r (.DELTA.f.sub.r0) in Eq. (9a) has the effective width (.DELTA.f.sub.r /.DELTA.f.sub.r0). For .DELTA.f.sub.r0 larger than .DELTA.f.sub.r /2, the cross-correlation function becomes very small, and for a measurement of the characteristic function .eta..sub.H, it is necessary to measure R(.DELTA.f.sub.r0) for .DELTA.f.sub.r0 <.DELTA.f.sub.r /2. While the theoretical form of the pupil function .GAMMA..sub.r (f.sub.r) is a step function of bandwidth .DELTA.f.sub.r, practical radar systems do not realize this functional form because of the antenna pattern and nonlinearities in the data recording and processing operations.

The result in Eq. (9a) could have been obtained directly from simple physical considerations by adding up the phases of the electric field reflected by the terrain within a resolution cell area of the terrain. The pupil function factor .GAMMA..sub.r (.DELTA.f.sub.r0) is a direct consequence of the nonvertical angle of incidence of the radar radiation, and the characteristic function factor .eta..sub.H is a consequence of the phase variations for reflection from a rough surface.

In order to obtain an experimental evaluation of <I.sub.P.sbsb.1 (f.sub.r0)I.sub.P.sbsb.1 (f.sub.r0 + .DELTA.f.sub.r0)> in Eq. (9a) for a given ocean area, the optical processor of FIG. 3 has been used to obtain the successive images I(x, y, f.sub.r0), I(x, y, f.sub.r0 + .DELTA.f.sub.r01), . . . I(x, y, f.sub.r0 + .DELTA.f.sub.r0n) at the frequencies f.sub.r0, f.sub.r0 + .DELTA.f.sub.r01, . . . f.sub.r0 + .DELTA.f.sub.r0n but for the fixed processing bandwidth of .DELTA.f.sub.r. The image I(x, y, f.sub.r0) in the map film 23 is placed in plane I of the optical cross-correlator in FIG. 4, and using the time codes in these images as a reference for accurate superposition of the image I(x, y, f.sub.r0) on I(x, y, f.sub.r0 + .DELTA.f.sub.r0k) at plane II, the quantity W.sub..DELTA.f.sbsb.r (.DELTA. f.sub.r0k, 0, 0) is measured for each successive I(x, y, f.sub.r0 + .DELTA.f.sub.r0k). In order to account for the effects of laser power fluctuations and background exposures on the film, a measurement of W.sub..DELTA.f.sbsb.r (.DELTA.f.sub.r0k) is also made for large .DELTA.x, .DELTA.y, and this is subtracted from W.sub..DELTA.f.sbsb.r (.DELTA.f.sub.r0k, 0, 0). To normalize the effects of average intensity variations for different I(x, y, f.sub.r0 + .DELTA.f.sub.r0k) further, the difference

is divided by the correlation number when there is no frequency shift, i.e., by W.sub..DELTA.f.sbsb.r (0, 0, 0) - W.sub..DELTA.f.sbsb.r (0, .DELTA.x, .DELTA.y). The resulting curve R.sup.E (.DELTA.f.sub.r0), which is defined by ##EQU25## is plotted as a function of .DELTA.f.sub.r0 and corresponds to the normalized value of .vertline.<I.sub.P.sbsb.1 (f.sub.r0)I.sub.P.sbsb.1 (fr.sub.r0 + .DELTA.f.sub.r0)>.vertline..sup.2, the area on the image correlated being controlled by the pupil P. The quantity R.sup.E (.DELTA.f.sub.r0) measured by the cross-correlator is equivalent to R(.DELTA.f.sub.r0) defined in Equation (9a), the superscript E being used to differentiate between the experimental and the theoretical functions involved.

FIG. 5 is a plot of the function R.sup.E (.DELTA.f.sub.r0) from experimental data where the bandwidth used in processing the successive images was 2 MHz, the processing chirp frequency interval for each successive pair being 0.5 MHz. Since the curves are a product of the characteristic function for the ocean wave heights and the radar pupil function, and since the radar pupil function falls quite drastically for the frequency separations greater than half the chirp bandwidth for the radar processor, it is necessary to confine the examination of the experimental curves for .DELTA.f.sub.r0 < 1 MHz to obtain the wave height information. For a given ocean area, the curves were obtained by performing cross-correlations over different portions of that area, and the deviation of the values of R.sup.E (.DELTA.f.sub.r0) for a given .DELTA.f.sub.r0 was approximately .+-.0.05. As expected, no effect was noticed for the variation of R.sup.E (.DELTA.f.sub.r0) for different positions along the y-axis since the radar aperture function is independent of the look angle.

An examination of R.sup.E (0.5 MHz) shows a variation in the wavelength decorrelation with respect to the ocean wave heights involved, and the points appear to fit the curve 0.65 exp(-0.066.sigma..sub.H), where .sigma..sub.H is the significant wave height. The coefficient 0.066 in the exponential is due to the decorrelation factor from phase changes from crest to trough of a wave, and is close to the theoretical value of [16.pi.(.DELTA.f.sub.r0 /c)]. The factor of 0.65 is the value of the radar chirp spectrum function at 0.5 MHz.

The dependence of the decorrelation function R.sup.E (.DELTA.f.sub.r0) on the chirp spectrum bandwidth used was verified by obtaining the curves R.sup.E (.DELTA.f.sub.r0) for 6-, 4-, 2-, 1-, and 0.5-MHz bandwidths, and extensive decorrelation was noted for the frequency separations .DELTA.f.sub.r0 equal to half the bandwidth employed in the processing. For an infinite-bandwidth radar, the effect on the decorrelation should be nonexistent. By extrapolating the R.sup.E (0.5 MHz) for an infinite bandwidth and comparing this result for R.sup.E (0.5 MHz) at 2 MHz, the value of 0.65 for the chirp aperture effect was verified.

When the cross-correlation function, R.sup.E (.DELTA.f.sub.r0), is evaluated for the frequency separation of .DELTA.f.sub.r0 = 1 MHz for the different wave heights, the functional value obtained is 0.29 exp (0.12.sigma..sub.H). The doubling of the factor in the exponent is due to the doubling of .DELTA.f.sub.r0. The effect of the chirp spectrum aperture function is given by the factor 0.29 for this case.

The chirp aperture function is the Fourier transform of the radar resolution in range. This resolution can be readily measured in the optical processor of FIG. 4 by plotting R.sup.E (O) for different position separations of the two films being correlated, and an experimental examination of this resolution function shows that the width of the resolution cell is inversely proportional to the chirp bandwidth used.

In these results, attention has been focused on measuring the effect on the radar images due to different effective frequencies of the radar operation. A similar calculation and experimental procedure can be employed to exploit the effect on the radar images of different look angles for the radar, and this can be done by processing the radar data for different Doppler bandwidths and cross-correlating these images as a function of the center Doppler frequency separations for these bandwidths.

To summarize the operation of the digital processor shown in FIG. 1(a) in accordance with the system theory discussed in detail above, a synthetic-aperture radar return signal S(t) of bandwidth f.sub.r is separated into a plurality of frequency channels by functional blocks (match filters 1, 2 . . . N, each channel of predetermined width, .DELTA.f.sub.r, over the entire available bandwidth as shown in FIG. 1(b) for a chirp or Doppler spectrum. The chirp spectrum is used for images to one side of the radar track, and the Doppler spectrum is used for images along the radar track. The outputs of the match filters are electric field images which are squared by the squarers 11, 12 . . . NN to obtain intensity images I.sub.1 (x), I.sub.2 (x) . . . I.sub.n (x). These images are degraded since only a part of the chirp or Doppler spectrum is used for each, but height information is the objective, not a high resolution image.

These intensity images are cross-correlated by the circuit 20, but first the image I.sub.1 (x) is auto-correlated (i.e., cross-correlated with itself) to obtain a normalizing factor to be applied to subsequent cross-correlations. Assuming a normalizing factor of unity for simplicity of description, the image I.sub.1 (x) is then cross-correlated with each of the images I.sub.2 (x) through I.sub.n (x) over a selected area. These cross-correlation functions are calculated and fixed (e.g., zero) components of these cross-correlation values are plotted as a function of the spacings of the center frequencies for the channels to obtain the characteristic functions of the ocean wave heights as shown in FIG. 5. These plotted curves calibrate the processor for subsequent wave height measurements. Thus a measure of the value of these curves for a given frequency spacing, e.g. .DELTA.f.sub.r0 = 0.5, provides a direct indication of the ocean wave heights as a function of the slope of the curves. Consequently, the output of the slope detector 21 is a direct indication of ocean wave height. In other words, after once calibrating the system from data plotted as in FIG. 5, the system will provide wave height data directly from the ratio of the cross-correlation function, R.sup.E (.DELTA.f.sub.r0), to the frequency spacing, .DELTA.f.sub.r0.

The functions of the digital processor shown in FIG. 1(a) are implemented optically in an analogous manner. The optical processor shown in FIG. 3 performs the functions of the match filters 1, 2 . . . N of FIG. 1(a), and the map film, which records the electric field images thus derived from the radar signal recorded on a data film 22, inherently squares the electric field images thus performing the functions of the squarers 11, 12 . . . NN of FIG. 1(a). Two copies are made of one of the intensity images thus recorded, such as for channel f.sub.r01. One is placed in plane I in the optical processor shown in FIG. 4 for cross-correlation with the other intensity images recorded on the other map film placed at plane II of FIG. 4. The pupil P selects the area of the intensity image in plane I to be cross-correlated. The cross-correlation functions at plane III are detected by a detector 32 and displayed on a cathode ray tube 33 and/or recorder 34. Once the system has been calibrated, the output of the detector 32 provides a direct indication of wave height just as in the case of the digital processor of FIG. 1(a). There is thus shown an optical processor that is directly analogous to the digital processor.

Although particular embodiments of the invention have been described and illustrated herein, it is recognized that modifications and variations may readily occur to those skilled in the art. It is therefore intended that the claims be interpreted to cover such modifications and variations.