Research on Turning Motion Targets and Velocity Estimation in High Resolution Spaceborne SAR.

The development of high resolution SAR makes the influence of moving target more prominent, which results in defocusing and other unexplained phenomena. This paper focuses on the research of imaging signatures and velocity estimation of turning motion targets. In this paper, the turning motion is regarded as the straight line motion of continuous change of moving direction. Through the analysis of the straight line motion with constant velocity and the geometric modeling of the turning motion in spaceborne SAR, the imaging signatures of the turning motion target are obtained, such as the broken line phenomenon at the curve. Furthermore, a method for estimating the turning velocity is proposed here. The radial velocity is calculated by the azimuth offset of the turning motion target and the azimuth velocity is calculated by the phase error compensated in the refocusing process. The amplitude and direction of the velocity can be obtained by using both of them. The results of simulation and GF-3 data prove the accuracy of the analysis of turning motion imaging signatures, and they also show the accuracy and validity of the velocity estimation method in this paper.


Introduction
As an active remote sensing technology, SAR has the ability to work all day and all weather, being hardly affected by bad weather or light conditions. Spaceborne SAR has been widely concerned because of its fast response, wide coverage, and strong anti-strike ability [1]. With the development of high-resolution spaceborne SAR, the influence of moving target becomes more prominent, which becomes the focus of research [2,3]. Moving targets, such as vehicles, trains, ships, and some military targets are the key targets of SAR imaging, which has very important application value in military reconnaissance, ocean observation, and environmental detection. The development of high resolution makes the influence of moving target more prominent, resulting in defocusing and other unexplained phenomena. On the one hand, it affects the interpretation and detection of the target, on the other hand, it also brings the possibility of parameter estimation of the moving target. Generally speaking, it is very important to study the moving target of high-resolution SAR.
Presently, the research on high-resolution SAR moving target mainly focuses on four aspects. The first is the analysis of imaging signatures of moving target. Most of the literatures analyzed the imaging effect of the uniform moving target, and pointed out that the target movement would lead to the target azimuth position dislocation and target defocusing [4][5][6]. There are also literature analyses of the imaging effect of moving targets with acceleration [7], that is, the third phase error and the asymmetric distortion of the image. Another part of the literature analyzes the influence of the complex motion of the target, such as the rotation motion of the target itself [8], three-dimensional motion [9,10], vibration, and other The rest of this report is organized, as follows. Section 2 presents the geometry and signal models, along with the binary relationship between the velocity and residual slant range. The imaging signatures of moving target with constant velocity and turning motion target are presented here. Section 3 provides the algorithm for refocusing a moving target image and the method for estimating the velocity of the turning motion target. The results of imaging signatures of turning motion target, along with the image refocusing and velocity estimation, are presented in Section 4. Section 5 provides the conclusions.

Analysis of Moving Target with Constant Velocity in Spaceborne SAR
The instantaneous slant range varies with the azimuth time and then creates the synthetic aperture due to the relative motion between the sensor and a given target on the rotating Earth. In spaceborne SAR cases, vectors are commonly used in the calculation. Consider the geometric situation in Figure 1, which shows the satellite-Earth geometry for a given target (T) and satellite (S 1 ) at a certain azimuth time. The reference coordinate system used in the figure is the inertial coordinate system. S is the satellite position when azimuth time η=0. The slant range history of the static target T can be expressed as: where R s (η) and R t (η) denote the position vector of the satellite and the target, respectively, at a certain azimuth time. R, V st , and A st are, respectively, the relative range vector, velocity vector, and acceleration vector of the satellite and target at η = 0. Subsequently, we can obtain the value of the distance between the radar and the target, as: where R is the module of the vector R.
Sensors 2020, 20, x FOR PEER REVIEW 3 of 19 The rest of this report is organized, as follows. Section 2 presents the geometry and signal models, along with the binary relationship between the velocity and residual slant range. The imaging signatures of moving target with constant velocity and turning motion target are presented here. Section 3 provides the algorithm for refocusing a moving target image and the method for estimating the velocity of the turning motion target. The results of imaging signatures of turning motion target, along with the image refocusing and velocity estimation, are presented in Section 4. Section 5 provides the conclusions.

