遥感 - 变化监测实习报告 - 图文

更新时间:2023-12-30 09:13:02 阅读量: 教育文库 文档下载

说明:文章内容仅供预览,部分内容可能不全。下载后的文档,内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的,是否完整无缺。

遥感影像非监督分类实习报告

变化监测实习报告

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

遥感影像: LS5_TM_20100725_023435_023501_121040_FASTB_L2

LS5_TM_20081210_022812_022837_121040_FASTB_L2

简析:影像为江西省鄱阳湖地区,在影像生成时间内。时值夏/冬季,但江西地区植被多常绿。

实习内容:就所下载遥感影像,采用非监督分类的方法,对影像中所放映的信息进行分类。

实习目的:掌握非监督分类的方法与过程,加深对非监督分类方法的理解。 实习步骤:

第 1 页 共 1 页

遥感影像非监督分类实习报告

第一步:调出非监督分类对话框

在 ERDAS 图标面板工具条中点击 Classifier 图标

→C1assification →Unsupervised Classification → Unsupervised classification。对话框如下:

第二步:进行非监督分类

在Unsupervised classification对话框输入数据(如上图右所示)。 确定输出文件(Input Raster File):caijianhou_40.img(要被分类的图像) →确定输出文件(Output File):非监督分类_caijianhou 4001.img即将产生的分类图像)

第 2 页 共 2 页

遥感影像非监督分类实习报告

→选择生成分类摸板文件:非监督分类_caijianhou_4001.sig(将产生一个模板文件)

→确定分类摸板文件(Filename ):非监督分类_caijianhou_4001.sig →对Clustering options选择Initialize fr om Statistics单选框 Initialize f rom Statistics指由图像文件整体(或其AOI区域)的统计值产主自由聚类,分出类别的多少由自己决定。Use Signature Means 是基于选定的模板文件进行非监督分类,类别的数目由模板文件决定。 →确定初始分类数(Number of classes): 35分出35个类别 一般设定为最终所分类别的4—6倍

→定义最大循环次数(Maximum Iterations): 24

最大循环次数(Maximum Iterations)是指ISODATA重新聚类的最多次数,这是为了避免程序运行时间太长或由于没有达到聚类标准而导致的死循环。一般在应用中将循环次数都取6次以上。 →设置循环收敛阈值(Convergence Threshold):0.95

收敛阈值(Convergence Threshold)是指两次分类结果相比保持不变的像元所占最大百分之此值的设立可以避免ISODATA无限循环下去。

→点击OK按钮(关闭Unsupervised Classification对话框,执行非监督分类,获得一个初步的分类结果) 初步分类后,结果如下图:

第 3 页 共 3 页

遥感影像非监督分类实习报告

第三部 分类评价:

1,在视窗中同时显示 caijianhou_40.img(下)和 非监督分类_caijianhou 4001.img(上)。caijianhou_40.img 显示方式用红5、绿4蓝2。

2,打开分类图像属性并调整字段显示顺序。选择Raster菜单项,Rster---Attributes,

打开Raster Attribute Editor对话框(非监督分类_caijianhou 4001.img的属性表),并调整字段显示顺序: Raster Attribute Editor 对话框菜单条:Edit→ Column Properties →column properties对话框

第 4 页 共 4 页

遥感影像非监督分类实习报告

如下图所示:

在Columns中选择要调整显示顺序的字段,通过Up、 Down、Top、Bottom等几个按钮调整其合适的位置,通过选择Display Width调整其显示宽度,通过Alignment调整其对齐方式。如果选择Editable复选框,则可以在Title中修改各个字段的名字及其它内容。

3,给各个类别赋相应的颜色

Raster Attribute Edi tor对话框,点击一个类别的Row字段从而选

第 5 页 共 5 页

遥感影像非监督分类实习报告

择该类别,右键点击该类别的Color字段(颜色显示区),打开As Is菜单并选择一种颜色。重复以上步骤直到给所有类别赋予合适的颜色 。操作后结果如下图所示:

Viewer显示结果如下图:

第 6 页 共 6 页

遥感影像非监督分类实习报告

4,不透明度设置。由于分类图像覆盖在原图像上面,为了对单个类别的判别精度进行分析,首先要把其它所有类别的不透明程度(Opacity)值设为0 (即改为透明),而要分析的类别的透明度设为1(即不透明)。

5,确定类别专题意义及其准确程度 。

视窗菜单条:Utility→flicker →viewer Flicker 对话框→Auto Mode 本小步是设置分类图像在原图像修背景上闪烁, 观察它与背景图像之间的关系从而断定该类别的专题意义,并分析其分类准确与否。

第 7 页 共 7 页

遥感影像非监督分类实习报告

6,标注类别的名称和相应颜色。Raster Attribute Editor对话框(打开非监督分类_caijianhou 4001.img的属性表)

→点击刚才分析类别的ROW字段从而选择该类别 →点击该类别的class Names字段从而进入输入状态

→在该类别的Class Names字段中输入其专题意义(如居民区建筑),并按回车键

→右键点击该类别的Color字段(颜色显示区) →As Is菜单→选择一种合适的颜色

重复以上4、5、6三步直到对所有类别都进行了分析与处理。 非监督分类结果:

第 8 页 共 8 页

遥感影像非监督分类实习报告

7,分类后处理之聚类统计

ERDAS图标面板工具条:点击Interpreter图标->GIS Analysis Clump->Clump对话框(如下图)

第 9 页 共 9 页

遥感影像非监督分类实习报告

在Clump对话框中,需要确定下列参数: ① 确定输入文件(Input File):40_2010_uned.img

② 定义输出文件(Output File): 40_2010_uned_clump .img ③文件坐标类型((Coordinate Type) : Map

④处理范围确定(Subset Definition) : ULX / Y, LRX / Y (缺省状态为整个图像范围,可以应用Inquire Box定义子区)

⑤确定聚类统计邻域大小(Connect Neighbors): 8,(统计分析将对每个像元四周的8个相邻像元进行)

⑥OK(关闭Clump对话框,执行聚类统计分析) 在virwer中打开40_2010_uned_clump .img 如下图

第 10 页 共 10 页

遥感影像非监督分类实习报告

8,分类后处理之去除处理

ERDAS图标面板工具条:点击Interpreter图标->GIS Analysis->Eliminate->Eliminate对话框

(如下图所示)

第 11 页 共 11 页

遥感影像非监督分类实习报告

在Eliminate对话框中,需要确定下列参数:

① 确定输入文件(Input File): 40_2010_uned_clump .img ② 定义输出文件(Output File): 40_2010_uned_eliminate .img ③文件坐标类型(Coordinate Tvpe ): Map

④处理范围确定(Subset Definition): ULX / Y, LRX / Y,(缺省状态为整个图像范围,可以应用Inquire'Box定义子区)

⑤确定最小图斑大小(Minimum): 16 pixels ⑥确定输出数据类型(Output): Unsigned 8 Bit ⑦OK(关闭Eliminate对话框,执行去除分析)

在viewer中打开40_2010_uned_eliminate .img,如下图所示

8,分类后处理之分类重编码

ERDAS图标面板工具条:点击Interpreter图标->GIS Analysis-> Recode->Recode对话框。如下图所示

第 12 页 共 12 页

遥感影像非监督分类实习报告

① 确定输入文件(Input File):40_2010_uned_eliminate.img ② 定义输出文件(Output File): 40_2010_uned_recode.img ③ 设置新的分类编码((Setup Recode ):点击Setup Recode按钮 ④ 打开Thematic Recode表格

⑤ 根据需要改变New Value字段的取值(直接输入)

⑥ 在本例中将原来的35类依次两两合并,形成五类(图4.20) ⑦OK(关闭Thematic Recode表格,完成新编码输入) ⑧确定输出数据类型(Output): Unsigned 4 Bit

⑨OK(关闭Recode对话框,执行图像重编码。输出图像将按照New Value变换专题分类图像属性,产生新的专题分类图像)

第 13 页 共 13 页

遥感影像非监督分类实习报告

在视窗中打开重编码后的专题分类图像,查看起分类属性表:

改变颜色

第 14 页 共 14 页

遥感影像非监督分类实习报告

9,精度评估

点击Classifier→

弹出如下窗口:

点击

,打开重分类之后的图像。点击Edit→Create/Random Point

第 15 页 共 15 页

遥感影像非监督分类实习报告

设置参数。其中Number of Points 数值要大于256,勾选

,其中Minimum Points 要大于20或30,最好设置

为30,点击OK,在弹出窗口中一直点击是。 在弹出窗口中点击View→Select Viewer 之后点击一幅图像,重复操作一次,第二次点击另外一幅图像,使之两幅图像都显示待验证的点。之后点击View→Show All 显示待验证的点。继续点击Edit→Show Class Values 显示原来的类别号。显示如下窗口:

之后找出每一个待验证的点进行验证,在Reference中输入确定后的类别号,直至验证完所有的点。

第 16 页 共 16 页

遥感影像非监督分类实习报告

我在精度评估的时候由于验证一些点之后ERDAS软件总是崩溃,所以最后没有验证完成。点还没有校正完,软件就停在这里了。

第 17 页 共 17 页

遥感影像非监督分类实习报告

10,做缓冲区

ERDAS图标面板工具条:点击Interpreter图标->GIS Analysis->Search->Search对话框。如下图所示

第 18 页 共 18 页

遥感影像非监督分类实习报告

① 确定输入文件(Input File):40_2010_uned_eliminate.img ② 定义输出文件(Output File): 40_2010_uned_search.img ③ Distance to search :>=10 ④ OK

⑤ ERDAS图标面板工具条:点击Modeler图标->Modeler Maker打开对话框; 如下图所示

第 19 页 共 19 页

遥感影像非监督分类实习报告

⑥ 建立模型

第 20 页 共 20 页

遥感影像非监督分类实习报告

⑦ 点OK

44

第 21 页 共 21 页

遥感影像非监督分类实习报告

FileType选择Thematic ⑧ OK

⑨ 在视窗中打开生成的search2.img

⑩ 分类改变颜色(建立缓冲区前后对比,左为建立缓冲区后影像) 按相同方法对 LS5_TM_20081210_022812_022837_121040_FASTB_L2

第 22 页 共 22 页

遥感影像非监督分类实习报告

进行系列处理,用以进行变化检测。 11,变化检测

① ERDAS图标面板工具条:点击Modeler图标->Modeler Maker打开对话框; 如下图所示

② 建立模型

第 23 页 共 23 页

遥感影像非监督分类实习报告

③ OK

④ 在视窗中打开变化检测后图像

第 24 页 共 24 页

遥感影像非监督分类实习报告

⑤ 查看属性表

第 25 页 共 25 页

遥感影像非监督分类实习报告

⑥ 变化分析 1 2009年数据 2 3 4 5

1 2584 2705 623 732 655 2 1236 1402 1086 1182 854 2008年数据 3 1012 1158 906 1010 771 4 409 474 451 529 277 5 544 314 626 279 407

第 26 页 共 26 页

遥感影像非监督分类实习报告

⑥ 变化分析 1 2009年数据 2 3 4 5

1 2584 2705 623 732 655 2 1236 1402 1086 1182 854 2008年数据 3 1012 1158 906 1010 771 4 409 474 451 529 277 5 544 314 626 279 407

第 26 页 共 26 页

本文来源:https://www.bwwdw.com/article/gp4x.html

Top