实习序号及题 目 实习 2. 影像几何精纠正 实习人姓名 专业班级及编号 09 地基 01 号 任课教师姓名 实习指导教师姓名 实习地点 榆 中 校 区 实 验 楼 A209 实习日期时间 2011 年 11 月 1 日 实习目的 （1）深入理解影像几何精纠正的原理 （2）学会使用影像对影像的几何精纠正方法和具体操作步骤 实习内容 （1）了解和熟悉控制点的选择方法，进行控制点的选择 （2）利用纠正好的全色波段高分辨率影像完成同景多光谱影像的几何精纠正，并进行 精度评价和对比 基本原理 几何精纠正的基本原理是回避成像的空间几何过程，直接利用地面控制点数据对遥 感图像的几何畸变本身进行数学模拟，并且认为遥感图像的总体畸变可以看作是挤压、 扭曲、缩放、偏移以及更高次的基本变形的综合作用的结果因此，校正前后图像相应 点的坐标关系可以用一个适当的数学模型来表示 几何精纠正的过程就是：利用地面控制点数据确定一个模拟几何变形畸变的数学模 型，以此建立原始图像空间与标准空间的某种对应关系，然后利用这种对应关系，把畸 变图像空间中的全部像素转换到标准空间中，从而实现图像的几何精纠正。
ENVI 几何校正模型分为以下三类：仿射变换（RST） ，多项式模型（Polynomial） ， 局部三角（Delaunay Triangulation） 下面着重介绍多项式模型（Polynomial） ： 几何精纠正的关键步骤如下： 1、控制点的选取 ①原则： 控制点的空间分布应尽量做到在影像工作区范围内尽量做到均匀分布； 在影像上具有明显、清晰的地物特征的标志点，如道路交叉点、河流汇合口，田间 地块的转折点； 地面控制点对应地物要相对稳定，不能随时间、季节而发生变化，以保证当两幅不 同时相或相邻的影像做地图几何精纠正时，均可以被正确定位； 在大多为山区的影像上选取控制点时，应尽量选择同高程的控制点（通常是在平原、 低地地区） ，不要从山地顶部选取控制点（除非采用 DEM 支持的正射影像纠正法） ； 图像的边缘一定要选，防止外推，特征变化大的地方多选； 进行最后的插值和重采样前对所有 GCP 进行筛选，剔除其中误差较大的点； 单点定位精度一般应小于 0.5 个像元大小，至少不能超过一个像元 ②控制点的数目 ： 最少控制点的个数：L=(n+1)(n+2)/2，n 为多项式阶数。
2、建立坐标变换函数 建立两图像像元点之间的对应关系； x,y: 变换前的图像坐标； u,v： 变换后的图像坐标； 通常数学关系 f 表示为二元 n 次多项式： 运用最小二乘法进行拟合求系数（其中 L 为计算样本点，或者 GCP 的个数） ： 运用已知 GCP 点对数据进行矩阵运算求解多项式系数 3、空间插值 ①前向映射插值 方法：计算输入影像中各个像元在变换后输出影像坐标系上相应的位置，把各个像 元的数据投影到该位置上 存在问题：输出影像上各个栅格点一般不会与这些投影点在位置上相对应，由这些 投影点解算各个栅格的灰度值仍存在算法实现上的复杂性 ②逆向映射插值 利用反变换关系式计算输出影像上各个像元在输入影像中相应的坐标位置进而采 用空间插值法计算出该点的灰度值，即为输出影像上对应点的灰度值此系目前商用软 件中的通用方法 4、重采样（灰度插值） 纠正后的新图像的每一个像元，根据变换函数，可以得到它在原始图像上的位置 如果求得的位置为整数（在原图的像元中心） ，则该位置处的像元灰度就是新图像的 灰度值多数位置不为整数，需根据（x,y） （变化前的图像坐标）与周围像素之间的关系 内插产生新的像素值，然后把像素值写到（u,v） （变化后的图像坐标） 。
),( ),(vufx vufyxy 21002100minminLnn i ij xlijll lijLnn i ij ylijll lijXa x yYb x y 1001002020njLn ijstx ijlllll lijstnjLn yijst ijlllll lijsta x yXx yab x yXx yb10011001njLnL ijstst ijlllllll lijlnjLnL ijstst ijlllllll lijla x yx yX x yb x yx yY x y 常用的重采样方法有： ①最近邻法（Nearest Neighbor） ②双线性内插法（Bilinear Interpolation） ③三次卷积法（Cubic Convolution） 数据准备 （1）已经纠正的 Spot 全色波段影像：bldr_sp.img （2）待纠正的 TM 多光谱影像：bldr_tm.img 以距内插点最近的栅格点灰度值为对应像元 的灰度值。
x=Integer(x+0.5)，y=Integer(y+0.5) 优点：简单易用，计算量小，几何位置上 的精度±0.5 像元，不损害原来的像元值 缺点：纠正后的图像亮度不连续，会影响 制图效果，图斑边缘会出现锯齿现象 任意位于 4 个像素 fi,j，fi,j+1，fi+1,j，fi+1,j+1之间 的值 p(x,y)，由双线性插值得到为 p=(1-dx)(1-dy) fi,j+dx(1-y)fi,j+1 +dy(1-x)fi+1,j+dxdyfi+1,j+1 式中，dx=x-int(x)，dy=y-int(y) 优点：该方法改变了原影像的数据值，与 最近邻法相比较，虽计算量增加，但精度明显 提高，特别是对亮度不连续现象或线状特征的 块状化现象有明显的改善 缺点：此方法具有低通滤波的性质，会损 失图像中的一些边缘或线性信息，导致图像模 糊 tttwsin)(使用内插点周围 16 个观测点的像元值， 采用 3 次卷积函数（辛克函数）的离散化形式对所求 像元值进行内插 优点：该方法改变了原来的数据，具有影 像的均衡化和清晰化的效果，可得到平滑程度 较高的影像，适合遥感制图。
缺点：计算量大 操作方法及过程 （1）在 Display#1 视窗中导入待纠正的 TM 多光谱影像：bldr_tm.img，以 4、3、2 合 成显示；在 Display#2 视窗中导入已经纠正的 Spot 全色波段影像：bldr_sp.img，以灰 度值显示 （2）在菜单栏中点击 Map/Registration/Select GCPs:Image to Image 弹出对话框，在 Base Image 中选择 Display#2，Warp Image 中选择 Display#1，后点击 OK （3）弹出对话框，点击 Options/Set Point Colors，弹出对话框， 可以对控制点的颜色进行修改， 以便于在影像图上辨认 （4）利用放大、缩小及漫游工具,在 Display#1 和 Display#2 中选择控制点（按照实 习原理中的控制点选择原则） ， Display#1 和 Display#2 中一对控制点必须是同一种地物， 当选择好一对控制点时，点击 add （5）刚开始时可以在影像工作区范围的四个拐角进行选择，这样有利于控制误差；当 控制点的个数足够多的时候，可以将 Degree 改为 2，在 Display#2 中选择一个控制点， 利用 Predict,可以在 Display#1 中大概预测出该控制点对应控制点，然后进行适当的修改即可；点击 Show List 即可弹出对话框，里面显示控制点的坐标及单点的误差，可以将误差较大的控制点删除。
（6）Display#1 中一个控制点对应 Display#2 中一个控制点，共选择 25 对控制点调 整控制点位置，使得单点误差小于 1 个像元，总体误差小于 0.5 个像元 （7）保存控制点设置：点击对话框下的 File/Save GCPs to ASCII，弹出对话框，利用 Choose 选择储存路径和文件名，点击OK （8） 重采样并导出几何精纠正后的遥感影像图： 点击对话框下的 Options/Warp File(as Image to Map)弹出对话框，从中选择bldr_tm.img，点击 OK,弹出对话框如下图所示，在 X Pixel Size 填 30m，Y Pixel Size 填 30m（原因是待纠正的 TM 影像分辨率是 30*30m），在 Method 中选择 Polynomial（多 项式纠正模型），Degree 中选择 2（二阶），Resampling 中选择 Bilinear（双线性内插 法），利用 Choose 选择输出路径和文件名，点击 OK 即可 （9）将几何精纠正后的影像与参考影像 bldr_sp.img 进行卷帘操作，验证几何精纠正 的效果：打开 ENVI Zoom，点击 File/Open，选择几何精纠正后的影像和参考影像bldr_sp.img，点击“打开”，点击工具栏中的卷帘工具，进行验证。
结果与分析 （1）利用纠正好的全色波段高分辨率影像完成同景多光谱影像的几何精纠正： 图 1 图 2 图 3 图 4 图 1 为需几何精纠正的图像，图 2 为同一地区参考图像，两图像中对应的各均匀分 布 25 个控制点，图 3 为 GCP 表，从该表中我们可以清楚地了解 25 个控制点中单点误差 最小的是 GCP#14，值为 0.0705，单点误差最大的是 GCP#3，值为 0.8305，剩下各点中大 部分点的单点误差在 0.1 到 0.5 之间，控制点总体误差为 0.371200，误差较小 （2）操作结果： 图 5 用几何精纠正后的图像与参考图像做卷帘操作，图像在居民区无明显的跳跃，而在水库的边上有稍微的跳跃，在山区也稍微有点跳跃，尤其在山中的路表现得较为突出 分析原因：首先在山区没有高程数据，在 Google 地球中发现该影像的山区高差较大，这 使得选择控制点时很难保证在同一高程，纠正后山区稍有一些变形；其次两个图像的分 辨率不一样，这使得控制点的选择比较困难，尤其 TM 影像的分辨率为 30*30m 识别地物 较为困难；水库的跳跃可能是季节的变化导致水库边界不一致，出现了跳跃。
（3）采用一阶多项式和二阶多项式进行影像几何精纠正的差异： 在几何精纠正中，最少控制点的个数：L=(n+1)(n+2)/2，n 为多项式阶数采用一阶 多项式进行影像几何精纠正，至少需要 3 个控制点，采用二阶多项式进行影像几何精纠 正，至少需要 6 个控制点在适当增加控制点数量时，采用二阶多项式进行影像几何精 纠正比一阶效果好 （4）使用 Google 地球找到亚特兰大市（Atlanta） ，用校正过的影像与 Google 地球 上的地图相比： 发现在居民区两者地物位置基本一致，但是在山区的路中，总是稍有一点偏差，如 图所示，这也印证了步骤 2 所出现的现象 存在问题与解决办法 （ 1 ） 在 实 习 的 过 程 中 遇 到 有 的 同 学 在 进 行 显 示 图 像 时 ， 如 图 所 示 选 择然后 Load Band 后，能显示影像；但是当控制点选完后，再进行几何精纠正时，在 Options 下没有 Warp File(as Image to Map)选项，不知道原因， 查资料未能找到解释，在下次实习课上，向老师请教 （2）选择的控制点单点误差太大 首先应尽量寻找明显的地物特征点，如地形地物交叉点、河流弯曲或分叉处、泊边 缘、城廓边缘、水坝和交叉路口等能准确定位的特征点,而且在特征变化大的地区应该多 选取，并且将参考影。