Ground Moving Target Imaging for Highly Squint SAR by Modified Minimum Entropy Algorithm and Spectrum Rotation

Ground moving target (GMT) is displaced and defocused in conventional synthetic aperture radar (SAR) image due to the residual phase error of non-cooperative GMT motion. In this study, a GMT imaging (GMTIm) method is proposed for highly squint SAR. As the squint angle become large, the displace and defocus effect of the GMT image become severe and the geometry distortion of the GMT image cannot be ignored. The proposed method first deduced the two-dimensional (2-D) frequency domain signal of the GMT and the bulk compression function of the Range Migration Algorithm (RMA) in highly squint SAR. Then GMT ROI data are extracted and a modified minimum entropy algorithm (MMEA) is proposed to refocus the GMT image. MMEA introduces the idea of bisection into the iteration process to converge more efficiently than the previous minimum entropy method. To overcome the geometry distortion of the GMT image, an equivalent squint angle spectrum rotation method is proposed. Finally, to suppress the GMT image sidelobe, the sparse characteristic of GMT is considered and a sparse enhancement method is adopted. The proposed method can realize GMTIm in highly squint SAR where the squint angle reaches to 75 degrees. The PSNR and ISLR of point target in highly squint SAR is close to that in side-looking SAR. The simulated point target data and ship data are used to validate the effectiveness of the proposed method.


Introduction
Ground moving target imaging (GMTIm) is an important application of synthetic aperture radar (SAR) [1][2][3][4][5][6]. The motion of ground moving target (GMT) causes the phase error which makes the SAR image of GMT defocused and smeared. The conventional SAR imaging algorithm, which is designed for stationary scene, is not capable of compensating the phase error. The movement of GMT would cause the Doppler frequency shift and the change of the azimuth frequency modulation, which are closely related to the image quality of GMT [7][8][9][10]. Highly quality image of GMT contributes to further SAR application, such as moving targets identification, classification etc. Therefore, it is worthwhile to study the refocusing of GMT in SAR images. The problem of GMTIm can be decomposed into clutter suppression and refocusing.
In most practical situations, GMT has low signal-to-clutter ratio in the received echo, so the clutter suppression is needed. Along-track interferometry (ATI) and displaced phase center antenna (DPCA) are two typical technologies of clutter suppression [11][12][13]. To achieve DPCA and ATI, receiving channels redundancy is needed, which may increase system complexity. Region of interest (ROI) [14][15][16] is another clutter suppression method. Zhang et al. [16] proposed a GMTIm method based on (ROI) with single channel SAR system, which uses ROI data containing defocused GMT complex image to achieve the parameter. Secondly, to increase the spectrum selection region and overcome the geometry distortion of GMT image, we utilize the equivalent squint angle spectrum rotation method to convert the 2-D spectrum signal of the GMT ROI data from the squint-looking to side-looking form. Finally, we use the sparse enhancement algorithm to suppress the cross sidelobe of the GMT image [39]. Experiments on simulated GMT data and ship data demonstrate the effectiveness of the proposed novel method.
This study is organized as follows. In Section 2, the signal model is constructed. Section 3 is the detail of the proposed MMEA, the equivalent squint angle spectrum rotation and sparsity enhancement. In Section 4, simulated data experiment is provided to evaluate the performance of the proposed method. Section 5 provides the conclusion.

Signal Model
Consider the squint-looking SAR geometry as illustrated in Figure 1, where the xaxis is the azimuth direction and r-axis are the range direction. It is assumed that the SAR platform moves at the speed of v along the azimuth direction, the wave beam squint angle is θ . The range fast time and azimuth slow time are denoted as t and s t , respectively.
The n-th pulse is transmitted at slow time Assume that the radar transmits the linear frequency modulation (LFM) signal. Then the base band echo can be written as  The instantaneous slant range between SAR platform and the GMT is expressed as Assume that the radar transmits the linear frequency modulation (LFM) signal. Then the base band echo can be written as where the rect(·) is the rectangular function, T p is the pulse width, T a is the synthetic aperture time, f c is the carrier frequency, γ is the chirp rate, and σ is the target scatter coefficient that assumed to be a constant. c is the light speed. To derive the 2-D spectrum formular of the GMT echo, firstly obtain the range frequency domain and azimuth time domain signal by taking the range Fourier transform (FT) of s r (t, t s ) the azimuth FT of the signal s 1 ( f r , t s ) can be expressed as where f r and f a are range and azimuth frequency, respectively. According to the principle of stationary phase (POSP), the 2-D spectrum expression can be written as where To simplify the expression, envelope functions and scattering coefficient are ignored, which does not influence the following procedure. According to the RMA, the bulk focusing reference function of squint-looking SAR can be constructed as where R re f is the reference distance. Multiplied with Equation (9), Equation (5) becomes It can be seen from Equation (10) that after the 2-D inverse Fourier transform (IFT) of Equation (10), the defocused GMT image can be obtained because of the existing of the high-order residual phase items. The level of defocus varies among the GMT which are in different position and with different velocity. For the stationary targets, the distance between the target point and the reference point becomes larger, the residual phase value is larger and the phenomenon of defocus becomes more severe. For GMT, in addition to the distance between the target point and the reference point, the velocity of GMT also affects the refocusing effect. The defocus caused by distance can be eliminated by accurate Stolt interpolation, while the defocus caused by velocity need to be solved by compensating for the residual phase related to velocity. The Stolt interpolation is expressed as after adopting the Stolt interpolation, Equation (10) becomes If the residual phase containing the unknown velocity is compensated, the refocused GMT image can be obtained immediately. However, due to the squint angle of SAR, the preliminary refocused GMT image still exists geometry distortion, thus geometry correction is needed. Next, sparsity enhancement is implemented to suppress the cross sidelobe of the GMT image. Geometry correction and sparsity enhancement is in the next section. For the problem of Doppler ambiguity which may occur in practice, detailed analysis has been stated in [38] and effective methods of removing Doppler ambiguity have been proposed by series of literatures [17,32,40]. Increasing the pulse repetition frequency (PRF) can also avoid the Doppler ambiguity. In this study, it is assumed that the Doppler ambiguity is removed by mature techniques.

Proposed Method for Highly Squint SAR GMTIm
As is illustrated in Figure 2, the proposed method for GMTIm can be divided into three main steps: ROI data extraction, parameter estimation and refocusing. The main operation of ROI data extraction is extracting the GMT region within the preliminary image result after turning (12) into time domain. The key process of parameter estimation is the proposed MMEA method, which introduces the idea of bisection into the iteration minimum entropy algorithm to increase the iteration speed. MMEA firstly sets the parameter value range calculated by possible GMT velocity, and then shrinks the value range according to the first-order derivative of entropy until the value range is close to zero. Compensate the residual phase error with the phase compensation function containing the estimated parameter. Then, come to the refocusing step which has two main sub-steps: spectrum rotation and sparse enhancement. Spectrum rotation and sparse enhancement is used to eliminate the geometry distortion and suppress the sidelobe, respectively.
The first step of the proposed method is extracting ROI data, which are the defocused complex sub-images containing the GMT. ROI data can be easily extracted from the whole SAR image and its 2-D frequency domain is expressed as where f r and f a are range and azimuth frequencies of ROI data. Assume that the size of original data is N r × N a , and the size of ROI data is n r × n a , where N r , n r denote the range samples and N a , n a denote the azimuth samples. The ROI data usually occupy a small part of the original image, which means the volume of data is reduced. The next of the proposed method is on the basis of ROI data.

MMEA for Parameter Estimation
It can be seen from Equation (13) that the defocus effect of GMT mainly comes from v e which contains the unknown velocity of GMT and causes the phase error. When the residual phase error is compensated by the phase compensation function with the correct velocity and 2D IFT is conducted, the preliminary refocused GMT image is obtained. 1/v 2 e can be considered as a single parameter α. In this section, MMEA is proposed to estimate the unknown parameter α. The residual phase compensation function can be constructed as where α is the estimated phase compensation parameter and R ref is the reference distance. Ideally, the estimated parameter α is equivalent to the actual parameter. The residual phase compensation function H 2(α) f r , f a can be written as a matrix H(α) when the parameter α is given, and the preliminary refocusing procedure can be formulated through a series of transformation: where F −1 a and F −1 r is the azimuth and range inverse Fourier transform matrix, respectively, S ROI is the matrix form of 2-D frequency domain ROI data and ⊗ denote the Hadamard product. Z(α) is the refocusing GMT image which is the function of parameter α. The Z(α) becomes the refocused image when the parameter α is the correct one.
(28) The first step of the proposed method is extracting ROI data, which are the defocused complex sub-images containing the GMT. ROI data can be easily extracted from the whole SAR image and its 2-D frequency domain is expressed as where r f and a f are range and azimuth frequencies of ROI data. Assume that the size of original data is × r a N N , and the size of ROI data is × r a n n , where r N , r n denote the range samples and a N , a n denote the azimuth samples. The ROI data usually occupy a small part of the original image, which means the volume of data is reduced. The next of The key idea of minimum entropy algorithm is that the refocused GMT image has the minimum image entropy. Iterate the parameter α until the GMT image reaches the minimum image entropy and the parameter of this time is the correct one. The GMT image entropy function is formulated as where k is the base of the logarithm, n r and n a are the numbers of range and azimuth samples, respectively. z (m,n) (α) is the pixel of n-th column and m-th row of the GMT image Z(α). It can be seen from Equations (15) and (16)  only determined by parameter α, and the refocusing procedure can be also transformed to a optimization problem as It is not easy to solve this problem directly. To convert this problem into a series of easier sub-problems, the IMEA constructs a surrogate function of (17). Different from IMEA, MMEA uses the idea of bisection to effectively obtain the optimal parameter α. The parameter interval and the first-order derivative of E(α) at the bisection point are two key points of MMEA. MMEA has two main steps at every iteration. The first step is obtaining the parameter interval [a, b] according to the radar platform velocity v and possible GMT velocity [v x ,v r ]. Then calculate the first-order derivative of E(α) with respect to parameter α at the bisection point (a + b)/2. The first-order derivative of E(α) can be written as where z i (α) and z r (α) denote the imaginary and real parts of z (m,n) (α), respectively. The specific expression and derivation process can be found in Appendix A. The second step is updating the parameter interval according to the derivative sign of E(α) at the bisection point (a + b)/2. Repeat the step 1 and step 2 until iterative converges. Convergence criterion is b − a < ε and the ε is the threshold value close to zero, which can be set according to the desired precision. However, to directly calculate the first-order derivative of E(α) may be time consuming because of the complexity of the formula. Notice that the sign of the first-order derivative of E(α) is enough for the parameter interval updating.
To get this sign, set a differential increment ∆α on both sides of the bisection point, and calculate the E(α) at the point α + ∆α and α − ∆α. Then the derivative sign of E(α) at the bisection point can be obtained by a subtraction between E(α + ∆α) and E(α − ∆α). The details of the proposed MMEA for GMTIm are summarized in Figure 3.

Spectrum Rotation for Geometry Correction
After the parameterα is estimated by MMEA, the corresponding phase compensation function can be constructed as where the estimatedα is equal to actual α in theory. Multiply S 1_ROI f r , f a with Equation (19), the compensated signal becomes Equation (20) The refocused GMT image is obtained after adopting Stolt interpolation and 2-D IFT to Equation (20). However, Stolt interpolation need to select spectrum support region which is rectangle, and with the increase of the squint angle, the size of spectrum support region decreases. As is shown in Figure 4a, the spectrum support region reaches the minimum size when the squint angle is 45 degrees, and the size of spectrum support region increases when the squint angle continues to increase. In squint mode SAR, the decrease of spectrum support region may influence the image quality. Besides, the image exists geometry distortion after directly adopting Stolt interpolation and 2-D IFT to Equation (20) in squint mode. gence criterion is ε − < b a and the ε is the threshold value close to zero, which can be set according to the desired precision. However, to directly calculate the first-order derivative of ( ) α E may be time consuming because of the complexity of the formula. Notice that the sign of the first-order derivative of ( ) α E is enough for the parameter interval updating. To get this sign, set a differential increment α Δ on both sides of the bisection point, and calculate the ( ) α E at the point α α + Δ and α α − Δ . Then the derivative sign of ( ) α E at the bisection point can be obtained by a subtraction between ( ) α α . The details of the proposed MMEA for GMTIm are summarized in Figure   3. ?

