首页 / 专利库 / 信号处理 / 信噪比 / 峰值信噪比 / 地波雷达海杂波中目标的三维检测方法

地波雷达海杂波中目标的三维检测方法

阅读:403发布:2020-05-13

专利汇可以提供地波雷达海杂波中目标的三维检测方法专利检索,专利查询,专利分析的服务。并且地波雷达海杂波中目标的三维检测方法,包括构造地波雷达通道-距离-多普勒回波谱,提取海杂波区域,以斜投影技术构建多级 空域 滤波器 ;将地波雷达通道-距离-多普勒回波谱馈送给每个空域滤波器,各滤波器依次对地波雷达 覆盖 的每个方位作为潜在目标方向,依次进行滤波;对各空域滤波后的谱数据进行取最小值的合成操作,得到合成谱,对合成谱的距离-多普勒-方位谱进行三维检测,以抑制非均匀海杂波,搜寻峰值,最终输出舰船目标的方位、距离和多普勒参数。本 发明 在传统距离-多普勒维目标二维检测 基础 上拓展方位维检测,进而形成方位-距离-多普勒维目标检测,通过联合海杂波和目标在这三个联合维度的回波差异,提高了海杂波中的目标检测能 力 。,下面是地波雷达海杂波中目标的三维检测方法专利的具体信息内容。

1.地波雷达海杂波中目标的三维检测方法,其特征在于包括以下步骤:
步骤1:构造地波雷达通道-距离-多普勒回波谱,提取海杂波区域:
对地波雷达接收通道采集的下变频时域采样信号,先后进行冲压缩和相干积累,得到地波雷达通道-距离-多普勒回波谱,采用二维信噪比法提取海杂波区域,记为x(c,r,f),其物理量为幅度,单位是dB;c代表接收通道,r代表距离,f代表多普勒;
步骤2:构建多级空域滤波器
采用斜投影技术构建多级空域滤波器;
步骤3:频谱合成:
将地波雷达通道-距离-多普勒回波谱馈送给每一个空域滤波器,各个滤波器依次对地波雷达覆盖的每一个方位作为潜在目标方向,并依次进行滤波,在对每一个方向进行滤波后,保留了该方向的谱信号,记为 其物理量为幅值,单位是dB;a代表方位;
对各个空域滤波后的谱数据进行取最小值的合成操作,即在每个空域滤波器的滤波结果中只保留其中的最小值,所得结果称为合成谱;
经过谱合成后,来自目标方位的期望信号在不引入幅值和相位畸变的情况下得到保留,除与目标方位相同的杂波外,其余大部分杂波被消除;
当将一个方位设置为目标的潜在方向时,只保留了该方向的距离多普勒谱和无方向性噪声,而抑制了其他方向的谱信号;
步骤4:目标检测及参数估计:
对步骤3产生的合成谱,对其距离维、多普勒维、方位维分别设置阈值,每个阈值由均值类恒虚警检测器确定,对合成谱的距离-多普勒-方位谱进行三维检测,以抑制非均匀海杂波,然后是峰值搜寻过程,最终输出舰船目标的方位、距离和多普勒参数。
2.如权利要求1所述的地波雷达海杂波中目标的三维检测方法,其特征在于所述步骤2采用斜投影技术构建多级空域滤波器如下:
当目标方向的空域导向矢量为A,且其他观测方向为潜在海杂波空域导向矢量B,则定义沿B投影到A的单个斜投影空域滤波EAB为
其中上标-1是矩阵的广义逆,上标H是矩阵的厄米共轭转置, 是矩阵B的正交补投影,定义为
该斜投影滤波EAB具有以下特性:
上述特性为构造多级空域滤波器奠定了基础,实现该多级滤波器的具体步骤如下:
1)确定使用的多级滤波器的数量,其数量P由地波雷达覆盖方位的空间范围决定,空域采样间隔计算为
P=int(|θL-θT+Δθ|+|θT+Δθ-θR|)/Δθ  (4)
其中int(.)表示取整数,θL和θR为地波雷达覆盖的左右方位角边界,θT是潜在的目标方向,Δθ为滤波时所选的凹口间隔。
2)在地波雷达覆盖的方位角范围内,选择杂波角度或干扰方位角进行空间滤波,一次一个角度,以Δθ为间隔且杂波角度或干扰方位角满足
θC∈[θL,θT-Δθ]∪[θT+Δθ,θR]  (5)
3)对于步骤2选择的杂波角度或干扰方位角,记其目标子空间和杂波子空间分别为及 然后根据式(1)构造各空域滤波器EAB;
4)重复步骤2)和3)P次,从而得到P个基于斜投影的滤波器,并将每个基于斜投影滤波器相乘,形成一个多级空域滤波器。
其中p代表第p个滤波器。
3.如权利要求1所述的地波雷达海杂波中目标的三维检测方法,其特征在于所述步骤3将通道-距离-多普勒谱馈送给多级空域滤波器的具体方式如下:利用逻辑积对多个空域滤波后的谱数据进行合成,最终三维目标谱为
其中a,r,f表示方位、距离和多普勒,min(.)表示取括号内的最小值,s(r,fj)和i(r,fk)分别代表船只目标和海杂波距离-多普勒谱;K代表海杂波的多普勒频点数,k代表第k个海杂波的多普勒频点;J代表目标的多普勒频点数,且K>J暗含船只回波频点小于海杂波频点;n(r,f)是高斯白噪声;Aj代表船只的导向矢量,Bk代表海杂波的导向矢量。
船只和海杂波的导向矢量分别为
Bk=[a(θk)…a(θK)]T  (9)
其中θ是方位角,且θj表示第j个目标多普勒频点的方位角,j为取值1到J的自然数;θk表T
示第k个海杂波多普勒频点的方位角,k为取值1至K的自然数,当k取值K时,公式9即[a(θK)] ;
d是阵元间距,M是阵元数,f0是雷达载波频率,c是光速,[.]T为矩阵转置。
4.如权利要求1所述的地波雷达海杂波中目标的三维检测方法,其特征在于所述步骤
4:对步骤3的合成谱,通过均值类恒虚警检测器在距离、方位以及多普勒域分别设置阈值,对方位-距离-多普勒谱进行三维检测,
其中采用峰值检测的三维复合目标检测表示为
D=peak(Dr||Da||Df)  (10)
其中D存储被探测船只的方位、距离和多普勒位置参数,peak(.)表示在括号中搜索幅度峰值,Dr,Da和Df分别存储在距离、方位和多普勒维的谱幅度值,这些幅度值均大于各自距离、方位以及多普勒域的检测阈值Tr,Ta和Tf,||表示逻辑和运算, 是特定多普勒为f和方位为a处的距离维幅度值, 是多普勒为f和距离为r处的方位维幅度值,是特定方位为a和距离为r处的多普勒维幅度值;
通过式(10)的峰值检测,最终输出舰船目标的方位、距离和多普勒参数。

说明书全文

地波雷达海杂波中目标的三维检测方法

技术领域

[0001] 本发明涉及一种地波雷达海杂波中目标的检测方法,具体涉及一种地波雷达海杂波中目标的三维检测方法。

背景技术

[0002] 地波雷达是海上船只目标大范围、超视距、连续探测的主要手段,它利用垂直极化高频电磁波(3-30MHz)沿海面绕射的特性,实现海上目标的超视距探测。地波雷达海杂波是该雷达发射的高频无线电波与海浪相互作用而散射回来的信号,分为一阶、二阶及高阶海杂波。其中,一阶海杂波表现为一对展宽的双峰,能量很强,造成船只检测盲区[1],降低了探测率。因此需要通过发展海杂波内的船只目标检测方法,突破“检测盲区”的限制,提高地波雷达海上目标的探测性能。
[0003] 常规的地波雷达目标检测方法通常采用的是二维阈值检测方法,即在距离维和多普勒维分别检测一次再将两次的检测结果相“与”,进而输出目标检测结果,同时输出目标的距离及多普勒参数。后期在对应的目标距离-多普勒谱点位置取通道数据,采用比幅或者查找表等波束形成方法测向,得到目标的方位参数。在二维阈值检测的过程中,多采用幅值检测。方法有均值类恒虚警检测器、排序统计类恒虚警检测器和自适应回归恒虚警检测器[2]。其中,均值类恒虚警检测器检测方法适合均匀背景,当参考窗滑窗范围越大,恒虚警检测损失越小,但无法形成有效的局部估计。排序统计类恒虚警检测器检测方法对弱杂波具有一定的抑制能,在强杂波区检测性能下降。自适应回归恒虚警检测方法对电离层及低速区杂波具有抑制效果,但在均匀检测背景区,检测性能不如均值类恒虚警检测器和自适应回归恒虚警检测器。该类二维检测方法通常在实际应用中需要通过区域分割,然后针对不同杂波区及类型,采用不同的阈值检测方法来达到较好的检测性能,但海杂波内的目标检测效果仍不理想。
[0004] 海杂波是来自随机方位的海面雷达回波,在通道-距离-多普勒谱中表现为在不同的距离及多普勒位置暗含不同的空域回波信息。考虑到产生谐振的海浪可能分布于大片的海洋表面,因而会有多个且随机的潜在海杂波回波方向,同时目标则来自于某一特定方位。海杂波及船只目标的这种空域回波差异体现在地波雷达接收阵列获取的通道数据中。因而,在目标检测的过程中应联合距离及多普勒域信息利用该空域差异,而这是现有技术所缺乏的。
[0005] 相关的参考文献如下:
[0006] [1]Leong H,Ponsford A.The effects of sea clutter on the performance of HF surface wave radar in ship detection.IEEE Radar Conference.Rome:IEEE,2008:1-6.
[0007] [2]Ji Yonggang,Xu Leda,Wang Yiming,and Chu Xiaoliang.Ship Detection in Strong Clutter Environment Based on Adaptive Regression Thresholding for HFSWR,2014International Conference on Computer Science and Electronic Technology(ICCSET 2014),Atlantis Press,pp.352-355.

发明内容

[0008] 本发明目的是提供一种地波雷达海杂波中目标的三维检测方法,该方法提出的三维目标检测方法不同于传统的距离-多普勒二维目标检测,增加了方位维检测,进而形成方位-距离-多普勒维目标检测。充分利用目标与海杂波在空间、距离及频率联合维的不同,实现海杂波中的目标检测。
[0009] 为了实现上述目的,本发明给出了地波雷达海杂波中目标三维检测方法的技术方案,其特征在于包括以下步骤:
[0010] 步骤1:构造地波雷达通道-距离-多普勒回波谱,提取海杂波区域:
[0011] 对地波雷达接收通道采集的下变频时域采样信号,先后进行冲压缩和相干积累,得到地波雷达通道-距离-多普勒回波谱,采用二维信噪比法提取海杂波区域,记为x(c,r,f),其物理量为幅度,单位是dB;c代表接收通道,r代表距离,f代表多普勒;
[0012] 步骤2:构建多级空域滤波器
[0013] 在地波雷达的舰船探测场景中,海杂波可能会具有各种多普勒和方位,甚至可能每个多普勒都有多个方位。但是每艘船只从特定的不同的多普勒和方位产生回波。因此,需要设置多个凹口来消除目标方向角以外的所有方向的杂波干扰。
[0014] 为了设置这些方向上的凹口,需要构造多个空域滤波器,每个滤波器都有一个相应的凹口,从而构成多级空域滤波器。斜投影空域滤波器级联相乘时具有等效的滤波效果,即来自目标方位角的信号不失真,而海杂波则会被连续的凹口剔除。本发明的各个空域滤波器采用斜投影技术构建,保留目标信号,抑除杂波干扰。
[0015] 步骤3:频谱合成:
[0016] 将地波雷达通道-距离-多普勒回波谱馈送给每一个空域滤波器,各个滤波器依次对地波雷达覆盖的每一个方位作为潜在目标方向,并依次进行滤波,在对每一个方向进行滤波后,保留了该方向的谱信号,记为 其物理量为幅值,单位是dB;a代表方位;对各个空域滤波后的谱数据进行取最小值的合成操作,即在每个空域滤波器的滤波结果中只保留其中的最小值,所得结果称为合成谱。经过谱合成后,来自目标方位的期望信号在不引入幅值和相位畸变的情况下得到保留,除与目标方位相同的杂波外,其余大部分杂波被消除。应该注意的是,无方向的噪声是不能被抑制的。由于噪声是高斯的,且相对较弱,其对目标谱的影响可以忽略不计。当将一个方位设置为目标的潜在方向时,只保留了该方向的距离多普勒谱和无方向性噪声,而抑制了其他方向的谱信号;
[0017] 步骤4:目标检测及参数估计:
[0018] 对步骤3产生的合成谱,对其距离维、多普勒维、方位维分别设置阈值,每个阈值由均值类恒虚警检测器确定,对合成谱的距离-多普勒-方位谱进行三维检测,以抑制非均匀海杂波,然后是峰值搜寻过程,最终输出舰船目标的方位、距离和多普勒参数。
[0019] 与现有技术相比,本发明的创新之处体现在以下方面:
[0020] 本专利通过发展海杂波内的船只目标检测方法,突破了传统意义上海杂波导致的“检测盲区”的限制,提高了地波雷达海上目标的探测性能。
[0021] 不同于传统的距离-多普勒二维目标检测,增加了方位维检测,进而形成方位-距离-多普勒维目标检测。充分利用目标与海杂波在空间、距离及频率联合维的不同,实现海杂波中的目标检测。
[0022] 海杂波及船只目标的这种空域回波差异体现在地波雷达接收阵列获取的通道数据中。因而,在目标检测的过程中应联合距离及多普勒域信息利用该空域差异。但地波雷达的接收整列孔径相对于其波长过小,导致即使采用常规的波束形成方法,其空域滤波效果也有限,无法形成有效的三维检测。本发明采用多级斜投影空域滤波的方式,达到了更好的空域滤波效果,使得三维联合检测成为可能。附图说明
[0023] 图1本发明的目标检测流程。
[0024] 图2一个接收通道的地波雷达距离-多普勒谱。
[0025] 图3不同空间主角下的信干噪比改善和干噪比。
[0026] 图4不同滤波凹口间隔下的空间主角和方位。
[0027] 图5实测地波雷达处理结果,
[0028] 其中,(a)常规DBF后的方位-距离-多普勒谱,(b)本专利的方位-距离-多普勒谱,(c)多普勒维结果对比,(d)距离维结果对比,(e)方位维结果对比。

具体实施方式

[0029] 下面结合公式和附图,对本发明的方法做进一步说明:
[0030] 如图1所示,利用高频地波雷达多域信息的强干扰环境下的一阶回波谱提取方法,主要包括回波谱构造及海杂波区域提取,多级空域滤波和频谱合成及目标检测及参数估计,最终输出检测到的目标距离、多普勒及方位信息。其具体步骤如下:
[0031] 步骤1:利用周期为T秒(本发明以291s为例)的地波雷达多通道下变频信号,经匹配滤波和相干积累处理后,形成通道-距离-多普勒谱。采用二维信噪比法提取海杂波区域,记为x(c,r,f),其物理量为幅度,单位是dB;c代表接收通道,r代表距离,f代表多普勒。其形式如图2所示。
[0032] 步骤2:该步骤采用斜投影方法构造多级空域滤波器,当目标方向的空域导向矢量为A,且其他观测方向为潜在海杂波空域导向矢量B,则定义沿B投影到A的单个斜投影空域滤波EAB为
[0033]
[0034] 其中上标-1是矩阵的广义逆,上标H是矩阵的厄米共轭转置, 是矩阵B的正交补投影,定义为
[0035]
[0036] 该斜投影滤波EAB具有以下特性:
[0037]
[0038] 上述特性为构造多级空域滤波器奠定了基础,实现该多级滤波器的具体步骤如下:
[0039] 1)确定使用的多级滤波器的数量,其数量P由地波雷达覆盖方位角的空间范围决定,
[0040] 空域采样间隔计算为
[0041] P=int(|θL-θT+Δθ|+|θT+Δθ-θR|)/Δθ  (4)
[0042] 其中int(.)表示取整数,θL和θR为地波雷达覆盖的左右方位角边界,θT是潜在的目标方向,Δθ为滤波时所选的凹口间隔。
[0043] 在实际应用中,目标方位与滤波凹口位置之间的空间主角对信干噪比有很大的影响。因此,在选择斜投影空域滤波参数之前,应确定不同干扰噪声比下的信干噪比改进量,这与空间主角密切相关。该空间主角由滤波方位的凹口间隔Δθ决定的。
[0044] 根据数值分析得到的图3可以看出,当主角很小为1°且信杂比SIR为-10dB时,信干噪比SINR很难得到改善。而当主角为90°且干扰大于噪声时,可以保证SINR得到改善。当主角为13°且INR大于15B时,SINR将得到改善,符合地波雷达的实际应用情况。在这种情况下,从图4可以看出,空间主角随着方位角间距的增大和目标方位角的接近而增大。因此,为了在[-40°,40°]的方位角覆盖范围内保持13°以上的主角,应选择至少2°的空间滤波凹口间隔。
[0045] 2)在地波雷达覆盖的方位角范围内,选择杂波角度或干扰方位角进行空间滤波,一次
[0046] 一个角度,以Δθ为间隔且杂波角度或干扰方位角满足
[0047] θC∈[θL,θT-Δθ]∪[θT+Δθ,θR]  (5)
[0048] 3)对于步骤2选择的杂波角度或干扰方位角,记其目标子空间和杂波子空间分别为 及 然后根据式(1)构造各空域滤波器EAB;
[0049] 4)重复步骤2)和3)P次,从而得到P个基于斜投影的滤波器,并将每个基于斜投影滤波器相乘,形成一个多级空域滤波器。
[0050]
[0051] 其中p代表第p个滤波器。
[0052] 步骤3:将步骤1提取的海杂波通道-距离-多普勒谱馈送到步骤3形成的多级空域滤波器Ecas,得到地波雷达的距离-多普勒-方位谱,见图5。
[0053] 将通道-距离-多普勒谱馈送给多级空域滤波器的具体方式如下:利用逻辑积对多个空域滤波后的谱数据进行合成,最终三维目标谱为
[0054]
[0055] 其中a,r,f表示方位、距离和多普勒,min(.)表示取括号内的最小值,s(r,fj)和i(r,fk)分别代表船只目标和海杂波距离-多普勒谱;K代表海杂波的多普勒频点数,k代表第k个海杂波的多普勒频点;J代表目标的多普勒频点数,且K>J暗含船只回波频点小于海杂波频点;n(r,f)是高斯白噪声;Aj代表船只的导向矢量,Bk代表海杂波的导向矢量。
[0056] 船只和海杂波的导向矢量分别为
[0057]
[0058] Bk=[a(θk) … a(θK)]T  (9)
[0059] 其中θ是方位角,且θj表示第j个目标多普勒频点的方位角,j为取值1到J的自然数;θk表示第k个海杂波多普勒频点的方位角,k为取值1至K的自然数,当k取值K时,公式9即[a(θK)]T;d是阵元间距,M是阵元数,f0是雷达载波频率,c是光速,[.]T为矩阵转置。
[0060] 由式(7),经过谱合成后,来自目标方位的期望信号在不引入幅值和相位畸变的情况下得到恢复,除与目标方位相同的杂波外,其余大部分杂波被消除。应该注意的是,无方向的噪声是不能被抑制的。由于噪声是高斯的,且相对较弱,其对目标谱的影响可以忽略不计。当将一个方位设置为目标的潜在方向时,只保留了该方向的距离多普勒谱和无方向性噪声,而抑制了其他方向的谱信号。滤波依次选择地波雷达覆盖的每一个方位作为潜在目标方向,并保留该方向的谱信号。
[0061] 步骤4:步骤3的合成谱,通过均值类恒虚警检测器在距离、方位以及多普勒域分别设置阈值,对方位-距离-多普勒谱进行三维检测。采用峰值检测的三维复合目标检测表示为
[0062] D=peak(Dr||Da||Df)  (10)
[0063]
[0064]
[0065]
[0066] 其中D存储被探测船只的方位、距离和多普勒位置参数,peak(.)表示在括号中搜索幅度峰值,Dr,Da和Df分别存储在距离、方位和多普勒维的谱幅度值,这些幅度值均大于各自距离、方位以及多普勒域的检测阈值Tr,Ta和Tf,||表示逻辑和运算, 是特定多普勒为f和方位为a处的距离维幅度值, 是多普勒为f和距离为r处的方位维幅度值,是特定方位为a和距离为r处的多普勒维幅度值。
[0067] 通过式(10)的峰值检测,最终输出舰船目标的方位、距离和多普勒参数。
[0068] 图5给出了本专利方法的实测地波雷达处理结果以及与常规DBF方法的对比。由图5(a)常规DBF后的方位-距离-多普勒谱,海杂波的条带特征仍然非常明显,难以发现其中的点状目标(目标所在的方位为10°,距离82.5~85Km,多普勒0.2198~0.2266Hz)。图5(b)本专利的方位-距离-多普勒谱,消除了海杂波的条带特征,点状目标显现。本专利对目标的信干噪比提高接近10dB,高于采用DBF方法的5dB。检测到的目标方位为8°,距离为84Km,多普勒为0.2266Hz。图5(c)、(d)、(e)给出了多普勒维、距离维、方位维的结果对比,进一步表明本专利方法的有效性。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