Analysis of Moving Target with Constant Velocity in Spaceborne SAR
The instantaneous slant range varies with the azimuth time and then creates the synthetic aperture due to the relative motion between the sensor and a given target on the rotating Earth. In spaceborne SAR cases, vectors are commonly used in the calculation. Consider the geometric situation in Figure 1, which shows the satellite-Earth geometry for a given target (T) and satellite (S1) at a certain azimuth time. The reference coordinate system used in the figure is the inertial coordinate system. S is the satellite position when azimuth time η=0. The slant range history of the static target T can be expressed as: where Rs(η) and Rt(η) denote the position vector of the satellite and the target, respectively, at a certain azimuth time. R, Vst, and Ast are, respectively, the relative range vector, velocity vector, and acceleration vector of the satellite and target at η=0. Subsequently, we can obtain the value of the distance between the radar and the target, as: where R is the module of the vector R. Assuming that T is moving along a straight line with constant velocity vector Vm, we can obtain the slant range of the moving target to the radar neglecting the high order items, as: Assuming that T is moving along a straight line with constant velocity vector V m , we can obtain the slant range of the moving target to the radar neglecting the high order items, as: Sensors 2020, 20, 2201 The moving target mainly affects the history of the slant range between the target and the radar. The slant range error that is introduced by target motion is the basis for SAR moving target imaging analysis. Through the previous analysis, we have obtained the slant range history of the stationary and moving targets, respectively, in Equation (2) and Equation (3). Subsequently, the slant range error of the moving target can be calculated, as: neglecting higher-order terms. Decomposing V m orthogonally into V x and V y , where V y directed along the radar boresight, and V x parallel to V st at azimuth time η = 0. We can obtain the following expressions: where γ is the look angle.
According to the relationships in Equation (5), we can rewrite Equation (4) as: In terms of the residual range expression, the phase error and Doppler frequency error can be written as: By analyzing the Doppler frequency error in Equation (8), we can get the imaging signatures of moving target SAR image, which are azimuth delocalization and azimuth defocusing phenomenon described in detail in the following paragraph.

Target Azimuth Delocalization
From the Doppler frequency error expression in Equation (8), we can obtain the Doppler centroid frequency error, as: According to the SAR imaging theory, the azimuth position of the target is related to the Doppler centroid frequency The Doppler centroid frequency error that was caused by the moving target will result in azimuth delocalization, which is related to Doppler modulation frequency and beam scanning speed, namely ground speed.
Let f r denotes the Doppler modulation frequency and V g denotes the ground speed. The azimuth position offset is shown, as: It can be seen from Equation (10) and Equation (9) that the factors affecting the azimuth position offset include the target radial velocity, the wavelength, the ground beam scanning speed, the Doppler Sensors 2020, 20, 2201 5 of 19 modulation frequency, and the look angle. This kind of azimuth delocalization of the target is usually called location error.
Many factors affect the location error of the target, among which the radial velocity of the target relative to the radar platform is an important aspect. It is the radial velocity of the target that causes the change of Doppler centroid frequency. The Doppler centroid frequency error is reflected as target azimuth delocalization in the SAR image. There is an approximately linear relationship between the target's azimuth position offset and radial velocity.

Target Azimuth Defocusing
From the Doppler frequency error expression in Equation (11), we can get the Doppler modulation frequency error, as: A matched filter mainly realizes the SAR azimuth imaging process. The design of matched filter is related to the Doppler information of signal, and the Doppler frequency modulation information is the main one. The target motion causes the change of the Doppler frequency modulation. The azimuth processing of the echo signal of the moving target with the Doppler frequency modulation rate of the stationary target will make it difficult to obtain the focused image due to the mismatch of the matching filter. Equation (11) shows that the most important factor affecting the Doppler frequency modulation is the azimuth velocity. Apparently, the azimuth velocity will generate a loss of focus in SAR image.
Through the above analysis, we can conclude that the signatures of moving target with constant velocity in SAR image are defocusing phenomenon mainly caused by azimuth velocity and azimuth delocalization mainly caused by radial velocity.

Analysis of Turning Motion Target
In this section, the turning motion model of multi-point targets, which is derived from the turning motion of train, will be constructed and its imaging signatures will be analyzed.
The geometric model of turning motion target is built in the ground plane in order to facilitate imaging analysis, which is shown in Figure 2. For the convenience of subsequent expression, the plane rectangular coordinate system is established in the ground plane, with the center point O of the scene as the origin of the coordinate system, the distance to the ground as the x-axis, and the azimuth position as the y-axis.
Generally, the curve is mainly composed of three parts. They are the straight part before entering the curve, the arc part of the curve, and the straight part leaving the curve. The two straight line parts are tangent to the arc part at the connection point. According to the geometry knowledge of a curve, four factors that determine the specific shape of a curve are: centering angle of arc part of curve, mark as β; 2.
radius of the arc part of the curve, mark as R l ; 3.
slope of line l1 at the beginning of the curve, mark as θ1; and, 4. Section AB is the straight track before entering the curve, OC is the straight track after leaving the curve, and BO is the curve between the two straight tracks, as shown in Figure 2. According to the curve geometry knowledge, the rotation center of the curve is the O' point marked in Figure 2 In fact, the part of the arc can be regarded as a small line segment with continuous change of slope. Subsequently, the turning movement is actually a "straight line movement" with continuous change of moving direction. Equation (6) shows that the slant range history error is mainly determined by the azimuth velocity component and the radial velocity component of the moving target. Taking a scattering point in any part of the curve as an example, where the target position is (rt,at) and the corresponding tangent slope is θt, then we can get: In the process of uniform turning movement, the component of velocity in azimuth and radial direction is changing due to the change of movement direction. The slant range history error relative to radar platform can be deduced from the range history error of linear moving target in Equation (6), as: where Vy(η) and Vx(η) represent the azimuth velocity component and radial velocity component, respectively, which are determined by the direction of motion: where Vm is the amplitude of target turning velocity and ϑ(η) is the target turning direction. In the curve part, the instantaneous motion direction of the target is the tangent direction of the target position. It is assumed that the target direction of motion at η=0 is θ0. According to the curve arc length and the amplitude of target velocity, the instantaneous direction of the target's motion can be calculated, as: Section AB is the straight track before entering the curve, OC is the straight track after leaving the curve, and BO is the curve between the two straight tracks, as shown in Figure 2. According to the curve geometry knowledge, the rotation center of the curve is the O' point marked in Figure 2 In fact, the part of the arc can be regarded as a small line segment with continuous change of slope. Subsequently, the turning movement is actually a "straight line movement" with continuous change of moving direction. Equation (6) shows that the slant range history error is mainly determined by the azimuth velocity component and the radial velocity component of the moving target. Taking a scattering point in any part of the curve as an example, where the target position is (r t ,a t ) and the corresponding tangent slope is θ t , then we can get: In the process of uniform turning movement, the component of velocity in azimuth and radial direction is changing due to the change of movement direction. The slant range history error relative to radar platform can be deduced from the range history error of linear moving target in Equation (6), as: where V y (η) and V x (η) represent the azimuth velocity component and radial velocity component, respectively, which are determined by the direction of motion: where V m is the amplitude of target turning velocity and ϑ(η) is the target turning direction. In the curve part, the instantaneous motion direction of the target is the tangent direction of the target position.
It is assumed that the target direction of motion at η=0 is θ 0 . According to the curve arc length and the amplitude of target velocity, the instantaneous direction of the target's motion can be calculated, as: Take Equation (15) into Equation (14), and rewrite Equation (13) as: The sine function and cosine function in Equation (16) are expanded by Taylor at η=0 to obtain: The polynomial representation of Equation (16) can be obtained by taking Equation (17) into the Equation (16): This calculation process is relatively tedious, therefore only the specific expression of the coefficient of the primary term is given here. The coefficients of other secondary terms are replaced by p n .
In terms of the residual range expression, the phase error, and Doppler frequency error can be written as: The errors in polynomial form are mainly the target azimuth position offset caused by primary term and the target azimuth defocusing phenomenon caused by higher term. Azimuth defocusing is a qualitative description, which can not be clearly matched in the results of multi-point target motion imaging. Therefore, this section focuses on the quantitative analysis of multi-point target azimuth location error.
According to Equation (10) and Equation (20), the azimuth position offset that is caused by turning movement of point target with position (r t ,a t ) can be obtained, as follows: It can be seen from Equation (21) that the azimuth position offset of the turning moving target is related to the slope of the tangent corresponding to the position at the time of azimuth center. θ 0 is a variable that changes with the target point, which can be expressed as: where (r',a') represents the stationary position of the target at any point of the turning movement. Therefore, the azimuth position offset of certain target can be obtained as: In fact, the product of V m and direction function cosθ(r',a') is the radial velocity of each target on the curve. The amplitude of the whole velocity V m is constant, so the different direction θ(r',a') of each target will lead to different radial velocity, which results in a different azimuth position offset. This different offset will make the targets no longer hold the shape of the curve, but the broken line of different shapes which depends on the magnitude of the velocity.
The specific imaging effects will be detailed in the following simulation results in Section 4.1. Figure 2 shows the turning motion model in this paper. In this paper, the turning motion is divided into three parts, two of which are uniform linear motion, and one is uniform turning motion. The speed of turning is equal to that of straight-line motion because the targets of turning motion are a whole (such as a train). Based on this assumption, we can determine the speed of turning by estimating the speed of a certain straight-line movement.

