摄影测量学空间前方后方交会.docxVIP

地球科修鸟探境工程修陀 实验报告书 一实习任务 在LPS中采集4个控制点及两个检查点的像平面坐标及其对应物方坐标；编写空间后方前方 交会的程序，利用该程序计算出相片的外方位元素，并II利用内外方位元素解算出两个检查 点的物方坐标，并与LPS工作站上的对应坐标相比较。 二实验原理 前方交会数学模型及公式 共线方程 x x _ f 3（*人-+ 一 片丿 + 勺匕人-Zj 一厂一勺（X「Xs） + ￡（N—￡） + C3（Z「Zs） a2 （Xa - X$ ） + 乞（Ya 一 5 ） + C? （ZA — Zs ） a3（X「Xs） + daA-3） + C3（ZA-Zs） 旋转矩阵R q = cos 0 cos k — sin 0 sin co sin k ci2 = - cos siii x; - sin（/）sill co cos k cl =_sin0cosy l\ = cos（o sin k = cos co cos k 代=-sin。 q = sin（））cos k + cos（/）sin co siii k c） = 一 sin 0 sin k + cos（/）sin co cos k c3 = cos 0 cos 0 （3）经过线性化，得到兀，y的误差方程式 + ^）2 +。\4卜0 + 勺厶厶。+ QeK + （工）—* vy = a2iA^, + 幺卫△气 + 幺匕门‘ + + %氏 + （x） _ y 矩阵形式如下： V = AX 系数方程 其中： 勺=K/ +勺刃/乙 13=[uJ 4-C^Xj/Z a21 =[a2f+ a3 刃/Z a22=lb2f+b3y]/Z a” = y sin 0 一 {x/ f [兀cosk 一 y sin k] + f cosk}cos? ay =一 / sink— 兀/ f [xsinK + ycosK] ai6=J a24 = 一兀sinfi?-[y/ f （^cos/c-ysin r）- t sin k]cos（o a2s =- / cos^-7 /（xsinxs-ycosA；） a26 = 一兀 —[（g,（y）-歼 x=a1（x_xs）+b1（y_y5）+cI（z_z5） y=a2（x_x5）+b2（Y_r5）+c2（z_z5） Z =勺（X 一 X$ ） + % （F 一岭）+ C3 （Z 一 Z$ ） 近似值计算公式如下 = JX [Z y-y0 = -fY/Z 由最小二乘间接平差原理可得： X = [△“ s b^s Afi? △?] xs = x/+aa：s1+a^2+- 3=F+?+△+??? Zs = Zs。+ +aZ$2 + … 0 = 0°+ △才 + + … 6? = d° + Ad?1 + Ad)2 + … K = K° +△/? +△? +… 后方交会数学模型，公式 在空间后方交会屮运用的数学模型为共线方程 〃i(x-x$)+zd(y-ys)+dw) -a3(X-Xs) + b3(y-?) + c3(Z-Zs) a2(X - Xs) + b2(Y-Ys) + c2(Z-Zs) -a3(X-Xs) + /?3(y-W) + e3(Z-Z$) 计算时使用迭代计算 附源代码 三实验思路及步骤 利用后方交会得出两张像片各自的外方位元素 1） 获取已知数据：从摄影资料中插曲像片比例尺、平均航高、内方位元素以及控制点的地 面摄影测量坐标及对应的像点坐标。 2） 确定未知数的初始值：在竖直摄彫的情况下，胶原素的初始值 为0,线元素其中Zs二m* 3） 计算旋转矩阵R。 4） 逐点计算像点坐标的近似值：利用共线方稈。 5） 组成误差方程并法化。 6） 解求外方位元素。 7） 检查计算是否收敛。 利用解求出的外方位元素进行前方交会 1） 用各自像片的角元素计算汕左右像片的旋转矩阵R1和R2. 2） 根据左右像片的外方位元素计算摄影基线分量Bx,By, Bzo 3） 逐点计算像点的空间辅助坐标。 4） 计算投影系数。 5） 计算未知点的地面摄影测量坐标。 6） 重复以上步骤完成所有点的地面坐标的计算。 程序框图 后方交会程序框图 五计算成果 C:\Users\Administrat0r\De5ld0p\ ?h.exe 序号 1 X -86.15 -53.4 -14.78 10.46 序号 1 X 36589.4 37631.139101 40426.5 像点坐标一 y -68.99 82.21-76.63 64.43 地窗座标一 V 25273.3 31324.5 24935 30319.8 Z 2195.17 728.69 2386.5 757.31 第/迭代 脳向顷角q改正数：-170.552秒 鲁向他角?戎正数：50 ? 952秒 橡H筵角戚正数：-