实验二、数据预处理

1、实习二、数据预处理一、预处理简介ERDAS IMAGING数据预处理模块是由一组实用的图像数据处理工具构成，包括生成单值图像（Creat New Image）、三维地形表面（Create Surface）、图像分幅裁剪（Subset Image）、图像几何校正（Image Geometric Correction）、图像拼接处理（Mosaic Images）、非监督分类（Unsupervised Classification）、以及图像投影变换（Reprojection Images）等，主要是根据工作区域的地理特征和专题信息提取的客观需要，对数据输入模块中获取的IMG图像文件进行范围调整、误

2、差校正、坐标转换等处理，以便进一步开展图像解译、专题分类等分析研究。数据预处理模块简称Data Preparation或DataPrep，可以通过两种途径启动：ERDAS图标面板菜单条：MainData PreparationData Preparation菜单（图2.1）ERDAS图标面板工具条：点击DataPrep图标Data Preparation菜单（图2.1）图2.1 Data Preparation菜单从图2.1可以看出，ERDAS IMAGIMG数据预处理模块包括了7项主要功能，其中第一项功能（生成单值图像）比较简单，第六项功能（非监督分类）将在图像分类中进行说明。下面将主要介绍

3、其余五项功能，重点是图像几何校正和图像拼接处理，因为这两项操作是从事遥感应用研究必须开展的基本工作过程。二、三维地形表面（3D Surfacing）三维地形表面工具允许用户在不规则空间点的基础上产生三维地形表面，所支持的输入数据类型包括：ASCII码点文件、ArcInfo的Coverage点文件和线文件，ERDAS IMAGING的注记数据层，以及栅格图像文件IMG。所有输入数据必须具有X、Y、Z值，三维地形表面工具所应用的TIN插值方法，所输出的是一个连续的栅格图像文件。每一个已知的空间点在输出的地形表面上保持Z值不变，而没有Z值的空间点，其输出表面的Z值是基于TIN其周围的已知点插值计算获

4、得的。在三维地形表面工具中提供了两种TIN插值方法：线性插值（Linera）与非线性插值（non-linear）。线性插值方法是应用一次多项式方程进行计算，输出的TIN三角面是一些有棱角的平面；非线性插值方法应用五次多项式方程进行计算，输出的是平滑表面，这种情况下，TIN三角面不是一个平面，而是具有弹性的曲面。线性插值方法速度快但结果简单，而非线性插值方法产生基于不规则分布数据集的非常连续的、圆滑的表面结果。1.启动三维地形表面（Surfacing Create Surface）ERDAS 图标面板菜单条：MainData Preparation菜单选择Create Surface打开3D S

5、urfacing对话框（图2.2）ERDAS 图标面板工具条：点击Data Pre图标，打开Data Preparation菜单选择Create Surface打开3D Surfacing对话框（图2.2）3D Surfacing对话框由菜单条（Menu Bar）、工具条（Tool Bar）和数据表格（Data CellArray）组成。菜单条主要由文件操作（File）、数据表操作（Row）和表面生成（Surface）菜单组成，而工具条则由读取数据（Read Point）文件、保存数据文件（Save Point）和生成地形表面（Perform Surfacing）图标组成。图2.2 3D Su

6、rfacing对话框（读入数据之后）2.定义地形表面参数（Defining Surface Parameters）3D Surfacing对话框菜单条：File Read 打开Read Point对话框（图2.3）图2.3 Read Point对话框在Read Point 对话框中需要定义下列参数：数据源文件类型（Source File Type）：ASCII File数据源文件名称（Source File Name）：Inpts.datOK（关闭Read Points对话框）打开Import Options对话框（图2.4）在Import Options栏目，需要定义下列参数：选择字段类型（

7、Field Type）：Delimited by Separation（分割字符）选择分割字符（Separation Charactor）：Comma（逗号分割）图2.4 Import Options 对话框（Field Definition栏目）每行结束字符（Row Terminator Charactor）：Return New Line（DOS）确定跳过几行（Number of Rows to Skip）：0（从头读）点击Input Preview标签，进入Input Preview栏目（图2.5）从Import Options对话框显示的原始数据可知，数据文件中的数据记录方式是一行一个

