一种交叉口密集区域的骨架线提取方法技术

本发明专利技术公开了一种交叉口密集区域的骨架线提取方法，包括如下步骤：S1、建立多边形的约束Delaunay三角网，对多边形中的长边进行加密，并更新约束Delaunay三角网；S2、基于更新后的约束Delaunay三角网，提取多边形骨架线；S3、基于提取的骨架线，对交叉口密集区进行识别；S4、对识别出的交叉口密集区进行骨架线调整优化。优点是：通过该方法提取骨架线后能够保持原有区域的整体特征，自然延展且符合人类视觉要求；该方法在提取骨架线的过程中不受分支复杂边界及排列结构的干扰，能够处理此区域的骨架线抖动，所提骨架线不会存在较大程度的扭曲，也不会丢失了此区域的整体结构。

A Skeleton Line Extraction Method for Dense Intersection Areas

【技术实现步骤摘要】
一种交叉口密集区域的骨架线提取方法
本专利技术涉及地图制图
，尤其涉及一种交叉口密集区域的骨架线提取方法。
技术介绍
骨架线提取是实现多边形降维、分裂融解等地图综合操作的关键步骤，骨架线的提取要顾及多边形的形状特征，能够概括多边形的主体结构及延展特性，在符合制图规范的同时要满足人的视觉认知。目前骨架线提取最常用的方法是DeLucia等在1987年提出了一种基于边界约束Delaunay三角网的骨架线提取方法，后续大量研究人员基于该方法进行改进。如陈涛等人的论文“多边形骨架线与形心自动搜寻算法研究”中利用Delaunay三角网分析多边形的形态结构，并在此基础上根据分支部分的面积，提取主骨架线和多边形的形心；专利“一种狭长图斑分裂线的处理方法”(申请号：CN201710973188.1)公开了一种狭长图斑分裂线的处理方法，其通过Delaunay三角网狭长图斑进行三角形划分,并提取所述狭长图斑中的骨架线，同时对骨架线上的抖动和去悬挂线处理，从而获取与原图形的形状更加一致的自然光滑的图斑分裂线。现有基于Delaunay三角网骨架线提取方法用于道路或水网的交叉口的骨架线提取时，是以每交叉口为处理单元，通过设定方向阈值及分支结点距离阈值，来获得交叉口区域的骨架线，当然而，当交叉口情况较为简单，如单一T字形或者十字型时，提取效果较好，如专利“一种狭长图斑分裂线的处理方法”即如此。但当交叉口复杂汇聚时，现有方法会因无法准确剥离单个分支汇聚区，导致所提骨架线不够准确。此外，交叉口密集区骨架线的形状复杂，以每一个分支汇聚区为单元进行处理，无法顾及各个分支之间相互关联形成的区域整体特征，导致区域整体结构被破坏。
技术实现思路
本专利技术的目的在于提供一种交叉口密集区域的骨架线提取方法，从而解决现有技术中存在的前述问题。为了实现上述目的，本专利技术采用的技术方案如下：一种交叉口密集区域的骨架线提取方法，包括如下步骤，S1、建立多边形的约束Delaunay三角网，对多边形中的长边进行加密，并更新约束Delaunay三角网；S2、基于更新后的约束Delaunay三角网，提取多边形骨架线；S3、基于提取的骨架线，对交叉口密集区进行识别；S4、对识别出的交叉口密集区进行骨架线调整优化。2、根据权利要求1所述的交叉口密集区域的骨架线提取方法，其特征在于：步骤S1还包括如下内容，S101、识别约束Delaunay三角网中III类三角形中的钝角三角形，并设定最小角阈值，将最小角小于最小角阈值的钝角三角形进行标记、并放入三角形集合；S102、识别三角形集合中所有三角形的最长边，并查找与这些最长边邻接、且有一个边是边界的II类三角形；S103、将所有最长边为边界边的II三角形挑选出来、标记为待加密三角形，且将其相应的最长边标记为待加密局部长边；S104、将与待加密三角形关联的三角形集合中III类三角形的最短边设置为加密步长，并将待加密三角形和所述待加密局部长边按加密步长进行加密。优选的，步骤S3包括如下内容，S301、为初始骨架线构建点、线拓扑，记录骨架线拓扑中的各个结点所关联的关联弧段数量；S302、选取关联弧段数量为三的结点，该结点为三叉口结点，其所在区域为三叉口区域；S303、计算三叉口区域的局部近似宽度；S304、将首尾端点均为三叉口结点的弧段定义为连接弧段；分别计算连接弧段首结点和尾结点所在区域的局部近似宽度，并将其中的较大值作为连接弧段所在区域局部近似宽度；S305、将连接弧段的首结点和尾结点间去除III类三角形内部长度的部分记为连接弧段的有效长度；并判断连接弧段的有效长度是否小于连接弧段所在区域局部近似宽度，若是，则所述连接弧段的首结点和尾结点相互关联，该连接弧段为关联弧段；若否，则所述连接弧段的首结点和尾结点不关联，该连接弧段为关联弧段；S306、将步骤S305中所有的关联弧段放入关联弧段集合中；并依次从关联弧段集合中选取一个关联弧段，记为第一关联弧段，并将第一关联弧段的首结点和尾结点作为追踪结点，将第一关联弧段的首结点或尾结点存在的关联弧段记为第二关联弧段，并将第二关联弧段记录到关联弧段集合中；待各个关联弧段都检测完毕，对各关联弧段进行聚类扩展，得到交叉口聚集结果。优选的，步骤S4包括如下内容，S401、利用CRITIC方法对各连接弧段的长度、连通度、接近度以及中介度进行加权，得到连接弧段的重要度；S402、选取一个仅有一侧存在关联弧段的三叉口结点作为stroke连接追踪起始结点，并将该三叉口结点的连接弧段作为追踪弧段，将该追踪弧段另一侧的结点作为追踪结点；S403、将追踪结点的连接弧段作为stroke连接候选集，并计算各个连接弧段的重要性；S404、将重要性较大的连接弧段优先与上一连接弧段连接构成stroke；S405、以步骤S403和步骤S404为基础，继续向上追踪计算stroke连接，直至不存在追踪结点的连接弧段，则单条stroke连接结束；S406、探索现有stroke连接的分支连接弧段，直至所有交叉口密集区内所有的连接弧段均已计算，stroke连接计算结束；S407、基于计算得到的交叉密集区stroke连接，以stroke连接上的两个连接弧段为基准弧段，取III类三角形边上的两点连线的中点作为三叉口结点的调整位置，将三叉口结点关联的三个连接弧段与该中点进行连接，完成各个三叉口区域的骨架线调整。优选的所述交叉口结点为III类三角形的中心点。优选的，所述三叉口区域的局部近似宽度为三叉口对应的III类三角形的内部骨架线中最大值的两倍。本专利技术的有益效果是：1、区域通过该方法提取骨架线后能够保持原有区域的整体特征，自然延展且符合人类视觉要求。2、实现交叉口复杂密集区域的识别判定和骨架线提取，同时依据stroke连接具备的良好延续性特征进行提取骨架线的优化，使最终骨架线结果更加科学准确。3、该方法在提取骨架线的过程中不受分支复杂边界及排列结构的干扰，能够处理此区域的骨架线抖动，所提骨架线不会存在较大程度的扭曲，也不会丢失了此区域的整体结构。附图说明图1是本专利技术实施例中骨架线提取方法的流程图；图2是本专利技术实施例中步骤S1的流程图；图3是本专利技术实施例中拉长的III类三角形的示意图；图4是本专利技术实施例中一边为边界的II类三角形的示意图；图5是本专利技术实施例中具有待加密局部长边三角形的示意图；图6是本专利技术实施例中加密步长进行加密的三角形的示意图；图7是本专利技术实施例中步骤S3的流程图；图8是本专利技术实施例中三叉口结点的示意图；图9是本专利技术实施例中连接弧段的示意图；图10是本专利技术实施例中三叉口聚集结果示意图；图11是本专利技术实施例中步骤S4的流程图；图12是本专利技术实施例中三叉口区域骨架线调整前的示意图；图13是本专利技术实施例中三叉口区域骨架线调整后的示意图；图14是本专利技术实施例中的原始骨架线的示意图；图15是本专利技术实施例中使用现有技术方法提取骨架线的示意图；图16是本专利技术实施例中使用本专利技术方法提取骨架线的示意图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本专利技术，并不用于限定本专利技术。如图1所示，本专利技术提供了一种交叉口密集区域的骨架线提取方法，本文档来自技高网...

【技术保护点】
1.一种交叉口密集区域的骨架线提取方法，其特征在于：包括如下步骤，S1、建立多边形的约束Delaunay三角网，对多边形中的长边进行加密，并更新约束Delaunay三角网；S2、基于更新后的约束Delaunay三角网，提取多边形骨架线；S3、基于提取的骨架线，对交叉口密集区进行识别；S4、对识别出的交叉口密集区进行骨架线调整优化。

【技术特征摘要】
1.一种交叉口密集区域的骨架线提取方法，其特征在于：包括如下步骤，S1、建立多边形的约束Delaunay三角网，对多边形中的长边进行加密，并更新约束Delaunay三角网；S2、基于更新后的约束Delaunay三角网，提取多边形骨架线；S3、基于提取的骨架线，对交叉口密集区进行识别；S4、对识别出的交叉口密集区进行骨架线调整优化。2.根据权利要求1所述的交叉口密集区域的骨架线提取方法，其特征在于：步骤S1还包括如下内容，S101、识别约束Delaunay三角网中III类三角形中的钝角三角形，并设定最小角阈值，将最小角小于最小角阈值的钝角三角形进行标记、并放入三角形集合；S102、识别三角形集合中所有三角形的最长边，并查找与这些最长边邻接、且有一个边是边界的II类三角形；S103、将所有最长边为边界边的II三角形挑选出来、标记为待加密三角形，且将其相应的最长边标记为待加密局部长边；S104、将与待加密三角形关联的三角形集合中III类三角形的最短边设置为加密步长，并将待加密三角形和所述待加密局部长边按加密步长进行加密。3.根据权利要求1所述的交叉口密集区域的骨架线提取方法，其特征在于：步骤S3包括如下内容，S301、为初始骨架线构建点、线拓扑，记录骨架线拓扑中的各个结点所关联的关联弧段数量；S302、选取关联弧段数量为三的结点，该结点为三叉口结点，其所在区域为三叉口区域；S303、计算三叉口区域的局部近似宽度；S304、将首尾端点均为三叉口结点的弧段定义为连接弧段；分别计算连接弧段首结点和尾结点所在区域的局部近似宽度，并将其中的较大值作为连接弧段所在区域局部近似宽度；S305、将连接弧段的首结点和尾结点间去除III类三角形内部长度的部分记为连接弧段的有效长度；并判断连接弧段的有效长度是否小于连接弧段所在区域局部近似宽度，若是，则所述连接弧段的首结点和尾结点相互关联，该连接弧段为关联弧段；若否，则所述连接弧段的首结点和尾结...

【专利技术属性】
技术研发人员：李成名，吴伟，殷勇，武鹏达，郭沛沛，刘晓丽，
申请(专利权)人：中国测绘科学研究院，
类型：发明
国别省市：北京,11