Robust GMTI Scheme for Highly Squinted Hypersonic Vehicle-Borne Multichannel SAR in Dive Mode

: While there are recent researches on hypersonic vehicle-borne multichannel synthetic aperture radar in ground moving target indication (HSV-MC-SAR/GMTI), this article, which speciﬁcally explores a robust GMTI scheme for the highly squinted HSV-MC-SAR in dive mode, is novel. First, an improved equivalent range model (IERM) for stationary targets and GMTs is explored, which enjoys a concise expression and therefore offers the potential to simplify the GMTI process. Then, based on the proposed model, a robust GMTI scheme is derived in detail, paying particular attention to Doppler ambiguity arising from the high-speed and high-resolution wide-swath. Furthermore, it retrieves the accurate two-dimensional speeds of GMTs and realizes the satisfactory performance of clutter rejection and GMT imaging, generating the matched beamforming and enhancing the GMT energy. Finally, it applies the inverse projection to revise the geometry shift induced by the vertical speed. Simulation examples are used to verify the proposed GMTI scheme.


Introduction
Flight-borne synthetic aperture radars (SAR) provide tremendous potentials to generate microwave imageries of ground moving targets (GMT) [1] and stationary targets. However, due to the flight characteristics of the platforms, traditional SAR processing exhibits some limitations [2]. To be specific, an air-borne SAR possesses a small detection range, as it flies at a low altitude and a low velocity. Likewise, a space-borne SAR [3], due to the high-altitude and static orbit, has a fixed observed area and high requirements of transmitted power [4,5]. On the other hand, hypersonic vehicle-borne SAR (HSV-SAR), which operates at the altitude of 20 to 100 km and speeds in excess of 1700 m/s, could not be processed using traditional SAR algorithms [6,7]. Fortunately, its processing is also not overwhelming since its skipping orbit can be simplified to a dive orbit, as five-sixths of the entire skipping orbit has no acceleration [8,9].
The GMT indication (GMTI) plays an important role in both civilian and military applications [10][11][12], especially for the vehicle detection. There are recent researches on HSV-SAR/GMT indication (GMTI) in the literature. For side-looking HSV-SAR, refs. [6,7] reject the stationary clutters and recover the GMT imageries, and ref. [8] investigates a SAR/GMTI with a dive orbit. Finally, the stationary target [13] and GMT [14,15] imaging of HSV-SAR are studied, in combination with squint-looking and a. horizontal orbit. On the other hand, there is no GMTI scheme for the HSV-borne multichannel SAR (HSV-MC-SAR) investigating region of interest on the high-squint side of a dive orbit.

1.
Traditionally, the impacts of high-order phase terms are negligible and may be disregarded in the GMT focusing stage, especially for the case with a side-looking and horizontal orbiting SAR [14,15,20]. However, when dealing with a high-squint, a highspeed, and a dive orbit, the range model of HSV-SAR/GMTI has non-negligible high-order or coupled phase terms and may induce distortion on the GMT envelopes [6][7][8]14,15].

