尧图网站建设 
简介这份雅江流域矢量数据集面向地理信息、水文生态与遥感测绘方向的研究人员和学生提供雅鲁藏布江流域的空间边界与要素表达可用于流域管理、洪水风险评估、水资源配置及生态保育规划等GIS分析场景。资源包共8个文件压缩后约1.31MB包含shp主文件、shx索引、dbf属性表、prj投影定义以及sbn、sbx空间索引和xml元数据、png预览图覆盖矢量数据读取、坐标定义与属性查询的完整环节。目前已有332人学习下载。数据参考中国三级流域范围空间分布矢量资料及相关研究论文整理能清晰呈现河流位置、流向与流域边界便于直接导入ArcGIS、QGIS等平台开展空间分析与模型运算为西藏地区地理、水文与生态研究提供实证基础。1. 雅江流域矢量数据集从 shp 文件到流域分析的落地路径雅江流域矢量数据集 shp 数据说白了就是把雅鲁藏布江流域的边界、水系、子流域划分这些地理要素以 Shapefile 格式打包好让你能直接拖进 QGIS 或 ArcGIS 里做空间分析。很多人第一次搜“雅鲁藏布江数据下载”往往是因为手头有个流域水文、生态评估或者遥感影像裁剪的需求结果发现公开渠道要么只有整块青藏高原的粗边界要么分辨率低到没法用。这个数据集解决的就是“精确到子流域级别”的问题——它通常包含流域面、干流线、支流线、湖泊面等图层适合做水文建模、淹没分析、土地利用统计的从业者。如果你正在找一份能直接参与计算、而不是只能看个大概的矢量底图那接下来的内容就是帮你把这份数据从下载到用起来的路走通。2. 雅江流域 shp 数据到底包含什么图层拆解与选型判断2.1 常见图层构成与属性字段一份可用的雅江流域矢量数据集通常不会只有一个boundary.shp。我拿到过的版本里至少包含以下四类图层每类都有明确的用途图层文件名示例几何类型关键属性字段典型用途yajiang_basin.shp面BASIN_ID,AREA_KM2,NAME整体裁剪、统计总面积yajiang_sub_basins.shp面SUB_ID,UP_AREA,ORDER子流域划分、水文响应单元yajiang_rivers.shp线RIVER_NAME,STRAHLER,LENGTH_KM河网提取、流向分析yajiang_lakes.shp面LAKE_NAME,AREA_KM2湖泊调蓄、水面变化监测属性字段的命名在不同来源里会有差异但核心逻辑一致AREA_KM2用来快速筛选面积阈值STRAHLER是 Strahler 分级数值越大代表主流级别越高。如果你要做水文站点的上游汇水区提取UP_AREA这个字段比你自己用 DEM 算要省事得多。2.2 坐标系与分辨率为什么直接拖进去会偏雅江流域横跨多个 UTM 带常见做法是提供 WGS84 地理坐标系EPSG:4326作为交换格式但真正做面积量算和距离分析时必须投影到适合的平面坐标系。我一般会先检查.prj文件里写的是什么# 查看 shp 文件的坐标系信息 ogrinfo -al -so yajiang_basin.shp | grep -A 5 Coordinate System如果输出是GCS_WGS_1984说明是地理坐标系单位是度。这时候你直接算面积得到的是平方度毫无意义。正确的做法是投影到EPSG:32645UTM Zone 45N或EPSG:32646具体看流域主体落在哪个带。雅江大部分在 Zone 45N 和 46N 之间我通常选 45N 做整体分析误差在可接受范围内。注意不要用 Web MercatorEPSG:3857做面积统计高纬度地区变形极大雅江流域用它会让你算出的面积比实际大出百分之十几。2.3 数据来源与版本判断怎么确认拿到的是能用的公开渠道能搜到的雅江矢量数据质量参差不齐。我判断一份数据能不能用的标准有三条第一子流域边界是否闭合有没有悬挂节点第二干流线是否连续有没有断成好几截第三属性表里有没有AREA_KM2这种能直接参与计算的字段。如果拿到的是从 DEM 自动提取的河网线要素往往会有大量毛刺需要先做平滑和拓扑检查。import geopandas as gpd # 读取子流域面 sub gpd.read_file(yajiang_sub_basins.shp) # 检查几何有效性 invalid sub[~sub.is_valid] print(f无效几何数量: {len(invalid)}) # 检查面积字段是否合理 if AREA_KM2 in sub.columns: print(sub[AREA_KM2].describe()) else: # 没有面积字段就自己算但必须先投影 sub_proj sub.to_crs(epsg32645) sub[AREA_KM2] sub_proj.area / 1e6 print(sub[AREA_KM2].describe())这段代码先做几何有效性检查再判断有没有现成的面积字段。如果没有就投影后自己算。to_crs(epsg32645)这一步是关键不投影直接.area得到的是平方度数值会小得离谱。/1e6是把平方米转成平方公里。跑完看describe()的输出如果最大面积和最小面积差了几个数量级说明子流域划分粒度不均匀用的时候要留意。3. 从下载到入库雅江 shp 数据的获取与预处理流程3.1 下载渠道的甄别与文件完整性检查搜“雅鲁藏布江数据下载”能出来一堆结果但很多是整块青藏高原的边界真正细化到雅江流域子流域级别的并不多。我一般会优先找附带元数据说明的压缩包里面至少要有.shp、.shx、.dbf、.prj四个文件。缺.prj的坐标系靠猜后期对不上就麻烦了。拿到压缩包后先解压再检查文件完整性# 列出压缩包内容确认 shp 相关文件齐全 unzip -l yajiang_vector.zip # 解压后检查每个图层的文件配套情况 for f in *.shp; do base${f%.shp} for ext in shp shx dbf prj; do if [ ! -f ${base}.${ext} ]; then echo 缺失: ${base}.${ext} fi done done第一行列出压缩包内容确认没有漏掉.prj或.dbf。第二段循环检查每个 shp 文件是否配套完整。.shx是索引文件缺了读不了.dbf存属性表缺了只剩空几何.prj定义坐标系缺了就得手动指定。这三个缺任何一个后期都会翻车。3.2 用 QGIS 做拓扑检查与修复QGIS 是检查 shp 数据最直观的工具。把子流域面拖进去先看有没有明显的缝隙或重叠。我通常用“拓扑检查器”跑一遍规则选“不能有重叠”和“不能有缝隙”。如果发现重叠用“融合”工具合并如果有缝隙用“捕捉几何”修。import geopandas as gpd from shapely.ops import unary_union # 读取子流域 sub gpd.read_file(yajiang_sub_basins.shp) # 检查重叠两两求交看有没有面积大于阈值的交集 overlaps [] for i, geom_i in enumerate(sub.geometry): for j, geom_j in enumerate(sub.geometry): if i j and geom_i.intersects(geom_j): inter geom_i.intersection(geom_j) if inter.area 1e-8: # 忽略极小交集 overlaps.append((i, j, inter.area)) print(f存在重叠的图斑对: {len(overlaps)}) # 如果有重叠用 unary_union 融合后再拆 if overlaps: merged unary_union(sub.geometry) sub_fixed gpd.GeoDataFrame(geometry[merged], crssub.crs) sub_fixed.to_file(yajiang_sub_basins_fixed.shp)这段代码做两两相交检查inter.area 1e-8是为了过滤掉边界接触产生的零面积交集。如果确实有重叠unary_union把所有面融成一个整体再根据需要重新拆分。注意unary_union之后属性表会丢失所以修复前先备份原始属性修复后用空间连接把属性挂回去。3.3 裁剪与投影把雅江数据对齐到你的研究区如果你只关心雅江的某个子区域比如拉萨河下游就需要用裁剪工具。我一般用geopandas的clip函数比 QGIS 里点来点去更可控。import geopandas as gpd # 读取雅江全流域和你的研究区边界 yajiang gpd.read_file(yajiang_basin.shp) study_area gpd.read_file(lhasa_downstream.shp) # 统一坐标系 study_area study_area.to_crs(yajiang.crs) # 裁剪 clipped gpd.clip(yajiang, study_area) # 投影到 UTM 45N 后保存 clipped_proj clipped.to_crs(epsg32645) clipped_proj.to_file(yajiang_clipped_utm45.shp)to_crs(yajiang.crs)确保两个图层坐标系一致否则clip会报错或结果错位。gpd.clip保留的是第一个图层中落在第二个图层范围内的部分。最后投影到 UTM 45N 再保存方便后续做面积和距离计算。如果你要裁剪的是线图层同样逻辑把面换成线即可。4. 雅江 shp 数据避坑指南五个让我返工的血泪教训4.1 坑一坐标系缺失导致面积算错十倍现象用geopandas读入 shp 后直接.area得到的数值和实际面积对不上差了一个数量级甚至更多。原因shp 文件没有.prj或者.prj写的是地理坐标系度但代码里没做投影转换就直接算面积。解决先gpd.read_file()后检查.crs如果是None或EPSG:4326必须.to_crs(epsg32645)再算面积。我现在的习惯是读入后第一件事就打印gdf.crs确认无误再往下走。4.2 坑二子流域边界不闭合水文分析直接断流现象做流向分析或上游汇水区提取时河网在某个子流域边界处突然断开下游拿不到上游的数据。原因子流域面之间有缝隙或者面边界和河网线没有严格对齐导致空间连接时匹配失败。解决用unary_union把所有子流域面融成一个整体再和河网做空间连接。或者用 QGIS 的“捕捉几何”工具把面边界顶点捕捉到河网线上。我一般会在融合后重新用polygonize拆分确保边界共享。4.3 坑三属性表编码乱码中文河名变问号现象.dbf文件里的中文字段在 QGIS 里显示正常但用geopandas读出来全是乱码或问号。原因.dbf的默认编码是GBK或Latin-1而geopandas默认按UTF-8读。解决读入时指定编码gdf gpd.read_file(yajiang_rivers.shp, encodinggbk)如果gbk不行就试latin-1再不行用chardet检测。我通常会在读入后把字段名和值都打印一遍确认没有乱码再继续。4.4 坑四线图层毛刺太多平滑后拓扑又断了现象从 DEM 提取的河网线锯齿严重做缓冲区分析时边界不平滑但用平滑算法处理后线又断成几截。原因平滑算法如simplify会改变线的几何形状如果容差设得太大相邻线段可能不再相交。解决先用小容差simplify(0.0001)做轻微平滑再用linemerge把断开的线合并。如果还不行就回到原始线用snap把端点捕捉到最近的其他线上。我一般会保留一份原始线作为备份平滑后的线只用于可视化分析还是用原始数据。4.5 坑五下载的 shp 缺少.shxQGIS 能开但代码读不了现象QGIS 里能正常加载 shp但geopandas.read_file()报错Unable to open ... .shx。原因.shx是索引文件QGIS 有时能容错但GDAL库严格要求它存在。解决如果压缩包里没有.shx用ogr2ogr重新生成ogr2ogr -f ESRI Shapefile output.shp input.shp这条命令会重新创建完整的 shp 文件集包括.shx。注意output.shp和input.shp不能同名否则会覆盖。5. 进阶技巧用雅江 shp 数据做子流域批量统计与自动化出图5.1 批量统计每个子流域的土地利用面积占比拿到雅江子流域 shp 后最常见的进阶需求是每个子流域里草地、林地、水体各占多少。手动一个个算太慢用geopandas的overlay加groupby可以批量出结果。import geopandas as gpd import pandas as pd # 读取子流域和土地利用数据 sub gpd.read_file(yajiang_sub_basins.shp).to_crs(epsg32645) landuse gpd.read_file(landuse_2020.shp).to_crs(epsg32645) # 叠加分析把土地利用切到每个子流域里 overlay gpd.overlay(landuse, sub, howintersection) # 计算每个子流域内各地类的面积 overlay[area_km2] overlay.area / 1e6 stats overlay.groupby([SUB_ID, LU_TYPE])[area_km2].sum().reset_index() # 透视成宽表每行一个子流域每列一个地类 pivot stats.pivot(indexSUB_ID, columnsLU_TYPE, valuesarea_km2).fillna(0) # 算占比 pivot_pct pivot.div(pivot.sum(axis1), axis0) * 100 pivot_pct.to_csv(sub_basin_landuse_pct.csv)gpd.overlay的howintersection保留两个图层的交集部分土地利用面被切碎到每个子流域里。groupby按子流域 ID 和地类汇总面积pivot把长表转宽表最后div算占比。输出的 CSV 可以直接拿去画堆叠柱状图。注意to_crs必须在overlay之前做否则面积计算还是错的。5.2 自动化出图用 matplotlib 批量渲染子流域专题图如果要把每个子流域单独出一张图用循环加matplotlib可以省掉大量重复劳动。import matplotlib.pyplot as plt fig, axes plt.subplots(3, 4, figsize(16, 12)) axes axes.flatten() for idx, (sub_id, group) in enumerate(sub.groupby(SUB_ID)): ax axes[idx] group.plot(axax, colorlightblue, edgecolorblack) ax.set_title(fSub-basin {sub_id}) ax.set_axis_off() plt.tight_layout() plt.savefig(sub_basins_grid.png, dpi150)这段代码把子流域按SUB_ID分组每组画在一个子图里。3x4的网格适合 12 个子流域多了就调整行列数。set_axis_off()去掉坐标轴出图更干净。保存时dpi150兼顾清晰度和文件大小。如果子流域数量不固定可以用math.ceil动态算行列。5.3 一个我常用的检查习惯每次拿到新的雅江 shp 数据我都会先跑一遍“三查”查坐标系、查几何有效性、查属性表字段。这三步花不了五分钟但能避免后面几个小时的返工。具体就是def quick_check(path): gdf gpd.read_file(path) print(fCRS: {gdf.crs}) print(f无效几何: {len(gdf[~gdf.is_valid])}) print(f字段: {list(gdf.columns)}) print(f要素数: {len(gdf)}) quick_check(yajiang_basin.shp)这个函数打印坐标系、无效几何数量、字段列表和要素总数。坐标系不对就投影几何无效就修复字段缺失就补算。养成这个习惯后我很少再因为数据本身的问题卡住。希望帮到你。本文还有配套的精品资源点击获取 雅江流域矢量数据集shp：从下载到流域分析的完整避坑指南
简介这份雅江流域矢量数据集面向地理信息、水文生态与遥感测绘方向的研究人员和学生提供雅鲁藏布江流域的空间边界与要素表达可用于流域管理、洪水风险评估、水资源配置及生态保育规划等GIS分析场景。资源包共8个文件压缩后约1.31MB包含shp主文件、shx索引、dbf属性表、prj投影定义以及sbn、sbx空间索引和xml元数据、png预览图覆盖矢量数据读取、坐标定义与属性查询的完整环节。目前已有332人学习下载。数据参考中国三级流域范围空间分布矢量资料及相关研究论文整理能清晰呈现河流位置、流向与流域边界便于直接导入ArcGIS、QGIS等平台开展空间分析与模型运算为西藏地区地理、水文与生态研究提供实证基础。1. 雅江流域矢量数据集从 shp 文件到流域分析的落地路径雅江流域矢量数据集 shp 数据说白了就是把雅鲁藏布江流域的边界、水系、子流域划分这些地理要素以 Shapefile 格式打包好让你能直接拖进 QGIS 或 ArcGIS 里做空间分析。很多人第一次搜“雅鲁藏布江数据下载”往往是因为手头有个流域水文、生态评估或者遥感影像裁剪的需求结果发现公开渠道要么只有整块青藏高原的粗边界要么分辨率低到没法用。这个数据集解决的就是“精确到子流域级别”的问题——它通常包含流域面、干流线、支流线、湖泊面等图层适合做水文建模、淹没分析、土地利用统计的从业者。如果你正在找一份能直接参与计算、而不是只能看个大概的矢量底图那接下来的内容就是帮你把这份数据从下载到用起来的路走通。2. 雅江流域 shp 数据到底包含什么图层拆解与选型判断2.1 常见图层构成与属性字段一份可用的雅江流域矢量数据集通常不会只有一个boundary.shp。我拿到过的版本里至少包含以下四类图层每类都有明确的用途图层文件名示例几何类型关键属性字段典型用途yajiang_basin.shp面BASIN_ID,AREA_KM2,NAME整体裁剪、统计总面积yajiang_sub_basins.shp面SUB_ID,UP_AREA,ORDER子流域划分、水文响应单元yajiang_rivers.shp线RIVER_NAME,STRAHLER,LENGTH_KM河网提取、流向分析yajiang_lakes.shp面LAKE_NAME,AREA_KM2湖泊调蓄、水面变化监测属性字段的命名在不同来源里会有差异但核心逻辑一致AREA_KM2用来快速筛选面积阈值STRAHLER是 Strahler 分级数值越大代表主流级别越高。如果你要做水文站点的上游汇水区提取UP_AREA这个字段比你自己用 DEM 算要省事得多。2.2 坐标系与分辨率为什么直接拖进去会偏雅江流域横跨多个 UTM 带常见做法是提供 WGS84 地理坐标系EPSG:4326作为交换格式但真正做面积量算和距离分析时必须投影到适合的平面坐标系。我一般会先检查.prj文件里写的是什么# 查看 shp 文件的坐标系信息 ogrinfo -al -so yajiang_basin.shp | grep -A 5 Coordinate System如果输出是GCS_WGS_1984说明是地理坐标系单位是度。这时候你直接算面积得到的是平方度毫无意义。正确的做法是投影到EPSG:32645UTM Zone 45N或EPSG:32646具体看流域主体落在哪个带。雅江大部分在 Zone 45N 和 46N 之间我通常选 45N 做整体分析误差在可接受范围内。注意不要用 Web MercatorEPSG:3857做面积统计高纬度地区变形极大雅江流域用它会让你算出的面积比实际大出百分之十几。2.3 数据来源与版本判断怎么确认拿到的是能用的公开渠道能搜到的雅江矢量数据质量参差不齐。我判断一份数据能不能用的标准有三条第一子流域边界是否闭合有没有悬挂节点第二干流线是否连续有没有断成好几截第三属性表里有没有AREA_KM2这种能直接参与计算的字段。如果拿到的是从 DEM 自动提取的河网线要素往往会有大量毛刺需要先做平滑和拓扑检查。import geopandas as gpd # 读取子流域面 sub gpd.read_file(yajiang_sub_basins.shp) # 检查几何有效性 invalid sub[~sub.is_valid] print(f无效几何数量: {len(invalid)}) # 检查面积字段是否合理 if AREA_KM2 in sub.columns: print(sub[AREA_KM2].describe()) else: # 没有面积字段就自己算但必须先投影 sub_proj sub.to_crs(epsg32645) sub[AREA_KM2] sub_proj.area / 1e6 print(sub[AREA_KM2].describe())这段代码先做几何有效性检查再判断有没有现成的面积字段。如果没有就投影后自己算。to_crs(epsg32645)这一步是关键不投影直接.area得到的是平方度数值会小得离谱。/1e6是把平方米转成平方公里。跑完看describe()的输出如果最大面积和最小面积差了几个数量级说明子流域划分粒度不均匀用的时候要留意。3. 从下载到入库雅江 shp 数据的获取与预处理流程3.1 下载渠道的甄别与文件完整性检查搜“雅鲁藏布江数据下载”能出来一堆结果但很多是整块青藏高原的边界真正细化到雅江流域子流域级别的并不多。我一般会优先找附带元数据说明的压缩包里面至少要有.shp、.shx、.dbf、.prj四个文件。缺.prj的坐标系靠猜后期对不上就麻烦了。拿到压缩包后先解压再检查文件完整性# 列出压缩包内容确认 shp 相关文件齐全 unzip -l yajiang_vector.zip # 解压后检查每个图层的文件配套情况 for f in *.shp; do base${f%.shp} for ext in shp shx dbf prj; do if [ ! -f ${base}.${ext} ]; then echo 缺失: ${base}.${ext} fi done done第一行列出压缩包内容确认没有漏掉.prj或.dbf。第二段循环检查每个 shp 文件是否配套完整。.shx是索引文件缺了读不了.dbf存属性表缺了只剩空几何.prj定义坐标系缺了就得手动指定。这三个缺任何一个后期都会翻车。3.2 用 QGIS 做拓扑检查与修复QGIS 是检查 shp 数据最直观的工具。把子流域面拖进去先看有没有明显的缝隙或重叠。我通常用“拓扑检查器”跑一遍规则选“不能有重叠”和“不能有缝隙”。如果发现重叠用“融合”工具合并如果有缝隙用“捕捉几何”修。import geopandas as gpd from shapely.ops import unary_union # 读取子流域 sub gpd.read_file(yajiang_sub_basins.shp) # 检查重叠两两求交看有没有面积大于阈值的交集 overlaps [] for i, geom_i in enumerate(sub.geometry): for j, geom_j in enumerate(sub.geometry): if i j and geom_i.intersects(geom_j): inter geom_i.intersection(geom_j) if inter.area 1e-8: # 忽略极小交集 overlaps.append((i, j, inter.area)) print(f存在重叠的图斑对: {len(overlaps)}) # 如果有重叠用 unary_union 融合后再拆 if overlaps: merged unary_union(sub.geometry) sub_fixed gpd.GeoDataFrame(geometry[merged], crssub.crs) sub_fixed.to_file(yajiang_sub_basins_fixed.shp)这段代码做两两相交检查inter.area 1e-8是为了过滤掉边界接触产生的零面积交集。如果确实有重叠unary_union把所有面融成一个整体再根据需要重新拆分。注意unary_union之后属性表会丢失所以修复前先备份原始属性修复后用空间连接把属性挂回去。3.3 裁剪与投影把雅江数据对齐到你的研究区如果你只关心雅江的某个子区域比如拉萨河下游就需要用裁剪工具。我一般用geopandas的clip函数比 QGIS 里点来点去更可控。import geopandas as gpd # 读取雅江全流域和你的研究区边界 yajiang gpd.read_file(yajiang_basin.shp) study_area gpd.read_file(lhasa_downstream.shp) # 统一坐标系 study_area study_area.to_crs(yajiang.crs) # 裁剪 clipped gpd.clip(yajiang, study_area) # 投影到 UTM 45N 后保存 clipped_proj clipped.to_crs(epsg32645) clipped_proj.to_file(yajiang_clipped_utm45.shp)to_crs(yajiang.crs)确保两个图层坐标系一致否则clip会报错或结果错位。gpd.clip保留的是第一个图层中落在第二个图层范围内的部分。最后投影到 UTM 45N 再保存方便后续做面积和距离计算。如果你要裁剪的是线图层同样逻辑把面换成线即可。4. 雅江 shp 数据避坑指南五个让我返工的血泪教训4.1 坑一坐标系缺失导致面积算错十倍现象用geopandas读入 shp 后直接.area得到的数值和实际面积对不上差了一个数量级甚至更多。原因shp 文件没有.prj或者.prj写的是地理坐标系度但代码里没做投影转换就直接算面积。解决先gpd.read_file()后检查.crs如果是None或EPSG:4326必须.to_crs(epsg32645)再算面积。我现在的习惯是读入后第一件事就打印gdf.crs确认无误再往下走。4.2 坑二子流域边界不闭合水文分析直接断流现象做流向分析或上游汇水区提取时河网在某个子流域边界处突然断开下游拿不到上游的数据。原因子流域面之间有缝隙或者面边界和河网线没有严格对齐导致空间连接时匹配失败。解决用unary_union把所有子流域面融成一个整体再和河网做空间连接。或者用 QGIS 的“捕捉几何”工具把面边界顶点捕捉到河网线上。我一般会在融合后重新用polygonize拆分确保边界共享。4.3 坑三属性表编码乱码中文河名变问号现象.dbf文件里的中文字段在 QGIS 里显示正常但用geopandas读出来全是乱码或问号。原因.dbf的默认编码是GBK或Latin-1而geopandas默认按UTF-8读。解决读入时指定编码gdf gpd.read_file(yajiang_rivers.shp, encodinggbk)如果gbk不行就试latin-1再不行用chardet检测。我通常会在读入后把字段名和值都打印一遍确认没有乱码再继续。4.4 坑四线图层毛刺太多平滑后拓扑又断了现象从 DEM 提取的河网线锯齿严重做缓冲区分析时边界不平滑但用平滑算法处理后线又断成几截。原因平滑算法如simplify会改变线的几何形状如果容差设得太大相邻线段可能不再相交。解决先用小容差simplify(0.0001)做轻微平滑再用linemerge把断开的线合并。如果还不行就回到原始线用snap把端点捕捉到最近的其他线上。我一般会保留一份原始线作为备份平滑后的线只用于可视化分析还是用原始数据。4.5 坑五下载的 shp 缺少.shxQGIS 能开但代码读不了现象QGIS 里能正常加载 shp但geopandas.read_file()报错Unable to open ... .shx。原因.shx是索引文件QGIS 有时能容错但GDAL库严格要求它存在。解决如果压缩包里没有.shx用ogr2ogr重新生成ogr2ogr -f ESRI Shapefile output.shp input.shp这条命令会重新创建完整的 shp 文件集包括.shx。注意output.shp和input.shp不能同名否则会覆盖。5. 进阶技巧用雅江 shp 数据做子流域批量统计与自动化出图5.1 批量统计每个子流域的土地利用面积占比拿到雅江子流域 shp 后最常见的进阶需求是每个子流域里草地、林地、水体各占多少。手动一个个算太慢用geopandas的overlay加groupby可以批量出结果。import geopandas as gpd import pandas as pd # 读取子流域和土地利用数据 sub gpd.read_file(yajiang_sub_basins.shp).to_crs(epsg32645) landuse gpd.read_file(landuse_2020.shp).to_crs(epsg32645) # 叠加分析把土地利用切到每个子流域里 overlay gpd.overlay(landuse, sub, howintersection) # 计算每个子流域内各地类的面积 overlay[area_km2] overlay.area / 1e6 stats overlay.groupby([SUB_ID, LU_TYPE])[area_km2].sum().reset_index() # 透视成宽表每行一个子流域每列一个地类 pivot stats.pivot(indexSUB_ID, columnsLU_TYPE, valuesarea_km2).fillna(0) # 算占比 pivot_pct pivot.div(pivot.sum(axis1), axis0) * 100 pivot_pct.to_csv(sub_basin_landuse_pct.csv)gpd.overlay的howintersection保留两个图层的交集部分土地利用面被切碎到每个子流域里。groupby按子流域 ID 和地类汇总面积pivot把长表转宽表最后div算占比。输出的 CSV 可以直接拿去画堆叠柱状图。注意to_crs必须在overlay之前做否则面积计算还是错的。5.2 自动化出图用 matplotlib 批量渲染子流域专题图如果要把每个子流域单独出一张图用循环加matplotlib可以省掉大量重复劳动。import matplotlib.pyplot as plt fig, axes plt.subplots(3, 4, figsize(16, 12)) axes axes.flatten() for idx, (sub_id, group) in enumerate(sub.groupby(SUB_ID)): ax axes[idx] group.plot(axax, colorlightblue, edgecolorblack) ax.set_title(fSub-basin {sub_id}) ax.set_axis_off() plt.tight_layout() plt.savefig(sub_basins_grid.png, dpi150)这段代码把子流域按SUB_ID分组每组画在一个子图里。3x4的网格适合 12 个子流域多了就调整行列数。set_axis_off()去掉坐标轴出图更干净。保存时dpi150兼顾清晰度和文件大小。如果子流域数量不固定可以用math.ceil动态算行列。5.3 一个我常用的检查习惯每次拿到新的雅江 shp 数据我都会先跑一遍“三查”查坐标系、查几何有效性、查属性表字段。这三步花不了五分钟但能避免后面几个小时的返工。具体就是def quick_check(path): gdf gpd.read_file(path) print(fCRS: {gdf.crs}) print(f无效几何: {len(gdf[~gdf.is_valid])}) print(f字段: {list(gdf.columns)}) print(f要素数: {len(gdf)}) quick_check(yajiang_basin.shp)这个函数打印坐标系、无效几何数量、字段列表和要素总数。坐标系不对就投影几何无效就修复字段缺失就补算。养成这个习惯后我很少再因为数据本身的问题卡住。希望帮到你。本文还有配套的精品资源点击获取