An Efficient Method for Ground Maneuvering Target Refocusing and Motion Parameter Estimation Based on DPT–KT–MFP

: The image of ground maneuvering targets may be defocused due to the Doppler ambiguity, high-order range migration (RM), and Doppler frequency migration (DFM) caused by the target’s complex motions in a synthetic aperture radar (SAR) system. To settle these problems, an efficient algorithm based on discrete polynomial-phase transform (DPT), keystone transform (KT), and matched filtering processing (MFP) is presented for ground maneuvering target refocusing and motion parameter estimation in this paper. Firstly, the DPT is applied to transform the cubic phase into the quadratic phase and simultaneously eliminate the quadratic RM, cubic RM, and quadratic DFM. Furthermore, the Doppler ambiguity containing Doppler center blur and Doppler spectrum ambiguity is also dealt with effectively by introducing a very small fixed lay time after DPT operation. Then, the KT is performed to correct the linear RM. After that, the matched filtering function related to the target’s equivalent third-order coefficient is constructed to compensate for the residual linear DFM in the range-time and slow-time domain. Lastly, a well-refocused image of the maneuvering target can be acquired, and the target’s motion parameters can be estimated effectively. The proposed algorithm has high computational efficiency and possesses favorable refocusing performance and motion parameter estimation precision. Simulation and real data processing results prove the effectiveness of the presented algorithm.


