本发明涉及一种降低目标虚警的滤波技术,属于雷达信号处理技术领域。
背景技术:
机场监视雷达(asr:airportsurveillanceradar)属于航空管制一次雷达,通常设置在航管控制中心或相应的航路节点上,作用距离在200公里内,主要为塔台管制员提供机场附近飞机目标的进出港信息。
随着生活水平的日益提高,人们使用飞机出行的频率也在逐步上升。为了提高航班进出港的效率,需要机场监视雷达提供更为精细的飞机目标位置点迹信息。然而,传统的点迹滤波技术难以实现机场监视雷达对复杂多变的空间环境中飞机目标点迹的精细化处理。
为此,有学者提出了航管雷达的精细化点迹处理方法,例如,《雷达科学与技术》近期公开发表了一篇名为“航管一次雷达精细化点迹处理设计与实现”的文章,其主要方法是对雷达回波信号进行脉压、动目标检测(mtd)滤波、恒虚警率检测,然后对检测出的可疑点迹进行点迹凝聚和四级点迹滤波,最后将满足条件的点迹送给终端。该方法的亮点在于点迹滤波部分采用了分级有针对性地清除不同类型的虚假点迹,获得一定的点迹滤除效果,但是从最终的处理结果来看,剩余的点迹质量还有待进一步提高,同时这种多级滤波的策略增加了真实目标点迹被判为虚假点迹的可能性。
技术实现要素:
针对以上问题,本发明从图像处理的角度提供了一种全新的点迹滤波方法。该方法是基于多帧cfar检测的二值图像中飞机目标和干扰信号的图像特征具有明显差异来设计的,其中飞机目标的连通域表现为比较规整的点阵状,而干扰信号的连通域表现为无规则的孤立散点,因此可以通过判断多帧cfar检测的二值图像中连通域的点阵面积大小来实现飞机目标和干扰信号的分离的目的。为此,本发明创新地提出了基于多帧cfar检测二值图像的连通域面积测量的点迹滤波技术,解决了输出点迹质量不佳和多级判决容易造成真实目标点迹丢失的问题,达到了航管一次雷达精细化点迹滤波的目的。
为了解决以上问题,本发明采用了如下技术方案:一种基于二值图像连通域面积测量的点迹滤波技术,其特征在于,包括以下步骤:
步骤1:对单帧多脉冲回波信号进行脉压、动目标显示(mti)处理;
步骤2:对脉压、mti后的单帧距离向信号进行cfar检测;
步骤3:对取模值后的单帧信号进行滑窗cfar检测;
步骤4:对多帧cfar检测的二值图像进行连通域标记及面积测量;
步骤5:对满足条件的连通域进行重心位置计算;
步骤6:将检出目标的重心位置转化为雷达坐标系下的距离-方位信息;
步骤7:得到雷达坐标系下的点迹输出。
所述步骤1具体包括,对单帧多脉冲回波信号进行脉压、动目标显示(mti)处理:
设脉压后的第i帧多脉冲信号xi为
其中,
为了抑制零频附近的杂波,对脉压后的单帧多脉冲信号xi进行mti处理得单帧距离信号
其中,wm为mti加权向量w中第m个加权因子,0<wm<1,m=1,2,…,m;上标h表示共轭转置算子;上标*表示共轭算子。
所述步骤2具体包括,对脉压、mti后的单帧距离向信号进行cfar检测:
对mti处理后的单帧距离信号yi进行取模值操作,如下式所示:
其中,
所述步骤3具体包括,对取模值后的单帧信号进行滑窗cfar检测,具体过程可分解为如下四个步骤:
a)设置检测希望的虚警概率pf,检测窗口长度nw(奇数)和保护窗口长度np(单边);
b)估计当前滑窗内的噪声功率σ2(n),n=1,2,…,n,下面分三种情况计算噪声功率:
(i)如果
(ii)如果1≤n≤(nw-1)/2,则
(iii)如果n-(nw-1)/2+1≤n≤n,则
c)根据所设定的虚警概率计算检测门限,如下式所示
式中,ln(pf)表示对pf取自然对数。
d)根据所得检测门限dth(n),利用平方律检波器hp(i,n)来判断待检测单元
所述步骤4具体包括,对连续ibw(2≤ibw≤i)帧检测结果的i×n矩阵hp进行连通域标记并统计每个连通域的面积大小,具体过程如下:
a)对hp进行逐行逐列扫描,并对非零连通区域进行标记,设求得连通区域个数记为znum,这里0<znum<500;
b)逐个连通区域统计其面积大小,即计算该连通域内像素的总个数,假定第l个连通区的面积为zarea(l),其中zarea为znum×1维的列向量,l=1,2,…,znum,这里0<zarea(l)<1000;
所述步骤5具体包括,判断第l个连通区的面积zarea(l)是否大于阈值tarea,这里取tarea=4。若不大于,则不处理,否则检出目标数增加1,并计算该连通域的重心位置(xc(l),yc(l)),即
式中,xc,yc分别为znum×1维列向量,xl,yl分别表示第l个连通区内所有像素位置坐标的横坐标和纵坐标列向量,这里0<xc(l),xl(u)<6000,0<yc(l),yl(u)<484。
所述步骤6具体包括,将检出目标的重心位置转化为雷达坐标系下的距离-方位信息,即
其中,r,a分别为znum×1维距离和方位列向量,r0和a0分别为距离和方位的起始位置,△r和△a分别为距离和方位的间隔长度,这里r0=3km,a0=0°,△r=15m,△a=0.74°。
为了验证本发明的实际工程应用价值和有益效果,我们在同等条件下比较了本发明方法和传统方法处理某航管一次雷达11圈扫描数据的点迹输出效果,见附图8和图9。
数据处理中选用的参数如下:(1)恒虚警率为10-6;(2)cfar检测窗口长度为21;(3)保护单元长度为1;(4)cfar检测结果二值图大小为484x6000,其中484为帧数,6000为距离单元数;(5)连通域标记选用8邻接方式;(6)连通域面积阈值为4。
处理结果表明:本发明方法能有效滤除虚假目标点迹,虚假点迹抑制率约为96.5%(其中传统方法共计检出5574个点迹,本发明方法共计检出197个点迹,这里定义虚假点迹抑制率为:1-197/5574=0.9647),能够满足航管一次雷达精细化点迹滤波的目的。
附图说明
图1(a)为mti滤波过程示意图(mti滤波前频谱及mti频率响应);
图1(b)为mti滤波过程示意图(mti滤波后频谱);
图2为cfar滑窗检测示意图;
图3(a)为4邻接连通域示意图;
图3(b)为8邻接连通域示意图;
图4为某二值图像(8x11);
图5为连通域标记结果;
图6为本发明技术方案流程图;
图7为雷达天线波束中心扫描示意图;
图8为传统方法点迹输出结果;
图9为本发明方法点迹输出结果;
图10为连通域面积统计结果。
具体实施方式
下面结合对本发明做进一步阐述。
为了更好地说明本发明,先简单介绍一下本发明的方案设计理论。
1、运动目标指示
运动目标指示(mti)是多普勒处理方法中较为常见的一种简单快速的滤波手段,其主要思想是通过构造一个高通滤波器,将零频附近的杂波能量抑制在一定水平内,从而提高运动目标的信杂比,有利于目标的检测。mti滤波过程示意图见图1(a)、图1(b)。
需要指出的是,mti技术一般用于脉冲数较少的场合,即滤波器的阶数受限,这直接导致了滤波器频率响应的凹口较宽,在抑制杂波的同时,也会对低速运动目标产生一定衰减。
常用的m脉冲mti滤波器传递函数可由二脉冲对消器级联获得,即
对应的冲击响应系数为
式中,(·)!表示排列算子。
2、恒虚警率检测
恒虚警率(cfar)检测是一种在复杂背景环境中提供可预测虚警率的检测技术。假设我们这里讨论的干扰为高斯白噪声,对于仅存在一个采样且采样值是归一化的,则其虚警率可由下式给出
其中,dth为检测门限。由上式我们不难可求得在给定虚警率pf下所需的检测门限值,即
dth=-lnpf(14)
对于非归一化的采样情况,检测门限的计算需要考虑噪声功率σ2的影响,即
dth=-σ2lnpf(15)
由上式不难看出,准确计算检测门限除了给定的虚警率pf外,还需要知道噪声功率σ2。对于大多数雷达系统而言,该参数是未知的,或者说是会随着时间和场景的变化而发生变化的,因此需要对其进行估计。而噪声估计的关键在于噪声样本的挑选。下图给出了一种常用于对距离单元进行处理的cfar滑窗检测示意图见图2,其中待检测单元位于检测窗口中央,其相邻两侧有一定长度的保护单元,剩余部分为用于噪声统计的参考样本数据。增加保护单元是因为目标在距离向上可能跨多个单元,此时待检测单元临近的样本内不仅含有噪声,还有目标的部分能量。
3、二值图像连通域标记及面积测量
连通域标记是指从二值图像(0和1)中将邻接的“1”集合提取出来,并为其标记相应的数字序号,并将该数字序号作为相应连通域集合内所有元素的新像素值。在图像处理领域,该操作是最为基本的,也是非常关键的一步,对图像处理的最终结果具有重要影响。
连通域面积测量是指在完成二值图像连通域标记后,对连通域内所含元素对应尺度面积的总和,对于尺度均匀且归一化的二值图像,连通域面积等于所含元素个数总和。
连通域标记的结果跟所选用的邻接方式有关。常见的邻接方式有两种:4邻接和8邻接。4邻接一共有4个点,即上下左右,见图3(a)。8邻接的点一共有8个,包括了对角线位置的点,见图3(b)。
下面我们举个二值图像连通域标记过程的例子。假设某二值图像(8x11),见图4,图中含有5个连通域。下面采用8邻接方式对该二值图像进行从左至右、从上之下逐个像素进行扫描,最终得到5个连通域标记结果,见图5。假设该二值图像的尺度是均匀且归一化的,则连通域面积统计结果见图10。
结合上述理论,如图6所示,本发明提供了一种基于二值图像连通域面积测量的点迹滤波技术,其特征在于,包括以下步骤:
步骤1:对单帧多脉冲回波信号进行脉压、动目标显示(mti)处理;
步骤2:对脉压、mti后的单帧距离向信号进行cfar检测;
步骤3:对取模值后的单帧信号进行滑窗cfar检测;
步骤4:对多帧cfar检测的二值图像进行连通域标记及面积测量;
步骤5:对满足条件的连通域进行重心位置计算;
步骤6:将检出目标的重心位置转化为雷达坐标系下的距离-方位信息;
步骤7:得到雷达坐标系下的点迹输出。
所述步骤1具体包括,对单帧多脉冲回波信号进行脉压、动目标显示(mti)处理:
设脉压后的第i帧多脉冲信号xi为
其中,
为了抑制零频附近的杂波,对脉压后的单帧多脉冲信号xi进行mti处理得单帧距离信号
其中,wm为mti加权向量w中第m个加权因子,0<wm<1,m=1,2,…,m;上标h表示共轭转置算子;上标*表示共轭算子。
所述步骤2具体包括,对脉压、mti后的单帧距离向信号进行cfar检测:
对mti处理后的单帧距离信号yi进行取模值操作,如下式所示:
其中,
所述步骤3具体包括,对取模值后的单帧信号进行滑窗cfar检测,具体过程可分解为如下四个步骤:
a)设置检测希望的虚警概率pf,检测窗口长度nw(奇数)和保护窗口长度np(单边),本发明取pf=10-6,nw=21,np=1;
b)估计当前滑窗内的噪声功率σ2(n),n=1,2,…,n,下面分三种情况计算噪声功率:
(i)如果
(ii)如果1≤n≤(nw-1)/2,则
(iii)如果n-(nw-1)/2+1≤n≤n,则
c)根据所设定的虚警概率计算检测门限,如下式所示
式中,ln(pf)表示对pf取自然对数。
d)根据所得检测门限dth(n),利用平方律检波器hp(i,n)来判断待检测单元
所述步骤4具体包括,对连续ibw(2≤ibw≤i)帧检测结果的i×n矩阵hp进行连通域标记并统计每个连通域的面积大小,具体过程如下:
a)对hp进行逐行逐列扫描,并对非零连通区域进行标记,假定求得连通区域个数记为znum,这里0<znum<500;
b)逐个连通区域统计其面积大小,即计算该连通域内像素的总个数,假定第l个连通区的面积为zarea(l),其中zarea为znum×1维的列向量,l=1,2,…,znum,这里0<zarea(l)<1000;
所述步骤5具体包括,判断第l个连通区的面积zarea(l)是否大于阈值tarea,这里取tarea=4。若不大于,则不处理,否则检出目标数增加1,并计算该连通域的重心位置(xc(l),yc(l)),即
式中,xc,yc分别为znum×1维列向量,xl,yl分别表示第l个连通区内所有像素位置坐标的横坐标和纵坐标列向量,这里0<xc(l),xl(u)<6000,0<yc(l),yl(u)<484。
所述步骤6具体包括,将检出目标的重心位置转化为雷达坐标系下的距离-方位信息,即
其中,r,a分别为znum×1维距离和方位列向量,r0和a0分别为距离和方位的起始位置,△r和△a分别为距离和方位的间隔长度,这里r0=3km,a0=0°,△r=15m,△a=0.74°。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。