遥感原理综合性实验报告(erdas)
- 格式:doc
- 大小:2.22 MB
- 文档页数:16
1 实验编号: 四川师范大学 实验报告 2014 年 8 月 30 日 地理与资源科学学院 2013级 1 班 实验名称:遥感数字处理图像实验 姓 名: 成 绩: 学 号: 指导教师: 程武学
1.实验目的和要求 帮助学生在了解RS基本组成与数据结构模型的基础上,重点学习使用ERDAS IMAGINE9.1软件进行视窗操作、数据数据预处理、图像解译、图像分类和矢量功能;掌握遥感图像前后处理和解译、分类地理的技术流程和方法,数据库建设以及地理数据的编辑和管理;配合具体实例运用GIS空间分析工具。通过系列实习过程,重点培养学习者掌握RS提取信息的基本过程和技巧,并可初步用来解决运用遥感提取信息的问题。 2.主要仪器与试剂 主要仪器:局域网系统;高级PC 计算机;Windows 2000,ERDAS IMAGINE8.5 3.实验步骤: 一、图像几何校正过程:
(1)数据预处理:ERDAS图标面板菜单:Main->Data Preparation->Image Geometric
Correction,打开Set Geo-Correction Input File对话框(如右图1) (2)显示图像文件:ERDAS图标面板菜单条:Session->Tile Viewers;然后,在Viewer #1中打开需要校正的Landsat TM图像:01.img;在Viewer #2中打开作为地理参考的校正过的SPOT图像:00.img(如图2)
图1 Set Geo-Connection Input File对话框 mask.img话框 2 图5 Geo Correction Tools对话框
图2打开两幅图像后的窗口 (3)启动几何校正模块:在Viewer #1菜单条:Raster->Geometric Correction,打开Set
Geometric Model对话框(图3)选择多项式几何校正计算模型:Polynomial – OK,同时打开Polynomial Model Properties对话框(图4)和Geo Correction Tools对话框(图5)
在Polynomial Model Properties对话框中,定义多项式模型参数及投影参数:定义多项式次方(Polynomial ); 再定义投影参数。 (4)Apply->Close,打开GCP Tool Reference Setup对话框(图5)
图3 Set Geometric Model对话框 图4 Polynomial ModelProperd对话框 3
(5)启动控制点工具:首先在GCP Toots Reference Setup对话框(图5)中选择采点模式:选择视窗采点模式:Existing Viewer,单击OK(关闭GCP Tools Reference Setup对话框),自动打开Viewer Selection Instiuctions指示器(图6) (6)在显示作为地理参考图像panAtlanta.img的Viewer #2中点击左键,打开Reference Map Information提示框(图7)(显示参考图像的投影信息),OK(关闭Reference Map Information提示框)
(7)采集地面控制点:整个屏幕将自动变化为如图8所示的状态:其中包含两个主视窗、两个放大窗口、两个关联方框(分别位于两个视窗中,指示放大窗口与主视窗的关系)、控制点工具对话框、几何校正工具等。
图6 Viewer Selection Instiuctions指示器 图7 Reference Map Information提示框 4
图8 Reference Map Information 窗口 (8)采集地面检查点:以上所采集的GCP的类型为Control Point(控制点),用于控制计算,建立转换模型及多项式方程。下面所要采集的GCP的类型均是Check(检查点),用于检验所建立的转换方程的精度和实用性。依然在GCP Tool对话框状态下: (a)在GCP Tool菜单条中确定GCP类型:Edit->Set Point Type->Check (b)在GCP Tool菜单条中确定GCP匹配参数(Matching Parameter):Edit-> Point Matching->打开GCP Matching对话框。 在GCP Matching对话框中,需要定义下列参数:①匹配参数(Matching Parameters);②最大搜索半径(Max. Search Radius)为3;③搜索窗口大小(Search Window Size): X:5 Y:5;④约束参数(Threshold Parameters ):设相关阈值(Correlation Threshold):0.8;⑤删除不匹配的点(Discard Unmatched Point):Active (c)在匹配所有/选择点(Match All / Selected Point)选项组中设从输入到参考(Reference from Input)或从参考到输入(Input from Reference) (d)Close(关闭GCP Matching对话框) 5
(e)确定地面检查点:在GCP Tool工具条中选择Create GCP图标,并将Lock图标打开,锁住Create GCP功能,如同选择控制点一样,分别在Viewer #l和Viewer #2中定义5个检查点,定义完毕后点击Unlock图标,解除Create GCP功能。 (f)计算检查点误差:在GCP Tool工具条中点击Compute Error图标,检查点的误差就会显示在GCP Tool的上方,只有所有检查点的误差均小于一个像元(Pixel),才能继续进行合理的重采样。一般来说,如果你的控制点(GCP)定位选择比较准确的话,检查点匹配会比较好,误差会在限差范围内。否则,若控制点定义不精确,检查点就无法匹配,误差会超标。 (9)检验校正结果:(a)打开两个平铺视窗(open and Tile two Viewer),视窗菜单条:File->Open->RasterOption->图像文件,ERDAS图标面板:Session->Tile Viewer->平铺视窗 (b)建立视窗地理连接关系(Geo Link two Viewer),在Viewer #1中:按住右键->快捷菜单->Geo Link /Unlink,在Viewer #2中:点击左键一建立与Viewer #1的连接 (c)通过查询光标进行检验(Check with Inquire Cursor),在Viewer #1中:按住右键->快捷菜单->Inquire Cursor->打开光标查询对话框,在Viewer #1中:移动查询光标,观测其在两屏幕中的位置及匹配程度,并注意光标查询对话框中数据的变化。:
图10 校正前TM影像 6 图11 校正后TM影像
二、图像去霾过程: (1)ERDAS 图标面板菜单条.: 单击Main->Image Interpreter->Radiometric Enhancement
打开 Haze Reduction 对话框(图12),在对话框中,需要设置下列参数: ①确定输入文件(Input File):klon_tm.img; ②定义输出文件(Ouput File):haze.img; ③定义坐标类型为(Coordinate Type):Map; ④处理范围确定(Subset Definition):在ULX/Y、LRX/Y微调框中输入需要的数值(缺省状态为整个图像范围,可以应用InquireBox定义子区); ⑤处理方法(Method):选择Landsat 5 TM (2)单击OK按钮,关闭 Haze Reduction对话框,执行去霾处理。去霾处理的前后对比效果如下图13和14所示:
图12 Haze Reduction 对话框 7
图13 去霾处理前影像 图14 去霾处理后影像 三、非监督分类过程:
(1)调出非监督分类对话框:
在ERDAS图标面板工具条中点击DataPrep图标Data Preparation菜单->Unsupervised 8
Classification菜单项Unsupervised Classification对话框(图17);或在ERDAS图标面板工具条中点击Classifier图标,Classification菜单,Unsupervised Classification菜单项,Unsupervised Classification对话框(图18)。
图19设置参数后的Unsupervised Classification对话框 (2)进行非监督分类
图17方法一打开的Unsupervised Classification对话框 图18方法二打开的Unsupervised
Classification对话框 9
Unsupervised Classification对话框(图19): ①确定输出文件: 11.img(要被分类的图像)②确定输出文件: 12.img〔即将产生的分类图像) ③选择生成分类模板文件: Output Signature Set)④确定分类模板文件(FileName): 13.sig ⑤对Clustering Options选择Initialize from Statistic单选框 Initialize from Statistics指由图像文件整体(或其AOl区域)的统计值产生自田聚类,分出类别的多少由用户自己决定。Use Signature Means是基于选定的模板文件进行非监督分类,类别的数目由模板文件决定。 ⑥确定初始分类数(Number of classes): 10(分出10个类别)。实际工作中一般将分类数取为最终分类数的2倍以上。 ⑦点击Initializing Options按钮可以调出File Statistics Options对话框以设置1SODATA的一些统计参数。 ⑧点击Color Scheme Options按钮可以调出‘Output Color Scheme Options对话框以决定输出的分类图像是彩色的还是黑白的。这两个设置项使用缺省值。 ⑨定义最大循环次数(Maximum Iterations): 24;最大循环次数(Maximum Iterations)是指ISODATA重新聚类的最多次数,这是为了避免程序运行时间太长或由于没有达到聚类标准而导致的死循环。一般在应用中将循环次数都取6次以上。 ⑩设置循环收敛阅值(Convergence Threshold): 0.95 收敛阂值(Convergence Threshold)是指两次分类结果相比保持不变的像元所占最大百分比。此值的设立可以避免ISODATA无限循环下去。 (3)点击OK按钮(关闭Unsupervised Classification对话框, 执行非监督分类,获得一个初步的分类结果,如图20)。