概述
本文只是一篇实用教程,不是专门的技术文档,更不是什么学术论文,主要是便于初学者能看懂,本教程只适合ArcGIS基础薄弱但又需要在规划项目中通过ArcGIS做初步简单的地形处理和分析的童鞋,当然最好有一定ArcGIS基础,起码要知道ArcGIS是干嘛的吧,本文不负责介绍ArcGIS软件的理论知识。
本教程尽量结合一点经验,把可能遇到的问题或难点以及相应的解决办法也写清楚,也是自认为区别于网上其他一些教程的特点吧。
ArcGIS做地形分析的基本步骤其实是非常简单的,网上也有非常多的教程进行阐述,其实大家按这个步骤也基本能做出来,但问题就出现在每个人所使用的基础数据都多多少少有不同的问题,因此在操作过程中就不同程度的受阻,因此本教程先列出最基本的操作步骤,再对所使用的数据及各操作步骤可能出现的问题进行分析并逐一列出解决方法。
ArcGIS地形分析的基本步骤如下:
- 整理并添加带有高程属性的CAD等高线或者高程点;
- 将导入的等高线或高程点数据转换为ArcGIS可编辑的Shapefile文件;
- 进一步整理、检查、修正Shapefile文件中的数据错误;
- 使用整理好的Shapefile生产数字高程三角模型TIN文件;
- 将TIN文件转换为栅格数字高程模型DEM文件;
- 基于DEM即可进一步进行基于地形的坡度、坡向、起伏度、汇水特征等专项分析了。
ArcGIS地形分析的核心是数字高程模型DEM数据,基于地形的高程、坡度、坡向、起伏度、汇水分析等都是基于DEM数据的,所以本教程首先重点介绍如何生成获得DEM数据。
DEM的全称是Digital Elevation Model,非常直接的翻译,就是“数字”、“高程”、“模型”三个词的直接拼接,可谓是简单粗暴好理解,可到底什么是“数字高程模型”这个概念可能是经常困扰大家的问题,貌似很高大上,尤其是一些学术论文,经常只写其英文名,还是简称——DEM,可谓是故作高深,搞得不明真相的童鞋越发糊涂。其实很简单,高程这个概念不用说,“数字模型”就是把高程数字化后的模型嘛,数字化就是在电脑里面显示就是数字化咯,数字高程模型就是在电脑里面进行显示的高程数据嘛。
搞建筑规划的童鞋肯定都是知道现状地形图里面的高程点、等高线的,从某种程度上来说高程点、等高线就是高程数据的一种数字模型,DEM只是把高程进行连续面状数字化了,那对应的等高线就是线状高程数字化模型,高程点就是点状高程数字化,这样就明白为什么我们可以用等高线或者高程点生成DEM了吧。当然DEM不只是用高程点或者等高线生产的,还有卫星遥感数据等渠道可以获取,这个后面会讲到。
在ArcGIS里面用高程点或等高线生成DEM也不是直接生成,需要通过先生产TIN文件,在将TIN文件装换为DEM。
数字高程模型DEM灰度可视化示意图
有童鞋可能要问TIN文件又是什么鬼,为什么不能直接生成DEM,而非要先生成TIN文件呢?
TIN的全称是Triangulated Irregular Network,翻译过来就是它的中文名——“不规则三角网”,其实也是一种数字模型,当然,TIN文件可以是表达高程的数字模型,也可以是表达其他数据的数字模型本文里面的TIN文件当然就是表达高程的数字模型了,它和DEM的区别就在于他是矢量的,而DEM是栅格的!为什么会这样呢?
因为我们知道CAD里面的高程点和等高线都是矢量的,我们把它们导入ArcGIS后依然还是矢量的,我们要把矢量的点状高程点或者线状等高线高程数据转换为连续面状的DEM高程数据模型是需要插值的(什么是插值?简单理解就是在空白地方增加数据,点状、线状转换为面状就是要增加数据的过程!),而这个插值过程是不能直接将矢量插值为栅格数据的,所以要先插值生成TIN格式的矢量数字三个网络模型,再将其转化为栅格数字高程模型DEM文件。
不规则三角网模型TIN文件的可视化效果
好了,关于ArcGIS地形分析里面两个核心数据TIN和DEM已经说清楚了,再来看看我们熟悉又让人琢磨不透的CAD高程数据吧。
我们常见的CAD高程数据基本就两种:高程点、等高线。做ArcGIS地形分析对CAD高程点和等高线的要求都是必须带有高程属性!
高程点的高程属性比较简单,如图3.3所示,我们在CAD现状地形图中的大部分高程点都是如图这种高程点图块,图块带有一个高程点及其高程值标注,此时我们双击这个高程点标注文字是不能像传统文字那样编辑的。
ArcGIS也是不能直接识别CAD中这样的图块,所以我们不能直接将这一的高程点导入ArcGIS,需要我们进行处理——炸块!
但要注意的是此处的炸块不是用CAD中的X命令,如果用X命令对这个块进行炸裂,会发现整个图块就消失了,这是很多童鞋在处理这样高程点时会遇到的问题,正确的做法是通过湘源控规(请自行安装)里面的“工具-图块-炸属性块”对高程点图块进行炸块处理。输入这个工具命令后,会提示是删除还是保留隐藏属性,这里可以默认即可。
炸属性块后根据不同情况会有不同显示,但标注高程点的文字肯定是会单独显示的,我们通过双击也发现是可以像传统文字那样进行编辑修改了,这正是我们需要用到的数据。
常见CAD高程点示意图
湘源控规炸属性块工具示意图
等高线的属性数据相对复杂一些,图中绿色的线即是我们需要提取使用的地形等高线,其中比较细的等高线叫做基本等高线,也叫首曲线,是一副地形图中按基本等高距绘制出的等高线;其中比较粗的等高线叫加粗等高线,也叫计曲线,是按每隔五条基本等高线加粗的一条等高线,主要是为了便于判读和计算高程使用。
两根计曲线的高程差为2.5米,则基本等高线的高程差为0.5米。
CAD等高线数据示意图
各条等高线的高程查看方法为,在CAD中选中需要查看高程的等高线打开进入其属性查看Elevation属性值即为其高程。
CAD查看图元属性可以通过PR(Properties)命令、Ctrl+1打开属性对话框或者通过LI(List)命令调取图元属性列表等方法。
CAD等高线高程属性查询
如果我们所拿到的CAD现状图中具有以上特征的高程点或等高线,那我们基本就能导入ArcGIS正确生成地形分析所需要的TIN文件和DEM数字高程模型,并做进一步的相关分析了。
1.CAD高程点和等高线的添加
如前文所述,需要在CAD中准备好带高程的高程点和等高线数据。通常我们的现状地形图是多种图层叠加的,需要首先将这两个数据所在图层单独分离出来(也有是所有图元被合并在一个图层的,这种比较麻烦些,后面会讲怎么操作),图层的分离可以在CAD里面提前做好了再导入ArcGIS,也可以导入ArcGIS转换为Shapefile后再做分层处理,此处主要讲一下如何导入ArcGIS后做分层处理(CAD比较简单,直接单独按图层单独选择并写块即可,但是如果数据量大的话在CAD里面会比较卡顿,所以建议先导入ArcGIS再做分离)。
首先需要在CAD里面明确高程点和等高线各自的图层名,如本教程所采用数据对应的图层名分别为GCD、DGX。打开ArcGIS for Desktop中的其中一个软件ArcMAP,添加准备好的CAD现状地形图。
步骤为:
方法一:点击菜单栏“文件-添加数据-添加数据”;
ArcMAP中添加数据方法一
方法二:在菜单栏空白处点击鼠标右键确认标准工具栏打开,点击标准工具栏里的数据添加工具进行添加;
ArcMAP中添加数据方法二
方法三:打开ArcCatalog(在ArcGIS启动目录下)或者在标准工具栏里打开目录工具,找到对应的地形CAD文件直接拖入即可。
ArcMAP中添加数据方法三
2.CAD高程点和等高线的提取
CAD添加进ArcGIS后的数据列表,我们可以看到一个CAD文件添加进ArcGIS后会出现5个文件,分别是Annontation、Point、Polyline、Polygon、Multipatch(具体是什么后面会另外再讲,这里不展开,有兴趣的小伙伴也可以先自行百度),我们分别对高程点和等高线所属的Annotation和Polyline文件进行分离提取并生成对应的Shapefile文件。(注意:这一步是非常重要的一步)
CAD导入ArcGIS后的数据列表
操作步骤:
(1)通过标准工具栏打开工具箱(ArcToolbox)窗口,并依此展开至“分析工具-提取分析-筛选”。
CAD导入ArcGIS后的数据列表
(2)双击“筛选”工具,在弹出的筛选工具窗口的输入栏选中Annotation数据;输出要素类位置可默认,也可按自己需要输入对应文件位置;表达式一栏点击后侧SQL表达式图标弹出表达式输入框, 依此双击各参数并单击相应按钮,此处的表达式含义就是将CAD中图层(Layer)名为“GCD”的高程点数据提取出来。
(前面说了要明确高程点和等高线的图层名就是此处要用到!)完成后点击确定即可生成对应的高程点Shapefile文件Annotation_Select1。
通过属性筛选工具提取GCD图层数据示意
按相同的方法将等高线对应的Polyline文件进行提取分离,这次的表达式为"Layer"="DGX",对应生成等高线Shapefile文件Polyline_Select1。
3.高程数据的初步修正
通过上述1、2步的操作,我们已经成功的将CAD格式的高程点和等高线数据添加到ArcGIS中并将其各自转换为ArcGIS可编辑的Shapefile格式,对应的文件分别为Aannotation_Select1、Polyline_Select1。
按一般的正常步骤到了这一步即可开始生成TIN文件了,但是通常我们导入的数据或多或少都有些问题,尤其是高程值不正常会导致生成的TIN文件出现错误,下面就介绍一下如何检查高程值的错误以及修复方法。
(1)高程点错误数据的修正
在内容列表里右键点击Aannotation_Select1图层,选择“打开属性表”,在属性表里我们找到Elevation、Text两列属性值,其中Elevation即为高程点对应的高程值,Text即为CAD里面标注的高程值,此处的Elevation值是真实值我们可以用,但是如果此处的Elevation值为空或者全部为零我们就不能用这个值来计算TIN文件了,就需要用到Text来间接计算。
我们实验数据里的Elevation值非空,可以使用,但要检查数据的准确性,一是要检验其数据的正确性,检验方法是通过和Text数据列对比,如果基本是能一一对应的,那数据也就基本是准确的;另外是要检查是否有错误值,检查方法是在表头处双击(啥是表头?就是显示Elevation这一列表的名字的地方!),切换数据按高低顺序显示,我们可以看到当Elevation数据列按从高到低排列后,出现了部分数据出现严重偏差,和Text数据列所标注高程不对等的情况,数据值为0甚至是负数,这是明显的错误高程值,需要对错误的高程点筛选出来进行修正,方法如下:
检查高程点数据中Elevation和Text属性值的差异性
Elevation和Text属性值存在偏差的部分
①点击属性表-表选项-按属性选择,打开属性选择对话框。
在属性表中打开属性选择工具
4.生成TIN文件
5.生成DEM文件