METHOD AND DEVICE FOR PROCESSING SAR RAW DATA

Abstract

A method according to the present invention comprises the steps of: dividing SAR raw data into one or more sub-aperture data by a predetermined number in an azimuth direction; performing a spectral length extension FFT on the sub-aperture data in the azimuth direction; multiplying the sub-aperture data by a chirp scaling function; performing a range FFT on the sub-aperture data; performing range compression, SRC, and a bulk RCMC on the sub-aperture data; performing an inverse chirp-z transform on the sub-aperture data in a range direction; multiplying the divided sub-aperture data by a predetermined first function; performing an IFFT on the sub-aperture data in the azimuth direction; recombining the sub-aperture data; multiplying the recombined data by a second function and deramping same; performing an azimuth FFT on the recombined data; performing an azimuth IFFT on the recombined data; multiplying the recombined data by a third function and deramping same; performing the azimuth FFT on the recombined data; performing azimuth compression by multiplying the recombined data by a fourth function; performing an azimuth inverse chirp-z transform on the recombined data; and multiplying the recombined data by a fifth function for phase preservation.

Claims

1. A method for processing synthetic aperture radar image (SAR) raw data, the method comprising: dividing the SAR raw data into one or more sub-aperture data by a predetermined number in an azimuth direction; performing a spectral length extension fast Fourier transform (FFT) on the divided sub-aperture data in the azimuth direction; multiplying the divided sub-aperture data by a chirp scaling function; performing the FFT on the divided sub-aperture data in a range direction; performing a range compression, a secondary range compression (SRC), and a bulk range cell migration correction (RCMC) on the divided sub-aperture data; performing an inverse chirp-z transform on the divided sub-aperture data in the range direction; multiplying the divided sub-aperture data by a first function predetermined for residual phase correction and azimuth scaling; performing an inverse fast Fourier transform (IFFT) on the divided sub-aperture data in the azimuth direction; recombining the divided sub-aperture data; multiplying the recombined data by a second function to perform deramping; performing the FFT on the recombined data in the azimuth direction; performing the IFFT on the recombined data in the azimuth direction; multiplying the recombined data by a third function to perform deramping; performing the FFT on the recombined data in the azimuth direction; performing an azimuth compression by multiplying the recombined data by a fourth function; performing the inverse chirp-z transform in an azimuth direction on the recombined data; and multiplying the recombined data by a fifth function for phase preservation.

2. (canceled)

3. (canceled)

4. (canceled)

5. (canceled)

6. The method of claim 1, wherein the inverse chirp-z transform is performed by Equation 1 below: $x_k=\text{?}.Math.\left[\mathrm{FFT}\left\{\mathrm{IFFT}\left\{Y_n\right\}.Math.\mathrm{IFFT}\left\{\text{?}\right\}\right\}\right],k=0,1,\mathrm{.Math.},M-1$ $Y_n=X\left(2_n\right).Math.B^n.Math.\text{?}.Math.n=\frac{F_n}{\Delta F}$ $B=B_0.Math.\exp \left(j2{\mathrm{\pi \theta }}_0\right).Math.{\theta }_0=\Delta F.Math.t_0$ $W=W_0.Math.\exp \left(\mathrm{j2}{\mathrm{\pi \varphi }}_0\right).Math.{\varphi }_0=\Delta F.Math.\Delta t$ $\text{?}\text{indicates text missing or illegible when filed}$ where, X(z.sub.n) is a signal in the input frequency domain, M is the number of output sample signals, ΔF is frequency spacings of the input spectrum signals, and B.sub.0 and W.sub.0 are amplitude constants, and the start time of the signal x.sub.k on the output time is set to t.sub.0, and the time spacing of the samples is set to Δt.

7. The method of claim 1, wherein the first function is defined by Equation 2, the second function by Equation 3, the third function by Equation 4, the fourth function by Equation 5, and the fifth function by Equation 6 as follows: $\begin{array}{cc}{H}_4\left(f_a,r\right)=M_1\left(w_{\eta }\right).Math.\exp \left[-j\frac{\pi }{K_{\mathrm{scl}}\left(r\right)}{f}_a^2\right]& \left[\mathrm{Equation}2\right]\end{array}$ where, $K_{\mathrm{scl}}\left(r\right)=-\frac{2v_{\mathrm{eff}}^2\left(r\right)}{\lambda r_{\mathrm{scl}}\left(r\right)},$
r.sub.scl(r)=r,
M.sub.1(w.sub.η)=exp[j{2(2π/λ+w.sub.r/c)R.sub.r2(η*)+w.sub.ηη*}],
R.sub.r2(η)=c.sub.4η.sup.4+(c.sub.3+4c.sub.4t.sub.3)η.sup.3+(c.sub.2+3c.sub.3t.sub.1+6c.sub.4t.sub.1.sup.2)η.sup.2 where, λ is the wavelength with respect to the center frequency of the transmission signal forming the beam, r is the closest approach range, f.sub.a is the azimuth Frequency (Doppler Frequency shift), v.sub.eff is the effective velocity, K.sub.scl(r) is a scaling Doppler rate, and c.sub.2, c.sub.3 and c.sub.4 are coefficients obtained from geometry including orbit information and attitude information of a satellite;
H.sub.5(t.sub.a,r)=exp[−jπK.sub.rot_geometry.Math.(t.sub.a−t.sub.mid).sup.2]  [Equation 3] where, t.sub.a is the azimuth time, t.sub.mid is the azimuth time of the selected derotation center, and r.sub.rot_geometry is the distance to the beam rotation center given by geometry, $K_{\mathrm{rot}1}=K_{\mathrm{rot}\_\mathrm{geometry}}=-\frac{2v_{\mathrm{eff}}^2\left(r\right)}{\lambda r_{\mathrm{ro}t\_\mathrm{geomtry}}}$ is an azimuth derotation Doppler rate;
H.sub.6(t.sub.a,r)=exp[−jπ(K.sub.rot2(r)−K.sub.rot_geometry).Math.(t.sub.a−t.sub.mid).sup.2]  [Equation 4] where $K_{\mathrm{rot}2}\left(r\right)=-\frac{2v_{\mathrm{eff}}^2\left(r\right)}{\lambda r_{\mathrm{rot}2}\left(r\right)}$ is an azimuth deramping Doppler rate, and
r.sub.rot2(r)=r.Math.ε; $\begin{array}{cc}{H}_7\left(f_a,r\right)=W\left(f_a\right).Math.\exp \left[j\frac{\pi }{K_{\mathrm{eff}}\left(r\right)}{f}_a^2\right]\mathrm{where},K_{\mathrm{eff}}\left(r\right)=K_{\mathrm{act}}\left(r\right)-K_{\mathrm{rot}2}\left(r\right);\mathrm{and}& \left[\mathrm{Equation}5\right]\\ H_8\left(t_a,r\right)=\exp \left[j\pi K_t\left(r\right).Math.{\left(1-\frac{1}{\mathrm{.Math.}}\right)}^2.Math.{\left(t_a-t_{\mathrm{mid}}\right)}^2\right]\mathrm{where},K_t\left(r\right)=-\frac{2v_{\mathrm{eff}}^2\left(r\right)}{\lambda .Math.r\left(\mathrm{.Math.}-1\right)}.& \left[\mathrm{Equation}6\right]\end{array}$

8. The method of claim 7, wherein, when the SAR operational mode is Stripmap, ScanSAR and TOPS mode, $\mathrm{.Math.}=\frac{r_{\mathrm{rot}\_\mathrm{geometry}}}{r_{\mathrm{mid}}}$ is set, and, when the SAR operational mode is Stripmap or ScanSAR, if the value of r.sub.rot_geometry cannot be calculated numerically, r.sub.rot_geometry=1000.Math.r.sub.mid, if |r.sub.rot_geometry|>1000.Math.r.sub.mid and r.sub.rot_geometry>0, r.sub.rot_geometry=1000.Math.r.sub.mid, if |r.sub.rot_geometry|>1000.Math.r.sub.mid and r.sub.rot_geometry<0, r.sub.rot_geometry=−1000 r.sub.mid, and if |r.sub.rot_geometry|≤1000.Math.r.sub.mid, the calculated value of r.sub.rot_geometry is applied as it is, where r.sub.mid is the closest distance to the center of a scene.

9. The method of claim 7, when a SAR operational mode is a spotlight operational mode including a sliding spotlight and a staring spotlight, ε is set by the following equation: ${\mathrm{.Math.}}_{\mathrm{optimized}}=\frac{{\mathrm{.Math.}}_{{\min \mathrm{\_\gamma }}_1,{\gamma }_2}+{\mathrm{.Math.}}_{{\max \mathrm{\_\gamma }}_1,{\gamma }_2}}{2}$ where ε.sub.min_γ.sub.1.sub.,γ.sub.2, ε.sub.max_γ.sub.1.sub.,γ.sub.2 are defined as an upper limit and a lower limit of ε that satisfy both Conditional Equation 1 and Conditional Equation 2 below: [Conditional Equation 1] ε.sub.min_γ.sub.1<ε<ε.sub.max_γ.sub.1, wherein, Conditional Equation 1 is calculated as ε range that is calculated for values of γ.sub.1 range such as $\frac{\Delta t_{a0}}{T_a}<{\gamma }_1<1.25\mathrm{in}\mathrm{.Math.}=1+\frac{\Delta t_{\mathrm{ao}}}{{\gamma }_1{T}_a-\Delta t_{\mathrm{ao}}},{\gamma }_1>0;$ and [Conditional Equation 2] ε.sub.min_γ.sub.1<ε<ε.sub.max_γ.sub.2, wherein Conditional Equation 2 is calculated as ε range that is calculated for values of γ.sub.2 range such as $1<{\gamma }_2<0.75.Math.\frac{B_{a\_\mathrm{Targe}t}}{B_{\mathrm{FOV}}}\mathrm{in}\mathrm{.Math.}=\frac{v_{\mathrm{eff}}\left(r_{\mathrm{mid}}\right).Math.T_{\mathrm{obs}}}{v_{\mathrm{eff}}\left(r_{\mathrm{mid}}\right).Math.T_{\mathrm{obs}}-{\gamma }_2.Math.{\theta }_{\mathrm{az}}.Math.r_{\mathrm{mid}}},{\gamma }_2>0,\mathrm{where},B_{a\_\mathrm{Target}}=\frac{2v_{\mathrm{eff}}^2\left(r_{\mathrm{mid}}\right)}{\lambda r_{\mathrm{mid}}}.Math.T_{\mathrm{obs}}$ and T.sub.obs is the target observation duration (Target dwell time).

10. A computer-readable recording medium storing a program for performing a method for processing a synthetic aperture radar image (SAR) raw data, the program comprising: an instruction set for dividing the SAR raw data into one or more sub-aperture data by a predetermined number in an azimuth direction; an instruction set for performing a spectral length extension fast Fourier transform (FFT) on the divided sub-aperture data in the azimuth direction; an instruction set for multiplying the divided sub-aperture data by a chirp scaling function; an instruction set for performing the FFT on the divided sub-aperture data in a range direction; an instruction set for performing a range compression, a secondary range compression (SRC), and a bulk range cell migration correction (RCMC) on the divided sub-aperture data; an instruction set for performing an inverse chirp-z transform on the divided sub-aperture data in the range direction; an instruction set for multiplying the divided sub-aperture data by a first function predetermined for residual phase correction and azimuth scaling; an instruction set for performing an inverse fast Fourier transform (IFFT) on the divided sub-aperture data in the azimuth direction; an instruction set for recombining the divided sub-aperture data; an instruction set for multiplying the recombined data by a second function to perform deramping; an instruction set for performing the FFT on the recombined data in the azimuth direction; an instruction set for performing the IFFT on the recombined data in the azimuth direction; an instruction set for multiplying the recombined data by a third function to perform deramping; an instruction set for performing the FFT on the recombined data in the azimuth direction; an instruction set for performing an azimuth compression by multiplying the recombined data by a fourth function; an instruction set for performing the inverse chirp-z transform in an azimuth direction on the recombined data; and an instruction set for multiplying the recombined data by a fifth function for phase preservation.

Description

BRIEF DESCRIPTION OF THE DRAWING

[0053] FIG. 1 is a diagram illustrating an image formation process of a SARP core algorithm.

[0054] FIG. 2 is an example of representing SARP focusing by mathematical modeling.

[0055] FIG. 3 is a diagram illustrating a focusing process of the SARP core algorithm.

[0056] FIG. 4 is a diagram showing a process of Sliding spotlight and TOPS SAR data processing using the conventional BAS.

[0057] FIG. 5 is a diagram provided to explain a SARP core algorithm according to an embodiment of the present invention.

DETAILED DESCRIPTION

[0058] Hereinafter, the exemplary embodiments of the present invention will be described in detail with reference to the accompanying drawings for those with ordinary knowledge in the art to be able to easily achieve the present invention.

[0059] FIG. 5 is a diagram provided to explain a method for processing SAR raw data by a SARP core algorithm according to an embodiment of the present invention.

