USH1181H

USH1181H - Method for high resolution radar imagery and accurate dimensional measurements

Info

Publication number: USH1181H
Application number: US07/120,216
Authority: US
Inventors: August W. Rihaczek
Original assignee: US Air Force
Current assignee: US Air Force
Priority date: 1985-05-16
Filing date: 1987-11-06
Publication date: 1993-05-04

Abstract

This invention involves a method of processing radar returns to form two-dimensional images of targets such as ground vehicles, aircraft, ships, and so forth. Resolution in one dimension is provided by range resolution, and resolution in the other dimension is provided by Doppler resolution. Man-made targets such as ground vehicles or ships consist of scattering centers which are extended smooth plates, long edges, rods and the like. These scattering units typically are so large compared with the radar wavelength that they have a highly lobed backscattering pattern. When the processing time is increased in order to achieve higher Doppler resolution, it is possible to observe over the processing time a full or even more than one sidelobe of the backscattering pattern. Doppler resolution theory shows that under these circumstances the responses of scatterers are smeared out into meaningless background interference, so that image detail is lost. The present invention avoids this problem by analyzing the signal in segments much shorter than used with conventional processing. This is done despite the fact that the required image detail appears to demand long processing times. Processing a scatterer only over the time interval in which it is dominant (it governs the behavior of the phase then) prevents the smearing of scatterers and there is no need to resolve one from another. The result is an improved radar image of actual targets.

Description

RIGHTS OF THE GOVERNMENT

The invention described herein may be manufactured and used by or for the Government of the United States for all governmental purposes without the payment of any royalty.

This application is a continuation-in-part of application Ser. No. 06/754,895, filed May 16, 1985, now abandoned.

BACKGROUND OF THE INVENTION

This invention relates to a processing method for producing high resolution radar imagery and accurate dimensional measurements using synthetic or inverse synthetic aperature radar. The data is from an airborne radar using a "stretch" format, that is, a long linear FM waveform is used for each individual radar pulse, and the returned waveforms are then mixed on reception with a reference chirp.

In many radar application it is necessary to form two-dimensional images of targets such as ground vehicles, aircraft, ships, and so forth. Resolution in one dimension is provided by range resolution, and resolution in the other dimension is provided by Doppler resolution. The principles are widely applicable and widely applied. Synthetic aperture radar forms two-dimensional images in range and cross range on this basis by utilizing the motion of the platform for Doppler resolution. Inverse synthetic aperture radar is accomplishing the same objective by utilizing the motion of the target to be imaged. In other applications the same principle of Doppler resolution is used, even though no specific name has been given to the process. The very same principles also if, instead of forming images, the processor uses range and Doppler resolution to resolve specific scatterers on the target, and then measures the separation of these scatterers in order to obtain target dimensions.

A serious problem appears with Doppler resolution for man-made targets. Underlying the principles of Doppler resolution is the assumption that a target can be modeled by a set of fixed point scatterers. The emphasis here is placed on "point" scatterers, which term implies that the backscattering behavior is not aspect-angle dependent. In other words, the backscattering is isotropic. In this case as the processing time is increased in order to achieve better Doppler resolution, it is indeed possible to resolve scatterers which are more closely spaced, thus obtaining an image with more detail. Unfortunately, however, man-made targets do not correspond to this model, and Doppler resolution is not working as desired.

Man-made targets such as ground vehicles or ships consist of scattering centers which are extended smooth plates, long edges rods, and the like. These scattering units typically are so large compared with the radar wavelength that they have a highly lobed backscattering pattern. In low-resolution applications coherent processing will extend over relatively small changes in the aspect angle of the target, so that the radar may stay within one and the same sidelobe of the backscattering pattern. However, when the processing time is increased in order to achieve higher Doppler resolution, there will be observed over the processing time a full or even more than one sidelobe of the backscattering pattern. Doppler resolution theory shows that under these circumstances the responses of scatterers are smeared out into meaningless background interference, so that image detail is lost. It is a serious effect which drastically degrades image quality. It is equally detrimental for dimensional measurements, since most of the scatterers will not be observable at all.

As an illustration of the problem, in FIG. 1 shows a two-dimensional image of a tank generated by conventional processing methods. The positions of responses, which then were plotted by the computer as circles. The larger circles mark the positions of stronger responses, whereas the small circles indicate responses so weak that they are more likely to be background interference. A diagram of the tank is overlaid. Note that the stronger responses form a type of L, which is characteristic for conventional images of stationary tanks. Evidently, much detail is lost. It takes a trained operator to recognize that this is the image of a stationary tank, and even for such an operator it is impossible to recognize what tank it is. Conventional imaging does not work at all when the vehicles are moving, so that no sample "image" is shown there.

