A Sparse SAR Imaging Method for Low-Oversampled Staggered Mode via Compound Regularization

: High-resolution wide-swath (HRWS) imaging is the research focus of the modern space-borne synthetic-aperture radar (SAR) imaging field, with significant relevance and vast application potential. Staggered SAR, as an innovative imaging system, mitigates blind areas across the entire swath by periodically altering the radar pulse repetition interval (PRI), thereby extending the swath width to multiples of that achievable by conventional systems. However, the staggered mode introduces inherent challenges, such as nonuniform azimuth sampling and echo data loss, leading to azimuth ambiguities and substantially impacting image quality. This paper proposes a sparse SAR imaging method for the low-oversampled staggered mode via compound regularization. The proposed method not only effectively suppresses azimuth ambiguities arising from nonuniform sampling without necessitating the restoration of missing echo data, but also incorporates total variation (TV) regularization into the sparse reconstruction model. This enhances the accurate reconstruction of distributed targets within the scene. The efficacy of the proposed method is substantiated through simulations and real data experiments from spaceborne missions.


Introduction
The synthetic aperture radar (SAR) is an important remote sensing technology capable of all-weather and all-time imaging, widely utilized in fields including natural environment monitoring, agricultural mapping, and change detection [1][2][3][4][5].High-resolution wide-swath (HRWS) imaging has emerged as a critical developmental direction for spaceborne SAR systems [6].However, traditional SAR systems are hampered by constraints such as the minimum antenna area [7,8] and noise-equivalent sigma zero (NESZ), among other parameters, preventing the achievement of HRWS imaging.To overcome these limitations, researchers globally have embarked on extensive research, leading to the proposition of innovative systems.These include the azimuth multichannel SAR system [9][10][11], multiple-input multiple-output (MIMO) SAR technology [12][13][14], and digital beamforming (DBF) technology [15,16].
In spaceborne SAR systems, to protect the receiving subsystem from being damaged by high-power transmitted pulses, signals cannot be transmitted and received simultaneously, leading to swath-width blind areas.In traditional SAR systems, these blind areas are fixedly located along the azimuth direction due to the constant pulse-repetition frequency (PRF).Staggered SAR, an innovative imaging system, introduces variable PRF to mitigate the issue of consistently positioned blind areas [17][18][19][20][21].By periodically altering the radar pulse repetition interval (PRI), staggered SAR disperses these blind areas across various slant range gates throughout the swath.Consequently, in staggered SAR, blind areas are not merely parallel strips along the track axis, but are evenly distributed throughout the swath in a pattern dictated by the sequence of PRIs [22].This strategic dispersion allows the staggered SAR to significantly widen the swath, achieving HRWS imaging-a capability exemplified by its application in the Tandem-L system [23][24][25].
The unique variable PRF mode inherent in staggered SAR systems introduces two principal challenges: nonuniform sampling and echo data loss, around which existing signal-processing methodologies are centered [26].The German Aerospace Center (DLR) has developed an innovative interpolation method utilizing best linear unbiased (BLU) estimation techniques [22].This method initially employs BLU interpolation to convert nonuniform echo data into a uniformly sampled dataset.Following this resampling, traditional SAR frequency domain imaging algorithms are applied to derive the final imagery.Nevertheless, it is important to note that the fidelity of BLU interpolation diminishes with a reduction in the oversampling rate, resulting in this algorithm only being applicable to highly oversampled staggered SAR.Although high oversampling rates play a critical role in minimizing azimuth ambiguity within staggered SAR frameworks, it is not desirable for spaceborne SAR systems due to expansive data volumes.
Recent studies have increasingly focused on low-oversampled staggered SAR, which presents a promising approach to substantially diminish the volume of echo data requiring storage [27,28].A notable advancement in this domain is the introduction of a twostep algorithm designed for handling low-oversampled staggered SAR data, as detailed in [29].This process begins with the application of the missing-data iterative adaptive algorithm (MIAA) [30], grounded in spectrum estimation, to infer the complete echo spectrum and recuperate absent echo data.Subsequently, a multichannel reconstruction algorithm is employed to reconstruct the echo signal.Despite its innovative approach, this algorithm encounters two significant obstacles: the high computational demand associated with the iterative recovery of missing data and the diminished accuracy of spectral estimation methods in reconstructing distributed targets.In an effort to refine this approach, [31] introduced an enhancement to the aforementioned algorithm, integrating it with the BLU interpolation technique.This enhanced method incorporates a preliminary step to determine if the scene of interest encompasses distributed targets.Depending on the outcome, the process either proceeds with BLU interpolation for distributed scenes or applies the original two-step algorithm otherwise.While the strategies proposed in [29] and [31] offer viable solutions for low-oversampled staggered SAR, they both hinge on the precise restoration of missing echo data to effectively mitigate azimuth ambiguities.
In the past decade, sparse signal processing algorithms based on compressive sensing (CS) theory [32,33] have been successfully introduced to SAR signal processing and markedly enhanced SAR imaging quality [34,35].Studies indicate that, under the condition where the observation matrix adheres to the restricted isometric property (RIP) [36], it is feasible to recover the original sparse signal from a smaller set of samples than what is mandated by the traditional Shannon-Nyquist sampling theorem [37,38].The pioneering work of Çetin et al. in 2001 [39] introduced the concept of applying regularization theory to SAR imaging, aiming at the enhancement of image features.Following this, Patel et al. in 2010 [40] formulated a more comprehensive CS-based SAR model, facilitating scene reconstruction by tackling the L1 regularization problem.This model has since been widely adopted in CS-based SAR imaging [41].Subsequently, Zhang et al. [42] and Çetin et al. [43] furthered this line of inquiry by integrating sparse signal processing with SAR imaging, coining the concept of "sparse SAR imaging".This approach leverages solving the   ( 0 <  ≤ 1) regularization problem to recover scenes of interest.Compared with matched filtering (MF), sparse SAR imaging exhibits superior capabilities in noise and sidelobe suppression, as well as in mitigating azimuth ambiguities attributed to undersampling.The introduction of the approximate observation concept based on azimuthrange decoupling [44] has enabled the application of sparse SAR imaging to large-scale scene reconstructions.This technique has found applications across various constant PRF modes, including strip-map [45], sliding spotlight [46], terrain observation by progressive scans (TOPS) [47], and variable PRF modes [48].Nevertheless, few studies have applied sparse signal processing to the staggered SAR system.In addition, SAR imaging methods that utilize the internal information of images have been focused on research.A novel SAR imaging strategy was proposed in [49].This method segregates the SAR image into sparse and low-rank matrices, reflecting the image's redundancy characteristics, thereby framing the SAR imaging process as a joint sparse and low-rank matrix recovery problem.In [50], a new structural sparse representation-based SAR imaging approach is proposed to effectively depict SAR image structures.This approach establishes an adaptive sparse space to accurately represent the varying local structures of images.
The L1 regularization problem, when translated into an equivalent convex quadratic optimization issue, can be resolved with notable efficiency.However, L1 regularization often introduces additional bias into estimations and affects the reconstruction accuracy [51].Recent studies in CS have highlighted the unique advantages of L1/2 regularization [52].As a nonconvex penalty, L1/2 regularization is lauded for its unbiasedness, capacity to enforce sparsity, and oracle properties, delivering solutions that are notably sparser than those yielded by L1 regularization [53].Moreover, region-based features play a crucial role in applications such as target classification and image segmentation.To this end, the total variation (TV) norm of image magnitude has been integrated as a constraint within the SAR imaging model [54], facilitating the maintenance of continuity in the backscattering coefficient across distributed targets within specified areas [55,56].Obviously, combining L1/2 regularization with TV regularization can significantly enhance the quality of reconstructed SAR images.In [57], the authors proposed an SAR imaging method based on   and TV composite norm regularization.However, this method will inevitably result in huge computational and memory costs.In recent years, deep learning methodologies, particularly those involving deep convolutional neural networks (CNNs), have demonstrated exceptional prowess across various domains, including image restoration [58] and speech signal processing [59].The fusion of deep learning techniques with sparse SAR imaging will emerge as a promising avenue for future research.
In this paper, we propose a sparse SAR imaging method for low-oversampled staggered mode via compound regularization.This method integrates L1/2-regularizationbased sparse signal processing to mitigate azimuth ambiguities, while the incorporation of the TV regularization term boosts the reconstruction accuracy of distributed targets.Our proposed method uniquely accounts for the positioning of blind areas within the sparse SAR imaging model, constructing a blind-area index matrix that signifies echo data loss.This allows for the suppression of azimuth ambiguities without necessitating the recovery of missing data.Furthermore, acknowledging the nonuniform sampling inherent in the staggered mode, our approach employs nonuniform Fourier-transform techniques in the formation of imaging and echo simulation operators, diverging from the traditional azimuth-range decouple operators used in conventional SAR systems.Additionally, we incorporate the TV regularization term into our sparse reconstruction model, facilitating the precise reconstruction of distributed targets.The effectiveness of our method is demonstrated through both simulated data and real spaceborne SAR data experiments.
The remainder of this paper is structured as follows.Section 2 delves into the sparse imaging and reconstruction models specific to staggered SAR.Section 3 outlines the reconstruction process employing L1/2&TV regularization, detailing the construction of imaging and echo simulation operators for staggered SAR and introducing our proposed method.Section 4 is dedicated to numerical simulations and experiments with real data.Section 5 thoroughly analyzes the experimental outcomes.Finally, Section 6 concludes the paper with a succinct summary.

Observation Matrix-Based Sparse Imaging Model of Staggered SAR
In strip-map SAR imaging, we assume that the observed scene is rectangular.Let  denote the backscattering coefficient matrix, whose (, ) entry is the backscattering coefficient, (, ), of the observed scene.The fully sampled echo data of all targets in the observed scene can be expressed as Let   ∈ ℂ ×  denote the 2-D fully sampled echo data matrix, where  is the total number of pulses transmitted in the azimuth direction and   is the number of samples in the range direction.Let  ∈ ℂ   ×  denote the 2-D backscattering coefficient matrix; then, Equation (2) can be rewritten as where   = (  ) and  = () are the vectorized versions of   and  , respectively, the vectorization operation (⋅) reshapes the matrix to vector in column order;  is the additive noise vector; and  is the observation matrix of staggered SAR, defined as It should be noted that, since the PRF of the staggered SAR is no longer constant, the azimuth sampling time   in Equation ( 5) takes values according to the periodically varying linear PRI sequence, which can be expressed as ,  = 2,3, . . ., , where  1 is the azimuth sampling time when the first pulse is transmitted;  is the number of PRIs within one linear variation period;   is the corresponding PRI sequence in this period, with  = 0,1, . . .,  − 1; (⋅) denotes the floor function; and (⋅) denotes the remainder function.
Owing to the particular variable PRF mode, the staggered SAR has the problem of echo data loss, so it is necessary to construct the sparse SAR imaging model in the case of under-sampling.Constructing the sparse imaging model of staggered SAR needs to generate the blind area index matrix indicating the echo data loss, whose constituent elements can be represented as where   is the pulse width.
Let  ∈ ℂ ×  denote the 2-D blind area index matrix; then, the observation matrixbased sparse imaging model of staggered SAR can be expressed as where  ∈ ℂ ×  is the 2-D under-sampled echo data matrix; ⊙ is the Hadamard product operator;  is the radar system observation matrix; and  is the 2-D noise matrix.

Sparse Reconstruction Model of Staggered SAR
In recent years, TV regularization has been widely applied in image processing.It can not only effectively suppress noise, but also restore the edge characteristics of distributed images well.Therefore, in order to maintain the continuity of the backscattering coefficient of distributed targets in a specific observation scene, it can be realized by introducing the TV norm [56].Let the SAR image  be a two-dimensional matrix of size  × .For a certain two-dimensional observation scene, the discrete TV norm of image  can be defined as where ∇ represents the gradient of the image: In Equation ( 10), (∇) , 1 and (∇) , 2 can be specifically represented as where  , represents the pixel value of image  at the -ℎ row and -ℎ column, with ,  = 1,2, … , .
For the observation-matrix-based sparse imaging model represented by Equation ( 8), by combining the   regularization term and TV regularization term into a compound regularization term, the sparse reconstruction model of staggered SAR can be expressed as where  1 is the regularization parameter of the   -norm term and  2 is the regularization parameter of the TV-norm term.In this paper, we take the value of  to be 1/2.In-troducing the TV regularization term into the reconstruction model can enhance the region-based features of the SAR image and reconstruct the distributed targets in the scene more accurately.

Reconstruction Process of L1/2&TV Regularization
It can be seen from Equation ( 13) that the reconstruction model contains two regularization terms, which are combined into a compound regularization term, so we need to find a solution to the linear inverse problem that includes the compound regularization term.For the 1-D linear inverse problem with compound regularizations: the unconstrained linear inverse problem can be transformed into the following equivalent constrained optimization problem according to the variable separation method [60]  ̂= where  1 and  2 are two auxiliary variables.According to the Lagrange multiplier method theorem, the optimization problem represented by Equation ( 15) can be solved by minimizing the following equation: where  1 and  2 are Lagrange multipliers.
In order to minimize (,  1 ,  2 ,  1 ,  2 ) in Equation ( 16), the alternating minimization method is used for the three variables ,  1 , and  2 , respectively.The specific iterative process is shown below.
For the minimization problem represented by Equation ( 17), as the objective function is a quadratic function, a linear equation with a solution can be generated: where  represents the identity matrix.
For the optimization problem represented by Equation (18), it is actually the L1/2 regularization problem [53].The corresponding complex soft thresholding function can be represented as where (  (+1) ) can be denoted as The specific solution steps for Equation (18) will be described in Section 3.3.
For the TV regularization problem represented by Equation ( 19), the solution can be obtained by solving the dual problem of TV regularization by utilizing the Chambolle algorithm [55].After the parameter conversion operation, Equation ( 19) is equivalent to the following equation For the convenience of expression, the TV-norm term regularization parameter  2 /2 2 of the dual problem is denoted as   .
For the optimization problem represented by Equation ( 23), its optimal solution can be expressed as where (⋅) denotes the sign function and  = ( 1 ,  2 ) is the dual variable of  2 , which can be solved by the gradient projection descent method: where ∇(⋅) denotes the gradient operator; (⋅) = −∇ * (⋅) denotes the divergence operator; and  denotes the iteration step size.

Imaging Operator and Echo Simulation Operator of Staggered SAR
In Section 2, a sparse imaging model of the staggered SAR is established, as shown in Equation ( 8).However, it should be noted that, due to the presence of azimuth-range coupling in the 2-D echo data domain, we cannot directly construct the observation matrix  .Therefore, the recovery of the observation scene cannot be achieved by solving the sparse reconstruction model in Equation (13).If we choose to vectorize the echo data matrix to achieve azimuth-range decoupling, perform matrix-vector multiplication, and transform the vector into the form of a matrix after 1-D sparse reconstruction, this will result in huge computational and memory costs.Such a large amount of computation is unacceptable for the large-scale scene.Inspired by [42], according to the idea of approximated observation, an azimuth-range decouple-based sparse SAR imaging method has been proposed [44].The coupling of the 2-D echo data can be removed by constructing the echo simulation operator to replace the observation matrix, which can reduce the computational complexity significantly [45,61].This section will construct the imaging operator and echo simulation operator of staggered SAR.
The chirp scaling algorithm (CSA) is a widely used SAR imaging algorithm [62].Since CSA completes range cell migration correction (RCMC), azimuth compression, and other operations through phase multiplication in the frequency domain, thus avoiding interpolation operations, it has the characteristics of simplicity and high efficiency.In this paper, the imaging operator is constructed based on CSA.
CSA can be summarized into the following three steps: (1) chirp scaling operation; (2) bulk RCMC, range compression, and secondary range compression (SRC) operation; and (3) azimuth focusing and residual phase compensation operation.The process of CSA can be expressed as: where  ̃ is the reconstructed SAR image matrix; ℛ(⋅) is the imaging operator of the staggered SAR system;   is the fully sampled echo data matrix of the staggered mode, as mentioned above; FFT  is the fast Fourier transform (FFT) operator in the range direction; IFFT  and IFFT  are the inverse fast Fourier transform (IFFT) operators in the range and azimuth directions, respectively;   is the nonuniform discrete Fourier transform (NUDFT) matrix, which will be defined in the following text; and  1 ∈ ℂ ×  ,  2 ∈ ℂ ×  , and  3 ∈ ℂ ×  are the chirp scaling operation matrix, the bulk RCMC, range compression and SRC operation matrix, and the azimuth focusing and residual phase compensation operation matrix, respectively.It should be noted that, due to the nonuniform sampling in the azimuth direction caused by the variable PRF mode, the conventional FFT operators can no longer be used in the processing of the azimuth direction in a staggered SAR system.For contrast, we used the NUDFT matrix represented by   in Equation ( 26) to transform the echo data from the time domain to the range-Doppler domain so that the imaging operator ℛ(⋅) of the staggered SAR system is different from that of the constant PRF-mode SAR system.Specifically,   = [ 1 ,  2 , ⋯ ,   ] consists of  column vectors, whose constituent elements are defined as follows: where   is the -ℎ sampling time in the azimuth direction of the staggered SAR system, given by Equation ( 6);   () = ( − 1)/(    ) represents the -ℎ value of the azimuth Doppler frequency, with  = 1,2, … ,   ;   is the number of uniform frequency grids in the azimuth Doppler domain; and   is the average PRI value of the varying PRI sequence.By performing the imaging operator, i.e., the MF operation, on the fully sampled echo data, we can obtain the equivalence relation shown below: However, owing to the presence of noise and sidelobes, the value of  ̃ is always approximately equal to ; that is,  ̃≈ .
For the equivalence relation expressed in Equation ( 28), if ℛ ≈ , then we can use ℛ −1 , the inverse process of the imaging operator of staggered SAR, to approximately replace the observation matrix .It can be seen from Equation ( 26) that the imaging operator ℛ(⋅) is reversible.Therefore, the approximate value of the fully sampled echo data   can be obtained by performing the inverse operation of the imaging operator on the real backscattering coefficient  of the observation scene.The process of inverse CSA can be expressed as where  ̃ is the approximate echo data matrix; ℳ(⋅) is the echo simulation operator of the staggered SAR system; FFT  is the FFT operator in the azimuth direction;   is the nonuniform inverse discrete Fourier transform (NUIDFT) matrix, which will be defined below; and  1 * ,  2 * , and  3 * are the conjugate transpositions of the above three operation matrices in CSA, respectively.Similar to the imaging operator, during the construction of the echo simulation operator, we use the NUIDFT matrix represented by   in Equation ( 29) to transform the data from the range-Doppler domain to the time domain.Therefore, the echo simulation operator ℳ(⋅) of the staggered SAR system is also different from that of a constant PRF mode SAR system, such as the random down-sampling SAR system.  = [ 1 ,  2 , ⋯ ,    ] consists of   column vectors, whose constituent elements are defined as follows: In the staggered SAR system, the imaging operator and echo simulation operator are inverse processes of each other, and their relationship can be expressed as As mentioned above, by approximately replacing the observation matrix  with the echo simulation operator ℳ(⋅), the observation matrix-based model denoted in Equation ( 8) can be rewritten as the azimuth-range decouple-based sparse imaging model of staggered SAR, represented as  ≈  ⊙ ℳ() + . (32)

Sparse Reconstruction Method of Staggered SAR
The iterative process in Section 3.1 (Equations ( 15)-( 25)) is the step of solving Equation (14).As introduced in the previous section, this iterative process needs to store the observation matrix  and perform matrix-vector multiplication operations, which will bring huge memory requirements and computational costs.Therefore, in this section, we need to introduce the idea of approximated observation and use the azimuth-range decouple operators to approximately replace the multiplication of the observation matrixbackscattering coefficient vector.Then, the azimuth-range decouple-based sparse reconstruction model of staggered SAR can be expressed as This section combines L1/2 regularization, TV regularization, and azimuth-range decouple operators effectively, and proposes a sparse SAR imaging method for low-oversampled staggered mode via compound regularization.This method can improve the reconstruction accuracy of distributed targets and maintain the continuity of the backscattering coefficients of distributed targets in a certain area.Table 1 summarizes the iterative procedures of the proposed method.In Table 1, σ represents the noise variance,  () denotes the recovered image, and the expression of the divergence operator is There are several algorithms to achieve sparse reconstruction, such as iterative soft thresholding (IST) [63][64][65] and complex approximated message passing (CAMP) [66,67].In this paper, we chose CAMP as the sparse reconstruction algorithm to solve the L1/2 regularization problem, as shown in Table 1 (step 2-step 6), where   is the real part of the complex soft thresholding function ;   is the imaginary part of ; and   /  and   /  are the partial derivatives of   and   with respect to the real and imaginary parts of the input elements, respectively.

