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

实验三、遥感图像的几何校正与裁剪实验内容:1.图像分幅裁剪(SubsetImage2.图像几何校正(GeometricCorrection3.图像拼接处理(MosaicImgaes4.生成三维地形表面(3DSurfacing1.图像分幅裁剪在实际工作中,经常需要根据研究工作范围对图像进行分幅裁剪,按照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对话框场协调控制电梯和解救被困人员,同时将被困电梯的具体方位、电梯编号、停留楼层通知工程部和电梯公司进行解救被困人员。4.2.监控中心通过电视监控屏观察电梯内人员情况,使用电梯应急电话与被困人员联系,做好解释安慰工作。4.3.如遇特殊情况无法消除故障或被困人员中有身体严重不适时,应立即报消防部门或求助有关部门解决。4.4.当值领班管理员主管填写特发事件报告,做好详细记录故障情况及处理经过。5.水浸处理程序5.1.当值各岗管理员如发现大厦范围内水浸,应立即将出事地点和情况报告当值领班、管理员主管,同时尽快采用就近的防水设施保护好受浸楼层各电梯口,以免电梯受损。5.2.当值领班、管理员主管接报后立即赶往现场查看情况,组织抢险并及时向上级领导(物业主任、物业经理和公司领导及工程部汇报情况。5.3.设法查明浸水原因,并采取有效阻截措施,如水浸是来自大厦外的暴雨洪水,应当在各低于水位的出入口使用备用拦水闸板和沙包;如水浸是来自市政地下水反溢,应当暂时将反溢的地下水道能往大厦的入口封闭,并用排水水泵将大厦的积水抽排出大厦外;如水浸来自大厦机管设施的损坏或故障,应当先关闭控制有关故障部位的水闸或供水泵。5.4.组织当值人员根据水浸情况,协同工程部采取有效措施,如将电梯开高离开受浸范围、关闭受浸区域的电闸、在水漫延的通道上摆放拦水沙包、疏通排水沟、开启排水泵、用吸水机等,尽可能减少水浸所致的损失。5.5.水浸排除后,立即通知清洁公司人员清除积水并清理现场环境,通知工程部查明原因,修复受损的设施,尽快恢复大厦的正常运作。5.6.当值领班、管理员主管填写特发事件报告,写出事件详细经过。6.停电事故的处理程序6.1.接通知大厦将在短时间内停电,应在停电前10分钟,用紧急广播通知用户,工程部安排电梯工,提前将电梯停在大厦首层;物业部派管理员提前10分钟到达电梯厅维持秩序,阻止客人使用电梯,以免停电期间有客人困梯。6.2.未预知的情况下大厦突然发生停电,立即联系工程部,尽快采取措施恢复供电。6.3.使用紧急照明,保证公共地方及主要通道的照明。6.4.派管理员到各主要电梯厅及大厦出入口维持秩序。6.5.消防监控中心和巡楼保安员密切注视大厦各楼层,以防有人趁机制造治安问题。6.6.当值领班、管理员主管填写特发事件报告,写出事情详细经过。7.业主/租户单元内发生刑事和治安灾害事故的处理程序。7.1.当值管理员应迅速向当值领班、管理员主管、物业主任、物业经理汇报或向公安部门报案,业主/租户受侵害的财物投购保险的,由被保人通知承保的保险公司。7.2.当值领班、管理员主管应根据具体情况,采取适当的方法把整个现场保护起来,禁止无关人员进入现场,以免破坏现场遗留的痕迹、物证、影响证据的收集。7.3.抓紧时机向发现人或周围群众了解案件、事故发生发现经过,收集群众的反映和议论,了解更多的情况并认真记录。7.4.向到达现场的公安人员认真汇报案件发生情况,协助破案。7.5.当值领班、管理员主管填写特发事件报告,写出事情详细经过。8.抢劫的处理程序在当值中如遇有公开使用暴力,威迫或其他手段如打、砸、抢、强行掠取公司或业主/租户财物的犯罪行为时,必须立即按如下程序处理。8.1.迅速制止犯罪,设法擒获或召集组织其他人员擒获罪犯,同时向当值领班、管理员主管、物业主任、物业经理报告或公安部门报案。8.2.如在楼内发生劫案,应通知大堂管理员留意,严格检查,组织人员进行清查,如罪犯逃走,追不上时,应看清楚犯罪分子的人数、衣着、面貌明显特征、身体特点、所用的交通工具及其型号、品牌、特征等,及时报告上级领导、公安部门。8.3.保护好现场,劫匪遗留凶器、作案工具等不要用手触摸,不要让无关人员进入现场,如无固定的现场的、无法将劫匪遗留的物品留在原处,要一一收拾起来,交给公安部门处理。8.4.访问目击群众收集发生劫案情况,提供给公安部门,同时公安部门未勘查现场或未勘查完毕不能离开。8.5.事主或在场如有受伤的,要立即设法将伤者送医院抢救医治,并迅速报告公安部门。8.6.当值领班、管理员主管填写特发事件报告,写出事件详细经过。9.对不执行规定不听劝阻的处理程序在当值过程中如遇到不登记不出示证件强行进入,强行在不准停车的地方停车,在不该堆放垃圾的地方堆放垃圾等,可按如下程序处理:9.1.处理问题时要有礼貌、态度和蔼、说话和气、以理服人。9.2.对不停劝阻者,尽量想办法查清其姓名或单位,如实记录下来以便向领导汇报。9.3.发生纠纷,管理人员一定要保持冷静,避免争吵。9.4.若遇到个别人员蛮横无理,打骂管理人员,可通知上级领导同此人共同协商妥善处理,若情节严重,报告公安部门依法处理。9.5.当值领班、管理员主管填写特发事件报告,写出事件详细经过。10.对醉酒闹事或精神病人的处理程序醉酒者或精神病人失去了正常的理智,有些处于不能自控的状态,很容易造成伤害,当值人员应按如下程序处理:10.1.进行劝阻或阻拦,让其离开大厦管辖范围。10.2.如果熟悉或认识醉酒者或精神病人,应设法立即通知其家人或工作单位,请他们派人领回,并采取监控和监护措施。10.3.若醉酒者或精神病人有危害社会安全的行为,可将其强制送交公安部门,处理过程中要注意保护自己。11.在当值中拾到遗失物品的处理程序11.1.群众交来或拾获群众遗失或遗留的物品,都要一一登记清楚,并报告当值保安主任或上级领导。11.2.处理拾遗物品,拾遗物品中有证件提供失主的详细地址或电话的,可通知失主前来认领,但要问明情况,如是贵重物品,必须交由部门领导处理。11.3.如拾遗物品是属于违禁品的登记后上报部门、公司领导与公安部门协调处理。11.4.对拾遗物品,必须报告部门领导,不得拾物不报、少报或私自挪用,违者重处,情节严重的以贪污论处。(九、常见的几种可疑情况及处理措施:加强巡逻检查是治安防范工作中预防、发现和打击违法犯罪分子的一项行之有效的措施。而对可疑情况视级别及采取相应的措施是巡逻工作的重点,以下几种常见可疑情况及处理方法:1.可疑情况:在大厦内游荡或借口找人说不出被访人姓名及所在单元和公司名称,或接二连三的进去单元内的。处理方法:密切注意其举动,必要时劝其离开。2.可疑情况:发现其身上带有铁棒、刀具、钳子、螺丝刀等工具的。处理方法:检查其扔带工具的用途,如用途不明的,先约束起来,送公安部门。3.可疑情况:携带物品繁多,如电脑、打印机等贵重物品,又无任何证明的(搬家有管理处放行条的除外。处理方法:暂时将人、物扣留,待其出具可靠证明后放行,如无任何证明,即通知部门领导或公安部门审查。4.可疑情况:在偏避、隐蔽处清理皮包或钱包的处理方法:立即设法拦截,询问证实,如属盗窃、抢劫财物的,报告部门领导或送交公安部门处理。5.可疑情况:单车、摩托车无牌、无证行驶、无钢印、有耗损痕迹的,或将未开锁的单车背走或提走的处理方法:当即扣留人、车,待查明后方行。6.可疑情况:遇到保安即转身远离或逃跑的人。处理问题:设法阻截(用对讲机向各岗位通告及当值保安主任报告擒获。7.可疑情况:低价出售物品与实际价值相差较大的。视窗菜单条:Raster→GeometricCorrection→打开SetGeometricModel对话框资源卫星图像校正具体过程:①视窗采点模式:如果已经拥有需要校正图像区域的数字图像、或经过校正的图像、或注记图层,就可以应用视窗采点模式,直接以数字地图、或经过校正的图像、或注记图层作为地理参考,在视窗中打开该图层,从中采集控制点。(以c:\ProgramFile\IMAGINE8.4\examples下的panAtlanta.img和tmAtlanta.img为例,前者是作为地理参考的图像,后者为需要校正的图像第一步:显示图像文件→在ERDASIMAGINE8.4图标面板工具条点击“Viewer”图标两次,打开两个视窗→ERDASIMAGINE8.4图标面板菜单条:Session→TileViewer然后,在Viewer#1中打开需要校正的LandsatTM图像:tmAtlanta.img在Viewer#2中打开作为地理参考的校正过的SPOT图像:panAtlanta.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:LoadInputSaveInputSaveInputAsLoadReferenceSaveReferenceSaveReferenceASClose文件操作:调用输入控制点文件(*.gcc将输入控制点保存在图像中保存输入控制点文件(*.gcc调用参考控制点文件(*.gcc将参考控制点保存到图像中保存参考控制点文件(*.gcc关闭控制点工具View:显示操作:ViewOnlySelectedGCP’sShowSelectedGCPinTableArrangeFramesonScreenToolsStartChipViewer视窗仅显示所选择的控制点在表格显示所选择的控制点重新排列屏幕中的组成要素调出控制点工具图标面板重新打开放大窗口Edit:SetPointType(Control/CheckResetReferenceSourceReferenceProjectionPointPredictionPointMatching编辑操作:设置采集点的类型(控制/检查改变参考控制点源文件改变参考文件的投影参数按照转换方程计算下一点位置匹配两幅图像像元的灰度值Help:HelpforGCPTool联机帮助:关于GCP工具的联机帮助GCPTool窗口工具图标命令及其功能(按图标从左到右的顺序命令功能Toggle自动GCP编辑模式开关键Calculate依据控制点求解几何校正模型Automatic设置自动转换计算开关ComputeError计算检查点的误差,更新RMS误差SelectGCP激活GCP选择工具,在视窗中选择GCPCreateGCP在视窗中选择定义GCPLock/Unlock锁住/释放当前命令,以便重复使用FindinInput选择寻找输入图像中的GCPFindinRefer选择寻找参考文件中的GCPZValue计算更新所选GCP的Z值Auto-ZValue自动更新所有GCP的Z值GCP窗口数据字段及其含义字段含义Point#GCP顺序号,系统自动产生PointIDGCP标识码,用户可以定义>GCP当前选择状态提示符号Color输入GCP显示颜色XInput输入GCP的X坐标YInput输入GCP的Y坐标Color参考GCP显示颜色XReference参考GCP的X坐标YReference参考GCP的Y坐标TypeGCP的类型(控制点/检查点XResidual单个GCP的X残差YResidual单个GCP的Y残差RMSError单个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中移动关联方框位置,寻找对应的地物特征点,作为参考GCP→在GCPTool窗口中点击CreateGCP图标,并在Viewer#4中点击左键定点,系统将自动把参考点的坐标(XReference,YReference显示在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工具条中点击“ComputeError”图标,检查点的误差值会显示在GCPTool的上方,只有所有检查点的误差均小于一个像元(即<30米,才能继续进行合理的重采样,如果误差大,则必须重复以上操作直到满足要求为止。→单击GeoCorrectionTools对话框的“Resample”图标→打开Resample对话框在Resample对话框中,定义重采样参数:→OutputFile(输出图像文件名:nj321resample.img→ResampleMethod(重采样方法:NearestNeighbor(邻近点插值法→选定IgnoreZeroinStats→其它设置均取系统默认值→单击OK(执行重采样3.图像拼接处理(以c:\ProgramFile\IMAGINE8.4\examples下的wasia1_mss