基于DEM的青龙河流域数字河网提取.docVIP

基于DEM的青龙河流域数字河网提取 　　摘 要：简单叙述了流域河网提取的两种办法，并以ArcGIS为操作平台，介绍了基于坡面径流漫流模型的数字河网提取原理和方法，其流程可简单分为：DEM的预处理、水流流向的确定、汇流累积的计算以及河网生成。该文基于SRTM DEM数据以青龙河水域为研究对象进行流域特征信息提取，提取结果与国家河网数据基本吻合，并就不同阈值划分下的河网进行比较。 　　关键词：河网 SRTM DEM ArcGIS 青龙河流域 　　中图分类号：P33 文献标识码：A 文章编号：1672-3791（2015）05（b）-0024-02 　　Abastract：The paper briefly describes two methods of extracting river basins.With the Hydrology Analysis tools of ArcGIS， the principle and method of ground water overflow modelare well introduced，the process of which includes the flow direction， flow accumulation， drainage network andthedrainagenetworkextraction. Based on the Qinglong River watershed SRTM DEM data， the fact that the river networksturn out with different threshold values is discussed. And the extraction results were verified and in accordance with the basic national river network data. 　　Key Words：drainagenetwork； SRTM DEM； ArcGIS； Qinglong River 　　流域河网水系是描述一个地区地理地貌和水文特征基本的重要条件，对流域整体规划、水资源优化配置、水利水电开发利用以及防洪抗旱、航运、养殖、灌溉和旅游等都具有十分重大的意义[1]。目前，利用DEM数据提取流域数字河网主要有谷点提取模型和基于流向提取模型两种[2]。谷点提取模型最初是由Greysukh，Pekuker， Douglas[3]等人提出，其原理是通过对比所研究网格的高程与相邻网格的高程找出DEM中的谷点单元，再将谷点单元连接后形成河网[4]。该方法的问题是利用谷点单元形成的河网往往不连续，且在地形起伏较小的平原地区难以识别谷点单元[5]。而目前广泛应用的算法是基于O’Callaghan和 Mark提出的坡面径流模拟算法，其原理是先确定流域各栅格间的流向，再根据每个栅格点上游的计算得到汇流累积量，最后按不同阈值提取出河网。其中确定流向的方法多采用D8的算法，曹玲玲[6]、吴冰[7]、周国园[8]、沈中原[9]、热紫燕[10]等人都在不同流域基于坡面径流模拟算法进行了河网提取。该文主要介绍利用DEM数据基于坡面径流漫流模型的数字河网提取方法，并在ArcGIS环境下通过青龙河流域的实例进行分析。 　　1 研究区背景概况 　　1.1 流域概况 　　青龙河是秦皇岛市的重要水源地，青龙河流域是秦皇岛市限制纳污控制考核的重要水功能区。青龙河流经滦河流域东侧，地处东经118°37′～119°37′，北纬39°51′ ～41°07′之间，流域面积6340 km2，占滦河流域面积的14%，是滦河水系年平均来水量最大的支流。该文拟以桃林口水坝址以上流域为例，进行流域数字河网的提取研究。 　　1.2 数据来源 　　DEM的数据来源是SRTM DEM（Shuttle Radar Topography Mission）数据，是2000年2月11日美国发射的“奋进”号航天飞机上搭载的SRTM系统采集的地表面数据。SRTM DEM数据采用WGS1984投影参数，通过中国科学院计算机网络信息中心国际科学数据镜像网站可下载空间分辨率为 3″（约90m）的中国区域数据[11]。由于下载的图像数据不一定刚好在一个图幅范围内，需对下载的数据进行剪裁及修正，从而得到青龙河流域部分。 　　2 河网提取方法 　　2.1 DEM预处理 　　在利用SRTM数据提取流域的数字河网之前，首先需要对SRTM数据进行预处理，以消除SRTM中存在的无效数据区域。此外，由于原始的DEM图像中数据缺陷或是一些真实地形的存在，使得DEM表面存在着一些凹陷的区域，使水流无法按照其本身流向流出