一种基于浅水多波束水深数据的水下目标实时分割方法及应用与流程

1.本发明涉及水底地形观测和水下目标探测技术领域，具体涉及一种基于浅水多波束水深数据的水下目标实时分割方法及应用。


背景技术：

2.随着长江流域经济的快速发展，长江已成为世界上通货量最大、运输最为繁忙的通航河流，为了维护长江航道通航能力，如何快速高效地进行水下碍航物探测和清理是实现长江航道运力最大化的保障之一。以kongsberg em2040、reson t50p、sonic2024为代表的第五代浅水多波束测深系统普遍具有宽覆盖、高精度、高分辨率等特点，不仅能获得精细的水下地形，还能实现航行障碍物探测、水下建筑物监测、水下物体搜寻、水下清淤扫测等目的，多波束水下测量技术已成为长江航道“数字航道”建设的核心测量技术。
3.目前，基于多波束的水下目标自动化探测大多利用后向散射图像和水体图像数据在后处理时进行，且识别率低，还没有多波束硬件和软件利用水深点数据进行实时目标探测的应用。基于三维激光点云的目标识别和重建技术在地面测量中已获得广泛应用，而浅水多波束系统同样能获得水下目标的高密度点云且分割目标通常较为单一，分割场景简单，但是目前仅用于定性分析，主要通过测量员利用显示软件目视找出目标位置。
4.点云分割是根据点云的空间、几何和纹理等特征对点云进行划分，使得同一划分点云具有相似特征，目前点云的分割方法主要有基于边的方法、基于面的方法、基于聚类的方法和混合分割方法，其中基于边的分割方法优点是速度快，对明显边界识别能力强，但是容易受测量噪声影响，因此不适用于多波束的点云分割；基于面的分割方法主要有区域生长法和ransac随机抽样一致性算法，这两种算法都比较依赖给定的平面分割阈值，区域生长法需要设定法线角度差和曲率两个阈值，而ransac算法需要设定离群点到拟合平面的容忍距离；区域生长法能一次性分割出多个平面特征，对规则平面如四四方方的建筑物等分割较好，但是难以处理噪点较多且凹凸不平的地面，因此不适用于多波束的点云分割。基于面的分割方法的缺陷是无法分割出非平面特征，对于多波束数据，平面分割能区分水下的平坦地形背景和混合目标，但是无法将多个目标分割开来。基于聚类的方法利用点云的局部几何特征参数，例如高斯曲率、平均曲率、法矢量、坐标、距离等进行聚类，为了加速邻域点的搜索，通常采用kd
‑
tree进行索引，但是聚类算法对噪声敏感。
5.通过上述分析，现有技术存在的问题及缺陷为：
6.1、目前，使用多波束扫海技术自动寻找水下目标的方式主要是利用二维图像数据，还没有利用三维水深点数据进行水下目标分割和探测的应用。
7.2、多波束水深数据大多存在噪点，给水下目标的分割造成了极大阻碍。


技术实现要素：