U.S. patents of interest include Tricoles et al U.S. Pat. No. 4,385,301, which discloses a system providing an image of the location of emitters of electromagnetic radiation behind enemy lines on a battlefield. The patented system comprises an array of antennas for receiving electromagnetic radiation from emitters, and for providing a received signal from each emitter in response to the received radiation; a receiver system coupled to each antenna of the array for measuring the phase and intensity of each received signal, and for providing separate coherent phase signals and amplitude signals that respectively indicate the measured phase and intensity; and a signal processor coupled to the receiver system for processing the coherent phase signals and amplitude signals to provide an image signal for generating an image of the emitters. U.S. Pat. No. 4,387,373 to Longuemare Jr. in column 4, lines 3-6 speaks of improving the quality of a synthetic monopulse radar image by increasing the processing rate. A high speed ambiguity function evaluation processor is taught by Tamukra in U.S. Pat. No. 4,389,373. In U.S. Pat. No. 4,209,853, Hyatt combines a high resolution narrow field-of-view array with a low resolution wide field-of-view array in a holographic processor. In U.S. Pat. No. 4,068,234, O'Meara processes imaging information with a computer programmed for inverse Fourier transform computations, and in U.S. Pat. No. 4,204,262 Fitelson et al perform optical signal processing with a direct electronic Fourier transform.

SUMMARY OF THE INVENTION

An object of the invention is to provide a radar data processing technique to measure accurately the positions of scattering centers on a target, and from the configurations of there scattering centers to obtain target dimensions and shape.

REPORT

The invention is disclosed in a technical report AFWAL-TR-81-1151, titled "Ground Moving Target Identification Via Target Motion Resolution" for Air Force contract F33615-80-C-1065. The report is part of the Defense Technical Information Center collection of documents with the accession number AD C027234. A microfiche copy is attached hereto as an appendix, and is hereby incorporated by reference.

BRIEF DESCRIPTION OF THE DRAWING

FIG. 1 is a diagram showing an image of a stationary tank with conventional radar signal processing;

FIG. 2 is a graph giving an example of the amplitude phase curve for a stationary tank;

FIG. 3 is a diagram showing the image of a stationary tank with the new processing method;

FIG. 4 is a graph showing a doppler plot for gate 4.877 m (16.0 ft), outside the range interval of the tank;

FIG. 5 is a graph showing a doppler plot for gate 5.18 m, first gate with a tank return;

FIG. 6 is a graph showing a doppler plot for gate 6.4 m, (tank);

FIG. 7 is a graph showing a doppler plot for gate 9.753 m, (tank);

FIG. 8 is a graph showing a doppler plot for gate 10.97 m, (tank);

FIG. 9 is a graph showing amplitude and phase for gate 5.18 m, (tank);

FIG. 10 is a graph showing amplitude and phase for gate 6.4 m, (tank);

FIG. 11 is a graph showing amplitude and phase for gate 9.753 m, (tank);

FIG. 12 is a graph showing range rate for response in gate 4.877 m, (tank);

FIG. 13 is a graph showing range rate for response in gate 5.18 m, (tank);

FIG. 14 is a graph showing range rate for response in gate 6.4 m, (tank); and

FIG. 15 is a graph showing range rate for response in gate 9.753 m, (tank).

DETAILED DESCRIPTION

As stated above, it is the highly lobed backscattering pattern of the scattering centers on man-made targets which causes smearing and loss of image detail. This sidelobe pattern implies that when the radar views several scattering centers in one range cell, their backscattering at any one time will be very different. For example, at a given time the radar may see the peak of the backscattering sidelobe from one scatterer, but the null (no backscattering) from another. This means that over some observation time scatterers will become stronger and weaker, alternately effectively appearing and disappearing. Over the processing time needed to form a high-resolution image, many such cycles may occur.

The problem is illustrated by FIG. 2, which shows the amplitude and phase functions for the return from a single range gate for the tank data which led to FIG. 1. Of interest is the phase function. If the range cell contained a single scatterer with isotropic behavior, the phase would increase or decrease linearly, the slope of the phase function giving the Doppler of the scatterer. If the range cell contained two such scatterers, the phase function would again be linearly increasing or decreasing, as for a single scatterer, but a roughly sinusoidal modulation would be impressed by the second scatterer.

