基于积分图像的二维CFAR门限计算方法技术

本发明专利技术涉及一种基于积分图像的二维恒虚警门限计算方法，属于雷达技术领域。包括构造积分图像矩阵；根据给定的参考窗和保护窗计算出参考窗与保护窗对应顶点坐标位置，由顶点坐标位置可知积分图像矩阵对应点的元素值；计算参考窗和保护窗中所有元素的累加和；计算参考窗和保护窗中所含元素的总个数；然后计算回形窗区域杂波幅度均值和；根据回形窗区域杂波幅度均值和计算CFAR门限阈值。与传统方法相比较，计算次数与给定的参考窗、保护窗大小无关，且总次数少于传统二维恒虚警门限计算方法，大大缩减了算法耗时，可以满足工程实时性要求，应用性更强。应用性更强。应用性更强。

【技术实现步骤摘要】
基于积分图像的二维CFAR门限计算方法


[0001]本专利技术属于雷达
，涉及雷达的恒虚警检测算法，具体涉及一种新的基于积分图像的二维恒虚警门限阈值计算方法。

技术介绍

[0002]恒虚警目标检测技术就是当杂波的幅值发生变化时，自适应调节雷达的检测门限，使虚警概率保持恒定。采用恒虚警检测技术能够保证信号处理机不因干扰太强而过载，从而使雷达目标检测系统正常工作。因此采用恒虚警检测技术意义重大，在雷达目标检测中，必须采用恒虚警检测技术。CFAR检测技术的具体实现形式是CFAR检测器。传统的雷达目标检测工程实现时，多采用均值类CFAR检测，包括CA
‑
CFAR、GO
‑
CFAR、SO
‑
CFAR等。图1为均值类CFAR检测的示意图。
[0003]研究发现，雷达背景中的噪声和杂波是在距离和多普勒域二维同时存在，并且局部噪声在剔除异常值后，可以认为是平稳的高斯随机过程。显然，二维CFAR检测利用待检单元附近距离
‑
多普勒二维空间上的单元作为参考滑窗要比传统的一维距离域或多普勒频率域参考滑窗使CFAR检测处理更趋合理有效。图2为二维CFAR检测的示意图。
[0004]积分图像是Viola等人于2001年提出的，他指出通过构造积分图像，矩形区域内的像素和的计算可通过三次加法完成，利用这一点可加快检测的速度。
[0005]对于一个数组对于一个数组a(x)，其前缀加法输出结果b(x)可表示为
[0006][0007]设有一幅图像I(x,y)，那么积分图像可以表示为
[0008][0009]如图3所示，积分图像在(x,y)处的值为原始图像在由(0,0)与(x,y)确定的矩形区域内的像素和。

技术实现思路