8.本发明的技术任务是针对现有技术的不足，提供一种基于浅水多波束水深数据的
水下目标实时分割方法及应用，能够解决利用多波束进行水下目标扫测依赖人工目视检测的问题。
9.本发明解决其技术问题所采用的技术方案是：
10.1、本发明提供一种基于浅水多波束水深数据的水下目标实时分割方法及应用，包括如下步骤：
11.s1、实时接收30～50ping多波束原始数据至缓冲区并进行波束点归位计算；
12.s2、采用半径滤波器，以平均水深的3％～5％作为滤波半径对缓冲区内的水深点进行大尺度噪点滤波；
13.s3、采用平均水深的8％～10％作为平面距离阈值对滤波后的水深点进行ransac平面分割，对分割后的局外点和局内点进行平均深度比较，判断局外点所在目标特征的凹凸性，若为凹陷则交换局外点和局内点；
14.s4、将分割出来的局外点进行欧式距离聚类得到单个聚类目标，将单个目标聚类点数小于5的聚类目标剔除；
15.s5、计算各个聚类目标的体积、中心位置等属性并三维显示，清空缓冲区并重复步骤s1～s4。
16.进一步地，所述步骤s2中，采用半径滤波器进行滤波时，若该半径范围内的其余波束点个数小于2则该点视为离群点；若采集的多波束数据质量较差，可适当放宽离群点数值限制。
17.进一步地，所述步骤s3中，平面距离分割阈值设置为缓冲区内所有经过滤波后水深点的平均水深的8％～10％，若局外点占比超过50％，则采用最小中值方差采样一致性方法进行平面分割，平面分割完成后将背景地形的水深点剔除。
18.进一步地，所述步骤s3中，采用ransac平面分割算法解算平面参数的具体公式为：
[0019][0020]
式中，a、b、c为平面方程ax+by+cz+1＝0的参数，x、y、z为水深点的坐标。
[0021]
进一步地，所述步骤s4中，聚类距离阈值设置为平均水深的8％～10％；若聚类目标含有较多噪点，需要设置最小聚类点数，一般设置为5即可，将低于此聚类点数的聚类目标删除。
[0022]
进一步地，所述步骤s4中，欧式距离聚类是一种最近邻聚类算法，通过计算相邻点之间的欧式距离：
[0023][0024]
将离散的点集按照一定的距离尺度划分为一系列的点簇，给定聚类距离d(水深的8％～10％)，如果两个相邻点之间的欧式距离小于或等于d，那么这两个点将被聚为同一类；否则为不同类，具体实现流程为：
[0025]
建立一个kd树存储点云数据p，设置一个空的点簇链表c和空队列q；对每个点p
i
∈p，添加p
i
到当前队列q，对每个点p
i
∈q：以d为半径，搜索点p
i
的邻域点(符合欧式距离小于
或等于d的原则)存入p
ik
，对每个检查该点是否被处理过，如果没有，则加入q，当链表中所有在q中的点都被处理过时，把q加入到链表c中，q置空；当所有p中的点都被处理过时，终止。
[0026]
2、本发明另提供一种基于浅水多波束水深数据的水下目标实时分割方法的应用，其特征在于，基于权利要求1所述的浅水多波束水深数据的水下目标实时分割方法，进行水下碍航物探测和清理。
[0027]
本发明的一种基于浅水多波束水深数据的水下目标实时分割方法及应用与现有技术相比，所产生的有益效果是，
[0028]
本发明针对高密度多波束水深点数据，通过归位计算、自动滤波等预处理，对处理后的多波束水深点数据采用基于ransac平面分割和欧式距离聚类的混合分割方法，实现了水下目标的实时快速分割以及物理属性和地理属性标定。该方法的实现无需用户干预，全自动化处理，目标分割效果良好。
[0029]
采用ransac平面分割去掉地形背景再结合欧式距离聚类能够有效分割如石块、沉船、木头等水下目标，在分割前和分割后都需要进行滤波处理，根据水深结合多波束系统参数合理地设置阈值是分割成功的关键。本发明实现了利用走航多波束快速检测水下障碍物的目的。第五代浅水多波束系统具有高分辨率、高精度等特性，能够探测厘米级的水下目标，目前关于声呐的探测基本上都是基于图像的检测，本发明直接基于水深点数据进行检测并三维重构，在显示上较二维图像更直观。
[0030]
说明书附图
[0031]
图1为较佳实施例的流程示意图；
[0032]
图2为本发明实施例提供的半径滤波示意图；
[0033]
图3为本发明实施例提供的沉船分割结果示意图；
[0034]
图4为本发明实施例提供的多组水下目标分割结果示意图；
[0035]
图5为本发明实施例提供的浅水多波束水下目标检测系统流程图；
[0036]
图6为本发明实施例提供的浅水多波束水下目标实时分割系统的示意图。
具体实施方式
[0037]
下面结合附图1
‑
6，对本发明的一种基于浅水多波束水深数据的水下目标实时分割方法及应用作以下详细地说明。
[0038]
如附1
‑
6所示，本发明的一种基于浅水多波束水深数据的水下目标实时分割方法及应用，包括如下步骤：
[0039]
s1、实时接收30～50ping多波束原始数据至缓冲区并进行波束点归位计算；
[0040]
s2、采用半径滤波器，以平均水深的3％～5％(em2040标称测深精度的10倍)作为滤波半径对缓冲区内的水深点进行大尺度噪点滤波；例如300khz短脉宽发射模式在20m水深精度为3cm，滤波半径设为0.1m。
[0041]
s3、采用平均水深的8％～10％作为平面距离阈值对滤波后的水深点进行ransac平面分割，对分割后的局外点和局内点进行平均深度比较，判断局外点所在目标特征的凹凸性，若为凹陷则交换局外点和局内点；具体过程如下：
[0042]
s31：在局部范围内，可以假设水底为一个平面：
[0043]
ax+by+cz+1＝0
ꢀꢀꢀ
(公式1)
[0044]
s32：分割的核心在于确定平面模型系数a、b、c，当有n个拟合点时，求解过程表示为矩阵的形式如下：
[0045][0046]
通过变换得到：
[0047][0048]
s33：通过反复的随机抽样子集确定平面参数模型，用得到的模型去测试所有的其它数据，如果某个点适用于估计的模型，则认为它也是局内点，如果有足够多的点被归类为假设的局内点，那么估计的模型就足够合理；最后通过所有的局内点重新计算模型并估计局内点与模型的错误率来评估模型。
[0049]
s4、将分割出来的局外点进行欧式距离聚类得到单个聚类目标，将单个目标聚类点数小于5的聚类目标剔除；
[0050]
s41：将离散的点集按照一定的距离尺度划分为一系列的点簇，给定聚类距离d(水深的8％～10％)，如果两个相邻点之间的欧式距离小于或等于d，那么这两个点将被聚为同一类；否则为不同类，具体实现流程为：
[0051]
建立一个kd树存储点云数据p，设置一个空的点簇链表c和空队列q；对每个点p
i
∈p，添加pi到当前队列q，对每个点p
i
∈q：以d为半径，搜索点pi的邻域点(符合欧式距离小于或等于d的原则)存入p
ik
，对每个检查该点是否被处理过，如果没有，则加入q，当链表中所有在q中的点都被处理过时，把q加入到链表c中，q置空；当所有p中的点都被处理过时，终止。
[0052]
s5、计算各个聚类目标的体积、中心位置等属性并三维显示，清空缓冲区并重复步骤s1～s4，直到完成所有数据的分割。
[0053]
表1为本发明实施例提供的不同样本目标的分割效果、阈值设置和耗时对比结果
[0054][0055][0056]
本说明书所公开的实施例只是对本发明单方面特征的一个例证，本发明的保护范围不限于此实施例，其他任何功能等效的实施例均落入本发明的保护范围内。对于本领域的技术人员来说，可根据以上描述的技术方案以及构思，做出其它各种相应的改变以及变形，而所有的这些改变以及变形都因该属于本发明权利要求的保护范围之内。