空間后方交會報告.doc

《空間后方交會報告.doc》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、任務已知f=153.24mm，m=10000，限差0.1’各點坐標點號像點坐標地面坐標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求近似垂直攝影情況下后方交會解設計任務1、確定未知數(shù)的初始值：Φ0=ω0=К0=0,內(nèi)方位元素x0=0.0，y0=0.0，f=153.24mm。Zs0=H+14

2、i=14Zi=mf+14i=14Zi=3049.3m；Xs0=14i=14Xi=38437m；Ys0=14i=14Yi=27963.16m2、計算旋轉(zhuǎn)矩陣R利用角元素的近似值計算方向余弦值，組成R陣根據(jù)《攝影測量學》P32中的公式（3-9），初步計算R陣R[0][0]=cos(Φ)*cos(K)-sin(Φ)*sin(W)*sin(K);R[0][1]=-cos(Φ)*sin(K)-sin(Φ)*sin(W)*cos(K);R[0][2]=-sin(Φ)*cos(W);R[1][0]=cos(W)*sin(K);R[1][1]=cos(W)*cos(K);R[1

3、][2]=-sin(W);R[2][0]=sin(Φ)*cos(K)+cos(Φ)*sin(W)*sin(K);R[2][1]=-sin(Φ)*sin(K)+cos(Φ)*sin(W)*cos(K);R[2][2]=cos(Φ)*cos(W);得初始R陣1000100013、逐點計算近似值（x），（y）：帶入《攝影測量學》P61的公式（5-1）；得(x1)1(y1)1(x2)1(y2)1(x3)1(y3)1(x4)1(y4)1-0.33147-0.48258-0.053220.221960.15351-0.700980.133020.157564、組成誤差方程式

4、：按（5-8）；（5-9b）、（5-4）式逐點計算誤差方程式的系數(shù)和常數(shù)項根據(jù)Lx=x-(x)；Ly=y-(y)得Lx1Ly1Lx2Ly2Lx3Ly3Lx4Ly40.245320.41359-0.00018-0.13975-0.168290.62347-0.12256-0.09313解得A陣為-0.000100.000056219-0.201678-0.038785-0.068990-0.00010.0000450209-0.00000387855-0.18430.08615-0.000100.0000348473-0.171848-0.28648-0.0822

5、10-0.0001-0.0000536479-0.0000028648-0.1973440.0534534-0.000100.000009645-0.154666-0.00739096-0.076630-0.00010.0000500065-0.000000739096-0.191560.01478-0.00010-0.00000682589-0.153954-0.004397920.064430-0.0001-0.0000420452-0.000000439792-0.18033-0.01046根據(jù)公式X=（ATA）-1ATL解除X陣，檢驗Φ1、ω1、К1并不

6、在限差0.1’之內(nèi);利用解得的Φ1、ω1、К1，帶入2中循環(huán)；知道解出合格結(jié)果！1、結(jié)果，通過C++程序解得結(jié)果2、源程序代碼#include#include#include#include#includeusingnamespacestd;constintn=6;voidinverse(doublec[n][n]);templatevoidtranspose(T1*mat1,T2*mat2,inta,intb);template<

7、typenameT1,typenameT2>voidmulti(T1*mat1,T2*mat2,T2*result,inta,intb,intc);templatevoidinput(T*mat,inta,intb);templatevoidoutput(T*mat,char*s,inta,intb);intmain(){ofstreamoutFile;cout.precision(5);doublex0=0.0,y0=0.0;doublefk=0.15324;//內(nèi)方位元素doublem=10000;//估算比例尺d

8、oubleB[4][5]