简介面向城市规划、GIS 分析与交通研究场景的武汉市路网矢量数据包适合需要本地化基础地理底图的学生、研究人员与制图从业者使用。数据以 shp 格式组织涵盖武汉市道路矢量图层以及各行政区边界和武汉市整体边界道路部分仅保留陆路数据不含地铁、铁路、水路与航空要素便于直接叠加分析或作为专题图底图。压缩包共 21 个文件由 shp、dbf、shx、prj、sbn、sbx、xml 等组成分别承载几何图形、属性表、空间索引、坐标系定义与元数据说明整体约 3.45MB体量轻便、结构规整导入 ArcGIS 等软件即可读取。已有 2928 人学习下载说明其在同类数据中具备一定参考价值。读者可获得一套可直接用于行政区划裁剪、路网密度统计与地图配准的矢量底图省去自行采集与整理边界、道路图层的时间成本。1. 拿到一份城市路网 shp先别急着往 ArcGIS 里拖上周有个做交通仿真的小哥找我说他手上有一份武汉市路网矢量数据 shp拖进 ArcGIS 之后地图一片空白属性表倒是满的坐标死活对不上。这不是数据坏了是坐标系没认出来。城市路网 shp 这类资源本质是一套带属性表的线要素集合字段里通常有道路等级、名称、长度这类信息能不能用起来八成取决于你对坐标系和字段结构的判断而不是软件操作熟不熟。这份数据适合做交通可达性分析、路径规划底图、城市肌理研究也适合当成 GIS 入门练手的第一份真实城市数据。下面按我实际拆包的顺序讲从结构、坐标系、字段清洗到空间分析把能抄的步骤和会翻车的地方都摆出来。2. 拆开 shp 包四个必需文件和字段结构先摸清2.1 shp 不是单个文件少一个都打不开很多人以为 shp 就是一个文件解压完只把.shp拖进软件结果报错或者属性丢失。shapefile 是一组文件的统称缺一不可的四个是扩展名作用缺失后果.shp存几何坐标直接打不开.shx几何索引能打开但无法空间查询.dbf属性表能显示图形但没有字段.prj坐标系定义坐标数值在但位置错乱除此之外常见还有.cpg字符编码声明和.sbn/.sbx空间索引。武汉市路网这种中文属性数据.cpg特别关键它决定了 dbf 里的中文字段是按 UTF-8 还是 GBK 读缺了它 ArcGIS 里道路名称大概率变乱码。我一般拿到包先做一件事确认这四个文件同名同目录再检查有没有.cpg。# 列出目录下所有同名文件确认 shp 家族完整 ls -l wuhan_road.* # 期望看到 wuhan_road.shp / .shx / .dbf / .prj 至少四个这段命令就是核对文件家族是否齐全。如果只有.shp别急着用先找原始提供方要齐或者用ogr2ogr从其他格式转出来补齐硬用残缺文件后面做空间连接一定出问题。2.2 用 ogrinfo 先读元数据比直接开软件快在往 ArcGIS 里拖之前我习惯先用 GDAL 的命令行工具把底细摸一遍不依赖图形界面信息还全。# 查看数据基本信息要素数、几何类型、坐标系、字段 ogrinfo -so -al wuhan_road.shp-so是 summary only只输出摘要不列每个要素-al是对所有图层操作。输出里重点看三行Geometry: Line String确认是线要素Feature Count看路网规模Layer SRS WKT看坐标系是不是投影坐标系。如果 SRS 显示(unknown)说明.prj缺失或损坏这就是第 1 章那个小哥坐标对不上的根因。2.3 字段结构决定这份数据能干什么路网 shp 的价值一半在几何一半在属性表。用下面命令把字段名和类型列出来# 只输出字段定义不读几何速度快 ogrinfo -so wuhan_road.shp wuhan_road -sql SELECT * FROM wuhan_road LIMIT 0输出里会看到类似NAME: String、CLASS: Integer、LENGTH: Real这样的字段。城市路网数据常见的字段语义是道路名称、道路等级快速路/主干路/次干路/支路、长度、方向。等级字段是后面做分层渲染和可达性加权的关键如果只有名称没有等级就得靠名称规则或人工补这一点在选数据时要心里有数。字段名如果是拼音缩写如DLMC、DLDJ别猜对照数据说明或者抽样看值域再定。3. 坐标系这关武汉路网为什么总“飘”到海里3.1 地理坐标系和投影坐标系差的是单位不是数字武汉位于东经 113°41′115°05′北纬 29°58′31°22′。如果一份路网 shp 的坐标值是114.3, 30.6这种量级那是地理坐标系WGS84 或 CGCS2000单位是度如果是500000, 3370000这种量级那是投影坐标系单位是米。两者混用图形要么缩成一个点要么飞到几千公里外。ArcGIS 里“飘到海里”的经典现象就是数据本身是投影坐标但被当成地理坐标显示或者.prj写的是 WGS84 而实际是别的基准。判断方法很简单看ogrinfo输出里的坐标范围或者直接在 ArcGIS 里看右下角坐标读数。度级数值带小数点后好几位米级数值是六到七位整数。3.2 用 ogr2ogr 重投影一步到位确认原始坐标系后如果要做长度计算、缓冲区、路网密度这类以米为单位的分析必须转到投影坐标系。国内常用 CGCS2000 高斯-克吕格投影武汉落在 3 度带第 38 带中央经线 114°E。# 从地理坐标系转到 CGCS2000 3度带38带投影 ogr2ogr -f ESRI Shapefile \ -t_srs EPSG:4547 \ wuhan_road_proj.shp wuhan_road.shp-t_srs是目标坐标系EPSG:4547 对应 CGCS2000 / 3-degree Gauss-Kruger CM 114E。如果你的数据本来就是投影坐标但带号不对用-s_srs显式指定源坐标系再转避免 GDAL 猜错。转完再用ogrinfo -so看一眼坐标范围数值应该落在几十万到几百万米之间这才对。提示重投影只改坐标表达不改几何形状。如果原始数据基准就是错的比如把西安80当CGCS2000重投影救不回来只能找控制点做配准这是另一套流程。3.3 ArcGIS 里的对应操作和常见误判在 ArcGIS Pro 里用“投影”工具Project而不是“定义投影”Define Projection。这两个工具的区别是血泪经验定义投影是“告诉软件这份数据是什么坐标系”不改坐标值投影是“把坐标值真正换算过去”。很多人坐标对不上就狂点定义投影结果越点越乱。正确顺序是先确认.prj内容对不对不对才用定义投影修正元数据再决定要不要投影转换。4. 字段清洗与拓扑检查路网能不能算看这一步4.1 属性表清洗空值、重复、等级缺失真实城市路网数据几乎不可能干净。常见问题是道路名称空值、同一路段重复录入、等级字段缺失。用 SQL 在 ogr2ogr 或 ArcGIS 里先做统计-- 统计各道路等级数量识别缺失和异常值 SELECT CLASS, COUNT(*) AS cnt FROM wuhan_road GROUP BY CLASS ORDER BY cnt DESC;如果CLASS出现空值或0、-1这类占位值说明等级没填全。处理策略分两种做可视化可以保留做加权分析比如按等级设通行速度就必须补全或剔除。我一般把缺失等级的路段单独导出成一个图层人工核对后再合并不直接删因为删了路网就断了。4.2 拓扑检查悬挂点和自相交是路径分析的隐形杀手路网做路径规划最怕两种几何问题悬挂点断头路没接上和自相交一条线自己绕回来。ArcGIS 里用“拓扑”工具建规则常见规则是“不能有悬挂点”“不能自相交”。在开源侧可以用shapely快速扫一遍import geopandas as gpd from shapely.validation import explain_validity gdf gpd.read_file(wuhan_road_proj.shp) # 检查无效几何 invalid gdf[~gdf.geometry.is_valid] for idx, row in invalid.iterrows(): print(idx, explain_validity(row.geometry))is_valid判断几何是否合法explain_validity给出具体原因比如Self-intersection。无效几何在构建网络数据集时会直接报错必须先buffer(0)修复或手工编辑。注意buffer(0)对线要素可能改变端点修复后要重新检查连通性。4.3 构建网络数据集前的连通性验证路径分析依赖路网在交叉口真正“接上”而不是视觉上挨着。判断方法是看交叉口处是否有共享节点。ArcGIS 的“整合”Integrate工具可以把容差内的端点吸附到一起但容差设大了会把平行道路误合并。我一般先用 1 米容差试跑完抽样看交叉口再决定是否调整。这一步没有后悔药整合前一定备份原始 shp。5. 避坑排查这五条我全踩过5.1 中文属性乱码现象ArcGIS 打开后道路名称显示为问号或方块。原因dbf 编码与系统默认编码不一致缺.cpg或.cpg写错。解决用ogrinfo确认编码或在 QGIS 里指定编码重新导出补一个内容为UTF-8的.cpg文件。5.2 坐标显示正常但长度算出来是 0.00x现象字段计算器算长度结果全是零点几。原因数据是地理坐标系单位是度长度按度算自然极小。解决先重投影到米制坐标系再算别在度坐标上直接算长度。5.3 路网看着连、路径分析却断开现象交叉口视觉上相接网络分析提示不连通。原因端点坐标有微小差异没有共享节点。解决用整合工具按小容差吸附或构建网络时设连通策略为“端点连通”并给容差。5.4 投影转换后位置整体偏移几百米现象转完投影和底图对不上整体平移。原因源坐标系基准判断错误或带号选错3 度带和 6 度带混用。解决核对中央经线武汉 3 度带是 114°E6 度带是 117°E别搞混。5.5 数据量大导致 ArcGIS 卡死现象拖入后软件无响应。原因路网要素数大且未建空间索引。解决先用ogrinfo看要素数超过十万级先建.sbn索引或切分图层别一次性全量渲染。6. 进阶把路网 shp 变成可计算的路网图6.1 用 networkx 把 shp 转成图结构做算法验证时我习惯把 shp 转成图脱离 GIS 软件跑。思路是把每条道路当边交叉口当节点。geopandas读几何networkx建图边权用长度或通行时间。import geopandas as gpd import networkx as nx from shapely.geometry import Point gdf gpd.read_file(wuhan_road_proj.shp) G nx.Graph() for idx, row in gdf.iterrows(): line row.geometry # 取线段首尾作为节点坐标四舍五入到 0.1 米避免浮点误差 start (round(line.coords[0][0], 1), round(line.coords[0][1], 1)) end (round(line.coords[-1][0], 1), round(line.coords[-1][1], 1)) length line.length # 投影坐标系下单位为米 G.add_edge(start, end, weightlength, namerow.get(NAME, )) print(nx.number_of_nodes(G), nx.number_of_edges(G))这里的关键是节点坐标取整浮点误差会让本该相同的交叉口变成两个节点图就断了。weight用几何长度后续跑nx.shortest_path就是最短路径。如果要做时间最短把 weight 换成length / speedspeed 按道路等级赋值。6.2 验证抽样对比 ArcGIS 网络分析结果图建好后别直接信抽几条起终点用nx.shortest_path_length算距离再和 ArcGIS 网络分析的最短路径对比误差应在米级以内。差异大就回头查连通性和节点取整精度。我一般抽 10 条跨区路径做校验全部吻合才敢往下用。6.3 一个具体技巧按等级分层抽稀做可视化全量路网渲染又慢又糊按等级分层是常规操作。快速路和主干路保留全部次干路和支路按长度阈值抽稀出图层次立刻清楚。用geopandas按CLASS过滤后分别导出再在制图软件里叠。这个技巧不改变分析结果只服务出图但能省掉大量渲染等待。从那以后我每次拿到新的城市路网 shp都强制先跑一遍ogrinfo -so看坐标系和字段再决定要不要投影、要不要清洗绝不直接拖进软件开干。这套顺序帮我省下的返工时间比任何快捷键都值。希望帮到你。本文还有配套的精品资源点击获取