基于统计熵的海面目标检测方法

1.本发明涉及一种本发明属于雷达信号处理技术领域，具体涉及一种基于统计熵的海面目标检测方法。


背景技术：

2.海杂波是雷达照射下海面产生的后向散射回波，海杂波的存在会严重影响雷达对海面目标的检测能力。海杂波的形成机理十分复杂，其特性不仅取决于当前海面环境参数，而且还受到雷达工作参数的影响。
3.传统恒虚警率(cfar)目标检测算法基于海杂波的统计特性，利用数理统计模型对海杂波进行统计建模，并计算出相应的检测门限实现目标检测。对于海杂波模型，瑞利分布、对数正态分布和韦布尔分布首先被用来描述海杂波的幅度统计特性，之后考虑到海杂波的时空相关性和拖尾特性，又提出了用k分布和pareto分布去拟合海杂波的统计特征。然而，当海况较高时，大量海尖峰和异常散射单元的出现可能会导致海杂波模型的失配，从而导致目标检测性能的下降。


技术实现要素：

4.要解决的技术问题
5.本发明克服了传统cfar目标检测算法在海杂波幅度模型不匹配时，目标检测能力下降的不足，利用熵的特征来描述目标单元和海杂波单元之间海杂波幅度直方图扩展程度不同的特征差异，提出一种基于海杂波幅度统计熵的目标检测算法，提高了海杂波背景下的目标检测性能。
6.技术方案
7.一种基于统计熵的海面目标检测方法，其特征在于步骤如下：
8.步骤1：海杂波序列归一化处理
9.设海杂波幅度序列为：x＝{xk,k＝1,2,3..n}，其中n表示海杂波幅度序列的长度，利用海杂波幅度序列的标准差进行归一化处理，得到新的序列x'：
10.x'＝x/std(x)
11.其中，std(x)为海杂波幅度序列的标准差；
12.步骤2：海杂波幅度直方图计算
13.步骤2.1：计算海杂波幅度直方图中统计量的数目
14.针对步骤1中得到的归一化海杂波序列x'，利用如下公式计算海杂波幅度直方图中所有统计量的数目：
[0015][0016]
其中lx表示所有统计量的数目，fix(
·
)表示向下取整，m、n表示任意正整数；
[0017]
步骤2.2：计算各统计量的幅度
[0018]
针对步骤2.1中计算的海杂波幅度直方图中统计量的数目，根据海杂波序列的幅
度，将海杂波序列放置在每一个所属的统计分量内，计算每个统计量内海杂波序列的数目nn，即各统计量的幅度被存储在向量nn中，计算步骤如下：
[0019]
d＝(max(x')-min(x'))/lx
[0020]
(id,(i 1)d),i＝0,1...,lx-1
[0021]
其中d表示各统计量的宽度，即各统计量中所包含海杂波序列幅度的最大值减最小值，(id,(i 1)d),i＝0,1
…
,lx-1表示第i个统计量中所能包含海杂波序列的幅度区间；针对上述计算结果绘制海杂波幅度统计直方图；
[0022]
步骤3：计算海杂波模型幅度统计熵
[0023]
针对步骤2中得到的统计直方图，计算海杂波模型幅度统计熵；构造海杂波模型幅度统计熵来描述目标距离单元和海杂波距离单元之间统计直方图的差异；所述海杂波模型幅度统计熵的计算过程如下：
[0024]
f＝nn/n
[0025][0026]
其中f＝nn/n表示将频数nn转换成频率f，e表示海杂波模型幅度统计熵的计算结果，f(i)表示第i个频率分量；
[0027]
步骤4：目标探测
[0028]
针对步骤3中计算得到的幅度统计熵e的结果作为特征检测器的输入进行目标检测，特征检测器的结构如下所示：
[0029][0030]
其中z表示判决门限，h0表示海杂波距离单元，h1表示目标距离单元。
[0031]
一种计算机系统，其特征在于包括：一个或多个处理器，计算机可读存储介质，用于存储一个或多个程序，其中，当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现上述的方法。
[0032]
一种计算机可读存储介质，其特征在于存储有计算机可执行指令，所述指令在被执行时用于实现上述的方法。
[0033]
有益效果
[0034]
本发明提出的一种基于统计熵的海面目标检测方法，该方法可应用于岸基/弹载/机载对海探测雷达，通过定义幅度统计熵来描述目标与杂波距离单元回波信号统计直方图的特征差异，并利用此特征差异进一步实现目标检测，提升了雷达对海面小目标的检测能力。具有以下有益效果：
[0035]
1、本发明通过计算海杂波序列幅度统计直方图，通过分析目标距离单元与海杂波距离单元回波信号统计直方图，提取出幅度统计熵特征作为目标检测器的输入。由于幅度统计熵特征能够充分描述海杂波与目标的差异，因此，提出基于幅度统计熵的检测方法能够提升目标检测性能。
[0036]
2、由于基于幅度统计熵的目标检测算法不依赖于海杂波幅度模型，因此克服了传统cfar目标检测方法因模型失配所导致的目标检测能力下降的不足，提高了雷达对海上小
目标的检测性能。
附图说明
[0037]
附图仅用于示出具体实施例的目的，而并不认为是对本发明的限制，在整个附图中，相同的参考符号表示相同的部件。
[0038]
图1是本发明的算法实现流程图；
[0039]
图2是x波段雷达回波的距离-时间-强度图像；
[0040]
图3是目标距离单元雷达回波幅度统计直方图和纯海杂波距离单元海杂波幅度：a、海杂波距离单元；b、目标距离单元；
[0041]
图4是基于x波段雷达实测回波和海杂波幅度统计熵的目标检测结果。
具体实施方式
[0042]
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图和实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。此外，下面描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0043]
参照图1，本发明的具体实现过程如下：
[0044]
步骤1、海杂波幅度序列归一化处理。设海杂波幅度序列为：x＝{xk,k＝1,2,3..n}利用海杂波幅度序列的标准差进行归一化处理，得到新的序列x'。
[0045]
x'＝x/std(x)
[0046]
步骤2、海杂波幅度直方图计算。
[0047]
a)计算海杂波幅度直方图各统计量的数目
[0048]
根据步骤1中得到的归一化海杂波序列x'，利用如下公式计算海杂波幅度直方图各统计量的数目。
[0049][0050]
其中lx表示所有统计量数目的和，fix(
·
)表示向下取整，m、n表示任意正整数，但m和n的选择要尽可能的适应海杂波序列的长度与幅度，使绘制出的统计直方图更加美观。
[0051]
b)计算各统计量的幅度
[0052]
根据a中计算的海杂波幅度直方图中各统计量的数目，根据海杂波序列的幅度，将海杂波序列放置在每一个所属的统计单元内，统计单元内海杂波序列的数目nn，即每根柱子的幅度被存储在向量nn中，计算步骤如下：
[0053]
d＝(max(x')-min(x'))/lx
[0054]
(id,(i 1)d),i＝0,1
…
,lx-1
[0055]
其中d表示任意一个统计量的宽度，即任意一个统计量中所包含海杂波序列幅度的最大值减最小值，(id,(i 1)d),i＝0,1
…
,lx-1表示第i个统计量所能包含海杂波序列的幅度区间。针对上述计算结果绘制海杂波幅度统计直方图。
[0056]
步骤3、根据步骤2中得到的统计直方图，计算海杂波模型幅度统计熵。海杂波模型幅度统计熵的计算过程如下：
[0057]
f＝nn/n
[0058][0059]
其中f＝nn/n表示将频数nn转换成频率f，e表示海杂波模型幅度统计熵的计算结果，f(i)表示第i个频率分量。
[0060]
步骤4、目标探测
[0061]
根据步骤3中计算得到的幅度统计熵e的结果作为特征检测器的输入进行目标检测，特征检测器的结构如下所示：
[0062][0063]
其中z表示判决门限，由于不同距离单元幅度统计熵e的值是独立同分布的,所以判决门限z的值可以由非参数cfar方法决定，h0表示在该距离单元不存在目标，h1表示在该距离单元存在目标。这样设置判决门限的原因是：当出现目标时，回波信号的幅度会相对增加，导致统计直方图中柱子的分布扩展，从而使海杂波幅度统计熵e的值增加。
[0064]
本发明的效果可以通过以下实测数据试验进一步说明：
[0065]
实验条件：本实验采用x波段雷达实测进行海杂波幅度统计直方图分析和海杂波幅度熵的计算，x波段雷达的高度为30m，本实验采用雷达极化方式在hh下采集到的数据，scr约为0～6db。该x波段雷达的主要参数指标如表1所示，图2显示了x波段雷达回波的距离-时间-强度的图像。
[0066]
表1测试x波段雷达参数指标表
[0067][0068]
试验内容与结果：
[0069]
本发明结合信息论理论中利用熵来描述数据分散程度的特征，重点关注目标距离单元和海杂波距离单元海杂波幅度统计直方图之间的差异，并分别计算二者的幅度熵，利用非参数化cfar进行目标检测。
[0070]
图3分别是目标距离单元和海杂波距离单元的幅度统计直方图，其中图3.a是海杂波距离单元的统计直方图，图3.b是目标距离单元的统计直方图。结合两图，从对比中可以发现，图3.b中各统计量的分布情况比图3.a更加分散，这也印证了当目标出现时幅度熵会增加的结论。
[0071]
图4显示了不同距离单元幅度熵的计算结果，其中虚线代表非参数化cfar下计算出的门限值，圆圈代表不同距离单元的幅度熵。从图中可以发现，只有在第七距离单元的幅度熵的值大于门限值，其他周围的距离单元的幅度熵均小于门限值，这说明目标只出现在第七距离单元。
[0072]
因此可以得到结论：利用目标距离单元和海杂波距离单元在幅度熵数值上的差异可以在海上目标检测中区分海杂波与目标。
[0073]
目标检测性能分析：在这一部分，分析了本发明所提出方法的检测性能。表2展示了本发提出的目标检测方法与传统的cfar目标检测算法在虚警率为千分之一时的检测概率对比。
[0074]
表2不同方法的检测性能对比
[0075][0076]
从表2可以看出，本发明的检测性能优于传统的cfar目标检测算法，提高了在低信杂比情况下对海上小目标的探测性能。
[0077]
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明公开的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。