遥感数据图像处理实验三、遥感图像的几何校正与裁剪

实验三、遥感图像的几何校正与裁剪实验内容：1.图像分幅裁剪（SubsetImage）2.图像几何校正（GeometricCorrection）3.图像拼接处理（MosaicImgaes）4.生成三维地形表面（3DSurfacing）1.图像分幅裁剪在实际工作中，经常需要根据研究工作范围对图像进行分幅裁剪，按照ERDASIMAGINE8.4实现图像分幅裁剪的过程，可以将图像分幅裁剪为两类型：规则分幅裁剪，不规则分幅裁剪。1.1规则分幅裁剪(以c:\ProgramFile\IMAGINE8.4\examples\lanier.img为例)规则分幅裁剪是指裁剪图像的范围是一个矩形，通过左上角和右上角两点的坐标可以确定图像的裁剪位置，过程如下：方法一:→ERDASIMAGINE8.4图标面板菜单条：Main→DataPreparation(或单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标)→打开DataPreparation对话框→单击SubsetImage按钮，打开Subset对话框在Subset对话框中需要设置下列参数：→输入文件名（InputFile）:lanier.img→输出文件名（OutputFile）:lanier_sub.img→坐标类型（CoordinateType）:Map→裁剪范围（SubsetDefinition）:ULX、ULY、LRX、LRY（注：ULX，ULY是指左上角的坐标，LRX,LRY是指右上角的坐标，缺省状态为整个图像范围）→输出数据类型（OutputDataType）:Unsigned8Bit→输出文件类型（OutputLayerType）:Continuous→输出统计忽略零值:IgnoreZeroInOutputStats→输出像元波段（SelectLayers）:2,3,4→OK(关闭Subset对话框，执行图像裁剪)方法二:→ERDASIMAGINE8.4图标面板菜单条:Main→StartIMAGINEViewer(或单击RDASIMAGINE8.4图标面板工具条“Viewer”图标)→打开一个二维视窗→单击视窗工具条最左端的“打开文件”图标→打开SelectLayerToAdd对话框在SelectLayerToAdd对话框完成以下设置：→LookIn:examples→FileName:lanier.img→Filesoftype:IMAGINEImage→双击OK按钮→在二维视窗中打开lanier.img文件→单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标→打开DataPreparation对话框→单击SubsetImage按钮→打开Subset对话框在Subset对话框中需要设置下列参数：→输入文件名（InputFile）:lanier.img→输出文件名（OutputFile）:lanier_sub.img→坐标类型（CoordinateType）:Map→输出数据类型（OutputDataType）:Unsigned8Bit→输出文件类型（OutputLayerType）:Continuous→输出统计忽略零值:IgnoreZeroInOutputStats→输出像元波段（SelectLayers）:2,3,4→单击FromInquireBox按钮→打开InvalidCoordinateType对话框→单击Continue→在显示图像文件lanier.img视窗中单击工具条的“+”按钮,打开InquireCursor对话框,在视窗中移动十字光标,确定裁剪范围左上角和右下角,读取其坐标分别填入SubsetImage对话框的ULX,ULY中和LRX,LRY中→单击OK按钮(关闭Subset对话框，执行图像裁剪)方法三:首先在视窗中打开lanier.img文件→AOI→Tools打开AOI工具面板→单击矩形框确定裁剪范围→File→Save→AOILayerAs→打开SaveAOIAs对话框,输入文件名:2→单击OK(退出SaveAOIAs对话框)→单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标→打开DataPreparation对话框→单击SubsetImage按钮→打开Subset对话框在Subset对话框中需要设置下列参数：→输入文件名（InputFile）:lanier.img→输出文件名（OutputFile）:lanier_sub.img→坐标类型（CoordinateType）:Map→输出数据类型（OutputDataType）:Unsigned8Bit→输出文件类型（OutputLayerType）:Continuous→输出统计忽略零值:IgnoreZeroInOutputStats→输出像元波段（SelectLayers）:2,3,4→单击AOI按钮→打开ChooseAOI对话框→在ChooseAOI对话框作如下设置：→AOISource:File→AOIFile:2→单击OK(退出ChooseAOI对话框)→单击OK(退出Subset对话框,执行图像裁剪)→单击OK(退出Modeler对话框,完成图像裁剪)1.2不规则分幅裁剪不规则分幅裁剪是指裁剪图像的边界范围是个任意多边形，无法通过左上角和右下角两点的坐标确定图像的裁剪位置，而必须事先生成一个完整的闭合多边形区域，可以是一个AOI多边形，也可以是ArcInfo的一个PolygonCoverage，针对不同的情况采用不同的裁剪过程。（一）AOI多边形裁剪(以c:\ProgramFile\IMAGINE8.4\examples\lanier.img为例)→在视窗中打开需要裁剪的图像→应用AOI工具绘制多边形可以将多边形AOI文件保存在文件中（＊.aoi）,也可以暂时不退出视窗，将图像与AOI多边形保留在视窗中，然后进行以下操作：→ERDASIMAGINE8.4图标面板菜单条：Main→DataPreparation(或单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标)，打开DataPreparation对话框→单击SubsetImage按钮，打开Subset对话框在Subset对话框中需要设置下列参数：→输入文件名（InputFile）:lanier.img→输出文件名（OutputFile）:lanier_sub.img→CoordinateType:Map→Output:Unsigned8bit→Oupput:Continuous→选定IgnoreZeroinOutputStats→SelectLayers:2,3,4→单击AOI按钮，打开ChooseAOI对话框→选定File→选择文件：2.aoi→单击OK(退出ChooseAOI对话框OK（退出Subset对话框→单击Modeler对话框中的OK按钮完成图像裁剪（二）ArcInfo多边形裁剪)→单击,执行图像裁剪）(以c:\ProgramFile\IMAGINE8.4\examples\gy_resample.img为例)如果是按照行政区划边界或自然区划边界进行图像的分幅裁剪，往往是首先利用ArcInfo或者ERDAS的Vector模块绘制精确的边界多边形（Polygon），然后以ArcInfo的Polygon为边界进行图像裁剪。对于这种情况，需要调用ERDAS其它模块的功能分两步完成。第一步：将ArcInfo多边形转换为栅格图像文件打开VectorToRaster对话框有两种方法：方法一：ERDASIMAGINE8.4图标面板菜单条：Main→ImageInterpreter(或单击ERDASIMAGINE8.4图标面板工具条“Interpreter”图标)→打开ImageInterpreter对话框→单击Utilities→打开Utilities对话框→选择VectorToRaster→打开VectorToRaster对话框方法二：ERDASIMAGINE8.4图标面板菜单条：Main→Vector(或单击ERDASIMAGINE8.4图标面板工具条“Vector”图标)→打开VectorUtilities对话框→选择VectortoRaster→打开VectortoRaster对话框本例中选择前一种方法：在VectortoRaster对话框中设置下列参数：→InputVectorFile（输入矢量文件名称）：boundgy.shp→VectorType(确定矢量文件类型)：polygon→UseAttributeasValue(使用矢量属性值)：GYBOUND_ID→OutputFile（输出栅格文件名称）：raster.img→DataType(栅格数据类型)：Unsigned8bit→LayerType(栅格文件类型)：Thematic→SizeDifinition(转换范围大小)：ULX,ULY,LRX，LRY（缺省条件下为整个图像范围）→Units（坐标单位）：Meters→CellSize（输出像元大小）：X:30/Y:30→选择SquareCell（正方形像元）→单击OK（关闭VectorToRaster对话框，执行矢量到栅格的转换）矢量图：栅格图：第二步：通过掩模运算（Mask）实现图像的不规则裁剪ERDASIMAGINE8.4图标面板菜单条：Main→ImageInterpreter(或单击ERDASIMAGINE8.4图标面板工具条“Interpreter”图标)→ImageInterpreter→Utilities→Mask→打开Mask对话框在Mask对话框中设置下列参数：→InputFile(输入需要裁剪的图像文件名称):gy_resample.img→InputMaskFile(输入掩模文件名称)：raster.img→点击SetupRecode按钮设置裁剪区域内新值（NewValue）为1，区域外取0值→选择Intersection（确定掩模区域作交集运算）→OutputFile（输出图像文件名称）：mask.img→Output(输出数据类型)：Unsigned8bit→单击OK（关闭Mask对话框，执行掩模运算）2.图像几何校正(GeometricCorrection)几何校正就是将图像数据投影到平面上，使其符合地图投影系统的过程；而将地图坐标系统赋予图像数据的过程，称为地理参考（Georeferencing）。由于所有地图投影系统都遵包含了地理参考过程。在ERDASIMAGINE8.4中进行图像几何校正,通常有两种途径启动几何校正模块：（1）数据预处理途径从于一定的地图坐标系统，所以几何校正过程(以c:\ProgramFile\IMAGINE8.4\examples\lanier.img为例)方法一：首先在视窗中打开lanier.img文件→ERDASIMAGINE8.4图标面板菜单条：Main→DataPreparation(或单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标)，打开DataPreparation对话框→选定ImageGeometricCorrection选项→打开SetGeoCorrectionInputFile对话框→选定FromViewer→单击SelectViewer按钮→打开ViewerSelectionInstructions对话框→在显示lanier.img图像的视窗中单击鼠标左键→打开SetGeometricModel对话框方法二：→ERDASIMAGINE8.4图标面板菜单条：Main→DataPreparation(或单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标)，打开DataPreparation对话框→选定ImageGeometricCorrection选项→打开SetGeoCorrectionInputFile对话框→选择FromImageFile→在InputImageFile中选择需要进行几何校正的图像：lanier.img→单击OK→同时打开一个二维视窗（lanier.img图像文件显示在其中）和SetGeometricModel对话框（2）窗栅格操作途径(以c:\ProgramFile\IMAGINE8.4\examples\lanier.img为例)首先在视窗中打开lanier.img文件视窗菜单条：Raster→GeometricCorrection→打开SetGeometricModel对话框资源卫星图像校正具体过程：①视窗采点模式：如果已经拥有需要校正图像区域的数字图像、或经过校正的图像、或注记图层，就可以应用视窗采点模式，直接以数字地图、或经过校正的图像、或注记图层作为地理参考，在视窗中打开该图层，从中采集控制点。（以c:\ProgramFile\IMAGINE8.4\examples下的panAtlanta.img和tmAtlanta.img为例，前者是作为地理参考的图像，后者为需要校正的图像）第一步：显示图像文件→在ERDASIMAGINE8.4图标面板工具条点击“Viewer”图标两次，打开两个视窗→ERDASIMAGINE8.4图标面板菜单条：Session→TileViewer然后，在Viewer#1中打开需要校正的LandsatTM图像：在Viewer#2中打开作为地理参考的校正过的SPOT图像：tmAtlanta.imgpanAtlanta.img第二步：启动几何校正模块Viewer#1菜单条：Raster→GeometricCorrection→打开SetGeometricModel对话框→选择多项式几何校正计算模型：Polynormial→单击OK按钮→同时打开GeoCorrectionTools对话框和PolynormialModelProperties对话框在PolynormialModelProperties对话框中，定义多项式模型参数及投影参数→定义多项式次方(Ploynomial)：2，其余取缺省值→单击Apply按钮→单击Close按钮→打开GCPToolReferenceSetup对话框说明：该实例是采用视窗采点模式，作为地理参考的SPOT图像已经含有投影信息，所以这里不需要定义投影参数。如果不是采用视窗采点模式，或者参考图像没有包含投影信息，则必须在这里定义投影信息，包括投影类型及其对应的投影参数。第三步：启动控制点工具在GCPToolReferenceSetup对话框中→选择视窗采点模式：ExistingViewer→单击OK按钮（关闭GCPToolReferenceSetup对话框）→打开ViewerSelectionInstruction指示器→在显示作为地理参考图像panAtlanta.img的Viewer#2中点击左键→打开ReferenceMapInformation提示框→单击OK按钮（关闭ReferenceMapInformation提示框）此时整个屏幕将显示：两个主视窗、两个放大窗口、两个关联方框（分别位于两个视窗中，指示放大窗口与主视窗的关系）、控制点工具窗口、几何校正工具。表明控制点工具启动，进入控制点采集状态第四步：采集地面控制点控制点工具窗口简介：控制点工具窗口（GCPTool）有菜单条、工具条、控制点数据表和状态条四个部分组成。GCPTool窗口菜单命令及其功能命令功能File:文件操作：LoadInput调用输入控制点文件（*.gcc）SaveInput将输入控制点保存在图像中SaveInputAs保存输入控制点文件（*.gcc）LoadReferenceSaveReferenceSaveReferenceASClose调用参考控制点文件（*.gcc）将参考控制点保存到图像中保存参考控制点文件（*.gcc）关闭控制点工具View:显示操作：ViewOnlySelectedGCP’sShowSelectedGCPinTableArrangeFramesonScreenTools视窗仅显示所选择的控制点在表格显示所选择的控制点重新排列屏幕中的组成要素调出控制点工具图标面板重新打开放大窗口StartChipViewerEdit:编辑操作：SetPointType(Control/Check)ResetReferenceSourceReferenceProjectionPointPrediction设置采集点的类型（控制／检查）改变参考控制点源文件改变参考文件的投影参数按照转换方程计算下一点位置匹配两幅图像像元的灰度值PointMatchingHelp:联机帮助：HelpforGCPTool关于GCP工具的联机帮助GCPTool窗口工具图标命令及其功能(按图标从左到右的顺序)命令功能Toggle自动GCP编辑模式开关键Calculate依据控制点求解几何校正模型设置自动转换计算开关AutomaticComputeErrorSelectGCPCreateGCPLock/UnlockFindinInputFindinReferZValue计算检查点的误差，更新RMS误差激活GCP选择工具，在视窗中选择GCP在视窗中选择定义GCP锁住／释放当前命令，以便重复使用选择寻找输入图像中的GCP选择寻找参考文件中的GCP计算更新所选GCP的Z值Auto-ZValue自动更新所有GCP的Z值GCP窗口数据字段及其含义字段含义GCP顺序号,系统自动产生GCP标识码,用户可以定义GCP当前选择状态提示符号输入GCP显示颜色输入GCP的X坐标输入GCP的Y坐标参考GCP显示颜色参考GCP的X坐标参考GCP的Y坐标GCP的类型(控制点／检查点)单个GCP的X残差Point#PointID>ColorXInputYInputColorXReferenceYReferenceTypeXResidualYResidualRMSError单个GCP的Y残差单个GCP的RMS误差Contribution单个GCP的贡献率Match两幅图像像元灰度值的匹配程度关于GCPTool窗口的说明:1.输入控制点(InputGCP)是在原始文件视窗中采集的，具有原文件的坐标系统；而参考控制点(ReferenceGCP)是在参考文件视窗中采集的，具有已知的参考坐标系统，GCP工具将根据对应点的坐标值自动生成转换模型。2.在GCP数据中，残差(Residuals)、中误差(RMS)、贡献率(Contribution)及匹配程度(Match)等参数，是在编辑GCP的过程中自动计算更新的，用户是不能任意改变的，但可以通过精确GCP位置来调整3.每个IMG文件都可以有一个GCP数据集与之相关联，GCP数据集保存在一个栅格层数据文件中；如果IMG文件有一个GCP数据集存在的话，只要打开GCPTool，GCP点都会出现在视窗中4.所有的输入控制点(InputGCP)都可以保存在图像文件中(SaveInput)，也可以保存在控制点文件中(SaveInputAs)。如果保存在图像文件中，调用的方法如3所述，如果是保存在控制点文件中，可以通过加载调用(LoadInput)。5.参考控制点(ReferenceInput)也可以类似地保存在参考图像中(SaveReference)或控制点文件中(SaveReferenceAs)，可以通过加载调用(LoadReference)采集地面控制点具体过程如下：→在GCPTool窗口中点击SelectGCP图标，进入GCP选择状态→在GCP数据表中将输入GCP的颜色（Color）设置为比较明显的黄色→在Viewer#1中移动关联方框位置，寻找明显的地物特征点，作为输入GCP→在GCPTool窗口中点击CreateGCP图标，并在Viewer#3中点击左键定点，GCP数据表将记录一个输入GCP，包括其编码、标识码。X坐标、Y坐标→在GCPTool窗口中点击SelectGCP图标，重新进入GCP选择状态→在GCP数据表中将参考GCP的颜色（Color）设置为比较明显的红色→在Viewer#2中移动关联方框位置，寻找对应的地物特→在GCPTool窗口中点击CreateGCP图标，并在Viewer#4中点击左键定点，系统将自动把参考点的坐标（XReference，YReference）显示在GCP数据表中征点，作为参考GCP→在GCPTool窗口中点击SelectGCP图标，重新进入GCP选择状态，并将光标移回到Viewer#1，准备采集另一个输入控制点→不断重复以上步骤，采集若干GCP，直到满足所选定的几何校正模型为止（多项式变换的次方数为二时，最少需要六个控制点），而后每采集一个InputGCP，系统就自动产生一个ReferenceGCP。第五步：采集地面检查点以上所采集的GCP的类型(Type)均为Control(控制点),用于控制计算,建立转换模型及多项式方程.下面所要采集的GCP的类型均是Check(检查点),用于检验所建立的转换方程的精度和适用性。在GCPTool窗口中→在GCPTool窗口菜单条中确定GCP类型：Edit→SetPointType→Check→在GCPTool窗口菜单条中确定GCP匹配参数（MatchingParameter）：Edit→PointMatching→打开GCPMatching对话框在GCPMatching对话框，需要定义下列参数：MatchingParameter（匹配参数）:→Max.SearchRadius（最大搜索半径）：3→SearchWindow（搜索窗口大小）：X:5；Y:5ThresholdParameters（约束参数）：→CorrectionThreshold（相关阈值）:0.8→选定DiscardUnmatchedPoint（删除不匹配的点）→单击Close（关闭GCPMatching对话框）→确定地面检查点：在GCPTool工具条中选择CreateGCP图标，并将Lock图标打开锁住CreateGCP功能，象选择控制点一样，分别在Viewer#1和Viewer#2中定义检查点，定义完毕后点击Unlock图标，解除CreateGCP功能→计算检查点误差：在GCPTool工具条中点击“ComputeError”图标，检查点的误差值会显示在GCPTool的上方，只有所有检查点的误差均小于一个像元，才能继续进行合理的重采样第六步：图像重采样重采样过程就是依据未校正图像像元值计算生成一幅校正图像的过程，原图像中所有栅格数据层都将进行重采样。ERDASIMAGINE8.4提供三种最常用的重采样方法：NearestNeighbor：邻近点插值法，将最邻近像元值直接赋予输出像元BilinearInterpolation:双线性插值法，用双线性方程和2×2窗口计算输出像元值CubicConvolution:立方卷积插值法,用立方方程和4×4窗口计算输出像元值具体操作如下:在GeoCorrectionTools对话框中选择“ImageResample”图标→打开ImageResample对话框在Resample对话框中,定义重采样参数:→OutputFile(输出图像文件):rectify.img→ResampleMethod(选择重采样方法)：NearestNeighbor→OutputCellSize(定义输出图像范围)：ULX、ULY；LRX、LRY→选定IgnoreZeroinStats→RecalculateOutputDefault(设置重新计算缺省值)：SkipFactor:X:1;Y:1→单击Ok（关闭Resample对话框）第八步：保存几何校正模型在GeoCorrectionTools对话框中点击Exit按钮，退出图像几何校正过程。按照系统提示将图像几何校正模型保存为模式文件resample.gms.第九步：检验校正结果检验校正结果的基本方法是：同时在两个视窗中打开两幅图像，其中一幅是校正后的图像，一幅是当时的参考图像，通过视窗地理连接（GeoLink／Unlink）功能及查询光标功能进行目视定性检验，具体过程如下：→在ERDASIMAGINE8.4图标面板具工条点击“Viewer”图标两次，打开两个视窗→ERDASIMAGINE8.4图标面板菜单条：Session→TileViewers然后，在Viewer#1中打开校正后的图像：rectify.img在Viewer#2中打开参考图像：panAtlanta.img→在Viewer#1中：按住右键→快捷菜单→GeoLink/Unlink→在Viewer#2中：点击左键→建立与Viewer#1的连接→在Viewer#1中：按住右键→快捷菜单→InquireCursor→打开光标查询对话框→在Viewer#1中：移动查询光标,观测其在两屏幕中的位置及匹配程度,并注意光标查询对话框中数据的变化。如果满意的话，关闭光标查询窗口。②地图采点模式：如果只有硬拷贝的地图或坐标纸作为参考，可以采用地图采点模式。首先在地图上选点并量算其坐标，在地图上寻找六个以上明显的地物特征点，并量算起坐标（多项式变换的次方数为二时，最少需要六个控制点）,然后通过键盘输入坐标数据。(以c:\ProgramFile\IMAGINE8.4\examples\nj321.img为例)首先在一个二维视窗中打开nj321.img图像文件

视窗菜单条：Raster→GeometricCorrection→打开SetGeometricModel对话框→选择多项式几何校正计算模型：Polynormial(注：卫星图像校正过程中多采用Polynormial（多项式变换）→单击OK按钮→同时打开GeoCorrectionTools对话框和PolynormialModelProperties对话框在PolynormialModelProperties对话框中，定义多项式模型参数及投影参数:→定义多项式次方(Ploynomial)：2，其余取缺省值→单击Apply按钮→单击Close按钮→打开GCPToolReferenceSetup对话框→选择KeyboardOnly（通过键盘输入控制点）→单击Ok→打开ReferenceMapInformation对话框→在MapUnits（地图单位）中选择：Other（说明：因为本例中采用的地图单位为km,在该下拉菜单中没有列出）→单击OK→打开GCPTool对话框在显示nj321.img图像的视窗Viewer#1中产生一个关联方框，同时产生一个放大窗口Viewer#2，放大窗口显示关联方框中图像的放大效果。→在Viewer#1中用关联方框查找与地图对应的一个控制点，单击GCPTool对话框工具条的“CreateGCP”图标，在放大窗口Viewer#2中点击左键定点，XInput和YInput值会自动产生，在XRef.和YRef.中输入控制点在地图上的坐标→按照类似方五个控制点的坐标→在GCPTool窗口菜单条中确定GCP类型：Edit→SetPointType→Check→当在放大窗口Viewer#2中点击第七个控制点时，XRef.和YRef.值会自动产生。法输入其余的左键定点→继续在Viewer#1中寻找四个明显的点，在放大窗口Viewer#2中点击左键定点，这五个点作为检查点→计算检查点误差：在GCPTool工具条中点击显示在GCPTool的上方，均小于一个像元（即<30米），才能继续进直到满足要求为止。“ComputeError”图标，检查点的误差值会只有所有检查点的误差行合理的重采样，如果误差大，则必须重复以上操作→单击GeoCorrectionTools对话框的“Resample”图标→打开Resample对话框在Resample对话框中,定义重采样参数:→OutputFile(输出图像文件名):nj321resample.img→ResampleMethod(重采样方法)：NearestNeighbor（邻近点插值法）→选定IgnoreZeroinStats→其它设置均取系统默认值→单击OK（执行重采样）3.图像拼接处理（以c:\ProgramFile\IMAGINE8.4\examples下的wasia1_mss.img、wasia2_mss.img、wasia3_mss.img）图像拼接处理是将具有地理参考的若干相邻图像合并成一幅图像或一图组像，需要拼接的输入图像必须含有地图投影信息，或者说输入图像必须经过几何校正处理或进行过校正标定。具体操作过程如下：→ERDASIMAGINE8.4图标面板菜单条：Main→DataPreparation(或单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标)→打开DataPreparation对话框→MosaicTool菜单→打开MosaicTool对话框MosaicTool视窗菜单条：Edit→AddImage→打开AddImageforMosaic对话框在AddImageforMosaic对话框中，需要设置下列参数：→ImageFilename（拼接图像文件）：→ImageAreaOptions(图像拼接区域)：ComputeActiveAreawasia1_mss.img→单击Add（图像wasia1_mss.img被加载到Mosaic视窗中）→按照同样的过程加载wasia2_mss.img和wasia3_mss.img→单击Close（关闭AddImageforMosaic对话框）→MosaicTool视窗菜单条：Edit→ImageMatching→打开MatchingOptions对话框→在MatchingOptions对话框中作如下设置：→在MatchingMethod中选择OverlapArea→单击OK按钮（关闭MatchingOptions对话框）→MosaicTool视窗菜单条：Edit→SetOverlapFunction→打开SetOverlapFunction对话框在SetOverlapFunction对话框，设置以下参数：→在IntersectionType中设置NoCutlineExits（没有裁切线）→在SelectFunction中设置Average（均值）→单击Apply→单击Close（关闭SetOverlapFunction对话框）→MosaicTool视窗菜单条：Process→RunMosaic→打开RunMosaic对话框在RunMosaic对话框中设置下列参数：→OutputFileName（确定输出文件名）：wasia_mosaic.img→WitchOutputs（确定输出图形区域）：All→在IgnoreInputValue中设置0→在OutputBackgroundValue中设置0→选定StatsIgnoreValue，并设置为0→单击OK（关闭RunMosaic对话框，运行图像拼接）→MosaicTool视窗菜单条：File→Close→打开VerifySaveonClose对话框→单击“否”按钮(关闭MosaicTool视窗,退出Mosaic工具)4.生成三维地形表面三维地形表面工具允许在不规则空间点的基础上产生三维地形表面，括：ASCII码点文件、ArcInfo的Coverage点文件和线文件、ERDASIMAGINE的注记数据层以及栅格图像文件IMG。所支持的类型包①以ASCII码点文件为例具体操作过程如下：(数据源以c:\ProgramFile\IMAGINE8.4\examples\lnpts.dat为例)第一步：启动三维地形表面ERDASIMAGINE8.4图标面板菜单条：Main→DataPreparation(或单击ERDASIMAGINE8.4图标面板工具条“DataPrep”图标)，→打开DataPreparation对话框→单击CreateSurface选项→打开3DSurfacing对话框第二步：定义地形表面参数→3DSurfacing对话框菜单条：File→Read→打开ReadPoints对话框→在SourceFileType中设置：ASCIIFile→在SourceFileName中设置：inpts.dat→单击OK（关闭ReadPoints对话框）→打开ImportOptions对话框在ImportOptions对话框中，需要设置如下参数：→FieldType（字段类型）：选择DelimitedbySeparator(分割字符)→SeparatorCharacter(分割字符)：Comma（逗号分割）→RowTerminatorCharacter(每行结束字符)：ReturnNewLine（DOS）→NumberofToSkip(跳过行数)：0→点击InputPreview选项，进入InputPreview栏从ImportOptions对话框（InputPreview栏）显示的原始数据可以看到数据的记录方式是一行一个点，每一行数据包括点号、X坐标、Y坐标、高程值四个字段，每两个字段之间的间隔号在此处读入时不需要，因此，必须在ImportOptions对话框的ColumnMappingZ与数据文件中字段的对应关系：→OutputColumnName：X对应InputFieldNumber:2→OutputColumnName：Y对应InputFieldNumber:3→OutputColumnName：Z对应InputFieldNumber:4→单击Ok（关闭ImportOptions对话框，读入数据）中确定X、Y、→返回3DSurfacing对话框第三步：生成三维地形表面→3DSurfacing对话框菜单条：Surface→Surfacing→打开Surfacing对话框在Surfacing对话框中，设置如下参数：→OutputFile（输出文件名）：surface.img→SurfacingMethod(表面插值方法)：LinearRubberSheeting（线性插值）→OutputCorners（输