尧图网页设计 
简介这份雅江流域矢量数据集面向地理信息、水文生态与遥感测绘方向的研究人员和学生提供雅鲁藏布江流域的空间边界与要素表达可用于流域管理、洪水风险评估、水资源配置及生态保育规划等空间分析与制图场景。资源包共8个文件压缩后约1.31MB包含shp主文件、shx索引、dbf属性表、prj投影定义以及sbn、sbx空间索引和xml元数据、png预览图完整覆盖GIS数据读取、坐标定义与属性查询所需的核心组成。雅鲁藏布江是世界海拔最高的大河之一流域面积约93.5万平方公里流经日喀则、拉萨、山南与林芝孕育了拉萨河、年楚河、多雄藏布等支流和雅鲁藏布大峡谷数据据此整理具备明确的地理与科研价值。目前已有332人学习下载适合需要快速获取流域矢量底图、开展空间建模或作为论文实证基础的用户参考使用。1. 雅江流域矢量数据集从 shp 文件到能用的边界数据做雅鲁藏布江流域的水文、生态或者地质灾害分析第一步往往不是跑模型而是找边界。雅江流域矢量数据集 shp 数据说白了就是一套把雅鲁藏布江干流、主要支流、流域分水界线用点线面矢量要素描述清楚的空间数据常见格式是 Esri Shapefile配套 .shp、.shx、.dbf、.prj 四件套。它解决的核心问题是你手上有一堆站点、栅格、遥感影像但没有一个统一的空间范围去裁剪、统计和叠加。适合谁做流域水文模拟的、算冰川融水补给量的、做土地利用变化统计的以及需要把雅江流域作为研究区写进论文的人。很多人卡在“数据下载”这一步下到的要么是栅格边界要么坐标系对不上要么属性表里连河名都没有。这篇就按我实际处理雅江数据的顺序把选型、下载后的整理、坐标系处理、裁剪统计和踩过的坑讲清楚让你拿到 shp 之后能直接进 QGIS 或 Python 里干活。2. 雅江流域 shp 数据拿到手先做什么图层结构与坐标系核对2.1 雅江矢量数据集通常包含哪些图层一套完整的雅江流域矢量数据按我经手的常见组织方式至少会拆成这么几类流域边界面watershed boundary、干流中心线river centerline、支流线tributaries、湖泊水面lakes、可能还有行政县界叠加层。为什么强调“图层结构”因为很多下载站把整个流域打包成一个 shp里面只有面没有河网你后面想做河网密度或者缓冲区分析就还得再找线数据。拿到数据第一件事不是急着打开看而是先列清单。图层名常见几何类型关键字段用途basin_boundaryPolygonNAME, AREA裁剪、统计范围main_riverLineStringNAME, ORDER干流长度、缓冲区tributariesLineStringNAME, ORDER河网密度lakesPolygonNAME, AREA水面变化admin_overlayPolygonADCODE, NAME分区统计这张表不是让你照抄字段名而是提醒你下载下来的数据字段名可能是拼音、英文缩写甚至乱码先看清楚再写代码。如果只有面没有线那这套数据只能做范围裁剪做不了河网分析得补数据。2.2 坐标系不核对后面全白干雅江流域横跨多个经度带常见坐标系有 WGS84 地理坐标EPSG:4326和 CGCS2000 投影坐标。我见过最典型的翻车现场下载的 shp 是地理坐标单位是度你直接拿它去算面积结果出来一堆小数还以为是数据错了。其实是因为没投影。核对方法很简单在 QGIS 里右键图层属性看 CRS或者用 Python 读 .prj 文件。下面这段代码就是干这个的import geopandas as gpd # 读取雅江流域边界 shp路径按你本地实际改 gdf gpd.read_file(./yajiang/basin_boundary.shp) # 打印坐标系信息重点看 EPSG 编号 print(CRS:, gdf.crs) print(几何类型:, gdf.geom_type.unique()) print(字段列表:, list(gdf.columns)) print(要素数量:, len(gdf)) # 如果是地理坐标单位度面积计算前必须投影 if gdf.crs and gdf.crs.is_geographic: print(当前是地理坐标系面积计算需先投影)逻辑说明gpd.read_file直接读 shp不需要额外装驱动gdf.crs返回坐标系对象is_geographic判断是不是经纬度。参数上路径一定要用你解压后的真实路径Windows 下注意反斜杠转义或者用正斜杠。如果打印出来 CRS 是 None说明 .prj 文件丢了这时候别硬算先去补 .prj 或者用已知坐标系重新定义。提示shp 四件套缺一不可尤其是 .prj。只拷贝 .shp 给别人等于把坐标系信息扔了。2.3 属性表乱码与字段类型检查dbf 文件里的中文属性在 QGIS 里可能显示正常在 Python 里读出来是乱码。原因是编码。常见做法是读的时候指定 encoding# 尝试用 gbk 读取属性解决中文乱码 gdf gpd.read_file(./yajiang/basin_boundary.shp, encodinggbk) print(gdf.head())如果 gbk 还不行就换 utf-8 或者 latin1 试。字段类型也要看AREA 字段如果是字符串你后面求和会报错得先astype(float)。这一步花五分钟能省后面两小时排查。3. 雅鲁藏布江数据下载后的整理从压缩包到可用的 GeoPackage3.1 下载来源与文件组织习惯雅鲁藏布江数据下载的渠道不少有科研机构共享的、有地理信息平台提供的也有从论文附件里扒出来的。我不建议你把下载的压缩包直接扔进项目文件夹而是按“原始数据 / 处理中间件 / 成果”三层来放。原始数据解压后先别改复制一份到工作目录再动手。常见做法是# 建立项目目录结构 mkdir -p yajiang_project/{raw,work,output} # 解压下载的 zip 到 raw unzip yajiang_shp.zip -d yajiang_project/raw/ # 查看解压后文件 ls -lh yajiang_project/raw/这样做的原因是shp 一旦被覆盖或者字段被改原始数据就没了后悔药没处买。raw 目录只读work 目录随便折腾。3.2 用 GeoPandas 合并多图层并转 GeoPackageshp 有个烦人的限制一个文件只能存一种几何类型字段名最长 10 个字符。雅江数据如果分了好几个 shp管理起来很碎。我一般会转成 GeoPackage.gpkg一个文件装多个图层字段名也不受限制。import geopandas as gpd from pathlib import Path raw_dir Path(./yajiang_project/raw) out_gpkg ./yajiang_project/work/yajiang.gpkg # 定义要合并的图层和对应文件名 layers { basin: basin_boundary.shp, main_river: main_river.shp, tributaries: tributaries.shp, } for layer_name, shp_file in layers.items(): shp_path raw_dir / shp_file if not shp_path.exists(): print(f跳过不存在的文件: {shp_path}) continue gdf gpd.read_file(shp_path, encodinggbk) # 统一写入 gpkgmodea 表示追加图层 gdf.to_file(out_gpkg, layerlayer_name, driverGPKG) print(f已写入图层 {layer_name}要素数 {len(gdf)})逻辑说明to_file指定 driverGPKG 就能写 GeoPackagelayer参数决定图层名。参数上如果 gpkg 已存在追加模式不会覆盖其他图层但同名图层会被替换。这一步做完你就有了一份干净的、多图层的雅江矢量库。3.3 几何有效性修复下载的 shp 经常有自相交、空几何、重复节点。这些在裁剪时会导致报错或者结果缺失。用 Shapely 检查并修复from shapely.validation import make_valid # 检查无效几何 invalid gdf[~gdf.is_valid] print(f无效几何数量: {len(invalid)}) # 修复无效几何 gdf[geometry] gdf[geometry].apply( lambda geom: make_valid(geom) if geom and not geom.is_valid else geom )make_valid是 Shapely 提供的修复函数比 buffer(0) 更稳。修复后最好再跑一次is_valid确认。如果无效几何特别多说明数据源质量差得考虑换一套。4. 雅江流域矢量数据的裁剪、统计与出图避坑4.1 用流域边界裁剪栅格与矢量拿到流域边界后最常见的操作是裁剪。矢量裁矢量用gpd.clip裁栅格用rasterio.mask。先讲矢量裁剪# 读取流域边界和待裁剪的支流 basin gpd.read_file(./yajiang_project/work/yajiang.gpkg, layerbasin) trib gpd.read_file(./yajiang_project/work/yajiang.gpkg, layertributaries) # 确保坐标系一致 trib trib.to_crs(basin.crs) # 裁剪 trib_clipped gpd.clip(trib, basin) print(f裁剪前 {len(trib)}裁剪后 {len(trib_clipped)}) trib_clipped.to_file(./yajiang_project/output/trib_clipped.gpkg, layertrib_clipped, driverGPKG)关键参数是to_crs如果两个图层坐标系不一致clip 会报错或者结果错位。裁剪前一定先统一 CRS。裁栅格的话用 basin 的几何做 maskcropTrue表示只保留边界内部分。4.2 面积与长度统计算法选择面积统计有个坑地理坐标系下算出来的“度平方”没有意义。正确做法是先投影到等面积投影比如 Albers 或者 UTM 对应带。雅江流域跨度大我一般用 Albers 自定义参数# 定义 Albers 等面积投影适合雅江流域 albers projaea lat_127 lat_233 lat_030 lon_088 datumWGS84 unitsm basin_proj basin.to_crs(albers) basin_proj[area_km2] basin_proj.area / 1e6 print(basin_proj[[NAME, area_km2]])lat_1、lat_2是标准纬线取流域南北范围附近lon_0是中央经线。这样算出来的面积单位是平方米除以 1e6 得平方公里。长度统计同理投影后.length才是米。4.3 出图时的标注与边界压盖用 QGIS 出图雅江流域边界线宽设 0.5 到 0.8 比较合适太粗会压盖内部河网。标注河流名称时如果线要素太长标注会重复可以在图层属性里设置“标注位置”为沿线并限制重复距离。导出 PDF 时勾选“裁剪到地图范围”避免白边。这些细节不影响分析但影响你给合作者看的观感。5. 雅江 shp 数据处理常见问题排查坐标系、编码与几何5.1 现象裁剪结果为空或只剩一小块原因两个图层坐标系不一致或者边界几何无效。解决先print(gdf.crs)对比再make_valid修复边界最后统一to_crs再裁剪。我遇到过边界是地理坐标、支流是投影坐标clip 出来直接空排查了半天。5.2 现象属性表中文变问号原因dbf 编码不是 utf-8。解决读的时候试encodinggbk如果还乱用 QGIS 打开另存为 utf-8 编码的 gpkg。别在 Excel 里直接改 dbf会破坏字段结构。5.3 现象面积算出来是 0.000x原因地理坐标系下直接算面积。解决投影到等面积坐标系再算参考 4.2 的 Albers 参数。这是新手最容易踩的坑没有之一。5.4 现象shp 文件打开报“缺少 .prj”原因下载或拷贝时丢了 .prj。解决如果知道数据来源坐标系用gdf.set_crs(EPSG:4326)补上再转出不知道就联系数据提供方别瞎猜。5.5 现象GeoPackage 写入后图层名重复原因多次运行脚本同名图层被覆盖或追加混乱。解决写入前判断图层是否存在或者每次输出带时间戳的文件名。我一般会在脚本里加if out_gpkg.exists(): out_gpkg.unlink()做全量重建避免脏数据累积。6. 让雅江矢量数据可复现一个检查脚本和我的固定习惯最后一章不讲大道理给一个我每次拿到新 shp 都会跑的检查脚本以及几个固定习惯。这个脚本把坐标系、几何有效性、字段类型、范围框一次性打出来三十秒判断这套雅江数据能不能直接用。import geopandas as gpd def inspect_shp(path, encodinggbk): gdf gpd.read_file(path, encodingencoding) print( * 40) print(f文件: {path}) print(fCRS: {gdf.crs}) print(f几何类型: {gdf.geom_type.unique()}) print(f要素数: {len(gdf)}) print(f字段: {list(gdf.columns)}) print(f无效几何: {(~gdf.is_valid).sum()}) print(f范围: {gdf.total_bounds}) # 检查是否有空几何 empty gdf[gdf.is_empty] print(f空几何: {len(empty)}) return gdf # 用法 gdf inspect_shp(./yajiang_project/raw/basin_boundary.shp)total_bounds返回 [minx, miny, maxx, maxy]你可以快速判断经度是不是在 80 到 100 之间、纬度是不是在 27 到 33 之间如果差太远说明坐标系或者数据范围有问题。is_empty检查空几何空几何在后续操作里会静默失败必须提前清掉。我的固定习惯有这么几条第一raw 目录永不写入所有修改在 work 目录做第二每次处理前先跑一遍 inspect把 CRS 和 bounds 记在脚本注释里第三输出统一用 GeoPackage不用 shp省得四件套丢文件第四面积长度统计前先投影投影参数写在脚本顶部常量区方便换流域时改。这些习惯不高级但能让你在换一台机器或者隔三个月再回头看时五分钟内复现整个流程。雅江流域矢量数据集 shp 数据这件事难点从来不是下载而是下载之后坐标系、几何、编码这三关。把这三关过了后面裁剪统计就是常规操作。希望帮到你。本文还有配套的精品资源点击获取 雅江流域矢量数据集shp处理全流程：坐标系、几何修复与裁剪统计
简介这份雅江流域矢量数据集面向地理信息、水文生态与遥感测绘方向的研究人员和学生提供雅鲁藏布江流域的空间边界与要素表达可用于流域管理、洪水风险评估、水资源配置及生态保育规划等空间分析与制图场景。资源包共8个文件压缩后约1.31MB包含shp主文件、shx索引、dbf属性表、prj投影定义以及sbn、sbx空间索引和xml元数据、png预览图完整覆盖GIS数据读取、坐标定义与属性查询所需的核心组成。雅鲁藏布江是世界海拔最高的大河之一流域面积约93.5万平方公里流经日喀则、拉萨、山南与林芝孕育了拉萨河、年楚河、多雄藏布等支流和雅鲁藏布大峡谷数据据此整理具备明确的地理与科研价值。目前已有332人学习下载适合需要快速获取流域矢量底图、开展空间建模或作为论文实证基础的用户参考使用。1. 雅江流域矢量数据集从 shp 文件到能用的边界数据做雅鲁藏布江流域的水文、生态或者地质灾害分析第一步往往不是跑模型而是找边界。雅江流域矢量数据集 shp 数据说白了就是一套把雅鲁藏布江干流、主要支流、流域分水界线用点线面矢量要素描述清楚的空间数据常见格式是 Esri Shapefile配套 .shp、.shx、.dbf、.prj 四件套。它解决的核心问题是你手上有一堆站点、栅格、遥感影像但没有一个统一的空间范围去裁剪、统计和叠加。适合谁做流域水文模拟的、算冰川融水补给量的、做土地利用变化统计的以及需要把雅江流域作为研究区写进论文的人。很多人卡在“数据下载”这一步下到的要么是栅格边界要么坐标系对不上要么属性表里连河名都没有。这篇就按我实际处理雅江数据的顺序把选型、下载后的整理、坐标系处理、裁剪统计和踩过的坑讲清楚让你拿到 shp 之后能直接进 QGIS 或 Python 里干活。2. 雅江流域 shp 数据拿到手先做什么图层结构与坐标系核对2.1 雅江矢量数据集通常包含哪些图层一套完整的雅江流域矢量数据按我经手的常见组织方式至少会拆成这么几类流域边界面watershed boundary、干流中心线river centerline、支流线tributaries、湖泊水面lakes、可能还有行政县界叠加层。为什么强调“图层结构”因为很多下载站把整个流域打包成一个 shp里面只有面没有河网你后面想做河网密度或者缓冲区分析就还得再找线数据。拿到数据第一件事不是急着打开看而是先列清单。图层名常见几何类型关键字段用途basin_boundaryPolygonNAME, AREA裁剪、统计范围main_riverLineStringNAME, ORDER干流长度、缓冲区tributariesLineStringNAME, ORDER河网密度lakesPolygonNAME, AREA水面变化admin_overlayPolygonADCODE, NAME分区统计这张表不是让你照抄字段名而是提醒你下载下来的数据字段名可能是拼音、英文缩写甚至乱码先看清楚再写代码。如果只有面没有线那这套数据只能做范围裁剪做不了河网分析得补数据。2.2 坐标系不核对后面全白干雅江流域横跨多个经度带常见坐标系有 WGS84 地理坐标EPSG:4326和 CGCS2000 投影坐标。我见过最典型的翻车现场下载的 shp 是地理坐标单位是度你直接拿它去算面积结果出来一堆小数还以为是数据错了。其实是因为没投影。核对方法很简单在 QGIS 里右键图层属性看 CRS或者用 Python 读 .prj 文件。下面这段代码就是干这个的import geopandas as gpd # 读取雅江流域边界 shp路径按你本地实际改 gdf gpd.read_file(./yajiang/basin_boundary.shp) # 打印坐标系信息重点看 EPSG 编号 print(CRS:, gdf.crs) print(几何类型:, gdf.geom_type.unique()) print(字段列表:, list(gdf.columns)) print(要素数量:, len(gdf)) # 如果是地理坐标单位度面积计算前必须投影 if gdf.crs and gdf.crs.is_geographic: print(当前是地理坐标系面积计算需先投影)逻辑说明gpd.read_file直接读 shp不需要额外装驱动gdf.crs返回坐标系对象is_geographic判断是不是经纬度。参数上路径一定要用你解压后的真实路径Windows 下注意反斜杠转义或者用正斜杠。如果打印出来 CRS 是 None说明 .prj 文件丢了这时候别硬算先去补 .prj 或者用已知坐标系重新定义。提示shp 四件套缺一不可尤其是 .prj。只拷贝 .shp 给别人等于把坐标系信息扔了。2.3 属性表乱码与字段类型检查dbf 文件里的中文属性在 QGIS 里可能显示正常在 Python 里读出来是乱码。原因是编码。常见做法是读的时候指定 encoding# 尝试用 gbk 读取属性解决中文乱码 gdf gpd.read_file(./yajiang/basin_boundary.shp, encodinggbk) print(gdf.head())如果 gbk 还不行就换 utf-8 或者 latin1 试。字段类型也要看AREA 字段如果是字符串你后面求和会报错得先astype(float)。这一步花五分钟能省后面两小时排查。3. 雅鲁藏布江数据下载后的整理从压缩包到可用的 GeoPackage3.1 下载来源与文件组织习惯雅鲁藏布江数据下载的渠道不少有科研机构共享的、有地理信息平台提供的也有从论文附件里扒出来的。我不建议你把下载的压缩包直接扔进项目文件夹而是按“原始数据 / 处理中间件 / 成果”三层来放。原始数据解压后先别改复制一份到工作目录再动手。常见做法是# 建立项目目录结构 mkdir -p yajiang_project/{raw,work,output} # 解压下载的 zip 到 raw unzip yajiang_shp.zip -d yajiang_project/raw/ # 查看解压后文件 ls -lh yajiang_project/raw/这样做的原因是shp 一旦被覆盖或者字段被改原始数据就没了后悔药没处买。raw 目录只读work 目录随便折腾。3.2 用 GeoPandas 合并多图层并转 GeoPackageshp 有个烦人的限制一个文件只能存一种几何类型字段名最长 10 个字符。雅江数据如果分了好几个 shp管理起来很碎。我一般会转成 GeoPackage.gpkg一个文件装多个图层字段名也不受限制。import geopandas as gpd from pathlib import Path raw_dir Path(./yajiang_project/raw) out_gpkg ./yajiang_project/work/yajiang.gpkg # 定义要合并的图层和对应文件名 layers { basin: basin_boundary.shp, main_river: main_river.shp, tributaries: tributaries.shp, } for layer_name, shp_file in layers.items(): shp_path raw_dir / shp_file if not shp_path.exists(): print(f跳过不存在的文件: {shp_path}) continue gdf gpd.read_file(shp_path, encodinggbk) # 统一写入 gpkgmodea 表示追加图层 gdf.to_file(out_gpkg, layerlayer_name, driverGPKG) print(f已写入图层 {layer_name}要素数 {len(gdf)})逻辑说明to_file指定 driverGPKG 就能写 GeoPackagelayer参数决定图层名。参数上如果 gpkg 已存在追加模式不会覆盖其他图层但同名图层会被替换。这一步做完你就有了一份干净的、多图层的雅江矢量库。3.3 几何有效性修复下载的 shp 经常有自相交、空几何、重复节点。这些在裁剪时会导致报错或者结果缺失。用 Shapely 检查并修复from shapely.validation import make_valid # 检查无效几何 invalid gdf[~gdf.is_valid] print(f无效几何数量: {len(invalid)}) # 修复无效几何 gdf[geometry] gdf[geometry].apply( lambda geom: make_valid(geom) if geom and not geom.is_valid else geom )make_valid是 Shapely 提供的修复函数比 buffer(0) 更稳。修复后最好再跑一次is_valid确认。如果无效几何特别多说明数据源质量差得考虑换一套。4. 雅江流域矢量数据的裁剪、统计与出图避坑4.1 用流域边界裁剪栅格与矢量拿到流域边界后最常见的操作是裁剪。矢量裁矢量用gpd.clip裁栅格用rasterio.mask。先讲矢量裁剪# 读取流域边界和待裁剪的支流 basin gpd.read_file(./yajiang_project/work/yajiang.gpkg, layerbasin) trib gpd.read_file(./yajiang_project/work/yajiang.gpkg, layertributaries) # 确保坐标系一致 trib trib.to_crs(basin.crs) # 裁剪 trib_clipped gpd.clip(trib, basin) print(f裁剪前 {len(trib)}裁剪后 {len(trib_clipped)}) trib_clipped.to_file(./yajiang_project/output/trib_clipped.gpkg, layertrib_clipped, driverGPKG)关键参数是to_crs如果两个图层坐标系不一致clip 会报错或者结果错位。裁剪前一定先统一 CRS。裁栅格的话用 basin 的几何做 maskcropTrue表示只保留边界内部分。4.2 面积与长度统计算法选择面积统计有个坑地理坐标系下算出来的“度平方”没有意义。正确做法是先投影到等面积投影比如 Albers 或者 UTM 对应带。雅江流域跨度大我一般用 Albers 自定义参数# 定义 Albers 等面积投影适合雅江流域 albers projaea lat_127 lat_233 lat_030 lon_088 datumWGS84 unitsm basin_proj basin.to_crs(albers) basin_proj[area_km2] basin_proj.area / 1e6 print(basin_proj[[NAME, area_km2]])lat_1、lat_2是标准纬线取流域南北范围附近lon_0是中央经线。这样算出来的面积单位是平方米除以 1e6 得平方公里。长度统计同理投影后.length才是米。4.3 出图时的标注与边界压盖用 QGIS 出图雅江流域边界线宽设 0.5 到 0.8 比较合适太粗会压盖内部河网。标注河流名称时如果线要素太长标注会重复可以在图层属性里设置“标注位置”为沿线并限制重复距离。导出 PDF 时勾选“裁剪到地图范围”避免白边。这些细节不影响分析但影响你给合作者看的观感。5. 雅江 shp 数据处理常见问题排查坐标系、编码与几何5.1 现象裁剪结果为空或只剩一小块原因两个图层坐标系不一致或者边界几何无效。解决先print(gdf.crs)对比再make_valid修复边界最后统一to_crs再裁剪。我遇到过边界是地理坐标、支流是投影坐标clip 出来直接空排查了半天。5.2 现象属性表中文变问号原因dbf 编码不是 utf-8。解决读的时候试encodinggbk如果还乱用 QGIS 打开另存为 utf-8 编码的 gpkg。别在 Excel 里直接改 dbf会破坏字段结构。5.3 现象面积算出来是 0.000x原因地理坐标系下直接算面积。解决投影到等面积坐标系再算参考 4.2 的 Albers 参数。这是新手最容易踩的坑没有之一。5.4 现象shp 文件打开报“缺少 .prj”原因下载或拷贝时丢了 .prj。解决如果知道数据来源坐标系用gdf.set_crs(EPSG:4326)补上再转出不知道就联系数据提供方别瞎猜。5.5 现象GeoPackage 写入后图层名重复原因多次运行脚本同名图层被覆盖或追加混乱。解决写入前判断图层是否存在或者每次输出带时间戳的文件名。我一般会在脚本里加if out_gpkg.exists(): out_gpkg.unlink()做全量重建避免脏数据累积。6. 让雅江矢量数据可复现一个检查脚本和我的固定习惯最后一章不讲大道理给一个我每次拿到新 shp 都会跑的检查脚本以及几个固定习惯。这个脚本把坐标系、几何有效性、字段类型、范围框一次性打出来三十秒判断这套雅江数据能不能直接用。import geopandas as gpd def inspect_shp(path, encodinggbk): gdf gpd.read_file(path, encodingencoding) print( * 40) print(f文件: {path}) print(fCRS: {gdf.crs}) print(f几何类型: {gdf.geom_type.unique()}) print(f要素数: {len(gdf)}) print(f字段: {list(gdf.columns)}) print(f无效几何: {(~gdf.is_valid).sum()}) print(f范围: {gdf.total_bounds}) # 检查是否有空几何 empty gdf[gdf.is_empty] print(f空几何: {len(empty)}) return gdf # 用法 gdf inspect_shp(./yajiang_project/raw/basin_boundary.shp)total_bounds返回 [minx, miny, maxx, maxy]你可以快速判断经度是不是在 80 到 100 之间、纬度是不是在 27 到 33 之间如果差太远说明坐标系或者数据范围有问题。is_empty检查空几何空几何在后续操作里会静默失败必须提前清掉。我的固定习惯有这么几条第一raw 目录永不写入所有修改在 work 目录做第二每次处理前先跑一遍 inspect把 CRS 和 bounds 记在脚本注释里第三输出统一用 GeoPackage不用 shp省得四件套丢文件第四面积长度统计前先投影投影参数写在脚本顶部常量区方便换流域时改。这些习惯不高级但能让你在换一台机器或者隔三个月再回头看时五分钟内复现整个流程。雅江流域矢量数据集 shp 数据这件事难点从来不是下载而是下载之后坐标系、几何、编码这三关。把这三关过了后面裁剪统计就是常规操作。希望帮到你。本文还有配套的精品资源点击获取