Next Article in Journal
Low-Latency and Energy-Efficient Data Preservation Mechanism in Low-Duty-Cycle Sensor Networks
Previous Article in Journal
Smartphone-Based Escalator Recognition for the Visually Impaired
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Generalized Chirp-Scaling Algorithm for Geosynchronous Orbit SAR Staring Observations

1
School of Electronics and Information, Northwestern Polytechnical University, N 1, Dongxiang Rd, Xi’an 710129, China
2
China Academy of Space Technology (Xi’an), Weiqu Street, Xi’an 710100, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(5), 1058; https://doi.org/10.3390/s17051058
Submission received: 19 March 2017 / Revised: 30 April 2017 / Accepted: 2 May 2017 / Published: 6 May 2017
(This article belongs to the Section Remote Sensors)

Abstract

:
Geosynchronous Orbit Synthetic Aperture Radar (GEO SAR) has recently received increasing attention due to its ability of performing staring observations of ground targets. However, GEO SAR staring observation has an ultra-long integration time that conventional frequency domain algorithms cannot handle because of the inaccurately assumed slant range model and existing azimuth aliasing. To overcome this problem, this paper proposes an improved chirp-scaling algorithm that uses a fifth-order slant range model where considering the impact of the “stop and go” assumption to overcome the inaccuracy of the conventional slant model and a two-step processing method to remove azimuth aliasing. Furthermore, the expression of two-dimensional spectrum is deduced based on a series of reversion methods, leading to an improved chirp-scaling algorithm including a high-order-phase coupling function compensation, range and azimuth compression. The important innovations of this algorithm are implementation of a fifth-order order slant range model and removal of azimuth aliasing for GEO SAR staring observations. A simulation of an L-band GEO SAR with 1800 s integration time is used to demonstrate the validity and accuracy of this algorithm.

1. Introduction

Owing to its ability of all-weather and all-day Earth observation, synthetic aperture radar (SAR) has been widely applied in a variety of microwave remote sensing fields, such as disaster monitoring [1], topographic mapping [2], soil moisture estimation [3], oil spill observation [4], crop growth monitoring [5], forest cover measurement [6], wetland mapping [7], maritime surveillance [8,9,10], and so on. However, with increasing demands of wider imaging swath and shorter revisit period, it is difficult for low Earth orbit (LEO) SAR to satisfy the application requirements. To solve this problem, the Geosynchronous Orbit SAR (GEO SAR) concept has been put forward.
Compared with LEO SAR, GEO SAR not only has more advantages, such as wider swath, and shorter revisiting period [11,12,13,14,15,16], but it also can realize long time observation of ground targets even with staring, which will be beneficial to its applications in microwave remote sensing fields. Staring observation is the classical spotlight mode with a steering to a rotation center within the imaged scene [17]. GEO SAR’s staring observation can be obtained using a certain beam forming strategy that could be implemented by adjusting the antenna beam direction, or through satellite attitude steering [18].
The first challenge in GEO SAR staring observation image formation is the integration time which is extended to thousands of seconds in comparison with the hundreds of seconds in GEO SAR conventional observation using strip-map mode. In this case, conventional GEO SAR frequency domain algorithms, such as the modified chirp scaling (CS) algorithm [19,20], improved frequency domain focusing method [21,22], generalized Omega-K algorithm [23], and sub-aperture processing methods [24], cannot be directly used because of the inaccuracy of the corresponding slant models. In above algorithms, the slant range model based on the fourth model fails to precisely describe the range history of the ultra-long integration time in the staring observation. Zhao et al. [25] have proposed an improved RD imaging algorithm based on a fifth Doppler range model that is mainly used for strip-map SAR. Moreover, the Doppler range model in [25] does not consider the impact of “stop and go” assumptions. Yu et al. [26] have proposed a time-frequency scaling algorithm to correct the high-order azimuth variance in GEO SAR, but the imaging mode is still in strip-map mode. Conventional time domain processing algorithms, such as the back projection algorithm [27], are not restricted by the slant range model. However, due to the prohibitively heavy computation load and inefficient processing, they are not practical to use. Besides, in the thousand seconds of integration time, many factors would affect the imaging focus, such as atmospheric refractive index, orbit perturbation, etc. Several methods have been proposed to solve these problems. Ruiz et al. [28] proposed the method of retrieving the atmospheric phase screen (APS) and compensating its effects on GEO SAR focusing. To improve GEO SAR focusing influenced by perturbations, Dong et al. [29] recommended the use of accurate orbit measurements or some signal processing methods. The autofocus technique for phase error correction has been proven to be a robust technique that can compensate for the phase error caused by atmosphere or orbit perturbations [30]. Although the need for deeper understanding of the effect on GEO SAR imaging, this paper will focus on imaging processes.
The other challenge encountered in GEO SAR staring observation of image formation is the azimuth aliasing as the spotlight SAR. Therefore, an imaging algorithm for GEO SAR staring observation should address the issues. This paper aims to study the main problems posed by above challenges.
The main contributions of this paper are as follows: first, a fifth-order slant range model considering the impact of the “stop and go” assumption is developed for GEO SAR staring observations. Then, an accurate 2-D spectral based on series inversions is deduced. And a two-step processing method is proposed to remove spectrum aliasing. At last, an improved chirp scaling algorithm based on fifth-order slant range is obtained to perform high-order-phase coupling function compensation, range and azimuth compression.
This paper is organized as follows: the signal model based on the fifth-order slant range for GEO SAR is investigated in Section 2. In Section 3, a two-dimensional spectrum of the echo is derived. In Section 4, the imaging algorithm is deduced in detail and the algorithm flow chart is given. Section 5 gives the simulation results to validate the improved imaging algorithm. Finally, conclusions are summarized in Section 6.

