High Speed Maneuvering Platform Squint TOPS SAR Imaging Based on Local Polar Coordinate and Angular Division

: This paper proposes an imaging algorithm for synthetic aperture radar (SAR) mounted on a high-speed maneuvering platform with squint terrain observation by progressive scan mode. To overcome the mismatch between range model and the signal after range walk correction, the range history is calculated in local polar format. The Doppler ambiguity is resolved by nonlinear derotation and zero-padding. The recovered signal is divided into several blocks in Doppler according to the angular division. Keystone transform is used to remove the space-variant range cell migration (RCM) components. Thus, the residual RCM terms can be compensated by a uniﬁed phase function. Frequency domain perturbation terms are introduced to correct the space-variant Doppler chirp rate term. The focusing parameters are calculated according to the scene center of each angular block and the signal of each block can be processed in parallel. The image of each block is focused in range-Doppler domain. After the geometric correction, the ﬁnal focused image can be obtained by directly combined the images of all angular blocks. Simulated SAR data has veriﬁed the effectiveness of the proposed algorithm.


Introduction
In the past a few decades, the technology of synthetic aperture radar (SAR) has been widely used in the field of geoscience and climate change research, environmental and earth system monitoring [1]. As an active sensor, radar transmits wideband signal for target detection and imaging [1][2][3][4]. It has an ability to get high resolution images in day and night and all weather conditions [1]. SAR has been mounted on numerous platforms, such as satellite, airplane, missiles, drones, and others. To get a well-focused SAR image, the corresponding imaging algorithms need to be developed according to the platform motion characteristics and the radar working mode. In most cases, the radar platform moves with a straight flight track and the radar works at side-looking mode. For airborne platforms, there are motion errors during a synthetic aperture which is usually caused by the air turbulence. Some motion compensation and autofocusing methods are proposed to handle this problem [5][6][7]. Recently, radar mounted on high-speed maneuvering platform becomes a research hotspot [8][9][10]. It has several features that are different from the common platform SAR imaging. One is the recorded data is usually part of the full synthetic aperture. The speed of platform is over Mach 1. For real-time imaging, the synthetic aperture time is usually less than a second. Thus, the SAR data is focused in range-Doppler domain for more efficiency. Another feature is the nonlinear flight In this paper, an imaging algorithm for high-speed maneuvering platform SAR with squint TOPS mode is proposed. The range history is calculated in local polar format coordinate. It is accurate and exactly matched with the signal after RWC without any variable substitution. To make the signal unambiguous in Doppler, the signal is first multiplied by the nonlinear derotation function. Then, the signal is recovered by zero-padding in Doppler and time domain phase compensation. The recovered signal is unambiguous both in time and Doppler domain. To weaken the influence of space-variation in third order term, the signal is divided into several blocks according to the angular division. Different from the existed methods, the data of each angular block can be processed in parallel. There is no need to combine the signal before the final focusing. The space-variant range envelop is removed by keystone transform (KT). Thus, the RCM can be corrected by unified phase function with reference to the scene center of each angular block. Because the space-variation of third order term can be neglected in a block, only the space-variant components in Doppler chirp rate term need to be considered. By introducing perturbation terms in Doppler, the space-variant components can be compensated in time domain. Then, the signal can be focused in Doppler after time domain deramp. The focusing parameters are all calculated with reference to the scene center of each angular block. The final focused image can be obtained by directly combining all the angular block images.
This paper is organized as follows. Section 2 gives the expressions of range history in local polar format. The accuracy of the proposed range model and space-variation properties are also analyzed. Section 3 is the proposed imaging algorithm based on angular division. Section 4 shows the experiments of simulated SAR data. The main processing steps have all been verified. Section 5 provides some discussions on the advantages of the proposed method. Section 6 conclude the paper.

Range Model in Local Polar Format Coordinate
The imaging geometry of maneuvering platform SAR with squint TOPS mode is shown in Figure 1. The platform moves along the curve ACD with velocity vector V = v x , v y , v z while the beam is scanning from Q s to Q e with a constant angular speed ω. B is an arbitrary point on the flight track corresponding to azimuth tine t a . C(0, 0, H) is the platform position at azimuth time t a = 0. During the following imaging processing, the signal spectrum is divided into several blocks corresponding to different angular areas. T is an arbitrary target within the angular area i whose center yaw angle is θ i . According to the local polar format definition in [14], the reference range of T is R i . So, the slant range of T can be expressed as where k n,i (R i , ∆θ) = (1/n!)(d n R T,i (t a ; R i , ∆θ)/dt n a )| t a =0 is the nth-order deviation at t a = 0 corresponding to angular area i. In each angular area i, an arbitrary target can be expressed with the reference range R i and an angular deviation ∆θ. The slant range terms k n,i are space-variant with R i and ∆θ which are the range and azimuth space-variation, respectively. The range-variation can be overcome by range blocking and chirp scaling methods [19]. For azimuth processing, more accurate phase compensations are needed, especially for the Doppler parameters.

Space-Variation of Doppler Parematers
The azimuthal space-variation of slant range terms in (1) need to be considered in the following processing.
is the linear term which will affect the azimuth focusing position. It is also the main component of RCM.
Some simulations are made to calculate the approximation errors with the system parameters in Table 1.  Figure 1. Imaging geometry of maneuvering platform SAR with TOPS mode.

Space-Variation of Doppler Parematers
The azimuthal space-variation of slant range terms in (1) need to be considered in the following processing. k n,i (R i , ∆θ)(n = 1, 2, 3, 4) are expanded at ∆θ = 0 as where k nm,i (R i ) = (1/m!)[d n k n,i (R i , ∆θ)/d∆θ]| ∆θ=0 . k 1,i (R i , ∆θ) is the linear term which will affect the azimuth focusing position. It is also the main component of RCM.
is the Doppler chirp rate term which directly affect the azimuth focusing quality. k 3,i (R i , ∆θ) and k 4,i (R i , ∆θ) are higher order terms. They usually affect the sidelobe level of a focused point. There are approximation errors in (2) when calculating the slant range terms. The corresponding phase errors can be expressed as Some simulations are made to calculate the approximation errors with the system parameters in Table 1. As the simulation results shown in Figure 2, the approximation errors in k 3,i (R i , ∆θ) and k 4,i (R i , ∆θ) can be neglected, for the maximum phase error is less than π/4. However, the expressions for k 1,i (R i , ∆θ) and k 2,i (R i , ∆θ) are not accurate enough, as shown in Figure 2a,b. In order to overcome this problem, angular division is used during the range cell migration correction (RCMC) and azimuth focusing. By dividing the whole steering angle into three parts, the approximation errors can be significantly reduced. As the simulation results shown in Figure 3, the approximation errors of k 1,i (R i , ∆θ) and k 2,i (R i , ∆θ) have been decreased to a value much less than π/4. So, the expressions in (2) are accurate enough for the following signal processing based on the angular division method.
Another advantage of the angular division method is the space-variation of k 3,i (R i , ∆θ) can also be neglected. For the existed space-variation correction methods [10,14], the spacevariation of k 2,i (R i , ∆θ) and k 3,i (R i , ∆θ) cannot be removed at the same time. Usually only the space-variation of k 2,i (R i , ∆θ) is compensated. However, to get a well-focused image, the space-variation of k 3,i (R i , ∆θ) also needs to be considered. As the simulation results shown in Figure 4a, the variation cannot be neglected, for the maximum phase value is greater than π/4. After angular division method, the space-variation can be neglected as the results shown in Figure 4b. Furthermore, there is no need for the azimuth signal combination. Therefore, the data of different angular blocks can be processed in parallel. In order to overcome this problem, angular division is used during the ran cell migration correction (RCMC) and azimuth focusing. By dividing the whole steer angle into three parts, the approximation errors can be significantly reduced. As the s ulation results shown in Figure 3, the approximation errors of ( ) have been decreased to a value much less than 4 π . So, the expressions in (2) are accur enough for the following signal processing based on the angular division method. Another advantage of the angular division method is the space-variation can also be neglected. For the existed space-variation correction metho [10,14], the space-variation of ( ) also needs to be considered. As simulation results shown in Figure 4a, the variation cannot be neglected, for the maxim phase value is greater than 4 π . After angular division method, the space-variation c be neglected as the results shown in Figure 4b. Furthermore, there is no need for the a muth signal combination. Therefore, the data of different angular blocks can be proces in parallel.

Imaging Algorithm
In this section, an imaging algorithm based on angular division is proposed. Figure  5 shows the flow charts of the algorithm. The SAR data is first multiplied by the unified RWC and nonlinear derotation function. Then, the signal is unambiguous in Doppler. After zero-padding and phase compensation, the signal is recovered in both time and Doppler domain. The signal is divided into several blocks in Doppler according to the angular division. For angular data block i, the data is compensated by the residual RWC function. The space-variant RCM components are removed by KT. The RCM is then corrected by a unified function. For the azimuth focusing, frequency perturbation terms are introduced to correct the space-variant Doppler chirp rate. All the focusing parameters are calculated according to the scene center of each angular block and the data can be processed in parallel. Each angular block i is focused in range-Doppler domain. The final focused SAR image can be obtained by the combination of all the angular block images after geometric correction.

Imaging Algorithm
In this section, an imaging algorithm based on angular division is proposed. Figure  5 shows the flow charts of the algorithm. The SAR data is first multiplied by the unified RWC and nonlinear derotation function. Then, the signal is unambiguous in Doppler. After zero-padding and phase compensation, the signal is recovered in both time and Doppler domain. The signal is divided into several blocks in Doppler according to the angular division. For angular data block i, the data is compensated by the residual RWC function. The space-variant RCM components are removed by KT. The RCM is then corrected by a unified function. For the azimuth focusing, frequency perturbation terms are introduced to correct the space-variant Doppler chirp rate. All the focusing parameters are calculated according to the scene center of each angular block and the data can be processed in parallel. Each angular block i is focused in range-Doppler domain. The final focused SAR image can be obtained by the combination of all the angular block images after geometric correction.

Imaging Algorithm
In this section, an imaging algorithm based on angular division is proposed. Figure 5 shows the flow charts of the algorithm. The SAR data is first multiplied by the unified RWC and nonlinear derotation function. Then, the signal is unambiguous in Doppler. After zero-padding and phase compensation, the signal is recovered in both time and Doppler domain. The signal is divided into several blocks in Doppler according to the angular division. For angular data block i, the data is compensated by the residual RWC function. The space-variant RCM components are removed by KT. The RCM is then corrected by a unified function. For the azimuth focusing, frequency perturbation terms are introduced to correct the space-variant Doppler chirp rate. All the focusing parameters are calculated according to the scene center of each angular block and the data can be processed in parallel. Each angular block i is focused in range-Doppler domain. The final focused SAR image can be obtained by the combination of all the angular block images after geometric correction.

Spectrum Recovering and Angular Division
The spectrum bandwidth of the TOPS SAR echo signal is constructed by two parts the resolution bandwidth and the beam steering bandwidth. As the time-frequency dia gram shown in Figure 6a, the resolution bandwidth is much smaller than the bandwidt caused by beam steering. The total Doppler bandwidth is much larger than the system PRF which will lead to Doppler ambiguity. The time-frequency line of beam steering ca be calculated by where ( )( )

Spectrum Recovering and Angular Division
By a range Fourier transform, the signal in range frequency domain can be expressed as where W r ( f r ) is the range envelop function in frequency domain. w a (t a ) is the azimuth beam window. φ cp ( f r , t a ; R, ∆θ) is the range-azimuth coupling phase term which can be expressed as where R T (t a ; R, ∆θ) is the slant range in (1) without angular division. The corresponding coefficients are k n (R, ∆θ) and k nm (R). The spectrum bandwidth of the TOPS SAR echo signal is constructed by two parts, the resolution bandwidth and the beam steering bandwidth. As the time-frequency diagram shown in Figure 6a, the resolution bandwidth is much smaller than the bandwidth caused by beam steering. The total Doppler bandwidth is much larger than the system PRF which will lead to Doppler ambiguity. The time-frequency line of beam steering can be calculated by where l n (R)(n = 0, 1, 2, 3) is corresponding to the slope of time-frequency line. l 0 (R) is the Doppler deviation caused by squint mode. It can be removed by unified RWC and range compression. The compensation function can be derived as l 1 (R) is the linear component caused by beam steering. Usually, this term can be compensated by linear derotation [17], as shown in Figure 6b. However, by nonlinear derotation, the required system PRF can be further reduced. The compensation function can be formed as (7) and (8) are both calculated with reference to the scene center corresponding to R re f . The time-frequency diagram of signal after nonlinear derotation is shown in Figure 6c. The time-frequency line in (6) becomes where l n (R) = l n (R) − l n R re f is the residual component after nonlinear derotation.  [17], as shown in Figure 6b. However, by nonlinear derotation, the required system PRF can be further reduced. The compensation function can be formed as (7) and (8)  The time-frequency diagram of signal after nonlinear derotation is shown in Figure 6c.
The time-frequency line in (6) becomes is the residual component after nonlinear derotation. It can be noticed that the signal is unambiguous in Doppler. Thus, the Doppler interval can be extended by zero-padding after an azimuth Fourier transform. In order to get the numbers of zero-padding, the Doppler bandwidth of echo signal and signal after nonlinear derotation need to be calculated. According to (6), the Doppler bandwidth can be derived as the beam steering bandwidth: Similarly, the signal Doppler bandwidth in Figure 6c can be expressed as It can be noticed that the signal is unambiguous in Doppler. Thus, the Doppler interval can be extended by zero-padding after an azimuth Fourier transform. In order to get the numbers of zero-padding, the Doppler bandwidth of echo signal and signal after nonlinear derotation need to be calculated. According to (6), the Doppler bandwidth can be derived as the beam steering bandwidth: Similarly, the signal Doppler bandwidth in Figure 6c can be expressed as So, the number of zero-padding can be derived with considering the sampling redundancy: Furthermore, the azimuth sampling frequency has been extended as The new azimuth time and frequency have changed to t a and f a , respectively. Then, the signal is transformed into azimuth time domain. A phase function is needed to compensate the influence caused by the nonlinear derotation in (8). The function can be formed as After phase compensation, the signal is unambiguous in Doppler, as shown in Figure 6d. Then, the signal can be divided into several blocks in Doppler according to the angular division. As shown in Figure 7a, the observation area is divided uniformly into n blocks. Each block has a same yaw angle. There is a mapping relation between the yaw angle and the Doppler, as shown in Figure 7b. For block i, the corresponding support region in Doppler can be easily found. To avoid resolution loss in the margin area, the selected data support region in Doppler is a bit larger than the calculated one. So, there is no need for the signal to be recombined before azimuth focusing. The signal after angular division can be processed in parallel without resolution loss. So, the number of zero-padding can be derived with considering the sampling redundancy: Furthermore, the azimuth sampling frequency has been extended as Then, the signal is transformed into azimuth time domain. A phase function is needed to compensate the influence caused by the nonlinear derotation in (8). The function can be formed as After phase compensation, the signal is unambiguous in Doppler, as shown in Figure 6d. Then, the signal can be divided into several blocks in Doppler according to the angular division. As shown in Figure 7a, the observation area is divided uniformly into n blocks. Each block has a same yaw angle. There is a mapping relation between the yaw angle and the Doppler, as shown in Figure 7b. For block i, the corresponding support region in Doppler can be easily found. To avoid resolution loss in the margin area, the selected data support region in Doppler is a bit larger than the calculated one. So, there is no need for the signal to be recombined before azimuth focusing. The signal after angular division can be processed in parallel without resolution loss.

Residual RWC and Space-Variant RCMC
After angular division and aforementioned processing steps, the phase in (5) becomes According to the simulation results in Figures 2 and 4, the third and the fourth order

Residual RWC and Space-Variant RCMC
After angular division and aforementioned processing steps, the phase in (5) becomes According to the simulation results in Figures 2 and 4, the third and the fourth order terms are processed as constant components. k 1,i (R i , ∆θ) = k 1,i (R i , ∆θ) − k 10 R re f is the residual linear component after the unified RWC in Equation (7). Figure 8 shows the RCM lines of targets T A , T B and T C which share a same reference range. There are still residual RWC components, as shown in Figure 8a. The signal after angular division needs to be processed with reference to the scene center of each angular block. So, the residual RWC need to be compensated by the function: Then, the linear term becomes k 1,i (R i , ∆θ) = k 1,i (R i , ∆θ) − k 10 R re f ,i . It can be noticed that there are still space-variant terms remained. As shown in Figure 8b, the RCM lines of T A and T C cannot be corrected by a same function with reference to T B . In order to remove the space-variant components, KT is applied by interpolating the signal into a new time interval. The mapping relation can be expressed as where t a is the azimuth time after KT. The corresponding azimuth frequency is f a . After KT, the phase terms in (15) becomes where φ cp,i f r , t a ; R i , ∆θ is the coupling phase term which can be decomposed into three terms. φ rcm,i f r , t a ; R i , ∆θ is the RCM component. φ src,i f r , t a ; R i , ∆θ is the secondary range compression (SRC) term which is caused by the KT. φ az,i t a ; R i , ∆θ is the azimuth phase term. It can be noticed that the linear RCM term including the space-variant components have been removed. The remained RCM and SRC terms can be compensated with reference to the scene center of each angular block. The compensation function can be formed as exp 4π c f c f 2 r k 20,i R re f ,i t a 2 + 3k 30,i R re f ,i t a 3 + 6k 40,i R re f ,i t a

(22)
After RCM and SRC compensation and a range inverse Fourier transform, the signal can be expressed as It can be seen in (23) that the targets sharing a same reference range have been focused into a same range cell. There is no mismatch between the range focused position and the azimuth phase [10] which is an advantage over the conventional RWC based methods.

Aizmuth Fousing
The phase terms in (21) are space-variant with the θ Δ .
( ) is the Doppler chirp rate term which will directly affect the azimuth focusing quality. Figure 9 shows the time-frequency lines of targets A T , B T and C T which are in a same range cell.
By introducing the frequency perturbation terms, the space-variation of ( )

Aizmuth Fousing
The phase terms in (21) are space-variant with the ∆θ. k 2,i (R i , ∆θ) is the Doppler chirp rate term which will directly affect the azimuth focusing quality. Figure 9 shows the timefrequency lines of targets T A , T B and T C which are in a same range cell. The slope of the time-frequency line mainly corresponding to the Doppler chirp rate which are γ A , γ B and γ C . In order to remove the space-variation of Doppler chirp rate, frequency perturbation terms are introduced. By an azimuth Fourier transform, the signal is transformed into Doppler domain: (24) where f dc = 2k 1,i (R i , ∆θ)/λ is the residual Doppler center after residual RWC of block i. The constant phase term in (24) is neglected, for it has no influence on the azimuth focusing. κ n,i (n = 2, 3, 4) can be expressed as By introducing the frequency perturbation terms, the space-variation of k 2,i (R i , ∆θ) can be removed in time domain. The perturbation function can be formed as (26) where α i and β i are coefficients to be determined. Multiplying (26) with (24), the signal is then transformed into azimuth time domain: where D n (n = 1, 2, 3, 4) is the coefficients after frequency perturbation which can be expressed as: The space-variation of D 3,i (R i , ∆θ) and D 4,i (R i , ∆θ) are neglected for the perturbation terms are small. For further analysis, D 2,i (R i , ∆θ) is expanded as Taylor series at ∆θ = 0. The first and second order coefficients can be derived as

Geometric Correction
The image of each angular block is focused in range-Doppler domain, as derived in (32). There is geometric distortion in azimuth. According to (28), the Doppler focused position is modulated by θ Δ , α and β . Assuming an arbitrary target in an angular block with coordinate ( ) The signal is divided in Doppler according to the angular division. Thus, there are some overlaps of different data blocks in Doppler. The final SAR image can be obtained by the direct combination of all angular block images.

Experimental Results
In this part, some experiments are made based on the parameters in Table 1 and the geometry in Figure 10. The key processing steps have been verified as follows.

Geometric Correction
The image of each angular block is focused in range-Doppler domain, as derived in (32). There is geometric distortion in azimuth. According to (28), the Doppler focused position is modulated by ∆θ, α and β. Assuming an arbitrary target in an angular block with coordinate (x, y). The corresponding yaw angle can be calculated as θ = atan(y/x). Comparing θ with θ i of each block, the block number of the target can be found. So, the reference range and angular deviation of target (x, y) in block i can be derived as The signal is divided in Doppler according to the angular division. Thus, there are some overlaps of different data blocks in Doppler. The final SAR image can be obtained by the direct combination of all angular block images.

