Next Article in Journal
Regional Vicarious Calibration of the SWIR-Based Atmospheric Correction Approach for MODIS-Aqua Measurements of Highly Turbid Inland Water
Previous Article in Journal
Fault Slip Model of the 2018 Mw 6.6 Hokkaido Eastern Iburi, Japan, Earthquake Estimated from Satellite Radar and GPS Measurements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Fast High-Resolution Imaging Algorithm for Helicopter-Borne Rotating Array SAR Based on 2-D Chirp-Z Transform

1
School of Microelectronics and Communication Engineering, Chongqing University, Chongqing 400044, China
2
Center of Communication and Tracking Telemetering Command, Chongqing University, Chongqing 400044, China
3
National Laboratory of Radar Signal Processing, Xidian University, Xian 710071, China
4
Shanghai Satellite Engineering Research Institute, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(14), 1669; https://doi.org/10.3390/rs11141669
Submission received: 18 April 2019 / Revised: 4 July 2019 / Accepted: 11 July 2019 / Published: 13 July 2019
(This article belongs to the Section Engineering Remote Sensing)

Abstract

:
In view of the azimuth resolution of the helicopter-borne rotating array synthetic aperture radar (RoASAR) depending on the azimuth reconstruction angle and the sector distortion caused by the azimuth Deramp processing, this paper proposes an efficient helicopter-borne RoASAR high-resolution imaging algorithm based on two-dimensional (2-D) chirp-z transform (CZT). First, the high-order Taylor series expansion is performed on the slant range, and the accurate 2-D spectral expression of the point target is obtained by using the method of series reversion (MSR). Based on that, the space-variant characteristics of the range cell migration (RCM) terms are analyzed. After that, the space-variant RCM and the sector distortion effect caused by the azimuth Deramp processing are removed by using efficient 2-D CZT, thereby increasing the azimuth reconstruction angle and improving the azimuth resolution. The proposed algorithm is efficient without the interpolation operation, and it is easy to implement in real-time. Finally, the simulations are provided to verify the effectiveness of the proposed algorithm.

1. Introduction

Helicopter-borne rotating array synthetic aperture radar (RoASAR) is a novel radar imaging model and system, which is installed and operated on a helicopter platform, where the radar antenna is set on a rigid support at the tip of the helicopter rotor. With the rotation of the rotor, the 360° area synthetic aperture radar (SAR) image around the helicopter platform can be obtained without the movement of the helicopter platform [1,2]. This advantage breaks through the shortcomings that traditional SAR imaging is restricted by the blind area at the front of the platform, and it has become a popular choice for helicopter-borne anti-collision systems.
In a complex low-altitude environment, the RoASAR system provides high-resolution SAR images in front of the helicopter platform, for full day, all-weather and real-time monitoring of potential threats [3,4,5]. In addition to that, the RoASAR system can perform continuous data acquisition for fixed targets. In recent years, some imaging methods for RoASAR have been proposed to obtain the SAR images. The time-domain correlation imaging method, such as the back-projection (BP) algorithm in [6], can be applied to any geometric configuration, but it is not suitable for real-time imaging processing due to its high computational complexity. In recent years, some improved fast time-domain algorithms have also been proposed to reduce the computational complexity [7,8,9]. The sub-aperture imaging method proposed in [10] is of a low computational complexity but suffers from low azimuth resolution. To overcome the low-resolution shortcoming, the Omega-k imaging method is proposed in [11,12,13], but it also has the great computation burden caused by the Stolt interpolation operation, and it is only applicable to the small imaging area near the center of the scene due to the range of spatial-variant characteristics of the tangential velocity. To reduce the computational complexity, a novel imaging method based on the chirp-z transform (CZT) is proposed in [14,15], in which an efficient CZT is utilized to remove the range-azimuth coupling. However, this method ignores the azimuth sector distortion caused by azimuth Deramp operation. With the requirements of a real-time, high-resolution and imaging scene scope requirements of the helicopter system, a fast high-resolution RoASAR imaging algorithm still needs to be studied.
To overcome the aforementioned issues, an efficient high resolution RoASAR imaging algorithm based on the two-dimensional (2-D) CZT is proposed in this work. First, for the imaging geometric configuration, the system parameters design and azimuth resolution dependence characteristics on the beam width and rotor rigid length are analyzed. With the high-order Taylor expansion of the instantaneous slant range, the precise 2-D frequency analytic expression is derived by method of series reversion (MSR) [16]. Furthermore, according to the analysis on the range-dependent range cell migration (RCM), an efficient inverse CZT (ICZT) is developed to remove the spatial-variant features of RCM. SAR image distortion in azimuth direction caused by the azimuth Deramp operation is analyzed and compensated by the improved CZT, and a well-focused high-resolution SAR image is obtained. Compared with the existing methods, the proposed method is computationally efficient thanks to avoiding the interpolation operation and is of high azimuth resolution because the large azimuth reconstruction angle is utilized. Finally, simulation results and analysis verify the correctness and effectiveness of the proposed algorithm.
The rest of the paper is organized as follows. In Section 2, the geometric and signal model of airborne RoASAR system are introduced, and then, the system parameters design, azimuth angular resolution and slant range approximation are analyzed. Section 3 describes the proposed 2-D CZT imaging algorithm. In Section 4, simulation results and analysis are provided to verify the feasibility of the algorithm. Section 5 concludes the paper briefly.

2. Overall Analysis of the Helicopter-Borne RoASAR System

2.1. Geometric Model of the Helicopter-Borne RoASAR System

Figure 1a,b show the basic imaging geometric model of the RoASAR system, in which the antenna is integrated at the end of the inelastic rotary arm and mounted on the helicopter rotor shaft. In this configuration, the helicopter platform hovers at a height of H above the ground, the length of the rotating arm is r a , and the angular velocity is ω . The irradiation direction of the radar is always perpendicular to the motion direction of the rotor end. Since the hovering state of the helicopter is regarded as a fixed platform, the lower angle θ is constant during the operation of the SAR.
Let t a represent the azimuth slow time variable and assume that the beam is illuminated on the positive half of the X-axis at time t a = 0 . The coordinate of the reference point scatterer P in the scene can be expressed as ( r p c o s θ p ,   r p s i n θ p ,   0 ) , where r p is the ground distance from the phase center of the antenna to the ground scatter P and ( r a c o s θ t ,   r a s i n θ t ,   H ) is the instantaneous coordinate of the radar antenna phase center at the time t a . The relationship between the angle and radius is
{ θ t = θ 0 + ω t a r p = r a + H t a n θ θ p = ω t p .
In the helicopter-borne RoASAR system, the positive half of the X-axis is selected as t a = 0, and therefore, the corresponding initial angle is θ 0 = 0 . The instantaneous slant range from the radar antenna to the scatter P is expressed as
R ( t a ; r p ) = ( r a c o s θ t r p c o s θ p ) 2 + ( r a s i n θ t r p s i n θ p ) 2 + H 2 = r a 2 + r p 2 + H 2 2 r a r p c o s ( ω t a ω t p ) .
Equation (2) is an exact expression of the slant range of RoASAR. Because of the cosine term in the slant range, it is difficult to obtain the accurate 2-D spectrum of the point target by using the principle of stationary phase (POSP), which brings difficulties to the design of an efficient frequency domain imaging algorithm. Therefore, some approximations must be made and the specific approximation process is given after the resolution and parameter analysis.

2.2. Resolution Analysis of the Helicopter-Borne RoASAR System

The analysis of RoASAR resolution and related system parameters selection are discussed now. The range resolution ρ r is determined by the bandwidth B r of the transmitted signal, i.e., ρ r = c / ( 2 B r ) , and the corresponding ground resolution is Δ R g = c / ( 2 B r s i n θ i n ) , where θ i n is the incident angle. While the azimuth resolution of the helicopter-borne RoASAR is different from the traditional stripmap SAR, it depends on the size of the azimuth reconstruction angle, the helicopter rotor length and the wavelength. According to the matched filtering theory, the azimuth temporal resolution depends on the azimuth Doppler bandwidth B a , i.e., Δ t a = 1 / B a . In RoASAR, bandwidth B a is
B a = 2 · 2 v λ c o s ( φ m a x ) ,
where φ m a x is the maximal angle between the antenna and the center of the aperture, and v is the tangential velocity of antenna corresponding to the uniform rotation ω. From the diagram of the beam radiation trajectory for the target P in Figure 2, it is calculated that φ m a x = c o s 1 ( r p α / 2 R a m a x ) , where α is the azimuth reconstruction angle and R a m a x is the maximum slant range corresponding to the point target in the reconstruction angle range. Therefore, the corresponding azimuth angular resolution is
Δ α = ω · Δ t = λ 4 r a c o s ( φ m a x ) = λ 4 r a 2 H 2 + r p 2 + r a 2 2 r a r p c o s ( α / 2 ) r p · α .
It can be seen that the azimuth angle resolution is related not only to the ground distance of the target r p , but also to the rotor length r a , helicopter platform height H and azimuth reconstruction angle α . In Figure 3a, the result of the angular resolution at the reconstructed angle α = 70 ° and target ground distance r p = 2000   m against the lengths of the rotor and the heights of the aircraft, is depicted. Figure 3b shows the azimuth angular resolution versus the reconstructed angle when the length of the rotor, the height of the aircraft and the target ground distance are fixed. Figure 3c shows the azimuth angular resolution versus the observation angle when the length of the rotor, the height of the aircraft and the target ground distance are fixed.
From Figure 3a, the azimuth angular resolution decreases with the increase in the height of the helicopter platform H . This is mainly because the increase in the height of the aircraft causes the component of the antenna linear velocity on the radar and target line of sight (LOS) to become smaller, decreasing the Doppler bandwidth. The longer the rotor length is, the higher the angular resolution is. The reason is that with the azimuth angle being fixed, the longer the rotor length is, the larger the synthetic aperture will be. Figure 3b shows that while other parameters are determined, the angular resolution of the target increases with an increase in reconstruction angle, which is consistent with the theoretical analysis.
Due to the physical limits of the helicopter’s structure, the rotor length cannot be too long. To obtain higher azimuth resolution capability, it is hoped that the azimuth reconstructed angle can be as large as possible according to Figure 3b. However, the instantaneous slant range expression contains the cosine term, shown in Equation (2), which prevents the efficient frequency domain imaging algorithm developed. Therefore, some approximations need to be developed. For this problem, it is well-known that when the phase error introduced by the slant range approximation error Δ R is less than π / 2 , the influence on the imaging result can be ignored, i.e., 4 π Δ R / λ < π / 2 . According to the aforementioned analysis, we set the system parameters of helicopter-borne RoASAR as the platform height H = 1000   m , rotor length r a = 2   m , rotation angular velocity ω = 2   π f = 15   r a d / s , radar wavelength λ = 0.03 and the point target coordinates ( 2000 , 0 , 0 ) .
The detailed equations of second-order and fourth-order Taylor series expansion of the instantaneous slant range are developed in Appendix A (A9) and (A10), respectively. Under this case, Figure 4a,b show the approximated errors introduced by the second-order and fourth-order Taylor series expansion. It can be seen from Figure 4 that, in order to produce a well-focused image, the azimuth reconstruction angle can only be taken in the range between 27 ° and + 27 ° , while the range of the reconstruction angle can be reached between 55 ° and + 55 ° , when the slant range is expanded as the fourth Taylor series. Therefore, in this work, slant range is expanded as a fourth-order Taylor series.
The slant range in Equation (2) is expanded as a fourth-order Taylor series, and the higher order terms are ignored, and two-way slant range now is
2 R ( t a ; r p ) 2 r a 2 + r p 2 + H 2 2 r a r p ( 1 ω 2 ( t a t p ) 2 2 + ω 4 ( t a t p ) 4 24 ) = 2 R p 2 + r a r p ω 2 ( t a t p ) 2 r a r p ω 4 12 ( t a t p ) 4 2 R p + r a r p ω 2 R p ( t a t p ) 2 ( r a r p ω 4 12 R p + r a 2 r p 2 ω 4 4 R p 3 ) ( t a t p ) 4 = R 0 + k 2 ( t a t p ) 2 + k 4 ( t a t p ) 4 ,
where R 0 , R p , k 2 and k 4 are expressed as
{ R 0 = 2 R p = 2 ( r p r a ) 2 + H 2 k 2 = r a r p ω 2 R p k 4 = ( r a r p ω 4 12 R p + r a 2 r p 2 ω 4 4 R p 3 ) ,
where the coefficients k 2 and k 4 denote the range curvature and higher order phase coupling term, and R 0 is the minimum slant range between the target and the antenna phase center.
Assume that the radar transmits a linear frequency modulated (LFM) signal, then the received baseband echo signal of the point target after demodulation is
S 1 ( t r ; t a ) = w r [ t r 2 R ( t a ; r p ) c ] w a ( t a ) e x p [ j π K r ( t r 2 R ( t a ; r p ) c ) 2 ] × e x p [ j 4 π λ R ( t a ; r p ) ]   ,
where w r ( · ) and w a ( · ) are range and azimuth envelopes, respectively, t r denotes the range fast time variable, λ is the wavelength, K r is the frequency modulation factor, and c is the speed of light.

3. Fast Imaging Algorithm for Helicopter-Borne RoASAR Based on 2-D CZT

In Section 2, the RoASAR geometry and signal model were presented, and some system parameters selection criteria were discussed. Based on that, the detailed derivations of the proposed fast RoASAR imaging algorithm using the 2-D CZT is now provided, which mainly contains the following three steps: (1) the precise 2-D frequency analytic expression, (2) spatial-variant RCM correction using an efficient, improved ICZT, and (3) SAR image azimuth distortion caused by the Deramp operation compensated by the CZT.

3.1. 2-D Spectrum Derivation Based on MSR

Due to the quartic item of ( t a t p ) in Equation (5), it is difficult to obtain the exact solution of the stationary phase point when the POSP is used to obtain the 2-D spectrum. In this paper, an MSR algorithm is introduced to derive the accurate 2-D spectrum analytic expression. After substituting Equations (5) and (6) into Equation (7) and using MSR (the detailed process of MSR is developed in Appendix B), the obtained accurate 2-D spectrum analytic expression is
S 2 ( f r ; f a ) = W r ( f r ) W a ( f a ) · e x p ( 2 π f a t p ) · e x p [ j Φ ( f r ; f a ) ] ,
where the phase term Φ ( f r ; f a ) is expressed as
Φ ( f r ; f a ) = 2 π ( f r + f c ) c R 0 π f r 2 K r + π c 2 k 2 ( f r + f c ) f a 2 π c 3 k 4 2 k 2 4 ( f r + f c ) 3 f a 4 ,
where f r is the range frequency, f a is the azimuth frequency and f c is the carrier frequency. According to Equation (7), the range frequency and the azimuth frequency domain are coupled together, and expanding the phase term in Equation (9) as a power series of f r yields
Φ ( f r ; f a ) = ϕ 0 ( f a ; R 0 ) + ϕ 1 ( f a ; R 0 ) f r + ϕ 2 ( f a ; R 0 ) f r 2 + ϕ 3 ( f a ; R 0 ) f r 3 ,
where
{ ϕ 0 ( f a ; R 0 ) = 2 π f c R 0 c + π c 2 k 2 f c f a 2 π c 3 k 4 8 k 2 4 f c 3 f a 4 ϕ 1 ( f a ; R 0 ) = 2 π R 0 c π c 2 k 2 f c 2 f a 2 + 3 π c 3 k 4 8 k 2 4 f c 4 f a 4 ϕ 2 ( f a ; R 0 ) = π K r + π c 2 k 2 f c 3 f a 2 3 π c 3 k 4 4 k 2 4 f c 5 f a 4 ϕ 3 ( f a ; R 0 ) = π c 2 k 2 f c 4 f a 2 + 5 π c 3 k 4 4 k 2 4 f c 6 f a 4 ,
where ϕ 0 ( f a ; R 0 ) is the range-independent azimuth modulation term, ϕ 1 ( f a ; R 0 ) is the RCM term, ϕ 2 ( f a ; R 0 ) is the second range compression (SRC) term, and ϕ 3 ( f a ; R 0 ) is the cubic cross-coupling term. Due to the range-dependent characteristic of the RCM term, it is impossible to perform a precise range compression for all of the targets located in different range cells. In order to obtain accurate range processing, an efficient improved ICZT algorithm is developed.

3.2. Range Processing

Because the image quality is seriously affected by the inaccurate range processing, it is necessary to analyze the characteristics of the range and azimuth coupling terms. In the 2-D frequency domain, it is inconvenient to analyze three variables at the same time. In order to better illustrate the range-dependent feature of SRC and cubic coupling terms, f a and f r in Equation (11) are replaced by B a / 2 and B r / 2 , respectively, and the approximated expression of which has only one variable R 0 that is derived as follows:
Δ ϕ 2 = | ϕ 2 ( f a ; R 0 ) f r 2 | ϕ 2 ( B a 2 ; R 0 ) · ( B r 2 ) 2 Δ ϕ 3 = | ϕ 3 ( f a ; R 0 ) f r 3 | ϕ 3 ( B a 2 ; R 0 ) · ( B r 2 ) 3
where B a is the Doppler bandwidth and B r is the transmitted signal bandwidth.
From Section 2, the RoASAR simulation parameters are listed in Table 1. Using Equation (12), the range-dependent characteristic curves of SRC and cubic terms are shown in Figure 5a,b. It is seen that the variation is far less than π / 4 and their range-dependent effect can be neglected in the range processing.
From Figure 5c, the RCM term varies greatly. The final imaging result will be affected if it is ignored. Therefore, in order to obtain accurate focusing result, an improved ICZT method is proposed to correct range spatial-variant RCM in this paper.
According to the aforementioned analysis, the SRC and cubic cross-coupling terms can be compensated with the reference distance, and one builds a compensate function as follows:
H 1 ( f r ; f a ; R S ) = e x p [ j ϕ 3 ( f a ; R S ) f r 3 j ϕ 2 ( f a ; R S ) f r 2 ] ,
where R S is the reference distance between the center of the scene and the phase center of the radar, and the 2-D signal expression after compensation is
S 3 , 1 ( f r ; f a ) = W r ( f r ) W a ( f a ) · e x p ( 2 π f a t p ) e x p [ j ϕ 0 ( f a ; R 0 ) + j ϕ 1 ( f a ; R 0 ) f r ] .
Since the RCM term cannot be ignored directly, the traditional Omega-k algorithm involves interpolation operations, but it increases the computational complexity. In this paper, an efficient improved ICZT is adopted to correct the RCM [17,18,19]. The RCM term ϕ 1 ( f a ; R 0 ) is decomposed into two parts: One related to the scene reference distance and another related to the variation distance. By introducing a range variation term Δ R = R 0 R s , the RCM term ϕ 1 ( f a ; R 0 ) can be rewritten as
ϕ 1 ( f a ; R 0 ) = R s ϕ 1 ( f a ; R s ) R 0 + Δ R ϕ 1 ( f a ; R 0 ) R 0 .
First, we compensate the bulk RCM term with the reference distance of the scene. The corresponding compensation function H b u l k ( f r ; f a ; R s ) is used to compensate the bulk RCM term. The expression of the compensation function is
H b u l k ( f r ; f a ; R s ) = e x p ( j R s ϕ 1 ( f a ; R s ) f r R 0 ) .
The compensated 2-D signal is expressed as
S 3 , 2 ( f r ; f a ) = W r ( f r ) W a ( f a ) · e x p ( 2 π f a t p ) e x p [ j ϕ 0 ( f a ; R 0 ) + j Δ R ϕ 1 ( f a ; R 0 ) f r R 0 ] .
According to 2 Δ R f r = Δ R 2 + f r 2 ( f r Δ R ) 2 , the first-order range term related to Δ R is decomposed as
e x p ( j Δ R ϕ 1 ( f a ; R 0 ) f r R 0 ) = e x p ( j ϕ 1 ( f a ; R 0 ) R 0 Δ R 2 2 ) e x p ( j ϕ 1 ( f a ; R 0 ) R 0 f r 2 2 ) × e x p ( j ϕ 1 ( f a ; R 0 ) R 0 ( f r Δ R ) 2 2 ) .
The echo signal obtained by range inverse Fourier transform along the range frequency f r is
S 3 ( f r ; f a ; R s ) = S 3 , 2 ( f r ; f a ) e x p [ j Δ R ϕ 1 ( f a ; R s ) f r R s ] d f r = e x p ( j ϕ 1 ( f a ; R s ) R s Δ R 2 2 ) × S 3 , 2 ( f r ; f a ) × e x p ( j ϕ 1 ( f a ; R s ) R s f r 2 2 ) e x p ( j ϕ 1 ( f a ; R s ) R s ( f r Δ R ) 2 2 ) d f r .
In Equation (19), the exponential term unrelated to f r is removed from the integral and R 0 is replaced by R s (because R s has little influence on the error), and it is simplified as
S 3 ( f r ; f a ; R S ) = e x p ( j ϕ 1 ( f a ; R s ) R S Δ R 2 2 ) × { [ S 3 , 2 ( f r ; f a ) e x p ( j ϕ 1 ( f a ; R s ) R s f r 2 2 ) ] e x p ( j ϕ 1 ( f a ; R s ) R s f r 2 2 ) } .
Obviously, Equation (20) is a process of ICZT after expansion. By use of twice phase multiplications and one convolution operation, the expression of the range Doppler domain of the signal is
S I C Z T ( t r ; f a ) = w r ( t r R 0 c ) W a ( f a ) · e x p ( 2 π f a t p ) e x p [ j ϕ 0 ( f a ; R 0 ) ] .
From Equation (21), we can see that the range migration factor is equalized to be in unity after the ICZT operation, which means that the range spatial RCM is completely removed. Then, the azimuth processing can be applied to compress the azimuth modulated term.

3.3. Azimuth Processing

The Doppler centroid frequency of the target varies with the azimuth position under the wide azimuth beam. Standard SPECAN algorithm cannot obtain good focusing images in azimuth direction through Deramp processing and fast Fourier transform (FFT) [20]. The output of direct FFT in a standard SPECAN algorithm after the Deramp processing is
S s p e c a n ( t r ; f a ) = w r ( t r R 0 c ) w a ( π F a N a ( f a 2 r a ( r p r a ) ω 2 λ R 0 t p ) ) · e x p ( j 2 π R 0 λ ) .
Assume that the number of azimuth samples is N a and the time domain output sample spacing is
Δ t = p r f N a · K a ,
where p r f is the pulse repetition frequency, then K a is the azimuth frequency modulation factor, which varies with range, which results in uneven spacing of output samples and distortion in final imaging results. Its expression is as follows
K a = r a ( r p r a ) ω 2 λ R 0 ,
In general, image distortion is tolerable if the distance r p varies within a small range. But for the RoASAR system which operates at low altitudes, the distance r p changes rapidly as the beam moves from the edge to the center. In addition, because Δ t is proportional to R 0 , target located at near range are elongated, and it is compressed at far range. Therefore, the non-uniform spacing sampling of azimuth will seriously affect the image quality.
The efficient CZT uses the scaling Fourier transform (SCFT), instead of the Fourier transform, to solve the non-uniform sampling problem [21,22,23,24]. The compression factor is
β = R s R 0 · ( r p r a ) r s ,
where r s = R s 2 H 2 represents the reference ground position.
The scaling factor is a function of range and requires performing a scaling operation at each range gate. After the scaling operation, the output sampling interval of each distance gate is uniform, and thus, an accurate focusing image can be obtained.
First, we remove the high-order phase term from the range processing result, and then replace the hyperbolic phase with the quadratic phase. The compensation function needed is
H 2 ( t r ; f a ; R 0 ) = e x p [ j Δ ϕ 0 ( f a ; R 0 ) ] ,
where
Δ ϕ 0 ( f a ; R 0 ) = π c 3 k 4 8 k 2 4 f c 3 f a 4 .
After multiplying the range processing results of Equation (21) by (26), we transform the remaining signals into the time domain, given by
S 4 ( t r ; t a ) = w r ( t r R 0 c ) w a ( t a ) e x p ( 2 π R 0 λ ) e x p [ j π 2 r a ( r p r a ) ω 2 λ R 0 ( t a t p ) 2 ] .
Based on the SPECAN algorithm, the Deramp processing is performed, and the azimuth-to-Deramp function is used to eliminate the quadratic phase. The Deramp function expression is
H 3 ( t r ; t a ; R 0 ) = e x p ( j π 2 r a ( r p r a ) ω 2 λ R 0 ( t a t p ) 2 ) .
Next, we propose to transform the signal into the range-Doppler domain by CZT to avoid interpolation. The kernel expression of CZT is
h C Z T = e x p ( j 2 π 1 β f a t a ) .
The final focusing image obtained is given by
S C Z T ( t r ; f a ) = S C F T { S 4 ( t r ; t a ) · H 4 } = w r ( t r R 0 c ) w a ( π F a N a ( f a β 2 r a r s ω 2 λ R s t p ) ) e x p ( 2 π R 0 λ ) .
The flowchart of the whole algorithm is shown in Figure 6.

3.4. Some Considerations of the Proposed Method in Applications

We now compare the computational complexity of our proposed method with the BP algorithm in [6] and Omega-k algorithm in [13]. The proposed algorithm avoids interpolation, and all the operations in the algorithm are realized by FFT and complex phase multiplications, which means it is of high efficiency. To be specific, we estimate the floating-point operations (FLOPs) of aforementioned algorithms. Generally speaking, an N -point FFT or inverse FFT (IFFT) need 5 N log 2 ( N ) FLOPs, and 6 N FLOPs are needed for one-time complex multiplication. Suppose that the range sample number is N r , the azimuth sample number is N a and the length of the interpolator kernel is N k e r . In this case, the computational complexity of BP algorithm in [6] is
C BP = 10 N a N r log 2 ( N r ) + 6 N a N r + N a N r ( 10 N a N r log 2 ( N r ) + 6 N a N r ) .
The computational load in terms of FLOPs for the Omega-k method in [13] is
C Omega k = 20 N a N r log 2 ( N r ) + 20 N a N r log 2 ( N a ) + 18 N a N r + 2 ( 2 N k e r 1 ) N a N r .
According to the flowchart of the proposed method in Figure 6, the computational load of the proposed method is mainly composed of the following steps. First, the 2-D frequency is obtained by 2D-FFT and the computational complexity is 5 N a N r log 2 ( N r ) + 5 N a N r log 2 ( N a ) FLOPs. In range processing, the bulk RCM, SRC and cubic cross-coupling terms are compensated with the reference distance, directly, and the computational complexity is 6 N a N r FLOPs. Then, the range-dependent RCM term is removed by using an ICZT operation, which contains twice phase multiplications and one convolution operation. The computational complexity is 10 N a N r log 2 ( N r ) + 18 N a N r FLOPs. In azimuth processing, the implemented procedure mainly includes twice phase multiplications, one FFT and one CZT. In these steps, the computational complexity is 15 N a N r log 2 ( N a ) + 30 N a N r FLOPs. To summarize, we combine the complex multiplication together, and the total computation load in terms of FLOPs for the proposed method is approximately expressed as
C proposed = 15 N a N r log 2 ( N r ) + 20 N a N r log 2 ( N a ) + 36 N a N r .
It is noted that the BP and Omega-k algorithms that contain the correlation operation point-by-point and interpolations operation, respectively, have a higher computational load than the proposed algorithm, which means the proposed algorithm is more suitable to real-time imaging processing.
The motion errors generated from the atmospheric turbulence and platform tremble are contributing factors in the most airborne platforms, particularly in helicopter-borne RoASAR systems. The two-step approach in [25] and [26] is a common airborne motion compensation algorithm and can be easily integrated with the proposed 2-D CZT algorithm. First, the phase error correction for a reference range is defined and then it is carried out directly with range uncompressed data. Second, the phase error correction of range dependence is defined, which is carried out after the range cell migration correction (RCMC). In our system, the motion compensation step is also fused in the proposed method that can reduce the effects of the caused errors, depicted in Figure 6. Moreover, the data-driven autofocusing techniques in [8,9] are also good candidate methods for the helicopter-borne RoASAR motion compensation, and its will be our future study.

4. Simulation Results and Analysis

In this section, simulation results are presented to demonstrate the effectiveness and performance of the proposed algorithm for high resolution helicopter-borne RoASAR, where the simulated data processing parameters are the same scene as the original data in Table 1. For fairness of comparison, no weighting window and sidelobe control approach are adopted here. A circular imaging region with twenty-five targets is adopted, as shown in Figure 7. To compare the performance, the results of the BP algorithm in [6], Omega-k method [13], and SPECAN algorithm in [20] are also provided.
For the edge points of the imaging scene, the results of the bulk range compensated echo with the phase of the reference point and the ICZT results proposed in this paper are respectively shown in Figure 8. It can be seen from Figure 8a,b that the RCM is well corrected by ICZT, and the main energy is well concentrated into one range gate.
Meanwhile, to demonstrate the effectiveness of the proposed algorithm, the comparisons with the BP algorithm in [6], Omega-k algorithm in [13] and SPECAN algorithm in [20] are also provided. The imaging result of the BP algorithm and the Omega-k algorithm is depicted in Figure 9a,b. Both of them obtained a well-focused image, but suffering from a great computational load due to the time-domain correlation point-by-point and Stolt interpolation operation, respectively, which are not suitable for real-time imaging processing. The standard SPECAN algorithm adopts the azimuth Deramp operation, which ignores the influence of the range coupling term on the position modulation frequency, and the imaging result is depicted in Figure 9c. The imaging result shows an obvious sector shaped distortion, which is inconsistent with the original actual scene. The image result obtained by the proposed algorithm is shown in Figure 9d. It can be seen that the sector distortion caused by the Deramp operation has been well corrected, verifying the correctness and effectiveness of the proposed algorithm.
To further demonstrate the performance of the proposed imaging algorithm, we present the magnified edge point target imaging results comparison in Figure 10. The range and azimuth profile of the BP algorithm and the Omega-k algorithm are depicted in Figure 10a–d, respectively. And the range and azimuth profile of the proposed method is provided in Figure 10e,f. It can be found that, from Figure 10, the well-focused SAR image can be obtained with the three methods.
The performance measures include displacement, integral sidelobe ratio (ISLR), peak sidelobe ratio (PSLR) and resolution, and are shown in Table 2 and Table 3. The true positions of targets are provided in Figure 7. It can be seen that the imaging accuracy of our algorithm is higher than the standard SPECAN algorithm, but slightly lower than the traditional time-domain algorithm and Omega-k algorithm. The comparison of imaging quality shows the effectiveness of proposed algorithm. However, the proposed method has the advantage in computational complexity without interpolation operation.
In Figure 11, the running times of different approaches versus imaging scene size are presented, where we used an Intel quad-core processor, CPU clocked frequency at 3.2 GHz, memory of 8 GB, and a Windows 10 operating system. It confirms that the computational times of the proposed method and the Omega-k algorithm are much less than the BP algorithm. And the computational cost of the proposed method is less than the Omega-k algorithm, since the interpolation operation is avoided, which is also consistent with the theoretical analysis. As the size of RoASAR imaging scene increases, the advantage of the proposed method in terms of running time is more obvious. To conclude, the proposed method can provide a high-resolution imaging result for helicopter-borne RoASAR, with a fast-computational speed in comparison with the BP algorithm in [6], fast factorized back-projection (FFBP) algorithm in [7], and Omega-k algorithm in [13]. This again demonstrates the proposed method is a promising candidate for practical applications.

5. Conclusions

Aiming at solving the problems of limited azimuth resolution and high computational complexity of a helicopter RoASAR system, an improved 2-D CZT imaging algorithm is proposed. Based on the fourth-order expansion of instantaneous slant range, the approximated phase equation of slant range tolerates a larger azimuth reconstruction angle, and thus, the azimuth resolution is improved. The improved ICZT is used to eliminate RCM, where FFT in the standard SPECAN algorithm is replaced, which avoids the interpolation operation to eliminate the sector distortion and reduces the computational complexity. Finally, simulation results are provided to verify the effectiveness of the algorithm.

Author Contributions

M.Q. and D.L. proposed the method, conceived and designed the experiments, and analyzed the data; X.T. and C.Z. performed the experiments and wrote the paper; and W.L. and L.X. revised the paper.

Funding

This work was funded by the National Key R&D Program of China (grant no. 2016YFE0200400), by the National Natural Science Foundation of China (grant no. 61771015 and 61801294), by the Graduate Scientific Research and Innovation Foundation of Chongqing (grant no. CYS19019), by the Key R&D Program of ShanXi Province (grant no. 2017KW-ZD-12), by the Innovative Research Group of National Natural Science Foundation of China (grant no. 61621005), by the Chongqing Research Program of Basic Research and Frontier Technology(cstc2018jcyjAX0351), by the Shanghai Aerospace Science and Technology Innovation Fund (grant no. SAST2017041), by the Guangxi Key Laboratory of Wireless Wideband Communication and Signal Processing, the 2017 Fund Project of the Director (grant no. GXKL06170202), and by the Fundamental Research Funds for the Central Universities (grant no. 2019CDQYTX012).

Acknowledgments

The authors would like to thank all the anonymous reviewers and editors for their useful comments and suggestions that greatly improved this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

We now give a detailed derivation process of the Taylor series expansion. The mathematical expression of the Taylor series expansion is as follows.
f ( x ) = n = 0 f ( n ) ( x ) n ! ( x x 0 ) n | x = x 0 .
Based on (A1), the instantaneous slant range can be expressed as
R ( t a ) R ( t P ) + R ( 1 ) ( t p ) ( t a t p ) 1 ! + R ( 2 ) ( t p ) ( t a t p ) 2 2 ! + R ( n ) ( t p ) ( t a t p ) n n ! | t a = t P .
where
R ( t a ) = r p 2 + r a 2 + H 2 2 r a r p cos ω ( t a t p ) .
The first derivative of instantaneous slant range is
R ( 1 ) ( t a ) = r a r p ω sin ω ( t a t p ) r p 2 + r a 2 + H 2 2 r a r p cos ω ( t a t p ) = ω r a r p sin ω ( t a t p ) R ( t a ) .
Similarly, the second derivative of instantaneous slant range is
R ( 2 ) ( t a ) = ω 2 r a r p cos ω ( t a t p ) R ( t a ) ω 2 r a 2 r p 2 sin 2 ω ( t a t p ) R 3 ( t a ) .
The third derivative of instantaneous slant range is
R ( 3 ) ( t a ) = ω 3 r a r p sin ω ( t a t p ) R ( t a ) 3 ω 3 r a 2 r p 2 sin [ ω ( t a t p ) ] cos [ ω ( t a t p ) ] R 3 ( t a ) + 3 ω 3 r a 3 r p 3 sin 3 ω ( t a t p ) R 5 ( t a ) .
Finally, the fourth derivative of instantaneous slant range is
R ( 4 ) ( t a ) = ω 4 r a r p cos ω ( t a t p ) R ( t a ) 3 ω 4 r a 2 r p 2 cos [ ω ( t a t p ) ] 4 ω 4 r a 2 r p 2 sin [ ω ( t a t p ) ] R 3 ( t a ) = 18 ω 4 r a 3 r p 3 sin 2 [ ω ( t a t p ) ] cos [ ω ( t a t p ) ] R 5 ( t a ) 15 ω 4 r a 4 r p 4 sin 4 ω ( t a t p ) R 7 ( t a ) .
When t a = t P , the instantaneous slant range and their derivatives are
R ( t P ) = r p 2 + r a 2 + H 2 R ( 1 ) ( t P ) = 0 R ( 2 ) ( t P ) = ω 2 r a r p R ( t P ) R ( 3 ) ( t P ) = 0 R ( 4 ) ( t P ) = ω 4 r a r p R ( t P ) 3 ω 4 r a 2 r p 2 R 3 ( t P ) .
Using the results obtained, the second-order approximation of instantaneous slant range is given in (A9), where the approximation error is corresponding to Figure 4a in the paper.
R ( t a ) R ( t P ) + 1 2 r a r p ω 2 R ( t P ) ( t a t p ) 2 .
And the fourth-order approximation of instantaneous slant range is given in (A10), where the approximation error is corresponding to Figure 4b in the paper.
R ( t a ) R ( t P ) + 1 2 r a r p ω 2 R ( t P ) ( t a t p ) 2 + 1 24 ( ω 4 r a r p R ( t P ) 3 ω 4 r a 2 r p 2 R 3 ( t P ) ) ( t a t p ) 4 = R ( t P ) + r a r p ω 2 2 R ( t P ) ( t a t p ) 2 ( ω 4 r a r p 24 R ( t P ) ω 4 r a 2 r p 2 8 R 3 ( t P ) ) ( t a t p ) 4 .
Using (A10), the fourth-order Taylor series expansion of two-way slant range can be expressed as Equation (5) in the article.

Appendix B

The MSR approach was originally developed in [16]. Based on the idea of MSR, the point target signal in the time domain after range compression is
S 1 ( t r ; t a ) = w r [ t r 2 R ( t a ; r p ) c ] w a ( t a ) e x p [ j 4 π λ R ( t a ; r p ) ] .
The fourth-order approximation of instantaneous slant range is
2 R ( t a ; r p ) = R 0 + k 2 ( t a t p ) 2 + k 4 ( t a t p ) 4 ,
where the derivation of the coefficients is given in the previous response
{ R 0 = 2 R p = 2 ( r p r a ) 2 + H 2 k 2 = r a r p ω 2 R p k 4 = ( r a r p ω 4 12 R p + r a 2 r p 2 ω 4 4 R p 3 ) .
Expanding the point target signal up to the fourth azimuth frequency terms, the 2-D point target spectrum is given by
S 2 ( f r ; t a ) = W r ( f r ) w a ( t a ) e x p [ j 2 π ( f r + f c ) 2 R ( t a ; r p ) c ] .
Next, we perform an azimuth FT. Using the method of stationary phase, azimuth frequency is related to azimuth time by
( c f r + f c ) f a = 2 k 2 ( t a t p ) + 4 k 4 ( t a t p ) 3 .
We can derive an expression of t a in terms of f a by using the series reversion. Using the results in the Appendix in reference [16], we arrive at
t a ( f a ) t p = A 1 ( c f r + f c f a ) + A 2 ( c f r + f c f a ) 2 + A 3 ( c f r + f c f a ) 3 ,
where
A 1 = ( 2 k 2 ) 1 A 2 = 0 A 3 = ( 2 k 2 ) 5 ( 0 2 k 2 · 4 k 4 ) .
Substituting (A16) to (A14), we obtain the 2-D spectrum of S 1 ( t r ; t a ) , given by
S 3 ( f r ; f a ) = W r ( f r ) W a ( f a ) · e x p ( 2 π f a t p ) · e x p [ j Φ ( f r ; f a ) ] ,
where the phase term Φ ( f r ; f a ) is expressed as
Φ ( f r ; f a ) = 2 π ( f r + f c ) c R 0 π f r 2 K r + π c 2 k 2 ( f r + f c ) f a 2 π c 3 k 4 2 k 2 4 ( f r + f c ) 3 f a 4 .

References

  1. Kreitmair-Steck, W.; Wolframm, A.P. Heliradar-a Rotating Antenna Synthetic Aperture Radar for Helicopter all Weather Operations. In Proceedings of the AGARD 59th Guidance and Control Panel, Pratica di Mare, Italy, 20–25 October 1994; pp. 56–62. [Google Scholar]
  2. Joelianto, E.; Budiyono, A.; Wijayanti, I.E.; Megawati, N.Y. Model Predictive Control for Obstacle Avoidance as Hybrid Systems of Small Scale Helicopter. In Proceedings of the 2013 3rd International Conference on Instrumentation Control and Automation, Ungasan, Indonesia, 28–30 August 2013. [Google Scholar]
  3. Kwag, Y.K.; Kang, J.W. Obstacle Awareness and Collision Avoidance Radar Sensor System for Low-Altitude Flying Smart UAV. In Proceedings of the 23rd Digital Avionics Systems Conference, Salt Lake City, UT, USA, 28–28 October 2004. [Google Scholar]
  4. Kreitmair, W. Heliradar: The Pilot’s Eye for Flights in Adverse Weather Conditions. In SPIE 1995 Aerosense Symposium; Enhanced and Synthetic Vision: Orlando, FL, USA, 1995. [Google Scholar]
  5. Guo, H.D.; Li, X.W. Characteristics and application expansion of new generation SAR earth observation technology. Chin. Sci. Bull. 2011, 56, 1155–1168. [Google Scholar] [CrossRef]
  6. Jakowatz, C.V., Jr.; Wahl, D.E.; Yocky, D.A. Beamforming as a Foundation for Spotlight-Mode SAR Image Formation by Backprojection. In Proceedings of the SPIE-The International Society for Optical Engineering, Orlando, FL, USA, 15 April 2008; Volume 6970. [Google Scholar]
  7. Ulander, L.M.; Hellsten, H.; Stenstrom, G. Synthetic-aperture radar processing using fast factorized back-projection. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 760–776. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, L.; Li, H.-L.; Qiao, Z.-J.; Xing, M.D.; Bao, Z. Integrating Autofocus Techniques with Fast Factorized Back-Projection for High-Resolution Spotlight SAR Imaging. IEEE Geosci. Remote Sens. Lett. 2013, 10, 1394–1398. [Google Scholar] [CrossRef]
  9. Zhou, S.; Yang, L.; Zhao, L.; Bi, G. Quasi-Polar-Based FFBP Algorithm for Miniature UAV SAR Imaging Without Navigational Data. IEEE Trans. Geosci. Remote Sens. 2017, 55, 7053–7065. [Google Scholar] [CrossRef]
  10. Chen, S.B.; Chun-sheng, J.L. Operation mode of circular trace scanning SAR for wide observation. J. Electron. Inf. Technol. 2008, 30, 2805–2808. [Google Scholar]
  11. Li, D.; Liao, G.S.; Xu, Q.; Wang, W. High resolution imaging algorithm for helicopter-borne rotating synthetic aperture radar. J. Syst. Eng. Electron. 2013, 35, 1389–1395. [Google Scholar]
  12. Zhang, J.; Liao, G.S.; Xu, J.W.; Yang, K.; Zhu, S. Research on high resolution wavenumber-domain imaging algorithm for ROSAR. J. Syst. Eng. Electron. 2017, 39, 1127–1234. [Google Scholar]
  13. Li, D.; Liao, G.; Yang, L. Modified Omega-K algorithm for processing helicopter-borne frequency modulated continuous waveform rotating synthetic aperture radar data. J. Syst. Eng. Electron. 2015, 26, 476–485. [Google Scholar] [CrossRef]
  14. Adams, J.W.; Bayma, R.W.; Lawrence, M.E.; Petrosian, L. On the modified chirp z transform for synthetic aperture radar systems. IEEE Trans. Signal Process. 1991, 39, 953–955. [Google Scholar] [CrossRef]
  15. Li, D.; Liu, H.; Liao, Y.; Gui, X. A Novel Helicopter-Borne Rotating SAR Imaging Model and Algorithm Based on Inverse Chirp-Z Transform Using Frequency-Modulated Continuous Wave. IEEE Geosci. Remote Sens. Lett. 2015, 12, 1–5. [Google Scholar]
  16. Neo, Y.L.; Wong, F.; Cumming, I.G. A two-dimensional spectrum for bistatic SAR processing using series reversion. IEEE Geosci. Remote Sens. Lett. 2007, 4, 93–96. [Google Scholar] [CrossRef]
  17. Lanari, R. A new method for the compensation of the SAR range cell migration based on the chirp-z transform. IEEE Trans. Geosci. Remote Sens. 1995, 33, 1296–1299. [Google Scholar] [CrossRef]
  18. Zhou, P.; Zhou, S.; Xiong, T.; Li, Y.C.; Xing, M.D. A chirp-z transform imaging algorithm for missile-borne SAR with diving maneuver based on the method of series reversion. J. Electron. Inf. Technol. 2010, 32, 2861–2867. [Google Scholar] [CrossRef]
  19. Li, C.; He, M. Imaging algorithm for geosynchronous orbit SAR based on chirp-z transform and azimuth scaling. J. Electron. Inf. Technol. 2015, 37, 1736–1742. [Google Scholar]
  20. Cumming, I.G.; Wong, F.H. Digital Processing of Synthetic Aperture Radar Data, 2nd ed.; Artech House: Norwood, MA, USA, 2005; pp. 249–260. [Google Scholar]
  21. Yang, W.; Chen, J.; Zeng, H.; Wang, P.; Liu, W. A wide-swath spaceborne TOPS SAR image formation algorithm based on chirp scaling and chirp-z transform. Sensors 2016, 16, 2095. [Google Scholar] [CrossRef] [PubMed]
  22. Zhong, H.; Zhang, Y.; Chang, Y.; Liu, E.; Tang, X.; Zhang, J. Focus High-Resolution Highly Squint SAR Data Using Azimuth-Variant Residual RCMC and Extended Nonlinear Chirp Scaling Based on a New Circle Model. IEEE Geosci. Remote Sens. Lett. 2018, 15, 547–551. [Google Scholar] [CrossRef]
  23. Bi, H.; Zhang, B.; Zhu, X.; Hong, W. Azimuth-range decouple-based L 1 regularization method for wide ScanSAR imaging via extended chirp scaling. J. Appl. Remote Sens. 2017, 11, 15007. [Google Scholar] [CrossRef]
  24. Wu, Y.; Sun, G.-C.; Xia, X.-G.; Xing, M.; Yang, J.; Bao, Z. An Azimuth Frequency Non-Linear Chirp Scaling (FNCS) Algorithm for TOPS SAR Imaging with High Squint Angle. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 213–221. [Google Scholar]
  25. Su, Y.X.; Li, D. Motion compensation and imaging algorithm for helicopter-borne rotating synthetic aperture radar. Radar Sci. Technol. 2018, 16, 383–390. [Google Scholar]
  26. Xing, M.; Jiang, X.; Wu, R.; Zhou, F.; Bao, Z. Motion Compensation for UAV SAR Based on Raw Radar Data. IEEE Trans. Geosci. Remote Sens. 2009, 47, 2870–2883. [Google Scholar] [CrossRef]
Figure 1. Rotating array synthetic aperture radar (RoASAR) imaging geometric configuration. (a) RoASAR 3-D imaging geometric model. (b) RoASAR front view image.
Figure 1. Rotating array synthetic aperture radar (RoASAR) imaging geometric configuration. (a) RoASAR 3-D imaging geometric model. (b) RoASAR front view image.
Remotesensing 11 01669 g001
Figure 2. Diagram of the beam radiation trajectory for the target P.
Figure 2. Diagram of the beam radiation trajectory for the target P.
Remotesensing 11 01669 g002
Figure 3. Variation curves of azimuth angular resolution. (a) Azimuth angular resolution varying with r a and H . (b) Azimuth angular resolution varying with reconstruction angle α . (c) Azimuth angular resolution versus observation angle θ .
Figure 3. Variation curves of azimuth angular resolution. (a) Azimuth angular resolution varying with r a and H . (b) Azimuth angular resolution varying with reconstruction angle α . (c) Azimuth angular resolution versus observation angle θ .
Remotesensing 11 01669 g003aRemotesensing 11 01669 g003b
Figure 4. Instantaneous slant range approximation error. (a) Second-order approximation. (b) Fourth-order approximation.
Figure 4. Instantaneous slant range approximation error. (a) Second-order approximation. (b) Fourth-order approximation.
Remotesensing 11 01669 g004
Figure 5. Range-dependent characteristic curves versus slant range. (a) Cubic cross-coupling term ϕ 3 . (b) SRC term ϕ 2 . (c) RCM term ϕ 1 .
Figure 5. Range-dependent characteristic curves versus slant range. (a) Cubic cross-coupling term ϕ 3 . (b) SRC term ϕ 2 . (c) RCM term ϕ 1 .
Remotesensing 11 01669 g005
Figure 6. The flowchart of the algorithm presented in this paper.
Figure 6. The flowchart of the algorithm presented in this paper.
Remotesensing 11 01669 g006
Figure 7. Imaging scene distribution.
Figure 7. Imaging scene distribution.
Remotesensing 11 01669 g007
Figure 8. (a) Edge point range cell migration (RCM) result by bulk range compensation. (b) Edge point RCM correction result by inverse chirp-z transform (ICZT).
Figure 8. (a) Edge point range cell migration (RCM) result by bulk range compensation. (b) Edge point RCM correction result by inverse chirp-z transform (ICZT).
Remotesensing 11 01669 g008
Figure 9. Comparison of the image results. (a) Back-projection (BP) algorithm in [6]. (b) Omega-k algorithm in [13]. (c) SPECAN algorithm in [20]. (d) Result by the proposed method.
Figure 9. Comparison of the image results. (a) Back-projection (BP) algorithm in [6]. (b) Omega-k algorithm in [13]. (c) SPECAN algorithm in [20]. (d) Result by the proposed method.
Remotesensing 11 01669 g009
Figure 10. Magnified point target imaging results comparison. (a) Range profile of the BP algorithm. (b) Azimuth profile of the BP algorithm. (c) Range profile of the Omega-k algorithm. (d) Azimuth profile of the Omega-k algorithm. (e) Range profile of our proposed method. (f) Azimuth profile of our proposed method.
Figure 10. Magnified point target imaging results comparison. (a) Range profile of the BP algorithm. (b) Azimuth profile of the BP algorithm. (c) Range profile of the Omega-k algorithm. (d) Azimuth profile of the Omega-k algorithm. (e) Range profile of our proposed method. (f) Azimuth profile of our proposed method.
Remotesensing 11 01669 g010aRemotesensing 11 01669 g010b
Figure 11. The running times versus RoASAR imaging scene size.
Figure 11. The running times versus RoASAR imaging scene size.
Remotesensing 11 01669 g011
Table 1. Simulation Parameters.
Table 1. Simulation Parameters.
Parameters NameValue
Rotator length2 m
Rotating angle speed15 rad/s
Platform altitude1000 m
Transmit bandwidth300 MHz
Carrier frequency10 GHz
Pulse repetition frequency10 kHz
Reference slant range2000 m
Azimuth beam width80°
Scene size (range*azimuth)700 m*70°
Table 2. Measurement of point displacement.
Table 2. Measurement of point displacement.
ApproachRange Displacement (m)Azimuth Displacement (°)
PT1PT2PT3PT4PT1PT2PT3PT4
BP−0.240+0.2500−0.0530+0.055
Omega-k−0.240+0.2500−0.0550+0.059
SPECAN−0.520+0.2300−0.1100+0.109
Proposed method−0.250+0.2500−0.0570+0.061
Table 3. Evaluation of the imaging performance.
Table 3. Evaluation of the imaging performance.
ApproachRangeAzimuth
PSLR (dB)ISLR (dB)Resolution (m)PSLR (dB)ISLR (dB)Resolution (°)
BP−13.29−9.760.55−13.22−9.720.1154
Omega-k−13.24−9.690.59−12.74−9.680.1263
Proposed method−13.22−9.680.59−12.96−9.640.1290
theoretical−13.4−9.800.50−13.4−9.800.1117

Share and Cite

MDPI and ACS Style

Qin, M.; Li, D.; Tang, X.; Zeng, C.; Li, W.; Xu, L. A Fast High-Resolution Imaging Algorithm for Helicopter-Borne Rotating Array SAR Based on 2-D Chirp-Z Transform. Remote Sens. 2019, 11, 1669. https://doi.org/10.3390/rs11141669

AMA Style

Qin M, Li D, Tang X, Zeng C, Li W, Xu L. A Fast High-Resolution Imaging Algorithm for Helicopter-Borne Rotating Array SAR Based on 2-D Chirp-Z Transform. Remote Sensing. 2019; 11(14):1669. https://doi.org/10.3390/rs11141669

Chicago/Turabian Style

Qin, Mengyi, Dong Li, Xianhui Tang, Cao Zeng, Wenjun Li, and Liying Xu. 2019. "A Fast High-Resolution Imaging Algorithm for Helicopter-Borne Rotating Array SAR Based on 2-D Chirp-Z Transform" Remote Sensing 11, no. 14: 1669. https://doi.org/10.3390/rs11141669

APA Style

Qin, M., Li, D., Tang, X., Zeng, C., Li, W., & Xu, L. (2019). A Fast High-Resolution Imaging Algorithm for Helicopter-Borne Rotating Array SAR Based on 2-D Chirp-Z Transform. Remote Sensing, 11(14), 1669. https://doi.org/10.3390/rs11141669

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