Introduction
Synthetic aperture radar (SAR), which can work on all time scales and in all weather conditions, is widely applied in the civil and military fields [1][2][3][4][5][6][7][8][9]. In practical applications, in order to improve the focusing performance and acquire a high-resolution SAR image, the synthetic aperture time is usually long. However, within the long integration time, the ground maneuvering target may be severely defocused because of the high-order range migration (RM) and Doppler frequency migration (DFM) induced by the target's complex motions. As a consequence, in order to obtain good maneuvering target imaging performance, the effects of RM and DFM should be first eliminated.
The maneuvering target's cross-track velocity may cause linear RM. The keystone transform (KT) [10][11][12][13][14] was presented to correct the linear RM by rescaling the azimuth slow-time axis. The axis rotation moving target detection (AR-MTD) [15] was developed to eliminate the linear RM by axis rotation. In [16][17][18], Xu et al. proposed the radon Fourier transform (RFT), which jointly searches the target's range and velocity to eliminate the linear RM. Furthermore, the frequency-domain deramp-keystone transform (FDDKT) [19] was presented to refocus the maneuvering target without any parameter searching procedures. Nevertheless, for a ground maneuvering target with along-track velocity and cross-track acceleration, the abovementioned methods are invalid since they cannot eliminate the quadratic RM and linear DFM caused by the along-track velocity and cross-track acceleration of maneuvering targets.
In order to deal with the quadratic RM and linear DFM, several valid algorithms were recently presented. The second-order KT (SOKT) [20][21][22] was proposed to remove quadratic RM via time-scale transform. The radon fractional Fourier transform (RFRFT) [23] was developed to correct the linear and quadratic RM by searching the velocity and acceleration of the maneuvering target and then integrating the target energy via FRFT. The second-order Wigner-Ville distribution (SOWVD) method [24] was proposed to correct the RM via KT and a compensation function related to SAR platform velocity, before focusing the maneuvering target via SOWVD. Tian et al. proposed an algorithm based on a symmetric autocorrelation function and window Lv transform (SAF-WLVT) [25] to simultaneously remove the linear and quadratic RM via SAF operation, before effectively integrating the target energy via WLVT. In addition, the symmetric autocorrelation function and scaled Fourier transform (SAF-SFT) method [26] and the deramp SOKT phase difference (deramp-SOKT-PD) method [27] were presented to eliminate the quadratic RM and linear DFM without any motion parameter searching procedures. However, when a ground target maneuvers with an along-track acceleration, cubic RM and quadratic DFM may occur in a high-resolution SAR system. Consequently, the focusing performance of the aforementioned algorithms may deteriorate dramatically because they ignore the effects of cubic RM and quadratic DFM. In addition, the focusing performance may further deteriorate due to the Doppler ambiguity containing the Doppler center blur and Doppler spectrum ambiguity caused by the target's fast maneuvering motions.
To address the abovementioned problems, the generalized radon Fourier transform (GRFT) [28] was proposed to eliminate the high-order RM and DFM, as well as deal with the Doppler ambiguity. Although this method has optimal coherent accumulation performance and parameter estimation precision, the computational complexity of GRFT on the basis of four-dimensional (4D) parameter search is so huge that it is difficult to implement in practice. The KT-based method was proposed in [29], which first employs the KT to eliminate the linear RM, and then eliminates the residual RM and DFM by jointly searching the acceleration, jerk, and Doppler ambiguity number. This method can acquire similar target detection capability and parameter estimation precision to the GRFT method, but the focusing performance of the proposed method may deteriorate severely when the Doppler spectrum of a maneuvering target, whose bandwidth is greater than half of the pulse repetition frequency (PRF), is not entirely in one PRF band. Moreover, the computational burden of the KT-based method proposed in [29] is still heavy due to a 3D parameters search. In [30], Li et al. presented the generalized KT and generalized dechirp process (GKTGDP) algorithm to achieve target detection, which corrected the RM by performing high-order KT three times and searching the Doppler ambiguity number. However, this algorithm may suffer from great performance losses because it ignores the effect of Doppler spectrum ambiguity. Huang et al. proposed an algorithm based on SOKT and Houghgeneralized high-order ambiguity function (SOKT-GHHAF) [31], which estimates the target's motion parameters by introducing high-order transform kernels and eliminating the Doppler spectrum ambiguity via a preprocessing step (shifting target Doppler center with half of the PRF). Nevertheless, the preprocessing is useless if the bandwidth of the target Doppler spectrum is greater than half of the PRF. The improved 2D matched filter algorithm [32] was presented to accumulate the weak target using a 2D grid search in the range-frequency and azimuth-Doppler domain. However, under the condition of a small time-bandwidth product, this algorithm may be inappropriate.
It is worth noting that the abovementioned methods [28][29][30][31][32] have a high computational burden due to multidimensional motion parameter search. To decrease the computational cost, several computationally efficient algorithms were recently presented. A method [33] based on time reversing transform, second-order keystone transform, and Lv distribution (TRT-SOKT-LVD) was developed to refocus the maneuvering only using a 1D parameter search. However, the bandwidth of the target spectrum is doubled, which may aggravate the Doppler spectrum ambiguity. The adjacent cross-correlation function iteratively (IACCF) method [34] was presented to accumulation the target energy via two ACCF operations. Although the IAACF can be quickly implemented via fast Fourier transform (FFT) and inverse FFT (IFFT), its integration performance is poor because the nonlinear transform is applied twice in the IACCF method. In [35], a fast algorithm based on auto-correlation function, polynomial-phase transform, and LVD (ACF-PPT-LVD) was developed to remove the RM and DFM by decreasing the phase order. Nevertheless, the application of a high-order nonlinear transform resulted in a serious deterioration of focusing performance.
Motivated by previous researches, we propose an efficient method for ground maneuvering target refocusing based on discrete polynomial-phase transform (DPT), KT, and matched filtering processing (MFP), i.e., DPT-KT-MFP. Firstly, the DPT is performed to simultaneously eliminate the quadratic RM, cubic RM, and quadratic DFM. Furthermore, the effects of Doppler ambiguity containing Doppler center blur and Doppler spectrum ambiguity are removed by introducing a very small fixed lay time after DPT operation. After that, the KT is applied to correct the linear RM. Then, a matched filtering function is constructed to compensate for the residual linear DFM caused by the equivalent third-order coefficient in the range-time and slow-time domain. After matching the proper equivalent third-order coefficient in the 1D parameter searching space and eliminating the effect of linear DFM, a well-refocused image of the maneuvering target can be acquired, and then the target's motion parameters can be effectively estimated. The contributions of this paper are as follows: 1. The presented DPT-KT-MFP method can effectively remove the high-order RM and DFM caused by a ground maneuvering target's complex motions. 2. The presented method can effectively deal with the Doppler ambiguity containing Doppler center blur and Doppler spectrum ambiguity. 3. The proposed algorithm has high compensational efficiency since only the target's equivalent third-order coefficient, which is usually very small for ground maneuvering targets, requires to be searched. 4. The proposed algorithm has high motion parameter estimation accuracy due to it avoiding the error propagation effect.
The remainder of this paper is organized as follows: in Section 2, the mathematical model of the received signal is established. In Section 3, the proposed refocusing algorithm is introduced. Section 4 presents some simulation experiments to evaluate the refocusing performance of the proposed method. At last, the work is concluded in Section 5.

Signal Model
The motion geometry relationship between the SAR platform with the constant speed v and the ground moving target on the 2D slant plane is given in Figure 1, where the radar works in the side-looking mode. The maneuvering target moves from a to b within the synthetic aperture time a T . c v , c a , a v , and a a are the target's cross-track velocity, cross-track acceleration, along-track velocity, and along-track acceleration, respectively, whereas 0 R is the target's nearest slant range. The instantaneous slant range between the radar and the maneuvering target is described as follows [31]: denote the equivalent first-order coefficient, second-order coefficient, and third-order coefficient of Assuming that the radar transmits a linear frequency modulated (LFM) signal, the received baseband signal can be described as follows [33,34]: where 1 A denotes the amplitude after range compression, whereas B and  indicate the bandwidth and wavelength of transmitted signal, respectively. In general, because of the maneuvering target's high-speed motion and low PRF, there is a velocity ambiguity problem. Therefore, the target's equivalent speed can be described as where 2 A denotes the echo amplitude after range FFT, and f denotes the range frequency variable corresponding to t .
Considering that (5) can be rewritten as From Equation (6), we can find that the term  Figure 2, the target spectrum distributions can be classified as the following three cases: (1) target spectrum is overall in a PRF band [13]; (2) target spectrum occupies two neighboring PRF bands (target spectrum bandwidth is less than one PRF) [36]; (3) target spectrum spans over several PRF bands (target spectrum bandwidth is larger than one PRF) [37]. Obviously, the target Doppler spectrum is split in cases 2 and 3 since it is not overall in a PRF band, thereby bringing about Doppler spectrum ambiguity. To solve these problems, we develop an efficient method for ground maneuvering target refocusing based on DPT-KT-MFP in this paper.

DPT-KT-MFP Algorithm with Single Target
The DPT [38,39] can be applied to decrease the target phase order, which can be stated as where  is the conjugate operation, and 0  is the fixed delay time, which is usually very small. Substituting Equation (6) into Equation (7) yields From Equation (8), we can find that the order of m t is reduced to second order, which means that the quadratic DFM and cubic RM caused by the third-order coefficient  are so small that the azimuth Doppler spectrum of the maneuvering target is entirely distributed in the non-ambiguous Doppler area. Thus, the problem of Doppler center blur and Doppler spectrum ambiguity is effectively addressed. However, the linear RM caused by 2 c and linear DFM induced by 3 c still exist after the DPT procedure, which can deteriorate the focusing performance.
In [14,24,39], several algorithms were proposed to correct the RM related to secondorder coefficient 2 c by constructing a compensation function. However, constructing the compensation function needs some prior knowledge, such as platform velocity or the target's nearest slant range, which cannot always be accurately obtained. Additionally, the compensation error of RM caused by these methods [14,24,39] may be larger than a range resolution cell in the high-resolution case, which may result in a serious refocusing performance loss.
Consequently, in this paper, the KT was applied to absolutely remove the linear RM without any prior knowledge. The KT can be stated as follows [10][11][12][13][14]24,29]: where n t indicates the new slow-time variable. Substituting Equation (9) into Equation am a n n a n c Under the narrowband condition, From Equation (11), one can see that the coupling between the term n t  and f is removed, which means that the linear RM is effectively corrected. Moreover, as for the ground maneuvering target, the   After performing the range IFFT on Equation (11), we obtain       where 3 A represents the echo amplitude after IFFT. From Equation (12), we can find that the RM containing linear, quadratic, and cubic RM is effectively removed via the previous operation. However, due to the existence of the linear DFM caused by third-order coefficient 3 c , the maneuvering target is still seriously defocused if we directly apply the azimuth FFT to accumulate the target energy. The previous methods proposed in [24,31,33,39] assumed that the target signal is an LFM signal after correcting the RM, and they utilized time-frequency analysis algorithms to remove the linear DFM. Nevertheless, these methods need to search the range cells for extracting the target energy, which greatly increases the computational complexity. Actually, the maximum offset of linear DFM induced by 3 c , which is equal to c only spans across several Doppler frequency resolution cells. Therefore, the searching third-order coefficient is much smaller than the searching range cells. Hence, we constructed a matched filtering function to remove the linear DFM, which can be expressed as  (14), we can find that the Doppler frequency related to 3 a is  is equal to a Doppler frequency resolution cell in one coherent process interval [23]. Therefore, the 3 a  can be described as When 3 a is equal to 3 c , multiplying Equation (13) by Equation (12) and performing the azimuth FFT yields       c where   FFT n t denotes the FFT operation, 4 A is the signal amplitude after FFT operation, and n f indicates the Doppler frequency variable corresponding to n t . From Equations (13) and (15), it can be seen that, when the third-order coefficient 3 c is matched exactly, the maneuvering target is well refocused at According to the aforementioned analyses, we know that the presented algorithm can effectively remove the effects of Doppler center blur, Doppler spectrum ambiguity, high-order RM, and DFM, as well as refocus the maneuvering target. Figure 3 gives the flowchart of the presented method.  Below, we give an example to demonstrate that the presented algorithm can effectively refocus a ground maneuvering target.
Example 1: Suppose a ground maneuvering target exists in the observation area. The radar system parameters are listed in Table 1 Figure 4a gives the range compression result. Obviously, serious RM containing linear, quadratic, and cubic RM occurs. The target Doppler spectrum is given in Figure 4b. It can be seen that the target spectrum occupies three PRF bands, as shown in Figure 2c, which means that serious Doppler spectrum ambiguity appears. The target trajectory after DPT operation is shown in Figure 4c, from which we can find that an oblique straight line is present, which means that the quadratic and cubic RM is effectively removed. The target spectrum after DPT operation is given in Figure 4d. We can observe that the target spectrum is overall distributed in one PRF band, which means that the problem of Doppler spectrum ambiguity is addressed. Additionally, the Doppler center blur is also eliminated after DPT operation, since the Doppler center is changed from 320 HZ ( 2 am n   ) to 166.2 HZ ( 0 am n  ). Figure 4e shows the result after KT operation. We can find that target trajectory becomes one straight line along the slow-time axis, which means that the residual linear RM is effectively corrected. Figure 4f shows the searching result of the third-order coefficient. Obviously, according to the peak position, 3 c can be estimated as 0.2624 m/s 3 . After correcting the linear DFM using the estimated third-order coefficient, we can observe that a distinct peak appears in Figure 4g, which indicates that the presented algorithm can effectively remove the Doppler center blur, Doppler spectrum ambiguity, highorder RM, and DFM, as well as refocus the maneuvering target. On the basis of peak detection technology [40], the first-and second-order coefficients can be estimated as follows: 1

DPT-KT-MFP Algorithm with Multiple Targets
With respect to refocusing multiple maneuvering targets, cross-term interference is inevitably introduced because of the existence of nonlinear transform in our proposed method. In this subsection, the effects of cross-terms are analyzed. Suppose that K maneuvering targets exist in the observation area; then, the range compression signal in m f t  domain can be stated as rect exp From Equations (19) and (20), it can be seen that, with respect to the auto-terms after DPT, the cubic RM and quadratic DFM are effectively removed. Furthermore, similar to refocusing the single target, the quadratic RM caused by are still coupled with f . Consequently, after carrying out the KT and MFP operations as described in Section 3.1, the energy of the auto-terms is well coherently accumulated, while the energy of the cross-terms cannot be integrated effectively because of the high-order RM and DFM. In addition, the proposed method can also refocus the maneuvering targets in the case of close multiple targets, since the 0,i R in auto-terms is effectively removed after DPT operation. On the basis of the aforementioned analyses, we can conclude that the presented method is able to remove the effects of cross-terms and refocus the ground maneuvering targets in the case of multiple targets. Next, several simulation experiments are presented to validate how the presented method suppresses the cross-terms and focuses the auto-terms. Example 2: Assume two maneuvering targets, represented by T1 and T2, exist in the observation area. The radar system parameters are the same as those in Example 1. In addition, the motion parameters of the two maneuvering targets are set as follows:  Figure 5 gives the simulation results. Figure 5a gives the target trajectories after range compression. Obviously, the targets' energy spans across several range cells during the coherent integration time. After DPT and KT operations, the RM correction result is given in Figure 5b, from which it can be seen that the energy of T1 and T2 is concentrated in the same range cells, while the energy of cross-terms still spans across several range cells. Figure 5c shows the searching result of the third-order coefficient. It can be seen that two distinct peaks exist. According to the peak positions, the third-order coefficients of T1 and T2 can be estimated as follows: . After compensating for the linear DFM using the estimated third-order coefficients, Figure 5d,e give the final focusing results of T1 and T2, respectively. Obviously, two well-focused peaks appear in Figure 5d,e, which demonstrates that the propose algorithm can effectively eliminate the effects of cross-terms and integrate the energy of auto-terms. According to the peak position, the remaining motion parameters of T1 and T2 can be estimated as follows: 1,1 31.9975 c   and 2,1 11.5162 c  for T1; 1,2 36.0019 c  and 2,2 9.9589 c  for T2.

Computational Complexity Analysis
Next, the computational cost of GRFT [28], SOKT-GHHAF [31], TRT-SOKT-LVD [33], and the proposed method was analyzed. Suppose N , M , and M  represent the numbers of range cells, echo pulses, and lag samples, respectively, and the number of searching first-, second-, and third-order coefficients is the same, denoted as 1 N . The computational complexity of GRFT based on 4D brute-force search is approximately   The procedures of the TRT-SOKT-LVD mainly include a TRT operation, whose computational cost is approxi- The computational complexity of the abovementioned four algorithms is listed in Table 2. SOKT-GHHAF, the presented method has a much lower computational burden. Compared with the TRT-SOKT-LVD, the proposed algorithm has a similar computational burden under the above-assumed condition. Actually, with respect to the ground maneuvering target, the number of searching range cells N is usually far greater than the number of searching third-order coefficients 1 N . Consequently, the presented algorithm has a higher computational efficiency than the TRT-SOKT-LVD algorithm.
Below, the computational time of these four methods is compared, where the target's motion parameters and radar parameters are the same as those in Example 1 (the processor of the computer was a Intel (R) Core (TM) i7-4790 central processing unit (CPU) @3.60 GHz; the random-access memory (RAM) was 8 GB; the software was MATLAB 2014a). Table 3 shows the computational time of these methods, from which it can be seen that the presented algorithm is more efficient in comparison with the other three algorithms. Table 2. Comparison of computational complexity. GRFT, generalized radon Fourier transform; SOKT-GHHAF, second-order KT and Hough-generalized high-order ambiguity function; TRT-SOKT-LVD, time reversing transform, second-order KT, and Lv distribution.

Methods
Computational Complexity

Analysis and Discussion of Simulation Results
In this subsection, we adopt a simulated ideal point target to validate the effectiveness of the presented DPT-KT-MFP algorithm under the Gaussian noise environment. The target's motion parameters and radar parameters are the same as those in Example 1. In order to better show the superior performance of the presented DPT-KT-MFP algorithm, GRFT [28], SOKT-GHHAF [31], TRT-SOKT-LVD [33], and moving target detection (MTD) are performed for comparison.

Maneuvering Target Integration Performance
In Figure 6, we first evaluate the integration performance of the presented algorithm for a maneuvering target. The signal-to-noise ratio (SNR) after range compression was set as 6dB . It is worth noting that, on the basis of the target's motion parameters and radar system parameters set above, Doppler ambiguity (containing Doppler center blur and Doppler spectrum ambiguity) occurs. The target trajectory after range compression is shown in Figure 6a, from which we can find that the target energy is distributed in multiple range cells. Figure 6b gives the direct KT result. One can see that the trajectory of the maneuvering target is split into three parts since its Doppler spectrum occupies three PRF bands. Figure 6c gives the result after DPT and KT, from which noise is removed to clearly show the target trajectory. We can observe that the target energy is distributed in one range cell, since the RM, Doppler center blur, and Doppler spectrum ambiguity are eliminated after DPT and KT operations. The refocusing result of the presented DPT-KT-MFP algorithm is given in Figure 6d. It can be seen that the target energy is well coherently integrated as a distinct peak. Figure 6e gives the MTD result. One can see that the MTD method cannot effectively accumulate the target energy because of the high-order RM and DFM caused by the complex motions of the maneuvering target. Figure 6f gives the focusing result of SOKT-GHHAF. Obviously, the SOKT-GHHAF method is invalid because of the Doppler spectrum ambiguity and the application of the incoherent integration algorithm. The result after TRT and SOKT is shown in Figure 6g. Similar to Figure 6c, the noise is removed to clearly show the target trajectory. We can observe that the target trajectory splits into five parts since the bandwidth of the target Doppler spectrum is doubled after TRT. Therefore, only a small fraction of target energy is accumulated by TRT-SOKT-LVD. Figure 6h gives the TRT-SOKT-LVD result. Obviously, the maneuvering target is still seriously defocused due to the effects of Doppler spectrum ambiguity. The GRFT result is shown in Figure 6i. We can observe that the integration performance of the GRFT algorithm is better than the presented algorithm. Nevertheless, the good accumulation performance is at the cost of a high computational burden.

Maneuvering Target Detection Ability
Below, we analyze the detection ability of MTD, SOKT-GHHAF, TRT-SOKT-LVD, GRFT, and the proposed method via Monte Carlo trials. The constant false alarm probability was set at 6 10  . After performing 500 Monte Carlo trials at every input SNR, Figure   7 shows the curves of detection probability varying with range compression SNRs. Obviously, the target detection performance of the presented algorithm was superior to that of the MTD, SOKT-GHHAF, and TRT-SOKT-LVD methods, since the effects of high-order RM and DFM, Doppler center blur, and Doppler spectrum ambiguity are effectively eliminated by the proposed algorithm. Compared with the GRFT method, the presented algorithm has a similar target detection ability but with a much lower computational cost when the input SNR is larger than 0dB . Unfortunately, the proposed method may suffer from a great detection performance loss in cases with low SNR, since the DPT is a nonlinear transform.

Motion Parameter Estimation Performance
Next, the motion parameter estimation performance of SOKT-GHHAF, TRT-SOKT-LVD, GRFT, and the presented algorithm was also compared via Monte Carlo trials, in which the input range-compressed SNRs varied from 15 B d  to 5 dB with a step size of 1dB . Moreover, at every input SNR, we carry out 500 Monte Carlo trials. Figure 8a-c show the root-mean-square errors (RMSEs) of the estimated first-, second-, and third-order coefficients varying with input SNRs, respectively. From Figure 8, one can see that the motion parameter estimation performance of SOKT-GHHAF method deteriorates severely, since the effects of Doppler ambiguity and cubic RM are ignored, and the HT, which is the incoherent integration operation, is used twice in the SOKT-GHHAF method. Moreover, SOKT-GHHAF has a huge computational burden due to the application of three 1D parameter searches. The TRT-SOKT-LVD method suffers from a great performance loss of parameter estimation, because it cannot remove the Doppler spectrum ambiguity. In comparison with the other three algorithms, the GRFT method can acquire the optimal motion parameter estimation performance, but it needs a 4D parameter search, which makes it difficult to implement in practice. When the input SNR is larger than 0dB , the motion parameter estimation performance of the presented algorithm is similar to that of the GRFT algorithm. Moreover, compared with the GRFT algorithm, the presented algorithm has a higher computational efficiency, since it only needs to search the third-order coefficient, whose value is usually very small for ground maneuvering targets.

Analysis and Discussion of Real Data Processing Results
In this subsection, we verify the effectiveness of the presented refocusing algorithm using real radar data recorded by an X-band airborne radar with a broadside mode. This radar has three apertures in the along-track direction, whereby one aperture transmits the signals and all three apertures receive the reflected echoes. The observed targets were cooperative moving trucks driving along the road. The SAR system parameters are listed in Table 4. Figure 9a gives the image of the ground maneuvering targets. We can observe that the maneuvering targets are submerged by strong ground clutter. In order to obtain a higher target signal-to-clutter-plus-noise-ratio (SCNR), we carried out the extended factored algorithm [41] to suppress the ground clutter. After clutter suppression, Figure 9b shows that several clear targets appear. However, the maneuvering targets are seriously smeared because of the RM and DFM. Below, we utilize a maneuvering target of interest marked by a red ellipse to demonstrate the effectiveness of the presented algorithm. The range compression result of the interested target is given in Figure 9c, from which we can find that the interested target energy is distributed in multiple range cells because of severe RM. The target trajectory after DPT and KT is shown in Figure 9d. it can be seen that the trajectory of the maneuvering target is parallel to the slow-time axis, which means that the RM is effectively eliminated. Figure 9e gives the focusing result of the presented algorithm. We can observe that a well-refocused peak appears since the high-order RM and DFM are absolutely removed, which proves the effectiveness of the presented DPT-KT-MFP algorithm.

Conclusions
In this paper, the DPT-KT-MFP algorithm was presented to address the refocusing and motion parameter estimation problems for ground maneuvering targets, involving Doppler ambiguity, high-order RM, and DFM during the coherent integration time. Obviously, the proposed algorithm can effectively remove the high-order RM and DFM induced by complex motions of ground maneuvering targets, as well as realize the target refocusing with low computational complexity. In addition, the Doppler ambiguity containing Doppler center blur and Doppler spectrum ambiguity can also be effectively dealt with due to the DPT operation. Furthermore, for refocusing of multiple targets, the presented algorithm is able to eliminate the effects of cross-terms in which high-order RM and DFM still exist. Compared with conventional algorithms which ignore the effects of Doppler ambiguity and the target's high-order motions, the presented algorithm can acquire better focusing performance and motion parameter estimation precision. Furthermore, compared with the GRFT method, the presented algorithm has higher computational efficiency and similar focusing performance in the case of a high SNR.  Data Availability Statement: Restrictions apply to the availability of these data. Data was obtained from China Academy of Space Technology (CAST) and are available from the authors with the permission of CAST.