一种二维缩比快速恒虚警检测方法与流程

本发明涉及一种二维缩比快速恒虚警检测方法。

背景技术：

恒虚警率处理技术是雷达信号处理的重要组成部分，在雷达检测中占有不可或缺的地位，从原始回波到经过预处理的数据，如多普勒频谱，甚至是合成孔径雷达图像都需要雷达检测。在实际工作中，雷达面临的大空域范围快速搜索目标、检测时间长等问题，使传统的恒虚警检测算法难以保持良好的检测性能。

目前主要应用的恒虚警检测包括单元平均(ca_cfar)恒虚警检测、单元平均选大(go_cfar)、单元平均选小(so_cfar)和有序统计量(os_cfar)恒虚警检测。可查到的相关专利有“应用在恒虚警检测及数据处理中的单元平均选小门限检测方法”(申请号：20160729115.3)、“基于单边检测单元累积平均的多目标恒虚警检测方法”(申请号：201610438762.9)、“二维自适应雷达恒虚警检测方法”(申请号：201310240053.6)。ca_cfar在非均匀杂波背景中性能严重下降，且对数据特性要求较高；go_cfar在多目标环境中的检测性能下降到非常严重；so_cfar在均匀杂波背景中的检测性能和边缘杂波中虚警性能都很差；os_cfar是在对所有参考单元的测量值进行有序统计排序的基础上，选取最合适的第k个参考单元进行背景杂波功率估计，它在抗脉冲干扰方面作用显著，但k值选择不合理时也会引起漏检或者虚警，同时需要进行排序处理。以上专利主要采取上述几种恒虚警检测方法，检测时间长，硬件开销大，不易于工程实现。

此外，相关公开文献包括《利用二维恒虚警进行非均匀噪声背景下的目标检测》是针对高频地波雷达中的非均匀噪声环境及各种干扰源时变的特点，先对非均匀噪声环境进行白化处理，然后在一维恒虚警的基础上，采用一种基于多普勒频率和距离的二维恒虚警处理算法同时进行检测；《雷达在不同杂波背景下的恒虚警率处理》主要论述符合高斯分布的噪声和杂波背景下，几种典型恒虚警检测器的优缺点以及适用环境，并结合工程实践经验给出各类恒虚警门限的优化设计原理，提出了在不同环境下的最优恒虚警检测处理流程。以上文章和专利都是针对现有恒虚警检测方法进行对比，均未提及降低雷达二维恒虚警检测时间。

技术实现要素：

本发明提供一种二维缩比快速恒虚警检测方法，既能在不同杂波背景中保持较好的检测性能，又能大大减少检测时间，实现了雷达大空域范围快速搜索目标的同时提高了工程实现的可能性。

为了达到上述目的，本发明提供一种二维缩比快速恒虚警检测方法，包含以下步骤：

s1、对假定的雷达二维检测平面进行等比例分块，得到多个子块：

假定雷达二维检测平面的距离维长度为r，频率维长度为f，对该雷达二维检测平面进行等比例分块，得到n个大小为r×f的子块，其中，β为缩比因子，r、f、r、f均不等于0。

s2、计算雷达二维检测平面各个子块的点数均值：

将雷达二维检测平面的n个子块中各点数分别相加，得到n1、n2、...、nn，并计算各个子块的点数均值，记为其中k＝1,2,...,n。

s3、计算雷达二维检测平面的二维恒虚警检测门限：

将雷达二维检测平面的距离维等分为四部分，分别为a1、a2、a3以及a4；

将雷达二维检测平面的频率维等分为四部分，分别为b1、b2、b3以及b4；

将得到的16块沿距离维依次命名t1、t2、...、t16，将t1、t4、t13以及t16四块内各点数相加分别得到t1、t4、t13以及t16；

比较t1、t4、t13以及t16各值的大小，将其中最大值记为tmax，并通过计算得到二维恒虚警检测门限t：

s4、利用二维恒虚警检测门限对雷达二维检测平面的各个子块进行遍历筛选，得到点数均值大于门限的子块记为有目标子块：

将各个子块的均值dk分别与二维恒虚警检测门限t进行比较，若均值dk大于检测门限t，则判断该子块为有目标子块，记为mk，其中k＝1,2,...,n。

s5、计算有目标子块的二维恒虚警检测门限：

将各个有目标子块mk的距离维等分为四部分，分别为c1、c2、c3以及c4；

将雷达二维检测平面的频率维等分为四部分，分别为h1、h2、h3以及h4；

将得到的16块沿距离维依次命名q1、q2、...、q16，将q1、q4、q13以及q16四块内各点数相加分别得到q1、q4、q13以及q16；