[0060] Referring to FIG. 5, first, a processor for processing the SAR raw data according to the present invention divides the synthetic aperture radar image (SAR) raw data into one or more sub-aperture data by a predetermined number in an azimuth direction, at S1.

[0061] When dividing and processing the SAR raw data based on a sub-aperture unit in the azimuth direction at S1, it is possible that the azimuth time length of one sub-aperture may be set as any value without limitations within the azimuth time length of the entire raw data. This is because the method for converting the data into the azimuth frequency domain while satisfying the Nyquist criteria is not a short azimuth FFT but a spectral length extension azimuth FFT. The number of sub-apertures may be set in consideration of the beam or operational mode of the designed SAR system. At S1, the number of sub-apertures may be set to a minimum value of 1 or 2 or more. The azimuth time length of the sub-aperture is determined according to the set number of sub-apertures.

[0062] When the short azimuth FFT is used, the processing accuracy of the SARP may be lowered for a certain SAR system. This is because the algorithm of the SARP is developed in the principle of stationary phase (POSP) manner and thus there is an error in approximation, and the error is increased as the time bandwidth product (TBP) value is decreased, and the short FFT decreases the TBP value. In addition, the short azimuth FFT method increases the number of sub-apertures, resulting in lowered quality of the image in a part where the sub-aperture images are recombined. On the other hand, the SARP core algorithm according to the present invention minimizes the number of sub-apertures while allowing adjustment, thereby ensuring maximum accuracy in the process of processing the signal for any SAR system. The shortcoming of this method can be that the processing time is increased due to the additional spectral length extension process, but this is acceptable when considering that the method can be applied to the entire operational modes.

[0063] Next, at S2, the spectral length extension azimuth FFT may be performed on the sub-aperture data divided at S1. Through S2, the sub-aperture data may be converted from the SAR signal domain into the azimuth frequency domain.

[0064] After performing S2, chirp scaling may be performed by multiplying the sub-aperture data converted into the frequency domain by a chirp scaling function (H.sub.cs), at S3.

[0065] After performing S3, a range FFT may be performed on the divided sub-aperture data in a range direction, at S4. The sub-aperture data divided through S4 may be converted into a two-dimensional frequency domain.

[0066] After performing S4, by multiplying the sub-aperture data converted into the two-dimensional frequency domain by a function H.sub.RC×H.sub.RC×H.sub.BV at S5, the range compression, the secondary range compression (SRC), and the bulk range cell migration correction (RCMC) may be performed on the divided sub-aperture data. H.sub.RC is a function to perform the range compression, H.sub.RC is a function to perform the SRC, and H.sub.BV is a function to perform the bulk RCMC. According to an embodiment, the range compression, the secondary range compression (SRC), and the bulk range cell migration correction (RCMC) may be performed.

[0067] First, H.sub.cs, H.sub.RC, H.sub.RC and H.sub.BV may use a function corresponding to any generally known SARP core algorithm. For example, the functions proposed in “Ian G. Cumming, Frank H. Wong, Digital Processing of Synthetic Aperture Radar Data, Artech House Inc., pp. 283-322, 2005.” may be used, the details of which are well known to a person skilled in the art and thus will not be redundantly described in detail herein.

[0068] Next, after performing S5, an inverse chirp-z transform may be performed on the divided sub-aperture data in the range direction, at S6. As a method for forming an image at S6, inverse chirp-z transformation (ICZT) may be used instead of the inverse FFT (IFFT). The method of the IFFT sets the sample spacings of an image by adjusting the number of samples of data in the frequency domain. However, in the IFFT method, since the number of samples is integer and can only be set discontinuously, there is a constraint that the sample spacings are also adjusted discontinuously. In the case of the SAR operational mode in which observation is performed with multiple beams, there is a high possibility that the sample spacings between formed images of the beams differs even by a very small value. When such images of multiple beams are simply mosaicked, the quality of the entire scene image is lowered. As the image size becomes larger, the distortion of the location information becomes greater. When the interpolation function is applied to solve this problem, the quality of the image is lowered and the processing time is increased due to the accuracy error of the interpolation itself. Therefore, according to an embodiment, by using the inverse chirp-z transform method capable of continuously adjusting the sample spacings, the accuracy of the positions or values of the pixels during image formation may be secured to the maximum for any SAR system.

[0069] At S6, the inverse chirp-z transform may be performed by Equation 1 below, and may perform a function of inverse transforming any input signal in the frequency domain into a signal in the time domain.

[00013]$\begin{array}{cc}{x}_k=W^{\frac{k^2}{2}}.Math.\left[\mathrm{FFT}\left\{\mathrm{IFFT}\left\{Y_n\right\}.Math.\mathrm{IFFT}\left\{W^{-\frac{n^2}{2}}\right\}\right\}\right],k=0,1,\mathrm{.Math.},M-1Y_a=X\left(z_n\right).Math.B^n.Math.W^{\frac{n^2}{2}},n=\frac{F_n}{\Delta F}B=B_0.Math.\exp \left(j2\pi {\theta }_0\right).Math.{\theta }_0=\Delta F.Math.t_{0}W=W_0.Math.\exp \left(j2{\mathrm{\pi \varphi }}_0\right).Math.{\varphi }_0=\Delta F.Math.\Delta t& \left[\mathrm{Equation}1\right]\end{array}$

[0070] where, X(z.sub.n) is a signal in the input frequency domain, M is the number of output sample signals, ΔF is frequency spacings of the input spectrum signals, and B.sub.0 and W.sub.0 are amplitude constants. The start time of the signal x.sub.k on the output time is set to t.sub.0, and the time spacing of the samples is set to Δt.

[0071] Next, at S7, with respect to the inverse chirp-z transformed sub-aperture data at S6, the predetermined first function H.sub.4 for the residual phase correction and azimuth scaling may be multiplied by the sub-aperture data that is inverse chirp-z transformed in the range direction at S6.

[0072] At S7, accurate azimuth matched filtering is performed on the sub-aperture data inverse chirp-z transformed in the range direction by the first function H.sub.4, and a quadratic phase signal is formed by using K.sub.scl(r) corresponding to the actual azimuth Doppler rate component.

[0073] Next, at S8, the inverse fast Fourier transform (IFFT) may be performed on the divided sub-aperture data in the azimuth direction. Through S8, the sub-aperture data may be converted from the range Doppler domain into the SAR signal domain.

[0074] After S1 to S8 are all performed on the sub-aperture data divided from the SAR raw data, the divided sub-aperture data may be recombined at S9. When the SAR raw data is processed as one sub-aperture data at S1, S1 to S8 may be performed only once.

[0075] Thereafter, the data recombined at S9 (hereinafter, “recombined data”) may be multiplied by a second function H.sub.5 to perform deramping at S10.

[0076] Next, at S11, FFT may be performed on the deramped recombined data at S10 in the azimuth direction, and azimuth antenna pattern compensation may be performed.

[0077] Next, at S12, IFFT may be performed in the azimuth direction on the FFT-processed recombined data in the azimuth direction at S11.

[0078] Thereafter, the IFFT-processed recombined data in the azimuth direction at S12 may be multiplied by a third function H.sub.6 to perform deramping at S13.

[0079] Then, at S14, the FFT may be performed in the azimuth direction on the deramped recombined data at S13.

[0080] Next, at S15, the azimuth compression (AC) may be performed by multiplying the FFT-processed recombined data in the azimuth direction at S14 by a fourth function H.sub.7.

[0081] Thereafter, at S16, the inverse chirp-z transform in the azimuth direction may be performed on the recombined data compressed in the azimuth direction at S16.

[0082] Finally, at S17, the inverse chirp-z transformed recombined data in the azimuth direction at S16 may be multiplied by a fifth function H.sub.8 for phase preservation to generate single look complex (SLC) data.

[0083] Hereinbelow, the first function H.sub.4, the second function H.sub.5, the third function H.sub.6, the fourth function H.sub.7, and the fifth function H.sub.8 used in the SARP core algorithm according to the present invention and the constraints to accurately process the SAR raw data in every SAR operational mode are described in detail.

[0084] The first function H.sub.4 is defined by Equation 2, the second function H.sub.5 is defined by Equation 3, the third function H.sub.6 is defined by Equation 4, the fourth function H.sub.7 is defined by Equation 5, and the fifth function H.sub.8 is defined by Equation 6.

[00014]$\begin{array}{cc}{H}_4\left(f_a,r\right)=M_1\left(w_{\eta }\right).Math.\exp \left[-j\frac{\pi }{K_{\mathrm{scl}}\left(r\right)}{f}_a^2\right]& \left[\mathrm{Equation}2\right]\end{array}$

[0085] where,

[00015]$K_{\mathrm{scl}}\left(r\right)=-\frac{2v_{\mathrm{eff}}^2\left(r\right)}{\lambda r_{\mathrm{scl}}\left(r\right)},$
r.sub.scl(r)=r

M.sub.1(w.sub.η)=exp[j{2(2π/λ+w.sub.r/c)R.sub.r2(η*)+w.sub.ηη*}],

R.sub.r2(η)=c.sub.4η.sup.4+(c.sub.3+4c.sub.4t.sub.3)η.sup.3+(c.sub.2+3c.sub.3t.sub.1+6c.sub.4t.sub.1.sup.2)η.sup.2

[0086] where, λ is the wavelength with respect to the center frequency of the transmission signal forming the beam, r is the closest approach range, f.sub.a is the azimuth Frequency (Doppler Frequency shift), v.sub.eff is the effective velocity, K.sub.scl(r) is a scaling Doppler rate, and c.sub.2, c.sub.3 and c.sub.4 are coefficients. According to an embodiment, c.sub.2, c.sub.3 and c.sub.4 are obtained by using geometry including orbit information and attitude information of a satellite.

H.sub.5(t.sub.a,r)=exp[−jπK.sub.rot_geometry.Math.(t.sub.a−t.sub.mid).sup.2]  [Equation 3]

[0087] where, t.sub.a is the azimuth time, t.sub.mid is the azimuth time of the selected derotation center, and r.sub.rot_geometry is the distance to the beam rotation center given by geometry.

[00016]$K_{\mathrm{rot}1}=K_{\mathrm{rot\_geometry}}=-\frac{2v_{\mathrm{eff}}^2\left(r\right)}{\lambda r_{\mathrm{rot\_geomtry}}}$

is an azimuth derotation Doppler rate.

H.sub.6(t.sub.a,r)=exp[−jπ(K.sub.rot2(r)−K.sub.rot_geometry).Math.(t.sub.a−t.sub.mid).sup.2]  [Equation 4]

[00017]$K_{\mathrm{rot}2}\left(r\right)=-\frac{2v_{\mathrm{eff}}^2\left(r\right)}{\lambda r_{\mathrm{rot}2}\left(r\right)}$

[0088] is the azimuth deramping Doppler rate, and r.sub.rot2(r)=r.Math.ε.

[00018]$\begin{array}{cc}{H}_7\left(f_a,r\right)=W\left(f_a\right).Math.\exp \left[j\frac{\pi }{K_{\mathrm{eff}}\left(r\right)}{f}_a^2\right]\mathrm{where},K_{\mathrm{eff}}\left(r\right)=K_{\mathrm{scl}}\left(r\right)-K_{\mathrm{rot}2}\left(r\right)& \left[\mathrm{Equation}5\right]\\ H_8\left(t_a,r\right)=\exp \left[j\pi K_t\left(r\right).Math.{\left(1-\frac{1}{\mathrm{.Math.}}\right)}^2.Math.{\left(t_a-t_{\mathrm{mid}}\right)}^2\right]\mathrm{where},K_t\left(r\right)=-\frac{2v_{\mathrm{eff}}^2\left(r\right)}{\lambda .Math.r\left(\mathrm{.Math.}-1\right)}& \left[\mathrm{Equation}6\right]\end{array}$

[0089] According to an embodiment, by adjusting ε, processing the SAR raw data using the first function H.sub.4, the second function H.sub.5, the third function H.sub.6, the fourth function H.sub.7, and the fifth function H.sub.8 can be applied to all modes of the SAR system, i.e., stripmap, ScanSAR, TOPS, and sliding spotlight, staring spotlight, and any operational mode between sliding spotlight and staring spotlight.

[0090] Hereinbelow, a method for adjusting E will be described in detail.

[0091] Among the five constraints described with reference to the conventional BAS system, some constraints are not applicable in the SARP core algorithm according to the present invention. Specifically, since the SARP core algorithm of the present invention performs the azimuth spectral length extension FFT for each sub-aperture, the condition that the azimuth bandwidth after derotation should satisfy the Nyquist criteria, which is the first constraint of the conventional BAS, is not applicable. The constraint in which the Doppler rate after the second derotation should have a proper large value other than ‘0’ is applicable. In addition, since the pixel spacings may be freely adjusted by using ICZT instead of IFFT for image formation in the SARP core algorithm according to the present invention, the condition that the azimuth pixel spacings of the processed image should not be too small compared to the azimuth resolution, which is the third constraint of the conventional BAS, is not applicable. The fourth constraint in which the azimuth time range of the azimuth scene should not be increased too much is not applicable. Since the SARP core algorithm of the present invention leaves a quadratic component of the actual azimuth signal when performing azimuth scaling, the condition that the time shift of the azimuth signal by H.sub.4 should not be too large, which is the fifth constraint of the BAS, is not applicable.

[0092] Therefore, the constraints of the present invention may be simplified as follows.

[0093] Firstly, the Doppler rate after the derotation should have a proper large value other than ‘0’. Secondly, the azimuth time range of the azimuth scene should not be increased too much.

[0094] The azimuth bandwidth after deramping may be represented by Equation 7 below.

[00019]$\begin{array}{cc}{B}_{\mathrm{a\_Total}\mathrm{\_Span}}=\frac{2v_{\mathrm{eff}}^2\left(r\right)}{r.Math.\lambda }.Math.\left\{.Math.\left(\frac{1}{\mathrm{.Math.}}-\frac{r}{r_{\mathrm{rot\_geometry}}}\right).Math.T_a.Math.+.Math.\left(1-\frac{1}{\mathrm{.Math.}}\right).Math.T_{\mathrm{ab}s}.Math.+B_{\mathrm{FOV}}\right\}& \left[\mathrm{Equation}7\right]\end{array}$

[0095] B.sub.FOV: Doppler frequency range of the data in instant field of view

[00020]$B_{\mathrm{FOV}}=\frac{2.Math.v_a.Math.{\theta }_{\mathrm{ax}}}{\lambda }$

[0096] The first constraint according to the present invention described above requires that the second component

[00021]$\frac{2v_{\mathrm{eff}}^2\left(r\right)}{r.Math.\lambda }.Math..Math.\left(1-\frac{1}{\mathrm{.Math.}}\right).Math.T_{\mathrm{ab}s}.Math.$

in Equation 7 have a proper value greater than 0.

[0097] Further, the second constraint requires that the following conditions be met.

[00022]$.Math.\frac{\mathrm{.Math.}}{\mathrm{.Math.}-1}.Math..Math.\Delta t_{a0}<{\gamma }_1.Math.T_a$

[0098] T.sub.a: total observation time

[0099] Δt.sub.a0: azimuth time length for scene size

[0100] γ.sub.1>0γ.sub.1: ratio of azimuth time range compared to Ta, on which scene would appear after application of H.sub.7 and azimuth ICZT.

[0101] For example, γ=1.25 may be set. However, in the case of ScanSAR mode and TOPS mode, the second constraint may not be applied for the operation purpose and characteristics.

[0102] Meanwhile, in order to satisfy the above two constraints, the SARP core algorithm according to the present invention may set c appropriately for each operational mode.

[0103] When the SAR operational mode is Stripmap, ScanSAR and TOPS mode,

[00023]$\mathrm{.Math.}=\frac{r_{\mathrm{rot\_geometry}}}{r_{\mathrm{mid}}}$

may be set.

[0104] However, when the SAR operational mode is Stripmap or ScanSAR, the following may be set.

[0105] If the value of r.sub.rot_geometry cannot be calculated numerically, r.sub.rot_geometry=1000.Math.r.sub.mid,

[0106] if |r.sub.rot_geometry|>1000.Math.r.sub.mid and r.sub.rot_geometry>0, r.sub.rot_geometry=1000.Math.r.sub.mid,

[0107] if |.sub.rot_geometry|>1000.Math.r.sub.mid and r.sub.rot_geometry<0, r.sub.rot_geometry=−1000.Math.r.sub.mid,

[0108] and if |r.sub.rot_geometry|≤1000.Math.r.sub.mid, the calculated r.sub.rot_geometry value can be applied as is. Here, r.sub.mid is the closest distance to the center of the scene.

[0109] Meanwhile, when the SAR operational mode is the spotlight operational mode including the sliding spotlight and the staring spotlight, c may be set by the following equation.

[00024]${\mathrm{.Math.}}_{\mathrm{optimized}}=\frac{{\mathrm{.Math.}}_{{\min \mathrm{\_\gamma }}_1,{\gamma }_2}+{\mathrm{.Math.}}_{{\max \mathrm{\_\gamma }}_1,{\gamma }_2}}{2}$

[0110] ε.sub.min_γ.sub.1.sub.,γ.sub.2, e.sub.max_γ.sub.1.sub.,γ.sub.2 may be defined as an upper limit and a lower limit of ε that satisfy both Conditional Equation 1 and Conditional Equation 2 below.

ε.sub.min_γ.sub.1<ε<ε.sub.max_γ.sub.1  [Conditional Equation 1]

[0111] For Conditional Equation 1, in

[00025]$\mathrm{.Math.}=1+\frac{\Delta t_{\mathrm{ao}}}{{\gamma }_1{T}_a-\Delta t_{\mathrm{ao}}},{\gamma }_1>0,$

[0112] it may be calculated as c range that is calculated for values of γ.sub.1 range such as

[00026]$\frac{\Delta t_{a0}}{T_a}<{\gamma }_1<1.25.$
ε.sub.min_γ.sub.2<ε<ε.sub.max_γ.sub.2  [Conditional Equation 2]

[0113] For Conditional Equation 2, in

[00027]$\mathrm{.Math.}=\frac{v_{\mathrm{eff}}\left(r_{\mathrm{mid}}\right).Math.T_{\mathrm{obs}}}{v_{\mathrm{eff}}\left(r_{\mathrm{mid}}\right).Math.T_{\mathrm{obs}}-{\gamma }_2.Math.{\theta }_{\mathrm{az}}.Math.r_{\mathrm{mid}}},{\gamma }_2>0,$

[0114] it may be calculated as ε range that is calculated for values of γ.sub.2 range such as

[00028]$1<{\gamma }_2<0.75.Math.\frac{B_{a\_\mathrm{Targe}t}}{B_{\mathrm{FOV}}}.$

Here,

[0115] [00029]$B_{a\_\mathrm{Target}}=\frac{2v_{\mathrm{eff}}^2\left(r_{\mathrm{mid}}\right)}{\lambda r_{\mathrm{mid}}}.Math.T_{\mathrm{obs}}$

and T.sub.obs is the target observation duration (Target dwell time).

[0116] The embodiments described above may be implemented as a hardware component, a software component, and/or a combination of a hardware component and a software component. For example, the devices, methods, and components described in the embodiments may be implemented by using one or more general computer or specific-purpose computer such as a processor, a controller, an arithmetic logic unit (ALU), a digital signal processor, a microcomputer, a field programmable gate array (FPGA), a programmable logic unit (PLU), a microprocessor, or any other device capable of executing instructions and responding thereto. The processing device may execute an operating system (OS) and one or more software applications executed on the operating system. Further, the processing device may access, store, operate, process, and generate data in response to the execution of software. For convenience of understanding, although it may be described that one processing device is used, one of ordinary skill in the art may understand that the processing device may include a plurality of processing elements and/or a plurality of types of processing elements. For example, the processing device may include a plurality of processors or one processor and one controller. In addition, other processing configurations such as a parallel processor are possible.

[0117] The software may include a computer program, code, instructions, or a combination of one or more of the above, and may configure the processing unit, or instruct the processing unit independently or collectively to operate as desired. Software and/or data may be interpreted by the processing device or, in order to provide instructions or data to the processing device, may be embodied in any type of machine, component, physical device, virtual equipment, computer storage medium or device, or signal wave transmission, permanently or temporarily. The software may be distributed over networked computer systems and stored or executed in a distributed manner. The software and data may be stored on one or more computer-readable recording media.

[0118] The method according to the embodiment may be implemented in the form of program instructions that can be executed through various computer means and recorded in a computer-readable medium. The computer readable medium may include program instructions, data files, data structures, and the like alone or in combination. The program instructions recorded on the medium may be those specially designed and configured for the purposes of the embodiments, or may be known and available to those skilled in computer software. Examples of computer readable recording medium include magnetic media such as hard disks, floppy disks, and magnetic tape, optical media such as CD-ROMs and DVDs, magneto-optical media such as floptical disks, and hardware devices specifically configured to store and execute program instructions such as ROM, RAM, flash memory, and the like. Examples of the program instructions include machine language codes such as those generated by a compiler, as well as high-level language codes that may be executed by a computer using an interpreter, and so on. The hardware device described above may be configured to operate as one or more software modules in order to perform the operations according to the embodiments, and vice versa.

[0119] As described above, although the embodiments have been described with reference to the limited drawings, a person of ordinary skill in the art can apply various technical modifications and variations based on the above. For example, even when the described techniques are performed in an order different from the described method, and/or even when the components of the system, structure, device, circuit, and the like are coupled or combined in a form different from the way described, or replaced or substituted by other components or equivalents, an appropriate result can be achieved.