地形分析实验报告.docx

实验项目（三）:地形分析一.实验目的：掌握基本地形因子计算，地形特征提取的基本方法和步骤，并进行不同类型地形因子的统计。二.实验器材：arcmap三.实验数据：dem数据四.实验步骤：1 .坡度的提取，不同类型区的坡度提取：D空间分析模块中的坡度提取：加载dem数据，和空间分析模块；SpatialAnalyst->surfaceanalysis->slope,打开相应的对话框；通过设置不同的输出像元大小，可以得到不同精度的坡度提取图。2）三维分析模块中的坡度提取：3DAnalyst->->surfaceanalysis->slope,打开相应的对话框；通过设置不同的输出像元大小，可以得到不同精度的坡度提取图。2 .坡向的提取，不同类型区的坡向提取：D空间分析模块中的坡向提取：加载dem数据，和空间分析模块；SpatialAnalyst->surfaceanalysis->asect,打开相应的对话框；通过设置不同的输出像元大小，可以得到不同精度的坡度提取图。2）三维分析模块中的坡向提取：3DAnalyst->->surfaceanalysis->aspect,打开相应的对话框；通过设置不同的输出像元大小，可以得到不同精度的坡度提取图.3 .山脊线,山谷线的提取：1）加载原始DEM数据。ArcToolbox,/AnalysisTools/Hydrology<.2）洼地填充：Hydrology/Fi11,进行原始DEM的洼地填充。在Inputsurfaceraster文本框中选择原始DEM数据dem,将输出数据命名为filldem,Zlimit取默认值，即填充所有洼地。3）无洼地的水流方向的计算：Hydrology：FlowDirection,在Inputsurfaceraster文本框中选择填充过的无洼地DEM数据filldem,设置输出文件名为flowdirfillo4）汇流累积量的计算：Hydrology/FlowAccumulationo选择flowdirfill作为输入的水流方向数据，输出数据命名为fIowacch5）汇流累积量为零值的提取：SpatialAnalyst/RasterCalculator,在文本框中填写公式：facc=（fIowacc=O）,单击evaluate进行计算。6）在ArCMaP中打开facc,会发现有很多的地方并不是山脊线，因此需对此数据做如下处理：利用邻域分析的方法（NeighborhoodStatiStiCS命令），对face。进行3X3邻域分析，求均值（mean）,使数据变得光滑，处理后的数据命名为neiborfaccO。7）SpatialAnalyst/SurferAnalysis:COntour和HiIIShade命令，分别生成原始DEM的等值线图ctour和晕渲图hi11shade。8）在neiborfaccO数据上单击右键，选择Properties命令，进入Symbology选项卡，进行重新分级，将数据分为两级，这时需要不断调整分级临界点，并以等值线图和晕渲图作为辅助判断。属性值越接近于1的栅格越有可能是山脊线的位置，最终确定的分界阈值为0.5541。9）将进行过二值化的neiborfaccO进行重分类为reneibor,将属性值接近1的那一类的属性值赋值为L其余的赋值为0。10）利用SpatialAnalyst菜单下的RasterCalculator将重分类后的reneibor数据与正地形数据Zhengdixing相乘，消除了那些存在于负地形区域中的错误的山脊线。然后将计算结果进行重分类，所有属性不为1的栅格值赋为NODATA,即得到了山脊线。4、计算表面阴影D选择用来计算阴影的表面2）设置太阳高度角与方位角;3）设定高程转换系数；4）指定输出栅格单元大小5）指定输出路径及文件名。