[发明专利]一种基于Delaunay三角网的栅格地图矢量化方法

权利要求书

1.一种基于Delaunay三角网的栅格地图矢量化方法，包括以下步骤：

第一步、对单要素栅格图进行黑白二值化；

第二步、寻找所述单要素栅格图中的线划要素；

第三步、提取各线划要素对应的栅格图的边界栅格，从而获得各线划要素对应的边界栅格图；

第四步、基于各边界栅格图的边界栅格中心点分别生成Delaunay三角网；

第五步、剔除不在所述线划要素范围内的干扰三角形；

第六步、利用Delaunay三角网三角形间的拓扑关系提取线划要素骨架线，即得到单要素栅格图的矢量化结果。

2.根据权利要求1所述的基于Delaunay三角网的栅格地图矢量化方法，其特征在于：如果待矢量化地图为多要素的彩色地图，则首先通过分色，得到多个单色的单要素栅格图，再分别对单要素栅格图执行第一步到第六步。

3.根据权利要求1所述的基于Delaunay三角网的栅格地图矢量化方法，其特征在于，第二步中，使用像元邻域分析的标记算法寻找所述单要素栅格图中的线划要素，具体方法如下：

1)顺序遍历单要素栅格图中的像元，对所有不为0的像元根据其扫描顺序依次赋值；

2)再次顺序遍历重新赋值后的单要素栅格图中的像元，对所有不为0的像素与其八邻域中不为0的像素进行比较，若其八邻域范围内所有标记值的最小值小于其本身的标记值，则修改其标记值为撰述最小值；

3)反向遍历步骤2)处理后的单要素栅格图中的像元，处理过程同步骤2)；

4)再次顺序遍历单要素栅格图中的像元，统计具有不同标记的像元种类数；

5)重复步骤2)-4)直到统计的像元种类数与上一次统计的像元种类数一致，即完成像元的分组，具有相同值的像元属于同一线划要素。

4.根据权利要求1所述的基于Delaunay三角网的栅格地图矢量化方法，其特征在于，第三步中，使用基于数学形态学的腐蚀算法提取栅格图中的边界栅格，具体方法如下：

对分组后的栅格图用3*3的腐蚀算子进行腐蚀运算，得到腐蚀后的栅格集合；用所述栅格图减去腐蚀后得到的栅格集合就得到了栅格图中的边界栅格。

5.根据权利要求1所述的基于Delaunay三角网的栅格地图矢量化方法，其特征在于，第五步中，剔除不在所述线划要素范围内的干扰三角形的方法为：判断各Delaunay三角网中各三角形，若三角形位于所述线划要素范围内，则保留该三角形，否则准备删除；判断所述准备删除的三角形各边，若所述边位于所述线划要素范围内，则保留该边，否则剔除。

6.根据权利要求5所述的基于Delaunay三角网的栅格地图矢量化方法，其特征在于，判断所述三角形位于线划要素范围内的具体方法如下：如果所述三角形内任一点位于线划要素范围内，则该三角形位于线划要素范围内。

7.根据权利要求5所述的基于Delaunay三角网的栅格地图矢量化方法，其特征在于，判断准备删除的三角形各边是否位于线划要素范围内的具体方法如下：若所述边的两端点的像元八邻域连通值与该线划要素相同，则该边位于线划要素范围内。