实验二、遥感图像的几何校正

《实验二、遥感图像的几何校正》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验二、遥感图像的几何校正实验目的：通过实习操作，掌握遥感图像儿何校正的基本方法和步骤，深刻理解遥感图像几何校正的意义。实验内容：ERDAS软件中图像预处理模块下的图像儿何校正。几何校正就是将图像数据投影到平面上，使其符合地图投影系统的过程。而将地图投影系统赋予图像数据的过程，称为地里参考(Geo-referencing)o由于所冇地图投影系统都遵循一定的地图处标系统，因此几何校正的过程包含了地理参考过程。1、图像几何校正的途径ERDAS图标血板I-具条：点击DataPrep图标，ImageGeometricConection

2、打开SetGeo-CorrectionInputFile对话框（图2・1）。ERDAS图标血板菜单条：MainfDataPreparation-*ImageGeometricCorrectionf打开SetGeo-CorrectionInputFile对话框（图2-1）o图2-1SetGeo-CorrectionInputFile对话框在SetGeo-CorrectionInputFile对话框(图1)中，需要确定校正图像,有两种选择情况：其一：首先确定來自视窗(FromViewer),然后选择显示图像视窗。其二：首先确定来自

3、文件(FromImageFile),然后选择输入图像。2、图像几何校正的计算模型(GeometricCorrectionModel)ERDAS提供的图像几何校正模型有7种，具体功能如下:表2・1几何校正计算模型与功能模型功能Affine图像仿射变换(不做投影变换)Polynomial多项式变换(同时作投影变换)Reproject投影变换(转换调用多项式变换)RubberSheeting非线性变换、非均匀变换Camera航空影像正射校正LandsatLantsat卫星图像正射校正SpotSpot卫星图像正射校正3、图像校正的具体

4、过程第一步：显示图像文件(DisplayImageFiles)首先，在ERDAS图标面板111点击Viewer图衣两次,打开两个视窗(Viewer1/Viewer2),并将两个视窗平铺放置，操作过程如下：ERDAS图农而板菜单条：Session-*TitleViewers然后，在Viewer1中打开盂要校-正的Landsat图像：vt87_sub2.img在Viewer2屮打开作为地理参考的校止过的Landsat图像：ws87rs.img第二步：启动几何校正模块(GeometricCorrectionTool)Viewer1

5、菜单条：Raster—GeometricCorrectionf打开SetGeometricModel对话框(2)选择多项式儿何校正模型：Polynomial-*OKf同时打开GeoCorrectionTools对话框(3)和PolynomialModelProperties对话框(4)0在PolynomialModelProperties对话框中，定义多项式模型参数以及.投影参数：定义多项式次方(PolynomialOrder):2f定义投影参数：(PROJECTION):略-*Apply-*Close打开GCPToolRef

6、erenseSetup对话框(5)WSetGeo>etricIodelSelectGeometricModel:AffinePolynomialRubberSheetingCameraLandsatSpotOpenExistingModel...UseExistingCalibrationCancel

7、Help图2-2SetGeometricModel对话框图2-3GeoCorrectionTools对话框图2-4PolynomialProperties对话框图2-5GCPToolReferenseSetup对话框第三步：启动

8、控制点工具(StartGCPTools)图2~6ViewerSelectionInstructions首先,在GCPToolReferenseSetup对话框（图5）中选择釆点模式:—选择视窗采点模式：ExistingViewer-*OKf打开ViewerSeiectionInstructions指示器（图2-6）f在显示作为地理参考图像ws87_rs.img的Viewer2屮点击左键f打开referenceMapTnformation提示框（图2-7）；-►OK-此时，整个屏幕将白动变化为如图7所示的状态，表明控制点工具被启

9、动，进入控制点采点状态。必.ReferenceMapInformation[凶CurrentReferenceMapProjection:Projection:UTMSpheroid:KrasovskyZoneNumber:50Datum:KrasovskyMapUnits: