实验二、数据预处理

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

2、围调整、误差校正、坐标转换等处理，以便进一步开展图像解译、专题分类等分析研究。数据预处理模块简称Data Prep aration或Data Prep,可以通过两种途径启动:ERDAS 图标面板菜单条：Main Data PreparationData Preparation 菜单（图 2.1）ERDAS图标面板工具条：点击Data Prep图标Data Prep aratio n菜单(图2.1)|呢血5 FrQrCreate New Image,.1|匚Surface.|Subset Image.IImage Geometric Correction.1Mosaic Images.|Uhsu

3、pervised 匚bssif»calion.|Reprojed Images.I1Close1Help|图 2.1 Data Preparation 菜单从图2.1可以看出，ERDAS IMAGIMG 数据预处理模块包括了 7项主要功能，其中第一项功能(生成单值图像)比较简单，第六项功能(非监督分类)将在图像分类中进行说明。F面将主要介绍其余五项功能，重点是图像几何校正和图像拼接处理,因为这两项操作是从事遥感应用研究必须开展的基本工作过程。二、三维地形表面(3D Surfacing)三维地形表面工具允许用户在不规则空间点的基础上产生三维地形表面,所支持的输入数据类型包括：ASCII

4、码点文件、Arcinfo的Coverage点文件和线文件,ERDAS IMAGING的注记数据层，以及栅格图像文件IMG。所有输入数据必须具有 X、Y、Z值，三维地形表面工具所应用的TIN插值方法，所输出的是一个连续的栅格图像文件。每一个已知的空间点在输出的地形表面上保持Z值不变,而没有Z值的空间点，其输出表面的Z值是基于TIN其周围的已知点插值计算获得的。在三维地形表面工具中提供了两种TIN插值方法：线性插值(Linera)与非线性插值(non-li near)。线性插值方法是应用一次多项式方程进行计算,输出的TIN三角面是一些有输出的是平滑表面，这种情况棱角的平面；非线性插值方法应用五次多

5、项式方程进行计算, 下，TIN三角面不是一个平面，而是具有弹性的曲面。线性插值方法速度快但结果简单，而 非线性插值方法产生基于不规则分布数据集的非常连续的、圆滑的表面结果。1.启动三维地形表面(Surfaci ng Create Surface)ERDAS 图标面板菜单条：Main 7 Data Preparation菜单Create Surface打开 3D Surfacing 对话框(图 2.2)7选择ERDAS图标面板工具条：点击3D Surfacing对话框由菜单条（Data Pre 图标，打开 Data Preparation 菜单Create Surface打开 3D Surfac

6、ing 对话框(图 2.2)Menu Bar ）、工具条（Tool Bar）和数据表格（DataCellArray ）组成。菜单条主要由文件操作 （File ）、数据表操作（Row）和表面生成（Surface） 菜单组成，而工具条则由读取数据（ Read Point）文件、保存数据文件（Save Point）和生成地形表面（Perform Surfacing ）图标组成。|说引 SurfacingFile Ko* Surface Help9 H ©Rm Use1XY1z1H1可5CO0G0.1252.0002X73457987549?3e6.E313.00034X723039.813

7、492246.1834.000X721715 2504S2331.7815.0005X731579 685438306.6566.000e7X729694.433498239.0947.000X729053.0004SG529.3448.0000X731671 625437557.5319.000910 =1X7329716.63S493978.00010.000X733253G304EG407.O9411.0004图2.2 3D Surfacing对话框（读入数据之后）2.定义地形表面参数(Defi ning Surface Parameters)Source Fib Type:| ASCI

8、I FileSource File Name: (" dat)OKCancel |图2.3 Read Point对话框在Read Point对话框中需要定义下列参数:睡二1117数据源文件类型(Source File Type) : ASCII File7 数据源文件名称(Source File Name): Inpts.dat7OK （关闭 Read Points对话框）7打开 import Options 对话框（图 2.4） 在Import Options栏目，需要定义下列参数：7选择字段类型（Field Type) : Delimited by Separation (分割字

9、符)7选择分割字符（Separation Charactor) : Comma (逗号分割)匸因d:/program files/innagine 3.5/ewamples/lnpts.(jatFile to Import From:FieldDeMon | Input ReviewField Type: 泾 S鱼血?妝邑剧2®Fixed WidthSeparator Character:I CommaRow Terminator Character:IR eturn N ewLine DOS)Comment Character:Number of Rows T Skip:OutDU

10、t Column NameInput Field NumberyX11Y2Z3斗11Column MappingCarcel IOKViewHelp图 2.4 Import Options 对话框(Field Definition 栏目)7 每行结束字符(Row Terminator Charactor): Return New Line ( DOS)7确定跳过几行（Number of Rows to Skip ）: 0 （从头读）7点击 Input Preview 标签，进入 Input Preview 栏目（图 2.5）从Import Options对话框显示的原始数据可知，数据文件中的数

11、据记录方式是一行一个点，每一行数据包括点号、X坐标、丫坐标、Z坐标（高程值）四个字段，其中点号在此处读入数据时不需要，因此，必须在Import Options对话框的Colum Mapping中确定X、Y、Z与数据文件中字段的对应关系:Out put Colu mn Name:X对应Input Field Number : 27 Out put Colu mn Name:Y对应Input Field Number : 37 Out put Colu mn Name:Z对应Input Field Number : 47OK （关闭Import Options对话框，读入数据）数据读到 3D Su

12、rfacing Data CellArray 中（图 2.2）Import OptionsFile to Import From:d:/program files/innagine 3.5/ewamples/lnpts.（jatFiekJDdinition Input Previei/JRow1Field 1Field 21gJ47591.3135,3-22,247591.5312533,247EG0.G09375,344,247670,2187555,247606.953125£elG,247834.953125,37411J47716 09375OutDUt Column Nam

13、eInput Field NumberLIX21Y3Z4J1Column MappingOKVi訓CancelHelp图 2.5 Import Options 对话框（Input Preview 栏目）如果需要，可以将读入的数据保存为ERDAS的Annotation Layer或Arcinfo的Coverage文件:3D Surfacing dialog 菜单：File 7 Save as 7 Save as对话框（图 2.6）7确定输出文件类型（Out put File Type ） : P oi nt Coverage确定输出文件名（Save as Arcinfo） : test poin

14、tOutput File Type:I Point CoverageAttribute for Z:ELEVATIONPrecisionI SingleOK口 ncel图2.6 Save as对话框3.生成三维地形表面( Creat ing 3D Surface)3D Surfacing 对话框菜单条：Surface 7 Surfacing 7打开 Surfacing 对话框(图 2.7)3D Surfacing对话框工具条：点击 Perform Surfacing图标7 打开Surfacing对话框(图2.7)ifl Sy facingSurfacing Method：tectsurface

15、.imgI Linear Rubber SheelirgColumns:511Rows:512Save Output Coords A$： (".arcHo)呂琳 poiritlBackground ValueOKBalchCarcell |HeIp I1233100 000000弓 lRX 124a405Le59375弓13307079 000000甘 lRY: 13791755 000000甘Output Cell Sizes;130.000000弓 Y： 130.000000甘R Square CellsOutput Corrers:ULX：ULY:X0.0000尸oinoin

16、oinoinoiiooiioiMoiioinoinoin“ro-5H ilgriore Zero b Output StatsiOutput Dta T>ipe:jUnsigned IE bit图2.7 Surfacing对话框在Surfacing对话框中需要设置下列参数:7输出文件名称(Output File ）： testsurface.img7 表面插值方法(Surfacing Method ）： Linear Rubber Sheeting7 输出文件范围(Output Corners ）：自动读取 ULX 、ULY、LRX 、LRY7 输出像元大小(Output Cell Si

17、ze ）：X：30/Y：307 输出像元形状：Square Cells7 输出图像背景值Background Value）： 07 输出统计忽略值Ignore Zero in Output Stats ）： 07 输出数据类型(Output Data Type）： Unsigned 16 bit7 OK (关闭 Surfacing 对话框，执行地形表面过程)2.显示三维地形表面( Display 3D Surface )由三维地形表面过程生成的图像文件， 就是 ERDAS 系统的 DEM 文件， DEM 文件常常用于三维图像( Image Drape )的显示或虚拟地理信息系统( Virtua

18、l GIS )操作。所以，可以在 Viewer 视窗中打开曲面所生成的地形表面文件 testsurface.img 文件，显示其二维平面效果或通过 Image Info 查看其定量信息，也可以在其上叠加对应的具有相同投影系统的图像文件，显示其三维立体效果。三、图像分幅裁剪( Subset Image)在实际工作中， 经常需要根据研究工作范围对图像进行分幅裁剪，按照 ERDAS 实现图像分幅裁剪的过程，可以将图像分幅裁剪分为两种类型：规则分幅和不规则分幅。1.规则分幅裁剪( Rectangle Subset Image) 规则分幅是指裁剪图像的边界范围是一个矩形，通过左上角和右下角两点的坐标，

19、就可以确定图像的裁剪位置，整个裁剪过程比较简单。ERDAS 图标面板菜单条： Main 7 Data Preparation 7 Data Preparation 菜单（图 2.1）7选择 Subset Image7 开 Subset Image 对话框（图 2.8）ERDAS图标面板工具条：点击Data Prep图标7打开Data Preparation菜单(图2.1)7 选择 Subset Image7 打开 Subset Image 对话框（图 2.8）iSSubsetIriput Fie： (".imgLlpanierimgOutput File： ('.img)la

20、nier_sub. imgCoordinate Type:Subset Definition:From Inquire BoxC MapLILX:0.00LR>i511 00a FileULY:0.00LR Y:511.00Data Type:Input： Unsigned 3 bitOutput: |Unsigred 3 bitOutput:ConrinuousNumter of Input layers:Output 0ptions:y ignore Zero in Output Slab.!17Select La卿說Uwe a comma for wepa同ed listi.e.

21、1,3,5 J or enter ranges u$inga''："(i-e. 2:5)KB吕tchI AOI .Cancel | Help |图 2.8 Subset Image 对话框在Subset Image对话框中需要设置下列参数:7输入文件名称(Input File ): Lanier.img7 输出文件名称(Out put File): Lan ier_sub.img7 坐标类型(Coordinate Type): File7 裁剪范围(Subset Definition )：输入 ULX、ULY、LRX、LRY7 输出数据类型(Out put Data

22、Type) : Un sig ned 8 bit7 输出统计忽略零值：ign ore Zero in Out put Stats输出像元波段（Select Layers）: 1:7 （表示选择从第一波段到第七波段）2,3,2 （表示选择2, 3, 2四个波段）7OK （关闭Subset Image对话框，执行图像裁剪）说明：在上述图像裁剪过程中，裁剪范围是通过直接输入左上角坐标和右下角坐标定In quire Box),然后义的。此外，还可以通过两种方式定义裁剪范围：其一是应用查询框（ 在Subset Image对话框中选择 From In quire Box功能；其二是应用 AOI ,然后在S

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

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

25、 Image对话框中需要设置下列参数:7输入文件名称(Input File ): Lanier.img7输出文件名称（Out put File ）:Lani er_sub.img7应用AOI确定裁剪范围：点击AOI按钮7 打开选择 AOI (Choose AOI)对话框（图2.9）NoneC Viewer_ & 01 File：orCancel I Help ISelect he AOI File： F am subset, aoi图2.9 Choose AOI对话框7在 Choose AOI 对话框中确定 AOI 的来源(AOI Source ): File 或 Viewer7 输出

26、数据类型(Out put Data Type) : Un sig ned 8 bit7 输出统计忽略零值：ign ore Zero in Out put Stats7输出像元波段（Select Layers）: 1:7 （表示选择从第一波段到第七波段）2,3,2 （表示选择2, 3, 2四个波段）OK （关闭Subset Image对话框，执行图像裁剪）Arcinfo 多边形裁剪（Polygon Coverage Subset Image）ArcI nfo如果按照行政区划边界或自然区划边界进行图像的分幅裁剪，往往是首先利用Polygon），然后以 ArcInfo 的 Polygon 为或ERD

27、AS的Vector模块绘制精确的边界多边形（边界条件进行图像裁剪。对于这种情况，需要调用ERDAS其他模块的操作分两步完成。第一步:将Arcinfo多边形转换成栅格图像文件ERDAS图标面板菜单条:Main 7 image Interpreter 7 Utilities 7 Vector to Raster7打开 Vector to Raster 对话框（图 2.10）ERDAS图标面板工具条:Mai n 宀 Vector 宀 Vector to Raster7打开 Vector to Raster 对话框（图 2.10）ERDAS图标面板工具条:'点击Interpreter 图标 7

28、 Utilities 7 Vector to Raster7打开 Vector to Raster 对话框（图 2.10）ERDAS图标面板工具条:点击 7 IVectot 图标 Vector to Raster7打开 Vector to Raster 对话框（图 2.10）f£i. Vector to BesterOutput Raster: （".imgzoneesT extures1 .imagineSSOergxlib1 Application Datagnutimelib=l Cookiesinfo1 Favoritesorthobase_j Local Sett

29、ingsvirtualgis二1 Mv DocumentszoneElS 1 NetHoodT |nj administratorexam ple$raster ingOKCancelJ -llHelpmmi-ni-i-, .isi图 2.10 Vector to Raster 对话框（a）£0 Convert Vector Layer to Raster LayerInput Coverage Type： 金 PolygonLine O PointWeight T able File:DaU Compression:Output Innage Type:OdadtAm lltefn

30、 as Pinell Value:Thematic Continuous 臣 Create FramkJ Layers|Z0NE3S-ID回ge全!2虫,© 乂 無虜(7OKCancelBatchHmlp ISubwet Definilwn in Mop Cowd's):ULX: 1709335.69ULY: I 501 S72.31:dLRX: 73995E.31d :dLRY:473099 59Output Options：Riwel Size | 30.00000Lookup Table File:Block Size:| 4Output Pixel Data Type

31、:I Unsigned 3ZbitH图 2.10 Vector to Raster 对话框(b)在Vector to Raster对话框中需要设置下列参数:7确定矢量文件名称(Input Vector File ): zone887 确定输出文件名称( Out put Raster File): raster.img7确定矢量文件类型(Vector Type): Polygon7转换范围大小(Subset Definition )： ULX、ULY、LRX、LRY7输出像元大小(Pixel Size ): 30栅格数据类型(Out put Pixel Data Type ): Un sig n

32、ed 8 bit7使用矢量属性值(An Item as Pixel Value ) : ZONE88-ID7 栅格文件类型(Out put Image Type) : Thematic7是否创建图层金字塔( Create Pyramid Layers)输出统计忽略零值：Ignore Zero in Out put StatsOK (关闭Vector to Raster对话框，执行失栅转换)第二步：通过掩模运算(Mask)实现图像不规则裁剪ERDAS 图标面板菜单条：Main 7 Image Interpreter 7 Utilities 7 Mask7打开Mask对话框(图2.11)ERDAS

33、 图标面板菜单条： Ma in 7 Data Prep aratio n7 Image Geometric Correction7打开Mask对话框(图2.11)% MaskInput File: img)I brier imgWindow;L UnJon 冷 lntef$ectioriOutput Fib: p.imgI mask.imgInput Mask File: (".img)I raster ing0 ignore2eio in Output Stdd nininiBinniBii-iiBivM -i-i'-i-i' n- rZero's Indi

34、cate encluded Area.D 旳 T ype：Input til:Unsigned S bitllinput W2:Unsigned S bitOutput：1 Unsigned 8 bit二|Setup Recode AOI .BatchVieiw.Help图2.11 Mask对话框在Mask对话框中需要设置下列参数:7 确定输入文件名称(Input Vector File ): Lanier.img7确定掩模文件名称(Input Mask File ): raster.img7点击Set Recode设置裁剪区域内新值(New Value)为1，区域外取0值7确定掩模区域做交集

35、运算：In tersection7输出图像文件名称(Out put File): mask.img7 输出数据类型(Out put Data Type) : Un sig ned 8 bit7 Ok (关闭Mask对话框，执行掩模运算)7实现图像不规则裁剪。四、图像几何校正(Geometric Correction)几何校正就是将图像数据投影到平面上，使其符合(Conform)地图投影系统的过程,而将地图坐标系统赋予图像数据的过程，称为地理参考(Geo-referencing )。由于所有地图投 影系统都遵从于一定的地图坐标系统，所以几何校正过程包含了地理参考过程。1. 图像几何校正概述(In

36、troduction to Geometric Correction )在正式开始介绍图像几何校正分幅和过程之前,首先对ERDAS图像几何校正过程中的几个普遍性的问题进行简单说明，以便于随后的操作。图像几何校正途径(Geometric Correction Process )数据预处理途径(Start from Data Preparation )7打开 Set Geo-Correction Input File 对话框（图 2.12）ERDAS 图标面板工具条：点击Data Prep 图标 7 Image Geometric Correction打开 Set Geo-Correction I

37、nput File 对话框（图 2.12）監 S电t G电Q Corfsetion Input FilsFrom VieuverFrom Image FifeCancel IHelp I图 2.12 Set Geo-Correction Input File 对话框在Set Geo-Correction Input File对话框中，需要确定校正图像，有两种选择情况其一：首先确定来自视窗(From Viewer)，然后选择显示图像视窗(Select Viewer)7打开 Set Geometric Model 对话框（图 2.13）7选择几何校正计算模型(Select Geometric Mo

38、del ) 7 OK7打开校正模型参数与投影参数设置对话框7定义校正模型参数与投影参数7Apply 7 Close凶丽ne1CameraLanchatPolvriomialReprojectRjbter SheetingSpotOpen Eh曲ifig Model.'匚伽詐湎罰帀1_M 1 Cancel 1Help |图 2.13 Set Geometric Model 对话框7打开 GCP Tool Referenee Setup 对话框（图 2.14）7确定米点模式，米点校正IHGCF Tool Reference S&tupExisting ViewerImage Lay

39、er New Viewer)Vector Layer (N ew ViewerArinotdtian Layer (N e鞫 VieifverGCP File .gccASCIIII FileDigitizing Tablet (Current ConfigurarionDigitizing T訪l&t (New ConfigurationKeyboard Only厂“EgCancelHelp图 2.14 GCP Tool Refere nee Set up 对话框其二：首先确定来自文件(From Image File )然后选择输入图像(Input Image File )打开 Se

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

41、rGeometric Correction7选择几何校正计算模型(Select Geometric Model ) 7 OK7打开校正模型参数与投影参数设置对话框7定义校正模型参数与投影参数7Apply 7 Close7打开 GCP Tool Referenee Setup 对话框(图 2.14)7确定米点模式，米点校正几何校正模块启动的两种途径相比，视窗栅格操作途径更为直观简便。2. 几何校正计算模型(Geometric Correction Model )ERDAS提供的图像几何校正计算模型有7种，具体功能如表 2.1所列:其中，多项式变换在卫星图像校正过程中应用较多,在调用多项式模型时，

42、需要确定多表2.1几何校正计算模型与功能模型功能Affine图像仿射变换（不做投影变换）P olyomial多项式变换（冋时做投影变换）Rep roject投影变换（转换调用多项式变换）Rubber Sheet ing非线性、非均匀变换Camera航空影像正射校正LandsatLa ndsat卫星图像正射校正SPOTSPOT卫星图像正射校正项式的次方数（Order），通常整景图像选择 3次方，次方数与所需要的最少控制点数是相关的，最少控制点计算公式为（t+1）（t+2）/2,式中t为次方数，即一次方程需要 3个控制点,二次方程需要6个控制点，三次方程需要 10个控制点，依此类推。3. 几何校正

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

44、:GCP File (*.gcc)ASCII File文件采点模式：在控制点文件中读点ASCII码文件中读点Map to Viewer :Digitiz ing Tablet (Curre nt)Digitizing Tablet ( New)Keyboard Only地图采点模式：在当前数字化仪上采点在新配置数字化仪上采点通过键盘输入控制点表2.2所列的三种几何校正采点模式，分别应用于不同的情况:如果已经拥有需要校正图像区域的数字地图，或经过校正的图像， 或注记图层的话， 就可以应用第一种模式（视窗采点模式），直接以数字地图、或经过校正的图像、或注记图层 作为地理参考，在另一种视窗中打开相应

45、的数据层，从中采集控制点。如果事先经过通过 GPS测量、或摄影测量、或其它途径获得了控制点的坐标数据，并保存为ERDAS的控制点文件格式或 ASCII数据文件的话，就应该调用第二种类型 （文件采 点模式），直接在数据文件中读取控制点坐标。如果前两种条件都不符合，只有硬拷贝的地图或坐标纸作为参考的话，则只好采用第三种类型（地图采点模式），要么首先在地图上选点并量算坐标，然后通过键盘输入坐标数据;要么在地图上选定后，借助数字化仪来采集控制点坐标。在实际工作中，这三种采点模式都是有可能遇到的。2.资源卫星图像校正（Rectify Landsat Image）图像校正的一般流程（General Rec

46、tification Workflow ）图2.15 图像校正的一般流程F面介绍基于SPOT图像和TM图像校正过程，其工作流程如图2.15所示:图像校正的具体过程(Concrete Rectification Process )第一步：显示图像文件（Dis play Image Files）首先，在ERDAS图面板中点击Viewer，打开两个视窗，并将两个视窗平铺放置，操作过程如下:REDAS图标面板菜单条:Sessi on 宀 Tile Viewers然后，在Viewer#1中打开需要校正的 La ndsat TM图像：tmAtlanta.img在 Viewer#2中打开作为地理参考的校正过

47、的SPOT 图像：panAtlanta.img第二步：启动几何校正模块（ Geometric Correction ToolViewer#1 菜单条： Raster 7 Geometric Correcti on2.13)7打开 Set Geometric Model 对话框(图7选择多项式几何校正计算模型Polynomial 7 OK7同时打开 Geo Correction Tools对话框（图2.16）和Poly no mial Model Prop erties 对话框(图 2.17)Exit匚arractigHelp图 2.16 Geo Correction Tools 对话框在Pol

48、ynomial Model Properties对话框中，定义多项式模型参数及投影参数:7定义多项式次方(Polynomial Order )： 27定义投影参数(Projection): 7 Apply 7 Closefff_ Fclynofhi al Ho 触 1 Proper tits （Ho File）Parameter | Transformation | Projection |Load LFF FileStatus:Model has no solution.图 2.17 Polynomial Model Properties 对话框说明：该实例是采用视窗采点模式，作为地理参考的

49、SPOT图像已经含有投影信息，所 以，这里不需要定义投影参数。如果不是采用视窗采点模式，或者参考图像没有包含投影信 息，则必须在这里定义投影信息，包括投影类型及其对应的投影参数。第三步：启动控制点工具(Start GCP Tools)首先在GCP Tools Referenee Setup对话框(图2.12)中选择采点模式:7选择视窗采点模式：Exiti ng Viewer7OK （关闭 GCP Tools Referenee Setup 对话框）7打开 Viewer Seleetion Intruetion 指示器（图 2.18）Vie*er Selection Instruct!onsfo

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

51、，进入控制点采集状态。Current Reference Map Projection：Pejection： State PlaneSpheroid:Clarke 1866Zone Number:367GDatum: NAO27UI-'Cancel Ii DKUmmwouf 121豳 y UffeHelp I图 2.19 Referenee Map Information 提示框第四步：采集地面控制点（Ground Control Point)控制点工具对话框简介（Introduction to GCP Tool)在正式开始采集控制点之前,首先对控制点工具对话框进行说明:GCP 对话框

52、（GCP Tool）由菜单条、工具条和控制点数据表、状态条四个部分组成,菜单条中命令及其功能如表2.3所列，工具条中的图标及其功能如表2.2所列，控制点坐标表的组成如表2.5所列。a落邈7 £ 迸I竺注-Inlxilb Vin融 03, 136905.23Ifa叭灯tCtate riaiie / Clarke 16bb) 扌f 凶 o + %陀 CCF 7 ooL : CCn-poj-t '如 dIIuiI, i 百)Qldftrbna*puiLtluit*. Lng)題 2塹緜lx©直糯册 2器r upe 311 D137.000|”9(7.00Ch 111 Conbol11 Conhcl图 2.20 Referenee Map Information 提示框表2.3 GCP菜单命令及其功能File :命令功 能Load InputSave In