Spectrum Rotation for Geometry Correction
After the parameter α is estimated by MMEA, the corresponding phase compensation function can be constructed as where the estimated α is equal to actual α in theory. Multiply The refocused GMT image is obtained after adopting Stolt interpolation and 2-D IFT to Equation (20). However, Stolt interpolation need to select spectrum support region which is rectangle, and with the increase of the squint angle, the size of spectrum support region decreases. As is shown in Figure 4a, the spectrum support region reaches the minimum size when the squint angle is 45 degrees, and the size of spectrum support region increases when the squint angle continues to increase. In squint mode SAR, the decrease of spectrum support region may influence the image quality. Besides, the image exists geometry distortion after directly adopting Stolt interpolation and 2-D IFT to Equation (20) in squint mode. Figure 4. (a) Spectrum support region selection before the spectrum rotation; (b) Spectrum support region selection after the spectrum rotation.
To get the maximum spectrum support region and remove the geometry distortion effect of GMT imaging, a spectrum rotation method is adopted. In this method, the rotation angle is the equivalent angle. As is shown in Figure 4b, the squinted 2-D spectrum is rotated into the side-looking form. For simplicity, set β α In ordinary spectrum rotation method aiming at the stationary scene, the rotation To get the maximum spectrum support region and remove the geometry distortion effect of GMT imaging, a spectrum rotation method is adopted. In this method, the rotation angle is the equivalent angle. As is shown in Figure 4b, the squinted 2-D spectrum is rotated into the side-looking form. For simplicity, set β =α − 1/v 2 . In Figure 4, F x = c β f a /2, . In ordinary spectrum rotation method aiming at the stationary scene, the rotation angle is the squint angle. However, for GMTIm application, the rotation angle is not simply equal to the squint angel because both of the GMT and SAR platform have motion. Different from the traditional squint angle spectrum rotation (SASR) method, this method is called equivalent squint angle spectrum rotation (ESASR) method. One key point of ESASR is that how to calculate the rotation angle which is relevant to the shift of the Doppler center. To obtain the Doppler center, the imaging geometry model of GMT is analyzed firstly, as is shown in Figure 5. In squinted SAR, the Doppler center shift f dc of the GMT echo consists of two parts, including the Doppler center shift f dc1 caused by the platform moving speed and the Doppler center shift f dc2 caused by the GMT velocity. Different from the traditional squint angle spectrum rotation (SASR) method, this method is called equivalent squint angle spectrum rotation (ESASR) method. One key point of ESASR is that how to calculate the rotation angle which is relevant to the shift of the Doppler center. To obtain the Doppler center, the imaging geometry model of GMT is analyzed firstly, as is shown in Figure 5. In squinted SAR, the Doppler center shift dc f of the GMT echo consists of two parts, including the Doppler center shift where θˆ is the angle shown in Figure 4b, which can be called equivalent squint angle. When the GTM velocity is zero, θ θ = , which means the rotation angle is equal to the squint angle. This is the condition of the stationary scene mentioned before. When the GMT has the velocity, the rotation angle is the equivalent squint angle θˆ. The Doppler center is obtained by Doppler center estimation method from the echo, rather than calculated by = − The spectrum rotation is achieved by the azimuth interpolation, the interpolation factor is ( ) where β α = − 2 1 v . Additionally, after the spectrum rotation, Equation (20) becomes From Figure 5 we know that f dc1 = 2v sin θ/λ and f dc2 = 2(v x sin θ + v r cos θ)/λ, the relationship of f dc1 , f dc2 and f dc is f dc = f dc1 − f dc2 . Then in analogy with the stationary scene rotation angle, the rotation angle of GMTIm iŝ whereθ is the angle shown in Figure 4b, which can be called equivalent squint angle. When the GTM velocity is zero,θ = θ, which means the rotation angle is equal to the squint angle. This is the condition of the stationary scene mentioned before. When the GMT has the velocity, the rotation angle is the equivalent squint angleθ. The Doppler center is obtained by Doppler center estimation method from the echo, rather than calculated by f dc = f dc1 − f dc2 because the GMT velocity is unknown. To sum up, the spectrum rotation factor is The spectrum rotation is achieved by the azimuth interpolation, the interpolation factor is where β =α − 1/v 2 . Additionally, after the spectrum rotation, Equation (20) becomes where After the spectrum rotation, adopt Stolt interpolation on Equation (24) with the β which containing the estimated and correct parameterα to eliminate the residual phase, the interpolation factor is after the interpolation, Equation (24) becomes Equation (28) S after taking 2-D IFT of Equation (26), the refocused GMT image with no geometry distortion is obtained.

Sparsity Enhancement for Sidelobe Suppression
Notice that the GMT image is usually sparse in space domain. If the GMT image is constrained with sparsity and the 2-D IFT is replaced by sparsity reconstruction, the sidelobe can be suppressed further. The process of reconstructing the GMT sparse result is represented as a L1 norm optimization problem where F a and F r are the azimuth and range FT matrix, respectively. X is the GMT sparse image. This optimization problem can be effectively solved by iterative threshold algorithm (ITA). The GMT sparse result at every iteration is where F −1 r and F −1 a are the range and azimuth IFT matrix, respectively. µ is the parameter related to iteration convergency. E denotes the thresholding operator For the full sampled data, the GMT sparse result is obtained through only one iteration. Suppose that the sparse result X is zero, then the algorithm of the sparsity enhancement is shown as Figure 6

Results
In this section, simulated data and measured ship data are used to give experiments of the proposed GMTIm algorithm. Section 4.1 is the experiment of a simulated scene. The simulated scene consists of two point targets. The process of GMTIm algorithm is displayed, including the conventional imaging results by RMA, focused GMT image without or with the spectrum rotation and their respective sparse enhancement image. Section 4.2 is the experiment of measured ship data collected by spaceborne SAR platform. The data cannot be directly used because it is side-looking mode and the proposed algorithm is aimed at squint mode. The side-looking SAR ship data are reused to form the squint SAR ship data. The quality of the refocused GMT image is analyzed in both simulated data and ship data experiment, and the effectiveness of the proposed GMTIm algorithm is demonstrated.

GMTIm Based on Simulated Data
In this section, simulated point targets data are used to analyze the performance of the proposed algorithm. As is shown in Figure 7, the point targets scene is comprised of two scatterers: moving target M1 and stationary target S1. The range and azimuth velocity of

Results
In this section, simulated data and measured ship data are used to give experiments of the proposed GMTIm algorithm. Section 4.1 is the experiment of a simulated scene. The simulated scene consists of two point targets. The process of GMTIm algorithm is displayed, including the conventional imaging results by RMA, focused GMT image without or with the spectrum rotation and their respective sparse enhancement image. Section 4.2 is the experiment of measured ship data collected by spaceborne SAR platform. The data cannot be directly used because it is side-looking mode and the proposed algorithm is aimed at squint mode. The side-looking SAR ship data are reused to form the squint SAR ship data. The quality of the refocused GMT image is analyzed in both simulated data and ship data experiment, and the effectiveness of the proposed GMTIm algorithm is demonstrated.

GMTIm Based on Simulated Data
In this section, simulated point targets data are used to analyze the performance of the proposed algorithm. As is shown in Figure 7, the point targets scene is comprised of two scatterers: moving target M1 and stationary target S1. The range and azimuth velocity of M1 are v r = 18 m/s and v x = 10 m/s, respectively. The radar platform flies at the velocity of 150 m/s and the other parameters of radar are: the signal center frequency is 10 GHz, the bandwidth is 75 MHz, the pulse duration is 2.2 µs, and the pulse repetition frequency (PRF) is 3000 Hz, the scene center range is 5000 m and the squint angle are 45, 60, 75 degrees.

Results
In this section, simulated data and measured ship data are used to give experiments of the proposed GMTIm algorithm. Section 4.1 is the experiment of a simulated scene. The simulated scene consists of two point targets. The process of GMTIm algorithm is displayed, including the conventional imaging results by RMA, focused GMT image without or with the spectrum rotation and their respective sparse enhancement image. Section 4.2 is the experiment of measured ship data collected by spaceborne SAR platform. The data cannot be directly used because it is side-looking mode and the proposed algorithm is aimed at squint mode. The side-looking SAR ship data are reused to form the squint SAR ship data. The quality of the refocused GMT image is analyzed in both simulated data and ship data experiment, and the effectiveness of the proposed GMTIm algorithm is demonstrated.

GMTIm Based on Simulated Data
In this section, simulated point targets data are used to analyze the performance of the proposed algorithm. As is shown in Figure 7, the point targets scene is comprised of two scatterers: moving target M1 and stationary target S1. The range and azimuth velocity of   Figure 8 is the conventional RMA imaging results of point targets. The moving target M1 is defocused compared to S1 before the residual phase is compensated. The following procedure is based on the GMT ROI data which are indicated by the red dash box. According to the proposed algorithm, the ROI data from conventional imaging result are firstly transformed into 2-D frequency domain and then turned into the MMEA parameter estimation step.
Remote Sens. 2021, 13, x FOR PEER REVIEW 13 of 20 Figure 8 is the conventional RMA imaging results of point targets. The moving target M1 is defocused compared to S1 before the residual phase is compensated. The following procedure is based on the GMT ROI data which are indicated by the red dash box. According to the proposed algorithm, the ROI data from conventional imaging result are firstly transformed into 2-D frequency domain and then turned into the MMEA parameter estimation step.  First step of MMEA is to obtain the initial parameter interval. In most of the practical situations, the GMT velocity is less than 30 m/s. As is shown in Figure 9a, according to The reason why MMEA cost such a few iterations is that the parameter interval shrinks to the half range in every iteration and the optimal parameter is guaranteed within the parameter interval all the time. If the initial parameter is set properly, IMEA can also converges fast. The problem for IMEA is that it is difficult to set the proper parameter due to the unknown GMT motion. The MMEA does not need to set an initial parameter but a parameter value range which is irrelevant to the actual GMT motion. Compensate the residual phase error of the GMT defocused image, we get the GMT focused image and sparse enhancement image with or without spectrum rotation, which First step of MMEA is to obtain the initial parameter interval. In most of the practical situations, the GMT velocity is less than 30 m/s. As is shown in Figure 9a, according to , the initial parameter interval is easily obtained, that is [3.003 × 10 −5 , 6.9444 × 10 −5 ]. The parameter iteration comparison graph of MMEA and IMEA is shown in Figure 9b, from which we can see MMEA and IMEA reach the same estimated parameter 5.35 × 10 −5 . However, IMEA iterated over 40 times to reach the final parameter, while the proposed MMEA iterated less than 10 times. The reason why MMEA cost such a few iterations is that the parameter interval shrinks to the half range in every iteration and the optimal parameter is guaranteed within the parameter interval all the time. If the initial parameter is set properly, IMEA can also converges fast. The problem for IMEA is that it is difficult to set the proper parameter due to the unknown GMT motion. The MMEA does not need to set an initial parameter but a parameter value range which is irrelevant to the actual GMT motion.
Remote Sens. 2021, 13, x FOR PEER REVIEW 13 of 20 Figure 8 is the conventional RMA imaging results of point targets. The moving target M1 is defocused compared to S1 before the residual phase is compensated. The following procedure is based on the GMT ROI data which are indicated by the red dash box. According to the proposed algorithm, the ROI data from conventional imaging result are firstly transformed into 2-D frequency domain and then turned into the MMEA parameter estimation step.  First step of MMEA is to obtain the initial parameter interval. In most of the practical situations, the GMT velocity is less than 30 m/s. As is shown in Figure 9a, according to The reason why MMEA cost such a few iterations is that the parameter interval shrinks to the half range in every iteration and the optimal parameter is guaranteed within the parameter interval all the time. If the initial parameter is set properly, IMEA can also converges fast. The problem for IMEA is that it is difficult to set the proper parameter due to the unknown GMT motion. The MMEA does not need to set an initial parameter but a parameter value range which is irrelevant to the actual GMT motion. Compensate the residual phase error of the GMT defocused image, we get the GMT focused image and sparse enhancement image with or without spectrum rotation, which Compensate the residual phase error of the GMT defocused image, we get the GMT focused image and sparse enhancement image with or without spectrum rotation, which is shown in Figure 10. In Figure 10a, the focused GMT image without spectrum rotation has severe geometry distortion that the point-like GMT turns into ellipse-like GMT. Accordingly, in Figure 10d, the sparse enhancement image of the line-like GMT image is also ellipse-like. As is shown in Figure 10b, although a part of the geometry distortion can be corrected by SASR, there still exists some distortion residual. As is shown in Figure 10c,f, the geometry distortion is almost corrected by ESASR and the corresponding sparse enhancement result has a low sidelobe.
ing. The spectrum rotation is needed to correct the geometry distortion. To furthe ate the performance of the spectrum rotation method, the imaging quality ind MMEA, MMEA+SASR, MMEA+ESASR are quantified as shown in Table 1. The indexes include peak side-lobe ratio (PSLR), integrated side-lobe (ISLR), and di angle (DA) which is denoted as θ DA . The θ DA is defined as the angle of the cro lobe, as is shown in Figure 10a-c. Evidently, compared with the MME MMEA+SASR, the result of MMEA+ESASR are more similar to the ideal value o ISLR, and especially the DA. It means MMEA+ESASR achieves precise focusing GMT and corrects a large proportion of the geometry distortion.

