遥感概论实验指导书

1、遥感概论实验指导书实验一 遥感图像的光学合成原理一、目的要求1、了解彩色的基本特性和相互关系；2、掌握三原色及其补色，掌握加色法；3、认识正负像片的产生过程。二、实验地点：机房三、实验步骤：本实习在CAI软件中进行。包括以下五部分：1、彩色基本特性及其相互关系（1） 从CAI软件主界面中进入“遥感光学合成原理”子目录；（2） 进入“彩色与非彩色”，对比彩色图像与非彩色图像对地物特征的表现；（3） 退回到“遥感光学合成原理”进入“彩色的特性”分别再进入“明度”、“色调”和“饱和度界面”，观察枫叶色彩的变化，分析对比彩色的三大特性在色彩中的影响；（4） 退回到“遥感光学合成原理”进入“颜色立体”，

2、观察颜色立体中彩色三大特性的表示方法，掌握明度、色调和饱和度的相互关系。2、三原色、补色和加色法（1）从 CAI 软件主界面进入“遥感光学合成原理”子目录；（2）进入“三原色”子目录；（3）当蓝、绿、红均为“0”时，屏幕中圆盘的颜色是 色。（4）分别拖动蓝、绿和红的标尺，使绿和红为0，将蓝色标尺从0向255拖动，观察屏幕圆盘颜色的变化情况为 。以同样的方式，分别观察绿和红随着标尺从0至255变化时，屏幕圆盘颜色的变化。（5） 将蓝、绿和红均设为255时，屏幕圆盘为 色。 （6） 使红为0，拖动标尺使蓝和绿均为255，屏幕圆盘为 色？同样，使绿为0，拖动红和蓝的标尺至255，屏幕中圆盘为 色？使

3、蓝为0，拖动红和绿的标尺至255，屏幕中圆盘为 色。这种颜色合成的方法为 。（7） 如何改变蓝、绿和红色，使屏幕圆盘为灰色？答：3、补色（1）从 CAI 软件主界面进入“遥感光学合成原理”子目录；（2）进入“互补色”子目录；（3）单击互补色按钮，圆盘中出现一对互补色，再单击融合按钮，观察圆盘中混合后的颜色为白色；重复上述过程，观察圆盘中出现的互补色，以及混合后圆盘中的颜色；（4） 击非互补色按钮，圆盘中出现一对非互补色，单击融合按钮，观察圆盘中混合后的颜色为各种彩色；重复上述操作，观察圆盘中出现的非互补色以及混合效果。5、彩色正负像片产生过程（1）从 CAI 软件主界面进入“遥感光学合成原理”

4、子目录；（2）分别进入“彩色负片的生成”、“彩色正片的生成”和“彩红外像片的生成”子目录；（3）观看系统演示过程，回答下面的问题。1）地物在彩色负片上的颜色与本身的颜色有何不同？2）彩色正片上地物的颜色与实际地物颜色相同吗？ 3）彩红外像片最适于记录哪类物体？实验二 ERDAS IMAGINE 8.4软件演示与介绍数据输入输出一、实验目的1.了解ERDAS IMAGINE 8.4 的启动与关闭方法 2.熟悉ERDAS IMAGINE 8.4界面 3.常用数据的输入与输出 二、实验数据p129r033_7t19990812_z48_nn10p129r033_7t19990812_z48_nn20

5、p129r033_7t19990812_z48_nn30p129r033_7t19990812_z48_nn40p129r033_7t19990812_z48_nn50p129r033_7t19990812_z48_nn61p129r033_7t19990812_z48_nn70三、实验步骤1. ERDAS IMAGINE 8.4的启动方法： 方法一：单击Windows任务栏的“开始”按钮，然后指向“程序”文件夹，再单击ERDAS IMAGINE 8.4程序的图标即可； 方法二：如果桌面上创建了ERDAS IMAGINE 8.4的快捷方式图标，双击快捷图标即可启动ERDAS IMAGINE 8

6、.4 2. ERDAS IMAGE 8.4的关闭方法： 方法一：单击ERDAS IMAGINE 8.4窗口右上角控制按钮的“×”符号 方法二：单击Session(综合菜单)中的Exit IMAGINE命令 方法三：单击ERDAS IMAGINE 8.4窗口左上角的图标，或右键单击控制按钮以左的区域打开控制菜单，选择“关闭”命令3. ERDAS IMAGINE 8.4的界面介绍：启动ERDAS IMAGINE 8.4以后，用户首先看到的是ERDAS IMAGINE 8.4的图标面板（ICON PANEL）,包括菜单条（Menu Bar）和工具条（Tool Bar）两部分，其中提供了启动

7、ERDAS IMAGINE 8.4软件模块的全部菜单和图标。菜单条共有五项：Session（综合菜单）；Main（主菜单）；Tools（工具菜单）；Utilities（实用菜单）；Help（帮助菜单），图标面板工具条共有十一项：ERDAS IMAGINE Credits；Viewer；Import；Dataprep；Composer；Interpreter；Catalog；Classifier；Modeler；Vector；Radar;Virtual GIS;OrthoBASEERDAS IMAGINE 8.4图标面板工具条 图标 命令 功能 ERDAS IMAGINE Credits IMA

8、GINE Credits 查阅ERDAS IMAGINE 8.4信用卡 Viewer Start IMAGINE Viewer 打开ERDAS IMAGINE 8.4视窗 Import Import/Export 启动数据输入输出模块 Dataprep Data Preparation 启动数据预处理模块 Composer Map Composer 启动专题制图模块 Interpreter Image Interpreter 启动图像解译模块 Catalog Image Catalog 启动图像库管理模块 Classifier Image Classification 启动图像分类模块 Mod

9、eler Spatial Modeler 启动空间建模工具 Vector Vector 启动矢量功能模块 Radar Radar 启动雷达图像处理模块 4. 数据输入输出 4.1 了解常用输入输出数据格式 ERDAS IMAGINE 8.4的数据输入输出功能（Import/Export），允许输入多种格式的数据供ERDAS IMAGINE 8.4使用，同时允许将ERDAS IMAGINE 8.4的文件转换成多种数据格式。 ERDAS IMAGINE 8.4允许输入输出的数据格式可以通过以下方法看到： ERDAS IMAGINE 8.4图标面板菜单条：MainImport/Export （或ER

10、DAS IMAGINE 8.4图标面板工具条：单击“Import”图标） Import/Export(数据输入输出)对话框在Import/Export(数据输入输出)对话框中，选定Import (输入数据) 时，单击“Type”右侧的下拉菜单，拖动垂直滚动条即可浏览ERDAS IMAGINE 8.4允许输入的所有数据格式；选定Export (输出数据) 时，单击“Type”右侧的下拉菜单即可浏览ERDAS IMAGINE 8.4允许输出的所有数据格式。4.2 TIFF图像数据输入输出 （以e:实验二etm下的文件为例） 在ERDAS IMAGINE 8.4中打开TIFF图像文件时，将文件类型指

11、定为TIFF格式就可以直接在视窗中显示TIFF图像，具体过程如下： ERDAS IMAGINE 8.4图标面板菜单条：MainStart IMAGINE Viewer（或ERDAS IMAGINE 8.4图标面板工具条：单击“Viewer”图标） 打开一个二维视窗Viewer #1 单击视窗工具条图标 打开Select Layer To Add对话框在File of type右侧的下拉菜单中选择:TIFF(*.tif) 在Look in右侧的下拉菜单中选择所要显示的文件的路径: e:imagecbers011056 在Look in下面方框中选择所要显示文件名: p129r033_7t1999

12、0812_z48_nn10.tif 单击OK按钮（退出Select Layer To Add对话框），完成文件的显示。 如果要在图像解译器（Interpreter）或其他模块下对图像作进一步处理操作，依然需要将TIFF文件转换为IMG文件 具体过程如下： 在二维视窗中打开TIFF格式文件p129r033_7t19990812_z48_nn10.tif视窗菜单条：FileSaveTop Layer As 打开Save L20000100215.tif As对话框单击Save L20000100215.tif as正下方方框右侧的文件夹 打开Save as对话框,切换到自己的文件夹。在File o

13、f type右侧的下拉菜单中选择:IMAGINE Image(*.img) 在File name右侧的方框中输入文件名：如tm1.img 双击OK按钮 返回Save L20000100215.tif As对话框 单击OK（退出Save 20000520477_2.tif As对话框） 打开Copying Image进度条 单击OK按钮完成文件的转换按照上述步骤完成其他6个波段的文件转换。4.3. 组合多波段数据 （应用前面生成的TM7个波段的图像） 在实际工作中，对遥感图像的处理和分析都是针对多波段图像进行的，所以还需要将若干波段图像文件组合成一个多波段图像文件，具体过程如下： ERDAS I

14、MAGINE 8.4图标面板菜单条：MainImage Interpreter(或ERDAS IMAGINE 8.4图标面板工具条：单击“Interpreter”图标) 打开Image Interpreter对话框 选择Utilities 打开Utilities对话框选择Layer Stack 打开Layer Selection and Stacking对话框在Layer Selection and Stacking对话框中依次Add（加载）单波段图像文件： Input File:(*.img)（输入单波段文件）:tm1.imgAdd Input File:(*.img)（输入单波段文件）:

15、tm 2.imgAdd Input File:(*.img)（输入单波段文件）: tm 3.imgAdd .Input File:(*.img)（输入单波段文件）: tm 7.imgAdd Output File:(*.img)（输出多波段文件）:layerstack.img Output (输出数据类型):Unsigned 8 Bit Output Options (波段组合选择):Union 选定Ignore Zero In Stats(输出统计忽略零值)单击OK(关闭Layer Selection and Stacking对话框,执行波段组合) 打开Modeler-running mod

16、el:layermerge.pmdl进度条 单击Modeler-running model:layermerge.pmdl进度条中的OK按钮完成波段组合实验三 图像的拼接与裁剪一、实验目的1.图像拼接处理（Mosaic Imgaes）2.图像分幅裁剪（Subset Image） 二、实验数据wasia1_mss.img、wasia2_mss.img、wasia3_mss.imglanier.img、layerstack.img、yinchuanshi.shp、yinchuanshi.prj、yinchuanshi.dbf、yinchuanshi.shx三、实验步骤1.图像拼接处理 （以c:Pr

17、ogram File IMAGINE 8.4examples下的wasia1_mss.img、wasia2_mss.img、wasia3_mss.img） 图像拼接处理是将具有地理参考的若干相邻图像合并成一幅图像或一组图像，需要拼接的输入图像必须含有地图投影信息，或者说输入图像必须经过几何校正处理或进行过校正标定。 具体操作过程如下： ERDAS IMAGINE 8.4 图标面板菜单条：MainData Preparation(或单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标) 打开Data Preparation对话框Mosaic Tool菜单 打开Mosai

18、c Tool对话框Mosaic Tool视窗菜单条：EditAdd Image 打开Add Image for Mosaic对话框在Add Image for Mosaic对话框中，需要设置下列参数： Image Filename（拼接图像文件）：wasia1_mss.img Image Area Options(图像拼接区域)：Compute Active Area单击Add（图像wasia1_mss.img被加载到Mosaic视窗中） 按照同样的过程加载wasia2_mss.img和wasia3_mss.img 单击Close（关闭Add Image for Mosaic对话框） Mosa

19、ic Tool视窗菜单条：EditImage Matching 打开Matching Options对话框 在Matching Options对话框中作如下设置： 在Matching Method中选择Overlap Area 单击OK按钮（关闭Matching Options对话框） Mosaic Tool视窗菜单条：EditSet Overlap Function 打开Set Overlap Function对话框 在Set Overlap Function对话框，设置以下参数： 在Intersection Type中设置No Cutline Exits（没有裁切线） 在Select Fu

20、nction中设置Average（均值）单击Apply 单击Close（关闭Set Overlap Function对话框） Mosaic Tool视窗菜单条：ProcessRun Mosaic 打开Run Mosaic对话框 在Run Mosaic对话框中设置下列参数： Output File Name（确定输出文件名）：wasia_mosaic.img Witch Outputs（确定输出图形区域）：All 在Ignore Input Value中设置0 在Output Background Value中设置0 选定Stats Ignore Value，并设置为0单击OK（关闭Run Mo

