DEM分析报告报告材料与可视化 下载本文

内容发布更新时间 : 2024/11/16 14:53:45星期一 下面是文章的全部内容请认真阅读。

实用标准文案

一.软件平台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对话框,在对话框中设置目

精彩文档