2.
The HSV-SAR/GMT range model, with a squint angle and a dive orbit, is complex. Some precise models comprise cumbersome expressions, thus limiting their application in subsequent processing of targets. Examples include the fourth-order range model (FORM4) [21] and the fourth-order Doppler range model (DRM4) [22]. While methods such as the advanced hyperbolic range equation (AHRE) [23] model and its modifications [24][25][26] have simplified the expressions, its accuracy is not sufficient to satisfy the high-resolution imaging requirement. In order to reduce the cumbersome range model and form a basis for stationary target processing, the wavenumberdomain imaging algorithm with modified equivalent range model (MERM) has been derived and discussed in [27], but the GMTI processing for HSV-SAR is still not available at present. 3.
The improved clutter rejection algorithms can enhance the effective accumulation for GMT and have extensive applications in the detection of faint GMT. These include the extended DPCA (EDPCA) [35][36][37] and the imaging STAP (ISTAP) [38]. Furthermore, the Deramp technique [39] enhances the ability to retrieve the ambiguity-free GMT, as it places the zeros in the DA clutter directions. However, these algorithms would suffer from GMT accumulation degradation when the limited channel number caused by the special aerodynamic characteristics of HSV-SAR are involved [6,7]. The chirp Fourier transform (CFT) and its modified algorithms enjoy a lower demand for the number of channels, because GMT is sparse in the surrounding clutter for coarse-focusing imageries. Examples are the CFT-center [20] and the CFT-two-step [6] algorithms. However, the beamforming mismatch caused by coarse cross-track velocity (CTV) may not be so readily eliminated and GMT energy decline is induced.
Furthermore, there have been some studies on CTV estimations. When it is combined with specific chirp-varying imageries, a parameter estimation algorithm with mono-channel is proposed in [40]. Inspecting the multichannel SAR (MC-SAR), refs. [41,42] calculate CTV by the along-track interferometry (ATI) operator. Furthermore, these algorithms with particular sweep function are explored that permit the CTV to be retrieved, yielding the subspace projection (SP) [43], adaptive matched filtering (AMF) [44], weighted AMF (WAMF) [45], and joint-pixel normalized sample covariance matrix (JPNSCM) [46]. Unfortunately, when faced with DA, these algorithms run into degradation of CTV estimation, envelope smearing, and decreased GMT energy. While the CFT-modified algorithm [15] rejects the DA clutter and performs the small-interval sweep on unknown CTV, it is computationally time-consuming.
To tackle these problems, this article presents a robust GMTI scheme for the highly squinted HSV-MC-SAR with a dive orbit. An improved equivalent range model (IERM) of stationary target and GMT is explored first. Then, incorporated with the proposed model, a robust clutter rejection and GMT imaging algorithm is derived in detail, especially for the case with DA. Finally, by using the inverse projection, the geometry shift of GMT can thence be revised. To summarize, this article has the following contributions:

1.
While there are recent researches on highly squinted HSV-MC-SAR/GMTI [14,15], this article, which specifically explores a GMTI scheme in dive mode, is novel. It presents the IERM of a stationary target and GMTs, performs the accurate clutter rejection with matched beamforming, and achieves the GMT imaging and location.

2.
Derivation of an IERM of stationary target and GMT. Due to its ability to transform the dive orbit into the horizontal orbit, it has a more concise expression compared to [8] and thus mitigating vertical speed impact and simplifying the GMTI processing.

3.
An improved clutter rejection with a two-step CTV sweep is established. Because of the smaller range of small-interval sweep, it has a shorter calculation time compared to [15]. Moreover, it performs a small-interval sweep of the CTV to form a matched beamforming, thus minimizing the effect of beamforming mismatch in [6,20].
This article is structured as follows. In Section 2, an IERM is introduced. A robust GMTI scheme for the highly squinted HSV-MC-SAR with a dive orbit is derived in detail in Section 3. The simulations and examples are presented in Section 4 and conclusions are drawn in Section 5.

Improved Equivalent Range Model
In order to form the basis for processing the GMTI, the IERMs of stationary targets and GMTs are presented.

IERM for Stationary Targets
For an HSV, five-sixths of the entire skipping orbit has no acceleration, so we can just focus on a uniform dive orbit [8,9]. Considering a region of interest on the squint side of HSV-SAR with a dive angle α, taking a stationary target P as an example, the range model is illustrated in Figure 1. θ 0 and R 0 represent squint angle and slant range. An HSV position is defined with altitude H, synthesized speed v, horizontal speed v x , and vertical speed v z . The instantaneous range model for a stationary target [27] can be modeled as where ξ a is the slow-time. Figure 1. Geometric model of highly squinted HSV-SAR with a dive orbit for a stationary target. Figure 2 shows the equivalent model for a stationary target. Equation (1) is simplified to a classic hyperbolic range equation (HRE), and the IERM for a stationary target can be described as Figure 2. Diagram of equivalent model for a stationary target.
Inspecting (3), the first equation denotes the synthesized speed of the platform, the second equation represents the speed projection operator in the beam center direction, and θ represents the equivalent squint angle. By utilizing the operator to transform the focusing plane onto Plane DEFP in Figure 2, one guarantees that the impacts of the dive angle can be mitigated. Although the operator simplifies the expression of the range history, there are side effects in the imaging geometry that will be analyzed later.
Applying Taylor series expansion on (2), we have

