ERDAS实验报告

遥感实验四图像增强和利用监督分类提

取地类信息

一、实验目的

根据任务分配表获取自己所需要的数据源，应用ERDAS和ArcGIS软件将本人所分配城市范围内城镇用地和水域面积进行统计并作对比, 并制作城市土地变化图。

二、实验数据

tm12103820041014.img、tm12103820090604.img、new_boundry.shp

三、实验内容及主要步骤

1.图像增强(主成分变换的方法)：打开Erdas软件，点击Interpreter->Spectral Enhancement->principal Comp，选择2004年数据为要进行增强的影像，选择输出文件路径和名称，并设置参数如下表，点击OK。

图1

图2

2.裁剪（AOI视图裁剪方法）：点击Data Prep->Subset Image，弹出Subset对话框，选择要输入和输出的文件路径和名称，打开视图窗口，加载new_boundry.shp，点击合肥市所在区域，点击AOI->Copy Selection to AOI，点击Subset对话框中的AOI按钮，选择Viewer，

点击OK。

图3 图4

图6

3.监督分类

i点击Classifier->Signature Editor，弹出Signature Editor对话框，点击视图窗口

AOI->Tools。利用按钮在裁剪的图上勾选出感兴趣区域，点击按钮添加到Signature Editor对话框中。添加9个水域的感兴趣区域后，选择所有添加的记录，点击Signature Editor

对话框将上面的记录归结为一类，并命名为water，点击Edit->Delete，删除选中的记录。同理，选择出绿地和城镇的分类信息。点击File->Save as，存储为*sig格式的分类信息文件。

ii点击Classifier->Supervised Classification，选择要进行分类的图像和分类后的图像名称和路径，选择之前做好的*sig格式分类信息文件导入，点击OK。

图7

4.合并地块

i点击Tnterpreter->GIS Analysis->Record，弹出Recode对话框。选择分类好的影像，选择输出的文件路径和名称，勾选lgnore Zero in Stats，点击Setup Recode，将同一地类的值设为相同值。点击Thematic Recode的OK按钮。点击Recode上的OK按钮。

图8

5.合并零碎图斑

i Clump聚合：点击Inerpreter->GIS Analysis->Clump，弹出Clump对话框，选择合并地类后的影像，选择输出的影响的路径和名称。勾选lgnore Zero in Output Stats。

图9

ii Eliminate融合：点击Inerpreter->GIS Analysis->Eliminate，在弹出的对话框选择输入的聚合后的影像和输出的影像路径和名称。在Minimum输入要融合的最小象元数，勾选lgnore Zero in Stats，点击OK。

图10

6.建模二值化

i 打开进行融合过的影像，点击Raster->Profile Tool，弹出Select Profile Tool对话框，选择Spatial，点击OK。利用SPATIAL PROFILE面板上的工具查看不同颜色的象元值。

其中白色区域为3，亮灰色区域为2，暗灰色区域为1。

图11

ii 点击Modeler->Model Maker，画好建模图，输入要建模的影像和建模二值化后的影像路径和名称。双击函数定义区域，弹出Function Definition，在Function下拉选择Conditional，具体选择EITHER ..IF.. OR..函数，在下面的函数表达式框中写入EITHER 1 IF ( $n1_quchu=3 ) OR 0 OTHERWISE，表示这个影像值为3的区域全为1，其他全为0，以此进行二值化。点击OK。点击按钮运行函数。这里已经将白色区域即城镇进行二值化。同理，可以将

水域和绿地进行二值化。

图12

图13城镇二值化图14水域二值化图15绿地二值化

7.栅格图转成矢量图

点击Vector->Raster to Vector，弹出Raster to Vector对话框，将已经二值化的影像输入，并输出转成矢量图的路径和名称。

8.分类变化成图

i 打开ArcGIS->ArcMap，点击工具栏的按钮添加之前转成矢量的数据。在数据管理窗口右击一个图层->Open Attribute Table，打开这个图层的属性表，在GRID-CODE一列右击Sort Ascending进行升序排列。点击Option->Select by Attribute，弹出Select by Attribute 对话框，双击GRID-CODE，点击“=”，在SELECT *FROM polygon WHERE框中的等号后面输入1。点击Apply。关闭对话框。右击属性表中GRID-CODE列，选择Field Calculate，弹出对话框，在GRID-CODE = 框中输入2。点击OK。同理更改另外两幅矢量图的GRID-CODE属性值为3和4。

图16 图17 ii 右击图层管理窗口的一个图层，选择Data->Export Data ，弹出对话框，选择输出文件的路径，这样就可以将*ArcInfo 格式转成*shp 格式。点击OK ，并在到地图窗口中。

图18

iii 点击ArcMap 工具栏上ArcToolbox 按钮，打开ArcToolbox ，选择Data Mangement Tools->General->Merge （合并工具），双击Merge ，弹出对话框。在Input Datasets 下拉框中选择三个*shp 格式的文件，合成一个shp 格式的文件。在Output Datasets 选择文件输出路径，点击OK 。合并后的结果为图20。