8、点，每一行数据包括点号、X坐标、Y坐标、Z坐标（高程值）四个字段，其中点号在此处读入数据时不需要，因此，必须在Import Options对话框的Colum Mapping中确定X、Y、Z与数据文件中字段的对应关系：Output Column Name：X对应Input Field Number：2Output Column Name：Y对应Input Field Number：3Output Column Name：Z对应Input Field Number：4OK（关闭Import Options对话框，读入数据）数据读到3D Surfacing Data CellArray中（图2.2）

9、图2.5 Import Options对话框（Input Preview栏目）如果需要，可以将读入的数据保存为ERDAS的Annotation Layer或ArcInfo的Coverage文件：3D Surfacing dialog菜单：File Save as Save as对话框（图2.6）确定输出文件类型（Output File Type）：Point Coverage确定输出文件名（Save as ArcInfo）：testpoint图2.6 Save as对话框3.生成三维地形表面（Creating 3D Surface）3D Surfacing对话框菜单条：Surface Surf

10、acing打开 Surfacing对话框（图2.7）3D Surfacing对话框工具条：点击Perform Surfacing图标 打开Surfacing对话框（图2.7）图2.7 Surfacing对话框在Surfacing对话框中需要设置下列参数：输出文件名称（Output File）：testsurface.img表面插值方法（Surfacing Method）：Linear Rubber Sheeting输出文件范围（Output Corners）：自动读取ULX、ULY、LRX、LRY输出像元大小（Output Cell Size）：X：30/Y：30输出像元形状：Square C

11、ells输出图像背景值（Background Value）：0输出统计忽略值（Ignore Zero in Output Stats）：0输出数据类型（Output Data Type）：Unsigned 16 bitOK（关闭Surfacing对话框，执行地形表面过程）2.显示三维地形表面（Display 3D Surface）由三维地形表面过程生成的图像文件，就是ERDAS系统的DEM文件，DEM文件常常用于三维图像（Image Drape）的显示或虚拟地理信息系统（Virtual GIS）操作。所以，可以在Viewer视窗中打开曲面所生成的地形表面文件testsurface.img文件，

12、显示其二维平面效果或通过Image Info查看其定量信息，也可以在其上叠加对应的具有相同投影系统的图像文件，显示其三维立体效果。三、图像分幅裁剪（Subset Image）在实际工作中，经常需要根据研究工作范围对图像进行分幅裁剪，按照ERDAS实现图像分幅裁剪的过程，可以将图像分幅裁剪分为两种类型：规则分幅和不规则分幅。1.规则分幅裁剪（Rectangle Subset Image）规则分幅是指裁剪图像的边界范围是一个矩形，通过左上角和右下角两点的坐标，就可以确定图像的裁剪位置，整个裁剪过程比较简单。ERDAS 图标面板菜单条：Main Data Preparation Data Prepa

13、ration菜单（图2.1）选择Subset Image开Subset Image对话框（图2.8）ERDAS 图标面板工具条：点击Data Prep图标打开 Data Preparation菜单（图2.1）选择Subset Image打开Subset Image对话框（图2.8）图2.8 Subset Image对话框在Subset Image对话框中需要设置下列参数：输入文件名称（Input File）：Lanier.img输出文件名称（Output File）：Lanier_sub.img坐标类型（Coordinate Type）：File裁剪范围（Subset Definition）：

14、输入ULX、ULY、LRX、LRY输出数据类型（Output Data Type）：Unsigned 8 bit输出统计忽略零值：Ignore Zero in Output Stats输出像元波段（Select Layers）：1:7（表示选择从第一波段到第七波段）2,3,2（表示选择2，3，2四个波段）OK（关闭Subset Image对话框，执行图像裁剪）说明：在上述图像裁剪过程中，裁剪范围是通过直接输入左上角坐标和右下角坐标定义的。此外，还可以通过两种方式定义裁剪范围：其一是应用查询框（Inquire Box），然后在Subset Image对话框中选择From Inquire Box功

15、能；其二是应用AOI，然后在Subset Image对话框中选择AOI功能，打开AOI对话框，并确定AOI区域来自图像视窗即可。输出波段的选择：一种是用冒号表示从第几波段到第几波段，另外一种是用逗号隔开表示仅仅选择这几个波段。2.不规则分幅裁剪（Polygon Subset Image）不规则分幅裁剪是指 裁剪图像的边界范围是个任意多边形，无法通过左上角和右下角两点的坐标确定图像的裁剪位置，而必须事先生成一个完整闭合多边形区域，可以是一个AOI多边形，也可以是ArcInfo的一个Polygon Coverage。针对不同的情况采用不同裁剪过程。AOI多边形裁剪（Polygon AOI Subs

16、et Image）首先在视窗中打开需要裁剪的图像，并应用AOI工具绘制多边形AOI，可以将多边形AOI保存在文件中（*.aoi），也可以暂时不退出视窗，将图像与AOI多边形保留在视窗中，然后：ERDAS 图标面板菜单条：Main Data Preparation Data Preparation菜单（图2.1）选择Subset Image开Subset Image对话框（图2.8）ERDAS 图标面板工具条：点击Data Prep图标打开 Data Preparation菜单（图2.1）选择Subset Image打开Subset Image对话框（图2.8）在Subset Image对话框中

17、需要设置下列参数：输入文件名称（Input File）：Lanier.img输出文件名称（Output File）：Lanier_sub.img应用AOI确定裁剪范围：点击AOI按钮打开选择AOI（Choose AOI）对话框（图2.9）图2.9 Choose AOI对话框在Choose AOI对话框中确定AOI的来源（AOI Source）：File或Viewer输出数据类型（Output Data Type）：Unsigned 8 bit输出统计忽略零值：Ignore Zero in Output Stats输出像元波段（Select Layers）：1:7（表示选择从第一波段到第七波段）

18、2,3,2（表示选择2，3，2四个波段）OK（关闭Subset Image对话框，执行图像裁剪）ArcInfo多边形裁剪（Polygon Coverage Subset Image）如果按照行政区划边界或自然区划边界进行图像的分幅裁剪，往往是首先利用ArcInfo或ERDAS的Vector模块绘制精确的边界多边形（Polygon），然后以ArcInfo的Polygon为边界条件进行图像裁剪。对于这种情况，需要调用ERDAS其他模块的操作分两步完成。第一步：将ArcInfo多边形转换成栅格图像文件ERDAS 图标面板菜单条：Main Image Interpreter Utilities Vec

19、tor to Raster 打开Vector to Raster对话框（图2.10）ERDAS 图标面板工具条：Main Vector Vector to Raster 打开Vector to Raster对话框（图2.10）ERDAS 图标面板工具条：点击Interpreter图标 Utilities Vector to Raster 打开Vector to Raster对话框（图2.10）ERDAS 图标面板工具条：点击IVectot图标 Vector to Raster 打开Vector to Raster对话框（图2.10）图2.10 Vector to Raster对话框（a）图2.

20、10 Vector to Raster对话框（b）在Vector to Raster对话框中需要设置下列参数：确定矢量文件名称（Input Vector File）：zone88确定输出文件名称（Output Raster File）：raster.img确定矢量文件类型（Vector Type）：Polygon转换范围大小（Subset Definition）：ULX、ULY、LRX、LRY输出像元大小（Pixel Size）：30栅格数据类型（Output Pixel Data Type）：Unsigned 8 bit使用矢量属性值（An Item as Pixel Value）：ZONE

21、88-ID栅格文件类型（Output Image Type）：Thematic是否创建图层金字塔（Create Pyramid Layers）输出统计忽略零值：Ignore Zero in Output StatsOK（关闭Vector to Raster对话框，执行失栅转换）第二步：通过掩模运算（Mask）实现图像不规则裁剪ERDAS 图标面板菜单条：Main Image Interpreter Utilities Mask 打开Mask对话框（图2.11）ERDAS 图标面板工具条：点击Interpreter图标 Utilities Maskr 打开Mask对话框（图2.11）图2.11

22、Mask对话框在Mask对话框中需要设置下列参数：确定输入文件名称（Input Vector File）：Lanier.img确定掩模文件名称（Input Mask File）：raster.img点击Set Recode设置裁剪区域内新值（New Value）为1，区域外取0值确定掩模区域做交集运算：Intersection输出图像文件名称（Output File）：mask.img输出数据类型（Output Data Type）：Unsigned 8 bitOk（关闭Mask对话框，执行掩模运算）实现图像不规则裁剪。四、图像几何校正（Geometric Correction）几何校正就是将

23、图像数据投影到平面上，使其符合（Conform）地图投影系统的过程，而将地图坐标系统赋予图像数据的过程，称为地理参考（Geo-referencing）。由于所有地图投影系统都遵从于一定的地图坐标系统，所以几何校正过程包含了地理参考过程。1.图像几何校正概述（Introduction to Geometric Correction）在正式开始介绍图像几何校正分幅和过程之前，首先对ERDAS图像几何校正过程中的几个普遍性的问题进行简单说明，以便于随后的操作。图像几何校正途径（Geometric Correction Process）数据预处理途径（Start from Data Preparati

24、on）ERDAS图标面板菜单条：Main Data Preparation Image Geometric Correction打开Set Geo-Correction Input File对话框（图2.12）ERDAS图标面板工具条：点击Data Prep图标 Image Geometric Correction打开Set Geo-Correction Input File对话框（图2.12）图2.12 Set Geo-Correction Input File对话框在Set Geo-Correction Input File对话框中，需要确定校正图像，有两种选择情况其一：首先确定来自视窗（

25、From Viewer），然后选择显示图像视窗（Select Viewer）打开Set Geometric Model对话框（图2.13）选择几何校正计算模型（Select Geometric Model）OK打开校正模型参数与投影参数设置对话框定义校正模型参数与投影参数Apply Close图2.13 Set Geometric Model对话框打开GCP Tool Reference Setup对话框（图2.14）确定采点模式，采点校正图2.14 GCP Tool Reference Setup对话框其二：首先确定来自文件（From Image File）然后选择输入图像（Input Im

26、age File）打开Set Geometric Model对话框（图2.13）选择几何校正计算模型（Select Geometric Model）OK打开校正模型参数与投影参数设置对话框定义校正模型参数与投影参数Apply Close打开GCP Tool Reference Setup对话框（图2.14）确定采点模式，采点校正对于第一种情况，必须事先在一个视窗中打开需要几何校正的图像，否则无法进行选择。视窗栅格操作途径（Start From Viewer Raster）这种途径是首先在一个视窗中打开需要校正的图像，然后在栅格操作菜单中启动几何校正模块，具体过程如下：视窗菜单条：RasterG

27、eometric Correction选择几何校正计算模型（Select Geometric Model）OK打开校正模型参数与投影参数设置对话框定义校正模型参数与投影参数Apply Close打开GCP Tool Reference Setup对话框（图2.14）确定采点模式，采点校正几何校正模块启动的两种途径相比，视窗栅格操作途径更为直观简便。2.几何校正计算模型（Geometric Correction Model）ERDAS提供的图像几何校正计算模型有7种，具体功能如表2.1所列：表2.1 几何校正计算模型与功能模型功能Affine图像仿射变换（不做投影变换）Polyomial多项式变

28、换（同时做投影变换）Reproject投影变换（转换调用多项式变换）Rubber Sheeting非线性、非均匀变换Camera航空影像正射校正LandsatLandsat卫星图像正射校正SPOTSPOT卫星图像正射校正其中，多项式变换在卫星图像校正过程中应用较多，在调用多项式模型时，需要确定多项式的次方数（Order），通常整景图像选择3次方，次方数与所需要的最少控制点数是相关的，最少控制点计算公式为（t+1）（t+2）/2，式中t为次方数，即一次方程需要3个控制点，二次方程需要6个控制点，三次方程需要10个控制点，依此类推。3.几何校正采点模式（Record Ground Control

29、Point）图2.12中列出了ERDAS系统提供的9种控制点采集模式，仔细分析这9种模式，可以归纳阿三大类，具体类型与含义如表2.2所列。表2.2 几何校正采点模式与含义模式含义Viewer to Viewer：Existing ViewerImage Layer（New Viewer）Vector Layer（New Viewer）Annotation Layer（New Viewer）视窗采点模式：在已经打开的视窗中采点在新打开的图像视窗中采点在新打开的矢量视窗中采点在新打开的注记视窗中采点File to Viewer：GCP File（*.gcc）ASCII File文件采点模式：在控制

30、点文件中读点ASCII码文件中读点Map to Viewer：Digitizing Tablet（Current）Digitizing Tablet（New）Keyboard Only地图采点模式：在当前数字化仪上采点在新配置数字化仪上采点通过键盘输入控制点表2.2所列的三种几何校正采点模式，分别应用于不同的情况：如果已经拥有需要校正图像区域的数字地图，或经过校正的图像，或注记图层的话，就可以应用第一种模式（视窗采点模式），直接以数字地图、或经过校正的图像、或注记图层作为地理参考，在另一种视窗中打开相应的数据层，从中采集控制点。如果事先经过通过GPS测量、或摄影测量、或其它途径获得了控制点的坐

31、标数据，并保存为ERDAS的控制点文件格式或ASCII数据文件的话，就应该调用第二种类型（文件采点模式），直接在数据文件中读取控制点坐标。如果前两种条件都不符合，只有硬拷贝的地图或坐标纸作为参考的话，则只好采用第三种类型（地图采点模式），要么首先在地图上选点并量算坐标，然后通过键盘输入坐标数据；要么在地图上选定后，借助数字化仪来采集控制点坐标。在实际工作中，这三种采点模式都是有可能遇到的。2.资源卫星图像校正（Rectify Landsat Image）图像校正的一般流程（General Rectification Workflow）下面介绍基于SPOT图像和TM图像校正过程，其工作流程如图2

32、.15所示：图2.15 图像校正的一般流程图像校正的具体过程（Concrete Rectification Process）第一步：显示图像文件（Display Image Files）首先，在ERDAS图 面板中点击Viewer，打开两个视窗，并将两个视窗平铺放置，操作过程如下：REDAS图标面板菜单条：Session Tile Viewers然后，在Viewer#1中打开需要校正的Landsat TM图像：tmAtlanta.img在Viewer#2中打开作为地理参考的校正过的SPOT图像：panAtlanta.img第二步：启动几何校正模块（Geometric Correction To

33、ol）Viewer#1菜单条：Raster Geometric Correction打开Set Geometric Model对话框（图2.13）选择多项式几何校正计算模型PolynomialOK同时打开Geo Correction Tools对话框（图2.16）和Polynomial Model Properties对话框（图2.17）图2.16 Geo Correction Tools对话框在Polynomial Model Properties对话框中，定义多项式模型参数及投影参数：定义多项式次方（Polynomial Order）：2定义投影参数（Projection）：Apply C

34、lose打开GCP Tool Reference Setup对话框图2.17 Polynomial Model Properties对话框说明：该实例是采用视窗采点模式，作为地理参考的SPOT图像已经含有投影信息，所以，这里不需要定义投影参数。如果不是采用视窗采点模式，或者参考图像没有包含投影信息，则必须在这里定义投影信息，包括投影类型及其对应的投影参数。第三步：启动控制点工具（Start GCP Tools）首先在GCP Tools Reference Setup对话框（图2.12）中选择采点模式：选择视窗采点模式：Exiting ViewerOK（关闭GCP Tools Reference

35、 Setup对话框）打开Viewer Selection Intruction指示器（图2.18）图2.18 Viewer Selection Instruction指示器在显示作为地理参考图像panAtlanta.img的Viewer#2中点击左键打开Reference Map Information提示框（图2.19）（显示参考图像的投影信息）OK（关闭Reference Map Information提示框）整个屏幕将自动变化为如图2.20所示的状态：其中包含两个主视窗，两个放大窗口、两个关联方框（分别用于两个视窗中，指示放大窗口与主视窗的关系）、控制点工具对话框、几何校正工具等。表明控

36、制点工具被启动，进入控制点采集状态。图2.19 Reference Map Information提示框第四步：采集地面控制点（Ground Control Point）控制点工具对话框简介（Introduction to GCP Tool）在正式开始采集控制点之前，首先对控制点工具对话框进行说明：GCP对话框（GCP Tool）由菜单条、工具条和控制点数据表、状态条四个部分组成，菜单条中命令及其功能如表2.3所列，工具条中的图标及其功能如表2.2所列，控制点坐标表的组成如表2.5所列。图2.20 Reference Map Information提示框表2.3 GCP菜单命令及其功能命 令功

37、 能File：Load InputSave InputSave Input AsLoad ReferenceSave ReferenceSave Reference AsClose文件操作：调用输入控制点文件（*.gcc）将输入控制点保存在图像中保存输入控制点文件（*.gcc）调用参考控制点文件（*.gcc）将参考控制点保存在图像中保存参考控制点文件（*.gcc）关闭控制点工具Viewer：View Only Selected GCPsShow Selected GCP in TableArrange Frames on ScreenToolsStart Chip Viewer显示操作：视窗仅

38、显示所选择的控制点在表格显示所选择的控制点重新排列屏幕中的组成要素调出控制点工具图标面板重新打开放大窗口Edit：Set Point Type（Control/Check）Reset Reference SourcePoint PredictionPoint Matching编辑操作：设置采集点的类型（控制/检查）改变参考控制点源文件改变参考文件的投影参数匹配两幅图像像元的灰度值Help：Help for GCP Tool联机帮助关于GCP工具的联机帮助表2.4 GCP工具图标及其功能图标命令功能Toggle自动GCP编辑模式开关键Calculate依据控制点求解几何校正模型Automatic

39、设置自动转换计算开关Compute Error计算检验点的误差，更新RMS误差Select GCP激活GCP选择工具，在视窗中选择GCPCreate GCP在视窗中选择定义GCPLock/Unlock锁住/释放当前命令，以便重复使用Find in Input选择寻找输入图像中的GCPFind in Refer选择寻找参考文件中的GCPZ Value计算更新所选的GCP的Z值Auto-Z Value自动更新所有的GCP的Z值表2.5 GCP数据表字段及其含义字 段含 义Point #GCP顺序号，系统自动产生Point IDGCP标识码，用户可以定义GCP当前选择状态提示符号Color输入GCP

40、显示颜色X Input输入GCP的X坐标Y Input输入GCP的Y坐标Color参考GCP显示颜色X Reference参考GCP的X坐标Y Reference参考GCP的Y坐标TypeGCP的类型（控制点/检查点）X Residual单个GCP的X残差Y Residual单个GCP的Y残差RMS Error单个GCP的RMS残差Contribution单个GCP的贡献率Match两幅图像像元灰度值的匹配程度关于GCP工具对话框，还需要说明几点：（a）输入控制点（Input GCP）的是在原始文件视窗中采集的，具有原文件的坐标系统；而参考控制点（Reference GCP）是在参考文件视窗中

41、采集的，具有已知的参考坐标系统，GCP工具将根据对应点的坐标值自动转换模型。（b）在GCP数据表中，残差（Residual）、中误差（RMS）、贡献率（Contribution）及匹配程度（Match）等参数，是在编辑GCP的过程中自动计算更新的，用户是不可以任意改变的，但可以通过精确GCP位置来调整。（c）每个IMG文件都可以有一个GCP数据集与之相关联，GCP数据保存在一个栅格层数据文件中；如果IMG文件有一个GCP数据集存在的话，只要打开GCP工具，GCP点就出现在视窗中。（d）所有的输入GCP都可以直接保存在图像文件中（Save Input），也可以保存在控制点文件中（Save Inp

