[发明专利]一种利用三角形进行地图网格化的划分算法

说明书

本发明公开了一种利用三角形进行地图网格化的划分算法，采用Douglas‑Peucker算法提取区域边界，用约束Delaunay三角剖分算法得到三角区域格网，所述采用Douglas‑Peucker算法提取区域边界，通过多边形来拟合曲线，具体的算法步骤如下：S1、连接曲线首尾两点A、B形成一条直线AB；S2、计算曲线上离该直线段距离最大的点C，计算其与AB的距离d；S3、比较该距离与预先给定的阈值threshold的大小，S4、当所有曲线都处理完毕后，依次连接各个分割点形成折线，作为原曲线的近似，本发明结构科学合理，使用安全方便，采用Douglas‑Peucker算法提取区域边界，用约束Delaunay三角剖分算法得到三角区域格网，可以快速准确的对地图区域进行合理的三角形网格划分。

技术领域

本发明涉及地理信息技术领域，具体为一种利用三角形进行地图网格化的划分算法。

背景技术

地理信息是地理数据所蕴含和表达的地理含义，是与地理环境要素有关的物质的数量、质量、性质、分布特征、联系和规律的数字、文字、图像和图形等的总称；

在地理信息数据处理过程中，经常会用到网格化处理的方法，如何通过合理的方法对给定的区域进行网格化，让网格能够充分的逼近区域地图，是在网格化处理中重要的处理方法。

发明内容

本发明提供一种利用三角形进行地图网格化的划分算法，可以有效解决上述背景技术中提出如何通过合理的方法对给定的区域进行网格化，让网格能够充分的逼近区域地图的问题。

为实现上述目的，本发明提供如下技术方案：一种利用三角形进行地图网格化的划分算法，采用Douglas-Peucker算法提取区域边界，用约束Delaunay三角剖分算法得到三角区域格网；

所述采用Douglas-Peucker算法提取区域边界，通过多边形来拟合曲线，具体的算法步骤如下：

S1、连接曲线首尾两点A、B形成一条直线AB；

S2、计算曲线上离该直线段距离最大的点C，计算其与AB的距离d；

S3、比较该距离与预先给定的阈值threshold的大小，

S4、当所有曲线都处理完毕后，依次连接各个分割点形成折线，作为原曲线的近似。

根据上述技术方案，所述S3中阈值threshold小于threshold，则以该直线作为曲线的近似，该段曲线处理完毕；

阈值threshold距离大于阈值，则用点C将曲线分为两段AC和BC，并分别对两段曲线进行步骤S1-S3的处理。

根据上述技术方案，所述约束Delaunay三角剖分算法得到三角区域格网，在完成了边界提取后，就可以开始向计算区域内插入点进行Delaunay三角剖分；

所述Delaunay三角剖分具体算法步骤如下：

A1、设置一个面积的控制量以便控制剖分的程度，它可以方便的控制三角网格的大小，计算三角形数组中的所有三角形面积大小；

A2、选择面积较大的三角形的最长边作为新的插入点，并将该点加到顶点数组后面；

A3、搜索整个三角形数组，找出所有插入点落在三角形外接圆中的三角形，计入优化三角形数组，并将三角形的边记入边数组，并且不能有重复边，若有重复边就要从数组中删除该边；

A4、连接插入点与边数组中每条边的两个顶点，形成新的三角形加入三角形数组；

A5、删除优化三角形数组中的三角形，更新整个三角形数组，完成一次点的插入；

A6、重复A2、A3、A4直到三角形的面积符合控制要求，结束点的插入，这样就完成了Delaunay剖分。