Moving Target Refocusing and Velocity Estimation of Turning Motion Target
The velocity of linear moving target can be divided into radial velocity and azimuth velocity. This paper deals with formed SLC SAR slices, that is, the region of interest (ROI). There are many methods of velocity estimation for this kind of moving target, among which the azimuthal velocity estimation is the most popular, while the radial velocity estimation is less studied and it has its application premise. In this paper, according to the signatures of the turning motion target, known is the original position of the target (for example, the rail in the image is the position where the train should be when it is stationary), the radial velocity of the target is calculated by the known target azimuth position offset. At the same time, the motion compensation algorithm that is based on autofocusing is used to refocus ROI target. The azimuth velocity of the target is calculated by the motion compensation phase and radial velocity of the target, combined with Equation (7). Figure 3 shows the overall algorithm flow chart. We number each computational block in order to distinguish the calculation steps of the method and expand the discussion of the method in detail. The italics in the Figure 3 represent the data dimension and domain output after the computational block, which is also the input of the next computational block. The process is briefly described, as follows: 1. intercept the target of the linear motion part as the ROI target, which size is M × N, M refers to azimuth direction while N refers to range direction; 2. obtain the equivalent raw data in Doppler domain sized M × N by applying inverse azimuth Fourier transform (IFFT) to ROI target; 3. autofocusing for the data obtained in the previous step through the minimum entropy phase adjustment algorithm and obtain the compensation phase curve in Doppler domain; 4. Obtain the refocused ROI target sized M × N in time domain by applying azimuth Fourier transform (FFT); 5. measure and calculate the displacement of azimuth position caused by target motion; 6. calculate the radial velocity V y according to Equation (9) and Equation; 7. The azimuth velocity V x is calculated by the compensation phase curve and the radial velocity estimation results combined with Equation (7); and, 8. calculate the amplitude V m and the direction θ m of the velocity through the radial velocity V y and the azimuthal velocity V x . 7. The azimuth velocity Vx is calculated by the compensation phase curve and the radial velocity estimation results combined with Equation (7); and, 8. calculate the amplitude Vm and the direction θm of the velocity through the radial velocity Vy and the azimuthal velocity Vx.  The motion compensation algorithm and velocity estimation algorithm based on autofocusing are described in detail below.

Motion Compensation Algorithm
This section, corresponding to steps 1-4 in Figure 3, is aimed to obtain the refocused moving target by autofocusing method, and to obtain the compensation phase curve that is caused by motion error, which is used to estimate the velocity of the target. The following are steps in detail.
Step 1. Obtain the ROI target We have a SAR complex image of a turning motion target. We intercept the linear motion part from the formed image as ROI target. Subsequently, the following technique used in the paper deals with the formed SAR ROI target. The size of the ROI target is M × N, M refers to azimuth direction, while N refers to range direction, and the initial state is in time domain.
Step 2. Azimuth IFFT The ROI target is inverted in order to obtain the equivalent raw data, which are the inputs of the motion compensation technique that were acquired by applying inverse Fourier transform (IFFT) and inverse compression in the azimuth direction. Thus, we can obtain the data in the Doppler frequency domain before the azimuth compression, as: where fη is the Doppler frequency, t is the range time, σ is the backscattering coefficient, λ is the radar wavelength, c is the light speed, and B is the bandwidth of the transmitted signal. The motion compensation algorithm and velocity estimation algorithm based on autofocusing are described in detail below.

Motion Compensation Algorithm
This section, corresponding to steps 1-4 in Figure 3, is aimed to obtain the refocused moving target by autofocusing method, and to obtain the compensation phase curve that is caused by motion error, which is used to estimate the velocity of the target. The following are steps in detail.
Step 1. Obtain the ROI target We have a SAR complex image of a turning motion target. We intercept the linear motion part from the formed image as ROI target. Subsequently, the following technique used in the paper deals with the formed SAR ROI target. The size of the ROI target is M × N, M refers to azimuth direction, while N refers to range direction, and the initial state is in time domain.
Step 2. Azimuth IFFT The ROI target is inverted in order to obtain the equivalent raw data, which are the inputs of the motion compensation technique that were acquired by applying inverse Fourier transform (IFFT) and inverse compression in the azimuth direction. Thus, we can obtain the data in the Doppler frequency domain before the azimuth compression, as: where f η is the Doppler frequency, t is the range time, σ is the backscattering coefficient, λ is the radar wavelength, c is the light speed, and B is the bandwidth of the transmitted signal.
The target is assumed to be a point-like scatterer, and the motion errors that are caused by the platform are assumed to have been perfectly compensated by chirp-scaling algorithm. Thus, we obtain the range migration, as: where R is the slant range at synthetic aperture time η=0, as mentioned in Section 2. ∆R(f η ) is the residual slant range expression that is caused only by the moving target in the Doppler domain.
Step 3. The minimum entropy phase adjustment autofocusing algorithm The range migration causes a variation in the phase terms. For the compensation of the time-varying Doppler phase, the minimum entropy phase adjustment based on the gradient descent algorithm is chosen [38]. This is found by iteratively solving an equation, which is derived by minimizing the entropy of the image and ultimately provides the compensating phase curve. It is worth mentioning that the compensated phase error curve is obtained in the azimuth Doppler domain, whose dimension is M × 1, that is, it changes along the azimuth direction.
Step 4. Azimuth FFT The data after motion compensation can be written, as: Finally, for the image formation, an FFT along the azimuth direction can be used to refocus the image.

Velocity Estimation Method
This section deals with the estimation velocity method that corresponds to steps 5-8 in Figure 3. The velocity estimation method in this section is based on the error analysis presented in Section 2 and the motion compensation technique in the previous section. The followings are steps in detail.
Step 5. Measure the displacement of azimuth position We have known the correct position of the target when it is stationary (for example, the rail is the position where the train should be when it is stationary). By comparing the position of the moving target with the position it should be when it is stationary along the azimuth direction, we can obtain the number of azimuth offset points as ∆B. It is assumed that the azimuth spatial resolution of the image is ρ a , and we can get the azimuth offset expressed in meters, as: Step 6. Calculate the radial velocity Vy The radial velocity of the target is calculated according to the offset of the turning moving target from the original position, which is the target azimuth position offset ∆a. We can obtain the V y according to Equation (9) and Equation (10): Step 7. Calculate the azimuth velocity Vx The basic principle of the velocity estimation method in this step is that, when the defocused ROI target is processed by autofocusing, the imaging quality is obviously improved and the focusing effect is good, the phase error caused by motion are completely compensated, that is: where ∆ϕ(f η ) represents the phase error of azimuth frequency domain that is caused by motion. ∆ϕ(m) represents the compensated phase obtained by the phase correction algorithm.
∆ϕ(f η ) contains parameters related to the velocity of motion. As long as the specific analytical expression of ∆ϕ(f η ) is obtained, the corresponding velocity estimates can be obtained.
The analytical expression for ∆ϕ(f η ) is derived below. In order to facilitate the analysis, the two-dimensional spectrum of the state target echo is rewritten, as follows: where G rcm (f η ,R) is the range migration term in azimuth frequency domain. G ac (f η ,R) is the azimuth modulation term. G o (f η ,R) is the second distance compression term. The expression of the three is: where f τ is the range frequency, V r is the equivalent velocity used in the imaging algorithm, and f 0 is the carrier frequency. The Doppler frequency error of moving target obtained from Section 2 is: In the basic theory of SAR imaging, the Doppler frequency information of the target represents the relationship between the azimuth time and the azimuth frequency of the target point. It can be considered that the azimuth frequency of moving target has changed as compared with that of stationary target, thus changing the expression of target echo in two-dimensional frequency domain. If the azimuth frequency of a moving target is f η ' and that of a stationary target is f η , the relationship between them can be approximately considered as: The azimuth modulation term of moving target in azimuth frequency domain are: After the matched filtering of the imaging algorithm, the residual phase error is: where G * ac (f η ,R) is the matched filter in the frequency domain of the azimuth of the stationary target. The results of Equation (35) are expanded by Taylor series on f η = 0 and retained to the quadratic term, as: The linear part of the phase error compensation term is not accurate due to the motion compensation algorithm, so only the quadratic term of the phase error compensation term can be used for velocity estimation.
Using f η as the independent variable to do quadratic polynomial fitting to ∆ϕ(m) to obtain the coefficient of quadratic term p 2 : The parameters in the Equation (37) are all known, and the result of V x can be obtained by solving the quadratic equation of one variable.
Step 8. Calculate the amplitude and direction of the estimated velocity The amplitude and direction of velocity are obtained, as follows, using V y and V x : where V m and θ m denote the amplitude and direction of the velocity, respectively.

Simulation Results and Actual Data of Imaging Signatures of Turning Motion
This section mainly shows the simulation results of turning movement at different turning speeds to prove the accuracy and rationality of the analysis in Section 2. Figure 2 shows the shape of simulation curve. Table 1 shows the simulation curve parameters. Table 2 shows the simulation radar parameters. In this simulation experiment, a total of twenty-three point targets are simulated, including seven straight lines entering the curve, nine curves, and seven straight lines out of the curve. Figure 4a shows the results of static curve target. It can be seen that the turning model is relatively accurate. In this experiment, three turning velocities are simulated, which are 20 km/h, 45 km/h, and 65 km/h. This kind of different azimuth offset makes the moving target of the curve show a line shape, rather than the original curved shape. When compared with the results of different velocities, it can be seen that the larger the amplitude of velocity is, the more obvious the deflection trace of the curve target is, and the more inconsistent the shape of the original curve is. When the amplitude of velocity reaches a certain value (45 km/h in the simulation), the moving target of the curve can even appear the phenomenon of approximately right angle bend. The broken line phenomenon proves that the analysis in Section 2 is correct.   In the 3m resolution image of GF-3, we also found that the turning movement, which caused the broken line phenomenon, as shown in Figure 5. It can be seen from the optical image that the red part in the Figure 5b is a railway track, so the brighter line in the Figure 5a is supposed to be a turning train. In the 3m resolution image of GF-3, we also found that the turning movement, which caused the broken line phenomenon, as shown in Figure 5. It can be seen from the optical image that the red part in the Figure 5b is a railway track, so the brighter line in the Figure 5a is supposed to be a turning train.  It can be seen from the figure that the train is separated from the rail track. The track is a normal curve model with two straight sections connecting with the curve section smoothly, while the train shows three broken lines. The straight parts of both ends of the train are consistent with the straight parts of the rail, and there is an identical azimuth position offset, which is the signatures of the uniform linear moving targets. However, the curve section of the train is completely different from the original rail track shape, which is caused by the continuous change of the radial velocity during the turning, resulting in the continuous change of the azimuth position offset. Different position targets have a different azimuth position offset, which will make the targets no longer hold the shape of the curve, but the broken line of different shapes, which depends on the magnitude of the velocity. Therefore, the imaging train's shape is different from the original shape, showing a broken line shape rather than an arc shape.
Through the analysis of Sections 2.1 and 2.2, it can be known that the radial velocity determines the target azimuth position offset, while the radial velocity of the turning moving target is related to the amplitude of the target velocity and the moving direction. The reason why the train in Figure 5 can have a broken line phenomenon of approximately right angle is that the continuous change of azimuth position offset that is caused by its moving velocity and direction is just matched, as good as the result of simulation v = 50km/h in Figure 4c. If the speed increases or decreases, it will also show the same phenomenon in the simulation, as in Figures 4b and 4d.
Generally speaking, the phenomenon of train broken line in the Figure 5 is caused by different azimuth position offset when turning. This is consistent with the simulation results. It shows that the analysis of imaging signatures of turning motion in this paper is correct.

Simulation Results and Actual Data Results of Velocity Estimation
This section mainly shows the velocity estimation results of turning movement to prove the accuracy and effectiveness of the method presented in Section 3. The experimental data are the moving target data of three turning velocities generated by the simulation in the previous section. In this experiment, the straight-line part of the entrance curve is used for velocity estimation.
We intercept the straight part of the turning motion for autofocusing, which has been circled in red in the Figures 4b-d. For convenience, only the result of V = 45km/ h is shown here as an example. Figures 6a-b shows the comparison before and after the focus of the moving target. It can be seen that the image quality of the point target after the autofocus has improved significantly. We extract this target for quality evaluation. The quality evaluation results are shown in Table 3 and the range and azimuth waveforms are shown in Figures 6c-d. It can be seen that the azimuth waveform and the range waveform meet the index, indicating that the phase error that is caused by the motion has been compensated. The compensation phase is shown in Figure 6e. At the same time, we measure the azimuth position offset between the moving point target and the original stationary target. Using the Equation (28) and Equation (37), we can calculate the velocity estimation result. Table 4 shows It can be seen from the figure that the train is separated from the rail track. The track is a normal curve model with two straight sections connecting with the curve section smoothly, while the train shows three broken lines. The straight parts of both ends of the train are consistent with the straight parts of the rail, and there is an identical azimuth position offset, which is the signatures of the uniform linear moving targets. However, the curve section of the train is completely different from the original rail track shape, which is caused by the continuous change of the radial velocity during the turning, resulting in the continuous change of the azimuth position offset. Different position targets have a different azimuth position offset, which will make the targets no longer hold the shape of the curve, but the broken line of different shapes, which depends on the magnitude of the velocity. Therefore, the imaging train's shape is different from the original shape, showing a broken line shape rather than an arc shape.
Through the analysis of Sections 2.1 and 2.2, it can be known that the radial velocity determines the target azimuth position offset, while the radial velocity of the turning moving target is related to the amplitude of the target velocity and the moving direction. The reason why the train in Figure 5 can have a broken line phenomenon of approximately right angle is that the continuous change of azimuth position offset that is caused by its moving velocity and direction is just matched, as good as the result of simulation v = 50km/h in Figure 4c. If the speed increases or decreases, it will also show the same phenomenon in the simulation, as in Figure 4b,d.
Generally speaking, the phenomenon of train broken line in the Figure 5 is caused by different azimuth position offset when turning. This is consistent with the simulation results. It shows that the analysis of imaging signatures of turning motion in this paper is correct.

