ENVI进行影像到影像的几何校正

卫星影像的几何校正以具有地理参考的SPOT 4 10m全色波段为基础，进行Landsat 5 TM 30m影像的几何校正过程，其流程如图1所示。图1 几何精校正流程目的：1、掌握利用地面控制点（GCP）进行影像到影像几何校正的方法2、影像上GCP的选取方法数据准备：bldr_tm.img 没有地理坐标的tm影像bldr_sp.img 带地理坐标的SPOT影像利用GCP进行几何校正的具体操作第一步 打开并显示影像文件（1）在#1窗口中打开bldr_tm.img作为待校正图像，在#2窗口中打开bldr_sp.img作为参考图像（图2）。图2 参考图像（左）与待校正图像（右）第二步 启动几何校正模块（1）一旦两幅图像都已经显示，选择主菜单MapRegistrationSelect GCPs: Image to Image，打开几何校正模块。（2）在Image to Image Registration对话框中，选择显示SPOT影像的Display作为基准图像（Base Image），显示TM影像的Display为待校正图像（Warp Image）（图3）。点击OK，进入采集地面控制点。图3 指定参考图像与待校正图像第三步 采集地面控制点（1）控制点工具对话框说明：图4 地面控制点工具对话框表1 菜单命令及功能菜单命令功能Filesave GCPs to ASCIIsave Coefficients to ASCIIrestore GCPs from ASCII保存GCP到ASCII文件保存多项式系数到ASCII文件从ASCII文件中打开GCPOptionWarp Displayed BandWarp FileWarp Displayed Band(as Image to Map)Warp File(as Image to Map)Reverse Base/Warp1st Degree(RST Only)Auto PredictLabel PointsOrder Point by ErrorClear All PointsSet Point ColorsAutomatically Generate Tie Points配准当前显示的波段配准整个文件校正当前显示的波段校正整个文件颠倒基准图像和被校正图像角色选择使用RST模型来计算误差打开/关闭自动预测点功能打开/关闭GCP标签打开/关闭根据误差从大到小对GCPs排序删除所有控制点设置控制点标示颜色启动自动寻找同名点（Tie）功能 当基准图像没有地理投影时选择这种配准命令；如果基准图像具有地理投影时选择此命令，得到的结果诸如投影参数、像元大小将与基准图像相同。 当基准图像有地理投影时，可以选择这种配准命令，在输出结果时候还可以更改校正图像的输出像元大小和投影参数表2其它功能按钮及功能选择文本框功能Base XBase YWarp XWarp YDegree基准图像上的Zoom显示窗口十字光标的X像素坐标（列数）基准图像上的Zoom显示窗口十字光标的Y像素坐标（行数）校正图像上的Zoom显示窗口十字光标的X像素坐标（列数）校正图像上的Zoom显示窗口十字光标的Y像素坐标（行数）预测控制点、计算误差（RMS）多项式次数按钮功能Add PointPredictShow/Hide ListDelete Last Point添加控制点预测点位置，当控制点数量达到多项式最少点要求是可用显示/关闭控制点列表删除最后一个收集的控制点标签意义Number of Selected PointRMS Error已收集的控制点个数累计误差（单位：像元） 当控制点数量达到一定数量时才能更改，如控制点数达到6，Degree值可以改为2，最大为3。（2）地面控制点采集在图像几何校正过程中，采集地面控制点是一项重要和繁重的工作，直接影响最后的校正结果，具体过程如下：1）在两个Display中移动方框位置，寻找明显的地物特征点作为输入GCP。2）在Zoom窗口中，通过将十字光标放置在两幅影像的相同地物点上。图5 采集地面控制点3）在Ground Control Points Selection上，单击Add Point按钮，将当前找到的地面控制点添加到地面控制点列表中。在Ground Control Points Selection上，单击Show List按钮，可以显示地面控制点列表（图6）图6 控制点列表4）用同样的方法继续寻找其余的点，当选择控制点的数量达到3时，RMS被自动计算。Ground Control Points Selection上的Predict按钮可用，这时在基准图像显示窗口上面定位一个特征点，单击Predict按钮，校正图像显示窗口上会自动预测区域，适当调整下位置，将控制点收集。随着控制点数量的增加，预测点的精度越来越精确。5）对于RMS过高的地面控制点，一般有两种处理方法，一是直接删除，选择要删除的控制点所在的行，单击Delete按钮；二是在两个图像的Zoom窗口上，将十字光标重新定位到正确的位置，单击Update按钮进行微调。5）当Ground Control Points Selection上的RMS值小于1个像素时，并且控制点的数量足够且分布均匀，完成控制点的选择（本次实验要求12个控制点）。控制点选择的示例如图7。6）保存控制点，Ground Control Points Selection菜单FileSave GCPs to ASCII。图7 控制点选择第四步 选择校正参数输出结果有两种校正输出方式：Warp File和Warp File (as Image Warp)1、Warp File（1）Ground Control Points Selection对话框上选择菜单OptionsWarp File ，选择校正图像（TM影像），打开校正参数设置对话框（图8）。图8 Warp File校正参数设置（2）在Registration Parameters对话框中：l 校正方法Method选择多项式Polynomiall 次数（Degree）选择2l 重采样（Resampling）选择Bilinear或者Cubic Convolutionl 背景值（Background）为0l Output Image Extent：默认是根据基准图像大小计算，可以进行适当调整l 选择输出路径和文件名单击OK，得到校正后影像。它的尺寸大小、投影参数和像元大小都与基准图像一致。2、Warp File (as Image Warp)（1）在Ground Control Points Selection对话框上，选择菜单OptionsWarp File(as Image Warp)，选择校正图像（TM影像），打开校正参数设置对话框（图8）。图9 Warp File (as Image Warp)校正参数设置（2）在Warp File (as Image Warp)校正参数对话框中，默认投影参数和像元大小与基准图像一致。各参数作如下设置：l 投影参数不变，在X和Y的像元大小输入30ml 校正方法Method选择多项式Polynomiall 次数（Degree）选择2l 重采样（Resampling）选择Bilinear或者Cubic Convolutionl 背景值（Background）为0l Output Image Extent：默认是根据基准图像大小计算，可以进行适当调整l 选择输出路径和文件名单击OK，得到校正后影