简介这份资源面向GIS从业者、地理研究人员及交通规划相关专业学生提供新疆维吾尔自治区各级公路与地级市行政区划的Shapefile矢量数据可用于空间分析、地图制图与交通网络研究。压缩包共35个文件约1.9MB以shp、shx、dbf、prj、sbn、sbx、xml等类型为主其中shp存储公路与行政区划几何要素shx与sbn/sbx支撑空间索引dbf保存属性字段prj定义坐标系统xml记录元数据构成完整可用的GIS数据集。数据按一级、二级、三级公路及高速公路分层组织并附各地级市行政边界便于按等级提取路网、叠加行政区划做统计与制图。已有1248人学习下载适合需要新疆路网与区划底图进行空间分析、交通流量模拟或城市规划研究的读者直接调用。1. 新疆各级公路shp文件从一份路网数据到能跑通的空间分析做交通可达性、物流路径规划或者区域经济研究的人大概率都绕不开一个问题手头缺一份靠谱的路网底图。尤其是新疆这种面积占全国六分之一、路网稀疏又跨度极大的区域随便找一份全国路网数据切到新疆往往要么属性字段残缺要么县道乡道全丢要么坐标系对不上。新疆各级公路shp文件说的就是按行政等级高速、国道、省道、县道、乡道等分层组织的矢量线数据通常以 Shapefile 格式分发包含几何线串和道路等级、编号、名称等属性字段。它能解决的核心问题是让你在 QGIS、ArcGIS 或 PostGIS 里直接做缓冲区、最短路径、服务区覆盖这类分析而不用从零数字化。适合做区域规划、GIS 二次开发、交通建模的从业者也适合需要本地化路网底图的数据工程师。下面按「数据拿到手怎么验、怎么用、坑在哪」的顺序讲透。2. 拿到新疆各级公路shp文件后先做这四步体检一份路网 shp 文件能不能用不取决于它有多少条记录而取决于坐标系、拓扑、属性和完整性这四关。很多人拿到文件直接拖进软件就画图结果一叠加底图就偏到几百米外或者做路径分析时发现断头路遍地。这一章把体检流程拆开讲每一步都给出可复现的操作。2.1 坐标系与投影先确认是地理坐标还是投影坐标Shapefile 本身可以带.prj文件也可以不带。新疆跨越东经 73° 到 96°如果数据是 WGS84 地理坐标EPSG:4326直接量算长度会得到「度」而不是「米」缓冲区分析也会严重变形。常见做法是先看.prj内容再用gdalinfo或 Python 读一下。# 查看 shp 的坐标系和范围gdalinfo 是 GDAL 自带工具 gdalinfo xinjiang_roads.shp | head -40输出里重点看两行Coordinate System is:后面跟的是不是GEOGCS[WGS 84...]以及Extent的经纬度范围是否落在新疆合理区间经度约 73–96纬度约 34–49。如果.prj缺失gdalinfo会显示Coordinate System is:为空这时候不能猜必须找数据提供方确认或者用已知控制点反推。# 用 geopandas 快速检查坐标系和几何类型 import geopandas as gpd gdf gpd.read_file(xinjiang_roads.shp) print(CRS:, gdf.crs) # 应为 EPSG:4326 或投影坐标系 print(几何类型:, gdf.geom_type.unique()) # 路网应为 LineString 或 MultiLineString print(记录数:, len(gdf)) print(字段:, list(gdf.columns))逻辑说明gdf.crs为None说明没有坐标系定义后续所有空间运算都不可信geom_type如果混入 Point 或 Polygon说明数据被污染。参数上如果要做长度量算建议统一转到适合新疆的投影坐标系比如 CGCS2000 的 3 度带或 Albers 等积投影而不是留在 4326 上硬算。2.2 属性字段道路等级和编号是分析的主键新疆各级公路 shp 文件的价值一大半在属性表里。常见字段包括NAME道路名称、CLASS或GRADE等级、CODE编号如 G30、S21、TYPE高速/国道/省道/县道。如果只有几何没有等级字段那这份数据只能当底图看做不了分层分析。# 统计各等级道路数量和总长度 gdf_proj gdf.to_crs(EPSG:3857) # 临时投影仅用于快速量算 gdf_proj[length_km] gdf_proj.length / 1000 summary gdf_proj.groupby(CLASS).agg( 条数(CLASS, count), 总长_km(length_km, sum) ).sort_values(总长_km, ascendingFalse) print(summary)逻辑说明to_crs(EPSG:3857)是 Web 墨卡托在新疆高纬度地区长度会偏大这里只用于快速对比各等级占比正式量算应换成 Albers。参数上groupby的字段名要按你数据里实际的等级字段改常见是CLASS、GRADE、ROAD_TYPE三种。如果发现县道、乡道记录数极少说明这份数据可能只收录了国省道需要补数据。2.3 拓扑检查断头路和重复线是路径分析的致命伤路网做最短路径前提是线在节点处真正相连。Shapefile 不存储拓扑关系两条线端点差几米就算断开。常见做法是用 QGIS 的「拓扑检查器」或 PostGIS 的ST_Equals、ST_Intersects排查。-- 在 PostGIS 中找出端点未与其他线相接的悬挂线 SELECT a.id FROM roads a WHERE NOT EXISTS ( SELECT 1 FROM roads b WHERE a.id b.id AND ST_DWithin(ST_StartPoint(a.geom), b.geom, 0.0001) );逻辑说明ST_DWithin的容差0.0001度约等于 10 米用于容忍数字化误差。如果悬挂线数量占比超过 5%直接做路径分析会大量失败。解决方式是用ST_Snap做节点捕捉或者用 QGIS 的「v.clean」工具做拓扑修复。注意修复前务必备份原始 shp拓扑修复不可逆。2.4 完整性核对和公开路网里程做交叉验证新疆公路总里程有公开统计口径你可以把 shp 里各等级道路长度加总和统计值做量级对比。如果高速里程只有几百公里而实际新疆高速早已过万公里说明数据严重不全。这一步不需要精确到公里但量级差一个数量级就说明数据不能用。检查项合格标准不合格时的处理坐标系有明确 EPSG 或投影定义找提供方确认不猜几何类型全部为 LineString/MultiLineString过滤或剔除异常几何等级字段至少含高速/国道/省道三级补数据或降级使用悬挂线占比低于 5%拓扑修复后再用里程量级与公开统计同量级判定为不完整数据3. 用 Python 把新疆各级公路shp文件跑成可用路网体检通过之后下一步是把它变成能算最短路径、能算服务范围的网络。这一章给出一套从读取、投影、拓扑构建到导出的完整流程代码可以直接抄参数按你的数据改。3.1 投影转换把经纬度转成米制坐标做距离和缓冲区分析必须用投影坐标系。新疆常用 Albers 等积投影参数如下import geopandas as gpd # 定义适合新疆的 Albers 投影 albers_proj projaea lat_134 lat_248 lat_041 lon_085 x_00 y_00 datumWGS84 unitsm no_defs gdf gpd.read_file(xinjiang_roads.shp) gdf_albers gdf.to_crs(albers_proj) gdf_albers.to_file(xinjiang_roads_albers.shp, encodingutf-8)逻辑说明lat_1和lat_2是双标准纬线取新疆南北跨度 34° 和 48°lon_085是中央经线。这样投影后长度变形在千分之一以内满足路网量算精度。参数上如果你的分析只涉及南疆可以把lat_1/lat_2收窄到 36/42变形更小。导出时encodingutf-8防止中文路名乱码。3.2 构建网络数据集节点和边的提取Shapefile 本身不是网络需要把线在交点处打断生成节点和边。常见做法是用networkx配合geopandas或者直接用 PostGIS 的pgr_createTopology。import networkx as nx from shapely.geometry import LineString # 简化示例把每条线作为一条边端点作为节点 G nx.Graph() for idx, row in gdf_albers.iterrows(): line row.geometry if not isinstance(line, LineString): continue start line.coords[0] end line.coords[-1] length line.length / 1000 # 公里 G.add_edge(start, end, weightlength, road_ididx, namerow.get(NAME)) print(节点数:, G.number_of_nodes()) print(边数:, G.number_of_edges())逻辑说明这里用端点坐标元组作为节点标识简单但有个前提——相交的线必须在端点处精确重合。如果数据里存在「十字交叉但未打断」的情况需要先用shapely.ops.unary_union或 PostGIS 的ST_Node打断。参数上weight用长度如果你要考虑通行时间可以改成length / speed速度按等级赋值。3.3 按等级赋权高速、国道、县道不能一个速度真实路径规划里高速 120 km/h、国道 80 km/h、县道 40 km/h权重差异巨大。下面按等级字段批量赋速度。speed_map { 高速: 120, 国道: 80, 省道: 60, 县道: 40, 乡道: 30 } gdf_albers[speed] gdf_albers[CLASS].map(speed_map).fillna(30) gdf_albers[travel_time_h] gdf_albers[length_km] / gdf_albers[speed]逻辑说明map按等级映射速度fillna(30)给未知等级一个保守值。参数上速度值要按你的分析目的调整做货运和做应急响应取值不同。注意这里没有考虑实际路况、限速和拥堵属于静态近似适合宏观规划不适合实时导航。3.4 导出 GeoPackage比 Shapefile 更适合后续分析Shapefile 有字段名长度限制10 字符、不支持 UTF-8 中文、单文件 2GB 上限。分析完成后建议导出 GeoPackage。gdf_albers.to_file(xinjiang_roads.gpkg, layerroads, driverGPKG)逻辑说明GeoPackage 是 OGC 标准单文件、支持长字段名和中文、支持空间索引。参数上layer是图层名后续在 QGIS 里可以直接加载。如果你要入库 PostGIS用ogr2ogr或geopandas.to_postgis均可。4. 新疆各级公路shp文件避坑五条血泪经验这一章集中讲我在用新疆路网数据时踩过的坑每条按「现象 → 原因 → 解决」写能帮你省掉大量返工时间。4.1 叠加底图偏移几百米现象把 shp 拖到在线底图上路网整体偏移高速跑到戈壁里。原因数据是 WGS84但底图是 GCJ02 或 CGCS2000坐标系不一致。解决先确认 shp 的.prj如果是在线地图截图配准来的数据大概率是 GCJ02需要做坐标纠偏后再用。不要直接用「动态投影」糊弄偏移不会消失。4.2 路径分析大量失败现象用 networkx 或 pgRouting 算最短路提示「节点不可达」。原因线在交叉处没有打断端点不重合网络图里根本没有连通。解决先用ST_Node或 QGIS「线相交」工具打断再做拓扑构建。容差设 0.0001 度左右太大误连太小漏连。4.3 中文路名乱码现象属性表里NAME字段显示为问号或方块。原因Shapefile 的.dbf默认编码是 GBK 或 Latin1不是 UTF-8。解决读取时指定encodinggbk或直接转成 GeoPackage。如果已经乱码且无备份基本不可逆只能重新获取数据。4.4 县道乡道大量缺失现象统计发现县道只有几十条和实际不符。原因很多公开路网数据只收录国省道县乡道靠众包补充覆盖不全。解决把这份数据定位为「干线骨架」县乡道用 OSM 或其他众包数据补补完后做拓扑合并。不要指望一份 shp 解决所有等级。4.5 长度量算偏大或偏小现象同一段路在 4326 下量是 0.5 度在 3857 下量是 60 公里实际约 45 公里。原因地理坐标量算单位是度Web 墨卡托在高纬度放大。解决统一用 Albers 等积投影量算或者用geopy的测地线距离。记住任何长度分析前先问一句「当前坐标系单位是什么」。5. 进阶用新疆各级公路shp文件做服务区覆盖分析前面讲的是「数据能用」这一章讲「数据用出价值」。服务区覆盖分析是路网数据最典型的落地场景之一给定高速服务区点位算 30 分钟车程能覆盖哪些区域。这个分析在新疆特别有意义因为服务区间距大覆盖盲区直接关系到应急保障。5.1 用等时圈替代简单缓冲区简单缓冲区是直线距离路网分析必须用等时圈。思路是以服务区为起点沿路网按速度累积时间超过阈值就停止。import networkx as nx def isochrone(G, source, max_time_h): 基于 Dijkstra 的等时圈节点提取 lengths nx.single_source_dijkstra_path_length( G, source, cutoffmax_time_h, weighttravel_time_h ) return lengths # 返回 {节点: 累计时间} # 假设已有一个服务区节点 source_node reachable isochrone(G, source_node, max_time_h0.5) print(30 分钟可达节点数:, len(reachable))逻辑说明single_source_dijkstra_path_length按travel_time_h累积cutoff截断。参数上max_time_h0.5就是 30 分钟。注意这里假设服务区节点已经吸附到路网最近节点上吸附容差建议不超过 500 米否则等时圈起点就偏了。5.2 把可达节点转成覆盖面节点可达不等于面覆盖。常见做法是对可达边做凸包或缓冲合并生成覆盖多边形。from shapely.ops import unary_union # 收集可达边对应的几何 reachable_edges [] for u, v, data in G.edges(dataTrue): if u in reachable and v in reachable: reachable_edges.append(data[geometry]) coverage unary_union(reachable_edges).buffer(5000) # 5 公里缓冲逻辑说明buffer(5000)是给路网两侧各 5 公里覆盖带模拟道路服务范围。参数上缓冲距离按你的分析尺度定宏观规划 5–10 公里精细分析 1–2 公里。注意unary_union前要确保几何有效否则会报错。5.3 覆盖盲区识别与验证把覆盖面和新疆行政区划叠加差集就是盲区。验证方法随机抽几个盲区内的点反查最近服务区车程是否确实超过 30 分钟。验证项方法合格标准等时圈边界抽查边界点车程误差小于 10%盲区面积与总面积比值符合预期无异常大块服务区吸附检查吸附距离小于 500 米我自己的习惯是每次做完等时圈一定随机抽 5 个点手工验算因为路网权重里一个速度值填错整张图就全偏了。这个后悔药没有只能靠验证提前吃。希望帮到你。本文还有配套的精品资源点击获取