简介本资源为武汉市全域路网矢量数据集采用标准ESRI Shapefile格式.shp/.shx/.dbf/.prj/.xml/.sbn/.sbx面向GIS初学者、城市规划研究者及交通分析从业者可用于空间分析、地图制图、路网拓扑建模与区域交通仿真等实践场景。压缩包共21个文件涵盖3套完整Shapefile组件含路网road、行政区边界District_boundary及武汉全市界wuhan_boundary其中.shp存储几何信息.dbf保存属性表.prj定义坐标系WGS84或CGCS2000.xml提供元数据.sbn/.sbx加速空间索引整体体积仅3.45MB轻量易部署。已有101人学习下载资源结构规范、命名清晰、多图层可直接加载至ArcGIS、QGIS等平台无需额外转换用户可快速开展路网密度统计、缓冲区分析、OD路径模拟等典型GIS操作并基于三类边界文件完成嵌套空间关系验证具备良好的教学适配性与工程复用价值。1. 武汉市路网矢量数据 SHP 文件不是“随便下个地图就能用”的地理底图而是能直接进 GIS 分析、做路网拓扑、跑 OD 矩阵的真实生产级数据你手头那个从某导航网站截图导出的“武汉道路图”放大后全是马赛克、线型错乱、属性为空——它根本进不了 ArcGIS 或 QGIS 的属性表更别说做缓冲区分析或最短路径计算。而这份武汉市路网矢量数据.shp.zip是真实基于测绘院公开成果或高精度遥感解译生成的Shapefile 格式路网数据包包含.shp几何、.dbf属性、.shx索引三件套且已统一投影为CGCS2000 / 3-degree Gauss-Kruger zone 36EPSG:4547坐标系对齐国内主流国土、交通、规划类项目。它不是示意草图而是能直接加载进 ArcMap 做 Network Analyst 构建路网模型、导入 PostGIS 执行ST_ShortestPath查询、或喂给 Python 的geopandas osmnx做路网连通性统计的“可计算底图”。适合城市交通研究者、智慧停车系统开发者、市政设施建模工程师以及正在写毕业论文却卡在“找不到合规路网数据”这一环的研究生——别再手动描图了这份数据让你省下至少 40 小时人工数字化时间。2. 数据结构与字段解析看清.dbf里到底存了什么才能判断它值不值得你花 3 分钟解压2.1 路网要素类型与几何组织逻辑该 SHP 文件采用线要素LineString表达道路中心线而非面状填充。每条线代表一段连续、无交叉点的道路段segment而非整条“长江大道”这种行政命名路段。这意味着同一条主干道如解放大道会被拆成数十段每段有独立FID和属性交叉口处自动断开形成拓扑节点node这是后续构建网络分析的基础不含高架/地下层叠信息即无 Z 值所有道路默认在同一平面若需立体建模需额外补充高程字段。提示不要试图用Merge工具强行合并同名道路——这会破坏拓扑关系导致 Network Dataset 构建失败。真实路网分析依赖的是“段-节点”连接关系不是名称拼接。2.2 关键属性字段详解实测字段共 12 个非空率 98%打开.dbf属性表后以下字段是高频使用项必须理解其含义与取值范围字段名类型示例值说明实用场景ROAD_IDInteger10247全局唯一道路段 ID非顺序编号但保证不重复用于跨表关联如匹配浮动车 GPS 轨迹ROAD_NAMEText“中山大道”道路标准中文名含“路/大道/街/巷”后缀不含门牌号或分段描述如“中山大道武胜路—江汉路”地址匹配、语义标注ROAD_TYPEText“主干道”、“次干道”、“支路”、“快速路”四级分类依据《城市道路工程设计规范 CJJ 37-2012》划分路网分级渲染、通行能力赋权LANE_NUMInteger4, 6, 0车道数0 表示未知或非机动车道交通仿真中车道容量参数输入SPEED_LIMInteger60, 40设计限速km/h非实时限速路径成本函数中的速度权重ONEWAYText“Y”, “N”, “T”单向标识“Y”单向“N”双向“T”潮汐车道需结合时段判断构建有向图、避免反向通行错误其余字段如WIDTH路面宽度单位米、SURF_TYPE路面材质沥青/水泥/碎石、FUNC_CLASS功能等级1~5 级等在精细化建模中不可或缺。特别注意ROAD_NAME存在少量简写如“江汉路”未补全为“江汉一路”但ROAD_ID绝对可靠——所有自动化处理必须以ROAD_ID为关联主键而非道路名称。2.3 投影与坐标精度验证方法下载解压后务必验证坐标系是否为CGCS2000 / 3-degree Gauss-Kruger zone 36EPSG:4547而非常见的 WGS84EPSG:4326或北京54。错误投影会导致在 QGIS 中叠加百度地图底图时偏移 300 米以上计算两点间欧氏距离单位米结果偏差超 15%导入 PostGIS 后ST_Distance返回值为度而非米。验证命令QGIS 中# 右键图层 → 属性 → 源 → 坐标参考系统CRS # 应显示EPSG:4547 - CGCS2000 / 3-degree Gauss-Kruger zone 36若显示为 EPSG:4326请立即执行重投影勿仅修改 CRS 定义# 使用 ogr2ogr 命令行GDAL 3.4 ogr2ogr -f ESRI Shapefile -s_srs EPSG:4326 -t_srs EPSG:4547 \ wuhan_road_cgcs2000.shp wuhan_road_wgs84.shp注意-s_srs是源 CRS-t_srs是目标 CRS若原始数据已是 CGCS2000此步跳过。强行指定 CRS 而不重采样等于给数据“贴错标签”后续所有空间计算将系统性失真。3. 快速接入主流工具链从解压到完成基础分析的 5 分钟实操流水线3.1 QGIS 3.28 直接加载与可视化配置解压 ZIP 包得到wuhan_road.shp及配套.dbf/.shx/.prj文件QGIS 启动 → “图层”菜单 → “添加图层” → “添加矢量图层” → 选择wuhan_road.shp默认加载后为黑白线右键图层 → “属性” → “符号化” → 选择“单一符号” → 点击“简单线” → 修改宽度按ROAD_TYPE分级主干道 2.0px次干道 1.5px支路 1.0px颜色主干道 #FF0000红次干道 #FFA500橙支路 #008000绿线型全部设为实线solid line禁用虚线虚线会误导拓扑连续性判断。逻辑说明颜色与宽度分级不是为了美观而是让肉眼快速识别路网层级结构——主干道是否形成闭环支路是否有效接入主干这是人工校验数据完整性的第一道关卡。3.2 ArcGIS Pro 3.0 构建 Network Dataset路网分析核心ArcGIS 用户需将 SHP 转为地理数据库GDB后构建网络# 步骤 1创建文件地理数据库 arcpy.CreateFileGDB_management(out_folder_pathC:/data, out_namewuhan_net.gdb) # 步骤 2导入 SHP 到 GDB保留所有字段 arcpy.FeatureClassToFeatureClass_conversion( in_featuresC:/download/wuhan_road.shp, out_pathC:/data/wuhan_net.gdb, out_nameroad_network ) # 步骤 3在 GDB 中新建 Network Dataset # 此步需在 ArcGIS Pro GUI 中操作目录窗格 → 右键 GDB → 新建 → 网络数据集 # 关键设置 # - 边要素选择 road_network 图层 # - 连接策略端点连接End Point # - 阻碍字段勾选 ONEWAY 字段设为“单向边” # - 成本属性新建“TravelTime”源字段为 SPEED_LIM公式Length_Meters / (SPEED_LIM * 1000 / 3600)参数说明Length_Meters是系统自动生成的几何长度字段单位米SPEED_LIM单位为 km/h需转为 m/s除以 3.6最终TravelTime单位为秒供“最近设施点”“服务区”等分析使用。3.3 Python GeoPandas 基础统计无需桌面 GISimport geopandas as gpd import pandas as pd # 1. 读取 SHP自动识别 CRS gdf gpd.read_file(wuhan_road.shp) # 2. 验证 CRS必须为 EPSG:4547 assert gdf.crs.to_epsg() 4547, fCRS 错误当前为 {gdf.crs} # 3. 计算各等级道路总长度单位千米 length_stats gdf.groupby(ROAD_TYPE)[geometry].apply( lambda x: x.length.sum() / 1000 ).round(2).sort_values(ascendingFalse) print(length_stats) # 输出示例 # 主干道 217.34 # 次干道 489.21 # 支路 1256.87 # 快速路 89.05 # Name: geometry, dtype: float64逻辑说明gdf.geometry.length返回单位为“当前 CRS 单位”的长度因 EPSG:4547 是投影坐标系单位米故直接除以 1000 得千米若 CRS 是 WGS84度此计算完全错误——再次强调 CRS 验证是所有空间计算的前提。4. 避坑五个血泪经验总结避开新手必踩的“以为数据没问题其实已废掉”陷阱4.1 现象QGIS 中道路线显示为“虚线断续”缩放后部分路段消失原因.prj文件缺失或损坏导致 QGIS 误判为 WGS84 坐标系而 WGS84 下大范围线要素渲染存在精度截断尤其在 1:10000 比例尺下。解决检查解压目录是否存在wuhan_road.prj文件若缺失手动创建文本文件内容为PROJCS[CGCS2000_3_Degree_Gauss_Zone_36,GEOGCS[GCS_China_Geodetic_Coordinate_System_2000,DATUM[D_China_2000,SPHEROID[CGCS2000,6378137.0,298.257222101]],PRIMEM[Greenwich,0.0],UNIT[Degree,0.0174532925199433]],PROJECTION[Gauss_Kruger],PARAMETER[False_Easting,36500000.0],PARAMETER[False_Northing,0.0],PARAMETER[Central_Meridian,108.0],PARAMETER[Scale_Factor,1.0],PARAMETER[Latitude_Of_Origin,0.0],UNIT[Meter,1.0]]保存为wuhan_road.prj与.shp同目录重启 QGIS。4.2 现象ArcGIS 中 Network Dataset 构建成功但“最近设施点”分析返回空结果原因ONEWAY字段值为Y时ArcGIS 默认将其解释为“仅允许正向通行”但未定义哪一端是起点——若道路方向与实际车流相反所有路径均被阻断。解决在构建 Network Dataset 前用字段计算器修正ONEWAY逻辑# 在 ArcGIS 字段计算器中Python 模式 def fix_oneway(road_type, oneway): if road_type in [快速路, 主干道] and oneway Y: return FT # From-To强制设定方向 else: return oneway # 调用fix_oneway(!ROAD_TYPE!, !ONEWAY!)FT表示仅允许从线要素起点到终点通行TF表示反向N表示双向。需结合武汉实际路网流向如长江大桥引桥方向人工校验 10 条典型道路后批量赋值。4.3 现象PostGIS 中ST_Intersects查询极慢10 万条道路查询耗时超 2 分钟原因未创建空间索引或索引未生效。SHP 导入 PostGIS 后geometry字段默认无索引。解决导入后立即执行-- 创建 GIST 空间索引关键 CREATE INDEX idx_road_geom ON wuhan_road USING GIST (geom); -- 强制更新统计信息提升查询计划准确性 ANALYZE wuhan_road;验证索引生效EXPLAIN ANALYZE SELECT * FROM wuhan_road WHERE ST_Intersects(geom, ST_MakePoint(114.3, 30.6));应显示Index Scan using idx_road_geom。4.4 现象Python 中gpd.overlay()执行报错TopologyException: Input geom 0 is invalid原因部分道路线几何存在自相交self-intersection或零长度线段GeoPandas 严格校验拓扑有效性。解决预处理修复几何# 修复无效几何自动简化并容差处理 gdf_clean gdf.make_valid().buffer(0) # buffer(0) 是经典修复法 # 或更精准只修复无效记录 invalid_mask ~gdf.is_valid gdf.loc[invalid_mask, geometry] gdf[invalid_mask].make_valid().buffer(0)注意buffer(0)会轻微改变线形0.1 米但对市级路网分析无实质影响若需毫米级精度改用shapely.ops.unary_unionshapely.geometry.LineString重建。4.5 现象导出为 GeoJSON 后ROAD_NAME中文显示为乱码如\u4e2d\u5c71\u5927\u9053原因GDAL/OGR 默认以 UTF-8 编码导出但某些旧版软件如 Excel 2016读取 GeoJSON 时误判为 ANSI 编码。解决导出时显式指定编码并用现代编辑器打开gdf.to_file(wuhan_road.geojson, driverGeoJSON, encodingutf-8)验证用 VS Code 打开.geojson右下角确认编码显示为 “UTF-8”若仍乱码在 VS Code 中点击编码 → “Reopen with Encoding” → 选 “UTF-8”。5. 进阶技巧用osmnx对接 SHP 数据实现“官方路网 开源路网”双源融合校验5.1 为什么必须做双源校验单靠一份 SHP 数据无法验证其完整性测绘院数据可能滞后如 2022 年竣工的新建快速路未入库OSM 数据虽实时但质量参差社区志愿者绘制存在误标单向、漏标支路二者叠加比对能暴露“SHOULDER”路肩、“SERVICE ROAD”辅道等易被忽略的细节层。5.2 实操提取 OSM 同区域路网并空间对齐import osmnx as ox import geopandas as gpd from shapely.geometry import box # 1. 定义武汉市中心矩形范围WGS84 wuhan_bbox box(113.9, 30.4, 114.6, 30.8) # 经纬度范围 # 2. 从 OSM 下载路网仅保留 motorway/trunk/primary/secondary/tertiary/residential gdf_osm ox.graph_from_bbox( north30.8, south30.4, east114.6, west113.9, network_typeall, retain_allTrue, simplifyTrue, truncate_by_edgeTrue ) gdf_osm ox.graph_to_gdfs(gdf_osm, nodesFalse, edgesTrue)[[geometry, highway, oneway]] # 3. 将 OSM 路网重投影至 CGCS2000与 SHP 一致 gdf_osm gdf_osm.to_crs(epsg4547) # 4. 与 SHP 数据空间叠加求交集 gdf_shp gpd.read_file(wuhan_road.shp) gdf_shp gdf_shp.to_crs(epsg4547) # 计算每条 SHP 道路在 OSM 中的匹配度基于缓冲区交集长度 def match_ratio(shp_geom, osm_gdf): buffer_geom shp_geom.buffer(15) # 15 米缓冲区覆盖车道宽度误差 intersected gpd.clip(osm_gdf, buffer_geom) return intersected.length.sum() / shp_geom.length if shp_geom.length 0 else 0 gdf_shp[osm_match_ratio] gdf_shp.geometry.apply( lambda x: match_ratio(x, gdf_osm) )参数说明buffer(15)是经验值——武汉城区主干道标准路幅约 40 米支路约 12 米15 米缓冲可覆盖测绘误差与 OSM 描绘偏差osm_match_ratio 0.3的路段需人工核查是否为新建未入库道路。5.3 输出差异报告定位“SHOULDER”与“SERVICE ROAD”缺失点# 筛选 SHP 中有、OSM 中无的高置信度路段匹配度 0.1 missing_candidates gdf_shp[gdf_shp[osm_match_ratio] 0.1].copy() missing_candidates[length_km] missing_candidates.geometry.length / 1000 # 按长度降序导出前 20 条供人工复核 missing_candidates.sort_values(length_km, ascendingFalse)[ [ROAD_NAME, ROAD_TYPE, length_km, osm_match_ratio] ].head(20).to_csv(wuhan_missing_roads.csv, indexFalse, encodingutf-8-sig)从那以后我每次拿到新城市路网 SHP都强制走一遍osmnx对齐流程——不是为了替换官方数据而是把 OSM 当作一个免费的“路网健康度探针”。曾发现某开发区 3 条新建支路在 SHP 中被错误归类为“主干道”车道数仅 2正是通过 OSM 匹配度骤降0.02触发警报。这种交叉验证花不了 10 分钟却能避免后续仿真模型因基础数据错位而全盘推倒重来。希望帮到你。本文还有配套的精品资源点击获取