GIS学习二:地形分析

1. TIN 及DEM 生成

1.1由高程点、等高线矢量数据生成TIN转为DEM

在ArcMap中新建一个地图文档。添加矢量数据:Elevpt_Clip、Elev_Clip、Boundary、Erhai。设置TIN的图层(符号)得到如下的效果。
在这里插入图片描述

高程+erhai
在这里插入图片描述

执行工具栏TIN转换到栅格,指定相关参数:属性:[高程],像素大小:[50],输出栅格的位置和名称: [TinGrid]。确定后得到DEM数据:TinGrid, 其中,每个栅格单元表示50m×50m的区域
在这里插入图片描述

1.2 TIN的显示及应用

关闭除[TIN]之外的所有图层的显示,编辑图层[tin]的属性,在图层属性对话框中,点击[符号] 选项页,将[ 边界类型 ] 和 [ 高程 ] 前面检查框中的勾去掉; 点击 [ 添加 ] 按钮
在 [添加渲染] 对话框中,将 [所有边用同一符号进行渲染] 和 [ 所有点用同一符号进行渲染 ] 这两项添加到TIN的显示列表中
在这里插入图片描述

将TIN图层局部放大
在这里插入图片描述

TIN 转换为 坡度多边形
新建地图文档,加载图层 [tin],参考上一步操作,将 [面坡度用颜色梯度表进行渲染] 和 [面坡向用颜色梯度进行渲染] 这两项添加到TIN的显示列表中,
在上面的对话框中,选中Slope,点击 [分类] 按钮,在下面的对框中,将[类] 指定为 5,然后在 [间隔值] 列表中输入间隔值:[ 8, 15,25, 35, 90]
点击两次 [确定] 后关闭图层属性对话框,图层[ tin ] 将根据指定的渲染方式进行渲染,效果如下图所示 :
在这里插入图片描述

执行[3D分析]工具栏中的命令[转换]>>[TIN转换到矢量]
得到多边形形图层:[ tinSlopef] ,它表示研究区内各类坡度的分布状况,结果是矢量格式,打开其属性表可以看到属性 [SlopeCode] 为 数值[1,2,3,4,5]
在这里插入图片描述

查看矢量图层:tinSlopef 中要素属性表,其中属性[SlopeCode]1,2,3,4,5分别表示坡度范围(0-8)、(8-15)、(15-25)、(25-35)、(>35)
在这里插入图片描述

Eliminate合并破碎多边形(选做,需要8-10分钟)
新建地图文档,加载坡度多边形图层:TinSlopef, 打开TinSlopef的属性表,添加一个字段Area(类型为Double),通过[计算值]操作,计算各个多边形的面积:
选中高级,输入VBA代码到[Pre-Logic VBA Script Code],输入变量[dblArea]到[Area=] 下的输入框中。[Pre-Logic VBA Script Code]框中代码如下:
dim dblArea as double
dim pArea as IArea
set pArea = [Shape]
dblArea = pArea.Area
以下的操作将会把面积小于10000平方米的多边形合并到周围与之有最长公共边的多边形中:
在这里插入图片描述

执行菜单命令[选择]>>[通过属性选择],查询”Area”<=10000 (平方米)的图斑

被选中的多边形以高亮方式显示,这些小的图斑将会被合并到与之相邻且有最大公共边的多边形。当然也可以选择合并到相邻的面积最大的多边形。
在这里插入图片描述

Eliminate(合并破碎多边形)操作原理
在这里插入图片描述

原始多边形 合并后多边形,选中的(面积<=10000m2 )多边形被合并到与之相邻的面积最大的多边形中

将地图适当放大,比较原始图层:tinSlopef 与合并后的图层:tinSlopef_Eliminate
在这里插入图片描述

TIN 转换为坡向多边形
得到坡向多边形图层后,得到的坡向多边形中属性AspectCode的数值(-1,1,2,3,4,5,6,7,8,9)分别表示当前图斑的坡向(平坦、北、东北、东、东南、南、西南、西、西北、北),其中1,9是相同的可以合并为1
在这里插入图片描述

2. DEM的应用

2.1坡度:Slope

新建地图文档,加载[1.1(6)]中得到的DEM数据:TINGrid。加载3D分析扩展模块,打开[3D分析]工具栏,执行菜单命令[3D分析]>>[表面分析]>>[坡度], 指定各参数。得到坡度栅格slope of TinGrid:
在这里插入图片描述

坡度栅格中,栅格单元的值在[ 0 -90 ] 度间变化
在这里插入图片描述

以下计算剖面曲率:
执行菜单命令:[3D分析]>>[表面分析]>>[坡度]。按如下所示,指定各参数:得到剖面曲率栅格:[Slope of Slope of tingrid]
在这里插入图片描述

2.2 坡向:Aspect

在上一步的基础上进行,关闭[Slope of tingrid]的显示。执行菜单命令:[3D分析]>>[表面分析]>>[坡向],按下图所示,指定各参数:得到坡向栅格:[Aspect of tingrid]
在这里插入图片描述

坡向栅格
以下计算平面曲率:执行菜单命令:[3D分析]>>[表面分析]>>[坡度],按下图所示指定各参数:生成平面曲率栅格:[Slope of Aspect of tingrid]:
在这里插入图片描述

2.3提取等高线

新建地图文档,加载DEM数据: [tingrid]。〔在执行以下操作时确保,3D分析扩展模块已激活〕,打开Arctoolbox,执行命令: [3D Analyst Tools]>>[ Raster Surface ]>> [ 等高线 ] ,按上图所示指定各参数。生成等高线矢量图层:Contour_tingrid:
在这里插入图片描述

2.4计算地形表面的阴影图

在上一步基础上进行,打开[3D 分析]工具栏,执行菜单命令:[ 3D分析 ]>>[表面分析]>>[ 山影 ],按下图所示指定各参数:生成地表阴影栅格:[ Hillshade of tinGrid ]:
在这里插入图片描述

DEM渲染:
如以下第2幅图所示,关闭除[tingrid] 和 [Hillshade of tingrid]以外所有图层的显示,并将[ tingrid ] 置于[ Hillshade of tirngrid] 之上,右键点击[ tingrid] ,在出现的右键菜单中执行[ 属性 ],在[图层属性]对话框中,参照下图所示设置[符号]选项页中颜色。
打开工具栏[效果],如下图所示,设置栅格图层[tingrid]的透明度为:[40%]左右。
在这里插入图片描述

2.5可视性分析

A.通视性分析
在上一步的基础上进行,打开[ 3D分析] 工具栏,从工具栏选择[ 通视线 ](Line of sight)工具:在出现的[ 通视线 ]Line of Sight对话框中输入[观察者偏移量] 和 [目标偏移量], 即距地面的距离,如图:
在地图显示区中从某点[A]沿不同方向绘制多条直线,可以得到观察点 [A] 到不同目标点 的通视性:
在这里插入图片描述

绿色线段表示可视的部分,红色线段表示不可见部分

B.可视区分析:移动发射基站信号覆盖分析
在上一步基础上进行,在内容列表区[TOC]中关闭除 [tingrid] 之外的所有图层,加载移动基站数据-矢量图层:[移动基站.shp]
在[3D 分析] 工具栏中,执行菜单命令:[3D 分析]>>[表面分析]>>[视域],按下图所示指定各参数:
生成可视区栅格:[ ViewShed of 移动基站 ]:
在这里插入图片描述

其中绿色表示现有发射基站信号已覆盖的区域,淡红色表示,无法接收到手机信号的区域

2.6地形剖面

在上一步基础上进行,打开 [ 3D分析 ] 工具栏,点击 [插入线] 工具,跟踪一条线段,这条线段可以从DEM:[TINGRID] 中得到高程值,
点击[ 创建剖面图 ] 按钮,得到上一步所生成的3D线段的剖面图:
在这里插入图片描述
额…都是自己写的哈,小学一下

  • 0
    点赞
  • 19
    收藏
    觉得还不错? 一键收藏
  • 5
    评论
评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值