10GIS_2010210618综合实验：基于ERDAS和ArcGIS的区域土地利用变更分析.doc

综合实验：基于ERDAS和ArcGIS的区域土地利用变更分析 实验任务：根据任务分配表获取自己所需要的数据源，应用ERDAS和ArcGIS软件将本人所分配城市范围内城镇用地和水域面积进行统计并作对比，如表1，并制作城市土地变化图，参见“实验四：保定市土地利用变化分析.pdf”文档30页。 表1 土地类型 2004（km2） 2004像元数 2009（km2） 2009像元数 建设用地 170.8883 189876 261.8406 290934 水域 65.0286 72254 65.0898 72322 其他 511.1253 567914 420.1128 466792 最终提交该实验的实验报告一份（word文档格式，要求图文并茂），其中包括该实验的操作流程和实验结果（实验结果主要是表1和城市土地变化图）。 实验报告的命名格式为：姓名_城市名 实验大致步骤： 图像增强处理 辐射增强处理： 辐射增强（Radiometric Enhancement） 辐射增强处理是对单个像元的灰度值进行变换达到图像增强的目的。 ERDAS IMAGINE供的辐射增强处理功能如表4.2所列。 表4.2 遥感图像辐射增强命令及其功能 直方图均衡化处理 该处理实质上是对图像进行非线性拉伸，重新分配图像像元值，使一定灰度范围内的像元的数量大致相等。 .在ERDAS图标面板工具条中，点击Interpreter/Radiometric Enhancement/Histogram Equalization—--打开Histogram Equalization对话框,并设置参数如下： 打开两个新窗口，分别将2004年原数据和已增强的数据打开，进行对比，如图所示： 如下图是2009年的原数据和已增强的数据对比： 图像裁切 AOI多边形裁剪（Polygon AOI Subset Image ） 首先在视窗中打开需要裁剪的图像，并应用AOI工具绘制多边形AOI：可以将多边形AOI保存在文件中（（*.aoi ），也可以暂时不退出视窗，将图像与AOI多边形保留在视窗中，然后 ERDAS图标面板菜单条Main-Data Preparation-Data Preparation菜单，选择Subset Image-打开Subset Image对话框ERDAS图标面板工具条:点击Data Prep图标-打开Data Preparation菜单 （3）选择Subset Image-打开Subset Image对话框，在Subset Image对话框中需要设置下列参数： 输入文件名称（Input File）；输出文件名称（Output File）；（4）应用AOI确定裁剪范围点击AOI按钮打开选择AOI （Choose AOI）对话框 （5）在Choose AOI对话框中确定AOI的来源（AOI Source） : File（或Viewer），如果选择了文件（File），则进一步确定AOI文件，否则，直接进入下一步 （6）输出数据类型（Output Data Type）:Unsigned 8 Bit （7）输出像元波段（Select Layers）:1:6（表示选择1,2, 3, 4, 5,6六个波段） （8）OK（关闭Subset Image对话框，执行图像裁剪）在Viewer视窗中打开剪裁后的图像如图所示 在Viewer视窗中打开剪裁后的图像如图所示 监督分类 3.1设置监督分类特征文件。 1.点击主菜单上的Classifer——Signature Editor，跳出Signature Editor面板。 2.点击主窗口菜单栏上的AOI——Tools，用多边形工具画出一类用地的训练区，点击加入特征文件编辑面板。 3.画一块，添加一块，每类用地在重点地域和全区域选5-10块。点击，可将选择的几块地归为一类。归并后，点击Edit——delete将原有的删除。并重命名地类。 4.将全区域土地分为几种地类。保存分类特征信息：在Signature Editor面板上，点击File——Save As…，在弹出的对话框中选择路径并为文件命名，点击OK。 3.2利用分类特征文件进行监督分类。点击Classifier——Supervised Classification，在弹出的对话框选择要分类的img图像（Input Raster File），选择特征文件（Input Signature File），选择分类后的文件路径，并命名（Classied File）。在Non-parametric Rule中选择Feature Space，其他按默认设置。 3.3合并地类。 1.打开经过监督分类解译后的2004年遥感影像：在主窗口上点击，打开。 2. 点击主菜单上的Inerpr