空间后方交会编程实习报告.doc

《空间后方交会编程实习报告.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、空间后方交会编程实习报告一实习目的用程序设计语言（VisualC++或者C语言）编写一个完整的单片空间后方交会程序，通过对提供的试验数据进行计算，输出像片的外方位元素并评定精度。本实验的目的在于让学生深入理解单片空间后方交会的原理，体会在有多余观测情况下，用最小二乘平差方法编程实现解求影像外方位元素的过程。通过上机调试程序加强动手能力的培养，通过对实验结果的分析，增强学生综合运用所学知识解决实际问题的能力。二实习内容利用一定数量的地面控制点，根据共线条件方程求解像片外方位元素。三实习数据已知航摄仪的内方位元素：fk=153.24mm，x0=y

2、0=0.0mm，摄影比例尺为1:50000；4个地面控制点的地面坐标及其对应像点的像片坐标：点号像片坐标(mm)地面点坐标(m)xyXYZ1-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四实习原理如果我们知道每幅影像的6个外方位元素，就能确定被摄物体与航摄影像的关系。因此，如何获取影像的外方位元素，一直是摄影测量工作

3、者所探讨的问题。可采取的方法有：利用雷达、全球定位系统（GPS）、惯性导航系统（INS）以及星相摄影机来获取影像的外方位元素；也可以利用影像覆盖范围内一定数量的控制点的空间坐标与摄影坐标，根据共线条件方程，反求该影像的外方位元素，这种方法称为单幅影像的空间后方交会。单像空间后方交会的基本思想是：以单幅影像为基础，从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发，根据共线条件方程，解求该影像在航空摄影时刻的外方位元素Xs，Ys，Zs，t，w，k。五实习流程（1）获取已知数据。从摄影资料中查取影像比例尺1/m，平均摄影距

4、离（航空摄影的航高、内方位元素x0，y0，f；获取控制点的空间坐标Xt，Yt，Zt。（2）量测控制点的像点坐标并进行必要的影像坐标系统误差改正，得到像点坐标。（3）确定未知数的初始值。单像空间后方交会必须给出待定参数的初始值，在竖直航空摄影且地面控制点大体对称分布的情况下，可按如下方法确定初始值：Z0s=H=m*f+ΣZi/4；X0s=Σxi/n;Y0s=ΣYi/n;t=ω=κ=0;式中：m为摄影比例尺分母；(4)计算旋转矩阵R。利用角元素的近似值按下式计算方向余弦值a1,a2,a3,b1,b2,b3,c1,c2,c3，组成R阵。（5）逐点计

5、算像点坐标的近似值。利用未知数的近似值按共线条件方程计算控制点像点坐标的近似值（x）、（y）。（6）按下式逐点计算误差方程式的系数和常数项，组成误差方程。(7)计算法方程的系数矩阵ATA与常数项ATL，组成法方程；（8）解求外方位元素。根据法方程，解求外方位元素的改正数，并与相应的近似值求和，得到外方位元素新的近似值。（9）检查计算是否收敛。将所求得的外方位元素的改正数与规定的限差比较，通常对t、ω、κ、Xs、Ys、Zs的改正数Δt，Δω，Δκ，ΔXs，ΔYs，ΔZs给予限差，当改正数小于限差时，迭代结束。否则用新的近似值重复（4）——（8）

6、步骤计算，直到满足要求为止。（10）空间后方交会的精度估计：按上述方法所求得的影像外方位元素的精度可以通过法方程式中未知数的系数矩阵的逆阵（ATA）-1来解求，此时视像点坐标为等精度不相关观测值。因为ATA）-1中第i个主对角线上的元素Qii就是法方程式中第i个未知数的权倒数，若单位权中误差为m0，则第i个未知数的中误差为：mi=0当参加空间后方交会的控制点有n个时，则单位权中误差可按下式计算：六主要代码与详解voidR(doublet,doublew,doublek,double*a,double*b,double*c){a[0]=cos(

7、t)*cos(k)-sin(t)*sin(w)*sin(k);a[1]=-cos(t)*sin(k)-sin(t)*sin(w)*cos(k);a[2]=-sin(t)*cos(w);b[0]=cos(w)*sin(k);b[1]=cos(w)*cos(k);b[2]=-sin(w);c[0]=sin(t)*cos(k)+cos(t)*sin(w)*sin(k);c[1]=-sin(t)*sin(k)+cos(t)*sin(w)*cos(k);c[2]=cos(t)*cos(w);}//子函数计算旋转矩阵R。利用角元素的近似值按下式计算方向余

8、弦值a1,a2,a3,b1,b2,b3,c1,c2,c3，组成R阵。voidmain(){inti,m,num;doublet,w,k,Xs,Ys,Zs,f;//六