空间后方交会程序.docVIP

空间后方交会程序 空间后方交会程序 PAGE 第 PAGE 8 页 共 NUMPAGES 8 页 空间后方交会程序 实验目的： 掌握摄影测量空间后方交会的原理，利用计算机编程语言实现空间后方交会外方位元素的解算。 仪器用具及已知数据文件： 计算机windows xp 系统，编程软件（VISUAL C++），地面控制点在摄影测量坐标系中的坐标及其像点坐标文件。 实验内容： 单张影像的空间后方交会：利用已知地面控制点数据及相应像点坐标根据共线方程反求影像的外方位元素。 数学模型：共线条件方程式： 求解过程： （1）获取已知数据。从航摄资料中查取平均航高与摄影机主距；获取控制点的地面测量坐标并转换为地面摄影测量坐标。 （2）量测控制点的像点坐标并做系统改正。 （3）确定未知数的初始值。在竖直摄影且地面控制点大致分布均匀的情况下，按如下方法确定初始值，即： ,, φ =ω=κ=0 式中；m为摄影比例尺分母；n为控制点个数。 （4）用三个角元素的初始值，计算个方向余弦，组成旋转矩阵R。 （5）逐点计算像点坐标的近似值。利用未知数的近似值和控制点的地面坐标代入共线方程式，逐点计算像点坐标的近似值（x）、（y）。 （6）逐点计算误差方程式的系数和常数项，组成误差方程式。 （7）计算法方程的系数矩阵和常数项，组成法方程式。 （8）解法方程，求得外方位元素的改正数,，，dφ,dω,dκ。 （9）用前次迭代取得的近似值，加本次迭代的改正数，计算外方位元素的新值。 （10）将求得的外方位元素改正数与规定的限差比较，若小于限差则迭代结束。否则用新的近似值重复（4）～（9），直到满足要求为止。 四．实验结果： 程序的源代码如下所示： #include #include #include #include #include #define N 4 void turn(double *A,double A2[],int m,int n) n); return; } for(i=0;iam;i++) for(j=0;jbn;j++) { u=i*bn+j; C[u]=; for(l=0;lan;l++) C[u]+=A[i*an+l]*B[l*bn+j]; } return; } double *inv(double *a,int n) f,M[i*n+j]); printf(\n); } printf(\n); } main() lf,Ys=%.5lf,Zs=%.5lf,t=%.5lf,w=%.5lf,k=%.5lf\n,Xs0,Ys0,Zs0,t,w,k); } printf(\n外方位元素为：\n); printf(Xs=%.5lf,Ys=%.5lf,Zs=%.5lf,t=%.5lf,w=%.5lf,k=%.5lf\n,Xs0,Ys0,Zs0,t,w,k); fclose(fp); } 程序运行的结果为： 五．实验总结： 通过这次的实验我学到了很多的东西，通过编程加深了对摄影测量空间后方交会相关知识的理解。在老师的严格要求下，我翻阅了很多的C语言的书籍，看了很多的算法，这个程序最难的地方是矩阵的求逆，以及如何将文本文档中的数据赋给程序中的数组变量，这些在以前的学习中都没有学到，感谢刘昭华老师的严格要求。