Numerical Simulations
In this part, the numerical simulation experiments will be exploited to compare the reconstruction performance of the matched filtering (MF) method, the BLU interpolation method, and the L1/2&TV-regularization-based method.The numerical simulation experiments were conducted under low oversampling conditions, and the major parameters are shown in Table 2.  PRI 21 In the simulation experiments, the fast-changing PRI sequence design scheme is adopted; that is, in the raw azimuth signal, two consecutive samples are never missed for all slant ranges of interest [22].The PRI sequence values corresponding to the simulation parameters in Table 2 and the blind areas' position of the sequence are shown in Figure 1, respectively.In order to verify the effectiveness of the proposed sparse SAR imaging method for the low-oversampled staggered mode in suppressing the azimuth ambiguities and sidelobes, we will conduct point target simulation experiments and present the quantitative results in this section.In the first simulation experiment, a one-dimensional scene is constructed along the azimuth direction, and one point target is set at the center of the scene.To compare the imaging performance at different slant ranges, that is, different loss rates of missing data, according to the diagram of the blind areas' position shown in Figure 1b, we selected three different slant ranges and set the range position of the point target to 956 km, 982 km, and 994 km, which means that the point target is located inside the blind areas, at the boundary of the blind areas and outside the blind areas, respectively.It should be explained that the inside, boundary, and outside of the blind areas respectively refer to the loss of two samples, the loss of one sample, and the loss of no samples within one PRI variation period.Figure 2 shows the imaging results of low oversampling staggered SAR with different methods, where the MF method refers to the direct imaging processing of nonuniformly sampled echo data using CSA.It can be seen from Figure 2a that the imaging result of the MF method includes multiple pairs of uniformly distributed weak azimuth ambiguities caused by nonuniform sampling and a pair of strong azimuth ambiguities caused by a nonideal azimuth antenna pattern (AAP), which are represented by the green rectangular dotted lines and the brown circular dotted lines, respectively.In order to better evaluate the azimuth ambiguity suppression ability of different methods, we select the integrated sidelobe ratio (ISLR) and the azimuth ambiguity-to-signal ratio (AASR) [68] as quantitative evaluation indicators in this paper, where ISLR only considers the ambiguity components caused by nonuniform sampling and echo data loss.The AASR value can be defined as where  is the main imaging area,   is the number of the pixel in ,  is the ambiguous area, and   is the number of the pixel in .
Table 3 lists the quantitative calculation results of three different imaging methods when the point target is located inside the blind areas (slant range = 956 km).In order to better illustrate the effectiveness of the proposed method, in the second simulation experiment, the ISLR values of different imaging methods were calculated at different oversampling rates.Figure 3 demonstrates the oversampling rate-ISLR relationship curves when the point target is separately located inside the blind areas, at the boundary of the blind areas, and outside the blind areas.Similar to the calculation approach in Table 3, the ambiguity components caused by nonideal AAP are not considered when calculating ISLR here.Next, in order to verify the effectiveness of the proposed sparse SAR imaging method in improving the reconstruction accuracy of distributed targets, distributed targets simulation experiments will be carried out and the corresponding quantitative calculation results will be provided.The distributed targets simulation experiments in this section were also conducted under low oversampling conditions, with the simulation parameters shown in Table 2. Similarly, in the third simulation experiment, a one-dimensional scene was constructed along the azimuth direction, and one distributed target was set at the center of the scene.The position of the distributed target in the range direction was set at the slant range of 956 km, which ensured that this distributed target was located inside the blind areas.Figure 4 shows the imaging results of the one-dimensional distributed target with different methods.
Figure 4a shows the imaging result of the ideal case, where a constant PRF is used to simulate the signal and the echo data are not lost.This constant PRF is selected as the mean PRF value of the fast-changing PRI sequence.It can be seen that the simulation imagery under uniform sampling conditions only contains azimuth ambiguities caused by nonideal AAP, located on both sides of the distributed target.Figure 4b-d are the imaging results of MF, BLU interpolation, and L1/2&TV regularization under the nonuniform sampling condition, respectively.In order to better evaluate the reconstruction accuracy of different methods for distributed targets, the normalized root-mean-square error (NRMSE) is chosen as the quantitative evaluation indicator, which can be defined as where  is the true value of the reference image and  ̂ is the actual value of the reconstructed image.Table 4 lists the quantitative calculation results of the three different imaging methods in Figure 4.In order to further demonstrate the effectiveness of the proposed method, in the fourth simulation experiment, white additive noise with different SNRs is added to the simulated data and the distributed target is placed inside the blind areas and at the boundary of the blind areas, respectively.We compared the reconstruction performance of the L1/2&TV-regularization-based method with BLU interpolation method under low oversampling conditions.For different SNRs, 200 Monte Carlo experiments were conducted, and the results are shown in Figure 5.In the fifth simulation experiment, the SNR was set to 10 dB.The NRMSE values of the BLU interpolation method and L1/2&TV-regularization-based method under different oversampling rates were calculated.To make the conclusions of the simulation experiment more reliable, 200 Monte Carlo experiments were equally used during this simulation process, and the results are shown in Figure 6. Figure 6 demonstrates the oversampling rate-NRMSE relationship curves when the distributed target is separately located inside the blind areas and at the boundary of the blind areas.

Real Data Experiments
In this part, the real data experiments will be carried out to verify the effectiveness of the proposed L1/2&TV-regularization-based method.Gaofen-3 single channel strip-map mode data are chosen as the real spaceborne SAR data.The major system parameters are shown in Table 5.As the original echo data were acquired under the condition of constant PRF, the staggered SAR echo data need to be generated from the uniformly sampled raw data.In order to obtain low-oversampling variable PRF-mode echo data, resampling operations and data loss processing operations in the blind areas were performed on the uniformly sampled raw data according to the parameters in Table 2.In order to explore the improvement effect of the proposed method on the reconstruction performance of distributed targets, two actual scenarios containing sea-land boundary targets were selected as the experimental research objects.Figure 7 shows the imaging results of different methods in the first sea-land boundary scenario.Figure 7a is the imaging result of the original echo data under the constant PRF condition, with only azimuth ambiguities caused by the nonideal AAP present in the imagery.Figure 7b is the imaging result of the direct MF operation on the missing echo data of the variable PRF mode.It can be seen that there were azimuth ambiguities caused by nonuniform sampling and echo data loss in the image.Figure 7c,d are the imaging results of BLU interpolation and L1/2&TV regularization under nonuniform sampling conditions, respectively.
Figure 8 compares the imaging results of the BLU interpolation method and the L1/2&TV regularization-based method, providing the one-dimensional profiles of the strong point target indicated by the blue arrows in Figure 7c,d along the azimuth direction.
In order to verify the effectiveness of the proposed method in improving the accuracy of distributed target reconstruction, three different regions in the first scenario were selected as the research objects, represented by red rectangular frames and named Region 1, Region 2, and Region 3, as shown in Figure 7.The islands where Region 1 and Region 2 are located are represented by green rectangular dashed frames, as shown in Figure 7b,d.The island areas contained in the dashed frames are enlarged and displayed, and the results are shown in Figure 9. Figure 10 shows the imaging results of different methods in the second sea-land boundary scenario.Similarly, two different regions in the second scenario are selected as the research objects, represented by red rectangular frames and named Region 4 and Region 5.
Next, quantitative calculations were conducted on the continuity of the reconstruction results for Region 1, Region 2, Region 3, Region 4, and Region 5.This paper selects the equivalent number of looks (ENL) as the indicator to measure the coherent speckle suppression effect of SAR images.ENL is an indicator that measures the continuity of a uniform region, and the larger the value of ENL, the better the smoothing effect of the image.The definition of ENL is as follows: where   represents the average value of the SAR image intensity and   2 represents the variance of the SAR image intensity.The ENL values of Region 1, Region 2, Region 3, Region 4, and Region 5 were calculated, and the results are shown in Table 6.

Discussion
In this section, we will discuss the experimental results of the previous section in detail.Firstly, the point target simulation experimental results will be discussed.As shown in Figure 2a, the 1-D imaging result of the MF method for the point target located inside the blind areas contains multiple pairs of weak azimuth ambiguities and one pair of strong azimuth ambiguities.Figure 2b presents the imaging result of the BLU interpolation method under the same conditions.Comparing Figure 2a with Figure 2b, it can be found that, under low-oversampling conditions, the results of the BLU interpolation method and the MF method were almost the same, and the azimuth ambiguities could not be suppressed.It can be seen from Figure 3 that, when the azimuth oversampling rate increases, the ISLR values of the three imaging methods all decreased.Moreover, when the oversampling rate is greater than 1.5, that is, in the case of high oversampling, the BLU interpolation method performs better in suppressing azimuth ambiguities caused by nonuniform sampling, but it is not suitable for lower oversampling rates (from 1.1 to 1.5), as shown in Figure 3b.It should be explained that the recovery accuracy of the BLU interpolation algorithm proposed in [22] depends largely on the PRF of the radar system [26].The accuracy of the interpolation decreases when the PRF is close to the Nyquist frequency and improves as the oversampling rate increases [27,29].Therefore, the imaging results in Figure 2 and the quantitative calculation results in Table 3 and Figure 3 are consistent with the theoretical analysis.For the point target, the L1/2&TV regularization actually degenerated into the L1/2 regularization.Figure 2c presents the imaging result of the L1/2-regularization-based method when the point target is located inside the blind areas.Combined with the quantitative results in Table 3, it indicates that the proposed method can effectively suppress the azimuth ambiguities caused by nonuniform sampling and echo data loss in the low-oversampled staggered mode, reducing the ISLR value by about 10 dB.In addition, the AASR values in Table 3 also show that the proposed method had a certain suppression effect on the strong azimuth ambiguities caused by the nonideal AAP.By comparing the results of Figure 2 vertically and combining them with Figure 3, we can find that the values of ISLR were different when the point target was located inside the blind areas, at the boundary of the blind areas, and outside the blind areas, respectively.This is because they lose different samples within one PRI variation period.It is obvious that the ambiguity of the point target was the lowest when it was located outside the blind areas.Comparing Figure 3c with Figure 3a, it can be found that the oversampling rate-ISLR relationship curve of the L1/2-regularization-based method was much lower than that of the MF method under various conditions.Comparing Figure 3c with Figure 3b, although the results of the L1/2-regularization-based method and the BLU interpolation method were similar under high-oversampling conditions, the former performed significantly better than the latter under low-oversampling conditions.Therefore, the proposed sparse SAR imaging method can effectively solve the problem of azimuth ambiguities caused by nonuniform sampling and echo data loss in the low-oversampled staggered mode.
Secondly, the results of the distributed-targets simulation experiments will be discussed.Figure 4b,c show the imaging results of the MF method and BLU interpolation method for the low-oversampled staggered mode, respectively.It can be seen from Figure 4b that, compared with the point target, for the distributed target, the azimuth ambiguities caused by nonuniform sampling were more dispersed.Although the BLU interpolation method could maintain the scattering characteristics of the distributed target, it could not effectively suppress the azimuth ambiguities under low-oversampling conditions, resulting in almost the same imaging results as the MF method, as shown in Figure 4c. Figure 4d shows the imaging result of the L1/2&TV regularization-based method proposed in this paper.It can be seen that the proposed method can effectively suppress the azimuth ambiguities caused by nonuniform sampling without recovering the missing echo data and reduce the reconstruction error of the distributed target.Furthermore, the quantitative calculation results in Table 4 indicate that, under the low-oversampling condition (oversampling rate 1.1), the proposed method can significantly reduce the NRMSE value, while the NRMSE value of the BLU interpolation method is not significantly different from that of the MF method.Combined with the imaging results in Figure 4, it can be concluded that the proposed method can reconstruct distributed targets more accurately.This can be explained by the TV norm characteristic [55,56], which can maintain the continuity of the backscattering coefficient of distributed targets.From Figure 5, it can be seen that, in the presence of noise, the proposed L1/2&TV-regularization-based method can obtain lower NRMSE values than the BLU interpolation method, verifying the accuracy of this method for distributed target reconstruction under low-oversampling conditions.In addition, by comparing Figure 5a,b, it can be seen that, when the distributed target is located at the boundary of the blind areas, the NRMSE values of both methods are lower than those when the target is located inside the blind areas.This indicates that the reconstruction accuracy is related to the number of samples lost within one PRI variation period, and the higher the data loss rate, the lower the reconstruction accuracy.From Figure 6, it can be seen that both the BLU interpolation method and the L1/2&TV-regularization-based method had reduced NRMSE values as the oversampling rate increased, regardless of whether the distributed target was located inside or at the boundary of the blind areas.This illustrates that the reconstruction accuracy was improved, and the larger the oversampling rate, the better the effect.When the oversampling rate was greater than 1.5, that is, in the case of high oversampling, the NRMSE values of the BLU method significantly decreased, being even lower than those of the L1/2&TV-regularization-based method.However, when the oversampling rate was less than 1.2, the reconstruction accuracy of the L1/2&TV-regularization-based method was superior to that of the BLU interpolation method, further verifying that the proposed method could improve the reconstruction accuracy of distributed targets under low-oversampling conditions.Figure 6 also illustrates the impact of the data loss rate on the accuracy of distributed target reconstruction, which was more evident in low-oversampling cases.The above discussion of the distributed target simulation experimental results prove that the proposed sparse SAR imaging method can improve the reconstruction accuracy of distributed targets and enhance the region-based features of SAR images in the low-oversampled staggered mode.
Thirdly, the results of real data experiments will be discussed.The effectiveness of the proposed method in suppressing azimuth ambiguities is discussed first.Figure 7c presents the imaging result of the BLU interpolation method for the first sea-land boundary scenario.It can be seen that azimuth ambiguities caused by nonuniform sampling still existed in the image, proving its limitations under low-oversampling conditions.Figure 7d is the imaging result of the L1/2&TV regularization-based method.It can be seen that, due to the azimuth-range decouple operators, the azimuth ambiguities and false targets in the image were effectively suppressed.Figure 8 further illustrates the conclusion.From the comparison results in Figure 8, it can be seen that the BLU interpolation method could not effectively suppress azimuth ambiguities, so it was not suitable for the low-oversampled staggered mode.Compared with the BLU interpolation method, the L1/2&TV regularization-based method could significantly improve the suppression effect of azimuth ambiguities caused by nonuniform sampling and echo data loss.Next, the effectiveness of the proposed method in improving the accuracy of the distributed target reconstruction is discussed.It can be seen from Figure 9 that there were lots of coherent speckles in the image reconstructed by MF, and the continuity of the island areas could not be satisfied, as shown in Figure 9a.The reconstructed image result using the proposed method was more uniform and continuous, and the coherent speckles and noise in the image were effectively suppressed, as shown in Figure 9b.This proves that the proposed method could maintain the continuity of the backscattering coefficient of distributed targets and improve the reconstruction accuracy of distributed targets for the low-oversampled staggered mode.The numerical calculation results in Table 6 are discussed at the end.As multi-look processing is not performed on real data and the BLU interpolation method actually utilizes matched filtering for imaging processing, the ENL values of the reconstruction result are approximately equal to 1, as shown in the first row of Table 6.In contrast, the L1/2&TV-regularization-based method can effectively improve the ENL values of the reconstruction result, as shown in the second row of Table 6, further demonstrating that the proposed method can enhance the region-based features of SAR images.
Finally, the comparison of simulation and real data experimental results will be discussed.For the point target simulation, we designed two experiments, as stated above.Both point target simulation experiments verified the effectiveness of the proposed method for suppressing azimuth ambiguities in low-oversampled staggered SAR.However, it should be noted that the point target simulation experiments only considered the L1/2 regularization term in compound regularization.Similarly, the distributed targets simulation verified the effectiveness of the proposed method for improving the reconstruction accuracy of distributed targets in low-oversampled staggered SAR.But these experiments mainly considered the TV regularization term in compound regularization.Therefore, it is necessary to select appropriate scenarios in real data experiments and verify the effects of L1/2 regularization TV regularization simultaneously, as shown in Figures 7 and 10.Since the selected scenarios included both strong point targets and distributed targets, such as island areas, compound regularization can be effectively verified.Comparing Figure 8 with Figure 2, it can be concluded that L1/2&TV regularization can indeed suppress the azimuth ambiguities caused by nonuniform sampling and echo data loss in the lowoversampled staggered mode.In addition, the ISLR indicators of the azimuth profiles of different imaging methods shown in Figure 8 were quantitatively calculated.The result was that the ISLR value of L1/2&TV regularization was 12 dB lower than that of BLU interpolation, which was consistent with the calculation results in Table 3. Comparing Table 6 with Table 4, it can be found that two different evaluation indicators (NRMSE and ENL) both demonstrated that L1/2&TV regularization had more significant advantages in distributed target reconstruction.In summary, the comparison between simulation and real data experiments fully demonstrates the effectiveness of compound regularization in both azimuth ambiguity suppression and accurate reconstruction of distributed targets.
It should be pointed out that, if we want to better suppress the azimuth ambiguities caused by nonideal AAP, the group sparsity property [69] should be considered when constructing the imaging model, and the nonuniform sampling problem should be solved during the group sparse reconstruction process [70].This is also work we will carry out in the future.

