基于海尖峰抑制和多帧联合的小目标检测方法技术

本发明专利技术公开了一种基于海尖峰抑制和多帧联合的小目标检测方法，主要解决现有技术不适用于在快速扫描模式下检测海面低、慢、小目标的问题。其实现步骤为：1)获取回波数据；2)将回波数据沿着脉冲维进行相干累计得到数据矩阵z；3)利用数据矩阵z，计算第一重检测门限T1并利用T1重新设置数据矩阵z；4)设计方向匹配滤波器，利用方向匹配滤波器和重新设置的数据矩阵z，计算待检测单元d的检验统计量ξd和第二重检测门限T；5)通过比较检验统计量ξd和第二重检测门限T的大小判断出目标是否存在。本发明专利技术提高了雷达检测性能，降低了虚警概率，可用于检测海面低速微弱目标。

【技术实现步骤摘要】
基于海尖峰抑制和多帧联合的小目标检测方法
本专利技术属于信号处理
，具体涉及一种小目标检测方法，可用于海面低、慢、小目标的识别与跟踪。
技术介绍
海杂波背景下的低速微弱目标检测方法研究在军事、民用方面均具有广阔的应用前景。岸基雷达和机载雷达通常工作在高分辨、快扫描模式下。其中快扫描模式下雷达在每个波位获得的脉冲数在10个以下。在距离分辨率较高、脉冲数较少的情况下，海面低速微弱运动目标检测的难度较高。而扫描内出现的海尖峰效应更加大了海面目标检测的难度。文献He,Y.,Guan,J.:Meng,X.W.etal.:‘RadartargetdetectionandCFARprocessing’,(TsinghuaUniversityPress,2011,2stedn.),pp.30-50和文献Watts,S.:‘Cell-averagingCFARgaininspatiallycorrelatedK-distributedclutter’,IEEERadarSonarNavig.,1996,143,pp.321-327提出的各类基于能量的恒虚警检测方法由于容易实现、计算速度快被广泛应用于雷达目标检测中。基于能量的恒虚警检测方法在海杂波统计特性未知、脉冲数较少的情况下仍可使用。但是，由于微弱目标的回波经常被淹没在强的海杂波中使得雷达获得的信杂比降低，在低信杂比情况下，基于能量的恒虚警检测方法的检测性能很差。而且，基于能量的恒虚警检测方法也不能减少海尖峰效应带来的大量虚警。
技术实现思路
本专利技术的目的在于提出一种基于海尖峰抑制和多帧联合的小目标检测方法，以提高在快扫描模式下对海面低、慢、小目标的检测性能，降低虚警概率。为实现上述目的，本专利技术的技术方案包括如下：(1)利用雷达发射机发射连续的脉冲信号，雷达接收机接收M×I×Q维的回波数据矩阵X，其中，M表示帧数即扫描数，I表示距离单元数，Q表示积累脉冲数；(2)利用回波数据矩阵X沿着脉冲维求非相干累积，计算数据矩阵z，其中，数据矩阵z的第m行i列的数据为z(m,i)：||表示取模，X(m,i,q)表示回波数据矩阵X的第m行i列q页的回波数据；(3)给定第一重检测的虚警概率p1和第二重检测的虚警概率p，通过蒙特卡罗实验计算第一重检测门限T1，其中，p1≥p；(4)根据第一重检测门限T1重新设置数据矩阵z，其中，重新设置的数据矩阵z的第m行i列的数据为：(5)给定窗长2L+1，设计方向匹配滤波器；(6)给定参考单元个数R，利用方向匹配滤波器和重新设置后的数据矩阵z，计算待检测单元d的检验统计量ξd；(7)根据第二重检测的虚警概率p，通过蒙特卡罗实验计算第二重检测门限T；(8)比较检验统计量ξd和第二重检测门限T的大小，判断出目标是否存在：如果ξd≥T，则表明待检测距离单元d有目标，如果ξd＜T，则表明待检测距离单元d没有目标。本专利技术与现有技术比较具有以下优点：1)本专利技术由于基于能量检测，相比已有的自适应检测方法更具有广泛性，能在海杂波统计特性未知、脉冲数较少的情况下对海面目标进行检测。2)本专利技术由于基于能量在方向上的累积且保证在所有方向均进行搜索，不仅利用了帧间的能量信息还利用了目标运动方向的信息，相比已有的依靠能量的恒虚警检测方法，更适合用于在快速扫描模式下对海面微弱目标检测。3)本专利技术由于使用双重检测，相比已有方法，能有效抑制海尖峰，提高了在快速扫描模式下雷达对海面低、慢、小目标的检测性能，降低了虚警概率。附图说明图1为本专利技术的实现流程图；图2为用本专利技术和现有方法在实测海杂波数据下进行目标检测的结果图；具体实施方式下面结合附图对本专利技术作进一步说明：参照图1，本专利技术的实现步骤如下：步骤1，获取回波数据。利用雷达发射机发射连续的脉冲信号，雷达接收机接收M×I×Q维的回波数据矩阵X，其中，M表示帧数即扫描数，I表示距离单元数，Q表示积累脉冲数。步骤2，利用回波数据矩阵X沿着脉冲维求非相干累积，计算数据矩阵z。求非相干累积计算数据矩阵z有两种方法：第一种是加性累积，其计算公式为:其中，z(m,i)表示数据矩阵z的第m行i列的数据，||表示取模，X(m,i,q)表示回波数据矩阵X的第m行i列q页的回波数据；第二种是乘性累积，其计算公式为：其中，z(m,i)表示数据矩阵z的第m行i列的数据，||表示取模，X(m,i,q)表示回波数据矩阵X的第m行i列q页的回波数据；本实例使用但不限于使用第二种方法计算数据矩阵z。步骤3，给定第一重检测的虚警概率p1和第二重检测的虚警概率p，通过蒙特卡罗实验计算第一重检测门限T1。(3.1)在数据矩阵z中选取V个纯杂波单元作为训练单元，V≥100/p，本实例取V＝100/p；(3.2)取每个纯杂波单元在数据矩阵z中对应的值作为检验统计量，将得到的V个检验统计量按降序排列，取排列后的第[Vp1]个检验统计量作为第一重检测门限T1，其中[Vp1]表示不超过实数Vp1的最大整数。步骤4，根据第一重检测门限T1，重新设置数据矩阵z，其中，重新设置的数据矩阵z的第m行i列的数据为：步骤5，设计方向匹配滤波器。(5.1)给定窗长2L+1，确定水平滤波器h，其中，L为正整数，水平滤波器h的第x行y列的值为：x＝-L,-L+1,...,L-1,L，y＝-L,-L+1,...,L-1,L，为获得更好的检测性能L≥4，本实例取L＝4；(5.2)利用水平滤波器h，保证所有方向都被搜索到，设计方向匹配滤波器：wj(θj)＝h(xcosθj-ysinθj,xsinθj+ycosθj)其中，wj(θj)表示第j个方向的方向匹配滤波器；θj表示第j个旋转角度，其为第j个方向的方向匹配滤波器的自变量；J表示方向匹配滤波器的个数，J≥2L，本实例取J＝16。步骤6，计算待检测单元d的检验统计量。(6.1)在重新设置后的数据矩阵z中，以待检测单元d为中心，选取(2L+1)×(2L+1)维数据为模版矩阵，将模版矩阵与方向匹配滤波器点乘并对点乘结果中的所有元素进行累积求和，得到待检测单元d在每个方向上的累积和；(6.2)给定参考单元个数R,选取待检测单元d周围的R个单元作为参考单元，利用R个参考单元分别替换步骤(5.1)中的待检测距离单元d，得到每个参考单元在每个方向上的累积和，再在每个方向上选取出R个累积和中的最大值；(6.3)利用待检测单元d在每个方向上的累积和与R个参考单元在每个方向上的累积和的最大值计算待检测单元d的检验统计量ξd：其中，max()表示取最大值，st表示待检测单元d在第t个方向上的累积和，nt表示R个参考单元在第t个方向上的累积和的最大值，t＝1,2,...,J。步骤7，根据系统给定的第二重检测虚警概率p，通过蒙特卡罗实验计算第二重检测门限T。(7.1)以步骤(3.1)中选取的V个纯杂波单元作为训练单元；(7.2)利用V个训练单元分别替换步骤6中的待检测距离单元d，重复步骤6得到每个训练单元的检验统计量；(7.3)将得到的V个检验统计量按降序排列，取排列后的第[Vp]个检验统计量作为检测门限T，其中[Vp]表示不超过实数Vp的最大整数。步骤8，通过比较检验统计量ξd和检测门限T的大小判断出目标是否存在：如果ξd≥T，则表明待检测距离单元d有目标，如本文档来自技高网...

