简介这份2024年全国地铁线路矢量数据面向GIS从业者、城市规划与交通研究人员以Shapefile标准格式提供全国各城市地铁线路的走向、站点分布与线路长度可直接用于线路规划、站点覆盖分析、人流模拟及应急疏散场景适合具备ArcGIS或QGIS基础的中高级用户。压缩包共9个文件除.shp几何数据、.dbf属性表、.shx索引、.prj坐标系定义外还包含.xml元数据、.md说明文档并附带.sbn/.sbx空间索引与.cpg编码文件整套约3.59MB结构紧凑加载后即可进行查询、统计与制图。已有682人学习下载数据基于2024年现状整理时效性与完整性较好。矢量格式可无限放大而不失真每条地铁线路作为独立图层呈现便于叠加路网、行政区划做缓冲区分析和客流模拟也可按线路筛选并制作专题地图为城市交通规划、线网评估和出行服务研究提供可靠的基础数据支撑。1. 一份全国地铁线路 shp拿到手先别急着拖进 ArcGIS做 WebGIS 或者做城市规划分析的人大概率都有过这种体验拿到一个“全国地铁线路矢量数据”的压缩包解压出来几个文件火急火燎拖进 ArcGIS结果要么线断成一截一截要么属性表里全是问号要么图直接跑到海上。真正的问题往往不是数据本身不对而是你不知道这份 shp 里装的是什么坐标系、什么拓扑结构、什么编码格式。这篇文章要拆的就是这类资源里最常见、也最实用的形态——一份整理过的全国地铁线路矢量数据shp 格式面向的是需要做地铁网络分析、线路可视化、站点周边分析或者地图发布的从业者。读完你能知道它内部结构长什么样怎么在各种 GIS 工具里正确加载以及真正干活时会在哪些地方翻车。2. 拆开这份矢量数据坐标系、字段与文件构成2.1 一份 shp 其实是一个“小文件柜”.shp / .shx / .dbf / .prj 各管什么很多人以为 shp 是一个文件其实它是一组文件。标准 ESRI Shapefile 由多个文件组成缺一个都会出问题。最基本的四个是文件后缀作用缺失后果.shp存储几何要素也就是线的坐标无几何图上没东西.shx几何索引用来快速定位要素部分软件打开会提示缺少索引QGIS 还能忍ArcGIS 很可能直接报错.dbf属性表存储线路名称、城市、开通年份等字段数据打不开或属性为空.prj坐标系描述文件记录投影和基准面软件按默认坐标系猜容易跑位另外还经常配套两个文件.cpg 用于声明属性表编码.sbn/.sbx 是空间索引。判断一份数据是否完整最简单的办法是看这四个主文件是否齐全。这一组文件的关系可以理解成.shp 存地图上的“形状”.dbf 存这张图对应的“表格”.shx 是查位置的“目录”.prj 是告诉你“这套坐标是怎么定义的”。四者必须配对出现单独拷贝一个 .shp 走到别的电脑上大概率打不开。在处理这类全国地铁线路数据时我一般会先用文件管理器把解压后的目录完整看一眼。光看名字不够重点看里面有没有 .prj没有 .prj 就要做好坐标系打架的心理准备。2.2 坐标系CGCS2000 与经纬度的关系国内公开发布的地理数据近几年绝大多数不再用北京 54、西安 80 那套老基准而是转向 CGCS2000China Geodetic Coordinate System 2000。地铁线路这类小比例尺全国数据常见情况是直接用经纬度坐标存储对应的坐标系记录通常是 GCS_China_Geodetic_Coordinate_System_2000 或 EPSG:4490CGCS2000 地理坐标系。这里要说明一个容易混淆的点如果一份数据是 CGCS2000 地理坐标系它的数值看起来就是经纬度范围大概在东经 73 到 135、北纬 3 到 54 之间。但如果你看到坐标值非常大比如横坐标几百万纵坐标几千万那很可能它被投影过采用的是高斯-克吕格投影或阿尔伯斯等积投影。对于全国地铁线路这种跨省跨市的数据我最推荐使用的是经纬度格式也就是 EPSG:4490 或 EPSG:4326WGS84。原因很简单和公开地图服务、在线底图对接时不需要额外投影转换导出 GeoJSON 或者叠加天地图都用得上。如果拿到手的数据没有 .prj 文件但坐标范围目测像经纬度我会先在软件里手动指定 EPSG:4490再检查位置是否合理。2.3 字段设计线路名、城市、开通年份是怎么落成属性的一份能直接干活的地铁线路数据属性表里通常会有这样几个字段字段典型含义示例值用处城市上海、广州按城市筛选、分色线路名称Line 1、1号线、M1标注、检索线路类型地铁、轻轨、磁浮区分制式开通年份1993、2024时序分析长度km41.7统计总里程需要提醒的是不同来源的数据字段命名差异很大有的叫 name有的叫 线路名称有的直接用拼音缩写。这也是为什么网上搜“省界线文件省 1 和省 2 的 shp 文件有什么区别”这类问题的人特别多——同一个名字的 shp字段设计、几何精度可能完全不同两个文件叠加到一起甚至对不上。拿到任何数据第一步不是画图而是把属性表完整翻一遍看字段含义、看统计范围。尤其不要迷信文件名。文件名只能说明主题说明不了坐标系、精度、更新日期。这份数据如果标题写的是 2024 年你需要重点确认属性表里是不是有 2024 年新开通的线路比如某些城市年底通车的段有没有补进去。3. 把 shp 从 zip 变成可用图层解压、校验、加载与坐标对齐3.1 解压前的三道检查文件列表、解压完整性、编码声明下载资源通常以 zip 压缩包形式提供。很多人拿到手直接双击解压解压失败之后就开始怀疑数据坏了。实际上 zip 解压出问题很大概率不是数据坏而是下载过程不完整。我一般先不解压而是在命令行里列压缩包内容unzip -l metro_lines_2024.zip这个命令会把 zip 内所有文件列出来我可以一眼看到里面是不是包含 .shp、.dbf、.prj、.shx以及有没有多余的文件。如果发现只有一个 .shp 孤零零躺在里面那就说明发布方打包时漏了配套文件后续要自己补坐标系定义。确认清单没问题之后再正式解压unzip -o metro_lines_2024.zip -d ./metro_data这里的-o表示覆盖已存在的同名文件-d指定解压目录。解压完成后不要急着关终端先看目录下有没有乱码文件。如果文件名出现一串无法识别的字符通常是 zip 打包时用了中文文件名且编码不是 UTF-8这在 Windows 上尤其常见。严格情况下解压后的 .dbf 字段也可能是 GBK 编码会在加载后表现为中文乱码这一点后面会专门讲。3.2 用 GDAL 读元数据加载之前先给数据“验明正身”每次拿到 shp我会先用 GDAL 看一眼它的元数据信息这是成本最低、也最有效的验货方式。GDAL 是开源 GIS 生态的事实标准工具集命令行工具ogrinfo可以直接读取 shp 的几何类型、要素数量、空间范围和坐标系ogrinfo -so -al metro_data/metro_lines.shps -o表示只输出概要信息索要元数据-al表示遍历所有图层。输出里会包含类似这样的信息INFO: Open of metro_lines.shp using driver ESRI Shapefile successful. Layer name: metro_lines Geometry: Line String Feature Count: 1200 Extent: (73.502355, 18.172138) - (134.4660 28, 54.3260115) Layer SRS: WGS_1984看到Geometry: Line String说明数据是线要素Feature Count代表总共有 1200 条记录Extent 如果落在我刚才说的经纬度合理区间说明坐标系大概率没有跑偏。这里有个细节值得说明如果输出显示Geometry: Unknown或Layer SRS: (unknown)那基本可以判定 .prj 缺失你需要之后手动定义坐标系。对于全国地级市级别的数据Feature Count 不应该只有几十条至少几百上千。如果明显过少说明它可能只包含了部分城市的地铁线不是完整全国版这对后续做分析有重大影响。3.3 在 QGIS 里加载尽力避免坐标系“猜错”导致的偏移QGIS 加载 shp 比较宽容即使缺少 .prj 也能打开但正因为它“好说话”很多人忽略了对坐标系的确认。正确做法是在加载后打开图层属性Layer Properties找到“信息”一栏看“CRS”确保它显示的是 EPSG:4490 或 EPSG:4326。如果打开后图形明显不对比如线路本该在广州结果跑到了内蒙古多半是因为软件用默认的 WGS84 去解释了一份投影坐标。这时候有两个工具要分清工具/操作作用使用场景设定图层坐标系Set Layer CRS强制让软件按指定坐标系读取这个图层数据本身坐标值没错只是缺 .prj需要告诉软件如何解释重新投影图层Reproject Layer把几何从当前坐标系转换成另一个坐标系数据坐标正确但你需要把它转到 Web 墨卡托或 WGS84很多人把这两个混淆是“玄学偏移”的高发原因。最简单的判断方法如果只是显示错误用“设定图层坐标系”修正如果坐标值本身不合预期用“重新投影”转换。在 ArcGIS 里对应的是“定义投影Define Projection”和“投影Project”两个工具概念一模一样。如果你是在 ArcGIS Pro 里操作路径是 Geoprocessing → Define Projection。3.4 和底图对齐把线路数据叠加到在线地图上的尝试加载完成后下一步通常是对齐到在线底图上。QGIS 里最常见的操作是使用 QuickMapServices 插件或 XYZ Tiles添加 OpenStreetMap 标准底图然后把地铁线路拖到上方。如果使用 CGCS2000 的地理坐标系叠加 OSM 时 QGIS 会自动进行动态投影对齐效果是正常的。但我遇到过一个特殊情况某些版本 QGIS 的动态投影运算较慢叠加大范围底图时刷新非常卡顿。这时我的习惯是直接把地铁线路数据转成 EPSG:3857Web 墨卡托让图层和在线底图保持同一坐标系性能会好很多ogr2ogr -t_srs EPSG:3857 metro_lines_3857.shp metro_data/metro_lines.shpogr2ogr的-t_srs参数指定目标坐标系。输出文件名不能和输入文件名相同否则会报错。转换完成后再用ogrinfo看一眼范围确认坐标已经变成以米为单位的 Web 墨卡托值。4. 属性与几何对齐编码、字段加工与空间配准4.1 属性表中文乱码GBK 与 UTF-8 的纠葛这是 shp 数据在国内最常见的坑。很多原始数据产自 ArcGIS 中文环境属性表被写成 GBK 编码而开源的 QGIS 和部分 Python 库默认按 UTF-8 读取结果就是打开属性表之后满屏问号。处理方式取决于你用什么工具。QGIS 里可以右键图层 → Layer Properties → Source修改“Data Source Encoding”把默认的 UTF-8 改成 GBK乱码通常会立刻恢复。ArcGIS Pro 对中文编码支持好一些但偶尔也会遇见字段名正常、字段值乱码的情况。更专业的处理是用 Python 做批量修复。这里给一段简易脚本思路是先用 pyshp 读取原始数据把字段值重新编码后写成一个新的 shpimport shapefile # 读取 GBK 编码的原始 shp 文件 reader shapefile.Reader(metro_lines.shp, encodinggbk) # 创建新的 shp 写入对象 writer shapefile.Writer(metro_lines_utf8.shp) # 复制字段定义 writer.fields reader.fields # 逐条读取几何和属性重新写入 for record in reader.records(): writer.record(*record) for shape in reader.shapes(): writer.shape(shape) writer.close()这段代码实际上不会做编码转换因为 pyshp 在读取时就通过encodinggbk参数解出了正确字符串写出时默认按 UTF-8 存。核心点在于告诉 reader 原始编码是 GBK输出就变成 UTF-8。字段名如果也是中文且被截断需要额外处理字段名长度问题下一节会说明。更省事的方案是使用 GDAL 的编码修复能力。常见的做法是找到图层根目录的.cpg文件直接把它的内容改成UTF-8。但要注意.cpg 文件只是在告诉软件“这个 .dbf 用 UTF-8 编码”如果实际内容本身是 GBK改了反而会乱得更厉害。正确的做法是先确认原始编码是什么再对症下药。4.2 字段名缩短DBF 的 10 字符限制DBF 格式对字段名长度有硬性限制通常是 10 个字节。如果一个字段名是中文比如“线路名称”在部分环境下会被截断成很短的乱码。这是业内流传已久的痛点尤其从 ArcGIS 导出到其他平台时尤其明显。我的习惯是在数据整理阶段就把字段名改成英文并且控制在 10 字节以内。例如原始含义字段名建议线路名称name所属城市city开通年份open_year线路长度length_km数据来源source这样做的代价是在 ArcGIS 里做地图标注时需要手动关联中文字段别名。但好处是所有下游工具都能稳定读取不会因为字段名长度被截断而丢属性。4.3 线和点的对齐方式从线路数据引出站点数据的使用注意如果你拿到的资源只有地铁线路Line没有站点Point而你的分析又需要站点数据这时别急着骂资源不全。线路数据中往往隐含了站点信息但需要你自行提取。常见做法是把线路数据在 QGIS 里使用“Extract Vertices”工具将折线上的所有节点提取成点。注意这个方法提取出来的点包括转弯点、线段端点并不等于真实站点。真实站点通常在线路的拐角或交叉处但也不完全是。要想准确获取站点还是得靠属性表里的站点名字段或另外下载站点数据。这里必须提醒一个限制地铁线路的实际走向可能非常弯曲提取出的节点数量会远大于站点数量。如果你直接用节点做缓冲区分析结果会失真。我一般会先用ogr2ogr或 QGIS 的简化工具Simplify对线路做轻度平滑再用节点提取获得相对合理的站点位置。简化程度需要反复试探过度简化会破坏线路走向。5. 避坑这份数据落地最容易翻车的五件事5.1 压缩包解压报“已损坏”现象从网盘或下载站拿到的 zip 解压到一半提示文件损坏或者解压出来只有部分文件。原因大概率是下载中途断流文件不完整或者浏览器下载时没有按“完整文件”保存。少部分情况下是发布方打包时把文件打进了多层目录解压软件没正确处理。解决先执行unzip -l 文件名.zip检查目录结构再用unzip -t 文件名.zip做完整性测试。如果提示 archive is not a valid zip重新下载一次换用支持断点续传的下载工具避免使用浏览器自带的临时下载机制。5.2 打开后线路跑到海上现象线路图形完全不在中国范围跑到非洲或者太平洋。原因最常见的是 .prj 文件缺失GIS 软件默认按 WGS84 解析而数据实际是 CGCS2000 投影坐标。少数情况下是 .prj 存在但描述不标准软件读错了。解决先运行ogrinfo看 Extent判断坐标值是否符合经纬度。如果是投影坐标用“定义投影”指定正确的坐标系例如 CGCS2000 / 3-degree Gauss-Kruger zone。拿不准时逐一带可能的投影定义看地图位置位置对了之后再用“投影”工具把它转成经纬度。5.3 属性表字段全是问号现象地图显示正常但打开属性表中文全部变成乱码。原因.dbf 文件是 GBK 编码而读取软件默认 UTF-8。这两个编码在中文环境下造成的乱码几乎是必然的。解决QGIS 里图层属性 → Source → 修改栏目编码为 GBK。ArcGIS Pro 里可以右键图层属性查看“代码页”将其从 UTF-8 改为 GBK。更一劳永逸的做法是用第 4.1 节的方法重新导出成 UTF-8 编码的 shp后续所有工具都省心。5.4 下载的文件数比预期少缺少几何文件现象解压后目录里只有 .dbf、.prj但找不到 .shp。原因某些网盘服务把多文件打包时过滤掉了一些后缀名或者杀毒软件把 .shp 当作未知文件隔离了。解决先去杀毒软件的隔离区恢复文件然后检查发布方的压缩包是不是用单独一层目录打包。如果确实缺少 .shp 文件那这个资源不完整只能重新下载或找替代源。不要试图用其他文件伪造替换shp 和 dbf 是严格匹配的图层明确不匹配会报字段数不匹配。5.5 线路断成很多短线无法做网络分析现象在地图上看起来一条线路被拆分成无数小线段每条线单独一个要素做缓冲区或网络分析时结果破碎。原因这个现象不一定是数据坏了。很多 shp 在输出时会按几何类型拆分比如把一条长线按路段拆分方便维护与更新。但如果分析需求是按照整条线路算长度、做连通性分析这种分段方式就不合适。解决QGIS 里使用“Dissolve”融合工具按城市字段或线路名字段融合线段。融合后一个要素代表一条完整线路。融合前需要确认字段值唯一比如线路名称是否存在同名问题否则会把不同城市的两条线路错误融合成一条。6. 让数据真正可用转 GeoJSON、压缩优化和一套自查习惯拿到这份全国地铁线路 shp 后真正的工程价值在把它喂给地图引擎或分析系统。这里分享一个我常用的收尾流程数据严格检查 → 转 GeoJSON → 压缩瘦身。GeoJSON 是目前 WebGIS 前端最友好的格式Leaflet、Mapbox、OpenLayers 直接加载不需要后端处理。用 ogr2ogr 转换一行命令搞定ogr2ogr -f GeoJSON -t_srs EPSG:4326 metro_lines_2024.json metro_lines.shp这一步把几何和属性完整保留坐标转成经纬度。如果不转投影直接输出前端加载后还需要额外转换容易引入偏移。GeoJSON 文件体积可能会比 shp 大得多尤其这种全国线要素字段多、坐标点密动不动几十上百兆。这时候可以用 topojson-cli 做拓扑压缩把重复坐标抽掉数据量经常能减少七成以上。常见的做法是先安装 TopoJSON 工具再用命令行压缩npx topojson -o metro_lines_topo.json -p metro_lines.json这里-p表示保留所有属性数据。压缩后的文件在 Mapbox 和 ECharts 里直接可用加载速度明显提升。注意TopoJSON 的坐标是量化过的精度适合可视化不适合做精确测量。如果下游分析要求厘米级精度继续用 GeoJSON 或原始 shp。这套流程我重复了无数次现在已经固定成习惯拿到任何一份国界、铁路、道路类的 shp第一步验证 zip 完整性第二步跑 ogrinfo 确认坐标系和要素数第三步加载后对底图看偏移第四步查属性表编码第五步再谈转格式和配图。坑踩得多了就会明白省掉任意一步都可能在后半程加倍补回来。希望这篇文章能让你拿到这份 2024 年全国地铁线路 shp 时少走一段弯路。本文还有配套的精品资源点击获取