IERM for GMT
Taking channel number N as an example, n = 1, 2, · · · , N, the data acquisition model of HSV-MC-SAR for a GMT is illustrated in Figure 3, where the data acquisition plane corresponds to the Plane DEFP in Figure 2. The data acquisition model is defined as follows: an arbitrary GMT is denoted by T(x T , R 0 ); the nearest distance between the radar carrier and T is represented by R b ; the speed vector along the X -axis and Y -axis are denoted by v x and v y ; along-track velocity (ATV) and CTV of a GMT are indicated by v a and v r ; the adjacent channel distance is represented by d; the distance between the first and the n-th channel is denoted by d n = (n − 1) × d.
The GMT 2D velocities are introduced as By extending (4), the range history in the n-th channel can be formulated as where ∆ξ n = d n v. From (7), the IERM for a GMT in the n-th channel can be modeled as and v R = v sin θ + v r represent the relative ATV and CTV between the radar carrier and T, respectively. Equation (8) enjoys a more concise expression compared to [8], because the dive orbit is transformed into the horizontal orbit. Moreover, its accuracy meets the requirements of high-resolution imageries [7,27].

Robust GMTI Scheme for Highly Squint-Looking HSV-MC-SAR in Dive Mode
From Figure 4, the whole flowchart can be grouped into four areas: (1) The coarse imaging using the cubic CFT (CCFT) function, alleviating the impacts of DA and cubic phase. (2) An improved clutter rejection algorithm with a two-step CTV sweep, which has a shorter calculation time because the range of the small-interval sweep is reduced.
(3) The GMT finer-imaging algorithm is presented, in combination with the clutter-free signal. (4) By using the inverse projection, the geometry shift of GMT is corrected.

Coarsely Imaging
To alleviate the RCM and DA influences, we provide a coarse imaging algorithm in this section with an IERM. After the range Fourier transform (FT) and range compression were accomplished, the received signal in terms of a GMT and the n-th channel can be introduced as where the coefficient σ T represents the complex amplitude of a GMT, c denotes the speed of transmitted linear frequency modulation (LFM) signal, f r and f c indicate the range and carrier frequencies, w r (·) denotes the range-frequency profile, and a a (·) represents the azimuth-profile. Substituting (8) into (10), yielding  Observing the exponential terms of (11), the second component leads to range walk, the third component results in range curvature, and the fourth component contains the information of third-order RCM.

Geometry correction
It is safe to assume that the radar carrier velocity is much greater than the GMT velocity for the HSV-MC-SAR, although the 2D velocities of GMT are typically unknown. The following equivalence relation is reasonable in the coarse imaging stage.
Inspecting (11), the function of coarse range walk correction (RWC) can be constructed as Multiplying (11) and (13), we have The influences of the DA and cubic phase make it difficult to focus the GMT. Inspired by square CFT [6,20], the coarse-imaging operator with the CCFT is expressed as where The procedure of coarse-imaging is as follows where w a (·) denotes the azimuth-frequency profile, the coefficient G az represents the GMT complex amplitude after azimuth integration, f a reflects the azimuth CFT frequency, The received signal in the azimuth CFT domain can be written as where the GMT complex amplitude after range integration is indicated by the coefficients G r , and the Doppler centroid is denoted by f dc = 2v r /λ. The number of DA is generally given by B a /PRF = 2L + 1, the Doppler bandwidth is represented by B a , the abbreviation of the pulse repetition frequency is represented as PRF, l is integral, and l ∈ [−L, L]. The ambiguity number of Doppler centroid is introduced as where · denotes the maximum integer. For the baseband received signal, the azimuth frequency, CTV, and Doppler centroid are modeled as where the first blind velocity is denoted by v PRF = PRF · λ/2. For a GMT, the baseband received signal of the coarse imageries is recovered from (18) and we define Inspecting (23), the last phase is connected to the channel number, and the GMT steering vector is formulated as From (23), by setting the motion parameter of a GMT to zero, the baseband received signal of stationary clutter in the coarse imageries can be modeled as Observing (25), the clutter steering vector is introduced as Inspecting (23) and (25), by sweeping and identifying the maximum output between different ambiguity regions, the ambiguity-free GMT can be retrieved. Since only a few GMTs are present in the region of interest, it can be considered as sparse [7,15,47]. Moreover, the effect of DA stationary clutter may be non-negligible, which will be discussed in the next subsection.

Improved Clutter Rejection Algorithm with Two-Step CTV Sweep
To reduce the small-interval sweep range of CTV and minimize the effect of DA stationary clutter, an improved clutter rejection algorithm based on a two-step CTV sweep is derived. Furthermore, the cross-track parameter estimation is divided into coarse and small-interval sweeps.
The signal in terms of single-pixel can be introduced as with S T (p, q) = [S T,1 (p, q), · · · , S T,n (p, q), · · · , S T,N (p, q)] T S C (p, q) = [S C,1 (p, q), · · · , S C,n (p, q), · · · , S C,N (p, q)] T where the range-cell and azimuth-cell indexes of an arbitrary pixel are indicated by p and q.
With the coarse imageries, the geometric model of the joint-pixel is summarized in Figure 5 [15,46]. Furthermore, the reconstructed signal of the joint-pixel can be formulated as follows Pixels used for covariance matrix recovery The components of the reconstructed signal in (30) are written as S T = S T,1 p JP , q JP , · · · , S T,n p JP , q JP , · · · , S T,N p JP , q JP T ∈ N C ×1 (31) S C = S C,1 p JP , q JP , · · · , S C,n p JP , q JP , · · · , S C,N p JP , q JP T ∈ N C ×1 with S T,n p JP , q JP = [S T,n (p − (N r − 1)/2, q − (N a − 1)/2), . . . , S T,n (p + (N r − 1)/2, q + (N a − 1)/2)] ∈ 1×N r ·N a (33) where the range-cell and azimuth-cell numbers in the pixel window are represented by N a and N r , the range-cell and azimuth-cell indexes in terms of the joint-pixel are indicated by p JP and q JP , the number of pixels in a pixel window is denoted by N C = N r · N a (N − 1) + 1.
The covariance matrix with joint-pixel information in (35) is formulated as [15,46] where N JP is the sample number, and N JP ≥ 2N C − 1. E[·] and [·] H indicate the mean and conjugate transpose operations.
From (35), the weight operator can be generated, viz., where the inverting operation is represented by [·] −1 .
The CTV sweep function can be modeled as v r ,l, where the finer CTV is indicated byv r , the GMT ambiguity-free region is denoted byl, and the optimized weight operator is represented byŴ opt . The clutter elimination processors can be introduced as The clutter rejection processors for different algorithms are shown in Figure 6, where the dotted line on the left denotes the general Doppler spectrum diagram (DSD) [7,49] of a GMT and clutter. T-point, C-point, and C -point denote a GMT, clutter, and a DA component of clutter. The green and orange parallelograms represent an arbitrary ambiguity region and the GMT ambiguity-free region, respectively. M-point, which is the midpoint of dashed line CC , represents the assumed GMT. A GMT can be obtained through placing the zeros in the DA clutter directions, because the GMT is sparse for the coarseimaging imageries [7,15]. The CFT-center algorithm locates the beamforming center in the assumed GMT direction, and the CFT-two-step algorithm steers the beamforming center as determined by the coarse CTV-they both have beamforming mismatch. The proposed algorithm, which has the matched beamforming in the GMT direction, enhances the SCNR and forms the basis for reducing the GMTI processing. Next, we shall introduce the CTV estimation based on the two-step sweep, and the process can be described as follows: 1.
CTV estimation with the large-interval sweep. (a) Coarse sweeping in an arbitrary region. Perform the large-interval sweep with the sweep function of (38) in the lth ambiguity region. (b). Repeat (a) for all ambiguity regions. By sweeping and identifying the maximum output between different ambiguity regions, the GMT true region can be retrieved. (c). Determine the small-interval sweep range. The smallinterval sweep range can be determined from the best and second best CTV values generated by large-interval sweep in the GMT ambiguity-free region.

2.
CTV estimation with the small-interval sweep. Sweeping in the small-interval sweep range, the finer CTV can be estimated from the output that maintains the maximum GMT power.

GMT Finer-Imaging
After the aforementioned processing, GMT together with the accurate CTV can be obtained. Here, a GMT finer-imaging algorithm for the highly squinted HSV-MC-SAR with a dive orbit is derived, ATV estimated and the GMT velocities revised. Compared to the CFT-two-step algorithm [6], the proposed GMT finer-imaging algorithm preprocesses the swept CTV and its RWC operator is reduced to a single-step operation.
Transforming the received signal into the range frequency domain, yielding where an arbitrary ξ a is denoted by ξ 0 , the aperture center time for the GMT is expressed by ξ c . Taking the second-order Keystone Transform (SOKT) [6][7][8]15,50,51], the coupling between range frequency and square term of slow-time can be minimized, and its transform operator is introduced as whereξ a represents the transformed slow-time.
In combination with transform operator of (42), (40) is rewritten as The range frequency is much smaller than carrier frequency, and thus the following expression is reasonable, viz., Substituting (44) into (43), yielding Observing (45), the coupling between range frequency and square term of slow-time is eliminated, and it is easy to obtain the finer RWC operator and yields With the finer RWC operator, the received signal can be generated as After that, we reconstruct the second-order phase and yield with Transforming (48) into the 2D time domains, that is, The processing of simplified fractional Fourier transform (SFrFT) [52,53] is described as where the kernel operator is introduced as Γ ρ ξ a , η = exp −j η −ξ a /2 sin ρ [54]. The rotation factor is recovered by sweeping the maximum power of (51), that is [55] [ρ e ] = arg max ρ,η The ATV sweep of GMT can hence be modeled as v a = cR 0 2 f c γ est (53) where γ est = PRF 2 cot ρ e /N sa reflects the chirp rate of the received signal, and N sa indicates the sampling number in the slow-time direction.
In order to realize the GMT finer-imaging, with the swept parameter in (38) and (53), the operators of azimuth matched filtering and third-order range cell migration correction (RCMC) can be generated.
wheref a illustrates the transformed version of frequency corresponding to the transformed slow-time.
Taking azimuth compression with (54) and revising the third-order phase with (55), the GMT finer-imaging can hence be realized, yielding where the coefficientσ T denotes the GMT complex amplitude after the finer-imaging.

Geometry Correction
With the aforementioned operation, we apply an IERM to achieve the GMTI for the highly squinted HSV-MC-SAR with a dive orbit. Although the speed projection operator (e.g., range model equivalence) simplifies the entire GMTI process, its side effect is geometry shift. Utilizing the operator to transform the focusing plane onto Plane DEFP, the geometry shift is more serious compared to the classic focusing model, as shown in Figure 7. The inverse projection is an effective tool to mitigate the shift induced by the vertical speed of HSV [27,56]. The first step of the tool is to place a uniform alignment grid on the ground, and the second step is to project a GMT to the Plane DEFP. From Figure 7, β is the azimuth angle, the following relation between T (x , y ) and T(x, y) can be modeled as where the Doppler frequency of the target on the ground is represented by f d . Figure 7. Formulation of geometric correction.

