遥感变化监测实习报告.docx

1、遥感变化监测实习报告变化监测实习报告实习名称变化监测实习课程遥感图像处理XX班级实习时间学号得分实习原理：非监督分类运用1SODATA算法，完全按照像元的光谱特性进行统计分类，常常用于对分类区没有什么了解的情况。使用该方法时。原始图像的所有波段都参于分类运算，分类结果往往是各类像元数大体等比例。由于人为干预较少，非监督分类过程的自动化程度较高。非监督分类一般要经过以下几个步骤：初始分类、专题判别、分类合并、色彩确定、分类后处理、色彩重定义、栅格矢量转换、统计分析。实习数据： 遥感影像： LS5_TM_20100725_023435_023501_121040_FASTB_L2LS5_TM_20

2、081210_022812_022837_121040_FASTB_L2 简析：影像为江西省鄱阳湖地区，在影像生成时间内。时值夏/冬季，但江XX区植被多常绿。实习内容：就所下载遥感影像，采用非监督分类的方法，对影像中所放映的信息进行分类。实习目的：掌握非监督分类的方法与过程，加深对非监督分类方法的理解。实习步骤：第一步：调出非监督分类对话框 在 ERDAS 图标面板工具条中点击 Classifier 图标 C1assification Unsupervised Classification Unsupervised classification。对话框如下：第二步：进行非监督分类在Unsupe

3、rvised classification对话框输入数据（如上图右所示）。确定输出文件（Input Raster File）:caijianhou_40.img（要被分类的图像） 确定输出文件（Output File）:非监督分类_caijianhou 4001.img即将产生的分类图像） 选择生成分类摸板文件:非监督分类_caijianhou_4001.sig（将产生一个模板文件） 对Clustering options选择Initialize fr om Statistics单选框 Initialize f rom Statistics指由图像文件整体（或其AOI区域）的统计值产主自由聚类

4、，分出类别的多少由自己决定。Use Signature Means 是基于选定的模板文件进行非监督分类，类别的数目由模板文件决定。确定初始分类数（Number of classes）： 35分出35个类别 一般设定为最终所分类别的46倍定义最大循环次数（Maximum Iterations）： 24 最大循环次数（Maximum Iterations）是指ISODATA重新聚类的最多次数，这是为了避免程序运行时间太长或由于没有达到聚类标准而导致的死循环。一般在应用中将循环次数都取6次以上。 设置循环收敛阈值（Convergence Threshold）：0.95 收敛阈值（Convergenc

5、e Threshold）是指两次分类结果相比保持不变的像元所占最大百分之此值的设立可以避免ISODATA无限循环下去。 点击OK按钮（关闭Unsupervised Classification对话框，执行非监督分类，获得一个初步的分类结果）初步分类后，结果如下图：第三部 分类评价： 1，在视窗中同时显示 caijianhou_40.img（下）和 非监督分类_caijianhou 4001.img（上）。caijianhou_40.img 显示方式用红5、绿4蓝2。 2，打开分类图像属性并调整字段显示顺序。选择Raster菜单项，Rster-Attributes，打开Raster Attrib

6、ute Editor对话框（非监督分类_caijianhou 4001.img的属性表），并调整字段显示顺序： Raster Attribute Editor 对话框菜单条：Edit Column Properties column properties对话框如下图所示：在Columns中选择要调整显示顺序的字段，通过Up、 Down、Top、Bottom等几个按钮调整其合适的位置，通过选择Display Width调整其显示宽度，通过Alignment调整其对齐方式。如果选择Editable复选框，那么可以在Title中修改各个字段的名字及其它内容。3，给各个类别赋相应的颜色Raster A

7、ttribute Edi tor对话框，点击一个类别的Row字段从而选择该类别，右键点击该类别的Color字段（颜色显示区），打开As Is菜单并选择一种颜色。重复以上步骤直到给所有类别赋予合适的颜色 。操作后结果如下图所示：Viewer显示结果如下图： 4，不透明度设置。由于分类图像覆盖在原图像上面，为了对单个类别的判别精度进行分析，首先要把其它所有类别的不透明程度（Opacity）值设为0 （即改为透明），而要分析的类别的透明度设为1（即不透明）。5，确定类别专题意义及其准确程度 。 视窗菜单条：Utilityflicker viewer Flicker 对话框Auto Mode 本小步是

8、设置分类图像在原图像修背景上闪烁， 观察它与背景图像之间的关系从而断定该类别的专题意义，并分析其分类准确与否。6，标注类别的名称和相应颜色。Raster Attribute Editor对话框（打开）点击刚才分析类别的ROW字段从而选择该类别 点击该类别的class Names字段从而进入输入状态 在该类别的Class Names字段中输入其专题意义（如居民区建筑），并按回车键 右键点击该类别的Color字段（颜色显示区） As Is菜单选择一种合适的颜色重复以上4、5、6三步直到对所有类别都进行了分析与处理。非监督分类结果：7，分类后处理之聚类统计ERDAS图标面板工具条:点击Interpr

9、eter图标-GIS Analysis Clump-Clump对话框（如下图）在Clump对话框中，需要确定下列参数:1确定输入文件（Input File）:40_2010_uned.img2定义输出文件（Output File）:40_2010_uned_clump .img文件坐标类型（Coordinate Type） : Map处理范围确定（Subset Definition） : ULX / Y, LRX / Y （缺省状态为整个图像范围，可以应用Inquire Box定义子区）确定聚类统计邻域大小（Connect Neighbors）: 8，（统计分析将对每个像元四周的8个相邻像元进

10、行）OK（关闭Clump对话框，执行聚类统计分析）在virwer中打开40_2010_uned_clump .img 如下图8，分类后处理之去除处理ERDAS图标面板工具条:点击Interpreter图标-GIS Analysis-Eliminate-Eliminate对话框（如下图所示）在Eliminate对话框中，需要确定下列参数:1确定输入文件（Input File）: 40_2010_uned_clump .img2定义输出文件（Output File）: 40_2010_uned_eliminate .img文件坐标类型（Coordinate Tvpe ）: Map处理范围确定（Su

11、bset Definition）: ULX / Y, LRX / Y，（缺省状态为整个图像范围，可以应用InquireBox定义子区）确定最小图斑大小（Minimum）: 16 pixels确定输出数据类型（Output）: Unsigned 8 BitOK（关闭Eliminate对话框，执行去除分析）在viewer中打开40_2010_uned_eliminate .img，如下图所示8，分类后处理之分类重编码ERDAS图标面板工具条:点击Interpreter图标-GIS Analysis-Recode-Recode对话框。如下图所示1确定输入文件（Input File）:40_2010_

12、uned_eliminate.img2定义输出文件（Output File）: 3设置新的分类编码（Setup Recode ）:点击Setup Recode按钮4打开Thematic Recode表格5根据需要改变New Value字段的取值（直接输入）6在本例中将原来的35类依次两两合并，形成五类（图4.20）OK（关闭Thematic Recode表格，完成新编码输入）确定输出数据类型（Output）: Unsigned 4 BitOK（关闭Recode对话框，执行图像重编码。输出图像将按照New Value变换专题分类图像属性，产生新的专题分类图像）在视窗中打开重编码后的专题分类图像，

13、查看起分类属性表：改变颜色 9，精度评估点击Classifier弹出如下窗口：点击，打开重分类之后的图像。点击EditCreate/Random Point 设置参数。其中Number of Points 数值要大于256，勾选，其中Minimum Points 要大于20或30，最好设置为30，点击OK，在弹出窗口中一直点击是。 在弹出窗口中点击ViewSelect Viewer 之后点击一幅图像，重复操作一次，第二次点击另外一幅图像，使之两幅图像都显示待验证的点。之后点击ViewShow All 显示待验证的点。继续点击EditShow Class Values 显示原来的类别号。显示如下

14、窗口：之后找出每一个待验证的点进行验证，在Reference中输入确定后的类别号，直至验证完所有的点。 我在精度评估的时候由于验证一些点之后ERDAS软件总是崩溃，所以最后没有验证完成。点还没有校正完，软件就停在这里了。10，做缓冲区ERDAS图标面板工具条:点击Interpreter图标-GIS Analysis-Search-Search对话框。如下图所示1确定输入文件（Input File）:40_2010_uned_eliminate.img2定义输出文件（Output File）: 3Distance to search :=104OK5ERDAS图标面板工具条:点击Modeler图

15、标-Modeler Maker打开对话框； 如下图所示6建立模型7点OK44FileType选择Thematic8OK910分类改变颜色（建立缓冲区前后对比，左为建立缓冲区后影像）按相同方法对 LS5_TM_20081210_022812_022837_121040_FASTB_L2进行系列处理，用以进行变化检测。11，变化检测1ERDAS图标面板工具条:点击Modeler图标-Modeler Maker打开对话框； 如下图所示2建立模型3OK4在视窗中打开变化检测后图像5查看属性表6变化分析2008年数据2009年数据1234512584123610124095442270514021158474314362310869064516264732118210105292795655854771277407