Robust Clutter Suppression and Ground Moving Target Imaging Method for a Multichannel SAR with High-Squint Angle Mounted on Hypersonic Vehicle

This paper investigates a robust clutter suppression and detection of ground moving target (GMT) imaging method for a multichannel synthetic aperture radar (MC-SAR) with high-squint angle mounted on hypersonic vehicle (HSV). A modified coarse-focused method with cubic chirp Fourier transform (CFT) is explored first that permits the coarsely focused imageries to be recovered, thus alleviated the impacts of GMT Doppler ambiguity and range cell migration (RCM). After that, in combination with joint-pixel model, a robust clutter suppression method which enhances the GMT integration, and improving the accuracy of radial speed (RS) recovery by modifying the matching between the beamformer center and GMT, is proposed. Due to that the first-order phase compensation and RS retrieval are predigested, the proposed algorithm has lower the algorithmic complexity. Finally, the feasibility of our proposed method are verified via experimental results based on simulated and real measured data.


Introduction
Recently, the hypersonic vehicle (HSV) which works in the near-space had received great interests, as it provides the tremendous potential to significantly bridge the gap between air-borne and space-borne remote sensing [1,2]. The HSV-borne synthetic aperture radar (HSV-SAR) possesses a speed advantage and a wider observed area, over the airborne SAR [3][4][5]. On the other hand, it has a flexible detection range and lower power requirements compared to the space-borne SAR [4]. Compare with the classical SAR, it exhibits the characteristics of fast responsiveness, multiple revisers, and persistent detections, due to its high-speed (e.g., 5 to 20 Mach), high-maneuverability, and global reaching [6][7][8][9].
The side-looking SAR, with a static beam direction, has narrow coverage [10][11][12][13][14]. Whereas the squint-looking SAR [15], due to its variable beam direction, achieves the flexibility of able to revisit interesting regions. More importantly, a high-squint SAR [16,17], with a squint angle of over 45°, is capable of detecting the ground moving targets (GMT) for the hidden areas that the side-looking SAR cannot match [18].
To the best of our knowledge, most existing research have been provided for HSV-SAR/GMT indication with side-looking or single-channel squint-looking. To be specific, Reference [3] focused on HSV-SAR imaging with squint angle. In [4], clutter cancellation and GMT imaging are investigated, in combination with the side-looking HSV-SAR, [5] deals with a HSV descending stage and also did some simulation analyses. Reference [6] primarily concentrated on GMT imaging of squint-looking HSV-SAR in a single-channel To reduce the DOF demanded for the unambiguous GMT extraction, the chirp Fourier transform (CFT) is a potential tool to minimize the Doppler ambiguity impacts [37]. Reference [19] requires approximately half the DOF of ISTAP method, compressing the Doppler spectrum, eliminating the clutter, and retaining the unambiguous GMT (we call it CFTbased method). Then, Reference [4] uses the coarse RS to optimize the beam-former and ameliorate the GMT integration (we call it CFT-improved method). However, the clutter cancellation beam-formers of CFT-based and CFT-improved methods are sensitive to the RS recovery accuracy process, which may lead to a mismatch between the beam-former center and GMT. Moreover, for a high-squint SAR mode, a side-looking SAR mode on which the GMT imaging of CFT-improved method is based is not valid.
To retrieve the target directions or RS, the wide-band polynomial-phase are investigated, in combination with the sensor arrays in [38]. Reference [39,40] focused on target direction estimations with narrow-band signal, Reference [41] deals with the clutter. However, combining these methods with SAR imaging models also may not be so easy. On the other hand, the adaptive matched filtering (AMF) [42] and the subspace projection (SP) [43] methods can be employed to implement the motion parameter recovery, but these methods would suffer from quality reduction when the coarsely focused imageries with co-registration error are involved. The joint-pixel model is developed for the cost function of RS search, especially for the case with a co-registration or channel phase errors. These methods such as the weighted AMF (WAMF) improves the ideal steering vector [44], while [45] concentrate on the heterogeneous clutter. However, the Doppler ambiguity or large relative horizontal speed (HS) may not be so easily removable and subsequently GMT profile distortion is induced.
To cope with these challenges, this paper proposes a robust clutter suppression and GMT imaging method for the high-squint HSV-MC-SAR. Coarse range cell migration correction (RCMC) and cubic CFT are conducted first, and the coarsely focused imageries are retrieved. After that, the recoveries of RS and unambiguous GMT for the coarsely focused imageries are achieved, in combination with the joint-pixel model. Finally, with the proposed CFT-improved technique, the GMT accurate-focused method for the high-squint is derived. Compare with the existing works, the proposed method has the following innovations: 1.
The modified coarse-focused method with cubic CFT would alleviate the impacts of GMT Doppler ambiguity and RCM, and form a basis for the subsequent motion parameter recovery, clutter suppression, and GMT imaging; 2.
The robust clutter suppression method, due to that the recovered RS is developed to modify the matching between the beam-former center and GMT, reduces the unwanted integration loss in the desired GMT direction and ameliorates the SCNR; 3.
By alleviating or inhibiting the effects of clutter Doppler ambiguity, channel phase mismatch, and co-registration error, the RS recovery with joint-pixel model is more precise; 4.
The proposed GMT imaging method, due to that the parameter recovery is preprocessed and the first-order phase correction is simplified to a single-step process, has lower computational burden compared to the CFT-improved method [4]. In the case of high-squint and high-speed, the impacts of GMT 2-D speeds are accurately compensated, and accurate-focused imagery of GMT can hence be recovered.
This paper is organized as follows. In Section 2, the high-squint HSV-MC-SAR is described, and the range histories are provided based on data acquisition geometry. A robust clutter suppression and GMT imaging method is derived in Section 3, and RCM effect is analyzed in Section 4. In Section 5, the proposed method is verified by simulation experiments and compared with existing methods. Finally, Section 6 concludes the whole article.

