尧图网站建设 
简介这份资源面向GIS初学者、水文地理研究者及高校师生提供一套可直接用于空间分析与制图练习的河流水系矢量数据。压缩包共46个文件约9.7MB以SHP格式为核心包含shp、shx、dbf、prj等配套文件分别承载几何图形、索引、属性表与坐标投影信息另有sbn、sbx空间索引及xml元数据、jpg预览图方便在ArcGIS、QGIS等软件中直接加载编辑。数据按河流等级与流域层级组织涵盖一级至五级河流及三级流域、三级以上河流等图层可用于流域水资源分布研究、洪水风险评估、河流生态保护规划等场景。目前已有138人学习下载读者可借助这些数据练习距离量测、缓冲区分析与空间查询并叠加人口、地形、气候等图层开展综合决策分析也可用于课堂教学帮助理解河流分级体系与流域管理逻辑。1. 拿到「河流及流域.zip」先别急着解压这份 SHP 学习资料到底装了什么如果你正在做流域边界提取、河网拓扑检查或者要给水文站做空间关联大概率会在某个资料包里翻到一个叫「河流及流域.zip」的东西。它不是一个可执行程序也不是某个软件的安装包而是一套以 SHPShapefile为核心的矢量地理数据学习资料。SHP 是 Esri 定义的一种经典矢量格式虽然老但在 QGIS、ArcGIS、GeoPandas、PostGIS 里都能直接读兼容性极好。这份资料的价值在于它把河流线要素和流域面要素放在同一个包里省去了你从不同来源分别找水系和汇水区、再手动对齐坐标系的麻烦。适合谁做水文分析的学生、需要快速验证流域划分算法的工程师、以及想拿真实矢量数据练手空间连接和叠加分析的人。但先别双击解压——SHP 不是单个文件解压后你会看到一堆同名不同后缀的文件漏掉一个就可能打不开。2. 拆开 SHP 的「黑匣子」从文件组成到坐标系判读2.1 为什么一个 SHP 会解压出五六个文件Shapefile 这个叫法其实有误导性它本质是一个文件集合。解压「河流及流域.zip」后你通常会看到.shp、.shx、.dbf、.prj有时还有.cpg和.sbn。.shp存几何形状.shx是索引.dbf是属性表.prj存坐标系定义.cpg指定字符编码。很多人只复制.shp去用结果属性表丢失或者中文乱码这就是血泪经验的来源。我一般会先看.prj里有没有GCS_WGS_1984或CGCS2000字样这决定了后面做面积量算时要不要投影。如果.prj缺失QGIS 会提示未知坐标系此时不要硬猜用ogrinfo先看一眼元数据。# 查看 SHP 的图层名、要素数量和坐标系 ogrinfo -so 河流.shp 河流 # 输出示例Layer name: 河流 Geometry: Line String Feature Count: 1247 # 以及 Extent 和 Layer SRS 信息-so表示 summary only只输出摘要不打印每个要素。河流是图层名通常和文件名一致。如果这里报Unable to open datasource先检查路径里有没有中文或空格GDAL 对中文路径的支持取决于编译选项稳妥做法是放到纯英文目录下再试。2.2 用 QGIS 快速判读河流与流域的拓扑关系打开 QGIS把河流线和流域面一起拖进去。先看流域面是不是闭合多边形再看河流线是否落在流域内部。常见问题是河流线超出了流域边界这在水文上意味着流域划分和河网提取用了不同精度的 DEM。你可以用「按位置选择」工具把河流线作为输入流域面作为选择条件选那些与流域相交的河流。如果数量对不上说明数据本身有拓扑错误。另一个检查点是属性表里的字段河流图层通常有NAME、LENGTH或ORDER流域图层可能有AREA、PERIMETER或BASIN_ID。这些字段名不是固定的但你可以通过字段类型判断——数值型字段适合做统计文本型适合做分类。import geopandas as gpd # 读取河流和流域注意指定编码避免中文乱码 rivers gpd.read_file(河流.shp, encodingutf-8) basins gpd.read_file(流域.shp, encodingutf-8) print(rivers.crs) # 查看坐标系 print(basins.columns) # 查看字段名 print(rivers.geom_type) # 应为 LineString 或 MultiLineString print(basins.geom_type) # 应为 Polygon 或 MultiPolygon # 检查河流是否全部落在流域范围内 mask rivers.within(basins.unary_union) print(f落在流域内的河流数: {mask.sum()} / {len(rivers)})encodingutf-8是为了让.dbf里的中文属性正确显示如果报错可以换成gbk。unary_union把多个流域面合并成一个整体几何within返回布尔序列。如果落在流域内的河流数明显少于总数说明有河流线在流域外需要检查数据来源是否一致。3. 把 SHP 用起来从坐标转换到空间连接的完整操作链3.1 坐标系不统一时先做投影别直接算面积很多新手拿到 SHP 后直接算流域面积结果单位是「平方度」数值小得离谱。这是因为地理坐标系经纬度下算面积没有物理意义。正确做法是先投影到合适的投影坐标系。中国常用 CGCS2000 高斯克吕格投影按经度带选择 EPSG 代码比如 EPSG:4547 对应 CGCS2000 / 3-degree Gauss-Kruger CM 114E。如果你不确定用哪个带可以用pyproj查。from pyproj import CRS # 查看当前坐标系 print(rivers.crs) # 如果输出 EPSG:4326说明是 WGS84 经纬度 # 投影到 CGCS2000 3度带中央经线114E rivers_proj rivers.to_crs(epsg4547) basins_proj basins.to_crs(epsg4547) # 现在算面积单位是平方米 basins_proj[area_km2] basins_proj.geometry.area / 1e6 print(basins_proj[[BASIN_ID, area_km2]].head())to_crs不改变几何形状只做坐标换算。epsg4547是投影坐标系的编号不同地区要换。area / 1e6把平方米转成平方公里。如果投影后面积仍然异常检查几何是否有效用basins_proj.is_valid看有没有自相交多边形。3.2 用空间连接把河流属性和流域属性挂到一起水文分析里经常需要知道每条河流属于哪个流域。GeoPandas 的sjoin可以按空间关系连接两个图层。默认是intersects即只要几何有交集就匹配。对于河流和流域更合理的是用within要求河流完全落在流域内。# 空间连接河流落在哪个流域内 joined gpd.sjoin(rivers_proj, basins_proj, howleft, predicatewithin) # 查看结果字段 print(joined.columns) # 现在每条河流都带上了流域的 BASIN_ID 和 area_km2 # 统计每个流域内的河流总长度 joined[length_km] joined.geometry.length / 1000 summary joined.groupby(BASIN_ID)[length_km].sum().reset_index() print(summary.head())howleft保留所有河流即使没匹配到流域也保留方便排查。predicatewithin比默认的intersects更严格。groupby按流域编号汇总河流长度这是计算河网密度的基础。如果结果里BASIN_ID出现空值说明有河流不在任何流域内需要回到 2.2 节检查拓扑。3.3 导出为 GeoJSON 或 PostGIS 供后续使用SHP 有 2GB 大小限制字段名也不能超过 10 个字符。如果数据量大或者字段名长建议转成 GeoJSON 或直接入库 PostGIS。GeoJSON 用文本存几何可读性好但文件体积大。PostGIS 适合做复杂查询和并发访问。# 用 ogr2ogr 转 GeoJSON ogr2ogr -f GeoJSON 河流.geojson 河流.shp -lco ENCODINGUTF-8 # 导入 PostGIS ogr2ogr -f PostgreSQL PG:dbnamehydrology userpostgres 流域.shp -nln basins -lco GEOMETRY_NAMEgeom-lco ENCODINGUTF-8确保中文属性不乱码。PostGIS 连接字符串里dbname、user按实际改。-nln basins指定表名。导入后可以用 SQL 做空间查询比如SELECT * FROM basins WHERE ST_Area(geom) 1000000;。4. 避坑与排查SHP 处理中最容易翻车的五个场景4.1 现象打开 SHP 后属性表全是乱码原因.dbf文件的字符编码和读取时指定的编码不一致。中文 Windows 下生成的 SHP 常用 GBK而 QGIS 或 GeoPandas 默认用 UTF-8。解决先看有没有.cpg文件里面会写编码。没有的话用gpd.read_file(..., encodinggbk)试或者用 QGIS 的「重新编码」功能转成 UTF-8。4.2 现象sjoin结果为空但两个图层明明有重叠原因坐标系不一致。一个图层是 EPSG:4326另一个是 EPSG:4547虽然视觉上重叠但数值范围差几个数量级。解决统一用to_crs转到同一坐标系再做空间操作。我一般会在sjoin前打印两个图层的crs和total_bounds确认范围接近。4.3 现象计算面积时结果为零或负数原因多边形几何无效比如自相交、环方向错误。解决用basins.geometry.is_valid检查对无效几何用basins.geometry.buffer(0)修复。buffer(0)有时能解决自相交但会轻微改变边界修复后要重新检查面积。4.4 现象ogrinfo报错Unable to open datasource原因路径含中文、空格或特殊字符或者.shp文件不完整缺少.shx或.dbf。解决把文件移到纯英文路径确保同名文件都在。如果是从压缩包直接打开先完整解压。4.5 现象导出 GeoJSON 后中文变成\uXXXX原因GeoJSON 规范默认用 Unicode 转义不是乱码。解决用ensure_asciiFalse参数或者在 QGIS 导出时勾选「不转义 Unicode」。如果是在 Web 端用浏览器能正确解析转义序列不影响显示。5. 进阶技巧用河流与流域数据做河网密度分级和验证河网密度是流域地貌的重要指标等于流域内河流总长度除以流域面积。用这份资料可以快速算出来并做分级渲染。先按 3.2 节把河流和流域连接再按BASIN_ID汇总长度最后除以面积。# 计算河网密度 density joined.groupby(BASIN_ID).agg( total_length_km(length_km, sum), area_km2(area_km2, first) ).reset_index() density[density] density[total_length_km] / density[area_km2] print(density.sort_values(density, ascendingFalse).head(10))agg里用first取面积因为同一流域内面积值相同。density单位是 km/km²数值越大说明河网越密集。你可以把结果导回 QGIS按density字段做渐变符号化直观看出哪些流域水系发育。验证方法拿一条已知长度的河流手动在 QGIS 里量一下和length_km对比。如果差异超过 5%检查投影是否合适。另外可以用rivers.geometry.length在投影坐标系下算不要在地理坐标系下算。# 验证单条河流长度 sample rivers_proj[rivers_proj[NAME] 长江] print(sample.geometry.length / 1000) # 应接近实际长度如果NAME字段不存在用FID或OBJECTID定位。长度单位是米除以 1000 得公里。这个习惯我每次拿到新 SHP 都会走一遍——先看坐标系再验一条已知要素确认无误才批量处理。从那以后我每次处理河流及流域数据都强制走一遍坐标系检查和单要素验证省去了很多返工。希望帮到你。本文还有配套的精品资源点击获取 河流及流域SHP数据实战：从坐标系判读到空间连接与河网密度计算
简介这份资源面向GIS初学者、水文地理研究者及高校师生提供一套可直接用于空间分析与制图练习的河流水系矢量数据。压缩包共46个文件约9.7MB以SHP格式为核心包含shp、shx、dbf、prj等配套文件分别承载几何图形、索引、属性表与坐标投影信息另有sbn、sbx空间索引及xml元数据、jpg预览图方便在ArcGIS、QGIS等软件中直接加载编辑。数据按河流等级与流域层级组织涵盖一级至五级河流及三级流域、三级以上河流等图层可用于流域水资源分布研究、洪水风险评估、河流生态保护规划等场景。目前已有138人学习下载读者可借助这些数据练习距离量测、缓冲区分析与空间查询并叠加人口、地形、气候等图层开展综合决策分析也可用于课堂教学帮助理解河流分级体系与流域管理逻辑。1. 拿到「河流及流域.zip」先别急着解压这份 SHP 学习资料到底装了什么如果你正在做流域边界提取、河网拓扑检查或者要给水文站做空间关联大概率会在某个资料包里翻到一个叫「河流及流域.zip」的东西。它不是一个可执行程序也不是某个软件的安装包而是一套以 SHPShapefile为核心的矢量地理数据学习资料。SHP 是 Esri 定义的一种经典矢量格式虽然老但在 QGIS、ArcGIS、GeoPandas、PostGIS 里都能直接读兼容性极好。这份资料的价值在于它把河流线要素和流域面要素放在同一个包里省去了你从不同来源分别找水系和汇水区、再手动对齐坐标系的麻烦。适合谁做水文分析的学生、需要快速验证流域划分算法的工程师、以及想拿真实矢量数据练手空间连接和叠加分析的人。但先别双击解压——SHP 不是单个文件解压后你会看到一堆同名不同后缀的文件漏掉一个就可能打不开。2. 拆开 SHP 的「黑匣子」从文件组成到坐标系判读2.1 为什么一个 SHP 会解压出五六个文件Shapefile 这个叫法其实有误导性它本质是一个文件集合。解压「河流及流域.zip」后你通常会看到.shp、.shx、.dbf、.prj有时还有.cpg和.sbn。.shp存几何形状.shx是索引.dbf是属性表.prj存坐标系定义.cpg指定字符编码。很多人只复制.shp去用结果属性表丢失或者中文乱码这就是血泪经验的来源。我一般会先看.prj里有没有GCS_WGS_1984或CGCS2000字样这决定了后面做面积量算时要不要投影。如果.prj缺失QGIS 会提示未知坐标系此时不要硬猜用ogrinfo先看一眼元数据。# 查看 SHP 的图层名、要素数量和坐标系 ogrinfo -so 河流.shp 河流 # 输出示例Layer name: 河流 Geometry: Line String Feature Count: 1247 # 以及 Extent 和 Layer SRS 信息-so表示 summary only只输出摘要不打印每个要素。河流是图层名通常和文件名一致。如果这里报Unable to open datasource先检查路径里有没有中文或空格GDAL 对中文路径的支持取决于编译选项稳妥做法是放到纯英文目录下再试。2.2 用 QGIS 快速判读河流与流域的拓扑关系打开 QGIS把河流线和流域面一起拖进去。先看流域面是不是闭合多边形再看河流线是否落在流域内部。常见问题是河流线超出了流域边界这在水文上意味着流域划分和河网提取用了不同精度的 DEM。你可以用「按位置选择」工具把河流线作为输入流域面作为选择条件选那些与流域相交的河流。如果数量对不上说明数据本身有拓扑错误。另一个检查点是属性表里的字段河流图层通常有NAME、LENGTH或ORDER流域图层可能有AREA、PERIMETER或BASIN_ID。这些字段名不是固定的但你可以通过字段类型判断——数值型字段适合做统计文本型适合做分类。import geopandas as gpd # 读取河流和流域注意指定编码避免中文乱码 rivers gpd.read_file(河流.shp, encodingutf-8) basins gpd.read_file(流域.shp, encodingutf-8) print(rivers.crs) # 查看坐标系 print(basins.columns) # 查看字段名 print(rivers.geom_type) # 应为 LineString 或 MultiLineString print(basins.geom_type) # 应为 Polygon 或 MultiPolygon # 检查河流是否全部落在流域范围内 mask rivers.within(basins.unary_union) print(f落在流域内的河流数: {mask.sum()} / {len(rivers)})encodingutf-8是为了让.dbf里的中文属性正确显示如果报错可以换成gbk。unary_union把多个流域面合并成一个整体几何within返回布尔序列。如果落在流域内的河流数明显少于总数说明有河流线在流域外需要检查数据来源是否一致。3. 把 SHP 用起来从坐标转换到空间连接的完整操作链3.1 坐标系不统一时先做投影别直接算面积很多新手拿到 SHP 后直接算流域面积结果单位是「平方度」数值小得离谱。这是因为地理坐标系经纬度下算面积没有物理意义。正确做法是先投影到合适的投影坐标系。中国常用 CGCS2000 高斯克吕格投影按经度带选择 EPSG 代码比如 EPSG:4547 对应 CGCS2000 / 3-degree Gauss-Kruger CM 114E。如果你不确定用哪个带可以用pyproj查。from pyproj import CRS # 查看当前坐标系 print(rivers.crs) # 如果输出 EPSG:4326说明是 WGS84 经纬度 # 投影到 CGCS2000 3度带中央经线114E rivers_proj rivers.to_crs(epsg4547) basins_proj basins.to_crs(epsg4547) # 现在算面积单位是平方米 basins_proj[area_km2] basins_proj.geometry.area / 1e6 print(basins_proj[[BASIN_ID, area_km2]].head())to_crs不改变几何形状只做坐标换算。epsg4547是投影坐标系的编号不同地区要换。area / 1e6把平方米转成平方公里。如果投影后面积仍然异常检查几何是否有效用basins_proj.is_valid看有没有自相交多边形。3.2 用空间连接把河流属性和流域属性挂到一起水文分析里经常需要知道每条河流属于哪个流域。GeoPandas 的sjoin可以按空间关系连接两个图层。默认是intersects即只要几何有交集就匹配。对于河流和流域更合理的是用within要求河流完全落在流域内。# 空间连接河流落在哪个流域内 joined gpd.sjoin(rivers_proj, basins_proj, howleft, predicatewithin) # 查看结果字段 print(joined.columns) # 现在每条河流都带上了流域的 BASIN_ID 和 area_km2 # 统计每个流域内的河流总长度 joined[length_km] joined.geometry.length / 1000 summary joined.groupby(BASIN_ID)[length_km].sum().reset_index() print(summary.head())howleft保留所有河流即使没匹配到流域也保留方便排查。predicatewithin比默认的intersects更严格。groupby按流域编号汇总河流长度这是计算河网密度的基础。如果结果里BASIN_ID出现空值说明有河流不在任何流域内需要回到 2.2 节检查拓扑。3.3 导出为 GeoJSON 或 PostGIS 供后续使用SHP 有 2GB 大小限制字段名也不能超过 10 个字符。如果数据量大或者字段名长建议转成 GeoJSON 或直接入库 PostGIS。GeoJSON 用文本存几何可读性好但文件体积大。PostGIS 适合做复杂查询和并发访问。# 用 ogr2ogr 转 GeoJSON ogr2ogr -f GeoJSON 河流.geojson 河流.shp -lco ENCODINGUTF-8 # 导入 PostGIS ogr2ogr -f PostgreSQL PG:dbnamehydrology userpostgres 流域.shp -nln basins -lco GEOMETRY_NAMEgeom-lco ENCODINGUTF-8确保中文属性不乱码。PostGIS 连接字符串里dbname、user按实际改。-nln basins指定表名。导入后可以用 SQL 做空间查询比如SELECT * FROM basins WHERE ST_Area(geom) 1000000;。4. 避坑与排查SHP 处理中最容易翻车的五个场景4.1 现象打开 SHP 后属性表全是乱码原因.dbf文件的字符编码和读取时指定的编码不一致。中文 Windows 下生成的 SHP 常用 GBK而 QGIS 或 GeoPandas 默认用 UTF-8。解决先看有没有.cpg文件里面会写编码。没有的话用gpd.read_file(..., encodinggbk)试或者用 QGIS 的「重新编码」功能转成 UTF-8。4.2 现象sjoin结果为空但两个图层明明有重叠原因坐标系不一致。一个图层是 EPSG:4326另一个是 EPSG:4547虽然视觉上重叠但数值范围差几个数量级。解决统一用to_crs转到同一坐标系再做空间操作。我一般会在sjoin前打印两个图层的crs和total_bounds确认范围接近。4.3 现象计算面积时结果为零或负数原因多边形几何无效比如自相交、环方向错误。解决用basins.geometry.is_valid检查对无效几何用basins.geometry.buffer(0)修复。buffer(0)有时能解决自相交但会轻微改变边界修复后要重新检查面积。4.4 现象ogrinfo报错Unable to open datasource原因路径含中文、空格或特殊字符或者.shp文件不完整缺少.shx或.dbf。解决把文件移到纯英文路径确保同名文件都在。如果是从压缩包直接打开先完整解压。4.5 现象导出 GeoJSON 后中文变成\uXXXX原因GeoJSON 规范默认用 Unicode 转义不是乱码。解决用ensure_asciiFalse参数或者在 QGIS 导出时勾选「不转义 Unicode」。如果是在 Web 端用浏览器能正确解析转义序列不影响显示。5. 进阶技巧用河流与流域数据做河网密度分级和验证河网密度是流域地貌的重要指标等于流域内河流总长度除以流域面积。用这份资料可以快速算出来并做分级渲染。先按 3.2 节把河流和流域连接再按BASIN_ID汇总长度最后除以面积。# 计算河网密度 density joined.groupby(BASIN_ID).agg( total_length_km(length_km, sum), area_km2(area_km2, first) ).reset_index() density[density] density[total_length_km] / density[area_km2] print(density.sort_values(density, ascendingFalse).head(10))agg里用first取面积因为同一流域内面积值相同。density单位是 km/km²数值越大说明河网越密集。你可以把结果导回 QGIS按density字段做渐变符号化直观看出哪些流域水系发育。验证方法拿一条已知长度的河流手动在 QGIS 里量一下和length_km对比。如果差异超过 5%检查投影是否合适。另外可以用rivers.geometry.length在投影坐标系下算不要在地理坐标系下算。# 验证单条河流长度 sample rivers_proj[rivers_proj[NAME] 长江] print(sample.geometry.length / 1000) # 应接近实际长度如果NAME字段不存在用FID或OBJECTID定位。长度单位是米除以 1000 得公里。这个习惯我每次拿到新 SHP 都会走一遍——先看坐标系再验一条已知要素确认无误才批量处理。从那以后我每次处理河流及流域数据都强制走一遍坐标系检查和单要素验证省去了很多返工。希望帮到你。本文还有配套的精品资源点击获取