Next Article in Journal
The Penetration Analysis of Airborne Ku-Band Radar Versus Satellite Infrared Lidar Based on the Height and Energy Percentiles in the Boreal Forest
Previous Article in Journal
Robust Controller for Pursuing Trajectory and Force Estimations of a Bilateral Tele-Operated Hydraulic Manipulator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Fast Back-Projection SAR Imaging Algorithm Based on Wavenumber Spectrum Fusion for High Maneuvering Platforms

1
National Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
2
Sichuan Aerospace Electronic Equipment Research Institute, Chengdu 610100, China
*
Author to whom correspondence should be addressed.
Submission received: 8 March 2021 / Revised: 14 April 2021 / Accepted: 19 April 2021 / Published: 23 April 2021

Abstract

:
In this paper, a fast back-projection (BP) synthetic aperture radar (SAR) imaging algorithm based on wavenumber spectrum fusion is proposed for high maneuvering platforms with high squint angles and high dive angles. When the platforms have nonlinear trajectory errors, the trajectory errors will introduce significant spectrum displacement, which varies with the target position. If the influence of trajectory error is not considered, the trajectory’s deviation from the ideal SAR imaging geometry will degrade the focusing performance. The dive angle will further enlarge the influence on SAR focusing. In this paper, the center function of the spatially variable sub-aperture wavenumber spectrum is derived, firstly, which is the basis for reasonably setting the oversampling rate. Then, the accurate spectrum fusion method is proposed, which is efficiently implemented by fast Fourier transform (FFT) and circular shifting. The proposed algorithm is called high-squint and high-dive accelerated factorized back-projection (HSHD-AFBP). Compared to the commonly used spectrum fusion method, HSHD-AFBP not only improves the imaging quality, but also expands the image width. Finally, the effectiveness of the proposed algorithm is verified by simulation data.

1. Introduction

As a powerful microwave imaging system, synthetic aperture radar (SAR) is widely concerned in the military and civil detection fields because of its advantages for all-time and all-weather working ability [1,2]. As is well known, SAR achieves range resolution by transmitting large time-bandwidth product signals, and azimuth resolution by the synthetic aperture, which is formulated by the relative motion of the platform and the target to obtain the azimuth-Doppler bandwidth. The high maneuvering platform has the characteristics of fast flight speed and high maneuverability, which can realize the requirement of advanced observation. With the rapid development of SAR applications, high maneuvering platforms SAR imaging technology has become a hot issue in SAR imaging application [3,4,5,6].
At present, SAR imaging algorithms are usually divided into two categories: frequency domain algorithms and time domain algorithms.
Traditional frequency domain algorithms, such as the range-Doppler algorithm (RDA) [7], chirp scaling algorithm (CSA) [8,9,10], nonlinear chirp scaling algorithm (NCSA) [11], many extensions of the NCSA (ENCSA) [12,13], etc., are widely used and deeply investigated. These methods use the azimuth translation-invariance of echo data to uniformly deal with the spatial variation of echo data in the range of the Doppler domain, so the efficiency of the algorithm is high. However, if the radar platform’s moving trajectory has a large dive angle and squint angle, simultaneously combining with high maneuvering ability, the radar echoes have inherent spatial variation, and the high-squint angle observation mode gives the echo data a serious coupling effect between range and azimuth, as well as Doppler folding. Therefore, the traditional frequency domain algorithms cannot be directly applied to high maneuvering platforms.
The widely used time domain algorithm is the back-projection algorithm (BPA) [14,15,16,17,18], which is suitable for arbitrary platforms and arbitrary trajectories, and adopts the method of point-by-point coherent accumulation to achieve high-precision SAR imaging with accurate motion compensation. The disadvantage of BPA is the intensive computational burden and that the motion trajectory is required to have high accuracy. To improve SAR imaging efficiency, two improved BP algorithms—the fast BP (FBP) [15] algorithm and the fast factorized BP (FFBP) algorithm [17]—have recently been developed. During sub-image fusion, the FFBP often needs two-dimensional image-domain complex interpolation for SAR imaging performance conservation, so it still produces a huge amount of calculation. Considering imaging efficiency and imaging performance, the accelerated fast BP algorithm (AFBP) is developed in [18]. Sub-aperture wavenumber spectrum fusion is carried out in the wavenumber spectrum domain by using a unified polar coordinate system. Compared with the conventional FBP algorithm, the AFBP can precisely focus more efficiently by means of the implementation of fast Fourier transform (FFT) and circular shifting for wavenumber domain SAR fusion. However, AFBP is suitable for a broadside mode, and it fails to correct the sub-image spectrum center accurately in high-squint SAR. Moreover, AFBP cannot deal with a nonlinear trajectory, which limits its practical application in the SAR field. The high-squint AFBP algorithm (HS-AFBP) for airborne SAR processing is proposed in [19], and two spectrum center functions are also given for linear and nonlinear trajectories. In [19], the radar platforms’ motion error is modeled as a linear function for two-dimensional spatial variance. The angular-variant coefficient, generally being a constant, is obtained by first-order fitting of the motion errors in the central range cell. However, the motion errors introduce the displacement of the sub-aperture spectrum centers of point targets for different angular positions even in the same range cell. That is to say, the angular-variant coefficient varies with the target position. Therefore, the linear function model of the two-dimensional spatial variance cannot be used accurately to represent the real motion errors. Moreover, this approximation error of the linear function model will increase with the motion errors and the dive angles’ increment.
In this paper, the high-squint and high-dive AFBP (HSHD-AFBP) is proposed for high maneuvering platforms SAR imaging. The main contributions can be specified as follows. Based on the high maneuvering platforms SAR imaging geometry, the relationship of the sub-aperture wavenumber spectrum centroid position and the target position is firstly derived. It is found that the center position of the sub-aperture wavenumber spectrum varies with the target position if the nonlinear motion errors occur for the high maneuvering platforms. The larger the motion errors and trajectory dive angles, the larger the center position variations. Secondly, an accurate sub-aperture spectrum fusion method is proposed, which eliminates spectrum folding by formulating the accurate spectrum shift function. The proposed spectrum fusion method can overcome the spatial variation influence of the center position of the sub-aperture spectrum, and then the SAR focusing performance can be guaranteed effectively.
This paper is organized as follows. In Section 2, the SAR signal model of high maneuvering platforms is established. In Section 3, the AFBP and HS-AFBP are briefly reviewed, followed by a detailed shortage description of the existed methods. In Section 4, the details of the proposed HSHD-AFBP algorithm are given. Firstly, the more accurate wavenumber spectrum center function of sub-aperture is derived and the constraint of the motion errors for the length of the sub-aperture is discussed. Then, the improved spectrum fusion method is proposed. Next, the integrated implementation process of the HSHD-AFBP algorithm is introduced. Finally, its computation load is discussed. In Section 5, the proposed algorithm is verified based on simulation experiments. Section 6 summarizes this paper.

