Home | About | Journals | Submit | Contact Us | Français |

**|**Sensors (Basel)**|**v.17(2); 2017 February**|**PMC5336078

Formats

Article sections

- Abstract
- 1. Introduction
- 2. Signal Characteristics Analysis
- 3. Precious Azimuth Weighting Method in a Two-Step Algorithm
- 4. Experimental Results
- 5. Conclusions
- References

Authors

Related links

Sensors (Basel). 2017 February; 17(2): 220.

Published online 2017 January 24. doi: 10.3390/s17020220

PMCID: PMC5336078

Jonathan Li, Academic Editor

Beijing Key Laboratory of Embedded Real-time Information Processing Technology, School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China; Email: nc.ude.tib@943010213 (F.X.); Email: moc.621@9040110211bx (B.X.); Email: nc.ude.tib@gnetgnol (T.L.)

Received 2016 November 7; Accepted 2017 January 20.

Copyright © 2017 by the authors.

Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

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.

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 (SPECAN), a highly accurate algorithm for spotlight data is proposed [8,21], whose azimuth weighting is not considered and needs sub-aperture processing. Back-projection (BP) and fast BP [22,23,24,25] are other ways to focus SAR echoes precisely without considering the work mode, but the sidelobe suppression of several work modes, such as the sliding spotlight mode and TOPS mode, need to be considered. Additionally, the computational load may be further improved. From former analysis, it is clearly that the azimuth weighting method in the sliding spotlight mode is not carefully discussed. Therefore, based on a two-step algorithm, an azimuth weighting method is proposed, which is suitable for high-resolution data processing.

The paper is organized as follows: In the second section, the signal model is established, which includes the geometry and echo model. In the third section, based on curved geometry, the algorithm of weighting windows in frequency is declared. In addition, the comparison of the former weighting method and the proposed method is illuminated. In the fourth section, the simulation results validate the presented analysis. Conclusions are drawn in the fifth section.

The planar curved imaging geometry of the sliding-spotlight mode is shown in Figure 1, which indicates that the azimuth beam is steered from fore to aft at a constant rotation rate as:

$${\omega}_{r}=\frac{d\theta}{dT},$$

(1)

where $\theta $ is the instantaneous squint angle and $T$ is the whole acquisition interval, and $R$ and $r$ are the slant ranges from the flight path to the scene center and the virtual rotation center to the scene center, respectively. $O$ is the geocenter, $\beta $ is the azimuth beam width, $X\approx \beta R$ is the width in azimuth direction, and ${X}_{f}$ represents the width of the scene. $\alpha $ is the rotation angle.

Space-borne sliding-spotlight mode planar imaging curved Earth geometry. A and B are the orbit start and end. The orange region represents the beam. The blue region represents the Earth’s interior.

For azimuth beam scanning at a constant rotation rate leading to ${\omega}_{r}=\theta /T$, the steering factor $A$ is defined as [14]:

$$A=\frac{r}{R+r},$$

(2)

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]:

$${\rho}_{AZ}=\frac{{D}_{a}}{2}\cdot \frac{{V}_{g}}{{V}_{s}}\cdot {\gamma}_{w,a}\cdot A=\frac{{D}_{a}}{2}\cdot A,$$

(3)

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 ${\gamma}_{w,a}$ is the impulse response width broadening factor due to an azimuth processing window. The approximation ${V}_{g}/{V}_{s}{\gamma}_{w,a}\approx 1$ is reasonable in the majority case.

If the transmitter illuminates the target scene with a baseband chirp signal $p(t)$ and a point target locates at $\left({x}_{0},R\right)$, the echo of transmitted pulse can be expressed as:

$$\begin{array}{cc}\hfill s\left(\tau ;{t}_{a},x\right)& =rect\left[\frac{\tau -\frac{2R\left({t}_{a}\right)}{c}}{{T}_{p}}\right]\cdot \mathrm{exp}\left[-j\frac{4\pi R\left({t}_{a}\right)}{\lambda}\right]\cdot \mathrm{exp}\left[j\pi {K}_{r}{\left(\tau -\frac{2R\left({t}_{a}\right)}{c}\right)}^{2}\right]\hfill \\ & \cdot \text{}rect\left[\frac{{V}_{f}{t}_{a}-{x}_{0}}{X}\right]\cdot rect\left[\frac{x}{{X}_{f}}\right]\hfill \end{array},$$

(4)

where $\lambda =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. $\tau $ 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\left({t}_{a}\right)$ can be written as:

$$R\left({t}_{a}\right)\approx \sqrt{{R}_{0}^{2}+{V}_{r}^{2}\cdot {{t}_{a}}^{2}},$$

(5)

where ${V}_{r}\approx \sqrt{{V}_{g}\cdot {V}_{s}}$. ${V}_{r}$ is usually used for the imaging focus and ${V}_{s}$ is used for calculating the Doppler bandwidth.

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.

C-band VV-polarized image of Hong Kong obtained by the GF-3 system in sliding spotlight mode in 2016. (**A**) Using a weighting operation, and (**B**) without a weighting operation. Azimuth resolution is about 1 m. The extension of the area is of about 8 km × **...**

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.

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}\left({t}_{a}\right)$ can be expressed as:

$$\begin{array}{cc}{R}_{c}\left({t}_{a}\right)\hfill & =\sqrt{{\left(R+r\right)}^{2}+{V}_{s}^{2}\cdot {{t}_{a}}^{2}}\hfill \\ & \approx \left(R+r\right)+\frac{{V}_{s}^{2}\cdot {{t}_{a}}^{2}}{2\left(R+r\right)}+\frac{{V}_{s}^{4}\cdot {{t}_{a}}^{4}}{8{\left(R+r\right)}^{3}}\hfill \end{array}.$$

(6)

Next, the ${f}_{dc}$ in the virtual point can be calculated via ${f}_{dc\_rot}\left({t}_{a}\right)=-2/\lambda \cdot d{R}_{c}\left({t}_{a}\right)/d{t}_{a}$.

Then the azimuth variation rate ${f}_{rot}$ causing by azimuth beam rotation can be expressed as:

$${f}_{rot}\left({t}_{a}\right)=\frac{d{f}_{dc\_rot}\left({t}_{a}\right)}{d{t}_{a}}=-\frac{2}{\lambda}\cdot \left[\frac{{V}_{s}^{2}}{\left(R+r\right)}+\frac{3{V}_{s}^{4}\cdot {{t}_{a}}^{2}}{2{\left(R+r\right)}^{3}}\right].$$

(7)

Under the typical low orbit parameter and carrier frequency is C band, ${f}_{rot}$ can be expressed approximately as:

$$\begin{array}{cc}\hfill {f}_{rot}& \approx -\frac{2}{\lambda}\cdot \frac{{V}_{s}^{2}}{\left(R+r\right)}\hfill \\ & =-\frac{2{V}_{s}}{\lambda}\cdot \frac{2\theta}{T}\hfill \\ & =-\frac{2{V}_{r}}{\lambda}\cdot \frac{2{\theta}_{r}}{T}\hfill \end{array},$$

(8)

where ${\theta}_{r}$ is equivalent instantaneous squint angle. Thus, the azimuth bandwidth caused by the beam rotation can be expressed as:

$$\begin{array}{cc}\hfill {B}_{rot}& =abs\left[{f}_{rot}\cdot T\right]\hfill \\ & =\frac{2{V}_{s}}{\lambda}\cdot 2\theta \hfill \\ & =\frac{2{V}_{r}}{\lambda}\cdot 2{\theta}_{r}\hfill \end{array}.$$

(9)

Compared with the strip-map mode, the azimuth original bandwidth caused by the beam width is:

$${B}_{beam}=\frac{2{V}_{s}}{\lambda}\cdot \beta .$$

(10)

Finally, the azimuth total Doppler bandwidth can be expressed as:

$${B}_{total}={B}_{rot}+{B}_{beam}.$$

(11)

Figure 3 shows the Doppler time-frequency relationship.

Where ${f}_{dr}=2{V}_{r}^{2}/\left(\lambda \cdot {R}_{0}\right)$ 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:

$$\begin{array}{cc}\hfill g\left({t}_{a}\right)& =\mathrm{exp}\left(-j\pi \cdot {f}_{rot}\cdot {{t}_{a}}^{2}\right)\hfill \\ & =exp\left(j\pi \cdot \frac{2{V}_{s}}{\lambda}\cdot \frac{2\theta}{T}\cdot {{t}_{a}}^{2}\right)\hfill \\ & =exp\left(j\pi \cdot \frac{2{V}_{s}^{2}}{\lambda \left(R+r\right)}\cdot {{t}_{a}}^{2}\right)\hfill \end{array},$$

(12)

While conducting the azimuth weighting processing in range time domain, the azimuth echo preprocessing can be accomplished by employing:

$$\begin{array}{c}c\left({t}_{\tau},{t}_{a}\right)\hfill \\ =S\left({t}_{\tau},{t}_{a}\right)\otimes g\left({t}_{a}\right)\hfill \\ ={\displaystyle \int S\left({t}_{\tau},z\right)\cdot g\left({t}_{a}-z\right)dz}\hfill \end{array}.$$

(13)

Substituting Equation (12) into Equation (13), Equation (13) can be rewritten as:

$$c\left({t}_{\tau},{t}_{a}\right)=\underset{\begin{array}{l}phase\hfill \\ compensation\hfill \end{array}}{\underbrace{\mathrm{exp}\left[j2\pi \cdot \frac{{V}_{s}^{2}{t}_{a}^{2}}{\lambda \left(R+r\right)}\right]}}\underset{FT}{\underbrace{\cdot {\displaystyle \int \underset{dechirp}{\underbrace{S\left({f}_{\tau},{t}_{a}\right)\cdot \mathrm{exp}\left[j2\pi \cdot \frac{{V}_{s}^{2}\cdot {z}^{2}}{\lambda \left(R+r\right)}\right]}}\cdot \mathrm{exp}\left[-j2\pi \cdot \frac{2{V}_{s}^{2}\cdot z\cdot {t}_{a}}{\lambda \left(R+r\right)}\right]dz}}}.$$

(14)

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.

The Doppler time-frequency relationship after dechirp processing, the signal totally overlaps in the support domain which can be weighted. The time-frequency relationship is under the condition of the narrowband signal.

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 $\lambda \approx 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:

$${\lambda}_{equal}\left({f}_{\tau}\right)=\frac{c}{{f}_{0}+{f}_{\tau}}$$

(15)

The parameter ${\lambda}_{equal}\left({f}_{\tau}\right)$ is defined in the range frequency domain, which is related with ${f}_{\tau}$. 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:

$$\begin{array}{cc}\hfill c\left({f}_{\tau},{t}_{a}\right)& =S\left({f}_{\tau},{t}_{a}\right)\otimes g\left({t}_{a}\right)\hfill \\ & ={\displaystyle \int S\left({f}_{\tau},z\right)\cdot g\left({t}_{a}-z\right)dz}\hfill \end{array}.$$

(16)

Combining Equations (14)–(16), then:

$$\begin{array}{cc}\hfill c\left({f}_{\tau},{t}_{a}\right)& =\underset{\begin{array}{l}phase\hfill \\ compensation\hfill \end{array}}{\underbrace{\mathrm{exp}\left[j2\pi \cdot \frac{{V}_{s}^{2}\cdot {t}_{a}^{2}}{\frac{c}{\left({f}_{0}+{f}_{\tau}\right)}\cdot \left(R+r\right)}\right]}}\hfill \\ & \underset{FT}{\underbrace{\cdot {\displaystyle \int \underset{dechirp}{\underbrace{S\left({f}_{\tau},{t}_{a}\right)\cdot \mathrm{exp}\left[j2\pi \cdot \frac{{V}_{s}^{2}\cdot {z}^{2}}{\frac{c}{\left({f}_{0}+{f}_{\tau}\right)}\cdot \left(R+r\right)}\right]}}\cdot \mathrm{exp}\left[-j2\pi \cdot \frac{2{V}_{s}^{2}\cdot z\cdot {t}_{a}}{\frac{c}{\left({f}_{0}+{f}_{\tau}\right)}\cdot \left(R+r\right)}\right]dz}}}\hfill \end{array}.$$

(17)

is derived. Implementing Equation (17) in the discrete domain, it can be written as:

