Abstract:
Distortions inherent in the formation of a range/doppler image by an airborne Synthetic Aperture Radar (SAR) of a ship under the influence of roll, pitch, and yaw motions characteristic of sea state conditions are removed in conjunction with a least squares linear regression solution to doppler processed interferometric azimuth angle data derived from ship radar reflections, resulting in a scaled high resolution range/doppler image representative of the true range/cross range (azimuth) distribution of the ship, so that continuous automatic tracking of a cursor imbedded in a single designated resolution cell of the ship&#39;s displayed image, essential to carrying out precision standoff command weapon guidance to that selected ship target cell, can be accomplished.

Description:
The Government has rights in this invention pursuant to contract N00019-79-C-0075 awarded by the Department of the Navy. 
    
    
     BACKGROUND OF THE INVENTION 
     The present invention relates to radar controlled weapons systems and, more particularly, to a method and apparatus for generating real-time high resolution Synthetic Aperture Radar (SAR) imagery from an airborne platform of a translating ship under the influence of roll, pitch, and yaw motions characteristic of sea state conditions. 
     A highly resolved ship image on an airborne display permits targeting to a particular part of the ship for the purpose of standoff command guidance weapon delivery. As such, this invention relates to airborne SAR systems used for generating real-time high resolution imagery of a ship target under the influence of sea state conditions and for accurately measuring and tracking the range and azimuth angle of a designated resolution cell within the aforesaid displayed target area so as to enable the accurate delivery of an air-to-ground missile or glide bomb from the SAR bearing aircraft to the ship target. The manner of weapon guidance depends upon reducing to zero the relative range and azimuth angle between weapon and designated target resolution cell. 
     Existing airborne synthetic aperture radar systems have the capability for producing a high resolution image of stationary ground targets by virtue of the rotational motion of the SAR with respect to the target field. High resolution in both range and cross-range (azimuth) directions are achieved with short-pulse radar techniques for range resolution, and by digital signal doppler processing for azimuth resolution, so as to produce a two-dimensional image. The application of such known SAR processing techniques to a ship at sea, however, can readily produce distorted and improperly focused ship imagery, precluding ship recognition, because of ship rotational motions brought about by sea state conditions. 
     A system for ameliorating some of the adverse effects brought about by ship rotational motions is described in a commonly assigned copending United States patent application of Sol Boles, Ser. No. 389,367, filed concurrently herewith, entitled &#34;Range/Azimuth Ship Imaging For Ordnance Control&#34;, the disclosure of which is incorporated herein by reference. The disclosed invention provides an undistorted two-dimensional image of the ship from a direct plot of range versus interferometrically determined azimuth angle of all essential scatterers comprising the ship. Notwithstanding the advantages achieved in conjunction with the disclosed invention, it does have certain operational limitations and drawbacks in certain applications. For example, the aforementioned techniques may not be capable of generating SAR imagery of high quality when the ship target is located in a high sea state environment because of the resulting degradation in the signal to noise ratio. 
     The invention disclosed herein is related to that described in the aforementioned application of Boles, entitled &#34;Range/Azimuth Ship Imaging For Ordnance Control&#34;, but by a unique processing implementation, is capable of achieving a significant advantage in ship image quality and signal to noise ratio. The latter advantage, in turn, permits operation to greater ranges than would otherwise be possible by the techniques described in &#34;Range/Azimuth Ship Imaging For Ordnance Control&#34; for given transmitter and antenna parameters. This advantage is brought about by processing and displaying a high resolution range/doppler image of the ship target, with the distortions inherent in such a presentation due to ship scatterer rotational motions brought about by sea state conditions removed in accordance with the teachings of this invention, rather than processing and displaying a range/interferometrically determined azimuth angle image. 
     As those skilled in the art will recognize, the range/doppler representation could have a very substantial signal to noise advantage over the range/azimuth approach to imagery because of the image quality limitations brought about by angle glint errors associated with the latter. 
     It is accordingly an object of the present invention to provide a method and apparatus for generating high resolution synthetic aperture radar displayed imagery from an airborne platform of a ship under the influence of sea state conditions. 
     It is another object of the present invention to provide a method and apparatus for determining the cross line-of-sight relative velocity associated with aircraft and ship rotational motions requisite to achieving prescribed cross-range (azimuth) resolution in the displayed image representation of the ship. 
     It is a still further object of the present invention to provide a method and apparatus for eliminating the distortions inherent in range/doppler ship imagery brought about by ship rotational motions accompanying sea state conditions by forming a scaled SAR high resolution range/doppler image projection of the ship on the basis of interferometric azimuth angle measurements derived from doppler processed range sampled data, so as to permit the cursoring and tracking of a particular range/doppler resolution cell of the ship target for the purpose of carrying out standoff command guidance weapon delivery to said target. 
     SUMMARY OF THE INVENTION 
     The foregoing and other objects and advantages which will be apparent in the following detailed description of the preferred embodiment, or in the practice of the invention, are achieved by the invention disclosed herein, which generally may be characterized as a method and apparatus for forming, in conjunction with an airborne synthetic aperture radar system having an interferometer antenna and a range versus azimuth angle display, a range/doppler image having improved azimuth resolution of a ship target under the influence of sea state conditions. 
     In accordance with the teachings of the invention, a range/doppler image having improved azimuth resolution of a ship target under the influence of sea state conditions is formed by processing the received signals from the scatterers comprising the ship target to obtain an estimate of the net doppler producing cross line-of-sight velocity of the radar bearing aircraft relative to the ship, and determining from the estimated cross line-of-sight relative velocity the values of predetermined system parameters to be used predictively in the succeeding integration interval in forming the range/doppler image. 
     The formed range/doppler image having azimuth resolution is displayed in a range versus azimuth angle coordinate frame, and continuous automatic tracking of a cursor imbedded in a single designated resolution cell of the ship&#39;s displayed image, essential to carrying out precision standoff command weapon guidance to that selected ship target cell, is thereby enhanced. 
     In addition, techniques are provided for removing the isodop distortion inherent in the formed range/doppler image, and for centering the displayed range/doppler image of the ship target. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     Serving to illustrate an exemplary embodiment of the invention are the drawings of which: 
     FIG. 1 illustrates aircraft and ground target geometrical relationships governing the formation of a synthetic aperture radar range/doppler image; 
     FIG. 2 illustrates aircraft and ship target geometrical and motional relationships governing the formation of a range/doppler ship image; 
     FIG. 3A depicts doppler versus azimuth angle functional relationships for assumed ship rotational rates of -W zz  and +W yy  ; 
     FIG. 3B illustrates a range/doppler image for the assumed rotational rates as represented in FIG. 3A; 
     FIG. 4A depicts doppler versus azimuth angle functional relationships for assumed ship rotational rates of +W zz  and +W yy  ; 
     FIG. 4B illustrates a range/doppler image for the assumed rotational rates as represented in FIG. 4A.; 
     FIG. 5 illustrates the determination of K R  and θ o  via a linear regression curve fit; and 
     FIG. 6 is a block diagram illustrating the physical embodiments of the present invention. 
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
     In order to provide for a better understanding and background of the present invention, it will be useful to describe and discuss briefly the contemporary art and some of its resulting problems. 
     Referring to FIG. 1 the aircraft and ground target geometrical relationships governing the formation of a synthetic aperture radar range/doppler image are illustrated. As shown therein, for an aircraft flying with respect to a fixed ground target field, the doppler shift, f d , associated with the signal reflected from a scatterer located at azimuth and elevation angle differences, Δθ and Δη, respectively, with respect to a boresight line drawn to focus point, F, is given by ##EQU1## where λ is the radiated wavelength, V xx  is the line-of-sight velocity component, V yy  is the cross line-of-sight velocity component, and V zz  forms an orthogonal set with V xx  and V yy  as shown. Also, Δθ is the azimuth angle in the slant plane containing V xx  and V yy  as shown. Now, from, ##EQU2## where, h=aircraft altitude, R o  is the slant range to the ground focus point, and η o  is the elevation angle to the focus point, after differentiation, ##EQU3## where ΔR is the slant range difference from the focus point to the scatterer in question. Substituting equation (2) into (1), yields for the doppler shift, ##EQU4## assuming a flat target surface. 
     According to contemporary art, a synthetic aperture radar image can be generated by sorting the signals received in a given range resolution increment located at ΔR from the range of the focus point, F, according to their doppler frequencies, using digital signal processing techniques. This is done after the application of a suitable phase correction to all incoming signals on a pulse-to-pulse basis, to remove the effects of the line-of-sight velocity component, V xx , on the measured doppler frequency, which, in turn, is accomplished on the basis of the system&#39;s estimate of line-of-sight velocity, V xx . The resulting doppler frequency is then proportional to azimuth angle from boresight, Δθ, (see equation (3)) plus an additive constant, -2V zz  (tan η o )ΔR/λR o . If this constant is removed in each range bin, by adding to each sorted doppler frequency an amount, 2V zz  (tan η o )ΔR/λR o , on the basis of the system&#39;s estimate of the pertinent variables, then a realistic range/azimuth angle map results, representative of the gridded ground matrix depicted in FIG. 1. It is also seen from equation (3) that the constant of proportionality between doppler and azimuth angle is 2V yy  /λ, so that a knowledge of aircraft cross line-of-sight velocity, V yy , is required to establish the angular resolution associated with a prescribed doppler filter spacing. 
     Referring now to FIG. 2, the aircraft and ship target geometrical and motional relationships governing the formation of a range/doppler ship image, the subject of the present invention, are illustrated. As shown therein, assuming that the ship is translating and undergoing roll, pitch, and yaw motions as a result of sea state conditions, the doppler-producing ship rotational motions can be represented by instantaneous axes of rotation, W zz  and W yy , both perpendicular to the line-of-sight vector from the aircraft to the center of rotation of the ship, and lying respectively in planes containing V xx  and V zz  and V xx  and V yy . It can be readily shown that the additional doppler shifts of the signals observed at the radar resulting from W zz  and W yy  rotations, for ship scatterers assumed to lie in the horizontal datum plane of FIG. 2, are given respectively by 2W zz  ΔθR o  /λ and -2W yy  (tan η o )ΔR/λ. Accordingly, the doppler shift expression, including both aircraft and ship motions, by analogy to equation (3), is given by, ##EQU5## where the following definitions are employed, 
     
         V&#39;.sub.xx ≡V.sub.xx                                  (5) 
    
     
         V&#39;.sub.yy ≡V.sub.yy +W.sub.zz R.sub.o                (6) 
    
     
         V&#39;.sub.zz ≡V.sub.zz +W.sub.yy R.sub.o                (7) 
    
     and where V xx , V yy  and V zz  are now the relative translational velocities between aircraft and ship, given by, ##EQU6## 
     As a prerequisite to doppler processing the date for the formation of a range/doppler image with contemporary art, the system estimate, V&#39; xx  of the true line-of-sight velocity, V&#39; xx  is used to motion compensate the data for focusing and image centering to correct for both aircraft and ship translational motions. Defining, 
     
         V&#39;.sub.xx =V&#39;.sub.xx +ε.sub.v&#39;xx                   (11) 
    
     where ε v&#39;xx  is the error in line-of-sight velocity estimate, equation (4) becomes, ##EQU7## Without further corrections, a doppler processed ship image would lack display centering in proportion to the error in line-of-sight velocity, ε v&#39;xx . Since ship velocity is normally unknown, this error source could ordinarily be of considerable magnitude. Furthermore, since for moderate to severe sea state conditions, (see equations (6) and (7)) 
     
         |W.sub.zz R.sub.o |&gt;&gt;|V.sub.yy | 
    
     and, 
     
         |W.sub.yy R.sub.o |&gt;&gt;|V.sub.zz |, 
    
     the doppler-producing cross line-of-sight velocities, V&#39; yy  and V&#39; zz  in equation (12) are largely determined by ship angular rotation rates, W zz  and W yy , which too, are unknown. Since azimuth resolutin for fixed filter bandwidth depends on V&#39; yy , large uncertainty in this value precludes the achievement of desired image azimuth resolution, which could thereby preclude image recognition. Scaling of the horizontal (azimuth) dimension of the resulting image is also correspondingly unknown. Uncertainty in V&#39; zz  prevents the proper removal of the image distortion (commonly referred to as &#34;isodop&#34; distortion) as a function of the range extent of the image, which could be very severe for the normally large values of V&#39; zz  encountered in a high sea. 
     FIG. 3B illustrates the &#34;isodop&#34; distortion inherent in a range/doppler image (shown dashed) associated with an actual range/azimuth angle representation (shown solid) for assumed ship rotational rates of -W zz  and +W yy . FIG. 3A depicts a functional plot of equation (12) for this condition (assuming ε v&#39;xx  =0) and is the basis for the range/doppler image of FIG. 3B. (The negative Δθ axis is taken to correspond to positive doppler which would be the case as depicted by FIG. 1, assuming aircraft motion only.) Similarly, FIGS. 4A and 4B illustrate the same for +W zz  and +W yy . It is to be noted that in addition to the isodop distortion exhibited by FIG. 4B, the image has also been inverted along the azimuth dimension due to the change in sign of the doppler gradient with respect to azimuth angle. 
     The objectives and goals of the present invention are achieved with a novel processing approach, wherein the velocity quantities, ε v&#39;xx , V&#39; yy , and V&#39; zz , and hence the ship translational and rotational rates, are determined from the body of data available to the SAR/Interferometer System. 
     A doppler cell index, i f , is defined, where ##EQU8## where Δf&#39; d  is the doppler filter bandwidth established in the digital signal processing. Also, doppler filter bandwidth relates to the synthetic aperture integration time, T, by 
     
         Δf&#39;.sub.d =1/T tm (14) 
    
     where 
     ti N P  =Tf r ,                                      (15) 
     N P  is the number of pulses integrated during the synthetic aperture and f r  is the pulse repetition frequency. Furthermore, a range cell index, i R , is defined, wherein, ##EQU9## in which RR=desired range resolution. Combining equations (13) through (16) with (12), and rearranging, results in, ##EQU10## which expresses azimuth angle in terms of doppler and range cell indices. Equation (17) is seen to be of the form, Δθ=θ o  +K f  i f  +K R  i R , and represents a family of straight lines, similar in form to FIGS. 3A and 4A. 
     By plotting measured interferometric azimuth angle, Δθ, for all range and doppler cell indices, i R  and i f , including ship returns, and using standard least squares linear regression techinques to fit straight lines to each range bin of data, a set of quantities are obtained from the resulting data, defined as follows: 
     θ o  ≡the averaged (over all range bins) azimuth angle intercept, corresponding to i R  =0 at i f  =0 
     K f  ≡the averaged (over all range bins) rate of change of Δθ with respect to i f  for constant i R   
     K r  ≡the averaged (over all range bins) rate of change of Δθ with respect to i R  for constant i f . 
     By analogy to equation (17), therefore, the following equalities are written: ##EQU11## and, 
     from which, the three velocity terms, ε v&#39;xx , V&#39; yy  and V&#39; zz  are derived: ##EQU12## 
     The regression constants, K f , K R , and θ o  are determined as follows: 
     The slope of the straight line fit (K f ) to doppler versus azimuth angle measurement data in the range bin denoted by i R  is found from standard least squares linear regression formulae, ##EQU13## where i f  and (Δθ)i f  are doppler cell index and interferometric azimuth angle, respectively, of that cell, and w i .sbsb.f is a relative weight proportional to doppler cell signal power. The summations are carried out over all doppler filters processed in each range bin, i R . An improved slope estimate, K f , for use in equations (21) through (23), is obtained by averaging over all range bins (R+1), ##EQU14## 
     Similarly, in each range bin, the azimuth angle intercept, (Δθ)i R  for i f  =0, is found from, ##EQU15## 
     As illustrated in FIG. 5, the values of θ o  and K R  for use in solving for the velocity terms ε v&#39;xx , V&#39; yy , and V&#39; zz  in equations (21) through (23), are found by plotting ##EQU16## versus i R  and determining slope and intercept values, using similar regression equations. Accordingly, ##EQU17## 
     In terms of the invention, the determined error in system line-of-sight velocity, ε v&#39;xx  permits a continual (aperture to aperture) update of the system V&#39; xx  estimate, thereby eliminating image azimuth centering errors by providing for an exact motion compensation correction for aircraft to ship net line-of-sight velocity. Furthermore, referring to equation (12), by adding to the doppler frequency, f d  of every resolution cell of each range bin, i R  an amount, ##EQU18## or equivalently, ##EQU19## rounded to the nearest unit, the isodop distortion inherent in the image is removed, so that the resulting range/doppler matrix truly represents the range/azimuth angle image of the ship. In this manner, the dashed image of FIG. 3B is restored to its true range/angle counterpart; the dashed image of FIG. 4B is transposed, after isodop correction, to the solid image shown. The processor recognizes, however, that the sign of W zz  is such as to have caused mirror imaging of the doppler image about the range axis, so that the isodop corrected image is accordingly reversed to agree with the true range/azimuth distribution of the ship target. 
     From equation (12) and the knowledge of V&#39; yy  determined from equation (22), the doppler filter bandwidth and filter frequency spacing, Δf&#39; d , is given by, ##EQU20## for presecribed azimuth resolution (Δθ) res  : Recognizing that linear azimuth resolution, d az , is 
     
         d.sub.az =R.sub.o (Δθ).sub.res.,               (30) 
    
     then from equations (29), (30), and (22), ##EQU21## and correspondingly, the synthetic aperture integration time, T, is given by T=1/Δf&#39; d . Also, the pulse repetition frequency, f r , for specified array length, N p , and azimuth resolution, d az , is given by, ##EQU22## Filter bandwidth, integration time, and pulse repetition frequency are used predictively for the next aperture, in terms of the present estimated value V&#39; yy . Furthermore, range performance can be extended appreciably by smoothing the derived velocity constants over many apertures using standard curve fitting or filtering techniques. 
     Ship translational motion must be tracked so that antenna boresignt and range swath start bear a constant relationship with respect to the ship. Interferometric azimuth angle data from each useful ship resolution cell are averaged on an array-to-array basis for antenna azimuth control. The updated line-of-sight velocity obtained from the regression solution is tracked so as to advance or retard the range swath start trigger in accordance with ship as well as aircraft motion. Range rate and azimuth rate so determined are applied to the system computer so as to also steer antenna boresight in both azimuth and elevation in accordance with both ship and aircraft translational motions. 
     Significant returns from ship scatterers at different elevation angles, such as from the mast of a ship, could produce data points which fall considerably off the straight line plot of azimuth angle versus doppler frequency due to their different velocities from those at similar azimuth angles but at a lower datum plane, thereby perturbing slope and intercept values used in the regression algorithms. Any such points can be eliminated by excluding any data points outside of a prescribed threshold band about the least squares fit formed from all available data points, and refitting the data after such exclusion. Should elevation angle measurement also be implemented in addition to azimuth, then elevation measurements can also be used as an added discriminant in the exclusion process. 
     For purposes of weapon guidance, a cursor is moved about the displayed image and placed over a particular range/doppler resolution cell, the designated target cell. The appropriate values of range and interferometrcially determined azimuth angle can be extracted from the system for relative range/azimuth angle weapon guidance of a weapon borne transponder to the designated target cell. This is done by recognizing from equation (17), that after isodop correction, the smoothed azimuth angle of any designated range/doppler cell of indices i R  and i f , can be expressed by, 
     
         Δθ=θ.sub.o +K.sub.f (i.sub.f)            (33) 
    
     where i f  is the doppler cell index of the cursor designated filter containing the isodop corrected target signal. The techniques embodied in this invention are seen to provide, in addition, a substantially improved azimuth angle accuracy of a designated target cell than would ordinarily be achieved on the basis of an inteferometric measurement to that cell alone, since the azimuth angle estimate denoted by equation (17), after isodop correction, Δθ=θ o  +K f  i f , derives benefit from all the azimuth angle measurement entering into the regression estimates. 
     As time progresses, cursor tracking of a designated target cell is accomplished in terms of its predicted range/azimuth location (suitably referred to the range/doppler coordinates of the display via available interferometric data) relative to the SAR bearing aircraft on the basis of the derived relative translation between aircraft and ship from aperture to aperture. 
     If cursoring is done to a resolution cell near the extremity of the ship, however, then cursor misalignment to that cell from frame to frame of from 1 to 2 resolution cell lengths could normally result principally due to ship yaw motions. This cursoring track error can be recovered, however, because of the fact that W zz , the essential contributor to such tracking error, is computable from equation (6), where V&#39; yy  is computed by equation (22) and V yy , the aircraft to ship relative cross line-of-signt velocity, is obtained from airborne inertial measurements and ship translational velocity estimates. From a knowledge of W zz , range and azimuth cursor track corrections, d (ΔR) and d (az), are shown in the aforementioned application of Boles, entitled &#34;Range/Azimuth Ship Imaging For Ordnance Control&#34;, to be 
     
         d(ΔR)=-i.sub.f (d.sub.az)∫W.sub.zz dt 
    
     and 
     
         d(az)=i.sub.R (RR)∫W.sub.zz dt 
    
     where i f  and i R  are initial cursor doppler and range cell indices, W zz  dt is the rotated angle about the zz axis, and where the integration is performed on a frame-to-frame basis from the time of initial cursor placement. 
     Referring now to FIG. 6, a block diagram of the preferred embodiment of the system used for practicing the present invention is illustrated. As shown therein, pulses of electro-magnetic energy generated in Coherent Transmitter 11 from reference signals derived in Exciter/Frequency Synthesizer 12 are radiated from Transmitting Antenna 9 so as to optimally illuminate a ship target under way on the surface of the sea. Signals reflected from the ship target are received by Interferometer Antenna 10 comprised of two separate receiving elements whose common boresight direction corresponds to that of Transmitting Antenna 9. Reference signals from Exciter/Frequency Synthesizer 12, coherent with the transmitted signal, are applied to Receivers 13 and 14 to demodulate the target signals received by the interferometer elements to quadrature related I and Q doppler shifted components, which represent the real and imaginary components of the instantaneous analog signal vectors out of the receivers. These analog signal components are digitized in the A/D Converters 15 and 16 at a sampling rate determined by system range resolution requirements. These digitized samples, representing separately the signals received by the two interferometer elements, are stored in Bulk Memory 19 on a pulse-to-pulse basis. Concurrently, motion compensation corrections for antenna phase center translational and rotational motions, as well as for ship translational motion, are computed and stored in Motion Compensation Unit 20 in a time sequenced order on the basis of computations performed in General Purpose Computer 17 of the two-way line-of-sight displacement change between antenna phase centers and the tracking point of the ship, predicated on the regression solution for line-of-sight velocity obtained from Velocity Computer 28 and from Inertial Navigation System 18 measurements. 
     At the end of the coherent integration interval for which data is being collected, corrections stored in Motion Compensation Unit 20 are applied to the time sequences stored in Bulk Memory 19 in the form of vector rotations reflecting the two-way motion compensation phase correction to each range sample of each pulse of the (2) sequences representing signals received at the two interferometer halves over the coherent integration interval. After motion compensation correction, data is read out of Bulk Memory 19 (as new data is being entered) for Fourier Transform digital signal processing to produce the necessary filtering so as to provide the desired azimuth resolution in each range bin. The filtering is performed in FFT Processors 21 and 22 which perform Fast Fourier Transform digital processing so as to produce azimuth resolved coherently integrated vector sums in each doppler filter of each range bin. The outputs of FFT Processors 21 and 22 are applied to Range/Doppler Map Generator 23 for noncoherent integration on a cell-by-cell basis and hence, to Scan Converter 25 for isodop distortion removal, for summation with cursor video signals generated in Scan Converter 25, and for image intensity scaling, prior to inputting to CRT Display 26 for image presentation of the ship target. 
     The range/doppler maps generated in FFT Processors 21 and 22 are also inputted to Doppler/Azimuth Angle Processor 24 where a phase comparison is made between corresponding range/doppler cell outputs in the two maps on a cell-by-cell basis, for each and every resolution cell of the two maps, to yield the interferometric azimuth angle of the signal received in each such resolution cell, in accordance with well known interferometric principles, i.e., Δθ=sin -1  [(λ/2πd)φ], where Δθ is the azimuth spatial angle, d is the interferometer baseline spacing, and φ  is the measured electrical phase. The doppler/azimuth angle data distributions in each range bin that are formed in Doppler/Azimuth Angle Processor 24 are applied to Velocity Computer 28 wherein the constants K f , K R , and θ o , are determined using linear regression techniques in accordance with the teachings of this invention and as exemplified by equations (24) through (28), from which the three velocity terms, ε v&#39;xx , V&#39; yy , and V&#39; zz  are computed as specified by equations (21) through (23). In these latter three equations, λ, f r , N p , R.sub. o, and RR are operating constants, and η o , the depression angle to the ship focus point, is derived from on-board altitude and slant range measurements, all of which are obtained from General Purpose Computer 17. 
     The value of V&#39; zz  so determined in Velocity Computer 28 is used to remove the isodop distortion inherent in the range/doppler image stored in Scan Converter 25, by shifting each doppler cell index of each range bin, i R , equally by an amount ##EQU23## rounded to the nearest doppler cell unit. The resulting range/doppler image in Scan Converter 25, which now represents the true range/azimuth distribution of the ship image, is inputted to CRT Display 26. Concurrently, the 3 velocity terms derived from the regression solutions in Velocity Computer 28, ε v&#39;xx , V&#39; yy  and V&#39; zz  are inputted to General Purpose Computer 17. 
     The value of line-of-sight velocity error, ε v&#39;xx  is used as a line-of-sight velocity correction, -ε v&#39;xx , in Motion Compensation Unit 20 for the purpose of focusing and image centering by the application of Motion Compensation phase corrections as described earlier. 
     The value of V&#39; yy  is used by General Purpose Computer 17 for the formation of doppler filter bandwidth and spacing, integration time, and pulse repetition frequency, as prescribed by equations (31) and (32), to be used in the formation of the next synthetic aperture. 
     Antenna azimuth steering commands are generated in General Purpose Computer 17 on the basis of aircraft navigational data inputs from the Inertial Navigation System 18 and averaged interferometric azimuth angle measurements from the various ship scatterers provided by Doppler/Azimuth Angle Processor 24, in such a manner as to maintain the average of all interferometric azimuth angles at zero (with respect to its electrical boresight) when smoothed with a tracking filter in General Purpose Computer 17 over successive apertures. The net pulse-to-pulse two-way line-of-sight displacement change due to relative translational motion between aircraft and ship computed in General Purpose Computer 17 for purposes of motion compensation, also serves as a basis for controlling the precise timing for the start of range sampling at A/D Converters 15 and 16, so that corresponding samples from pulse-to-pulse over the coherent integration interval represent the same range increment of the ship. The net line-of-sight displacement change with time, in conjunction with aircraft altitude available in General Purpose Computer 17 from other on-board sensors, is used to provide elevation steering to the Transmitting antenna 9 and receiving Interferometer Antenna 10. 
     Cursor placement to a desired resolution cell of the ship image is achieved by operator designation of video signals through General Purpose Computer 17 and applied to Scan Converter 25. Cursor tracking of the designated ship target resolution cell during the weapon delivery phase is similarly controlled by General Purpose Computer 17 on the basis of its computations of target cell range and azimuth angle change with time derived from its navigational solutions for translational motions between aircraft and ship, wherein said cursor video signals are injected into appropriate range/doppler resolution cells as time progresses. Second order corrections to cursor tracking of the designated ship target cell to eliminate potential tracking errors arising from ship rotational (yaw) motion are effected as described in the aforementioned application of Boles, entitled &#34;Range/Azimuth Ship Imaging For Ordnance Control&#34;. 
     Control signals for the pointing of Antennae 9 and 10 and for pulse repetition frequency (prf) control of coherent Transmitter 11 are obtained from General Purpose Computer 17. All task management such as data transmission and initiation of subroutine sequences, are performed by the Radar Data Processor 27. 
     Although separate transmitting and receiving antennae have been illustrated in the block diagram if FIG. 6, both transmitting and receiving functions could be combined into a single aperture. All elements in the block diagram can be either manufactured or purchased as separate commercial entities and incorporated into the embodiments of this invention by those skilled in the art. Also, many of the functional units shown separately in FIG. 6 can be combined for purposes of design economy. 
     Thus, although the invention has been described and illustrated in detail, it is to be understood that the same is exemplary and is not to be taken by way of limitation, the spirit and scope of this invention being limited only by the terms of the appended claims.