文档介绍：地球科修鸟探境工程修陀实验报告书一实****任务在LPS中采集4个控制点及两个检查点的像平面坐标及其对应物方坐标;编写空间后方前方交会的程序,利用该程序计算出相片的外方位元素,并II利用内外方位元素解算出两个检查点的物方坐标,并与LPS工作站上的对应坐标相比较。二实验原理前方交会数学模型及公式共线方程xx_f3(*人-+一片丿+勺匕人-Zj一厂一勺(X「Xs)+£(N—£)+C3(Z「Zs)a2(Xa-X$)+乞(Ya一5)+C?(ZA—Zs)a3(X「Xs)+daA-3)+C3(ZA-Zs)旋转矩阵Rq=cos0cosk—sin0sincosinkci2=-cossiiix;-sin(/)sillcocoskcl=_sin0cos<yl\=cos(osink=coscocosk代=-sin。q=sin())cosk+cos(/)sincosiiikc)=一sin0sink+cos(/)sincocoskc3=cos0cos0(3)经过线性化,得到兀,y的误差方程式+^)2 +。\4卜0+勺厶厶。+Q'e'K+(工)—*vy=a2iA^,+幺卫△气+幺匕门‘+ +%氏+(x)_y矩阵形式如下:V=AX系数方程其中:勺=K/+勺刃/乙"13=[uJ4-C^Xj/Za21=[a2f+a3刃/Za22=lb2f+b3y]/Za”=ysin0一{x/f[兀cosk一ysink]+fcosk}cos©ay=一/sink—兀/f[xsinK+ycosK]ai6=Ja24=一兀sinfi?-[y/f(^cos/c-ysinr)-tsink]cos(oa2s=-/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[Zy-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) ) 根据左右像片的外方位元素计算摄影基线分量Bx,By,Bzo3) 逐点计算像点的空间辅助坐标。4) 计算投影系数。5) 计算未知点的地面摄影测量坐标。6) 重复以上步骤完成所有点的地面坐标的计算。程序框图后方交会程序框图五计算成果C:\Users\Administrat0r\De5ld0p\©---.**********.5像点坐标一y--:-•戎正数:50•952秒橡H筵角戚正数:--5••0022625260538*-1冰鬣黏数数盘顷黑正正正2<向向舟改改第辱像xsvsZS秒秒秒■r-x—vr-x128471?211••-3-2^^.•徴落卽|:歙数数戏顷星正正正心向向事改改第B鬲畑lysES改正粼小于6秒,迭代完成!迭代伏数为:3遂弍外方位元素的解题向顷角q:-«::-?5?:-: