4数据的处理与变换.docxVIP

练习：数据更新变换 1． 背景：由于空间数据（包括地形图与 DEM）都是分幅存储的，造成某一特定研究区域跨越了不同的图幅。而当我们要获取有特定边界的研究区域时，就要对数据进行裁切、拼接(overlay/union)、提取（extract/select）等操作，有时还要进行相应的投影变换。 2． 目的：获取具有投影坐标系统的特定边界的 DEM 数据。 3． 要求：通过两幅给定的 DEM 数据，提取出白水县县界范围内的 DEM 数据，并将数据转换成高斯克吕格投影系统。通过练习，掌握数据提取、裁切、拼接及投影变换的方法。 4． 数据：1 幅 1：25 万矢量数据，为地理坐标系统，其中大地基准是D_North_American_1927,参考椭球体是 Clarke 1866，这是 Arcgis 为 Shapefile 类型的数据假设的地理坐标系统（实验数据在 Chapter4/Exercise1 中）。2 幅 1：25 万 DEM 数据，为地理坐标系统，大地基准是 D_Krasovsky_1940,参考椭球体是 Krasovsky_1940。 1：25万矢量数据 白水县县界 提取 裁切 2幅1：25万DEM数据 合并 1幅1：25万DEM数据 白水县DEM数据 定义投影 结果数据 图1 工作流程 5． 操作步骤： （1）白水县的行政范围的提取 1） 打开 1：25 万矢量数据（图 2）。 图 2 原始矢量数据 2） 利用 Analysis Tools 工具箱，Extract 工具集中的 Select 工具，依据“name”字段，即 SQL 表达式设置为“NAME = 白水县”，提取出白水县（图 3）。 展开Analysis Tools工具箱，打开Extract工具集，双击Select，打开Select 对话框。 在 Input Features 文本框中选择输入“E：/ChP4/Ex1/Vector”矢量数据。 在 Output Feature Class 文本框键入输出的数据的路径与名称“E：/ChP4/Ex1/vector_Select”。 单击Expression可选文本框旁边的按钮，打开Query Builder对话框，设置SQL表达式“NAME = 白水县”。 单击OK按钮，执行Select操作。 图 3 白水县矢量数据 （2）DEM 数据的拼接 1） 打开白水县横跨的两幅 DEM 数据，DEM1 和 DEM2（图 4）。 图 4 原始 DEM 数据 2） 利用 Data Management Tools 工具箱，Raster 工具集中的 Mosaic To New Raster 工具进行数据拼接，得到拼接结果 DEM（图 5）。其中，Cellsize 是默认状态，则结果数据与原始数据的栅格大小一致。因为原始数据类型为 16_BIT_UNSIGNED，所以在 Pixel type 窗口选择该类型。 展开 Data Management Tools 工具箱，打开Raster工具集，双击 Mosaic To New Raster，打开 Mosaic To New Raster 对话框。 在 Input Rasters 文本框中选择输入 DEM1 和 DEM2。 在 Output Location 文本框键入输出数据存储的位置“E：/ChP4/Ex1”。 在 Raster dataset name with extension 文本框设置输出数据的名称“DEM”。 E. 在 Pixel type 可选窗口，设置输出数据栅格的类型为 16_bit_UNSIGNED。 在 Mosaic Mothod 可选窗口，确定镶嵌重叠部分的方法，本次拼接方法选择MEAN，表示重叠部分的结果数据取重叠栅格的平均值。 单击OK按钮，执行 Mosaic To New Raster 操作。 图 5 DEM 拼接结果 （3）白水县 DEM 的裁切 利用 Spatial Analyst Tools 工具箱，Extraction 工具集中的 Extract by Mask 工具，以白水县矢量数据裁切拼接的 DEM 数据，获取白水县 DEM（图 6）。从图 6 下面状态栏的坐标可以看出，白水县的 DEM 是以地理坐标系统显示的，为了便于量算以及与其它数据叠合分析，应该把地理坐标系统转换为投影坐标系统。 展开 Spatial Analyst Tools 工具箱，打开 Extraction 工具集，双击 Extract by Mask，打开 Extract by Mask 对话框。 在 Input raster 文本框中选择输入需要裁切的栅格数据“E：/ChP4/Ex1/DEM”。 在 Input raster or feature mask d