地信遥感课设 - 图文

遥感技术课程设计报告

图2.5 点残差图 图2.4 部分几何校正图

9页 第

遥感技术课程设计报告

图2.7 变换函数图

（3）根据变换函数，取四个角求出整体图像的范围重复（2）输入新的行列数得到整体的矫正图像

图2.8 点残差图 图2.9 整体几何校正图

2.2.1 纠正变换函数的求取

先比较几个较可靠点的残差值，去掉残差值大的点，如果不可靠的点都去掉了，总的RMS值大于0.5，则开始比较可靠点的残差值并去掉残差值大的点，直到总的RMS值小于0.5。最后点确定扭则会出现纠正后的图像

说明控制点的取舍情况、写出最终纠正变换函数，及相应的均方误差。 旧坐标：左上角（0,742） 右上角（871,742） 左下角（0,0） 右下角(871,0)

ULX：232679.0 ULY：3394695.0 DRX：259319.0 DRY：3374367.0

列数=（Max X-Min X）/分辨率=（259319.0-232679.0）/24=1110 行数=（3394695.0-3374367.0）/24=847

a=12779.3366204466729,a1=0.0325947719709593,a2=-0.0060196114357260;

10页 第

遥感技术课程设计报告

b=-110487.1547799512628,b1=0.0063519315430515,b2=0.0322664407680667;

由x=fx(X，Y)=a0+a1X+a2Y y=fy(X，Y)=b0+b1X+b2Y得出

X=(a2*(y-b)-b2*(x-a))/(a2*b1-a1*b2); Y=(a1*(y-b)-b1*(x-a))/(a1*b2-a2*b1);

将上面四个点代进去求出新的四个点为左上角（235985,3400750） 右上角（261770,3395680） 左下角（231887,3378560） 右下角(257672,337349)

新坐标：MAXX=261770 MINX=231887 MAXY=3400750 MINY=3373490

列数=（Max X-Min X）/分辨率=（261770-231887）/24=1245

行数=（Max Y-Min Y）/分辨率=（3400750-3373490）/24=1135

2.2.2 纠正后图像边界范围的确定

在原始图像上测出左上角、右上角、左下角、右下角四个点的坐标，代入几何部分纠正的变换函数中，可求出新的X最大值、最小值，Y最大值、最小值再根据列数=（Max X-Min X）/分辨率， 行数=（Max Y-Min Y）/分辨率求出几何纠正图像的行列数，对原始图像进行整体的几何纠正得到整体几何纠正图像。 #include #include using namespace std; void main() {

double x,y; double

a=12779.3366204466729,a1=0.0325947719709593,a2=-0.0060196114357260; double

b=-110487.1547799512628,b1=0.0063519315430515,b2=0.0322664407680667; double X,Y;

cout<<\请输入原始坐标的x:\<>x;

cout<<\请输入原始坐标的y:\<>y;

X=(a2*(y-b)-b2*(x-a))/(a2*b1-a1*b2); Y=(a1*(y-b)-b1*(x-a))/(a1*b2-a2*b1); cout<<\<

}

11页 第

遥感技术课程设计报告

截图结果如下：

2.2.3纠正后图像与原始图像的对比分析

1.着重说明纠正前后地物的空间位置变化。

12页 第