envi中的几何校正.doc

几何校正几何校正1.遥感图像产生几何畸变的原因遥感图像产生几何畸变的原因地物目标发出的电磁波被卫星上所载传感器接收,这些电磁波上记录和传达了地物目 标的信息,这是遥感图像成像的过程也是它的内在规律在这个过程中图像的几何畸变也随 即产生了,其中原因很多,主要表现在以下几个方面: 1. 1 卫星位置和运动状态变化的影响 卫星围绕地球按椭圆轨道运动,引起卫星航高和飞行速度的变化,导致图像对应产生偏离与在 卫星前进方向上的位置错动另外,运动过程中卫星的偏航、 翻滚和俯仰变化也能引起图 像的畸变 以上误差总的来说,都是因为传感器相对于地物的位置、 姿态和运动速度变化产生的, 属于外部误差此外,由于传感器本身原因产生的误差,即内部误差,这类误差一般很小,通常 人们不作考虑 1. 2 地球自转的影响 大多数卫星都是在轨道运行的降段接收图像,即当地球自西向东自转时,卫星自北向南运 动这种相对运动的结果会使卫星的星下位置产生偏离,从而使所成图像产生畸变 1. 3 地球表面曲率的影响 地球表面是不规则的曲面,这使卫星影像成像时像点发生移动,像元对应于地面的宽度不 等特别是当传感器扫描角度较大时,影响更加突出。

1. 4 地形起伏的影响 当地形存在起伏时,使原来要反映的理想的地面点被垂直在其上的实际某高点所代替,引 起图像上像点也产生相应的偏离 1. 5 大气折射的影响 由于大气圈的密度是不均匀分布的,从下向上越来越小,使得整个大气圈的折射率不断变 化,当地物发出的电磁波穿越大气圈时,经折射后的传播路径不再是直线而是一条曲线,从而 导致传感器接收的像点发生位移 2.进行几何校正并保证精度的必要性进行几何校正并保证精度的必要性 遥感图像几何校正的精确与否直接关系到应用遥感信息反应地表地物的地理位置和面 积的精确度,关系到从图像上获取的信息准确与否,因此在选择控制点上要十分小心,尽可能 提高其精度,并且要对校正结果进行反复的分析比较,必要时还要进行多次校正几何校正让 图像上地物对应的像元出现在它应该在的地方,再通过辐射校正、 影像增强等遥感图像处 理技术,还图像以“ 本来面目” 然后通过对图像的识别、 分类、 解译处理实现地面空 间上各类资源信息的空间分析研究,使遥感技术投入到实际生产应用中 3.几何精校正几何精校正 遥感影像图的几何校正目前有 3 种方案,即系统校正、 利用控制点校正以及混合校正。

遥感数据接收后,首先由接收部门进行校正,这种校正叫系统校正 (又叫几何粗校正) ,即把遥 感传感器的校准数据、 传感器的位置、 卫星姿态等测量值代入理论校正公式进行几何畸 变校正;而用户拿到这种产品后,由于使用目的不同或投影及比例尺不同,仍旧需要做进一步 的几何校正,这就需要对其进行几何精校正即利用地面控制点 GCP ( GroundContr ol2Point, 遥感图像上易于识别,并可精确定位的点 )对因其他因素引起的遥感图像几何畸变进行纠正 混合校正则是由一般地面站提供的遥感 CCT 已经完成了第一阶段的几何粗校正,用户所要 完成的仅仅是对图像做进一步的几何精校正 几何精校正就是利用地面控制点 GCP 对各种因素引起的遥感图像几何畸变进行校正 从数学上说,其原理是通过一组 GCP 建立原始的畸变图像空间与校正空间的坐标变换关系,利用这种对应关系把畸变空间中全部元素变换到校正空间中去,从而实现几何精校正 在 ENVI 中对影像进行几何精校正步骤如下： 将原始数据 LS5_TM_20090324_032106_032131_129042_FASTB_L2.zp 和 LS5_TM_20090324_000000_000000_129042_GEOTIFF_L4.zip 分别解压到相应的文件夹， 为了便于识别，将 LS5_TM_20090324_032106_032131_129042_FASTB_L2 重命名为 TM2009_129042_L2，将 LS5_TM_20090324_000000_000000_129042_GEOTIFF_L4 重命名 为 TM_2009129042_L4。

打开 ENVI4.6 软件，并加载 TM2009_129042_L2 和 TM_2009129042_L4 两幅影像，加载 TM2009_129042_L2 步骤如下：打开 file-open external file-landsat-fast 此时弹出加载数据选 项窗口选择 header.dat，点击“打开”按钮，显示结果如图选择 RGB Color,选择 543 波段组合，并将其显示在 Display1 中加载 TM_2009129042_L4 步骤如下：打开 file-open image file 或者在 Available Bands list 窗 体中选择 file-open image file 此时弹出加载数据选项窗口单击“打开”按钮，此时所需的两幅影像波段信息显示在 Available Bands list 窗体中，如 图将 L5-TM-129-042-20090324-L 下的波段以 RGB 方式按 543 波段组合并显示于 Display2 中打开 Map-Registration-Select GCPs: Image to Image,弹出如下图所示窗口：在左边“base image”下提供可供选择的基准图像，在“warp image”下提供可供选择的待 校正图像。

由于加载于 Display1 中的图像是 TM2009_129042_L2 为二级产品是没有经过几何精校正的 图像，而加载于 Display2 中的图像 TM2009_129042_L4 为四级产品已经过几何精校正，可 以作为基准图像以校正未经过几何精校正的二级产品；因此在上图中选择 Display2（即基 准图像 TM2009_129042_L4） ，而 warp image 下应选择 display1（即待校正图像 TM2009_129042_L2） 待基准图像和待校正图像选择完毕后，弹出窗口此时可以在 Display1 和 Display2 中选择同名点并增加控制点，选择同名点时，应注意尽量 避免选择状态时常变动的，如河流等，因为河流由于流水水位的不确定性和不稳定性，很 容易造成人为控制误差；应尽量选择状态相对稳定的地物作为同名点，如道路交叉的地方， 小型建筑物（如桥梁） ，山峰等在两个显示窗口中确定好同名点后点击如上图列表中“Add point”按钮以确认控制点的选择，如此继续增加控制点；知道控制点误差不再显示 为零（在 envi4.6 中一般选择控制点数多达 4 个便可以识别误差） ，可以启动自动预测功能 （注：只能进行自动预测而不能自动加点，预测出的区域还需用鼠标重新定位进行微调， 并点击 Add point 增加控制点） ，自动预测功能的启动为：在 Ground Control Point Select 窗 口中，选择 options-auto predict。

查看误差列表可选择上图中“show list”按钮，此时弹出如下显示误差信息的窗口，可以 让列表中的误差（RMS）从大到小排列（点击如下图的 options，在下拉列表中选择 order point by error,也可在 Ground Control Point Select 窗口中选择 options，再勾选 options 下拉列 表中的 order point by error） ，控制点的误差（RMS）须控制在 0.5 以下，不满足要求的控制 点应定位（选择所要调整的控制点如下图蓝色阴影部分）到控制点所在位置进行微调，微 调完成后点击下图中”update”按钮以更新控制点,同样如要删除控制点也只需选中需要删除 的控制点，点击“Delete”按钮完成控制点选择后如下图所示：控制点选择完成后应对控制点进行保存，保存方法为：在 Ground Control Point Select 窗口 中选择 file-save GCPs to ASCII，将控制点保存到相应位置以便以后使用 确认控制点误差在 0.5 以下后便可以进行校正处理，在 Ground Control Point Select 窗口中 选择 options-warp files，弹出窗口选择待校正文件 header.dat（即 TM2009_129042_L2） ，点击“ok”弹出如下校正参数设置 窗口（Registration Parameters）上图中 method 为几何校正计算模型，RST 为仿射变换模型，为仿射变换模型，Polynomial 为为多项式模型，多项式模型，Triangulation 为局部三角网模型为局部三角网模型。

多项式模型在卫星图像校正过程中应用较多在应用此模型时，需要确定多项式的次方数，通常选择 2 次或 3 次 选择的次方数与所需要的最少控制点是相关的，最少控制点计算公式为（n+1）*(n+2)/2， 其中 n 为次方数，即选择多少次多项式 n 就是多少，而所需要的最少控制点数目即为 （n+1）*(n+2)/2，如利用二次多项式模型，所需要的最少控制点数目就是 6 个上图中degree 的数值就代表 n，resampling 表示重采样方式，在此提供了三种重表示重采样方式，在此提供了三种重采样方式：最邻近点法（采样方式：最邻近点法（nearest neighbor）） ，双线性内插法，双线性内插法（（bilinear）以及三次卷积内插法（）以及三次卷积内插法（cubic convolution））可根据需要选择不同的内插方法进行重采样，选择不同的采样方法效果不一样 其他参数为默认值，校正参数设置完成后，选择校正文件的输出路径，点击“ok”执行几 何校正 检验几何校正的结果可以检测校正效果，具体方法是：同时在两个窗口中打开图像，其中 一幅是校正后的图像，一幅是基准图像，通过视窗链接（link displays）及十字光标或者地 理链接（geographic link）进行关联。

在显示校正结果的 image 窗口中，从右键快捷菜单中选择 geographic link 命令，选择需要 链接的两个窗口，打开十字光标进行查看，如图所示：若校正结果不理想可打开控制点文件，对控制点进行调整，并重新执行几何校正。