文档介绍:摄影测量实习报告-摄影测量实习报告64
一、实习任务利用自己所熟悉的一种编程语言,实现单像空间后方交会,解求此张像片的6个外方位元素, , , ,ω,κ,摄影测量实习报告。二、实习目的1、深刻理解单张像片空间后方交会的原理与意义;2、在存在多余观测值时,利用最小二乘平差方法,经过迭代,求的外方位元素的最佳值;3、熟悉vc编程方法,利用编程实现计算。三、实习原理以单幅影像为基础,从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发,根据共线条件方程,求解该影象在航空摄影时刻的像片外方位元素, , , ,ω,κ共线条件方程如下:x-x0=-f*/y-y0=-f*/其中:x,y为像点的像平面坐标; x0,y0,f为影像的外方位元素;, ,为摄站点的物方空间坐标;x,y,z为物方点的物方空间坐标;旋转矩阵r为;由于此共线条件方程是非线性方程,先对其进行线性化,利用泰勒展开得:=(x)-x++++++++=(y)-y++++++++像点观测值一般视为等权,即p=i;矩阵形式:v=ax-l,p=i;通过间接平差,为提高精度,增加多余观测方程,根据最小二乘平差原理,可计算出外方位元素的改正数。
经过迭代计算,每次迭代用未知数的近似值与上次迭代计算的改正数之和作为新的近似值,重复计算,求出新的改正数,这样反复趋近,直到改正数小于某个限值为止。四、程序框图输入原始数据归算像点坐标x,y 计算并确定初值, , , ,组成旋转矩阵r 计算和逐点组成误差方程式并法化所有点完否? 解法方程,求未知数改正数计算改正后的外方位元素未知数改正数摄影测量实习报告》。经过比较发现,计算出的6个外方位元素与所给参考值相比,相差很小,计算结果符合要求:线元素误差小于米;角元素误差30秒。计算其精度,可以通过法方程式中未知数的系数矩阵的逆阵-1来求解,此时,视像点坐标为等精度不相关观测值。因为逆阵中第i个主对角线上元素qii就是法方程式中第i个未知数的权倒数,若单位权中误差为m0,则第i个未知数的中误差为:mi=当参加空间后方交会的控制点有n个时,则单位权中误差可按下式计算:m0=要求:线元素精度mx等,高于米;角元素精度高于弧度。
计算结果都达到标准。在此次计算中,我运用了所给的全部控制点,而空间后方交会所运用的控制点,应该避免位于一个圆柱面上,否则会出现解不唯一
的情况。选点时