21、saic对话框，运行图像拼接） Mosaic Tool视窗菜单条：FileClose 打开Verify Save on Close对话框单击“否”按钮(关闭Mosaic Tool视窗,退出Mosaic工具)在viewer窗口中分别打开wasia1_mss.img、wasia2_mss.img、wasia3_mss.img以及拼接后的图像wasia_mosaic.img，体会图像拼接的意义。2.图像分幅裁剪 在实际工作中，经常需要根据研究工作范围对图像进行分幅裁剪，按照ERDAS IMAGINE 8.4实现图像分幅裁剪的过程，可以将图像分幅裁剪为两类型：规则分幅裁剪，不规则分幅裁剪。 2.1规则

22、分幅裁剪 (以c:Program File IMAGINE 8.4exampleslanier.img为例) 规则分幅裁剪是指裁剪图像的范围是一个矩形，通过左上角和右下角两点的坐标可以确定图像的裁剪位置，过程如下： 方法一: ERDAS IMAGINE 8.4 图标面板菜单条：MainData Preparation(或单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标) 打开Data Preparation 对话框单击Subset Image按钮，打开Subset对话框在Subset对话框中需要设置下列参数： 输入文件名（Input File）:lanier.im

23、g 输出文件名（Output File）:lanier_subset1.img 坐标类型（Coordinate Type）:Map 裁剪范围（Subset Definition）:ULX、ULY、LRX、LRY （注：ULX，ULY是指左上角的坐标，LRX,LRY是指右下角的坐标，缺省状态为整个图像范围） 本次裁剪坐标已确定左上角ULX：236355.00 ULY：3800230.00 右下角LRX：244455.00 LRY：3793420.00将该坐标值对应输入ULX、ULY、LRX、LRY中输出数据类型（Output Data Type）:Unsigned 8 Bit 输出文件类型（Ou

24、tput Layer Type）:Continuous 输出统计忽略零值:Ignore Zero In Output Stats 输出像元波段（Select Layers）:2,3,4OK(关闭Subset对话框，执行图像裁剪)在viewer窗口中打开lanier_subset1.img，并与lanier.img对比。方法二: ERDAS IMAGINE 8.4图标面板菜单条:MainStart IMAGINE Viewer(或单击RDAS IMAGINE 8.4图标面板工具条“Viewer”图标) 打开一个二维视窗 单击视窗工具条图标 打开Select Layer To Add对话框 在Se

25、lect Layer To Add对话框完成以下设置： Look In:examples File Name:lanier.img Files of type:IMAGINE Image 双击OK按钮 在二维视窗中打开lanier.img文件单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标 打开Data Preparation对话框单击Subset Image按钮 打开Subset对话框在Subset对话框中需要设置下列参数： 输入文件名（Input File）:lanier.img 输出文件名（Output File）:lanier_subset2.img 坐标

26、类型（Coordinate Type）:Map 输出数据类型（Output Data Type）:Unsigned 8 Bit 输出文件类型（Output Layer Type）:Continuous 输出统计忽略零值:Ignore Zero In Output Stats 输出像元波段（Select Layers）:2,3,4单击From Inquire Box按钮 打开Invalid Coordinate Type对话框 单击Continue 在显示图像文件lanier.img视窗中单击工具条的按钮,打开Inquire Cursor对话框,在视窗中移动十字光标,确定裁剪范围左上角和右下角(

27、自己确定左上角和右下角的范围),读取其坐标分别填入Subset Image对话框的ULX,ULY中和LRX,LRY中 单击OK按钮(关闭Subset对话框，执行图像裁剪)在view中分别打开lanier_subset2.img和原始图像lanier.img，并进行对比。方法三: 首先在视窗中打开lanier.img文件AOITools打开AOI工具面板单击矩形框确定裁剪范围FileSaveAOI Layer As 打开Save AOI As对话框,输入文件名:lanier.aoi 单击OK(退出Save AOI As对话框)单击ERDAS IMAGINE 8.4 图标面板工具条“DataPre

28、p”图标 打开Data Preparation对话框 单击Subset Image按钮 打开Subset对话框 在Subset对话框中需要设置下列参数： 输入文件名（Input File）:lanier.img 输出文件名（Output File）:lanier_subset3.img 坐标类型（Coordinate Type）:Map 输出数据类型（Output Data Type）:Unsigned 8 Bit 输出文件类型（Output Layer Type）:Continuous 输出统计忽略零值:Ignore Zero In Output Stats 输出像元波段（Select La

29、yers）:2,3,4 单击AOI按钮 打开Choose AOI对话框在Choose AOI对话框作如下设置： AOI Source:File AOI File: lanier.aoi单击OK(退出Choose AOI对话框) 单击OK(退出Subset对话框,执行图像裁剪) 单击OK(退出Modeler对话框,完成图像裁剪)在viewer中打开lanier_subset3.img和lanier.img，并进行对比。1.2不规则分幅裁剪 不规则分幅裁剪是指裁剪图像的边界范围是个任意多边形，无法通过左上角和右下角两点的坐标确定图像的裁剪位置，而必须事先生成一个完整的闭合多边形区域，可以是一个AO

30、I多边形，也可以是ArcInfo的一个Polygon Coverage，针对不同的情况采用不同的裁剪过程。 （一） AOI多边形裁剪 (以c:Program File IMAGINE 8.4exampleslanier.img为例) 在视窗中打开需要裁剪的图像应用AOI工具绘制多边形FileSaveAOI Layer As 打开Save AOI As对话框,输入文件名:lanier2.aoi ERDAS IMAGINE 8.4 图标面板菜单条：MainData Preparation(或单击ERDAS IMAGINE 8.4 图标面板工具条“DataPrep”图标)，打开Data Prepar

31、ation 对话框 单击Subset Image按钮，打开Subset对话框在Subset对话框中需要设置下列参数： 输入文件名（Input File）:lanier.img 输出文件名（Output File）:lanier_subset4.img Coordinate Type:Map Output:Unsigned 8 bit Oupput:Continuous 选定Ignore Zero in Output Stats Select Layers:2,3,4单击AOI按钮，打开Choose AOI对话框选定File 选择文件：lanier2.aoi单击OK(退出Choose AOI对话

32、框) 单击OK（退出Subset对话框,执行图像裁剪） 单击Modeler对话框中的OK按钮完成图像裁剪.在viewer中打开lanier_subset4.img和lanier.img，并进行对比。（二） ArcInfo多边形裁剪 (以layerstack.img为例) 如果是按照行政区划边界或自然区划边界进行图像的分幅裁剪，往往是首先利用Arcinfo或者ERDAS的Vector模块绘制精确的边界多边形（Polygon），然后以Arcinfo的Polygon为边界进行图像裁剪。对于这种情况，需要调用ERDAS其它模块的功能分两步完成。 第一步：将Arcinfo多边形转换为栅格图像文件 ERD

33、AS IMAGINE 8.4 图标面板菜单条：MainImage Interpreter (或单击ERDAS IMAGINE 8.4 图标面板工具条“Interpreter”图标) 打开Image Interpreter对话框 单击Utilities打开Utilities对话框选择Vector To Raster 打开Vector To Raster对话框在Vector to Raster对话框中设置下列参数： Input Vector File（输入矢量文件名称）：yinchuanshi.shp Vector Type(确定矢量文件类型)：polygon Use Attribute as V

34、alue(使用矢量属性值)：idOutput File（输出栅格文件名称）：raster.img Data Type(栅格数据类型)：Unsigned 8 bit Layer Type(栅格文件类型)：Thematic Size Difinition(转换范围大小)：ULX,ULY,LRX，LRY（缺省条件下为整个图像范围） Units（坐标单位）：Meters Cell Size（输出像元大小）：X:30/Y:30 选择Square Cell（正方形像元）单击OK（关闭Vector To Raster对话框，执行矢量到栅格的转换）矢量图栅格图第二步：通过掩模运算（Mask）实现图像的不规则裁

35、剪 ERDAS IMAGINE 8.4 图标面板菜单条：MainImage Interpreter (或单击ERDAS IMAGINE 8.4 图标面板工具条“Interpreter”图标) Image Interpreter Utilities Mask 打开Mask对话框在Mask对话框中设置下列参数： Input File(输入需要裁剪的图像文件名称):layerstack.img Input Mask File(输入掩模文件名称)：raster.img 点击Setup Recode按钮设置裁剪区域内新值（New Value）为1，区域外取0值 选择Intersection（确定掩模区域

36、作交集运算） Output File（输出图像文件名称）：mask.img Output (输出数据类型)：Unsigned 8 bit单击OK（关闭Mask对话框，执行掩模运算）打开mask.img，如下图。实验四 遥感图像几何校正一、 实验目的掌握在Erdas imagine环境中进行图像几何校正的方法。几何校正就是将图像数据投影到平面上，使其符合地图投影系统的过程；而将地图坐标系统赋予图像数据的过程，称为地理参考（Georeferencing）。由于所有地图投影系统都遵从于一定的地图坐标系统，所以几何校正过程包含了地理参考过程。 二、 实验数据panAtlanta.img、tmAtlan

37、ta.img三、 实验步骤视窗采点模式： 如果已经拥有需要校正图像区域的数字图像、或经过校正的图像、或注记图层，就可以应用视窗采点模式，直接以数字地图、或经过校正的图像、或注记图层作为地理参考，在视窗中打开该图层，从中采集控制点。 （以c:Program File IMAGINE 8.4examples下的panAtlanta.img和tmAtlanta.img为例，前者是作为地理参考的图像，后者为需要校正的图像） 第一步：显示图像文件 在ERDAS IMAGINE 8.4 图标面板工具条点击“Viewer”图标两次，打开两个视窗 ERDAS IMAGINE 8.4 图标面板菜单条：Sessi

38、onTile Viewer 然后，在Viewer#1中打开需要校正的Landsat TM图像：tmAtlanta.img 在Viewer#2中打开作为地理参考的校正过的SPOT图像：panAtlanta.img第二步：启动几何校正模块 Viewer#1菜单条：RasterGeometric Correction 打开Set Geometric Model对话框 选择多项式几何校正计算模型：Polynormial单击OK按钮 同时打开Geo Correction Tools对话框和Polynormial Model Properties对话框在Polynormial Model Properti

39、es对话框中，定义多项式模型参数及投影参数 定义多项式次方(Ploynomial)：2，其余取缺省值单击Apply按钮 单击Close按钮 打开GCP Tool Reference Setup对话框说明：该实例是采用视窗采点模式，作为地理参考的SPOT图像已经含有投影信息，所以这里不需要定义投影参数。如果不是采用视窗采点模式，或者参考图像没有包含投影信息，则必须在这里定义投影信息，包括投影类型及其对应的投影参数。 第三步：启动控制点工具 在GCP Tool Reference Setup对话框中 选择视窗采点模式：Existing Viewer 单击OK按钮（关闭GCP Tool Refere

40、nce Setup对话框）打开Viewer Selection Instruction指示器 在显示作为地理参考图像panAtlanta.img的Viewer#2中点击左键 打开Reference Map Information提示框单击OK按钮（关闭Reference Map Information提示框） 此时整个屏幕将显示：两个主视窗、两个放大窗口、两个关联方框（分别位于两个视窗中，指示放大窗口与主视窗的关系）、控制点工具窗口、几何校正工具。表明控制点工具启动，进入控制点采集状态 第四步：采集地面控制点 控制点工具窗口简介： 控制点工具窗口（GCP Tool）有菜单条、工具条、控制点数据

41、表和状态条四个部分组成。GCP窗口数据字段及其含义 字段 含义 Point# GCP顺序号,系统自动产生 Point ID GCP标识码,用户可以定义 > GCP当前选择状态提示符号 Color 输入GCP显示颜色 X Input 输入GCP的X坐标 Y Input 输入GCP的Y坐标 Color 参考GCP显示颜色 X Reference 参考GCP的X坐标 Y Reference 参考GCP的Y坐标 Type GCP的类型(控制点检查点) X Residual 单个GCP的X残差 Y Residual 单个GCP的Y残差 RMS Error 单个GCP的RMS误差 关于GCP Too

42、l窗口的说明: 1.输入控制点(Input GCP)是在原始文件视窗中采集的，具有原文件的坐标系统；而参考控制点(Reference GCP)是在参考文件视窗中采集的，具有已知的参考坐标系统， GCP工具将根据对应点的坐标值自动生成转换模型。 2.在GCP数据中，残差(Residuals)、中误差(RMS)、贡献率(Contribution)及匹配程度(Match)等参数，是在编辑GCP的过程中自动计算更新的，用户是不能任意改变的，但可以通过精确GCP位置来调整 3.每个IMG文件都可以有一个GCP数据集与之相关联，GCP数据集保存在一个栅格层数据文件中；如果IMG文件有一个GCP数据集存在的

43、话，只要打开GCP Tool，GCP点都会出现在视窗中 4.所有的输入控制点(Input GCP)都可以保存在图像文件中(Save Input)，也可以保存在控制点文件中(Save Input As)。如果保存在图像文件中，调用的方法如3所述，如果是保存在控制点文件中，可以通过加载调用(Load Input)。 5.参考控制点(Reference Input)也可以类似地保存在参考图像中(Save Reference)或控制点文件中(Save Reference As)，可以通过加载调用(Load Reference) 采集地面控制点具体过程如下： 在GCP Tool窗口中点击Select G

44、CP图标，进入GCP选择状态 在GCP数据表中将输入GCP的颜色（Color）设置为比较明显的黄色 在Viewer#1中移动关联方框位置，寻找明显的地物特征点，作为输入GCP 在GCP Tool窗口中点击Create GCP图标，并在Viewer#3中点击左键定点，GCP数据表将记录一个输入GCP，包括其编码、标识码。X坐标、Y坐标 在GCP Tool窗口中点击Select GCP图标，重新进入GCP选择状态 在GCP数据表中将参考GCP的颜色（Color）设置为比较明显的红色 在Viewer#2中移动关联方框位置，寻找对应的地物特征点，作为参考GCP 在GCP Tool窗口中点击Create

45、 GCP图标，并在Viewer#4中点击左键定点，系统将自动把参考点的坐标（X Reference，Y Reference）显示在GCP数据表中 在GCP Tool窗口中点击Select GCP图标，重新进入GCP选择状态，并将光标移回到Viewer#1，准备采集另一个输入控制点 不断重复以上步骤，采集若干GCP，直到满足所选定的几何校正模型为止（多项式变换的次方数为二时，最少需要六个控制点），而后每采集一个Input GCP，系统就自动产生一个Reference GCP。 第五步：采集地面检查点 以上所采集的GCP的类型(Type)均为Control(控制点),用于控制计算,建立转换模型及多

46、项式方程.下面所要采集的GCP的类型均是Check(检查点),用于检验所建立的转换方程的精度和适用性。 在GCP Tool窗口中 在GCP Tool窗口菜单条中确定GCP类型：EditSet Point TypeCheck 在GCP Tool窗口菜单条中确定GCP匹配参数（Matching Parameter）：EditPoint Matching打开GCP Matching对话框 在GCP Matching对话框，需要定义下列参数： Matching Parameter（匹配参数）: Max.Search Radius（最大搜索半径）：3Search Window （搜索窗口大小）：X:5；

47、Y:5 Threshold Parameters（约束参数）： Correction Threshold（相关阈值）:0.8 选定 Discard Unmatched Point（删除不匹配的点）单击Close（关闭GCP Matching对话框） 确定地面检查点：在GCP Tool工具条中选择Create GCP图标，并将Lock图标打开锁住Create GCP功能，象选择控制点一样，分别在Viewer#1和 Viewer#2中定义检查点，定义完毕后点击Unlock图标，解除Create GCP功能 计算检查点误差：在GCP Tool工具条中点击“Compute Error”图标，检查点的误

48、差值会显示在GCP Tool的上方，只有所有检查点的误差均小于一个像元，才能继续进行合理的重采样 第六步：图像重采样 重采样过程就是依据未校正图像像元值计算生成一幅校正图像的过程，原图像中所有栅格数据层都将进行重采样。ERDAS IMAGINE 8.4提供三种最常用的重采样方法： Nearest Neighbor：邻近点插值法，将最邻近像元值直接赋予输出像元 Bilinear Interpolation:双线性插值法，用双线性方程和2×2窗口计算输出像元值Cubic Convolution:立方卷积插值法,用立方方程和4×4窗口计算输出像元值 具体操作如下: 在Geo Co

49、rrection Tools对话框中选择“Image Resample”图标 打开Image Resample对话框 在Resample对话框中,定义重采样参数: Output File(输出图像文件):rectify.img Resample Method(选择重采样方法)：Nearest Neighbor Output Cell Size(定义输出图像范围)：ULX、ULY；LRX、LRY 定义输出像元大小为30。选定Ignore Zero in Stats Recalculate Output Default(设置重新计算缺省值)：Skip Factor:X:1;Y:1 单击Ok（关闭R

50、esample对话框） 第七步：保存几何校正模型 在Geo Correction Tools对话框中点击Exit按钮，退出图像几何校正过程。按照系统提示将图像几何校正模型保存为模式文件resample.gms. 第八步：检验校正结果 检验校正结果的基本方法是：同时在两个视窗中打开两幅图像，其中一幅是校正后的图像，一幅是当时的参考图像，通过视窗地理连接（Geo LinkUnlink）功能及查询光标功能进行目视定性检验，具体过程如下： 在ERDAS IMAGINE 8.4 图标面板工具条点击“Viewer”图标两次，打开两个视窗 ERDAS IMAGINE 8.4 图标面板菜单条：SessionT

51、ile Viewers 然后，在Viewer#1中打开校正后的图像：rectify.img 在Viewer#2中打开参考图像：panAtlanta.img 在Viewer#1中：按住右键快捷菜单Geo Link/Unlink 在Viewer#2中：点击左键建立与Viewer#1的连接 在Viewer#1中：按住右键快捷菜单Inquire Cursor打开光标查询对话框 在Viewer#1中：移动查询光标,观测其在两屏幕中的位置及匹配程度,并注意光标查询对话框中数据的变化。如果满意的话，关闭光标查询窗口。 实验五 图像增强处理一、 实验目的掌握在Erads Imagine下进行直方图调整、彩色变

52、换、图像运算、多光谱变换、滤波运算、分辨率融合等图像增强处理方法。二、 实验数据lanier.img、wasia1_mss.img、wasia2_mss.img、dmtm.img、spots.img、tmatlanta.img三、 实验步骤（一）直方图均衡化 （以d:program files imagine 8.4exampleslanier.img为例） 直方图均衡化实质上是对图像进行非线性拉伸，重新分配图像像元值，使一定灰度范围内的像元的数量大致相等；这样，原来直方图中间的峰顶部分对比度得到增强，而两侧的谷底部分对比度降低，输出图像的直方图是一较平的分段直方图。 ERDAS IMAGIN

53、E 8.4图标面板菜单条:MainImage Interpreter(或单击ERDAS IMAGINE 8.4图标面板工具条“Interpreter”图标) 打开Image Interpreter对话框选择Radiometric Enhancement 打开Radiometric Enhancement对话框选择Histogram Equalization 打开Histogram Equalization对话框在Histogram Equalization对话框中，设置下列参数： Input File（确定输入文件）：lanier.img Output File（定义输出文件）：equaliz

54、ation.img Coordinate Type（文件坐标类型）：Map Subset Definition（处理范围确定）：ULX,ULY;LRX，LRY（缺省状态为整个图像范围） Number of Bins（输出数据类型）：256 选定Ignore Zero in Stats(输出数据统计时忽略零值)单击OK（关闭Histogram Equalization对话框，执行直方图均衡化）（二）直方图匹配 （以d:program files imagine 8.4exampleswasia1_mss.img和wasia2_mss.img为例） 直方图匹配是对图像查找表进行数学变换，使一幅图像的直方图与另一幅图像类似，直方图匹配经常作为相邻图像拼接或应用多时相遥感图像进行动态变化研究的预处理工作，通过直方图匹配可以消除由于太阳高度角或大气影响造成的相邻图像的效果差异。 ERDAS IMAGINE 8.4图标面板菜单条:MainImage I