Method of image generation by means of two-dimensional data processing in connection with a radar with synthetic aperture

• US 5,627,543 A
• Filed: 08/03/1995
• Issued: 05/06/1997
• Est. Priority Date: 08/05/1994
• Status: Expired due to Term

• Current Assignee: Deutsche Telekom AG
• Original Assignee: Deutsche Telekom AG
• Inventors: Moreira, Alberto
• Primary Examiner(s): Lobo, Ian J.
• Application Number:
  US08/510,698
• Time in Patent Office: 642 Days
• Field of Search: 342/25, 342/195, 342/196, 342/191
• US Class Current: 342/25.00A
• CPC Class Codes: G01S 13/9011 with frequency domain processing of the SAR signals in azimuth G01S13/9005 takes precedence
• CPC Class Codes: G01S 7/4004 of parts of a radar system

Abstract

In connection with a method for image generation by means of two-dimensional data processing, received SAR data are multiplied by a phase correction (H_mc) for a reference range (r_ref) for the insertion of a motion compensation and for processing at a high drift angle, and an additional cubic phase term is inserted for compensating a range migration. The entire range migration is then eliminated by means of an additional linear frequency displacement; subsequently the SAR data are transformed back into the "range-Doppler" domain. A remaining phase error, created by a "chirp scaling" correction, is corrected, the SAR data are transformed back into the time domain and a phase correction as a function of the range is performed by multiplication for the exact motion compensation in the time domain. The one-dimensional reference function is performed in the frequency domain for azimuth compression, by means of which two-dimensional SAR data are obtained.

Claims

1. A method for image generation by means of two-dimensional data processing in connection with a radar with synthetic aperture SAR, characterized in that, for the insertion of a motion compensation and for processing at a large drift angle, received SAR data (s(τ

, t;

r_o)),

wherein the time in a range direction is identified by τ

, the time in the azimuth direction by t and the shortest range between the target subject and the radar system by r_o,

are multiplied by a phase correction;

space="preserve" listing-type="equation">H.sub.mc (τ

, t;

r.sub.ref)=exp [j·

φ

.sub.mc (τ

, t;

r.sub.ref)]

where j describes a complex part and φ

_mc is a line-of-sight phase formation compensation for a reference range (r_ref) in a first multiplication unit;

an additional cubic phase term is inserted via a second multiplication unit following the performance of azimuth FFTs for compensating a range migration by means of a function;

##EQU12## wherein f_a is the azimuth frequency, k is a modulation rate as a function of the range and azimuth frequency, a'"'"' is a linear scaling factor, R'"'"' the range migration, k_r the modulation rate of a transmitted chirp signal, λ

the wave length, c the speed of light and β

² and β

³ are respectively a square and a cube of a parameter β

defined as ##EQU13## wherein v is a platform speed;

following a performance of range FFTs after compensation of the range migration by means of the function H₁, the range migration is entirely eliminated by means of an additional linear frequency displacement (f_r) in a function ##EQU14## by multiplication in a third multiplication unit;

subsequently the SAR data are transformed back into a two-dimensional range-Doppler domain having coordinates of range time and azimuth frequency, by means of range IFFTs;

a remaining phase flow created by a chirp scaling function is corrected by multiplication in a fourth multiplication unit by a function;

##EQU15## subsequently the SAR data are transformed back into a time domain by means of azimuth IFFTs;

then a range-dependent phase correction

space="preserve" listing-type="equation">H.sub.mc (τ

, t;

r.sub.o)=exp[j·

(φ

.sub.mc (τ

, t;

r.sub.o)-φ

.sub.mc (τ

, t;

r.sub.ref))]

is performed by multiplication in a fifth multiplication unit for the exact motion compensation in the time domain;

then, after performing further azimuth FFTs for azimuth compression, a one-dimensional reference function in the frequency domain ##EQU16## is performed by multiplication in a sixth multiplication unit, and following the performance of azimuth IFFTs, two-dimensional SAR data (corresponding to a function f(τ

, t)) are obtained.

2. The method in accordance with claim 1, characterized in that, for processing with a variable Doppler centroid in a range domain following the performance of the first azimuth FFTs,

a spectral extension of a azimuth spectrum, by means of which the point number after each azimuth FFT is artificially increased is performed in the azimuth direction by determining in accordance with ##EQU17## of an extended frequency variation ##EQU18## wherein a Doppler centroid is identified by f_dc and the pulse repetition frequency by PRF;

the frequency variation of the functions H₁, H₂ and H₃ is increased corresponding to the frequency extension in the azimuth direction;

following a multiplication in the fourth multiplication unit by the function H₃ spectral extension in the azimuth direction is removed, and

during azimuth processing the function H₄ is generated for each range position by multiplication in the sixth multiplication unit by an appropriated value of the Doppler centroid f_dc ##EQU19## wherein φ

_i is the angle of incidence of the radar beams on the ground, φ

_d is the drift angle of the platform and φ

_n the pitch angle of the platform.

3. The method in accordance with claim 1, characterized in that for processing with a variable Doppler centroid in the azimuth direction

a transformation of the SAR data (s(τ

, t;

r_o)) from the time domain to the range-Doppler domain is performed by means of short azimuth FFTs, so that a plurality of sub-apertures in the azimuth direction are created;

a point number in the azimuth direction of the functions H₁, H₂, and H₃ is reduced corresponding to the length of the azimuth FFTs, and

the sub-apertures are transformed back into the time domain by means of short azimuth IFFTs and are subsequently combined for forming the entire aperture.

4. The method in accordance with claim 3, characterized in that for inserting a motion extraction method following the phase correction by means of the function H₃

the azimuth spectra of each sub-aperture are evaluated via the fourth multiplication unit in accordance with the RDM (Reflectivity Displacement Method) which performs a cross correlation of adjacent azimuth spectra, and

the correction phase H_mc (τ

, t;

r_o) is multiplied in the fifth multiplication unit for the motion compensation with the SAR signal after the combination of the sub-apertures in the time domain.