【技术保护点】
一种基于海尖峰抑制和多帧联合的小目标检测方法，包括：(1)利用雷达发射机发射连续的脉冲信号，雷达接收机接收M×I×Q维的回波数据矩阵X，其中，M表示帧数即扫描数，I表示距离单元数，Q表示积累脉冲数；(2)利用回波数据矩阵X沿着脉冲维求非相干累积，计算数据矩阵z，其中，数据矩阵z的第m行i列的数据为z(m,i)：

【技术特征摘要】
1.一种基于海尖峰抑制和多帧联合的小目标检测方法，包括：(1)利用雷达发射机发射连续的脉冲信号，雷达接收机接收M×I×Q维的回波数据矩阵X，其中，M表示帧数即扫描数，I表示距离单元数，Q表示积累脉冲数；(2)利用回波数据矩阵X沿着脉冲维求非相干累积，计算数据矩阵z，其中，数据矩阵z的第m行i列的数据为z(m,i)：||表示取模，X(m,i,q)表示回波数据矩阵X的第m行i列q页的回波数据；(3)给定第一重检测的虚警概率p1和第二重检测的虚警概率p，通过蒙特卡罗实验计算第一重检测门限T1，其中，p1≥p；(4)根据第一重检测门限T1重新设置数据矩阵z，其中，重新设置的数据矩阵z的第m行i列的数据为：(5)给定窗长2L+1，设计方向匹配滤波器；(6)给定参考单元个数R，利用方向匹配滤波器和重新设置后的数据矩阵z，计算待检测单元d的检验统计量ξd；(7)根据第二重检测的虚警概率p，通过蒙特卡罗实验计算第二重检测门限T；(8)比较检验统计量ξd和第二重检测门限T的大小，判断出目标是否存在：如果ξd≥T，则表明待检测距离单元d有目标，如果ξd＜T，则表明待检测距离单元d没有目标。2.如权利要求1所述的方法，其特征在于,步骤(3)中计算第一重检测门限T1，按如下步骤进行：(3a)在数据矩阵z中选取V个纯杂波单元作为训练单元，V≥100/p；(3b)取每个纯杂波单元在数据矩阵z中对应的值作为检验统计量，将得到的V个检验统计量按降序排列，取排列后的第[Vp1]个检验统计量作为第一重检测门限T1，其中[Vp1]表示不超过实数Vp1的最大整数。3.如权利要求1所述的方法，其特征在于,步骤(5)中设计方向匹配滤波器，按如下步骤进行：(5a)给定窗长2L+1，确定水平滤波器h，其中，L为正整数，水平滤波器h的第x行y列的值为：(5b)利用水平滤波器h，保证所有方向都被搜索到，设计方向匹...

【专利技术属性】
技术研发人员：水鹏朗，杨春娇，施赛楠，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61