Experimental Results
In this part, some experiments are made based on the parameters in Table 1 and the geometry in Figure 10. The key processing steps have been verified as follows.

Geometric Correction
The image of each angular block is focused in range-Doppler domain, as derived in (32). There is geometric distortion in azimuth. According to (28), the Doppler focused position is modulated by θ Δ , α and β . Assuming an arbitrary target in an angular block with coordinate ( ) The signal is divided in Doppler according to the angular division. Thus, there are some overlaps of different data blocks in Doppler. The final SAR image can be obtained by the direct combination of all angular block images.

Experimental Results
In this part, some experiments are made based on the parameters in Table 1 and the geometry in Figure 10. The key processing steps have been verified as follows.

Correction of Space-Variant RCM
In order to verify the method of space-variant RCMC, several targets with different reference range are selected. The target distribution geometry is shown in Figure 10a. Target A, B and C are located at near, center and far reference ranges, respectively. Then, two different RCMC methods are used. The simulation results are shown in Figure 11. By time domain RCMC without KT, there are still space-variant RCM components remained. It can be noticed in Figure 11a,c that the residual RCM is much larger than a range cell. By the proposed RCMC method, the space-variant RCM can be fully removed. The RCM can be accurately compensated, as shown in Figure 11d

Correction of Space-Variant RCM
In order to verify the method of space-variant RCMC, several targets with different reference range are selected. The target distribution geometry is shown in Figure 10a. Target A , B and C are located at near, center and far reference ranges, respectively. Then, two different RCMC methods are used. The simulation results are shown in Figure 11. By time domain RCMC without KT, there are still space-variant RCM components remained. It can be noticed in Figure 11a,c that the residual RCM is much larger than a range cell. By the proposed RCMC method, the space-variant RCM can be fully removed. The RCM can be accurately compensated, as shown in Figure 11d,f. Figure 11. RCM lines of targets in Figure 10

Performance of Geometric Correction and Azimuth Focusing
To verify the effectiveness of geometric correction, 9 × 9-point targets matrix with a distribution area of 5 km × 5 km are used. Nine targets are selected and the distribution geometry is shown in Figure 10b. By dividing the targets into three angular blocks, the imaging results are shown in Figure 12. Because the azimuth focusing position is modulated by the frequency perturbation terms, there is obvious geometric distortion in each angular block image. The geometric correction results of angular block images are shown in Figure 13a-c. According to the coordinates of selected targets-A1 (26,26), C1(2526, 26), A2 (26,426), C2(2526, 426), A3 (26,826) and C3(2526, 826)-the geometry distortion within each block has been corrected. By directly combining the angular block images, the final focused image can be obtained, as shown in Figure 13d. To check the image distortion in the final image, targets A1, C1, B2, C1 and C3 are selected in Figure 12. Their locations in Figure 13d are (26,26), (2526,26), (1276, 1276), (26,2526) and (2526,2526). The resolution cell of the interpolated image in Figure 13 are 2 m × 2 m. It can be found that the geometric distortion is controlled within one resolution cell.