Now consider the actual example of FIG. 2 Ignoring the high-frequency noisy variations and considering only the trend indicated by the lines drawn in the figure, Line 1 shows that a particular scatterer can be observed at this time, but only for a very short interval. Then the scatterer associated with the segment of Line 2 becomes dominant. Later another scatterer is observed, associated with lines 3 and 3'. The two segments are associated with the same scatterer, because the lines have the same slopes. With the curve marked 4, again for a short period, a pair of scatterers is observed. This the fact is verified that scatterers appear and disappear.

Conventional processing typically extends over the total time interval shown in FIG. 2, so as to obtain high "resolution." However, with a backscattering behavior of the type shown in FIG. 2 the result will be smearing of responses and their disappearance, rather than high resolution. The actual targets simply do not conform to the model underlying conventional Doppler processing.

Based on these findings, the invention provides the following solution to the problem: If scattering centers do not persist over the times needed to achieve the desired Doppler resolution, then shorter processing times are required. In fact, the processing times cannot be longer than the time intervals over which the scatterers persist. By processing a scatterer only over the time interval in which it is dominant (it governs the behavior of the phase then), there is no need to resolve it from other scatterers. A measurement of its position can be made to a small fraction of the width of the resolution cell, obtaining an accuracy corresponding to long processing times without actually using long processing times. The same approach also prevents the smearing of scatterers, and hence their loss. The central idea of this invention is to use short processing times, despite the fact that the required image detail appears to demand long processing times. Typically, the processing time may be 10% of that used with conventional processing, but the number is not important.

There are various ways in which the processing principle can be implemented. Three methods successfully have been tried, but there might be other implementations of the same principle. The common feature is always that the signal is analyzed in segments much shorter than used with conventional processing.

As one possibility, such phase segments as marked in FIG. 2 can be identified. Where, aside from the high-frequency variations, the trend is linear, the slope can be measured, which is the Doppler and can be translated in to cross range. For each such slope a point can be marked at which the corresponding scatterer appears. If so desired, the strength of the scatterer can be obtained from the signal amplitude, shown in FIG. 2, at the time at which the scatterer appears. An analysis of pairs of scatterers, such as in Segment 4 of FIG. 2, follows similar rules.

As another possibility, short-term Fourier transforms can be used, rather than the much longer-term Fourier transforms of conventional processing, followed by analysis of the result to determine positions of scatterers. The third approach is to correct phase jumps, such as the one which occurs between Segments 3 and 3' in FIG. 2, so as to lengthen the available time interval. Although in this case the transform then is taken over a longer time, note that the signal first was analyzed on a pulse-by-pulse basis through its phase function. This a the extreme of a "short" processing time.

To indicate the success obtained with this new processing approach, FIG. 3 shows a reprocessing of the same image as in FIG. 1 with the novel method. The increase in the image detail is striking. Note that the overall dimensions of the tank are now accurately defined. Even the machine gun on top of the turret is visible. In this case, the image was obtained by using short-term Fourier transforms, with a duration about 10% of that which would have to be used with conventional processing. However, the other forms of implementation have been tested on a moving tank and a moving jeep. In either case, highly accurate positions of scatterers are obtained without losing most of the scatterers, using processing times far too short to provide high Doppler resolution with conventional processing methods. Moreover, the approach is robust because there is no attempt to utilize superresolution techniques, which work on paper but not in practice. Note that the images shown here as illustrations were derived from real radar data taken under practical conditions.

For a more detailed description of the invention, section 5 of the technical report AFWAL-TR-81-1151 is substantially reproduced below.

5. IMAGING OF THE STATIONARY TANK 5.1 The Problem With Conventional SAR Images

The objective is to obtain processing results superior to those achieved with conventional SAR processing. To this end, it is necessary first to determine the specific problem with SAR imaging. As already pointed out, the problem is caused by the backscattering behavior of man-made targets. However, the better the details of the problem are understood, the more successful should be the evolution of superior processing methods. It is desired to improve upon the image of a stationary tank, but it would also be desirable to improve upon the results achieved for the moving tank, because a stationary tank should be easier to analyze than a moving tank.

FIG. 4 shows a Doppler plot for a range cell just outside the boundaries of the tank. Each frame represents the spectrum within the range gate corresponding to the computer printout for the center of the ordinate, in meters. In the case of FIG. 4 this would be a range gate of 4.887 m or 16.00 ft. Otherwise, the ordinate is a true time axis, with the time separation of adjacent frames being 0.0375 sec. Since there are 120 frames plotted, the spectra are shown over a total time period of 120×0.0375=4.5 sec. The plot shows that the range gate contains only clutter; or at least that the clutter in the range gate completely overrides any return from the tank (the gate may contain the gun barrel).

FIG. 5 shows the same type of plot one range gate farther away, which is the gate at 5.18 m or 17.0 ft. The only change in the plot parameters is that the amplitude scale of FIG. 5 was decreased by the factor of 0.27 from that of FIG. 4, so that the amplitudes in FIG. 4 are correspondingly higher. The tank response now appears at about the zero-Doppler vertical line, where it should appear due to the motion compensation. FIG. 6 shows the same type of plot for Gate 6.4 m or 21.0 ft, FIG. 7 for Gate 9.753 m, and FIG. 8 for Gate 10.97 m.

The characteristics of these Doppler plots evidently are similar. The responses are fluctuating with time, and the Doppler positions of the peaks are varying slightly. (The distortions at the top of the plot are due to the phase jumps in the data and should be ignored). Under the circumstances the Doppler coordinate is translated into cross range in meters by multiplication with the factor of 0.0929. This, for a tank width of 3.6 m, the Doppler spread of the tank becomes 3.6/0.0929=38.8 Hz. A look at the Doppler scale of the plots shows how small this Doppler spread actually is, and how significant the "slight" variations in the peak positions are. With conventional SAR processing one simply increases the processing time for the spectrum over the entire 4.5 sec, obtaining only a single spectrum. It is necessary to determine why this leads to difficulties.

The situation is most easily understood by plotting the phase function of the signal in a given range gate. This phase function is shown for Gate 5.18 m in FIG. 9 together with the amplitude of the signal. Two segments of the phase function appear over which it has (aside from the noisy variations) essentially constant slope. Since Doppler is the time derivative of phase, the response has constant Doppler over these segments. The processing time could be increased to cover the individual segments without any detrimental consequences. However, if the processing time is extended over the entire observation time, as is done with SAR processing, evidently averaging occurs over various types of phase discontinuities. It would be an unlikely accident if this led to good results.

FIG. 10 is the phase plot for Gate 6.4 m, and it illustrates an additional problem. One now can readily identify several distinct phase slopes, which implies that the Doppler is changing in steps. The practical interpretation is the following: At any one time, a specific scattering center at a fixed location on the tank (and hence having a fixed Doppler) is predominant. With time, the return from this scatterer becomes weak, and another scatterer having a different Doppler becomes predominant, and so on over the observation time. This behavior is confirmed by the amplitude plot on top, which shows increases and decreases that correlate with the changes in the phase slope. Thus, each scatterer is observed strongly only over a short period, after which it is "replaced" by another scatterer. Evidently, if one simply extends the processing time over the entire observation interval, so as to increase Doppler resolution, the resulting "response" will be some meaningless average over the responses from the individual scatterers. This is why SAR processing fails.

Lastly in FIG. 11 is shown the phase and amplitude plots for Gate 9.753 m. The effects are the same, or perhaps somewhat worse. To summarize the situation, one can state the following: Doppler processing, and with it SAR processing, is based on the assumption that each scatterer is a point scatterer whose amplitude is constant over the entire processing time. In actual fact, the scatterers on a man-made target persist only over periods so short that they cannot be resolved in Doppler from each other if an equally short processing time is used. If the processing time is increased in order to achieve "Doppler resolution", the result will be a smearing of the responses into more or less meaningless clutter. Some strong responses may remain near their nominal positions. This is why SAR processing of a stationary tank shows a few responses which are roughly arranged in the from of an L. The two near edges of a tank are responsible for this behavior. However, a similar effect might occur with APCs or even the larger trucks. SAR imaging thus will not allow reliable classification, and evidently is unsuited for identification.

5.2 Choice of an Improved Processing Method

The question to be answered is: Given the peculiar backscattering behavior of man-made targets, what is the best processing method for pinpointing the location of individual scatterers in each range gate? Unfortunately, the investigation needed to answer this question goes much beyond the scope of this program. Hence, the question posed under this program is: What type of processing should be used in order to obtain in the simplest possible fashion a significant improvement over the results with conventional SAR processing?

Earlier the problems caused by amplitude modulation introduced by the backscattering pattern of individual scattering centers were discussed at great length. In order to obtain meaningful results for the moving tank, in this first phase of the program elaborate phase correction procedures were used that offset much of the effects from the pattern modulation. In this second phase, the need for phase corrections have been successfully avoided by using the following argument: If the phase function within a range gate is well behaved over extended sections, the data evaluation can be restricted to those sections and ignore the transitions where phase corrections would be needed. This approach may not allow extraction of all the information within the signals, but it would lead to significant improvement over conventional imagery.

There has not been the opportunity to check whether or not the same approach would also be successful for the moving tank processed during Phase 1 of this program. However, even if the method were restricted to stationary tanks (and other vehicles), it still would be valuable because even under operational conditions it is always known whether the vehicle is stationary before processing for identification. Thus one need not necessarily use the same identification procedure for both stationary and moving vehicles.

The first idea was tested and was to confine giving attention to the linear-phase segments in the phase plot for each range gate, and to determined the corresponding Doppler (which is then converted to cross range). This will yield as many scatterers per range gate as there are identifiable linear segments of different slopes in the phase function. The results were far superior to conventional imagery. However, one concern was that the phase slope could not be measured as accurately as desirable, in particular if the process were to be automated.

The next possibility investigated was that of overcoming this possible limitation by adding to the data in each range gate an artificial strong scatterer with a Doppler offset. This has the consequence of translating a linear phase slope into a sinusoidal variation, with the frequency of the sinusoid indicative of the phase slope, and hence of the Doppler position of the associated scatterer. This approach also was unsuccessful. However, a further investigation led to the conclusion that the same objective could be achieved in a simpler manner by the following processing method.

One treats the data as if conventional SAR imagery were to be performed. However, then one does not use the long processing times required to obtain the desired degree of Doppler, or cross range, resolution, because of the problems discussed earlier. Instead, one uses a rather short processing time, perhaps 10% of what would be used for SAR imagery. Whenever a linear phase segment occurs in a given range gate, the short processing time will lead to a clean response, whose Doppler position will indicate the short processing window, the Doppler response will be so broad that Doppler resolution is useless for the purpose on hand. However, since a linear phase segment occurs only when one scatterer becomes predominant, no Doppler resolution is needed. One still can determine the Doppler position of the response peak to a small fraction of the response width, and this provides the requisite measurement accuracy.

In practice the situation is more complex. Although one scatterer may be predominant, other scatterers will be interfering. If a second scatterer is present, it will introduce an essentially sinusoidal modulation into the linear phase of the dominant scatterer, and this will cause a meandering back and forth of the Doppler history of the dominant scatterer. The mean Doppler of the Doppler history will give the Doppler of the dominant scatterer and the frequency of the meandering will be the Doppler offset of the interference scatterer.

If a third scatterer is significant, it will likewise introduce a phase modulation corresponding to its Doppler offset, and will modulate the Doppler history of the dominant scatterer. There was no time available under this program to develop the theory fully. Intstead, at most a second interfering scatterer was considered. If the situation is more complicated, this will be evident from the processing results, in which case one simply does not perform the measurement. In other words, this program simply used the processor output only for those segments where its interpretation was evident. Although this may be considered only a first usage stage in the full utilization of the information in radar returns. it will be seen that this simple procedure was amazingly successful. Again, note that it was not investigated whether the same procedure also could be used on moving vehicles.

5.3 Processing Results for the Stationary Tank

One important question concerns the choice of the coherent processing time, or the window length. As was explained above, the window length cannot be too large, because one then would "average" over different scattering centers, rather than resolve the scattering centers. On the other hand, if the processing window is too short, the accuracy with which the Doppler position of a response peak can be measured will be poor, and so will the image quality.

As an example, consider the processing time of 0.3 sec used in FIGS. 4 to 8. With Hamming weighting for Doppler sidelobe suppression, the 3 dB response width will then be 4.5 Hz. It is reasonable to assume that the peak position can be measured to about one-thirtieth of the 3 dB, width, which would be 0.15 Hz. The width of the tank was stated to be 38.8 Hz in computer units, which would be 38.8/15.84=2.4 Hz in actual frequency units. A measurement error of 0.15 Hz thus would correspond to 0.15/2.4=6.0% of the tank width, or 22 cm. Since the measurement accuracy of one-thirtieth of the 3 dB width may be optimistic, and a 6% error is larger than desirable, we have used a processing time of 0.6 sec for forming the image. Note, however, that this processing time was selected without any analysis of the problem of optimizing the processing time.

With a processing time of 0.6 sec, there was then calculated in the computer the same type of Doppler plots as shown in FIGS. 4 to 8 and the computer plotted the position of the highest peak in each frame (the plot was done in range rate coordinates rather than Doppler, with range rate converted into Doppler by division by 9.42×10^-4). To convert range rate into actual cross range, one must use the multiplication of 94.6. The following examples illustrate how the two-dimensional image was formed.