比较q1、q4、q13以及q16各值的大小，将其中最大值记为qmax，并通过计算得到二维恒虚警检测门限pk：

s6、利用有目标子块的二维恒虚警检测门限对有目标子块进行遍历筛选，得到大于门限的点记为目标，并获得目标位置：

根据二维恒虚警检测门限pk，遍历有目标子块mk中各个点mk；

若有mk(x，y)大于等于检测门限pk，则判定该点为目标，并记录此时的目标位置(βx,βy)，其中1<＝x<＝r，1<＝y<＝f；

若有mk(x，y)小于检测门限pk，则判定该点不是目标。

本发明与背景技术相比具有如下优点：

1、从雷达检测时间出发，引入二维缩比恒虚警检测方法，大大减少雷达检测时间，实现了雷达大空域范围快速搜索目标；

2、在不同杂波背景中能保持较好的检测性能，提高了工程实现的可能性。

附图说明

图1是本发明提供的一种二维缩比快速恒虚警检测方法的流程图。

图2是恒虚警检测时间对比图。

具体实施方式

以下根据图1和图2具体说明本发明的较佳实施例。

如图1所示，本发明提供一种二维缩比快速恒虚警检测方法，包含以下步骤：

步骤s1、对假定的雷达二维检测平面进行等比例分块，得到多个子块；

步骤s2、计算雷达二维检测平面各个子块的点数均值；

步骤s3、计算雷达二维检测平面的二维恒虚警检测门限；

步骤s4、利用二维恒虚警检测门限对雷达二维检测平面的各个子块进行遍历筛选，得到点数均值大于门限的子块记为有目标子块；

步骤s5、计算有目标子块的二维恒虚警检测门限；

步骤s6、利用有目标子块的二维恒虚警检测门限对有目标子块进行遍历筛选，得到大于门限的点记为目标，并获得目标位置。

进一步，所述的步骤s1中，假定雷达二维检测平面的距离维长度为r，频率维长度为f，对该雷达二维检测平面进行等比例分块，得到n个大小为r×f的子块，其中，β为缩比因子，r、f、r、f均不等于0。

所述的步骤s2中，将雷达二维检测平面的n个子块中各点数分别相加，得到n1、n2、...、nn，并计算各个子块的点数均值，记为其中k＝1,2,...,n。

所述的步骤s3包含以下步骤：

步骤s3.1、将雷达二维检测平面的距离维等分为四部分，分别为a1、a2、a3以及a4；

步骤s3.2、将雷达二维检测平面的频率维等分为四部分，分别为b1、b2、b3以及b4；

步骤s3.3、将得到的16块沿距离维依次命名t1、t2、...、t16，将t1、t4、t13以及t16四块内各点数相加分别得到t1、t4、t13以及t16；

步骤s3.4、比较t1、t4、t13以及t16各值的大小，将其中最大值记为tmax，并通过计算得到二维恒虚警检测门限t：

所述的步骤s4中，将各个子块的均值dk分别与二维恒虚警检测门限t进行比较，若均值dk大于检测门限t，则判断该子块为有目标子块，记为mk，其中k＝1,2,...,n。

所述的步骤s5包含以下步骤：

步骤s5.1、将各个有目标子块mk的距离维等分为四部分，分别为c1、c2、c3以及c4；

步骤s5.2、将雷达二维检测平面的频率维等分为四部分，分别为h1、h2、h3以及h4；

步骤s5.3、将得到的16块沿距离维依次命名q1、q2、...、q16，将q1、q4、q13以及q16四块内各点数相加分别得到q1、q4、q13以及q16；

步骤s5.4、比较q1、q4、q13以及q16各值的大小，将其中最大值记为qmax，并通过计算得到二维恒虚警检测门限pk：

所述的步骤s6中，根据二维恒虚警检测门限pk，遍历有目标子块mk中各个点mk；

若有mk(x，y)大于等于检测门限pk，则判定该点为目标，并记录此时的目标位置(βx,βy)，其中1<＝x<＝r，1<＝y<＝f；

若有mk(x，y)小于检测门限pk，则判定该点不是目标。

图2为恒虚警检测时间对比图，数据来源于实测数据，选取截取后点数分别176，266，353，533，1100，1316，1433，1933，2600，可以看出二维缩比快速恒虚警检测方法可以大大减少雷达检测时间，实现了雷达大空域范围快速搜索目标。仿真结果验证了本发明的可行性。

本发明既能在不同杂波背景中保持较好的检测性能，又能大大减少检测时间，实现了雷达大空域范围快速搜索目标的同时提高了工程实现的可能性。

尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。