Range bins Range bins Range bins
Range bins Range bins Range bins   As is mentioned before, MMEA converges faster than IMEA in the aspect of refocusing. The spectrum rotation is needed to correct the geometry distortion. To further evaluate the performance of the spectrum rotation method, the imaging quality indexes of MMEA, MMEA+SASR, MMEA+ESASR are quantified as shown in Table 1. The quality indexes include peak side-lobe ratio (PSLR), integrated side-lobe (ISLR), and distortion angle (DA) which is denoted as θ DA . The θ DA is defined as the angle of the cross side-lobe, as is shown in Figure 10a-c. Evidently, compared with the MMEA and MMEA+SASR, the result of MMEA+ESASR are more similar to the ideal value of PSLR, ISLR, and especially the DA. It means MMEA+ESASR achieves precise focusing of the GMT and corrects a large proportion of the geometry distortion. The point target response in azimuth and range profile of the GMT image with ESASR and its sparse enhancement (SE) image are shown in Figure 11, where the red line denotes the point target response of GMT with ESASR and the green line denotes the ESASR-SE result. The ESASR-SE results have a narrower main lobe and lower side lobe than the refocused GMT image without sparse enhancement.
The point target response in azimuth and range profile of the GMT image with ESASR and its sparse enhancement (SE) image are shown in Figure 11, where the red line denotes the point target response of GMT with ESASR and the green line denotes the ESASR-SE result. The ESASR-SE results have a narrower main lobe and lower side lobe than the refocused GMT image without sparse enhancement.

Magnitude(dB)
Magnitude ( To demonstrate the proposed algorithm works in situation of higher squint angle, the following experiments is implemented when the squint angle are 60 or 75 degrees. As is shown in Figure 12, the focused GMT image is displayed after the MMEA. Before the spectrum rotation is used, the target geometry distortion becomes more severe with the increase of squint angle. The point-like target nearly turns into line-like. When the SASR is implemented, the target geometry distortion is relieved a little, while the ESASR can eliminates almost all the distortion. The following procedure is the sparce enhancement which is the same as part A. Table 2 is the comparison of distortion angle.  To demonstrate the proposed algorithm works in situation of higher squint angle, the following experiments is implemented when the squint angle are 60 or 75 degrees. As is shown in Figure 12, the focused GMT image is displayed after the MMEA. Before the spectrum rotation is used, the target geometry distortion becomes more severe with the increase of squint angle. The point-like target nearly turns into line-like. When the SASR is implemented, the target geometry distortion is relieved a little, while the ESASR can eliminates almost all the distortion. The following procedure is the sparce enhancement which is the same as part A. Table 2 is the comparison of distortion angle. The point target response in azimuth and range profile of the GMT image with ESASR and its sparse enhancement (SE) image are shown in Figure 11, where the red line denotes the point target response of GMT with ESASR and the green line denotes the ESASR-SE result. The ESASR-SE results have a narrower main lobe and lower side lobe than the refocused GMT image without sparse enhancement. To demonstrate the proposed algorithm works in situation of higher squint angle, the following experiments is implemented when the squint angle are 60 or 75 degrees. As is shown in Figure 12, the focused GMT image is displayed after the MMEA. Before the spectrum rotation is used, the target geometry distortion becomes more severe with the increase of squint angle. The point-like target nearly turns into line-like. When the SASR is implemented, the target geometry distortion is relieved a little, while the ESASR can eliminates almost all the distortion. The following procedure is the sparce enhancement which is the same as part A. Table 2 is the comparison of distortion angle.

GMTIm Based on GF-3 Data
The original SAR data of ship target are collected by a side-looking SAR platform GF-3 satellite [14]. To make the side-looking SAR data suitable for the squint SAR experiments, ship target data are used to calculate the scattering echoes of squint mode SAR. The motion parameter of the ship target is [v x , v r ] = [5 m/s,10 m/s]. The radar platform flies at a 10 km height with a speed of 150 m/s. The down-sight angle is 60 degrees. The other parameters of radar are: the signal center frequency is 10 GHz, the bandwidth is 150 MHz, the pulse duration is 2.2 µs, and the pulse repetition frequency (PRF) is 1275 Hz, the scene center range is 5000 m, the radar operates in squint-look mode with squint angle 50 degrees.
The ROI data extraction and parameter estimation process is shown in Figure 13. The GMT ROI data are indicated by the red dashed box in Figure 13a which is the conventional RMA imaging result. The parameter iteration process of MMEA and IMEA is made a contrast in Figure 13b, and the estimated parameter is 4.877 × 10 −5 . MMEA iterated only 9 times to reach the final parameter, which is less than that of IMEA. Figure 14a is the defocused ROI image. As is shown in Figure 14b, after the residual phase compensation operation, the focused image without spectrum rotation is obtained. There is some influence of the geometry distortion on the target recognition application. After the spectrum rotation, the geometry distortion of the GMT is corrected and the focused image is shown in Figure 14c. Even if there exists some orientation difference from the actual ship target, it is easier to accomplish the target recognition. The final result shown in Figure 14d is the sparse enhancement result whose sidelobe is suppressed and target characteristic is more distinct.

GMTIm Based on GF-3 Data
The original SAR data of ship target are collected by a side-looking SAR platform GF-3 satellite [14]. To make the side-looking SAR data suitable for the squint SAR experiments, ship target data are used to calculate the scattering echoes of squint mode SAR. The motion parameter of the ship target is [ x v , r v ] = [5 m/s,10 m/s]. The radar platform flies at a 10km height with a speed of 150 m/s. The down-sight angle is 60 degrees. The other parameters of radar are: the signal center frequency is 10GHz, the bandwidth is 150 MHz, the pulse duration is 2.2 μs , and the pulse repetition frequency (PRF) is 1275 Hz, the scene center range is 5000 m, the radar operates in squint-look mode with squint angle 50 degrees.
The ROI data extraction and parameter estimation process is shown in Figure 13. The GMT ROI data are indicated by the red dashed box in Figure 13a which is the conventional RMA imaging result. The parameter iteration process of MMEA and IMEA is made a contrast in Figure 13b, and the estimated parameter is − × 5 4.877 10 . MMEA iterated only 9 times to reach the final parameter, which is less than that of IMEA. Figure 14a is the defocused ROI image. As is shown in Figure 14b, after the residual phase compensation operation, the focused image without spectrum rotation is obtained. There is some influence of the geometry distortion on the target recognition application. After the spectrum rotation, the geometry distortion of the GMT is corrected and the focused image is shown in Figure 14c. Even if there exists some orientation difference from the actual ship target, it is easier to accomplish the target recognition. The final result shown in Figure 14d is the sparse enhancement result whose sidelobe is suppressed and target characteristic is more distinct.

Conclusions
In this paper, a GMTIm method for the squint mode SAR was proposed, which had the following steps: ROI data extraction, parameter estimation by MMEA, residual phase compensation, spectrum rotation, and sparse enhancement. The ROI data of the GMT were extracted from the preliminary SAR image of the RMA. In squint SAR, the RMA signal model of the GMT was derived, and the accurate residual phase compensation was achieved because of the accurate slant range model of the RMA. The ROI data were useful for decreasing the processing data volume and restraining the clutter at the same time. The parameter estimation step took advantage of the idea of bisection to boost up the convergence speed of the minimum entropy algorithm. The new efficient minimum entropy algorithm is called MMEA. After the residual phase compensation by the compensation function with the estimated parameter, the refocused GMT image was obtained. Spectrum rotation method ESASR corrected the geometry distortion of the GMT image, which was helpful for the target recognition application. The final sparse enhancement step fully made use of the GMT sparsity to further suppress the GMT image sidelobe. Experiment results validate the effectiveness of the proposed method on GMTIm in squint mode SAR. The simulated targets move with a constant velocity, which may not conform the physical truth that the targets have a complicated motion. This will be our future work.
Author Contributions: Conceptualization and methodology, S.X. and J.N.; software, S.X.; resources, Q.Z and Y.L.; writing-review and editing, S.X., J.N., Q.Z. and L.Y. All authors have read and agreed to the published version of the manuscript.  we can obtain the first-order derivative of the entropy function with respect to by substituting Equations (A5)-(A9) into Equation (A1).