遥感实习报告毕业论文.pdf

《遥感实习报告毕业论文.pdf》由会员分享，可在线阅读，更多相关《遥感实习报告毕业论文.pdf（22页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、-遥感实习报告 实验课程名称:ERDAS IMAGINE 8.6 系统的数据预处理 实 验 时 间:实验指导老师:学 号:班 级:姓 名:-目录 实习 1 ERDAS IMAGINE 8.6 系统简介与入门.2 练习 2 数据的输入/输出简介.6 实习 2 ERDAS IMAGINE 8.6 系统的数据预处理.10 练习 2 图像几何校正（GEOMETRIC CORRECTION）.12 练习 3 卫星图像拼接（SATELLITE IMAGE MOSAIC）.18 练习 4 图像投影变换（REPROJECT IMAGES）.20 实习总结：.22 -实习 1 ERDAS IMAGINE 8.6

2、 系统简介与入门 1 遥感的应用领域和 ERDAS IMAGINE 系统的基本概念和功能：遥感技术作为对地观测，提取地表最现势状况的最有利工具，被广泛应用在各行各业，包括测绘，自然资源管理，林业，水利，交通，环境保护，电力电信、防震减灾、城市规划、国防军事等，ERDAS IMAGINE 作为遥感界的排头兵，为大家提供了不仅仅是增强，滤波，纠正，融合等简单的基本应用，而是提供了强大的工具，使你在定量化的分析方面，系统功能的可扩充性方面使您更加得心应手，如专家分类，子象元分类（混合象元），三维可视化分析，数字摄影测量等，同时，还给您带来与 GIS 一化集成的解决方案，如查询检索编辑ArcInfo

3、的地理信息，建立矢量层后的人工解译，直接得到目标的矢量数椐，还可将分好类的专题影像转换成 ArcInfo 的矢量数据（Coverage,Shape File）,使分析的结果可以直接为地理信息系统管理与应用，从而发挥更大的作用。ERDAS IMAGINE 是美国 ERDAS 公司开发的专业遥感图像处理与地理信息系统软件。ERDAS IMAGINE 是以模块化的方式提供给用户的，可使用户根据自己的应用要求、资金情况合理地选择不同功能模块及其不同组合，对系统进行剪裁，充分利用软硬件资源，并最大限度地满足用户的专业应用要求。ERDAS IMAGINE 面向不同需求的用户，对于系统的扩展功能采用开放的体

4、系结构，以 IMAGINE Essentials、IMAGINE Advamage、IMAGINE Professional 的形式为用户提供了低、中、高三档产品架构，并有丰富的功能扩展模块供用户选择，使产品模块的组合具有极大的灵活性，1.1 IMAGINE Essentials 级 是一个花费极少的，包括制图和可视化核心功能的图像工具软件无论您是独立地从事工作或是处在企业协同计算的环境下，都可以借助 IMAGINE Essentials 完成二维三维显示、数据输入，排序与管理、地图配准，专题制图以及简单的分析。可以集成使用多种数据类型，井在保持相同的易于使用和易于剪裁的界面下升级到其它的 E

5、RDAS 公司产品。1.3 IMAGINE Professional 级 是面向从事复杂分析，需要最新和最全面处理工具，经验丰富的专业用户。Professional是功能完整丰富的图像地理信息系统。除了 Essentials 和 Advantage 中包含的功能以外，IMAGINE Professional 还提供轻松易用的空间建模工具（使用简单的图形化界面），高级的-参数/非参数分类器，知识工程师和专家分类器，分类优化和精度评定，以及雷达图像分析工具 2、ERDAS IMAGINE 软件的主要特点 为什么我们选择 ERDAS IMAGINE 作为学习遥感信息技术处理的操作软件？1 图像处理方

6、面 1)方便和直观的操作步骤使用户操作非常灵活：ERDAS IMAGINE 具有非常友好、方便地管理多窗口的功能。不论是几何校正还是航片、卫片区域正射矫正以及其它与多个窗口有关的功能，IMAGINE 都将相关的多个窗口非常方便地组织起来，免去了用户开关窗口、排列窗口、组织窗口的麻烦，应用方便因而加快了产品的生产速度。IMAGINE 的窗口提供了卷帘、闪烁、设置透明度以及根据坐标进行窗口联接的功能，为多个相关图像的比较提供了方便的工具。IMAGINE 的窗口还提供了整倍的放大缩小、任意矩形放大缩小、实时交互式放大缩小、虚拟及类似动画游戏式漫游等工具，方便对图像进行各种形式的观看与比较。2)ERD

7、AS IMAGINE 为不同的应用提供了 250 多种地图投影系统。支持用户添加自己定义的坐标系统。支持不同投影间的实时转换、不同投影图像的同时显示对不同投影图像直接进行操作等。支持相对坐标的应用。另外有非常方便的坐标转换工具，经纬度到大地坐标，反之亦然。3)常用的图像处理算法都可用图形菜单驱动，用户也可指定批处理方式（batch），使图像处理操作在用户指定的时刻开始执行；4)图像的处理过程可以由图像的属性信息控制，而上层属性信息可存在于本层或任何其他数据层次；5)图像处理过程可以用于具有不同分辨率的图像数据上，输出结果的分辨率可由用户指定；6)支持对不同图像数据源的交集、并集和补集的图像处理

8、；7)图解空间建模语言，EML 和 C 语言开发包的应用使得解决应用问题的客户化更加容易与简单。用户可以对 IMAGINE 本身应用的功能进行客户化的编辑，满足自己专业的独特需求。还可以将自己多年探索、研究的成果及工作流程以模型的形式表现出来。模型既可以单独运行也可以和界面结合象其它功能一样运行。更可以利用 C Toolkit 进行新型算法及功能的开发。8)独一无二的专家工程师及专家分类器工具，为高光谱、高分辨率图像的快速高精度分类提供了可能。此工具突破了传统分类只能利用光谱信息的局限，可以利用空间信息辅助分类。-此工具可以将我们所积累的几乎所有数字信息应用于分类，是分类应用的一大飞越。其功能

9、强大且应用方便，其提供的游标功能使知识库的优化成为轻而易举的操作。其知识库的可移动性为其它非专业人员进行分类工作提供了方便，为成熟知识库的推广应用提供了方便易行的途径。利用专家的知识还可以建立决策支持系统，为决策人提供工具。2 与地理信息系统的集成方面 ERDAS IMAGINE 系统已经内含了 ArcInfo Coverage 矢量数据模型，可以不经转换地读取、查询、检索其 coverage、GRID、SHAPEFILE、SDE 矢量数据，并可以直接编辑 coverage、SHAPEFILE数据。如果ERDAS IMAGINE再加上扩展功能，还可实现GIS的建立拓扑关系、图形拼接、专题分类图

10、与矢量二者相互转换。节省了工作流程中让人头疼、费时费力的数据转换工作，解决了信息丢失问题，可大大提高工作效率，使遥感定量化分析更完善。1、数据叠加显示 数据叠加显示（Blend,Swipe,Flicker）是针对具有相同的地理参考系统（地图投影和坐标系统）的两个文件进行操作的，所以，在进行数据叠加操作之前，首先需要按照在一个窗口中同时打开两个文件，需要说明的是：在打开第 2 个文件的时候，一定要在 Raster Options或者 Vector Options 中设置不清除窗口中已经打开的文件（取消选中 Clear Display 复选框）。见图 1.4 1 叠加数据准备（Prepare Ov

11、erlay Data）视窗菜单条:File-Open-Raster Layer-Select Layer To Add 对话框（图 1.3）视窗工具条:鼠标左键点击“打开文件”图标-Select Layer To Add 对话框（图 1.3）（1）确定文件路径（Look in）:%sampledata/openimage（2）确定文件类型（Files of type）:IMAGINE Image（*.img）（3）选择文件名称（File Name）:lanier.img（4）单击 OK（关闭 Select Layer to Add 对话框,打开下层图像文件）视窗工具条:鼠标左键点击“打开文件”

12、图标-Select Layer To Add 对话框（图 1.3）（1）确定文件路径（Look in）:%sampledata/openimage（2）确定文件类型（Files of Type）:IMAGINE Image（.img）（3）选择文件名称（File Name）:lnlandc.img（4）在 Raster Options 栏目中:不选择 Clear Display-（5）单击 OK（关闭 Select Layer to Add 对话框，打开上层图像文件）2 叠加显示操作（Operate Overlay Display）IMAGINE 系统所提供的数据叠加显示工具有三个，分别是混合

13、显示工具（Blend Tool）,卷帘显示工具（Swipe Tool）、和闪烁显示工具（Flicker Tool），都集成在实用菜单中。（1）混合显示工具（Blend Tool）本操作通过控制上层图像显示的透明度大小，使得上下两层图像混合显示.视窗菜单条:Utility-Blend-Viewer Blend/Fade 对话框（图 1.5）:在 Viewer Blend/Fade 对话框中，用户既可以通过设置 Blend/Fade Percentage（0-100）达到混合显示效果，也可以通过定义 Speed 和选择 Auto Mode 自动显示文件混合效果。2、文件显示顺序的问题 在实际工作中

14、，经常需要在同一个视窗中同时打开多个文件，包括图像文件、图形文件、AOI 文件、注记文件等，可以应用 Arrange Layers 命令调整文件显示顺序。为了说明文件显示顺序操作功能，需要首先在视窗中依次打开一组图像文件（lnlandc.img,Lanier.img,lndem.img,lnlakes.img），注意打开上层图像时，不要清除视窗 y 中已经打开的图像。依次打开%sampledata/openimage 下的以上文件 视窗菜单条:View-Arrange Layers-Arrange Layers Viewer 对话框（图 1.9）在 Arrange Layers Viewer

15、对话框中点击鼠标左键或拖动文件，达到调整文件顺序之目的，然后应用（Apple）显示顺序调整，并关闭（Close）Arrange Layers Viewer 对话框，结束文件显示顺序操作。显示比例操作（Display Scale）用于调整文件显示比例及其与视窗的对比关系，Scale 菜单对应的二级下拉菜单中包括四个，命令依次是：（1）Image to Window:按照视窗大小调整文件显示比例；（2）Window to Image:按照文件尺寸调整视窗大小；（3）Extent:显示文件整体范围；（4）Scale Tool:通过比例工具定义显示比例。练习 2 数据的输入/输出简介 ERDAS IM

16、AGINE 的数据输入输出功能（import/Export），允许您输入多种格式的数据供 IMAGINE 使用，同时允许您将 IMAGINE 的文件转换成多种数据格式。目前，IMAGINE-可以输入的数据格式达 70 多种，可以输出的数据格式近 34 中，几乎包括常用或常见的栅格数据和矢量数据格式，具体的数据格式都罗列在 IMAGINE 输入输出对话框中（图 1.10）。数据输入输出的一般操作过程如下:（在此我们以地区的 TM 多光谱卫星数据输入为例，演示如何操作数据的输入/输出工作。首先我们获得的数据是刻录在光盘上的，为了工作方便，我们先把数据复制到硬盘上你的工作目录%sampledata/

17、im-export/19860730，复制完成后我们发现，TM 数据是分波段存放的。根据需要，我们有输入单波段数据和组合多波段数据两种方法。1、输入单波段数据 ERDAS图 标 面 板 菜 单条:Main-Import/Export-数据输入输出对话框（图 1.10）在数据输入输出对话框中，通常需要设置下列参数信息:、对于这种数据，必须按照普通二进制来（Generic Binary）输入，而不能按照TM图像或SPOT图像来输入。同时，虽然数据文件是存储在只读光盘（CD-ROM）或磁带（Tape）中，但为了提高数据转换速度、并保证转换质量，最好是将数据文件直接复制到计算机硬盘中，而后选择文件（F

18、ile）作为输入媒体、而不要选择CD-ROM或Tape。2、组合多波段数据（LaverStack Bands Data）上面的数据输入只是将单波段的普通二进制数据文件转换成 ERDAS 自己的单波段IMG 文件，而在实际工作中，对遥感图像的处理和分析都是针对多波段图像进行的，所以，还需要将若干单波段图像文件组合（Laver Stack）成一个多波段图像文件，具体过程如下:ERDAS 图标面板菜单条:Main-Image Interpreter-Utilities-Layer Stack 打开 Layer Selection and Stacking 对话框（图 1.11）ERDAS 图标面板工

19、具条:点击 Interpreter图标-Utilities-Layer Stack 打开 Layer Selection and Stacking 对话框（图 1.11）在 Layer Selection and Stacking 对话框中，依次选择并加载（Add）单波段图像:图 1.10 Import/Export 图标数据输入输出对话框 -（1）输入单波段文件（InputFile:*.img）:Band1.img-Add（2）输入单波段文件（Input File:img）:Band2,img-Add（3）输入单波段文件（InputFile:*.img）:Band3,img-Add （4）输

20、出多波段文件（Output File:*.img）:bandstack,img（5）输出数据类型（Output Data Type）:Unsigned 8 Bit（6）波段组合选择（OutputOption）:Union（7）输出统计忽略零值:Ignore Zero In Stats（8）单击 OK（关闭 Layer Selection and Stacking对话框，执行波段组合）我们就得到一个具有 7 个波段的多光谱卫星影像，保存到工作目录下，文件名为：19860730tm.img。3.TIFF 图像数据输入输出（Import/Export TIFF Data）TIFF 图像数据是非常通用

21、的图像文件格式，从 ERDAS IMAGINE8.4 起增加了一个 TIFF DLL动态连接库，从而使ERDAS IMAGINE支持6.0版本的TIFF图像数据格式的直接读写，包括普通 TIFF 和 GeoTIFF。用户在使用TIFF图像数据时，不需要再象以前那样通过Import/Export来转换TIFF文件，而是只要在打开图像文件时，将文件类型指定为 TIFF 格式就可以直接在视窗中显示 TIFF图像。如果要在图像解译器（Interpreter）或其它模块下对图像做进一步的处理操作，依然需要将 TIFF 文件转换为 IMG 文件，这种转换非常简单，只要在打开 TIFF 的视窗中将 TIFF

22、 文件另存为（Save As）IMG 文件就可以了。同样，如果 ERDAS IMAGINE 的 IMG 文件需要转换为 GeoTIFF 文件，只要在打开 IMG 图像文件的视窗中将 IMG 文件另存为 TIFF 文件就可以了。4.输出JPEG图像数据（Export JPEG Data）JPEG 图像数据是一种通用的图像文件格式，ERDAS 可以将自己的 IMG 图像文件输出成 JPEG 图像文件，供其它图 1.11 Layer Selection and Stacking 对话框 图 1.12 Export JFIF Data 对话框-图像处理系统或办公软件使用，在 ERDAS 图标面板菜单条

23、:Main-Import/Export 输入输出对话框中（如图 1.12）（1）选择输出数据操作:Export（2）选择输出数据类型（Type）为 JPEG:JFIF（JPEG）（3）选择输出数据媒体（Media）为文件:File（4）确定输入文件路径和文件名为 19860730tm.img（5）确定输出文件路径和文件名为 19860730tm.jpg，单击 OK 按钮（关闭 Import/Export 对话框），打开 Export JFIF Data 对话框（如图 1.13）。在 Export JFIF Data 对话框中设置下列输出参数：设置图像对比度调整（Contrast Options

24、）为 Apply Standard Deviation Stretch；设标准差拉伸倍数（Standard Deviations）为 2；设图像转换质量（Quality）为 100；在 Export JFIF Data 中，定义下列参数：选择波段（Select Layers）为 4,3,2；坐标类型（Coordinate Type）为 Map；定 义 子 区（Subset Definition）为ULX,ULY,LRX,LRY；单击 OK 按钮（关闭 Export Options 对话框，结束输出参数定义）；返回 Export JFIF Data 对话框；（6）单击 OK 按钮（关闭 Expo

25、rt Options 对话框，执行 JPG 数据输出）。图 1.13 Export Options 对话框-实习 2 ERDAS IMAGINE 8.6 系统的数据预处理 1、规则分幅裁剪（Rectangle Subset Image）规则分幅裁剪是指裁剪图像的边界范围是一个矩形，通讨左下角和右下角两点的坐标，就可以确定图像的裁剪位置，整个裁剪过程比较简单。（1）ERDAS 图标面板菜单条:Main-Data Preparation-Data Preparation 菜单（图 2.1.）或 ERDAS 图标面板工具条:点击 Data Prep 图标-打开 Data Preparation 菜单

26、（图 2.1）选择 Subset Image-打开 Subset Image 对话框（图 2.1）在 Subset Image 对话框中需要设置下列参数:2、不规则分幅裁剪（Polygon Subset Image）在实际工作中，我们经常遇到的问题是不规则的边界剪裁问题。不规则分幅裁剪是指裁剪图像的边界范围是个任意多边形，无法通过左上角和右下角两点的坐标确定图像的裁剪位置，而必须事先生成一个完整的闭合多边形区域，可以是一个 AOI 多边形，也可以是 ArcInfo 的一个 Polygon Coverage，针对不同的情况采用不同裁剪过程。2、ArcInfo 多边形裁剪（Polygon Cove

27、rage Subset Image）如果是按照行政区划边界或自然区划边界进行图像的分幅裁剪，往往是首先利用ArcInfo 或 ERDAS 的 Vector 模块绘制精确的边界多边形（Polygon），然后以.ArcInfo 的Polygon 为边界条件进行图像裁剪。对于这种情况，需要调用 ERDAS 其它模块的功能分两步完成。（1）将 ArcInfo 多边形转换成栅格图像文件（有 4 种方式可以打开）ERDAS图标面板菜单条:Main-Image Interprete:-Utilities-Vector to Raste 图 2.1 Subset Image 对话框 图 2.2 绘制多边形 A

28、OI 对话图 2.6 Vector to Raster 对话框-打开 Vector to Raster 对话框（图 2.6）ERDAS 图标面板菜单条:Main-Vector-Vector to Raster，打开 Vector to Raster对话框（图 2.6）ERDAS 图标面板工具条:点击 Interpreter 图标-Utilities-Vector to Raster 打开 Vector to Raster 对话框（图 2.6）ERDAS 图标面板工具条:点击 Vector 图标-Vector to Raster（2）打开 Vector to Raster 对话框（图 2.6）在

29、 Vector to Raster 对话框中需要设置下列参数:输入矢量文件名称（Input Vector File）:zone88；确定矢量文件类型（Vector Type）:Polygon；输出栅格文件名称（Output Image File）:raster.img；栅格数据类型（Data Type）:Unsigned 8 bit；栅格文件类型（Layer Type）:Continuous；转换范围大小（Size Definition）:ULX,ULY,LRX,LRY；坐标单位（Units）:feet；输出像元大小（Cell Size）:X:30/Y:30；选择正方形像元:Squire Ce

