CN112614206B - 道路中心线提取方法和装置 - Google Patents

Abstract

本发明实施例提供一种道路中心线提取方法和装置，属于计算机技术领域。该方法包括：读取双线道路路网的数据；在所述双线道路路网的两侧形成缓冲区，以生成道路几何图形；根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段；基于均质化网格，对所述道路几何图形进行切分，以得到多个网格；针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线；对所述初版道路中心线进行简化处理，得到终版道路中心线。该道路中心线提取方法和装置可以使得提取道路中心线准确、高效。

Description

道路中心线提取方法和装置

技术领域

本发明涉及计算机技术领域，具体地涉及一种道路中心线提取方法和装置。

背景技术

现有技术中，城市规划中通常是以道路中心线路网为基础，进行城市单元地块划分以及空间句法计算。因此，在城市规划的大数据分析领域中，中心线的提取尤为重要。

对于较为简单的道路网，中心线的提取比较简单，提取效率也高。但是，对于大规模复杂的道路网，由于涉及的道路较多，几何形体相互嵌套，交叉路口多，中心线抽取场景相对复杂，导致计算量急剧增大。同时，传统的中心线提取方法无法将交叉路口简化，会出现多余节点，该多余节点对整体道路网络影响较大，难以用到真实项目中。因此，传统的道路中心线提取方式没有办法高效率地进行提取。

发明内容

本发明实施例的目的是提供一种道路中心线提取方法和装置，该道路中心线提取方法和装置可以使得提取道路中心线准确、高效。

为了实现上述目的，本发明实施例提供一种道路中心线提取方法，该方法包括：读取双线道路路网的数据；在所述双线道路路网的两侧形成缓冲区，以生成道路几何图形；根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段；基于均质化网格，对所述道路几何图形进行切分，以得到多个网格；针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线；对所述初版道路中心线进行简化处理，得到终版道路中心线。

进一步的，所述根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段包括：在所述道路几何图形的线框上，生成多个间距为预设距离的点；根据所述多个间距为预设距离的点，生成泰森多边形，以得到基于泰森多边形原理的多个线段。

进一步的，所述对所述初版道路中心线进行简化处理，得到终版道路中心线包括：删除所述初版道路中心线中的断头线；根据道格拉斯-普克算法，对删除断头线之后的初版道路中心线进行重组；获取重组后的初版道路中心线中长度小于预设长度的线段；针对每个长度小于预设长度的线段，将所述长度小于预设长度的线段的中点和与所述长度小于预设长度的线段连接的线段的外端点连接，所述外端点是不在所述长度小于预设长度的线段上的端点；删除所述长度小于预设长度的线段以及与所述长度小于预设长度的线段连接的线段，以得到所述终版道路中心线。

进一步的，在生成道路几何图形之后，该方法还包括：删除所述道路几何图形中，面积小于预设面积或宽度小于预设宽度的图块；根据删除图块后的道路几何图形，执行生成泰森多边形的步骤。

进一步的，针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线包括：针对所述多个网格中的每个网格，提取两个端点都在所述道路几何图形内的线段；对所提取的线段，删除与所述道路几何图形有交点的线段，并进行组合，以生成初版道路中心线。

本发明实施例还提供一种道路中心线提取装置，该装置包括处理单元，该处理单元用于：读取双线道路路网的数据；在所述双线道路路网的两侧形成缓冲区，以生成道路几何图形；根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段；基于均质化网格，对所述道路几何图形进行切分，以得到多个网格；针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线；对所述初版道路中心线进行简化处理，得到终版道路中心线。

进一步的，所述处理单元还用于：在所述道路几何图形的线框上，生成多个间距为预设距离的点；根据所述多个间距为预设距离的点，生成泰森多边形，以得到基于泰森多边形原理的多个线段。

进一步的，所述处理单元还用于：删除所述初版道路中心线中的断头线；根据道格拉斯-普克算法，对删除断头线之后的初版道路中心线进行重组；获取重组后的初版道路中心线中长度小于预设长度的线段；针对每个长度小于预设长度的线段，将所述长度小于预设长度的线段的中点和与所述长度小于预设长度的线段连接的线段的外端点连接，所述外端点是不在所述长度小于预设长度的线段上的端点；删除所述长度小于预设长度的线段以及与所述长度小于预设长度的线段连接的线段，以得到所述终版道路中心线。