2. Signal Model for GEO SAR Staring Observation

The geometry of satellite and ground targets for GEO SAR staring observation is shown in Figure 1.
In the case of staring observations, the antenna (moving from position A to position B) is constantly steered to the point target which results in a long illumination time. Assuming the position vector of the satellite and the target vector in the Earth fixed coordinate system to be R S ( t a ) and R T ( t a ) respectively, and then the slant range from the satellite to the target point can be expressed as R ( t a ) = R S ( t a ) R T ( t a ) , where t a is the azimuth time.
Taking a group of typical GEO SAR parameters, e.g., the semi-major axis of the orbit is 42,164 km and the orbit eccentricity is 0, the orbit inclination is 20° and the azimuth time is 1800s. From Figure 2, the maximal phase error caused by the conventional fourth-order slant range model is greater than π/4 whereas the fifth-order slant range model is far less than π/4. Nevertheless, the fourth-order slant range model cannot meet the requirement of the imaging algorithm. Consequently, the fifth-order range model is highly desired. Under the rule that phase error is less than π/4, we can get the relationship between different order slant range model and synthetic aperture time, which is shown in Figure 3.
As can be seen from the Figure 3, the third-order slant range model can be used when the synthetic aperture time does not exceed 680 s. The fourth-order slant distance model can be used when not more than 1580 s, and so on. However, when the synthetic aperture time exceeds 2830 s, the sixth-order slant range model must be used. In the 1800 s of synthetic aperture time, the fifth diagonal model can meet the requirements. Using Taylor series expansion in azimuth time, the slant range is extended to a fifth-order model by following equation:
R ( t a ) = R c + V c t a + 1 2 A c t a 2 + 1 6 B c t a 3 + 1 24 C c t a 4 + 1 120 D c t a 5 +
where, the expression of R c ,   V c ,   A c ,   B c ,   C c can be obtained in [31]. The value of D c can be further expressed as:
D c = ( 15 x 1 5 / 32 75 x 1 3 x 2 / 4 + 45 x 3 x 1 2 / 2 + 45 x 1 x 2 2 / 2 30 x 4 x 1 30 x 3 x 2 + 60 x 5 ) R S T
in which:
R S T = | R S T | ,
x 1 = 2 R S T , V S T R S T 2 , x 2 = R S T , A S T + V S T , V S T R S T 2 , x 3 = R S T , B S T / 3 + V S T , A S T R S T 2 , x 4 = R S T , C S T / 12 + V S T , B S T / 3 + A S T , A S T / 4 R S T 2 , x 5 = R S T , D S T / 60 + V S T , C S T / 12 + A S T , B S T / 6 R S T 2 , R S T = R S R T ,   V S T = V S V T ,   A S T = A S A T ,   B S T = B S B T ,   C S T = C S C T ,   D S T = D S D T
where R ,   V ,   A ,   B ,   C are all state vectors and represent the satellite and target positions vector, velocity vector, acceleration vector, rate of the acceleration vector, and the rate of the acceleration vector, respectively. The subscript S and T denote position and target, respectively. Operator , represents the inner product operation. The expression of D S and D T can be obtained with the following equations:
D S = μ 2 V S / | R S | 6   3   μ 2 R S , V S R S / | R S | 5 +   9   μ A S R S , V S / | R S | 5 +   9   μ V S , V S V S / | R S | 5 +   9   μ R S , A S V S / | R S | 5 +   9   μ V S , A S R S / | R S | 5   45 μ R S , V S 2 V S / | R S | 7   45   μ R S , V S V S , R S R S / | R S | 7   45   μ A S , R S R S , V S R S / | R S | 7   3   μ 2 R S , V S R S / | R S | 8 +   9   μ 2 R S , V S R S / | R S | 8 +   105   μ R S , V S 3 R S / | R S | 9
D T = W e C T
where W e is angular velocity vector of the Earth rotation, μ is the constant of gravitation and operator ⊗ stands for the outer product.
In the slant range, because of the ultra-long integration time, the effect of perturbing orbital elements on the satellite and target kinematic parameters should be taken into account [29]. As indicated in [31], the satellite and target kinematic parameters that can be obtained from the orbital elements include the following: a , i.e., the semimajor axis of the orbit, e , i.e., the eccentricity of the orbit, i , i.e., the inclination of the orbit, and Ω, i.e., the longitude of the ascending node, f , i.e., the true anomaly, w , i.e., the argument of perigee. Ignoring the slight stochastic perturbation errors, the osculating orbital elements can be expressed as:
{ a = a m + ( Δ a ) s e c + ( Δ a ) l p + ( Δ a ) s p e = e m + ( Δ e ) s e c + ( Δ e ) l p + ( Δ e ) s p i = i m + ( Δ i ) s e c + ( Δ i ) l p + ( Δ i ) s p Ω = Ω m + ( Δ Ω ) s e c + ( Δ Ω ) l p + ( Δ Ω ) s p f = f m + ( Δ f ) s e c + ( Δ f ) l p + ( Δ f ) s p w = w m + ( Δ w ) s e c + ( Δ w ) l p + ( Δ w ) s p
where a m ,   e m ,   i m ,   Ω m ,   f m ,   w m are the mean orbital elements; Δ a ,   Δ e ,   Δ i ,   Δ Ω ,   Δ f ,   Δ w are the perturbing orbital elements; and the subscripts “sec”, “lp”, and “sp” denote the secular perturbation term, the longperiod term, and the short-period term, respectively. Therefore, by computation of the orbit perturbation parameters and the accurate measurement of the orbit, we can obtain the slant distance accurately.
In the conventional SAR echo signal model, the transmission point and receiving point can be considered to be at the same position. However, for GEO SAR, the orbit height is about 36,000 km, so this assumption is invalid [22]. To facilitate deduction and use of the algorithm, the impact of “stop-and-go” assumption must be compensated. Subsequently the imaging algorithm can be expressed in the form of monostatic SAR.
The error caused by the “stop-and-go” assumption on the slant range can be expressed as follows [22]:
Δ R ( t a ) = Δ R c + Δ V c t a + 1 2 Δ A c t a 2 + 1 6 Δ B c t a 3 +
where Δ R c = V S T , R S T c , Δ A c = B S T , R S T + 3 A S T , V S T 2 c , Δ B c = C S T , R S T + 4 B S T , V S T 6 c + R S T , R S T 2 c , and c is the speed of light.
The reached range error using the same orbit parameters as in Section 2 and the third-order approximation is 0.416 mm, leading to a phase error far less than π/4. Thus, the error distance caused by the “stop and go” assumption can be fully derived from the third-order approximation. Taking the error caused by the “stop and go” assumption into account, the slant range in GEO SAR can be written as:
R ( t a ) R c + V c t a + 1 2 V c t a 2 + 1 6 B c t a 3 + 1 24 C t c a 4 + 1 120 D c t a 5
where R c = R c + Δ R c , V c = V c + Δ V c , A c = A c + Δ A c , B c = B c + Δ B c .
According to the SAR principle, the echo signal after demodulation can be expressed as:
s 0 ( t r , t a ) = w r ( t r 2 R ( t a ) c ) w a ( t a ) exp ( j 4 π f 0 R ( t a ) c ) exp ( j π K r ( t r 2 R ( t a ) c ) 2 )
where t r is the range time, t a is the azimuth time, w r and w a stand for range envelope and azimuth envelope, respectively; f 0 denotes carrier frequency, K r is the frequency modulate rate. R ( t a ) is the slant range depending on azimuth time.

3. Two-Dimensional Spectrum Analysis

The two-dimensional spectrum analysis plays an important role in frequency domain algorithms, and the improved imaging algorithm can be deduced from which. If the fast Fourier transform FFT is performed in the range domain, the echo signal can be written as:
s 1 ( f r , t a ) = W r ( f r ) w a ( t a ) exp ( j 4 π ( f 0 + f r ) R ( t a ) c ) exp ( j π f r 2 K r )
where f r is the range frequency.
Then, we perform the azimuth FFT and using the stationary phase principle and series inversion principle. The signal’s two-dimensional frequency domain expression can be obtained as:
Θ ( f r , f a ) π f r 2 K r 4 π ( f c + f r ) c { R c A 1 2 [ c f a 2 ( f c + f r ) V c ] 2 A 2 3 [ c f a 2 ( f c + f r ) V c ] 3 A 3 4 [ c f a 2 ( f c + f r ) V c ] 4 A 4 5 [ c f a 2 ( f c + f r ) V c ] 5 }
where f a is azimuth frequency, and:
{ A 1 = 1 A c A 2 = B c A c 3 A 3 = 3 B c 2 A c C c 6 A c 5 A 4 = B c 3 2 A c 7 + 7 B c C c 12 A c 6 D c A c 5
because f r f c is far less than 1, so the series expansion can be represented in terms of ( 1 + f r f c ) 1 , ( 1 + f r f c ) 2 and ( 1 + f r f c ) 3 . With this substitution, the two-dimensional spectral phase can be rewritten as follows:
Θ ( f r , f a ) φ 0 ( f a ) + φ 1 ( f a ) f r + φ 2 ( f a ) f r 2 + φ 3 ( f a ) f r 3
where φ 0 ( f a ) is the azimuth phase, which is irrelevant to range frequency, can be expressed as:
φ 0 ( f a ) = 4 π λ [ R c + λ 2 A 1 8 ( f a f d ) 2 λ 3 A 2 24 ( f a f d ) 3 + λ 4 A 3 64 ( f a f d ) 4 λ 5 A 4 160 ( f a f d ) 5 ]
where f d = 2 V c λ , φ 1 ( f a ) is the first order coupling phase describing RCM, which is relevant on range frequency can be expressed as:
φ 1 ( f a ) = 4 π c { R c + λ 2 A 1 8 [ ( f a f d ) 2 + 2 f d ( f a f d ) ] λ 3 A 2 24 [ 2 ( f a f d ) 3 + 3 f d ( f a f d ) 2 ] + λ 4 A 3 64 [ 3 ( f a f d ) 4 + 4 f d ( f a f d ) 3 ] λ 5 A 4 160 [ 4 ( f a f d ) 5 + 5 f d ( f a f d ) 4 ] }
φ 2 ( f a ) is the second-order coupling phase which contains the range modulation phase and the second-order cross-coupling, expressed as:
φ 2 ( f a ) = π f a 2 f c 2 { λ A 1 2 λ 2 A 2 6 [ 3 ( f a f d ) ] + λ 3 A 3 16 [ 6 ( f a f d ) 2 ] λ 4 A 4 40 [ 10 ( f a f d ) 3 ] } π K r
φ 3 ( f a ) is the third-order coupling phase which is relevant on the third- order range frequency. It can be defined in:
φ 3 ( f a ) = π f a 2 f c 3 { λ A 1 2 λ 2 A 2 6 [ 4 ( f a f d ) + f d ] + λ 3 A 3 16 [ 10 ( f a f d ) 2 + 4 f d ( f a f d ) ] λ 4 A 4 40 [ 20 ( f a f d ) 2 + 10 f d ( f a f d ) ] }

4. Imaging Algorithm

In this section, we concentrate on the corresponding imaging algorithm. The proposed algorithm includes two parts. The first part, the azimuth preprocessing, is used to remove the aliasing in the staring observation mode; the second one is the improved CS algorithm.

4.1. Azimuth Preprocessing

In the staring observation mode, the azimuth Doppler bandwidth is composed of two parts. One part is the beam coverage Doppler bandwidth of the radar, known as the instantaneous Doppler bandwidth as strip-map mode. The other part is steering of the beam around the observation point causing rotation of the Doppler center. It introduces an extra bandwidth called beam rotation bandwidth, which may be several times larger than the instantaneous Doppler bandwidth. In space-borne SAR, to obtain a larger imaging swath, pulse repetition frequency is usually designed by instantaneous Doppler bandwidth that will cause the azimuth spectrum aliasing in the staring imaging mode. To remove the aliasing, two methods are usually adopted. One is the sub-aperture method [32,33], where the data is divided into small segments in the azimuth and imaging respectively processing. However, the sub-aperture operation is inefficient and complex. The other one is azimuth preprocessing based on two-step. The two step processing method is illustrated in [34,35,36], usually used to remove the aliasing in spotlight SAR and sliding-spotlight SAR. We employ it to remove azimuth aliasing in staring observation. Subsequently, the signal aliasing is removed, resulting in imaging processing without ambiguity.

4.2. Improved CS Algorithm

Due to the curved trajectory and high order slant mode, classical imaging algorithms cannot be directly used in the case of GEO SAR. Therefore, an improved CS algorithm based on the fifth-order slant range is adapted to imaging processing. After the azimuth preprocessing, the signal third-order coupling phase can be compensated in frequency domain. The reference function is described as:
H 1 ( f r , f a , R r e f ) = exp ( j φ 3 ( f a , R r e f ) f r 3 )
where R r e f is the slant range in the beam center.
Then, the range inverse fast Fourier transform (IFFT) of the signal is computed, and the expression can be obtained through the following expression:
S ( t r , f a ; R c ) = W r ( t r + φ 1 ( f a ) π ) W a ( f a f d ) exp ( j φ 0 ( f a ; R c ) j π 2 φ 2 ( f a ) ( t r t d ( f a , R c ) ) 2 )
where, t d ( f a , R c ) = φ 1 ( f a , R c ) / ( 2 π ) .
The next step is to determine the scaling factor of the improved CS algorithm. Due to the high-order slant mode, the linear scaling factor is not easy to get and difficult to describe with the indicated equation. However, we can get it under the numerical polynomial fitting model as described in [37].
According to the numerical fitting, the linear scaling factor C s ( f a ) can be obtained, and the line module frequency scaling function is described as:
H 2 ( t r , f a , R r e f ) = exp ( j π ( C s ( f a ) 1 ) π φ 2 ( f a , R r e f ) ( t r t d ( f a , R r e f ) ) 2 )
After multiplying by the linear scaling function and computing the ranged FFT, the expression of the echo signals can be expressed as follows:
S ( t r , f a ; R c ) = W r ( π f r φ 2 ( f a , R c ) C s ( f a ) ) W a ( f a f d ) exp ( j φ 2 ( f a , R c ) f r 2 C s ( f a ) ) exp ( j 4 π f r c ( R c + R r e f ( C ( f a ) 1 ) ) ) exp ( j φ 0 ( f a , R c ) ) exp ( j 4 π 2 c 2 C s ( f a ) ( C s ( f a ) 1 ) 1 φ 2 ( f a ; R c ) ( R c R r e f ) 2 )
Then, the phase of the range cell migration correction phase and the secondary range compression phase are compensated in the 2-D frequency domain. The function can be expressed as:
H 3 ( t r , f a , R r e f ) = exp ( j π f r 2 C s ( f a ) φ 2 ( f a , R r e f ) ) exp ( j 4 π ( C s ( f a ) 1 ) R r e f c )
After applying range IFFT, the azimuth compression and residual phase functions are presented by a fifth-order polynomial in the frequency domain, which is given by:
H 4 ( t r , f a , R r e f ) = exp ( j φ 0 ( f a ) ) exp ( j 4 π 2 c 2 ( C s ( f a ) 1 ) C s ( f a ) φ 2 ( f a , R r e f ) ( R c R r e f ) 2 )
Finally, the echo data are transformed into the 2-D time domain, and the image is obtained. The flowchart of the imaging algorithm is shown in Figure 4.

5. Simulation Results

In order to validate the algorithm, simulation experiments are performed. The imaging quality of the proposed algorithm is evaluated. The parameters used in our experiments are described in Table 1.
The satellite orbit time is selected at perigee passage time 0 and the designed imaging scene with size of 80 km × 40 km is shown in Figure 5. The integration time is 1800 s for the targets in the scene.
The target T2’s scene bandwidth including the extra bandwidth introduced by the steering of the antenna and the instantaneous bandwidth of the beam is calculated, which is about 260 Hz. The instantaneous bandwidth is varied with the perigee passage time, as shown in Figure 6. It can be seen that the instantaneous bandwidth is about 20 Hz at perigee passage time 0 h, far less than scene bandwidth.
The pulse repetition frequency is selected to 90 Hz. Figure 7a shows that the aliasing is occurs in the echo spectrum. However, when the operation of the azimuth preprocessing is used, the aliasing in the azimuth is removed, as shown in Figure 7b.

5.1. Simulation Results of the Conventional Algorithm

To compare the imaging results of algorithm, the conventional algorithm based on the fourth-order slant range mode is used [19]. The imaging results are shown in Figure 8. It can be clearly seen that all the point targets are severely defocused in the scene.

5.2. Simulation Results of the Improved Algorithm

The improved scaling imaging algorithm based on the fifth-order slant model is used for imaging processing. The imaging results are shown in Figure 9.
The imaging performance including resolution (RZ), peak sidelobe rate (PSLR) and integrated sidelobe ratio (ISLR) of point targets are evaluated without windowing treatment as shown in Table 2. In order to compare the imaging results obtained by proposed algorithm and algorithm in [19], simulation results of the algorithm in [19] are also shown in Table 2.
Comparing the analytical results utilizing the improved algorithm of T1, T2 and T3, the PSLR in azimuth slightly decreases along the azimuth direction, which is mainly caused by the azimuth-variant compressed function. Similarly, the PLSR and ISLR in range are also slowly falling along the range direction which is mainly caused by the differential RCMC in range for the degradation, whereas from the table can be seen clearly that imaging performance close to the theoretical value which can meet the demand of imaging. However, azimuth and range defocusing occurs by using the algorithm [19], therefore, the azimuth resolution and ISLR has seriously deteriorated. By comparison, the proposed algorithm has better performance.
In order to further demonstrate the advantage of the proposed algorithm, different orbital positions such as at perigee passage time 2 h are randomly selected for imaging simulation. In the case of satellite yaw angle control, the imaging results obtained are shown in Figure 10.
The imaging performance at perigee passage time 2 h are evaluated without windowing treatment as shown in Table 3. As shown in Table 3, the resolution is slightly different from Table 2, which is caused by the orbital characteristics. The ideal PSLR should be −13.26 dB. The maximal loss of range and azimuth PSLRs is less than 0.23 dB, indicating good focusing quality. The range and azimuth ISLRs over the whole swath are around the standard value, also indicating good focusing quality.
In order to further verify the improved algorithm, the real image is used (Huashan Mountain, Shaanxi, China, obtained from sliding rail vehicle experiment) as the input radar cross section information in echo generation. The simulation parameters are the same as listed in Table 1. The synthetic aperture time is 1800 s. Both the conventional algorithm and the improved algorithm are used for focusing the SAR raw data. The imaging results of Huashan Mountain area are depicted.
As can be seen from Figure 11, it is obvious that the well-focused results are obtained by the improved scaling imaging algorithm based on the fifth-order slant model and as shown in Figure 11b, but the conventional algorithm results are defocusing as shown in Figure 11a. All the point target simulation results and real data simulation results prove the effectiveness of the improved algorithm.

6. Conclusions

GEO SAR has recently received increasing attention because of its great potential for future spaceborne SAR missions. An imaging algorithm for GEO SAR with staring observation is proposed. Based on the fifth-order slant range mode, the 2-D spectrum is deduced with consideration of the impact of “stop-and-go” assumption, and then a two step processing method is used to remove the signal spectrum aliasing. An improved scaling algorithm is presented for imaging processing. Finally, based on the orbit parameters and radar parameters, simulation experiments have demonstrated that the algorithm is effective for an L-band GEO SAR with 1800 s integration time. Besides, when it comes to longer synthetic aperture and wider swaths, the effects of spatial variation, atmospheric refractive index and orbit perturbation on image formation and the higher order slant model should be studied in depth. All of these will be major topics in our future work.

Acknowledgments

Financial support for this study was provided by the National Defense Key Laboratory Foundation of China under Grants 9140C5305020902 and National Natural Science Foundation of China (project number 61420106007).

Author Contributions

All authors have made great contributions to the work. Caipin Li and Mingyi He carried out the theoretical framework, conceived and designed the experiments; Caipin Li performed the experiments and wrote the manuscript; Mingyi He gave insightful suggestions for the work and revised the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Martino, G.D.; Iodice, A.; Riccio, D.; Ruello, G. Disaster Monitoring by Extracting Geophysical Parameters from SAR Data. In Proceedings of the IGARSS, Barcelona, Spain, 23–28 July 2007; pp. 4949–4952. [Google Scholar]
  2. Heygster, G.; Dannenberg, J.; Notholt, J. Topographic Mapping of the German Tidal Flats Analyzing SAR Images with the Waterline Method. IEEE Trans. Geosci. Remote Sens. 2010, 48, 1019–1030. [Google Scholar] [CrossRef]
  3. Srivastava, H.S.; Patel, P.; Sharma, Y.; Navalgund, R.R. Large-area soil moisture estimation using multi-incidence-angle RADARSAT-1 SAR data. IEEE Trans. Geosci. Remote Sens. 2009, 47, 2528–2535. [Google Scholar] [CrossRef]
  4. Brekke, C.; Solberg, A.H.S. Oil spill detection by satellite remote sensing. Remote Sens. Environ. 2005, 95, 1–13. [Google Scholar] [CrossRef]
  5. Kurosu, T.; Fujita, M.; Chiba, K. Monitoring of rice crop growth from space using the ERS-1 C-band SAR. IEEE Trans. Geosci. Remote Sens. 1995, 33, 1092–1096. [Google Scholar] [CrossRef]
  6. Le Toan, T.; Quegan, S.; Davidson, M.W.J.; Balzter, H.; Paillou, P.; Papathanassiou, K. The BIOMASS mission: Mapping global forest biomass to better understand the terrestrial carbon cycle. Remote Sens. Environ. 2011, 115, 2850–2860. [Google Scholar] [CrossRef]
  7. Henderson, F.M.; Lewis, A.J. Radar Detection of Wetland Ecosystems: A Review. Int. J. Remote Sens. 2008, 29, 5809–5835. [Google Scholar] [CrossRef]
  8. Sadjadi, F. New Comparative Experiments in Range Migration Mitigation Methods Using Polarimetric Inverse Synthetic Aperture Radar Signatures of Small Boats. In Proceedings of the IEEE Radar Conference, Cincinnati, OH, USA, 19–23 May 2014; pp. 0613–0616. [Google Scholar]
  9. Sadjadi, F. Radar beam sharpening using an optimum FIR filter. Circuits Syst. Signal Process. 2000, 19, 121–129. [Google Scholar] [CrossRef]
  10. Sadjadi, F.A. New Experiments in Inverse Synthetic Aperture Radar Image Exploitation for Maritime Surveillance. In Proceedings of the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Baltimore, MD, USA, 13 June 2014; p. 09. [Google Scholar]
  11. Hobbs, S.; Mitchell, C.; Forte, B.; Holley, R. System design for geosynchronous synthetic aperture radar missions. IEEE Trans. Geosci. Remote Sens. 2014, 52, 7750–7763. [Google Scholar] [CrossRef]
  12. Tomiyasu, K.; Pacelli, J.L. Synthetic aperture radar imaging from an inclined geosynchronous orbit. IEEE Trans. Geosci. Remote Sens. 1983, 21, 324–329. [Google Scholar] [CrossRef]
  13. Bruno, D.; Hobbs, S.E.; Ottavianelli, G. Geosynchronous synthetic aperture radar: Concept design, properties and possible applications. Acta Astron. 2006, 59, 149–156. [Google Scholar] [CrossRef]
  14. Bruno, D.; Hobbs, S.E. Radar imaging from geosynchronous orbit: Temporal decorrelation aspects. IEEE Trans. Geosci. Remote Sens. 2010, 48, 2924–2929. [Google Scholar] [CrossRef]
  15. Sheng, W.; Hobbs, S.E. Research on compensation of motion, Earthcurvature and tropospheric delay in GEO SAR. Acta Astron. 2011, 68, 2005–2011. [Google Scholar] [CrossRef]
  16. Ruiz-Rodon, J.; Broquetas, A.; Makhoul, E.; Monti Guarnieri, A. Nearly zero inclination geosynchronous SAR mission analysis with long integration time for earth observation. IEEE Trans. Geosci. Remote Sens. 2014, 52, 6379–6391. [Google Scholar] [CrossRef]
  17. Mittermayer, J.; Wollstadt, S.; Prats-Iraola, P.; Scheiber, R. The TerraSAR-X Staring Spotlight Mode Concept. IEEE Trans. Geosci. Remote Sens. 2014, 52, 3695–3706. [Google Scholar] [CrossRef]
  18. Li, C.P.; He, M.Y. A novel attitude steering strategy for GEO SAR staring imaging. In Proceedings of the IEEE China SIP, Chengdu, China, 12–15 July 2015; pp. 630–634. [Google Scholar]
  19. Bao, M.; Xing, M. Chirp scaling algorithm for GEO SAR based on fourth-order range equation. Electron. Lett. 2012, 48, 56–57. [Google Scholar] [CrossRef]
  20. Hu, C.; Liu, Z.P.; Long, T. An improved CS algorithm based on the curved trajectory in geosynchronous SAR. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2012, 5, 795–808. [Google Scholar] [CrossRef]
  21. Liu, Z.P.; Hu, C.; Zeng, T. Improved secondary range compression focusing method in GEO SAR’. In Proceedings of the IEEE ICASSP, Prague, Czech Republic, 22–27 May 2011; pp. 1373–1376. [Google Scholar]
  22. Hu, C.; Long, T.; Liu, Z.P.; Zeng, T. An improved frequency domain focusing method in geosynchronous SAR. IEEE Trans. Geosci. Remote Sens. 2014, 52, 5514–5528. [Google Scholar]
  23. Hu, B.; Jiang, Y.C.; Zhang, S.; Zhang, Y. Generalized Omega-K Algorithm for Geosynchronous SAR Image Formation. IEEE Trans. Geosci. Remote Sens. Lett. 2015, 12, 2286–2290. [Google Scholar] [CrossRef]
  24. Sun, G.C.; Xing, M.D.; Wang, Y.; Yang, J. A 2-D space-variant chirp scaling algorithm based on the RCM equalization and subband synthesis to process geosynchronous SAR data. IEEE Trans. Geosci. Remote Sens. 2014, 52, 4868–4880. [Google Scholar]
  25. Zhao, B.J.; Han, Y.Z.; Gao, W.J.; Luo, Y.; Han, X. A New Imaging Algorithm for Geosynchronous SAR Based on the Fifth-Order Doppler Parameters. Prog. Electromagn. Res. B 2013, 55, 195–215. [Google Scholar] [CrossRef]
  26. Yu, Z.; Lin, P.; Xiao, P.; Kang, L.; Li, C. Correcting Spatial Variance of RCM for GEO SAR Imaging Based on Time-Frequency Scaling. Sensors 2016, 16, 1091. [Google Scholar] [CrossRef] [PubMed]
  27. Li, Z.; Li, C.S.; Yu, Z.; Zhou, J.; Chen, J. Back projection algorithm for high resolution GEO-SAR image formation. In Proceedings of the IGARSS, Vancouver, BC, Canada, 24–29 July 2011; pp. 336–339. [Google Scholar]
  28. Ruiz Rodon, J.; Broquetas, A.; Monti Guarnieri, A.; Rocca, F. Geosynchronous SAR focusing with atmospheric phase screen retrieval and compensation. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4397–4404. [Google Scholar] [CrossRef]
  29. Dong, X.; Hu, C.; Long, T.; Li, Y. Numerical Analysis of Orbital Perturbation Effects on Inclined Geosynchronous SAR. Sensors 2016, 16, 1420. [Google Scholar] [CrossRef] [PubMed]
  30. Hu, C.; Tian, Y.; Zeng, T.; Long, T. Adaptive Secondary Range Compression Algorithm in Geosynchronous SAR. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 9, 1397–1413. [Google Scholar] [CrossRef]
  31. Jiang, M.; Hu, W.; Ding, C.; Liu, G. The effects of orbital perturbation on geosynchronous synthetic aperture radar imaging. IEEE Trans. Geosci. Remote Sens. 2015, 12, 1106–1110. [Google Scholar] [CrossRef]
  32. Mittermayer, J.; Moreira, A. Spotlight SAR processing using the extended chirp scaling algorithm. In Proceedings of the IEEE IGARSS, Singapore, 3–8 August 1997; pp. 2021–2023. [Google Scholar]
  33. Prats, P.; Scheiber, R.; Mittermayer, J.; Meta, A. Processing of sliding spotlight and TOPS SAR data using baseband azimuth scaling. IEEE Trans. Geosci. Remote Sens. 2010, 48, 70–780. [Google Scholar] [CrossRef]
  34. Lanari, R.; Tesauro, M.; Sansosti, E.; Fornaro, G. Spotlight SAR data focusing based on a two-step processing approach. IEEE Trans. Geosci. Remote Sens. 2011, 39, 1993–2004. [Google Scholar] [CrossRef]
  35. Liu, F.F.; Ding, Z.G.; Zeng, T.; Long, T. Performance analysis of two-step algorithm in sliding spotlights space-borne SAR. In Proceedings of the Radar Conference, Washington, DC, USA, 10–14 May 2011; pp. 965–968. [Google Scholar]
  36. Sun, G.C.; Xing, M.D.; Wang, Y.; Wu, Y. Sliding spotlight and TOPS SAR data processing without subaperture. IEEE Trans. Geosci. Remote Sens. Lett. 2011, 8, 1036–1040. [Google Scholar] [CrossRef]
  37. Huang, L.J.; Qiu, X.L.; Hu, D.H.; Ding, C.B. Focusing of medium-earth-orbit SAR with advanced nonlinear chirp scaling algorithm. IEEE Trans. Geosci. Remote Sens. 2011, 49, 500–508. [Google Scholar] [CrossRef]
Figure 1. The geometry of satellite and ground targets.
Figure 1. The geometry of satellite and ground targets.
Sensors 17 01058 g001
Figure 2. Phase errors caused by slant range mode. (a) Phase error of the fourth-order slant range model as a function of time and (b) Phase error of the fifth-order slant range model as a function of time.
Figure 2. Phase errors caused by slant range mode. (a) Phase error of the fourth-order slant range model as a function of time and (b) Phase error of the fifth-order slant range model as a function of time.
Sensors 17 01058 g002
Figure 3. Relationship between different order slant range model and synthetic aperture time.
Figure 3. Relationship between different order slant range model and synthetic aperture time.
Sensors 17 01058 g003
Figure 4. Flowchart of the imaging algorithm.
Figure 4. Flowchart of the imaging algorithm.
Sensors 17 01058 g004
Figure 5. Target positions in the scene.
Figure 5. Target positions in the scene.
Sensors 17 01058 g005
Figure 6. The instantaneous bandwidth.
Figure 6. The instantaneous bandwidth.
Sensors 17 01058 g006
Figure 7. Signal spectrum without and with azimuth pre-processing. (a) signal spectrum aliasing without azimuth pre-processing;(b) signal spectrum aliasing removed with azimuth pre-processing.
Figure 7. Signal spectrum without and with azimuth pre-processing. (a) signal spectrum aliasing without azimuth pre-processing;(b) signal spectrum aliasing removed with azimuth pre-processing.
Sensors 17 01058 g007
Figure 8. Imaging results by using conventional fourth-order algorithm (a) Target T1; (b) Target T2; (c) Target T3; (d) Range profile of point target T1; (e) Range profile of point target T2; (f) Range profile of point target T3; (g) Azimuth profile of point target T1; (h) Azimuth profile of point target T2; (i) Azimuth profile of point target T3.
Figure 8. Imaging results by using conventional fourth-order algorithm (a) Target T1; (b) Target T2; (c) Target T3; (d) Range profile of point target T1; (e) Range profile of point target T2; (f) Range profile of point target T3; (g) Azimuth profile of point target T1; (h) Azimuth profile of point target T2; (i) Azimuth profile of point target T3.
Sensors 17 01058 g008aSensors 17 01058 g008b
Figure 9. Imaging results at perigee passage time 0 h. (a) 2-D contour of point target T1; (b) 2-D contour of point target T2; (c) 2-D contour of point target T3; (d) Range profile of point target T1; (e) Range profile of point target T2; (f) Range profile of point target T3; (g) Azimuth profile of point target T1; (h) Azimuth profile of point target T2; (i) Azimuth profile of point target T3.
Figure 9. Imaging results at perigee passage time 0 h. (a) 2-D contour of point target T1; (b) 2-D contour of point target T2; (c) 2-D contour of point target T3; (d) Range profile of point target T1; (e) Range profile of point target T2; (f) Range profile of point target T3; (g) Azimuth profile of point target T1; (h) Azimuth profile of point target T2; (i) Azimuth profile of point target T3.
Sensors 17 01058 g009aSensors 17 01058 g009b
Figure 10. Imaging results at perigee passage time 2 h. (a) 2-D contour of point target T1; (b) 2-D contour of point target T2; (c) 2-D contour of point target T3; (d) Range profile of point target T1; (e) Range profile of point target T2; (f) Range profile of point target T3; (g) Azimuth profile of point target T1; (h) Azimuth profile of point target T2; (i) Azimuth profile of point target T3.
Figure 10. Imaging results at perigee passage time 2 h. (a) 2-D contour of point target T1; (b) 2-D contour of point target T2; (c) 2-D contour of point target T3; (d) Range profile of point target T1; (e) Range profile of point target T2; (f) Range profile of point target T3; (g) Azimuth profile of point target T1; (h) Azimuth profile of point target T2; (i) Azimuth profile of point target T3.
Sensors 17 01058 g010aSensors 17 01058 g010b
Figure 11. Comparison of imaging results of Huashan Mountain area. (a) conventional algorithm; (b) improved algorithm based on the fifth-order slant model.
Figure 11. Comparison of imaging results of Huashan Mountain area. (a) conventional algorithm; (b) improved algorithm based on the fifth-order slant model.
Sensors 17 01058 g011
Table 1. Simulation Parameters.
Table 1. Simulation Parameters.
ParametersValue
Orbit semi-major axis42,164 km
Orbit eccentricity0
Orbit inclination20°
Argument of perigee95°
Right ascension of ascending node97°
Radar frequency1.25 GHz
Pulse repetition frequency90 Hz
Bandwidth80 MHz
Antenna beam width0.5°
Off-nadir angle3.5°
Table 2. Imaging results of point targets at perigee passage time 0 h.
Table 2. Imaging results of point targets at perigee passage time 0 h.
Target PositionsProposed AlgorithmAlgorithm in [19]
RangeAzimuthRangeAzimuth
RZ (m)PSLR (dB)ISLR (dB)RZ (m)PSLR (dB)ISLR (dB)RZ (m)PSLR (dB)ISLR (dB)RZ (m)PSLR (dB)ISLR (dB)
T13.80−13.22−9.832.11−13.02−9.354.75−2.357−7.3546−4.21−2.52
T23.86−13.27−9.752.02−13.18−9.683.81−2.8658.8632−2.361.25
T33.92−13.20−9.912.13−13.05−9.574.63−1.08−1.0848−5.01−1.63
Table 3. Imaging results of point targets at perigee passage time 2 h.
Table 3. Imaging results of point targets at perigee passage time 2 h.
Target PositionsRangeAzimuth
RZ (m)PSLR (dB)ISLR (dB)RZ (m)PSLR (dB)ISLR (dB)
T13.82−13.08−9.522.35−13.20−9.45
T23.93−13.19−9.692.27−13.16−9.71
T34.01−13.06−9.652.47−13.03−9.35

Share and Cite

MDPI and ACS Style

Li, C.; He, M. A Generalized Chirp-Scaling Algorithm for Geosynchronous Orbit SAR Staring Observations. Sensors 2017, 17, 1058. https://doi.org/10.3390/s17051058

AMA Style

Li C, He M. A Generalized Chirp-Scaling Algorithm for Geosynchronous Orbit SAR Staring Observations. Sensors. 2017; 17(5):1058. https://doi.org/10.3390/s17051058

Chicago/Turabian Style

Li, Caipin, and Mingyi He. 2017. "A Generalized Chirp-Scaling Algorithm for Geosynchronous Orbit SAR Staring Observations" Sensors 17, no. 5: 1058. https://doi.org/10.3390/s17051058

APA Style

Li, C., & He, M. (2017). A Generalized Chirp-Scaling Algorithm for Geosynchronous Orbit SAR Staring Observations. Sensors, 17(5), 1058. https://doi.org/10.3390/s17051058

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop