Next Article in Journal
Lower Limb Motion Recognition with Improved SVM Based on Surface Electromyography
Previous Article in Journal
Portable Arduino-Based Multi-Sensor Device (SBEDAD): Measuring the Built Environment in Street Cycling Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Sparse Recovery Algorithm for Suppressing Multiple Linear Frequency Modulation Interference in the Synthetic Aperture Radar Image Domain

School of Electronic and Optical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(10), 3095; https://doi.org/10.3390/s24103095
Submission received: 22 January 2024 / Revised: 30 April 2024 / Accepted: 9 May 2024 / Published: 13 May 2024
(This article belongs to the Section Radar Sensors)

Abstract

:
In synthetic aperture radar (SAR) signal processing, compared with the raw data of level-0, level-1 SAR images are more readily accessible and available in larger quantities. However, an amount of level-1 images are affected by radio frequency interference (RFI), which typically originates from Linear Frequency Modulation (LFM) signals emitted by ground-based radars. Existing research on interference suppression in level-1 data has primarily focused on two methods: transforming SAR images into simulated echo data for interference suppression, or focusing interference in the frequency domain and applying notching filters to reduce interference energy. However, these methods overlook the effective utilization of the interference parameters or are confined to suppressing only one type of LFM interference at a time. In certain SAR images, multiple types of LFM interference manifest bright radiation artifacts that exhibit varying lengths along the range direction while remaining constant in the azimuth direction. It is necessary to suppress multiple LFM interference on SAR images when original echo data are unavailable. This article proposes a joint sparse recovery algorithm for interference suppression in the SAR image domain. In the SAR image domain, two-dimensional LFM interference typically exhibits differences in parameters such as frequency modulation rate and pulse width in the range direction, while maintaining consistency in the azimuth direction. Based on this observation, this article constructs a series of focusing operators for LFM interference in SAR images. These operators enable the sparse representation of dispersed LFM interference. Subsequently, an optimization model is developed that can effectively suppress multi-LFM interference and reduce image loss with the assistance of a regularization term in the image domain. Simulation experiments conducted in various scenarios validate the superior performance of the proposed method.

1. Introduction

Synthetic aperture radar (SAR) is an instrumental remote sensing technology renowned for its ability to provide high-resolution imaging capabilities across the globe, regardless of weather conditions. However, SAR imaging quality is susceptible to the adverse effects of diverse and intricate radio frequency interference (RFI). The interferences typically originate from non-intentional sources of electromagnetic devices operating within a similar frequency band, including neighboring satellites in adjacent orbits and ground broadcasting stations, etc. [1,2,3,4,5,6,7]. RFI suppression methods can be divided into two categories: RFI suppression methods based on level-0 products and those based on level-1 products [8].
Parametric, nonparametric, and semi-parametric methods are used for RFI suppression based on level-0 products [9]. The maximum likelihood estimation (MLE) algorithm [10], parameter maximum likelihood (PML) algorithm [11], least mean square (LMS) algorithm [12], and others estimate the parameters of RFI, reconstruct it, and subtract it from the SAR echo data. However, accurate estimation is required for effective RFI suppression. Nonparametric methods, such as the notched filter method, discard the frequency spectrum contaminated by narrowband RFI [13]. The eigen-subspace projection (ESP) method also removes the signal of interest (SOI) that overlaps with RFI in the frequency domain or eigen-subspace. In order to effectively preserve the SOI, semi-parametric methods based on sparse recovery (SR) and low-rank representation (LRR) have been extensively utilized for RFI suppression. In these methods, the sparsity or low-rank characteristics of RFI are incorporated as regularization terms and minimized through an optimization problem [14,15]. The optimization problem also models a regularization term for the SOI to prevent overfitting. Most methods based on SR or LRR focus on accurately capturing the characteristics of SOI through regularization terms. In Ultra-Wideband (UWB) SAR systems based on the pulse, the L 1 norm of the SOI in the time domain of raw data is utilized as a regularization term to ensure the imaging of strong scattering points [14,16,17,18,19,20]. In pulse-compression-based SAR systems, the L 1 norm of the distance profile is incorporated as a regularization term in the optimization model [21,22,23]. For semi-parametric methods, the exclusive reliance on the sparsity of the distance profile is inadequate in ensuring the precise recovery of the SOI.
Traditional algorithms mostly utilize level-0 SAR data, specifically in the raw data domain, for interference suppression, without extensively focusing on the characterization of RFI signals. Access to level-0 SAR data is often limited for users, leading them to primarily work with level-1 data [24,25,26,27,28,29]. As a result, it becomes imperative to prioritize the development of interference suppression techniques that are specifically designed for level-1 data.
In SAR systems, the primary source of RFI often stems from ground-based radars. These ground-based radars typically employ the emission of Linear Frequency Modulation (LFM) signals. For level-1 data affected by LFM interference, there are two categories of methods for interference suppression: those based on the simulation of SAR raw data using SAR images and those based on the analysis of SAR images themselves. In the study by Reigber et al. [30], the transformation of synthesized SAR images into a representation where common raw-data interference filtering methods can be applied is discussed. Building upon this, the work of [9] proposes a defocusing procedure to obtain pseudo-echo data and uses an improved ESP method and traditional frequency domain notching filter (FNF) for mitigating different scenarios of RFI. The research published in [31] introduces a robust principal component analysis (RPCA) method to directly remove RFI artifacts caused by LFM signals in single-look complex (SLC) SAR images. The work of [32] proposes a two-dimensional SPECtral ANalysis (2-D SPECAN) for removing LFM radar interference in spaceborne SAR images [33].
For methods based on the simulation of SAR raw data, the defocusing procedure and subsequent RFI suppression involve a complex process. It may lead to errors in polarimetric calibration due to the potential misestimation of calibration parameters in RFI-affected regions. Methods based on the analysis of SAR images adjust thresholds and discard real echoes at the same locations as the interference, resulting in the loss of target scenes. Interference suppression methods based on both categories are typically designed to address LFM interference from a single parameter category. However, in reality, SAR imaging covers a wide area, and there may be multiple ground-based radars emitting LFM interference with multiple parameters. Additionally, considering the long synthetic aperture of SAR, even if there is only one ground-based radar, its operational mode and waveform parameters may switch within the time span of the beam coverage [34,35]. Therefore, it is crucial to research and develop methods that can effectively suppress LFM interference with multiple parameters. These methods should be capable of mitigating the effects of various types of LFM interference to improve the overall quality of SAR imaging.
The work of [15] introduces a RFI suppression model that combines multiple image domain regularizations. RFI suppression is formulated as an optimization problem with regularization terms to characterize both RFI and the SOI. It aims to preserve different types of targets by incorporating various image domain regularizations. The generalized regularization-based RFI suppression model and solver serve as inspiration for the research presented in this paper.
In this paper, by simultaneously focusing on different parameters of interference in the level-1 SAR images, the method effectively suppresses multiple parameters of ground-based LFM interference. Furthermore, with the assistance of a terrain regularization term, it minimizes the loss of terrain information. To begin, the time–frequency characteristics of SAR images in the range and azimuth directions are observed using Short-Time Fourier Transform (STFT) [36,37,38]. It is found that LFM signals are still present in the SAR images, albeit with different parameters. Next, an analysis is conducted to determine the number of LFM interference types contained in SAR images and to identify the parameters for each LFM interference. Matching signals are then constructed in the range and azimuth directions for each LFM interference. Subsequently, a focusing operator is designed for LFM interference in the level-1 image domain, and a joint optimization model is formulated. This model simultaneously optimizes the degree of interference focusing and the loss of terrain information, enabling the recovery of interference signals and terrain images from contaminated images affected by interference.
The rest of this paper is organized as follows. Section 2 introduces the signal model in the echo domain and image domain. Section 3 shows the proposed the algorithm based on ADMM cyclic iteration. Section 4 introduces the experimental results as well as an analysis of different RFI suppression algorithms for point targets and surface targets. Finally, Section 5 concludes this paper.

2. Signal Model

2.1. Signal Model of LFM Interference in Echo Domain

In a SAR system, it is assumed that the SAR satellite travels along the azimuth direction with a velocity of v. Nearby interfering satellites transmit pulse LFM signals with a similar carrier frequency but with different parameters in multiple segments randomly. During the SAR motion, the received SAR signal is contaminated by LFM interference from other satellites. Assuming there are n ( n = 1 N ) types of interference signals, and each type consists of m ( m = 1 M ) pulses, after the carrier frequency removal, the model for the interference signal can be written as
s rfi ( η ; t ) = n = 1 N m = 1 M A m , n δ m , n ( η ) rect t τ m , n T n · exp j π k n ( t τ m , n ) 2 + j 2 π Δ f n ( t τ m , n )
where η denotes slow time in the azimuth dimension and t donates fast time in the range dimension. A m , n and τ m , n denote the amplitude and time delay of the m-th pulse of the n-th kind of LFM interference, respectively. τ m , n is calculated by dividing the slant range between the source of interference and the SAR antenna phase center by the speed of light. T n and k n are the pulse width and frequency modulation rate, respectively. Δ f n is the carrier frequency difference between the n-th kind of LFM interference signal and the SAR signal. The RFI exhibits random and short-duration characteristics. Considering that SAR operates in a step-by-step manner, in the azimuth direction, the interference signal can be represented as an impulse function. Therefore, it is represented using the Dirac function δ m , n .
The echo received after illuminating the targets with the transmitted signal from the SAR satellite can be written as
s 0 ( η ; t ) = i = 1 K A i a a ( η ) rect t 2 R i / c T 0 · exp j π k 0 ( t 2 R i c ) 2 exp j 4 π λ R i
where A i ( i = 1 K ) denotes the radar cross section (RCS) of each target. a a donates the azimuth window function. T 0 and k 0 donate the pulse width and frequency modulation rate of the transmitted signal, respectively. λ denotes the wavelength of the carrier signal. R i donates the slant range between the i-th target and the instantaneous position of the radar platform.
The model of the total signal can be written as
s ( η ; t ) = s 0 ( η ; t ) + s rfi ( η ; t ) + s z
where s z denotes the Gaussian white noise.

2.2. Signal Model of LFM Interference in Image Domain

The SAR received signal, after carrier removal and pulse compression processing, forms an L1 level Single-Look Complex (SLC) image. The reference signal for matched filtering can be written as
s ref ( η ; t ) = rect ( t T 0 ) exp ( j π k 0 t 2 ) a a ( η ) exp ( j π k a η 2 )
where k a denotes the frequency modulation rate in the azimuth dimension. After matched filtering, a target is represented by a sinc-like shape. The interference in the azimuth direction is modulated from a Dirac function to the LFM signal. The LFM signal has a frequency modulation opposite to the azimuth signal of s 0 ( η ; t ) . In the range direction, the interference and the reference signal of matched filtering do not match, resulting in changes in the frequency modulation rate and pulse width of the LFM interference. The LFM interference in the image domain can be written as
s rfi ( η ; t ) = n = 1 N m = 1 M A m , n rect η η m , n T a exp j π k a ( η η m , n ) 2 · rect t τ m , n T n exp j π k n ( t τ m , n ) 2 + j 2 π Δ f n ( t τ m , n )
where η and t are the azimuth time and range time in the image domain. T a denotes the synthetic aperture time, while k a is opposite to the azimuth frequency modulation rate of the SAR system. According to the work in [32], k n k 0 , k n = ( k n k 0 ) / ( k n k 0 ) and T n = T 0 | k n k 0 | / k 0 .
Casting the above model in matrix form allows the problem to be analyzed more clearly. Assuming there are two types of interference, denoted as R 1 and R 2 , the model for the SAR received signal in Formula (3) can be written as
S = S 0 + R 1 + R 2 + Z
S 0 = G 0 ( Y 0 )
where Y 0 represents the location of targets. Any non-zero element at a certain position in the matrix Y 0 indicates the presence of a target at that location. G 0 ( · ) denotes the transmission operator of the signal. The process of matched filtering can be represented by the operator G 0 1 ( · ) , and their specific forms can be written as
G 0 ( · ) = F a 1 F r 1 F a F r ( · ) Θ rc Θ rcmc Θ ac G 0 1 ( · ) = F r 1 F a 1 F r F a ( · ) Θ ac Θ rcmc Θ rc
where the operators F a ( · ) and F r ( · ) represent the Fourier transforms in the azimuth and range dimensions, respectively. The superscript −1 denotes the inverse transform. The range compression, azimuth compression, and bulk range cell migration correction (RCMC) operators are denoted as Θ rc , Θ ac , and Θ rcmc , respectively. The Hadamard product and division operators are represented by ⊙ and ⊘, respectively.
The received signal in Formula (6) is processed through matched filtering to obtain the level-1 SAR image, written as
Y = G 0 1 ( G 0 ( Y 0 ) + R 1 + R 2 + Z ) = Y 0 + G 0 1 ( R 1 ) + G 0 1 ( R 2 ) + G 0 1 ( Z )
Due to the effects of sidelobes, noise, and interference, Y is only an approximation of Y 0 . Similarly, G i 1 ( · ) ( i = 1 , 2 ) denotes the focusing process on the interference signal R i . Y i denotes the result of focusing on the interference of the SAR image Y . Taking i = 1 as an example, it yields the interference signal focused in the image domain, written as
Y 1 G 1 1 ( Y ) G 1 1 ( G 0 1 ( R 1 ) )
Under the assumption of sparsity in SAR images, when interference is focused in the image domain, its distribution also satisfies the sparsity condition. Therefore, the optimization problem of interference suppression can be written as
min Y 0 , Y 1 , Y 2 λ 0 Y 0 1 + λ 1 Y 1 1 + λ 2 Y 2 1 , s . t .     Y Y 0 G 1 ( Y 1 ) G 2 ( Y 2 ) F 2 < ϵ
where λ i ( i = 0 , 1 , 2 ) is the scalar parameters of Y i regularizations. L 1 norm ‘ · 1 ’ minimizes the sparsity of regular terms. The model’s objective is to minimize the sparsity of the images subject to the constraint that the Frobenius norm ‘ · F ’ of the difference between the original image and the reconstructed image is less than ϵ . The objective of this optimization problem is to suppress the impact of interference by minimizing the sum of sparse terms. The solution to the optimization problem will separate the focused interference ‘image’ from the reconstructed image with a lessened interference impact, thereby achieving interference suppression.

3. Sparse Recovery Algorithm Based on ADMM Optimization

This paper presents an algorithm for suppressing LFM interference in SAR images named LFMIS-JointSR. The algorithm utilizes an ADMM-based iterative approach to solve the joint sparse optimization problem, aiming to achieve more accurate estimation results for SAR images and focused images of different interference. The fundamental concept behind ADMM is to reframe the original problem into a sequence of subproblems by introducing Lagrange multipliers. In the case of the joint sparse optimization problem depicted by Formula (11), the Lagrangian function L is denoted as
L ( Y , Y 0 , Y 1 , Y 2 , U 0 , U 1 , U 2 , Φ 0 , Φ 1 , Φ 2 ) = 1 2 Y Y 0 G 1 ( Y 1 ) G 2 ( Y 2 ) F 2 + λ 0 U 0 1 + λ 1 U 1 1 + λ 2 U 2 1 + Φ 0 , U 0 Y 0 + Φ 1 , U 1 Y 1 + Φ 2 , U 2 Y 2 + μ 2 U 0 Y 0 F 2 + μ 2 U 1 Y 1 F 2 + μ 2 U 2 Y 2 F 2
where U i ( i = 0 , 1 , 2 ) denote the augmented variables. μ > 0 is a penalty parameter. Φ i = μ ( U i Y i ) denote the Lagrangian multipliers. By rewriting U i Y i in the penalty term U i Y i F 2 as U i Y i + Φ i μ Φ i μ , and then splitting it, combining the inner product terms, and dividing the Lagrangian function L by μ , we can obtain the new Lagrangian function, denoted as
L ( Y , Y 0 , Y 1 , Y 2 , U 0 , U 1 , U 2 , Φ 0 , Φ 1 , Φ 2 ) = 1 2 X G 1 ( Y 1 ) G 2 ( Y 2 ) Y 0 F 2 + λ 0 U 0 1 + λ 1 U 1 1 + λ 2 U 2 1 + 1 2 U 0 Y 0 + Φ 0 μ F 2 Φ 0 μ F 2 + 1 2 U 1 Y 1 + Φ 1 μ F 2 Φ 1 μ F 2 + 1 2 U 2 Y 2 + Φ 2 μ F 2 Φ 2 μ F 2
Next, the solution to the subproblems will be presented through the following steps of updating variables.
  • Optimize the augmented variables U i :
    By holding the other variables constant and assuming that they have been pre-estimated, the optimization problem derived from the Lagrangian function, as indicated in Formula (13), can be simplified as
    U 0 ( k + 1 ) = arg min U 0 λ 0 μ ( k ) U 0 ( k ) 1 + 1 2 U 0 ( k ) Y 0 ( k ) + Φ 0 ( k ) μ ( k ) F 2 U 1 ( k + 1 ) = arg min U 1 λ 1 μ ( k ) U 1 ( k ) 1 + 1 2 U 1 ( k ) Y 1 ( k ) + Φ 1 ( k ) μ ( k ) F 2 U 2 ( k + 1 ) = arg min U 2 λ 2 μ ( k ) U 2 ( k ) 1 + 1 2 U 2 ( k ) Y 2 ( k ) + Φ 2 ( k ) μ ( k ) F 2
    where the superscript k represents the iteration counter for the variables. By utilizing the soft threshold algorithm to solve the aforementioned optimization problems, we can obtain
    U 0 ( k + 1 ) = S T Y 0 ( k ) Φ 0 ( k ) μ ( k ) , λ 0 μ ( k ) U 1 ( k + 1 ) = S T Y 1 ( k ) Φ 1 ( k ) μ ( k ) , λ 1 μ ( k ) U 2 ( k + 1 ) = S T Y 2 ( k ) Φ 2 ( k ) μ ( k ) , λ 2 μ ( k )
    where
    S T A , b i , j = A i , j , if A i , j > b 0 , otherwise
  • Optimize Y i :
    Similarly, by holding the other variables, including the previously estimated U i , constant, the optimization problems for Y i can be written as
    Y 0 ( k + 1 ) = arg min Y 0 1 2 μ ( k ) Y Y 0 ( k ) G 1 ( Y 1 ( k ) ) G 2 ( Y 2 ( k ) ) F 2 + 1 2 U 0 ( k + 1 ) Y 0 ( k ) + Φ 0 ( k ) μ ( k ) F 2 Y 1 ( k + 1 ) = arg min Y 1 1 2 μ ( k ) Y Y 0 ( k + 1 ) G 1 ( Y 1 ( k ) ) G 2 ( Y 2 ( k ) ) F 2 + 1 2 U 1 ( k + 1 ) Y 1 ( k ) + Φ 1 ( k ) μ ( k ) F 2 Y 2 ( k + 1 ) = arg min Y 2 1 2 μ ( k ) Y Y 0 ( k + 1 ) G 1 ( Y 1 ( k + 1 ) ) G 2 ( Y 2 ( k ) ) F 2 + 1 2 U 2 ( k + 1 ) Y 2 ( k ) + Φ 2 ( k ) μ ( k ) F 2
    Based on the derivation in [39], it is demonstrated that the application of an operator on a matrix leaves its Frobenius norm unaffected. Setting the first-order derivative to zero yields the minimizers of Y i . The results are as follows:
    Y 0 ( k + 1 ) = 1 μ ( k ) + 1 μ ( k ) U 0 ( k + 1 ) + Φ 0 ( k ) + Y G 1 ( Y 1 ( k ) ) G 2 ( Y 2 ( k ) ) Y 1 ( k + 1 ) = 1 μ ( k ) + 1 μ ( k ) U 1 ( k + 1 ) + Φ 1 ( k ) + G 1 1 ( Y Y 0 ( k + 1 ) G 2 ( Y 2 ( k ) ) ) Y 2 ( k + 1 ) = 1 μ ( k ) + 1 μ ( k ) U 2 ( k + 1 ) + Φ 2 ( k ) + G 2 1 ( Y Y 0 ( k + 1 ) G 1 ( Y 1 ( k + 1 ) ) )
  • Update the Lagrange multipliers Φ i and penalty parameter μ :
    Respectively, the Lagrangian multipliers and the penalty parameter are updated via
    Φ 0 ( k + 1 ) = Φ 0 ( k ) + μ ( k ) ( U 0 ( k + 1 ) Y 0 ( k + 1 ) ) Φ 1 ( k + 1 ) = Φ 1 ( k ) + μ ( k ) ( U 1 ( k + 1 ) Y 1 ( k + 1 ) ) Φ 2 ( k + 1 ) = Φ 2 ( k ) + μ ( k ) ( U 2 ( k + 1 ) Y 2 ( k + 1 ) )
    and
    μ ( k + 1 ) = min ( μ max , ρ μ ( k ) )
    where the constant ρ > 1 is employed to facilitate the acceleration of convergence. μ max denotes the maximum permissible value of μ .
  • Termination criteria for iterations:
    The termination criteria for iterations include the following:
    • The maximum number of iterations, denoted as K max , is reached;
    • The objective function value between consecutive iterations falls below the predefined threshold, i.e., Y Y 0 G 1 ( Y 1 ) G 2 ( Y 2 ) F 2 / Y F 2 < ϵ .
The proposed ADMM algorithm can be summarized as follows:
U 0 ( k + 1 ) = S T Y 0 ( k ) Φ 0 ( k ) μ ( k ) , λ 0 μ ( k )
U 1 ( k + 1 ) = S T Y 1 ( k ) Φ 1 ( k ) μ ( k ) , λ 1 μ ( k )
U 2 ( k + 1 ) = S T Y 2 ( k ) Φ 2 ( k ) μ ( k ) , λ 2 μ ( k )
Y 0 ( k + 1 ) = 1 μ ( k ) + 1 μ ( k ) U 0 ( k + 1 ) + Φ 0 ( k ) + Y G 1 ( Y 1 ( k ) ) G 2 ( Y 2 ( k ) )
Y 1 ( k + 1 ) = 1 μ ( k ) + 1 μ ( k ) U 1 ( k + 1 ) + Φ 1 ( k ) + G 1 1 ( Y Y 0 ( k + 1 ) G 2 ( Y 2 ( k ) ) )
Y 2 ( k + 1 ) = 1 μ ( k ) + 1 μ ( k ) U 2 ( k + 1 ) + Φ 2 ( k ) + G 2 1 ( Y Y 0 ( k + 1 ) G 1 ( Y 1 ( k + 1 ) ) )
Φ 0 ( k + 1 ) = Φ 0 ( k ) + μ ( k ) ( U 0 ( k + 1 ) Y 0 ( k + 1 ) )
Φ 1 ( k + 1 ) = Φ 1 ( k ) + μ ( k ) ( U 1 ( k + 1 ) Y 1 ( k + 1 ) )
Φ 2 ( k + 1 ) = Φ 2 ( k ) + μ ( k ) ( U 2 ( k + 1 ) Y 2 ( k + 1 ) )
μ ( k + 1 ) = min ( μ max , ρ μ ( k ) )
The steps of the proposed algorithm are depicted in Algorithm 1.   
Algorithm 1: Multiple LFM Interference Suppression Algorithm Based on Joint Sparse Regularization
  Input: Complex-valued SAR image Y
  Output: Recovered complex-valued SAR image Y 0 , LFM interference image Y 1 and Y 2
1  Initialize k = 0 ;
2  Initialize U 0 ( 0 ) , U 1 ( 0 ) , U 2 ( 0 ) , Y 0 ( Y ) , Y 1 ( G 1 1 ( Y ) ) , Y 2 ( G 2 1 ( Y ) ) , Φ 0 ( 0 ) , Φ 1 ( 0 ) , Φ 2 ( 0 ) , μ ( μ 0 ) ;
3  Repeat
4      k = k + 1 ;
5      calculate U 0 ( k + 1 ) , U 1 ( k + 1 ) , U 2 ( k + 1 ) by (21)–(23);
6      calculate Y 0 ( k + 1 ) , Y 1 ( k + 1 ) , Y 2 ( k + 1 ) by (24)–(26);
7      calculate Φ 0 ( k + 1 ) , Φ 1 ( k + 1 ) , Φ 2 ( k + 1 ) by (27)–(29);
8      calculate μ by (30);
9  Until  k K max or Y Y 0 G 1 ( Y 1 ) G 2 ( Y 2 ) F 2 / Y F 2 < ϵ

4. Experiment and Analysis

In this section, the effectiveness of the proposed method will be validated using simulated SAR data. The simulation parameters are provided in Table 1. Several methods, namely FNF, PCA, RPCA, and LFMIS-JointSR, were utilized for conducting interference suppression in the simulation experiments. The objective was to compare their respective outcomes based on the target-to-background ratio (TBR) and mean square error (MSE). In Section 4.1, comparisons are conducted between various interference conditions and target scenarios, as shown in Table 2. The effectiveness of the proposed algorithms was demonstrated through the use of the TBR metric. In Section 4.2, the performance of different algorithms under different Signal-to-Interference Ratio (SIR) conditions was compared using the MSE metric, showcasing the robustness of the proposed algorithm. In Section 4.3, a quantitative comparison was carried out to better evaluate the proposed method. This comparison included a recently published similar work [32] as well as traditional algorithms.
In addition, it is worth noting that the signal and interference were processed in different domains. When LFM interference is subjected to matched filtering, the resulting signal exhibits a high sidelobe level when processed through a soft thresholding algorithm followed by an inverse-matched filtering transformation. Therefore, for the imaging operator of LFM interference, a Chirp Scaling (CS) imaging operator was experimentally adopted. Considering RFI 1 as an example, the imaging operator can be written as
G 1 ( · ) = F a 1 F r 1 F r [ F a ( · ) Θ rcmc ] Θ rc Θ ac Θ pc
G 1 ( · ) 1 = F a 1 F r 1 F r [ F a ( · ) Θ pc ] Θ ac Θ rc Θ rcmc
where Θ rcmc and Θ pc denote the RCMC and additional phase correction operators of the Chirp Scaling algorithm, respectively. Θ rc and Θ ac denote the range and azimuth compression operators, respectively.
MSE can be used to assess the discrepancy between the image obtained through the RFI suppression method and the original RFI-free image. It is defined as
MSE = 10 log 10 Y Y ^ F 2 Y F 2
where Y denotes the original RFI-free image. Y ^ denotes the recovered image after interference suppression. A smaller MSE serves as an indicator of superior performance in interference suppression.
TBR is commonly used to assess imaging quality [40]. It is defined as
TBR = 10 log 10 Y ^ ( T ) F 2 1 N B Y ^ ( B ) F 2
Respectively, T and B are the image patch containing a strong target and background, whose pixel values in the recovered image Y ^ are denoted as Y ^ ( T ) and Y ^ ( B ) . N B denotes the number of pixels in B . A higher TBR signifies enhanced performance in effectively recovering strong targets while minimizing the average sidelobe level.
In the first case, a point target and two LFM interferences, namely RFI 1 and RFI 2 , are present, as depicted in Figure 1a. The parameters of the transmitted signal, RFI 1 and RFI 2 , are provided in Table 1. The bandwidth and chirp duration set by RFI 1 and RFI 2 closely align with the transmitted signal, and the signal-to-noise-plus-interference ratio (SINR) is 2.5754 dB. The noise is Gaussian white noise. The processing results of the FNF, PCA, RPCA, and LFMIS-JointSR methods are displayed in Figure 1b–e. The FNF method effectively reduces the interference energy but also diminishes the energy of the target, resulting in a reduction in the target’s amplitude. The PCA algorithm is designed to retain the main signal components to suppress interference. It does not exhibit a satisfactory interference suppression effect in the presence of RFI 1 and RFI 2 . This is due to the similarity in parameters between the LFM interference and the transmitted signal. Both the RPCA and LFMIS-JointSR algorithms demonstrate interference suppression capabilities. Figure 1d–e reveal that the LFMIS-JointSR algorithm achieves superior interference suppression compared to other methods.
In the second case, three-point targets are present, as depicted in Figure 2a. Pulse LFM interferences with RFI 1 or RFI 2 ’s parameters, each consisting of three pulses, overlap with the echoes of the point targets. The SINR is 1.0234 dB. The FNF and PCA algorithms exhibit poor interference suppression effects, while the RPCA algorithm still exhibits residual interference energy. The LFMIS-JointSR algorithm consistently maintains interference suppression capabilities.

4.1. The Extended Target Case

The primary objective of simulation in the case of point targets is to validate the feasibility of the algorithm. In the case of interference suppression simulation with complex backgrounds, the main purpose is to evaluate the algorithm’s performance. In Figure 3a, a simulated SAR image is shown wherein the imaging of extended targets is overlapped with LFM interference. The SINR is −1.9719 dB. Figure 3b–e depict the results obtained by the FNF, PCA, RPCA, and LFMIS-JointSR methods, respectively. The results of the FNF and PCA algorithms indicate that, in the presence of complex backgrounds, neither in the frequency domain nor in the feature domain can interference and signal be effectively separated, resulting in suboptimal interference suppression.However, the RPCA and LFMIS-JointSR algorithms demonstrate the capability to suppress interference in complex backgrounds. Specifically, the LFMIS-JointSR algorithm exhibits a superior performance in separating signal and interference components accurately during interference suppression, thereby providing clearer and more reliable reconstructed images.

4.2. Robustness Analysis

The robustness of each algorithm was evaluated by comparing the resulting images after interference suppression under SIR. Figure 4a depicts a simulated SAR image with an SIR of −15 dB. In this image, there is interference present on both sides of the target point, labeled as RFI 1 or RFI 2 , respectively, significantly impacting the observation of the target.
Figure 4b–e show the results obtained using the FNF, PCA, RPCA, and LFMIS-JointSR methods, respectively. The FNF method partially removes the interference energy in the frequency domain. The interference and target signals are not completely separated in the frequency domain, leading to a reduction in target energy while residual interference energy remains. On the other hand, the PCA, RPCA, and LFMIS-JointSR methods all demonstrate interference suppression capabilities, indicating their robustness to interference. Among them, the LFMIS-JointSR algorithm exhibits the best interference suppression performance under low SIR conditions, highlighting its strong robustness.

4.3. Simulation with a Real SAR Image

