## 1.

## Introduction

Synthetic aperture radar (SAR) plays a significant role in remote sensing.^{1}2.^{–}^{3} Recently, airborne squint spotlight SAR has become one of the most important topics in microwave remote sensing, as it has the advantages of illuminating a particular area within a single pass of the platform by foreword/backward-squinted antenna and providing higher azimuth resolution by operating at a spotlight mode.^{4}^{,}^{5} Therefore, squint spotlight SAR currently has significant applications in intelligence surveillance, security monitoring, target detection, etc.^{6}^{,}^{7}

However, the state-of-the-art ultrahigh-resolution airborne squint spotlight SAR has to deal with an extremely huge volume of raw data due to the large range cell migration (RCM) that results from the highly squinted angle and high-signal sampling rate corresponding to a high-range resolution. As a result, it leads to higher computation loads and lower efficiency in SAR image formation processing. In this article, the adaptive sliding receive-window (ASRW) technique was introduced for effectively reducing RCM in a -squinted airborne spotlight SAR system. In this system, its receive-window starting time varies pulse-by-pulse as a function of range-walk, namely, the linear term of RCM. Consequently, both data volume and computation loads can be significantly reduced. As a consequence the squint SAR system makes real-time processing feasible, especially for the intelligence reconnaissance airborne SAR system. However, those classical squint SAR data image algorithms, such as the nonlinear chirp-scaling algorithm, extended nonlinear chirp-scaling algorithm, and wave number domain algorithm (WDA),^{8}9.10.^{–}^{11} cannot be employed to focus squint spotlight ASRW-SAR data directly for not considering the sliding receive-window effects.

In this article, a fast image-formation algorithm for ultrahigh-resolution airborne squint spotlight ASRW-SAR was proposed based on the principle of WDA and the azimuth deramping technique. Azimuth deramping preprocessing was implemented for eliminating azimuth spectrum aliasing.^{12} Moreover, WDA is a good method for focusing squint spotlight SAR data because of its high accuracy and efficiency.^{1}^{,}^{13} In the WDA part, bulk compression and modified Stolt mapping were utilized for implementing high-precision focusing. Geometric correction was also employed for compensating the geometric distortion that results from the ASRW technique.

In this article, Sec. 2 analyzes the airborne squint spotlight SAR in detail using the ASRW technique. Subsequently, a fast image-formation algorithm for airborne squint spotlight ASRW-SAR was proposed in Sec. 3. Then, the point targets simulation experiments were carried out in Sec. 4 to prove the validity of the proposed algorithm, and it was verified to have a high-processing accuracy and efficiency for squint spotlight ASRW-SAR data. Finally, conclusions are summarized in Sec. 5.

## 2.

## Airborne Squint Spotlight Synthetic Aperture Radar Using Adaptive Sliding Receive-Window Technique

Figure 1 shows the acquisition geometry of squint spotlight SAR. We define the rectangular coordinate system $OXYZ$: the origin is the nadir at azimuth beginning time, the $x$-axis is parallel to the flight path and points along the azimuth direction, the $y$-axis is perpendicular to the $x$-axis on the ground and points along the range direction, and the $z$-axis provides a right-hand Cartesian coordinate system. As shown in Fig. 1, radar S is at ($x$, 0, $H$), moving at speed $v$. $\phi $ is the squint angle, $r$ is the reference range, and $x$ is the function of the azimuth time $t$ which indicates the position of the airplane along the azimuth direction. C is the scene center and A is a certain target in the scene. Also, the azimuth distance between the scene center C and target A is ${x}_{A}$, where ${x}_{A}=v{t}_{A}$. $R(t)$ denotes the instantaneous slant-range distance from the antenna phase center to a certain target A.

Combining all the geometry, the instantaneous slant-range distance $R(t)$ can be written as^{1}^{,}^{9}

## (1)

