实验八:栅格数据处理
一、实验目的
1、掌握栅格数据进行裁切、拼接、提取等操作; 2、掌握投影变换方法。
二、实验准备
数据准备:
1幅1:25万矢量数据(Vector),为白水县的行政范围。地理坐标系统,其中大地基准是D_North_American_1927,参考椭球体是Clarke 1866。
2幅1:25万DEM数据(DEM1和DEM2)。地理坐标系统,其中大地基准是D_Krasovsky_1940,参考椭球体是Krasovsky_1940。
软件准备:
ArcGIS Desktop9.x,ArcCatalog
三、实验内容
白水县跨两个DEM图幅,提取出白水县的DEM数据,并将数据转换成高斯克吕格投影系统。获取具有投影坐标系统的特定边界DEM数据。工作流程如图1所示。
2幅1:25万DEM数据 合并 1幅1:25万DEM数据 1:25万矢量数据 提取 白水县县界 裁切 白水县DEM数据 定义投影 结果数据 图1 工作流程
四、实验步骤
(1) 白水县行政范围的提取
1) 打开1:25万矢量数据(图2)。
图2 原始矢量数据
2) 利用Analysis Tools工具箱,Extract工具集中的Select工具,依据“name”字段,
即SQL表达式设置为“\白水县'”,提取出白水县地图数据(图3)。 A B
C
展开Analysis Tools工具箱,打开Extract工具集,双击Select,打开Select对话框。
在Input Features文本框中选择输入“E:/ChP4/Ex1/Vector”矢量数据。 在Output Feature Class文本框键入输出的数据的路径与名称“E:/ChP4/Ex1/vector_Select”。 D E
单击Expression可选文本框旁边的SQL表达式“\ = '白水县'”。 单击OK按钮,完成操作。
按钮,打开Query Builder对话框,设置
图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窗口选择该类型。
A 展开Data Management Tools工具箱,打开Raster工具集,双击Mosaic To New B C
Raster,打开Mosaic To New Raster对话框。
在Input Rasters 文本框中选择DEM1和DEM2。
在Output Location文本框键入输出数据存储的位置“E:/ChP4/Ex1”。
D E F G
在Raster dataset name with extension文本框设置输出数据的名称“DEM”。 在Pixel type可选窗口,设置输出数据栅格的类型为16_bit_UNSIGNED。 在Mosaic Mothod可选窗口,确定镶嵌重叠部分的方法,本次拼接方法选择MEAN,表示重叠部分的结果数据取重叠栅格的平均值。 单击OK按钮,完成操作。
图5 DEM拼接结果
(3) 白水县DEM的裁切
利用Spatial Analyst Tools工具箱,Extraction工具集中的Extract by Mask工具,以白水县矢量数据裁切拼接的DEM数据,获取白水县DEM(图6)。从图6下面状态栏的坐标可以看出,白水县DEM是以地理坐标系统显示的,为了便于量算以及与其它数据叠合分析,应该把地理坐标系统转换为投影坐标系统。
1) 展开Spatial Analyst Tools工具箱,打开Extraction工具集,双击Extract by Mask,
打开Extract by Mask对话框。
2) 在Input raster文本框中选择需要裁切的栅格数据“E:/ChP4/Ex1/DEM”。
3) 在Input raster or feature mask data文本框定义进行裁切数据“E:/ChP4/Ex1/ vector_Select”。
4) 在Output raster文本框键入输出数据的路径与名称“E:/ChP4/Ex1/extract_dem”。 5) 单击OK按钮,完成操作。
图6 白水县DEM数据
(4) 白水县DEM的投影变换
利用Data Management Tools工具箱,Projections and Transformations中的Raster工具集,选择其中的Project Raster工具,进行白水县DEM的投影变换。我国大中比例尺地形图规定采用以克拉索夫斯基椭球体元素计算的高斯—克吕格投影。因此,投影方式选择Xian 1980 GK Zone 19.prj(图7),即为高斯—克吕格投影,西安1980大地基准,中央经线为111°。转换结果如图8所示,从下面状态栏显示的坐标可以看出,白水县的DEM已转换为平面直角坐标。
图7 投影坐标系统的选择
1) 展开Data Management Tools工具箱,打开 Projections and Transformations中的Raster工具集,双击Project Raster,打开Project Raster对话框。 2) 在Input raster文本框中选择进行投影变换的栅格数据“E:/ChP4/Ex1/extract_dem”。 3) 在Output raster文本框键入输出的栅格数据的路径与名称E:
/ChP4/Ex1/Result/project_dem”。 4) 单击Output coordinate system文本框旁边的
图标,打开Spatial Reference属性对
话框,单击Select按钮,打开Browe for Coordinate System对话框,选择Xian 1980 GK Zone 19.prj投影。
5) Resampling technique是选择栅格数据在新投影类型下的重采样方式,选择
NEAREST。
6) 单击OK按钮,完成操作。
图8 投影结果图
百度搜索“70edu”或“70教育网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,70教育网,提供经典综合文库Arcgis操作 实验八:栅格数据处理在线全文阅读。
相关推荐: