Synthetic aperture radar radio frequency interference suppression method, system, storage device and electronic equipment based on s transform

Abstract

The application discloses a synthetic aperture radar radio frequency interference suppression method and system based on S transform, a storage device and electronic equipment, and comprises the following steps: detecting interference-containing echo data by a method based on distance spectrum flatness and symmetry, and marking interference-containing pulses; calculating the amplitude spectrum of the interference-containing pulses, and performing smoothing and normalization processing; calculating the optimized window parameters after scaling the normalized amplitude spectrum; converting the interference-containing pulses to the time-frequency domain based on the optimized window parameters; positioning the interference position in the time-frequency domain through the OSTU algorithm, and eliminating the radio frequency interference through the wave-trap technology; inversely converting the processed interference-free pulses to the time domain based on the optimized window parameters; and combining all the processed interference-containing pulses and interference-free pulses to form new interference-free echo data. The application can effectively separate the radio frequency interference from the useful signal, thereby significantly reducing the influence of the radio frequency interference on the image quality of the synthetic aperture radar.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of signal processing, in particular to a method and system for suppressing radio frequency interference (RFI) of synthetic aperture radar (SAR) based on S transform, a computer readable storage medium and an electronic device. BACKGROUND

[0002] As an active microwave remote sensing device, synthetic aperture radar (SAR) can obtain high-precision ground target information at all times and in all weather. It realizes high resolution in azimuth and range directions through synthetic aperture and pulse compression technology. However, with the rapid development of modern communication technology, the spectrum resource is becoming increasingly scarce. As an open radar system, SAR has to share the spectrum with other devices in the same frequency band, which makes the electromagnetic environment more complex and the problem of radio frequency interference (RFI) more serious. These RFI problems pose a great threat to the SAR system and seriously affect the performance of SAR. On the one hand, the impulse response disturbance caused by RFI reduces the quality of SAR imaging, producing bright lines, fog-like artifacts and blurring in the image. On the other hand, the existence of RFI will lead to inaccurate spatial and radiometric measurements, and the resulting polarization and phase distortion will further hinder the business application of SAR images, such as agricultural investigation, land planning and ecological environment monitoring.

[0003] Under this background, time-frequency transform technology, as a method combining time and frequency information, shows its important role in interference feature extraction. Compared with traditional time-domain and frequency-domain suppression methods, time-frequency domain suppression methods can more accurately capture the interference characteristics in the time-frequency plane, thereby more effectively detecting and suppressing the interference. For this reason, time-frequency domain suppression has gradually become an important trend to deal with the problem of interference in complex electromagnetic environment.

[0004] Currently, the interference suppression method in time-frequency domain mainly converts the data to time-frequency domain by using time-frequency transform technology, and then detects and suppresses according to the characteristics of interference in time-frequency domain. For example, Zhang et al. designed an interference suppression filter based on the analysis of the time-frequency characteristics of NBI and WBI and the constant false alarm rate algorithm. Su et al. proposed an interference suppression algorithm based on time-frequency domain robust principal component analysis (RPCA), which first introduced the RPCA algorithm into time-frequency domain signal separation. Lyu et al. proposed a wideband complex parameter interference suppression algorithm based on signal time-frequency domain low-rank sparse matrix decomposition. Han et al. proposed a WBI suppression and single component extraction method combining instantaneous frequency estimation and regular time-frequency filtering. However, the above researches always choose simple and fast STFT when performing time-frequency transform, ignoring the influence of poor time-frequency resolution of STFT.

[0005] In summary, the existing time-frequency domain interference suppression method focuses on the research of interference suppression method after time-frequency transformation, ignoring the potential influence of time-frequency transformation itself on interference suppression effect. In order to ensure the quality of SAR image, the high-resolution time-frequency representation containing interference pulses can provide reliable technical support for SAR anti-jamming, and has important significance for improving the survival ability and practical efficiency of SAR system in complex electromagnetic environment. SUMMARY

[0006] The present application aims to provide a time-frequency domain synthetic aperture radar radio frequency interference suppression method, system, storage device and electronic equipment to solve or alleviate the problems existing in the prior art.

[0007] In order to achieve the above-mentioned purpose, the present application provides the following technical scheme:

[0008] The present application provides a time-frequency domain synthetic aperture radar radio frequency interference suppression method, comprising: step S101, applying a method based on distance spectrum flatness and symmetry to detect interference containing echo data, and marking out interference containing pulses; step S102, calculating the amplitude spectrum of interference containing pulses, and performing smoothing and normalization processing; step S103, after scaling the normalized amplitude spectrum, the optimized window parameter is calculated; step S104, based on the window parameter of step S103, the interference containing pulses are converted to time-frequency domain; step S105, the interference position in time-frequency domain is located by OSTU algorithm, and the radio frequency interference is eliminated by notch technology; step S106, the processed interference-free pulse is inverse transformed to time domain; step S107, all processed interference containing pulses and interference-free pulses are combined to form new interference-free echo data.

[0009] Preferably, in step S101, the method based on distance spectrum flatness and symmetry is applied to detect interference containing echo data, and the interference containing pulses are marked out, specifically:

[0010] The kurtosis and residual energy of each echo are calculated in a pulse-by-pulse manner to obtain a kurtosis value sequence and a residual energy sequence, and the K-means algorithm is used for classification.

[0011] Let K(η) and C(η) represent the radio frequency interference detection results based on flatness and symmetry indexes respectively, and multiply the classification results of flatness and symmetry. The result that meets flatness and symmetry is the radio frequency interference-free echo, and the rest is the radio frequency interference containing echo. This process can be represented as:

[0012]

[0013] Where Y(η) represents the radio frequency interference detection result, and the constant 1 represents that there is no interference in the echo, and the constant 0 represents that there is interference in the echo.

[0014] Preferably, in step S102, the amplitude spectrum of the interference pulse is calculated and smoothed and normalized, specifically:

[0015] The original echo signal x(t) is calculated, and the smoothed and normalized adjustment is made to obtain the normalized amplitude spectrum X(f) norm (f) is:

[0016] X(f) = smooth{abs[FT(x(t))]}

[0017]

[0018] Where smooth represents smoothing; FT represents Fourier transform.

[0019] Preferably, in step S103, the normalized amplitude spectrum is scaled to obtain the optimized window parameter, specifically:

[0020] X r (f) = [rX norm (f)] + 1

[0021] For SAR echo signals, the parameter r is usually 2 or 3.

[0022] Next, the window parameter p(f) varying with frequency is calculated to better match the signals of different frequency components:

[0023]

[0024] Where N r represents the number of distance sampling points.

[0025] Finally, the parameters α and β are introduced to further optimize the window parameter p(f), and the window parameter of the improved S transform is constructed:

[0026]

[0027] The improved window function w ′ (t,f) is:

[0028]

[0029] Preferably, in step S104, the interference pulse is converted to the time-frequency domain according to the window parameter of step S103, specifically:

[0030] For the radio frequency interference detection result Y(η) of step S101, set x η (t) as the ηth pulse data and Y(η) = 0, and its time-frequency domain representation is:

[0031]

[0032] where w(t) is the optimized window function in step S103.

[0033] Preferably, in step S105, the interference position in the time-frequency domain is located by the OSTU algorithm, and the radio frequency interference is eliminated by the notch technology, specifically:

[0034] For the time-frequency representation containing the interference pulse, the interference position matrix P i (t,f) is obtained by using the OSTU algorithm, and on this basis, the notch operation is performed, and the mathematical representation of this process is:

[0035]

[0036] that is, the result after interference suppression.

[0037] Preferably, in step S106, the processed interference-free pulse is inverse transformed to the time domain, specifically:

[0038]

[0039] The original signal is recovered using the parameters corresponding to the time-frequency transformation in step S103.

[0040] Preferably, in step S107, the interference-free echo data is recombined, specifically:

[0041] All processed interference-containing pulses and the pulses not marked in step S101 are combined to form new interference-free echo data.

[0042] The embodiment of the application also provides a method for suppressing radio frequency interference of a time-frequency domain synthetic aperture radar, comprising: an interference pulse detection and marking unit configured to detect and mark interference-containing pulses from synthetic aperture radar data input into a system; an amplitude spectrum calculation and preprocessing unit configured to obtain a normalized amplitude spectrum of original echo data; a window parameter optimization unit configured to calculate optimized window parameters after scaling the normalized amplitude spectrum; a time-frequency domain transformation unit configured to convert the interference-containing pulses to the time-frequency domain by using the optimized window parameters; an interference positioning and suppression unit configured to locate the interference position in the time-frequency domain and eliminate the radio frequency interference by using the notch technology; a time-frequency domain inverse transformation unit configured to inverse transform the processed interference-free pulses to the time domain by using the optimized window parameters; and an interference-free pulse recombination unit configured to recombine all processed interference-containing pulses with interference-free pulses, thereby obtaining new interference-free echo data.

[0043] ​The embodiment of the present application further provides a computer readable storage medium, which has a computer program stored thereon, and the program is the time-frequency domain synthetic aperture radar radio frequency interference suppression method as any of the above.

[0044] The embodiment of the present application further provides an electronic device, which comprises a memory, a processor, and a program stored in the memory and executable on the processor, and the processor implements the time-frequency domain synthetic aperture radar radio frequency interference suppression method as any of the above when executing the program.

[0045] Beneficial effects:

[0046] In the time-frequency domain synthetic aperture radar radio frequency interference suppression method provided by the present application, first, the interference pulse is detected and marked from the synthetic aperture radar data input into the system; the interference pulse amplitude spectrum is calculated and is subjected to smoothing and normalization processing; the normalized amplitude spectrum is scaled, and the optimized window parameter is calculated; the interference pulse is converted to the time-frequency domain; the interference position in the time-frequency domain is located, and the radio frequency interference is eliminated through the notch technology; the processed interference-free pulse is inversely converted to the time domain; finally, all the processed interference pulses and the interference-free pulses are combined to form new interference-free echo data. Through the present application, the limitations of low calculation efficiency and insufficient time-frequency focusing of the S transform in suppressing radio frequency interference can be overcome, and the present application can be widely applied to image processing in the synthetic aperture radar system. BRIEF DESCRIPTION OF DRAWINGS

[0047] The drawings constituting a part of the specification of the present application are used to provide further understanding of the present application, the illustrative embodiments of the present application and the description thereof are used to explain the present application, and do not constitute improper limitation on the present application. Among them:

[0048] Figure 1 It is a flowchart according to the present application;

[0049] Figure 2 It is an image feature of the synthetic aperture radar radio frequency interference signal according to some embodiments of the present application;

[0050] Figure 3 It is a different domain feature of the synthetic aperture radar radio frequency interference signal according to some embodiments of the present application;

[0051] Figure 4 It is a synthetic aperture radar radio frequency interference suppression result schematic diagram according to some embodiments of the present application;

[0052] Figure 5 It is a unit configuration diagram according to the present application. DETAILED DESCRIPTION

[0053] The present application will be described in detail below with reference to the attached drawings and embodiments. Various examples are provided by way of explanation of the present application but not to limit the present application. It will be apparent to those skilled in the art that modifications and variations can be made in the present application without departing from the scope or spirit of the present application. For example, features shown or described as part of one embodiment can be used in another embodiment to yield still a further embodiment. It is, therefore, desired that the present application be deemed as including all such modifications and variations as fall within the scope of the appended claims and their equivalents.

[0054] Exemplary method

[0055] As shown in the figure, the time-frequency domain synthetic aperture radar radio frequency interference suppression method comprises: Figure 1

[0056] Step S101, using a method based on distance spectrum flatness and symmetry, detecting the interference-containing echo data, and marking out the interference-containing pulses;

[0057] Specifically, the kurtosis and residual energy of each echo are calculated in a pulse-by-pulse manner to obtain a kurtosis value sequence and a residual energy sequence, and the kurtosis value sequence and the residual energy sequence are classified by K-means algorithm.

[0058] Let K(η) and C(η) represent the radio frequency interference detection results based on flatness and symmetry indicators respectively, and multiply the classification results of flatness and symmetry. The result that meets flatness and symmetry is the radio frequency interference-free echo, and the rest is the radio frequency interference-containing echo. This process can be represented as:

[0059]

[0060] Where Y(η) represents the radio frequency interference detection result, and the constant 1 represents that there is no interference in the echo, and the constant 0 represents that there is interference in the echo.

[0061] Step S102, calculating the interference-containing pulse amplitude spectrum, and performing smoothing and normalization processing;

[0062] The original echo signal x(t) is calculated, and the smoothing and normalization adjustment is performed to obtain the normalized amplitude spectrum X norm (f) is:

[0063] X(f)=smooth{abs[FT(x(t))]}

[0064]

[0065] In the formula, smooth represents smoothing processing; FT represents Fourier transform.