进一步的，所述处理单元还用于：删除所述道路几何图形中，面积小于预设面积或宽度小于预设宽度的图块；根据删除图块后的道路几何图形，执行生成泰森多边形的步骤。

进一步的，所述处理单元还用于：针对所述多个网格中的每个网格，提取两个端点都在所述道路几何图形内的线段；对所提取的线段，删除与所述道路几何图形有交点的线段，并进行组合，以生成初版道路中心线。

通过上述技术方案，读取双线道路路网的数据；在所述双线道路路网的两侧形成缓冲区，以生成道路几何图形；根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段；基于均质化网格，对所述道路几何图形进行切分，以得到多个网格；针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线；对所述初版道路中心线进行简化处理，得到终版道路中心线。该道路中心线提取方法和装置可以使得提取道路中心线准确、高效。

本发明实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。

附图说明

附图是用来提供对本发明实施例的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明实施例，但并不构成对本发明实施例的限制。在附图中：

图1是本发明一实施例提供的道路中心线提取方法的流程图；

图2是本发明一实施例提供的基于泰森多边形原理的多个线段获得方法的流程图；

图3是本发明一实施例提供的提取与道路几何图形无交点的线段的方法的流程图；

图4是本发明一实施例提供的简化处理初版道路中心线的方法的流程图；

图5A-5J是本发明一实施例提供的道路中心线提取过程中各步骤结果的可视化示意图。

具体实施方式

以下结合附图对本发明实施例的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明实施例，并不用于限制本发明实施例。

图1是本发明一实施例提供的道路中心线提取方法的流程图。如图1所示，该方法包括：

步骤S11，读取双线道路路网的数据；

例如，读取双线道路路网的SHP数据或Json数据，其内容主要包括双线道路路网的矢量数据。可以通过在Python中运行以下部分代码实现：

df＝readData(file_path,city)

df_geom＝multilines2list(list(df.geometry))

步骤S12，在所述双线道路路网的两侧形成缓冲区，以生成道路几何图形；

例如，在双线道路路网两侧生成宽为15米(不限于此)的缓冲区，从而得到几何图形的双线道路路网。可以通过在Python中运行以下部分代码实现：

为了保证中心线提取的高效率及准确性，可以删除所述道路几何图形中，面积小于预设面积(例如200m²，但不限于此)或宽度小于预设宽度(例如25m，但不限于此)的图块，然后再执行后续的生成泰森多边形的步骤。可以通过在Python中运行以下部分代码实现：

#Create the Voronoi diagram

vor＝Voronoi(all_points)

all_vertices＝vor.vertices

print("Finish creating Thiessen Polygon...")

calculate_running_time(start_time)

步骤S13，根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段；

例如，可以包括以下步骤：

步骤S21，在所述道路几何图形的线框上，生成多个间距为预设距离的点；

例如，可以通过在Python中运行以下部分代码实现：

all_points＝polygons_to_point_series(buffers,point_dist)

all_points＝np.array(all_points)

print("Finish extracting points from the buffer polygon...")

calculate_running_time(start_time)

步骤S22，根据所述多个间距为预设距离的点，生成泰森多边形，以得到基于泰森多边形原理的多个线段。

例如，可以通过在Python中运行以下部分代码实现：

#Create the Voronoi diagram

vor＝Voronoi(all_points)

all_vertices＝vor.vertices

print("Finish creating Thiessen Polygon...")

calculate_running_time(start_time)

步骤S14，基于均质化网格，对所述道路几何图形进行切分，以得到多个网格；

例如，对道路几何图形进行500米*500米(但不限于此)的正方形密铺。例如，分割成20个正方形。因此，基于不同正方形，对道路几何图形进行切割。同时，为处于正方形内部的线段进行编号。例如，处于正方形01的，则为01编号，可以通过在Python中运行以下部分代码实现：

步骤S15，针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线；

例如，可以包括以下步骤：

步骤S31，针对所述多个网格中的每个网格，提取两个端点都在所述道路几何图形内的线段；

步骤S32，对所提取的线段，删除与所述道路几何图形有交点的线段，并进行组合，以生成初版道路中心线。

例如，可以通过在Python中运行以下部分代码实现：

ridges_on_buffer＝get_ridges_inside_buffer(vor.ridge_vertices,all_vertices,all_grids)

print("Finish filtering ridges...")

calculate_running_time(start_time)

步骤S16，对所述初版道路中心线进行简化处理，得到终版道路中心线。

例如，可以包括以下步骤：

步骤S41，删除所述初版道路中心线中的断头线；

例如，初版道路中心线会留有一部分的位于道路几何图形内部的线段，但本身不可作为路网，因此需要利用删除断头线的算法，对其进行删除，可以通过在Python中运行以下代码实现：

#Clean cul-de-sacs among the ridges

o_list＝[x[0]for x in ridges_on_buffer]

d_list＝[x[1]for x in ridges_on_buffer]

non_dangles_id＝get_non_dangles_id(o_list,d_list)

ridges_on_buffer＝[ridges_on_buffer[i]for i in non_dangles_id]

print("Finish deleting dangles...")

calculate_running_time(start_time)

步骤S42，根据道格拉斯-普克算法，对删除断头线之后的初版道路中心线进行重组；

例如，根据道格拉斯-普克算法，可以在不改变主要方向的前提下，对原本过于曲折的线段进行重组，可以通过在Python中运行以下部分代码实现：

#Simplify roads

roads_seperated＝[np.asarray(x)for x in allRoadList]

roads_single＝[rdp(x,epsilon＝rdp_eps)for x in roads_seperated]

print("Finish simplifying roads...")

calculate_running_time(start_time)

步骤S43，获取重组后的初版道路中心线中长度小于预设长度的线段；

步骤S44，针对每个长度小于预设长度的线段，将所述长度小于预设长度的线段的中点和与所述长度小于预设长度的线段连接的线段的外端点连接，所述外端点是不在所述长度小于预设长度的线段上的端点；

步骤S45，删除所述长度小于预设长度的线段以及与所述长度小于预设长度的线段连接的线段，以得到所述终版道路中心线。

例如，初版道路中心线中，在交叉路口中，会出现冗余线段，不可直接使用，因此可以对所有线段进行判断，选出长度小于预设长度(例如20米，但不限于此)的线段，记录与该线段相互连接线段的端点(不在该线段上)的序号，同时求出该线段中点。将记录后的序号全部赋值于该线段中点，以连接中点与周围端点，形成新的线段，并删除原有线段，可以通过在Python中运行以下部分代码实现：

#Clean too-short roads

roads_exploded＝explodeLineList(roads_single)

final_roads＝delete_short_roads(roads_exploded,threshold,margin_pos)

print("Finish cleaning short roads...")

calculate_running_time(start_time)

针对上述每一步，都可以使用例如QGIS进行可视化，下文展示每一步执行后的道路路网和中心线的变化：

如图5A所示，示出了执行读取双线道路路网的数据的步骤后的结果；

如图5B所示，示出了执行建立缓冲区的步骤后的结果；

如图5C所示，示出了执行剔除小图块的步骤后的结果；

如图5D所示，示出了执行按照特定距离生成点的步骤后的结果；

如图5E所示，示出了执行生成泰森多边形的步骤后的结果；

如图5F所示，示出了执行道路几何图形切分的步骤后的结果；

如图5G所示，示出了执行提取与所述道路几何图形无交点的线段的步骤后的结果；

如图5H所示，示出了执行剔除断头线的步骤后的结果；

如图5I所示，示出了执行初版道路中心线进行重组的步骤后的结果；

如图5J所示，示出了执行删除冗余线段的步骤后的结果。

在测试案例中，有7000条路网，基于现有的方法，需要算4.5个小时左右。但在用了本发明的方法，可将计算时间缩短至2分钟。

本发明还提供一种道路中心线提取装置，该装置包括处理单元，该处理单元用于：读取双线道路路网的数据；在所述双线道路路网的两侧形成缓冲区，以生成道路几何图形；根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段；基于均质化网格，对所述道路几何图形进行切分，以得到多个网格；针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线；对所述初版道路中心线进行简化处理，得到终版道路中心线。

进一步的，所述处理单元还用于：在所述道路几何图形的线框上，生成多个间距为预设距离的点；根据所述多个间距为预设距离的点，生成泰森多边形，以得到基于泰森多边形原理的多个线段。

进一步的，所述处理单元还用于：删除所述初版道路中心线中的断头线；根据道格拉斯-普克算法，对删除断头线之后的初版道路中心线进行重组；获取重组后的初版道路中心线中长度小于预设长度的线段；针对每个长度小于预设长度的线段，将所述长度小于预设长度的线段的中点和与所述长度小于预设长度的线段连接的线段的外端点连接，所述外端点是不在所述长度小于预设长度的线段上的端点；删除所述长度小于预设长度的线段以及与所述长度小于预设长度的线段连接的线段，以得到所述终版道路中心线。

进一步的，所述处理单元还用于：删除所述道路几何图形中，面积小于预设面积或宽度小于预设宽度的图块；根据删除图块后的道路几何图形，执行生成泰森多边形的步骤。

进一步的，所述处理单元还用于：针对所述多个网格中的每个网格，提取两个端点都在所述道路几何图形内的线段；对所提取的线段，删除与所述道路几何图形有交点的线段，并进行组合，以生成初版道路中心线。

上文所述的道路中心线提取装置的实施例与上文所述的道路中心线提取方法的实施例类似，在此不再赘述。

所述道路中心线提取装置包括处理器和存储器，上述处理单元等均作为程序单元存储在存储器中，由处理器执行存储在存储器中的上述程序单元来实现相应的功能。

处理器中包含内核，由内核去存储器中调取相应的程序单元。内核可以设置一个或以上，通过调整内核参数来提取道路中心线。

存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)，存储器包括至少一个存储芯片。

本发明实施例提供了一种存储介质，其上存储有程序，该程序被处理器执行时实现所述道路中心线提取方法。

本发明实施例提供了一种处理器，所述处理器用于运行程序，其中，所述程序运行时执行所述道路中心线提取方法。

本发明实施例提供了一种设备，设备包括处理器、存储器及存储在存储器上并可在处理器上运行的程序，处理器执行程序时实现以下步骤：

读取双线道路路网的数据；在所述双线道路路网的两侧形成缓冲区，以生成道路几何图形；根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段；基于均质化网格，对所述道路几何图形进行切分，以得到多个网格；针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线；对所述初版道路中心线进行简化处理，得到终版道路中心线。

进一步的，所述根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段包括：在所述道路几何图形的线框上，生成多个间距为预设距离的点；根据所述多个间距为预设距离的点，生成泰森多边形，以得到基于泰森多边形原理的多个线段。

进一步的，所述对所述初版道路中心线进行简化处理，得到终版道路中心线包括：删除所述初版道路中心线中的断头线；根据道格拉斯-普克算法，对删除断头线之后的初版道路中心线进行重组；获取重组后的初版道路中心线中长度小于预设长度的线段；针对每个长度小于预设长度的线段，将所述长度小于预设长度的线段的中点和与所述长度小于预设长度的线段连接的线段的外端点连接，所述外端点是不在所述长度小于预设长度的线段上的端点；删除所述长度小于预设长度的线段以及与所述长度小于预设长度的线段连接的线段，以得到所述终版道路中心线。

进一步的，在生成道路几何图形之后，该方法还包括：删除所述道路几何图形中，面积小于预设面积或宽度小于预设宽度的图块；根据删除图块后的道路几何图形，执行生成泰森多边形的步骤。

进一步的，针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线包括：针对所述多个网格中的每个网格，提取两个端点都在所述道路几何图形内的线段；对所提取的线段，删除与所述道路几何图形有交点的线段，并进行组合，以生成初版道路中心线。

本文中的设备可以是服务器、PC、PAD、手机等。

本申请还提供了一种计算机程序产品，当在数据处理设备上执行时，适于执行初始化有如下方法步骤的程序：

读取双线道路路网的数据；在所述双线道路路网的两侧形成缓冲区，以生成道路几何图形；根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段；基于均质化网格，对所述道路几何图形进行切分，以得到多个网格；针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线；对所述初版道路中心线进行简化处理，得到终版道路中心线。

进一步的，所述根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段包括：在所述道路几何图形的线框上，生成多个间距为预设距离的点；根据所述多个间距为预设距离的点，生成泰森多边形，以得到基于泰森多边形原理的多个线段。

进一步的，所述对所述初版道路中心线进行简化处理，得到终版道路中心线包括：删除所述初版道路中心线中的断头线；根据道格拉斯-普克算法，对删除断头线之后的初版道路中心线进行重组；获取重组后的初版道路中心线中长度小于预设长度的线段；针对每个长度小于预设长度的线段，将所述长度小于预设长度的线段的中点和与所述长度小于预设长度的线段连接的线段的外端点连接，所述外端点是不在所述长度小于预设长度的线段上的端点；删除所述长度小于预设长度的线段以及与所述长度小于预设长度的线段连接的线段，以得到所述终版道路中心线。

