Investigation of Wavenumber Domain Imaging Algorithm for Ground-Based Arc Array SAR

Ground-based synthetic aperture radar (GB-SAR) has become an important technique for remote sensing deformation monitoring. However, most of the existing GB-SAR systems realize synthetic aperture by exploiting two closely spaced horn antennas to move along a linear rail. In order to obtain higher data acquisition efficiency and a wider view angle, we introduce arc antenna array technology into the GB-SAR system, which realizes a novel kind of system: ground-based arc array SAR (GB-AA-SAR). In this paper, we analyze arc observation geometry and derive analytic expressions of sampling criteria. Then, we propose a novel wavenumber domain imaging algorithm for GB-AA-SAR, which can achieve high image reconstruction precision through numerical solutions in the wavenumber domain. The proposed algorithm can be applied in wide azimuth view angle scenarios, and the problem of azimuth mismatch caused by distance approximation in arc geometric efficient omega-k imaging can be solved successfully. Finally, we analyze the two-dimensional (2D) spatial resolution of GB-AA-SAR, and verify the effectiveness of the proposed algorithm through numerical simulation experiments.


Introduction
Landslides, mudslides, and other deformation disasters seriously endanger people's personal safety, as well as the safety of property. As a new type of technology for microwave remote sensing deformation measurement, ground-based synthetic aperture radar (GB-SAR) has been widely regarded for its extensive application prospects in deformation monitoring. Many research institutions have successfully developed GB-SAR systems, which also have been applied to deformation monitoring and disaster prevention [1][2][3][4][5][6][7]. The GB-SAR interferometry exploits non-contact microwave remote sensing measurement, which can achieve millimeter or even sub-millimeter level deformation monitoring accuracy [1,2]. The conventional GB-SAR system obtains high-range resolution by controlling antennas to transmit and receive wideband step frequency continue wave (SFCW) signals in the range direction, and the radar moves along a linear or arc mechanical rail to realize synthetic aperture. Then, the two-dimensional (2D) resolving of the observation area can be achieved. For catering to the development trends of miniaturization, lightness, and portability, the frequency modulation continuous wave (FMCW) signal has been applied into the design of the GB-SAR system. In order to further upgrade and optimize the existing GB-SAR system, this paper applies the arc antenna array technology to the GB-SAR system, which can obtain higher data acquisition efficiency and a wider view angle than conventional GB-SAR.
Arc antenna array technology is the technology that combines the arc observation geometry with an antenna array configuration. In [8], Klausing Helmut and Keydel Wolfgang firstly stated that a synthetic aperture can be obtained via the rotational movement, then subsequently proposed the corresponding imaging algorithms. Arc GB-SAR has attracted a great amount of interest in recent years, and some institutions have already developed their own prototype systems [9][10][11]. However, most of the current arc GB-SAR systems exploit the movement of a single pair of antennas rather than utilize an antenna array. The Joint Research Centre (JRC) has developed an array antenna system called MELISSA, but the observation geometry is equivalent to a liner array [12,13]. In this paper, an arc antenna array configuration is applied to the design of the GB-SAR system, which is referred to as ground-based arc array SAR (GB-AA-SAR). High-range resolution can be achieved for a FMCW signal with a large time-width product, whereas in the azimuth direction, the switched array antenna channels scan quickly over an arc aperture, and finally, high-speed and wide view angle data acquisition is realized.
This paper is organized as follows. In Section 2, the configuration of the GB-AA-SAR system and arc observation geometry are firstly introduced. Then, the corresponding arc array sampling criteria is derived, and the comparison between GB-AA-SAR and ground-based linear array SAR (GB-LA-SAR) is given. In Section 3, a novel wavenumber domain imaging algorithm for GB-AA-SAR is proposed, and the attainable spatial resolutions are analyzed. Finally, in Section 4, the GB-AA-SAR concept and the proposed imaging algorithm are validated by numerical simulations, and this paper is concluded in Section 5. Figure 1a shows the designed arc antenna array, where "AC" represents the transmitting antennas, and "BD" represents the receiving antennas. The transmitting and receiving antennas distribute uniformly along an arc, where the angle interval between adjacent antenna elements is fixed to be ∆θ Interval , and the radius of the arc is R arc . It is worth noting that the horn antenna element is adopted in this paper; on one hand, the horn antenna usually has a higher gain than that of other antenna forms, on the other hand, antennas such as microstrips cannot always support a wide bandwidth as much as the one used in this paper, which has reached up to 1 GHz. According to the equivalent phase center principle, a pair consisting of a transmitting antenna and a receiving antenna working in the quasi-stationary mode can be equivalent to a co-located virtual phase center. Thus, a uniform sampling of the arc array can be achieved. The transmitting and receiving antennas can be arranged in a staggered manner, which can realize single-transmitting and double-receiving data acquisition, and bring the benefits of improving data acquisition efficiency and the azimuth sampling rate [14]. As shown in Figure 1b, our GB-AA-SAR system is composed of three sub-blocks: an arc antenna array, a microwave switch network, and an FMCW transmitting/receiving (T/R) module. A microwave switch network can sequentially choose the antenna elements of arc array, and the FMCW T/R module also needs to transmit and receive the electromagnetic signal synchronously to finish the data acquisition. Compared with the mechanical movement of conventional GB-SAR systems, the electronic scanning mode of arc antenna array has outstanding advantages, such as for example, efficient data acquisition, a wide view angle, and a stable azimuth-sampling position in azimuth. Therefore, the GB-AA-SAR has great application potential in remote sensing deformation monitoring. GB-AA-SAR is usually applied in the deformation monitoring of slopes, dams, bridges, etc., which often lie in the far field of the antenna. The observation geometry of the GB-AA-SAR system is shown in Figure 2, where O is the center of the arc array, θn is the angle between the nth antenna element and coordinate axis OP, R0 is the distance from target to center O, Rn is the distance from target to the nth antenna element, and θt is the angle between the target and coordinate axis OP.

Arc Antenna Array and System Configuration
According to the geometric model, the distance from a point target located at (R0, θt) to the nth antenna element can be expressed as: (1) For the FMCW GB-AA-SAR system, the transmit signal is: GB-AA-SAR is usually applied in the deformation monitoring of slopes, dams, bridges, etc., which often lie in the far field of the antenna. The observation geometry of the GB-AA-SAR system is shown in Figure 2, where O is the center of the arc array, θ n is the angle between the nth antenna element and coordinate axis OP, R 0 is the distance from target to center O, R n is the distance from target to the nth antenna element, and θ t is the angle between the target and coordinate axis OP. According to the geometric model, the distance from a point target located at (R 0 , θ t ) to the nth antenna element can be expressed as: GB-AA-SAR is usually applied in the deformation monitoring of slopes, dams, bridges, etc., which often lie in the far field of the antenna. The observation geometry of the GB-AA-SAR system is shown in Figure 2, where O is the center of the arc array, θn is the angle between the nth antenna element and coordinate axis OP, R0 is the distance from target to center O, Rn is the distance from target to the nth antenna element, and θt is the angle between the target and coordinate axis OP.
According to the geometric model, the distance from a point target located at (R0, θt) to the nth antenna element can be expressed as: (1) For the FMCW GB-AA-SAR system, the transmit signal is: For the FMCW GB-AA-SAR system, the transmit signal is: where rect · T represents the rectangular envelope with duration T, t is the fast time, T is the pulse duration, f c is the center frequency, K is the frequency modulation slope, and B = KT is the signal bandwidth. If the reflectivity of the point target is σ, the echo signal can be expressed as: where p r (t) is the point spread function in the range direction. Then, the RVP can be compensated by multiplying it with the following spatial-varying compensation terms [14]: The spatial domain signal after compensation is:

Sampling Criteria of Ground-Based Arc Array SAR
In order to avoid aliasing during the imaging, the GB-AA-SAR system must obey the Nyquist sampling criterion. For fast-time domain sampling, the frequency domain support band of the dechirped echo signal is determined by: where k r = 2π f c is the range wavenumber, k rmin = min{k r } is the minimum value of the range wavenumber, k rmax = max{k r } is the maximum value of the range wavenumber. In the range direction, when Fourier analysis is performed, the range sampling rate needs to satisfy the following condition [14]: where ∆R is the range resolution cell. For the angle sampling interval in the azimuth direction, the instantaneous angular frequency can be expressed as: where θ d is the maximum angle difference between all of the antenna elements and the zero Doppler center of the observed target. Then, the angle-sampling interval in azimuth can be derived as follows: where k θmax is the maximum value of k θ , and k θmin is the minimum value of k θ . The azimuth angle-sampling interval in Equation (11) depends on the system frequency, the target distance, and the beamwidth of the antenna element. However, the arrangement of the arc antenna array requires a specific sampling interval, so we will combine antenna beamwidth with arc observation geometry for further analysis. Figure 3 shows the geometric model of GB-AA-SAR for targets at different distances, which takes into account the effect of antenna beamwidth on the observed target. The antenna beamwidth is θ s . The farther the target is, the more antenna elements will illuminate the target, and vice versa. Firstly, taking the far target T2 as an example, according to Equation (11), the angle-sampling interval can be expressed as follows: The following relationship can be obtained, according to the observation geometry in Figure 3.
where R  is the range resolution cell.
For the angle sampling interval in the azimuth direction, the instantaneous angular frequency can be expressed as: where  d is the maximum angle difference between all of the antenna elements and the zero Doppler center of the observed target. Then, the angle-sampling interval in azimuth can be derived as follows: where ' max k  is the maximum value of k '  , and ' min k  is the minimum value of k '  . The azimuth angle-sampling interval in Equation (11) depends on the system frequency, the target distance, and the beamwidth of the antenna element. However, the arrangement of the arc antenna array requires a specific sampling interval, so we will combine antenna beamwidth with arc observation geometry for further analysis. Figure 3 shows the geometric model of GB-AA-SAR for targets at different distances, which takes into account the effect of antenna beamwidth on the observed target. The antenna beamwidth is s  . The farther the target is, the more antenna elements will illuminate the target, and vice versa.
Firstly, taking the far target T2 as an example, according to Equation (11), the angle-sampling interval can be expressed as follows: The following relationship can be obtained, according to the observation geometry in Figure 3.   Substituting Equation (13) into Equation (12), the angle-sampling interval is available: The same conclusion can be obtained for the near target T1, so the azimuth angle-sampling interval of the arc array is independent of the target distance. Further to analysis, taking the maximum frequency f max into Equation (14), the sampling interval in the azimuth direction can be expressed as: where λ min = c/ f max is the minimum wavelength. It can be seen that the azimuth angle-sampling interval of the arc array is only related to shortest wavelength, arc radius, and antenna beamwidth, which are known in the system design process. Therefore, the arrangement of the antenna elements is available.

Comparison with GB-LA-SAR
In order to illustrate the advantages of GB-AA-SAR compared with GB-LA-SAR, a comparison of both configurations are demonstrated in Table 1. Single-transmitting and double-receiving data acquisition are performed in both configurations. It can be seen when the rest of the parameters are given and similar in GB-AA-SAR and GB-LA-SAR, the view angle range of GB-AA-SAR is wider than that of GB-LA-SAR. Since the view angle range of GB-LA-SAR is limited by the antenna beamwidth, the view angle range of GB-AA-SAR is dependent on the length of the arc, which can be increased easily, even to the length of the circle. So, the GB-AA-SAR has the capability to achieve a wide view angle or even 360 • deformation monitoring.

Wavenumber Domain Imaging Algorithm for GB-AA-SAR
The wavenumber domain imaging is performed in a two-dimensional frequency domain. The spatial signal can be transformed into the frequency domain by fast Fourier transform (FFT), and the typical imaging algorithm is the omega-k algorithm [15,16]. However, for the arc observation geometry, the data acquisition of the SAR signal is realized in polar coordinate format. So, the two-dimensional FFT is not applicable, and the conventional omega-k imaging cannot be directly applied to arc observation geometry. Of course, the spatial interpolation can realize the transformation from polar coordinates to Cartesian coordinates. However, the interpolation requires some computation. More importantly, a large number of zero-padding operations in azimuth are required, which will bring a great burden on the memory and calculation [15]. Therefore, this scheme is not suitable for the GB-AA-SAR imaging with a wide view angle. Dallinger has proposed an efficient omega-k algorithm in polar coordinates [16]. However, due to the Taylor approximation in the efficient omega-k algorithm, the mismatch problem in azimuth is brought about, which is unacceptable for GB-AA-SAR wide view angle imaging. The parameter k y in Dallinger's algorithm is related to the target distance. For airborne and spaceborne SAR imaging, since the observation scene is much smaller than the distance from the observation scene to radar, the reference distance can take the center of the observation scene. However, the approximation is not precise in a ground-based platform, so the algorithm is not suitable for GB-AA-SAR imaging.
For the signal S RVP (t, θ n ) shown in Equation (7), range Fourier transform is applied, and the frequency domain signal can be expressed as follows: where f = f c + Kt is the signal frequency. After Fourier transform in azimuth, the signal in the two-dimensional spatial frequency domain could be described as follows: where k θ = 2π f η w s is the spatial angular frequency in azimuth direction, the phase in Equation (17) is: In order to convert the spatial variant properties into spatial invariant properties, and realize the efficient omega-k imaging in polar coordinates for arc observation geometry, it is necessary to rewrite Equation (18) as follows: where P(k, k θ ) represents a function of k, k θ , Q(R 0 ) represents a function of R 0 , H(k θ ) represents a function of k θ , and G(θ c ) represents a function of θ c . Therefore, the stationary phase point θ pos n in Equation (19) could be solved through the Taylor expansion approximation and principle of stationary phase (POSP): where α = R 0 − R arc and β = R 0 R arc , taking Equation (20) into Equation (17), the two-dimensional spectrum in polar coordinates can be rewritten as follows: Then, a nonlinear mapping referred to as the STOLT interpolation [15] is performed on the two-dimensional spectrum, where the mapping relationship is: Finally, the efficient omega-k imaging is finished by two-dimensional inverse fast Fourier transform (IFFT).
Due to the efficient omega-k algorithm not being suitable for GB-AA-SAR, this paper presents a novel wavenumber domain imaging algorithm for GB-AA-SAR. The focus process in the range direction is realized by numerical solution in the wavenumber domain, and then, azimuth FFT without zero padding is applied to obtain the accurate SAR image, specifically: Firstly, let v = θ n − θ t , then Equation (17) can be rewritten as follows: The phase in Equation (23) is: Then the stationary phase can be obtained by solving Equation (24): The result is derived as follows: where α 0 = R 0 2 + R arc 2 , taking the results in Equation (26) into Equation (23), the result could be computed as follows: For the target in distance R i , let the compensation phase be: Then, the range accumulation is applied in the wavenumber domain: For the observation scene with N r sampling points, the range profile can be obtained by performing the operation to all of the sampling distances in the observation scene, and then the accurate SAR image is obtained through FFT operation in the azimuth direction.
In summary, the wavenumber domain imaging procedure of GB-AA-SAR is illustrated in Figure 4, which can be split into the following steps: Step 1 Perform an IFFT on the measured echo data S IF (t, θ n ) with respect to the fast time t in the range direction.
Step 2 Compensate the RVP using phase factor H, and obtain the result S RVP (t, θ n ).
Step 3 Perform the 2D FFT of the S RVP (t, θ n ) with respect to t and θ n to obtain the 2D frequency domain signal S kr, kθ . Step 4 For R = R i , realize a numerical solution in the wavenumber domain to acquire the stationary phase Ψ i .
Step 5 Compensate the phase Ψ i ; then, perform range accumulation in order to accomplish the focus procedure at R = R i . Step 6 Repeat Steps 4 to 5 for all of the available R i , which can be rapidly implemented by parallel processing, then perform azimuth IFFT to the accumulated result S R, kθ to accomplish the focus procedure in the azimuth direction.
Step 7 Finally, display the SAR image. range direction.
Step 2 Compensate the RVP using phase factor H, and obtain the result SRVP(t, θn).

