资源描述:
《摄影测量学单像空间后方交会编程实习报告》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、摄影测量学单像空间后方交会编程实习报告班级:130x姓名:xx学号:2013302590xxx指导老师:李欣一、实习目的通过对提供的数据进行计算,输出像片的外方位元素并评定精度。深入理解单像空间后方交会的思想,体会在有多余观测情况下,用最小二乘平差方法编程实现解求影像外方位元素的过程。通过尝试编程实现加强编程处理问题的能力和对实习内容的理解,通过对实验结果的分析,增强综合运用所学知识解决实际问题的能力。了解摄影测量平差的基本过程,掌握空间后方交会的定义和实现算法。二、实习内容根据学习的单像空间后方交会的知识,用程序设计语言(C++或C语言)编写一个完整的单像空间后方
2、交会程序,通过对提供的数据进行计算,输出像片的外方位元素并评定精度。三、实习数据已知航摄仪的内方位元素:fk=153.24mm,x0=y0=0,摄影比例尺为1:15000;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四、实习原理
3、如果我们知道每幅影像的6个外方位元素,就能确定被摄物体与航摄影像的关系。因此,如何获取影像的外方位元素,一直是摄影测量工作者所探讨的问题。可采取的方法有:利用雷达、全球定位系统(GPS)、惯性导航系统(INS)以及星相摄影机来获取影像的外方位元素;也可以利用影像覆盖范围内一定数量的控制点的空间坐标与摄影坐标,根据共线条件方程,反求该影像的外方位元素,这种方法称为单幅影像的空间后方交会。 单像空间后方交会的基本思想是:以单幅影像为基础,从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发,根据共线条件方程,解求该影像在航空摄影时刻的外方位元素Xs
4、,Ys,Zs,ϕ,ω,κ。一、实习流程1.获取已知数据。从摄影资料中查取影像比例尺1/m,平均摄影距离(航空摄影的航高、内方位元素x0,y0,f;获取控制点的空间坐标Xt,Yt,Zt。 2.量测控制点的像点坐标并进行必要的影像坐标系统误差改正,得到像点坐标。3.确定未知数的初始值。单像空间后方交会必须给出待定参数的初始值,在竖直航空摄影且地面控制点大体对称分布的情况下,可按如下方法确定初始值:Zs0=H=m∙fXs0=1ni-1nXtiYs0=1ni-1nYtiφ0=ω0=0式中:m为摄影比例尺分母,n为控制点个数; 4.计算旋转矩阵R。利用角元素的近似值按公式计算
5、方向余弦值a1,a2,a3,b1,b2,b3,c1,c2,c3,组成R阵。 5.逐点计算像点坐标的近似值。利用未知数的近似值按共线条件方程计算控制点像点坐标的近似值(x),(y)。 6.逐点计算误差方程式的系数和常数项,组成误差方程。7.计算法方程的系数矩阵ATA与常数项ATL,组成法方程。 8.解求外方位元素。根据法方程,解求外方位元素的改正数,并与相应的近似值求和,得到外方位元素新的近似值。9.检查计算是否收敛。将所求得的外方位元素的改正数与规定的限差比较,通常对ϕ,ω,κ的改正数△φ,△ω,△κ给予限差,当改正数小于限差时,迭代结束。否则用新的近似值重复(4)
6、——(8)步骤计算,直到满足要求为止。 1.空间后方交会的精度估计: 按上述方法所求得的影像外方位元素的精度可以通过法方程式中未知数的系数矩阵的逆阵(ATA)-1来解求,此时视像点坐标为等精度不相关观测值。因为(ATA)-1中第i个主对角线上的元素Qii就是法方程式中第i个未知数的权倒数,若单位权中误差为m0,则第i个未知数的中误差为:mi=Qiim0当参加空间后方交会的控制点有n个时,则单位权中误差可按下式计算:m0=±vv2n-6流程图如下:一、程序代码#include#include"Matrix.h"//矩阵运算头文件,来自网络,已上传vo
7、idmain(){doubleXs,Ys,Zs,q,w,k;//外方位元素doublex0,y0,f;//内方位元素doublex[4],y[4];//影像坐标doubleX0[4],Y0[4],Z0[4];//地面坐标doublem;//比例尺doublea1,a2,a3,b1,b2,b3,c1,c2,c3;//旋转矩阵intn=0;//迭代次数//地面坐标X0[0]=36589.41;X0[1]=37631.08;X0[2]=39100.97;X0[3]=40426.54;Y0[0]=25273.32;Y0[1]=31324.51;Y0[2]=24934.9