地理信息系统课程报告 林力.doc

《地理信息系统》课程设计报告 学号 1106601-30 姓名 林振坤 设计时间 2013年12月23日——2014年1月3日 设计地点 湖南城市学院信息楼406机房 设计的目的 理解基于DEM结合水文分析的方法提取研究区域和山谷线的原理；掌握水流方向.汇流累积量的提取方法以及他们的提取原理；能将水文分析的方法和其他的空间分析方法结合以及应用问题。 设计小组成员 1106601-24田玉姣 1106601-30林振坤 1106601-32阳星 1106601-31陈敏 指导教师 刘雯 设计课题 利用水文分析方法提取山脊、山谷线 基本思想及关键问题得解决方法 对于水文物理过程研究而言，由于山脊，山谷分别表示分水性与汇水性，山脊线和山谷线的提取实质上也是分水线与汇水线的提取。因此，对于山脊线和山谷线就可以利用水文分析的方法进行提取。 基于DEM的这种地形表面流水物理模拟分析的原理是：对于山脊线而言，由于它同时也是分水线，那么对于分水线上的那些栅格，由于分水线的性质是水流的起源点，通过地表径流模拟计算值有的这些山歌的水流方向都应该只具有流出方向而不存在流入方向，也就是其山歌的汇流累积量为零。通过对零值的汇流累积值的栅格的提取，就可以得到分水线，也就得到了山脊线；对于山谷线而言，由于其具有汇水的性质，那么对于山谷线的提取，可以利用反地形的特点，即是利用一个较大的数据减去原始的DEM数据，而得到了与原始地形完全相反的地形数据，也就是原始的DEM中的山脊变成负地形的山谷，而原始DEM中的山谷在负地形中就变成了山脊，那么，山谷线的提取就可以在负地形中利用提取山脊现代方法进行提取。 技术路线图 实验步骤 正负地形的提取 启动ArcToolbox，展开Analysis Tools 工具箱，打开hydrology 工具集。在图层管理器中加在研究区域的原始DEM数据。 加载Spatial Analyst 模块，点击Spatial Analyst 模块的下拉箭头，点击neighborhood statistics 菜单工具，利用领域分析的方法以11*11的窗口计算平均值。分析结构命名为meandem 。 点击spatial analyst 中的raster calculator 菜单工具，对原始DEM数据与领域分析之后的数据meandem 做减法运算，并将运算结果重分为两级，分级界线为0，那么大于0的区域在原始DEM上就是正地形区域，小于0的区域在原始DEM上就使负地形区域。 对于上一部得到的二值化数据惊醒两次重分类，一次将正地形区域属性值为1，负地形区域属性赋值为0，命名为zhengdixing；另一次将正地形区域属性赋值为0，负地形区域属性赋值为1，命名为fudixing。 山脊线的提取 在ArcMap 中加载研究区域的原始DEM数据。 洼地填充：双击hudrology 工具集中的fill工具，进行原始DEM的洼地点进行填充。在Input surface raster文本框中选择原始的DEM数据dem，将输出数据命名为filldem,因为选择的是将所有洼地全部填充，所有在填充容限Z limit 为默认值。 基于无洼地的水流方向的计算：双击hydrology 工具集中的flow direction 工具，在Input surface raster 文本框中选择填充过的无洼地DEM数据filldem，将输出的水流方向数据命名为flowdirfill。 汇流累积量的计算：双击hydrology 工具集中的flow accumulation 工具。悬着flowditfill作为输入的水流方向数据；输出数据命名为flowaccl。 回流累积量零值的提取：加载Spatial Analyst 模块，点击Spatial Analyst 模块的下拉箭头，然后单击raster calculator 菜单，打开栅格计算对话框，在文本框中填写汇流累积量零值的提取公式:facc0=(flowacc=0),然后点击evaluate 进行计算。计算结果为所有的汇流累积量为0的栅格。 在ArcMap中打开facc0，发现所有提取出的栅格很乱，有很多的地方并不是山脊线的位置，因此应对这个数据进行处理。处理过程可以利用领域分析的方法，对提取出的汇流累积量等于零值的数据进行3*3领域分析进行光滑处理，处理后的数据命名为neiborfacc0。 单击spatial analyst模块中的surfer analyst 中的countline和hillshade菜单命令，分别生成原始DEM的等值线图ctour和晕渲图hillshade。 打开neiborfacc0数据的属性信息，进行重新分类，将分类级别设置为两类，不断的调整分界数据大小，并以由DEM生成的等值线