一种基于背景感知的超视距雷达多策略CFAR检测方法

一种基于背景感知的超视距雷达多策略cfar检测方法
技术领域
1.本发明属于雷达回波目标检测领域，具体涉及一种基于背景感知的超视距雷达多策略cfar检测方法。


背景技术：

2.超视距雷达工作在高频波段，利用该频段的电磁波沿电离层返回散射传播机理以及可以沿海面低衰减绕射的特点，来实现超视距探测，具有探测距离远、安全系数高的优点，但同时，超视距雷达回波会受到电离层杂波和较严重的海杂波的影响，其检测背景往往不是均匀背景、杂波边缘背景、多目标背景这三类背景中的一种，而是多种背景混杂，并且杂波的分布往往是时变的(尤其是海杂波)，这些因素给目标检测造成了困难。在超视距雷达回波的目标检测中，传统的cfar(恒虚警)检测器往往难以发挥出较好的性能。均值类检测器中，ca-cfar(单元平均恒虚警检测器)在均匀背景下具有最优检测性能，但在超视距雷达的复杂背景下检测性能会下降；go-cfar(最大选择恒虚警检测器)和so-cfar(最小选择恒虚警检测器)分别改善了ca-cfar在杂波边缘和多目标环境下的检测性能，但在复杂、时变的检测背景下检测性能也会降低。有序统计类检测器，例如os-cfar(有序统计恒虚警检测器)、cmld-cfar(筛选平均恒虚警检测器)等，通过删除可能是干扰信号的参考单元，相比于均值类检测器，在多目标环境下检测性能更加优秀，但在复杂背景以及干扰目标数超过限度时检测性能会快速下降。自适应检测器，例如vi-cfar(变化指数恒虚警检测器)、hce-cfar(自适应杂波边缘恒虚警检测器)、e-cfar(删除强干扰恒虚警检测器)等，可以在特定环境下获得较强的检测性能，vi-cfar综合了ca-cfar、go-cfar、so-cfar的优点，在均匀背景下具有较强的检测性能，hce-cfar在杂波边缘背景下检测性能较强，e-cfar通过删除超过删除阈值的强干扰采样的方式，在多目标背景下有较优秀的检测性能，但在超视距雷达复杂、时变的背景下，它们的检测性能均会下降。超视距雷达复杂、时变的检测背景使得均值类、有序统计类、自适应等cfar检测器的目标检测性能下降，无法达到预期的目标检测效果。


技术实现要素：

3.本发明为了解决现有的cfar检测器在超视距雷达回波复杂、时变的检测背景下检测性能下降的问题，进而提供一种基于背景感知的超视距雷达多策略cfar检测方法。
4.一种基于背景感知的超视距雷达多策略cfar检测方法，所述方法是通过以下步骤实现的：
5.步骤一：获取超视距雷达回波ard(方位距离多普勒)谱数据；
6.步骤二：对雷达回波ard谱进行平方律检波；
7.步骤三：对要进行检测的rd(距离多普勒)谱，分别计算整个rd谱的标准差以及rd谱上每个数据点的标准差；
8.步骤四：计算rd谱上各个数据点与整个rd谱之间的kl散度值，将该kl散度值作为
数据点与整个rd谱之间的分布差异，得到kl散度矩阵；
9.步骤五：根据得到的kl散度矩阵，使用类间方差作为依据来得到kl散度矩阵的分割阈值；
10.步骤六：使用分割阈值kld将rd谱分为两个区域：噪声区域和杂波区域，kl散度值大于kld的数据点认为其分布与整体rd谱分布差距较大，更趋向于噪声，将其划分为噪声区域的数据点，kl散度值小于kld的数据点认为其分布与整体rd谱分布差距较小，更趋向于杂波，将其划分为杂波区域的数据点；
11.步骤七：对噪声区域和杂波区域的数据点采用不同的检测策略：对杂波区域的数据点，认为其为非目标点，对噪声区域的数据点，认为其可能是目标点，将其传入下一步检测步骤中；
12.步骤八：对属于噪声区域的数据点，对其进行峰值检测，如果在距离、多普勒、方位向(距离、多普勒、方位向参考单元分别取检测单元左右各1、2、1个数据点)上均为极大值点，则认为可能是目标点，将其传入下一步检测步骤中，否则认为其为非目标点；
13.步骤九：对经过峰值检测为极大值的数据点，使用gosso-cfar(广义有序统计最小选择恒虚警)检测器对其进行目标检测；
14.步骤十：根据gosso-cfar的判定结果得到最终的检测结果，如果gosso-cfar判定为目标点，则最终判定检测单元为目标点，否则，判定为非目标点；
15.进一步地，所述步骤三中整个rd谱的标准差计算公式如下：
[0016][0017]
其中，σ为整个rd谱的标准差，m是整个rd谱内数据点的个数，xi为rd谱内第i个数据点的值，μ为整个rd谱的均值；
[0018]
rd谱上各个数据点的标准差计算公式如下：
[0019][0020]
其中σk代表rd谱上第k个数据点的标准差，n为计算该标准差时使用的计算范围的总点数，由于只使用单个数据点无法计算标准差，所以计算数据点的标准差时使用以该数据点为中心的大小为3
×
3的矩阵作为计算范围，即n＝9，xi为计算范围内第i个数据点的值，μk代表计算范围内数据的均值，s为rd谱上要计算标准差的数据点的个数；
[0021]
进一步地，所述步骤四中kl散度值的计算公式如下：
[0022][0023]
其中，klk代表rd谱上第k个数据点与整个rd谱之间的kl散度值，σ为整个rd谱的标准差，σk为rd谱上第k个数据点的标准差；
[0024]
进一步地，所述步骤五的具体步骤如下：
[0025]
步骤五一：假定一个kl散度值，以该值为分割阈值将kl散度矩阵分割为两部分：区域1和区域2，设区域1内有n1个数据点，区域2内有n2个数据点，令p1、p2分别表示区域1和区域2占总区域的比重，则p1、p2的计算公式如下：
[0026][0027][0028]
步骤五二：令m1、m2分别表示区域1和区域2内kl散度值的均值，则m1、m2的计算公式如下：
[0029][0030][0031]
其中，kli代表区域1中第i个数据点的kl散度值，klj代表区域2中第j个数据点的kl散度值；
[0032]
步骤五三：计算区域1和区域2之间的类间方差σ2，其计算公式如下：
[0033]
σ2＝p1p2(m
1-m2)2；
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0034]
步骤五四：改变设定的kl散度阈值，再计算类间方差，使类间方差最大的kl散度值即为最终的kl散度分割阈值kld，其公式如下：
[0035]
kld＝arg
max
σ2。
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)；
[0036]
进一步地，所述步骤九中gosso-cfar检测流程如下：
[0037]
步骤九一：载入检测单元及其左右滑窗内的经过平方律检波后的数据，设检测单元的值为d，其左侧滑窗内的参考单元数据为x1,x2,x3......xn，右侧滑窗内的参考单元数据为y1,y2,y3......yn，其中n为设定的参考单元数量；
[0038]
步骤九二：对左侧的参考单元按由小到大的顺序排序，选出第k1个最小单元x，对右侧的参考单元按由小到大的顺序排序，选出第k2个最小单元y；
[0039]
步骤九三：根据步骤九一和步骤九二所得参数确定背景噪声的估计值z为：
[0040]
z＝min(x,y)；
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)
[0041]
步骤九四：检测阈值s＝tz，其中t为标称因子，其可由下式计算出来：
[0042][0043]
其中，p
fa
为虚警率，γ()为γ函数，对于整数x，γ(x)＝(x-1)！；
[0044]
步骤九五：根据计算得到的检测阈值，做出判定结果，设为h1有目标，h2为无目标，则判定策略如下：
[0045][0046]
即如果检测单元的值大于阈值，则判定为目标点，反之，则判定为非目标点；
[0047]
本发明与现有技术相比具有以下有益效果：
[0048]
本发明针对超视距雷达检测背景复杂、时变的问题设计了一种基于背景感知的多策略cfar检测器，该检测器根据kl散度值，使用类间方差最大化的方法确定阈值，将回波rd谱分为噪声、杂波两个区域，对两个区域内的检测点采用不同的检测策略进行检测，从而达到可以根据检测背景变化自适应地确定检测策略，进行目标检测的效果，在超视距雷达回波复杂、时变的检测背景下具有良好的检测性能。
附图说明
[0049]
图1为本发明的流程框图
[0050]
图2为本发明中gosso-cfar检测流程图
[0051]
图3为本发明中实测背景注入目标下各方法的检测性能曲线图
[0052]
图4为本发明中实测背景注入目标下各方法的虚警率曲线图；
具体实施方式
[0053]
具体实施方式一：参照图1至图4说明本实施方式，本实施方式提供一种基于背景感知的超视距雷达多策略cfar检测方法，其特征在于：所述方法是通过以下步骤实现的：
[0054]
步骤一：获取超视距雷达回波ard(方位距离多普勒)谱数据；
[0055]
步骤二：对雷达回波ard谱进行平方律检波；
[0056]
步骤三：对要进行检测的rd(距离多普勒)谱，分别计算整个rd谱的标准差以及rd谱上每个数据点的标准差；
[0057]
步骤四：计算rd谱上各个数据点与整个rd谱之间的kl散度值，将该kl散度值作为数据点与整个rd谱之间的分布差异，得到kl散度矩阵；
[0058]
步骤五：根据得到的kl散度矩阵，使用类间方差作为依据来得到kl散度矩阵的分割阈值；
[0059]
步骤六：使用分割阈值kld将rd谱分为两个区域：噪声区域和杂波区域，kl散度值大于kld的数据点认为其分布与整体rd谱分布差距较大，更趋向于噪声，将其划分为噪声区域的数据点，kl散度值小于kld的数据点认为其分布与整体rd谱分布差距较小，更趋向于杂波，将其划分为杂波区域的数据点；
[0060]
步骤七：对噪声区域和杂波区域的数据点采用不同的检测策略：对杂波区域的数据点，认为其为非目标点，对噪声区域的数据点，认为其可能是目标点，将其传入下一步检测步骤中；
[0061]
步骤八：对属于噪声区域的数据点，对其进行峰值检测，如果在距离、多普勒、方位向(距离、多普勒、方位向参考单元分别取检测单元左右各1、2、1个数据点)上均为极大值点，则认为可能是目标点，将其传入下一步检测步骤中，否则认为其为非目标点；
[0062]
步骤九：对经过峰值检测为极大值的数据点，使用gosso-cfar(广义有序统计最小选择恒虚警)检测器对其进行目标检测；
[0063]
步骤十：根据gosso-cfar的判定结果得到最终的检测结果，如果gosso-cfar判定为目标点，则最终判定检测单元为目标点，否则，判定为非目标点。
[0064]
具体实施方式二：参照图1至图4说明本实施方式，本实施方式是对具体实施方式
一所述的步骤三作进一步限定，本实施方式中所述步骤三中整个rd谱的标准差计算公式如下：
[0065][0066]
其中，σ为整个rd谱的标准差，m是整个rd谱内数据点的个数，xi为rd谱内第i个数据点的值，μ为整个rd谱的均值；
[0067]
rd谱上各个数据点的标准差计算公式如下：
[0068][0069]
其中σk代表rd谱上第k个数据点的标准差，n为计算该标准差时使用的计算范围的总点数，由于只使用单个数据点无法计算标准差，所以计算数据点的标准差时使用以该数据点为中心的大小为3
×
3的矩阵作为计算范围，即n＝9，xi为计算范围内第i个数据点的值，μk代表计算范围内数据的均值，s为rd谱上要计算标准差的数据点的个数。其它组成及连接方式与具体实施方式一相同。
[0070]
具体实施方式三：参照图1至图4说明本实施方式，本实施方式是对具体实施方式二所述的步骤四作进一步限定，本实施方式中所述步骤四中kl散度值的计算公式如下：
[0071][0072]
其中，klk代表rd谱上第k个数据点与整个rd谱之间的kl散度值，σ为整个rd谱的标准差，σk为rd谱上第k个数据点的标准差。其它组成及连接方式与具体实施方式二相同。
[0073]
具体实施方式四：参照图1至图4说明本实施方式，本实施方式是对具体实施方式三所述的步骤五作进一步限定，本实施方式中所述步骤五的具体步骤如下：
[0074]
步骤五一：假定一个kl散度值，以该值为分割阈值将kl散度矩阵分割为两部分：区域1和区域2，设区域1内有n1个数据点，区域2内有n2个数据点，令p1、p2分别表示区域1和区域2占总区域的比重，则p1、p2的计算公式如下：
[0075][0076][0077]
步骤五二：令m1、m2分别表示区域1和区域2内kl散度值的均值，则m1、m2的计算公式如下：
[0078][0079][0080]
其中，kli代表区域1中第i个数据点的kl散度值，klj代表区域2中第j个数据点的kl散度值；
[0081]
步骤五三：计算区域1和区域2之间的类间方差σ2，其计算公式如下：
[0082]
σ2＝p1p2(m
1-m2)2；
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0083]
步骤五四：改变设定的kl散度阈值，再计算类间方差，使类间方差最大的kl散度值即为最终的kl散度分割阈值kld，其公式如下：
[0084]
kld＝arg
max
σ2。
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0085]
其它组成及连接方式与具体实施方式三相同。
[0086]
具体实施方式五：参照图1至图4说明本实施方式，本实施方式是对具体实施方式一所述的步骤九中心定位轴组件8作进一步限定，本实施方式中所述步骤九中gosso-cfar检测流程如下：
[0087]
步骤九一：载入检测单元及其左右滑窗内的经过平方律检波后的数据，设检测单元的值为d，其左侧滑窗内的参考单元数据为x1,x2,x3......xn，右侧滑窗内的参考单元数据为y1,y2,y3......yn，其中n为设定的参考单元数量；
[0088]
步骤九二：对左侧的参考单元按由小到大的顺序排序，选出第k1个最小单元x，对右侧的参考单元按由小到大的顺序排序，选出第k2个最小单元y；
[0089]
步骤九三：根据步骤九一和步骤九二所得参数确定背景噪声的估计值z为：
[0090]
z＝min(x,y)；
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)
[0091]
步骤九四：检测阈值s＝tz，其中t为标称因子，其可由下式计算出来：
[0092][0093]
其中，p
fa
为虚警率，γ()为γ函数，对于整数x，γ(x)＝(x-1)！；
[0094]
步骤九五：根据计算得到的检测阈值，做出判定结果，设为h1有目标，h2为无目标，则判定策略如下：
[0095][0096]
即如果检测单元的值大于阈值，则判定为目标点，反之，则判定为非目标点。其它组成及连接方式与具体实施方式四相同。
[0097]
本发明已以较佳实施案例揭示如上，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可以利用上述揭示的结构及技术内容做出些许的更动或修饰为等同变化的等效实施案例，但是凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施案例所做的任何简单修改、等同变化与修饰，均仍属本发明技术方案范围。