Step 3
Perform the 2D FFT of the SRVP(t, θn) with respect to t and θn to obtain the 2D frequency domain signal Skr, kθ.

Step 4
For R = Ri, realize a numerical solution in the wavenumber domain to acquire the stationary phase Ѱi.
Step 5 Compensate the phase Ѱi; then, perform range accumulation in order to accomplish the focus procedure at R = Ri.
Step 6 Repeat Steps 4 to 5 for all of the available Ri, which can be rapidly implemented by parallel processing, then perform azimuth IFFT to the accumulated result SR, kθ to accomplish the focus procedure in the azimuth direction.
Step 7 Finally, display the SAR image.

Range Resolution and Azimuth Angular Resolution
For GB-AA-SAR, the point-spread function (PSF) of a point target is: where sin c(x) = sin(πx)/πx, B kr is the wavenumber bandwidth in the range direction, and B kθ is the wavenumber bandwidth in the azimuth direction. The resolution is usually defined as −3 dB bandwidth of the PSF main lobe. Then, according to the discussion in Section 2.2, the range resolution ρ r and the angular resolution ρ θ can be expressed as follows: It can be seen that, if the system parameters are given, the azimuth angle resolution of GB-AA-SAR is a constant, which is independent of the target distance. Figure 5 shows the azimuth angle resolution curve of GB-AA-SAR for different antenna beamwidths and arc radii.
where       c x = x x sin sin / , kr B is the wavenumber bandwidth in the range direction, and k B  is the wavenumber bandwidth in the azimuth direction. The resolution is usually defined as −3 dB bandwidth of the PSF main lobe. Then, according to the discussion in Section 2.2, the range resolution  r and the angular resolution   can be expressed as follows: It can be seen that, if the system parameters are given, the azimuth angle resolution of GB-AA-SAR is a constant, which is independent of the target distance. Figure 5 shows the azimuth angle resolution curve of GB-AA-SAR for different antenna beamwidths and arc radii.

Numerical Simulation Experiments
In this section, numerical simulation experiments for GB-AA-SAR are performed. Firstly, point target imaging simulation is performed to verify the effectiveness of the proposed algorithm through comparing it with the back projection algorithm (BP) and efficient omega-k algorithm. Then, multipoint target imaging is carried out with different distances and azimuth angles, and the 2D resolution of GB-AA-SAR is analyzed. The spatial invariant property of the azimuth angle resolution and the correctness of its theoretical value are verified. Simulation parameters are shown in Table 2.

Numerical Simulation Experiments
In this section, numerical simulation experiments for GB-AA-SAR are performed. Firstly, point target imaging simulation is performed to verify the effectiveness of the proposed algorithm through comparing it with the back projection algorithm (BP) and efficient omega-k algorithm. Then, multi-point target imaging is carried out with different distances and azimuth angles, and the 2D resolution of GB-AA-SAR is analyzed. The spatial invariant property of the azimuth angle resolution and the correctness of its theoretical value are verified. Simulation parameters are shown in Table 2. The point target imaging simulation using the proposed wavenumber imaging algorithm is carried out. The imaging results are compared with the results of the BP algorithm and the efficient omega-k algorithm in order to verify the imaging performance of the proposed algorithm. Figure 6 shows the imaging results of the three algorithms. Figure 6a is the point target imaging result using the proposed wavenumber domain imaging algorithm, Figure 6b is the point target imaging result of the BP algorithm, and Figure 6c gives the result of the efficient omega-k algorithm. Correspondingly, Figure 6d-f show the azimuth profile of the three imaging algorithms. It can be seen that all of the algorithms can finish the point target focus. In order to quantitatively compare the imaging performance of the algorithms, we measure the image quality parameters: impulse response width (IRW), peak sidelobe ratio (PSLR), and integrated sidelobe ratio (ISLR). IRW is the main lobe width of impulse response value above −3 dB, which is also known as resolution [17]. The point target imaging simulation using the proposed wavenumber imaging algorithm is carried out. The imaging results are compared with the results of the BP algorithm and the efficient omega-k algorithm in order to verify the imaging performance of the proposed algorithm. Figure 6 shows the imaging results of the three algorithms. Figure 6a is the point target imaging result using the proposed wavenumber domain imaging algorithm, Figure 6b is the point target imaging result of the BP algorithm, and Figure 6c gives the result of the efficient omega-k algorithm. Correspondingly, Figure 6d-f show the azimuth profile of the three imaging algorithms. It can be seen that all of the algorithms can finish the point target focus. In order to quantitatively compare the imaging performance of the algorithms, we measure the image quality parameters: impulse response width (IRW), peak sidelobe ratio (PSLR), and integrated sidelobe ratio (ISLR). IRW is the main lobe width of impulse response value above −3 dB, which is also known as resolution [17].  Table 3. It can be seen that the performance of the proposed algorithm is similar to the BP algorithm, but the proposed algorithm could be more efficient than the BP imaging algorithm. The proposed algorithm is inefficient compared with the omega-k algorithm, but can solve the problem of azimuth mismatch. (c,f) omega-k algorithm.
All of the imaging processes are performed on the same computer with the same imaging range, which is [570 m, 630 m]. The analysis results are shown in Table 3. It can be seen that the performance of the proposed algorithm is similar to the BP algorithm, but the proposed algorithm could be more efficient than the BP imaging algorithm. The proposed algorithm is inefficient compared with the omega-k algorithm, but can solve the problem of azimuth mismatch. In Section 3.2, we analyzed the 2D resolution of the GB-AA-SAR system, and derived the analytic expression of the range resolution and the azimuth resolution. The conclusion showed that the azimuth angle resolution is independent of the target distance and only depends on the wavelength, the radius of the arc array, and the antenna beamwidth. So, the spatial resolution of GB-AA-SAR in polar coordinates has spatially invariant properties. For further study, an analysis of the targets with different observation distances and different azimuth angles should be performed, and two-dimensional IRW given to verify the above conclusions.
In order to verify the 2D resolution sufficiently, point targets at (600 m, 0 • ), (10 m, 0 • ), (600 m, 30 • ), and (600 m, 45 • ) are analyzed respectively. The distribution of the targets is shown in Figure 7. As shown in Figure 8a,b, it can be seen that the resolutions are the same for the targets with different distances. Figure 8a,c shows the focus results of the targets with the same observation distance, but different azimuth angles. The resolutions are consistent with the conclusion and the theoretical values previously analyzed. However, in Figure 8d, when the azimuth angle of the target expanded to 45 • , the azimuth resolution decreased. The reason is that the angle range of the arc array is designed as [−60 • , 60 • ], when the target is at 45 • in azimuth, there are not enough array antennas to illuminate the target [15,18]. Correspondingly, the azimuth resolution decreases with the reducing of the azimuth synthesis aperture, which is not contradictory to the conclusion that the azimuth resolution has spatially invariant properties. This phenomenon also occurs in the linear array case; thus, it is not a new phenomenon. In practice, if it is necessary to carry out regional observation with a super-wide view angle, and the cost is allowed, the optimal solution is to increase the arc length and the number of antenna elements. The sub-optimal solution is to realize data acquisition by the mechanical rotating of the arc array, which is low-cost but inefficient. Both solutions can realize high-resolution and super-wide view angle imaging processing.
In Section 3.2, we analyzed the 2D resolution of the GB-AA-SAR system, and derived the analytic expression of the range resolution and the azimuth resolution. The conclusion showed that the azimuth angle resolution is independent of the target distance and only depends on the wavelength, the radius of the arc array, and the antenna beamwidth. So, the spatial resolution of GB-AA-SAR in polar coordinates has spatially invariant properties. For further study, an analysis of the targets with different observation distances and different azimuth angles should be performed, and twodimensional IRW given to verify the above conclusions.
In order to verify the 2D resolution sufficiently, point targets at (600 m, 0°), (10 m, 0°), (600 m, 30°), and (600 m, 45°) are analyzed respectively. The distribution of the targets is shown in Figure 7. As shown in Figure 8a,b, it can be seen that the resolutions are the same for the targets with different distances. Figure 8a,c shows the focus results of the targets with the same observation distance, but different azimuth angles. The resolutions are consistent with the conclusion and the theoretical values previously analyzed. However, in Figure 8d, when the azimuth angle of the target expanded to 45°, the azimuth resolution decreased. The reason is that the angle range of the arc array is designed as [−60°, 60°], when the target is at 45° in azimuth, there are not enough array antennas to illuminate the target [15,18]. Correspondingly, the azimuth resolution decreases with the reducing of the azimuth synthesis aperture, which is not contradictory to the conclusion that the azimuth resolution has spatially invariant properties. This phenomenon also occurs in the linear array case; thus, it is not a new phenomenon. In practice, if it is necessary to carry out regional observation with a super-wide view angle, and the cost is allowed, the optimal solution is to increase the arc length and the number of antenna elements. The sub-optimal solution is to realize data acquisition by the mechanical rotating of the arc array, which is low-cost but inefficient. Both solutions can realize high-resolution and super-wide view angle imaging processing.

Conclusions
In order to obtain efficient data acquisition and wide view angle imaging, the arc antenna array technology is applied into the GB-SAR system in this paper, which realizes a novel system: GB-AA-SAR. Then, we derived analytic expressions of sampling criteria, and proposed a novel wavenumber domain imaging algorithm suitable for GB-AA-SAR. Through comparing with the BP algorithm and efficient omega-k algorithm, it can be concluded that the proposed algorithm has the same accuracy as the BP algorithm, but the efficiency is higher than the BP algorithm. Further, the proposed imaging algorithm is inefficient compared with the omega-k algorithm, but can solve the problem of azimuth mismatch. Finally, we analyzed the resolutions of GB-AA-SAR and carried out multi-point targets numerical simulation, which not only verified the spatially invariant properties and the theoretical value of the resolutions, but also verified that the proposed imaging algorithm can finish highprecision imaging processing for wide azimuth view angle scenarios.

Conclusions
In order to obtain efficient data acquisition and wide view angle imaging, the arc antenna array technology is applied into the GB-SAR system in this paper, which realizes a novel system: GB-AA-SAR. Then, we derived analytic expressions of sampling criteria, and proposed a novel wavenumber domain imaging algorithm suitable for GB-AA-SAR. Through comparing with the BP algorithm and efficient omega-k algorithm, it can be concluded that the proposed algorithm has the same accuracy as the BP algorithm, but the efficiency is higher than the BP algorithm. Further, the proposed imaging algorithm is inefficient compared with the omega-k algorithm, but can solve the problem of azimuth mismatch. Finally, we analyzed the resolutions of GB-AA-SAR and carried out multi-point targets numerical simulation, which not only verified the spatially invariant properties and the theoretical value of the resolutions, but also verified that the proposed imaging algorithm can finish high-precision imaging processing for wide azimuth view angle scenarios.