Performance of Geometric Correction and Azimuth Focusing
To verify the effectiveness of geometric correction, 9 × 9-point targets matrix with a distribution area of 5 km × 5 km are used. Nine targets are selected and the distribution geometry is shown in Figure 10b. By dividing the targets into three angular blocks, the imaging results are shown in Figure 12. Because the azimuth focusing position is modulated by the frequency perturbation terms, there is obvious geometric distortion in each angular block image. The geometric correction results of angular block images are shown in Figure 13a A comparison is made to verify the effectiveness of the proposed method. The method in [10] is chosen as a reference. Considering the method is not designed to process TOPS SAR data. Processing steps of the proposed methods are used except the signal model. By the reference method, the signal model is low accurate after variable substitution. To evaluate the focusing performance, 2-D contour images of nine targets selected in Figure 10 are shown in Figures 14 and 15. Like the contour images in [28], the 2-D resolution is visible from the figure. The space-variant components of the third-order term are still remained. This will lead to azimuth defocusing in edge area, as shown in Figure 14. Compared with the reference method, the proposed range model can be directly matched with the signal after RWC which is more accurate than the reference one. The variation of the third-order phase term can be compensated as constant term within each angular block. So, the azimuth focusing performance can be improved, as shown in Figure 15. A comparison is made to verify the effectiveness of the proposed method. The method in [10] is chosen as a reference. Considering the method is not designed to process TOPS SAR data. Processing steps of the proposed methods are used except the signal model. By the reference method, the signal model is low accurate after variable substitution. To evaluate the focusing performance, 2-D contour images of nine targets selected in Figure 10 are shown in Figures 14 and 15. Like the contour images in [28], the 2-D resolution is visible from the figure. The space-variant components of the third-order term are still remained. This will lead to azimuth defocusing in edge area, as shown in Figure 14. Compared with the reference method, the proposed range model can be directly matched with the signal after RWC which is more accurate than the reference one. The variation of the third-order phase term can be compensated as constant term within each angular block. So, the azimuth focusing performance can be improved, as shown in Figure 15. A comparison is made to verify the effectiveness of the proposed method. The method in [10] is chosen as a reference. Considering the method is not designed to process TOPS SAR data. Processing steps of the proposed methods are used except the signal model. By the reference method, the signal model is low accurate after variable substitution. To evaluate the focusing performance, 2-D contour images of nine targets selected in Figure 10 are shown in Figures 14 and 15. Like the contour images in [28], the 2-D resolution is visible from the figure. The space-variant components of the third-order term are still remained. This will lead to azimuth defocusing in edge area, as shown in Figure 14. Compared with the reference method, the proposed range model can be directly matched with the signal after RWC which is more accurate than the reference one. The variation of the third-order phase term can be compensated as constant term within each angular block. So, the azimuth focusing performance can be improved, as shown in Figure 15. Furthermore, the 2-D PSLR, ISLR and IRW of these targets processed by the two methods are calculated, as listed in Tables 2 and 3 Furthermore, the 2-D PSLR, ISLR and IRW of these targets processed by the two methods are calculated, as listed in Tables 2 and 3.       Figure 10b by the reference method.  Table 3. PSLR, ISLR and IRW of selected targets in Figure 10b by the proposed method.    Figure 10b by the reference method. Table 3. PSLR, ISLR and IRW of selected targets in Figure 10b by the proposed method.

Discussion
This paper is based on previous works [14]. It provides another solution for highspeed maneuvering SAR imaging in TOPS mode. For the focusing performance, methods proposed in the two papers are in a same category. However, several main procedures are different: the preprocessing steps and the processing architecture. Signal ambiguity in Doppler is a feature of TOPS SAR data. The preprocessing steps in [14] is designed with reference to the methods in [16,17]. The key idea is to extend the data support region in Doppler by sacrificing that in azimuth time. It causes extra processing complexity when performing KT and other time domain compensations. For the method in this paper, the signal is unambiguous both in azimuth time and Doppler. So, the compensation functions can be easily applied without restrictions. A disadvantage is that the zeropadding increases the computational complexity. Thus, the parallel processing architecture is used to accelerate the data processing. It also has the potential ability to achieve higher focusing performance when the azimuth swath is larger.

Conclusions
This paper proposed an imaging algorithm for high-speed maneuvering platform SAR with squint TOPS mode. The range model is formed in local polar format which can overcome the mismatch between the range model and the signal after RWC. The Doppler folded signal is recovered by nonlinear derotation and zero-padding. The signal is then divided into several blocks according to the angular block division and can be processed in parallel. By KT, the space-variant RCM components are removed. Then, frequency perturbation terms are introduced to compensate the Doppler chirp rate variation. The focusing parameters are calculated according to the scene center of each angular block. Finally, the geometric distortion of each angular block is corrected by 2-D interpolation. The final focused image can be obtained by directly combining all the angular block images. Simulated SAR data has verified the effectiveness of the proposed algorithm.

Data Availability Statement:
The data is available to readers by contacting the corresponding author.