路网自动划分交通小区的系统、方法、设备及存储介质与流程

1.本技术涉及一种自动划分交通小区的系统、方法、设备及存储介质，尤其涉及一种基于图像技术的路网自动划分交通小区的系统、方法、设备及存储介质，属于智能交通技术领域。


背景技术：

2.交通模型技术是交通行业重要的理论支撑依据，可以为城市的规划建设提供重要的数据支撑。交通小区划分是交通模型构建过程中的重要一环，交通小区划分的合理与否，直接关系到了建成交通模型的准确性，然后交通小区划分是一个比较繁琐的过程，需要耗费大量的人力工作，为了简化交通小区划分的繁琐性，有不少自动划分交通小区的方法被工程师提出。一部分技术基于城市路网和用地性质，另一部分技术基于海量城市交通数据，从统计分析的角度进行小区划分的确定，前者基本都是基于delaunay三角网进行的分区划分，小区边界和和真实道路相差较多，且划分过程中参数意义不明朗，难以被使用者理解，后者对输入数据的需求要求过高，使用成本较高。
3.现有最接近的方案有两种，一是基于voronoi图的交通道路网络自动划分交通小区；第二种是基于空间统计分析方法的城市路网交通小区自动划分。
4.（1）基于voronoi图的交通道路网络自动划分交通小区voronoi图，又叫泰森多边形或dirichlet图，它是由一组由连接两邻点直线的垂直平分线组成的连续多边形组成。该方法通过提取所有道路节点，根据点集创建delaunay三角网，然后建立对偶图voronoi多边形，形成初始的交通小区，然后考虑水系绿地，根据面积阈值进行合并、叠加分析，生成最终的交通小区划分结果。
5.该方法的主要局限性体现在以下方面：没有考虑城市道路等级，因为主要道路网络才能成为交通小区的分界线，次要道路网络如支路等应该成为小区内部道路。
6.（2）基于空间统计分析方法的城市路网交通小区自动划分目前该类方法主要基于空间聚类分析方法开展交通小区的自动划分，通过大量实时数据和历史数据挖掘交通特征较为相似的交通小区。
7.该方法的主要不足在于：需要大量的、完整的交通数据，成本高、效率低，划分出的交通小区边界和真实道路相差甚远，对于用户的输入要求过高，不适合作为一个通用工具，而适合于专题研究。
8.综合上述可知，现有交通小区划分主要存在的缺点为：（1）没有考虑城市道路等级；（2）需要大量的、完整的交通数据，成本高、效率低，划分出的交通小区边界和真实道路相差甚远；不适合作为一个通用工具，而适合于专题研究；（3）对于用户的输入要求过高；（4）难以应用在实际模型搭建过程中，非标准化工具，仅适合于专题研究。


技术实现要素：

9.有鉴于此，本技术提供一种路网自动划分交通小区的系统、方法、设备及存储介质，以减少交通小区划分过程中的人工投入成本、提高小区划分效率、简化用户输入为目的，只要求用户给出城市路网即可自动划分，解决了目前已有的自动划分交通小区的技术方案，划分出的交通小区边界和真实道路相差甚远，要求提供过多的输入数据，难以应用在实际模型搭建过程中，因此基于本技术技术方案所构建的是一个可标准化的工具，可以被大规模应用在交通小区划分工作中。
10.本技术的技术方案是这样实现的：方案一：一种基于图像技术的路网自动划分交通小区的系统，包括：路网缓冲区创建模块，用于将复杂道路的路网进行合并，将线形矢量图形转化为面域矢量图形，存储为二值图像格式；坐标信息记录模块，用于计算二值图像单位像素的经纬度跨度；缓冲区中轴线骨架提取模块，依据二值图像单位像素的经纬度跨度，通过迭代计算提取二值图像的主线轮廓，得到骨架图像；拓扑网络生成模块，基于骨架图像获取拓扑网络；查找最小基本环模块，采用贪心算法对拓扑网络进行最小基本环的查找；交通小区划分模块，基于最小基本环面层对面域矢量图形（街道图层）进行空间切割，得到初步的交通小区面层文件；交通小区合并模块，空间切割后，在街道图层的边界处形成小面积的交通小区，将小面积的交通小区合并到其邻近的交通小区中，得到最终的交通小区面层文件。
11.方案二：一种基于图像技术的路网自动划分交通小区的方法，包括：步骤一，路网缓冲区创建，将复杂道路的路网进行合并，将线形矢量图形转化为面域矢量图形，存储为二值图像格式；步骤二，坐标信息记录，计算二值图像单位像素的经纬度跨度；步骤三，缓冲区骨架提取，依据二值图像单位像素的经纬度跨度，通过迭代计算提取二值图像的主线轮廓，得到骨架图像；步骤四，拓扑网络生成，基于骨架图像建立一圈0分析值像素作为保护层，根据标识单位像素的分析值进行等值填充，标记结点在二值图像中的位置，追踪结点之间的连接像素，记录线层信息，将得到的拓扑网络的像素坐标还原为经纬度坐标；步骤五，查找最小基本环，采用贪心算法对拓扑网络进行最小基本环的查找；步骤六，交通小区划分，基于最小基本环面层对面域矢量图形（街道图层）进行空间切割，得到初步的交通小区面层文件；步骤七，交通小区合并，
空间切割后，在街道图层的边界处形成小面积的交通小区，将小面积的交通小区合并到其邻近的交通小区中，得到最终的交通小区面层文件。
12.进一步地，所述步骤一，路网缓冲区创建，具体步骤为：设复杂道路的路网有n条交通路线link，基于输入的路网，对每一条link对象创建缓冲区，将线形矢量图形转化为多边形的面域矢量图形，进行合并后，得到一个面域对象；记录面域矢量图形中所有几何点的最大经度max_lon、最大纬度max_lat、最小纬度min_lat、最小经度min_lon，然后将面域矢量图形存储为二值图像。
13.进一步地，所述步骤二，坐标信息记录，具体步骤为：得到的二值图像为一个矩阵，其元素的值为0或者1，找出矩阵中1值的最小行索引min_row_index、最大行索引max_row_index、最小列索引min_col_index、最大列索引max_col_index，然后对其进行裁切，这一步裁切的目的就是去除矩阵的周边零元素，以便正确记录单位像素跨越的经纬度，裁切后记目前矩阵的行像素个数为row_pixel_num，列像素个数为col_pixel_num，经度方向上单位像素所跨越的经度为lng_para，纬度方向上单位像素所跨越的纬度为lat_para，按照公式一、公式二计算，得到的矩阵作为骨架提取的输入数据；
ꢀꢀꢀꢀ
公式一
ꢀꢀ
公式二。
14.进一步地，所述步骤三，缓冲区骨架提取，具体步骤为：二值图像中的每一个元素代表一个像素，像素的值为0或者1，值为1的像素构成了图像，值为0的像素为空白；首先，定义二值图像的一个像素点p1的八个方向上的邻域，北向为p2，东北为p3，东向为p4，东南向为p5，南向为p6，西南向为p7，西向为p8，西北向为p9；其次，进行迭代计算，每次迭代分为两个子项，在第一次子迭代中，若像素点p1满足公式3及公式4的条件，则，该像素点p1被删除；
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式三
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式四
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式五
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式六公式三中表示非0邻域的数量，公式四中是有序集合中01模式的数量，如果不满足公式三至公式六中任何一个条件，那么，像素点p1不被删除；在第二个子迭代中，公式三、公式四不变，公式五、公式六替换为公式七和公式八，
如下：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式七
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式八公式七和公式八的解是或者或者且，被移除的像素点p1在图像的东边或者南边或者西北角点，在第二个子迭代中删除的像素点p1在西北或者东南角点；经过多轮迭代后，获得骨架图像。
15.进一步地，所述步骤四，拓扑网络生成，具体步骤为：1) 在骨架图像的周围加一圈0值像素作为图像的保护层；2) 若某像素值原来为0，经过迭代后仍然为0；3) 如果有某非0像素点p
i
，且b(p
i
)=2，则令该像素点的分析值为1，代表该像素点为骨架图像的中间段的像素；4) 如果有某非0像素点p
i
，且b(p
i
)≠2，则令该像素点的分析值为2，代表该像素点为骨架图像的端点和交点；5) 然后依次遍历分析值为2的像素点，进行等值填充，并且标记端点和交点（端点和交点统称为结点）在骨架图像中的位置；6) 最后追踪结点之间的连接像素，依次记录其线层信息，得到的拓扑网络的坐标是像素坐标，将其还原为经纬度坐标。
16.进一步地，所述步骤五，在生成的拓扑网络中查找最小基本环，具体步骤为：1) 找到拓扑网络中每个无向图的所有顶点之间的最短路径；2) 对于每一个顶点v和每一条边(x，y)，创建一个回路集：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式九p(v，x)表示v,x之间的最短路径，p(v，y)表示v，y之间的最短路径；3) 对回路集中的回路按照权重排序；4) 采用贪心算法找到回路集中的最小基本环。
17.进一步地，所述步骤七，交通小区合并，具体步骤为：获得最小基本环的顶点序列（顶点序列为若干个顶点组成的序列）后，结合经纬度坐标可转化为面层数据：进一步地，所述步骤七，交通小区合并，具体步骤为：设街道图层中有k个面域，记为；基于最小基本环面层对基础面域矢量图形（街道图层）进行空间切割后，中有z
i
个交通小区，初步划分的交通小区集合为，表示由第i个基础面域划分出的第j个初始交通小区，设置面积阈值为area_threshold；合并规则为：for（重复）：
找出初步划分的交通小区，筛选面积小于等于area_threshold的交通小区，设筛选出的交通小区的索引构成的集合为s，此时的交通小区构成集合；for
ꢀꢀ
(此循环嵌套在上层循环中)：在中找到与当前交通小区相邻的小区，若有多个，则匹配面积最大的，然后将小面积的交通小区合并到其邻近的交通小区中，得到最终的交通小区面层文件。此种设计机制保证了小区只会在同一个基础面层中进行合并，即不会生成跨基础面层（街道图层）的交通小区。
18.所述方案二是基于方案一所述系统实现的。
19.方案三：一种电子设备，包括处理器和用于存储能够在处理器上运行的计算机程序的存储器，其中，所述处理器用于运行所述计算机程序时，执行方案二所述方法的步骤。
20.方案四： 一种存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现方案二所述方法的步骤。
21.本技术有益效果体现在：本技术仍然基于路网进行交通小区划分，为了保留“主要道路作为交通小区边界”这一特征，通过提取道路骨架形成新的拓扑网络，然后基于城市基础图层进行小区划分，计算速度快，输入数据简单，泛化能力好，可快速帮助交通模型师划分城市交通小区，可大大简化划分工作的人力投入。具有的主要效果有：1) 提出一套基于路网自动划分交通小区的可行方法；2) 整合了路网骨架提取、骨架转拓扑、最小基本环查找、面域空间分析算法，提供了一套完整的系统及方法；3) 支持不同的道路使用不同的半径来创建缓冲区，进而可以有效提取多线道路的骨架；4) 所需基础数据少、参数含义明确、运行速度快，可标准化为一个模型工具。
附图说明
22.通过阅读参照以下附图所作的对非限制性实施例所作的详细描述，本技术的其它特征、目的和优点将会变得更明显：图1为本技术实施例一提供的路网自动划分交通小区的系统框图；图2为本技术实施例二提供的路网自动划分交通小区的方法流程图；图3为本技术实施例二方法细化流程图；图4为本技术实施例二中获取边界经纬度示意图；图5为本技术实施例二中像素点p1的八邻域示意图；图6为本技术实施例二中形成的像素序列示意图；图7为本技术实施例二中基于最小基本环面层对（街道图层）进行空间切割前示意图；
图8为本技术实施例二中基于最小基本环面层对（街道图层）进行空间切割后示意图；图9为本技术实施例二中交通小区合并后示意图；图10为本技术实施例二中以深圳路网作为切分路网、使用街道图层作为基础面层示意图；图11为本技术实施例二中深圳路网骨架网络拓扑图；图12为本技术实施例二中深圳路网最终划分得到的交通小区图层示意图；图13为本技术的一种电子设备的结构示意图；图14为本技术实施例二中交通小区合并流程示意图；图15为本技术实施例二中图像骨架提取示例图。
具体实施方式
23.下面结合附图和实施例对本技术作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释相关申请，而非对该申请的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与申请相关的部分。
24.需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本技术。
25.实施例相关符号说明如下：
link
i
路网中的第i条link多段线对象link_buffer_polygon
i
由第i条link对象生成的buffer面域merged_polygon所有link_buffer_polygon合并后的面域max_lon、min_lon最大经度、最小经度max_lat、min_lat最大纬度、最小纬度max_row_index、min_row_index最大行索引、最小行索引max_col_index、min_col_index最大列索引、最小列索引row_pixel_num行像素个数col_pixel_num列像素个数lng_para经度方向上单位像素所跨越的经度lat_para纬度方向上单位像素所跨越的纬度p
i
第i个像素点的值（或第i个像素点对象）b(p
i
)第i个像素点非0邻域数量a(p
i
)有序集合中01模式的数量base_region
i
第i个基础面域对象initial_taz
ij
由第i个基础面域图层划分出来的第j个初始交通小区small_taz
ij
由第i个基础面域图层划分出来的第j个面积小于阈值的初始交通小区area_threshold面积阈值
实施例一本技术实施例一提供了一种基于图像技术的路网自动划分交通小区的系统（参见图1），包括：路网缓冲区创建模块，用于将复杂道路的路网进行合并，将线形矢量图形转化为面域矢量图形，存储为二值图像格式；坐标信息记录模块，用于计算二值图像单位像素的经纬度跨度；
缓冲区中轴线骨架提取模块，依据二值图像单位像素的经纬度跨度，通过迭代计算提取二值图像的主线轮廓，得到骨架图像；拓扑网络生成模块，基于骨架图像获取拓扑网络；查找最小基本环模块，采用贪心算法对拓扑网络进行最小基本环的查找；交通小区划分模块，基于最小基本环面层对街道图层进行空间切割，得到初步的交通小区面层文件；交通小区合并模块，空间切割后，在街道图层的边界处形成小面积的交通小区，将小面积的交通小区合并到其邻近的交通小区中，得到最终的交通小区面层文件。
26.实施例二本技术实施例二提供了一种基于图像技术的路网自动划分交通小区的方法（参见图2
‑
图9），该方法具体为：s1，路网缓冲区创建，缓冲区的创建，对于小区的划分合理性有较为重要的影响，该步骤的主要目的有：1) 将立交、多线快速路等复杂道路进行合并，方便提取道路中心线；2) 将线形文件转化为面域文件，以存储为图片格式，作为提取骨架图像的输入文件。
27.具体步骤为：设复杂道路的路网有n条交通路线link，基于输入的路网，路网信息只要求包含几何信息即可，对每一条link对象创建缓冲区，将线形矢量图形转化为多边形的面域矢量图形，进行合并后，得到一个面域对象；记录面域矢量图形中所有几何点的最大经度max_lon、最大纬度max_lat、最小纬度min_lat、最小经度min_lon，然后将面域矢量图形存储为二值图像。边界经纬度参见图4。
28.s2，坐标信息记录，矢量文件被二值化后，只有像素坐标，为了在后续步骤中能够还原地理坐标，需要计算单位像素的经纬度跨度。
29.计算二值图像单位像素的经纬度跨度，作为骨架图像提取的输入数据；具体步骤为：得到的二值图像为一个矩阵，其元素的值为0或者1，找出矩阵中1值的最小行索引min_row_index、最大行索引max_row_index、最小列索引min_col_index、最大列索引max_col_index，然后对其进行裁切，这一步裁切的目的就是去除矩阵的周边零元素，以便正确记录单位像素跨越的经纬度，裁切后记目前矩阵的行像素个数为row_pixel_num，列像素个数为col_pixel_num，经度方向上单位像素所跨越的经度为lng_para，纬度方向上单位像素所跨越的纬度为lat_para，按照公式一、公式二计算，得到的矩阵作为骨架提取的输入数据；
ꢀꢀ
公式一
ꢀꢀꢀ
公式二
s3，缓冲区骨架提取，骨架提取是一个图形学分支下的形态学问题，直观的理解，就是提取目标图像的主线轮廓，本质上是提取目标在图像上的中心像素轮廓，细化后的目标都是单层像素宽度。
30.依据二值图像单位像素的经纬度跨度，通过迭代计算提取二值图像的主线轮廓，得到骨架图像；具体步骤为：一副二值图可由矩阵m来定义，矩阵中的每一个元素m
ij
代表一个像素，像素的值为0或者1，值为1的像素构成了图像，黑色部分的路网像素的值为1，白色部分的像素值为0。在图像的骨架化过程中，就是要删除那些属于骨架之外的所有的像素点，为了完成这个过程，需要进行迭代计算；首先，定义二值图像的一个像素点p1的八个方向上的邻域（参见图5），北向为p2，东北为p3，东向为p4，东南向为p5，南向为p6，西南向为p7，西向为p8，西北向为p9；其次，进行迭代计算，每次迭代分为两个子项，在第一次子迭代中，若像素点p1满足公式3及公式4的条件，则，该像素点p1被删除；
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式三
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式四
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式五
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式六公式三中表示非0邻域的数量，公式四中是有序集合中01模式的数量，如果不满足公式三至公式六中任何一个条件，那么，像素点p1不被删除；关于01模式数量的说明（参见图6），形成的像素序列值为：0（p2）1（p3）0（p4）0（p5）1（p6）0（p7）0（p8）0（p9），此时=2。
31.在第二个子迭代中，公式三、公式四不变，公式五、公式六替换为公式七和公式八，如下：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式七
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式八公式七和公式八的解是或者或者且，被移除的像素点p1在图像的东边或者南边或者西北角点，在第二个子迭代中删除的像素点p1在西北或者东南角点；经过多轮迭代后，获得骨架图像。
32.s4，拓扑网络生成，基于骨架图像建立一圈0分析值像素作为保护层，根据标识单位像素的分析值进
行等值填充，标记结点在二值图像中的位置，追踪结点之间的连接像素，记录线层信息，将得到的拓扑网络的像素坐标还原为经纬度坐标；拓扑网络生成，具体步骤为：1) 在骨架图像的周围加一圈0值像素作为保护层；2) 像素值原来为0的像素点，仍然为0；3) 如果某非0像素，其八个邻域有两个非0像素，那么标识该像素为1，代表了骨架图像中的中间段的像素；4) 如果某非0像素，其八邻域中的非0像素个数不是2（比如是1、3等），那么标识该像素为分析值2，代表了骨架图像中端点和交点部分的像素；5) 然后依次遍历分析值为2的像素，进行等值填充，并且标记结点（结点为端点和交点的统称）在骨架图像中的位置；6) 最后追踪结点之间的连接像素，记录其线层信息，得到的拓扑网络的坐标是像素坐标，将其还原为经纬度坐标。
33.s5，查找最小基本环，提取出道路骨架的拓扑文件后，需要对网络进行最小基本环的查找，查找算法如下：对于一个无向图g(v，e)(边权值为正数)，设图中边的数量为m，节点的数量为n；查找最小基本环，具体步骤为：1) 找到拓扑网络中每个无向图的所有顶点之间的最短路径；2) 对于每一个顶点v和每一条边(x，y)，创建一个回路集：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
公式九p(v，x)表示v,x之间的最短路径，p(v，x)表示v，x之间的最短路径，p(v，y)表示v，y之间的最短路径；3) 对回路集中的回路按照权重排序；4) 采用贪心算法找到回路集中的最小基本环。
34.s6，交通小区划分，基于最小基本环面层对街道图层进行空间切割，得到初步的交通小区面层文件；s7，交通小区合并（参见图14），由于部分道路和基础街道图层的边界距离较近，在使用路网分割街道图层后，在街道图层的边界处会形成较多的面积很小的小区，这部分小区应该被合并到其邻近的交通小区中。交通小区合并，具体步骤为：设街道图层中有k个面域，记为；基于最小基本环面层对街道图层进行空间切割后（切割前参见图7，切割后参见图8，合并优化后见图9，为了方便理解，这里仅仅使用局部图层进行划分说明），中有z
i
个交通小区，初步划分的交通小区集合为，表示属于第i个基础面域的第j个初始交通小区，设置面积阈值为area_threshold；
合并规则为：for（重复）：找出初步划分的交通小区，筛选面积小于等于area_threshold的交通小区，设筛选出的交通小区的索引构成的集合为s，此时的交通小区构成集合；for
ꢀꢀ
：在中找到与当前交通小区相邻的小区，若有多个，则匹配面积最大的，然后将小面积的交通小区合并到其邻近的交通小区中（图8中的initial_taz1‑2和initial_taz2‑2面积小于阈值，将和其父图层合并），得到最终的交通小区面层文件。
35.合并后交通小区面层文件如图9所示。
36.本示例使用深圳路网作为切分路网、使用街道图层作为基础面层（见图10），通过本方法生成得到的骨架网络拓扑见图11，最终划分得到的交通小区图层见图12。
37.实施例三本技术实施例三提供一种电子设备，参见图13，电子设备以通用计算设备的形式表现。电子设备的组件可以包括但不限于：一个或者多个处理器或者处理单元，用于存储能够在处理器上运行的计算机程序的存储器，连接不同系统组件(包括存储器、一个或者多个处理器或者处理单元)的总线。
38.其中，所述一个或者多个处理器或者处理单元用于运行所述计算机程序时，执行实施例二所述方法的步骤。所述处理器所用类型包括中央处理器、通用处理器、数字信号处理器、专用集成电路、现场可编程门阵列或者其他可编程逻辑器件、晶体管逻辑器件、硬件部件或者其任意组合。
39.其中，总线表示几类总线结构中的一种或多种，包括存储器总线或者存储器控制器，外围总线，图形加速端口，处理器或者使用多种总线结构中的任意总线结构的局域总线。举例来说，这些体系结构包括但不限于工业标准体系结构(isa)总线，微通道体系结构(mac)总线，增强型isa总线、视频电子标准协会(vesa)局域总线以及外围组件互连(pci)总线。
40.实施例四本技术实施例四提供一种存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现实施例二所述方法的步骤。
41.需要说明的是，本技术所示的存储介质可以是计算机可读信号介质或者存储介质或者是上述两者的任意组合。存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。存储介质的更具体的例子可以包括但不限于：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(ram)、只读存储器(rom)、可擦式可编程只读存储器(eprom或闪存)、光纤、便携式紧凑磁盘只读存储器(cd
‑
rom)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本技术中，存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者
器件使用或者与其结合使用。而在本技术中，存储介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。存储介质还可以是存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于：无线、电线、光缆、rf等等，或者上述的任意合适的组合。
42.着重强调本技术的技术关键点：1、缓冲区中轴线骨架提取方法；2、拓扑网络生成方法；3、整合了路网骨架提取、骨架转拓扑、最小基本环查找、面域空间分析算法，提供了一套完整的系统及方法和python源代码实现。
43.本技术中缩略语和关键术语定义如下：（1）交通模型交通模型是反映交通系统内在规律的数学模型组合，并通过数字、图形、影像、视频等形式来描述，融合交通工程学与社会学、人口学、经济学、统计学、行为学、信息学等多种学科理论，运用数理方法和计算机软硬件设备，对交通政策与规划、建设与投资、运行与管理等各阶段决策提供支持的重要定量分析技术。
44.（2）城市交通模型城市交通模型是一系列对城市交通进行解析，以揭示和反映交通系统内在规律的数学模型。对城市交通解析的透彻程度将决定模型质量的优劣。现状交通模型是对城市现状的解释，规划交通模型则是对城市未来发展的判断，也即对城市规划的解读。
45.（3）交通小区在城市交通研究中，为了将交通需求的产生、吸引与一定区域的社会经济指标联系起来，将交通需求在空间上的流动用交通小区之间的交通分布图表现出来，同时为了便于用交通分配理论模拟道路网上的交通流，在不打破行政区划、保持交通区的用地性质、交通特性应尽量一致的情况下，将城市划分为多个交通小区。
46.（4）二值图像二值图像是指在图像中，灰度等级只有两种，也就是说，图像中的任何像素点的灰度值均为0或者1，分别代表黑色和白色。
47.（5）图像骨架提取图像骨架提取可直观的理解为：通过一定的算法规则识别数字图像中非骨架多余像素点后删除，仅仅保留骨架像素点，示例见图15 。在复杂多线道路场景下，骨架提取对于道路简化非常有用。
48.（6）面域对象和多段线对象面域对象（polygon）和多段线对象（linestring）用于描述空间地理信息元素，一个面域对象在可视化软件中表现为一个二维平面区域，一个多段线对象在可视化软件中表现为一条多折点折线，多用来模拟真实路网，这两种对象均存储着坐标信息，因此可以进行合并、打断、切分等空间运算。
49.以上所述的实施例，对本技术的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本技术的实施例而已，并不用于限定本技术的保护范围，凡在本技术的技术方案的基础之上，所做的任何修改、等同替换、改进等，均应包括在本技术的保护范围之内。