简介黄河流域2000年河网矢量shp格式数据面向地理信息系统、水文水资源与生态遥感方向的科研人员、工程师及高校师生主要用于流域河网结构分析、历史演变对比与专题制图。数据采用WGS-84坐标系以shp矢量格式存储可在ArcGIS、QGIS等主流GIS软件中直接打开且已经过遥感影像叠加校正位置准确、拓扑关系完整可直接用于实验或科研使用。压缩包共7个文件包含shp主文件、dbf属性表、shx索引、prj坐标系定义及sbn/sbx空间索引等要素整体仅6.5MB下载与解压都十分便捷。目前已有1463人学习使用用户既能查看2000年黄河流域河网的整体格局也可与其他年份数据叠加对比省去自行矢量化与坐标配准的繁琐步骤快速获得可用于分析的基础数据。对于需要快速获取历史河网底图的科研与教学任务这份数据提供了一个可靠、直接的起点。1. 为什么是2000年的河网shp一次流域数据落地的起点拿到“黄河流域2000年河网矢量shp格式数据”时第一反应不是拉进ArcGIS里调样式而是先回答三个问题几何是线还是面、坐标系是什么、2000年这个时间截面代表哪一版水系。做水文模拟、河湖制图和流域空间统计的人几乎都会在某一步碰上这个数据名。shp格式能跨软件交换但只有一个.shp文件并不完整丢了.shx或.dbf就打不开坐标系不写清楚叠加遥感影像后河流能偏出几公里。下面按一条完整的数据处理链展开先拆文件结构、再做坐标系统一、然后加载成可用图层最后完成河流分级和shp转txt把数据导出成模型能直接用的格式。2. 河网shp数据到手怎么检查文件结构、字段与数据源2.1 一份河网shp拆开看四个必需文件与属性字段shapefile不是单文件格式。主文件 .shp 存放几何坐标.shx 是几何索引.dbf 是属性表河名、级别、长度、流向这些字段都记录在里面.prj 写坐标参考。部分数据还会带 .cpg 编码说明文件早期作业成果里经常缺失。拿到数据第一步先检查整套文件是否完整ls -lh /data/huanghe_2000/huanghe_river_2000.*这段命令列出该目录下所有同名拓展名文件从文件大小和修改时间能快速判断 .shx、.dbf、.prj 是否齐全缺 .prj 会在后续叠加底图时给出坐标未知警告。接着用 GDAL 的 ogr 读取要素类型和字段定义from osgeo import ogr ds ogr.Open(/data/huanghe_2000/huanghe_river_2000.shp, 0) layer ds.GetLayer() print(几何类型:, layer.GetGeomType()) print(要素数量:, layer.GetFeatureCount()) layer_defn layer.GetLayerDefn() for i in range(layer_defn.GetFieldCount()): fd layer_defn.GetFieldDefn(i) print(fd.GetName(), fd.GetType())这里参数 0 表示只读方式打开文件GetGeomType() 返回的如果是线类型就符合预期如果显示点或面说明文件已经被处理过不能当作原始河网使用。字段名和类型打印出来以后后面写分级和导出脚本时按实际字段名写而不是想当然用 Name。河网shp属性表常见字段大致如下不同作业单位的命名差别很大Level 也可能叫 Grade 或 Order需要以上面脚本读到的为准。字段名常见类型内容说明Name字符串河流名称如“黄河”“渭河”Level整数河网级别Strahler 分级结果Length浮点河段长度需确认单位是米还是千米Width浮点平均河宽多数数据没有GridCode整数原始栅格河网的网格编码2.2 2000年这个时间点意味着什么2000年河网数据代表的是世纪之交的水系格局和现在的河道相比差异主要集中在下游引黄灌区人工渠道与调水工程。用“2020全国道路shp下载”那类矢量检索方式去找“黄河流域shp”容易先搜出近年遥感解译产物时间面不同做年际对比时会出现系统性偏差。常见使用场景是水文模型输入SWAT建模时把该shp作为河道文件2000年作为率定期起点另一个场景是河湖演变分析拿它与近年同级数据比较主槽摆动范围。别只看文件名里的“黄河”数据通常覆盖整个黄河流域需要提前准备好流域边界shp做裁剪否则统计时会把渭河、汾河之外的邻近水系也算进去。因此拿到数据后不能只认文件名里的“2000”还要从元数据里确认数据源比例尺1:25万到1:100万之间比较常见。2.3 拿到数据先做四件事到手后按顺序处理。第一确认 .shp/.shx/.dbf/.prj 全套文件都在缺 .prj 后续无法定位。第二用上面的 ogr 脚本读属性表中文乱码时优先怀疑 dbf 编码。第三比对要素数量与资料说明里的河段数数量差得多说明转换过程丢过要素。第四检查无效几何检查代码如下from osgeo import ogr src ogr.Open(/data/huanghe_2000/huanghe_river_2000.shp, 0) layer src.GetLayer() invalid 0 total 0 for feat in layer: total 1 geom feat.geometry() if geom is None or not geom.IsValid(): invalid 1 print(f共 {total} 条无效几何 {invalid} 条)IsValid() 按 OGC 简单要素规则判断自相交、重复点等问题。数量大的河网里有几条无效几何正常如果超过总量的百分之一就先做几何修复再做任何空间分析否则后面与DEM提取河网做重合率统计时无效要素会直接拉低匹配结果。3. 坐标系与投影把河网shp统一到WGS84或CGCS20003.1 先分清经纬度坐标与投影坐标河网数据进场时可能已经是WGS84或CGCS2000地理坐标也可能是高斯-克吕格投影的平面坐标。快速区分方法坐标值在0到180之间是经纬度坐标值动辄百万且不带负号大概率是高斯平面坐标。黄河流域跨度大叠加遥感影像和做制图统计时选坐标系策略不一样可以按用途直接对号入座用途推荐坐标系补充说明快速浏览、叠加在线底图EPSG:4326 或 EPSG:44904490 对应 CGCS2000在线底图做动态投影即可全流域面积与密度统计Albers 等积圆锥投影黄河流域常用中央经线 105°E标准纬线 25°N、47°N局部工程出图高斯-克吕格 3 度带带号按河段所在经度选择先读 .prj 再对照做面积统计却坚持用WGS84地理坐标量算结果会掺杂纬度变形反过来只是叠加底图却转成Albers桌面软件实时投影会正常工作但导出给其他系统的数据容易被误用因为对方拿到的是已经变了坐标值的数据。3.2 用ogr2ogr做坐标转换并修正编码一条命令同时完成重投影和输出dbf编码设置是最省事的做法ogr2ogr -t_srs EPSG:4490 \ -lco ENCODINGUTF-8 \ /data/huanghe_cgcs2000/huanghe_river_2000.shp \ /data/huanghe_2000/huanghe_river_2000.shp参数含义-t_srs 指定输出坐标系-lco 是创建输出图层的选项这里把dbf写成UTF-8避免中文属性在后续步骤里乱码输出目录必须提前创建输出文件名不能与输入相同。如果原文件没有 .prj这条命令会报“无法读取坐标系”错误需要先在QGIS或ArcGIS里给数据“定义投影”指定一个能确定数据真实位置的坐标系再执行转换。注意定义投影只修改元数据描述不重算坐标真正的坐标变换是这里 -t_srs 做的事或者GIS里的“投影和变换”工具。3.3 三个常见坐标坑坑一shp在QGIS里显示正常发布到GeoServer后要素跑到海里。这种多数是 .prj 里写的是QGIS自定义文本GeoServer解析不了统一转成EPSG:4490或4326再发布。坑二两个数据都带.prj但底图是Web墨卡托河网是WGS84叠加时软件行为不同。ArcMap默认按数据框坐标系做动态投影低版本对未定义坐标系的图层不会自动转换需要先把数据框坐标系设置为目标系再做“投影和变换”重算一份。坑三平面坐标状态下读出的Length字段直接当真实距离用。黄河流域中下游影响不大上游青海、甘肃段投影变形累计明显做径流长度统计前应转成等距或等积投影再量算反过来地理坐标系下Length读出来的是度不是米必须先投影再算长度。4. QGIS/ArcGIS加载河网shp编码、样式、修复与拓扑检查4.1 加载河网shp先解决乱码再叠加遥感影像QGIS加载shp可以直接拖拽但中文属性乱码时打开“数据源管理器 → 矢量”选择对应的shp在选项里把编码从Autodetect改为GBK或UTF-8重新加载即可。ArcMap在图层属性→源里查看坐标系如果显示“未知”先定义坐标再叠加遥感影像否则影像和河网对不上。目视检查方法不复杂放大到一段弯曲河段看线中心是否落在影像河流中线附近偏差超过一两个像元优先怀疑坐标而不是数据本身画错。4.2 按Level字段分级显示主河、支流一眼分开河网全画成同一粗细很难读。惯例是用Level字段做分类渲染打开图层属性→符号化→按字段分类分类字段选Level线宽用表达式控制常用的是 scale_exp(level, 1, 10, 0.4, 2.2, 0.5)意思是1级支流画0.4毫米10级干流画2.2毫米中间按对数过渡。设置完成后主流和支流层次清楚。需要留意如果后续要做河网网络分析线方向必须统一为上游到下游否则流向字段全部失去意义。检查方向可以给线图层加箭头符号方向混乱时ArcGIS用编辑工具“翻转线”QGIS用处理工具箱里的“反转线方向”算法跑一遍。如果数据里没有等级字段退而求其次的做法是按Shape_Length降序排列把前20%的河段当作干流候选再按连续性人工合并。4.3 修复几何与拓扑检查shapechecker之外的自动化方案河网数据里常见的几何错误类型和处理手段如下表几何错误表现处理手段自相交线自己穿过自己修复几何或 shapechecker 逐条修复悬垂端点支流末端没汇入干流拓扑检查“不能有悬垂端点”后编辑重复线段同一段河画了两遍删除重复几何空几何属性有记录但没有任何坐标删除对应要素批量修复用QGIS处理脚本比shapechecker逐条点选效率高import processing processing.run(native:fixgeometries, { INPUT: /data/huanghe_2000/huanghe_river_2000.shp, OUTPUT: /data/huanghe_2000/huanghe_river_fixed.shp })INPUT 指定待修复数据OUTPUT 指定新文件不覆盖原图。使用shapechecker修复shp步骤通常是加载shp、执行检查、逐条确认、另存适合几百条要素的小数据要素超过十万级时先跑上面的脚本再用拓扑检查器验证。QGIS的拓扑检查器可以新建规则“不能有悬垂端点”“不能有伪节点”和“不能重叠”错误列表可以导出成点shp后续继续用脚本处理。“只保留外边界线”的需求也常在这一步出现只留黄河干流时不要做凸包那会得到包含全部支流的轮廓正解是按Name字段过滤或者按Level字段选中最高级导出。字段语义不明确时先看属性表再操作。5. 从河网shp提取分析结果河流分级、shp转txt与渔网统计5.1 没有现成等级字段时用DEM重算Strahler分级有些河网shp属性表里没有Level字段。这种情况常见做法是基于DEM重建填洼→流向→汇流累积量→阈值提取河网→Stream Order→栅格转矢量。阈值的选取以“生成河网总长度与2000年shp接近”为目标试算黄河流域上游山区阈值可以调小下游平原区要调大避免把田间沟渠也提成河道。如果想验证已有shp用空间连接统计两套河网的重合率低于60%先查投影一致性和阈值参数。5.2 把河网shp转txt坐标点序列一次导出给数值模型提供河道断面或点位文件时常用ogr2ogr直接转成带WKT的CSVogr2ogr -f CSV /data/river/huanghe_river.csv \ /data/huanghe_2000/huanghe_river_2000.shp \ -lco GEOMETRYAS_WKT每行是一条河段的WKT几何加属性。若模型需要逐折点坐标用Python遍历线几何输出文本from osgeo import ogr src ogr.Open(/data/huanghe_2000/huanghe_river_2000.shp, 0) layer src.GetLayer() with open(/data/river/river_points.txt, w, encodingutf-8) as f: for feat in layer: geom feat.geometry() name feat.GetField(Name) if geom is None: continue lines [] gtype geom.GetGeometryType() if gtype ogr.wkbLineString: lines [geom] # 单线直接处理 elif gtype ogr.wkbMultiLineString: lines [geom.GetGeometryRef(i) for i in range(geom.GetGeometryCount())] for line in lines: for i in range(line.GetPointCount()): x, y, _ line.GetPoint(i) f.write(f{name},{x:.6f},{y:.6f}\n)代码先按几何类型展开MultiLineString再逐点输出横纵坐标写成6位小数约0.1米精度。字段名Name先用前面的ogr脚本确认防止实际数据里是RiverName导致KeyError。5.3 用渔网分割shp做批量统计做全流域河网密度对比时用行政边界容易掩盖上中下游差异均匀渔网更合适。QGIS里“处理工具箱→创建网格”网格类型选矩形间距先按2km试再用“相交”让河网与网格相交输出结果带grid_id按grid_id汇总线长即可。间距不是越小越好1:25万的河网数据用500m网格只会让相交计算变慢统计结果并不会因此有数量级变化。相交会把跨网格的河段切断汇总时用相交后的长度不要拿原始Length字段替代否则统计口径会前后矛盾。整个流程可以放进图形模型修复几何→相交→按属性汇总一次跑完全流域配好参数后这条流水线可以直接复用于其它年份同范围河网shp。本文还有配套的精品资源点击获取 
长河君 · 2026/9/14 16:18:54 · 阅读 1.2k 黄河流域2000年河网shp数据：从检查到坐标转换与分级导出
简介黄河流域2000年河网矢量shp格式数据面向地理信息系统、水文水资源与生态遥感方向的科研人员、工程师及高校师生主要用于流域河网结构分析、历史演变对比与专题制图。数据采用WGS-84坐标系以shp矢量格式存储可在ArcGIS、QGIS等主流GIS软件中直接打开且已经过遥感影像叠加校正位置准确、拓扑关系完整可直接用于实验或科研使用。压缩包共7个文件包含shp主文件、dbf属性表、shx索引、prj坐标系定义及sbn/sbx空间索引等要素整体仅6.5MB下载与解压都十分便捷。目前已有1463人学习使用用户既能查看2000年黄河流域河网的整体格局也可与其他年份数据叠加对比省去自行矢量化与坐标配准的繁琐步骤快速获得可用于分析的基础数据。对于需要快速获取历史河网底图的科研与教学任务这份数据提供了一个可靠、直接的起点。1. 为什么是2000年的河网shp一次流域数据落地的起点拿到“黄河流域2000年河网矢量shp格式数据”时第一反应不是拉进ArcGIS里调样式而是先回答三个问题几何是线还是面、坐标系是什么、2000年这个时间截面代表哪一版水系。做水文模拟、河湖制图和流域空间统计的人几乎都会在某一步碰上这个数据名。shp格式能跨软件交换但只有一个.shp文件并不完整丢了.shx或.dbf就打不开坐标系不写清楚叠加遥感影像后河流能偏出几公里。下面按一条完整的数据处理链展开先拆文件结构、再做坐标系统一、然后加载成可用图层最后完成河流分级和shp转txt把数据导出成模型能直接用的格式。2. 河网shp数据到手怎么检查文件结构、字段与数据源2.1 一份河网shp拆开看四个必需文件与属性字段shapefile不是单文件格式。主文件 .shp 存放几何坐标.shx 是几何索引.dbf 是属性表河名、级别、长度、流向这些字段都记录在里面.prj 写坐标参考。部分数据还会带 .cpg 编码说明文件早期作业成果里经常缺失。拿到数据第一步先检查整套文件是否完整ls -lh /data/huanghe_2000/huanghe_river_2000.*这段命令列出该目录下所有同名拓展名文件从文件大小和修改时间能快速判断 .shx、.dbf、.prj 是否齐全缺 .prj 会在后续叠加底图时给出坐标未知警告。接着用 GDAL 的 ogr 读取要素类型和字段定义from osgeo import ogr ds ogr.Open(/data/huanghe_2000/huanghe_river_2000.shp, 0) layer ds.GetLayer() print(几何类型:, layer.GetGeomType()) print(要素数量:, layer.GetFeatureCount()) layer_defn layer.GetLayerDefn() for i in range(layer_defn.GetFieldCount()): fd layer_defn.GetFieldDefn(i) print(fd.GetName(), fd.GetType())这里参数 0 表示只读方式打开文件GetGeomType() 返回的如果是线类型就符合预期如果显示点或面说明文件已经被处理过不能当作原始河网使用。字段名和类型打印出来以后后面写分级和导出脚本时按实际字段名写而不是想当然用 Name。河网shp属性表常见字段大致如下不同作业单位的命名差别很大Level 也可能叫 Grade 或 Order需要以上面脚本读到的为准。字段名常见类型内容说明Name字符串河流名称如“黄河”“渭河”Level整数河网级别Strahler 分级结果Length浮点河段长度需确认单位是米还是千米Width浮点平均河宽多数数据没有GridCode整数原始栅格河网的网格编码2.2 2000年这个时间点意味着什么2000年河网数据代表的是世纪之交的水系格局和现在的河道相比差异主要集中在下游引黄灌区人工渠道与调水工程。用“2020全国道路shp下载”那类矢量检索方式去找“黄河流域shp”容易先搜出近年遥感解译产物时间面不同做年际对比时会出现系统性偏差。常见使用场景是水文模型输入SWAT建模时把该shp作为河道文件2000年作为率定期起点另一个场景是河湖演变分析拿它与近年同级数据比较主槽摆动范围。别只看文件名里的“黄河”数据通常覆盖整个黄河流域需要提前准备好流域边界shp做裁剪否则统计时会把渭河、汾河之外的邻近水系也算进去。因此拿到数据后不能只认文件名里的“2000”还要从元数据里确认数据源比例尺1:25万到1:100万之间比较常见。2.3 拿到数据先做四件事到手后按顺序处理。第一确认 .shp/.shx/.dbf/.prj 全套文件都在缺 .prj 后续无法定位。第二用上面的 ogr 脚本读属性表中文乱码时优先怀疑 dbf 编码。第三比对要素数量与资料说明里的河段数数量差得多说明转换过程丢过要素。第四检查无效几何检查代码如下from osgeo import ogr src ogr.Open(/data/huanghe_2000/huanghe_river_2000.shp, 0) layer src.GetLayer() invalid 0 total 0 for feat in layer: total 1 geom feat.geometry() if geom is None or not geom.IsValid(): invalid 1 print(f共 {total} 条无效几何 {invalid} 条)IsValid() 按 OGC 简单要素规则判断自相交、重复点等问题。数量大的河网里有几条无效几何正常如果超过总量的百分之一就先做几何修复再做任何空间分析否则后面与DEM提取河网做重合率统计时无效要素会直接拉低匹配结果。3. 坐标系与投影把河网shp统一到WGS84或CGCS20003.1 先分清经纬度坐标与投影坐标河网数据进场时可能已经是WGS84或CGCS2000地理坐标也可能是高斯-克吕格投影的平面坐标。快速区分方法坐标值在0到180之间是经纬度坐标值动辄百万且不带负号大概率是高斯平面坐标。黄河流域跨度大叠加遥感影像和做制图统计时选坐标系策略不一样可以按用途直接对号入座用途推荐坐标系补充说明快速浏览、叠加在线底图EPSG:4326 或 EPSG:44904490 对应 CGCS2000在线底图做动态投影即可全流域面积与密度统计Albers 等积圆锥投影黄河流域常用中央经线 105°E标准纬线 25°N、47°N局部工程出图高斯-克吕格 3 度带带号按河段所在经度选择先读 .prj 再对照做面积统计却坚持用WGS84地理坐标量算结果会掺杂纬度变形反过来只是叠加底图却转成Albers桌面软件实时投影会正常工作但导出给其他系统的数据容易被误用因为对方拿到的是已经变了坐标值的数据。3.2 用ogr2ogr做坐标转换并修正编码一条命令同时完成重投影和输出dbf编码设置是最省事的做法ogr2ogr -t_srs EPSG:4490 \ -lco ENCODINGUTF-8 \ /data/huanghe_cgcs2000/huanghe_river_2000.shp \ /data/huanghe_2000/huanghe_river_2000.shp参数含义-t_srs 指定输出坐标系-lco 是创建输出图层的选项这里把dbf写成UTF-8避免中文属性在后续步骤里乱码输出目录必须提前创建输出文件名不能与输入相同。如果原文件没有 .prj这条命令会报“无法读取坐标系”错误需要先在QGIS或ArcGIS里给数据“定义投影”指定一个能确定数据真实位置的坐标系再执行转换。注意定义投影只修改元数据描述不重算坐标真正的坐标变换是这里 -t_srs 做的事或者GIS里的“投影和变换”工具。3.3 三个常见坐标坑坑一shp在QGIS里显示正常发布到GeoServer后要素跑到海里。这种多数是 .prj 里写的是QGIS自定义文本GeoServer解析不了统一转成EPSG:4490或4326再发布。坑二两个数据都带.prj但底图是Web墨卡托河网是WGS84叠加时软件行为不同。ArcMap默认按数据框坐标系做动态投影低版本对未定义坐标系的图层不会自动转换需要先把数据框坐标系设置为目标系再做“投影和变换”重算一份。坑三平面坐标状态下读出的Length字段直接当真实距离用。黄河流域中下游影响不大上游青海、甘肃段投影变形累计明显做径流长度统计前应转成等距或等积投影再量算反过来地理坐标系下Length读出来的是度不是米必须先投影再算长度。4. QGIS/ArcGIS加载河网shp编码、样式、修复与拓扑检查4.1 加载河网shp先解决乱码再叠加遥感影像QGIS加载shp可以直接拖拽但中文属性乱码时打开“数据源管理器 → 矢量”选择对应的shp在选项里把编码从Autodetect改为GBK或UTF-8重新加载即可。ArcMap在图层属性→源里查看坐标系如果显示“未知”先定义坐标再叠加遥感影像否则影像和河网对不上。目视检查方法不复杂放大到一段弯曲河段看线中心是否落在影像河流中线附近偏差超过一两个像元优先怀疑坐标而不是数据本身画错。4.2 按Level字段分级显示主河、支流一眼分开河网全画成同一粗细很难读。惯例是用Level字段做分类渲染打开图层属性→符号化→按字段分类分类字段选Level线宽用表达式控制常用的是 scale_exp(level, 1, 10, 0.4, 2.2, 0.5)意思是1级支流画0.4毫米10级干流画2.2毫米中间按对数过渡。设置完成后主流和支流层次清楚。需要留意如果后续要做河网网络分析线方向必须统一为上游到下游否则流向字段全部失去意义。检查方向可以给线图层加箭头符号方向混乱时ArcGIS用编辑工具“翻转线”QGIS用处理工具箱里的“反转线方向”算法跑一遍。如果数据里没有等级字段退而求其次的做法是按Shape_Length降序排列把前20%的河段当作干流候选再按连续性人工合并。4.3 修复几何与拓扑检查shapechecker之外的自动化方案河网数据里常见的几何错误类型和处理手段如下表几何错误表现处理手段自相交线自己穿过自己修复几何或 shapechecker 逐条修复悬垂端点支流末端没汇入干流拓扑检查“不能有悬垂端点”后编辑重复线段同一段河画了两遍删除重复几何空几何属性有记录但没有任何坐标删除对应要素批量修复用QGIS处理脚本比shapechecker逐条点选效率高import processing processing.run(native:fixgeometries, { INPUT: /data/huanghe_2000/huanghe_river_2000.shp, OUTPUT: /data/huanghe_2000/huanghe_river_fixed.shp })INPUT 指定待修复数据OUTPUT 指定新文件不覆盖原图。使用shapechecker修复shp步骤通常是加载shp、执行检查、逐条确认、另存适合几百条要素的小数据要素超过十万级时先跑上面的脚本再用拓扑检查器验证。QGIS的拓扑检查器可以新建规则“不能有悬垂端点”“不能有伪节点”和“不能重叠”错误列表可以导出成点shp后续继续用脚本处理。“只保留外边界线”的需求也常在这一步出现只留黄河干流时不要做凸包那会得到包含全部支流的轮廓正解是按Name字段过滤或者按Level字段选中最高级导出。字段语义不明确时先看属性表再操作。5. 从河网shp提取分析结果河流分级、shp转txt与渔网统计5.1 没有现成等级字段时用DEM重算Strahler分级有些河网shp属性表里没有Level字段。这种情况常见做法是基于DEM重建填洼→流向→汇流累积量→阈值提取河网→Stream Order→栅格转矢量。阈值的选取以“生成河网总长度与2000年shp接近”为目标试算黄河流域上游山区阈值可以调小下游平原区要调大避免把田间沟渠也提成河道。如果想验证已有shp用空间连接统计两套河网的重合率低于60%先查投影一致性和阈值参数。5.2 把河网shp转txt坐标点序列一次导出给数值模型提供河道断面或点位文件时常用ogr2ogr直接转成带WKT的CSVogr2ogr -f CSV /data/river/huanghe_river.csv \ /data/huanghe_2000/huanghe_river_2000.shp \ -lco GEOMETRYAS_WKT每行是一条河段的WKT几何加属性。若模型需要逐折点坐标用Python遍历线几何输出文本from osgeo import ogr src ogr.Open(/data/huanghe_2000/huanghe_river_2000.shp, 0) layer src.GetLayer() with open(/data/river/river_points.txt, w, encodingutf-8) as f: for feat in layer: geom feat.geometry() name feat.GetField(Name) if geom is None: continue lines [] gtype geom.GetGeometryType() if gtype ogr.wkbLineString: lines [geom] # 单线直接处理 elif gtype ogr.wkbMultiLineString: lines [geom.GetGeometryRef(i) for i in range(geom.GetGeometryCount())] for line in lines: for i in range(line.GetPointCount()): x, y, _ line.GetPoint(i) f.write(f{name},{x:.6f},{y:.6f}\n)代码先按几何类型展开MultiLineString再逐点输出横纵坐标写成6位小数约0.1米精度。字段名Name先用前面的ogr脚本确认防止实际数据里是RiverName导致KeyError。5.3 用渔网分割shp做批量统计做全流域河网密度对比时用行政边界容易掩盖上中下游差异均匀渔网更合适。QGIS里“处理工具箱→创建网格”网格类型选矩形间距先按2km试再用“相交”让河网与网格相交输出结果带grid_id按grid_id汇总线长即可。间距不是越小越好1:25万的河网数据用500m网格只会让相交计算变慢统计结果并不会因此有数量级变化。相交会把跨网格的河段切断汇总时用相交后的长度不要拿原始Length字段替代否则统计口径会前后矛盾。整个流程可以放进图形模型修复几何→相交→按属性汇总一次跑完全流域配好参数后这条流水线可以直接复用于其它年份同范围河网shp。本文还有配套的精品资源点击获取