Point Targets
In order to demonstrate the effectiveness of the proposed scheme for the highly squinted HSV-MC-SAR with a dive trajetory, the simulation results of point targets are provided and the simulation parameters are shown in Table 1. We consider a scattering point model comprising a 0 dB signal-to-clutter ratio (SCR) and 10 dB signal-to-noise ratio (SNR), and the geometric configuration is depicted in Figure 8. The ATV and CTV of the first GMT are, respectively, 8 and 8 m/s, and of the second GMT are 0 and 8 m/s.  First, we provide the results after range compression as shown in Figure 9a without the coarse RWC, and in Figure 9b with the coarse RWC. It is noted that the point target envelopes are basically straightened after the coarse RWC. Then, the coarse-focusing images of point targets are illustrated in Figure 10. The envelopes of point targets with azimuth FT are shown in Figure 10a, and one notices that the second GMT overlaps with the clutter. Figure 10b,c represents the coarse-imaging results, coarsely focusing image with the square CFT and coarsely focusing image with CCFT. One can notice that all point targets are roughly focused and the GMTs are now distinguished from the clutter. Moreover, the CCFT, which compensates for the high-order phase terms, outperforms the square CFT in the coarse-focusing of point targets. Applying inverse projection to revise the geometry shift in Figure 10c, the "square imagery" is yielded as shown in Figure 11, where the imagery corresponds to the distribution of point targets in Figure 8b.   Figure 12a illustrates the CTV sweep errors for different algorithms with respect to different CTV. One can see that the proposed CTV sweep algorithm obtains less error compared to the other algorithms. When faced with DA and phase errors, the other algorithms run into degradation of CTV sweep. The performance of the CTV sweep will drop when CTV is close to an integer multiple of the first blind velocity, because the GMT and clutter have close azimuth Doppler frequency and are difficult to distinguish. Figure 12b shows the errors of the CTV sweep with v r = 8 m/s, and the performance of the CTV sweep becomes better as the SCR increases. Figure 13 shows the beamforming for the proposed and the other algorithms, and the red arrow represents the true CTV of a GMT. The results show that the proposed algorithm has the highest amplitude in the red arrow direction compared with the other algorithms. Due to the lack of swept CTV for the CFT-center and CFT-two-step algorithms in this stage, they generate the beams that do not match a GMT. Fortunately, the proposed algorithm generates a matched beamforming by swept CTV.
After clutter rejection, the images are as shown in Figure 14, where the white rectangular areas represent the enlarged image of the GMT profiles. With the square CFT [6,20], the GMT images are smeared due to the inevitable third-order phase, as illustrated in Figure 14a. As shown in Figure 14b, the proposed scheme obtains better GMT envelopes because the third-order phase is corrected by CCFT and matched beamforming is generated from swept CTV [7,15].
After that, the finer-imaging processing for the first GMT is shown in Figure 15. Transforming the received signal into the time domain, its misaligned envelope is shown in Figure 15a. The image obtained is shown in Figure 15b with SOKT and shown in Figure 15c with finer RWC. One can see that the straightened envelope of the first GMT is retrieved. After third-order RCMC and geometry correction, the sidelobe of the first GMT and geometry shift are inhibited, and the finer-imaging GMT image is recovered, as shown in Figure 15d.   Finally, we further evaluated the finer-imaging result of the first GMT in Figure 15d, the contour plots are shown in Figure 16. The 2D contour plots are shown in Figure 16a, where the main-lobes and sidelobes are not coupled or "crossed". We consider the evaluation index comprising −13.27 dB at the ideal peak sidelobe ratio (PSLR) and −10.24 dB at the ideal integrated sidelobe ratio (ISLR). The azimuth slice of 2D contour plots is shown in Figure 16b, where the PSLR and ISLR are, respectively, −13.14 and −10.15 dB, which are close to the ideal values. In order to compare the proposed algorithm with existing algorithms, the results are shown in Table 2. It can be seen that the proposed algorithm has the highest SCNR compared to the other algorithms. The reason is that the ISTAP algorithm suffers degradation of energy accumulation due to the insufficient channel numbers of the HSV-SAR. The CFTcenter algorithm locates the beamforming center in the assumed GMT direction, and the CFT-two-step algorithm steers the beamforming center as determined by the coarse CTV, they both have beamforming mismatch and energy loss. On the other hand, the proposed algorithm has the best performance of GMT 2D velocity sweep compared to the other algorithms. During clutter suppression, the CFT-center algorithm does not sweep for CTV and the CFT-two-step algorithm achieves only a coarse CTV sweep. More importantly, the performance of parameter sweep becomes better as the SCNR increases. The algorithm complexities of primary steps for different algorithms are listed in Table 3, where the consumption time is mainly determined by clutter rejection. The ISTAP algorithm performs clutter rejection with 2D velocity sweep and is computationally timeconsuming. The proposed algorithm, due to the smaller range of small-interval sweep, has a shorter calculation time compared to the CFT-modified algorithm. While the complexities of the proposed algorithm are slightly higher than the CFT-center and CFT-two-step algorithms, the 2D velocity sweep is more accurate, beamforming is better matched, and extracted GMT energy is stronger.
In Table 3, the sweep numbers of ATV and CTV are represented by N na and N nr , the range pluses and the large-interval sweep number of CTV are represented by N sr and N coarse nr . The phase decoupling with SOKT consumes N sr N 2 sa multiplications and N sr N sa (N sa − 1) additions. The ATV sweep with SFrFT costs N sr N sa log 2 (N sa /2) multiplications and N sr N sa log 2 N sa additions. The RWC with the efficient Radon transform estimation (ERTE) requires N sr N 2 sa multiplications and N sr (N sa − 1) additions [6]. It is assumed that the above five algorithms have the same complexities for the GMT finerimaging, without counting the parameter estimation and phase decoupling.

Surface Targets
In this section, we evaluate the performance of the proposed algorithm by simulating the imaging of a surface moving target. The simulation parameters are listed in Table 1. The geometric configuration of a surface target is sketched in Figure 17. The surface moving target is a car that has v a = 5 m/s and v r = 5 m/s. The finer-imaging results after the proposed scheme was accomplished are shown in Figure 18a. Although the focused image of a vehicle is satisfactory, the geometry shift induced by the vertical speed is present. After geometry correction, the geometry shift is revised and the result is as illustrated in Figure 18b, very close to that sketched in Figure 17b. Therefore, the GMT finer-imaging and geometry correction of the proposed scheme are feasible.

Multiple Targets and Extended Scene
In this section, we evaluate the performance of the proposed algorithm by simulating the SAR imaging and localization of multiple targets. Figure 19a shows seven vehicles driving on a road in the region of interest. These simulated targets are then superimposed on the SAR image obtained by air-borne X-band multichannel SAR (MC-SAR), as illustrated in Figure 19b.  Figure 20a shows the finer-imaging results, one can see that the geometry shift induced by the vertical speed is present. Figure 20b illustrates the localization results after geometry correction, where the yellow and red marks represent the location and relocation results of the targets. We notice the scene is very close to that sketched in Figure 19b where the targets are marked on or beside the road. Therefore, the finer-imaging and target localization of the proposed scheme are effective.

Conclusions
Existing researches of HSV-SAR/GMT indication (GMTI) focus their region of interest on side-looking mode or squinted mode with a horizontal orbit. The proposed GMTI scheme for HSV-MC-SAR is new and focuses the region of interested on the high-squint side of a dive orbit. At first, an IERM, which provides a concise expression and simplifies the GMTI process, is explored. Then, a robust clutter rejection and GMT imaging algorithm for the highly squinted HSV-MC-SAR with a dive orbit is derived in detail, incorporated with the IERM model. Finally, the geometry shift of GMT is revised by inverse projection. Compared with the existing researches, the proposed scheme has accurate recoveries of GMT 2D speeds, matched beamforming, and satisfactory GMTI performance.