ENVI下数字图像几何校正

ENVI下数字图像几何校正

http://blog.sina.com.cn/s/blog_54388b830100gqsg.html

ENVI下数字图像几何校正(影像地理坐标定位和配准) (2010-01-06 09:41:51)转载▼ 标签： rs rnvi 杂谈 分类： ENVI/IDL


几何校正就是从具有几何变形的图像中消除变形的过程。而将地图投影系统赋予图像数据的过程，称为地理参考（Geo-referencing）。由于几何校正的遥感影像都需要遵循一定的地图坐标系统，因此几何校正的过程包含了地理参考过程。

多项式纠正法（Polynomial），利用已知象元在影像上的位置（行列坐标）及其相应地理坐标（可利用地图获取）建立校正变换多项式，进而对整幅图像进行校正，是最常用的几何校正方法。基于多项式的几何校正关键是控制点的选区。

多项式纠正法的精度与地面控制点（GCPs）的精度、分布、数量及纠正范围有关，GCPs的位置精度越高，则几何纠正的精度越高；GCPs分布应尽可能在整幅图像内均匀分布，否则在GCPs密集区精度较高，在GCPs分布稀疏区出现较大误差。

实验数据：…\data2\cztRef1 ； …\data2\czt_raw1

实验步骤：

1.重新启动ENVI。

2.打开文件cztRef1 和czt_raw1。

3.查看地图投影信息。记录：czt_raw1，投影( 未知 ),椭球体参数（未知），分辨率（未知）；cztRef1，投影Projection( UTM),椭球体参数Datum（WGS—84），分辨率（28.5）。

4.显示图像。Available Bands List窗口：点击Display->New Display，利用czt_raw1合成彩色图像（R-5、G-4、B-3波段合成），点击Load RGB在图像窗口#1显示图像；点击Display #1->New Display，利用cztRef1合成彩色图像（R-5、G-4、B-3波段合成），点击Load RGB在图像窗口#2显示图像。

5.打开控制点选择窗口（Ground Controls Points Selection）。主菜单：Map-> Registration->Select GCPs: Image to Map;窗口Image to Map Registration中设置: Select Registration Projection-- （ UTM），点击Datum…，选择Available Datums --（WGS—84），点击OK按钮，Zone--（49），X Pixel Size-（28.5 ），Y Pixel Size-（28.5 ）。

6.打开#2 Pixel Locator。主影像窗口#2菜单：Tools-> Pixel Locator…。拖动主影像窗口#2中的红色框，观察 #2 Zoom窗口中的变化和 #2 Pixel Locator中位置数据的变化。

7.选择控制点。将两个Zoom窗口中的“十”字定位到同一地物对应位置，点击#2 Pixel Locator中的Export按钮，点击窗口Ground Controls Points Selection的Add Point按钮，增加控制点。用同样方法选择10个以上控制点。如果刚输入的控制点有误，可点击Delete Last Point 按钮删除。也可点击Show List查看、修改控制点。

8.检查控制点。改变Degree为1，查看记录：选取到的控制点数量有（11）个，最大均方根误差RMS Error（0.4662），总体均方根误差RMS Error（0.3719）；改变Degree为2，查看记录：选取到的控制点数量有（11）个，最大均方根误差RMS Error（0.3684），总体均方根误差RMS Error（0.2349）。可见当Degree增大时RMS Error（减小）。（总体均方误差要自己计算的吧！）

9.保存控制点数据。窗口Ground Controls Points Selection菜单：File->Save GCPs w/ Coords…。指定文件名，保存控制点数据。拷贝反映控制点分布的Scroll图像，如图（）。拷贝反映控制点信息的Image to Map GCP List 窗口图像，如图（）。

10.校正图像。窗口Ground Controls Points Selection菜单：Options->Warp File…；选择文件czt_raw1，点击OK。窗口Registration Parameters:改变Resampling为三次卷积内插（Cubic Convolution），点击Choose按钮，输入文件名，点击打开。窗口Registration Parameters:点击OK。

11.查看校正结果图像的地图投影信息。记录：校正结果图像文件名（xiaozhengjieguo），投影(UTM ),椭球体参数（WGS-84 ），分辨率（28.5m ）。

12.显示校正结果图像。关闭除主菜单窗口和Available Bands List窗口以外的所有窗口，按步骤4中指定的方法，显示校正结果图像（窗口#1）和图像cztRef1（窗口#2）。拷贝校正结果图像的Scroll图像，如图（）。

检查校正结果精度。影像窗口#1：Tools->Link->Geographic Link…（on）（on）;在两个Zoom窗口中打开“十”字定位器；在影像的不同位置拖动主影像窗口#1中的红色框，观察 #1 Zoom窗口中“十”字定位器位置和 #2 Zoom中窗口“十”字定位器位置是否对应同一地物，如果偏差较大说明校正结果存在较大的误差，分析其原因，重新选取控制点进行校正。(貌似结果还可以吧，偏差不是很大。)