Next Article in Journal
Applicability of Precipitation Products in the Endorheic Basin of the Yellow River under Multi-Scale in Time and Modality
Previous Article in Journal
Hyperspectral Image Classification Based on Mutually Guided Image Filtering
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel SV-PRI Strategy and Signal Processing Approach for High-Squint Spotlight SAR

1
Department of Space Microwave Remote Sensing System, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100190, China
2
School of Electronic, Electrical and Communication Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Submission received: 26 January 2024 / Revised: 27 February 2024 / Accepted: 28 February 2024 / Published: 29 February 2024

Abstract

:
High-resolution and high-squint spaceborne spotlight synthetic aperture radar (SAR) has significant potential for extensive application in remote sensing, but its swath width effectiveness is constrained by a critical factor: severe range cell migration (RCM). To address this, pulse repetition interval (PRI) variation offers a practical scheme for raw data reception. However, the current designs for continuously varying PRI (CV-PRI) exhibit high complexity in engineering. In response to the issue, this paper proposes a novel strategy of stepwise varying PRI (SV-PRI), which demonstrates higher reconstruction accuracy compared with CV-PRI. Furthermore, confronting the azimuth non-uniform sampling characteristics induced by the PRI variation, this paper introduces a complete uniform reconstruction processing based on the azimuth partitioning methodology, which effectively alleviates the inherent contradiction between resolution and swath width. The processing flow, utilizing the temporal point remapping (TPR) concept, ensures the uniformity and coherence of dataset partitioning and reassembly in the context of the interpolation on non-uniform grids. Finally, according to the simulation results, the point target data, processed through the processing flow proposed in this study, have demonstrated effective focusing results.

1. Introduction

Synthetic aperture radar (SAR) is a powerful remote sensing technique, capable of providing high-resolution microwave images of the Earth’s surface regardless of weather conditions and sunlight illumination [1]. The growing demand for Earth observation necessitates improvements in radar technology to achieve more refined and informative imagery, making high-resolution imaging a critical criterion [2,3,4,5,6,7,8,9]. As a way to meet the requirement of a high-resolution azimuth, the spotlight mode has gained increasing prominence in recent years. In the spotlight mode, the spaceborne SAR system controls the direction of the radar antenna’s beam to continuously and extensively illuminate a specific observation area. This operational mode facilitates the acquisition of detailed terrain feature information and enables a more precise inversion of terrain parameters [2,8]. Based on the current in-orbit satellite data, spotlight SAR systems have demonstrated high-precision imaging capabilities with a sub-meter resolution, as exemplified by the prominent German TerraSAR-X satellite [3]. For specific regions of interest, the high-squint mode is imperative. The inclusion of a large angle significantly augments the acquisition of nuanced information regarding target characteristics and angle-dependent backscattering coefficients, enhancing the SAR system’s adaptability and observational capacity [9,10,11,12,13,14,15,16,17,18]. However, compared with the conventional mode, the imaging data reception in the high-squint spotlight SAR is constrained by severe range cell migration (RCM). This is a limitation arising from the inherent contradiction between resolution and swath width. Therefore, the acquisition of high-quality, focused images necessitates a trade-off, often resulting in reduced scanning width. In the case, the pulse repetition interval (PRI) variation, that is, PRI is dynamically adjusted in real time according to the changes in parameters such as slant range or angle from the reference target, emerges as an effective strategy for mitigating pulse interference caused by RCM [19].
In recent years, scholars have adopted the variable PRI technique. They proposed numerous methodologies involving continuously varying PRI (CV-PRI) to tackle blind zones resulting from fixed PRI sequences, as detailed in references [20,21,22,23,24,25,26]. However, due to the temporal design complexity, implementing CV-PRI sequences in engineering applications presents notable challenges. Theoretically, an alternative solution is to employ stepwise varying PRI (SV-PRI) sequences, which essentially strikes a balance between design complexity and data transmission rate.
The variation of PRI, equivalent to non-uniform sampling in an azimuth, introduces several challenges to subsequent imaging processing. Therefore, the processing of non-uniformly sampled data before imaging is a crucial aspect to consider. The Back-Projection Algorithm (BPA), as recommended in the reference [22], is well-suited for focusing non-uniformly sampled data in the azimuth-time domain. However, the principle of reconstructing a uniformly sampled signal opens up opportunities to employ both efficient and conventional frequency domain algorithms, as discussed in references [24,25,26]. These processing methods typically encompass two primary aspects: de-aliasing and interpolation techniques. The subsequent exposition will delve into these two aspects in a comprehensive manner.
Since the pulse repetition frequency (PRF) is lower than the Doppler bandwidth, the raw signal data are under-sampled and aliased, which implies that direct resampling will lead to unacceptable errors. Reference [13] advocates for de-ramping in two-step processing to remove aliasing and introduces a nonlinear shift processing method which removes the squint angle’s impact on focusing. Reference [9] proposes the extended two-step approach (TSA) to process the full-aperture and high-resolution data. However, reference [26] points out that even when utilizing the TSA, the spotlight SAR data, characterized by their high resolution, high squint angle, and large transmitted bandwidth, still suffer from severe aliasing. Furthermore, while the reference [26] recommends employing phase compensation based on the reference phase history [27], this method imposes an upper bound constraint on the swath width in the direction of the radar’s line of sight. This issue fundamentally stems from the inherent contradiction between resolution and swath width. Moreover, the significant Doppler effect, a consequence of the large transmission bandwidth, requires extensive zero-padding to maintain the effectiveness of the time-domain segmentation in the range direction. Undoubtedly, this method leads to a significant increase in data volume. Consequently, it becomes essential to explore alternative solutions that address the resolution-swath width contradiction.
In interpolation processing, the predominant methods include the non-uniform fast Fourier transform (NUFFT) and various interpolation algorithms. In reference [26], the analysis demonstrated that the weighted NUFFT is accurate and efficient in comparison with the best linear unbiased (BLU) interpolation and other NUFFT implementations. Furthermore, reference [28] proposes a novel Sinc-interpolation method based on the non-uniform discrete Fourier transform (NUDFT). This method eliminates ambiguous targets and offers focusing performance comparable to that of the NUDFT. Additionally, reference [12] describes performing baseband Lagrange interpolation to resample the data onto a uniform grid.
This study mainly analyzes the reception and processing of spotlight SAR raw data with a decimeter resolution (0.1 m) and large squint angle (40°). In terms of parameter design, the paper presents a specific strategic solution for the SV-PRI sequences with variable granularity, aimed at reducing engineering complexity. The phenomenon of spaceborne echoes spanning multiple pulse repetition intervals often triggers cascading effects in the design of timing sequences. Consequently, the sharp decrease characteristic of SV-PRI may result in the loss of echo data. To solve this issue, this study introduces constraints on abrupt reductions in SV-PRI, and theoretical deduction demonstrates that the proposed sequence effectively circumvents the problem. Regarding data processing, we propose a processing flow based on the azimuth partition to obtain uniformly sampled signal data, thereby alleviating the inherent contradiction between resolution and swath width. This study tackles the compatibility issues associated with azimuth data partitioning and reassembly in azimuth resampling, and introduces a novel partitioning approach based on temporal point remapping (TPR). The TPR strategy facilitates the effective reconstruction of a uniform grid, substantially simplifying the processing workflow.
The rest of this paper is organized as follows. Section 2 starts with a brief review of high-squint spotlight SAR imaging geometry, and then the analysis of the SV-PRI sequence design is described. Section 3 discusses the proposed processing of non-uniformly sampled signals. The experimental results of the proposed method are presented in Section 4. Discussion and remarks on possible further work are given finally in Section 5.

2. SV-PRI Sequences for High-Squint Spotlight SAR

2.1. Acquisition Geometry

As the orbital model is not the focal point of this paper, an equivalent simplified model of the satellite-ground geometry is constructed for the parameter design analysis. Under the assumptions of a linear satellite trajectory and a planar Earth surface, and ignoring Earth’s rotation, a Cartesian coordinate system is established for a high-resolution and high-squint spotlight SAR. This system takes the subsatellite point as the origin, where the y-axis aligns with the satellite’s velocity vector, the x-axis is orthogonal to the y-axis and extends along the ground range, and the z-axis is directed upward and perpendicular to the Earth’s surface. As illustrated in Figure 1, the imaging geometry considers an orbital altitude of H, an off-nadir angle of α 0 , an along-track velocity of V r , and a range swath width of W r . The reference target in the scene is denoted by the point P, and R 0 represents the minimum slant range to the radar platform in a positive side-looking configuration. The spotlight irradiation begins at time η s and ends at time η e , with corresponding squint angles marked as θ max and θ min . The central time‘s squint angle is θ s q , and the beam rotation angle is θ s y n .
During data acquisition in spotlight SAR, the antenna’s direction and velocity are meticulously controlled to facilitate prolonged illumination of the observed area. Concurrently, using high-squint angles broadens the variation range in slant range between a point target and the radar platform throughout the synthetic aperture period, consequently increasing the magnitude of RCM. The simulation analysis, depicted in Figure 2, investigates the maximal migration distance of the point target in spaceborne spotlight SAR when the nearest slant range is 1200 km, in relation to changes in resolution and squint angle. In Figure 2a, the central squint angle of the spaceborne SAR is maintained at 5°, whereas in Figure 2b, the azimuth resolution is fixed at 0.1 m. The analysis reveals that when azimuth resolution is 0.1m, the RCM can surpass 18 km at small squint angles. Notably, when the squint angle is increased to 45°, the RCM can escalate to as much as 300 km. Such a high magnitude of RCM significantly restricts the imaging swath and hinders the system’s ability to achieve the required azimuth resolution. Therefore, to enable high-resolution imaging within the effective swath width, it is critical to consider strategies for reducing RCM.

2.2. Design of the SV-PRI Sequence

PRI variation offers a practical scheme to significantly reduce RCM in echo data. In this section, we will detail the principles of pulse transmission and reception corresponding to fixed PRI and varying PRI, and introduce a novel iterative computation strategy for the SV-PRI sequence.
In single-base spaceborne SAR, pulse transmission and echo reception typically occur within a single PRI. Figure 3 illustrates the pulse timing of three different sequence types. Assuming a constant length and opening time of the receive echo window, denoted by T w and T s , the pulse duration is T p . P R I s in Figure 3a signifies a fixed interval, P R I n in Figure 3b represents the nth PRI value, and P R I k in Figure 3c represents the PRI value for the kth period. Figure 3 demonstrates that data blind zone can be mitigated by aligning echo data from the reference target with the receive echo window, using varying PRI calculated from variable parameters like slant range or squint angle. For the CV-PRI sequence, the interval between the arrival of the reflected pulse and the transmission of the previous pulse is kept constant at Δ η c . The SV-PRI sequence aims to keep the nth interval Δ η n n = 1 , 2 , 3 , , N within an optimal range, and N is the total number of transmitted pulses. These two PRI variations allow the radar system to continuously gather signal data from the target area without gaps. However, considering the design complexity of the CV-PRI sequence, the SV-PRI sequence is more advantageous. For the accuracy comparison in data processing between the two types of sequences, an in-depth analysis will be conducted in the subsequent simulations.
In the following, we will introduce a method for calculating the SV-PRI sequence. The discussion begins by focusing on the receive echo window. As illustrated in Figure 4, T w s represents the interval for safeguarding data, and T e c h o signifies the echo length, which can be estimated by
T e c h o = 2 ( R f a r R n e a r ) c + T p
where R f a r and R n e a r are the slant ranges of the farthest and nearest point targets in the scene at the central moment. To accommodate the specific transceiver mechanism, the receive echo window must meet two primary conditions, which are given by
T s > 0
T w P R I min T p 2 T s
Formula (2) stipulates that the window only activates after the transmitter is fully deactivated. Additionally, it follows from Formula (3) that the duration of reception must not exceed the minimum sampling interval P R I min in the variable PRI sequence.
Take the target at scene center as the reference target. First of all, the initial sampling time interval P R I s should be determined according to the acquisition geometric parameters and the beam Doppler bandwidth. Assume that the pulse transmission time vector is η u = [ η u 1 , η u 2 , η u 3 , , η u N ] , and η u n means the start time of the nth transmitted pulse. The echo delay time of nth transmitted pulse, defined as η d n , can be calculated using slant range models R r e f ( η u n ) with different accuracy. Based on start-stop model, η d n can be approximately equivalent to
η d n 2 R r e f ( η u n ) c
where c is the speed of light. To meet the high-precision timing design requirements for high-resolution imaging, η d n can be updated according to [26]:
η d n = R r e f ( η u n ) + R r e f ( η u n + η d n ) c
The echo delay time from the reference target at the central time is defined as
η d c = 2 R 0 c cos θ s q
Let M denote the number of complete P R I s spanned by the echo delay time, which is given by
M = η d c P R I s
Then, as shown in Figure 3a, Δ η n equals
Δ η n = η d n M · P R I s
With the fixed T w , T s , and T p , as shown in Figure 3b, CV-PRI sequences can reduce migration in a finite receive echo window by controlling Δ η n to approximate Δ η c , which can be fixed to
Δ η c = T w 2 + T p 2 + T s
In the given acquisition geometry, the slant range of the reference target gradually decreases, resulting in a consistent leftward shift of the echo data within the window when employing a fixed PRI sequence. For the SV-PRI sequence, as shown in Figure 4, to maximize the utilization rate of the receive echo window, the first echo data of each period should be positioned on the right side as much as possible, which requires updating Δ η c ; then,
Δ η c = T p + T s + T w 2 T e c h o 2 T w s
P R I s and squint angle θ s q in the azimuth central time determine the start time η u 1 of the first transmitting pulse in the synthetic aperture time. η d 1 can be calculated using Formulas (4) and (5), and then,
P R I 1 = η d 1 Δ η c M
To specify the period of the SV-PRI sequence, we define a granularity factor δ δ > 0 , and the SV-PRI sequence period is denoted by δ M = δ · M δ M 1 , which requires positive integer quantification. As previously described, Δ η n needs to be maintained within a certain range; that is,
T p + T s + T w s < Δ η n < Δ η c
To achieve the goal of complete echo capture within this range, δ M must meet the condition of
η d 1 η d δ M < T w T e c h o 2 · T w s
that is,
2 R r e f ( η u 1 ) R r e f ( η u 1 + ( δ M 1 ) · P R I 1 ) c < T w T e c h o 2 · T w s
The number of periods in the SV-PRI sequence is
K = N δ M
Using the iterative idea, subsequent P R I n and η u n can be calculated until the entire sequence is obtained; that is,
P R I n = ( η d n Δ η c ) ( η d n Δ η c ) M M , n = 1 , δ M + 1 , 2 δ M + 1 , , K · δ M + 1 P R I n 1 , e s l e
η u ( n + 1 ) = η u n + P R I n , n = 1 , 2 , 3 , , N 1

2.3. Analysis of Echo Data Loss

As mentioned in reference [24], due to the echo delay of M PRI before arrival in spaceborne SAR, there is a potential risk of echo data loss between two consecutive periods in the SV-PRI sequence. This paper conducts a thorough theoretical analysis to determine if the proposed SV-PRI sequence experiences this issue.
Assuming full capture of the last M pulse echoes in the kth period, the Formula (12) can be rewritten as
T p + T s + T w s < η d ( ( k + 1 ) · δ M M + l ) ( M l + 1 ) · P R I k ( l 1 ) · P R I k + 1 < Δ η c
where k = 1 , 2 , 3 , , K 1 , and 1 < l M . Define the abrupt decrement value between two consecutive periods as Δ α k = P R I k P R I k + 1 , and the range of Δ α k is given by
M · P R I k + ( T p + T s + T w s ) η d ( ( k + 1 ) · δ M M + l ) l 1 < Δ α k < M · P R I k + Δ η c η d ( ( k + 1 ) · δ M M + l ) l 1
Formula (19) delineates the specific range for abrupt decreases in SV-PRI. According to the imaging geometry model constructed in this study, data loss is related to the constraints on the right side of Formula (19). As a result, this study focuses solely on the right-hand side constraint. From Formula (16), the values of P R I k and P R I k + 1 are
P R I k = η d ( k · δ M + 1 ) Δ η c M
P R I k + 1 = η d ( ( k + 1 ) · δ M + 1 ) Δ η c M
Therefore, Δ α k can be updated to
Δ α k = η d ( k · δ M + 1 ) η d ( ( k + 1 ) · δ M + 1 ) M
And the right side of Formula (19) can be rewritten as
η d ( k · δ M + 1 ) η d ( ( k + 1 ) · δ M + 1 ) M < η d ( k · δ M + 1 ) η d ( ( k + 1 ) · δ M M + l ) l 1
that is,
δ · η d ( k · δ M + 1 ) η d ( ( k + 1 ) · δ M + 1 ) δ · M < δ M M + l 1 l 1 · η d ( k · δ M + 1 ) η d ( ( k + 1 ) · δ M M + l ) δ M M + l 1
In Figure 1, the backward difference in slant range of the reference target presents a decreasing trend, so Formula (24) can be updated to
δ 1 + δ 1 l 1 · M
that is
( 1 M l 1 ) · δ 1 M l 1
Formula (26) applies when δ > 0 and 1 < l M . Therefore, it can be inferred that the solution scheme proposed in this paper meets the necessary and sufficient conditions to prevent pulse echo loss, thereby confirming its effectiveness in preserving signal integrity within the SAR system.

3. Processing of Non-Uniformly Sampled Data

In this section, the processing of non-uniformly sampled data is discussed. As illustrated in Figure 5, the proposed flow primarily encompasses three key aspects: Doppler history recovery, azimuth partitioning and de-aliasing, and azimuth resampling and reassembly. The analysis will be sequentially conducted on these aspects.

3.1. Doppler History Recovery

For any single point target, the baseband signal in the range-frequency domain after non-uniform sampling and demodulation can be expressed as
S ( f τ , η u ) = W r ( f τ ) exp j 4 π ( f 0 + f τ ) R ( η u ) c exp j π f τ 2 K r exp j 2 π f τ Δ τ
where f τ denotes the range frequency, f 0 denotes the carrier frequency, W r ( f τ ) denotes the range spectral envelope, and K r denotes the rate of frequency modulation in the chirp signal. The modification of PRI results in a temporal displacement Δ τ during the initial phase of rapid time sampling. In this process, Δ τ leads to the archiving and documentation of the 2D signal at a specific position, distinct from the range gate determined by the fixed PRI. To recover the Doppler history information, the consistent starting sampling can be compensated by multiplying the SV-PRI signal with the phase function H d ( f τ , η u ) in the range-frequency domain. The expression of H d ( f τ , η u ) is
H d f τ , η u = exp j 2 π f τ Δ τ
However, this phase compensation results in significant signal folding in the range-time domain. In high-resolution and high-squint SAR systems, the conventional remediation, zero-padding, is no longer feasible. Consequently, the signal reconstruction methodology proposed in this paper is entirely based on the range-frequency domain, which effectively avoids the influence of this phenomenon.

3.2. TPR-Based Azimuth Partitioning and De-Aliasing

The partitioning approach presented in this paper primarily aims to surmount the limitations imposed by the range swath, thereby facilitating effective de-aliasing. As noted in the reference [26], the range swath along the radar’s line of sight is limited by
Δ R < α 1 α · P R I max R r e f ( η u c ) B d o p , r e f
where Δ R represents the maximum differential range of slant ranges in the direction of the radar’s line of sight at azimuth central moment η u c , and B d o p , r e f denotes the Doppler bandwidth of the reference target. Considering the impact of the signal folding phenomenon mentioned earlier, dividing the range swath in the range-time domain to meet the condition given by Formula (29) is not an optimal solution. Essentially, Formula (29) exemplifies the contradiction between resolution and swath width. A more effective approach is to shift directly from range partitioning in the range-time domain to azimuth partitioning in the range-frequency domain. The entirety of the data in the azimuth time can be partitioned into multiple single-segment datasets, each acquired over a short azimuth time η u ( k ) (k means the index of the kth dataset). In this way, the resultant residual bandwidth is compensated through an appropriate oversampling factor, and Formula (29) can be rewritten as
α k 1 α k · P R I max k R r e f ( η u c ( k ) ) B d o p , r e f ( k ) > Δ R
The concept of azimuth partitioning in varying-PRI spotlight SAR, as previously introduced in the literature [29], is primarily aimed at conforming to the stepwise characteristics of the PRI sequence. This method undergoes two stages of interpolation. Initially, each segment of uniformly sampled data is zero-padded in the azimuth frequency domain to achieve a similar sampling rate. Then, after reassembling the datasets, a secondary, more precise interpolation is performed to establish a complete dataset with uniform grid points. These two interpolation processes cumulatively contribute to a considerable computational burden. To further simplify this processing, this paper introduces a data partitioning and reassembly strategy based on TPR. The strategy allows each data segment to undergo a single interpolation, resulting in a strictly uniform sampling signal after reassembly. More importantly, this method reduces the number of interpolations required, significantly decreasing the computational load. Next, the details of this method will be elaborated upon.
Firstly, by analyzing the non-uniform transmission pulse vector η u and the number of sampling points N, the corresponding uniform transmission pulse vector η n and uniform interval P R I n e w for the reconstructed uniform grid can be determined; that is,
P R I n e w = η u N η u 1 N 1
η n = η n i η n i = η u 1 + ( i 1 ) · P R I n e w , i = 1 , 2 , , N
Determine the point number N b in each single-segment dataset that satisfies the condition of Formula (30), deriving the quantity of dataset K. The numbers of the kth non-uniformly and uniformly sampled points are denoted by N k and N k . The kth non-uniformly and uniformly sampled data are calibrated as P k and P k , where k = 1 , 2 , , K . Note that point number N K of datum P K is not strictly equal to N b . Then, the single-segment pulse transmission time vector is given by
η u ( k ) = η u i ( k ) i ( k 1 ) N b + 1 , k N b + 1 , k = 1 , , K 1
η u ( k ) = η u i ( k ) i ( k 1 ) N b + 1 , N , k = K
Define a set of temporal points T b , which is given by
T b = t b k t b k = η u ( k 1 ) N b + 1 ( k ) , k = 1 , 2 , , K
Define the remapping coefficient set M b , which is given by
M b = m b k m b k = t b k P R I n e w , k = 1 , 2 , K
The mapped uniform transmission pulse vector is denoted by η n ( k ) ; then,
η n ( k ) = m P R I n e w m m b k , m b k + 1 , k = 1
η n ( k ) = m P R I n e w m m b k 1 + 1 , m b k , k = 2 , , K
Following the computation of Formula (31) to Formula (38), the uniform grid η n ( k ) is remapped from the non-uniform grid η u ( k ) . Then, the de-aliasing operation is performed via a multiplication with the phase function [29]:
Q k ( f τ , η u ( k ) ) = exp j 4 π ( f 0 + f τ ) R r e f ( η u ( k ) ) c
Figure 6 presents a concrete example of the TPR-based partitioning strategy, depicting a situation where the data are divided into three segments. From the illustration, it is evident that the start points of each segment on the non-uniform axis form a set of temporal points T b . Utilizing the TPR approach, each start point on the non-uniform axis is mapped to the end point of the previous segment and the start point of the subsequent segment on the uniform axis. The interval between these two points is maintained at P R I n e w , ensuring a completely uniform sampled grid after data reassembly. The TPR strategy is based on data correlation for remapping. It does not require setting the number of divided datasets to be the same as the number of changes in the SV-PRI sequence. Thus, the processing method is also applicable to SV-PRI with smaller δ and CV-PRI.

3.3. Azimuth Resampling and Reassembly

After the processing of azimuth partitioning and de-aliasing, the sampling intervals of the single-segment data P k all conform to the Nyquist criterion. For the data P k , that is, s ( η u ( k ) ) , an efficient and precise interpolation function can be utilized to reconstruct the values s ( η n ( k ) ) in a uniform grid η n ( k ) . This paper selects a weighted-Sinc interpolation suitable for non-uniformly resampled signals [28], with the function expressed as
s ( η n ( k ) ) = 1 P R I n e w i = 0 N k 1 s ( η u i ( k ) ) Δ η u i ( k ) s i n c 1 P R I n e w ( η n ( k ) η u i ( k ) ) , k = 1 , 2 , , K
where Δ η u i ( k ) is the ith non-uniform sampling interval of datum P k . Then, the resampled datasets should be multiplied with the function
Q k ( f τ , η n ( k ) ) = exp j 4 π ( f 0 + f τ ) R r e f ( η n ( k ) ) c
After phase compensation, each dataset is reassembled based on its indices; that is,
η n = η n ( 1 ) , η n ( 2 ) , η n ( 3 ) , , η n ( K )
s ( η n ) = s ( η n ( 1 ) ) , s ( η n ( 2 ) ) , s ( η n ( 3 ) ) , , s ( η n ( K ) )
where s ( η n ) represents the signals after uniform reconstruction in the range-frequency and azimuth-time domain. By utilizing azimuth fast Fourier transform (FFT), the 2D spectrum that is equivalent to uniformly sampled data can be acquired.

4. Simulation Results

In this section, simulated spaceborne SAR data are used to analyze the performance of the proposed method. The simulation focuses on the reference target and utilizes its performance as a benchmark to evaluate the effectiveness of the proposed SV-PRI sequences and the uniform signal reconstruction process. The data geometry acquisition model in the simulation is shown in Figure 1, and specific system parameters are summarized in Table 1. All the simulation results in this paper are implemented in the MATLAB 2022b software.

4.1. Performance Analysis of the SV-PRI Sequence

4.1.1. Granularity Factor and Residual RCM

Firstly, we will discuss the inhibitory effects of the SV-PRI sequence on the RCM. The simulation provides two example results of SV-PRI sequences, as shown in Figure 7. Figure 7a,b illustrate the SV-PRI sequences with granularity factors of δ = 1 and δ = 10 and their residual RCM, respectively. From the left side of Figure 7 and with its magnified view, it is evident that these sequences exhibit distinct stepwise characteristics. Sequences with larger granularity factor δ possess a larger period of variation while also having lower temporal design complexity. The right side of Figure 7 indicates that both SV-PRI sequences reduce the residual RCM to within a controllable range. When the granularity factor δ = 10, the absolute value of the residual RCM is reduced to approximately 200 m. Meanwhile, when the granularity factor δ = 1, the absolute value of the residual RCM decreases to about 0.00025 m. The numerical results depicted in Figure 7 indicate that the finer the granularity factor of the sequence is, the stronger the inhibitory effect on RCM will be. Compared with the previously mentioned RCM of several hundred kilometers, the SV-PRI sequences proposed in this paper have a markedly significant effect on the suppression of RCM.
In summary, as the granularity factor increases, temporal design complexity is reduced, while the residual RCM escalates. In practical engineering, an appropriate granularity factor can be selected based on specific requirements, seeking a balance between temporal design complexity and residual RCM.

4.1.2. Resampling Relative Error

To evaluate the SV-PRI resampling performance, we define the relative error as the ratio of the error energy to the true signal energy [30], which is given by
ε = p = 1 N s ( p ) s ( p ) 2 p = 1 N s ( p ) 2
In Formula (44), s ( · ) is the signal sample containing the interpolation error, and s ( · ) is the true uniformly sampled signal. For comparison, this simulation selects the CV-PRI sequences in references [24,25], and the relative errors are as shown in Table 2. The results demonstrate that, compared with the CV-PRI, the SV-PRI sequences’ errors are smaller. And it is particularly noteworthy that the SV-PRI with the period M exhibits superior reconstruction accuracy, which is reduced by about 9 dB. Obviously, the proposed SV-PRI sequences show advantageous reconstruction performance with higher precision.

4.2. Analysis of Uniform Signal Reconstruction Processing

In this section, the approaches described in references [28,29] have been simulated and contrasted with the method proposed in this study. To quantitatively assess the uniform reconstruction performance of this paper’s processing flow compared with the other literature, the relative error described in Section 4.1.2 continues to be selected as the evaluation metric.
For generality, we selected two SV-PRI sequences with the highest reconstruction accuracy and the lowest temporal design complexity for simulation and comparative analysis. The resampling errors in the uniform signal reconstruction methods are presented in Table 3. The results indicate that the TSA method, combined with weighted-Sinc interpolation as referenced in [28], leads to severe deterioration errors. This phenomenon can be attributed to the fact that, as demonstrated in Figure 8a, the high-resolution, high-squint spotlight SAR data exhibit severe aliasing in the 2D spectrum even after TSA de-ramping. This issue is fundamentally rooted in the inherent contradiction between resolution and swath width. In Figure 8b, the TPR-based azimuth partitioning method alleviates this inherent contradiction, and as a result, the 2D spectrum of the reference target has been effectively compressed. Figure 8c demonstrates that the 2D spectrum achieves more accurate compression in the azimuth direction. Compared with the reference phase history method combined with the NUFFT in reference [29], the proposed method’s resampling error is similar when the granularity factor δ = 1. However, at a granularity factor of δ = 10, due to the NUFFT amplifying the periodic effects of the SV-PRI sequences, our method shows superior reconstruction accuracy.

4.3. Focusing Results

To analyze the impact of SV-PRI and the proposed method for uniform signal reconstruction on focused imaging, the simulation employs the BPA to focus on the reference point target. This study compares the focusing results of uniform signals with those of non-uniform signals processed using the proposed uniform reconstruction method. Table 4 presents the simulation comparison results for the Peak Sidelobe Ratio (PSLR), Integrated Sidelobe Level Ratio (ISLR), and Impulse Response Width (IRW). The consistency between these results underscores that the imaging quality of the non-uniform signal, when processed with the method proposed in this study, is on par with that of the uniform signal. As illustrated in Figure 9, when processed according to the method outlined in this paper, the imaging of the reference point target achieves favorable focusing effects.

5. Discussion

In summary, the innovative work proposed in this paper is aimed at raw unfocused data for high-resolution images, primarily involving two major aspects: the reception and processing of the raw signal. For the significant challenges of severe RCM in high-resolution and high-squint spaceborne spotlight SAR, this paper proposes a novel strategy involving the variable granularity SV-PRI sequence, which ensures the reception of raw signal data and lays the foundation for high-resolution imaging. Compared with the CV-PRI sequence, the SV-PRI sequence shows higher resampling accuracy. Additionally, its stepwise feature minimizes the timing design complexity, thus simplifying engineering implementation. Varying-PRI technology resolves the issue of RCM, but leads to a non-uniformly sampled signal. To address the problem, this paper proposes a comprehensive processing workflow with TPR-based azimuth partitioning and weight-Sinc interpolation for uniform reconstruction. The results show that the proposed method exhibits superior performance.
This method has achieved commendable results, but it also possesses certain limitations. The corresponding data on the non-uniform and uniform grid are not of exactly the same length, which has a certain scale of temporal scaling. This characteristic imposes higher demands on subsequent resampling techniques.
It is noteworthy that the application of certain special mathematical sequences in SV-PRI sequences may enhance reconstruction accuracy and RCM suppression capability. In the future, we will focus on the potential application of special mathematical sequences in varying PRI techniques, such as coprime sequences, the Chinese Remainder Theorem, prime gap sequences, and so on. On the other hand, research will be primarily focused on developing interpolation and resampling algorithms with higher precision, expanding the application of the proposed method in high-resolution and high-squint spotlight SAR.

6. Conclusions

This paper introduces a novel strategy involving the variable granularity SV-PRI sequence to address the challenges of severe RCM in high-resolution, high-squint spaceborne spotlight SAR. This sequence significantly improves the reception of raw data, thereby establishing a solid foundation for achieving high-resolution imaging. In comparison with the conventional CV-PRI sequences, the SV-PRI methodology enhances resampling accuracy and reduces the complexity of timing design, facilitating simpler engineering implementations. Additionally, a comprehensive processing workflow that includes TPR-based azimuth partitioning and weighted-Sinc interpolation is introduced to resolve non-uniform sampling issues, which achieves notable results. In the future, more efforts will be concentrated on developing more precise interpolation and resampling algorithms to extend the application of this method in high-resolution and high-squint spotlight SAR scenarios.

Author Contributions

Conceptualization, Y.H. and W.W.; methodology, Y.H. and W.W.; software, Y.H.; validation, Y.H., W.W. and X.W.; formal analysis, Y.H., W.W., X.W. and Y.D.; investigation, W.W., X.W., Y.D. and D.X.; resources, W.W. and Y.D.; data curation, X.W. and D.X.; writing—original draft preparation, Y.H.; writing—review and editing, W.W.; visualization, Y.H.; supervision, Y.D.; project administration, W.W. and Y.D.; funding acquisition, W.W., X.W., Y.D. and D.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Fund, grant number 61971401 and by the Youth Innovation Promotion Association, CAS.

Data Availability Statement

The data are not publicly available due to privacy.

Acknowledgments

The authors express their heartfelt gratitude to the editors and reviewers for their diligent efforts in facilitating the publication of this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Moreira, A.; Prats-Iraola, P.; Younis, M.; Krieger, G.; Hajnsek, I.; Papathanassiou, K.P. Papathanassiou A tutorial on synthetic aperture radar. IEEE Geosci. Remote Sens. Mag. 2013, 1, 6–43. [Google Scholar] [CrossRef]
  2. Jakowatz, C.V.; Wahl, D.E.; Eichel, P.H.; Ghiglia, D.C.; Thompson, P.A. Spotlight-Mode Synthetic Aperture Radar: A Signal Processing Approach, 1st ed.; Springer Science & Business Media: New York, NY, USA, 2012. [Google Scholar]
  3. Boerner, E.; Lord, R.; Mittermayer, J.; Bamler, R. Evaluation of TerraSAR-X spotlight processing accuracy based on a new spotlight raw data simulator. In Proceedings of the 2003 IEEE International Geoscience and Remote Sensing Symposium, Toulouse, France, 21–25 July 2003; Volume 2, pp. 1323–1325. [Google Scholar] [CrossRef]
  4. Reigber, A.; Scheiber, R.; Jager, M.; Prats-Iraola, P.; Hajnsek, I.; Jagdhuber, T.; Papathanassiou, K.P.; Nannini, M.; Aguilera, E.; Baumgartner, S.; et al. Very-high-resolution airborne synthetic aperture radar imaging: Signal processing and applications. Proc. IEEE 2012, 101, 759–783. [Google Scholar] [CrossRef]
  5. Prats-Iraola, P.; Scheiber, R.; Rodriguez-Cassola, M.; Mittermayer, J.; Wollstadt, S.; de Zen, F.; Bräutigam, B.; Schwerdt, M.; Reigber, A.; Moreira, A. On the processing of very high resolution spaceborne SAR data. IEEE Trans. Geosci. Remote Sens. 2014, 52, 6003–6016. [Google Scholar] [CrossRef]
  6. de Almeida, F.Q.; Younis, M.; Prats-Iraola, P.; Rodriguez-Cassola, M.; Krieger, G.; Moreira, A. Slow Pulse Repetition Interval Variation for High-Resolution Wide-Swath SAR Imaging. IEEE Trans. Geosci. Remote Sens. 2021, 59, 5665–5686. [Google Scholar] [CrossRef]
  7. Dong, B.; Li, G.; Zhang, Q. High-Resolution and Wide-Swath Imaging of Spaceborne SAR via Random PRF Variation Constrained by the Coverage Diagram. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5241016. [Google Scholar] [CrossRef]
  8. Wang, P.; Liu, W.; Chen, J.; Niu, M.; Yang, W. A High-Order Imaging Algorithm for High-Resolution Spaceborne SAR Based on a Modified Equivalent Squint Range Model. IEEE Trans. Geosci. Remote Sens. 2015, 53, 1225–1235. [Google Scholar] [CrossRef]
  9. Zhu, D.; Xiang, T.; Wei, W.; Ren, Z.; Yang, M.; Zhang, Y.; Zhu, Z. An Extended Two Step Approach to High-Resolution Airborne and Spaceborne SAR Full-Aperture Processing. IEEE Trans. Geosci. Remote Sens. 2021, 59, 8382–8397. [Google Scholar] [CrossRef]
  10. Sun, X.; Yeo, T.S.; Zhang, C.; Lu, Y.; Kooi, P.S. Time-varying step-transform algorithm for high squint SAR imaging. IEEE Trans. Geosci. Remote Sens. 1999, 6, 2668–2677. [Google Scholar] [CrossRef]
  11. Davidson, G.W.; Cumming, I. Signal properties of spaceborne squint-mode SAR. IEEE Trans. Geosci. Remote Sens. 1997, 3, 611–617. [Google Scholar] [CrossRef]
  12. Zeng, H.C.; Chen, J.; Liu, W.; Yang, W. Modified Omega-k algorithm for high-speed platform highly-squint staggered SAR based on azimuth non-uniform interpolation. Sensors 2015, 15, 3750–3765. [Google Scholar] [CrossRef] [PubMed]
  13. An, D.; Huang, X.; Jin, T.; Zhou, Z. Extended two-step focusing approach for squinted spotlight SAR imaging. IEEE Trans. Geosci. Remote Sens. 2011, 50, 2889–2900. [Google Scholar] [CrossRef]
  14. Fan, W.; Zhang, M.; Li, J.; Wei, P. Modified Range-Doppler Algorithm for High Squint SAR Echo Processing. IEEE Geosci. Remote Sens. Lett. 2019, 16, 422–426. [Google Scholar] [CrossRef]
  15. Zeng, L.; Liang, Y.; Xing, M.; Li, Z. A Novel Motion Compensation Approach for Airborne Spotlight SAR of High-Resolution and High-Squint Mode. IEEE Geosci. Remote Sens. Lett. 2016, 13, 429–433. [Google Scholar] [CrossRef]
  16. Xuying, X.; Gen, L.; Yanheng, M.; Lina, C. New slant range model and azimuth perturbation resampling based high-squint maneuvering platform SAR imaging. J. Syst. Eng. Electron. 2021, 32, 545–558. [Google Scholar] [CrossRef]
  17. Ding, Z.; Zheng, P.; Li, H.; Zhang, T.; Li, Z. Spaceborne High-Squint High-Resolution SAR Imaging Based on Two-Dimensional Spatial-Variant Range Cell Migration Correction. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5240114. [Google Scholar] [CrossRef]
  18. Chen, X.; Sun, G.-C.; Xing, M.; Li, B.; Yang, J.; Bao, Z. Ground Cartesian Back-Projection Algorithm for High Squint Diving TOPS SAR Imaging. IEEE Trans. Geosci. Remote Sens. 2021, 59, 5812–5827. [Google Scholar] [CrossRef]
  19. Zhang, Y.; Yu, Z.; Li, C. Effects of PRF variation on spaceborne SAR imaging. In Proceedings of the 2013 IEEE International Geoscience and Remote Sensing Symposium-IGARSS, Melbourne, VC, Australia, 21–26 July 2013; pp. 1336–1339. [Google Scholar] [CrossRef]
  20. Xu, H.; Gao, J.; Li, J. A variable PRF imaging method for high squint diving SAR. Prog. Electromagn. Res. 2013, 135, 215–229. [Google Scholar] [CrossRef]
  21. Yin, W.; Ding, Z.; Yang, S.; Li, Y.; Zeng, T.; Long, T. A continuous PRI variation method for geosynchronous SAR with elliptical orbit. In Proceedings of the 2015 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Milan, Italy, 26–31 July 2015; pp. 4582–4585. [Google Scholar] [CrossRef]
  22. Luo, X.; Xu, W.; Han, L. The Application of Staggered PRF Technology in Slant-Focus SAR. J. Radars 2015, 4, 70–77. [Google Scholar]
  23. Wang, Y.; Li, J.; Yang, J.; Sun, B. A novel spaceborne sliding spotlight range sweep synthetic aperture radar: System and imaging. Remote Sens. 2017, 9, 783. [Google Scholar] [CrossRef]
  24. Wang, P.; Xu, W.; Li, N.; Dong, W. Research on spaceborne high squint spotlight SAR imaging technology with variable PRI. J. Electron. Inf. Technol. 2018, 40, 2470–2477. [Google Scholar]
  25. Zou, H.; Zhao, F.; Jia, X.; Wang, W.; Zhang, H. Continuous pulse repetition interval variation for highly squint high-resolution spotlight SAR imaging. Remote Sens. Lett. 2021, 12, 209–218. [Google Scholar] [CrossRef]
  26. Zhou, Z.X.; Deng, Y.; Wang, W.; Jia, X.; Wang, R. Analysis of Varying-PRI Spotlight SAR Data. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5221020. [Google Scholar] [CrossRef]
  27. Meng, D.; Ding, C.; Hu, D.; Qiu, X.; Huang, L.; Han, B.; Liu, J.; Xu, N. On the processing of very high resolution spaceborne SAR data: A chirp-modulated back projection approach. IEEE Trans. Geosci. Remote Sens. 2017, 56, 191–201. [Google Scholar] [CrossRef]
  28. Chen, S.; Huang, L.; Yu, L. An Improved Sinc Interpolation-Based Staggered PRF Sampling for Beamforming SAR Imaging. J. Radars 2019, 8, 527–536. [Google Scholar]
  29. Zhang, Z.; Xu, W.; Huang, P.; Tan, W.; Qi, Y. Azimuth Preprocessing of Squinted Sliding Spotlight Synthetic Aperture Radar Data with Block Varying PRF. In Proceedings of the 2021 Photonics & Electromagnetics Research Symposium, Hangzhou, China, 21–25 November 2021. [Google Scholar]
  30. Villano, M.; Krieger, G.; Jäger, M.; Moreira, A. Staggered SAR: Performance analysis and experiments with real data. IEEE Trans. Geosci. Remote Sens. 2017, 55, 6617–6638. [Google Scholar] [CrossRef]
Figure 1. Acquisition geometry of high-resolution and high-squint spaceborne spotlight synthetic aperture radar (SAR).
Figure 1. Acquisition geometry of high-resolution and high-squint spaceborne spotlight synthetic aperture radar (SAR).
Remotesensing 16 00871 g001
Figure 2. Variation of range cell migration with resolution (a) and squint angle (b) at a nearest slant range of 1200 km. (a) Central squint angle is 5°. (b) Azimuth resolution is 0.1 m.
Figure 2. Variation of range cell migration with resolution (a) and squint angle (b) at a nearest slant range of 1200 km. (a) Central squint angle is 5°. (b) Azimuth resolution is 0.1 m.
Remotesensing 16 00871 g002
Figure 3. The transmission and reception processes with different sequences. (a) Fixed pulse repetition interval (PRI) sequence. (b) Continuously varying PRI (CV-PRI) sequence. (c) Stepwise varying PRI (SV-PRI) sequence with period M.
Figure 3. The transmission and reception processes with different sequences. (a) Fixed pulse repetition interval (PRI) sequence. (b) Continuously varying PRI (CV-PRI) sequence. (c) Stepwise varying PRI (SV-PRI) sequence with period M.
Remotesensing 16 00871 g003
Figure 4. The receive echo window.
Figure 4. The receive echo window.
Remotesensing 16 00871 g004
Figure 5. Flowchart of the uniform signal reconstruction.
Figure 5. Flowchart of the uniform signal reconstruction.
Remotesensing 16 00871 g005
Figure 6. Schematic diagram of TPR-based azimuth partitioning.
Figure 6. Schematic diagram of TPR-based azimuth partitioning.
Remotesensing 16 00871 g006
Figure 7. SV-PRI sequences and residual range cell migration (RCM) with different granularity factors δ . (a) δ = 1. (b) δ = 10.
Figure 7. SV-PRI sequences and residual range cell migration (RCM) with different granularity factors δ . (a) δ = 1. (b) δ = 10.
Remotesensing 16 00871 g007
Figure 8. Resampling results after de-aliasing. (a) 2D spectrum after using the method in [28]. (b) 2D spectrum after using TPR-based partitioning and the method in [28]. (c) 2D spectrum after using proposed method.
Figure 8. Resampling results after de-aliasing. (a) 2D spectrum after using the method in [28]. (b) 2D spectrum after using TPR-based partitioning and the method in [28]. (c) 2D spectrum after using proposed method.
Remotesensing 16 00871 g008
Figure 9. Focusing performances of the reference target. (a) Focusing result of data resampled with proposed method. (b) Cross range profile.
Figure 9. Focusing performances of the reference target. (a) Focusing result of data resampled with proposed method. (b) Cross range profile.
Remotesensing 16 00871 g009
Table 1. SAR system parameters used in the simulation.
Table 1. SAR system parameters used in the simulation.
ParameterValue
Orbit height1000 km
Antenna length12 m
Carrier frequency9.6 GHz
Pulse bandwidth3 GHz
Off-nadir angle30°
Squint angle40°
Cross-range resolution0.1 m
Initialized PRF1018.1 Hz
Table 2. Resampling errors with different types of PRI.
Table 2. Resampling errors with different types of PRI.
Types of PRIResampling Error (dB)
CV-PRI in  [24]−33.73
CV-PRI in  [25]−32.96
SV-PRI ( δ = 1)−42.77
SV-PRI ( δ = 2)−35.01
SV-PRI ( δ = 5)−34.98
SV-PRI ( δ = 10)−35.03
Table 3. Resampling errors in uniform signal reconstruction methods.
Table 3. Resampling errors in uniform signal reconstruction methods.
Method of Uniform ReconstructionResampling Error (dB)
δ = 1 δ = 10
Method in [9]2.972.98
Method in [29]−42.76−34.33
Proposed method−42.77−35.03
Table 4. Focusing results of the reference target.
Table 4. Focusing results of the reference target.
PSLR (dB)ISLR (dB)IRW
Uniformly sampling−13.2546−10.04660.8696
Proposed method−13.2549−10.04680.8696
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hu, Y.; Wang, W.; Wu, X.; Deng, Y.; Xiao, D. A Novel SV-PRI Strategy and Signal Processing Approach for High-Squint Spotlight SAR. Remote Sens. 2024, 16, 871. https://0-doi-org.brum.beds.ac.uk/10.3390/rs16050871

AMA Style

Hu Y, Wang W, Wu X, Deng Y, Xiao D. A Novel SV-PRI Strategy and Signal Processing Approach for High-Squint Spotlight SAR. Remote Sensing. 2024; 16(5):871. https://0-doi-org.brum.beds.ac.uk/10.3390/rs16050871

Chicago/Turabian Style

Hu, Yuzhi, Wei Wang, Xiayi Wu, Yunkai Deng, and Dengjun Xiao. 2024. "A Novel SV-PRI Strategy and Signal Processing Approach for High-Squint Spotlight SAR" Remote Sensing 16, no. 5: 871. https://0-doi-org.brum.beds.ac.uk/10.3390/rs16050871

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop