CN105974376A - 一种sar射频干扰抑制方法 - Google Patents
一种sar射频干扰抑制方法 Download PDFInfo
- Publication number
- CN105974376A CN105974376A CN201610298996.8A CN201610298996A CN105974376A CN 105974376 A CN105974376 A CN 105974376A CN 201610298996 A CN201610298996 A CN 201610298996A CN 105974376 A CN105974376 A CN 105974376A
- Authority
- CN
- China
- Prior art keywords
- signal
- echo
- reference signal
- sar
- interference
- 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.)
- Granted
Links
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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种SAR射频干扰抑制方法。技术方案是:以SAR数据的一条距离线回波信号为输入,根据回波信号构建相空间矩阵,利用子空间投影方法得到参考信号的相空间矩阵;采用对角平均化方法,从参考信号相空间矩阵中恢复出一维的参考信号;根据参考信号,确定干扰信号的判别门限;根据确定的干扰判别门限,确定回波中存在的干扰频点,并用参考信号的幅度均值作为替换,完成该条距离线的射频干扰抑制。本发明可以有效提升干扰抑制后的SAR图像质量。
Description
技术领域
本发明属于航天和微波遥感的交叉技术领域,特别涉及一种SAR(SyntheticAperture Radar,合成孔径雷达)数据中RFI(Radio Frequency Interference,射频干扰)抑制方法。
背景技术
SAR(Synthetic Aperture Radar,合成孔径雷达)通过主动发射电磁波来获取目标信息,是一种全天候、全天时的高分辨率成像系统,在目标侦察和地球遥感应用中具有巨大的潜力。然而,SAR(尤其是低波段SAR)在复杂的电磁环境中容易受到RFI的影响,进而降低图像质量,影响图像识别和应用。
SAR系统由于其良好的穿透特性,在隐蔽目标侦察以及森林植被测绘方面有着独特的优势,但其工作频段内存在大量的广播、电视信号以及地面警戒雷达发射的信号。这些干扰信号会导致SAR图像性能严重恶化。
陷波法和子空间投影法是当前广泛使用的RFI抑制方法,但传统的陷波法容易造成有用信息的损失,而子空间投影法的局限性在于过分依赖于子空间的正交性,当干扰信号功率较大时,会造成子空间误判,在图像中未受干扰区域引入虚警。针对这一现状,提出了一种基于联合滤波的RFI抑制方法,该方法综合了陷波法和子空间投影法的优点,在保证干扰抑制效果的前提下最大限度地保留了有用信号,且性能稳健。
发明内容
本发明的目的是:提出了一种基于联合滤波的RFI抑制方法,该方法是一种非参数方法,充分结合了子空间投影法和陷波法在干扰抑制方面的优点,具有更好的干扰抑制效果,且性能稳健。
本发明技术方案的思路是:首先构建回波的相空间矩阵,利用子空间投影法将回波的相空间矩阵投影到目标信号子空间,获得参考信号,并依据该参考信号,确定干扰信号的判别门限,然后,进行干扰检测、陷波和赋值处理,即可将干扰信号抑制掉。
本发明技术方案是:
假设x(n),1≤n≤Nr表示某一慢时刻SAR接收的回波信号,其中Nr表示距离向采样点数。
第一步:提取参考信号
(1)将雷达接收的回波信号x(n)按照时延嵌套构造成相空间矩阵S。
式中,L的大小取决于射频干扰的个数,且M=Nr+1-L。L过小则特征分解不充分,L过大则计算量过大,并且成像质量也不一定得到提升,通常根据实际情况确定。本发明的具体实施例中L取32。令协方差矩阵R=SSH。
(2)对协方差矩阵R进行特征值分解
R=UΛUH
其中Λ=diag[λ1,λ2,…,λL],λ1≥λ2≥…≥λL是一个L×L的对角矩阵,λi,1≤i≤L为特征值,对应的特征矢量为U=[u1,u2,…,uL]。
将前k个特征值所对应的特征矢量提取出,组成参考子空间Uk=[u1,u2,…,uk],根据实际情况确定,或根据文献“郑慧芳,杨淋,冯锦.SAR窄带干扰抑制的子带子空间滤波技术研究[J].电子与信息学报,2013,35(12):2836-2842”。以表示参考信号相空间矩阵。
(3)采用对角平均化方法从参考信号相空间矩阵Sk中恢复出1维参考信号序列xk(n),1≤n≤Nr。该过程可以看作是嵌套的逆过程,只需要计算Sk的每条与逆对角线平行的副逆对角线的均值,该Nr个均值就构成了1×Nr维参考信号xk(n),1≤n≤Nr。其公式为
第二步,对回波信号进行滤波
将回波信号x(n)和参考信号xk(n)分别进行FFT变换到频域
X=[X(1),X(2),…,X(Nr)]
Xk=[Xk(1),Xk(2),…,Xk(Nr)]
利用下式计算判别门限γ:
利用下式对回波信号进行滤波,得到射频干扰抑制后的回波信号
针对SAR的每个慢时刻的回波信号重复上述操作,即可达到抑制干扰的目的。
采用本发明可取得以下技术效果:
本发明提出的基于联合滤波的SAR射频干扰抑制方法,利用获得的参考信号确定干扰判别门限,并通过滤波去除干扰,根据参考信号对回波信号重新赋值。相比现有的方法,本方法在确定参考信号时利用子空间投影法可以更好地去除干扰的影响,得到更优的干扰判别门限,从而更准确地找到干扰的位置。根据参考信号对回波信号重新赋值时,对干扰频点的赋值更准确,进而实现更好的干扰抑制效果,并且性能稳健。
附图说明
图1为本发明提供的星载SAR射频干扰抑制方法的原理流程图;
图2是仿真实验中使用的点目标的成像结果;
图3是对图2添加射频干扰后点目标的成像结果;
图4为采用陷波法对图3进行干扰抑制后成像结果;
图5为采用子空间投影法对图3进行干扰抑制后成像结果;
图6为采用本文发明提供的方法对图3进行干扰抑制后成像结果;
图7是仿真实验中使用的实测数据的成像结果;
图8是对图7添加射频干扰后数据的成像结果;
图9是利用陷波法对图8进行干扰抑制后的成像结果;
图10是利用子空间投影法对图8进行干扰抑制后的成像结果;
图11是利用本发明提供的方法对图8进行干扰抑制后的成像结果。
具体实施方式
图1为本发明提供的星载SAR射频干扰抑制方法的原理流程图。整个流程包括:以SAR数据的一条距离线回波信号为输入,根据回波信号构建相空间矩阵,利用子空间投影方法得到参考信号的相空间矩阵;采用对角平均化方法,从参考信号相空间矩阵中恢复出一维的参考信号;根据参考信号,确定干扰信号的判别门限;根据确定的干扰判别门限,确定回波中存在的干扰频点,并用参考信号的幅度均值作为替换,完成该条距离线的射频干扰抑制。
图2~图6是进行点目标仿真实验结果。仿真了单个点目标回波,在此回波中添加了带宽为0.5MHz,1.0MHz和1.5MHz的射频干扰信号,干信比为20dB。图2为原始点目标回波成像结果,从图像中可以看出,点目标聚焦效果良好。图3是在图2中添加射频干扰后回波的成像结果,可以看到点目标已经被干扰信号湮没。图4是采用陷波法对图3进行干扰抑制后的成像结果,从图像中可以看出,干扰聚焦良好,但仍有较多的干扰残留。图5是采用子空间投影法对图3进行干扰抑制后的成像结果,点目标不仅不能良好聚焦,图像还残留较多的干扰。图6是采用本发明提供的方法进行干扰抑制后的成像结果,可以看到,点目标聚焦良好,且干扰残留较少。
图7~图11是进行面目标仿真实验的结果。仿真实验选取了一幅无干扰的RadarSat-1回波实测数据,在此数据中添加了带宽为0.5MHz,1.0MHz和1.5MHz的射频干扰信号,干信比为20dB。
图7是仿真实验中使用的原始数据的成像结果。从图中可以看出,图像清晰,对比度高。图8是对图7添加射频干扰后数据的成像结果。从图中可以看出,由于干扰的存在,图像目标信息被覆盖,图像对比度下降。图9和图10分别是利用陷波法和子空间投影法对图8进行干扰抑制后的成像结果,干扰抑制后,图像质量得到提升,不过仍有较多干扰残留,细节信息不明显。图11是利用本发明提供的方法对图8进行干扰抑制后的成像结果,干扰抑制后,图像质量明显提升,细节信息恢复较为明显,图像对比度有明显提高。
以上仿真实验说明了本发明方法可以有效地弥补现有方法在干扰抑制方面的不足,有效提升干扰抑制后的SAR图像质量。
Claims (1)
1.一种SAR射频干扰抑制方法,SAR是指合成孔径雷达,假设x(n),1≤n≤Nr,表示某一慢时刻SAR接收的回波信号,其中Nr表示距离向采样点数,其特征在于,包括下述步骤:
第一步:提取参考信号:
将雷达接收的回波信号x(n)按照时延嵌套构造成相空间矩阵S:
上式中M=Nr+1-L,L根据实际情况确定;
令协方差矩阵R=SSH,对协方差矩阵R进行特征值分解:
R=UΛUH
其中Λ=diag[λ1,λ2,…,λL]是一个L×L的对角矩阵,对应的特征矢量矩阵为U=[u1,u2,…,uL],并且λ1≥λ2≥…≥λL,λi为特征值,i=1,2,…L;上标H表示求矩阵的共轭转置;
将前k个特征值所对应的特征矢量提取出,组成参考子空间Uk=[u1,u2,…,uk],k的取值根据实际情况确定;
令参考信号相空间矩阵采用对角平均化方法从参考信号相空间矩阵Sk中恢复出1维参考信号序列xk(n);
第二步,对回波信号进行滤波
将回波信号x(n)和参考信号xk(n)分别进行快速傅里叶变换得到X和Xk:
X=[X(1),X(2),…,X(Nr)]
Xk=[Xk(1),Xk(2),…,Xk(Nr)]
利用下式计算判别门限γ:
利用下式对回波信号进行滤波,得到射频干扰抑制后的回波信号
针对SAR的每个慢时刻回波信号重复上述两个步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610298996.8A CN105974376B (zh) | 2016-05-09 | 2016-05-09 | 一种sar射频干扰抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610298996.8A CN105974376B (zh) | 2016-05-09 | 2016-05-09 | 一种sar射频干扰抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105974376A true CN105974376A (zh) | 2016-09-28 |
CN105974376B CN105974376B (zh) | 2018-05-25 |
Family
ID=56992145
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610298996.8A Active CN105974376B (zh) | 2016-05-09 | 2016-05-09 | 一种sar射频干扰抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105974376B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106501781A (zh) * | 2016-11-04 | 2017-03-15 | 中国人民解放军国防科学技术大学 | 基于斜投影滤波的sar射频干扰抑制方法 |
CN106597443A (zh) * | 2016-12-20 | 2017-04-26 | 中国人民解放军国防科学技术大学 | 低频合成孔径雷达图像对比相消射频干扰抑制方法 |
CN106680790A (zh) * | 2016-12-29 | 2017-05-17 | 中国电子科技集团公司第五十四研究所 | 一种伪码调相连续波雷达干扰检测方法 |
CN110221256A (zh) * | 2019-06-21 | 2019-09-10 | 西安电子科技大学 | 基于深度残差网络的sar干扰抑制方法 |
CN110531327A (zh) * | 2019-10-11 | 2019-12-03 | 上海无线电设备研究所 | 一种基于自动增益控制的雷达抗干扰探测方法 |
CN110554367A (zh) * | 2018-06-04 | 2019-12-10 | 重庆测威科技有限公司 | 一种基于压缩感知的目标散射特性测量干扰去除方法 |
CN110618428A (zh) * | 2019-09-19 | 2019-12-27 | 上海海洋大学 | 一种用于多普勒三维数据处理方法及软件 |
CN111580107A (zh) * | 2020-07-06 | 2020-08-25 | 内蒙古工业大学 | 一种sar回波信号的射频干扰抑制方法、装置及成像方法 |
CN113064122A (zh) * | 2021-02-25 | 2021-07-02 | 上海卫星工程研究所 | P波段sar干扰抑制算法的性能评价方法、系统及介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090091492A1 (en) * | 2007-10-09 | 2009-04-09 | The Mitre Corporation | Detection and mitigation radio frequency memory (DRFM)-based interference in synthetic aperture radar (SAR) images |
CN102243300A (zh) * | 2011-04-15 | 2011-11-16 | 中国人民解放军国防科学技术大学 | 低频合成孔径雷达射频干扰抑制及误差校正方法 |
CN102608586A (zh) * | 2012-02-11 | 2012-07-25 | 中国人民解放军陆军军官学院 | 一种基于神经网络的合成孔径雷达射频干扰抑制方法 |
CN103091665A (zh) * | 2013-02-01 | 2013-05-08 | 北京航空航天大学 | 一种压缩感知合成孔径雷达射频干扰抑制处理方法 |
CN103176184A (zh) * | 2013-02-06 | 2013-06-26 | 中国科学院电子学研究所 | 一种结合干扰抑制方法的p波段sar成像处理方法 |
CN103675768A (zh) * | 2013-12-17 | 2014-03-26 | 中国人民解放军国防科学技术大学 | 一种星载sar射频干扰抑制方法 |
-
2016
- 2016-05-09 CN CN201610298996.8A patent/CN105974376B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090091492A1 (en) * | 2007-10-09 | 2009-04-09 | The Mitre Corporation | Detection and mitigation radio frequency memory (DRFM)-based interference in synthetic aperture radar (SAR) images |
CN102243300A (zh) * | 2011-04-15 | 2011-11-16 | 中国人民解放军国防科学技术大学 | 低频合成孔径雷达射频干扰抑制及误差校正方法 |
CN102608586A (zh) * | 2012-02-11 | 2012-07-25 | 中国人民解放军陆军军官学院 | 一种基于神经网络的合成孔径雷达射频干扰抑制方法 |
CN103091665A (zh) * | 2013-02-01 | 2013-05-08 | 北京航空航天大学 | 一种压缩感知合成孔径雷达射频干扰抑制处理方法 |
CN103176184A (zh) * | 2013-02-06 | 2013-06-26 | 中国科学院电子学研究所 | 一种结合干扰抑制方法的p波段sar成像处理方法 |
CN103675768A (zh) * | 2013-12-17 | 2014-03-26 | 中国人民解放军国防科学技术大学 | 一种星载sar射频干扰抑制方法 |
Non-Patent Citations (1)
Title |
---|
郑慧芳 等: "SAR 窄带干扰抑制的子带子空间滤波技术研究", 《电子与信息学报》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106501781B (zh) * | 2016-11-04 | 2019-05-31 | 中国人民解放军国防科学技术大学 | 基于斜投影滤波的sar射频干扰抑制方法 |
CN106501781A (zh) * | 2016-11-04 | 2017-03-15 | 中国人民解放军国防科学技术大学 | 基于斜投影滤波的sar射频干扰抑制方法 |
CN106597443A (zh) * | 2016-12-20 | 2017-04-26 | 中国人民解放军国防科学技术大学 | 低频合成孔径雷达图像对比相消射频干扰抑制方法 |
CN106680790A (zh) * | 2016-12-29 | 2017-05-17 | 中国电子科技集团公司第五十四研究所 | 一种伪码调相连续波雷达干扰检测方法 |
CN106680790B (zh) * | 2016-12-29 | 2019-05-17 | 中国电子科技集团公司第五十四研究所 | 一种伪码调相连续波雷达干扰检测方法 |
CN110554367A (zh) * | 2018-06-04 | 2019-12-10 | 重庆测威科技有限公司 | 一种基于压缩感知的目标散射特性测量干扰去除方法 |
CN110554367B (zh) * | 2018-06-04 | 2023-06-06 | 重庆测威科技有限公司 | 一种基于压缩感知的目标散射特性测量干扰去除方法 |
CN110221256B (zh) * | 2019-06-21 | 2023-01-17 | 西安电子科技大学 | 基于深度残差网络的sar干扰抑制方法 |
CN110221256A (zh) * | 2019-06-21 | 2019-09-10 | 西安电子科技大学 | 基于深度残差网络的sar干扰抑制方法 |
CN110618428A (zh) * | 2019-09-19 | 2019-12-27 | 上海海洋大学 | 一种用于多普勒三维数据处理方法及软件 |
CN110531327B (zh) * | 2019-10-11 | 2021-08-31 | 上海无线电设备研究所 | 一种基于自动增益控制的雷达抗干扰探测方法 |
CN110531327A (zh) * | 2019-10-11 | 2019-12-03 | 上海无线电设备研究所 | 一种基于自动增益控制的雷达抗干扰探测方法 |
CN111580107A (zh) * | 2020-07-06 | 2020-08-25 | 内蒙古工业大学 | 一种sar回波信号的射频干扰抑制方法、装置及成像方法 |
CN111580107B (zh) * | 2020-07-06 | 2023-03-31 | 内蒙古工业大学 | 一种sar回波信号的射频干扰抑制方法、装置及成像方法 |
CN113064122A (zh) * | 2021-02-25 | 2021-07-02 | 上海卫星工程研究所 | P波段sar干扰抑制算法的性能评价方法、系统及介质 |
CN113064122B (zh) * | 2021-02-25 | 2023-02-17 | 上海卫星工程研究所 | P波段sar干扰抑制算法的性能评价方法、系统及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105974376B (zh) | 2018-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105974376A (zh) | 一种sar射频干扰抑制方法 | |
CN104635221B (zh) | 一种基于内定标数据的子带拼接方法 | |
CN103278820B (zh) | 临近空间慢速平台sar动目标检测方法及成像方法 | |
CN102520396A (zh) | 基于复数经验模态分解的时变窄带干扰抑制方法 | |
CN102608586B (zh) | 一种基于神经网络的合成孔径雷达射频干扰抑制方法 | |
CN104698459B (zh) | 一种应用于缺失数据的条带sar压缩感知成像方法 | |
CN106707258A (zh) | 一种非高斯背景下微动目标多参数估计方法 | |
CN104849708B (zh) | 基于频域多项式相位变换的高速机动目标参数估计方法 | |
CN106597443B (zh) | 低频合成孔径雷达图像对比相消射频干扰抑制方法 | |
CN106842148A (zh) | 基于frft的线性调频雷达干扰快速抑制方法 | |
CN102721952A (zh) | 基于超分辨多普勒谱的地面运动目标稳健分类方法 | |
CN103675768B (zh) | 一种星载sar射频干扰抑制方法 | |
CN109324315A (zh) | 基于双层次块稀疏性的空时自适应处理雷达杂波抑制方法 | |
CN103674001A (zh) | 一种基于增强自适应时频峰值滤波的光纤陀螺去噪方法 | |
CN103901416A (zh) | 一种基于稳健主成分分析法的多通道杂波抑制方法 | |
CN109031299B (zh) | 低信噪比条件下基于相位差分的isar平动补偿方法 | |
CN104111454B (zh) | 一种扫描雷达角超分辨率方法 | |
Yang et al. | Two-dimensional spectral analysis filter for removal of LFM radar interference in spaceborne SAR imagery | |
CN102901966B (zh) | 基于解斜和楔石处理的动目标瞬时距离多普勒成像方法 | |
Chen et al. | Suppressive interference suppression for airborne SAR using BSS for singular value and eigenvalue decomposition based on information entropy | |
CN110865345A (zh) | 一种快速自适应脉冲压缩方法 | |
CN106569188A (zh) | 基于改进pga的电离层相位污染校正算法 | |
CN104101868B (zh) | 基于干扰子空间重构的雷达多假目标干扰抑制方法 | |
Nguyen et al. | Radio-frequency interference separation and suppression from ultrawideband radar data via low-rank modeling | |
CN106772373B (zh) | 针对任意地面运动目标的sar成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |