一种基于多波束测深数据的Delaunay三角网快速生成方法技术

本发明专利技术公开了一种基于多波束测深数据的Delaunay三角网快速生成方法，针对目前不规则三角网生成算法中定位待插入点所在三角形及规格化三角网效率不高的问题，本发明专利技术提出一种基于多波束测深数据快速定位的方法和一种通过对角的正切值计算来判断三角网是否需要优化的方法。根据多波束测深数据每一PING数据在一条直线上的特性，对每一PING数据进行关联插入，相比于传统的逐个检测的方法，大大缩短了寻找时间；在进行优化判断方面，创新性提出了用正切值代替角度计算来判断，节省了判断时间。

【技术实现步骤摘要】
一种基于多波束测深数据的Delaunay三角网快速生成方法
本专利技术涉及是多数据波束测深的方法，具体涉及是一种基于多波束测深数据的Delaunay三角网快速生成方法。
技术介绍
不规则三角网作为对二维平面域内任意离散点集构网的一种图形表现手段和分析工具，具有几何结构良好、数据结构简单的特点。和规则格网相比，不规则区域数据点分布适应能力更强，因此在地学、有限元分析等领域得到了广泛的应用。而应用的关键就在于搜索定位三角形和重构方法上，三角网因其构网的唯一性和合理性使得其被应用在构网的生成上。随着三角剖分法的广泛应用，近几十年来该算法得到了广大学者的关注和研究，根据Delaunay三角网构建过程的不同，把生成Delaunay三角网的算法分为逐点插入法、三角网生长法和分治法。分治法通过将点集区域递归划分构建子网，而后进行子网合并生成最终的三角网，构网效率较高但是算法设计复杂且在递归过程中占用较多的内存。逐点插入算法主要在生成初始三角网的基础上，以逐点插入的形式构建三角网，具有算法实现简单，占用内存少的特点，但是随着三角网的增多判断插入点所在的三角形就会占据大量的时间，因此构网效率并不高。三角网生长算法主要是基于基边寻找第三点构建三角网的过程，和分治算法相比实现简单，不用进行大量子网的递归生成与合并，和逐点插入法相比有着同样的缺点就是算法的实现效率不高特别是在离散点较多的时候，需要花费大量的时间检索合适的第三点。在Delaunay三角网的构网过程中需要大量的影响域判定的计算，因此对判定过程的简化会使整个构网效率大大提高。综上所述，为了提高构网作业的工作效率和保证构网的正确性，迫切需要开发一种改进的三角网逐点插入算法，以满足实际使用的需要。因此，这里根据多波束测深数据的特点采用的逐点插入算法考虑到搜索定位三角形和影响域判定的时间优化，以提高逐点插入算法构网的效率同时保证构网的正确性。
技术实现思路
专利技术的目的在于提供一种基于多波束测深数据的Delaunay三角网快速生成方法，已解决在Delaunay三角网的构网过程中需要大量的影响域判定的计算，从而减少了判定时间，提高构网作业的工作效率和保证构网的正确性。本专利技术一种基于多波束测深数据的Delaunay三角网快速生成方法，根据多波束测深数据每一PING数据在一条直线上的特性，每一PING数据进行关联插入，并用正切值代替角度计算来判断三角网优化过程，算法具体实现过程步骤如下：(1)建立超级三角形；建立一个超级三角形使得所有数据点都在该超级三角形内部遍历所有需要插入的点，找到Xmin，Xmax，Ymin，Ymax；即所有点都在这个矩形内部，所以超级三角形只需包含这个矩形即可；(2)从所有数据点中取出一PING；采用传统方法在整个三角网中搜寻包含第一个点的三角形，即对所有三角形逐个进行检测；(3)按照从里到外的顺序进行该点影响域的判断；从P所在的三角形开始，搜索该三角形的邻接三角形，并进行空外接圆检测；找到外接圆包含点P的所有三角形并删除这些三角形，形成一个包含P的多边形空腔，我们称之为Delaunay空腔；然后连接P与Delaunay腔的每一个顶点，形成新的Delaunay三角网格，完成局部优化；(4)沿着PING线的方向搜寻包含下一点的三角形；重复步骤(3)，完成该点的局部三角网更新；(5)重复步骤(4)，直至这一PING数据点处理完毕；(6)重复步骤(2)到步骤(5)，直到所有数据处理完毕；(7)待所有点都插入完毕后，将含有超级三角形顶点的三角形删除；(8)对整体三角形网络进行凹凸边缘检测，将所有边缘凹角通过补三角方法变成凸角。所述的步骤(3)所执行结构，按照从里到外的顺序进行该点影响域的判断，从P所在的三角形开始，搜索该三角形的邻接三角形，并进行空外接圆检测；找到外接圆包含点P的所有三角形并删除这些三角形，形成一个包含P的多边形空腔，称其为Delaunay空腔；连接P与Delaunay空腔的每一个顶点，形成新的Delaunay三角网格。所述的步骤(3)中通过对角正切值来计算判断该点的影响域采用的方法如下：(3.1)当两个角都是锐角或都是钝角时即Tanα*Tanβ>0，只需要判断Tanα的正负即可判断出α+β是否大于180°；(3.2)当一个锐角一个钝角时即Tanα*Tanβ<0，取临界条件，即两角和等于180度情况，假设α为锐角，β为钝角，由正切角函数关系Tanα+Tanβ＝0此时α+β＝180°；又已知正切函数是单调递增的，所以Tanα+Tanβ<0此时α+β<180°或是Tanα+Tanβ>0此时α+β>180°；(3.3)所以只需计算两个角的正切值即可判断出两角和是否大于180°，进而确定影响域。所述的步骤(4)所执行结构，首先判断上一点所在的Delaunay空腔中与PING线相交的边，具体方法如下：(4.1)首先选取Delaunay空腔相邻的两个顶点，进行边侧检验；(4.2)然后判断若在PING线同一侧，则进行下一组检测；若在PING线两侧，判断该顶点与PING线的夹角，从而确定是与PING线相交还是与PING线的反向延长线相交，直到找到与PING线相交的边；(4.3)最后找到与PING线相交的边后便可以沿着PING线方向顺序进行逐个检测，直到找到该点所在的三角形。本专利技术的有益效果在于：本专利技术一种基于多波束测深数据的Delaunay三角网快速生成方法，根据多数据波束测深数据的特点，每一PING数据进行关联插入相比于传统的逐个检测的方法，大大缩短了寻找时间；在进行优化判断方面，创新性提出了用正切值代替角度计算来判断，节省了判断时间。附图说明图1本专利技术中的数据插入顺序流程图；图2本专利技术中的三角形优化图；图3本专利技术中的影响域判定过程耗时对比表；图4本专利技术中的一点所在三角形判断顺序；图5本专利技术中的一点所在三角形判断流程图。具体实施方式结合附图和实施例对本专利技术进一步说明。本专利技术的目的在于提供一种三角网快速插入算法。经过对三角网逐点插入算法的研究，借鉴多种结构和思想的优点，结合多波束测深系统数据的特点，提出了本专利技术的方法，因为多波束测深系统的同一PING波束的数据点在一条直线上，取首末位置的两个点构成有向线段，记为PING线，将一PING的第一个点用传统的方法插入到已存在的三角网络中，在第一点的Delaunay腔中，找到与PING线相交的边，对第二个点进行边侧检验，顺着PING线方向，直到找到第二点所在三角形，后面的点以此类推，直到PING线上所有点都插入到三角网络中。三角形优化过程采用Delaunay腔法，本专利技术通过利用对角的正切值计算代替求出对角角度实现插入点影响域判定的过程，减少了判定时间。本专利技术采用逐点插入方式生成Delaunay三角网的算法的步骤如下：(1)建立初始三角网格。针对给定的点集V，找到一个包含该点集超级三角形R，称R为辅助窗口。(2)逐点插入。假设目前已经有一个Delaunay三角网格T，现在在它里面再插入一个点P，需要先找到点P所在的三角形。连接P点和P点所在三角形的三个顶点形成三个三角形，搜索该三角形的邻接三角形，并进行空外接圆检测。进行三角网络局部优本文档来自技高网...

【技术保护点】
1.一种基于多波束测深数据的Delaunay三角网快速生成方法，其特征在于：根据多波束测深数据每一PING数据在一条直线上的特性，每一PING数据进行关联插入，并用正切值代替角度计算来判断三角网优化过程，算法具体实现过程步骤如下：(1)建立超级三角形；建立一个超级三角形使得所有数据点都在该超级三角形内部遍历所有需要插入的点，找到Xmin，Xmax，Ymin，Ymax；即所有点都在这个矩形内部，所以超级三角形只需包含这个矩形即可；(2)从所有数据点中取出一PING；采用传统方法在整个三角网中搜寻包含第一个点的三角形，即对所有三角形逐个进行检测；(3)按照从里到外的顺序进行该点影响域的判断；从P所在的三角形开始，搜索该三角形的邻接三角形，并进行空外接圆检测；找到外接圆包含点P的所有三角形并删除这些三角形，形成一个包含P的多边形空腔，我们称之为Delaunay空腔；然后连接P与Delaunay腔的每一个顶点，形成新的Delaunay三角网格，完成局部优化；(4)沿着PING线的方向搜寻包含下一点的三角形；重复步骤(3)，完成该点的局部三角网更新；(5)重复步骤(4)，直至这一PING数据点处理完毕；(6)重复步骤(2)到步骤(5)，直到所有数据处理完毕；(7)待所有点都插入完毕后，将含有超级三角形顶点的三角形删除；(8)对整体三角形网络进行凹凸边缘检测，将所有边缘凹角通过补三角方法变成凸角。...

【技术特征摘要】
1.一种基于多波束测深数据的Delaunay三角网快速生成方法，其特征在于：根据多波束测深数据每一PING数据在一条直线上的特性，每一PING数据进行关联插入，并用正切值代替角度计算来判断三角网优化过程，算法具体实现过程步骤如下：(1)建立超级三角形；建立一个超级三角形使得所有数据点都在该超级三角形内部遍历所有需要插入的点，找到Xmin，Xmax，Ymin，Ymax；即所有点都在这个矩形内部，所以超级三角形只需包含这个矩形即可；(2)从所有数据点中取出一PING；采用传统方法在整个三角网中搜寻包含第一个点的三角形，即对所有三角形逐个进行检测；(3)按照从里到外的顺序进行该点影响域的判断；从P所在的三角形开始，搜索该三角形的邻接三角形，并进行空外接圆检测；找到外接圆包含点P的所有三角形并删除这些三角形，形成一个包含P的多边形空腔，我们称之为Delaunay空腔；然后连接P与Delaunay腔的每一个顶点，形成新的Delaunay三角网格，完成局部优化；(4)沿着PING线的方向搜寻包含下一点的三角形；重复步骤(3)，完成该点的局部三角网更新；(5)重复步骤(4)，直至这一PING数据点处理完毕；(6)重复步骤(2)到步骤(5)，直到所有数据处理完毕；(7)待所有点都插入完毕后，将含有超级三角形顶点的三角形删除；(8)对整体三角形网络进行凹凸边缘检测，将所有边缘凹角通过补三角方法变成凸角。2.根据权利要求1所述的一种基于多波束测深数据的Delaunay三角网快速生成方法，其特征在于：所述的步骤(3)所执行结构，按照从里到外的顺序进行该点影响域的判断，从P所在的三角形开始，搜索该三角形的邻接三角形，并进行空外接圆检测；找到外接圆包含点P的所有三角形并删除这些...

【专利技术属性】
技术研发人员：李树声，刘影，杜伟东，徐超，王珏，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙江,23