In this section, we also conducted a quantitative comparison with a recently published similar work [32]. The work in [32] proposes a 2D algorithm that aims to suppress LFM interference in the SAR image domain. However, this algorithm is only effective for LFM interference with same parameters. To achieve optimal performance in this simulation, we manually executed the algorithm twice. As shown in Figure 5, the SAR image used in the experiment is the SLC image of the Sentinel-1 satellite. The image contains approximately 30 dBw of noise and numerous strong scattering points. We added RFI 1 and RFI 2 with a power of 40 dBw to the image. Figure 5c–g show the imaging results of FNF, PCA, RPCA, 2-D SPECAN, and LFMIS-JointSR algorithms, respectively.
Figure 5d–g show that the LFMIS-JointSR, RPCA, PCA, and 2-D SPECAN algorithms achieve good interference suppression effects, with LFMIS-JointSR and RPCA algorithms showing some noise suppression capabilities. The images’ TBRs are quantitatively compared and ranked in descending order as follows: LFMIS-JointSR achieves the highest TBR at 76.2689 dB, followed by RPCA at 73.6323 dB, PCA at 65.2427 dB, 2-D SPECAN at 65.2286 dB, and FNF at 60.3189 dB. The simulation results based on the real SAR image demonstrate that the proposed algorithm can effectively suppress multiple LFM interferences and noise in SAR image domain, and it outperforms recently published similar work.

5. Conclusions

In this paper, we focused on mitigating LFM interference in SLC images at the L1 level of SAR satellites. We commenced by establishing an echo signal model to concentrate the LFM interference on the SAR image. We then matrixized the process and replaced the dictionary design and sampling process of sparse signals with imaging operators. By incorporating augmented variables into the LFMIS-JointSR algorithm, we simplified the optimization problem of multiple variables into the task of finding a single variable multiple times. We provided a comprehensive description of the algorithm’s operation. Through simulations, we demonstrated that the proposed algorithm effectively focuses on various LFM interference and achieves interference suppression by concentrating and isolating the interference energy. Moreover, when compared with other interference suppression algorithms, the proposed algorithm exhibited superior interference suppression capabilities and robustness.

Author Contributions

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

Funding