Signal Model
A GMT range history of HSV-MC-SAR is derived in high-squint mode. With reference to the geometric configuration summarized in Figure 1, taking five antennas as an example (e.g., N = 5). H and v represent altitude and speed of HSV, θ (θ > 45 • ), and ϕ indicate the squint and pitch angles. W r denotes the swath width in range, R 0 indicates the slant range of a GMT at the initial position, and R b indicates the nearest vertical distance between the radar platform and a GMT. The effective phase center (EPC) model [46] in slant range plane as illustrated in Figure 1b. d and d n denote the distances between adjacent channels and from the 1st to the n-th receiving channels, and d n = (n − 1) × d, n = 1, 2, · · · , N. Suppose the initial position of a GMT is in the scene center, and located at T(x T , R 0 ). t m is the azimuth time along the HSV flight direction, and the initial position time is defined as t m = 0. In the data acquisition plane, the speeds are along the trajectory parallel to the X-axis and Y -axis, namely, v x and v y , also, v a and v r indicate the HS and RS of a GMT, respectively.
The instantaneous range history between the n-th channel and a GMT can be defined as For the HSV platform with high-squint and high-speed, the second-order Taylor expansions of instantaneous range history is no longer adequate. By expanding (1) into third-order Taylor series, one can yield where v a_rel = −v cos θ + v a and v r_rel = v sin θ + v r reflect the relative HS and RS between a GMT and the radar system. HS and RS are expressed as v a = v y sin θ + v x cos θ and v r = v y cos θ − v x sin θ.

Robust Clutter Suppression and GMT Imaging Method for the High-Squint HSV-MC-SAR
In this section, a robust clutter suppression and GMT imaging method for the highsquint HSV-MC-SAR is introduced in detail, and the flowchart summarized the processing steps as in Figure 2. (1) In the coarse-focusing stage, coarse range cell migration correction (RCMC), together with the Doppler compression are conducted. (2) The robust clutter suppression method with joint-pixel model retrieves precise RS, rejects the static clutter plus its ambiguous components, and effectively alleviates the unwanted GMT integration loss.
(3) The derivation of GMT accurate-focusing is performed under a high-squint condition. Note that the three parts of the flowchart are correlated with Sections 3.1-3.3.

Modified Coarse-Focusing of GMT and Clutter
In order to mitigate or minimize the impacts of RCM and GMT Doppler ambiguity, we describe here a modified coarse-focused method with cubic CFT for the GMT and clutter.
The echoes of GMT in the n-th channel can be defined as where τ n = 2R n (t m , R 0 )/c represents the time-delay for the n-th channel, σ T is the GMT magnitude,t indicates distance-time, and γ denotes the slope of linear frequency modulation (LFM) signal. a r (·) and a a (·) represent the window functions of range and azimuth in the time domain. λ and c reflect the transmitted carrier wavelength and propagation speed. After the range Fourier transform (FT) and range compression, (3) is rewritten as where the range and carrier frequencies indicated by f r and f c , and the distance-envelope represented by w r (·) in the Doppler domain.
For the HSV-SAR, GMT speeds are far less than the platform speed, and the following relationships are approximated in the coarse-focusing stage.
Substituting (2) and (5) into (4), we have It should be noted that the third, fourth, and fifth exponential terms in (6) are associated with range walk, range curvature, and cubic RCM terms, respectively. Note that the impacts of quadratic and cubic phases need to be calibrated, and the fourth-order phase can be disregarded that will be discussed later.
The coarse range walk correction (RWC) function can be derived from (6), and this yields The first and second phase terms situated on the right side of (7) can be employed to calibrate the range walk and channel phase errors. Multiplying (6) and (7), the range envelope alignment is basically achieved and yields The HSV-SAR systems, with their severe Doppler broadening or ambiguity echoes, have difficulty in accurately indicate the GMTs. CFT, as an extension of FT, and with azimuth pulse compression capabilities, can mitigate or minimize the Doppler ambiguity that the traditional FT cannot [19]. However, the impact of high-order phase errors is inevitable for the high-squint HSV-MC-SAR due to its high maneuverability. Thus the second-order Taylor series of range history on which the traditional dechirp and CFT-based plus CFT-improved methods are based is no longer valid. By extending the secondorder CFT [4,19] to the cubic CFT, the compression function in the n-th channel can be formulated as To alleviate the RCM and GMT Doppler ambiguity impacts, the coarse-focusing processors can be expressed as where w a (·) represents the azimuth-envelope in the Doppler domain. G az is the signal magnitude after azimuth accumulation, f a reflects the azimuth CFT frequency, and f d = 2v r ( f r + f c )/c. Inspecting (11), we can find that the cubic CFT simultaneously compresses the echoes in azimuth, compensates the cubic RCM, and calibrates the envelope offset brought by the channel errors.
Then, after the range inverse FT (IFT), the coarse-focused GMT signal in the n-th channel can be obtained as where G r denotes the signal magnitude after range accumulation, and f dc = 2v r /λ indicates the Doppler center of GMT. The Doppler ambiguity number can be defined as B a /PRF = 2L + 1, and PRF represents the pulse repetition frequency. the original Doppler spectrum is divided into 2L + 1 ambiguity areas, and B a indicates the azimuth Doppler bandwidth. L is a positive integer, and l ∈ [−L, L]. The GMT Doppler center and the CFT frequency spectrum both have Doppler ambiguity, and the ambiguity number of Doppler center defined as where · reflects the ceiling operator. According to the Doppler spectrum diagram [47], baseband frequency is given by For a baseband signal, RS and Doppler center of GMT can be calculated by where v PRF = PRF · λ/2 reflects the first blind speed.
Assuming that a GMT is located in the l-th ambiguity area, the coarsely focused imagery of GMT for the Nyquist band is formulated as One notices that only the last exponential term contains the space information of channel, and the steering vector of GMT can hence be defined as Zeroing the GMT speeds in (17), the expression of coarsely focused clutter imagery is retrieved and we define Similarly, the steering vector of stationary clutter can be written as From (17) and (19), one can see that clutter plus its ambiguous components are located in same channel. In fact, only few GMTs are located in covering areas and coarsely-focused GMTs are sparse in the surrounding stationary clutter, and thus the Doppler ambiguity impact of GMT has been mitigated [48]. Besides, it is necessary to reject the clutter Doppler ambiguity to indicate the GMT, this will be discussed in detail below.

Robust Clutter Suppression Method with Joint-Pixel Model
In order to recover the RS so as to achieve a satisfactory quality of beam-former and extract the GMT from the Doppler ambiguity echoes, a robust clutter suppression method which modifies the beam-former through precise RS recovery, eliminates the clutter plus its ambiguous components and maintains the GMT profiles is next introduced.
Here, we utilize the joint-pixel model to effectively retrieve RS and eliminate clutter, especially when the co-registration error and channel phase mismatch exist [44,49]. Figure 3 provides the joint-pixel model for the coarsely focused imagery. For the channel 1, image 1 with single-pixel model is given by Pixels used for covariance matrix recovery The joint-pixel model is conducted for motion parameter recovery, which is implemented through the detected and adjacent pixels of coarsely focused imagery, and the hypotheses can be defined as follow Inspecting (22), the GMT and clutter signals of coarsely focused imagery with jointpixel model can be formulated as S C = [S C,1 (p, q), · · · , S C,n (p, q), · · · , S C,N (p, q)] T ∈ C N C ×1 where S T,n (p, q) = [S T,n (p − (N r − 1)/2, q − (N a − 1)/2), . . . , N a and N r (N a = N r = 3) indicate the unit number in the azimuth and range directions of the selected window, p and q denote the range and azimuth pixels, [·] T indicates transpose operator, and N C = N r · N a (N − 1) + 1 represents pixel number of a joint-pixel data. The covariance matrix based on joint-pixel model and pixel offset vector [45] are formulated as where N JP is the sample number of joint pixel, and N JP 2N C − 1. E[·] and [·] H indicate the mean and conjugate transpose operators. ξ τ e jϕ τ (τ = 1, 2, ..., N C − 1) denotes the signal off-set between the auxiliary and reference channels, which are affected by the co-registration or channel phase errors. After that, using optimal pixel offset vector in (27) to alleviate GMT steering vector [45,50], yielding a whereâ T,l (v r ) = [1, a T,l (v r ) ⊗ ones(N r · N a , 1)], and ⊗ reflect the Hadamard and Kronecker operators.
Inspecting (30), the expression of weight vector can be given bŷ where [·] −1 reflects the inverse operator. Thereby, the RS search has the following form wherev r0 denotes the precise RS, which corresponds to the unambiguity areal of GMT, W opt reflects the optimal weight vector. Finally, the results of clutter suppression can be expressed as follows Figure 4 shows the beam-former of clutter cancellation for different methods, where the red and blue lines represent the ambiguous Doppler spectrums of GMT and clutter, and the profiles are shown by red and blue dots. To be specific, the CFT-based method [19] assumes that the beam-former is aligned to the middle of unambiguous area (e.g., corresponding to the green dot), which should consider the impacts caused by the beam-former mismatch. The CFT-improved method [4] which retrieves a rough RS via a large interval retrieval, relieves the beam-former mismatch, but the unwanted integration loss in a desired GMT direction still exists to vary degrees. In contrast, the beam-former center of the proposed method matches a GMT, which means that the unwanted integration loss can be excellently avoided.
The RS recovery can be divided into two steps:

GMT Accurate-Focusing
After the above processes, the ambiguity-free GMT, together with the precise RS, can be retrieved. In this subsection, we propose an accurate GMT imaging method for the high-squint HSV-MC-SAR, which estimates HS and inhibits the HS and RS impacts. Compared with the CFT-improved method, the proposed imaging method has fewer steps, preprocessing the RS retrieval and simplifying the first-order phase correction.
After range FT and azimuth IFT were accomplished, the GMT signal can be formulated as where where t 0 and t c denote arbitrary azimuth-time and synthetic aperture center time. A HSV-MC-SAR exhibits a speed advantage compared to the air-borne SAR. The arbitrary azimuthtime is close to the synthetic aperture center time, and, thus, the derivation of (34) are valid.
To decouple the f r and t 2 m , the second-order keystone transform (SOKT) [4,52,53] has proven to be a reliable tool, and it can be expressed as wheret m is called new azimuth-time after keystone transform. Substituting (40) into (34), it yields Clearly, f r f c is valid, and the approximate formulation can hence be obtained as Applying (42) in (41), we have Obviously, the decoupling of f r andt m is achieved. By using the precise RS to revise the RWC function, its accurate expression can be generated.
After an accurate RWC was performed, the echo signal can be expressed as follows To recover the HS, the phase reconstruction processors are derived as where After range IFT, viz., The simplified fractional Fourier transform (SFrFT) is a reliable way for HS recovery, especially when the impact of Doppler ambiguity is mitigated [54,55]. We apply the SFrFT processors and yield where Γ ρ (t m , η) = exp(−j(η −t m )/2 sin ρ) denotes the kernel function [56]. The rotationangle operator denoted by is estimated by searching the maximum output of the following formulation [57], that is, Then, the HS calculation can be conducted as follows: v a = cR 0 2 f c γ est (51) where γ est = PRF 2 cot ρ e /N nan reflects the Doppler rate, and N nan indicates the number of processing data in azimuth. To achieve the accurate focusing, the second-order azimuth compression function and the correction function of cubic RCM are constructed from the recovered values of GMT's motion parameters, we have wheref a indicates the new Doppler frequency fort m . After applying the azimuth compression and cubic RCMC, the accurately-focused of GMT in the 2-D time domains can be expressed as whereσ T illustrates the magnitude of accurate-focused GMT.

RCM Analysis
Here, the RCM impacts of HSV-SAR for different platform speeds or squint angles are discussed in detail. Figures 5 and 6 show the simulation results of RCM for different squint angles or platform speeds, where simulation experiments are conducted with specific parameters listed in Table 1. It can be seen from Figures 5a and 6a that with the increase in the squint angles or platform speeds, the range walk becomes larger. In Figure 5b, as the squint angle increases, the range curvature decreases. In the case of θ = 50 • , the range walk is much greater than one range resolution unit, and the range curvature is less than 1 m. Thereby, the coarse RWC can be performed without considering the range curvature effect.   For the HSV, taking the platform speed v = 2380 m/s, for example, the relationship between the neighboring channel errors and squint angles is as illustrated in Figure 7. In the high-squint mode, the envelope error of neighboring channel after range compression has exceeded three range resolution units. To calibrate the imagery, it is necessary to compensate the channel errors.

High-Order Phases
With the simulation parameters in Table 1 and the platform speed v = 2380 m/s, the high-order phases for different squint angles are summarized in Figures 8 and 9. The spatial variations of cubic RCM are shown in Figure 8a. Additionally, for a squint angle over 45 • , the cubic RCM has exceeded one hundredth of azimuth resolution, and its cubic phase error in azimuth more than 1 rad as shown in Figure 8b. From Figure 9, the phase errors of fourth-order are less than π/4 in a high-squint case, and its impact on focusing quality is negligible and may be ignored. Therefore, in a high-squint HSV-SAR mode, the third-order Taylor expansions of instantaneous range history can satisfy the accuracy of high-resolution imagery.

Simulation Results
In this section, the proposed method was demonstrated to be effective on clutter suppression and GMT imaging in simulation data of high-squint HSV-MC-SAR, where the specific parameters of simulation examples are listed in Table 2. The distribution diagram of point scatterers are sketched in Figure 10, and the signal-to-clutter ratio (SCR) and signal-to-noise ratio (SNR) were preset to be 0 and 10 dB, respectively. Additionally, for the first GMT, HS and RS are noted to be 0 and 14 m/s, and the 2-D motion parameters of second GMT were both set at v a = v r = 14 m/s.  First, the processing results of two GMTs plus clutter are summarized in Figure 11. Figure 11a illustrates the range compression results without coarse RWC, one notices that these profiles are severely oblique, which indicates that the effect of range walk is not calibrated. The processing results after coarse RWC are shown in Figure 11b, and the profiles of two GMTs and clutter are both roughly aligned. After that, the signal processing results of two GMTs and some clutter in the range compression and azimuth Doppler frequency domain are summarized in Figure 12. From Figure 12a, the point scatterer profiles after azimuth FT are shown first, and the local enlarge results are illustrated in Figure 12b. We notice that the profile aliasing of two GMTs exists, and the profiles of first GMT and clutter are garbled. The classical FT operation, due to that the phenomenon of Doppler ambiguity is present and GMT trajectory are submerged in the surrounding stationary clutter, brings tremendous challenge to the subsequent clutter suppression and GMT imaging. Figure 12c illustrates the profiles of point scatterers after azimuth second-order CFT operation [4,19], where we can see that two GMTs are now separated from the surrounding stationary clutter. By using the operation to compress the Doppler spectrum, the smeared coarsely-focused imagery of two GMTs and clutter is recovered. Because the second-order CFT would suffer from quality degradation in a highsquint case without high-order RCMC. After cubic CFT in azimuth, the coarsely focused imagery of two GMTs and clutter as shown in Figure 12d. One notices that the smeared imagery of two GMTs does not occur, which means the impacts of Doppler ambiguity and RCM are alleviated. Next, Figure 13a shows the quality assessment results of RS recovery with respect to RS from 9 to 15 m/s. We notice that the proposed RS recovery method yields somewhat smaller absolute error as compared to the other methods. The performance of the AMF and SP methods declines when the co-registration and channel phase errors are involved. The WAMF method, due to the power heterogeneous sample pixels, suffers from performance loss of RS recovery [45]. With the searching that maximizes the magnitude, the parameter recovery quality will decrease, even when RS approaches an integer multiple of the first blind speed (a GMT and clutter in close Doppler frequency unit).
In v r = 12 m/s case, the absolute errors of RS recovery for different SCR are summarized in Figure 13b. One notices that as the SCR increases, the absolute errors of parameter recovery become smaller. For the different clutter suppression methods, the magnitude of beam-former is plotted as a function of the RS, as shown in Figure 14, where the green arrow reflects real RS direction. One notice that the proposed method outperforms the other two methods in the desired GMT direction, and the quality of GMT integration is remarkably increased. The reason is that the RS is unknown in clutter suppression stage for the CFT-based and CFT-improved methods, which induces a mismatch between the beam-former center and GMT, whereas the beam-former of the proposed method is optimized by precise RS. After clutter suppression, the results as summarized in Figure 15, and the enlarged view of the GMT envelopes corresponds to the white rectangle. Figure 15a shows the GMT profiles after second-order CFT and robust clutter suppression, we notice that the envelopes of two GMTs are distorted, due to the presence of cubic phase errors. Figure 15b illustrates the GMT extraction results of our proposed method, one notices that the GMT profiles can be recovered extremely well, because the coarse RCMC is conducted before clutter suppression and an accurate beam-former is formed.
Then, the accurate focusing processors of second GMT as summarized in Figure 16. The responses in the 2-D time domain are shown in Figure 16a, and the profiles is seen to be misaligned. Figure 16b,c show the processing results after decoupling and accurate RWC, one notices that second GMT profiles is satisfactory aligned. To aid further rejection of the GMT sidelobe, we correct the cubic RCM and finally yield the point-like GMT, the accurate-focused imagery as illustrated in Figure 16d.
Furthermore, Figure 17 shows the imagery quality assessment of second GMT after accurate focusing. From Figure 17b,c, the peak sidelobe ratio (PSLR), whose ideal value of −13.27 dB, are seen to be −13.17 and −13.20 dB, in azimuth and range, respectively. The integrated sidelobe ratio (ISLR) of azimuth and range are increased by 0.16 and 0.09 dB compared to its ideal value of −10.24 dB [58]. Therefore, the proposed GMT accurate-focused method is capable of producing the satisfactory imagery results, and its effectiveness has been demonstrated. To assess the performance for ISTAP and CFT plus its modified methods, the primary quality assessment indexes are summarized in Table 3. Among the four methods above, the proposed method produces the highest SCNR, precise GMT 2-D speeds, and in comparatively short calculation time. Note that the results in Table 3 are the averages of 500 Monte-Carlo experiments tested on a personal PC with intel i7 CPU@3.6 GHz and 16 GB RAM, the codes are run by MATLAB 2016a. Further observations from Table 3 include the followings: First, the SCNR of the proposed method outperforms the other three methods. To be specific, the ISTAP method, due to that the radar system parameters used in this simulation cannot provide sufficient DOF, and, thus, is unable to effectively carry out GMT integration. On the other hand, during clutter suppression, the CFT-based method assumes that RS is half the first blind speed, and that the coarse RS of the CFT-improved method is obtained through the large interval retrieval. Therefore, the CFT-based and CFT-improved methods suffer from the problems of beam-former mismatch and unwanted integration loss.
Second, for the proposed method, the absolute error of GMT 2-D speeds is relatively small compared to the other methods. The reason is that the performance of motion parameter recovery is sensitive to SCNR and the proposed method is able to achieve the highest SCNR among all methods. Moreover, GMT parameter recovery accuracy of the CFT-improved method is slightly inferior, because it takes the finer retrieval of GMT 2-D speeds during GMT accurate-focusing and exhibits a higher false alarm probability [59].
Finally, the calculation times of the proposed and the other methods are quantified. In the clutter suppression step, the computation complexity of the ISTAP method is N nr N na × O N sr N 3 sa N 3 , where N sa and N sr represent the sampling number of simulation data in the azimuth and range directions, the searching numbers of HS and RS are denoted by N na and N nr . The computation complexities of the other three methods are both O N sr N 3 sa N 3 . During the GMT focusing, the precise RS recovery of CFT-improved method burdens the radar systems, since it adds N sr N 2 sa multiplications and N sr (N sa − 1) additions. For all the above algorithms, the computation complexities of clutter suppression are much greater than that of GMT focusing, and, thus, the calculation times of the ISTAP method is the highest. Since the RS estimation is preprocessed and the RWC is simplified, the proposed method consumes less computation time than the CFT-improved method. Although the computation complexity of our proposed method is not better than that of the CFT-based method, the recovery of the GMT speeds is better and the results of GMT imaging is also better in a high-squint case.

Real Data Results
Here, the proposed method was verified to be effective in relocation of real measured data. The parameters of the SAR systems are shown in Table 4. The targets are vehicles driving on a road in the coverage region, the distribution diagram of the scene is as shown in Figure 18 [50].  Figure 18. Sketch of the experimental scene. Figure 19 summarizes the relocation results of four different methods, where the blue and red circles represent the location and relocation positions of the moving targets. In Figure 19a,b, the relocation results show that most of the moving targets are not on the road. The reason is that the phenomenon of co-registration and channel phase errors occur, resulting in the failure of AMF and SP methods. In contrast, from the relocation positions of WAMF and our proposed methods as illustrated in Figure 19c,d, one notices that the moving targets are indicated on or beside the road. More importantly, the proposed method, which calibrates the impact of 2-D spends during GMT accurate-focusing, has a better relocation results compared to the WAMF method.

Conclusions
The existing GMT indication studies for SAR are usually based on air-borne or spaceborne systems, while most researches of HSV-SAR GMT indication work in a side-looking or single-channel squint-looking mode. This paper explores a robust clutter suppression and GMT imaging method for the high-squint HSV-MC-SAR. These problems (high-order phase error, motion parameter coupling, Doppler ambiguity, antenna size limitations) brought by the high-speed and high-squint are analyzed and solved, and the mismatch between the beam-former and GMT is also minimized. At first, the RCM impact is analyzed and corrected, and the point scatterers are coarsely focused. Then, robust clutter suppression method based on joint-pixel model is proposed, which has the ability to eliminate the static clutter plus its ambiguous components and mitigate the unwanted GMT integration loss. Finally, accurately-focused GMT imagery is obtained under a high-squint condition. Compared to the existing methods, the proposed method has the highest SCNR and motion parameter recovery accuracy, and requires shorter calculation times.
Future works include extending the method in this article to GMT location, and skipping trajectory should be considered to improve this method and enhance the ability to adapt to more complex flight trajectories.