erdas 几何校正学习

erdas 几何校正学习

本人将TM（分辨率30米）影像GM07.img与1：10000深圳独立坐标的地形图校正，选择多项式几何校正模型，校正前在地形图上选择了18个参考点，均是道路交叉点，单独生成cad文件后用GIS生成coverage文件，在erdas中将coverage文件转成Raster文件“SZreference.img”，校正后的文件为“REgm07.img”。粗略步骤：使用Image Geometric Correction →打开Set Geometric Model对话框→选择多项式几何校正模型：Polynomial→OK →同时打开Geo Correction Tools 对话框和Polynomial Model Properties对话框。在Polynomial Model Properties对话框中，定义多项式模型参数以及投影参数：→定义多项式次方（Polynomial Order）:4 →定义投影参数：（PROJECTION）:略→Apply→Close →打开GCP Tool Referense Setup 对话框，在需校正图像和参考图像上采点，不断重复，→采好点→重采样，生成校正后的文件。

具体做法第一步：显示图像文件（Display Image Files）在ERDAS图标面板中点击Viewer图表两次，打开两个视窗（Viewer1/Viewer2），并将两个视窗平铺放置，操作过程如下：ERDAS图表面板菜单条：Session→Title Viewers 在Viewer1中打开需要校正的图像：GM07.img；在Viewer2中打开作为地理参考的Raster图像：SZreference.img

第二步：启动几何校正模块（Geometric Correction Tool）Viewer1菜单条：Raster→Geometric Correction →打开Set Geometric Model对话框）→选择多项式几何校正模型：Polynomial→OK →同时打开Geo Correction Tools对话框和Polynomial Model Properties对话框；在Polynomial Model Properties对话框中，定义多项式模型参数以及投影参数：→定义多项式次方（Polynomial Order）:4（若此处定义的次方数为T，则需配准的点数为（T+1）*（T+2）/2，若为4，应该配置15个点或更多）→定义投影参数：（PROJECTION）:略→Close →打开GCP Tool Referense Setup 对话框

第三步：启动控制点工具（Start GCP Tools）首先，在GCP Tool Referense Setup 对话框中选择采点模式：→选择视窗采点模式：Existing Viewer→OK →打开Viewer Selection Instructions指示器→在显示作为地理参考图像SZreference.img 的Viewer2中点击左键→打开reference Map Information 提示框；→OK→表明控制点工具被启动，进入控制点采点状态。

第四步：采集地面控制点（Ground Control Point）GCP的具体采集过程：在图像几何校正过程中，采集控制点是一项非常重要和繁重的工作，1、在GCP 工具对话框中，点击Select GCP图表，进入GCP选择状态；2、在Viewer1中移动关联方框位置，寻找明显的地物特征点，作为输入GCP。3、在GCP工具对话框中，点击Create GCP图标，并在Viewer3中点击左键定点，GCP数据表将记录一个输入GCP，包括其编号、标识码、X坐标和Y坐标。4、在GCP 对话框中，点击Select GCP图标，重新进入GCP选择状态。5、在Viewer2中，移动关联方框位置，寻找对应的地物特征点，作为参考GCP。6、在GCP工具对话框中，点击Create GCP图标，并在Viewer4中点击对应点，系统将自动将参考点的坐标（X、Y）显示在GCP数据表中。7、在GCP对话框中，点击Select GCP图标，重新进入GCP选择状态，并将光标移回到Viewer1中，准备采集另一个输入控制点。8、不断重复1-7(切记一定要先选择输入GCP，输入

对应的参考GCP后，再采集另外一对输入－参考GCP)，采集若干控制点GCP，直到满足所选定的几何模型为止。

第五步：采集地面检查点（Ground Check Point）以上采集的GCP的类型均为控制点，用于控制计算，建立转换模型及多项式方程。

第六步：计算转换模型（Compute Transformation）在控制点采集过程中，一般是设置为自动转换计算模型。所以随着控制点采集过程的完成，转换模型就自动计算生成（这也要求输入－参考GCP要一一对应）。在Geo-Correction Tools 对话框中，点击Display Model Properties 图表，可以查阅模型。

第七步：图像重采样（Resample the Image）重采样过程就是依据未校正图像的像元值，计算生成一幅校正图像的过程。原图像中所有删格数据层都要进行重采样。ERDAS IMAGE 提供了三种最常用的重采样方法。图像重采样的过程：在Geo-Correction Tools对话框中选择Image Resample 图标。在Image Resample 对话框中，定义重采样参数；→输出图像文件明（OutputFile）：REgm07.img →选择重采样方法（Resample Method）:Nearest Neighbor →定义输出图像范围：采用缺省值→定义输出像元的大小：采用缺省值→设置输出统计中忽略零值: →定义重新计算输出采用缺省值。

第八步：保存几何校正模式（Save rectification Model）在Geo-Correction Tools 对话框中点击Exit按钮，推出几何校正过程，按照系统提示，选择保存图像几何校正模式，并定义模式文件，下一次直接利用。

第九步：检验校正结果（Verify rectification Result）基本方法：同时在两个视窗中打开两幅图像，一幅是矫正以后的图像，一幅是当时的参考图像，通过视窗地理连接功能，及查询光标功能进行目视定性检验。我的方法是把行政区边界生成.aoi文件，加入校正后的图像进行比对，比对后顺便进行不规则裁剪了，具体见我的另外一个帖子“erdas中进行AOI不规则多边形裁剪――aoi，Subset”