FIG. 12 shows the range rate curve for Gate 4.877 m, which is the gate of FIG. 4. The discontinuous nature of the curve implies that there is no usable scatterer. The gate is thus considered empty.

Next, in FIG. 13 is shown the corresponding range rate curve for Gate 5.18 m, the gate of FIG. 5 and the first one containing a measurable return from the tank. Only two segments of smooth behavior are seen in the plot, as marked in the figure. One simply reads the coordinates of the points halfway between maxima and minima, as marked. By multiplying the three coordinates of 0.018 and 0.015, and 0.02 m/sec by the conversion factor of 94.6 we obtain three cross-range coordinates of scatterers: 1.70, 1.42, and 1.89 m. These scatterer positions then are plotted for Gate 5.18 m. There will be no attempt to utilize the modulation in the curve of FIG. 13 for locating additional scatterers because for the first range gate containing any tank return these additional weak scatterers are likely to be clutter.

FIG. 14 gives the plot for Gate 6.4 m (also of FIG. 6). The curve is smooth throughout the entire observation interval. It should not be used toward the end, where the data have the phase jump problem discussed earlier. One may also not want to use the curve where it changes over a large interval, because there may be an accuracy problem. Otherwise one simply marks the halfway points and converts range rate to cross range, as before. In the figure there is also marked a modulation period of 0.034 sec, which corresponds to a Doppler difference of 1/0.034=29.4 Hz. By multiplication with 0.0929 this is converted into a cross-range difference of 2.73 m. The average between the two range-rate readings over this period is 0.01/m/sec, which by multiplication with 94.6 converts into a cross range of 0.95 m. The direction of the second scatterer, whether it is to the left or right of the dominant scatterer, could be determined by analyzing the amplitude/phase plot of the gate. For simplicity, one simply argues here that if it were counted to the right, it would fall by a large margin outside the boundaries of the tank. Thus, one takes the displacement negative at -2.73 m, which places the scatterer at -2.73+0.95=-1.78 m. If the procedure were repeated with the modulation period to the right of the one marked in FIG. 14, another scatterer would be found, and it also would be within the boundaries of the tank. We have attempted to map all scatterers.

The last example is given by FIG. 15 for Gate 9.753 m, the same gate as for FIG. 7. The same rules are followed for deriving the position of scatterers as illustrated above. The rules evidently are rather crude, and in all likelihood could be refined with further study. However, the coarseness of the rules also implies that automation of the analysis would be relatively simple.

The result of processing all range gates and performing the indicated measurements is given by FIG. 3. It is evident that the method has led to a relatively accurate definition of the overall size and shape of the vehicle (the object next to the tank again is the corner reflector, which appears to be resolved by the radar). The measured scatterer locations often coincide with actual discontinuities on the tank. However, the processing method would be expected to show the shifting scattering center on smooth surfaces as the aspect angle is changing because of the platform motion. No matter how much the processing is improved, one will always expect the speckled image familiar from coherent light. It is interesting to note that the image appears to reproduce correctly the offset machine gun on top of the turret, whereas the main gun has not appeared with this type of processing. This may be due to differences in the surface roughness. The machine gun might not have as smooth a barrel as the main gun. This point has not been investigated further.

For comparison, note the image that would be obtained with conventional SAR processing. For a better comparison with FIG. 3, it is desirable to change the method of presentation somewhat from that used earlier. First, if a response peak extends over several range gates, there will be plotted in the image a response in every range gate, and not just in the gate where the response is strongest. Second, in order to distinguish between meaningful responses and spurious responses introduced by the smearing of the major responses, bigger circles are used for the major responses and smaller circles for minor responses or, perhaps, spurious responses.

The result is shown in FIG. 1. If one considers for the moment only the major responses indicated by the larger circles, one observes the typical L-image of stationary tanks. The smaller circles help somewhat to fill out the image, but not much. A coparison with FIG. 3 shows that the stationary tank is much less well defined than with out modified processing approach. One cannot obtain the overall dimensions with acceptable accuracy from FIG. 1 and identification of the type of tank appears impossible. It would require the comparison with images from trucks to determine whether the L-shape characteristics of the image even permit relitable discrimination between trucks and tanks.

Section 5 of the technical report ends here.

COMPUTER PROGRAM

A copy of a "FORTRAN Listing of Computer Code for Algorithim Execution" (MRI 186-21A) is attached hereto as an appendix, and an explanation thereof follows.