进一步的，在生成道路几何图形之后，该方法还包括：删除所述道路几何图形中，面积小于预设面积或宽度小于预设宽度的图块；根据删除图块后的道路几何图形，执行生成泰森多边形的步骤。

进一步的，针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线包括：针对所述多个网格中的每个网格，提取两个端点都在所述道路几何图形内的线段；对所提取的线段，删除与所述道路几何图形有交点的线段，并进行组合，以生成初版道路中心线。

本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。

存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。

计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。

还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。

以上仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (8)

1.一种道路中心线提取方法，其特征在于，该方法包括：

读取双线道路路网的数据；

在所述双线道路路网的两侧形成缓冲区，以生成道路几何图形；

根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段；

基于均质化网格，对所述道路几何图形进行切分，以得到多个网格；

针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线；

对所述初版道路中心线进行简化处理，得到终版道路中心线；

所述对所述初版道路中心线进行简化处理，得到终版道路中心线包括：

删除所述初版道路中心线中的断头线；

根据道格拉斯-普克算法，对删除断头线之后的初版道路中心线进行重组；

获取重组后的初版道路中心线中长度小于预设长度的线段；

针对每个长度小于预设长度的线段，将所述长度小于预设长度的线段的中点和与所述长度小于预设长度的线段连接的线段的外端点连接，所述外端点是不在所述长度小于预设长度的线段上的端点；

删除所述长度小于预设长度的线段以及与所述长度小于预设长度的线段连接的线段，以得到所述终版道路中心线。

2.根据权利要求1所述的道路中心线提取方法，其特征在于，所述根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段包括：

在所述道路几何图形的线框上，生成多个间距为预设距离的点；

根据所述多个间距为预设距离的点，生成泰森多边形，以得到基于泰森多边形原理的多个线段。

3.根据权利要求1所述的道路中心线提取方法，其特征在于，在生成道路几何图形之后，该方法还包括：

删除所述道路几何图形中，面积小于预设面积或宽度小于预设宽度的图块；

根据删除图块后的道路几何图形，执行生成泰森多边形的步骤。

4.根据权利要求1所述的道路中心线提取方法，其特征在于，针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线包括：

针对所述多个网格中的每个网格，提取两个端点都在所述道路几何图形内的线段；

对所提取的线段，删除与所述道路几何图形有交点的线段，并进行组合，以生成初版道路中心线。

5.一种道路中心线提取装置，其特征在于，该装置包括处理单元，该处理单元用于：

读取双线道路路网的数据；

在所述双线道路路网的两侧形成缓冲区，以生成道路几何图形；

根据所述道路几何图形，生成泰森多边形，以得到基于泰森多边形原理的多个线段；

基于均质化网格，对所述道路几何图形进行切分，以得到多个网格；

针对所述多个网格中的每个网格，提取与所述道路几何图形无交点的线段并进行组合，以生成初版道路中心线；

对所述初版道路中心线进行简化处理，得到终版道路中心线；

所述处理单元还用于：

删除所述初版道路中心线中的断头线；

根据道格拉斯-普克算法，对删除断头线之后的初版道路中心线进行重组；

获取重组后的初版道路中心线中长度小于预设长度的线段；

针对每个长度小于预设长度的线段，将所述长度小于预设长度的线段的中点和与所述长度小于预设长度的线段连接的线段的外端点连接，所述外端点是不在所述长度小于预设长度的线段上的端点；

删除所述长度小于预设长度的线段以及与所述长度小于预设长度的线段连接的线段，以得到所述终版道路中心线。

6.根据权利要求5所述的道路中心线提取装置，其特征在于，所述处理单元还用于：

在所述道路几何图形的线框上，生成多个间距为预设距离的点；

根据所述多个间距为预设距离的点，生成泰森多边形，以得到基于泰森多边形原理的多个线段。

7.根据权利要求5所述的道路中心线提取装置，其特征在于，所述处理单元还用于：

删除所述道路几何图形中，面积小于预设面积或宽度小于预设宽度的图块；

根据删除图块后的道路几何图形，执行生成泰森多边形的步骤。

8.根据权利要求5所述的道路中心线提取装置，其特征在于，所述处理单元还用于：

针对所述多个网格中的每个网格，提取两个端点都在所述道路几何图形内的线段；

对所提取的线段，删除与所述道路几何图形有交点的线段，并进行组合，以生成初版道路中心线。

Images