图3.6 选择多项式次数对话框
Polynomial Order 表示多项式纠正采用的次数,一般选择2次多项式。
在Projection标签页下单击Add/Change Projection按钮,按如下图设置投影:
点击OK,关闭多项式属性对话框。出现地面控制点(GCP)工具参考设置对话框 选择Existing Viewer,点击OK选择参考图像,点击OK,启动GCP TOOL。屏幕显示包括两个放大的视窗Viewer3,Viewer4,分别是两个主视窗Viewer1,Viewer2中两个关联方框的放大框。以及GCP TOOL:包含菜单和工具图标,图标下面显示两幅图像坐标数据,开始选取控制点。连接框对应的放大框。
现在Viewer1中选择明显点 ,然后Viewer2中的选择同名点 。
图3.11 Viewer1、Viewer2同名点
根据多项式纠正的要求,同名点对至少要在7对以上,重复1,2两个步骤,直到选择的控制点数满足纠正的数量要求。当选择了六对同名点后,第七个及以上同名点的选择:只需在Viewer1中选择一个明显点,其同名点会在Viewer1中自动显示,如果位置有误差,可以用鼠标左键调整点位,确保选择的是同名点。
系统自动计算得到单点误差和其贡献。
全部控制点平差后的精度显示在右上角,只有当Total小于1,才满足纠正精度要求,继续后面的重采样。
多项式纠正模型参数对话框选择Transformation则可以看到计算的变换参数
图3.17 显示多项式中a,b的值
保存相应的input和reference控制点文件,进行灰度重采样 在几何纠正工具(Geo Correction Tools)中
图3.18 几何纠正工具图标框
选择灰度重采样模块,出现重采样对话框(图2.19),输入纠正后的文件名。
图3.19 重采样对话框
在对话框里(Output File)输入纠正以后的图像文件名,该对话框显示了输出图像的大小行数(Number rows)列数(Number Columns),输出图像左上角的地面坐标(LUX,ULY)以及右下角的地面坐标(LRX,LRY)。输出像元的地面大小(Output Cell Sizes),确定参数后点击OK,重采样完成。
六、实验结果分析
在计算参数时已经利用检查点做了验证,现将纠正后的TM图像与作为参考图像的SPOT图像关联起来,通过人工比较来验证结果。在两个窗口中分别打开SPOT 和TM图像(图2.13)。
图3.13 已校正的TM和SPOT图像
然后点击鼠标右键,弹出菜单,选择Geo. Link/unlink,表示连接两幅图像,提示在要连接的窗口中点击鼠标,将两个窗口关联起来。点击视窗按扭条中的+,移动光标检查两边的位置是否对应。如果两边对应,说明纠正图像满足要求。
对纠正可能出现的结果要能正确分析。
如果总体精度不能满足纠正要求,其原因可能有: (i) (ii) (iii) (iv) (v)
多项式次数不能满足纠正要求;
控制点选择不符合要求:分布不均匀,或者选点精度不符合要求,甚至选错; 地形起伏引起的投影差超限; 参考坐标精度有问题;
模型有问题,需要用共线方程纠正。
四、辐射校正
大气校正 1.打开待校正图,进行辐射定标。
首先,点击Raster->Attributes在打开的对话框中统计出6个波段每个的到达1000像元以上的最小灰度值。
TM大气辐射校正G、B参数值 Band1 G 1.2688 B -0.01 Band2 2.98126 -0.0232 Band3 1.76186 -0.0078 Band4 2.81771 -0.0193 Band5 0.65227 -0.008 Band7 0.44375 -0.004 打开经过大气校正后的图像,与原图对比,点击图标格式如下。
打开直方图,单击选择植被,然后调整