CN111580107B - Radio frequency interference suppression method and device for SAR echo signal and imaging method - Google Patents
Radio frequency interference suppression method and device for SAR echo signal and imaging method Download PDFInfo
- Publication number
- CN111580107B CN111580107B CN202010640070.9A CN202010640070A CN111580107B CN 111580107 B CN111580107 B CN 111580107B CN 202010640070 A CN202010640070 A CN 202010640070A CN 111580107 B CN111580107 B CN 111580107B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- notch
- frequency spectrum
- processing
- echo signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000003384 imaging method Methods 0.000 title claims abstract description 17
- 230000001629 suppression Effects 0.000 title claims abstract description 14
- 238000001228 spectrum Methods 0.000 claims abstract description 264
- 238000012545 processing Methods 0.000 claims abstract description 80
- 238000001514 detection method Methods 0.000 claims abstract description 33
- 238000013213 extrapolation Methods 0.000 claims abstract description 22
- 238000001914 filtration Methods 0.000 claims abstract description 18
- 238000005070 sampling Methods 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000009466 transformation Effects 0.000 claims description 4
- 238000011084 recovery Methods 0.000 claims description 3
- 238000012217 deletion Methods 0.000 claims 1
- 230000037430 deletion Effects 0.000 claims 1
- 230000003595 spectral effect Effects 0.000 abstract 4
- 230000008569 process Effects 0.000 description 11
- 238000010276 construction Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 1
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 1
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种SAR回波信号的射频干扰抑制方法、装置及成像方法,SAR回波信号的射频干扰抑制方法包括:对接收到的初始回波信号进行距离向傅里叶变换,获得初始频谱;基于预设的陷波干扰检测门限值对初始频谱进行陷波处理,获得陷波处理后的频谱;对陷波处理后的频谱进行滤波处理,获得滤波后的频谱;基于预设的预测模型对滤波处理后的频谱进行频谱外推处理,获得前向外推频谱数据以及后向外推频谱数据;基于前向外推频谱数据、后向外推频谱数据以及滤波后的频谱进行重建处理,获得无缺失的目标频谱;基于无缺失的目标频谱进行距离向信号的调频信息的恢复,获得抑制干扰且频谱完整的目标回波信号。本发明能够准确的对射频干扰进行抑制。
The invention discloses a radio frequency interference suppression method, device and imaging method of a SAR echo signal. The radio frequency interference suppression method of a SAR echo signal includes: performing range-to-Fourier transform on the received initial echo signal to obtain an initial Spectrum; performing notch processing on the initial spectrum based on the preset notch interference detection threshold to obtain a notch-processed spectrum; filtering the notch-processed spectrum to obtain a filtered spectrum; based on the preset The prediction model performs spectrum extrapolation processing on the filtered spectrum to obtain forward extrapolated spectral data and backward extrapolated spectral data; reconstructs based on forward extrapolated spectral data, backward extrapolated spectral data and filtered spectrum processing to obtain the target spectrum without loss; based on the target spectrum without loss, the frequency modulation information of the range signal is recovered, and the target echo signal with suppressed interference and complete spectrum is obtained. The invention can accurately suppress the radio frequency interference.
Description
技术领域Technical Field
本发明涉及航天和微波遥感的交叉技术领域,特别涉及一种SAR回波信号的射频干扰抑制方法、装置及成像方法。The present invention relates to the cross-technical field of aerospace and microwave remote sensing, and in particular to a method and device for suppressing radio frequency interference of SAR echo signals and an imaging method.
背景技术Background Art
SAR(Synthetic Aperture Radar,合成孔径雷达)技术是一种主动的航天、航空遥感手段,微波成像技术具有全天时、全天候工作特点,在环境保护、灾害监测、海洋观测、资源勘察、精细农业、地质测绘等方面有着广泛的应用,目前已成为高分对地观测和全球资源管理的最重要手段之一。然而,SAR在复杂电磁环境中容易受到RFI(Radio FrequencyInterference,射频干扰)的影响,尤其是低波段(L波段、P波段),进而降低图像质量,影响目标检测和图像识别。SAR (Synthetic Aperture Radar) technology is an active aerospace and aviation remote sensing method. Microwave imaging technology has the characteristics of working all day and all weather. It has been widely used in environmental protection, disaster monitoring, ocean observation, resource exploration, precision agriculture, geological mapping, etc., and has become one of the most important means of high-resolution earth observation and global resource management. However, SAR is easily affected by RFI (Radio Frequency Interference) in complex electromagnetic environments, especially in low bands (L band, P band), which reduces image quality and affects target detection and image recognition.
低频SAR的工作频率为VHF(Very High Frequency,甚高频)或者UHF(Ultra HighFrequency,超高频)频段,由于该频段被ITU(International Telecommunications Union,国际电信联盟)首先分配给地面广播电视系统、通信系统、GPS(Global PositioningSystem,全球定位系统)和地面警戒雷达,即使在SAR系统设计时考虑了诸多抗干扰技术,但是仍有大量的干扰信号进入到SAR系统接受机,最终导致SAR系统成像质量下降。The operating frequency of low-frequency SAR is VHF (Very High Frequency) or UHF (Ultra High Frequency) band. Since this band was first allocated by ITU (International Telecommunications Union) to terrestrial broadcasting and television systems, communication systems, GPS (Global Positioning System) and ground warning radar, even if many anti-interference technologies are considered in the design of the SAR system, there are still a large number of interference signals entering the SAR system receiver, which ultimately leads to a decrease in the imaging quality of the SAR system.
现有的RFI抑制方法主要有:频域陷波法、子空间分解法、干扰提取法等。由于传统的频域陷波法容易造成有用信号的缺失,而子空间投影法过分依赖于子空间的正交性,当干扰信号功率较大时,会造成子空间误判。并且干扰提取法计算复杂度过高,射频干扰的数学模型相当复杂,而且会出现模型失配,产生模型参数估计误差,进而导致RFI估计不准确。The existing RFI suppression methods mainly include: frequency domain notch method, subspace decomposition method, interference extraction method, etc. Since the traditional frequency domain notch method easily causes the loss of useful signals, and the subspace projection method is overly dependent on the orthogonality of the subspace, when the interference signal power is large, it will cause subspace misjudgment. In addition, the calculation complexity of the interference extraction method is too high, the mathematical model of radio frequency interference is quite complex, and there will be model mismatch, resulting in model parameter estimation errors, which will lead to inaccurate RFI estimation.
发明内容Summary of the invention
本发明实施例的目的在于提供一种SAR回波信号的射频干扰抑制方法、装置及成像方法,用于解决现有技术中射频干扰抑制不够准确的问题。The purpose of the embodiments of the present invention is to provide a method, device and imaging method for radio frequency interference suppression of SAR echo signals, so as to solve the problem that radio frequency interference suppression is not accurate enough in the prior art.
为了解决上述技术问题,本申请的实施例采用了如下技术方案:一种SAR回波信号的射频干扰抑制方法,包括如下步骤:In order to solve the above technical problems, the embodiments of the present application adopt the following technical solutions: a method for suppressing radio frequency interference of SAR echo signals, comprising the following steps:
对接收到的初始回波信号进行距离向傅里叶变换,获得与所述初始回波信号对应的初始频谱;Performing a range-direction Fourier transform on the received initial echo signal to obtain an initial frequency spectrum corresponding to the initial echo signal;
基于预设的陷波干扰检测门限值对所述初始频谱进行陷波处理,获得陷波处理后的频谱;Performing notch processing on the initial spectrum based on a preset notch interference detection threshold value to obtain a spectrum after notch processing;
对所述陷波处理后的频谱进行滤波处理,获得滤波后的频谱;Performing filtering on the frequency spectrum after the notch processing to obtain a filtered frequency spectrum;
基于预设的预测模型对所述滤波处理后的频谱进行频谱外推处理,获得前向外推频谱数据以及后向外推频谱数据;Performing spectrum extrapolation processing on the filtered spectrum based on a preset prediction model to obtain forward extrapolated spectrum data and backward extrapolated spectrum data;
基于所述前向外推频谱数据、所述后向外推频谱数据以及所述滤波后的频谱进行重建处理,获得无缺失的目标频谱;Performing reconstruction processing based on the forward extrapolated spectrum data, the backward extrapolated spectrum data and the filtered spectrum to obtain a target spectrum without missing components;
基于所述无缺失的目标频谱进行距离向信号的调频信息的恢复,获得抑制干扰且频谱完整的目标回波信号。The frequency modulation information of the range signal is restored based on the target spectrum without missing any, so as to obtain a target echo signal with suppressed interference and complete spectrum.
可选的,所述基于预设的陷波干扰检测门限值对所述初始频谱进行陷波处理,获得陷波处理后的频谱,具体包括:Optionally, performing notch processing on the initial spectrum based on a preset notch interference detection threshold value to obtain a spectrum after notch processing specifically includes:
基于预设的第一陷波干扰检测门限值对所述初始频谱进行陷波处理,获得第一频谱;Performing notch processing on the initial spectrum based on a preset first notch interference detection threshold value to obtain a first spectrum;
基于预设的第二陷波干扰检测门限值对所述第一频谱进行陷波处理,获得所述陷波处理后的频谱。The first frequency spectrum is notched based on a preset second notch interference detection threshold value to obtain the notched frequency spectrum.
可选的,所述基于所述前向外推频谱数据、所述后向外推频谱数据以及所述滤波后的频谱进行重建处理,获得无缺失的目标频谱,具体包括:Optionally, the reconstructing based on the forward extrapolated spectrum data, the backward extrapolated spectrum data and the filtered spectrum to obtain a target spectrum without missing specifically includes:
基于所述前向外推频谱以及所述后向外推频谱进行缺失频谱重建,获得缺失频谱;Reconstructing the missing spectrum based on the forward extrapolated spectrum and the backward extrapolated spectrum to obtain the missing spectrum;
基于所述滤波后的频谱以及所述缺失频谱进行拼接处理,获得所述无缺失的目标频谱。A splicing process is performed based on the filtered spectrum and the missing spectrum to obtain the target spectrum without missing anything.
可选的,所述预测模型包括AR线性预测模型;所述方法还包括构建所述AR线性预测模型,具体包括:Optionally, the prediction model includes an AR linear prediction model; the method further includes constructing the AR linear prediction model, specifically including:
基于距离向采样点的个数确定AR线性预测模型的阶数;Determine the order of the AR linear prediction model based on the number of distance sampling points;
基于预设的计算方法计算所述AR线性预测模型的系数;Calculating the coefficients of the AR linear prediction model based on a preset calculation method;
基于所述阶数和系数构建所述AR线性预测模型。The AR linear prediction model is constructed based on the order and coefficients.
可选的,所述预设方法包括如下一种或几种:伯格算法和改进的协方差算法。Optionally, the preset method includes one or more of the following: a Berg algorithm and an improved covariance algorithm.
本申请的实施例提供一种SAR回波信号的射频干扰抑制装置,包括:An embodiment of the present application provides a SAR echo signal radio frequency interference suppression device, comprising:
变换模块,用于对接收到的初始回波信号进行距离向傅里叶变换,获得与所述初始回波信号对应的初始频谱;A transformation module, used for performing a range-oriented Fourier transform on the received initial echo signal to obtain an initial spectrum corresponding to the initial echo signal;
陷波处理模块,用于基于预设的陷波干扰检测门限值对所述初始频谱进行陷波处理,获得陷波处理后的频谱;A notch processing module, used to perform notch processing on the initial spectrum based on a preset notch interference detection threshold value to obtain a spectrum after notch processing;
滤波处理模块,用于对所述陷波处理后的频谱进行滤波处理,获得滤波后的频谱;A filtering processing module, used for filtering the spectrum after the notch processing to obtain a filtered spectrum;
频谱外推处理模块,用于基于预设的预测模型对所述滤波处理后的频谱进行频谱外推处理,获得前向外推频谱数据以及后向外推频谱数据;A spectrum extrapolation processing module, used for performing spectrum extrapolation processing on the spectrum after filtering based on a preset prediction model to obtain forward extrapolated spectrum data and backward extrapolated spectrum data;
重建模块,用于基于所述前向外推频谱数据、所述后向外推频谱数据以及所述滤波后的频谱进行重建处理,获得无缺失的目标频谱;A reconstruction module, used for performing reconstruction processing based on the forward extrapolated spectrum data, the backward extrapolated spectrum data and the filtered spectrum to obtain a target spectrum without missing components;
恢复模块,用于基于所述无缺失的目标频谱进行距离向信号的调频信息的恢复,获得抑制干扰且频谱完整的目标回波信号。The recovery module is used to recover the frequency modulation information of the range signal based on the target spectrum without missing any, so as to obtain the target echo signal with suppressed interference and complete spectrum.
可选的,所述陷波处理模块具体用于:Optionally, the notch processing module is specifically used for:
基于预设的第一陷波干扰检测门限值对所述初始频谱进行陷波处理,获得第一频谱;Performing notch processing on the initial spectrum based on a preset first notch interference detection threshold value to obtain a first spectrum;
基于预设的第二陷波干扰检测门限值对所述第一频谱进行陷波处理,获得所述陷波处理后的频谱。The first frequency spectrum is notched based on a preset second notch interference detection threshold value to obtain the notched frequency spectrum.
可选的,所述重建模块具体用于:Optionally, the reconstruction module is specifically used for:
基于所述前向外推频谱以及所述后向外推频谱进行缺失频谱重建,获得缺失频谱;Reconstructing the missing spectrum based on the forward extrapolated spectrum and the backward extrapolated spectrum to obtain the missing spectrum;
基于所述滤波后的频谱以及所述缺失频谱进行拼接处理,获得所述无缺失的目标频谱。A splicing process is performed based on the filtered spectrum and the missing spectrum to obtain the target spectrum without missing anything.
可选的,所述预测模型包括AR线性预测模型;所述装置还包括用于构建所述AR线性预测模型的构建模块,所述构建模块具体用于:Optionally, the prediction model includes an AR linear prediction model; the device further includes a construction module for constructing the AR linear prediction model, and the construction module is specifically used to:
基于距离向采样点的个数确定AR线性预测模型的阶数;Determine the order of the AR linear prediction model based on the number of distance sampling points;
基于预设的计算方法计算所述AR线性预测模型的系数;Calculating the coefficients of the AR linear prediction model based on a preset calculation method;
基于所述阶数和系数构建所述AR线性预测模型。The AR linear prediction model is constructed based on the order and coefficients.
本申请的实施例提供一种成像方法,包括如下步骤:The embodiment of the present application provides an imaging method, comprising the following steps:
基于上述任一项所述的SAR回波信号的射频干扰抑制方法获得的目标回波信号,获得相应的二维SAR数据;Obtain corresponding two-dimensional SAR data based on a target echo signal obtained by the radio frequency interference suppression method for SAR echo signals described in any one of the above items;
利用所述二维SAR数据对目标对象进行成像处理。The two-dimensional SAR data is used to perform imaging processing on the target object.
本发明通过对回波信号进行陷波处理以及对陷波后频谱进行频谱外推处理,能够准确的对射频干扰进行抑制,解决了现有技术中射频干扰抑制不够准确的问题,同时本申请中的方法克服了传统陷波法的细节丢失的技术缺陷,有利于后续基于SAR数据进行成像。The present invention can accurately suppress radio frequency interference by performing notch processing on the echo signal and performing spectrum extrapolation processing on the spectrum after notching, thereby solving the problem of inaccurate radio frequency interference suppression in the prior art. At the same time, the method in the present application overcomes the technical defect of detail loss in the traditional notch method, which is beneficial to subsequent imaging based on SAR data.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1为本发明一实施例SAR回波信号的射频干扰抑制方法的流程图;FIG1 is a flow chart of a method for suppressing radio frequency interference of a SAR echo signal according to an embodiment of the present invention;
图2为本发明又一实施例SAR回波信号的射频干扰抑制装置的结构框图;FIG2 is a structural block diagram of a device for suppressing radio frequency interference of SAR echo signals according to another embodiment of the present invention;
图3为本发明另一实施一种成像方法的方法流程图。FIG. 3 is a flow chart of another embodiment of an imaging method according to the present invention.
具体实施方式DETAILED DESCRIPTION
此处参考附图描述本申请的各种方案以及特征。Various aspects and features of the present application are described herein with reference to the accompanying drawings.
应理解的是,可以对此处申请的实施例做出各种修改。因此,上述说明书不应该视为限制,而仅是作为实施例的范例。本领域的技术人员将想到在本申请的范围和精神内的其他修改。It should be understood that various modifications may be made to the embodiments of the present application. Therefore, the above description should not be considered as limiting, but only as an example of the embodiments. Other modifications within the scope and spirit of the present application will occur to those skilled in the art.
包含在说明书中并构成说明书的一部分的附图示出了本申请的实施例,并且与上面给出的对本申请的大致描述以及下面给出的对实施例的详细描述一起用于解释本申请的原理。The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate embodiments of the present application and, together with the general description of the present application given above and the detailed description of the embodiments given below, serve to explain the principles of the present application.
通过下面参照附图对给定为非限制性实例的实施例的优选形式的描述,本申请的这些和其它特性将会变得显而易见。These and other characteristics of the present application will become apparent from the following description of a preferred form of embodiment given as a non-limiting example with reference to the accompanying drawings.
还应当理解,尽管已经参照一些具体实例对本申请进行了描述,但本领域技术人员能够确定地实现本申请的很多其它等效形式,它们具有如权利要求所述的特征并因此都位于借此所限定的保护范围内。It should also be understood that, although the present application has been described with reference to some specific examples, those skilled in the art will be able to implement many other equivalent forms of the present application that have the features described in the claims and are therefore within the scope of protection defined thereby.
当结合附图时,鉴于以下详细说明,本申请的上述和其他方面、特征和优势将变得更为显而易见。The above and other aspects, features and advantages of the present application will become more apparent in view of the following detailed description when taken in conjunction with the accompanying drawings.
此后参照附图描述本申请的具体实施例;然而,应当理解,所申请的实施例仅仅是本申请的实例,其可采用多种方式实施。熟知和/或重复的功能和机构并未详细描述以避免不必要或多余的细节使得本申请模糊不清。因此,本文所申请的具体的机构性和功能性细节并非意在限定,而是仅仅作为权利要求的基础和代表性基础用于教导本领域技术人员以实质上任意合适的详细机构多样地使用本申请。Specific embodiments of the present application are described hereinafter with reference to the accompanying drawings; however, it should be understood that the embodiments applied for are merely examples of the present application, which may be implemented in a variety of ways. Well-known and/or repeated functions and mechanisms are not described in detail to avoid unnecessary or redundant details that obscure the present application. Therefore, the specific structural and functional details applied for herein are not intended to be limiting, but merely serve as the basis and representative basis for the claims to teach those skilled in the art to use the present application in a variety of ways with substantially any suitable detailed mechanism.
本说明书可使用词组“在一种实施例中”、“在另一个实施例中”、“在又一实施例中”或“在其他实施例中”,其均可指代根据本申请的相同或不同实施例中的一个或多个。This specification may use the phrases "in one embodiment," "in another embodiment," "in yet another embodiment," or "in other embodiments," all of which may refer to one or more of the same or different embodiments according to the present application.
本发明实施例提供一种SAR回波信号的射频干扰抑制方法,包括如下步骤:The embodiment of the present invention provides a method for suppressing radio frequency interference of a SAR echo signal, comprising the following steps:
步骤S101,对接收到的初始回波信号进行距离向傅里叶变换,获得与所述初始回波信号对应的初始频谱;Step S101, performing a range Fourier transform on the received initial echo signal to obtain an initial frequency spectrum corresponding to the initial echo signal;
本步骤在具体实施过程中是要将接收的初始回波信号即SAR回波信号x(n)进行距离向傅里叶变换,得到方位时域距离频域的回波信号X(n),即初始频谱,也可以称为回波距离线频谱,其中-Nr/2≤n≤Nr/2-1,Nr表示距离向采样点数。In the specific implementation process, this step is to perform a range Fourier transform on the received initial echo signal, i.e., the SAR echo signal x(n), to obtain an echo signal X(n) in the azimuth time domain and range frequency domain, i.e., the initial spectrum, which can also be called the echo range line spectrum, where -Nr / 2≤n≤Nr /2-1, and Nr represents the number of sampling points in the range direction.
步骤S102,基于预设的陷波干扰检测门限值对所述初始频谱进行陷波处理,获得陷波处理后的频谱;Step S102, performing notch processing on the initial spectrum based on a preset notch interference detection threshold value to obtain a spectrum after notch processing;
本步骤在具体实施过程中是需要先基于预设的第一陷波干扰检测门限值对所述初始频谱进行陷波处理,获得第一频谱;然后再基于预设的第二陷波干扰检测门限值对所述第一频谱进行陷波处理,获得所述陷波处理后的频谱。即对初始频谱进行两次陷波处理。In the specific implementation process of this step, it is necessary to first perform notch processing on the initial spectrum based on the preset first notch interference detection threshold value to obtain the first spectrum; and then perform notch processing on the first spectrum based on the preset second notch interference detection threshold value to obtain the notch-processed spectrum. That is, the initial spectrum is notched twice.
更加具体的,首先对初始频谱进行第一次门限陷波,具体步骤如下:计算第一次陷波干扰检测门限值λ1,计算公式为:More specifically, the initial spectrum is firstly notched by the first threshold, and the specific steps are as follows: Calculate the first notch interference detection threshold value λ 1 , and the calculation formula is:
λ1=μ1+σ1 (1)λ 1 =μ 1 +σ 1 (1)
其中,μ1表示回波距离线频谱的均值;X(n)表示回波距离线信号的频谱即初始频谱;σ1表示回波距离线频谱的标准差;Nr表示距离向的采样点数;N0表示过渡带点数。Wherein, μ 1 represents the mean of the echo range line spectrum; X(n) represents the spectrum of the echo range line signal, i.e., the initial spectrum; σ 1 represents the standard deviation of the echo range line spectrum; N r represents the number of sampling points in the range direction; and N 0 represents the number of transition zone points.
在利用上述计算公式计算获得λ1后,就可以基于λ1对初始频谱X(n)进行第一次陷波,即将满足|X(n)|>λ1的信号采样点置零,得到一次陷波后回波信号的第一频谱X1(n):After obtaining λ 1 by using the above calculation formula, the initial spectrum X(n) can be notched for the first time based on λ 1 , that is, the signal sampling points satisfying |X(n)|>λ 1 are set to zero, and the first spectrum X 1 (n) of the echo signal after the first notch is obtained:
其中,X(n)表示回波距离线的频谱,λ1表示第一次陷波的门限值,|·|表示信号的绝对值。Where X(n) represents the frequency spectrum of the echo range line, λ 1 represents the threshold value of the first notch, and |·| represents the absolute value of the signal.
在获得了第一频谱X1(n)后,就可以对回波信号进行第二次门限陷波,具体步骤如下:首先计算第二陷波干扰检测门限值λ2,计算公式为:After obtaining the first spectrum X 1 (n), the echo signal can be subjected to a second threshold notch. The specific steps are as follows: First, the second notch interference detection threshold λ 2 is calculated. The calculation formula is:
λ2=μ2+σ2 (5)λ 2 =μ 2 +σ 2 (5)
其中,μ2表示第一次陷波后非零频谱的均值;X1(n)表示回波信号一次陷波后的频谱即第一频谱;σ2表示第一次陷波后非零频谱的标准差;N1为第一次陷波后置零的干扰频点数;Nr表示距离向的采样点数;N0表示过渡带点数。Among them, μ 2 represents the mean of the non-zero spectrum after the first notch; X 1 (n) represents the spectrum of the echo signal after one notch, that is, the first spectrum; σ 2 represents the standard deviation of the non-zero spectrum after the first notch; N 1 is the number of interference frequency points set to zero after the first notch; N r represents the number of sampling points in the range direction; N 0 represents the number of transition band points.
对回波信号频谱进行第二次陷波,即将满足|X1(n)|>λ2的信号采样点置零,得到二次陷波后回波信号的频谱X2(n):Perform a second notch on the echo signal spectrum, that is, set the signal sampling points satisfying |X 1 (n)|>λ 2 to zero, and obtain the spectrum X 2 (n) of the echo signal after the second notch:
其中,X1(n)表示回波信号一次陷波后的距离线频谱,X(n)表示回波距离线的频谱,λ2表示第二次陷波的门限值,|·|表示信号的绝对值。Wherein, X 1 (n) represents the range line spectrum of the echo signal after the first notch, X(n) represents the spectrum of the echo range line, λ 2 represents the threshold value of the second notch, and |·| represents the absolute value of the signal.
在计算获得第二陷波干扰检测门限值λ2后,就可以基于λ2对第一频谱X1(n)进行第二次陷波,即将满足|X1(n)|>λ2的信号采样点置零,得到二次陷波处理后的回波信号的频谱X2(n)。After the second notch interference detection threshold λ 2 is calculated, the first spectrum X 1 (n) can be notched for the second time based on λ 2 , that is, signal sampling points satisfying |X 1 (n)|>λ 2 are set to zero to obtain the spectrum X 2 (n) of the echo signal after the second notch processing.
步骤S103,对所述陷波处理后的频谱进行滤波处理,获得滤波后的频谱;Step S103, filtering the spectrum after the notch processing to obtain a filtered spectrum;
本步骤在具体实施过程中,在获得了陷波处理后的频谱X2(n)后,就可以对X2(n)进行距离向匹配滤波,得到匹配滤波后的频谱X3(n),公式为:In the specific implementation process of this step, after obtaining the spectrum X 2 (n) after notch processing, X 2 (n) can be subjected to range matched filtering to obtain the spectrum X 3 (n) after matched filtering, and the formula is:
X3(n)=X2(n)·H(f) (9)X 3 (n) = X 2 (n)·H (f) (9)
其中,匹配滤波器H(f)为:Among them, the matched filter H(f) is:
其中,rect(·)表示矩形窗函数,X2(n)为二次陷波后回波信号的频谱,f为距离向频率,Kr为距离向调频率,Br为信号带宽。Wherein, rect(·) represents the rectangular window function, X 2 (n) is the frequency spectrum of the echo signal after the secondary notch, f is the range frequency, K r is the range modulation frequency, and Br is the signal bandwidth.
步骤S104,基于预设的预测模型对所述滤波处理后的频谱进行频谱外推处理,获得前向外推频谱数据以及后向外推频谱数据;Step S104, performing spectrum extrapolation processing on the filtered spectrum based on a preset prediction model to obtain forward extrapolated spectrum data and backward extrapolated spectrum data;
本步骤中的预测模块包括AR线性预测模型。AR线性预测模型是预先构建获得,构建步骤包括:基于距离向采样点的个数确定AR线性预测模型的阶数;基于预设的计算方法计算所述AR线性预测模型的系数;基于所述阶数和系数构建所述AR线性预测模型。The prediction module in this step includes an AR linear prediction model. The AR linear prediction model is pre-constructed, and the construction steps include: determining the order of the AR linear prediction model based on the number of distance sampling points; calculating the coefficients of the AR linear prediction model based on a preset calculation method; and constructing the AR linear prediction model based on the order and coefficients.
更加具体的,本实施例中在进行频谱外推处理之前,需要先确定AR线性预测模型阶数M和系数a(m),以此来构建AR线性预测模型。即基于距离向的采样点的个数来确定AR线性预测模型阶数,距离向的采样点的个数为Nr,则AR线性预测模型阶数M一般取Nr/3~Nr/2之间的值。AR线性预测模型的系数a(m)可以用Burg(伯格)算法、改进的协方差等方法计算。More specifically, in this embodiment, before performing spectrum extrapolation processing, it is necessary to first determine the AR linear prediction model order M and coefficient a(m) to construct the AR linear prediction model. That is, the AR linear prediction model order is determined based on the number of sampling points in the distance direction. The number of sampling points in the distance direction is N r , and the AR linear prediction model order M generally takes a value between N r /3 and N r /2. The coefficient a(m) of the AR linear prediction model can be calculated using methods such as the Burg algorithm and the improved covariance.
在确定了AR线性预测模型阶数M和系数a(m)之后,就可以基于滤波后的频谱X3(n)利用该预测模块来外推缺失频谱,分别得到前向外推频谱X4(n)和后向外推频谱X5(n)。前向外推频谱X4(n)和后向外推频谱X5(n)具体是利用AR线性预测模的两个公式来计算获得的,公式如下:After determining the AR linear prediction model order M and coefficient a(m), the prediction module can be used to extrapolate the missing spectrum based on the filtered spectrum X 3 (n) to obtain the forward extrapolated spectrum X 4 (n) and the backward extrapolated spectrum X 5 (n). The forward extrapolated spectrum X 4 (n) and the backward extrapolated spectrum X 5 (n) are specifically calculated using two formulas of the AR linear prediction module, as follows:
其中,X3(n)表示距离向匹配滤波后的频谱;M表示预测模型阶数;a(m)表示预测模型系数;a*(m)表示预测模型系数的共轭复数;fLow,l和fHig,l分别表示第l个缺失频谱间隙的最低频率和最高频率;Δf表示频率采样间隔。Wherein, X 3 (n) represents the spectrum after range matched filtering; M represents the prediction model order; a(m) represents the prediction model coefficient; a * (m) represents the conjugate complex number of the prediction model coefficient; f Low,l and f Hig,l represent the lowest frequency and the highest frequency of the lth missing spectrum gap respectively; Δf represents the frequency sampling interval.
步骤S105,基于所述前向外推频谱数据、所述后向外推频谱数据以及所述滤波后的频谱进行重建处理,获得无缺失的目标频谱;Step S105, performing reconstruction processing based on the forward extrapolated spectrum data, the backward extrapolated spectrum data and the filtered spectrum to obtain a target spectrum without missing components;
本步骤在具体实施过程中,需要先基于所述前向外推频谱以及所述后向外推频谱进行缺失频谱重建,获得缺失频谱;然后再基于所述滤波后的频谱以及所述缺失频谱进行拼接处理,获得所述无缺失的目标频谱。In the specific implementation process of this step, it is necessary to first reconstruct the missing spectrum based on the forward extrapolated spectrum and the backward extrapolated spectrum to obtain the missing spectrum; and then perform splicing processing based on the filtered spectrum and the missing spectrum to obtain the target spectrum without missing spectrum.
更加具体的,在获得前向外推频谱X4(n)和后向外推频谱X5(n)后,就可以将得到的前向外推频谱数据X4(n)和后向外推频谱数据X5(n)按照如下计算公式进行加权,得到缺失数据重建的最终结果计算公式为:More specifically, after obtaining the forward extrapolated spectrum X 4 (n) and the backward extrapolated spectrum X 5 (n), the forward extrapolated spectrum data X 4 (n) and the backward extrapolated spectrum data X 5 (n) can be weighted according to the following calculation formula to obtain the final result of missing data reconstruction: The calculation formula is:
其中,表示缺失频谱;NLow,l和NHig,l分别表示第l个缺失频谱间隙的最低频率和最高频率对应的频点;X3(n)表示距离向匹配滤波后的频谱;M表示预测模型阶数;a(m)表示预测模型系数;a*(m)表示预测模型系数的共轭复数;fLow,l和fHig,l分别表示第l个缺失频谱间隙的最低和最高频率;Δf表示频率采样间隔。in, represents the missing spectrum; N Low,l and N Hig,l represent the frequency points corresponding to the lowest frequency and the highest frequency of the lth missing spectrum gap respectively; X 3 (n) represents the spectrum after distance matching filtering; M represents the prediction model order; a(m) represents the prediction model coefficient; a * (m) represents the conjugate complex number of the prediction model coefficient; f Low,l and f Hig,l represent the lowest and highest frequencies of the lth missing spectrum gap respectively; Δf represents the frequency sampling interval.
在获得缺失频谱后,就可以基于该获得缺失频谱以及滤波后的频谱X3(n)利用如下公式来进行拼接处理,获得无缺失的目标频谱X6(n),公式为:In obtaining the missing spectrum Then, the missing spectrum can be obtained based on The filtered spectrum X 3 (n) is concatenated using the following formula to obtain a target spectrum X 6 (n) without missing components. The formula is:
其中,H*(f)为H(f)的复共轭,X6(n)为无缺失的匹配滤波后的信号频谱。Wherein, H * (f) is the complex conjugate of H(f), and X 6 (n) is the signal spectrum after matched filtering without missing parts.
步骤S106,基于所述无缺失的目标频谱进行距离向信号的调频信息的恢复,获得抑制干扰且频谱完整的目标回波信号。Step S106, recovering the frequency modulation information of the range signal based on the target spectrum without missing anything, to obtain a target echo signal with suppressed interference and complete spectrum.
本步骤在在获得了目标频谱X6(n)后,就可以通过如下公式来恢复SAR数据距离向信号的调频信息,得到抑制干扰且频谱完整的目标回波信号X7(n),公式为:After obtaining the target spectrum X 6 (n) in this step, the frequency modulation information of the SAR data range signal can be restored by the following formula to obtain the target echo signal X 7 (n) with suppressed interference and complete spectrum. The formula is:
X7(n)=X6(n)·H*(f) (15)X 7 (n)=X 6 (n)·H * (f) (15)
其中,H*(f)为H(f)的复共轭,X6(n)为无缺失的匹配滤波后的信号频谱即目标频谱。Wherein, H * (f) is the complex conjugate of H(f), and X 6 (n) is the signal spectrum after matched filtering without missing, that is, the target spectrum.
本发明实施通过在二次陷波后添加了频谱外推,使陷波后的频谱重新恢复,能够准确的对射频干扰进行抑制,解决了现有技术中射频干扰抑制不够准确的问题。本发明中的方法不仅具有传统方法的优点,而且弥补了传统方法的不足,有利于改善SAR图像的质量,提高了图像中细节部分的识别度。与此同时本发明中的方法降低了计算复杂度,提高了运算效率。The present invention implements the method of adding spectrum extrapolation after the secondary notch to restore the spectrum after the notch, which can accurately suppress radio frequency interference and solve the problem of inaccurate radio frequency interference suppression in the prior art. The method of the present invention not only has the advantages of the traditional method, but also makes up for the shortcomings of the traditional method, which is conducive to improving the quality of SAR images and improving the recognition of details in the image. At the same time, the method of the present invention reduces the computational complexity and improves the computational efficiency.
本发明另一实施例提供一种SAR回波信号的射频干扰抑制装置,包括:Another embodiment of the present invention provides a SAR echo signal radio frequency interference suppression device, comprising:
变换模块1,用于对接收到的初始回波信号进行距离向傅里叶变换,获得与所述初始回波信号对应的初始频谱;Transformation module 1, used for performing range Fourier transform on the received initial echo signal to obtain an initial spectrum corresponding to the initial echo signal;
陷波处理模块2,用于基于预设的陷波干扰检测门限值对所述初始频谱进行陷波处理,获得陷波处理后的频谱;A notch processing module 2, configured to perform notch processing on the initial spectrum based on a preset notch interference detection threshold value to obtain a spectrum after notch processing;
滤波处理模块3,用于对所述陷波处理后的频谱进行滤波处理,获得滤波后的频谱;A filtering processing module 3 is used to filter the spectrum after the notch processing to obtain a filtered spectrum;
频谱外推处理模块4,用于基于预设的预测模型对所述滤波处理后的频谱进行频谱外推处理,获得前向外推频谱数据以及后向外推频谱数据;The spectrum extrapolation processing module 4 is used to perform spectrum extrapolation processing on the spectrum after filtering based on a preset prediction model to obtain forward extrapolated spectrum data and backward extrapolated spectrum data;
重建模块5,用于基于所述前向外推频谱数据、所述后向外推频谱数据以及所述滤波后的频谱进行重建处理,获得无缺失的目标频谱;A reconstruction module 5, configured to perform reconstruction processing based on the forward extrapolated spectrum data, the backward extrapolated spectrum data and the filtered spectrum to obtain a target spectrum without missing components;
恢复模块6,用于基于所述无缺失的目标频谱进行距离向信号的调频信息的恢复,获得抑制干扰且频谱完整的目标回波信号。The recovery module 6 is used to recover the frequency modulation information of the range signal based on the target spectrum without missing any, so as to obtain the target echo signal with suppressed interference and complete spectrum.
本实施例在具体实施过程中,所述陷波处理模块具体用于:基于预设的第一陷波干扰检测门限值对所述初始频谱进行陷波处理,获得第一频谱;基于预设的第二陷波干扰检测门限值对所述第一频谱进行陷波处理,获得所述陷波处理后的频谱。In the specific implementation process of this embodiment, the notch processing module is specifically used to: perform notch processing on the initial spectrum based on a preset first notch interference detection threshold value to obtain a first spectrum; perform notch processing on the first spectrum based on a preset second notch interference detection threshold value to obtain the notch-processed spectrum.
本实施例在具体实施过程中,所述重建模块具体用于:基于所述前向外推频谱以及所述后向外推频谱进行缺失频谱重建,获得缺失频谱;In the specific implementation process of this embodiment, the reconstruction module is specifically used to: reconstruct the missing spectrum based on the forward extrapolated spectrum and the backward extrapolated spectrum to obtain the missing spectrum;
基于所述滤波后的频谱以及所述缺失频谱进行拼接处理,获得所述无缺失的目标频谱。A splicing process is performed based on the filtered spectrum and the missing spectrum to obtain the target spectrum without missing anything.
本实施例在具体实施过程中,所述预测模型包括AR线性预测模型;所述装置还包括用于构建所述AR线性预测模型的构建模块,所述构建模块具体用于:基于距离向采样点的个数确定AR线性预测模型的阶数;基于预设的计算方法计算所述AR线性预测模型的系数;基于所述阶数和系数构建所述AR线性预测模型。本实施例中预设的计算方法包括:伯格算法和改进的协方差算法。In the specific implementation process of this embodiment, the prediction model includes an AR linear prediction model; the device also includes a construction module for constructing the AR linear prediction model, and the construction module is specifically used to: determine the order of the AR linear prediction model based on the number of distance sampling points; calculate the coefficients of the AR linear prediction model based on a preset calculation method; and construct the AR linear prediction model based on the order and coefficients. The preset calculation methods in this embodiment include: the Berg algorithm and the improved covariance algorithm.
本发明实施通过利用陷波处理模块来对初始频谱进行二次陷波后处理利用频谱外推模块来对陷波后的频谱重新恢复,能够准确的对射频干扰进行抑制。本发明中的方法不仅具有传统方法的优点,而且弥补了传统方法的不足,有利于改善SAR图像的质量,提高了图像中细节部分的识别度。与此同时本发明中的方法降低了计算复杂度,提高了运算效率。The present invention implements the method of performing secondary notching post-processing on the initial spectrum by using a notch processing module and restoring the notched spectrum by using a spectrum extrapolation module, so as to accurately suppress radio frequency interference. The method of the present invention not only has the advantages of the traditional method, but also makes up for the shortcomings of the traditional method, is conducive to improving the quality of SAR images, and improves the recognition of details in the image. At the same time, the method of the present invention reduces the computational complexity and improves the computational efficiency.
本发明另一实施例提一种成像方法,包括:基于上述的SAR回波信号的射频干扰抑制方法获得的目标回波信号,获得相应的二维SAR数据;利用所述二维SAR数据对目标对象进行成像处理。具体包括如下步骤:Another embodiment of the present invention provides an imaging method, comprising: obtaining corresponding two-dimensional SAR data based on the target echo signal obtained by the above-mentioned SAR echo signal radio frequency interference suppression method; and performing imaging processing on the target object using the two-dimensional SAR data. Specifically, the method comprises the following steps:
步骤S301,基于采样时间对目标对象的各方位向进行回波数据采样,获得与各方位向对应的初始回波信号;Step S301, sampling echo data in each direction of the target object based on the sampling time to obtain initial echo signals corresponding to each direction;
步骤S302,对各所述初始回波信号进行距离向傅里叶变换,获得与各所述初始回波信号对应的初始频谱;Step S302, performing range Fourier transform on each of the initial echo signals to obtain an initial frequency spectrum corresponding to each of the initial echo signals;
步骤S303,基于预设的陷波干扰检测门限值对各所述初始频谱进行陷波处理,获得与各所述初始频谱对应的陷波处理后的频谱;Step S303, performing notch processing on each of the initial spectra based on a preset notch interference detection threshold value to obtain a notch-processed spectrum corresponding to each of the initial spectra;
步骤S304,对各所述陷波处理后的频谱进行滤波处理,获得与各所述陷波处理后的频谱对应的滤波后的频谱;Step S304, filtering each of the frequency spectra after notch processing to obtain a filtered frequency spectrum corresponding to each of the frequency spectra after notch processing;
步骤S305,基于预设的预测模型对各所述滤波处理后的频谱进行频谱外推处理,获得各所述滤波处理后的频谱对应的前向外推频谱数据以及后向外推频谱数据;Step S305, performing spectrum extrapolation processing on each of the filtered spectrums based on a preset prediction model to obtain forward extrapolated spectrum data and backward extrapolated spectrum data corresponding to each of the filtered spectrums;
步骤S306,基于各所述滤波后的频谱、与所述滤波后的频谱对应的前向外推频谱数据以及后向外推频谱数据进行重建处理,获得与各所述初始回波信号对应的无缺失的目标频谱;Step S306, performing reconstruction processing based on each of the filtered spectra, the forward extrapolated spectrum data corresponding to the filtered spectrum, and the backward extrapolated spectrum data to obtain a target spectrum without missing parts corresponding to each of the initial echo signals;
步骤S307,基于各所述无缺失的目标频谱进行距离向信号的调频信息的恢复,获得与各所述初始回波信号对应的目标回波信号。Step S307 , recovering the frequency modulation information of the range signal based on each of the target frequency spectra without missing items, and obtaining the target echo signal corresponding to each of the initial echo signals.
步骤S308,基于各所述目标回波信号获得与目标对象对应的二维SAR数据;Step S308, obtaining two-dimensional SAR data corresponding to the target object based on each of the target echo signals;
步骤S309,基于所述二维SAR数据进行成像处理,获得所述目标对象的SAR图像。Step S309: performing imaging processing based on the two-dimensional SAR data to obtain a SAR image of the target object.
本发明实施例在实施过程中,步骤S302-步骤S307的具体实施过程可以包括上述实施例提供的SAR回波信号的射频干扰抑制方法中的步骤,在次不再赘述。During the implementation of the embodiment of the present invention, the specific implementation process of step S302 to step S307 may include the steps in the method for suppressing radio frequency interference of SAR echo signals provided in the above embodiment, which will not be described in detail here.
本发明实施通过对二次陷波后的频谱进行频谱外推处理,能够使陷波后的频谱重新恢复,准确的对射频干扰进行了抑制,有利于改善SAR图像的质量,提高了图像中细节部分的识别度。与此同时本发明中的方法降低了计算复杂度,提高了运算效率。可以在更短时间内得到更好成像效果,即同时提高了成像效果和运算效率。The present invention implements spectrum extrapolation processing on the spectrum after the secondary notch, so that the spectrum after the notch can be restored, and the radio frequency interference is accurately suppressed, which is beneficial to improving the quality of the SAR image and improving the recognition of the details in the image. At the same time, the method of the present invention reduces the computational complexity and improves the operation efficiency. A better imaging effect can be obtained in a shorter time, that is, the imaging effect and the operation efficiency are improved at the same time.
以上实施例仅为本发明的示例性实施例,不用于限制本发明,本发明的保护范围由权利要求书限定。本领域技术人员可以在本发明的实质和保护范围内,对本发明做出各种修改或等同替换,这种修改或等同替换也应视为落在本发明的保护范围内。The above embodiments are only exemplary embodiments of the present invention and are not intended to limit the present invention. The protection scope of the present invention is defined by the claims. Those skilled in the art may make various modifications or equivalent substitutions to the present invention within the essence and protection scope of the present invention, and such modifications or equivalent substitutions shall also be deemed to fall within the protection scope of the present invention.
Claims (8)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010640070.9A CN111580107B (en) | 2020-07-06 | 2020-07-06 | Radio frequency interference suppression method and device for SAR echo signal and imaging method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010640070.9A CN111580107B (en) | 2020-07-06 | 2020-07-06 | Radio frequency interference suppression method and device for SAR echo signal and imaging method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111580107A CN111580107A (en) | 2020-08-25 |
CN111580107B true CN111580107B (en) | 2023-03-31 |
Family
ID=72127530
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010640070.9A Active CN111580107B (en) | 2020-07-06 | 2020-07-06 | Radio frequency interference suppression method and device for SAR echo signal and imaging method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111580107B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111913162B (en) * | 2020-09-04 | 2023-02-24 | 内蒙古工业大学 | Radio frequency interference detection method and device and radio frequency interference suppression method and device |
CN113064122B (en) * | 2021-02-25 | 2023-02-17 | 上海卫星工程研究所 | Performance evaluation method, system and medium of P-band SAR (synthetic Aperture Radar) interference suppression algorithm |
CN114200410B (en) * | 2021-11-12 | 2025-02-07 | 中山大学 | Radio Frequency Interference Suppression Method for Low-Frequency Ultra-Wideband SAR Based on Two-Dimensional Frequency Domain Notching |
CN116299465B (en) * | 2023-05-18 | 2023-08-01 | 西安电子科技大学 | A Backprojection Imaging Method for Bistatic SAR Based on Subspace Time-Frequency Mapping |
CN118625282B (en) * | 2024-08-12 | 2024-10-15 | 中国科学院空天信息创新研究院 | Distributed spaceborne SAR missing data recovery method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103176184A (en) * | 2013-02-06 | 2013-06-26 | 中国科学院电子学研究所 | P band SAR (synthetic aperture radar) imaging processing method combining interference suppression |
CN105974376A (en) * | 2016-05-09 | 2016-09-28 | 中国人民解放军国防科学技术大学 | SAR radio frequency interference suppressing method |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2342402C (en) * | 1999-02-02 | 2003-09-02 | Illinois Superconductor Canada Corporation | A method for detecting and eliminating narrowband channel interference |
SE0301231D0 (en) * | 2003-04-25 | 2003-04-25 | Ericsson Telefon Ab L M | SAR interference suppression |
US8385483B2 (en) * | 2008-11-11 | 2013-02-26 | Isco International, Llc | Self-adaptive digital RF bandpass and bandstop filter architecture |
CN102590829B (en) * | 2012-03-14 | 2013-09-25 | 西安电子科技大学 | Complete self-adaptive notch filter for satellite navigation system and notch filtering method of same |
CN103383454A (en) * | 2012-05-03 | 2013-11-06 | 中国科学院电子学研究所 | Method for restraining analog television signal in P waveband SAR echo signal |
CN103226325A (en) * | 2013-04-16 | 2013-07-31 | 深圳市汇川技术股份有限公司 | Method and device for eliminating resonance of servomotor |
CN106597443B (en) * | 2016-12-20 | 2019-02-26 | 中国人民解放军国防科学技术大学 | Low-frequency synthetic aperture radar image contrast destructive radio frequency interference suppression method |
-
2020
- 2020-07-06 CN CN202010640070.9A patent/CN111580107B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103176184A (en) * | 2013-02-06 | 2013-06-26 | 中国科学院电子学研究所 | P band SAR (synthetic aperture radar) imaging processing method combining interference suppression |
CN105974376A (en) * | 2016-05-09 | 2016-09-28 | 中国人民解放军国防科学技术大学 | SAR radio frequency interference suppressing method |
Also Published As
Publication number | Publication date |
---|---|
CN111580107A (en) | 2020-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111580107B (en) | Radio frequency interference suppression method and device for SAR echo signal and imaging method | |
CN111142105B (en) | ISAR imaging method for complex moving target | |
CN112269168B (en) | SAR broadband interference suppression method based on Bayesian theory and low-rank decomposition | |
EP3229038A1 (en) | Wavelet domain insar interferometric phase filtering method in combination with local frequency estimation | |
CN104133198B (en) | Ionospheric interference suppressing method in a kind of high-frequency ground wave radar | |
CN109324315B (en) | Space-time adaptive radar clutter suppression method based on double-layer block sparsity | |
CN105974376B (en) | A kind of SAR radio frequency interferences suppressing method | |
CN107942327A (en) | Single channel HRWS SAR imaging methods based on impulse phase coding | |
CN110488285B (en) | Non-interrupt bistatic SAR (synthetic aperture radar) phase synchronization signal processing method based on coding signal | |
CN103675768B (en) | A kind of suppression method for satellite-borne SAR radio-frequency interference | |
CN109143231B (en) | Digital television passive bistatic radar target detection method based on cyclic cancellation | |
CN112327259B (en) | Method and device for eliminating interference signals in SAR image | |
CN115494470A (en) | Geological imaging radar integrated signal processing method, system, terminal and medium | |
CN113238193A (en) | Multi-component combined reconstruction SAR echo broadband interference suppression method | |
CN113176543A (en) | Radio frequency interference suppression method and system based on multi-dimensional information combination | |
CN116148802B (en) | A Phase Synchronous Anti-jamming Method for Bistatic SAR | |
CN113848532A (en) | FMCW radar signal noise reduction system and method based on noise reduction model | |
CN106842204A (en) | A kind of high-resolution Spotlight SAR Imaging two dimension missing data image processing method | |
CN112882026B (en) | Tree height inversion method and device, computer storage medium and terminal | |
CN108008378A (en) | A kind of SAR false targets energy ratio determines method and device | |
Dutta et al. | SAR image segmentation using wavelets and Gaussian mixture model | |
CN115097391A (en) | Synthetic aperture radar interference suppression method based on automatic encoder | |
CN118759487B (en) | A method for suppressing SAR pulsed radio frequency interference | |
Zhang et al. | An Interference Suppression Method For Spaceborne Sar Image Via Space-Channel Attention Network | |
CN116973849A (en) | SAR image radio frequency interference suppression method based on spectrum outlier detection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |