一种sar射频干扰抑制方法

文档序号:10611767阅读:789来源:国知局
一种sar射频干扰抑制方法
【专利摘要】本发明提供一种SAR射频干扰抑制方法。技术方案是:以SAR数据的一条距离线回波信号为输入,根据回波信号构建相空间矩阵,利用子空间投影方法得到参考信号的相空间矩阵;采用对角平均化方法,从参考信号相空间矩阵中恢复出一维的参考信号;根据参考信号,确定干扰信号的判别门限;根据确定的干扰判别门限,确定回波中存在的干扰频点,并用参考信号的幅度均值作为替换,完成该条距离线的射频干扰抑制。本发明可以有效提升干扰抑制后的SAR图像质量。
【专利说明】
一种SAR射频干扰抑制方法
技术领域
[0001] 本发明属于航天和微波遥感的交叉技术领域,特别涉及一种SAR(Synthetic Aperture Radar,合成孔径雷达)数据中RFI(Radio Frequency Interference,射频干扰) 抑制方法。
【背景技术】
[0002] SAR(Synthetic Aperture Radar,合成孔径雷达)通过主动发射电磁波来获取目 标信息,是一种全天候、全天时的高分辨率成像系统,在目标侦察和地球遥感应用中具有巨 大的潜力。然而,SAR(尤其是低波段SAR)在复杂的电磁环境中容易受到RFI的影响,进而降 低图像质量,影响图像识别和应用。
[0003] SAR系统由于其良好的穿透特性,在隐蔽目标侦察以及森林植被测绘方面有着独 特的优势,但其工作频段内存在大量的广播、电视信号以及地面警戒雷达发射的信号。这些 干扰信号会导致SAR图像性能严重恶化。
[0004] 陷波法和子空间投影法是当前广泛使用的RFI抑制方法,但传统的陷波法容易造 成有用信息的损失,而子空间投影法的局限性在于过分依赖于子空间的正交性,当干扰信 号功率较大时,会造成子空间误判,在图像中未受干扰区域引入虚警。针对这一现状,提出 了一种基于联合滤波的RFI抑制方法,该方法综合了陷波法和子空间投影法的优点,在保证 干扰抑制效果的前提下最大限度地保留了有用信号,且性能稳健。

【发明内容】

[0005] 本发明的目的是:提出了一种基于联合滤波的RFI抑制方法,该方法是一种非参数 方法,充分结合了子空间投影法和陷波法在干扰抑制方面的优点,具有更好的干扰抑制效 果,且性能稳健。
[0006] 本发明技术方案的思路是:首先构建回波的相空间矩阵,利用子空间投影法将回 波的相空间矩阵投影到目标信号子空间,获得参考信号,并依据该参考信号,确定干扰信号 的判别门限,然后,进行干扰检测、陷波和赋值处理,即可将干扰信号抑制掉。
[0007] 本发明技术方案是:
[0008] 假设X(n),l彡η彡Nr表示某一慢时亥ijSAR接收的回波信号,其中Nr表示距离向采样 点数。
[0009] 第一步:提取参考信号
[0010] (1)将雷达接收的回波信号x(n)按照时延嵌套构造成相空间矩阵S。
[0011]
[0012]式中,L的大小取决于射频干扰的个数,且MiNr+l-Ll过小则特征分解不充分,L 过大则计算量过大,并且成像质量也不一定得到提升,通常根据实际情况确定。本发明的具 体实施例中L取32。令协方差矩阵R = SSh。
[0013] (2)对协方差矩阵R进行特征值分解
[0014] R = UAUh
[0015] 其中 Λ =diag[A1,A2,.",AL],A1>A 2多…彡是一个LXL的对角矩阵,Ai,l彡i$L 为特征值,对应的特征矢量为U= [U1,U2,…,UL]。
[0016]将前1^个特征值所对应的特征矢量提取出,组成参考子空间1^=[111,112,"_,1110,根 据实际情况确定,或根据文献"郑慧芳,杨淋,冯锦.SAR窄带干扰抑制的子带子空间滤波技 术研究[J].电子与信息学报,2013,35( 12): 2836-2842"。以S = S-UtUfS表示参考信号相空 间矩阵。
[0017] (3)采用对角平均化方法从参考信号相空间矩阵Sk中恢复出1维参考信号序列x k (n),l<n<Nr。该过程可以看作是嵌套的逆过程,只需要计算Sk的每条与逆对角线平行的副 逆对角线的均值,该N r个均值就构成了 lXNr维参考信号Xk(n),l<n<Nr。其公式为
[0018]
[0019]第二步,对回波信号进行滤波
[0020] 将回波信号X(n)和参考信号a(n)分别进行FFT变换到频域
[0021] X=[X(l),X(2),.",X(Nr)]
[0022] Xk=[Xk(l),Xk(2),---,Xk(Nr)]
[0023] 利用下式计算判别门限γ :
[0024]
[0025]
[0026]
[0027] 利用下式对回波信号进行滤波,得到射频干扰抑制后的回波信号
[0028]
[0029] 针对SAR的每个慢时刻的回波信号重复上述操作,即可达到抑制干扰的目的。
[0030] 采用本发明可取得以下技术效果:
[0031] 本发明提出的基于联合滤波的SAR射频干扰抑制方法,利用获得的参考信号确定 干扰判别门限,并通过滤波去除干扰,根据参考信号对回波信号重新赋值。相比现有的方 法,本方法在确定参考信号时利用子空间投影法可以更好地去除干扰的影响,得到更优的 干扰判别门限,从而更准确地找到干扰的位置。根据参考信号对回波信号重新赋值时,对干 扰频点的赋值更准确,进而实现更好的干扰抑制效果,并且性能稳健。
【附图说明】
[0032]图1为本发明提供的星载SAR射频干扰抑制方法的原理流程图;
[0033]图2是仿真实验中使用的点目标的成像结果;
[0034]图3是对图2添加射频干扰后点目标的成像结果;
[0035] 图4为采用陷波法对图3进行干扰抑制后成像结果;
[0036] 图5为采用子空间投影法对图3进行干扰抑制后成像结果;
[0037] 图6为采用本文发明提供的方法对图3进行干扰抑制后成像结果;
[0038]图7是仿真实验中使用的实测数据的成像结果;
[0039]图8是对图7添加射频干扰后数据的成像结果;
[0040] 图9是利用陷波法对图8进行干扰抑制后的成像结果;
[0041] 图10是利用子空间投影法对图8进行干扰抑制后的成像结果;
[0042] 图11是利用本发明提供的方法对图8进行干扰抑制后的成像结果。
【具体实施方式】
[0043]图1为本发明提供的星载SAR射频干扰抑制方法的原理流程图。整个流程包括:以 SAR数据的一条距离线回波信号为输入,根据回波信号构建相空间矩阵,利用子空间投影方 法得到参考信号的相空间矩阵;采用对角平均化方法,从参考信号相空间矩阵中恢复出一 维的参考信号;根据参考信号,确定干扰信号的判别门限;根据确定的干扰判别门限,确定 回波中存在的干扰频点,并用参考信号的幅度均值作为替换,完成该条距离线的射频干扰 抑制。
[0044] 图2~图6是进行点目标仿真实验结果。仿真了单个点目标回波,在此回波中添加 了带宽为0.5MHz,1.0MHz和1.5MHz的射频干扰信号,干信比为20dB。图2为原始点目标回波 成像结果,从图像中可以看出,点目标聚焦效果良好。图3是在图2中添加射频干扰后回波的 成像结果,可以看到点目标已经被干扰信号湮没。图4是采用陷波法对图3进行干扰抑制后 的成像结果,从图像中可以看出,干扰聚焦良好,但仍有较多的干扰残留。图5是采用子空间 投影法对图3进行干扰抑制后的成像结果,点目标不仅不能良好聚焦,图像还残留较多的干 扰。图6是采用本发明提供的方法进行干扰抑制后的成像结果,可以看到,点目标聚焦良好, 且干扰残留较少。
[0045] 图7~图11是进行面目标仿真实验的结果。仿真实验选取了 一幅无干扰的 RadarSat-Ι回波实测数据,在此数据中添加了带宽为0.5MHz,1.0MHz和1.5MHz的射频干扰 信号,干信比为20dB。
[0046] 图7是仿真实验中使用的原始数据的成像结果。从图中可以看出,图像清晰,对比 度高。图8是对图7添加射频干扰后数据的成像结果。从图中可以看出,由于干扰的存在,图 像目标信息被覆盖,图像对比度下降。图9和图10分别是利用陷波法和子空间投影法对图8 进行干扰抑制后的成像结果,干扰抑制后,图像质量得到提升,不过仍有较多干扰残留,细 节信息不明显。图11是利用本发明提供的方法对图8进行干扰抑制后的成像结果,干扰抑制 后,图像质量明显提升,细节信息恢复较为明显,图像对比度有明显提高。
[0047]以上仿真实验说明了本发明方法可以有效地弥补现有方法在干扰抑制方面的不 足,有效提升干扰抑制后的SAR图像质量。
【主权项】
1. 一种SAR射频干扰抑制方法,SAR是指合成孔径雷达,假设X(η),1《n《Nr,表示某一慢 时刻SAR接收的回波信号,其中Nr表示距离向采样点数,其特征在于,包括下述步骤: 第一步:提取参考信号: 将雷达接收的回波信号x(n)按照时延嵌套构造成相空间矩阵S:上式中M=Nr+l-L,L根据实际情况确定; 令协方差矩阵R= SSH,对协方差矩阵R进行特征值分解: R=U 八 IJH 其中Λ=diag[λl,λ2,…,λL]是一个LXL的对角矩阵,对应的特征矢量矩阵为U=[山, U2,…,UL],并且λ冷λ2>…,λ功特征值,i = 1,2,…レ上标Η表示求矩阵的共辆转置; 将前k个特征值所对应的特征矢量提取出,组成参考子空间化=[111,112,-,,1^]^的取值 根据实际情况确定; 令参考信号相空间矩阵=S-U;UfS,采用对角平均化方法从参考信号相空间矩阵Sk 中恢复出1维参考信号序列xk(n); 第二步,对回波信号进行滤波 将回波信号x(n)和参考信号xk(n)分别进行快速傅里叶变换得到X和Xk: Χ=[Χ(1),Χ(2),···,Χ(Νγ)] Xk=[Xk(l),Xk(2),---,Xk(Nr)] 利用下式计算判别口限丫:利用下式对回波信号进行滤波,得到射频干扰抑制后的回波信号^的);针对SAR的每个慢时刻回波信号重复上述两个步骤。
【文档编号】G01S7/36GK105974376SQ201610298996
【公开日】2016年9月28日
【申请日】2016年5月9日
【发明人】张永胜, 吴鹏, 张启雷, 董臻, 黄海风, 余安喜, 何峰, 孙造宇, 金光虎, 何志华, 杜湘瑜
【申请人】中国人民解放军国防科学技术大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1