DEM上机0.docVIP

数字高程模型上机实验报告 专业年级：08测绘工程 姓名：李旭 学号：0808030101 1. 地形指标提取 （1）坡度变率 1) 激活 DEM 主题，选择Spatial Analysis － Surface Analysis － Slope 命令，提取DEM 主题的坡度，得到主题Slope of DEM。 2) 激活主题Slope of DEM，再对其用上述的方法提取坡度，得到DEM 主题坡度的坡度，即坡度变率主题。 （2）坡向变率 1) 求取原始DEM数据层的最大高程值，记为H；通过Spatial Analysis 下的栅格计算器Calculator，公式为（H－DEM），得到与原来地形相反的DEM数据层，即反地形DEM 数据； 2) 基于反地形DEM 数据求算坡向值； 3) 利用 SOA 方法求算反地形的坡向变率，记为SOA2，由原始DEM 数据求算出的坡向变率值为SOA1； 4) 在 Spatial Analysis 下使用栅格计算器Calculator，公式为SOA ＝（（[SOA1]+[SOA2]）-Abs（[SOA1]-[SOA2]））/ 2，即可求出没有误差的DEM 的坡向变率。 （3）地形起伏度 1) 激活 DEM 数据，在Spatial Analysis 下使用栅格邻域计算工具Neighborhood Statistics。设置Statistic type 为最大值，邻域的类型为矩形（也可以为圆），邻域的大小为11×11（这个值也可以根据自己的需要进行改变），则可得到一个邻域为11×11 的矩形的最大值层面，记为A； 2) 重复 1、2，只是把Statistic type 值设置为最小值，即可得到DEM 数据的最小值层面，记为B； 3) 在 Spatial Analysis 下使用栅格计算器Calculator，公式为[A]-[B]，即可得到一个新层面，其每个栅格的值是以这个栅格为中心的确定邻域的地形起伏值。 （4）地面粗糙度： 1) 激活 DEM 主题，选择Spatial Analysis － Surface Analysis － Slope 命令，提取DEM 主题的坡度，得到主题Slope of DEM； 2) 激活主题Slope of DEM，在Spatial Analysis 下使用栅格计算器Calculator，公式为:1 / Cos（[Slope of DEM]*3.14159/180）即可得到地面粗糙度的层面。 2. 地形特征信息提取 1) 激活 DEM 数据，在Spatial Analysis 下使用surface 菜单下的Derive Aspect 命令，提取DEM 坡向层面记为A； 2) 激活 A 层面，在Spatial Analysis 下使用surface 菜单下的Derive Slope 命令，提取A 层面的坡度信息，记为SOA1； 3) 求取原始 DEM 数据层的最大高程值，记为H；通过Spatial Analysis 下的栅格计算器Calculator，公式为（H－DEM），得到与原来地形相反的DEM 数据层，即反地形DEM 数据； 4) 基于反地形DEM 数据求算坡向值； 5) 利用 SOA 方法求算反地形的坡向变率，记为SOA2； 6) 在 Spatial Analysis 下使用栅格计算器Calculator，公式为SOA ＝（（[SOA1]+[SOA2]）-Abs（[SOA1]-[SOA2]））/ 2，即可求出没有误差的DEM 的坡向变率SOA； 7) 激活原始DEM 数据，在Spatial Analysis 下使用栅格邻域计算工具Neighborhood Statistics；设置Statistic type 为平均值，邻域的类型为矩形（也可以为圆），邻域的大小为11×11（这个值也可以根据自己的需要进行改变），则可得到一个邻域为11×11 的矩形的平均值层面，记为B； 8) 在 Spatial Analysis 下使用栅格计算器Calculator，公式为C ＝[DEM]-[B]，即可求出正负地形分布区域。 9) 在 Spatial Analysis 下使用栅格计算器Calculator，公式为D ＝[C] 0 SOA 70,即可求出山脊线。 10) 同理，在栅格计算器 Calculator 中，修改公式为D ＝[C] 0 SOA 70,即可求出山谷线。 3. 表面创建及景观图制作 实验方法及其制作步骤： （1）启动ArcScene，打开场景文件Exercise3.sxd（光盘中..\Chp9\Ex3\Exercise3.sxd），其中已添加以下数据层：等高线数据层Arc_Clip，道路数据层Arc_Clip