内容发布更新时间 : 2024/7/5 5:06:06星期一 下面是文章的全部内容请认真阅读。

实用标准文案

一.软件平台ArcGIS或MapGIS（软件测试部分）：

（1）数据处理：拓扑构建、误差校正、地图投影

（2）数据管理：属性表创建、属性表关联、图形与属性数据挂接、属性表导出

（3）空间分析：查询检索、叠加分析、缓冲区分析

（4）数字高程模型：GRID及TIN模型创建，DEM分析（包括坡度、坡向、粗糙度、可视性、洪水淹没、流域地貌等分析）

（5）数据转换：ArcGIS、MapGIS、MapInfo、AutoCAD等数据间格式转换

实验四 基于ArcGIS的DEM分析与可视化 一、实验目的

1、掌握利用ArcGIS三维分析模块进行创建表面的基本方法 2、掌握地形特征信息的提取方法，能利用ArcGIS软件基于DEM对山脊线和山谷线的提取，显示粗糙度

3、掌握三维场景中表面及矢量要素的立体显示其原理与方法，熟练掌握ArcGIS软件表面及矢量要素杂场景中的三维显示及其叠加显示

4、熟练掌握ArcScene三维场景中要素、表面的多种可视化方法。 二、主要实验器材（软硬件、实验数据等）

计算机硬件：性能较高的PC；计算机软件：ArcGIS9.3软件；实验数据：《ArcGIS地理信息系统空间分析实验教程》随书光盘或其他中

三、实验内容与要求 1、地形特征信息提取 实验数据：dem

要求：利用所给区域DEM数据，提取该区域山脊线、山谷线栅格数据层。 具体操作： 1.打开arcmap，添加dem数据，点击DEM数据，打开Arctoolbox，使用Spatial Analysis tools\\Surface Analysis\\Aspect工具，提取DEM的坡向数据层，命名为A。

2.点击数据层A，使用Spatial Analysis tools\\Surface Analysis\\Slope工具，提取数据层A的坡度数据，命名为SOA1。（地面坡向变率，是指在地表的坡向提取基础之上，进行对坡向变化率值的二次提取，亦即坡向之坡度（Slope of Aspect, SOA）。它可以很好的反映等高线弯曲程度。）

3.求取原始DEM数据层的最大高程值，记为H；使用空间分析工具集中的栅格计算器（Raster Calculator），公式为（H—DEM），得到与原来地形相反的数据层，即反地形DEM数据。记为“-DEM”。

4.基于“-DEM”数据，使用Spatial Analysis tools\\Surface Analysis\\Aspect工具，提取-DEM的坡向数据层，命名为-A。。

5. 点击数据层-A，使用Spatial Analysis tools\\Surface Analysis\\Slope工具，提取反地形的坡向变率，记为SOA2。

6.使用空间分析工具集中的栅格计算器（Raster Calculator），公式为SOA=(([SOA1]+[SOA2])-Abs([SOA1]+[SOA2]))/2，这样就可以求出没有误差的DEM的坡向变

精彩文档

实用标准文案

率SOA。

7.再次点击初始DEM数据，点击Spatial Analyst---->Neighborhood Statistics（栅格邻域计算工具），打开对话框，设置统计类型（Statistic type）为平均值（Mean），邻域的类型为矩形（也可以为圆），领域的大小为11×11（这个值也可以根据自己的需要进行改变），则可得到一个邻域为11×11的矩形的平均值数据层，记为B，各项设置如图4-1所示。

8.使用空间分析工具集中的栅格计算器（Raster Calculator），公式为C=[DEM]-[B]，即可求出正负地形分布区域，如图4-2所示。

图4-1 栅格邻域计算对话框 图4-2正负地形分别示意图 9.使用空间分析工具集中的栅格计算器（Raster Calculator），公式为shanji=[C]>0&SOA>70，即可求出山脊线，如图4-3所示。

10.同理，在栅格计算器（Raster Calculator）中，键入公式为shangu=[C]<0&SOA>70，即可求出山谷线，如图4-4所示。

精彩文档

实用标准文案

图4-3 山脊线的界面 图4-4 山谷线的界面

11.地面粗糙度：

地面粗糙度是指在一个特定的区域内，地球表面积与其投影面积之比。它也是反映地表形态的一个宏观指标。根据地面粗糙度的定义，求每个栅格单元的表面积与其投影面积之比，

可以用如下方法来完成。

如图 8，假如 ABC 是一个栅格单元的纵剖面，α为此栅格单元的坡度，则 AB 面的面积为

此栅格的表面积，AC 面为此栅格的投影面积（也既是此栅格的面积），根据公式：Cosα=AC/AB

则可得出此栅格单元的地面粗糙度 M 为：M =“AB 面的面积”/“AC 栅格单元的面积”= （AC * AB ）/ （AC * AC） = 1 / Cosα

地面粗糙度的提取步骤如下：

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、表面创建及景观图制作

实验数据具体见《ArcGIS地理信息系统空间分析实验教程》随书光盘\\ch9\\EX3。 要求：

（1）利用所给等高线数据建立景区栅格表面；

（2）在ArcScene三维场景中，实现表面与其它要素叠加三维显示； （3）设计各要素如道路、水系等的符号化显示；

（4）综合考虑表面及各要素，生成美观大方的区域景观图； 具体操作：

思路是利用原有图层创建区域TIN表面，然后建立地表的平面景观，在生成三维景观图。 第一，创建区域TIN表面

1、 启动ArcScence,打开数据,在3D Analyst工具栏中选择Create/Modify TIN中的Create TIN From Features...打开Create TIN From Features对话框，在对话框中设置目

精彩文档