Fast image-formation algorithm for ultrahigh-resolution airborne squint spotlight synthetic aperture radar based on adaptive sliding receive-window technique

Wei Yang; Hong-cheng Zeng; Jie Chen; Peng-bo Wang

doi:10.1117/1.JRS.8.083620

4 June 2014 Fast image-formation algorithm for ultrahigh-resolution airborne squint spotlight synthetic aperture radar based on adaptive sliding receive-window technique

Wei Yang, Hong-cheng Zeng, Jie Chen, Peng-bo Wang

Author Affiliations +

Funded by: National Natural Science Foundation of China (NSFC), National Natural Science Foundation of China, Natural Science Foundation of China, Key National Natural Science Foundation of China

Journal of Applied Remote Sensing, Vol. 8, Issue 1, 083620 (June 2014). https://doi.org/10.1117/1.JRS.8.083620

Abstract

Adaptive sliding receive-window (ASRW) technique was usually introduced in airborne squint synthetic aperture radar (SAR) systems. Airborne squint spotlight SAR varies its receive-window starting time pulse-by-pulse as a function of range-walk, namely, the linear term of range cell migration (RCM). As a result, a huge data volume of the highly squint spotlight SAR echo signal can be significantly reduced. Because the ASRW technique changes the echo-receive starting time and Doppler history, the conventional image algorithm cannot be employed to directly focus airborne squint spotlight ASRW-SAR data. Therefore, a fast image-formation algorithm, based on the principle of the wave number domain algorithm (WDA) and azimuth deramping processing, was proposed for accurately and efficiently focusing the squint spotlight ASRW-SAR data. Azimuth deramping preprocessing was implemented for eliminating azimuth spectrum aliasing. Moreover, bulk compression and modified Stolt mapping were utilized for high-precision focusing. Additionally, geometric correction was employed for compensating the image distortion resulting from the ASRW technique. The proposed algorithm was verified by evaluating the image performance of point targets in different squint angles. In addition, a detailed analysis of computation loads in the appendix indicates that the processing efficiency can be greatly improved, e.g., the processing efficiency could be improved by 17 times in the 70-deg squint angle by applying the proposed image algorithm to the squint spotlight ASRW-SAR data.

1. Introduction

Synthetic aperture radar (SAR) plays a significant role in remote sensing.¹^–³ 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.⁴^,⁵ Therefore, squint spotlight SAR currently has significant applications in intelligence surveillance, security monitoring, target detection, etc.⁶^,⁷

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),⁸^–¹¹ 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.¹² Moreover, WDA is a good method for focusing squint spotlight SAR data because of its high accuracy and efficiency.¹^,¹³ 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 $O X Y Z$ : 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$ . $φ$ 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.

Fig. 1

Acquisition geometry of airborne squint spotlight synthetic aperture radar (SAR) mode.

Combining all the geometry, the instantaneous slant-range distance $R (t)$ can be written as¹^,⁹

Eq. (1)

R (t) = \sqrt{r^{2} + v^{2} {(t - t_{A})}^{2} - 2 r v (t - t_{A}) \sin φ} = r - v (t - t_{A}) \sin φ + \frac{v^{2} {(t - t_{A})}^{2} \cos^{2} φ}{2 r} + \dots = r + \binom{\underset{︸}{\frac{λ}{2} f_{d} (t - t_{A})}}{range - walk} + \binom{\underset{︸}{\frac{λ}{4} f_{r} {(t - t_{A})}^{2}}}{range - curve} + \dots,

where

λ

is the signal wavelength and

f_{d} = - 2 v \sin φ / λ

and

f_{r} = 2 v^{2} \cos^{2} φ / (λ r)

are referred to as the Doppler centroid and Doppler rate, respectively.⁹ 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:

Eq. (2)

Δ T (t) = 2 \cdot \frac{\frac{λ}{2} f_{d} (t - t_{A})}{c} = \frac{2 v (t - t_{A}) \sin φ}{c} = \frac{2 (x - x_{A}) \sin φ}{c},

where

c

is the velocity of light and

x = v t

.

Fig. 2

Range cell migration (RCM) variation as a function of squint angle.

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
$f_{s}$	1.7 GHz	$r$	30.9 Km
$φ$	$40 \deg / 70 \deg$	$D$	0.6 m
$τ$	2 μs	$f_{c}$	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.

Fig. 3

Operation mechanism for airborne squint spotlight SAR using adaptive sliding receive window (ASRW).

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

Eq. (3)

Δ T_{dig} (t) = ⌊ \frac{Δ T (t)}{1 / f_{s}} ⌋ \cdot f_{s},

where

⌊ \cdot ⌋

is an integral function. Also, a time-shift of the starting time results from

Δ T_{dig} (t)

is given by

Eq. (4)

Δ T_{shift} (t) = Δ T (t) - Δ T_{dig} (t) .

So, the phase error, caused by $Δ T_{shift} (t)$ , can be compensated for by Eq. (5) before the imaging focusing process. Then, the continuous variable $Δ T (t)$ will still be used in the subsequent mathematical derivation

Eq. (5)

Ω_{1} (f_{τ}, t) = \exp {j 2 π f_{τ} \cdot Δ T_{shift} (t)},

where

f_{τ}

is range frequency.

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

Eq. (6)

S (τ, x) = p [τ - \frac{2 R (x)}{c} - \frac{2 (x - x_{A}) \sin φ}{c}] \cdot \exp {- j π b {[τ - \frac{2 R (x)}{c} - \frac{2 (x - x_{A}) \sin φ}{c}]}^{2}} \cdot \exp {- j \frac{4 π R (x)}{λ}},

where

p [\cdot]

is the signal envelope,

b

is the range chirp FM rate, and

τ

is the range time. Simply put, after range compression, the 2-D matched filtering function signal

h_{r}

(

τ, x

) is given by adopting the impulse response function instead of the sinc function as follows:

Eq. (7)

h_{r} (τ, x) = δ (τ - \frac{2 Δ R (x)}{c} - \frac{2 (x - x_{A}) \sin φ}{c}) \cdot \exp (- j \frac{4 π f_{c} Δ R (x)}{c}),

where

Δ R (x) = R (x)

–

r

and

f_{c}

denotes the carrier frequency. Because the position of the range-compressed signal is changed with

Δ T (t)

, the wave number distribution varied accordingly. By applying the principle of stationary phase (POSP) to Eq. (7),

h_{r}

(

τ, x

) is changed from the range time domain to the range frequency domain

Eq. (8)

h_{r} (f_{τ}, x) = \exp {- j 4 π \frac{f_{c} + f_{τ}}{c} (Δ R (x) + x \sin φ)} \cdot \exp {j 4 π \frac{f_{c}}{c} x \sin φ} \cdot \exp {j Φ_{1} (f_{τ}, x_{A})},

where

Φ_{1} (f_{τ}, x_{A}) = 4 π f_{τ} x_{A} \sin φ / c

. As

Φ_{1} (f_{τ}, x_{A})

is independent of

x

, it can be compensated for by a matched filter in the range frequency and azimuth time domain. Therefore,

h_{r}

(f_{τ}, x)

can be rewritten after compensation as follows:

Eq. (9)

h_{r} (f_{τ}, x) = \exp {- j 4 π \frac{f_{c} + f_{τ}}{c} (Δ R (x) + x \sin φ)} \cdot \exp {j 4 π \frac{f_{c}}{c} \sin φ},

where the azimuth wave number

k_{x} = 2 π f_{a} / v

,

k_{0} = 2 π f_{d} / v

, the range wave number

k_{rc} = 4 π (f_{c} + f_{τ}) / c

, and

f_{a}

is the Doppler frequency. Then, the signal

h_{r}

(

f_{τ}, x

) is converted to the range frequency and azimuth wave number domain

Eq. (10)

H_{r} (f_{τ}, k_{x}) = \int \exp {- j k_{rc} (Δ R (x) + x \sin φ)} \cdot \exp {- j 2 π f_{d} \frac{x}{v}} \exp {- j k_{x} x} d x .

Substituting $k_{0} = 2 π f_{d} / v$ into Eq. (10),

Eq. (11)

H_{r} (f_{τ}, k_{x} + k_{0}) = \int \exp {- j k_{rc} (Δ R (x) + x \sin φ)} \cdot \exp {- j (k_{x} + k_{0}) x} d x .

In squint mode, the center of the azimuth wave number is – $k_{0} x$ , and $k_{x} \in$ [ $- π (f_{PRF} - f_{d}) / v$ and $π (f_{PRF} - f_{d}) / v$ ] in Eq. (11), where $f_{PRF}$ is equal to PRF. By using variable substitution, Eq. (11) can be written as follows:

Eq. (12)

H_{r} (f_{τ}, k_{x}) = \int \exp {- j k_{r c} (Δ R (x) + x \sin φ)} \cdot \exp {- j k_{x} x} d x,

where

k_{x} \in

(

- π f_{PRF} / v

,

π f_{PRF} / v

). Then, by using POSP, the stationary phase point

x_{0}

is given by

Eq. (13)

x_{0} = r \sin φ - \frac{(k_{rc} \sin φ + k_{x}) r \cos φ}{\sqrt{k_{rc}^{2} - {(k_{rc} \sin φ + 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

Eq. (14)

H_{r} (f_{τ}, k_{x}) = \exp {- j r [\cos φ \sqrt{k_{rc}^{2} - {(k_{rc} \sin φ + k_{x})}^{2}} + \sin φ (k_{rc} \sin φ + k_{x}) - k_{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

Eq. (15)

H_{r} (k_{r}, k_{x}) = δ (k_{r} + [\cos φ \sqrt{k_{rc}^{2} - {(k_{rc} \sin φ + k_{x})}^{2}} + \sin φ (k_{rc} \sin φ + k_{x}) - k_{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_{τ}, k_{x})$ , an IFFT was implemented with respect to $f_{τ}$ after the matched filtering operation

Eq. (16)

P (r, k_{x}) = \int p (f_{τ}, k_{x}) H_{r}^{*} (f_{τ}, k_{x}) \exp (j \frac{4 π f_{τ}}{c} r) d f_{τ},

where

H_{r}^{*} (f_{τ}, k_{x})

is the complex conjugate of

H_{r} (f_{τ}, k_{x})

. Then, by applying FFT with respect to

r

Eq. (17)

P (k_{r}, k_{x}) = \int \int p (f_{τ}, k_{x}) \cdot H_{r}^{*} (f_{τ}, k_{x}) \cdot \exp (j \frac{4 π f_{τ}}{c} r) \cdot \exp (- j k_{r} r) d r d f_{τ} = \int p (f_{τ}, k_{x}) H_{r}^{*} (\frac{4 π f_{τ}}{c} - k_{r}, k_{x}) d f_{τ} .

Substituting Eq. (15) into Eq. (17),

Eq. (18)

P (k_{r}, k_{x}) = \int δ (k_{r} - [\cos φ \sqrt{k_{rc}^{2} - {(k_{rc} \sin φ + k_{x})}^{2}} + \sin φ (k_{rc} \sin φ + k_{x}) - \frac{4 π f_{c}}{c}]) \cdot p (f_{τ}, k_{x}) d f_{τ} .

So, a new mapping relationship is illustrated as follows

Eq. (19)

k_{r} = \cos φ \sqrt{k_{rc}^{2} - {(k_{rc} \sin φ + k_{x})}^{2}} + \sin φ (k_{rc} \sin φ + k_{x}) - \frac{4 π f_{c}}{c} .

According to the new mapping relationship, it is easy to convert $p (f_{τ}, 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.

Fig. 4

Flowchart of the proposed wave number domain algorithm (WDA).

As analyzed in Sec. 2, the time-shift $Δ T_{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

Eq. (20)

Ω_{2} (f_{τ}, x) = \exp {- j \frac{4 π f_{τ} x_{A} \sin φ}{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.

Fig. 5

Flowchart of block processing.

Moreover, after performing the azimuth FFT, the residual phase caused by azimuth deramping processing¹⁶ should be compensated for by

Eq. (21)

Ω_{3} (f_{a}) = \exp {j π \frac{f_{a}^{2}}{f_{r}}} .

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

Eq. (22)

Ω_{4} (f_{τ}, k_{x}) = \exp {j r \cdot \cos φ \sqrt{k_{rc}^{2} - {(k_{rc} \sin φ + k_{x})}^{2}} + \sin φ (k_{rc} \sin φ + k_{x})} \cdot \exp {- j π \frac{f_{τ}^{2}}{b}} \cdot \exp {- j \frac{4 π f_{τ} r_{\min}}{c}},

where

r_{\min}

is minimum slant range,

r_{\min} = r

–

c \cdot N_{r} / (2 f_{s})

, and

N_{r}

denotes the range resample point. In Eq. (22), the first exponential term compensates for the phase at the reference slant range as the 2-D matched filtering function

H_{r}

(

f_{τ}, k_{x}

) shows in Eq. (14). The second exponential term finishes the range compression and the third exponential term allows for redundant phase compensation, which are all caused by digital processing. After the bulk compression, the target at the reference range is properly focused, but a residual phase still exists for the targets at other ranges.

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

Eq. (23)

Ω_{5} (f_{τ}, x) = \exp {j 4 π f_{τ} \frac{\sin φ}{c} x} .

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 $ρ_{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 m^{2}$ , resulting in $5000 \times 5000$ pixels in the SAR image) is illustrated by Fig. 6.

Fig. 6

Simulation scene.

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 ( $φ = 40 \deg / 70 \deg$ ). Also, the range sample number in a conventional-SAR and ASRW-SAR are $N_{r}$ and $N_{r}^{'}$ , 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.

Fig. 7

Comparing the echo data of ASRW-SAR and conventional SAR systems at different squint angles. (a) conventional SAR, $φ = 40 \deg$ , Nr = 32768, (b) ASRW SAR, $φ = 40 \deg$ , Nr =8192, (c) conventional SAR, $φ = 70 \deg$ , Nr = 262144, and (d) ASRW SAR, $φ = 70 \deg$ , Nr = 16384.

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 $ρ_{a}$ , range resolution $ρ_{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 $γ = 4.1 / 17.0$ in $φ = 40 \deg / 70 \deg$ . The advantage of less-computation load indicates its potential suitability in a real-time imaging system.

Fig. 8

Comparing imaging results of conventional WDA and improved algorithm. (a) Conventional WDA. (b) Proposed imaging algorithm.

Fig. 9

Imaging results of the improved algorithm in $φ = 40 \deg$ . (a) Target A. (b) Target C. (c) Target E.

Fig. 10

Imaging results of the improved algorithm in $φ = 70 \deg$ . (a) Target A. (b) Target C. (c) Target E.

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$

$ρ_{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_{ken}$ . Also, the azimuth sample numbers in conventional SAR and ASRW-SAR are $N_{a}$ and $N_{a}^{'}$ , 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	$8 N_{a} N_{r} + 2 N_{a} N_{r} \log_{2} N_{a}$	$8 N_{a}^{'} N_{r}^{'} + 2 N_{a}^{'} N_{r}^{'} \log_{2} N_{a}^{'}$
Range FFT	$2 N_{a} N_{r} \log_{2} N_{r}$	$2 N_{a}^{'} N_{r}^{'} \log_{2} N_{r}^{'}$
Azimuth FFT	$2 N_{a} N_{r} \log_{2} N_{a}$	$2 N_{a}^{'} N_{r}^{'} \log_{2} N_{a}^{'}$
Bulk compression	$4 N_{a} N_{r}$	$4 N_{a}^{'} N_{r}^{'}$
Stolt mapping	$2 M_{ken} N_{a} N_{r}$	$2 M_{ken} N_{a}^{'} N_{r}^{'}$
Azimuth IFFT	$2 N_{a} N_{r} \log_{2} N_{a}$	$2 N_{a}^{'} N_{r}^{'} \log_{2} N_{a}^{'}$
Geometric correction	0	$4 N_{a}^{'} N_{r}^{'}$
Range IFFT	$2 N_{a} N_{r} \log_{2} N_{r}$	$2 N_{a}^{'} N_{r}^{'} \log_{2} N_{r}^{'}$

$C_{con}$ denotes the entire computational load for conventional-SAR.

Eq. (24)

C_{con} = (12 + 2 M_{ken}) N_{a} N_{r} + 6 N_{a} N_{r} \log_{2} N_{a} + 4 N_{a} N_{r} \log_{2} N_{r} .

Also, the entire computational load $C_{ASRW}$ for ASRW-SAR is expressed as

Eq. (25)

C_{ASRW} = (16 + 2 M_{ken}) N_{a}^{'} N_{r}^{'} + 6 N_{a}^{'} N_{r}^{'} \log_{2} N_{a}^{'} + 4 N_{a}^{'} N_{r}^{'} \log_{2} N_{r}^{'} .

Assume that the computational efficiency $γ$ is defined as

Eq. (26)

γ = \frac{C_{con}}{C_{ASRW}} .

With the condition of the same parameters as listed in Table 1, $N_{a}$ and $N_{a}^{'}$ are equivalent, $N_{a} = N_{a}^{'} = [(v \cdot PRF) / (ρ_{a} \cdot f_{r})]$ assuming an 8-point interpolation is adopted ( $M_{ken} = 8$ ). And in different squint angle $φ = 40 \deg / 70 \deg$ , $N_{a}$ (or $N_{a}^{'}$ ) is equal to 17,500/87,792. Also, the $N_{r}$ and $N_{r}^{'}$ are listed in Fig. 6 in different squint angles. Therefore, the computation efficiency $γ$ is calculated by Eq. (26) in different squint angles, and the final value $γ$ 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

1.

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

2.

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). http://dx.doi.org/10.1587/transcom.E92.B.3883 ITRCEC 0916-8516 Google Scholar

3.

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). http://dx.doi.org/10.1587/transcom.E94.B.2983 ITRCEC 0916-8516 Google Scholar

4.

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

5.

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

6.

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

7.

D. Realeet al., “Advanced techniques and new high resolution SAR sensors for monitoring urban areas,” 1800 –1803 2011). Google Scholar

8.

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

9.

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). http://dx.doi.org/10.1109/TGRS.2012.2183606 IGRSD2 0196-2892 Google Scholar

10.

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

11.

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

12.

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

13.

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

14.

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). http://dx.doi.org/10.1155/2007/59704 1110-8657 Google Scholar

15.

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

16.

D. GuoH. XuJ. Li, “Extended wavenumber domain algorithm for highly squinted sliding spotlight SAR data processing,” Prog. Electromagn. Res., 114 17 –32 (2011). PELREX 1043-626X Google 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.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

Wei Yang, Hong-cheng Zeng, Jie Chen, and Peng-bo Wang "Fast image-formation algorithm for ultrahigh-resolution airborne squint spotlight synthetic aperture radar based on adaptive sliding receive-window technique," Journal of Applied Remote Sensing 8(1), 083620 (4 June 2014). https://doi.org/10.1117/1.JRS.8.083620

Published: 4 June 2014

Access the abstract

JOURNAL ARTICLE
13 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 5 scholarly publications.

Explore citations on Lens.org

KEYWORDS

Synthetic aperture radar

Detection and tracking algorithms

Image processing

Doppler effect

Signal processing

Distortion

Antennas

1.

Introduction

2.

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

Fig. 1

Eq. (1)

Eq. (2)

Fig. 2

Table 1

Fig. 3

Eq. (3)

Eq. (4)

Eq. (5)

3.

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

3.1.

Two-Dimensional Matched Filtering Function

Eq. (6)

Eq. (7)

Eq. (8)

Eq. (9)

Eq. (10)

Eq. (11)

Eq. (12)

Eq. (13)

Eq. (14)

Eq. (15)

3.2.

New Stolt-Mapping Relationship

Eq. (16)

Eq. (17)

Eq. (18)

Eq. (19)

3.3.

Processing Scheme of Fast Image-Formation Algorithm

Fig. 4

Eq. (20)

Fig. 5

Eq. (21)

Eq. (22)

Eq. (23)

4.

Experimental Results

Fig. 6

Fig. 7

Fig. 8

Fig. 9

Fig. 10

Table 2

5.

Conclusion

Appendices

Appendix:

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

Table 3

Eq. (24)

Eq. (25)

Eq. (26)

Acknowledgments

References

Biography

Show All Keywords

Keywords/Phrases

Search In:

Publication Years