简介武汉市路网矢量数据是一套面向GIS用户的陆路交通与行政边界数据包适用于城市规划、交通空间分析、专题制图等场景。这个zip压缩包共包含21个文件核心内容为三组shp矢量图层分别对应武汉市道路路网、区级行政边界与全市行政边界同时配有dbf属性表、prj投影坐标、shx空间索引及xml元数据便于在ArcGIS等桌面软件中直接加载与二次编辑。数据只保留公路等陆路要素不含地铁、铁路、水路或航空线路适合需要纯净路网底图的入门及进阶使用者且整个压缩包仅3.45MB目前已有2928人浏览学习。解压后可以快速获取道路线图层与两级行政区划边界用来做范围裁剪、缓冲区分析、路网密度统计或基础底图叠加输出足以支撑武汉市相关的空间分析与可视化课题。1. 武汉市路网矢量数据shp一份能直接进ArcGIS的路网基底做城市分析、交通可达性测算或者片区规划底图的时候最头疼的往往不是算法而是手里没有一套干净、连续、属性完整的路网矢量数据。武汉市路网矢量数据shp这套文件解决的就是这个基础问题。它不像网上那些零散的OSM导出块而是把武汉市域范围的道路按等级分层整理好道路中心线、路名、等级、长度这些关键字段基本齐全拿到手可以直接挂到ArcGIS里做显示、查询、缓冲区分析或者网络分析。这套数据适合谁用做城市规划的、搞交通研究的、做GIS课程设计的以及任何需要武汉市路网做底图或分析建模的开发者和学生。它不包含建筑物轮廓、POI或者地块边界专注把路网结构做扎实。2. 数据结构和坐标系先看懂shp里到底装了什么2.1 shp文件家族的构成和打开方式Shp格式和普通Excel不一样一个完整的路网数据不是单个文件而是由至少三个基础文件组成。.shp存储几何坐标.shp的兄弟文件.dbf存储属性表.shx做几何索引。很多新手拿到文件后只看到一个.shp就双击打开结果ArcGIS提示缺少文件或者显示空白就是因为漏掉了.dbf和.shx。在ArcMap或ArcGIS Pro里添加数据时应该直接选择.shp文件但要把三个文件放在同一个文件夹里路径最好不要带中文和空格否则容易出现读取异常。武汉市路网矢量数据shp的几何类型是Polyline也就是道路中心线。每个要素代表一条独立的路段路段之间通过节点衔接形成完整的路网拓扑。用ArcGIS打开后打开属性表能看到字段一般包括道路等级、道路名称、长度、方向等字段。2.2 坐标系CGCS2000还是Web墨卡托直接影响分析结果武汉路网shp的坐标系是一个关键参数。通常来说国内的路网shp分两类一类是WGS84经纬度坐标一类是CGCS2000投影坐标或Web墨卡托投影坐标。拿到数据的第一步在ArcGIS里通过右键图层属性查看源选项卡确认坐标系类型。如果数据是经纬度坐标直接做缓冲区分析或者长度计算结果会是度而不是米数值看起来很小用起来全是坑。我一般拿到这套武汉路网shp后会先做一次投影转换。武汉地处东经114度附近适合使用CGCS2000 3度分带的中央经线114度投影带这样路网长度计算出来的单位是米做缓冲区、网络分析都能对上物理距离。操作上在ArcToolbox里选择数据管理工具投影和变换要素投影选择一个本地投影坐标系输出即可。# 投影转换示例使用ArcGIS Pro的Python窗口或独立环境 import arcpy input_fc rD:\wuhan_road\road_network.shp output_fc rD:\wuhan_road\road_network_prj.shp # 定义CGCS2000 3度分带中央经线114度 out_coordinate_system arcpy.SpatialReference(CGCS2000 3-degree Gauss-Kruger zone 37) arcpy.management.Project(input_fc, output_fc, out_coordinate_system) print(投影转换完成)逻辑说明先读取原始坐标系再指定目标坐标系。CGCS2000 3度分带中央经线114度对应的zone参数是37如果不确定带号直接用arcpy.Describe()读取原始数据的extent来判断经纬度范围也能反推带号。参数说明out_coordinate_system可以换成其他投影比如Web墨卡托WGS 1984 Web Mercator Auxiliary Sphere但是如果做面积和长度分析Web墨卡托的形变会比较大不建议在精度要求高的场景下使用。2.3 道路分层和属性字段怎么区分快速路、主干道、次干道和支路武汉市路网矢量数据shp在制作时一般会按照城市道路标准进行分级。快速路、主干路、次干路、支路在字段里用数字编码区分通常1表示快速路2表示主干路3表示次干路4表示支路。没有这个字段的数据道路就是一堆没有语义的线段只能看形状不能做任何等级相关的分析。打开属性表后用符号系统按等级字段做唯一值渲染快速路用红色主干道用橙色次干路用黄色支路用灰色一张标准的城市路网图就出来了。如果发现等级字段缺失可以基于道路宽度或车道数人工补录但这属于比较大的工作量建议优先筛选已有字段里是否有名称、等级相关的字段。属性表中还有一个实用的字段是道路名称比如“中山大道”“解放大道”“珞喻路”等这个字段可以直接用于地址匹配和路径规划的显示。需要注意的是有些道路名称在不同区域有重复分析时要注意加上行政区代码做联合主键。3. 路网数据预处理与可用性校验从加载到能用的完整流程3.1 字段清理与坐标系探查的标准化操作拿到武汉市路网矢量数据shp文件的第一天我建议按一套固定流程走一遍第一步加载数据查看几何形状和坐标系第二步打开属性表查看字段类型和空值情况第三步做一次几何修复和投影转换第四步做一次拓扑检查确认道路之间没有断裂和悬挂节点。import arcpy # 设置工作空间 arcpy.env.workspace rD:\wuhan_road arcpy.env.overwriteOutput True # 检查几何是否存在错误 repair_result arcpy.management.RepairGeometry(road_network.shp) if repair_result.status Succeeded: print(几何修复完成) # 统计字段的空值数量 field_names [ROAD_NAME, ROAD_CLASS, LENGTH] for field in field_names: null_count arcpy.management.GetCount( arcpy.analysis.Select(road_network.shp, in_memory/selected, f{field} IS NULL) ) print(f{field} 空值数量: {null_count})逻辑说明RepairGeometry解决要素几何对象损坏的问题比如自相交、环方向异常、坐标非法等。空值统计则是判断路名和等级字段是否完整如果一个路网数据有30%以上的道路名称缺失那么后续路径规划时无法给出有意义的方向指引这种情况下要做字段补充或者放弃名称依赖。参数说明{field} IS NULL语法是ArcGIS的标准SQL查询表达式注意字段名要用包围。3.2 在ArcGIS中做路网显示优化和符号化设置武汉市路网数据在ArcGIS里加载后用符号系统做分层渲染。按道路等级字段设置线宽和颜色快速路5个点宽、主干路3个点、次干路2个、支路1个。做这一步的目的不只是好看是为了在制图输出和检查道路连通性时能直观分辨断头路和等级突变。显示优化之后做一次标注标注字段选道路名称设置合适的显示级别让缩放到1:10000时主要道路名称显示缩放到1:5000时支路名称显示。# 使用arcpy.mp设置图层符号系统 import arcpy aprx arcpy.mp.ArcGISProject(rD:\wuhan_road\wuhan.aprx) map aprx.listMaps(Map)[0] layer map.listLayers(road_network_prj)[0] # 统一符号化之后导出制图 aprx.save() print(制图工程保存完成)逻辑说明这一步核心是把数据的视觉语义调整到适合审阅的状态符号化本身不影响数据内容但影响后续人工检查效率。符号设置的具体数值以图层文件的saveLayerFile为准在ArcGIS Pro里修改完符号后保存为.lyrx文件便于后续复用。3.3 按行政区切片导出武汉各区路网武汉市地域范围大在做片区分析时通常只需要某一个区的路网。按行政区边界做裁剪导出一份分区的路网数据能减少后续分析的计算量。用武汉市某行政区的面要素作为裁剪范围逐区导出。import arcpy borough_fc rD:\wuhan_road\boroughs.shp road_fc rD:\wuhan_road\road_network_prj.shp output_dir rD:\wuhan_road\by_borough # 遍历每个行政区要素 with arcpy.da.SearchCursor(borough_fc, [NAME]) as cursor: for row in cursor: borough_name row[0] where_clause fNAME {borough_name} borough_temp in_memory/borough_filter road_clipped arcpy.analysis.Select(borough_fc, borough_temp, where_clause) out_name f{output_dir}\\{borough_name}_road.shp arcpy.analysis.Clip(road_fc, road_clipped, out_name) print(f{borough_name} 路网已导出)逻辑说明遍历行政区名称用Select提取当前行政区边界再以这个边界做Clip裁剪路网。in_memory是内存工作区效率比磁盘读写快。参数说明行政区名称字段要和属性表里的实际字段名一致比如“NAME”或者“区名全称”如果字段名是中文或者带空格需要用字段名转义写法否则查询会报语法错误。3.4 路网长度字段的计算与单位校正路网分析经常要统计每个行政区的道路总长度、路网密度。坐标系未投影的路网数据计算结果没有物理意义所以长度计算必须在投影坐标系下执行。import arcpy # 在投影后的路网中添加长度字段 road_fc rD:\wuhan_road\road_network_prj.shp arcpy.management.AddField(road_fc, LEN_M, DOUBLE, field_alias长度(米)) # 几何计算长度 arcpy.management.CalculateGeometryAttributes( road_fc, [[LEN_M, LENGTH]], coordinate_systemarcpy.SpatialReference(CGCS2000 3-degree Gauss-Kruger zone 37) ) # 分行政区统计长度 stat_result arcpy.analysis.Statistics( road_fc, rD:\wuhan_road\road_stats.dbf, [[LEN_M, SUM]], BOROUGH ) print(路网长度统计完成)逻辑说明CalculateGeometryAttributes在计算长度时指定坐标系是避免单位错误的关键。之后再按行政区字段做Statistics汇总得到每个区域的路网总长度。参数说明BOROUGH字段是行政区名称字段如果没有这个字段就做空间连接先把行政区名称附加到每条道路要素上再做统计。4. 避坑与常见问题排查武汉路网数据处理中的五个实战坑4.1 属性表中文乱码现象在ArcMap里打开武汉路网shp发现中文字段名和道路名称显示成乱码类似“µÀ·”“???”这样的符号。原因shp的.dbf文件属性编码格式是GBK或者ANSI而ArcGIS默认读取的是UTF-8系统区域设置不同导致乱码。解决打开ArcMap前先设置环境变量或者直接在ArcGIS Pro里对dbf文件执行Toolbox → Conversion Tools → Excel相关转码操作。最稳妥的做法是使用Python按GBK读取字段内容后重新生成一个UTF-8编码的dbf文件还有一种就是在ArcCatalog里选中dbf文件右键属性里修改代码页为GBK或936再加载到地图。4.2 道路显示偏移或坐标飞点现象将武汉市路网shp加载到底图或与其他数据叠加时路网与影像底图错位有的甚至出现在地图窗口之外看不到。原因坐标系混乱。常见情况是数据的坐标系是CGCS2000投影坐标但.prj文件缺失或者写错成WGS84导致ArcGIS按错误方式动态投影。解决先看.prj文件是否存在缺失就按数据采样点的经纬度范围反推坐标系使用Define Projection工具手动指定如果数据本身偏移距离固定用Shift工具整体平移校正。4.3 道路在路口处断裂现象放大路网后发现交叉路口处有细小的空隙道路线没有真正相连路径规划时路线无法转弯或者连通性很差。原因原始路网数据采集自不同来源拼接时节点没有做捕捉(Snap)线段端点悬浮在交叉点附近但不重合。解决用拓扑工具集的Integrate工具设置合适的聚类容差把距离小于容差范围内的节点合并。注意容差不能设置太大否则会把平行道路粘到一起。import arcpy road_fc rD:\wuhan_road\road_network_prj.shp # 设置容差0.5米把端点捕获到一起 arcpy.management.Integrate([road_fc], 0.5 meters) print(节点捕捉完成)逻辑说明Integrate直接修改原始要素的几何执行前要备份数据。容差0.5米适用于城市路网如果是郊区快速路间距大容差可以适当调整。参数说明0.5 meters是带单位的字符串不写单位时默认使用数据框单位容易出错建议写完整。4.4 道路等级字段值全是0或者空值现象按等级字段做符号化时所有道路都是同一个颜色检查属性表发现等级字段全是0或者空着。原因数据生产时等级赋值过程遗漏或者字段映射失败。解决基于道路名称或者道路宽度做推断比如宽度在50米以上的通常为主干路20米以下的通常是支路。也可以用空间关系匹配把已知等级的道路做缓冲区覆盖到的未知道路自动继承等级。4.5 导出SHP后大小异常或要素丢失现象从ArcGIS导出武汉路网为新的shp文件后文件数量比预想中少打开后发现部分要素缺失。原因导出时设置了范围约束或者数据中存在空几何要素被过滤掉。还有一种情况是字段名过长导出时被截断导致覆盖。解决导出时确认环境设置里输出坐标系和范围设置明确为“数据的范围”而不是“当前视图范围”。用Repair Geometry和Check Geometry两个工具先跑一遍确保不存在Null几何再执行导出。5. 基于武汉路网做网络分析从拓扑检查到OD成本矩阵5.1 构建网络数据集的准备工作路网数据要用于路径规划需要先构建网络数据集。武汉市路网shp构建网络数据集之前需要确保每个交叉口节点是连续的。经过前一步的Integrate捕捉处理后路网拓扑已经基本干净接下来需要为网络分析设置连通性策略。在ArcCatalog里右键路网要素集选择新建网络数据集连通性策略选择端点连通也就是说道路只在端点处交汇。武汉的道路绝大多数是平面交叉这个策略最符合实际。如果某个立交桥需要跨线不连通可以通过添加高程字段模拟立交关系但默认不设置。5.2 单点路径规划和时间成本设置import arcpy network_dataset rD:\wuhan_road\road_network_nd.nd arcpy.env.workspace rD:\wuhan_road arcpy.env.overwriteOutput True # 创建路径分析图层 arcpy.na.MakeRouteAnalysisLayer( network_dataset, route_analysis, Driving Time, 0.1 Meters, ALLOW_UTURNS, ONE_WAY_AGAINST_FLOW, NO_HIERARCHY, NO_LINES, NO_SIMPLIFY ) # 添加起点和终点位置 stops_fc rD:\wuhan_road\route_stops.shp arcpy.na.AddLocations(route_analysis, Stops, stops_fc, , 5000 METERS) # 求解并输出路径 arcpy.na.Solve(route_analysis, SKIP) print(路径求解完成)逻辑说明MakeRouteAnalysisLayer指定了出行模型为“Driving Time”这样道路等级、长度、方向限制都会参与阻抗计算。武汉道路中单向交通的比例不低ONE_WAY_AGAINST_FLOW参数确保路线不会违反单行线规则。参数说明5000 METERS是搜索距离当起点或终点不在路网节点上时ArcGIS在5000米范围内寻找最近的可通行节点。如果数据经过仔细预处理路网点位准确这个值可以缩小到1000米减少匹配到错误路段的风险。5.3 OD成本矩阵批量计算与路网连通性检验做武汉市多个行政区之间的通勤成本分析用OD成本矩阵批量计算比逐条做路径规划效率高得多。起点用各区中心点终点同样输出的是两两之间的行驶时间矩阵。import arcpy network_dataset rD:\wuhan_road\road_network_nd.nd arcpy.na.MakeODCostMatrixLayer( network_dataset, od_matrix, Driving Time, 2, 5000 METERS, ONE_WAY_AGAINST_FLOW, NO_HIERARCHY ) origins rD:\wuhan_road\origins.shp destinations rD:\wuhan_road\destinations.shp arcpy.na.AddLocations(od_matrix, Origins, origins, , 5000 METERS) arcpy.na.AddLocations(od_matrix, Destinations, destinations, , 5000 METERS) arcpy.na.Solve(od_matrix, SKIP) # 导出OD结果表 od_table rD:\wuhan_road\od_results.dbf arcpy.management.CopyRows(od_matrix, od_table) print(OD成本矩阵已生成)逻辑说明MakeODCostMatrixLayer的2表示同时计算起点到终点的去程和返程成本对于没有走行方向属性的路网来说去返成本通常一致但武汉的单行线会导致差异保留这个设置能检查方向不对称路段。参数说明CopyRows输出的是标准dbf表包含OriginID、DestinationID、Total_DrivingTime字段。这张表可以直接用于做行政区可达性排序和可视化。5.4 路网连通性验证的Quick检查方式网络数据集构建完成之后可以用一个简单方法验证路网的连通性。随机抽取地图上的几个远离干道的点位做任意两点间路径规划看是否能计算出路线。如果某个区域路网断裂求解会返回空结果或路线为0。批量做这个方法可以用循环遍历所有行政区中心点两两匹配。从那以后我在每套武汉市路网数据交付给分析流程之前都会强制走一遍这个验证方式构建网络数据集、跑一次OD矩阵、抽查若干个跨区域路线结果确认没有空值。这个过程能够把大量隐蔽的数据问题在正式分析之前暴露出来希望帮到你——拿到这套武汉市路网矢量数据shp之后先按第3章的流程做投影和拓扑修复再做网络分析能省掉一大批回头改数据的麻烦。本文还有配套的精品资源点击获取
