httpreadpudnc0md0wnl0ads74s0urcec0dewind0wscsharp269472wdqkjhfjh空间后方交会实习报告d0c.docVIP

遥感信息工程学院 04023班 王丹青 200432590117 空间后方交会实习报告 1 实习目的： 用Visual C++编写一个完整的单片空间后方交会程序，通过对提供的试验数据进行计算，输出像片的外方位元素并评定精度。深入理解单片空间后方交会的原理，体会在有多余观测情况下，用最小二乘平差方法编程实现解求影像外方位元素的过程。通过上机调试程序加强动手能力的培养，通过对实验结果的分析，增强综合运用所学知识解决实际问题的能力。 2实习环境： 2.1硬件环境：window操作系统 2.2软件环境：vc++6.0 3实习内容： 利用一定数量的地面控制点，根据共线条件方程求解像片外方位元素： 4 数据准备： 4.1 已知航摄仪内方位元素f＝153.24mm，Xo＝Yo＝0。 4.2 已知4对点的影像坐标和地面坐标：  影像坐标地面坐标x（mm）y（mm）X(m)Y（m）Z(m)1-86.15-68.9936589.4125273.322195.172-53.4082.2137631.0831324.51728.693-14.78-76.6339100.9724934.982386.50410.4664.4340426.5430319.81757.31 5实习过程： 5.1学习单张像片空间后方交会的基本理论，掌握其基本思想。 如果我们知道每幅影像的6个外方位元素，就能确定被摄物体与航摄影像的关系。而单像空间后方交会就是用于测定像片的外方位元素的，它的基本思想是：以单幅影像为基础，从影像所覆盖的地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发，根据共线方程，解求该影像在航空摄影时刻的外方位元素Xs,Ys,Zs,p,w,k.由于空间后方交会所采用的数学模型共线方程是非线性函数，为了便于外方位元素的解求，首先将其线性化。 5.2在纸上绘出空间后方交会的计算机程序框图。 为了能够在宏观上指导我们编写程序，我们需要在草稿纸上绘出程序框图。框图如下： 输入原始数据 ↓ 归算像点坐标x,y ↓ 计算和确定初值Xs0,Ys0,Zs0,p0,w0,k0 ↓ ︱ 组成旋转矩阵R ︱ ↓ ︱ 计算(x),(y)和lx,ly 迭 ↓ 迭 逐点组成误差方程式并法化 次 ↓ 数 否 所有点完否？ 小 ↓完 于 解法方程，求未知数改正数 限 ↓ 差 计算改正后的外方位元素 否 否 ↓ ？ 未知数改正数限差否？ 否︱ ↓是 ︱ 整理并输出计算结果 输出中间结果和出错信息 ↓ ︱ 正常结束 非正常结束 6. 按照程序框图编写程序。 程序代码如下： //#include stdafx.h #include math.h #include stdlib.h #include stdio.h #define q 206265 //矩阵相乘 void mult(double *m1,double *m2,double *result,int m,int p,int n) { int i,j,k; for(i=0;im;i++) for(j=0;jn;j++) { result[i*n+j]=0.0; for(k=0;kp;k++) result[i*n+j]+=m1[i*p+k]*m2[j+k*n]; }