一种多波束测深数据化简方法及装置与流程

本发明涉及一种多波束测深数据化简方法及装置，属于海底地形综合技术领域。

背景技术：

点要素化简是制图综合中的重要研究内容和经典研究问题，当前点要素化简方法有很多，主要包括基于曲率参数设计的点化简方法、基于距离及高程差指标的点化简方法、基于不规则三角网的点化简方法和基于douglas-peucker算法的点化简方法等。如发表于《西安交通大学学报》2004年第7期上一篇名为《逆向工程中的测量数据精简技术研究》的论文通过分析角度—弦高联合准则法和包围盒法化简点云方法的使用性特点，提出了同时基于角度—弦高法和包围盒法的改进型数据的直接精简方法，实现了处理曲率变化大、附加特征多的点云数据化简。通用的点化简方法将关键特征保留作为方法设计的基础，由于多波束测深数据的一大主要用途是作为海图制图的源数据，而航海图是使用最为广泛的一类海图，因此针对多波束测深数据的化简，需以满足航海图制图应用为重要约束条件，即化简过程应保留局部区域深度值最小点以确保能够指导舰船安全航行，然而常用的点化简方法没有考虑此类应用需求。

面对上述需求，有人提出了新的化简方法，例如，名为《顾忌地形复杂度的多波束测深数据抽稀算法》的论文公开了一种多波束测深数据化简方法，该方法采用k-means算法将原始多波束测深点集进行聚类分组，然后以不同的分组为基本处理单元，实现特征点的选取与点集的化简工作。该方法虽然在顾及航行安全“取浅”原则的同时能够保留海底真实地形特征，但是该抽稀方法是基于对原始多波束测深点进行k-means聚类的基础上进行的，初始聚类中心的选择和非地形特征点中的离群点均易对聚类结果产生较大影响，，若聚类结果不恰当，将会降低点化简结果的适用性。

技术实现要素：

本发明的目的是提供一种多波束测深数据化简方法及装置，以解决目前点化简方法中存在的因聚类结果不恰当导致的化简结果适用性差的问题。

本发明为解决上述技术问题而提供了一种多波束测深数据化简方法，该化简方法包括以下步骤：

1)计算待化简区域的地形复杂度，根据计算得到的地形复杂度将待化简区域划分成若干具有不同复杂程度的子区域；

2)对整个待化简区域进行规则格网划分，规则格网尺寸由待化简区域中各不同地形复杂程度类别区域的面积和化简后保留点的个数确定；

3)在各子区域的格网内采用相应的水深点化简规则，使得化简后各子区域水深点密度与该子区域的地形复杂度正相关。

本发明还提供了一种多波束测深数据化简装置，该化简装置包括存储器和处理器，以及存储在所述存储器上并在所述处理器上运行的计算机程序，所述处理器与所述存储器相耦合，所述处理器执行所述计算机程序时实现本发明的多波束测深数据的化简方法。

本发明采用稳定性更高的格网分割法实现测深数据点集的分组，避免了出现因聚类结果不正确而导致的点化简效果较差的情况，提高了化简后数据的适用性；同时，根据各局部区域地形复杂程度的不同，采用在复杂度较高区域内保留点密度较大的点化简策略，保证了化简后数据表达地形的准确性，使得化简结果不仅满足航海图制图应用，还可用于构建海底地形模型和进行地形分析。

进一步地，为保证化简过程顾及航行安全性准则，在对各子区域进行化简时需保留该区域各个规则格网内水深最浅点。

进一步地，所述步骤1)按复杂程度分成了平坦区域、一般区域和复杂区域，从舰船航行安全方面考虑，可将地形平坦区域划分为复杂区域，而不能反向划分，因此在对各类别地形的复杂度范围进行边界确定时采用缩平坦扩复杂的原则。

进一步地，本发明给出了各区域化简的具体化简策略，所述平坦区域各格网仅保留该格网内水深值最小点；一般区域各格网保留该格网内水深值最小点和水深值最大点；复杂区域各格网保留该格网内水深值最小点、该格网内水深值最大点和该格网内与平均水深值之差最小点。

进一步地，本发明给出了具体的复杂度确定方式，所述步骤1)采用多种地形因子的复杂度评价指标计算方法确定待化简区域的地形复杂度，地形因子包括坡度、局部地形粗糙度、局部地表起伏度、局部标准差和局部地表切割深度。

进一步地，本发明给出了复杂度的计算公式，所述多种地形因子的复杂度计算方法采用的计算公式为：

c为复杂度评价指标，coei为第i个地形因子，n为地形因子的个数，ti为第i个地形因子的贡献系数，通过因子分析法得到。

进一步地，本发明给出了化简后保留点数与化简规则的具体约束关系，化简后保留点数和各子区域包含的规则格网数量之间的关系为：

na＝np+2ny+3nf

其中na为化简后保留点个数，np为地形平坦子区域包含的规则格网数量，ny为地形一般子区域包含的规则格网数量，nf为地形复杂子区域包含的规则格网数量。

进一步地，本发明给出了化简后保留点数和规则格网尺寸的具体关系，化简后保留点数和各子区域面积及划分的规则格网尺寸之间的关系为：

其中na为化简后保留点个数，sp为地形平坦子区域的面积，sy为地形一般子区域的面积，sf为地形复杂子区域的面，d为规则格网边长。

附图说明

图1为本发明实施例中各类别区域化简规则示意图；

图2为本发明实施例中某海区原始海底地形三维模拟示意图；

图3为本发明实施例中某海区地形复杂程度分类图；

图4为本发明实施例中化简率为70％时某海区水深点化简结果构建的地形模型图；

图5为本发明实施例中化简率为60％时某海区水深点化简结果构建的地形模型图；

图6为本发明实施例中化简率为50％时某海区水深点化简结果构建的地形模型图；

图7为本发明实施例中化简率为40％时某海区水深点化简结果构建的地形模型图。

具体实施方式

下面结合附图对本发明的具体实施步骤作进一步地说明。

方法实施例

本发明首先利用地形分类方法将待化简海区划分成若干具有不同复杂程度类别的子区域，然后对待化简海区进行规则格网划分，规则格网的尺寸可由各类别地形的面积与整个实验海区化简后保留点的个数确定，最后在不同地形类别区域的格网内采用不同的水深点化简规则，使得化简后各地形类别区域水深点密度与该类别的地形复杂程度呈正相关。

下面以某海区为例对本发明具体实现过程进行详细说明，其中该海区的原始多波束测深数据对应的海底地形三维模拟示意图如图2所示。

1.地形复杂度计算。

本发明利用顾及多种单一地形因子的复杂度评价方法对真实海底地形的复杂度进行量化描述。具体而言，选用的单一地形因子包括坡度、局部地形粗糙度、局部地表起伏度、局部标准差和局部地表切割深度，每种地形因子的贡献系数可采用因子分析法进行求解。本实施例采用的地形复杂度的具体计算公式为：

c为地形复杂度，coei为第i个地形因子，n为地形因子的个数，ti为第i个地形因子的贡献系数。

其中坡度因子的计算可参考发表于《四川农业大学学报》2013年第4期上一篇名为《基于arcgis计算srtm数据dem坡度的最优方法探讨》论文中提到的三阶不带权差分法；局部地形粗糙度、局部地表起伏度和局部标准差三个因子的计算方法可参考发表于《山地学报》2012年第5期上的一篇名为《地形复杂度的多因子综合评价方法》论文；局部地表切割深度因子的计算方法可参考发表于《重庆大学学报》2012年第11期上的一篇名为《基于gis的地表破碎指数计算》论文。

2.根据复杂度进行分类。

通过选取若干地形复杂程度为平坦、一般和复杂三类典型地形区域，采用步骤1中地形复杂度量化表达方法分别计算各类地形的复杂度值并统计复杂度范围，根据统计结果确定各类别地形对应的地形复杂度量化范围，定义统计得到的三类典型地形区域的复杂度范围分别为平坦：(0,c1)、一般：(c2,c3)、复杂：(c4,1)。

由于各类地形之间划分的模糊性和主观性，得到的各类地形区域的复杂度范围并非严格连续且不重叠，如：c1<c2或c1>c2，基于此分类范围对地形进行分类时有可能无法确定区域所属地形类别或导致地形分类出现错误，因此需明确各类别地形的复杂度范围为严格连续且无重叠。从舰船航行安全方面考虑，可将地形平坦区域划分为复杂区域，而不能反向划分，因此在对各类别地形的复杂度范围进行边界确定时采用缩平坦扩复杂的原则。因此本发明定义各类别地形对应的地形复杂度量化范围为：平坦(0，min(c1,c2))、一般(min(c1,c2)，min(c3，c4))、复杂(min(c3,c4)，1)。

在上述思想基础上，对本实施例中的某海区进行地形复杂程度分类时，只需采用步骤1中的方法计算地形复杂度，根据以上确定的各类别地形的复杂度量化范围对其进行地形类别的划分，分成平坦区域、一般区域和复杂区域，其结果如图3所示。

3.确定格网尺寸，对各区域进行规则格网划分。

按照步骤2的地形划分结果，确定各类区域地形面积，按照设定的化简率和各类区域的化简策略确定格网尺寸，基于确定的规则格网尺寸对各类区域进行规则格网划分。化简后保留点数量na与设定的化简率r有关，

na＝n×r

其中n为原始测深点数量。

4.各类别区域化简规则制定

综合考虑保证舰船安全航行的“取浅原则”和构建真实准确的海底地形模型两个重要约束因素的影响，在不同类别区域制定不同的化简策略，使得化简结果满足上述两个领域的应用需求。

本实施例采用的化简规则如图1所示，针对不同的化简规则分别对每个地形类别中各规则格网内的测深点进行化简，其中平坦区域各格网仅保留该格网内水深值最小点，一般区域各格网保留该格网内水深值最小点和最大点，复杂区域各格网保留该格网内水深值最小点、该格网内水深值最大点和该格网内与平均水深值之差最小点。通过上述化简规则，可实现对各类区域的测深数据进行化简。

化简后保留点数和各子区域包含的规则格网数量之间的关系为：

na＝np+2ny+3nf

其中na为化简后保留点个数，np为地形平坦子区域包含的规则格网数量，ny为地形一般子区域包含的规则格网数量，nf为地形复杂子区域包含的规则格网数量。

化简后保留点数和各子区域面积及划分的规则格网尺寸之间的关系为：

其中na为化简后保留点个数，sp为地形平坦子区域的面积，sy为地形一般子区域的面积，sf为地形复杂子区域的面，d为规则格网边长。

在不同化简率时，利用本发明所提出的化简方法对某海区的原始水深点按照不同的化简率进行化简，化简结果构建的海底地形模型分别为图4—图7所示。

装置实施例

本发明的多波束测深数据的化简装置，包括存储器和处理器，以及存储在所述存储器上并在处理器上运行的计算机程序，处理器与所述存储器相耦合，处理器执行计算机程序时实现方法实施例中的多波束测深数据的化简方法，具体过程已在方法实施例中进行了详述，这里不再赘述。

根据本发明的各类别区域点化简方法可知，本发明能够确保点化简结果保留局部区域水深最浅点的同时能够保持区域地形形态的变化特征，保证了化简后数据表达地形的准确性，使得化简结果不仅满足航海图制图应用，还可用于构建海底地形模型和进行地形分析。同时本发明采用稳定性更高的格网分割法实现点集的分组，不会出现因聚类结果不正确而导致的点化简效果较差的情况，在实际应用过程中更加易用。