# A Bistatic Imaging Method for GEOSAR in the Strip Mode and UAVSAR in the Steering Beam Mode

Authors: W. Zhuo-qun; L. Ya-jun; S. Sheng; L. Shuang-shuang; X. Jin-guo; W. Jun-qiang; Z. Jian-hua
Journal: International Journal of Antennas and Propagation
 (2018)

---

## Abstract

The bistatic configuration with a geosynchronous orbital SAR (GEOSAR) transmitter and unmanned aerial vehicle SAR (UAVSAR) receiver can continuously image in any dangerous and interesting district. In this paper, the new imaging method in the case with the smaller orbital inclination of geosynchronous earth orbit and the steering beam working mode of UAVSAR was mainly studied and analyzed. GEOSAR can be approximately expressed as a static state, and only the receiver provides all the Doppler information. UAVSAR works in the steering beam modes, such as spotlight, sliding spotlight, and TOPS (Terrain Observation by Progressive Scan) mode. The azimuth bandwidth increased by the steering beam causes an aliasing situation in the azimuth frequency domain. To solve this problem, the proposed imaging method corrects the azimuth frequency aliasing using the scaling transform and the bulk azimuth compression. Compared with the traditional imaging method, the simulation validates perfectly the effectiveness of the bistatic imaging algorithm.


---

## Body

## 1. Introduction

Geosynchronous synthetic aperture radar (GEOSAR) is suitable for the continuous imaging in the specific partial region within the short revisit period of nearly 24 h [1–3]. However, the main difficulty for its hardware implementation is the large antenna and enough power of GEOSAR. The bistatic SAR configuration with a GEOSAR transmitter and unmanned aerial vehicle SAR receiver (GEO-UAV BiSAR) can reduce well the transmitted power and realize the imaging widely in any interesting and dangerous region [2].For the short synthetic aperture time of this bistatic radar, GEOSAR on the small orbital inclination (≤5°) may be approximately considered a static situation, and only the receiver contributes to the all azimuth modulations [4, 5]. The classical bistatic imaging algorithm cannot be employed to carry out this radar echo signal [6–9]. The imaging algorithm deals well with the data from the stationary transmitter and airborne receiver in the various working modes [10]. However, it is not considered that the azimuth bandwidth increased by antenna steering causes the aliasing in the azimuth frequency domain [11, 12]. To solve this problem, one way is to increase the pulse repetition frequency (PRF). However, the high PRF leads to range ambiguity and limits the processing speed of data acquired by SAR. The other way using a subaperture approach [13, 14] is proposed, but the general section principle cannot be confirmed well. Hence, the scaling transform and the bulk azimuth compression are introduced to handle the echoes of GEO-UAV BiSAR in this paper.This paper is organized as follows. In Section2, the characteristic of GEO-UAV BiSAR configuration and the unified signal model are described. Section 3 presents the bistatic imaging method which eliminates the azimuth aliasing through the scaling transform and the bulk azimuth compression. Section 4 shows the simulations which indicate the correctness and the effectiveness of the proposed imaging algorithm.


## 2. Uniform Signal Model of GEO-UAV BiSAR

When the beam steering of the receiver is in the working mode, the time-frequency relationship and the bistatic geometrical configuration are shown in Table1. In different working modes, the azimuth resolution ability is listed in the order TOPS < sliding spotlight < spotlight, whereas the surface coverage is listed in the reverse order. The particular bistatic geometric construction in the TOPS working mode is illustrated in Figure 1.Table 1
The geometrical configuration of the GEO-UAV BiSAR. The time-frequency of received signals is plotted from three-point targets in near range (a), center range (b), and far range (c).Bi represents the instantaneous signal bandwidth.
Mode & description
Time-frequency relationship
Geometric construction
Spotlight (VRg=0)
Azimuth resolution:λRcR/2LfR
Merits: this radar has the higher azimuth resolution than the other three modes
Defects: the surface swath is also smallest [15]
Sliding spotlight (0<VRg<VR)
Azimuth resolution:D/2μ
Merits: the receiver operated in a sliding spotlight mode achieves a tradeoff between the azimuth scene size and the azimuth resolution
Defects: as far as going studies in a large number of literatures are concerned, the less imaging method is appropriate [16]
TOPS (0<VR<VRg)
Azimuth resolution:D/2μ
Merits: the receiver in this type can switch the antenna to targets pertaining to different subswaths
Defects: this radar has the low azimuth resolution. A less imaging scheme can be applied [17]Figure 1
Geometrical configuration of GEO-UAV BiSAR (UAVSAR works in the TOPS mode).UAVSAR generates the virtual rotation point by antenna steering. The slant range between the receiver and the virtual rotation point isRRrot. VRg shows the ground velocity of the receiver’s antenna footprint. The beam velocity on the ground can be shown as VRg=μVR, where in the different working mode, μspot=0, μslid=RRrot−R0R/RRrot, and μtops=RRrot+R0R/RRrot. The received signal from the target Pm on the earth’s surface after demodulation are expressed as
(1)sstr,ta=rectta−ta0R/μLfR/VRrecttr−RTta0R,Pm+RRta,Pm/CTrexpjπKrtr−RTta0R,Pm+RRta,PmC2exp−j2πλRTta0R,Pm+RRta,Pm,where the first item rect⋅ represents the composite antenna pattern and shows the uniform illumination over the ground in any working mode of the receiver. The second item rect⋅ represents the pulse envelope. tr and ta are the fast time and the slow time, respectively. LfR is the length of the receiver’s flight path in the accumulated time. ta0R is the zero Doppler time of the receiver relative to the azimuth time origin ta=0. VR and C represent the platform velocity of the receiver and the speed of radio, respectively. λ and Fc are the wavelength and the carrier frequency of the transmitted LFM, respectively. Tr and Br show, respectively, the pulse duration and the bandwidth of transmitted LFM. Kr denotes the FM rate.R
T
t
a
,
P
m is the distance between the transmitter GEOSAR and the target Pm.
(2)RTta,Pm=R20T+VTsVTgta−ta0R2,where R0T denotes the slant range of the closest approach of the target to a virtual transmitter path which is assumed to be parallel with the receiver’s trajectory.When the orbital inclination is equal to 5°, the maximum of the approximated velocity is 22.07 m/s in the orbital period. The variation ofRTta,Pm calculated by (2) is less than λ/4 [18]. Hence, RTta,Pm is simplified to R0T.R
R
t
a
,
P
m represents the range from the receiver to the target.
(3)RRta,Pm=R0R2+VR2ta−ta0R2,where R0R is the slant range of the closest approach of the target to the receiver track.Using the principle of stationary phase (POSP), the echoes are redescribed by Fourier transform in range as(4)Ssfr,ta=rectta−ta0R/μLfR/VRrectfrBrexp−jπfr2Krexp−j2πFc+frCRTta0R,Pm+RRta,Pm.It is important to note that the transmitter-related term in (4) does not depend on the azimuth slow time ta but only on the fixed position of the point target. And all the Doppler information comes from the receiver platform.


## 3. Imaging Method for GEO-UAV BiSAR

Beam steering increases the azimuth bandwidth during the data acquisition. Thus, the classical methods are not suitable to manage the signal processing of GEO-UAV BiSAR. In this section, we show the new imaging technology which consists of three parts including resolving aliasing, the range compression and the RFM processing, and the azimuth deramping operation. The basic operations of the new algorithm in Figure2 are further illustrated in the following.Figure 2
Diagram of the GEO-UAV BiSAR imaging method.


### 3.1. Block 1: Resolving Aliasing in Azimuth Frequency

The azimuth bandwidth increased with antenna steering. Thus, the aliasing problem in the azimuth frequency domain needs to be urgently solved. In the first step, a convolution between the echo and the reference function is executed. The reference functionsrefta is given as
(5)srefta=expjΩπta2,Ω=1−μVR2λR0R.We simplify the calculation by the extension of the Bluenstein formula [19, 20].
(6)Ss′fr,ta′=Ssfr,ta⊗srefta=expjπΩta′2×∫Ssfr,zexpjπΩz2⋅exp−j2πΩtazdz.Equation (6) mainly contains two signal multiplications and a convolution. This course should be seen as a scaled inverse Fourier transform (SCFT−1) from the original coordinate ta to the new one ta′. And the relation between ta′ and fa is ta′=1/Ωfa. The signal Ss′fr,ta′ in the two-dimensional frequency domain can be redescribed as
(7)SS′fr,fa′=FFTSs′fr,ta′=SSfr,fa′Sreffa′.To obtain the echoes of the no-aliasing in the azimuth frequency domain, the signalSS′fr,fa′ is compensated by the function H1fa′ in (8).
(8)SS1′fr,fa′=SS′fr,fa′H1fa′,where H1fa′=expjπ1/Ωfa′2.The signalSS1′fr,fa′ without aliasing is obtained in (9), but the original PRF is increased through the above processing flow.
(9)SS1′fr,fa′=rectfa′−fdcBarectfrBrexp−jπfr2γexp2πfa′ta0Rexp−j2πFc+frCR0Texp−j2πR0RλFc+frFc2−λfa′VR2,where Br and Ba show, respectively, the range bandwidth and the Doppler bandwidth.


### 3.2. Block 2: Range Compression and RFM Processing

Range compression is performed for the signalSS1′fr,fa′, and the matched function H2fr is written as
(10)H2fr=expjπfr2γ.For the existence of a space variant phase (it stands for azimuth compression, RCMC, and secondary range compression), we need to implement the RFM processing to focus the data at a reference slant range. And the matched filter is expressed as(11)H3fr,fa′=expj2πRcRλFc+frFc2−λfa′VR2,where RcR is the slant range from the scene center to the receiver. RFM filtering correctly focuses the data at this reference slant range, partially compensating the phase of the target at other ranges. Through the steps H1:H3, the echoes can be described as
(12)SS2′fr,fa′=rectfa′−fdcBarectfrBrexp2πfa′ta0Rexp−j2πFc+frCRcTexp−j2πR0R−RcRCFc+fr2−Cfa′VR2.Here,fdc is the Doppler centroid. Fc+fr2−Cfa′/VR2 can be expanded by the Taylor expansion in the four-order form.
(13)Fc+fr2−Cfa′VR2≈Fcγfa′+frγfa′−1−γ2fa′2Fcγ3fa′⋅fr2+1−γ2fa′2Fc2γ5fa′⋅fr3,where γfa′=1−λfa′/VR2.On the right side of (13), the third term and the fourth term, respectively, stand for the high-order coupling of the two-dimensional frequency. These can be compensated fully by (14) and (15).
(14)H4fr,fa′=expjπR0R−RcR1−γ2fa′Cfc2γ5fa′⋅fr3,(15)H5fr,fa′=exp−jπR0R−RcR1−γ2fa′Cfcγ3fa′⋅fr2.Then, the echoes can be written as(16)SS3′fr,fa′=rectfa′−fdcBarectfrBrexp2πfa′ta0Rexp−j2πFc+frCR0Texp−j2πR0R−RcRCFcγfa′+frγfa′.


### 3.3. Block 3: Azimuth Deramping Operation

Based on the aforementioned imaging processing, the partial echoes are compensated at the no-reference slant range.SS3′fr,fa′ without azimuth frequency aliasing is gained. However, there still exists the aliasing in the azimuth time domain. Therefore, Block 3 mainly consisted of three components that include secondary RFM, the scaled transform, and the deramping operation.In the large coverage (especially TOPS mode), RFM processing can only make the most part of echoes at the reference slant range, and the rest of the echoes are compressed at the other range. Hence, in order to check all the echoes in the reference slant range, we further carry out the secondary RFM processing implemented by dividing the data into the small range blocks. The secondary RFM function is expressed as(17)H6fr,fa=expj2πRm−RcRCFcγfa′+frγfa′,where the subscript m stands for the index of the range blocks across the whole swath. And Rm is referred to as the reference slant range of the mth block (the midswath range in the block). Due to RFM errors, the azimuth broadening is less than 2%, and the focusing imaging results are obtained. The signal in the range time and azimuth frequency domain is written as
(18)SS4′tr,fa′=rectfa′−fdcBaexp2πfa′ta0Rwrtr−R0T+R0R−RmCexp−j2πλR0T+R0R−Rmγfa′.Then, the scaled transform is used to obtain the signal of no-aliasing in azimuth, including the rotation processingH7fa′ and the rotation transform H8fa′.
(19)H7fa′=expj2πλR0R−Rmγfa′,H8fa′=expjπ1−μΩfa′2exp−jπ1Ωfa′2.After azimuth IFFT, we obtain the signalss5′tr,ta′ without the aliasing in the two-dimensional time domain.
(20)ss5′tr,ta′=rectta′−κta0RLfR/VRexp−j2πλR0Twrtr−R0T+R0R−RmCexpjπΩVR2VR−VRgta′−ta0R2.Then, the deramping processing is implemented in this step, and the compensated function is defined as(21)H9ta′=exp−jπΩμta′2.Finally, after the azimuth Fourier transform, the echoes can be redescribed as(22)ss6′tr,ta′=wrtr−R0T+RcR−RmCwata′−ta0Rexp−j2πλR0T.


## 3.1. Block 1: Resolving Aliasing in Azimuth Frequency

The azimuth bandwidth increased with antenna steering. Thus, the aliasing problem in the azimuth frequency domain needs to be urgently solved. In the first step, a convolution between the echo and the reference function is executed. The reference functionsrefta is given as
(5)srefta=expjΩπta2,Ω=1−μVR2λR0R.We simplify the calculation by the extension of the Bluenstein formula [19, 20].
(6)Ss′fr,ta′=Ssfr,ta⊗srefta=expjπΩta′2×∫Ssfr,zexpjπΩz2⋅exp−j2πΩtazdz.Equation (6) mainly contains two signal multiplications and a convolution. This course should be seen as a scaled inverse Fourier transform (SCFT−1) from the original coordinate ta to the new one ta′. And the relation between ta′ and fa is ta′=1/Ωfa. The signal Ss′fr,ta′ in the two-dimensional frequency domain can be redescribed as
(7)SS′fr,fa′=FFTSs′fr,ta′=SSfr,fa′Sreffa′.To obtain the echoes of the no-aliasing in the azimuth frequency domain, the signalSS′fr,fa′ is compensated by the function H1fa′ in (8).
(8)SS1′fr,fa′=SS′fr,fa′H1fa′,where H1fa′=expjπ1/Ωfa′2.The signalSS1′fr,fa′ without aliasing is obtained in (9), but the original PRF is increased through the above processing flow.
(9)SS1′fr,fa′=rectfa′−fdcBarectfrBrexp−jπfr2γexp2πfa′ta0Rexp−j2πFc+frCR0Texp−j2πR0RλFc+frFc2−λfa′VR2,where Br and Ba show, respectively, the range bandwidth and the Doppler bandwidth.


## 3.2. Block 2: Range Compression and RFM Processing

Range compression is performed for the signalSS1′fr,fa′, and the matched function H2fr is written as
(10)H2fr=expjπfr2γ.For the existence of a space variant phase (it stands for azimuth compression, RCMC, and secondary range compression), we need to implement the RFM processing to focus the data at a reference slant range. And the matched filter is expressed as(11)H3fr,fa′=expj2πRcRλFc+frFc2−λfa′VR2,where RcR is the slant range from the scene center to the receiver. RFM filtering correctly focuses the data at this reference slant range, partially compensating the phase of the target at other ranges. Through the steps H1:H3, the echoes can be described as
(12)SS2′fr,fa′=rectfa′−fdcBarectfrBrexp2πfa′ta0Rexp−j2πFc+frCRcTexp−j2πR0R−RcRCFc+fr2−Cfa′VR2.Here,fdc is the Doppler centroid. Fc+fr2−Cfa′/VR2 can be expanded by the Taylor expansion in the four-order form.
(13)Fc+fr2−Cfa′VR2≈Fcγfa′+frγfa′−1−γ2fa′2Fcγ3fa′⋅fr2+1−γ2fa′2Fc2γ5fa′⋅fr3,where γfa′=1−λfa′/VR2.On the right side of (13), the third term and the fourth term, respectively, stand for the high-order coupling of the two-dimensional frequency. These can be compensated fully by (14) and (15).
(14)H4fr,fa′=expjπR0R−RcR1−γ2fa′Cfc2γ5fa′⋅fr3,(15)H5fr,fa′=exp−jπR0R−RcR1−γ2fa′Cfcγ3fa′⋅fr2.Then, the echoes can be written as(16)SS3′fr,fa′=rectfa′−fdcBarectfrBrexp2πfa′ta0Rexp−j2πFc+frCR0Texp−j2πR0R−RcRCFcγfa′+frγfa′.


## 3.3. Block 3: Azimuth Deramping Operation

Based on the aforementioned imaging processing, the partial echoes are compensated at the no-reference slant range.SS3′fr,fa′ without azimuth frequency aliasing is gained. However, there still exists the aliasing in the azimuth time domain. Therefore, Block 3 mainly consisted of three components that include secondary RFM, the scaled transform, and the deramping operation.In the large coverage (especially TOPS mode), RFM processing can only make the most part of echoes at the reference slant range, and the rest of the echoes are compressed at the other range. Hence, in order to check all the echoes in the reference slant range, we further carry out the secondary RFM processing implemented by dividing the data into the small range blocks. The secondary RFM function is expressed as(17)H6fr,fa=expj2πRm−RcRCFcγfa′+frγfa′,where the subscript m stands for the index of the range blocks across the whole swath. And Rm is referred to as the reference slant range of the mth block (the midswath range in the block). Due to RFM errors, the azimuth broadening is less than 2%, and the focusing imaging results are obtained. The signal in the range time and azimuth frequency domain is written as
(18)SS4′tr,fa′=rectfa′−fdcBaexp2πfa′ta0Rwrtr−R0T+R0R−RmCexp−j2πλR0T+R0R−Rmγfa′.Then, the scaled transform is used to obtain the signal of no-aliasing in azimuth, including the rotation processingH7fa′ and the rotation transform H8fa′.
(19)H7fa′=expj2πλR0R−Rmγfa′,H8fa′=expjπ1−μΩfa′2exp−jπ1Ωfa′2.After azimuth IFFT, we obtain the signalss5′tr,ta′ without the aliasing in the two-dimensional time domain.
(20)ss5′tr,ta′=rectta′−κta0RLfR/VRexp−j2πλR0Twrtr−R0T+R0R−RmCexpjπΩVR2VR−VRgta′−ta0R2.Then, the deramping processing is implemented in this step, and the compensated function is defined as(21)H9ta′=exp−jπΩμta′2.Finally, after the azimuth Fourier transform, the echoes can be redescribed as(22)ss6′tr,ta′=wrtr−R0T+RcR−RmCwata′−ta0Rexp−j2πλR0T.


## 4. Simulation Analysis

To verify the correctness of the proposed imaging method, we use the data collection geometry illustrated in Figure3 to perform the simulation.Figure 3
The composite scene with 9-point targets. Here,DX and DY represent the azimuth and the range resolution, respectively.


### 4.1. Imaging Simulation for the Single-Point Target

We use the proposed method and the old imaging method [10] to handle the echoes from the point target P5 in the spotlight mode of the receiver, respectively.The simulation parameters are listed in Table2. GEOSAR works at X-band. Under the condition of the LFM bandwidth = 80 MHz, the theoretical range resolution attains 1.86 m. PRF is set to 120% of the instantaneous bandwidth, i.e., PRF = 1555 Hz. And UAV works with the velocity 200 m/s in the spotlight mode. We, respectively, use the proposed method and the method referenced in [10] to handle the echoes from the point target P5 located on the position coordinate (0 m, 0 m, and 0 m). In the receiver’s spotlight mode, the azimuth resolution reaches 0.15 m.Table 2
Simulation parameters.
Parameter
Value
GEOSAR altitude
36,000 km
Wavelength
0.03 m (X-band)
Bandwidth
80 MHz
Pulse duration
20μs
UAVSAR altitude
8000 m
Antenna size of UAVSAR
2 m
UAVSAR velocity
200 m/s
Azimuth aperture angle in the case of the receiver’s spotlight mode
12°
The beam velocity of UAVSAR in the sliding spotlight mode
140 m/s
The beam velocity of UAVSAR in the TOPS mode
270 m/sThe echo signal of the single-point target is managed, respectively, using the imaging method referenced in [10] and the proposed imaging algorithm. The imaging results are shown in Figures 4 and 5. Using the imaging algorithm referenced in [10], the aliasing situation in the azimuth frequency spectrum is more obvious. The imaging results show the three-point targets in the scene. However, in fact, all echoes only come from a real point target P5, and the two points on the edge of Figure 4(b) are illusory. Hence, the real location of the point target cannot be clearly known by this method.Figure 4
Imaging results by the traditional method in literature [10].
(a)
Two-dimensional frequency spectrum of the point target
(b)
Imaging results of the point targetFigure 5
Imaging result by the proposed imaging method.
(a)
Two-dimensional frequency spectrum of the point target
(b)
Imaging results of the point targetThen, using the proposed imaging algorithm in this paper, the two-dimensional frequency spectrum with no aliasing is obtained. The results indicate a clear point target in the center of the scene. It is the reason why the new imaging method perfectly solves the azimuth frequency aliasing for antenna steering. Compared with the imaging algorithm referenced in [10], the proposed approach reduces the instantaneous azimuth bandwidth that increased with antenna steering and avoids the azimuth aliasing issue well. To further confirm the effectiveness of the new method quantitatively, imaging results of the point target will be evaluated by three indexes including the peak sidelobe ratio (PSLR), the integration sidelobe ratio (ISLR), and the impulse response width (IRW). In the range direction, PSLR = −13.25 dB, ISLR = −9.45 dB, and IRW = 1.95 m, and in azimuth, PSLR = −13.27 dB, ISLR = −9.48 dB, and IRW = 0.15 m. It can be clearly found that the imaging results are close to the theoretical values. For the receiver’s various working mode, the radar in the spotlight mode has the larger increments of the azimuth bandwidth than that in the other mode. Hence, this new method is appropriate for GEO-UAV BiSAR in the receiver’s various working mode.


### 4.2. Imaging Simulation for the Multipoint Targets

To present the validity of the proposed imaging method, we perform the imaging simulation of the multipoint targets in the various modes of the receiver. PRF is 20% greater than the instantaneous azimuth bandwidth in the simulation, andPRFspot=1555Hz, PRFsliding=171Hz, and PRFstrip=88Hz. Based on the simulated parameters listed in Table 2, we obtain the ideal range resolution to 1.87 m. In UAVSAR’s various working modes, the theoretical azimuth resolution follows DXspot=0.15m, DXsliding=1.40m, and DXTOPS=2.70m.The received signals from the multipoint targets are operated by the proposed imaging algorithm. Figures6–8(a) mainly show the imaging results of 9-point targets in the steering antenna working mode. And Figures 6, 7, and 8(b)–8(d), respectively, describe the two-dimensional profiles of the single-point P3. Subsequently, the imaging qualities of the partial point targets, i.e., PSLR, ISLR, and IRW, are revealed in Table 3. When UAVSAR works in the spotlight mode, the targets of P3 and P7 had a maximum expansion of 4.28% in range and 6.67% in azimuth. In the sliding spotlight working mode of the receiver, P5, located at the screen center, has the good imaging quality with 1.87 m in range and 1.40 m in azimuth. And in the receiver’s TOPS mode, P5 possesses an ideal range impulse response and the azimuth impulse response expands to 2.96%. Then, it is seen that, in the various working modes of the receiver, the reference target and its neighbors are well focused in the imaging results of the scene. The experiment demonstrates that the proposed algorithm can be perfectly applied to processing the GEO-UAV BiSAR echoes.Figure 6
Imaging results in the spotlight mode of the receiver.
(a)
Multipoint imaging results
(b)
Imaging result of P3
(c)
Range response profile of P3
(d)
Azimuth response profile of P3Figure 7
Imaging results in the sliding spotlight mode of the receiver.
(a)
Multipoint imaging results
(b)
Imaging result of P3
(c)
Range response profile of P3
(d)
Azimuth response profile of P3Figure 8
Imaging results in the TOPS mode of the receiver.
(a)
Multipoint imaging results
(b)
Imaging result of P3
(c)
Range response profile of P3
(d)
Azimuth response profile of P3Table 3
Imaging qualities of the multipoint targets (P1, P5, and P9).
Receiver mode
Location of point targets
Range
Azimuth
PSLR (dB)
ISLR (dB)
Resolution (m)
PSLR (dB)
ISLR (dB)
Resolution (m)
Spotlight
P3
−13.27
−9.26
1.95
−13.79
−9.38
0.16
P5
−13.25
−9.45
1.95
−13.27
−9.52
0.15
P7
−13.32
−9.47
1.95
−13.76
−9.54
0.16
Sliding spotlight
P3
−13.15
−9.52
1.87
−13.54
−9.45
1.54
P5
−13.23
−9.51
1.87
−13.26
−9.63
1.40
P7
−12.89
−9.33
1.93
−12.85
−9.42
1.55
TOPS
P3
−13.22
−9.64
1.87
−13.34
−9.45
2.81
P5
−13.01
−9.26
1.87
−13.21
−9.43
2.78
P7
−12.84
−9.62
1.87
−13.10
−9.75
2.80


## 4.1. Imaging Simulation for the Single-Point Target

We use the proposed method and the old imaging method [10] to handle the echoes from the point target P5 in the spotlight mode of the receiver, respectively.The simulation parameters are listed in Table2. GEOSAR works at X-band. Under the condition of the LFM bandwidth = 80 MHz, the theoretical range resolution attains 1.86 m. PRF is set to 120% of the instantaneous bandwidth, i.e., PRF = 1555 Hz. And UAV works with the velocity 200 m/s in the spotlight mode. We, respectively, use the proposed method and the method referenced in [10] to handle the echoes from the point target P5 located on the position coordinate (0 m, 0 m, and 0 m). In the receiver’s spotlight mode, the azimuth resolution reaches 0.15 m.Table 2
Simulation parameters.
Parameter
Value
GEOSAR altitude
36,000 km
Wavelength
0.03 m (X-band)
Bandwidth
80 MHz
Pulse duration
20μs
UAVSAR altitude
8000 m
Antenna size of UAVSAR
2 m
UAVSAR velocity
200 m/s
Azimuth aperture angle in the case of the receiver’s spotlight mode
12°
The beam velocity of UAVSAR in the sliding spotlight mode
140 m/s
The beam velocity of UAVSAR in the TOPS mode
270 m/sThe echo signal of the single-point target is managed, respectively, using the imaging method referenced in [10] and the proposed imaging algorithm. The imaging results are shown in Figures 4 and 5. Using the imaging algorithm referenced in [10], the aliasing situation in the azimuth frequency spectrum is more obvious. The imaging results show the three-point targets in the scene. However, in fact, all echoes only come from a real point target P5, and the two points on the edge of Figure 4(b) are illusory. Hence, the real location of the point target cannot be clearly known by this method.Figure 4
Imaging results by the traditional method in literature [10].
(a)
Two-dimensional frequency spectrum of the point target
(b)
Imaging results of the point targetFigure 5
Imaging result by the proposed imaging method.
(a)
Two-dimensional frequency spectrum of the point target
(b)
Imaging results of the point targetThen, using the proposed imaging algorithm in this paper, the two-dimensional frequency spectrum with no aliasing is obtained. The results indicate a clear point target in the center of the scene. It is the reason why the new imaging method perfectly solves the azimuth frequency aliasing for antenna steering. Compared with the imaging algorithm referenced in [10], the proposed approach reduces the instantaneous azimuth bandwidth that increased with antenna steering and avoids the azimuth aliasing issue well. To further confirm the effectiveness of the new method quantitatively, imaging results of the point target will be evaluated by three indexes including the peak sidelobe ratio (PSLR), the integration sidelobe ratio (ISLR), and the impulse response width (IRW). In the range direction, PSLR = −13.25 dB, ISLR = −9.45 dB, and IRW = 1.95 m, and in azimuth, PSLR = −13.27 dB, ISLR = −9.48 dB, and IRW = 0.15 m. It can be clearly found that the imaging results are close to the theoretical values. For the receiver’s various working mode, the radar in the spotlight mode has the larger increments of the azimuth bandwidth than that in the other mode. Hence, this new method is appropriate for GEO-UAV BiSAR in the receiver’s various working mode.


## 4.2. Imaging Simulation for the Multipoint Targets

To present the validity of the proposed imaging method, we perform the imaging simulation of the multipoint targets in the various modes of the receiver. PRF is 20% greater than the instantaneous azimuth bandwidth in the simulation, andPRFspot=1555Hz, PRFsliding=171Hz, and PRFstrip=88Hz. Based on the simulated parameters listed in Table 2, we obtain the ideal range resolution to 1.87 m. In UAVSAR’s various working modes, the theoretical azimuth resolution follows DXspot=0.15m, DXsliding=1.40m, and DXTOPS=2.70m.The received signals from the multipoint targets are operated by the proposed imaging algorithm. Figures6–8(a) mainly show the imaging results of 9-point targets in the steering antenna working mode. And Figures 6, 7, and 8(b)–8(d), respectively, describe the two-dimensional profiles of the single-point P3. Subsequently, the imaging qualities of the partial point targets, i.e., PSLR, ISLR, and IRW, are revealed in Table 3. When UAVSAR works in the spotlight mode, the targets of P3 and P7 had a maximum expansion of 4.28% in range and 6.67% in azimuth. In the sliding spotlight working mode of the receiver, P5, located at the screen center, has the good imaging quality with 1.87 m in range and 1.40 m in azimuth. And in the receiver’s TOPS mode, P5 possesses an ideal range impulse response and the azimuth impulse response expands to 2.96%. Then, it is seen that, in the various working modes of the receiver, the reference target and its neighbors are well focused in the imaging results of the scene. The experiment demonstrates that the proposed algorithm can be perfectly applied to processing the GEO-UAV BiSAR echoes.Figure 6
Imaging results in the spotlight mode of the receiver.
(a)
Multipoint imaging results
(b)
Imaging result of P3
(c)
Range response profile of P3
(d)
Azimuth response profile of P3Figure 7
Imaging results in the sliding spotlight mode of the receiver.
(a)
Multipoint imaging results
(b)
Imaging result of P3
(c)
Range response profile of P3
(d)
Azimuth response profile of P3Figure 8
Imaging results in the TOPS mode of the receiver.
(a)
Multipoint imaging results
(b)
Imaging result of P3
(c)
Range response profile of P3
(d)
Azimuth response profile of P3Table 3
Imaging qualities of the multipoint targets (P1, P5, and P9).
Receiver mode
Location of point targets
Range
Azimuth
PSLR (dB)
ISLR (dB)
Resolution (m)
PSLR (dB)
ISLR (dB)
Resolution (m)
Spotlight
P3
−13.27
−9.26
1.95
−13.79
−9.38
0.16
P5
−13.25
−9.45
1.95
−13.27
−9.52
0.15
P7
−13.32
−9.47
1.95
−13.76
−9.54
0.16
Sliding spotlight
P3
−13.15
−9.52
1.87
−13.54
−9.45
1.54
P5
−13.23
−9.51
1.87
−13.26
−9.63
1.40
P7
−12.89
−9.33
1.93
−12.85
−9.42
1.55
TOPS
P3
−13.22
−9.64
1.87
−13.34
−9.45
2.81
P5
−13.01
−9.26
1.87
−13.21
−9.43
2.78
P7
−12.84
−9.62
1.87
−13.10
−9.75
2.80


## 5. Conclusions

We have introduced a uniform bistatic imaging algorithm for processing the echoes from the GEO-UAV BiSAR in the antenna steering working mode of the receiver (i.e., spotlight, sliding spotlight, and TOPS). The transmitter with GEOSAR on the small inclination is confirmed as a quiescent state, and UAVSAR offers all Doppler information. That Doppler bandwidth increased by antenna steering generates the aliasing problem in the azimuth frequency domain making the classical bistatic imaging methods inappropriate for this radar. The proposed imaging algorithm mainly includes the scale transform and bulk azimuth compression, the range compression and the RFM processing, the scaling transform and deramping approaches. In comparison with the imaging method referenced in [10], the proposed algorithm verifies its correctness and effectiveness via simulations. It solves the aliasing frequency problem and well focuses the bistatic data on the real position of the point targets in the large scene.In our future work, we will further study the bistatic imaging method for GEOSAR on the large orbital inclination.



---

*Source: 6138709-2018-10-31.xml*