$$\begin{array}{cc}\hfill c\left({f}_{\tau},\Delta {x}_{1}\right)& =\mathrm{exp}\left[j2\pi \cdot \frac{{\left(n\cdot \Delta {x}_{2}\right)}^{2}}{\frac{c}{\left({f}_{0}+{f}_{\tau}\right)}\cdot \left(R+r\right)}\right]\hfill \\ & \cdot {\displaystyle \sum _{i=-\frac{{N}_{a}}{2}}^{i=\frac{{N}_{a}}{2}}S\left({f}_{\tau},\Delta {x}_{1}\right)}\cdot \mathrm{exp}\left[j2\pi \cdot \frac{{\left(i\cdot \Delta {x}_{3}\right)}^{2}}{\frac{c}{\left({f}_{0}+{f}_{\tau}\right)}\cdot \left(R+r\right)}\right]\cdot \mathrm{exp}\left[-j2\pi \frac{2\cdot \Delta {x}_{2}\cdot \Delta {x}_{3}}{\frac{c}{\left({f}_{0}+{f}_{\tau}\right)}\cdot \left(R+r\right)}\right]\hfill \end{array},$$

(18)

where:

$$\begin{array}{l}\Delta {x}_{1}=\frac{{V}_{r}}{PRF}\hfill \\ \Delta {x}_{2}=\frac{{V}_{s}}{PR{F}_{up}}\hfill \\ \Delta {x}_{3}=\frac{{V}_{s}}{PRF}\hfill \end{array},$$

(19)

then Equation (18) can be expressed as:

$$\begin{array}{cc}\hfill c\left(n\Delta {x}^{\u2033}\right)& =\mathrm{exp}\left[j2\pi \cdot \frac{{\left(n\cdot \Delta {x}_{2}\right)}^{2}}{\frac{c}{\left({f}_{0}+{f}_{\tau}\right)}\cdot \left(R+r\right)}\right]\hfill \\ & \cdot {\displaystyle \sum _{i=-\frac{P}{2}}^{i=\frac{P}{2}}S\left({f}_{\tau},\Delta {x}_{1}\right)}\cdot \mathrm{exp}\left[j2\pi \cdot \frac{{\left(i\cdot \Delta {x}_{3}\right)}^{2}}{\frac{c}{\left({f}_{0}+{f}_{\tau}\right)}\cdot \left(R+r\right)}\right]\cdot \mathrm{exp}\left[-j2\pi \cdot \frac{2\cdot \Delta {x}_{2}\cdot \Delta {x}_{3}}{\frac{c}{\left({f}_{0}+{f}_{\tau}\right)}\cdot \left(R+r\right)}\right]\hfill \\ & \approx \mathrm{exp}\left[j2\pi \cdot \frac{{\left(n\cdot \Delta {x}_{2}\right)}^{2}}{\frac{c}{\left({f}_{0}+{f}_{\tau}\right)}\cdot \left(R+r\right)}\right]\cdot DFT\left\{S\left({f}_{\tau},\Delta {x}_{1}\right)\cdot \mathrm{exp}\left[j2\pi \cdot \frac{{\left(i\cdot \Delta {x}_{3}\right)}^{2}}{\frac{c}{\left({f}_{0}+{f}_{\tau}\right)}\cdot \left(R+r\right)}\right]\right\}\hfill \end{array},$$

(20)

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:

$${B}_{ar}\approx \frac{2{V}_{s}\cdot \beta}{c}\cdot \left({f}_{0}+{f}_{r}\right).$$

(21)

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 $PRF>2{V}_{s}/\left(c/{f}_{0}\right)\cdot \beta +{V}_{s}/{\rho}_{r}\cdot \beta $, where ${\rho}_{r}=c/\left(2{B}_{r}\right)$ is the range 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.

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:

• 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.

• 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.

Different errors in the weighting operation. The purple line shows the window and the yellow region is the effective bandwidth that needs to be weighted. (**a**) The length of the window is shorter than the effective bandwidth; (**b**) The length of the 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 ${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.

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}\approx 2{V}_{s}\cdot \beta /\lambda $ and new bandwidth ${B}_{ar}\left({f}_{r}\right)\approx 2{V}_{s}\cdot \beta /c\cdot \left({f}_{0}+{f}_{r}\right)$.

• 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}$.

Conducting the weighting operation in the range frequency domain. The figure shows the spectrum in the two-dimensional frequency domain after a deramping operation. (**a**) Using the former weighting bandwidth${B}_{ar}\approx 2{V}_{s}\cdot \beta /\lambda $, the **...**

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 $\alpha $ is introduced, which is defined as:

$$\alpha =\frac{\raisebox{1ex}{${B}_{r}$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}{{f}_{0}}=\frac{{f}_{rot\_error}}{{f}_{rot\_ref}}=\frac{\frac{2{V}_{s}}{\raisebox{1ex}{$2c$}\!\left/ \!\raisebox{-1ex}{${B}_{r}$}\right.}\cdot \frac{2\theta}{T}}{\frac{2{V}_{s}}{\raisebox{1ex}{$c$}\!\left/ \!\raisebox{-1ex}{${f}_{0}$}\right.}\cdot \frac{2\theta}{T}}.$$

(22)

where $\alpha $ 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\_ref}$. 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:

$$mean(\mathrm{error})=\frac{\frac{1}{2}\alpha {\displaystyle \underset{0}{\overset{\frac{1}{2}{B}_{r}}{\int}}d{f}_{r}}}{{B}_{r}}=\frac{1}{4}\alpha .$$

(23)

After dechirp processing the Doppler time-frequency relationship in the range frequency domain, which shows the error caused by the wideband signal without updating ${f}_{r}$, the max error is caused by ${f}_{r}=\pm {B}_{r}/2$. The red line shows the Doppler time-frequency **...**

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:

${f}_{rot}\left({f}_{r}\right)-{f}_{rot}\left(0\right)=-2{V}_{s}\cdot \left({f}_{0}+{f}_{r}\right)/c\cdot 2\theta /T+2{V}_{s}\cdot \left({f}_{0}\right)/c\cdot 2\theta /T>0$, 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:

${f}_{rot}\left({f}_{r}\right)-{f}_{rot}\left(0\right)=-2{V}_{s}\cdot \left({f}_{0}+{f}_{r}\right)/c\cdot 2\theta /T+2{V}_{s}\cdot \left({f}_{0}\right)/c\cdot 2\theta /T<0$, 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. $\alpha $ 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.

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:

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 ${f}_{0}=5.4\text{}\mathrm{GHz}$. The result of the simulations are shown in Figure 10 and Figure 11.

The relationship between PSLR loss and range bandwidth weighting in the range time domain. The blue line shows the PSLR loss in the azimuth direction and the red line shows the PSLR loss in the range direction. The numbers beside the line show the two-dimensional **...**

The relationship between resolution loss and range bandwidth weighting in the range time domain. The blue line shows the resolution loss in the azimuth direction and the red line shows the resolution loss in the range direction. The numbers beside the **...**

In Figure 10 and Figure 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.

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 Figure 12 and Figure 13.

The relationship between PSLR loss and $\alpha $. The blue line shows the PSLR loss in the azimuth direction without considering the variation of ${f}_{r}$, the black line shows the PSLR loss in the azimuth direction with considering the variation of ${f}_{r}$. The **...**

The relationship between resolution loss and $\alpha $. The blue line shows the resolution loss in the azimuth direction without considering the variation of ${f}_{r}$. The red line shows the theoretic resolution loss without considering the variation of ${f}_{r}$ **...**

In Figure 12 and Figure 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 $\alpha $, 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 $\alpha $, 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.

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 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.

The comparison of two dimensional spectrum of a single target at 0.3 m resolution using different azimuth weighting methods. (**a**) Conducting azimuth weighting in the range time domain; (**b**) Conducting azimuth weighting in the range frequency domain without **...**

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.

Simulated image using different window-adding methods. Contour lines are just an example to show the result. (**a**) Conducting azimuth weighting in the range time domain; (**b**) Conducting azimuth weighting in the range frequency domain without updating ${f}_{r}$ **...**

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.

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.

A 0.3 m resolution sliding spotlight image. (**A**) Using the proposed weighting method; (**B**) Using the azimuth weighting method in the range time domain. The extension of the area is of about 1 km × 1 km.

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.

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.

This work was supported by the National Natural Science Foundation of China under Grant Nos. 61370017, 61225005, 61120106004 and the project of China high resolution earth observation system under Grant Nos. 41-Y20A13-9001-15/16, 12-Y20A15-9001-15/16.

Author Contributions

Feng Xiao conceived the idea and designed the experiments; Zegang Ding analyzed the data; Bin Xiong wrote the paper; Teng Long supervised the research, including the experiments and development.

Conflicts of Interest

The authors declare no conflict of interest.

1. Lanari R., Zoffoli S., Sansosti E., Fornaro G., Serafino F. New Approach for Hybrid Strip-map/Spotlight SAR Data Focusing. IEEE Proc. Radar Sonar Navig. 2001;148:363–372. doi: 10.1049/ip-rsn:20010662. [Cross Ref]

2. Lanari R., Franceschetti P., Tesauro M., Sansosti E. Spotlight SAR Image Generation Based on Strip Mode Focusing Techniques; Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS); New York, NY, USA. 28 June–2 July 1999; pp. 1761–1763.

3. Belcher D.P., Baker C.J. High Resolution Processing of Hybrid Stripmap/Spotlight Mode SAR. Proc. Inst. Elect. Eng. Radar Sonar Navig. 1996;143:366–374. doi: 10.1049/ip-rsn:19960790. [Cross Ref]

4. Prats-Iraola P., Scheiber R., Rodriguez-Cassola M., Mittermayer J. On the Processing of Very High Resolution Spaceborne SAR Data. IEEE Trans. Geosci. Remote Sens. 2014;52:6003–6016. doi: 10.1109/TGRS.2013.2294353. [Cross Ref]

5. Carrara W.G., Majewski R.M., Goodman R.S. Spotlight Synthetic Apecture Radar: Signal Processing Algorithms. Artech House; Norwood, MA, USA: 1995.

6. Mittermayer J., Lord R., Borner E. Sliding Spotlight SAR Processing for TerraSAR-X Using a New Formulation of the Extended Chirp Scaling Algorithm; Proceedings of the IEEE International Geoscience and Remote Sensing Symposium(IGARSS); Toulouse, France. 21–25 July 2003; pp. 1462–1464.

7. Lanari R., Tesauro M., .Sansosti E., Fornaro G. Spotlight SAR Data Focusing Based on a Two-Step Processing Approach. IEEE Trans. Geosci. Remote Sens. 2001;39:1993–2003. doi: 10.1109/36.951090. [Cross Ref]

8. Mittermayer J., Moreira A., Loffeld O. Spotlight SAR Data Processing Using the Frequency Scaling Algorithm. IEEE Trans. Geosci. Remote Sens. 1999;37:2198–2214. doi: 10.1109/36.789617. [Cross Ref]

9. Wu Y., Sun G.C., Yang C., Yang J., Xing M.D., Bao Z. Processing of Very High Resolution Spaceborne Sliding Spotlight SAR Data Using Velocity Scaling. IEEE Trans. Geosci. Remote Sens. 2016;54:1505–1518. doi: 10.1109/TGRS.2015.2481923. [Cross Ref]

10. Fornaro G., Lanari R., .Sansosti E., Tesauro M. A Two-step Spotlight SAR Data Focusing Approach; Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS); Honolulu, HI, USA. 24–28 July 2000; pp. 84–86.

11. Liu F.F., Ding Z.G., Zeng T., Long T. Performance Analysis of Two-Step Algorithm in Sliding Spotlight Space-borne SAR; Proceedings of the IEEE Radar Conference; Arlington, TX, USA. 10–14 May 2010; pp. 965–968.

12. Kuang H., Chen J., Yang W. A refined two-step algorithm for high resolution spaceborne SAR with squinted sliding spotlight mode; Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS); Milano, Italy. 26–31 July 2015; pp. 4476–4479.

13. Stangl M., Werninghaus R., Schweizer B., Fischer C. TerraSAR-X Technologies and First Results. IEEE Proc. Radar Sonar Navig. 2006;153:86–95. doi: 10.1049/ip-rsn:20045119. [Cross Ref]

14. Wang P.B., Zhou Y.Q., Chen J., Li C.S., Yu Z., Min H. A Deramp Frequency Scaling Algorithm for Processing Space-Borne Spotlight SAR Data; Proceedings of the IEEE International Symposium on Geoscience and Remote Sensing; Denver, CO, USA. 31 July–4 August 2006; pp. 3148–3151.

15. Yang W., Li C.S., Chen J., Wang P.B. Extended three-step focusing algorithm for sliding spotlight and tops data image formation; Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS); Vancouver, BC, Canada. 24–29 July 2011; pp. 479–482.

16. Buckreuss S., Werninghaus R., Pitz W. The German Satellite Mission TerraSAR-X. IEEE Aerosp. Electron. Syst. Mag. 2009;24:4–9. doi: 10.1109/MAES.2009.5344175. [Cross Ref]

17. Mora O., Granda J., Biescas E., Urdiroz A. TerraSAR-X High Resolution SAR Data: Ground Motion and Mapping Applications for Infrastructure, Oil & Gas and Public Health Domain; Proceedings of the European Conference on Synthetic Aperture Radar (EUSAR); Friedrichshafen, Germany. 2–5 June 2008; pp. 1–4.

18. Breit H., Fritz T., Balss U., Lachaise M. TerraSAR-X SAR Processing and Products. IEEE Trans. Geosci. Remote Sens. 2010;48:727–740. doi: 10.1109/TGRS.2009.2035497. [Cross Ref]

19. Boerner E., Lord R., Mittermayer J., Bamler R. Evaluation of TerraSAR-X Spotlight Processing Accuracy Based on a New Spotlight Raw Data Simulator; Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS); Toulouse, France. 21–25 July 2003; pp. 1323–1325.

20. Prats P., Scheiber R., Mittermayer J., Meta A., Moreira A. Processing of Sliding Spotlight and TOPS SAR Data Using Baseband Azimuth Scaling. IEEE Trans. Geosci. Remote Sens. 2010;48:770–780. doi: 10.1109/TGRS.2009.2027701. [Cross Ref]

21. Mittermayer J., Moreira A., Loffeld O. Spotlight Processing of Wide-Beam Stripmap SAR Data Using the Frequency Scaling Algorithm; Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS); Seattle, WA, USA. 6–10 July 1998; pp. 1177–1179.

22. Desai M.D., Jenkins W.K. Convolution Backprojection Image Reconstruction for Spotlight Mode Synthetic Aperture Radar. IEEE Trans. Image Proc. 1992;1:505–517. doi: 10.1109/83.199920. [PubMed] [Cross Ref]

23. Yegulalp A.F. Fast Backprojection Algorithm for Synthetic Aperture Radar; Proceedings of the Record of the 1999 IEEE Radar Conference; Waltham, MA, USA. 20–24 April 1999; pp. 60–65.

24. Zhang L., Li H.-L., Qiao Z.-J., Xu Z.-W. A Fast BP Algorithm with Wavenumber Spectrum Fusion for High-Resolution Spotlight SAR Imaging. IEEE Geosci. Remote Sens. Lett. 2014;11:1460–1464. doi: 10.1109/LGRS.2013.2295326. [Cross Ref]

25. Gaibel A., Boag A. Back-projection SAR imaging using FFT; Proceedings of the European Radar Conference (EuRAD); London, UK. 5–7 October 2016; pp. 69–72.

26. Shu B.Z., Xu Y.H., Ding Z.G., Zhu Y. An Azimuth Weighting Method of Slide-spotlight SAR Based on Subaperture Chirp Scaling Algorithm; Proceedings of the IET International Radar Conference; Hangzhou, China. 14–16 October 2015; pp. 1–5.

27. Cumming I.G., Wong F.H. Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation. Artech House; Norwood, MA, USA: 2005.

28. Raney K., Runge H., Bamler R., Cumming I., Wong F. Precision SAR Processing Using Chirp Scaling. IEEE Trans. Geosci. Remote Sens. 1994;32:786–799. doi: 10.1109/36.298008. [Cross Ref]

29. Davidson G.W., Cumming I.G., Ito M.R. A Chirp Scaling Approach for Processing Squint Mode SAR Data. IEEE Trans. Aerosp. Electron. Syst. 1996;32:121–133. doi: 10.1109/7.481254. [Cross Ref]

30. An D.X., Huang X.T., Jin T., Zhou Z.M. Extended Nonlinear Chirp Scaling Algorithm for High-Resolution Highly Squint SAR Data Focusing. IEEE Trans. Geosci. Remote Sens. 2012;50:3595–3609. doi: 10.1109/TGRS.2012.2183606. [Cross Ref]

31. Sun G.C., Jiang X.W., Xing M.D., Qiao Z.J., Wu Y.R., Bao Z. Focus Improvement of Highly Squinted Data Based on Azimuth Nonlinear Scaling. IEEE Trans. Geosci. Remote Sens. 2011;49:2308–2322. doi: 10.1109/TGRS.2010.2102040. [Cross Ref]

Articles from Sensors (Basel, Switzerland) are provided here courtesy of **Multidisciplinary Digital Publishing Institute (MDPI)**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |