A Modified Azimuth Weighting Method in a Two-Step Process Approach for Sliding Spotlight Data Processing

Low sidelobes are important and essential in all SAR (Synthetic Aperture Radar) images, regardless of the imaging mode, for fewer artificial targets. For strip-map mode all targets overlap in frequency, which is convenient to suppress sidelobes. However, weighting requires total overlap in the time or frequency domain, which a sliding spotlight signal could not satisfy. Furthermore, the wavelength cannot be regarded as a constant value under the condition of a wideband chirp signal, which leads to the variation of the Doppler bandwidth along with the range frequency. In this article, an azimuth weighting method is proposed that considers the influence of a wideband based on a two-step algorithm. The computer simulation is given to verify the presented method.


Introduction
By means of properly steering a radar antenna in the along-track direction, SAR enlarges its applications in both civil and military domains, especially in the sliding spotlight mode [1][2][3][4][5][6]. Compared to the traditional spotlight mode, the sliding spotlight mode could enable wider coverage in the azimuth direction. Just as in the spotlight case, the key problem with the sliding spotlight focus processing is in the aliasing of the azimuth echo.
However, the shift of the Doppler centroid introduced by beam rotation will cause the Doppler bandwidth of the scenario to be larger than the pulse repetition frequency (PRF) [7][8][9]. Additionally, under the condition of a wideband signal, the variation of residual Doppler bandwidth along the range frequency must be taken into consideration while conducting the weighting operation.
The classical two-step imaging algorithm is introduced to overcome the azimuth spectrum aliasing [7,[10][11][12], but the azimuth residual Doppler bandwidth is not carefully analyzed and the sidelobe suppression is not considered, especially for high-resolution imaging. Sub-aperture algorithms are another way to deal with the azimuth spectrum aliasing [6,13], but it is difficult to conduct an azimuth weighting operation since the different azimuth targets share different time support domains and different frequency support domains. Based on the azimuth frequency de-ramping principle, a novel processing approach is proposed [14,15], with the disadvantage of more FFT (Fast Fourier Transform) and IFFT (Inverse Fast Fourier Transform) operations and lower efficiency. Furthermore, a new algorithm named band azimuth scaling (BAS) is proposed for both TOPS (Terrain Observation by Progressive Scans) and sliding spotlight data imaging [16][17][18][19][20]. However, the BAS algorithm needs a sub-aperture combination. Based on a frequency scaling algorithm (FSA) and spectral analysis where  is the instantaneous squint angle and T is the whole acquisition interval, and For azimuth beam scanning at a constant rotation rate leading to r T  

, the steering factor
A is defined as [14]: Assuming that the azimuth beam steering from fore to aft is positive, the azimuth resolution of the sliding-spotlight mode can be approximately computed as [13]: For azimuth beam scanning at a constant rotation rate leading to ω r = θ/T, the steering factor A is defined as [14]: Sensors 2017, 17, 220 3 of 17 Assuming that the azimuth beam steering from fore to aft is positive, the azimuth resolution of the sliding-spotlight mode can be approximately computed as [13]: where V s and V g are the physical velocity of the SAR sensor and the footprint velocity, respectively without taking into account the azimuth beam steering. D a is the antenna length and γ w,a is the impulse response width broadening factor due to an azimuth processing window. The approximation V g /V s γ w,a ≈ 1 is reasonable in the majority case.

Properties of the Echo Signal
If the transmitter illuminates the target scene with a baseband chirp signal p(t) and a point target locates at (x 0 , R), the echo of transmitted pulse can be expressed as: where λ = c/ f 0 is the wavelength, K r = B/T p is the FM (frequency modulation) rate. V f is the footprint velocity taking into account azimuth beam steering. τ and t a are fast time and slow time variables, respectively. The azimuth signal of the echo is aliasing which cannot be weighted now. Based on rectilinear imaging geometry of sliding-spotlight mode, R(t a ) can be written as: where V r ≈ V g · V s . V r is usually used for the imaging focus and V s is used for calculating the Doppler bandwidth.

Precious Azimuth Weighting Method in a Two-Step Algorithm
Generally speaking, the sidelobes after pulse compression are not as low as expected, which leads to shade weak objects. Thus, a weighting technique is necessary in the image processing. The comparison between employing a weighting method and without weighting processing in the sliding-spotlight mode is shown in Figure 2.
Sensors 2017, 17, 220 3 of 17 where s V and g V are the physical velocity of the SAR sensor and the footprint velocity, respectively without taking into account the azimuth beam steering. a D is the antenna length and , wa  is the impulse response width broadening factor due to an azimuth processing window. The approximation VV   is reasonable in the majority case.

Properties of the Echo Signal
If the transmitter illuminates the target scene with a baseband chirp signal () pt and a point target locates at   0 , xR, the echo of transmitted pulse can be expressed as: can be written as: where r V is usually used for the imaging focus and s V is used for calculating the Doppler bandwidth.

Precious Azimuth Weighting Method in a Two-Step Algorithm
Generally speaking, the sidelobes after pulse compression are not as low as expected, which leads to shade weak objects. Thus, a weighting technique is necessary in the image processing. The comparison between employing a weighting method and without weighting processing in the sliding-spotlight mode is shown in Figure 2.  Real data from GF-3 is shown as Figure 2. The quantization methods are same. It clearly shows the weighting effects in the SAR image. The sidelobes of the strong target in the right picture are stronger than the surrounding targets and the weak targets are buried. Compared to the right figure, which is not using a weighting operation, the sidelobes of the left picture are lower.
In this section, based on a deramping method, the total weighting method is proposed in the first part, the comparison of the azimuth weighting processing in the range time domain [26] and proposed method are analyzed in the second part.

Azimuth Preprocessing and Weighting Method
Considering the rotation of the azimuth beam under the curved Earth geometry, the instantaneous range from the satellite to the virtual rotation center R c (t a ) can be expressed as: Next, the f dc in the virtual point can be calculated via f dc_rot (t a ) = −2/λ · dR c (t a )/dt a . Then the azimuth variation rate f rot causing by azimuth beam rotation can be expressed as: Under the typical low orbit parameter and carrier frequency is C band, f rot can be expressed approximately as: where θ r is equivalent instantaneous squint angle. Thus, the azimuth bandwidth caused by the beam rotation can be expressed as: Compared with the strip-map mode, the azimuth original bandwidth caused by the beam width is: Finally, the azimuth total Doppler bandwidth can be expressed as: In general, the total bandwidth is larger than the PRF in the sliding spotlight mode. In order to overcome the aliasing of the azimuth echo in the frequency domain, azimuth convolution processing is conducted, which is the key point of azimuth preprocessing. The quadratic phase signal is expressed as: While conducting the azimuth weighting processing in range time domain, the azimuth echo preprocessing can be accomplished by employing: Substituting Equation (12) into Equation (13), Equation (13) can be rewritten as: From Equation (14), the azimuth convolution processing includes three parts: dechirp processing, Fourier transform, and phase compensation processing. Figure 4 shows the azimuth timefrequency relationship after dechirp processing. Where f dr = 2V 2 r /(λ · R 0 ) is the azimuth FM rate. T s represents synthetic aperture time and B o is the Doppler course of a single point, respectively.
In general, the total bandwidth is larger than the PRF in the sliding spotlight mode. In order to overcome the aliasing of the azimuth echo in the frequency domain, azimuth convolution processing is conducted, which is the key point of azimuth preprocessing. The quadratic phase signal is expressed as: While conducting the azimuth weighting processing in range time domain, the azimuth echo preprocessing can be accomplished by employing: Substituting Equation (12) into Equation (13), Equation (13) can be rewritten as: From Equation (14), the azimuth convolution processing includes three parts: dechirp processing, Fourier transform, and phase compensation processing. Figure 4 shows the azimuth time-frequency relationship after dechirp processing. Note that the wavelength is defined for a single-frequency signal but, for the chirp signal, the wavelength is not a constant value, so it can be expressed approximately as 0 cf   , and the azimuth convolution can be carried out both in the range time domain and range frequency domain under the condition of the narrowband signal, while under the condition of the wideband signal, the convolution must be carried out in range frequency domain [7] and: The parameter   equal f   is defined in the range frequency domain, which is related with f  .
This parameter indicates all wavelengths in the chirp signal and can be used in the range frequency domain to describe the wavelength for the chirp signal precisely instead of in the range time domain. The error analysis of the azimuth weighting method in [26] is analyzed in the next part.
Conducting the azimuth weighting processing in the range frequency domain, the azimuth echo preprocessing by convolution can be written as: Combining Equations (14)-(16), then: is derived. Implementing Equation (17) in the discrete domain, it can be written as: Note that the wavelength is defined for a single-frequency signal but, for the chirp signal, the wavelength is not a constant value, so it can be expressed approximately as λ ≈ c/ f 0 , and the azimuth convolution can be carried out both in the range time domain and range frequency domain under the condition of the narrowband signal, while under the condition of the wideband signal, the convolution must be carried out in range frequency domain [7] and: The parameter λ equal ( f τ ) is defined in the range frequency domain, which is related with f τ . This parameter indicates all wavelengths in the chirp signal and can be used in the range frequency domain to describe the wavelength for the chirp signal precisely instead of in the range time domain. The error analysis of the azimuth weighting method in [26] is analyzed in the next part.
Conducting the azimuth weighting processing in the range frequency domain, the azimuth echo preprocessing by convolution can be written as: Combining Equations (14)- (16), then: Sensors 2017, 17, 220 7 of 17 is derived. Implementing Equation (17) in the discrete domain, it can be written as: where: then Equation (18) can be expressed as: which means using FFT is a good choice to improve algorithm efficiency. After dechirp processing the signal totally overlaps, so the weighting operation can be carried out at this time. The total residual bandwidth can be written as: The azimuth weighting method in [26] is carried out in the range time domain, which is accurate in the high-resolution imaging mode, but this method causes significant errors while the bandwidth is about an order of the carrier frequency. It is impossible to consider the effect by the range frequency if the azimuth preprocessing is conducted in the range time domain. Only by updating the f r in the range frequency domain is the azimuth window is precisely produced. A detailed analysis and more experiments illustrate this below.
We further remark that the parameter design of PRF must satisfy The complete algorithm is shown in Figure 5. Firstly, a dechirp operation is carried out for removing the azimuth beam rotation; here, a weighting operation is applied. Secondly, the imaging algorithm, i.e., the chirp scaling algorithm [27] (pp. 206-210) [28] or nonlinear chirp scaling algorithm [29][30][31] is applied to focus the SAR image. The main difference between the presented weighting method in this article and the weighting method in [26] is in which domain the weighting operation is performed, which influences the sidelobe suppression performance, especially under the condition of high resolution. The complete algorithm is shown in Figure 5. Firstly, a dechirp operation is carried out for removing the azimuth beam rotation; here, a weighting operation is applied. Secondly, the imaging algorithm, i.e., the chirp scaling algorithm [27] (pp. 206-210) [28] or nonlinear chirp scaling algorithm [29][30][31] is applied to focus the SAR image. The main difference between the presented weighting method in this article and the weighting method in [26] is in which domain the weighting operation is performed, which influences the sidelobe suppression performance, especially under the condition of high resolution.

The Comparison of the Former Weighting Method and Proposed Weighting Method
Up to now, a weighting method in azimuth preprocessing has been proposed, since the error of the proposed method is small, but the error of the former weighting method is great. This article is supposed to give an analysis of the error of the weighting and conduct a deramping operation in the range time domain [26], and perform weighting and conduct a deramping operation in the range frequency domain.
Actually, the deterioration of the azimuth resolution and PSLR (Peak Side Lobe Ratio) indeed changes along the varying resolution; the reason is complex and a detailed discussion is illustrated in this part.
Consider two different processing methods: Due to the range cell migration (RCM), the amplitude of a single target in the range-Doppler

The Comparison of the Former Weighting Method and Proposed Weighting Method
Up to now, a weighting method in azimuth preprocessing has been proposed, since the error of the proposed method is small, but the error of the former weighting method is great. This article is supposed to give an analysis of the error of the weighting and conduct a deramping operation in the range time domain [26], and perform weighting and conduct a deramping operation in the range frequency domain.
Actually, the deterioration of the azimuth resolution and PSLR (Peak Side Lobe Ratio) indeed changes along the varying resolution; the reason is complex and a detailed discussion is illustrated in this part.
Consider two different processing methods:

Weighting and Conducting a Deramping Operation in the Range Time Domain
• The loss of resolution.
Due to the range cell migration (RCM), the amplitude of a single target in the range-Doppler domain is shown in Figure 6.
In Figure 6, the range-Doppler spectrum is lost due to weighting operation. However, notice that the updating of f r along the range direction is impossible, which is due to the aliasing of the range signal in the time domain. Due to the RCM, after the weighting operation [26] the residual range-Doppler spectrum is curved, as the red line shows. Finally, resolution will be lost both along the range direction and azimuth direction because of the two-dimensional spectrum losses. Sensors 2017, 17, 220 9 of 17

Azimuth(Hz )
Range(km) Figure 6. Simulated image using the former weighting method [26] with the deramping operation in the range time domain. The red line shows the residual spectrum in the range-Doppler domain after a deremping operation while using the former weighting method [26], and the dotted line shows the spectrum in the range-Doppler domain after a deremping operation without a weighting operation.
In Figure 6, the range-Doppler spectrum is lost due to weighting operation. However, notice that the updating of r f along the range direction is impossible, which is due to the aliasing of the range signal in the time domain. Due to the RCM, after the weighting operation [26] the residual range-Doppler spectrum is curved, as the red line shows. Finally, resolution will be lost both along the range direction and azimuth direction because of the two-dimensional spectrum losses.


The loss of PSLR.
The azimuth PSLR and ISLR (Integrated Side Lobe Ratio) will be lost due to the wrong boundary of the weighting operation. The effect of different errors in the weighting operation is shown in Figure 7; if the length of the window is shorter than the effective bandwidth, resolution loss occurs. If the length of the window is larger than the effective bandwidth, PSLR loss occurs due to the effective window being just a part of the total window.
To sum up, the loss will be found both in the range direction and the azimuth direction, which includes resolution loss and PSLR loss. If the resolution is higher, the resolution loss and PSLR loss is higher, which is unacceptable. The quantitative analysis is too complicated and it is related with r B , 0 f , and azimuth resolution. Additionally, the necessity of the quantitative error analysis is not enough due to the former weighting method being easy to improve by weighting in the range frequency domain. Thus, the simulation results are illuminated in Section 4. Figure 6. Simulated image using the former weighting method [26] with the deramping operation in the range time domain. The red line shows the residual spectrum in the range-Doppler domain after a deremping operation while using the former weighting method [26], and the dotted line shows the spectrum in the range-Doppler domain after a deremping operation without a weighting operation.
• The loss of PSLR.
The azimuth PSLR and ISLR (Integrated Side Lobe Ratio) will be lost due to the wrong boundary of the weighting operation.
The effect of different errors in the weighting operation is shown in Figure 7; if the length of the window is shorter than the effective bandwidth, resolution loss occurs. If the length of the window is larger than the effective bandwidth, PSLR loss occurs due to the effective window being just a part of the total window. Range(km) Figure 6. Simulated image using the former weighting method [26] with the deramping operation in the range time domain. The red line shows the residual spectrum in the range-Doppler domain after a deremping operation while using the former weighting method [26], and the dotted line shows the spectrum in the range-Doppler domain after a deremping operation without a weighting operation.
In Figure 6, the range-Doppler spectrum is lost due to weighting operation. However, notice that the updating of r f along the range direction is impossible, which is due to the aliasing of the range signal in the time domain. Due to the RCM, after the weighting operation [26] the residual range-Doppler spectrum is curved, as the red line shows. Finally, resolution will be lost both along the range direction and azimuth direction because of the two-dimensional spectrum losses.


The loss of PSLR.
The azimuth PSLR and ISLR (Integrated Side Lobe Ratio) will be lost due to the wrong boundary of the weighting operation. The effect of different errors in the weighting operation is shown in Figure 7; if the length of the window is shorter than the effective bandwidth, resolution loss occurs. If the length of the window is larger than the effective bandwidth, PSLR loss occurs due to the effective window being just a part of the total window.
To sum up, the loss will be found both in the range direction and the azimuth direction, which includes resolution loss and PSLR loss. If the resolution is higher, the resolution loss and PSLR loss is higher, which is unacceptable. The quantitative analysis is too complicated and it is related with  To sum up, the loss will be found both in the range direction and the azimuth direction, which includes resolution loss and PSLR loss. If the resolution is higher, the resolution loss and PSLR loss is higher, which is unacceptable. The quantitative analysis is too complicated and it is related with B r , f 0 , and azimuth resolution. Additionally, the necessity of the quantitative error analysis is not enough due to the former weighting method being easy to improve by weighting in the range frequency domain. Thus, the simulation results are illuminated in Section 4.

Weighting and Conducting a Deramping Operation in the Range Frequency Domain
Conducting the weighting operation in the range frequency domain is reasonable and easy to achieve. The purpose of this part is to analyze the influence of the former bandwidth B ar ≈ 2V s · β/λ and new bandwidth B ar ( f r ) ≈ 2V s · β/c · ( f 0 + f r ).
• The loss of resolution.
In the two-dimensional frequency domains, as Figure 8 shows, RCM does not exist, while the frequency spectrum is trapezoidal. The two-dimensional frequency spectra are lost due to the weighting operation without updating f r along the range. However, the proposed weighting method can match the frequency spectrum perfectly and lead to no loss. Therefore, the following error analysis is about conducting the weighting operation in the range frequency domain without updating f r .

Weighting and Conducting a Deramping Operation in the Range Frequency Domain
Conducting the weighting operation in the range frequency domain is reasonable and easy to achieve. The purpose of this part is to analyze the influence of the former bandwidth 2 The loss of resolution.
In the two-dimensional frequency domains, as Figure 8 shows, RCM does not exist, while the frequency spectrum is trapezoidal. The two-dimensional frequency spectra are lost due to the weighting operation without updating r f along the range. However, the proposed weighting method can match the frequency spectrum perfectly and lead to no loss. Therefore, the following error analysis is about conducting the weighting operation in the range frequency domain without updating Similarly, the loss of resolution is because of the loss of the spectrum in the two-dimensional frequency domains. Since the spectrum is trapezoidal, calculating the loss of the azimuth resolution is achievable. The parameter  is introduced, which is defined as: where  represents the ratio of the range bandwidth to the carrier frequency; this parameter can reflect the rot f error in the whole deramping operation and finally determines the loss of the azimuth resolution and PSLR.
The error caused by the wideband signal without updating r f is shown in Figure 9, the loss of resolution is shown with dash line, and the loss of azimuth bandwidth is proportionate to Similarly, the loss of resolution is because of the loss of the spectrum in the two-dimensional frequency domains. Since the spectrum is trapezoidal, calculating the loss of the azimuth resolution is achievable. The parameter α is introduced, which is defined as: where α represents the ratio of the range bandwidth to the carrier frequency; this parameter can reflect the f rot error in the whole deramping operation and finally determines the loss of the azimuth resolution and PSLR. The error caused by the wideband signal without updating f r is shown in Figure 9, the loss of resolution is shown with dash line, and the loss of azimuth bandwidth is proportionate to f rot_error / f rot_re f . Moreover, the total error is the integration along the range frequency f r , the range frequency distribution is uniform, and the error integral is expressed as: Sensors 2017, 17, 220 11 of 17 frequency distribution is uniform, and the error integral is expressed as: In Figure 9, the error caused by the wideband signal means every point is inaccurate if This error leads to two other kinds of errors: 1.
If 0 r f  , as the red line shows: , the dechirp operation exceeds the original ratio, leading to over-dechirp, and finally causes the azimuth PSLR loss, which is because the azimuth weighting window is larger than the residual bandwidth.

2.
If 0 r f  , as the purple line shows: , the dechirp operation is less than the original ratio, leading to the lack of dechirp, and finally causes the azimuth resolution loss, which is because the azimuth weighting window is shorter than the residual bandwidth.
Only by updating r f in the weighting operation could the length of the azimuth window be controlled and, finally, can precise azimuth weighting be achieved.


The loss of PSLR Similarly, the azimuth PSLR and ISLR will be lost due to the wrong boundary of the weighting operation, as Figure 7 shows. The simulation results are illuminated in Section 4.
To sum up, the loss will be found both in the range direction and the azimuth direction, which includes resolution loss and PSLR loss.
 is directly proportional to the resolution loss and shows a positive correlation with the PSLR loss. We further mention the weighting result of our method is close to the theoretical value, as presented in the experimental part. In Figure 9, the error caused by the wideband signal means every point is inaccurate if f r = 0. This error leads to two other kinds of errors: 1. If f r < 0, as the red line shows: the dechirp operation exceeds the original ratio, leading to over-dechirp, and finally causes the azimuth PSLR loss, which is because the azimuth weighting window is larger than the residual bandwidth 2. If f r > 0, as the purple line shows: the dechirp operation is less than the original ratio, leading to the lack of dechirp, and finally causes the azimuth resolution loss, which is because the azimuth weighting window is shorter than the residual bandwidth.
Only by updating f r in the weighting operation could the length of the azimuth window be controlled and, finally, can precise azimuth weighting be achieved.
• The loss of PSLR Similarly, the azimuth PSLR and ISLR will be lost due to the wrong boundary of the weighting operation, as Figure 7 shows. The simulation results are illuminated in Section 4.
To sum up, the loss will be found both in the range direction and the azimuth direction, which includes resolution loss and PSLR loss. α is directly proportional to the resolution loss and shows a positive correlation with the PSLR loss. We further mention the weighting result of our method is close to the theoretical value, as presented in the experimental part.

Experimental Results
A number of simulations have been carried out with the resolution varying from 0.2 m to 1 m, which are illustrated as follows. The imaging algorithms are the same and a −25 dB Taylor window is adopted in the following simulations.
Consider two different processing methods:

Weighting and Conducting a Deramping Operation in the Range Time Domain
To show the loss of PSLR and resolution in sub-meter resolution data processing in the range time domain [26], we simulated a single target with the resolution varying from 0. Consider two different processing methods:

Weighting and Conducting a Deramping Operation in the Range Time Domain
To show the loss of PSLR and resolution in sub-meter resolution data processing in the range time domain [26], we simulated a single target with the resolution varying from 0.2 m to 1 m. The carrier frequency 0 5.4 GHz f  . The result of the simulations are shown in Figures 10 and 11. In Figures 10 and 11, the loss of resolution and PSLR are both in the range and azimuth directions because of two-dimensional coupling. As we can see, the loss of resolution and PSLR are both intolerable with the increase of resolution.

Weighting and Conducting a Deramping Operation in the Range Time Domain
To show the loss of PSLR and resolution in sub-meter resolution data processing in the range time domain [26], we simulated a single target with the resolution varying from 0.2 m to 1 m. The carrier frequency 0 5.4 GHz f  . The result of the simulations are shown in Figures 10 and 11. In Figures 10 and 11, the loss of resolution and PSLR are both in the range and azimuth directions because of two-dimensional coupling. As we can see, the loss of resolution and PSLR are both intolerable with the increase of resolution. In Figures 10 and 11, the loss of resolution and PSLR are both in the range and azimuth directions because of two-dimensional coupling. As we can see, the loss of resolution and PSLR are both intolerable with the increase of resolution.
We further note if the carrier frequency is increased, the loss of resolution and PSLR are both decreased.

Weighting and Conducting a Deramping Operation in the Range Frequency Domain
To show the loss of PSLR and resolution in sub-meter resolution data processing in the range frequency domain, many experiments were carried out. We simulated a single target with the resolution varying from 0.2 m to 1 m. The result of the simulations are illustrated in Figures 12 and 13. decreased.

Weighting and Conducting a Deramping Operation in the Range Frequency Domain
To show the loss of PSLR and resolution in sub-meter resolution data processing in the range frequency domain, many experiments were carried out. We simulated a single target with the resolution varying from 0. In Figures 12 and 13, the loss of resolution and PSLR are given in the azimuth direction. The loss in the range direction is so small that we can neglect it. As we can see, with the increase of  , the loss of resolution and PSLR are both intolerable. Additionally the theoretical value is shown by the red line, which is close to the simulated values. However, the performance of the proposed weighting decreased.

Weighting and Conducting a Deramping Operation in the Range Frequency Domain
To show the loss of PSLR and resolution in sub-meter resolution data processing in the range frequency domain, many experiments were carried out. We simulated a single target with the resolution varying from 0.2 m to 1 m. The result of the simulations are illustrated in Figures 12 and 13. In Figures 12 and 13, the loss of resolution and PSLR are given in the azimuth direction. The loss in the range direction is so small that we can neglect it. As we can see, with the increase of  , the loss of resolution and PSLR are both intolerable. Additionally the theoretical value is shown by the red line, which is close to the simulated values. However, the performance of the proposed weighting In Figures 12 and 13, the loss of resolution and PSLR are given in the azimuth direction. The loss in the range direction is so small that we can neglect it. As we can see, with the increase of α, the loss of resolution and PSLR are both intolerable. Additionally the theoretical value is shown by the red line, which is close to the simulated values. However, the performance of the proposed weighting method is hardly influenced by α, which means the proposed method is suitable for the high-resolution spotlight sliding method azimuth weighting.
We further mention if the carrier frequency is increased, the loss of resolution and PSLR are both decreased. Namely, the loss of resolution is inversely proportional to the carrier frequency.

The Simulation at 0.3 m Resolution
In order to verify the weighting algorithm performance at high-resolution, a 0.3 m resolution sliding spotlight mode simulation is conducted as a specific example. Three weighting methods are Sensors 2017, 17, 220 14 of 17 adopted in the imaging simulation with the same system parameters as in Table 1. The imaging algorithms are all NCSA (Nonlinear Chirp Scaling Algorithm) while the difference between a, b, and c procedures in Figure 14 are the weighting method. algorithms are all NCSA (Nonlinear Chirp Scaling Algorithm) while the difference between a, b, and c procedures in Figure 14 are the weighting method.   Now the two-dimensional spectrum losses are shown in Figure 14. Weighting in the range time domain losses a lot of spectrum, which leads to great resolution loss. Weighting in the range frequency domain without updating f r also leads to some resolution losses. The proposed weighting method in the range frequency domain with updating f r performs well. The comparison of the three azimuth weighting method is shown in Table 2 and the contour of the three imaged point targets is shown in Figure 15. In Figure 15 and Table 2, the losses exist both in weighting methods a and b, which are unacceptable, while the proposed weighting method performs well. The data are in the two-dimensional frequency domain after a precise deramping method and followed by 2-dimensional compensation naturally, which means the computational load of the proposed weighting method is not increased.  In Figure 15 and Table 2, the losses exist both in weighting methods a and b, which are unacceptable, while the proposed weighting method performs well. The data are in the twodimensional frequency domain after a precise deramping method and followed by 2-dimensional compensation naturally, which means the computational load of the proposed weighting method is not increased.
To compare the influence of the area target, the simulation with the same parameters in Table 1 is conducted, and the result is shown in Figure 16. The boundary in the left picture is clearer than the right picture in Figure 16. Improving the quality of the focusing images without any extra computation burden is significant.

Conclusions
The spaceborne sliding-spotlight mode could be widely adopted for future spaceborne remote sensing. However, the wideband signal will influence its weighting performance when a two-step algorithm is adopted in the system. Based on classical azimuth preprocessing methods, carefully analysis of the signal character, and proposing the azimuth weighting method in the range frequency To compare the influence of the area target, the simulation with the same parameters in Table 1 is conducted, and the result is shown in Figure 16.  In Figure 15 and Table 2, the losses exist both in weighting methods a and b, which are unacceptable, while the proposed weighting method performs well. The data are in the twodimensional frequency domain after a precise deramping method and followed by 2-dimensional compensation naturally, which means the computational load of the proposed weighting method is not increased.
To compare the influence of the area target, the simulation with the same parameters in Table 1 is conducted, and the result is shown in Figure 16. The boundary in the left picture is clearer than the right picture in Figure 16. Improving the quality of the focusing images without any extra computation burden is significant.

Conclusions
The spaceborne sliding-spotlight mode could be widely adopted for future spaceborne remote sensing. However, the wideband signal will influence its weighting performance when a two-step algorithm is adopted in the system. Based on classical azimuth preprocessing methods, carefully analysis of the signal character, and proposing the azimuth weighting method in the range frequency The boundary in the left picture is clearer than the right picture in Figure 16. Improving the quality of the focusing images without any extra computation burden is significant.

Conclusions
The spaceborne sliding-spotlight mode could be widely adopted for future spaceborne remote sensing. However, the wideband signal will influence its weighting performance when a two-step algorithm is adopted in the system. Based on classical azimuth preprocessing methods, carefully analysis of the signal character, and proposing the azimuth weighting method in the range frequency domain, updating f r , then based on the weighting requirement we compare the former weighting method and the proposed weighting method with the varying resolution. Finally, the simulation results validate the presented analysis in the paper.