基于双谱三特征的海面目标检测方法及系统与流程

1.本发明涉及信号处理技术领域，尤其涉及一种基于双谱三特征的海面目标检测方法及系统。


背景技术：

2.在雷达探测领域，对海面慢速小目标的检测一般通过能量检测方法和特征检测方法两种方式实现。其中，能量检测方法主要依据海杂波的局部幅度或功率水平信息构造似然比，根据门限因子形成检测门限并对目标存在与否做出判决。该方法所需的雷达驻留时间较短，但对目标的信杂比要求较高，且在海尖峰密集场景下易引起大量虚警。
3.特征检测方法主要通过挖掘海杂波与目标之间的差异特征，将二者从高重叠的观测空间转换到低重叠的特征空间，在特征空间内实现目标检测。特征检测方法所需的雷达驻留时间较长，需要对海面进行长时间的观测积累，通常可达几百毫秒甚至几秒。然而，在雷达的扫描观测模式下，目标驻留时间通常难以达到上述量级，随着积累时间的降低，特征检测方法性能下滑严重。
4.不难发现，上述两种方式得到的目标检测结果均不够准确，难以满足实际应用需求。
5.因此，现在亟需一种海面目标检测方法来解决上述问题。


技术实现要素：

6.本发明提供一种基于双谱三特征的海面目标检测方法及系统，用以解决现有技术中海面目标检测结果不准确的技术问题。
7.第一方面，本发明提供一种基于双谱三特征的海面目标检测方法，包括：
8.获取待检测信号，构造所述待检测信号的对角积分双谱；
9.从所述对角积分双谱中提取谱特征并进行累积得到双谱三特征，构建三特征空间，得到所述待检测信号的特征点；其中，所述双谱三特征包括累积峰值、累积谱宽以及全变差；
10.对所述待检测信号的特征点进行检测分类，得到目标检测结果。
11.根据本发明提供的一种基于双谱三特征的海面目标检测方法，所述构造所述待检测信号的对角积分双谱，包括：
12.计算所述待检测信号的二维双谱；
13.从所述二维双谱中提取得到对角积分双谱。
14.根据本发明提供的一种基于双谱三特征的海面目标检测方法，所述二维双谱由所述待检测信号的三阶累积量做二维傅里叶变换得到。
15.根据本发明提供的一种基于双谱三特征的海面目标检测方法，所述从所述对角积分双谱中提取谱特征并进行累积得到双谱三特征，构建三特征空间，得到所述待检测信号的特征点，包括：
16.从所述对角积分双谱中提取峰值、质心频率和谱宽；
17.分别将所述峰值和所述谱宽沿帧时间维度进行累积，得到累积峰值和累积谱宽；
18.提取用于衡量帧时间维度上所述质心频率特征稳定性的全变差；
19.以累积峰值特征为x轴、累积谱宽特征为y轴、全变差特征为z轴，构建得到三特征空间，得到所述待检测信号的特征点。
20.根据本发明提供的一种基于双谱三特征的海面目标检测方法，对所述待检测信号的特征点进行检测分类，得到目标检测结果，包括：
21.采集海杂波数据，通过凸包算法对所述海杂波数据进行特征训练，得到决策区域；
22.判断所述待检测信号的特征点是否落入所述决策区域，若所述待检测信号的特征点落入所述决策区域，则所述待检测信号为杂波单元；若所述待检测信号的特征点未落入所述决策区域，则所述待检测信号为目标单元。
23.根据本发明提供的一种基于双谱三特征的海面目标检测方法，对所述待检测信号的特征点进行检测分类，得到目标检测结果，还包括：
24.当所述海杂波数据的特征因海况或雷达视角的改变而发生变化时，重新采集海杂波数据以确定新的决策区域。
25.第二方面，本发明还提供一种基于双谱三特征的海面目标检测系统，包括：
26.对角积分双谱构造模块，用于获取待检测信号，构造所述待检测信号的对角积分双谱；
27.三特征提取模块，用于从所述对角积分双谱中提取谱特征并进行累积得到双谱三特征，构建三特征空间，得到所述待检测信号的特征点；其中，所述双谱三特征包括累积峰值、累积谱宽以及全变差；
28.目标检测模块，用于对所述待检测信号的特征点进行检测分类，得到目标检测结果。
29.根据本发明提供的一种基于双谱三特征的海面目标检测系统，所述三特征提取模块包括：
30.谱特征提取单元，用于从所述对角积分双谱中提取峰值、质心频率和谱宽；
31.第一特征提取单元，用于分别将所述峰值和所述谱宽沿帧时间维度进行累积，得到累积峰值和累积谱宽；
32.第二特征提取单元，用于提取用于衡量帧时间维度上质心频率特征稳定性的全变差；
33.特征点获取单元，用于以累积峰值特征为x轴、累积谱宽特征为y轴、全变差特征为z轴，构建得到三特征空间，得到待检测信号的特征点。
34.第三方面，本发明还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述任一种所述基于双谱三特征的海面目标检测方法的步骤。
35.第四方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述任一种所述基于双谱三特征的海面目标检测方法的步骤。
36.本发明提供的基于双谱三特征的海面目标检测方法及系统，通过构造待检测信号
的对角积分双谱，便于提取累积峰值、累积谱宽以及全变差这三个特征，利用融合上述三个特征的特征点进行目标检测，可以更好的区分目标单元与杂波单元，得到的目标检测结果更加准确可靠。
附图说明
37.为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
38.图1是本发明提供的基于双谱三特征的海面目标检测方法的流程示意图；
39.图2(a)是应用本发明提供的基于双谱三特征的海面目标检测方法对实测数据进行目标检测得到的检测结果示意图；
40.图2(b)是应用基于时频三特征的检测器对实测数据进行目标检测得到的检测结果示意图；
41.图2(c)是应用幅度、多普勒三特征检测器对实测数据进行目标检测得到的检测结果示意图；
42.图2(d)是应用分形检测器对实测数据进行目标检测得到的检测结果示意图；
43.图3为本发明提供的基于双谱三特征的海面目标检测系统的结构架构示意图；
44.图4为基于双谱三特征的海面目标检测系统中三特征提取模块的结构架构示意图；
45.图5为本发明实施例提供的电子设备的结构架构示意图。
具体实施方式
46.为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
47.图1示出了本发明实施例提供的基于双谱三特征的海面目标检测方法，包括：
48.s110：获取待检测信号，构造待检测信号的对角积分双谱。
49.在进行目标检测前，本发明实施例假设对海雷达在一个波束内发射一系列相干脉冲，并在每个距离单元接收到长度为k的回波时间序列x(k)。若该回波信息未受到目标影响，则待检测单元中应只包含海杂波及噪声信号c(k)，否则还应混有目标信号s(k)。基于此，海杂波中的目标检测可以归结为的二元假设检验问题，进而可以表示为：
[0050][0051]
可以理解的是，本步骤中提到的待检测信号，本质上是一串有限长的时间序列。
[0052]
首先计算待检测信号{x(t),t＝1,2,
…
,t}的二维双谱b(ω1,ω2)，它被定义为信号x(t)的三阶累积量c3(τ1,τ2)的二维傅里叶变换，其表达式如下：
[0053][0054]
式中，c3(τ1,τ2)需满足绝对可和条件，且存在：
[0055]
c3(τ1,τ2)＝e[x
h
(t)x(t τ1)x(t τ2)]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0056]
进一步，提取其中的对角积分双谱dib，其表达式如下：
[0057][0058]
式中，坐标系(ω1′
,ω2′
)由原双谱坐标系(ω1,ω2)逆时针旋转45度得到。
[0059]
s120：从对角积分双谱中提取谱特征并进行累积得到双谱三特征，构建三特征空间，得到待检测信号的特征点；其中，双谱三特征包括累积峰值、累积谱宽以及全变差。
[0060]
首先，根据海杂波单元与目标单元之间的非线性耦合差异性，从对角积分双谱中提取峰值pv、质心频率fc、谱宽sw三种特征，三个谱特征的表达式如下：
[0061]
pv＝max{dib(ω)}
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0062]
fc＝∫ω
·
dib(ω)dω/∫dib(ω)dω
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0063][0064]
考虑到雷达观测时间较短时，海杂波与目标之间的可分性不足。为此，本发明通过多帧扫描历史数据和当前帧数据的综合应用，将峰值、谱宽这两个特征分别沿帧时间维度进行累积，得到累积峰值cpv、累积谱宽csw两个特征，同时提取全变差tv这一特征，用于衡量帧时间维度上质心频率特征的稳定性，累积峰值cpv、累积谱宽csw以及全变差tv三个累积特征的表达式如下：
[0065][0066][0067][0068]
其中，s表示谱特征的累积起点，l表示谱特征的累积帧数。
[0069]
据此，可以构造一个三特征空间，其中，x轴为累积峰值特征cpv，y轴为累积谱宽特征csw，z轴为全变差特征tv。从而得到待检测信号的特征点，该特征点包含目标单元与海杂波单元之间的三个差异特征，通过该特征点，可以更加准确的区分目标单元与海杂波单元。同时，本发明提供的方法在三特征空间进行融合检测，能够避免单特征检测的局限性。
[0070]
s130：对待检测信号的特征点进行检测分类，得到目标检测结果。
[0071]
在本实施例中，前期需要采集海杂波数据，通过凸包算法对海杂波数据进行特征训练，得到决策区域。
[0072]
采用凸包算法构建决策区域的过程，具体如下：
[0073]
(1)初始化：设置海杂波数据的特征训练集为h，其中包含w个数据，计算虚警数为l＝wp
f
，其中p
f
为设定的虚警率。令起始变量l＝0。
[0074]
(2)计算当前数据点的凸包ch(h)，其中，凸包ch(h)的顶点为{v1,v2,
…
,v
r
}。然后
统计落入凸包ch(h)中的数据点数量，设为n
all
。
[0075]
(3)设置一个循环变量q，令其从1到r，计算凸包ch(h
‑
{v
q
})，即从集合h中去除顶点v
q
，然后计算处于新凸包ch(h
‑
{v
q
})中的特征点数量，设为n
q
。
[0076]
(4)比较集合{n
all
‑
n1,
…
,n
all
‑
n
r
}，去掉该集合中最大数值所对应的顶点v
i
。
[0077]
(5)令h
‑
{v
i
}＝h,l 1＝l。
[0078]
(6)若l＜l，返回步骤(2)继续下一个顶点的移除。否则，终止移除过程，输出最终的决策区域ω＝ch(h)。
[0079]
由于该决策区域是在满足虚警率要求的前提下，由凸包算法逐步缩小凸包划定范围后得到的，可以更好的区分当前特征点是杂波单元还是目标单元，从而更加准确的实现目标检测。
[0080]
之后，判断待检测信号的特征点是否落入上述决策区域ω内，若待检测信号的特征点落入决策区域ω，则待检测信号为杂波单元；若待检测信号的特征点未落入决策区域ω，则待检测信号为目标单元。
[0081]
在本步骤中，采集的海杂波数据需要通过三特征空间获得相应的特征点，由海杂波数据对应的特征点构建得到特征训练集。
[0082]
可以理解的是，当海杂波特征因海况或雷达视角的改变而发生明显变化时，需要重新采集海杂波数据以确定新的决策区域。
[0083]
考虑到高阶谱作为信号分析的有效工具，可以直观展现各谐波分量之间的频率、相位耦合现象，近年来在图像和物体识别等领域获得了广泛应用。其中，双谱作为阶数最低的高阶谱，在具有时移不变性、相位保持性、对称性等优良特性的同时依然能够保持高阶谱的所有特征，因此高阶谱分析时通常选择双谱作为研究对象。本发明从待检测信号的双谱特性入手，分析目标单元与海杂波单元之间的差异特征，并从中提取三个特征用于海面小目标检测，避免了单特征检测的局限性，从而有效提高了目标检测的准确率。
[0084]
由于现有技术中存在已知的基于时频三特征的检测器，该检测器融合应用归一化平滑伪维格纳
‑
威利分布中的时频脊累积量、连通区域个数和最大连通区域尺寸信息获得。同时，还存在幅度、多普勒三特征检测器，该检测器从雷达回波中提取了相对平均幅度、相对多普勒峰高以及相对多普勒熵三个特征，并在三维空间中利用凸包算法实现检测分类。
[0085]
为展示本发明提供的目标检测方法的优势，本实施例采用实测数据进行分析与验证，具体通过本发明提供的方法、上述基于时频三特征的检测器和幅度、多普勒三特征检测器以及公知的分形检测器分别对同一实测数据进行目标检测，虚警率均设置为10
‑3。其中，本发明实施例得到的检测结果如图2(a)所示，图2(b)、图2(c)、图2(d)分别示出了已有方法中基于时频三特征检测器、幅度、多普勒三特征检测器以及分形检测器的检测结果，从各自的检测结果图中可以直观的看出，相比已有的特征检测器，本发明实施例给出的检测方法在实际虚警率控制和检测准确率提升方面均具有明显优势。
[0086]
由此可见，本发明实施例提供的上述基于双谱三特征的海面目标检测方法，通过构造待检测信号的对角积分双谱，便于获取累积峰值、累积谱宽以及全变差这三个差异特征，进而利用融合上述三个特征的特征点进行目标检测，可以更好的区分目标单元与杂波单元，得到的目标检测结果更加准确可靠。
[0087]
图3示出了本发明实施例提供的基于双谱三特征的海面目标检测系统，包括：
[0088]
对角积分双谱构造模块310，用于获取待检测信号，构造待检测信号的对角积分双谱；
[0089]
三特征提取模块320，用于从对角积分双谱中提取谱特征并进行累积得到双谱三特征，构建三特征空间，得到待检测信号的特征点；其中，双谱三特征包括累积峰值、累积谱宽以及全变差；
[0090]
目标检测模块330，用于对待检测信号的特征点进行检测分类，得到目标检测结果。
[0091]
本实施例中对角积分双谱构造模块310包括获取单元、二维双谱计算单元以及对角积分双谱提取单元，具体地，对角积分双谱构造模块通过获取单元获取待检测信号，获取到待检测信号后，进一步通过二维双谱计算单元计算待检测信号的二维双谱；最后通过对角积分双谱提取单元从二维双谱中提取得到对角积分双谱。
[0092]
具体地，参见附图4，三特征提取模块320包括：
[0093]
谱特征提取单元321，用于从对角积分双谱中提取峰值、质心频率和谱宽；
[0094]
第一特征提取单元322，用于分别将峰值和谱宽沿帧时间维度进行累积，得到累积峰值和累积谱宽；
[0095]
第二特征提取单元323，用于提取用于衡量帧时间维度上质心频率特征稳定性的全变差；
[0096]
特征点获取单元324，用于以累积峰值特征为x轴、累积谱宽特征为y轴、全变差特征为z轴，构建得到三特征空间，得到待检测信号的特征点。
[0097]
本实施例中目标检测模块330包括决策区域生成单元和判断单元，具体地，目标检测模块首先通过决策区域生成单元采集海杂波数据，利用凸包算法对海杂波数据进行特征训练，得到决策区域；然后通过判断单元判断待检测信号的特征点是否落入决策区域，若待检测信号的特征点落入决策区域，则待检测信号为杂波单元；若待检测信号的特征点未落入决策区域，则待检测信号为目标单元。
[0098]
更优地，为了保证系统数据处理的实时性和可靠性，目标检测模块330还可以包括决策区域更新单元，当海杂波数据的特征因海况或雷达视角的改变而发生变化时，通过决策区域更新单元可以重新采集海杂波数据以确定新的决策区域。
[0099]
由此可见，本发明提供的基于双谱三特征的海面目标检测系统，通过构造待检测信号的对角积分双谱，便于提取累积峰值、累积谱宽以及全变差这三个特征，通过目标检测模块利用融合上述三个特征的特征点进行目标检测，可以更准确的区分目标单元与杂波单元，得到的目标检测结果更加准确可靠。
[0100]
图5示例了一种电子设备的实体结构示意图，如图5所示，该电子设备可以包括：处理器(processor)510、通信接口(communications interface)520、存储器(memory)530和通信总线540，其中，处理器510，通信接口520，存储器530通过通信总线540完成相互间的通信。处理器510可以调用存储器530中的逻辑指令，以执行基于双谱三特征的海面目标检测方法，该方法包括：获取待检测信号，构造待检测信号的对角积分双谱；从对角积分双谱中提取谱特征并进行累积得到双谱三特征，构建三特征空间，得到待检测信号的特征点；其中，双谱三特征包括累积峰值、累积谱宽以及全变差；对待检测信号的特征点进行检测分类，得到目标检测结果。
[0101]
此外，上述的存储器530中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read
‑
only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0102]
另一方面，本发明还提供一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法所提供的基于双谱三特征的海面目标检测方法，该方法包括：获取待检测信号，构造待检测信号的对角积分双谱；从对角积分双谱中提取谱特征并进行累积得到双谱三特征，构建三特征空间，得到待检测信号的特征点；其中，双谱三特征包括累积峰值、累积谱宽以及全变差；对待检测信号的特征点进行检测分类，得到目标检测结果。
[0103]
又一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各提供的基于双谱三特征的海面目标检测方法，该方法包括：获取待检测信号，构造待检测信号的对角积分双谱；从对角积分双谱中提取谱特征并进行累积得到双谱三特征，构建三特征空间，得到待检测信号的特征点；其中，双谱三特征包括累积峰值、累积谱宽以及全变差；对待检测信号的特征点进行检测分类，得到目标检测结果。
[0104]
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
[0105]
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如rom/ram、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
[0106]
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。