地理信息系统基础技术-显示和分析DEM.docVIP

显示和分析DEM 在已经生成的DEM基础上，我们来进行一些具体应用和分析。因为我们制作出的DEM的坐标单位为度（DD），而在应用一些算法时，如计算山影和坡度时，算法要求地图单位为米，所以我们首先将DEM转换为UTM投影系。 添加地理投影属性 打开一个浏览窗口（Viewer），显示DEM，如下图。在浏览窗口的工具栏中点击图层信息（Layer Info）的按纽。 将打开一个显示图层信息的窗体，从中我们可以看到许多有用的信息。目前在地理投影信息中，还几乎是空白。 从菜单编辑（Edit）中选择改变地图模型（Change Map Model）。 在弹出的窗口中的投影（Projection）后指明为经纬度（Geographic Lat/Lon），单击 OK回到前一屏。 在菜单编辑（Edit）中选择添加/改变投影系（Add/Change Projection） 在新弹出的窗体中指明投影系为经纬度（Geographic Lat/Lon），椭球体和光滑化的名称都选WGS 84。单击OK回到上一屏。 现在我们可以从图层信息（Layer Info）窗体中看出DEM已经有了地理投影属性。 进行坐标转换 现在我们将经纬度的地理坐标转为UTM的。 单击功能面板上的数据前处理（DataPrep）按纽，在新出现的功能列表中单击按纽图象重投影（Reproject Image）。 在新出现的坐标转换（Reproject Images）窗口中指定以下参数。指明输入文件名（Input File,），输出文件名（Output File），改变目标投影为UTM WGS 84 North，如图。输出栅格大小为30M，也可以做的更小一点，其他参数可以不变，单击OK，开始进行坐标转换。 将出现一个工作进度指示器，计算完成后单击OK退出。 打开一个浏览窗口显示进行过坐标变换的DEM，可以在浏览窗口下面的状态栏中看到坐标已经发生了改变。同时图象有一些倾斜。 在ArcView中显示DEM 打开ArcView，从菜单文件（File）选择扩展模块（Extensions）。 在新出现的加载扩展模块窗体中浏览到空间分析模块（Spatial Analyst），单击前面的小方块，添加一个对勾，单击OK 从视图菜单文件（File），选择设置工作目录（Set Working Directory） 在新出现的对话框中设置工作目录。 将我们前面生成的DEM文件，“class_dem_utm.img”，添加到视图（View）中 选中该图层，从菜单图层（Theme）中选择转换为GRID格式（Convert to Grid） 在随后出现的窗体中浏览到存放目标文件的文件夹，指定文件名，并将新生成的GRID文件加入视图中。GRID格式是ArcView 的空间分析模块可以处理的格式。 选中Utm-dem的GRID图层，从菜单表面（Surface）选择计算山影（Compute Hillshade ） 在随后弹出的山影计算参数选单中可以填写一些计算中用到的参数。这两个参数指明了光源方位角（Azimuth）和高度角（Altitude）。 计算完成后生成了一个临时的GRID文件存放山影数据。现在我们需要永久的存储它。选中该图层，从视图菜单图层（Theme）中选择转换为GRID格式（Convert to Grid）。 在随后出现的文件存储对话框中指明新文件的位置和名称。 然后我们需要清除一些已经用不到的图层。选中这些图层，从菜单编辑（Edit）选删除图层（Delete Themes），删除它们，只留下DEM和山影图层。该步骤不会永久的删除文件，而只是将图层从显示的列表中删除。 现在我们改变DEM的显示。双击DEM图层，在随后弹出的图例编辑器（Legend Editor）中点击分类（Classify）。 在随后出现的分类参数窗口中指定分类种类（Type）为自然间隔（Natural Breaks），分类数（Number of classes）为15。单击OK回到图例编辑器。 在图例编辑器的色带（Color Ramps）下拉选单中选择海拔带谱1（Elevation #1）。 点击在图例编辑器左下角的高级选项（Advanced），在随后弹出的对话框中指定亮度图层（Brightness Theme）为我们刚刚生成的山影图层。这一操作将在显示DEM图层时参考山影图层，对每一个点显示的亮度做调整，使得整个图看上去有立体感。 指定最小亮度（Minimum Cell Brightness）为50，最大亮度（Maximum Cell Brightness）为100。单击OK。关闭图例编辑器。 可以看到DEM的显示方式发生了变化。 目前每个图层都有许多图例数据，我们可以隐藏它们。选中需要隐藏图例的图层，从菜单图层（Theme）中选择隐藏/显示图