2. Signal Model

The geometry of high maneuvering platforms SAR is shown in Figure 1. All of the parameters in Figure 1a are defined in the local ground coordinate system. The origin o of the coordinate system is the projection point of the aperture center on the ground. The x o y plane is parallel to the surface of the Earth. The positive direction of the x -axis is the projection direction of the platform’s speed in the x o y plane. The solid line with arrow represents a true trajectory, and dashed lines with arrows indicate a nominal trajectory. The ideal trajectory is obtained by one-dimensional fitting of the actual trajectory. The motion speed of the ideal trajectory is v . The nominal trajectory intersects the z -axis at M and the x -axis at N . γ represents the angle between the speed direction and x axis, i.e., the dive angle. α p represents the angle between the speed direction and the straight line, MP . θ p represents the angle between the x -axis and the straight line, OP . The platform’s height is H .
We denote τ as the slow time. At an arbitrary azimuth time, the ideal antenna phase center position of uniform linear motion can be expressed as [ X ( τ ) , 0 , Z ( τ ) ] with X ( τ ) = v cos γ τ , Z ( τ ) = H v sin γ τ . The actual antenna phase center position can be expressed as [ X ( τ ) + Δ X ( τ ) , Δ Y ( τ ) , Z ( τ ) + Δ Z ( τ ) ] . For the convenience of formula derivation, X ( τ ) , Z ( τ ) , Δ X ( τ ) , Δ Y ( τ ) , Δ Z ( τ ) will be abbreviated as X , Z , Δ X , Δ Y , ΔZ, in turn, in the following. The actual instantaneous slant range of an arbitrary point target P ( x p , y p , 0 ) in the scene can be expressed as:
R ( τ ; x p , y p ) = ( X + Δ X x p ) 2 + ( Δ Y y p ) 2 + ( Z + Δ Z ) 2
The MNP plane is the imaging slant plane of target P . Figure 1b shows the polar coordinate system established on the MNP plane. The whole-aperture length is L . The center of the aperture and the origin of the polar coordinate system is M . Point target P can be expressed as ( r p , Θ p ) in the polar coordinate system, denoted as Θ p = cos α p . Then, the conversion relationship between the point target P in the Cartesian coordinate system and the polar coordinate system can be established:
{ x p = ( r p Θ p H sin γ ) cos γ y p = r p 2 H 2 x p 2
Substituting Equation (2) into Equation (1), we can get:
R ( τ ; x p , r p ) = ( X + Δ X x p ) 2 + ( Δ Y r p 2 H 2 x p 2 ) 2 + ( Z + Δ Z ) 2
Ideally, the platform flies along the linear trajectory with a constant speed, v . In the polar coordinate system, the instantaneous slant range between the target P and the ideal trajectory can be expressed as [18]:
R ( r p , Θ p ; D ) = r p 2 + D 2 2 r p D Θ p ,   L 2 < D L 2
where D = v τ .
Assumed that the SAR system transmits a linear frequency modulation (LFM) signal with a bandwidth of B and a wavelength of λ , the range compressed signal of the target P can be expressed as:
S r c ( t , D ) = r e c t ( D L ) sin c [ B ( t 2 R ( r p , Θ p ; D ) c ) ] exp [ j 4 π R ( r p , Θ p ; D ) λ ]
where
sin c ( x ) = sin ( π x ) π x
r e c t ( D L ) = { 1 ,    L 2 < D L 2 0 ,               else

3. Problem Description

The AFBP algorithm was proposed in [18]. Before that, all of the FBP imaging algorithms adopted an interpolation method for sub-aperture wavenumber spectrum fusion, and the echo data focusing effect was greatly affected by the accuracy of the interpolation function. The AFBP algorithm is realized under the unified polar coordinate system instead of the traditional interpolation fusion method by the wavenumber spectrum fusion. The core of the AFBP algorithm is to establish the Fourier transform relationship between the angular-domain impulse response function and the angular-domain wavenumber spectrum function. The concrete derivation process is as follows.
Since the main energy of point target P is concentrated in the small angle neighborhood with P as the center, the response function of point target P can be obtained by Equation (8). Then, integrating along the range history of R ( r p , Θ ; D ) to S r c ( t , D ) within the small angle range Θ , where Θ Θ p + Δ Θ with Θ p as the center. The impulse response function of point target P in the u t h sub-aperture can be expressed as:
I u ( r p , Θ ) = l 2 + x u l 2 + x u S r c ( t = 2 R ( r p , Θ p ; D ) c , D ) · e x p ( j 4 π 2 R ( r p , Θ ; D ) λ ) d D
where D x + x u ,     l / 2 x < l / 2 , x u represents the center of the u t h sub-aperture. Substituting Equation (5) into Equation (8), we obtain the following expression:
I u ( r p , Θ ) = l 2 + x u l 2 + x u e x p ( j 4 π Δ R ( r p , Θ ; D ) λ ) d D
where
Δ R ( r p , Θ ; D ) = R ( r p , Θ p ; D ) R ( r p , Θ ; D ) ( Θ Θ p ) · D
In Equation (10), both the quadratic term and higher-order term can be ignored for the broadside mode [18]. Substituting Equation (10) into Equation (9), the following can be obtained:
I u ( r p , Θ ) = l 2 + x u l 2 + x u e x p [ j 4 π ( Θ Θ p ) · D λ ] d D = Δ K α / 2 + K a u Δ K α / 2 + K a u e x p [ j ( K α + K a u ) ( Θ Θ p ) ] d K α = l · s i n c [ 2 l λ ( Θ Θ p ) ] · e x p [ j K a u ( Θ Θ p ) ] ,   Θ [ Δ A u 2 , Δ A u 2 ]
where Δ A u is the angular domain range of the pre-set sub-aperture SAR image, and K α and K a u represent the shape and center of the wavenumber spectrum in the angular domain, respectively.
{ K α = 4 π λ x K a u = 4 π λ x u
It can be clearly seen from Equation (11) that K α and Θ are Fourier transform pairs, and the impulse response function I u ( r p , Θ ) of the u th sub-aperture image can be abbreviated as I u ( Θ ) .
I u ( K α ) = Δ A u 2 + A c u Δ A u 2 + A c u I u ( Θ ) · e x p ( j K a Θ ) d Θ = r e c t ( K α K a u Δ K α ) · e x p ( j ( K a K a u ) Θ p ) K α ( Δ K α 2 , Δ K α 2 ] + K a u
where A c u is the initialized angular center of the unified polar coordinate system and Δ K α = 4 π / λ · l is the wavenumber width.
The AFBP algorithm is derived based on a linear trajectory with a constant speed assumption. When there are nonlinear motion errors in the trajectory of the platforms, spectra centroid displacement will be introduced for different ground targets. When fusing sub-aperture wavenumber spectra, the fusing error will degrade the SAR focusing performance if the influence of trajectory error is not considered. Moreover, the AFBP algorithm is only suitable for the broadside mode, not for the high-squint mode.
To deal with these constraints, the HS-AFBP method is proposed [19]. The linear function is used in the HS-AFBP method to fit the motion errors of the two-dimensional spatial variation, which extends the application range. For point target P ( r p , Θ p ) with HS-AFBP, the motion errors can be modeled as:
R e r r o r ( r p , Θ p ; D ) = c 0 + c 1 ( D ) · r p + c 2 ( D ) · Θ p
where c 0 is a constant, c 1 denotes the range-variant coefficient, and c 2 represents the angular-variant coefficient. The R e r r o r ( r p , Θ p ; D ) represents the difference between the ideal slant range and the actual slant range of the point target P. Then, the impulse response function of the point target under the nonlinear trajectory is derived as [19]:
I u ( r , Θ ) = l 2 + x u l 2 + x u e x p [ j ( K α + K a u ) ( Θ Θ p ) ] d x
where
{ K α = 4 π λ · ( 1 + 2 Θ p x u r p ) · x K a u = 4 π λ · ( x u + Θ p r p x u 2 c 2 ( x u ) )
The previous HS-AFBP method is only suitable for a small trajectory error and small dive angle, and its imaging performance will deteriorate greatly, even becoming invalid, for a large trajectory error and high dive angle. The specific analysis will be introduced in detail in the next section.

4. HSHD-AFBP Algorithm

4.1. Derivation of HSHD-AFBP Algorithm

According to Equation (3), the slant range difference between point target P and its adjacent angular-domain targets can be expressed as Δ R ( τ ; x p , r p ) .
Δ R ( τ ; x p , r p ) = R ( τ ; x p , r p ) R ( τ ; x p + Δ x , r p ) = ( X + Δ X x p ) 2 + ( Δ Y r p 2 H 2 x p 2 ) 2 + ( Z + Δ Z ) 2                   ( X + Δ X ( x p + Δ x ) ) 2 + ( Δ Y r p 2 H 2 ( x p + Δ x ) 2 ) 2 + ( Z + Δ Z ) 2
The second-order Taylor expansion of Equation (17) with respect to Δ x can be obtained as follows:
Δ R ( τ ; x p , r p ) = [ X Δ Y x p y p + Δ X ] · 1 r p · Δ x + 1 2 · Δ Y r p · r p 2 H 2 y p 3 · Δ x 2 + o ( Δ x 3 )
According to Equation (2), we obtain
Δ x = r p c o s γ ( Θ Θ p ) = r p c o s γ · Δ Θ
Substituting Equation (19) into Equation (18), Δ R ( τ ; x p , r p ) can be re-expressed as:
Δ R ( τ ; x p , r p ) [ X Δ Y x p y p + Δ X ] · 1 c o s γ · Δ Θ + 1 2 · Δ Y r p · r p 2 H 2 y p 3 · ( r p c o s γ ) 2 Δ Θ 2 + o ( Δ Θ 3 ) D · Δ Θ + [ Δ X Δ Y x p y p ] · 1 c o s γ · Δ Θ + 1 2 · Δ Y ( r p 2 H 2 ) y p 3 · r p ( c o s γ ) 2 Δ Θ 2 = D · Δ Θ + [ Δ X Δ Y x p y p ] · 1 c o s γ · Δ Θ + Δ r
where Δ r = 1 2 · Δ Y ( r p 2 H 2 ) y p 3 · r p ( c o s γ ) 2 Δ Θ 2 .
The term Δ r causes quadratic phase error (QPE) in BP integration, which is expressed by
φ = 4 π λ · 1 2 · Δ Y ( r p 2 H 2 ) y p 3 · r p ( c o s γ ) 2 · ( Δ Θ ) 2 = 2 π λ · Δ Y ( c o s r ) 2 · s i n α p · ( s i n θ p ) 2 · ( Δ Θ ) 2
For the sub-aperture length l u   , the angular resolution can be expressed as Δ Θ = λ / ( 2 l u ) [20]. If φ can be constrained within π / 8 , its effect on BP integration can be neglected [20]. Substituting Δ Θ = λ / ( 2 l u ) into Equation (21), we can get the constraint of sub-aperture length l u :
l u 1 c o s r · s i n θ p 4 λ · Δ Y s i n α p
For typical high maneuvering platform parameters, the constraint is easy to satisfy. Therefore, Equation (20) is rewritten as
Δ R ( r p , Θ ; D ) D · Δ Θ + [ Δ X Δ Y x p y p ] · 1 c o s γ · Δ Θ
For an ideal linear trajectory, Δ X = Δ Y = 0 . In this case, Equation (20) can be simplified as
Δ R ( r p , Θ ; D ) D · Δ Θ
Comparing Equations (23) and (24), the slant range difference between point target P and its adjacent angular-domain targets caused by motion error can be expressed as:
Δ R e r r o r ( r p , Θ p ; D ) = R e r r o r ( r p , Θ p ; D ) R e r r o r ( r p , Θ ; D ) = [ Δ X Δ Y x p y p ] · 1 c o s γ · Δ Θ
Equation (16) can be rewritten as:
{ K α = 4 π λ · ( 1 + 2 Θ p x u r p ) · x K a u = 4 π λ · ( x u + Θ p r p x u 2 + [ Δ X Δ Y x p y p ] · 1 c o s γ )
It can be seen from Equation (26) that the error component of the nonlinear motion trajectory causes the center position deviation of the target azimuth wavenumber spectrum; the displacement is related to motion error, the target position and the dive angle.
The impulse response function of the u t h sub-aperture can be expressed as:
I u ( r p , Θ ) = l u 2 + x u l u 2 + x u e x p [ j ( K α + K a u ) ( Θ Θ p ) ] d D = l u · s i n c [ 2 l u λ ( Θ Θ p ) ] · e x p [ j K a u ( Θ Θ p ) ]
Since K a u is related to the position of target point P , I u ( r p , Θ ) cannot be abbreviated as I u ( Θ ) . By executing an inverse Fourier transform of I u ( r p , Θ ) , the wavenumber spectrum function of the sub-aperture image can be obtained.
I u ( r p , K α ) = Δ A u 2 + A c u Δ A u 2 + A c u I u ( r p , Θ ) · e x p ( j K a Θ ) d Θ               = r e c t ( K α K a u Δ K α ) · e x p ( j ( K a K a u ) Θ p ) , K α ( Δ K α 2 , Δ K α 2 ] + K a u
where Δ K α = 4 π / λ · l u .
From Equations (26) and (28), it is found that the motion error along the x-axis causes a consistent spectral displacement, K C , and the motion error along the y-axis causes a spectral displacement, K V , that varies with the target’s position.
{ K C = 4 π λ · Δ X c o s γ K V = 4 π λ · x p y p · Δ Y c o s γ

4.2. Sub-Aperture Spectrum Fusion

To achieve accurate spectrum fusion, we propose an improved sub-aperture wavenumber spectrum fusion method.
Step 1: Sub-aperture spectrum shift calculation
Calculate the integral number of the central position shift of the sub-aperture wavenumber spectrum caused by motion error. The error component of the u t h sub-aperture center along the x axis is denoted as Δ X u , the error component along the y axis is denoted as Δ Y u , and the scene center point is marked as ( x 0 , y 0 , 0 ) in the Cartesian coordinate system and ( r 0 , Θ 0 ) in the polar coordinate system. The whole-aperture data after range compression are divided into sub-aperture data with the same azimuth sampling number, and each sub-aperture contains N 0 pulses. Taking the center point target of the angular domain as a reference, the integral number of the u t h sub-aperture wavenumber spectrum center position shift caused by motion error is calculated as follows:
Δ u = floor ( ( x u + Θ 0 r 0 x u 2 + [ Δ X u Δ Y u x p y p ] · 1 c o s γ ) / l u N 0 )
where floor ( · ) means the operation of taking the nearest integer toward negative infinity.
Assuming that G ( r p , Θ p ) = x p / y p , Δ G ( r p , Θ p ) = x p / y p x 0 / y 0 , the minimum oversampling rate k 0 is expressed as:
k 0 2 · Δ Y u c o s γ · l u · [ Δ G ( r p , Θ p ) ] max + 1
To improve computational efficiency, k 0 generally takes the power of 2.
Step 2: Wavenumber spectrum shifting
The wavenumber spectrum center of the u t h sub-aperture is shifted to the zero frequency, and the shift function H1 can be formulated as:
H 1 = e x p ( j 4 π λ · ( Δ u · l u N 0 ) ( Θ p Θ 0 ) )
Then the u t h sub-aperture SAR image after centroid shifting denoted as
I u ( r p , Θ ) = I u ( r p , Θ ) · H 1
Step 3: Sub-aperture spectrum fusion
Since the sub-aperture wavenumber spectrum center is shifted by the integral number, no interpolation operation is additionally needed in the sub-aperture wavenumber spectrum fusion procedure. The actual position of each sub-aperture center can be calculated according to Δ u , and then a simple sub-aperture spectrum superposition can be performed. Figure 2 and Figure 3 show the traditional and improved spectra fusion methods, respectively. In Figure 2 and Figure 3, each sub-aperture wavenumber spectrum has been moved to the zero frequency. Assuming that the oversampling factor of the sub-aperture wavenumber spectrum is k , the traditional method of spectra fusion is to simply overlap the adjacent sub-aperture wavenumber spectra by ( k 1 ) / k to obtain the whole-aperture wavenumber spectrum. When the platform moves along a uniform straight line, the correct wavenumber spectrum with whole-aperture can be obtained in this way. However, when the platform moves at a varying speed and along a curved trajectory, the sub-aperture width is different and the corresponding sub-aperture wavenumber spectrum width is also different. If the spectrum fusion is performed in the traditional way, gaps and overlaps between adjacent sub-apertures cannot be avoided. Considering the variation of spectrum width, the improved spectra fusion method effectively avoids the situation of spectrum gaps and overlaps.

4.3. Algorithm Implementation

Figure 4 shows the processing flow chart of the proposed method. The implementation process of the HSHD-AFBP algorithm is the same as that for the HS-AFBP algorithm [19], but their sub-aperture wavenumber spectrum center correction and sub-aperture spectrum fusion are different, as marked in Figure 4 with the blue color. The implementation process of the HSHD-AFBP algorithm is divided into two stages: (1) the prepare phase and (2) the process phase.
The prepare phase stage includes three steps:
(1)
Range interpolation. Interpolation operation is carried out on the range compressed data, so that the echo data can be coherently accumulated.
(2)
Sub-aperture division. The whole-aperture data after range compression are divided into Ns sub-apertures. Each sub-aperture contains the same number of pulses.
(3)
SAR imaging grid construction. The angle interval ΔΘ, established in the unified polar coordinate system, depends on the length of the sub-aperture [20].
Δ Θ λ 2 l u
The process phase stage includes four steps:
(1)
Sub-image generation. Sub-images I u ( r p , Θ ) , u = 1 , 2 , , N s , where N s is the number of sub-apertures, will be focused by the sub-apertures’ radar echoes with low resolution by using the BP integral.
(2)
Sub-aperture wavenumber spectrum center correction. Firstly, the correction of each sub-aperture wavenumber spectrum is performed by Equation (33) to obtain I u ( r p , Θ ) . Then, the angular inverse Fourier transform is applied to I u ( r p , Θ ) to obtain the wavenumber spectrum I u ( r p , K α ) of each sub-image.
(3)
Spectrum fusion. The whole-aperture wavenumber spectrum can be obtained by a simple sub-aperture spectrum superposition.
(4)
Full-resolution image generation. The angular Fourier transform of the fusion wavenumber spectrum I ( r p , K α ) is applied to obtain the full-resolution SAR image.

4.4. Computational Complexity Analysis

For simplicity, it is assumed that the size of the imaging scene is N × N and the whole aperture contains N pulses. The computational burden of the proposed algorithm mostly includes four aspects: sub-image generation ( N 3 / N s ), sub-aperture wavenumber spectrum center correction ( N 2 / N s ), sub-image angular inverse Fourier transform ( N 2 · l o g 2 N / N s ), and the whole-aperture wavenumber spectrum angular Fourier transform ( N 2 · l o g 2 N ). Hence, why the computational complexity is approximated as ( N / N s + 1 / N s + l o g 2 N / N s + l o g 2 N ) · N 2 . Comparing HSHD-AFBP and HS-AFBP, the main difference in computational complexity is in the processing steps of sub-aperture wavenumber spectrum center correction and sub-aperture spectrum fusion. The computational complexity of these two steps is much less than the complexity of sub-image generation. Therefore, HSHD-AFBP and HS-AFBP have almost the same computational complexity.

5. Experiments and Analysis

In this section, we will use the point targets and simulated nature scene to verify the effectiveness of the proposed method, as well as by comparing the focused SAR image performance of the HS-AFBP method.

5.1. Point Targets’ Simulation Results

Table 1 lists the main simulation parameters. Nine point targets are placed in the scene, as is shown in Figure 5. The coordinates of the nine point targets are shown in Table 2. The number of azimuth pulses is 2048. The whole aperture length is 446.9 m. The synthetic aperture is divided into 64 sub-apertures. Therefore, each sub-aperture contains 32 pulses. Figure 6 shows us the trajectory deviations in the corresponding coordinate-axis between the actual trajectory and the ideal trajectory, i.e., the fitting trajectory.
According to Equation (31) and the simulation parameters, we obtain k 0 = 2 if the first sub-aperture is taken as an example. Figure 7 shows the first sub-aperture two-dimensional wavenumber spectrum of three point targets in sequence. It can be seen that the spectra displacement of targets at different angles in the same slant range unit is different in the polar coordinate system. Table 3 lists the first sub-aperture wavenumber spectrum center positions of the three point targets obtained by the method of HS-AFBP, HSHD-AFBP, and the true values of the center positions of the sub-aperture wavenumber spectrum of the three point targets. HS-AFBP estimates the first sub-aperture wavenumber spectrum center positions according to Equation (16). HSHD-AFBP estimates the sub-aperture wavenumber spectrum center positions according to Equation (26). The true values of the center positions of the sub-aperture wavenumber spectrum of the three point targets can be obtained according to Figure 7. The values in Table 3 indicate the azimuth units where the center of the spectrum is folded multiple times. It can be seen from Table 3 that, whether the target is in the center of the angular-dimensional or on the edge of the angular-dimensional, HSHD-AFBP can accurately calculate wavenumber spectra displacement.
Taking the target P 4 as an example, the whole-aperture wavenumber spectrum is obtained by fusing the wavenumber spectra of 64 sub-apertures, as is shown in Figure 8. Figure 8a is the result obtained by HSHD-AFBP. Figure 8b is the result obtained by HS-AFBP. Figure 8a,b are whole-aperture wavenumber spectra. Figure 8c,d are the results of local position amplification of the whole-aperture wavenumber spectra. HS-AFBP first uses (16) to correct the center of the sub-aperture spectrum, and then uses the traditional spectra fusion method to obtain the whole-aperture wavenumber spectrum. Because the influence of trajectory error on the spectral width is not considered, there will be gaps and overlaps in the spectrum fusion result. It can be seen that there are obvious spectra gaps between adjacent sub-aperture spectra, as is shown in Figure 8d. The spectra gaps and overlaps introduce additional phase errors, which will lead to target defocusing. HSHD-AFBP uses (23) to correct the center of the sub-aperture spectrum, and then uses an accurate sub-aperture spectra fusion method to obtain the whole-aperture wavenumber spectrum. As can be seen from Figure 8c, the whole-aperture wavenumber spectrum is well fused.
Figure 9 shows the imaging results of the point targets. To prove the imaging efficiency of the proposed algorithm, the processing time of the two algorithms is measured in the same condition. They are programmed in the Matlab on a computer with a 2.80 GHz Central Processing Unit and 24 GB Random Access Memory. The processing times of HSHD-AFBP and HS-AFBP are 47.40 s and 46.26 s, respectively.
To compare the 2D focusing performance of HSHD-AFBP and HS-AFBP, the contour images of targets are shown in Figure 10. Figure 10a is the result obtained by HSHD-AFBP. Figure 10b is the result obtained by HS-AFBP. The contour image of each target in Figure 10 is centered on its own location. It can be clearly seen that, for the targets on the edge of the angular-dimensional, HS-AFBP can focus well, but HS-AFBP will cause serious defocusing. The resolution, peak side lobe radio (PSLR), and integral side lobe radio (ISLR) in the range and azimuth directions are listed in Table 4. It can be clearly seen that HSHD-AFBP has the ability to obtain well-focused SAR images for all of the nine point targets.
The simulation results reflect the following facts: (1) the trajectory errors cause different spectrum center shifts with different target locations, which will lead to SAR images’ focusing performance deterioration. These offsets must be calculated accurately to set the oversampling rate reasonably and ensure that spectral ambiguity does not occur. (2) Sub-aperture spectra fusion should consider the effect of trajectory errors on spectrum width, otherwise, spectrum gaps and overlaps will occur inevitably, which will affect the amplitude and phase of the whole-aperture spectrum, and then the SAR images will be defocused.

5.2. Nature Scene Simulation Results

To future-proof the feasibility of the proposed algorithm, a set of simulated nature scene data are processed by HSHD-AFBP and HS-AFBP. As is shown in [21], according to the point scattering model, the single look data of the SAR image can be regarded as nature scene targets’ electromagnetic backscattering coefficient in substance. Here, we use a high-resolution real airborne SAR image in X-band to simulate the high maneuvering platforms’ echo signal of the natural scene. The size of the natural scene is 820 m × 820 m (range × azimuth). The system parameters are shown in Table 5. The slant range of the scene’s center is 24.9 km. The plane flies at a speed of 1000 m/s. The pulse repetition frequency (PRF) is 4582 Hz. Figure 11 shows the raw signal of the nature scene.
Figure 12 shows the reconstructed images of the simulated nature scene by HSHD-AFBP and HS-AFBP. From Figure 12, we can see that the natural scene reconstructed by HSHD-AFBP demonstrates a better focusing effect than that of HS-AFBP. The processing times of HSHD-AFBP and HS-AFBP are 62.07 s and 60.36 s, respectively. For better comparison, we give the range and azimuth profiles of a point-like scatterer in Figure 13. The scatterer is circled in red in Figure 12. From Figure 13, one can see that, compared to those of HS-AFBP, the profiles of HSHD-AFBP have higher resolutions and lower sidelobes. All of the above results effectively prove that the proposed method can focus well on the SAR image of a high maneuvering platform.

6. Conclusions

In this paper, the HSHD-AFBP method was proposed for high squint angle and high dive angle SAR imaging, which is suitable for the nonlinear motion trajectory of high maneuvering platforms. The HSHD-AFBP can calculate the wavenumber spectra displacement of targets at different angles, while the HS-AFBP can only estimate the wavenumber spectrum displacement of the target located at the angle center. Therefore, the HSHD-AFBP can accurately calculate the oversampling rate required for sub-aperture SAR focusing and avoid sub-aperture spectrum blurring. In addition, this paper proposed an accurate sub-aperture spectra fusion method to deal with the influence of the nonlinear trajectory on spectrum width. Compared to the traditional spectrum fusion method, the proposed method improves the SAR’s focusing quality and enlarges the imaging scene’s width. Simulation experiments verified the effectiveness of the proposed method. However, a SAR image application, such as image matching and target recognition, requires the SAR image in the ground plane. Therefore, we will focus on how to quickly and accurately achieve SAR image focusing in the ground plane in future research, to extend the scope of the proposed HSHD-AFBP.

Author Contributions

Conceptualization, Y.G. and Z.S.; methodology, Y.G.; validation, Y.G., Z.S., and P.J.; formal analysis, H.L.; investigation, H.L.; resources, P.J.; data curation, Z.S.; writing—original draft preparation, Y.G.; writing—review and editing, Z.S.; visualization, H.L.; supervision, P.J.; funding acquisition, Z.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China under grant number 61671355.

Institutional Review Board Statement

Not applicable as this study did not involve human or animal subjects.

Informed Consent Statement

Not applicable as this study did not involve human or animal subjects.

Data Availability Statement

The data presented in this study are available on request from the corresponding authors.

Acknowledgments

We thank the National Natural Science Foundation of China. And we also thank the anonymous reviewers for their valuable comments to improve the paper’s quality.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cumming, I.G.; Wong, F.H. Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation; Artech House: Norwood, MA, USA, 2005. [Google Scholar]
  2. Soumekh, M. Synthetic Aperture Radar Signal Processing with MATLAB Algorithms; John Wiley: New York, NY, USA, 1999. [Google Scholar]
  3. Zeng, T.; Li, Y.; Ding, Z.; Long, T.; Yao, D.; Sun, Y. Subaperture Approach Based on Azimuth-Dependent Range Cell Migration Correction and Azimuth Focusing Parameter Equalization for Maneuvering High-Squint-Mode SAR. IEEE Trans. Geosci. Remote Sens. 2015, 53, 6718–6734. [Google Scholar] [CrossRef]
  4. Li, Z.; Xing, M.; Liang, Y.; Gao, Y.; Chen, J.; Huai, Y.; Zeng, L.; Sun, G.-C.; Bao, Z. A Frequency-Domain Imaging Algorithm for Highly Squinted SAR Mounted on Maneuvering Platforms With Nonlinear Trajectory. IEEE Trans. Geosci. Remote Sens. 2016, 54, 4023–4038. [Google Scholar] [CrossRef]
  5. Bie, B.; Sun, G.-C.; Xia, X.-G.; Xing, M.; Guo, L.; Bao, Z. High-Speed Maneuvering Platforms Squint Beam-Steering SAR Imaging Without Subaperture. IEEE Trans. Geosci. Remote Sens. 2019, 57, 6974–6985. [Google Scholar] [CrossRef]
  6. Mei, H.; Meng, Z.; Liu, M.; Li, Y.; Quan, Y.; Zhu, S.; Xing, M. Thorough Understanding Property of Bistatic Forward-Looking High-Speed Maneuvering Platform SAR. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 1826–1845. [Google Scholar] [CrossRef]
  7. Smith, A.M. A new approach to range-Doppler SAR processing. Int. J. Remote Sens. 1991, 12, 235–251. [Google Scholar] [CrossRef]
  8. Raney, R.; Runge, H.; Bamler, R.; Cumming, I.; Wong, F. Precision SAR processing using chirp scaling. IEEE Trans. Geosci. Remote Sens. 1994, 32, 786–799. [Google Scholar] [CrossRef]
  9. Davidson, G.; Cumming, I.; Ito, M. A chirp scaling approach for processing squint mode SAR data. IEEE Trans. Aerosp. Electron. Syst. 1996, 32, 121–133. [Google Scholar] [CrossRef]
  10. Moreira, A.; Huang, Y. Airborne SAR processing of highly squinted data using a chirp scaling approach with integrated motion compensation. IEEE Trans. Geosci. Remote Sens. 1994, 32, 1029–1040. [Google Scholar] [CrossRef]
  11. Wong, F.; Yeo, T. New applications of nonlinear chirp scaling in SAR data processing. IEEE Trans. Geosci. Remote Sens. 2001, 39, 946–953. [Google Scholar] [CrossRef]
  12. An, D.; Huang, X.; Jin, T.; Zhou, Z. Extended Nonlinear Chirp Scaling Algorithm for High-Resolution Highly Squint SAR Data Focusing. IEEE Trans. Geosci. Remote Sens. 2012, 50, 3595–3609. [Google Scholar] [CrossRef]
  13. Sun, G.; Xing, M.; Liu, Y.; Sun, L.; Bao, Z.; Wu, Y. Extended NCS Based on Method of Series Reversion for Imaging of Highly Squinted SAR. IEEE Geosci. Remote Sens. Lett. 2010, 8, 446–450. [Google Scholar] [CrossRef]
  14. Vu, V.T.; Sjogren, T.K.; Petersson, M.I. A comparison between fast farctorized backprojection and frequency-domain algo-rithms in UWB low frequency SAR. In Proceedings of the IEEE International Geoscience and Remote Sensing, Boston, MA, USA, 7–11 July 2008; pp. 1293–1296. [Google Scholar]
  15. Yegulalp, A. Fast backprojection algorithm for synthetic aperture radar. In Proceedings of the IEEE Radar Conference, Waltham, MA, USA, 22–22 April 1999; pp. 60–65. [Google Scholar]
  16. Wahl, D.E.; Yocky, D.A.; Jakowatz, J.C.V. An implementation of a fast backprojection image formation algorithm for spotlight-mode SAR. In Proceedings of the SPIE—Algorithms for Synthetic Aperture Radar Imagery XV, Orlando, FL, USA, 17–18 March 2008. [Google Scholar] [CrossRef]
  17. Ulander, L.; 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]
  18. Zhang, L.; Li, H.L.; Qiao, Z.J.; Xu, Z.W. A fast BP algorithm with wavenumber spectrum fusion for high resolution spotlight SAR imagery. IEEE Geosci. Remote Sens. Lett. 2014, 11, 1460–1464. [Google Scholar] [CrossRef]
  19. Ran, L.; Liu, Z.; Xie, R.; Zhang, L. Focusing High-Squint Synthetic Aperture Radar Data Based on Factorized Back-Projection and Precise Spectrum Fusion. Remote Sens. 2019, 11, 2885. [Google Scholar] [CrossRef] [Green Version]
  20. Carrara, W.G.; Goodman, R.S.; Majewski, R.M. Spotlight Synthetic Aperture Radar:Signal Processing Algorithm; Artech House: Boston, MA, USA, 1995; pp. 245–254. [Google Scholar]
  21. Zhang, S.S.; Long, T.; Zeng, T.; Ding, Z.G. Space-borne synthetic aperture radar received data simulation based on air-borne SAR image data. Adv. Space Res. 2008, 41, 1818–1821. [Google Scholar] [CrossRef]
Figure 1. SAR imaging geometric model of high maneuvering platforms; (a) Cartesian coordinate system, (b) polar coordinate system.
Figure 1. SAR imaging geometric model of high maneuvering platforms; (a) Cartesian coordinate system, (b) polar coordinate system.
Remotesensing 13 01649 g001
Figure 2. Schematic diagram of traditional spectrum fusion method.
Figure 2. Schematic diagram of traditional spectrum fusion method.
Remotesensing 13 01649 g002
Figure 3. Schematic diagram of improved spectrum fusion method.
Figure 3. Schematic diagram of improved spectrum fusion method.
Remotesensing 13 01649 g003
Figure 4. Flowchart of HSHD-AFBP.
Figure 4. Flowchart of HSHD-AFBP.
Remotesensing 13 01649 g004
Figure 5. Point targets’ positions.
Figure 5. Point targets’ positions.
Remotesensing 13 01649 g005
Figure 6. Trajectory deviations.
Figure 6. Trajectory deviations.
Remotesensing 13 01649 g006
Figure 7. The first sub-aperture two-dimensional wavenumber spectrum of three point targets; (a)   P 4 , (b)   P 5 , (c)   P 6 .
Figure 7. The first sub-aperture two-dimensional wavenumber spectrum of three point targets; (a)   P 4 , (b)   P 5 , (c)   P 6 .
Remotesensing 13 01649 g007
Figure 8. Spectrum fusion results of P 4 ; (a) the whole-aperture wavenumber spectrum in HSHD-AFBP, (b) the whole-aperture wavenumber spectrum in HS-AFBP, (c) the result of local position amplification in (a), (d) the result of local position amplification in (b).
Figure 8. Spectrum fusion results of P 4 ; (a) the whole-aperture wavenumber spectrum in HSHD-AFBP, (b) the whole-aperture wavenumber spectrum in HS-AFBP, (c) the result of local position amplification in (a), (d) the result of local position amplification in (b).
Remotesensing 13 01649 g008
Figure 9. Image results obtained by using (a) HSHD-AFBP and (b) HS-AFBP.
Figure 9. Image results obtained by using (a) HSHD-AFBP and (b) HS-AFBP.
Remotesensing 13 01649 g009
Figure 10. The contour images of the targets P 1 , P 5 , and P 9 , in turn; (a) HSHD-AFBP, (b) HS-AFBP.
Figure 10. The contour images of the targets P 1 , P 5 , and P 9 , in turn; (a) HSHD-AFBP, (b) HS-AFBP.
Remotesensing 13 01649 g010
Figure 11. Generated raw signals of the simulated nature scene; (a) amplitude, (b) phase.
Figure 11. Generated raw signals of the simulated nature scene; (a) amplitude, (b) phase.
Remotesensing 13 01649 g011
Figure 12. Image results obtained using (a) HSHD-AFBP and (b) HS-AFBP.
Figure 12. Image results obtained using (a) HSHD-AFBP and (b) HS-AFBP.
Remotesensing 13 01649 g012
Figure 13. Range and azimuth profile comparison of HSHD-AFBP and HS-AFBP; (a) range profiles, (b) azimuth profiles.
Figure 13. Range and azimuth profile comparison of HSHD-AFBP and HS-AFBP; (a) range profiles, (b) azimuth profiles.
Remotesensing 13 01649 g013
Table 1. Main system parameters.
Table 1. Main system parameters.
ParameterValue
Wave bandKu
Bandwidth100 M
Center line of scene20.89 km
Squint angle 70
Dive angle 60
Whole-aperture length446.9 m
Table 2. Point targets’ coordinates.
Table 2. Point targets’ coordinates.
Target(x,y,z)/km (r/km, Θ)
P 1 (10.22, 7.90, 0)(20.56, 0.9224)
P 2 (10.93, 6.88, 0)(20.56, 0.9397)
P 3 (11.64, 5.60, 0)(20.56, 0.9570)
P 4 (10.82, 7.95, 0)(20.89, 0.9224)
P 5 (11.54, 6.86, 0)(20.89, 0.9397)
P 6 (12.26, 5.47, 0)(20.89, 0.9570)
P 7 (11.42, 7.97, 0)(21.22, 0.9224)
P 8 (12.16, 6.80, 0)(21.22, 0.9397)
P 9 (12.89, 5.28, 0)(21.22, 0.9570)
Table 3. The spectrum center azimuth unit.
Table 3. The spectrum center azimuth unit.
MethodSpectrum Center Azimuth Unit
P 4 P 5   P 6  
HS-AFBP0.66 10.660.66
HSHD-AFBP60.46 63.855.78
Truth value60.440.185.28
1 The values indicate the azimuth units where the center of the spectrum is folded multiple times.
Table 4. Focusing performance evaluation of the nine point targets when using HSHD-AFBP.
Table 4. Focusing performance evaluation of the nine point targets when using HSHD-AFBP.
TargetResolution (m)PSLR (dB)ISLR (dB)
RangeAzimuthRangeAzimuthRangeAzimuth
P 1 1.320.98−13.32−13.19−11.07−9.99
P 2 1.321.10−13.29−13,23−11.05−10.12
P 3 1.321.35−13.27−13.30−11.00−10.32
P 4 1.320.99−13.29−13.16−11.04−9.76
P 5 1.321.12−12.89−13.16−11.05−9.73
P 6 1.321.32−13.29−13.21−11.03−9.77
P 7 1.321.03−13.31−13.24−11.04−10.21
P 8 1.321.18−13.28−13.27−11.02−10.22
P 9 1.321.39−13.30−13.28−11.03−10.20
Table 5. Main system parameters.
Table 5. Main system parameters.
ParameterValue
Wave bandKu
Bandwidth100 M
Center line of scene24.9 km
Squint angle 70
Dive angle 25
Whole-aperture length446.9 m
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, Y.; Suo, Z.; Jiang, P.; Li, H. A Fast Back-Projection SAR Imaging Algorithm Based on Wavenumber Spectrum Fusion for High Maneuvering Platforms. Remote Sens. 2021, 13, 1649. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13091649

AMA Style

Guo Y, Suo Z, Jiang P, Li H. A Fast Back-Projection SAR Imaging Algorithm Based on Wavenumber Spectrum Fusion for High Maneuvering Platforms. Remote Sensing. 2021; 13(9):1649. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13091649

Chicago/Turabian Style

Guo, Yuan, Zhiyong Suo, Penghui Jiang, and Han Li. 2021. "A Fast Back-Projection SAR Imaging Algorithm Based on Wavenumber Spectrum Fusion for High Maneuvering Platforms" Remote Sensing 13, no. 9: 1649. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13091649

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