Conclusions
This paper proposes a sparse SAR imaging method for the low-oversampled staggered mode via compound regularization.By incorporating the blind area index matrix to represent echo data loss within the sparse imaging model and integrating the L1/2 regularization term into the sparse reconstruction model, our method effectively mitigates azimuth ambiguities in staggered SAR systems without the need to recover missing data.Furthermore, by applying both the imaging and echo simulation operators specific to staggered SAR systems for azimuth-range decoupling and introducing the TV regularization term into the sparse reconstruction model, our method significantly enhances the reconstruction accuracy of distributed targets.The method's effectiveness is substantiated through experimental results derived from both simulated data and real spaceborne SAR data.

Figure 1 .
Figure 1.PRI sequence values and blind areas' position in the simulation experiments.(a) The fastchanging PRI sequence values correspond to the parameters in Table 2, where the x-axis 'm' is the sequence number; (b) The blind areas' position of the fast-changing PRI sequence.

Figure 2 .
Figure 2. The 1-D imaging results with different methods.(a-c) are the imaging results of the point target located inside the blind areas, respectively (slant range = 956 km); (d-f) are the imaging results of the point target located at the boundary of the blind areas, respectively (slant range = 982 km); (g-i) are the imaging results of the point target located outside the blind areas, respectively (slant range = 994 km).

Figure 4 .
Figure 4.The imaging results of a one-dimensional distributed target with different methods, where the distributed target is located inside the blind areas (slant range = 956 km).(a) The ideal case; (b) Using the MF method; (c) Using the BLU interpolation method; and (d) Using the L1/2&TV-regularization-based method.

Figure 5 .
Figure 5.Comparison curves of the NRMSE reconstructed by BLU interpolation and L1/2&TV regularization under different SNRs.(a) Inside the blind areas (slant range = 956 km); (b) At the boundary of the blind areas (slant range = 982 km).

Figure 7 .Figure 8 .Figure 9 .
Figure 7.The imaging results of different methods in the first actual scenario.(a) Original echo data imaging processing under constant PRF condition; (b) Direct imaging on staggered-mode SAR data (MF); (c) Using the BLU interpolation method; (d) Using the L1/2&TV-regularization-based method.

Figure 10 .
Figure 10.The imaging results of different methods in the second scenario.(a) Original echo data imaging processing under constant PRF condition; (b) Direct imaging on staggered mode SAR data (MF); (c) Using the BLU interpolation method; (d) Using the L1/2&TV regularization-based method.
is the observed scene;  and  are the slow time in azimuth and fast time in range, respectively;  and  are the azimuth and range coordinate indexes of the target in the observation scene, respectively;   (⋅) is the azimuth antenna weight;  is the speed of light;  is the platform velocity;  is the wavelength of the transmitted signal; (, , ) is the slant range from the transmitting antenna to the target; and (⋅) is the transmitted pulse signal.Let   and   denote the number of points in the azimuth direction and range direction of the discrete grid, respectively.By discretizing the observed scene, the discrete model of (1) can be expressed as   (   ,    ) = ∑ ∑ (   ,    ) ⋅ [   ,    ,    ,    ]    ,    ) represents the discrete sampling of the   -th range cell of the original echo signal   (, ) at the   -th pulse; (   ,    ) represents the discrete value of the backscattering coefficient (, ) at the   -th azimuth grid and the   -th range grid; and [   ,    ,    ,    ] represents the value of the discretized observation matrix, reflecting the imaging geometric relationship between the radar and the observed scene:

Table 1 .
Pseudocode of the proposed sparse SAR imaging method.

Table 2 .
Major parameters in numerical simulation experiments.

Table
Quantitative calculation results of ambiguities suppression with different imaging methods.

Table 4 .
Quantitative calculation results of NRMSE with different imaging methods.

Table 5 .
Major system parameters of Gaofen-3 single-channel strip-map mode.

Table 6 .
ENL values of Region 1 to Region 5 with different imaging methods.