42、ut As）。如果是保存在文件中，调用的方法如（c）所述，如果是保存在GCP文件中，可以通过加载调用（Load Input）。（e）参考GCP也可以类似的保存在参考图像中（Save Reference）或GCP文件中（Save Reference As），以便以后调用。GCP的具体采集过程（Steps of GCP Selection）在图像几何校正过程中，采集控制点是一项非常重要和相当繁重的工作，具体过程如下：a.在GCP工具对话框中点击Select GCP图标，进入GCP选择状态；b.在GCP数据表中将输入GCP的颜色（Color）设置为比较明显的黄色；c.在Viewer#1中移动关联方框

43、位置，寻找明显的地物特征点，作为输入GCP；d.在GCP工具对话框中点击Create GCP图标，并在Viewer#3中点击左键定点，GCP数据表将记录一个输入GCP，包括其编号、标识码、X坐标、Y坐标；e.在GCP工具对话框中点击Select GCP图标，重新进入GCP选择状态；f.在GCP数据表中将参考GCP的颜色设置为比较明显的红色；g. 在Viewer#2中移动关联方框位置，寻找对应的地物特征点，作为参考GCP；h. 在GCP工具对话框中点击Create GCP图标，并在Viewer#2中点击左键定点，系统将自动把参考点的坐标（X Reference、Y Reference）显示在GC

44、P数据表中；i.在GCP工具对话框中点击Select GCP图标，重新进入GCP选择状态，并将光标移回到Viewer #1中，准备采集另一个控制点；j.不断重复上述步骤（ai），采集若干GCP，知道满足所选定的几何校正模型为止；而后，每采集一个Input GCP，系统就自动产生一个Ref.GCP，通过移动Ref.GCP可以逐步优化校正模型，采集GCP后，GCP数据表如图2.21所示。图2.21 GCP工具对话框第五步：采集地面检查点（Ground Check Point）以上所采集的GCP的类型均为Control（控制点），用于控制计算，建立转换模型及多项式方程。下面所要采集的GCP的类型均为

45、Check（检查点），用于检验所建立的转换方程的精度和实用性。依然在GCP工具对话框状态下：a. 在GCP工具对话框菜单条中确定GCP类型：Edit Set Point Type Checkb. 在GCP Tool菜单条中确定GCP匹配参数（Matching Parameter）Edit Point Matching 打开GCP Matching对话框GCP Matching对话框在GCP Matching对话框中，需要定义下列参数：匹配参数（Matching Parameter）：最大搜索半径（Max.Search Radius）：3搜索窗口大小（Search Window Size）：X:

46、5 Y:5约束参数（Threshod Parameter）相关阈值（Correlation Threshold）：0.8删除不匹配的点（Discard Unmatched Point）：Active匹配所有/选择点（Match All/Selected Point）：从输入到参考（Reference from Input）从参考到输入（Input from Reference）Close（关闭GCP Matching对话框）c. 确定地面检查点，在GCP工具条中选择Create GCP图标，并将Lock图标打开，锁住Create GCP功能，如同选择控制点一样，在Viewer#1和Viewer#2中定义5个检查点，定义完毕后，点击Unlock图标，解除Create GCP功能d. 计算检查点误差，在GCP工具条中点击Compute Error图标，检查点的误差就会显示在GCP Tool的上方，只有所有检查点的误差均小于一个像元，才能继续进行合理的重采样。一般来说，如果你的控制点定位选择比较准确，检查点匹配会比较好，误差会在限差范围内；否则，若控制点定义不准确，检查点就无法匹配，误差会超标。第六步：计算转换模型（Compute Transformation）在控制点采集过程中，一般是设置为自动转换计算模式，所