This research was funded by the National Natural Science Foundation of China under Grant 62001229, 62101260 and 62101264, the Nature Science Foundation of Jiangsu Province under Grant BK20210334, BK20230915, the China Postdoctoral Science Foundation under Grant 2020M681604, and the Jiangsu Province Postdoctoral Science Foundation under Grant 2020Z441.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The real SAR image presented in the study are openly available in the Earth Science Data Systems (ESDS) Program at https://www.earthdata.nasa.gov/.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carrara, W.G.; Goodman, R.S.; Majewski, R.M. Spotlight Synthetic Aperture Radar: Signal Processing Algorithms; Artech House: London, UK, 1995. [Google Scholar]
  2. Shimada, M. L-band radio interferences observed by the JERS-1 SAR and its global distribution. In Proceedings of the 2005 IEEE International Geoscience and Remote Sensing Symposium, 2005. IGARSS’05, Seoul, Republic of Korea, 29 July 2005; IEEE: Piscataway, NJ, USA, 2005; Volume 4, pp. 2752–2755. [Google Scholar]
  3. Rosen, P.A.; Hensley, S.; Le, C. Observations and mitigation of RFI in ALOS PALSAR SAR data: Implications for the DESDynI mission. In Proceedings of the 2008 IEEE Radar Conference, Rome, Italy, 26–30 May 2008; IEEE: Piscataway, NJ, USA, 2008; pp. 1–6. [Google Scholar]
  4. Zhou, F.; Tao, M. Research on Methods for Narrow-Band Interference Suppression in Synthetic Aperture Radar Data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 3476–3485. [Google Scholar] [CrossRef]
  5. Yang, L.; Zheng, H.; Feng, J.; Li, N.; Chen, J. Detection and suppression of narrow band RFI for synthetic aperture radar imaging. Chin. J. Aeronaut. 2015, 28, 1189–1198. [Google Scholar] [CrossRef]
  6. Griffiths, H.; Cohen, L.; Watts, S.; Mokole, E.; Baker, C.; Wicks, M.; Blunt, S. Radar Spectrum Engineering and Management: Technical and Regulatory Issues. Proc. IEEE 2015, 103, 85–102. [Google Scholar] [CrossRef]
  7. Tao, M.; Su, J.; Huang, Y.; Wang, L. Mitigation of Radio Frequency Interference in Synthetic Aperture Radar Data: Current Status and Future Trends. Remote Sens. 2019, 11, 2438. [Google Scholar] [CrossRef]
  8. Li, N.; Zhang, H.; Lv, Z.; Min, L.; Guo, Z. Simultaneous Screening and Detection of RFI From Massive SAR Images: A Case Study on European Sentinel-1. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5231917. [Google Scholar] [CrossRef]
  9. Li, N.; Lv, Z.; Guo, Z. Observation and Mitigation of Mutual RFI Between SAR Satellites: A Case Study Between Chinese GaoFen-3 and European Sentinel-1A. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5112819. [Google Scholar] [CrossRef]
  10. Ulug, B. An Algorithm for Sinusoidal Interference Reduction Using Iterative Maximum Likelihood Estimation Techniques. Ph.D. Thesis, Ohio State University, Columbus, OH, USA, 1992. [Google Scholar]
  11. Vu, V.T.; Sjögren, T.K.; Pettersson, M.I.; Håkansson, L.; Gustavsson, A.; Ulander, L.M. RFI suppression in ultrawideband SAR using an adaptive line enhancer. IEEE Geosci. Remote Sens. Lett. 2010, 7, 694–698. [Google Scholar] [CrossRef]
  12. Lord, R.T.; Inggs, M.R. Efficient RFI suppression in SAR using LMS adaptive filter integrated with range/Doppler algorithm. Electron. Lett.-IEE 1999, 35, 629–630. [Google Scholar] [CrossRef]
  13. Mao, Y.; Huang, Y.; Yu, X.; Xin, Y.; Wang, Y.; Hong, W. An Radio Frequency Interference Mitigation Approach for Spaceborne SAR System in Low SINR Condition. IEEE Trans. Geosci. Remote Sens. 2023, 61, 5217414. [Google Scholar] [CrossRef]
  14. Huang, Y.; Liao, G.; Zhang, Z.; Xiang, Y.; Li, J.; Nehorai, A. Fast Narrowband RFI Suppression Algorithms for SAR Systems via Matrix-Factorization Techniques. IEEE Trans. Geosci. Remote Sens. 2019, 57, 250–262. [Google Scholar] [CrossRef]
  15. Lu, X.; Yang, J.; Yu, W.; Su, W.; Gu, H.; Yeo, T.S. Enhanced LRR-Based RFI Suppression for SAR Imaging Using the Common Sparsity of Range Profiles for Accurate Signal Recovery. IEEE Trans. Geosci. Remote Sens. 2021, 59, 1302–1318. [Google Scholar] [CrossRef]
  16. Huang, Y.; Liao, G.; Li, J.; Xu, J. Narrowband RFI suppression for SAR system via fast implementation of joint sparsity and low-rank property. IEEE Trans. Geosci. Remote Sens. 2018, 56, 2748–2761. [Google Scholar] [CrossRef]
  17. Nguyen, L.H.; Tran, T.D. Interference separation for UWB radar signals from entropy-driven robust PCA. In Proceedings of the 2017 IEEE Radar Conference (RadarConf), Seattle, WA, USA, 8–12 May 2017; IEEE: Piscataway, NJ, USA, 2017; pp. 0389–0393. [Google Scholar]
  18. Nguyen, L.H.; Tran, T.D. RFI-radar signal separation via simultaneous low-rank and sparse recovery. In Proceedings of the 2016 IEEE Radar Conference (RadarConf), Philadelphia, PA, USA, 2–6 May 2016; IEEE: Piscataway, NJ, USA, 2016; pp. 1–5. [Google Scholar]
  19. Huang, Y.; Liao, G.; Xu, J.; Li, J. Narrowband RFI suppression for SAR system via efficient parameter-free decomposition algorithm. IEEE Trans. Geosci. Remote Sens. 2018, 56, 3311–3322. [Google Scholar] [CrossRef]
  20. Ren, J.; Zhang, T.; Li, J.; Nguyen, L.H.; Stoica, P. RFI mitigation for UWB radar via hyperparameter-free sparse SPICE methods. IEEE Trans. Geosci. Remote Sens. 2018, 57, 3105–3118. [Google Scholar] [CrossRef]
  21. Lu, X.; Su, W.; Yang, J.; Gu, H.; Zhang, H.; Yu, W.; Yeo, T.S. Radio frequency interference suppression for SAR via block sparse Bayesian learning. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 4835–4847. [Google Scholar] [CrossRef]
  22. Liu, H.; Li, D.; Zhou, Y.; Truong, T.K. Joint wideband interference suppression and SAR signal recovery based on sparse representations. IEEE Geosci. Remote Sens. Lett. 2017, 14, 1542–1546. [Google Scholar] [CrossRef]
  23. Nguyen, L.H.; Dao, M.D.; Tran, T.D. Joint sparse and low-rank model for radio-frequency interference suppression in ultra-wideband radar applications. In Proceedings of the 2014 48th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 2–5 November 2014; IEEE: Piscataway, NJ, USA, 2014; pp. 864–868. [Google Scholar]
  24. Ali, I.; Cao, S.; Naeimi, V.; Paulik, C.; Wagner, W. Methods to Remove the Border Noise From Sentinel-1 Synthetic Aperture Radar Data: Implications and Importance For Time-Series Analysis. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 777–786. [Google Scholar] [CrossRef]
  25. Hodgson, M.E.; Jensen, J.R.; Schmidt, L.; Schill, S.; Davis, B. An evaluation of LIDAR- and IFSAR-derived digital elevation models in leaf-on conditions with USGS Level 1 and Level 2 DEMs. Remote Sens. Environ. 2003, 84, 295–308. [Google Scholar] [CrossRef]
  26. Filipponi, F. Sentinel-1 GRD Preprocessing Workflow. Proceedings 2019, 18, 11. [Google Scholar] [CrossRef]
  27. Chojka, A.; Artiemjew, P.; Rapiński, J. RFI Artefacts Detection in Sentinel-1 Level-1 SLC Data Based On Image Processing Techniques. Sensors 2020, 20, 2919. [Google Scholar] [CrossRef] [PubMed]
  28. Kellndorfer, J.; Pierce, L.; Dobson, M.; Ulaby, F. Toward consistent regional-to-global-scale vegetation characterization using orbital SAR systems. IEEE Trans. Geosci. Remote Sens. 1998, 36, 1396–1411. [Google Scholar] [CrossRef]
  29. Tao, M.; Lai, S.; Li, J.; Su, J.; Fan, Y.; Wang, L. Extraction and Mitigation of Radio Frequency Interference Artifacts Based on Time-Series Sentinel-1 SAR Data. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5217211. [Google Scholar] [CrossRef]
  30. Reigber, A.; Ferro-Famil, L. Interference suppression in synthesized SAR images. IEEE Geosci. Remote Sens. Lett. 2005, 2, 45–49. [Google Scholar] [CrossRef]
  31. Yang, H.; Tao, M.; Chen, S.; Xi, F.; Liu, Z. On the mutual interference between spaceborne SARs: Modeling, characterization, and mitigation. IEEE Trans. Geosci. Remote Sens. 2020, 59, 8470–8485. [Google Scholar] [CrossRef]
  32. Yang, H.; He, Y.; Du, Y.; Zhang, T.; Yin, J.; Yang, J. Two-dimensional spectral analysis filter for removal of LFM radar interference in spaceborne SAR imagery. IEEE Trans. Geosci. Remote Sens. 2021, 60, 5219016. [Google Scholar] [CrossRef]
  33. Cumming, I.G.; Wong, F.H. Digital processing of synthetic aperture radar data. Artech House 2005, 1, 108–110. [Google Scholar]
  34. Kari, S.S.; Raj, A.A.B.; Balasubramanian, K. Evolutionary Developments of Today’s Remote Sensing Radar Technology—Right From the Telemobiloscope: A review. IEEE Geosci. Remote Sens. Mag. 2023, 12, 67–107. [Google Scholar] [CrossRef]
  35. Zhang, S.; Lu, X.; Yu, J.; Dai, Z.; Su, W.; Gu, H. Clutter Suppression for Radar via Deep Joint Sparse Recovery Network. IEEE Geoscience and Remote Sensing Letters 2024, 21, 3332035. [Google Scholar] [CrossRef]
  36. Sturmel, N.; Daudet, L. Signal reconstruction from STFT magnitude: A state of the art. In Proceedings of the International Conference on Digital Audio Effects (DAFx), Paris, France, 19–23 September 2011; pp. 375–386. [Google Scholar]
  37. Průša, Z.; Balazs, P.; Søndergaard, P.L. A Noniterative Method for Reconstruction of Phase From STFT Magnitude. IEEE/ACM Trans. Audio Speech Lang. Process. 2017, 25, 1154–1164. [Google Scholar] [CrossRef]
  38. Kim, B.; Kong, S.H.; Kim, S. Low Computational Enhancement of STFT-Based Parameter Estimation. IEEE J. Sel. Top. Signal Process. 2015, 9, 1610–1619. [Google Scholar] [CrossRef]
  39. Lu, X.; Yang, J.; Yeo, T.S.; Su, W.; Gu, H.; Yu, W. Accurate SAR Image Recovery From RFI Contaminated Raw Data by Using Image Domain Mixed Regularizations. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5102813. [Google Scholar] [CrossRef]
  40. Cetin, M.; Karl, W.C.; Castanon, D.A. Evaluation of a regularized SAR imaging technique based on recognition-oriented features. In Proceedings of the Algorithms for Synthetic Aperture Radar Imagery VII. SPIE, Orlando, FL, USA, 24–28 April 2000; Volume 4053, pp. 40–51. [Google Scholar]