[0066] Step S103, after scaling the normalized amplitude spectrum, the optimized window parameters are calculated; specifically,​

[0067] X r (f) = [rX norm (f)] + 1

[0068] For SAR echo signal, parameter r usually takes 2 or 3.

[0069] Next, the window parameter p(f) varying with frequency is calculated to better match the signal of different frequency components:

[0070]

[0071] In the formula, N r represents the number of sampling points in the distance direction.

[0072] Finally, parameters a and β are introduced to further optimize the window parameter p(f), and the window parameter of the improved S transform is constructed:

[0073]

[0074] The improved window function w ′ (t,f) is:

[0075]

[0076] Step S104, convert the interference pulse to the time-frequency domain according to the window parameter of step S103; for the radio frequency interference detection result Y(η) of step S101, set x η (t) as the ηth pulse data and Y(η) = 0, and its time-frequency domain representation is:

[0077]

[0078] Wherein, w(t) is the window function optimized in step S103.

[0079] Step S105, locate the interference position in the time-frequency domain by the OSTU algorithm, and eliminate the radio frequency interference by the notch technology;

[0080] For the time-frequency representation of the interference pulse, the OSTU algorithm is used to obtain the interference position matrix P i (t,f), and on this basis, the notch operation is performed, and the mathematical representation of this process is:

[0081]

[0082] That is the result after interference suppression.

[0083] Step S106, inverse transform the processed interference-free pulse to the time domain.

[0084]

[0085] The original signal is recovered using parameters corresponding to the frequency transformation in step S103.

[0086] Step S107, all processed interference pulse-containing data A new interference-free echo data is formed with the pulses not marked in step S101.

[0087] The present application can distinguish different frequency components according to the amplitude of the actual signal by optimizing the window function, significantly improve the time-frequency resolution of S transform, improve the focusing of interference signals, and more accurately locate the interference signals, providing strong support for subsequent radio frequency interference suppression.

[0088] Figure 2 The imaging results of C-band Sentinel-1A raw data containing radio frequency interference are shown, which were collected on February 18, 2022, and the artifacts caused by radio frequency interference can be seen throughout the image.

[0089] Figure 3 The different domain features of the data are shown. Figure 2 The upper graph of figure (a) represents a two-dimensional time-domain amplitude image, in which the bright lines formed by radio frequency interference densely cover the SAR echo data. Figure 3 The lower graph of figure (a) represents a range-frequency domain-azimuth time domain amplitude image. Figure 3

[0090] The imaging results after using different interference suppression methods are shown. Figure 4 (a) shows the SAR image processed by the traditional method based on STFT, although most of the interference signals are suppressed, but there are still some interference residues. In addition, the details of the roads, buildings and other non-interference areas become blurred, indicating that the loss of useful signals is more serious. Figure 4 (b) is the result processed by the traditional method based on S transform, the interference residues are reduced, the details are relatively well preserved, but there is still some blurring phenomenon, and the overall performance is better than that of the STFT method. Figure 4 (c) shows the effect of the method proposed in the present application, the white interference stripes are almost completely suppressed, and the details such as roads and buildings in the non-interference area remain clear, indicating that the proposed method has more advantages in interference suppression and signal preservation. Figure 4 (d)~(f) are local magnified images of the region of interest in the red box, which further intuitively show the differences between the methods in suppressing interference and protecting useful signals.

[0091] Exemplary system

[0092] Figure 5 The application further provides a time-frequency domain synthetic aperture radar radio frequency interference suppression system, comprising: an interference pulse detection and marking unit configured to detect and mark interference pulses from synthetic aperture radar data input into the system; an amplitude spectrum calculation and preprocessing unit configured to obtain normalized amplitude spectrum of original echo data; a window parameter optimization unit configured to calculate optimized window parameters after scaling the normalized amplitude spectrum; a time-frequency domain conversion unit configured to convert the interference pulses to the time-frequency domain using the optimized window parameters; an interference positioning and suppression unit configured to position the interference in the time-frequency domain and eliminate the radio frequency interference through a notch technology; a time-frequency domain inverse conversion unit configured to inverse convert the processed interference-free pulses to the time domain using the optimized window parameters; and an interference-free echo recombination unit configured to recombine all the processed interference pulses and interference-free pulses to obtain new interference-free echo data.

[0093] The time-frequency domain synthetic aperture radar radio frequency interference suppression system provided by the application can implement any of the above time-frequency domain synthetic aperture radar radio frequency interference suppression steps and processes and achieve the same technical effects, which will not be described in detail here.

[0094] Exemplary apparatus

[0095] The application provides an electronic device, comprising a storage device, a processor, and a memory, wherein the processor is adapted to execute programs, the memory is used for storing programs, and the memory implements the time-frequency domain synthetic aperture radar radio frequency interference suppression method according to any one of claims 1-7 when executing the programs on the processor.

[0096] Since the time-frequency domain synthetic aperture radar radio frequency interference suppression steps have been described in detail in the specific implementation method examples, they will not be described in detail here.

[0097] The processor can include a central processing unit (CPU), a network processor (NP), and the like, and can also be a digital signal processor, an application-specific integrated circuit, a ready programmable gate array, or other programmable logic device, a discrete gate or transistor logic device, a discrete hardware component. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor.

[0098] The processor can be specifically configured to: detect and mark the interference pulse from the input system synthetic aperture radar data; calculate a normalized amplitude spectrum of the original echo data; calculate an optimized window parameter after scaling the normalized amplitude spectrum; convert the interference pulse to a time-frequency domain by using the optimized window parameter; locate the interference position in the time-frequency domain and eliminate the radio frequency interference by using a notch technology; inverse transform the processed interference-free pulse to the time domain by using the optimized window parameter; and recombine all the processed interference pulses and the interference-free pulses to obtain new interference-free echo data.

[0099] It should be noted that, according to the needs of implementation, each component / step described in the embodiments of the present application can be split into more components / steps, or two or more components / steps or part of the operation of the components / steps can be combined into a new component / step, to achieve the purpose of the embodiments of the present application.

[0100] The above method according to the embodiments of the present application can be implemented in hardware, firmware, or as software or computer code that can be stored in a recording medium such as a CD ROM, a RAM, a floppy disk, a hard disk or an optical disk, or be downloaded from a network and stored in a remote recording medium or a non-transitory machine storage medium and then stored in a local recording medium, so that the method described herein can be processed by such software on a recording medium using a general computer, a special processor or programmable or special hardware such as an ASIC or an FPGA. It can be understood that the computer, the processor, the microprocessor controller or the programmable hardware includes a storage component (for example, RAM, ROM, flash memory, etc.) that can store or receive software or computer code, when the software or computer code is accessed and executed by the computer, the processor or the hardware, the time-frequency domain synthetic aperture radar interference suppression method described herein is implemented. In addition, when the general computer accesses the code for implementing the method shown herein, the execution of the code will convert the general computer into a special computer for executing the method shown herein.

[0101] Those of ordinary skill in the art can realize that the units and method steps of the examples described in combination with the embodiments disclosed herein can be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether the functions are implemented in hardware or software depends on the specific application and involved constraints. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the embodiments of the present application.

[0102] It should be noted that each of the embodiments described in the specification of the present application adopts a progressive manner, and the same or similar parts between the embodiments can be mutually referred to. Each embodiment focuses on the differences from other embodiments. In particular, the device and system embodiments are described more simply because they are basically similar to the method embodiments, and the relevant parts can be referred to the part of the method embodiments.

[0103] The device and system embodiments described above are only schematic, and the units not shown in the description can or can not be physically separated, and the units prompted can or can not be physical units, i.e. they can be located in one place or distributed on multiple network units. Part or all of the modules can be selected to achieve the purpose of the embodiment according to actual needs. Those skilled in the art can understand and implement it without creative labor.

[0104] The above is only the preferred embodiment of the present application, and is not intended to limit the present application. Those skilled in the art can make various modifications and changes to the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for suppressing radio frequency interference in synthetic aperture radar based on S-transform, characterized in that, The method comprises the following steps: Step S101, detecting the interference-containing echo data based on the flatness and symmetry indicators of the distance spectrum, and marking the interference-containing pulses; Step S102, calculating the amplitude spectrum of the interference-containing pulses, and performing smoothing and normalization processing on the amplitude spectrum; Step S103, calculating the optimized window parameters after scaling the normalized amplitude spectrum; Step S104, converting the interference-containing pulses to the time-frequency domain based on the window parameters in step S103; Step S105, positioning the interference positions in the time-frequency domain through the OSTU algorithm, and eliminating the radio frequency interference through the notch technology; Step S106, inversely transforming the processed interference-free pulses to the time domain; Step S107, combining all the processed interference-containing pulses and the interference-free pulses to form new interference-free echo data.

2. The S-transform based synthetic aperture radar radio frequency interference mitigation method of claim 1, wherein, Step S101 specifically comprises the following steps: The kurtosis and residual energy of each echo are calculated in a pulse-by-pulse manner to obtain a kurtosis value sequence and a residual energy sequence, and the sequences are classified through the K-means algorithm; Let K(η) and C(η) represent the radio frequency interference detection results based on the flatness and symmetry indicators respectively, and multiply the classification results of the flatness and symmetry; the result that meets the flatness and symmetry is a radio frequency interference-free echo, and the rest is a radio frequency interference-containing echo, and the process can be represented as: Where Y(η) represents the radio frequency interference detection result, and the constant 1 indicates that there is no interference in the echo, and the constant 0 indicates that there is interference in the echo.

3. The S-transform based synthetic aperture radar radio frequency interference mitigation method of claim 1, wherein, Step S102 specifically comprises the following steps: The original echo signal x(t) is calculated by smoothing and normalizing adjustment to obtain the normalized amplitude spectrum X norm (f) that is: X(f)=smooth{abs[FT(x(t))]} In the formula, smooth represents smoothing processing; and FT represents Fourier transform.

4. The S-transform based synthetic aperture radar radio frequency interference mitigation method of claim 3, wherein, Step S103 specifically comprises the following steps: X r (f) = [rX norm (f)] + 1 For SAR echo signals, the parameter r is usually 2 or 3; Next, the window parameter p(f) varying with the frequency is calculated to better match the signals of different frequency components: In the formula, N r represents the number of sampling points in the distance direction; Finally, the parameters α and β are introduced to further optimize the window parameter p(f), and the window parameter of the improved S transform is constructed: The improved window function w ′ (t,f) is: 。 5. The S-transform based synthetic aperture radar radio frequency interference mitigation method of claim 1, wherein, Step S104 specifically comprises: For the radio frequency interference detection result Y(η) of step S101, set x η (τ) is the ηth pulse data and Y(η) = 0, then its time-frequency domain representation is: Where w(t-τ) is the optimized window function in step S103.

6. The S-transform based synthetic aperture radar radio frequency interference mitigation method of claim 5, wherein, Step S105 specifically comprises: For the time-frequency representation containing interference pulses, the OSTU algorithm is used to obtain the interference position matrix P i (t,f), on this basis, the notch operation is carried out, and the mathematical representation of the process is: i.e. non-interfering pulses.

7. The S-transform based synthetic aperture radar radio frequency interference mitigation method of claim 1, wherein, In step S106, the processed interference-free pulse inverse transform to time domain, specifically: The parameters corresponding to the time-frequency transformation in step S103 are used to restore the original signal.

8. A S-transform based synthetic aperture radar radio frequency interference suppression system, characterized in that, Comprise: An interference pulse detection marking unit configured to detect interference-containing echo data based on the flatness and symmetry indicators of the distance spectrum, and mark the interference-containing pulses; An amplitude spectrum calculation and preprocessing unit configured to obtain the normalized amplitude spectrum of the original echo data, calculate the amplitude spectrum of the interference-containing pulses, and perform smoothing and normalization processing on the amplitude spectrum; A window parameter optimization unit configured to calculate the optimized window parameters after scaling the normalized amplitude spectrum; A time-frequency domain transformation unit configured to convert the interference-containing pulses to the time-frequency domain by using the optimized window parameters; An interference positioning and suppression unit configured to position the interference positions in the time-frequency domain through the OSTU algorithm, and eliminate the radio frequency interference through the notch technology; A time-frequency domain inverse transformation unit configured to inversely transform the processed interference-free pulses to the time domain by using the optimized window parameters; The interference-free echo recombination unit is configured to recombine all the processed interference-containing pulses and interference-free pulses to obtain new interference-free echo data.

9. A storage device in which a plurality of programs are stored, characterized by The program is loaded and executed by the processor to implement the S transform based synthetic aperture radar radio frequency interference suppression method according to any one of claims 1-7.

10. An electronic device comprising a storage device, a processor; the processor is adapted to execute each program; the memory is used to store a plurality of programs; characterized in that, The memory implements the S transform based synthetic aperture radar radio frequency interference suppression method according to any one of claims 1-7 when the program on the processor is executed.

Images