When scatterers have highly lobed backscattering patters a particular scatterer may dominate the radar signal from a target for only a short time. During that interval the slope of the phase function indicates the cross-range position of the dominant scatterer. Below is described an automated algorithm for finding and measuring phase slopes corresponding to scatterer locations.

The program begins with the unwrapped phase function of the return signal for a specific range gate. This is the phase function generated by taking the phase measured for each pulse and, if necessary, adding or subtracting a full phase cycle in order to minimize the phase change from one pulse to the next. The duration of the signal is normalized to one second. The program then locates jumps in the phase and uses them to define breakpoints, the limits of time intervals which correspond to dominance of particular scatterers.

The program searches the phase function for points of extremal curvature (there will be an extremum on each side of a phase jump) with curvatures of magnitude greater than SCHNG (default value 0.3 cycles/sec/sec; although default values given throughout this description have been derived through applicatoon of this algorithm to a variety of data sets, optimal values for any particular data may vary). If the phase difference at the extremal points is greater than PJMP (default value 0.23 cycles) and the corresponding amplitude function has a minimum between the times of the extremal phase curvatures then a phase jump is said to have occurred. An amplitude minimum is accepted only if the the amplitude at the minimum is less than AMPIN (default value 0.15) times the amplitude of the first amplitude peak after the minimum under consideration. When a phase jump occurs the extremal points surrounding it are taken as breakpoints.

Each segment resulting from the above step is considered for further segmentation. A linear least-squares fit is made to the phase of each segment. No further segmentation is done to segments of duration less than SEGMIN (default 0.00093 sec), segments of duration less than PJSEGLN (default 0.04 sec) with fewer than NMBREXT (default 3) residual extrema in the difference between the fit line and the phase, and segments of duration less than PJSEGLN with fewer than two amplitude extrema and a standard deviation about the fit less than SDMIN (default 0.015 cycles).

Remaining segments are searched for deviation from linearity. Additional breakpoints are added at points of maximal deviation. The first additional breakpoints are assigned through a multiple (per segment) breakpoint assignment algorithm, as follows. The intersection points of the fit line and the phase are found. Breakpoints are assigned between these intersection points at the times of maximal deviation between the line and the phase, if this deviation is greater than 0.25 cycles.

Linear least-squares fits are made to the phases of the resultant segments. If the maximum deviation between the fit line and the phase occurs at eh end of a segment, and is greater than ED (default value 0.14 cycles), the end point is dropped from consideration for phase slope measurement and a breakpoint is assigned at the new segment endpoint. If the maximum deviation occurs away from the segment endpoints, and is greater than CD (default value 0.17 cycles), an additional breakpoint is assigned at the point of maximum deviation. Whenever a new segment is defined a linear least-squares fit is made and deviations are examined. Breakpoints are assigned and the process iterated until all deviations are less than ED and CD.

The resultant segments are examined for compliance with criteria which indicate a single dominant scatterer, measured with sufficient accuracy over a long enough duration so that the phase slope corresponds to the scatterer cross-range location. The first criterion is that an acceptable segment have a phase slope between MNSLOPE and MXSLOPE (settings supplied be the user should correspond to the image size unless target limts are known, or devived by a related algorithm). An acceptable segment must also have minimum length MSL (Default 012 sec). Any segment with length greater than SLL (default value 0.33 sec) is accepted without further testing.

Candidate segments shorter than PJSEGLN with either a phase jump break at an endpoint or with NRE (default 3) or fewer residual extrema in the difference between the phase and the fit line are rejected. Segments with a standard deviation of the phase about the fit line greater than MAXSD (default 0.07 cycles) are rejected, as are segments with a standard deviation less than MINSD (default 0.01 cycles) and fewer that two amplitude maxima.

Segments which have not been eliminated by the above tests are accepted. Their phase slopes indicate scatterer crossrange locations. The above procedure also may be applied to the return spectrum for a specific cross-range gate. In that case the measured phase slope indicates the range location of a scatterer dominant for part of the radar bandwidth.

It is understood that certain modifications to the invention as described may be made, as might occur to one with skill in the field of the invention, within the scope of the appended claims. Therefore, all embodiments contemplentated hereunder which achieve the object of the present invention have not been shown in complete detail. Other embodiments may be developed without departing from the scope of the appended claims.

The computer program listing follows. ##SPC1##

Claims

What is claimed is:

1. A signal processing method for data collected by airborne radar apparatus, the signal processing method being used to form two-dimensional images of ground targets, in which resolution in one dimension is provided by resolution of range cells, and resolution in the other dimension is provided by resolution of Doppler cells;

said ground targets being of types which have scattering units which are large compared with the radar wavelength and therefor have a highly lobed backscattering pattern, so that over a given observation time comprising a range cell the scattering units produce responses which become stronger and weaker, alternately effectively appearing and disappearing, the response for each scattering unit having a phase function, each scattering unit being identified by the slope of the phase function during intervals in which there is a dominant response from only one scattering unit;

wherein the signal processing method includes the step of processing the responses over time intervals which are short compared to said given observation time, so that measurements of the positions of the scattering units are made with an accuracy to a small fraction of the cross-range width of the range cell being processed, the scattering units being identified by the slopes of the responses.

2. The signal processing method according to claim 1, wherein during intervals in which there are responses from two scattering units, the phase function has a slope for one scattering unit, plus an approximately sinusoidal modulation at some frequency impressed by the response for the other scattering unit, and wherein the signal processing step over a short time interval in which the sinusoidal modulation appears includes the step of identifying the two responses, so that the positions of these scattering units are also made with an accuracy which is a small fraction of the width of said range cell.

3. The signal processing method according to claim 2, wherein the step of processing the responses over short time intervals comprises identifying phase segments where, aside from high-frequency variations, the trend is linear, and measuring the slope, which is the Doppler, and translating the slip to a measurement of the cross range position, so that for each slope a point is marked at which a corresponding scattering unit appears, with a similar analysis for intervals in which there are responses from two scattering units.

4. The signal processing method according to claim 1, wherein the step of processing the responses over short intervals comprises finding and measuring phase slopes corresponding to scatterer locations as follows:

beginning with the unwrapped phase function of the return signal for a specific range gate, which is the phase function minimize the phase change from one pulse to the next, the duration of the signal being normalized to one second, then locating jumps in the phase and using them to define breakpoints (the limits of time intervals which correspond to dominance of particular scatterers);

searching the phase function for points of extremal curvature (there will be an extremum on each side of a phase jump) with curvatures of magnitude greater than a value SCHNG (if the phase difference at the extremal points is greater than a value PJUMP and the corresponding amplitude function has a minimum between the times of the extremal phase curvatures then a phase jump is said to have occurred), an amplitude minimum being accepted only if the the amplitude at the minimum is less than a value AMPMIN times the amplitude of the first amplitude peak after the minimum under consideration (when a phase jump occurs the extremal points surrounding it are taken as breakpoints);

each segment resulting from the above step being considered for further segmentation, making a linear least-squares fit to the phase of each segment, no further segmentation being done to segments of duration less than a value SEGMIN, segments of duration less than a value PJSEGLN with fewer than a value NMBREXT residual extrema in the difference between the fit line and the phase, and segments of duration less than a value PJSEGLN with fewer than two amplitude extrema and a standard deviation about the fit less than a value SDMIN;

searching remaining segments for deviation from linearity, additional breakpoints being added at points of maximal deviation, the first additional breakpoints being assigned through a multiple (per segment) breakpoint assignment algorithm, as follows, the intersection points of the fit line and the phase being found, breakpoints being assigned between these intersection points at the times of maximal deviation between the line and the phase, if this deviation is greater than 0.25 cycles;

making linear least-squares fits to the phases of the resultant segments, if the maximum deviation between the fit line and the phase occurs at the end of a segment, and is greater than a value ED, the end point being dropped from consideration for phase slope measurement and a breakpoint being assigned at the new segment endpoint, if the maximum deviation occurs away from the segment endpoints, and is greater than a value CD, an additional breakpoint is assigned at the point of maximum deviation, whenever a new segment is defined a linear least-squares fit being made and deviations examined, breakpoints being assigned and the process iterated unit all deviations are less than the values ED and CD;

examining the resultant segments for compliance with criteria which indicate a single dominant scatterer, measured with sufficient accuracy over a long enough duration so that the phase slope corresponds to the scatterer cross-range location, the first criterion being that an acceptable segment have a phase slope between values MNSLOPE and MXSLOPE wherein an acceptable segment must also have minimum length MSL, any segment with length greater than a value SLL being accepted without further testing;

rejecting candidate segments shorter than a value PJSEGLN with either a phase jump break at an endpoint or with a value NRE or fewer residual extrema in the difference between the phase and the fit line, segments with a standard deviation of the phase about the fit line greater than a value MAXSD being rejected, as are segments with a standard deviation less than a value MINSD and fewer than two amplitude maxima;

segments which have not been eliminated by the above tests being accepted, their phase slopes indicating scatterer cross-range locations.