[0010]要解决的技术问题
[0011]本专利技术要解决的技术问题是：如图2所示，传统的二维CFAR门限阈值计算方法是给定合理的参考窗和保护窗后，利用累加法分别求出待测目标对应的参考窗和保护窗中所有元素的和，再用参考窗中元素和减去保护窗中元素和的差除以对应回形窗区域内元素个数，可以得到回形窗区域杂波幅度均值，均值乘以门限因子即为二维CFAR门限阈值。用这种方法求解门限阈值要进行的加法次数与所选参考窗、保护窗大小相关，且重复计算次数多、计算量大、耗时长，工程实用性不强。
[0012]技术方案
[0013]一种基于积分图像的二维CFAR门限计算方法，其特征在于步骤如下：
[0014]步骤1：已知原幅度模值矩阵为I，矩阵行数为dis，列数为dpl，选定距离维参考单元数为ref_dis，保护单元数为prt_dis，方位维参考单元数为ref_dpl，保护单元数为prt_dpl；构建全0矩阵Amp、Num，矩阵行数均为dis+2*(ref_dis+prt_dis)，列数均为dpl+2*(ref_dpl+prt_dpl)，为矩阵赋值如下：
[0015]Amp(ref_dis+prt_dis+1:dis+ref_dis+prt_dis,ref_dpl+prt_dpl+1:dpl+ref_dpl+prt_dpl)＝I
[0016]Num(ref_dis+prt_dis+1:dis+ref_dis+prt_dis,ref_dpl+prt_dpl+1:dpl+ref_dpl+prt_dpl)＝1
[0017]步骤2：对上面构建的两个矩阵分别按行、按列进行求和操作，分别得到幅度和矩阵SumAmp、点数和矩阵SumNum如下：
[0018][0019][0020]其中，0＜m≤dis+2*(ref_dis+prt_dis),0＜n≤dpl+2*(ref_dpl+prt_dpl)；
[0021]步骤3：设待检测点O坐标为(x,y)，O点对应参考窗为ABCD，保护窗为abcd，求出A、B、C、D，a、b、c、d坐标A(x+ref_dis+prt_dis,y+ref_dpl+prt_dpl)，B(x+ref_dis+prt_dis,y
‑
ref_dpl
‑
prt_dpl)，C(x
‑
ref_dis
‑
prt_dis,y
‑
ref_dpl
‑
prt_dpl)，D(x
‑
ref_dis
‑
prt_dis,y+ref_dpl+prt_dpl)，a(x+prt_dis,y+prt_dpl)，b(x+prt_dis,y
‑
prt_dpl)，c(x
‑
prt_dis,y
‑
prt_dpl)，d(x
‑
prt_dis,y+prt_dpl)，分别求参考窗、保护窗内所有元素的累加和Sum
ABCD
、Sum
abcd
及元素个数Sum
NABCD
、Sum
Nabcd
：
[0022]Sum
ABCD
＝SumAmp(m
C
,n
C
)
‑
SumAmp(m
B
,n
B
)
‑
SumAmp(m
D
,n
D
)+SumAmp(m
A
,n
A
)
[0023]Sum
abcd
＝SumAmp(m
c
,n
c
)
‑
SumAmp(m
b
,n
b
)
‑
SumAmp(m
d
,n
d
)+SumAmp(m
a
,n
a
)
[0024]Sum
NABCD
＝SumNum(m
C
,n
C
)
‑
SumNum(m
B
,n
B
)
‑
SumNum(m
D
,n
D
)+SumNum(m
A
,n
A
)
[0025]Sum
Nabcd
＝SumNum(m
c
,n
c
)
‑
SumNum(m
b
,n
b
)
‑
SumNum(m
d
,n
d
)+SumNum(m
a
,n
a
)
[0026]步骤4：计算由ABCD和abcd构成的回形窗区域杂波幅度均值为：
[0027][0028]步骤5：设恒虚警的门限阀值为k，门限因子为factor，可得
[0029][0030]一种计算机系统，其本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于积分图像的二维CFAR门限计算方法，其特征在于步骤如下：步骤1：已知原幅度模值矩阵为I，矩阵行数为dis，列数为dpl，选定距离维参考单元数为ref_dis，保护单元数为prt_dis，方位维参考单元数为ref_dpl，保护单元数为prt_dpl；构建全0矩阵Amp、Num，矩阵行数均为dis+2*(ref_dis+prt_dis)，列数均为dpl+2*(ref_dpl+prt_dpl)，为矩阵赋值如下：Amp(ref_dis+prt_dis+1：dis+ref_dis+prt_dis，ref_dpl+prt_dpl+1：dpl+ref_dpl+prt_dpl)＝INum(ref_dis+prt_dis+1：dis+ref_dis+prt_dis，ref_dpl+prt_dpl+1：dpl+ref_dpl+prt_dpl)＝1步骤2：对上面构建的两个矩阵分别按行、按列进行求和操作，分别得到幅度和矩阵SumAmp、点数和矩阵SumNum如下：点数和矩阵SumNum如下：其中，0＜m≤dis+2*(ref_dis+prt_dis)，0＜n≤dpl+2*(ref_dpl+prt_dpl)；步骤3：设待检测点O坐标为(x，y)，O点对应参考窗为ABCD，保护窗为abcd，求出A、B、C、D，a、b、c、d坐标A(x+ref_dis+prt_dis，y+ref_dpl+prt_dpl)，B(x+ref_dis+prt_dis，y
‑
ref_dpl
‑
prt_dpl)，C(x
‑
ref_dis
‑
prt_dis，y
‑
ref_dpl
‑
prt_dpl)，D(x
‑
ref_dis
‑
prt_dis，y+ref_dpl+prt_dpl)，a(x+prt_dis，y+prt_dpl)，b(x+prt_dis，y
‑
prt_dpl)，c(x
‑
prt_dis，y
‑
prt_dpl)，d(x
‑
prt_dis，t+prt_dpl)，分别求参考窗、保护窗内所有元素的累加和Sum
ABCD
、Sum
abcd
及元素个数Sum
NABCD
、Sum
Nabcd
：Sum
ABCD
＝SumA...

【专利技术属性】
技术研发人员：汪慧，冯超，任泽宇，戴巧娜，吕腾蛟，李静静，
申请(专利权)人：西安电子工程研究所，
类型：发明
国别省市：