Figure 1. Comparison of different algorithms for suppressing LFM interference in single-point target case. (a) The original SAR image. (be) Resulting images after interference suppression using FNF, PCA, RPCA, and LFMIS-JointSR algorithms, respectively.
Figure 1. Comparison of different algorithms for suppressing LFM interference in single-point target case. (a) The original SAR image. (be) Resulting images after interference suppression using FNF, PCA, RPCA, and LFMIS-JointSR algorithms, respectively.
Sensors 24 03095 g001
Figure 2. Comparison of different algorithms for suppressing LFM interference in multiple point targets and interference case. (a) The origrinal SAR image. (be) Resulting images after interference suppression using FNF, PCA, RPCA, and LFMIS-JointSR algorithms, respectively.
Figure 2. Comparison of different algorithms for suppressing LFM interference in multiple point targets and interference case. (a) The origrinal SAR image. (be) Resulting images after interference suppression using FNF, PCA, RPCA, and LFMIS-JointSR algorithms, respectively.
Sensors 24 03095 g002
Figure 3. Comparison of different algorithms for suppressing LFM interference in extended target and multiple interference case. (a) The original SAR image. (be) Resulting images after interference suppression using FNF, PCA, RPCA, and LFMIS-JointSR algorithms, respectively.
Figure 3. Comparison of different algorithms for suppressing LFM interference in extended target and multiple interference case. (a) The original SAR image. (be) Resulting images after interference suppression using FNF, PCA, RPCA, and LFMIS-JointSR algorithms, respectively.
Sensors 24 03095 g003
Figure 4. Comparison of different algorithms for suppressing LFM interference in low-SIR case. (a) The original SAR image. (be) Resulting images after interference suppression using FNF, PCA, RPCA, and LFMIS-JointSR algorithms, respectively.
Figure 4. Comparison of different algorithms for suppressing LFM interference in low-SIR case. (a) The original SAR image. (be) Resulting images after interference suppression using FNF, PCA, RPCA, and LFMIS-JointSR algorithms, respectively.
Sensors 24 03095 g004
Figure 5. Comparison of different algorithms for suppressing LFM interference in real SAR image. (a) The original SAR image. (b) The SAR image with RFI. (cg) Resulting images after interference suppression using FNF, PCA, RPCA, 2-D SPECAN, and LFMIS-JointSR algorithms, respectively.
Figure 5. Comparison of different algorithms for suppressing LFM interference in real SAR image. (a) The original SAR image. (b) The SAR image with RFI. (cg) Resulting images after interference suppression using FNF, PCA, RPCA, 2-D SPECAN, and LFMIS-JointSR algorithms, respectively.
Sensors 24 03095 g005
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParameterValueParameterValue
Transmitted waveformLFM signalCarrier frequency5.405 GHz
Signal bandwidth1 MHzPulse repetition frequency2500 KHz
Chirp duration5 μsRange between centre of the flight track and scene center887.23 km
Platform velocity7249 m s−1Synthetic aperture time0.1359 s
Theoretical range resolution150 mTheoretical azimuth resolution25 m
Signal bandwidth of RFI 1 1.7 MHz *Chirp duration of RFI 1 10 μs
Signal bandwidth of RFI 2 1.7 MHz *Chirp duration of RFI 2 10 μs
* The frequency of RFI 1 increases, while the frequency of RFI 2 decreases.
Table 2. TBRs of the FNF, PCA, RPCA, and LFMIS-JointSR for interference suppression results in simulation experiments.
Table 2. TBRs of the FNF, PCA, RPCA, and LFMIS-JointSR for interference suppression results in simulation experiments.
CaseMethodTBR (dB)
The first-point target caseFNF51.9474
PCA51.4356
RPCA69.1947
LFMIS-JointSR154.0224
The second-point target caseFNF43.8465
PCA42.0622
RPCA56.8874
LFMIS-JointSR154.4594
The extended target caseFNF50.0983
PCA54.4416
RPCA70.6983
LFMIS-JointSR144.1300
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

Tong, G.; Lu, X.; Yang, J.; Yu, W.; Gu, H.; Su, W. A Sparse Recovery Algorithm for Suppressing Multiple Linear Frequency Modulation Interference in the Synthetic Aperture Radar Image Domain. Sensors 2024, 24, 3095. https://doi.org/10.3390/s24103095

AMA Style

Tong G, Lu X, Yang J, Yu W, Gu H, Su W. A Sparse Recovery Algorithm for Suppressing Multiple Linear Frequency Modulation Interference in the Synthetic Aperture Radar Image Domain. Sensors. 2024; 24(10):3095. https://doi.org/10.3390/s24103095

Chicago/Turabian Style

Tong, Guanqi, Xingyu Lu, Jianchao Yang, Wenchao Yu, Hong Gu, and Weimin Su. 2024. "A Sparse Recovery Algorithm for Suppressing Multiple Linear Frequency Modulation Interference in the Synthetic Aperture Radar Image Domain" Sensors 24, no. 10: 3095. https://doi.org/10.3390/s24103095

APA Style

Tong, G., Lu, X., Yang, J., Yu, W., Gu, H., & Su, W. (2024). A Sparse Recovery Algorithm for Suppressing Multiple Linear Frequency Modulation Interference in the Synthetic Aperture Radar Image Domain. Sensors, 24(10), 3095. https://doi.org/10.3390/s24103095

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