Simulation Results and Actual Data Results of Velocity Estimation
This section mainly shows the velocity estimation results of turning movement to prove the accuracy and effectiveness of the method presented in Section 3. The experimental data are the moving target data of three turning velocities generated by the simulation in the previous section. In this experiment, the straight-line part of the entrance curve is used for velocity estimation.
We intercept the straight part of the turning motion for autofocusing, which has been circled in red in the Figure 4b-d. For convenience, only the result of V = 45km/ h is shown here as an example. Figure 6a-b shows the comparison before and after the focus of the moving target. It can be seen that the image quality of the point target after the autofocus has improved significantly. We extract this target for quality evaluation. The quality evaluation results are shown in Table 3 and the range and azimuth waveforms are shown in Figure 6c-d. It can be seen that the azimuth waveform and the range waveform meet the index, indicating that the phase error that is caused by the motion has been compensated. The compensation phase is shown in Figure 6e. At the same time, we measure the azimuth position offset between the moving point target and the original stationary target. Using the Equation (28) and Equation (37), we can calculate the velocity estimation result. Table 4 shows the velocity estimation results of simulation. The estimation results are basically consistent with the actual simulation velocity, which proves the correctness of the proposed method.
The following is a brief analysis and assessment of the simulation results. Firstly, the target motion compensation algorithm that is based on autofocusing is exact and effective. The point target result is evaluated to be in accord with the index, as shown in Table 3 and Figure 6, which illustrates that the algorithm performs well. Next, the velocity estimation method is effective. Although there are small errors between the velocity estimation results and the real target velocity and direction, they are generally consistent with each other, which illustrates that the velocity estimation method performs well. Finally, we discuss the reason of the velocity estimation error. It can be seen from the results that the error of velocity estimation is 0.15-0.25 m/s (we converted the unit km/h in Table 4 into m/s) and 0.5-3 degrees. In fact, the error mainly comes from two ways, which are the error caused by the azimuth and radial velocity estimation method. On the one hand, the azimuth position offset is obtained by measuring the offset points and azimuth resolution, and its accuracy is limited by the resolution. The error offset of 3m calculated from simulation parameters is enough to produce the error of 0.05-0.1m/s; on the other hand, the azimuth velocity is calculated by the quadratic term of the phase error compensation curve fitting. There are small errors in fitting process and approximate omissions in equation calculation, which leads to the estimation error of azimuth velocity. We will continue to explore in-depth in the future research.
The actual data still select the train target in Figure 5a. Table 5 shows the parameters of the data. The straight-line part before entering the curve is autofocused, and the results are shown in Figure 7. It can be seen that the quality of train target has been improved after autofocusing as the enlarged part shown by the red line in Figure 7a-b. Figure 7c shows the phase error curve. We measured the displacement of the train and the rail track along the azimuth direction that is shown in Figure 7d. Using the Equation (28) and Equation (37), we can obtain the velocity estimation results that are shown in Table 6. the velocity estimation results of simulation. The estimation results are basically consistent with the actual simulation velocity, which proves the correctness of the proposed method.
The following is a brief analysis and assessment of the simulation results. Firstly, the target motion compensation algorithm that is based on autofocusing is exact and effective. The point target result is evaluated to be in accord with the index, as shown in Table 3 and Figure 6, which illustrates that the algorithm performs well. Next, the velocity estimation method is effective. Although there are small errors between the velocity estimation results and the real target velocity and direction, they are generally consistent with each other, which illustrates that the velocity estimation method performs well. Finally, we discuss the reason of the velocity estimation error. It can be seen from the results that the error of velocity estimation is 0.15-0.25 m/s (we converted the unit km/h in Table 4 into m/s) and 0.5-3 degrees. In fact, the error mainly comes from two ways, which are the error caused by the azimuth and radial velocity estimation method. On the one hand, the azimuth position offset is obtained by measuring the offset points and azimuth resolution, and its accuracy is limited by the resolution. The error offset of 3m calculated from simulation parameters is enough to produce the error of 0.05-0.1m/s; on the other hand, the azimuth velocity is calculated by the quadratic term of the phase error compensation curve fitting. There are small errors in fitting process and approximate omissions in equation calculation, which leads to the estimation error of azimuth velocity. We will continue to explore in-depth in the future research.
The actual data still select the train target in Figure 5a. Table 5 shows the parameters of the data. The straight-line part before entering the curve is autofocused, and the results are shown in Figure 7. It can be seen that the quality of train target has been improved after autofocusing as the enlarged part shown by the red line in Figures 7a-b. Figure 7c shows the phase error curve. We measured the displacement of the train and the rail track along the azimuth direction that is shown in Figure 7d. Using the Equation (28) and Equation (37), we can obtain the velocity estimation results that are shown in Table 6.           The rail track is located near Guanting Reservoir in Beijing with longitude E115.61 and latitude N40.22, which can be found in the optical map. The direction of the track is approximately 40-50 degrees west north, which is not much different from our estimation result. The estimated amplitude of velocity of the train is 65.7 km/h, which is a reasonable value for the freight train. It can show the correctness and validity of the method.

Conclusions
This paper mainly studies the imaging signatures of SAR image of turning motion, and it points out that the turning motion will produce the phenomenon of broken line. In this paper, the turning velocity is estimated by using the displacement of the turning movement, which deviates from the original position and the phase error obtained by autofocusing. The results of the simulation and actual data prove the accuracy of the analysis of turning motion imaging signatures, and also show the accuracy and validity of the velocity estimation method in this paper.