DEM数据处理与分析Word格式文档下载.docx

1、其次，怎么去除高程影像中的空值（如-32767），让它在一个正常的范围内显示呢？小编这里以TIF格式的DEM高程影像为例（IMG的处理方式同样），一步步带大家来操作。1. 在ARCMAP里打开一幅DEM高程数据（ADD DATA），可以从左边看到其显示的数据范围是-32767-32726，右侧为灰色影像。2. 在ArcMap里打开Spatial analyze工具，选择Raster Calculator，设置DEM高程数据值为0并进行计算（点击Evaluate按钮），页面如下：3. 计算之后结果显示为Calculation，页面如下图：4. 再次打开Spatial analyze工具，选择Re

2、classify，选中chang missing values to Nodata打钩, 点击OK。5. 然后打开Raster Calculator面板，将上一步计算出来的Reclass of Calculation与原始的DEM 高程数据相乘, 点击Evaluate，界面如下：6. 计算结果，如图：7. 鼠标右键点击Calculation2，打开Properties, 选择Symbology选项，在左侧选中Stretched, 点击确定。如下图：8. 在这里，我们看到，经过处理，高程影像的值域范围已经变为正常的海拔范围（105-4336）了：9. 我们还可以给它设一个漂亮的显示颜色，简单的双

3、击左侧的灰度条即可，到这一步，DEM数据的预处理就做完了。（二）其他处理1. 首先添加DEM数据：2. 如下图显示，全黑，是由于存在坏点，可使用Calculate Statistics处理一下。3. 如下图依次展开Arc ToolboxData Management ToolsRasterRaster PropertiesCalculate Statistics4. 点击Calculate Statistics，打开出现下图选择所要处理的DEM数据。5. 点击OK，进行计算，计算完成如下图。（三）坐标转换（计算坡度之前的预处理）下载的DEM数据是WGS-1984坐标系，属于地心坐标系，XY坐标

4、单位是度分秒，必须将其单位转换为米，否则将导致所得坡度不正确，坡度平均会在80-90度（如下图）。故将其转换成以米为单位的坐标系，北京1954坐标系。1. 输入Project Raster，点击后出现下图：2. 在Out Coordinate System中点击后面的按钮，选择Gauss Kruger的Beijing 1954坐标系。根据数据所在六度分带选择投影方式。点击确定。3. 在Geographic Transformation中选择一种。（此处不确定都是什么意思）4. 点击OK开始处理，处理完成坐标转换成功。三、DEM数据拼接DEM数字高程包括SRTM90米分辨率原始高程数据和30米分

5、辨率数字高程数据产品。SRTM90米分辨率原始高程数据、30米分辨率数字高程数据产品。其中，SRTM由美国太空总署（NASA）和国防部国家测绘局（NIMA）联合测量。2000年2月11日，美国发射的“奋进”号航天飞机上搭载SRTM系统，SRTM系统获取的雷达影像的数据量约9.8万亿字节，经过两年多的数据处理，制成了数字地形高程模型（DEM），即现在的SRTM地形产品数据。如果下载30米的DEM数据，其文件为压缩文件，每个压缩包内包含三个文件，即数字高程模型（_dem）文件、质量评估（_num）文件和高程快视图（_jpg）文件。（一）获取打开ArcMap，将所有涉及到的数字高程模型文件加载到Ar

6、cMap中。（婺城区共两个文件，下图以四川省为例，共六个文件）。问题：为什么9个DEM数据的相接处的显示有很明显的分界？（二）镶嵌镶嵌是指两个或多个图像的组合或合并。在ArcGIS中，您可以通过将多个栅格数据集镶嵌到一起来创建一个单个栅格数据集。此外，还可以通过一系列栅格数据集创建镶嵌数据集和虚拟镶嵌。使用ArcToolbox中Data Management Tools（数据管理工具）/Raster（栅格）/Raster Dataset（栅格数据集）/Mosaic（镶嵌）或Mosaic To New Raster（镶嵌至新栅格）进行拼接。如运行“镶嵌至新栅格”，输入下载的多个DEM数据、输出的

7、位置（已经存在的目录）和合并后的栅格数据文件名，波段数据设置为1。如果选择“镶嵌”，则合并到第一个栅格文件中。镶嵌完成，得到如下结果，并把其它的图层移除镶嵌后是什么坐标系统？共有多少行和列？（三）裁剪合并后的DEM数据可能与研究区域不重合，需要根据边界对DEM数据进行裁剪，如将上面拼接的结果，按四川省边界图层进行提取。1. 加入四川省边界图层。2. 使用ArcToolbox/Spatial Analyst Tools（空间分析工具）/Extraction（提取分析）/Extract By Mask（按掩膜提取）提取工具。提取的结果文件命名为四川省90米DEM。（数据将用于后续实验任务）。结果如

8、下：ArcGis中除了可以按掩膜提取外，还可以按什么方式提取？四、地形属性提取（一）坡度提取ArcGis中坡度工具用于为每个像元计算值在从该像元到与其相邻的像元方向上的最大变化率。实际上，高程随着像元与其相邻的八个像元之间距离的变化而产生的最大变化率可用来识别自该像元开始的最陡坡降。该工具会将一个平面与要处理的像元或中心像元周围一个3x3的像元邻域的z值进行拟合。该平面的坡度值通过最大平均值法来计算（请参阅参考书目）。该平面的朝向就是待处理像元的坡向。坡度值越小，地势越平坦；坡度值越大，地势越陡峭。1.启动ArcMap，添加数据框，并更名为“任务2”，将四川省90米DEM（任务1得到的结果栅格

9、图层）加入。2.使用ArcToolbox中Spatial Analyst Tools（空间分析工具）/Surface（表面分析）工具集中Slope（坡度）工具提取坡度。怎样把坡度25度以上的地方查找出来？（二）坡向提取坡向用于识别出从每个像元到其相邻像元方向上值的变化率最大的下坡方向。坡向可以被视为坡度方向。输出栅格中各像元的值可指示出各像元位置处表面的朝向的罗盘方向。将按照顺时针方向进行测量，角度范围介于0（正北）到360（仍是正北）之间，即完整的圆。使用ArcToolbox中Spatial Analyst Tools/Surface工具集中Aspect（坡向）工具提取坡向。提取的结果如下：

10、（三）表面曲率提取曲率工具会逐个像元地计算输入表面的二阶导数值。曲率工具的输出是该表面的二阶导数（例如，坡度的坡度），从应用的角度看，该工具的输出可用于描述流域盆地的物理特征，从而便于理解侵蚀过程和径流形成过程。坡度会影响下坡时的总体移动速率。坡向将决定流向。剖面曲率将影响流动的加速和减速，进而将影响到侵蚀和沉积。平面曲率会影响流动的汇聚和分散。使用ArcToolbox中Spatial Analyst Tools空间分析工具/Surface表面分析工具集中Curvature工具提取表面曲率。通过点击图层右键，选择属性中的符号系统，显示方式设置为已分类，并设置分类中的类别数量及中断值，将其显示为

11、0.5总4个级次。五、透视图建立（一）设置抬升高度打开ArcScene，将场景名改为“任务3”，并将四川省90米DEM数据加载到视图中。并通过Properties（右键属性）中BaseHeights（基本高度）栏中，选择 Obtain heights for layer from surface，选择DEM数据为地表抬升高度数据。转换系数设置为8，表示抬升高度（米）是实际高度的8倍。（二）修改显示符号系统在符号系统中，选择拉伸方式显示，并设置色带为Elevation #1。（三）设置渲染在Rendering（渲染）栏中进行渲染，设置光照效果。（四）其它图层（栅格或矢量）数据按地形高度进行抬升将

12、四川省地州界及四川省主要公路图层加入场景。并设置四川省地州界及四川省主要公路图层的基本高度参数。附加：将四川省气象站点的温度插值栅格图层（实验五的结果）加在场景中，并以DEM为基础高程抬升（结果见下图）。六、建立和显示TINArcGis中栅格转 TIN 工具的用途是创建表面偏离输入栅格不超过指定 Z 容差的不规则三角网 （TIN）。 （一）TIN转换1. 启动ArcMap，添加数据框，并更名为“任务4”，将四川省90米DEM（任务1得到的结果栅格图层）加入。2. 使用ArcToolbox中3D Analyst工具/转换工具/栅格转TIN进行转换。（二）TIN属性描述1. z_tolerance

13、（可选），输入栅格与输出 TIN 之间所允许的最大高度差（z 单位）。默认情况下，z容差是输入栅格 z 范围的 1/10。如果指定的 Z 值容差较小或者栅格表面在地形构造方面很粗糙复杂，则需要更多的点来构建 TIN。2. max_points（可选），将在处理过程终止前添加到 TIN 的最大点数。默认情况下，该过程将一直持续到所有点被添加完。如果已指定，则点的最大数量将作为 TIN 的大小限制。如果在 TIN 结点计数达到此限制时，仍不符合 Z 容差，则该工具将会停止，并返回操作失败。3. z_factor（可选），输入 TIN 的高度转换为输出 TIN 的高度时所乘的系数。用于将 z 单位转

14、换为 x 和 y 单位。（三）TIN渲染修改Z容差值及最多点数，查看转换后的TIN组成。七、创建等高线（一）创建等高线1. 启动ArcMap，添加数据框，并更名为“任务5”，将四川省90米DEM（任务1得到的结果栅格图层）加入。2. 使用ArcToolbox中Spatial Analyst Tools空间分析工具/Surface表面分析工具集中等值线Contour工具提取表面等值线。也可以使用ArcToolbox中3D Analyst工具/栅格表面/等值线工具。注意：路径、文件名最好用英文字母，不要太深、不要有特殊字符，比如括号么、“-”等。起始等值线可用作起点或参考，其他所有等值线（包括正等

15、值线和负等值线）都将从该条等值线派生出来。起始等值线的默认值为 0.0。如果输入 5.0 作为等值线间距，则可生成 0、5、10 等正值，以及 0、-5、-10 等负值。实际生成的等值线将取决于此间距模式与输入 terrain 数据集或 TIN 的实际 z 范围间距模式的叠加。如果将起始等值线值指定为 0.1，使用相同的间距 5.0，则等值线模式将变为 0.1、5.1、10.1 等正值以及 -4.9、-9.9、-14.9 等负值。3. 标注等高线值。使用在属性选项中标注等高线的值，选择Contour为标注字段。（二）创建垂直剖面1. 启动ArcMap，添加数据框，并更名为“任务6”，将四川省9

16、0米DEM（任务1得到的结果栅格图层）加入。加入四川省主要公路图层。2. 打开3D Analyst工具栏。选择3D Analyst工具栏中的“插入线Interpolate Line”功能，使用鼠标沿一段主要公路划出创建垂直剖面图经过的点的线。 3. 使用3D Analyst工具栏中的创建剖面图Create Profile Graph工具，创建沿划线经过点的剖面图。通过属性和高级属性修改显示样式。（三）坡度分级分析四川省坡度小于等于10、10至20、大于20度各类点的百分之，并估算其面积。1. 启动ArcMap，添加数据框，并更名为“任务7”，将四川省坡度图层（任务2提取的坡度图层）加入。2.

17、使用重分类功能，分为三类，将中断值改为10、20。3. 通过属性表进行计算三类中的点数。面积怎样估算？七、DEM相关应用DEM应用之坡度：Slope新建地图文档，加载【ArcGIS地形分析-TIN及DEM的生成，TIN的显示】经验教程中得到的DEM数据：TINGrid。在【ArcToolbox】中，执行命令3D Analyst工具栅格表面坡度， 参照下图所示，指定各参数：执行后，得到坡度栅格Slope_tingri1：坡度栅格中，栅格单元的值在0 -82 度间变化。右键点击图层Slope_TinGrid，执行属性命令，设置图层符号系统，重新调整坡度分级。将类别调整为5，点分类按钮，用手动分级法

18、，将中断值调整为：8，15，25，35，90。【下面计算剖面曲率】在【ArcToolbox】中，执行命令3D Analyst工具栅格表面坡度。按如下所示，指定各参数。得到剖面曲率栅格：Slope _SlopeDEM应用之坡向：Aspect打开【ArcToolbox】，执行命令3D Analyst工具栅格表面坡向，按下图所示，指定各参数：执行结果为（得到坡向栅格：Aspect_tgrid）：【以下计算平面曲率】：在【ArcToolbox】中，执行命令3D Analyst工具栅格表面坡度，按下图所示，指定各参数，按下图所示指定各参数：执行后生成平面曲率栅格Slope_Aspect：DEM应用之提取

19、等高线新建地图文档，加载DEM数据：tingrid。（在执行以下操作时确保，3D Analyst扩展模块已激活，以后不说表示本人已经激活，各位在操作的时候这点应该也懂的吧！）然后在【ArcToolbox】中，执行命令3DAnalyst工具等值线，按下图所示指定各参数：执行后生成等高线矢量图层：Contour_tingrid:DEM应用之计算地形表面的阴影图在【ArcToolbox】工具箱中，执行命令3D Analyst工具 山体阴影，按下图所示指定各参数：执行后生成地表阴影栅格：Hillsha_ting1:下面进行【DEM渲染】：关闭除tingrid 和 Hillsha_ting以外所有图层的

20、显示，并将tingrid 置于Hillsha_ting 之上，右键点击tingrid ，在出现的右键菜单中执行属性，在图层属性对话框中，参照下图所示设置符号系统选项页中颜色。如图：在工具栏空白处右键，打开工具栏效果，如下图所示，设置栅格图层tingrid的透明度为：45%左右，结果如图：DEM应用之可视性分析（一）【通视性分析】数据使用上一步骤的结果！在上一步的基础上进行。打开3D Analyst 工具栏，从工具栏选择创建通视线（Line of sight）工具，如图：在出现的通视分析对话框中输入观察者偏移量和目标偏移量, 即距地面的距离，如图：在地图显示区中从某点A沿不同方向绘制多条直线，可

21、以得到观察点A到不同目标点的通视性：（二）【可视区分析：移动发射基站信号覆盖分析】（1）在上一步基础上进行，在内容列表区TOC中关闭除tingrid之外的所有图层，加载移动基站数据矢量图层：移动基站.shp（数据后边会提高给大家下载） （2）在【ArcToolbox】中，执行命令3D Analyst工具视域，按下图所示指定各参数：（其中绿色表示现有发射基站信号已覆盖的区域，淡红色表示，无法接收到手机信号的区域）DEM应用之地形剖面本步骤同样是在上一步骤的结果下进行的：在上一步基础上进行，打开3D Analyst工具栏，点击线插值工具，跟踪一条线段，这条线段可以从DEM:TINGRID 中得到高

22、程值点击剖面图 按钮，得到上一步所生成的3D线段的剖面图：到此，DEM应用除了【水文分析】外，主要的应用基本讲完。由于每个经验教程只能上传20张图片，故【水文分析】将会在下一个教程讲授！八、说明本文主要通过收集、整理并汇总网络上关于DEM数据处理与分析的各种实验方案及教程，形成了较为完整、细致的文档。汇总的内容主要来源有：1.如何在地理空间数据下载免费DEM数据？2.DEM高程数据处理方法吴飞3.DEM做坡度（详细步骤）及常见问题解决方法4.成都信息工程学院 地理信息系统原理实验指导 实验六、DEM的建立与应用王增武5.ArcGIS教程之DEM的应用（坡度坡向、提等高线）XX经验其他相关深入实习：1.基于ArcGIS的DEM信息提取及应用分析2.实习七、DEM建立与应用