41用已矢量化的等高线在arcgis里做坡度、坡向、起伏度的分析

在arcgis中中，进行如下操作：

1、创建TIN

打开3d analyst模块，利用creat /modify TIN---creat TIN from features命令（height source 选择高程字段），先将等高线转为TIN；

2、从TIN中创建栅格表面

打开3d analyst模块，利用convert---TIN to raster命令（attribute选择elevation，cell size自定义，若为大比例尺数据可以选择5或10，可以参考相关研究文献），生成栅格表面，即DEM；

（备注：矢量化的等高线必须比研究区的范围大些，创建TIN并生成Raster后，再用研究区边界来裁切，这样的DEM数据才能满足精度要求）

3、地形因子分析

打开3d analyst模块，利用surface analysis---slope命令，生成坡度数据；

打开3d analyst模块，利用surface analysis---aspect命令，生成坡向数据；

打spatial analyst模块，利用neighborhood tatistics命令进行邻域分析，先将statistic type设为最大值，输出栅格为A，再将statistic type设为最小值，输出栅格为B，利用raster calculator 生成地形起伏度数据，公式为[A]-[B]；

以上的地形数据，可以根据需要进行reclassfy重分类处理，分类标准参考相关文献，就可以获取所需的地形因子统计数据。

制图时，用view---layout view，添加比例尺、指北针、图例，就可以整饰出图

用arcgis绘制等高线

用arcgis绘制等高线 目录 1.获得高程数据 (2) 1.1空间地理数据云 (2) 1.2Google地球 (2) 2.导入高程数据 (3) 2.1空间地理数据云DEM数据的导入 (3) 2.2 Google地球取点数据的导入 (3) 3.插值 (3) 4.裁剪与相交 (4) 4.1裁剪 (4) 4.2相交 (4) 5.绘制等高线 (5) 6.处理等高线 (5) 6.1拓扑 (5) 6.2掩膜 (6) 6.3平滑处理 (6) 7.制图效果 (7) 8.制图心得 (7) 9.出图 (8)

关于高程数据的来源主要有两种，分别是空间地理数据云和Google地球。下面分别介绍获取方法。 1.1空间地理数据云 1)注册空间地理数据云账号并登录 2)在主页面点击“模型服务”，选择“DEM高程数据切割”并运 行 3)在新页面地图视口选择需要的区域，并单击下载，将其保存 于制图文件夹中 1.2Google地球 1)下载Google地球取点工具，打开谷歌地图 2)将目标区域选入Google地球的地图视口中，运行取点工具， 选择需要的取点精度，点击开始取点 3)取点结束后保存，并默认文件名为：采集文本1

基于不同高程数据来源，储存高程数据的格式不同，其也有不同的导入方式。下面介绍两种高程数据的导入。 2.1空间地理数据云DEM数据的导入 1)在ArcMap标准工具条中点击【导入数据】，在保存DEM数据 的文件夹中选中下载的DEM数据，点击确定 2.2 Google地球取点数据的导入 1)打开取点工具所保存的文本文档数据文件，在最前面加一行， 并输入““N”，“Y”，“X”，“Z””，其代表四列数据。N列表示 点号，Y列表示对应点的纬度，X列表示对应点的经度，Z列 表示对应点的高程。 2)在Arcmap主菜单中单击【文件】→【添加数据】→【添加 X,Y数据】，并在对话框中X字段、Y字段、Z字段中分别输 入X、Y、Z，点击确定 3.插值 在空间地理数据云上下载的高程数据可以忽略此步，无需插值。用Google地球取点下载的高程数据相当于点要素，需要插值转化为删格数据，然后才能绘制等高线。具体步骤如下： 1)在Arcmap标准工具条中单击【ArcToolbox窗口】

地形因子计算详解

第七章1、本章主题编号 2、本章内容概述 （1）概述 ●坡面因子的分类及提取方法 ●确定坡面因子提取的算法基础 ●提取坡面因子的常用分析窗口 （2）坡度、坡向 ●坡度的提取 ●坡向的提取

（3）坡形 ●宏观坡形因子 ●地面曲率因子 ●地面变率因子 （4）坡长 （5）坡位 （6）坡面复杂度因子 3、本章内容 3.1 概述 （1）坡面因子的分类及提取方法 ●坡面因子的分类 按照坡面因子所描述的空间区域范围，可以将坡面因子划分为微观坡面因子与宏观坡面因子两种基本类型。常用的微观坡面因子主要有：坡度、坡向、坡长、坡度变率、坡向变率、平面曲率、剖面曲率等。常用的宏观坡面因子主要有：地形粗糙度、地形起伏度、高程变异系数、地表切割深度，以及宏观坡形因子（直线形斜坡、凸形斜坡、凹形斜坡、台阶形斜坡）等。

按照提取坡面因子差分计算的阶数，可以将坡面因子分为一阶坡面因子、二阶坡面因子和高阶坡面因子。一阶坡面地形因子主要有坡度和坡向因子。二阶坡面因子主要有坡度变率、坡向变率、平面曲率、剖面曲率等因子。复合坡面因子有坡长、坡形因子、地形粗糙度、地形起伏度、高程变异系数和地表切割深度等。 按照坡面的形态特征，可将坡面因子进一步划分为：坡面姿态因子，坡形因子，坡位因子，坡长因子以及坡面复杂度因子五大类。 ●提取坡面因子的基本方法 首先将坡面的形态特征或各个坡面因子进行定量化描述，完成求导的数学模型，在此基础上，建立其以DEM为基本信息源进行提取的技术路线，并通过软件实现形成一套易于计算机操作的方法。 （2）确定坡面因子提取的算法基础 ●DEM格网数据的空间矢量表达（如图7.1） 图7.1 DEM格网数据的空间矢量模型