$$R(t)=\sqrt{{r}^{2}+{v}^{2}{(t-{t}_{A})}^{2}-2rv(t-{t}_{A})\mathrm{sin}\text{\hspace{0.17em}}\phi}=r-v(t-{t}_{A})\mathrm{sin}\text{\hspace{0.17em}}\phi +\frac{{v}^{2}{(t-{t}_{A})}^{2}{\mathrm{cos}}^{2}\text{\hspace{0.17em}}\phi}{2r}+\cdots \phantom{\rule{0ex}{0ex}}=r+\genfrac{}{}{0ex}{}{\underbrace{\frac{\lambda}{2}{f}_{d}(t-{t}_{A})}}{\text{range}-\text{walk}}+\genfrac{}{}{0ex}{}{\underbrace{\frac{\lambda}{4}{f}_{r}{(t-{t}_{A})}^{2}}}{\text{range}-\text{curve}}+\cdots ,$$^{9}Figure 2 shows range-walk and range-curve variations as a function of the squint angle using the SAR system parameters listed in Table 1 (in Sec. 4). As range-walk is the prominent part of RCM in squint SAR, the removal of range-walk can be considered as the law of the ASRW technique. Also, the receive-window starting time has been changed as follows:

## (2)

$$\mathrm{\Delta}T(t)=2\xb7\frac{\frac{\lambda}{2}{f}_{d}(t-{t}_{A})}{c}=\frac{2v(t-{t}_{A})\mathrm{sin}\text{\hspace{0.17em}}\phi}{c}=\frac{2(x-{x}_{A})\mathrm{sin}\text{\hspace{0.17em}}\phi}{c},$$## Table 1

Synthetic aperture radar (SAR) parameters.

Parameter | Value | Parameter | Value |
---|---|---|---|

λ | 1.0 cm | H | 8000 m |

v | 150 m/s | Pulse repetition frequency (PRF) | 600 Hz |

fs | 1.7 GHz | r | 30.9 Km |

φ | 40 deg/70 deg | D | 0.6 m |

τ | 2 μs | fc | 30.0 GHz |

Figure 3 shows the operation mechanism for airborne squint spotlight SAR using the ASRW technique. The geometry of the squint spotlight SAR operating mode is demonstrated in Fig. 3(a). Also, the echo from the target can be received by sliding the receive window as the red lines show in Fig. 3(b). However, if the receive window is fixed, as the green lines in Fig. 3(c) show, the RCM closes to the receive window and the receive window must be long enough to receive all the RCM. As the RCM increases with high resolution and squint angle, the RCM may be outside the receive window. Figure 3(d) shows both of the two modes’ SAR echoes in two-dimensional (2-D) signal processor memory, and the red lines illustrate that the RCM is reduced. As the range-walk is removed, the RCM is reduced significantly. This implies that the ASRW technique reduces the received echo data volume and saves the storage. Besides, most of the reduced echo data are invalid for the RCM, which would not affect the final imaging results. From the view of information entropy, the ASRW technique actually performs a lossless compression in the data acquisition stage. Therefore, the ASRW technique has no adverse effect on the final imaging.

However, the receive window should start at the specified range of the sampling time due to the radar hardware design principle. This implies that the receive-window starting time is discrete

## (3)

$$\mathrm{\Delta}{T}_{\text{dig}}(t)=\lfloor \frac{\mathrm{\Delta}T(t)}{1/{f}_{s}}\rfloor \xb7{f}_{s},$$So, the phase error, caused by $\mathrm{\Delta}{T}_{\text{shift}}(t)$, can be compensated for by Eq. (5) before the imaging focusing process. Then, the continuous variable $\mathrm{\Delta}T(t)$ will still be used in the subsequent mathematical derivation

## (5)

$${\mathrm{\Omega}}_{1}({f}_{\tau},t)=\mathrm{exp}\{j2\pi {f}_{\tau}\xb7\mathrm{\Delta}{T}_{\text{shift}}(t)\},$$## 3.

## Adaptive Sliding Receive-Window Synthetic Aperture Radar Fast Image-Formation Algorithm

Three components were proposed in this section: First, the expression of 2-D matched filtering function in the wave number domain was derived in detail. Then, based on the principle of WDA, a new Stolt-mapping relationship was illustrated according to the 2-D matched filtering function. Finally, a presentation of the process flow of the fast image-formation algorithm was given, which was proven to focus the ultrahigh-resolution airborne squint spotlight ASRW-SAR data accurately and efficiently.

## 3.1.

### Two-Dimensional Matched Filtering Function

In the subsequent analysis, the antenna beam pattern, backscatter coefficient, and other nonessential amplitude factors have been corrected. So, the point-target spectrum of the innovative mode can be given by

## (6)

$$S(\tau ,x)=p[\tau -\frac{2R(x)}{c}-\frac{2(x-{x}_{A})\mathrm{sin}\text{\hspace{0.17em}}\phi}{c}]\xb7\mathrm{exp}\{-j\pi b{[\tau -\frac{2R(x)}{c}-\frac{2(x-{x}_{A})\mathrm{sin}\text{\hspace{0.17em}}\phi}{c}]}^{2}\}\phantom{\rule{0ex}{0ex}}\xb7\mathrm{exp}\{-j\frac{4\pi R(x)}{\lambda}\},$$## (7)

$${h}_{r}(\tau ,x)=\delta (\tau -\frac{2\mathrm{\Delta}R(x)}{c}-\frac{2(x-{x}_{A})\mathrm{sin}\text{\hspace{0.17em}}\phi}{c})\xb7\mathrm{exp}(-j\frac{4\pi {f}_{c}\mathrm{\Delta}R(x)}{c}),$$## (8)

$${h}_{r}({f}_{\tau},x)=\mathrm{exp}\{-j4\pi \frac{{f}_{c}+{f}_{\tau}}{c}(\mathrm{\Delta}R(x)+x\mathrm{sin}\text{\hspace{0.17em}}\phi )\}\xb7\mathrm{exp}\left\{j4\pi \frac{{f}_{c}}{c}x\mathrm{sin}\text{\hspace{0.17em}}\phi \right\}\phantom{\rule{0ex}{0ex}}\xb7\mathrm{exp}\{j{\mathrm{\Phi}}_{1}({f}_{\tau},{x}_{A})\},$$## (9)

$${h}_{r}({f}_{\tau},x)=\mathrm{exp}\{-j4\pi \frac{{f}_{c}+{f}_{\tau}}{c}(\mathrm{\Delta}R(x)+x\mathrm{sin}\text{\hspace{0.17em}}\phi )\}\xb7\mathrm{exp}\left\{j4\pi \frac{{f}_{c}}{c}\mathrm{sin}\text{\hspace{0.17em}}\phi \right\},$$## (10)

$${H}_{r}({f}_{\tau},{k}_{x})=\int \mathrm{exp}\{-j{k}_{\mathrm{rc}}(\mathrm{\Delta}R(x)+x\mathrm{sin}\text{\hspace{0.17em}}\phi )\}\xb7\mathrm{exp}\{-j2\pi {f}_{d}\frac{x}{v}\}\mathrm{exp}\{-j{k}_{x}x\}\mathrm{d}x.$$Substituting ${k}_{0}=2\pi {f}_{d}/v$ into Eq. (10),

## (11)

$${H}_{r}({f}_{\tau},{k}_{x}+{k}_{0})=\int \mathrm{exp}\{-j{k}_{\mathrm{rc}}(\mathrm{\Delta}R(x)+x\mathrm{sin}\text{\hspace{0.17em}}\phi )\}\xb7\mathrm{exp}\{-j({k}_{x}+{k}_{0})x\}\mathrm{d}x.$$In squint mode, the center of the azimuth wave number is –${k}_{0}x$, and ${k}_{x}\in $ [$-\pi ({f}_{\text{PRF}}-{f}_{d})/v$ and $\pi ({f}_{\text{PRF}}-{f}_{d})/v$ ] in Eq. (11), where ${f}_{\text{PRF}}$ is equal to PRF. By using variable substitution, Eq. (11) can be written as follows:

## (12)

$${H}_{r}({f}_{\tau},{k}_{x})=\int \mathrm{exp}\{-j{k}_{rc}(\mathrm{\Delta}R(x)+x\mathrm{sin}\text{\hspace{0.17em}}\phi )\}\xb7\mathrm{exp}\{-j{k}_{x}x\}\mathrm{d}x,$$## (13)

$${x}_{0}=r\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi -\frac{({k}_{\mathrm{rc}}\mathrm{sin}\text{\hspace{0.17em}}\phi +{k}_{x})r\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\phi}{\sqrt{{k}_{\mathrm{rc}}^{2}-{({k}_{\mathrm{rc}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi +{k}_{x})}^{2}}}-{x}_{A}.$$So, the expression of the 2-D matched filtering function in the range frequency domain and azimuth wave number domain is

## (14)

$${H}_{r}({f}_{\tau},{k}_{x})=\mathrm{exp}\{-jr[\mathrm{cos}\text{\hspace{0.17em}}\phi \sqrt{{k}_{\mathrm{rc}}^{2}-{({k}_{\mathrm{rc}}\text{\hspace{0.17em}}\mathrm{sin}\phi +{k}_{x})}^{2}}+\mathrm{sin}\text{\hspace{0.17em}}\phi ({k}_{\mathrm{rc}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi +{k}_{x})-{k}_{\mathrm{rc}}]\}.$$Then, taking FFT of Eq. (14) with respect to $r$, we can get the expression of the 2-D matched filtering function in the wave number domain as

## (15)

$${H}_{r}({k}_{r},{k}_{x})=\delta ({k}_{r}+[\mathrm{cos}\text{\hspace{0.17em}}\phi \sqrt{{k}_{\mathrm{rc}}^{2}-{({k}_{\mathrm{rc}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi +{k}_{x})}^{2}}+\mathrm{sin}\text{\hspace{0.17em}}\phi ({k}_{\mathrm{rc}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi +{k}_{x})-{k}_{\mathrm{rc}}]).$$## 3.2.

### New Stolt-Mapping Relationship

Several modified Stolt-mapping relationships have been proposed in Refs. 14 and 15, which rewrite the expression of RCM. But they do not essentially change the mapping relationship between the frequency domain and wave number domain. For the image-formation algorithm of the squint spotlight ASRW-SAR, a new mapping relationship should be found immediately due to the sliding receive window. Assuming the range compress signal is $p({f}_{\tau},{k}_{x})$, an IFFT was implemented with respect to ${f}_{\tau}$ after the matched filtering operation

## (16)

$$P(r,{k}_{x})=\int p({f}_{\tau},{k}_{x}){H}_{r}^{*}({f}_{\tau},{k}_{x})\mathrm{exp}\left(j\frac{4\pi {f}_{\tau}}{c}r\right)\mathrm{d}{f}_{\tau},$$## (17)

$$P({k}_{r},{k}_{x})=\int \int p({f}_{\tau},{k}_{x})\xb7{H}_{r}^{*}({f}_{\tau},{k}_{x})\xb7\mathrm{exp}\left(j\frac{4\pi {f}_{\tau}}{c}r\right)\xb7\mathrm{exp}(-j{k}_{r}r)\mathrm{d}r\mathrm{d}{f}_{\tau}\phantom{\rule{0ex}{0ex}}=\int p({f}_{\tau},{k}_{x}){H}_{r}^{*}(\frac{4\pi {f}_{\tau}}{c}-{k}_{r},{k}_{x})\mathrm{d}{f}_{\tau}.$$Substituting Eq. (15) into Eq. (17),

## (18)

$$P({k}_{r},{k}_{x})=\int \delta ({k}_{r}-[\mathrm{cos}\text{\hspace{0.17em}}\phi \sqrt{{k}_{\mathrm{rc}}^{2}-{({k}_{\mathrm{rc}}\mathrm{sin}\text{\hspace{0.17em}}\phi +{k}_{x})}^{2}}+\mathrm{sin}\text{\hspace{0.17em}}\phi ({k}_{\mathrm{rc}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi +{k}_{x})-\frac{4\pi {f}_{c}}{c}\left]\right)\phantom{\rule{0ex}{0ex}}\xb7p({f}_{\tau},{k}_{x})\mathrm{d}{f}_{\tau}.$$So, a new mapping relationship is illustrated as follows

## (19)

$${k}_{r}=\mathrm{cos}\text{\hspace{0.17em}}\phi \sqrt{{k}_{\mathrm{rc}}^{2}-{({k}_{\mathrm{rc}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi +{k}_{x})}^{2}}+\mathrm{sin}\text{\hspace{0.17em}}\phi ({k}_{\mathrm{rc}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi +{k}_{x})-\frac{4\pi {f}_{c}}{c}.$$According to the new mapping relationship, it is easy to convert $p({f}_{\tau},{k}_{x})$ from the range frequency domain to the 2-D wave number domain $p({k}_{r},{k}_{x})$ by interpolation. Finally, fine-focusing the image can be done using the azimuth IFFT, geometric correction, and range IFFT operation.

## 3.3.

### Processing Scheme of Fast Image-Formation Algorithm

As Fig. 4 shows, the processing steps of the fast image-formation algorithm are divided into six stages as follows: time-shift compensation, azimuth deramping, data transform, bulk compression, modified Stolt mapping, and geometric correction.

As analyzed in Sec. 2, the time-shift $\mathrm{\Delta}{T}_{\text{shift}}$ ($t$) causes a phase error. So, the first step of the proposed algorithm is to compensate for the phase error as in Eq. (5).

The overlapped azimuth spectrum limits the direct application of the available stripmap imaging algorithm in the spotlight SAR mode. To resolve the problem, azimuth deramping is adopted to eliminate the azimuth spectrum aliasing. Azimuth deramping shares the same processing flow as in Ref. 16.

The data transform stage includes four processing steps, namely, range FFT, phase compression, azimuth FFT, and residual phase compensation. For ASRW-SAR, the receive-window starting time is changed pulse-by-pulse in the azimuth direction. So, the Doppler history is variable in azimuth, which means that the target located only in the scene center is precisely focused. As analyzed in Sec. 3.1, in order to guarantee the focusing accuracy of all the targets located in different azimuth positions, the third exponential phase term in Eq. (8) should be compensated for, which is given by

## (20)

$${\mathrm{\Omega}}_{2}({f}_{\tau},x)=\mathrm{exp}\{-j\frac{4\pi {f}_{\tau}{x}_{A}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi}{c}\}.$$Only the targets apart from the scene center at the distance ${x}_{A}$ in azimuth direction will be precisely focused. To ensure the imaging quality over the whole scene, a method based on block processing was proposed in Fig. 5. Echo data can be divided into $n$ blocks along the azimuth direction, and each block can be compensated for by using different ${x}_{A}$ values.

Moreover, after performing the azimuth FFT, the residual phase caused by azimuth deramping processing^{16} should be compensated for by

In the bulk compensation stage, the reference function is different from the traditional function in Ref. 1. Therefore, in a practical procedure, the reference function in bulk compensation is given by

## (22)

$${\mathrm{\Omega}}_{4}({f}_{\tau},{k}_{x})=\mathrm{exp}\{jr\xb7\mathrm{cos}\text{\hspace{0.17em}}\phi \sqrt{{k}_{\mathrm{rc}}^{2}-{({k}_{\mathrm{rc}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi +{k}_{x})}^{2}}+\mathrm{sin}\text{\hspace{0.17em}}\phi ({k}_{\mathrm{rc}}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\phi +{k}_{x})\}\phantom{\rule{0ex}{0ex}}\xb7\mathrm{exp}\{-j\pi \frac{{f}_{\tau}^{2}}{b}\}\xb7\mathrm{exp}\{-j\frac{4\pi {f}_{\tau}{r}_{\mathrm{min}}}{c}\},$$Different from the conventional Stolt mapping, the proposed Stolt mapping is modified as in Eq. (19) due to the change in the data-mapping relationship caused by the ASRW technique. The modified Stolt mapping performs the differential RCMC, differential SRC, and differential azimuth compression. After the modified Stolt mapping processing, the whole scene is acurately focused.

The last processing part is the geometric correction, which compensates for the change in the receive-window starting time caused by the ASRW technique. After the modified Stolt mapping, the azimuth IFFT is performed to transform the data into the azimuth time and range frequency domain. After that, geometric distortion is corrected by

## (23)

$${\mathrm{\Omega}}_{5}({f}_{\tau},x)=\mathrm{exp}\left\{j4\pi {f}_{\tau}\frac{\mathrm{sin}\text{\hspace{0.17em}}\phi}{c}x\right\}.$$Finally, a range IFFT is applied to get the fine-focusing and a nondistorted SAR image.

## 4.

## Experimental Results

In this section, a Ka-band linear frequency modulation transmission signal is adopted for simulation, and the azimuth resolution ${\rho}_{a}$ is 0.1 m. First, to prove the effectiveness of the ASRW technique, the experiments comparing the raw echo between the ASRW-SAR and conventional-SAR are carried out. Then, to verify the validity of the proposed image-formation algorithm, the experiments with point targets imaging experiments are carried out. The specific simulation parameters are listed in Table 1, and the simulation scene ($500\times 500\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{m}}^{2}$, resulting in $5000\times 5000$ pixels in the SAR image) is illustrated by Fig. 6.

In ASRW-SAR, the ASRW technique changes the starting time of the receive window in each azimuth direction by removing the range-walk. Also, the length of the range dimension is dependent on the RCM, whose main part is range-walk. Figure 7 shows the raw echo data of ASRW-SAR and conventional-SAR systems at different squint angles ($\phi =\phantom{\rule{0ex}{0ex}}40\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}/70\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$). Also, the range sample number in a conventional-SAR and ASRW-SAR are ${N}_{r}$ and ${N}_{r}^{\prime}$, respectively. As shown in Fig. 7, the RCM is reduced by ASRW, and the range sampling number decreases significantly. Further, the higher the squinted angle, the bigger the range sampling number scaling and the more immediate the demand of ASRW.

In the imaging experiments part, the imaging results of a conventional WDA and the proposed imaging algorithm are shown in Fig. 8. The comparison results of the proposed imaging algorithm with a conventional WDA shows that the proposed method is effective. Also, Figs. 9 and 10 show the interpolated contour plots of the corresponding focusing targets A, C, and E in different squint angles. The spatial resolution (azimuth resolution ${\rho}_{a}$, range resolution ${\rho}_{r}$), peak side lobe ratio, and integrated side lobe ratio of the three point targets processing results are listed in Table 2. Both the contour plots and image evaluation results demonstrate that the fine-focusing performance of the proposed improved image-formation algorithm and all the evaluation results indicate that the image results are slightly affected by the spatial variant slant range. Besides, the comparison of the computation load of conventional-SAR and ASRW-SAR is represented in the Appendix, and the computation efficiency is $\gamma =4.1/17.0$ in $\phi =40\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}/70\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$. The advantage of less-computation load indicates its potential suitability in a real-time imaging system.

## Table 2

Targets A, B, and C evaluation results in different squint angles.

Azimuth | Range | ||||||
---|---|---|---|---|---|---|---|

ρa (m) | Peak side lobe ratio (PSLR) (dB) | Integrated side lobe ratio (ISLR) (dB) | ρr* (m) | PSLR (dB) | ISLR (dB) | ||

φ=40 deg | A | 0.100 | −13.28 | −10.28 | 0.089 | −13.30 | −10.12 |

B | 0.100 | −13.26 | −10.27 | 0.089 | −13.30 | −10.12 | |

C | 0.100 | −13.26 | −10.27 | 0.089 | −13.30 | −10.12 | |

D | 0.100 | −13.26 | −10.27 | 0.089 | −13.30 | −10.12 | |

E | 0.100 | −13.27 | −10.28 | 0.089 | −13.30 | −10.12 | |

φ=70 deg | A | 0.102 | −12.88 | −9.62 | 0.088 | −13.32 | −10.39 |

B | 0.101 | −12.86 | −9.63 | 0.087 | −13.31 | −10.39 | |

C | 0.101 | −12.85 | −9.64 | 0.087 | −13.31 | −10.39 | |

D | 0.101 | −12.87 | −9.62 | 0.087 | −13.31 | −10.39 | |

E | 0.103 | −12.89 | −9.61 | 0.088 | −13.32 | −10.38 |

${\rho}_{r}^{*}$ denotes the slant range resolution.

## 5.

## Conclusion

In this article, an ASRW technique was introduced for efficiently acquiring the airborne highly squint spotlight SAR echo. The airborne squint spotlight ASRW-SAR varies its receive-window starting time pulse-by-pulse as a function of range-walk, namely, the linear term of RCM. So, the RCM of ASRW-SAR is reduced dramatically and the range sampling number is decreased, which means the data volume and the computation loads are reduced. Therefore, the ASRW technique decreases the SAR system complexity, and has great potential value in real-time imaging systems.

In order to deal with the varying receive-window starting times in airborne squint spotlight ASRW-SAR, a fast image-formation algorithm is presented for highly precise focusing of the ultrahigh-resolution (e.g., 0.1 m) airborne squint spotlight ASRW-SAR data. First, azimuth deramping preprocessing is performed to eliminate the azimuth spectrum aliasing. Then, bulk compression and modified Stolt mapping are utilized for precise focusing. At last, geometric correction is implemented to remove the effect of the varying receive-window starting time by the ASRW technique. The validity and accuracy of the proposed algorithm have been demonstrated by point targets simulation data. Also, the processing results of point targets at different squint angles indicate that the proposed algorithm is suitable for high-resolution and high squint angle airborne spotlight SAR image formation. Further, the computation load is reduced, which would be more markedly decreased in high squint angle, e.g., 17 times improvement of computation efficiency at a 70-deg squint angle.

## Appendices

## Appendix:

### Comparison of Computation Load of Conventional-SAR and ASRW-SAR

The number of multiplications is adopted as a computation load indicator. Let the interpolation kernel length be ${M}_{\text{ken}}$. Also, the azimuth sample numbers in conventional SAR and ASRW-SAR are ${N}_{a}$ and ${N}_{a}^{\prime}$, respectively. The particular computational load in each imaging stage is listed in Table 3.

## Table 3

Computation load comparison.

Conventional SAR | ASRW-SAR | |
---|---|---|

Azimuth deramping | 8NaNr+2NaNr log2Na | 8Na′Nr′+2Na′Nr′ log2Na′ |

Range FFT | 2NaNr log2Nr | 2Na′Nr′ log2Nr′ |

Azimuth FFT | 2NaNr log2Na | 2Na′Nr′ log2Na′ |

Bulk compression | 4NaNr | 4Na′Nr′ |

Stolt mapping | 2MkenNaNr | 2MkenNa′Nr′ |

Azimuth IFFT | 2NaNr log2Na | 2Na′Nr′ log2Na′ |

Geometric correction | 0 | 4Na′Nr′ |

Range IFFT | 2NaNr log2Nr | 2Na′Nr′ log2Nr′ |

${C}_{\text{con}}$ denotes the entire computational load for conventional-SAR.

## (24)

$${C}_{\text{con}}=(12+2{M}_{\text{ken}}){N}_{a}{N}_{r}+6{N}_{a}{N}_{r}{\mathrm{log}}_{2}\text{\hspace{0.17em}}{N}_{a}+4{N}_{a}{N}_{r}{\mathrm{log}}_{2}\text{\hspace{0.17em}}{N}_{r}.$$Also, the entire computational load ${C}_{\mathrm{ASRW}}$ for ASRW-SAR is expressed as

## (25)

$${C}_{\mathrm{ASRW}}=(16+2{M}_{\text{ken}}){N}_{a}^{\prime}{N}_{r}^{\prime}+6{N}_{a}^{\prime}{N}_{r}^{\prime}{\mathrm{log}}_{2}\text{\hspace{0.17em}}{N}_{a}^{\prime}+4{N}_{a}^{\prime}{N}_{r}^{\prime}{\mathrm{log}}_{2}\text{\hspace{0.17em}}{N}_{r}^{\prime}.$$Assume that the computational efficiency $\gamma $ is defined as

With the condition of the same parameters as listed in Table 1, ${N}_{a}$ and ${N}_{a}^{\prime}$ are equivalent, ${N}_{a}={N}_{a}^{\prime}=[(v\xb7\mathrm{PRF})/({\rho}_{a}\xb7{f}_{r})]$ assuming an 8-point interpolation is adopted (${M}_{\text{ken}}=8$). And in different squint angle $\phi =40\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}/70\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$, ${N}_{a}$ (or ${N}_{a}^{\prime}$) is equal to 17,500/87,792. Also, the ${N}_{r}$ and ${N}_{r}^{\text{'}}$ are listed in Fig. 6 in different squint angles. Therefore, the computation efficiency $\gamma $ is calculated by Eq. (26) in different squint angles, and the final value $\gamma $ is listed in Table 2.

## Acknowledgments

The authors would like to thank the anonymous reviewers for their valuable comments and useful suggestions. This work was supported in part by National Natural Science Foundation of China (NSFC) under Grant No. 61132006, and, in part, by National Natural Science Foundation of China (NSFC) under Grant No. 61171123.

## References

I. G. CummingF. H. Wong, Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation, Artech House, Norwood, Massachusetts (2005).Google Scholar

S.-I. HwangH. WangK. Ouchi, “Comparision and evaluation of ship detection and identification alogrithms using small boats and ALOS-PALSAR,” IEICE Trans. Commun. E92-B(12), 3883–3892 (2009).ITRCEC0916-8516http://dx.doi.org/10.1587/transcom.E92.B.3883Google Scholar

H. Anget al., “Assessment of building damage in 2008 wenchuan earthquake from multi-temporal SAR images using getis statistic,” IEICE Trans. Commun. E94-B(11), 2983–2986 (2011).ITRCEC0916-8516http://dx.doi.org/10.1587/transcom.E94.B.2983Google Scholar

S. Cimminoet al., “Efficient spotlight SAR raw signal simulation of extended scenes,” IEEE Trans. Geosci. Remote Sens. 41, 2329–2337 (2003).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2003.815239Google Scholar

G.W. DavidsonI. Cumming, “Signal properties of spaceborne squint-mode SAR,” IEEE Trans. Geosci. Remote Sens. 35(3), 611–617 (1997).IGRSD20196-2892http://dx.doi.org/10.1109/36.581976Google Scholar

V. C. Kooet al., “A new unmanned aerial vehicle synthetic aperture radar for environmental monitoring,” Prog. Electromagn. Res. 125, 527–543 (2012).PELREX1043-626Xhttp://dx.doi.org/10.2528/PIER11121101Google Scholar

D. Realeet al., “Advanced techniques and new high resolution SAR sensors for monitoring urban areas,” in Geoscience and Remote Sensing Symposium (IGARSS), 2010 IEEE International, Honolulu, HI, pp. 1800–1803 (2011).Google Scholar

H. ChengT. LongY. Tian, “An improved nonlinear chirp scaling algorithm based on curved trajectory in geosynchronous SAR,” Prog. Electromagn. Res. 135, 481–513 (2013).PELREX1043-626Xhttp://dx.doi.org/10.2528/PIER12092603Google Scholar

D. X. Anet al., “Extended nonlinear chirp scaling algorithm for high-resolution highly squint SAR data focusing,” IEEE Trans. Geosci. Remote Sens. 50(9), 3595–3609 (2012).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2012.2183606Google Scholar

Y. WangO. LoffeldS. Knedlik, “Spotlight-mode SAR data focusing using a modified wavenumber domain algorithm,” in Geoscience and Remote Sensing Symposium, IGARSS, pp. 567–570 (2007).Google Scholar

H.-S. ShinJ.-T. Lim, “Omega-k algorithm for spaceborne spotlight SAR imaging,” IEEE Geosci. Remote Sens. Lett. 9(3), 343–347 (2012).IGRSBY1545-598Xhttp://dx.doi.org/10.1109/LGRS.2011.2168380Google Scholar

L. Riccardoet al., “Spotlight SAR data focusing based on a two-step processing approach,” IEEE Trans. Geosci. Remote Sens. 39(9), 1993–2004 (2001).IGRSD20196-2892http://dx.doi.org/10.1109/36.951090Google Scholar

R. Bamler, “A comparision of range-Doppler and wavenumber domain SAR focusing algorithm,” IEEE Geosci. Remote Sens. 30(4), 706–713 (1992).IGRSD20196-2892http://dx.doi.org/10.1109/36.158864Google Scholar

M. VandewalR. SpeckH. Sub, “Efficient and precise processing for squinted spotlight SAR through a modified stolt mapping,” EURSAIP J. Adv. Signal Process. 2007(1), 33–39 (2007).1110-8657http://dx.doi.org/10.1155/2007/59704Google Scholar

A. ReigberE. AlivizatosA. Potsis, “Extended wavenumber domain synthetic apeture radar focusing with integrated motion compensation,” IEEE Proc., Sonar Navig. 153(3), 301–310 (2006).IRSNE21350-2395http://dx.doi.org/10.1049/ip-rsn:20045087Google Scholar

D. GuoH. XuJ. Li, “Extended wavenumber domain algorithm for highly squinted sliding spotlight SAR data processing,” Prog. Electromagn. Res. 114, 17–32 (2011).PELREX1043-626XGoogle Scholar

## Biography

**Wei Yang** received the MS and PhD degrees from Beihang University, China, in 2008 and 2011, respectively. He has been a lecturer with the School of Electronics and Information Engineering, Beihang University, since 2013. His current research interests include high-resolution spaceborne SAR image formation and modeling and simulation for spaceborne SAR systems.

**Hong-cheng Zeng** received the BS degrees in college of Information and Electrical Engineering from China Agriculture University in 2011. He is currently working towards his PhD at Beihang University. His research interest is high-resolution spaceborne SAR image formation.

**Jie Chen** received the BS and PhD degrees from Beihang University, China, in 1996 and 2002, respectively. He is a professor with the School of Electronics and Information Engineering, Beihang University (BUAA) from 2004. His current research interests include remote sensing information acquisition and signal processing; topside ionosphere exploration based on spaceborne HF-SAR; high-resolution spaceborne SAR image formation; modeling and simulation for spaceborne SAR systems.

**Peng-bo Wang** received the PhD degrees from Beihang University (BUAA), China, in 2007. He has been a lecturer with the School of Electronics and Information Engineering, Beihang University, since 2010. His current research interests include remote sensing information acquisition and signal processing; high-resolution spaceborne SAR image formation; modeling and simulation for spaceborne SAR systems.