实验七、DEM的生成及应用

  • 格式:doc
  • 大小:820.00 KB
  • 文档页数:7

《地理信息系统软件应用》实验报告

实验项目名称:DEM的生成及应用

一、实验目的

DEM是对地形地貌的一种离散的数字表达,是对地面特性进行空间描述的一种数字方法、途径,它的应用可遍及整个地学领域。通过对本次实习的学习,我们应:

a) 加深对TIN建立过程的原理、方法的认识;

b) 熟练掌握ArcGIS中建立DEM、TIN的技术方法。

c) 掌握根据DEM或TIN 计算坡度、坡向的方法。

d) 结合实际,掌握应用DEM解决地学空间分析问题的能力。

二、实验内容、实验要求及相应理论基础

要求:根据上述所给数据完成TIN及DEM的生成及应用

相关理论:在ArcGIS中支持的3D数据包括3类:栅格表面数据(Raster surfaces)、TIN表面数据(TIN surfaces)、含有Z值的特征要素(3D Feature)。

在ArcGIS中TIN的建立是基于一系列含有X、Y、Z值坐标的点,这些点就构成了TIN中三角形的顶点(Nodes),连接这些节点间的线构成了TIN中三角形的边(Edges)。TIN一旦建立,TIN表面上的任何一个点的高程值都可以通过数学方法估算出来或是利用该点所在三角形三顶点的高程值内插得到,TIN中的每个面的坡度、坡向等也都可以计算出来。

将TIN转换成Raster的目的一般有:⑴表面建模的需要;⑵从TIN中提取出表面的坡度、坡向以及高程信息。需要将Raster转换成TIN的情况有:(1)表面建模的需要;(2)由于可视化的需要,简化表面模型

三、使用仪器、材料

(1)软件准备:ArcGIS Desktop 9.3--ArcMap(3D分析模块)

(2)实验数据:top point.shp、river arc.shp、contour.shp, rivepoly polygon.shp、border

polygon.shp

四、实验步骤

(一):由高程点、等高线矢量数据生成TIN转为DEM

(1) 添加矢量数据top point.shp、river arc.shp、contour.shp, rivepoly polygon.shp、border polygon.shp(同时选中:在点击的同时按住Shift)

(2) 激活“3D Analyst”扩展模块(执行菜单命令 [工具]>>[扩展],在出现的对话框中选中3D分析模块),在工具栏空白区域点右键打开[3D分析] 工具栏

(3)执行工具栏[3D分析]中的菜单命令[3D分析]>>[创建/修改TIN]>>[从要素生成TIN];

(4)在对话框[从要素生成TIN中]中定义每个图层的数据使用方式;

在[从要素生成TIN中]对话框中,在需要参与构造TIN的图层名称前的检查框上打上勾,指定每个图层中的一个字段作为高度源(Height Source),设定三角网特征输入(Input as)方式。可以选定某一个值的字段作为属性信息(可以为None)。在这里指定图层[counter] 的参数:high[三角网作为:]指定为hard line,其它图层参数使用默认值即可。

(5)确定生成文件的名称及其路径,生成新的图层tin,在TOC(内容列表)中关闭除[TIN]和[counter]之外的其它图层的显示,设置TIN的图层(符号)得到如下的效果。

(6)执行工具栏[3D分析]中的命令[转换]>>[TIN转换到栅格],指定相关参数:属性:[高程],像素大小:保留默认值4,输出栅格的位置和名称: [TinGrid]

确定后得到DEM数据:TinGrid, 其中,每个栅格单元表示47.12m×47.12m的区域

(二):TIN的显示及应用

(1) 在上一步操作的基础上进行,关闭除[TIN]之外的所有图层的显示,编辑图层[tin]的属性,选择tin图层右击打开“properties”在图层属性对话框中,点击[符号] 选项页,将[ 边界类型 ]

和 [ 高程 ] 前面检查框中的勾去掉; 点击 [ 添加 ] 按钮

(2) 在 [添加渲染] 对话框中,将 [所有边用同一符号进行渲染] 和 [ 所有点用同一符号进行渲染 ] 这两项添加么TIN的显示列表中.

(3) 按上述操作后将TIN图层局部放大,认真理解TIN的存储模式及显示方式,得到如下简图:

(4) TIN 转换为 坡度多边形

新建地图文档,加载图层 [tin],参考上一步操作,将 [面坡度用颜色梯度表进行渲染]

和 [面坡向用颜色梯度进行渲染] 这两项添加到TIN的显示列表中,然后选择确定。

1、 在上面的对话框中,选中Slope,点击 [分类] 按钮,在下面的对框中,将[类] 指定为 5,然后在 [间隔值] 列表中输入间隔值: [ 8, 18,27, 47, 90] ,如下图所示

2、 点击两次 [确定] 后关闭图层属性对话框,图层[ tin ] 将根据指定的渲染方式进行渲染,效果如下图所示 :

3、 执行[3D分析]工具栏中的命令[转换]>>[TIN转换到矢量],按下图所示指定各参数:

4、

得到多边形形图层:[ tinSlope] ,它表示研究区内各类坡度的分布状况,结果是矢量格式,打开其属性表可以看到属性 [SlopeCode] 为 数值[1,2,3,4,5]

4、查看矢量图层:tinSlopef 中要素属性表,其中属性[SlopeCode]1,2,3,4,5分别表示坡度范围(0-8)、(8-17)、(17-27、(27-47)、(>475)

(5)Eliminate合并破碎多边形

1、计算多边形面积

A:新建地图文档,加载坡度多边形图层:TinSlope 打开TinSlopef的属性表,添加一个字段Area(类型为Double).

B:选择Area字段,右击打开选择File cacular,打开对话框后把其中的高级设置勾选上。

C:右击“help”,选择面积的计算方式:To calculate area:(Dim Output as double、Dim pArea

as Iarea、Set pArea = [shape]、Output = pArea.area),

D;输入参数后单击确定即计算面积,如下图所示:

2、选择面积小于8000平方米的多边形并把其合并到周围面积最大的多边形中:

A:右击打开tinslope图层的属性表,选择Option.右击选择 “select by attribute”

B:执行菜单命令[选择]>>[通过属性选择],查询”Area”<=8000 (平方米)的图,建立关系式“Area<=8000”,单击OK。

C:选中的区域属性表和图中对应的部分都显示成蓝色,如下图所示:

被选中的多边形以高亮方式显示,这些小的图斑将会被合并到与之相邻且面积最大的多边形

3、打开“Arctoolbox”——“Data management tools——“Genenalization”——“Eliminate”

然后择输入图层为tinSlope,指定输出路径和名字tinSlope Eliminate,即可将面积<=8000平方米的多边形合并到大地多边形中。

(6)TIN 转换为坡向多边形,方法参照以上第(4)步,得到坡向多边形图层得到如下图:

得到的坡向多边形中属性AspectCode的数值(-1,1,2,3,4,5,6,7,8,9)分别表示当前图斑的坡向(平坦、北、东北、东、东南、南、西南、西、西北、北),其中1,9是相同的可以合并为1.

五、实验过程原始记录(数据、图表、计算等)

(1)计算tinslope图层的面积(截取部分)

(2)查找面积小于8000平方米的多边形:建立关系式“Area<=8000”

(3)更多图表见制图成果。

六、实验结果及分析

第一次做这个实验的时候,在由高程点、等高线矢量数据生成TIN的时候,因为没能弄清建模的实质,在操作过程中系统总是报错,而且建模时所选择的字段不是高程获只有部分高程被选上,导致做出来的实验效果远不是预期所期待的那样,也导致在后续步骤无法进行。后来我通过上网查资料才弄清建模原理,再次实验时才做成功。

七、实验的心得体会

通过这次实验我更加熟悉了Arcgis的相关操作,也了解了部分3D分析功能的使用,在由高程点、等高线矢量数据生成TIN时,我发现找对高程点和高程来源是关键,否则TIN无法建成,或者出错。栅格转 TIN 工具的用途是创建表面偏离输入栅格不超过指定 Z 容差的不规则三角网

(TIN)。栅格转 TIN 常用于将通过数字高程模型 (DEM) 获得的栅格转换为 TIN 表面模型。