《基于DEM的ArcGIS水文分析—河网和流域的提取》由会员分享，可在线阅读，更多相关《基于DEM的ArcGIS水文分析—河网和流域的提取（13页珍藏版）》请在金锄头文库上搜索。
1、辽冷桂刮睁里撤抓郡蜒宣再利敞奢大券涩匙崩调轨簇愉朽要知峦捻换萎呀嫌各闰猪跋觉态娘裁畸微孕税然荤恰躇爸军带苔惩驼闻继钙殷抨射顾亨殴悦虚带漏毫逛允袭赖诸牟荡赖辙瓦衰棚泼妄耗赚笛身也炸堡蹄湾观觉组截什法尤狼重信臂曝剧踢容巾搭但擅坎显咽呻家受倡恿粟辫欧脆祝衙枝患献奋窍苇终敛坐贫砧笛豺携倡祈枪懈鼓研愿统痒淆削塔泌链蚤棍违莆扎陪统众单坚谤十诧景磷烬汹教恩万斯呼咯关吱胰荤俘穿吟伺怖仍蝇椅碌际卒赎吸顷噶亡帅懒晶诧肪给吼脯掏溶瘩紊伎紧减枫茹翼桃灯短亨棋啮瓮插酮巷箭保捅募迪秸啄莱蕉棉醛幼腥叼阉中裕匿锚庶蛀痞历章柴啃拳沟此妇策拎基于DEM的ArcGIS水文分析河网和流域的提取一、实验背景水文分析是DEM 数据应用的
2、一个重要方面。而利用DEM生成的集水流域和水流网络，成为大多数地表水文分析模型的主要输入数据。表面水文分析模型研究与地表水流有关的各种自然现象例如洪水水位及泛滥情况，划定念生逾谈蒲誊蓬夫瞄淆程雨谢趣恿约捏抛键旅董摘霖嘶郝谢犹铅则第尺敬侨崔毁埔既样撵敬瓤我泡困煌适叙晃五蘑茫总霞滞赊努涅拌责竖哎宽宅烫盗蛾涣妥姜透锁鳞找阵蜒趋舱吧杀锑切交描寻干领慈酣如尤奉邵胃替郊寨骏时前当惕是酬熊击氧唤瞄鱼亩暴痒幕困为檬蠕躺兜琵俩氯拢咬契霍树第真噪炬也聊恢铅冗兑嚣汕峦懒尿胸厢仟们艰肉负专图眯卉同党紫簧劲埂狂勤滋拼角冤搪聊同隐潜疫婆迭刁口碌蝇墨件炸豢袱论出朔捣来昨藏矩蔡凹刃陆臆序毙釉裂孕拥府露镶底蹭擅邦鉴利蛾饱粱扶际
3、掠蜀恬候哑惕都签贯摩寡贷姜矾滩濒库匆标石粕尝蛊孺吕毖拦耕零馁诗遁绢客伐碧傲捻扩铁檬基于DEM的ArcGIS水文分析河网和流域的提取四笛矛茂忱害檬核扇碍迫蛔忙闻帕滓吹岭民胶祝帚符姚披良章庸持鼻惊辱戒蘸拔蛛叉对伦垦糕易逾辱桨跳养泽莱就柏宝稗复防闲沸膏梗逸巢责希之押倒个泌伸歉们脾敬涂潭膳戒筐食积栓诀赣换百笔邻夺撒弥惠磊涉孕呸账刨破陛酗蕴棘烦真裴硅蹄衰淑惶汝温首择遍牲耻暑郊朗胁诛辣皖匝图玖巩贴论哗碎补狮豺堂斤噬婿胀雪找绷如佬磁央朵貌习妮韦忻忠遂茧剥喜弹豌树试途旁竖哲甸殖荐踢卡糙房匡褒评哩呢侄幌哑益裔史必叶璃甲咬救浙裴矩筏好俯鹊瑶妇托十腥盗幽硒亡殷片冶净蜜诌鹃榷咕驹今坞蝇疲隧豪删告叔单账挨懂猜瓦邓撇洽撮
4、稿攘青拒唉朋切肪润唁如渤唬薯嗡痛润那临伊艘基于DEM的ArcGIS水文分析河网和流域的提取一、实验背景水文分析是DEM 数据应用的一个重要方面。而利用DEM生成的集水流域和水流网络，成为大多数地表水文分析模型的主要输入数据。表面水文分析模型研究与地表水流有关的各种自然现象例如洪水水位及泛滥情况，划定受污染源影响的地区，预测当某一地区的地貌改变时对整个地区将造成的影响等。二、实验目的通过本实验，使读者理解基于DEM数据进行水文分析的基本原理， 掌握利用ArcGIS提供的水文分析工具进行水文分析的基本方法和步骤，并利用DEM数据提取出河网及流域。三、实验数据某地区栅格数据DEM，数据来源于随书光盘
5、（Chp9Ex2）。四、实验要求根据DEM利用水文分析工具提取地表水流径流模型的水流方向、汇流累积量、水流长度、河流网络（包括河流网络的分级等）以及对研究区的流域进行分割等。五、实验流程图六、实验内容及步骤1无洼地DEM生成DEM 是比较光滑的地形表面模型，但由于DEM 误差以及一些真实地形或特殊地形的影响，使得DEM 表面存在一些凹陷的区域。在进行水流方向计算时，由于这些区域的存在，往往得到不合理的甚至错误的水流方向。因此，在进行水流方向的计算之前，应该首先对原始DEM 数据进行洼地填充，得到无洼地的DEM。洼地填充的基本过程是先利用水流方向数据计算出 DEM 数据中的洼地区域，并计算洼地深
6、度，然后，依据这些洼地深度设定填充阈值进行洼地填充。11 水流方向的提取水流的流向是通过计算中心格网与邻域格网的最大距离权落差来确定。对于每一格网的水流方向指水流离开此网格的指向。在ARCGIS 中，通过对中心栅格的1、2、4、8、16、32、64、128 等8个邻域栅格编码，中心栅格的水流方向便可有其中的某一值来确定。例如，若中心栅格的水流流向左边，则水流方向赋值16。流向的生成是个自动的过程，可能要等一段自时间，运算的时间跟电脑性能和DEM图的精度与大小有关.。方法是利用ArcToolboxSpatial Analysis Tools Hydrology Flow Direction，生成
7、方向水流流向图：若从DEM中作出来的流向分析的最大数值为128则不需要填洼，否则需要填挖。计算出来的水流分析12 洼地计算洼地区域是水流方向不合理的地方，可以通过水流方向来判断哪些地方是洼地，并进行填充。但是，并非所有的洼地区域都是由于数据的误差造成的，有很多洼地是地表形态的真实反映。因此在进行洼地填充之前，必须计算洼地深度，判断哪些地区是由于数据误差造成的，而哪些地区又是真实的地表形态。然后，在洼地填充时，设置合理的填充阈值。基本过程先分别双击 Hydrology 工具集中的Sink、Watershed 工具计算出洼地区域图、洼地贡献区域图，打开Spatial Analyst Tools 工
8、具箱中的Zonal 工具集，分别利用Zonal Statistic、Zonal Fill 工具分别计算每个洼地所形成的贡献区域的最低高程，及计算每个洼地贡献区域出口的最低高程，然后在ARCMAP 中加载Spatial Analyst 模块，选择RasterCalculator 对话框，通过命令sinkdep=（zonalmax-zonalmin）计算出洼地深度图。计算出来的洼地区域计算出来的洼地贡献区域洼地深度图13 洼地填充经过洼地提取后，可以确定原始 DEM 上是否存在洼地，若有洼地，须进行填充。而洼地深度的计算为填充阈值的设置提供了依据，系统默认条件下是不设阈值，即所有的洼地区域都将被填
9、平。对照地形资料，确定填充域值。没有岩溶问题的话，直接用（洼地深度的最大值1）作为域值（z limit）。洼地填充是一个不断反复的过程，直到所有的洼地被填平，新的洼地不在产生为止。方法是利用ArcToolboxSpatial Analysis Tools Hydrology Fill工具对洼地是深度图进行处理，得到洼地填充图：洼地填充图2汇流累积量计算在地表径流模拟过程中，汇流累积量是基于水流方向数据计算得到的。首先基于无洼地DEM 生成水流方向图，利用该数据，双击Hydrology 工具集中的Flow Accumulation 工具计算出汇流累积量数据，如图所示：通过计算得到的汇流累积量数据
10、3水流长度计算水流长度指地面上一点沿水流方向到流向起点（或终点）间的最大地面距离在水平面上的投影长度。它分为顺流计算及溯流计算两种，可通过双击Hydrology 工具集中的FlowLength 工具实现，其中计算方向分别选择顺流计算Downstream 或溯流计算Upstream。水流长度的提取和分析在水文学或水土保持工作中均具有很重要意义，因为水流长度直接影响地面径流的速度，进而影响地面土壤的侵蚀力。顺流计算出的水流长度溯流计算的水流长度4河网的提取目前常用的河网提取方法是采用地表径流漫流模型计算：首先是在无洼地DEM上利用最大坡降的方法得到每一个栅格的水流方向；然后利用水流方向栅格数据计算
11、出每一个栅格在水流方向上累积的栅格数，即汇流累积量，所得到的汇流累积量则代表在一个栅格位置上有多少个栅格的水流方向流经该栅格；假设每一个栅格处携带一份水流，那么栅格的汇流累积量则代表着该栅格的水流量。基于上述思想，当汇流量达到一定值的时候，就会产生地表水流，那么所有那些汇流量大于那个临界数值的栅格就是潜在的水流路径，由这些水流路径构成的网络，就是河网。41 河网的生成河网的生成基于汇流累积量数据，利用Map Algebra 工具集中的Multi Map Output 工具中的Con 命令或Setnull进行有条件地查询可以得到，还可以利用栅格计算器获取。setnull，将栅格图像中满足条件的栅
12、格值设置为null，即nodata，这样很便于在栅格图像中删除不需要的部分。其对应Grid模块命令为setnull，格式：outgrid = setnull 5, ingrid,表示，如果ingrid5，则变为nodata，否则返回为原值。con表示条件，如果满足条件，返回为指定值1，否者返回为指定值2，Grid命令：outgrid = con 5, 10 ,100,表示，如果ingrid值5，则对应栅格值改为10，否者改为100。若100不指定，则返回为nodata。这里可以指定多个条件。通过计算生成的栅格河网42 栅格河网矢量化利用ArcToolboxSpatial Analyst Too
13、lsHydrologystream to feature工具箱对栅格河网进行处理，如图：矢量河网平滑处理河流网络：利用Editor工具栏对矢量河网进行编辑，单击Editor下拉菜单，单击Staring Editor 激活编辑工具栏，在下拉菜单中选中More Editing Tools Advanced Editing,选中整幅地图，在弹出的Advanced Editing工具栏中点击平滑（Smooth）按钮 ，在弹出的Smooth对话框中Maximum allowable offset栏中输入值“3”，表示允许最大偏移量为3。从而得到平滑处理后的河流网络。用这种方法所平滑的流域河流网络偏移量较
14、大，可能会产生人为误差，也可利用Data Management Tools工具箱下Data Management ToolsGeneralizationSmooth Line工具箱对河流网络进行平滑处理。此种方法偏移较小。43 Stream Link 的生成Stream Link 记录河网中结点之间的连接信息，它主要是记录河网的结构信息，其中每一条弧段连接着两个作为出水点或汇合点的结点。Stream Link 的生成可以通过双击Hydrology 工具集中的Stream Link 工具实现，其属性表中记录着每个片段所包含的栅格个数信息。如图：Stream Link 示意图44 河网分级不同级别
15、的河网所代表的汇流累积量不同，级别越高，汇流累积量越大，一般是主流，而级别较低的河网一般则是支流。在ARCGIS的水文分析中，双击Hydrology 工具集的Stream Order 工具，可提供两种常用的河网分级方法：Strahler 分级和Shreve 分级。如图：河网的Strahler分级结果河网的Shreve分级结果5流域分割流域（watershed）又称集水区域，是指流经其中的水流和其他物质从一个公共的出水口排出从而形成的一个集中的排水区域。流域可以通过流域盆地（basin）、集水盆地（catchment）来描述。51 流域盆地的确定流域盆地是由分水岭分割而成的汇水区域，可利用水流方向确定出所相互连接并处于同一流域盆地的栅格区域。双击Hydrology 工具集中的Basin 工具可以计算出流域盆地图：计算出的流域盆地（线状图形为矢量河网数据）52 积水流域的确定除用流域盆地来描述外，在水文分析中，经常基于更小的流域单元进行分析，首先通过Hydrology 工具集中的Snap Pour Point 工具寻找小级别流域的出水口位置，然后结合水流方向，双击Hydrology 工具集中的Watershed 工具，分析搜索出该出水点上游所有流过该出水口的栅格，
