ENVI实习之几何精纠正.docxVIP

内容 以具有地理参考的全色波段影像为基础，对没有地理坐标的多光谱影像进行几何校正。 步骤与要求 打开并显示图像文件，启动ENVI几何校正模块进行几何精纠正操作 采用二阶多项式的影像对影像纠正法，从纠正好的全色波段影像中读取对应点的坐标控制点。 控制点数不少于20个，并均匀分布在图幅的各个方向上。要求单点定位误差不得超过0.5个像元大小（RMS error ≤0.5）。 记录并分析GCP表中的数据，给出精度评价结果。 选择options菜单中的Warp File，重采样生成纠正多光谱影像文件 目视检测纠正精度。在全色影像主窗口中选择Tools | swipe打开卷帘窗口，查看纠正后影像与纠正好的影像，结合效果，检验精纠正效果。 基本原理由于非系统性因素（传感器姿态，高度，地球曲率，地形等因素）造成影像发生单元的畸变，整景影像的畸变，使用几何精纠正可以消除这种畸变 几何精纠正就是利用一系列的GCP和由其所构建的几何纠正模型完成对整景影像的纠正，使其单元格大小一致等，此外，还添加了投影信息（和参考影像投影信息一致） 几何精纠正精度保证是：几何精纠正模型，这个模型一般选取多项式模型，根据GCP来确定多项式模型的参数，完成模型构建。 之后，空间插值：利用逆向映射插值完成空间插值 之后，灰度插值，这里有三种方案：最邻近法，线性内插法，三次卷激法。这三种均可以完成灰度重采样，但是结果并不相同 最邻近法：选取据内插点最近的栅格点的灰度值作为内插点的灰度值，结果往往不连续 线性内插法和三次卷积法 采用内插点领域的结果作为内插点的灰度值 数据准备实习提供数据 某一地区的1994年的spot影像，分辨率为10米，utm投影 该地区2009年的TM影像，分辨率30米，无投影信息 操作方法及过程打开两幅影像（bldr-ap 和bldr-tm） open image map，使两幅影像显示 观察：TM影像没有投影信息 SP有投影信息（UTM） 打开 map—registration—select GCPS image to image模块，设置矫正模型。TM显示窗口为 warp map ，SP显示窗口为base map 进入了GCP选择对话框，既可以输入矫正坐标，也可以在两窗口中点取 先在四个角选取四个控制点，便于构建最简单的模型（3个点既可，但是为了保证精度，选取四个点） add point 完成三个点采集时，模型即可建立，此时就可以使用predict工具了。根据所建立的模型，可以预测base中点对应的warp中的大概位置。通过调节，便可以选取同名地。添加到列表中 添加第五个点后，ERROR x，y和rms便已经出现，说明模型基本建立 继续添加控制点，保证RMS控制在0.5内 完成20个点采集 但是，RMS=0.95 通过 UPDATE工具调节，直至单点RMS0.5 通过调节，单点RMS均小于0.5 完成模型构建 选择 warp fiel（对tm影像进行几何精纠正） 修改分辨率为30米（因为tm影像的分辨率为30m，纠正为10 没有意义，只是徒增了大小） Method 选取 polynomial 多项式，选取2项，gcp大于6 输出位置选择 在ENVI zoom中加载两幅影像使用SWIP进行对比 结果：在分界线右侧基本吻合，在右侧略有变形 结果与分析纠正后的影像具有了投影信息 utm 矫正后的影像左上角坐标和SP影像的左上角坐标不同，这是因为两幅影像的范围不同所致 矫正后的影像和TM原始影像对比发现，处在山区略有不同外，其他区域基本吻合 在山区，路是吻合的，只有山的角略有偏移 存在问题???解决办法 但完成三个点的拾取时，模型一阶多项式模型建立，这是就可以使用predict功能了，并且直接添加predict的点，这时，误差很小。 其实，这个模型是逐步建立的，首先三个点建立一个简单模型，使用predict预测这个模型中，点对应的内插点，我们必须利用这个预测的点，找到base中的点的同名点，找到后的这组点，参与模型构建 为什么会在第五个点添加时，才显示RMS，而不是第四个点？？ 总结几何精纠正 消除非系统几何畸变 给影像投影信息，将影像与实际地物对应起来 关于update X 负 向左移 正右移 Y 负 向上移 正下移 同坐标系 按方向键，移动距离0.5单元格