mapgis绘制等高线方法

用Mapgis程序绘制物化探等值线图及点位数据图方法说明 赵明亮（李金忠整理） 一、建立数据表 1、用电子表格建立成图元素表（不要表头，建议将平面坐标按Y、X、数据…列次序移到表前1、2…列，Y、X坐标分别为7位、6位整数）； 2、另存为文本文件（制表符分隔）类型的新文件并取名(如文件名为茨竹箐原始数据表)。 二、Mapgiss投影转换 1、打开Mapgiss程序主菜单→点击″实用服务″→点击″投影变换″→点击″P投影转换″→点击″U 用户文件投影转换″→打开文本文件(茨竹箐原始数据表)； 2、选″指定分隔符″→指定坐标X、Y列（注意X、Y列可能倒置）； 3、点击″设置分隔符″→选″Tab″及″空格″→点击″数据类型″→将所有的字符串点为″双精度″→在″属性名称所在行″选框中选择数据（注意数据不能选″无″）→在″加入″列选定所需的坐标及元素值→确定； 4、点击″点图元参数″→选择子图参数（建议选编号35的圆点的子图号，高度宽度一般均为1，颜色为1）→确定； 5、选″不需要投影″→点击″数据生成″→确定→关闭； 6、另取点文件名（如茨竹箐数据投影)→保存→退出。 三、等值线图的Mapgiss的DTM分析 1、打开Mapgiss程序主菜单→点击″空间分析″→点击″DTM分析″→点击″文件″→点击″打开数据文件″→点击″点数据文件″→打开点文件（如茨竹箐数据投影）； 2、点击″处理点线（P）″→点击″点数据高程提取″→在″高程属性项″中点击要成图的元素的值→确定； 3、点击″Tin模型″→点击″快速生成三角剖分网″→在是否保存选项中选择″否″； 4、点击″Tin模型″→点击″删除三角剖分网边″→手工删除不合理的剖分网边； 5、点击″Tin模型″→点击″追踪剖分等值线″； ①选″等值线套区″→选″绘制色阶″→可选或不选″保留边界″→选″等值线光滑处理″→在″光滑度″选项中选择″高程度″； ②点击″等值层值″→对等值线层分段参数中的″起始Z″、″终止Z″、″步长增″及″起始色″、″终止色″进行设置→点击″更新当前分段″→确定；

ArcGIS支持的栅格数据集文件格式

ArcGIS支持的栅格数据集文件格式 在ArcGIS 中，有两种使用栅格数据的方法：作为栅格数据集和作为栅格类型。栅格数据集用于定义像素的存储方式，例如，行数和列数、波段数、实际像素值，以及其他栅格 格式特定参数。栅格类型特定于某些产品，例如卫星传感器或文件格式范围之外的商业数据格式。栅格类型用于与栅格格式一起识别元数据，例如地理配准、采集日 期和传感器类型。 大多数使用栅格数据集的情况是将其添加到显示和执行地理处理等操作。但在向镶嵌数据集添加数据时，需要识别是将其作为栅格数据集添加，还是使用特定栅格类型进行添加。要确定对镶嵌数据集使用添加栅格至镶嵌数据集工具时支持哪种数据形式（栅格类型或栅格数据集），请参阅栅格数据和栅格类型支持列表。了解有关栅格类型的信息地 理数据库是ArcGIS 中的原生数据模型，可用于储存包括栅格数据集、镶嵌数据集和栅格目录在内的地理信息；然而，还有很多可能用到的文件格式保留在地理数据库之外。下表描述了 受支持的栅格格式（栅格数据集）及其扩展名，并识别出了它们是否为只读格式或是否还可以由ArcGIS 写入。

注：可以指定需要ArcGIS 识别哪些产品（自定义> ArcMap 选项> 栅格> 文件格式）；关闭不使用的产品可以提高性能。您也可在数据加载过程中过滤数据。有关详细信息，请参阅显示特定栅格格式。 格式描述扩展模块读/写 ARC 数字化栅格图形(ADRG) 由美国国家地理空间情报局(NGA) 以CD-ROM 的形式发布。使用等弧秒栅格图/地图(ARC) 系统（将地球划分为18 个纬度带或区域）对ADRG 进行地理配准。该数据由通过扫描源文档所生成的栅格影像及其他图形组成。 多个文件数据文件- 扩展名*.img 或*.ovr图例文件- 扩展名*.lgg 只读 ArcSDE 栅格

ArcGIS ArcMap中修改栅格图的像元值

在处理图象数据时，我们经常会碰到要求修改栅格图象象元值的问题，比如说DEM图的部分数据错误，我们要进行修改；再比如说栅格图象中有些与周围均匀色彩不一致的错误斑块要更正等等，那我们如何来处理这一类问题呢？ 现我以一DEM栅格图（名字为eldodem）为例，现在我要修改它的部分象元值，总结出以下三种方法，大家可以参考一下。第一种方法不大实用，但可借鉴，第二三种方法针对的条件不一致，大家可以在具体情况下进行选择。 一、直接运用转换，思路简单，易操作，但实用性 1、栅格——ascii文件——栅格这种方法是先将栅格图用工具 直接转成ascii文件，然后在ascii文本文件中直接修改需要修改的象元的值，修改好后又用工具 转换成栅格图。这种方法可行，但是不实用，因为我们要搜索到指定的象元好像不是那么简单。那么有没有别的比较好的方法呢？ 二、在栅格计算器中操作，方法灵活，可操作性强，实用性强

准备工作 先要在option中设置保留的栅格范围，通常情况下默认的为相交后的部分，这里我们要保留整个DEM，所以要改为以下设置： 2、通过点的位置修改点象元值 2．1* 问题一、要是我想修改图象中的指定的行列的点的象元值（比如说把第100行，200列的点的值修改为0），那该如何操作呢？ 这时候我们可以在栅格计算器中输入以下公式： 1.con(($$rowmap == 99 & $$colmap == 199),0,[eldodem.img]) (栅格图的编号是从0行0列开始的)

2．2* 问题二、那若是要修改指定行列范围内的栅格的象元值，比如说把第101行，251列到第401行，301列的部分的象元值改为0，又改如何进行呢？ 1.con(($$rowmap < 400 & $$rowmap > 100 & $$colmap < 300 & $$colmap > 250),0,[eldodem.img]) 结果如下

地形因子

第七章 1、本章主题编号 2、本章内容概述 （1）概述 ● 坡面因子的分类及提取方法 ● 确定坡面因子提取的算法基础 ● 提取坡面因子的常用分析窗口 （2）坡度、坡向 ● 坡度的提取 ● 坡向的提取 （3）坡形 ● 宏观坡形因子 ● 地面曲率因子 ● 地面变率因子 （4）坡长 （5）坡位 （6）坡面复杂度因子 3、本章内容 3.1 概述 （1）坡面因子的分类及提取方法 ● 坡面因子的分类 按照坡面因子所描述的空间区域范围，可以将坡面因子划分为微观坡面因子与宏观坡面因子两种基本类型。常用的微观坡面因子主要有：坡度、坡向、坡长、坡度变率、坡向变率、平面曲率、剖面曲率等。常用的宏观坡面因子主要有：地形粗糙度、地形起伏度、高程变异系数、地表切割深度，以及宏观坡形因子（直线形斜坡、凸形斜坡、凹形斜坡、台阶形斜坡）等。 按照提取坡面因子差分计算的阶数，可以将坡面因子分为一阶坡面因子、二

阶坡面因子和高阶坡面因子。一阶坡面地形因子主要有坡度和坡向因子。二阶坡面因子主要有坡度变率、坡向变率、平面曲率、剖面曲率等因子。复合坡面因子有坡长、坡形因子、地形粗糙度、地形起伏度、高程变异系数和地表切割深度等。 按照坡面的形态特征，可将坡面因子进一步划分为：坡面姿态因子，坡形因子，坡位因子，坡长因子以及坡面复杂度因子五大类。 ● 提取坡面因子的基本方法 首先将坡面的形态特征或各个坡面因子进行定量化描述，完成求导的数学模型，在此基础上，建立其以DEM为基本信息源进行提取的技术路线，并通过软件实现形成一套易于计算机操作的方法。 （2）确定坡面因子提取的算法基础 ● DEM格网数据的空间矢量表达（如图7.1） 图7.1 DEM格网数据的空间矢量模型 ● 基于空间矢量模型的差分计算 算法主要有数值分析方法、局部曲面拟合算法、空间矢量法、快速傅立叶变换等。其中数值分析方法包含有简单差分算法、二阶差分、三阶差分（带权或不带权）和Frame差分；局部曲面拟合又有线性回归平面、二次曲面和不完全四次曲面（据刘学军，2002）。 （3）提取坡面因子的常用分析窗口 ● 窗口分析（领域分析）的基本原理是：对栅格数据系统中的一个、多个栅格点或全部数据，开辟一个有固定分析半径的分析窗口，并在该窗口内进行诸如极值、均值、标准差等一系列统计计算，或进行差分及与其它层面信息的复合分析等，实现栅格数据有效的水平方向扩展分析。 ● 在坡面信息提取中，按照分析窗口的形状，可以将分析窗口划分为以下几类： 矩形窗口：以目标栅格为中心，分别向周围八个方向扩展一层或多层栅格。 圆形窗口：以目标栅格为中心，向周围作一等距离搜索区，构成一圆形分析窗口。

ArcGIS教程：栅格计算器的工作原理

栅格计算器用于在工具中创建和执行“地图代数”表达式。与其他地理处理工具一样，栅格计算器可以在“模型构建器”中使用，从而可以更轻松地在工作中发挥“地图代数”的功用。 栅格计算器具有如下的独有优势： 1、执行单行代数表达式。 2、使用“模型构建器”时，支持在“地图代数”中使用变量。 3、为一个表达式的三个或更多的输入应用 Spatial Analyst 运算符。 4、在一个表达式中使用多个 Spatial Analyst 工具。 栅格计算器用于通过简单的、类似计算器的工具界面，使用多种工具和运算符来执行单行代数表达式。当在一个表达式中使用多种工具或运算符时，该公式的执行速度一般会比分别执行每个运算符或工具的速度要快。 使用栅格计算器工具 工具对话框中有四个主要区域用于创建“地图代数”表达式： 图层和变量 输入图层和变量列表识别了可以在“地图代数”表达式中使用的输入。单独使用此工具时，该列表将包含内容表中的图层。在“模型构建器”中使用此工具时，该列表将包含内容列表中的图层、模型中创建的输出以及其他模型变量。 计算器和运算符按钮

使用计算器按钮可以在表达式中输入数值。使用运算符按钮可以将数学运算符(加、除等)和逻辑运算符(大于、等于等)输入表达式。通过单击这些按钮，可以将数字或运算符输入到表达式中光标当前所在的位置。注意确保表达式遵循正确的规则地图代数语法规则。 工具 在该工具列表中，可以方便地选择能在“地图代数”表达式中使用的工具。通过单击列表中的工具，该工具的名称以及括号的左右半边 [()] 将被放置在表达式中光标当前所在的位置。然后必须输入该工具所需的其他输入内容。虽然工具可以放置在表达式的任何位置，但是它应放置在可以生成有效的“地图代数”语法的位置。 表达式 “地图代数”将要执行的表达式。此表达式必须使用有效的语法输入。请参阅以下部分，了解有关“地图代数”语法规则的详细信息。 地图代数语言 “地图代数”是一种简单而强大的代数语言，利用它您可以操作所有 Spatial Analyst 工具、运算符及函数以执行各种地理分析。 在栅格计算器中使用的“地图代数”含有语法或规则集，“地图代数”必须遵循这些语法或规则集才能创建有效的表达式。如果不遵循这些规则，创建的表达式可能会无效并且无法执行，或者得不到预期的结果。 “地图代数”不仅可以通过栅格计算器工具来访问，还可以通过使用 Spatial Analyst ArcPy 模块的Python 脚本来访问。在此工具中使用相同的“地图代数”语法，但以下情况除外： · 不必将输出栅格名称或等号 (=) 填入表达式，因为该输出名称已在输出栅格参数中指定。 · 在使用运算符时，不必将输入栅格转换为栅格对象。

arcgis栅格计算器

ARCGIS中的栅格计算 栅格计算是栅格数数据空间分析中数据处理和分析中最为常用的方法，应用非常广泛，能够解决各种类型的问题，尤其重要的是，它是建立复杂的应用数学模型的基本模块。ArcGIS 9 提供了非常友好的图形化栅格计算器，利用栅格计算器，不仅可以方便的完成基于数学运算符的栅格运算，以及基于数学函数的栅格运算，而且它还支持直接调用ArcGIS 自带的栅格数据空间分析函数，并且可以方便的实现多条语句的同时输入和运行。 一数学运算 数学运算主要是针对具有相同输入单元的两个或多个栅格数据逐网格进行计算的。主要包括三组数学运算符：算术运算符，布尔运算符和关系运算符。 1. 算术运算 算术运算主要包括加、减、乘、除四种。可以完成两个或多个栅格数据相对应单元之间直接的加、减、乘、除运算。 例如，以今年与去年的降水量数据为基础，用公式（今年降水量-去年降水量）/去年降水量，可以计算出去年降水量的变化程度，如图8.65。（单位：毫米）

图1 算术运算示意图 2. 布尔运算 布尔运算主要包括：和（And）、或（Or）、异或（Xor）、非（Not）。它是基于布尔运算来对栅格数据进行判断的。经判断后，如果为“真”，则输出结果为1，如果为“假”，则输出结果为0。 （1）和（&）：比较两个或两个以上栅格数据层，如果对应的栅格值均为非0 值，则输出结果为真（赋值为1），否则输出结果为假（赋值为0）。 （2）或（|）：比较两个或两个以上栅格数据层，对应的栅格值中只要有一个或一个以上为非0 值，则输出结果为真（赋值为1），否则输出结果为假（赋值为0）。 （3）异或(!)：比较两个或两个以上栅格数据层，如果对应的栅格值在逻辑真假互不相同（一个为0，一个必为非0 值），则输出结果为真（赋值为1），否则输出结果为假（赋值为0）。 （4）非(^)：对一个栅格数据层进行逻辑“非”运算。如果栅格值为0 ，则输出结果为1；如果栅格值非0，则输出结果为0。

ArcGIS实验-Ex22-坡向变率(SOA)

第十二章空间分析建模 练习2：坡向变率（SOA） 一、背景 平面曲率即地面坡向变率，是指在地表的坡向提取基础之上，进行对坡向变化率值的二次提取，亦即坡向之坡度（Slope of Aspect，SOA）。地面坡向变率是一个反映等高线弯曲程度的指标，可以反映出地表所有的山脊线、山谷线。 值得注意的是：SOA在提取过程中在不同的坡面上将会有误差的产生，即在坡面的南北两侧，北面坡由于在坡向算法将会有误差产生，所以要对北坡的SOA结果进行纠正，因为从理论上讲SOA在地表北坡上将产生误差，北坡上坡向值范围为0-90°和270°-360°，在正北方向附近，15°和345°之间坡向差值只是30°，而在计算中却是差了330°，所以要利用反地形将北坡地区的坡向变率误差进行纠正。 二、目的 通过纠正平面曲率的例子，使读者了解如何在模型中调用用户自定义的模型，明晰模型嵌套的过程。 三、要求 为了得到纠正后的SOA，必须得到正地形的SOA1和反地形的SOA2，修正的公式为((SOA1+ SOA2) - (SOA1-SOA2).abs)/2。其数据流如图7所示。 四、数据 DEM栅格数据。 五、操作步骤 1、建立生成SOA的模型 (1) 打开ArcMap，启动ArcToolbox。 (2) 右键ArcToolbox，选择new toolbox，生成toolbox2，如图1。 (3) 右键toolbox2，在new中选择model，，则生成model1，如图2。 图1 选择new toolbox 图2 选择model (4) 右键图标，选择rename命令，输入generate SOA。 (5) 在模型窗口右键，选择create variable命令，在数据类型选择框中选中raster dataset，如图3所示。

用ARCGIS 10 CAD等高线转为DEM的教程

1.启动Spatial ETL Tool工具： 在“Untitled - ArcMap - ArcInfo”中左键单击 打开catalog 在“Toolbox (框线项目)”(位于“Untitled - ArcMap - ArcInfo”中)上用户左键单击选中toolbox 在“Toolbox (框线项目)”(位于“Untitled - ArcMap - ArcInfo”中)上用户右键单击展开菜单

在“Spatial ETL Tool (菜单项目)”上左键单击 在弹出的菜单中选中Spatial ETL Tool 2.设置Spatial ETL Tool工具 在“Create Translation Workspace Wizard (应用程序)”(位于“Create Translation Workspace Wizard”中)上左键单击打开输入文件的类型窗口

在“FME Reader Gallery (应用程序)”(位于“FME Reader Gallery”中)上左键单击 选中AutoCAD DWG/DXF，以打开dxf和dwg文件 在“Create Translation Workspace Wizard (应用程序)”(位于“Create Translation Workspace Wizard”中)上用户左键单击进入下一步，选择导入的CAD文件

在“Create Translation Workspace Wizard (应用程序)”(位于“Create Translation Workspace Wizard”中)上左键单击打开CAD文件所在位置 在“处理后地形图.dwg (列表项目)”(位于“Select File”中)上左键单击 选中CAD文件

arcgis实之栅格运算

栅格分析 从DEM上提取山顶点，凹陷点 数据 DEM数据，地理信息：KRASOVSKY-1940-Transverse-Mercator Cellsize 5×5 单位：meter 学校选址 数据 DEM LANDUSE 栅格数据无坐标信息 School rec_sites 矢量数据点类型无坐标信息 设置栅格分析环境 Cellsize设置，输出单元格大小设置 Mask设置，分析范围设置 从Dem上提取山顶点，凹陷点 原理 DEM （数字高程模型）属于场模型（连续的栅格），表达和存储时，使用单元格中间值代替整个栅格的值，可视化时，采用插值（如线性）生成连续的栅格 采用领域的方法，可以用领域中的最高值或最低值代替该单元格的值，当该单元格是山顶点时（相对领域而言），领域计算后的值（取最高值）与原值相同，相反，凹陷点也是如

此 操作 1.加载DEM add date 或直接从catalog 中拖入 使用IDENTIFY工具查看每个单元格时，每个单元格只有一个值，对单元格而言是离散的，但对DEM而言，是连续的 2.环境设置 在10.1中，设置在ARCTOOLS的环境参数设置中完成 常规设置 输出工作空间设置 栅格分析设置 栅格存储设置 关于构建金字塔等，此处默认 领域运算 在10.1中，这个工具是focal Ststistics 在这个工具中，由两部分 NEIGHBERHOOD部分，设置领域类型，大小。类型有：矩形，园，环，楔形，和用户自定义的，大小用户自己定义 Stastistics type 就是领域的功能选择，此处选择max，既将领域最大值赋予单元格 结果预计：计算后的高值和DEM高值相等，低值应该大于等于DEM的低值

DEM内插算法对坡度坡向的影响

第15卷第6期2008年12月 水土保持研究 Research of Soil and Water Conservation Vol.15,No.6 Dec.,2008 　 D EM内插算法对坡度坡向的影响3 陈吉龙1,武伟2,刘洪斌1 (1.西南大学资源环境学院,重庆400716;2.西南大学计算机与信息科学学院,重庆400716) 摘　要:虽然DEM的应用越来越广泛,但是在实际研究运用中获取DEM的方式主要是利用各种矢量高程信息通过插值而来,不同的插值算法必然会影响DEM分析的结果;以南方丘陵地区为研究区域,利用样区的1∶10000的矢量地形图为数据源,通过对不同算法生成不同分辨率的DEM坡度、坡向的对比研究发现:在坡度小于30°时,不同的算法对坡度具有显著的影响,而大于30°时分辨率和算法均不会对坡度产生显著影响;研究表明:在具有大量平坦区域的地区,尽管ANUDEM算法生成的DEM总体上的精度最高,但是不会产生平坦区域,所以并不适合于这类地区的地形分析。 关键词:DEM;坡度;坡向;地形分析 中图分类号:S157;TP79 文献标识码:A 文章编号:100523409(2008)0620014204 E ffect of DEM Interpolation on the Slope and Aspect C H EN Ji2long1,WU Wei2,L IU Hong2bin1 (1.College of Resource and Envi ronment,S outhwest Universit y,Chongqing400716,Chi na;2.College of Com puter and I nf ormation Science,S outhw est Universit y,Chongqing400716,China) Abstract:In many cases of the research and application,digital elevation data may be the only source for DEM generation with algorithms,which are available or affordable.It is known that the analytic result f rom the DEM can vary in quality depending on their algorithms.Taking hilly region in southern China as research region,using the digital contour with the scale of1∶10000as the source,this paper compared the slope,aspect f rom the DEM generated with five algorithms,the results showed that there were significant effect on the slope among the algorithms when the slope<30°,but there were no significant effect on the slope among the algorithms and resolution when the slope>30°.Furthermore,the result also showed that DEM interpola2 ted with ANUDEM is the most precise,however,there was no flat area in the result f rom the terrain analysis in the region with extensive flat area,so it can’t be suitable for performing terrain in such region. K ey w ords:DEM;slope degree;slope direction;terrain analysis 1　引言 自从DEM理论形成以来,不同的学者就对DEM做过大量的研究,研究领域侧重于利用DEM来进行地形分析、水文分析、DEM精度评价、地形描述误差、数学模拟方法以及探索新的DEM算法[1];其中地形分析是DEM研究的重要内容,是地质地貌研究的重要手段;虽然地形分析中所采用的DEM建立方法多种多样,但是当前乃至今后一段时期内DEM的生产依然是利用各种矢量高程信息通过构建TIN建立DEM;然而该方法的缺点是仅利用等高线生成DEM不能很好地反映河流和流域边界等重要地貌特征[2]; Hutchinson于20世纪80年代提出的ANUDEM插值算法[3],大量的研究也表明ANUDEM算法生成的DEM能更加精确地描述地貌特征[4],其质量较现有基于TIN方法建立的DEM有明显提高[5],因而在全球范围内得到广泛应用[627],然而国内对于ANUDEMS算法的应用很少见;在地形分析中,虽然众多学者对相同数据源生成的DEM的不同分辨率对地形分析结果的影响进行了研究[729],然而其中很多的分辨率太低,忽略了地形描述精度,使分析结果受到严重的影响;本文对不同算法生成的不同分辨率的DEM提取的坡度、坡向进行对比分析,以期为地形分析中DEM的建立和分辨率的选择及相关的研究提供参考。 2　材料与方法 2.1　研究区概况 研究区地处长江上游,重庆市西南边缘,位于重庆江津市,介于东经106°10′45″-106°23′42″和北纬29°22′50″-29°31′17″,面积79.36km2;地势中高西低,中部低山呈东北走向,东西以丘陵为住,最高海拔765.58m,最低海拔255.34 m,相对高差达510.24m。 3收稿日期:2008204218 　基金项目:重庆市科委重点攻关项目(2006AB1015);重庆市教委科技资助项目(0182) 　作者简介:陈吉龙(1983-),男,重庆巫山人,硕士研究生,从事3S应用研究。E2mial:cjl47168@https://www.wendangku.net/doc/c52328058.html, 通信作者:刘洪斌(1966-),男,重庆梁平人,研究员,从事3S应用研究。E2mial:lqb2000@https://www.wendangku.net/doc/c52328058.html,

利用ArcGIS制作真实地形图的方法

多媒体教学 2015 年 第 14 期 在中学地理教材上有许多地理信息均以地图的形式表达出来，且部分地理信息只能通过地图的形式加以呈现。地图对于地理的学习非常重要，地图有地理学的“第二语言”之称。近年来，对于等高线地形图的判读与应用归纳和总结的文章较多，然而很少有地理教师关注如何从出题者的角度科学地制作出一幅地形真实且精度较高的地形图，进而使原创试题设计更加合理。 一、地形图的分类 地形图描述了地表形态高低起伏的变化特征，一般包括等高线地形图、分层设色地形图和地形剖面图。其中，把陆地表面海拔相等的点连接成的线叫等高线，把海洋或者湖泊中深度相同的各点连接成的线叫等深线。在绘有等高线或等深线的地形图上，按照不同的高度和深度，着上深浅不同的各种颜色，使得地表形态和海底起伏状况清晰可见，这种地形图叫分层设色地形图，一般而言，对于黑白灰着色没有明确的要求，彩色中，往往以绿色、黄色、棕色等表示平原、高原和山地，以深浅不同的蓝色表示海洋的深度。地形剖面图则是根据等高线地形图沿着某一方向的直线所作的垂直断面图，这种地形图可以直观地表示地面上沿着某一方向地形的起伏状况。下面将从这三个层次依次介绍地形图的制作过程。 二、数据来源 要想制作一幅真实且具有一定精度的地形图，前提是获取相应精度的地形高程数据。一般而言，高程数 彭云龙 （江苏省南师附中新城初中黄山路分校, 江苏 南京 210019） 利用制作真实地形图的方法 摘要：新课程改革使地理信息技术进入中学地理教材，这要求中学地理教师需加深对地理信息技术的了解和应用。GIS技术在中学地理教学中应用广泛，本文以地形图的制作为例，探讨一线地理教师地形图制作过程中的不足，并分析地形图的种类，最后依次介绍了GIS技术在不同种类地形图制作中的具体过程。 关键词：ArcGIS；地形图；等高线 据的来源有两个途径：一是携带专业设备进行野外测量，然后在室内利用一定的插值方法进行空间插值得到高程数据；二是可以从权威机构共享的数据库中下载得到。第一种方法获取的数据精度较高，但难度较大，数据范围较小。对于中学地理教师而言，第二种数据源是比较切合实际的，并且数据精度一般可以达到出题者的要求。常见数据类型有SRTM-DEM数据和ASTER-GDEM数据[1]。SRTM-DEM数据有两种精度类型，即30m的SRTM-1和90m的SRTM-3高程数据。90m 的高程数据全球共享，可以从NASA官网获取（数据：https://www.wendangku.net/doc/c52328058.html,/srtm/），但是30m的高程数据仅覆盖美国地区。如果想获取30m的DEM数据，可以从中国科学院网络信息中心国际科学数据镜像网站（数据：https://www.wendangku.net/doc/c52328058.html,）免费下载。本文以30m的庐山地区的ASTER-GDEM数据为例。 三、地形图的制作过程1.等高线地形图的制作 (1)加载高(1)加载高程数据程数据 打开ArcMap9.3，点击“工具栏”上的“Arc Catalog”图标，找到高程数据（已经裁剪好的DEM）直接拖入“ArcMap9.3”主窗口中即可加载数据。如图1，黑色表示的海拔较低，最低为-53米，因为该区域属于长江水系，河流侵蚀基准面高于平均海平面（0米），故最低海拔是数据误差造成的，白色表示的海拔较高[2]，该地区最高海拔为1446米。

arcgis中对矢量和栅格数据进行裁剪切割的方法

arcgis中对矢量和栅格数据进行裁剪切割的方法 [引用 2011-11-07 13:21:46] ?? 字号：大中小? ArcGIS9:?? ArcToolbox-->Spatail Analyst Tools-->Extraction : 在ArcMap --> Spatial Analyst扩展模块工具条--〉Raster? Calculator 或者直接使用Mask ? 1、对矢量数据进行裁减：Arctoolbox中，spatialtool>extract>clip 在InputFeatures中选择被裁剪的图层，在ClipFeatures中选择裁剪形状 2、利用矩形对栅格数据裁剪:Arctoolbox中，datamanagement>raster>clip InputRaster中选择被裁剪的栅格数据，设定好矩形四个顶点即可 3、在SpatialAnalystTools中提供了多种对栅格数据的提取方法，Arctoolbox 中，SpatialAnalyst Tools>extraction包括提取值到点，根据属性提取，用圆提取，用多边形提取，用掩模提取，用点提取，用矩形提取等，其中的用掩模提取功能可以让我们通过不规则边界来获取需要的栅格数据。 ? 可以先建一个polygon的图层，再用这个图层clip那个raster ? 加载空间分析模块， option选项中设置mask图层，范围，cell尺寸，

用其栅格计算器计算： ? 图层名.shape*1, 点击evaluation， ok ? 如果不要求一定要多变形裁剪的话，建议使用clip data management toolbox---clip 这个工具可以保持原有像元信息，但是这里arcmap只提供了矩形切割。 这个约束应该不会影响一般的应用，否则的话建议使用图像处理软件中的mask 裁剪。 另外上面大家说得mask什么的方法，只能出来灰度图，可能原先是设计成处理grid的数据吧 ? ? ? ? ? (2) 用任意多边形剪切栅格数据(矢量数据转换为栅格数据) 在ArcCatlog下新建一个要素类（要素类型为:多边形），命名为：

坡度,坡向,提取等高线

?工具/原料 DEM的应用DEM的应用包括：坡度：Slope、坡向：Aspect、提取等高线、算地形坡度：Slope、坡向：Aspect、提取等高线、算地形表面的阴影图、可视性分析、地形剖面、水文分析表面的阴影图、可视性分析、地形剖面、水文分析等，其中涉及的知识点有： a)对TIN建立过程的原理、方法的认识； b)掌握ArcGIS中建立DEM、TIN的技术方法。 （对于这两步的教程本人之前有做过，下面教程不会再重复） c)掌握根据DEM  计算坡度、坡向的方法。 d)理解基于DEM数据进行水文分析的基本原理。 e)利用ArcGIS的提供的水文分析工具进行水文分析的基本方法和步骤。下面开始教程： 软件准备：ArcGIS   Desktop  10.0---ArcMap(3D  Analyst模块和spatial  a nalyst模块) 数据：DEM和TIN（使用由本人前面的教程【ArcGIS地形分析--TIN及DE M的生成，TIN的显示】得到的结果数据。 原始数据下载：https://www.wendangku.net/doc/c52328058.html,/s/1GGzT2

1 1

执行后，得到坡度栅格Slope_tingri1：坡度栅格中，栅格单元的值在[ 2 0  -82]  度间变化

3

【下面计算剖面曲率】 4 [3D  Analyst工具]——[栅格表面]在【ArcToolbox】中，执行命令[3D  Analyst工具]——[栅格表面] 到剖面曲率栅格：[————[坡度]  。按如下所示，指定各参数。得到剖面曲率栅格：[————[坡度] Slope  _Slope] 如图

ARCGIS中矢量裁剪栅格图像

ARCGIS中矢量裁剪栅格图像 (1) 是否需要裁剪栅格图象区域通过一个面状的shapefile表达出来？ 如果可以，那么就很简单了。 在ArcMap中，调用空间分析扩展模块，将你感兴趣区的shapefile多边形图层设置为掩膜，然后在栅格计算器中重新计算一下你的图象，它就会沿掩膜裁出。 设置掩膜：空间分析工具条的下拉菜单>option里面设置 (2) 用任意多边形剪切栅格数据(矢量数据转换为栅格数据) 2.1在ArcCatlog下新建一个要素类（要素类型为:多边形），命名为：ClipPoly.shp 2.2在ArcMap中，加载栅格数据：例如kunming.img、和ClipPoly.shp 2.3打开编辑器工具栏，开始编辑ClipPoly ，根据要剪切的区域，绘制一个任意形状的多边形。打开属性表，修改多边形的字段“ID”的值为1，保存修改，停止编辑。 2.4打开空间分析工具栏 执行命令：<空间分析>－<转换>--<要素到栅格> 指定栅格大小：查询要剪切的栅格图层kunming的栅格大小，这里假设指定为1 指定输出栅格的名称为路径 2.5执行命令: <空间分析>-<栅格计算器> 2.6构造表达式：[kunming]*[polyClip4-polyclip4] ，执行栅格图层：kunming和用以剪切的栅格polyClip4 之间的相乘运算 (3) 1、对矢量数据进行裁减：Arctoolbox中，spatial tool > extract>clip 在InputFeatures中选择被裁剪的图层，在ClipFeatures中选择裁剪形状 2、利用矩形对栅格数据裁剪:Arctoolbox中，data management>raster>clip InputRaster中选择被裁剪的栅格数据，设定好矩形四个顶点即可 3、在Spatial Analyst Tools中提供了多种对栅格数据的提取方法，Arctoolbox中，Spatial Analyst Tools>extraction包括提取值到点，根据属性提取，用圆提取，用多边形提取，用掩模提取，用点提取，用矩形提取等，其中的用掩模提取功能可以让我们通过不规则边界来获取需要的栅格数据。

第七章 坡面地形因子提取

1、本章主题编号 2、本章内容概述 （1）概述 ● 坡面因子的分类及提取方法 ● 确定坡面因子提取的算法基础 ● 提取坡面因子的常用分析窗口 （2）坡度、坡向 ● 坡度的提取 ● 坡向的提取 （3）坡形 ● 宏观坡形因子 ● 地面曲率因子 ● 地面变率因子 （4）坡长 （5）坡位 （6）坡面复杂度因子 3、本章内容 3.1 概述 （1）坡面因子的分类及提取方法 ● 坡面因子的分类 按照坡面因子所描述的空间区域范围，可以将坡面因子划分为微观坡面因子与宏观坡面因子两种基本类型。常用的微观坡面因子主要有：坡度、坡向、坡长、坡度变率、坡向变率、平面曲率、剖面曲率等。常用的宏观坡面因子主要有：地形粗糙度、地形起伏度、高程变异系数、地表切割深度，以及宏观坡形因子（直线形斜坡、凸形斜坡、凹形斜坡、台阶形斜坡）等。 按照提取坡面因子差分计算的阶数，可以将坡面因子分为一阶坡面因子、二阶坡面因子和高阶坡面因子。一阶坡面地形因子主要有坡度和坡向因子。二阶

坡面因子主要有坡度变率、坡向变率、平面曲率、剖面曲率等因子。复合坡面因子有坡长、坡形因子、地形粗糙度、地形起伏度、高程变异系数和地表切割深度等。 按照坡面的形态特征，可将坡面因子进一步划分为：坡面姿态因子，坡形因子，坡位因子，坡长因子以及坡面复杂度因子五大类。 ● 提取坡面因子的基本方法 首先将坡面的形态特征或各个坡面因子进行定量化描述，完成求导的数学模型，在此基础上，建立其以DEM为基本信息源进行提取的技术路线，并通过软件实现形成一套易于计算机操作的方法。 （2）确定坡面因子提取的算法基础 ● DEM格网数据的空间矢量表达（如图7.1） 图7.1 DEM格网数据的空间矢量模型 ● 基于空间矢量模型的差分计算 算法主要有数值分析方法、局部曲面拟合算法、空间矢量法、快速傅立叶变换等。其中数值分析方法包含有简单差分算法、二阶差分、三阶差分（带权或不带权）和Frame差分；局部曲面拟合又有线性回归平面、二次曲面和不完全四次曲面（据刘学军，2002）。 （3）提取坡面因子的常用分析窗口 ● 窗口分析（领域分析）的基本原理是：对栅格数据系统中的一个、多个栅格点或全部数据，开辟一个有固定分析半径的分析窗口，并在该窗口内进行诸如极值、均值、标准差等一系列统计计算，或进行差分及与其它层面信息的复合分析等，实现栅格数据有效的水平方向扩展分析。 ● 在坡面信息提取中，按照分析窗口的形状，可以将分析窗口划分为以下几类： 矩形窗口：以目标栅格为中心，分别向周围八个方向扩展一层或多层栅格。 圆形窗口：以目标栅格为中心，向周围作一等距离搜索区，构成一圆形分析窗口。

在ArcGIS中由等高线生成dem的步骤

在ArcGIS中由等高线生成dem的步骤 在arcgis中由等高线生成dem的步骤 1．进入arcgis的workstation模块 2．在Arc命令行下输入下面的命令（等高线的各层文件存放在el5目录中） Arc：arctin d:\el5 d:\tin line elev (黑色为提示符，蓝色为输入的命令，下同) （即为对el5 建立tin ，elev代表等高线的高程值，并且只有line 参与运算），这样就由等高线生成了tin 注：可以在Arc命令行设置workspace路径，以后的操作不必每次都有写上绝对路径，相对路径就可以了。命令为：Arc：wakespace d:\el5 Arc：w 可以显示当前系统的worksapce目录。 3．由tin生成lattice，需要输入如下命令 Arc：tinlattice d:\tin d:\lat 这样就有tin生成了lattice，转化为了grid形式，分辨率设置为30米 Enter distance between lattice mesh points :后要出入分辨率，对于其它的设置取其默认值即可。 4．最有一步，由lattice生成dem，命令如下 Arc：latticedem d:\lat dem 这样便由lattice转换得到了dem，运行结果如下： 然后在arcmap中可以打开生成的dem，同时也可以显示生成的tin，我们已经在d盘根目录建立了tin 和lattice子目录，目录名字分别为tin 、lat，在arcmap中会有对应的选项 分别单击之，则会加入到arcmap中成为图层。对显示tin的结果：