黄河流域数据集：DEM、水系、界线一次配齐，水文分析不再东拼西凑

简介一套覆盖黄河流域多要素的地理数据集面向地理信息系统、水文学及区域规划方向的研究者与从业者可为流域地形分析、水资源管理和生态保护研究提供基础底图与矢量图层。内容包含数字高程模型DEM、干流与主要支流、湖泊分布、流域界线、黄土高原地区界以及城市坐标点等要素既能提取坡度、坡向、汇水路径等地形因子也便于开展洪水风险评价、水土流失敏感性分析和湿地生态监测。压缩包共178个文件以shp、dbf、tif等常用格式为主并包含gdb地理数据库和配套投影文件整体大小约36.07MB文件命名规范、目录层次清晰可直接导入ArcGIS、QGIS等GIS平台。目前已有1068人学习下载。借助该数据集可针对黄河流域开展空间叠加、制图建模与专项分析免去自行检索和预处理多源空间数据的繁琐流程尤其适用于科研教学、流域规划及工程前期研究场景。黄河流域数据集DEM、水系、界线一次配齐做水文分析终于不用再东拼西凑了做黄河流域相关研究的同学应该都有过这种经历开题报告交上去导师说“先把手头数据整理好”然后就开始了漫长的数据苦旅。矢量边界找一个网站DEM去另一个平台下水系数据好不容易找到了又是十年前的旧版本坐标系统一都花了三天。我前阵子拿到一套整理好的黄河流域数据集里面从DEM到干流、支流、湖泊、流域界线、黄土高原地区界全都有说实话第一反应是“这种东西居然真的有人愿意整理出来分享”。这篇文章就把这套数据的实际内容、使用方法和我在实操中踩过的坑详细梳理一遍给正在做黄河流域研究或者准备做水文分析的读者一个参考。这套数据的核心构成其实不复杂矢量部分覆盖了黄河流域城市点、干流、重点支流、湖泊、流域界线、黄土高原地区界栅格部分提供了DEM数据。但就是这几样东西恰恰是水文分析、水土保持研究、流域生态评价和制图出图最常用到的底图数据。无论你是做SWAT模型、HEC-HMS水文模拟还是单纯想画一张黄河流域专题图这套数据都能直接作为基础底图使用省掉大量数据收集和预处理的时间。1. 为什么黄河流域研究急需一套整合好的基础数据1.1 黄河流域的复杂性决定了数据准备是硬骨头黄河流域从青藏高原一直延伸到渤海之滨跨越了九个省区从上游的高原草甸到中游的黄土高原再到下游的华北平原地形地貌差异极大。黄河干流全长约5464公里沿途接纳了渭河、汾河、湟水、洮河等几十条重要支流流域面积约75万平方公里。这个尺度下任何一个研究项目都会面临一个共同问题数据源分散、格式多样、精度不一。我见过太多人卡在数据准备阶段。有的人从地理空间数据云下了DEM发现边界裁剪不干净有的人从网上找到一份黄河流域边界结果坐标系是老的北京54跟最新研究区的数据套不在一起还有人花了一个星期整理支流水系最后发现不同来源的河网对不上支流位置偏移好几公里。这些问题在黄河流域这种大尺度研究区域尤其明显——数据范围越大拼接、配准、统一的难度就越高。与其每次从头找数据不如直接用一套已经经过整合和坐标统一的数据集作为基底再根据具体研究需求做局部补充。这套数据的价值恰恰在于此它把最常用、最基础的图层集合到了一起拿到手就能直接用而不是再从零开始东拼西凑。1.2 一套数据在手能省掉多少重复劳动从实际操作角度算一笔账如果完全靠自己收集黄河流域DEM至少要花半天到一天下载和拼接城市点坐标需要从公开地理库整理干流和支流水系要么找旧版矢量图手动校正要么用DEM自己提取再跟真实河网对比验证流域界线更要翻文献找来源黄土高原地区界更是冷门数据不知道要翻多少篇论文才能找到可下载的版本。整套下来一个人少说也得两到三周才能整理出一份勉强能用的底图数据还得祈祷各图层之间坐标统一、边界吻合。而拿到这套整理好的数据之后我大致验证了一下各图层坐标体系一致边界重叠情况良好直接进入分析和制图阶段完全可行。对于做毕业设计的学生、刚接触流域研究的年轻研究者或者项目周期紧的从业人员来说这相当于把最枯燥、最磨人的环节直接跳过了。2. 这份数据集里到底装了什么好东西2.1 城市点、干流与湖泊最基础的空间骨架先说城市点。这个图层收集了黄河流域范围内的主要城市点位数据千万别小看这个点图层它在实际研究中的用处比想象中大得多。做空间插值的时候气象站、水文站的数据往往需要跟城市位置关联展示做经济地理分析的时候城市点是基础单元哪怕只是画一张研究区概况图标注城市位置也是最直观的空间参考。我拿到这套数据后最先干的事就是把城市点叠到DEM上做了一张地势图上中下游的关系一眼就清楚了。干流和湖泊数据是黄河流域水系的核心骨架。黄河干流在数据里是以线要素形式存储的沿着河道走势完整覆盖了从源头到河口的主河道。湖泊数据则覆盖了流域内主要湖泊和水库像扎陵湖、鄂陵湖这些源头湖泊以及乌梁素海等重点湖泊都有体现。这里要说一句题外话数据说明里写的“重点直流”应该是“重点支流”的笔误做GIS数据的人对这类错别字早已见怪不怪——之前还见过把“秦岭”写成“秦玲”的不影响用就行。水系要素在实际分析中可以用来做缓冲区分析、计算河流长度、提取河流纵剖面、做河流分级甚至结合DEM做距离分析模拟洪水淹没范围。它和黄河流域城区、湖泊、行政边界叠加在一起就是一张标准的流域地理底图。2.2 流域界线与重点支流水文分析的核心输入流域界线这个图层是整个数据集里含金量比较高的一个。它划定的是黄河流域的完整集水范围。水文分析的第一步——确定研究区域边界——就要用它。你下载的DEM原始影像往往覆盖范围远大于流域边界如果没有精确的流域界线做裁剪后续步所有的栅格计算都会包含大量无关注区域既增加计算量又影响结果精度。重点支流图层单独列出来也说明整理者对水文研究的实际需求理解到位。黄河流域的重点支流就那几条渭河、汾河、湟水、洮河、大黑河、窟野河、无定河、洛河、沁河等。这些支流的径流量、输沙量直接关系到黄河干流的水沙变化是水文模型率定和验证的关键参照。单独建一个重点支流图层做研究区示意图、支流对比分析的时候非常方便不用再手工从全量河网里裁剪提取了。2.3 黄土高原地区界一个容易被忽略但极其重要的图层看到数据清单里有黄土高原地区界我当时就觉得整理者肯定是真正做过黄河流域研究的人。黄土高原的边界没有深入研究的人往往想不起来要这个数据但对搞黄河泥沙、水土保持研究的人来说这个界线比很多“大路货”图层都重要。黄土高原位于太行山以西、日月山以东、长城以南、秦岭以北总面积约64万平方公里是黄河泥沙的主要来源区。黄河“一碗水半碗沙”的说法根源就在黄土高原。做水土流失治理效益评价、退耕还林还草空间选址、泥沙来源分析都需要以黄土高原地区界作为基本分析范围。单独提供这个图层意味着你可以直接把研究区裁剪到黄土高原内部不用自己从一堆省级行政区里手动合并拼出来。这个图层的使用场景还有一个而且我私下觉得特别实用在做土地利用分析或者植被覆盖度分析的时候黄土高原地区界可以作为空间分层依据——界内和界外的土地利用结构、植被覆盖变化规律完全不同分开统计才能看出真正的区域差异。3. DEM数据在水文分析中的实操要点3.1 分辨率怎么选12.5米还是30米这套数据集里的DEM是哪一格分辨率的我没细看但不管拿到什么分辨率都要结合自己的研究尺度来选。热词列表里“12.5米dem下载”、“alos 12.5米dem数据”被频繁搜索说明大家都很关注高分辨率DEM。目前公开免费且比较常用的DEM就那几款DEM数据源分辨率适合尺度获取平台ALOS PALSAR12.5米中小流域、精细地形分析ASF DAAC、OpenTopographySRTM30米原90米中尺度流域、区域分析地理空间数据云、USGSGDEM V230米大尺度区域、制图底图地理空间数据云GLO-3030米全球性研究Copernicus Data Space如果你做的是整个黄河流域尺度的大范围分析30米SRTM或者GDEM就够用了数据量适中处理速度能接受提取出来的河网和实际水系的吻合度也还可以。但如果你是做某个支流的子流域模拟比如无定河流域的水沙过程模拟强烈建议用12.5米ALOS数据它在地形细节上的刻画明显好于30米数据提取出来的河网更平滑、更贴近实际河道走势。需要注意的一点是ALOS 12.5米数据是分幅存储的覆盖黄河流域可能需要拼接几十景影像。拼接的时候要注意条带之间是否有接边误差如果直接做填洼和水流方向分析接边处的错误会造成伪河道。建议先拼接、投影、裁剪再做后续水文预处理。3.2 用DEM提取河网与计算流域面积的完整流程DEM拿到手里不是为了看的而是为了做分析。在ArcGIS里用DEM提取河网的流程已经是非常经典的操作了填洼Fill→ 计算流向Flow Direction→ 计算汇流累积量Flow Accumulation→ 设置阈值提取河网Raster Calculator→ 转为矢量Stream to Feature。填洼是整个流程最关键的步骤。真实地表存在很多人工建筑、噪声造成的伪洼地如果直接计算流向水流会被困在这些洼地里面形成错误的水系。但填洼也不是越狠越好过度填洼会抹平真实地形特征所以参数要试。实际操作中我一般先看原始DEM的最大洼地深度然后设置一个略高于噪声水平的限制值做条件填洼。提取河网的阈值怎么定这里有个经验规律阈值越低提取的河网越密集伪河道越多阈值越高河网越稀疏细小支流会被丢弃。一个比较笨但有效的方法是设置不同阈值反复试把提取结果和真实支流数据叠在一起看。比如你对照重点支流图层发现阈值10万时渭河、汾河都能完整呈现但无定河的细小上游河道丢了就适当调低阈值。大部分大尺度分析集水面积阈值在5万到20万平方米之间比较合适。计算流域面积就更直接了。在提取出的河网上选择你关心的出水口位置添加一个点要素然后用Spatial Analyst里的Watershed工具把流向栅格和出水口点一起输入就能生成这个出水口以上的汇水区再用Zonal Geometry或直接查看面要素的属性表面积字段就得到流域面积了。很多人做流域面积计算容易忽略一步——确保投影坐标系是等积投影如果用的是地理坐标系WGS84面积单位是度算出来的数字没有任何意义。4. 在ArcGIS里把数据集用起来的完整流程4.1 数据加载与坐标系检查拿到这套数据之后第一步不是急着做分析而是检查坐标系。这是所有GIS项目里最值得养成习惯的一步。我处理数据这些年碰到的各种“灵异事件”——要素凭空偏移、面积算错、叠加对不上——九成以上都是坐标系问题引起的。把数据加载进ArcMap或ArcGIS Pro之后右键图层打开属性在Source选项卡里查看坐标系信息。如果所有图层显示的都是同一套坐标系比如CGCS2000或者WGS84那恭喜你这套数据整理者确实用心了。如果发现坐标系不一致就要用Project工具统一转换后再进行任何分析操作。这里补充一个投影坐标系的建议做黄河流域全流域尺度的面积量算、密度分析用Albers等积圆锥投影比较合适它是国内大区域制图的标准投影面积计算不会产生明显误差。如果做的是局部精细分析比如某个支流子流域的坡度坡向计算就投影到该区域对应的UTM分区精度更高。4.2 一次典型的流域分析实战我拿到数据之后做了一次典型的小实验用DEM提取黄河中游某重点支流的子流域边界再跟数据里给的流域界线做对比。这里把完整操作记录下来供参考。先加载DEM数据用Fill工具填洼。填洼结束后用Flow Direction计算流向再运行Flow Accumulation得到汇流累积量栅格。因为只是测试我在Spatial Analyst工具里选了Con表达式做河网提取汇流累积量大于10万的栅格赋值为1否则为NoData。然后Stream to Feature转为矢量河网。把生成的河网跟数据里的重点支流图层做对比观察中游主河道走向基本吻合但提取结果在上游源区多出了一些细碎河道这是阈值设置偏低造成的不影响大方向判断。接着我用Watershed工具在重点支流汇入黄河干流的位置设置出水口提取了这一支流的整个汇水区跟数据集自带的流域界线做空间叠加两者的边界走向高度一致仅在局部山脊区域有细微偏差。最后用Zonal Geometry as Table统计提取出的子流域面积对比数据集流域界线对应的理论集水面积误差在3%以内。这说明这套数据集的DEM精度和流域界线质量都经得起验证直接用于研究基底没有问题。5. 常见问题与排查技巧实录5.1 坐标对不上、边界偏差过大很多人拿到多图层数据后做的第一件事就是把所有图层叠加在一起然后发现水系和边界对不上或者城市点位置飘到了河道中央。遇到这种情况先不要急着做地理配准也不要手动移动要素。第一步是确认所有图层的坐标系是否一致第二步确认范围是否合理打开各图层的属性查看Extent范围第三步才是考虑数据源本身的问题。这套数据各图层的吻合度比我预想的好但如果你拿到的版本在局部区域存在偏差大概率是原始数据源拼接时造成的。处理办法是寻找一个精度较高的底图作为参考比如高分影像或者实测控制点对偏差较大的局部要素做空间校正不建议全局盲目配准。5.2 DEM负值、数据空洞与填洼误区DEM数据常见的三大坑这套数据也不能完全避免。负值问题水面和特殊地形区域往往会出现负高程这是数据采集时水面反射导致的通常出现在湖泊和水库区域。处理方式是检查DEM的数值范围如果负值数量较少且集中在已知水体区域可以直接用周围有效高程插值替换。数据空洞在山区陡峭地形下雷达影像经常出现黑色空洞。如果空洞面积小用栅格插值或者邻域分析可以补上空洞面积大就找另一期的DEM数据互补填补千万不要直接忽略。填洼误区很多人拿到DEM先无脑Fill这是最大的错误。填洼是为了消除伪洼地但真实存在的洼地比如喀斯特地区的天坑、人工水库也会被一并填平导致后续提取地表径流路径失真。正确的做法是先计算洼地深度人工判断哪些是真洼地哪些是伪洼地再做相应处理。5.3 提取河网与真实水系对不上DEM提取的河网和真实水系对不上基本是两个原因一是DEM精度不够平坦地区的河道在DEM上没有足够的高差信号水流方向计算不出来二是人为改道和水利工程的影响黄河下游是典型的地上河河道比地面高好几米DEM上根本看不到正常的汇水关系。遇到这种情况我的建议是不要试图只用DEM提取全流域河网。下游平原区直接用手工数字化的真实水系数据这套数据集里的干流和支流图层就能用来替代提取结果中游丘陵山区用DEM提取然后两段拼接。以实际需求为导向不纠结于所有河段都来自DEM自动提取。整理自己的数据底图永远是最费时间也最不出成果的环节。这套黄河流域数据集最让我满意的地方不是单一图层的精度有多高而是所有数据能直接拼在一起用城市点、干流、支流、湖泊、流域界线、黄土高原界还有DEM互相之间的空间关系是正确的。对于着急做研究的人来说这比什么都重要。最后分享一个个人习惯拿到任何数据集第一件事不要急着做分析先加载所有图层把符号化调整到合适状态叠加检查一遍边界有没有漏洞、要素有没有明显错位然后截一张全览图存档。这样后面每一步分析出了异常都能快速判断是数据问题还是操作问题省下的时间绝对远超一开始花掉的十分钟。本文还有配套的精品资源点击获取