30、ll（3）OK（关闭 Vector to Raster 对话框，执行矢栅转换）注意：1、为保证掩膜制作的精确和准确，在Vector 转为 Raster之前，一定要把zone88添加投影信息使之与待掩膜的图象的投影信息保持一致。如在本例子中，projection type:State Plane,State plane zone:4126,Zone type:Old USGS（D0154）,NAD:NAD27。2、同时需要注意，由于germtm.img的单位是feet，所以需要在投影变换的时候把raster的单位也调整为feet。(在viewer中打开zone88,点击工具栏第三个按钮info工

31、具，在vectorinfo对话框中点edti，打开add projection).（4）通过掩膜运算（Mask）实现图像不规则裁剪（5）OK（关闭 Mask 对话框，执行掩膜运算）打开多边形剪切后的图像 mask.img 如图 2.8 所示。-练 习2 图 像 几 何 校 正（Geometric Correction）几何校正就是将图像数据投影到平面上，使其符合（Conform）地图投影系统的过程;而将地图坐标系 统 赋 予 图 像 数 据 的 过 程，称 为 地 理 参 考（Geo-referencing）。由于所有地图投影系统都遵从于一定的地图坐标系统，所以几何校正过程包含了地理参考过程。

32、数据存放在%sampledata/gc-image 下。1、图像几何校正途径（Geometric Correction Process）在 ERDAS IMAGINE 系统中进行图像几何校正，通常有两种途径启动几何校正模块。（1）数据预处理途径（Start from Data Preparation）ERDAS 图标面板菜单条:Main-Data Preparation-Image Geometric Correction,打开 Set Geo-Correction Input File 对话框（图 2.9）ERDAS 图标面板工具条:点击 Data Prep 图标-Image Geometr

33、ic Correction 打开 Set Geo-Correction Input File 对话框（图 2.9）2、几何校正计算模型（Geometric Correction Model）如图所示，ERDAS 提供的图像几何校正计算模型有 10 种，具体功能如表 4.1 所列:以资源卫星图像校正为例，介绍的是以已经具有地理参考的 SPOT 图像为基础，进行Landsat TM 图像校正过程。（1）显示图像文件 首先，在 ERDAS 图标面板中点击 Viewer 图标两次，打开两个视窗Viewer#I/viewer#2），并将两个视窗平铺放置，操作过程如下:ERDAS 图标面板菜单条:Sess

34、ion-Tile Viewers 然后，在 Viewer#1 中打开需要校正的 Landsat TM 图像:tmAtlanta.img 在Viewer 图2.9 Set Geo-Connection Input File 对话框 图 2.10 打开两幅图像后的窗口-#2 中打开作为地理参考的校正过的 SPOT 图像:panAtlanta.img:（2）启动几何校正模块（Geometric Correction Tool）Viewer#1 菜单条:Raster-Geometric Correction，打开 Set Geometric Model 对话框（图 2.11）选择多项式几何校正计算模型

35、:Polynomial OK，同时打开 Polynomial Model Properties 对话框（图 2.12）和Geo Correction Tools 对话框（图 2.13）在 Polynomial Model Properties 对话框中，定义多项式模型参数及投影参数:1）定义多项式次方（Polynomial Order）；2）定义投影参数（Projection）（略）；（3）Apply-Close，打开 GCP Tool Reference Setup 对话框（图 2.13）说明:该实例是采用视窗采点模式，作为地理参考的SPOT图像已经含有投影信息，所以，这里不需要定义投影参数

36、。如果不是采用视窗采点模式，或者参考图像没有包含投影信息，图 2.11 Set Geometric Model 对话框 图 2.12 Polynomial ModelProperd 对话框 图2.13 Geo Correction Tools对话框-则必须在这里定义投影信息，包括投影类型及其对应的投影参数。（4）启动控制点工具（Start GCP Tools）首先在 GCP Toots Reference Setup 对话框（图 2.13）中选择采点模式:选择视窗采点模式:Existing Viewer，单击 OK（关闭 GCP Tools Reference Setup 对话框），自动打开

37、Viewer Selection Instiuctions 指示器（图 2.14）（5）在显示作为地理参考图像 panAtlanta.img 的 Viewer#2中点击左键，打开 Reference Map Information 提示框（图 2.15）（显示参考图像的投影信息），OK（关闭 Reference Map Information 提示框）（6）整个屏幕将自动变化为如图 2.16 所示的状态:其中包含两个主视窗、两个放大窗口、两个关联方框（分别位于两个视窗中，指示放大窗口与主视窗的关系）、控制点工具对话框、几何校正工具等。表明控制点工具被启动，进入控制点采集状态。（7）采集地面控制

38、点（Ground Control Point）说明：控制点工具对话框简介（Introduction to GCP Tool）在正式开始采集控制点之前，首先对控制点工具对话框进行说明GCP工具对话框（GCP Tool）由菜单条、工具条、控制点数据表（GCP Cell Array）及状态条（Status Bar）4个部分组成，菜单条中菜单命令及其功能详见软件菜单说明。关于GCP工具对话框，还需要说明几点:（a）输入控制点（Input GCP）的是在原始文件视窗中采集的，具有原文件的坐标系统;而参考控制点（Reference GCP）是在参考文件视窗中采集的，具有己知的参考坐标系统，GCP图 2.1

39、4 Viewer Selection Instiuctions 指示器 图 2.16Reference Map Information 窗口 图 2.15Reference Map Information 提示框-工具将根据对应点的坐标值自动生成转换模型。（b）在GCP数据表中，残差（Residuals）、中误差（RMS）、贡献率（Contribution）及匹配程度（Match）等参数，是在编辑GCP的过程中自动计算更新的，用户是不可以任意改变的，但可以通过精确GCP位置来调整。（c）每个IMG文件都可以有一个GCP数据集与之相关联，GCP数据集保存在一个栅格层数据文件中;如果IMG文件有一

40、个GCP数据集存在的话，只要打开GCP工具，GCP点就会出现在视窗中。（d）所有的输入GCP都可以直接保存在图像文件中（Save Ii）put），也可以保存在控制点文件中（Save Input As）。如果是保存在文件中，调用的方法如（c）所述，如果是保存在GCP文件中，可以通过加载调用（Load Input）。（e）参考GCP也可以类似地保存在参考图像中（Save Reference）或GCP文件中（Save Reference As），便于以后调用。（8）GCP 的具体采集过程（Steps of GCP Selection）在图像几何校正过程中，采集控制点是一项非常重要和相当繁重的工作，具

41、体过程如下:（a）在 GCP 工具对话框中点击 Select GCP 图标，进入 GCP 选择状态;（b）在 GCP 数据表中将输入 GCP 的颜色（Color）设置为比较明显的黄色，-（c）在 Viewer#1 中移动关联方框位置，寻找明显的地物特征点，作为输入 GCP:（d）在 GCP 工具对话框中点击 Create GCP 图标，并在 Viewer#3 中点击左键定点，GCP.数据表将记录一个输入 GCP，包括其编号、标识码、X 坐标、Y 坐标:（e）在 GCP 工具对话框中点击 Select GCP 图标，重新进入 GCP 选择状态;（f）在 GCP 数据表中将参考 GCP 的颜色（C

42、olor）设置为比较明显的红色:（g）在 Viewer#2 中移动关联方框位置，寻找对应的地物特征点，作为参考 GCP:（h）在 GCP 工具对话框中点击 Create GCP 图标，并在 Viewer#4 中点击左键定点，系统将自动把参考点的坐标X Reference,Y Reference显示在 GCP 表中，（i）在 GCP 工具对话框中点击 Select GCP 图标，重新进入 GCP 选择状态:并将光标移回到Viewer#1，准备采集另一个输入控制点。（j）不断重复（a）-（i），采集若干 GCP,直到满足所选定的几何校正模型为止；而后，每采集一个Input GCP,系统自动产生 一

43、个Ref.GCP,通过移动Ref.GCP 可以逐步优化校正模型。采集 GCP 以后，GCP 数据表如图 2.17 而示:（实际数据以所采集的GSP的点位不同而不同。）-（9）采集地面检查点 以上所采集的 GCP 的类型为 Control Point（控制点），用于控制计算，建立转换模型及多项式方程。下面所要采集的 GCP 的类型均是 Check（检查点），用于检验所建立的转换方程的精度和实用性。依然在 GCP Tool 对话框状态下:（a）在 GCP Tool 菜单条中确定 GCP 类型:Edit-Set Point Type-Check（b）在 GCP Tool 菜单条中确定 GCP 匹配参

44、数（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

45、 All/Selected Point）选项组中设从输入到参考（Reference from Input）或从参考到输入（Input from Reference）（d）Close（关闭 GCP Matching 对话框）（e）确定地面检查点:在GCP Tool工具条中选择 Create GCP 图标，并将 Lock 图标打开，锁住 Create GCP 功能，如同选择控制点一样，分别在 Viewer#l 和 Viewer#2 中定义 5 个检查点，定义完毕后点击 Unlock 图标，解除 Create GCP 功能。（f）计算检查点误差:在 GCP Tool 工具条中点击 Compute E

46、rror图标，检查点的误差就会显示在 GCP Tool 的上方，只有所有检查点的误差均小于一个像元（Pixel），才能继续进行合理的重采样。一般来说，如果你的控制点（GCP）定位选择比较准确的话,检查点匹配会比较好，误差会在限差范围内。否则，若图 2.17GCP Tool 对话框与 GCP 数据表 图2.18 GCP Matching对话框-控制点定义不精确，检查点就无法匹配，误差会超标。（10）计算转换模型（Compute Transformation）:在控制点采集过程中，一般是设置为自动转换计算模式，所以，随着控制点采集过程的完成，转换模型就自动计算生成，下面是转换模型的查阅过程:（a）

47、在 Geo-Correction Tools 对话框中点击 Display Model Properties 图标 （b）打开 Polynomial Model Properties（多项式模型参数）对话框（图 2.12）（c）在多项式模型参数对话框中查阅模型参数，并记录转换模型 （d）Close（关闭模型特性对话框，进入图像重采样阶段）（11）图像重采样（Resample the Image）图像重采样简介（Introduction to image Resample）重采样过程就是依据未校正图像像元值计算生成一幅校正图像的过程，原图像中所有栅格数据层都将进行重采样。ERDAS IMAGIN

48、E提供三种最常用的重采样方法:（1）Neatest Neighbor:邻近点插值法，将最邻近像元值直接赋予输出像元;（2）Bilinear Interpolation:双线性插值法，用双线性方程和2X2窗口计算输出像元值;（3）Cubic Convolution:立方卷积插值法，用立方方程和4X4窗口计算输出像元值;图像重采样过程（Process of Image Resample）:（a）首先，在 Geo Correction Too、对话框中选择 Image Resampl。图标（b）打开 Image Resample（图像重采样）对话框（c）然后，在 Image Resample 对话框

49、中，定义重采样参数:输出图像文件名（Output File）:rectify.img;选择重采样方法（Resample Method）:Nearest Neighbor;定义输出图像范围（0output Corners）:ULX,ULY,LRX,LRY;定义输出像元大小（Output Cell Size）:X:30.Y:30;设置输出统计中忽略零值:Ignore Zero in Stats;设置重新计算输出缺省值（Recalculate Output Default）:Skip Factor.10 单击 OK（关闭 Image Resample 对话框，启动重采样进程）（12）保存几何校正模式

50、（Save Rectification Mode）在 Geo Correction Tools 对话框中点击 Exit 按钮，退出图像几何校正过程，按照系统提示选择保存图像几何校正模式，并定义模式文件（*.gm），以便下次直接使用。（13）检验校正结果（Verify Rectification Result）检验校正结果的基本方法是:同时在两个视窗中打开两幅图像，其中一幅是校正以后的-图像，一幅是当时的参考图像，通过视窗地理连接（Geo Link/Unlink）功能及查询光标（Inquirecursor）功能进行目视定性检脸，具体过程如下:（a）打开两个平铺视窗（open and Tile t