空間后方交會(huì)程序

空間后方交會(huì)程序空間后方交會(huì)程序空間后方交會(huì)程序xxx公司空間后方交會(huì)程序文件編號(hào)：文件日期：修訂次數(shù)：第1.0次更改批準(zhǔn)審核制定方案設(shè)計(jì)，管理制度實(shí)驗(yàn)?zāi)康模赫莆諗z影測(cè)量空間后方交會(huì)的原理，利用計(jì)算機(jī)編程語(yǔ)言實(shí)現(xiàn)空間后方交會(huì)外方位元素的解算。儀器用具及已知數(shù)據(jù)文件：計(jì)算機(jī)windowsxp系統(tǒng)，編程軟件（VISUALC++），地面控制點(diǎn)在攝影測(cè)量坐標(biāo)系中的坐標(biāo)及其像點(diǎn)坐標(biāo)文件。實(shí)驗(yàn)內(nèi)容：?jiǎn)螐堄跋竦目臻g后方交會(huì)：利用已知地面控制點(diǎn)數(shù)據(jù)及相應(yīng)像點(diǎn)坐標(biāo)根據(jù)共線方程反求影像的外方位元素。數(shù)學(xué)模型：共線條件方程式：求解過(guò)程：（1）獲取已知數(shù)據(jù)。從航攝資料中查取平均航高與攝影機(jī)主距；獲取控制點(diǎn)的地面測(cè)量坐標(biāo)并轉(zhuǎn)換為地面攝影測(cè)量坐標(biāo)。（2）量測(cè)控制點(diǎn)的像點(diǎn)坐標(biāo)并做系統(tǒng)改正。（3）確定未知數(shù)的初始值。在豎直攝影且地面控制點(diǎn)大致分布均勻的情況下，按如下方法確定初始值，即：,,φ=ω=κ=0式中；m為攝影比例尺分母；n為控制點(diǎn)個(gè)數(shù)。（4）用三個(gè)角元素的初始值，計(jì)算個(gè)方向余弦，組成旋轉(zhuǎn)矩陣R。（5）逐點(diǎn)計(jì)算像點(diǎn)坐標(biāo)的近似值。利用未知數(shù)的近似值和控制點(diǎn)的地面坐標(biāo)代入共線方程式，逐點(diǎn)計(jì)算像點(diǎn)坐標(biāo)的近似值（x）、（y）。（6）逐點(diǎn)計(jì)算誤差方程式的系數(shù)和常數(shù)項(xiàng)，組成誤差方程式。（7）計(jì)算法方程的系數(shù)矩陣和常數(shù)項(xiàng)，組成法方程式。（8）解法方程，求得外方位元素的改正數(shù),，，dφ,dω,dκ。（9）用前次迭代取得的近似值，加本次迭代的改正數(shù)，計(jì)算外方位元素的新值。（10）將求得的外方位元素改正數(shù)與規(guī)定的限差比較，若小于限差則迭代結(jié)束。否則用新的近似值重復(fù)（4）～（9），直到滿(mǎn)足要求為止。四．實(shí)驗(yàn)結(jié)果：程序的源代碼如下所示：#include<>#include<>#include<>#include<>#include<>#defineN4voidturn(double*A,doubleA2[],intm,intn)n"); return; } for(i=0;i<am;i++) for(j=0;j<bn;j++) { u=i*bn+j; C[u]=; for(l=0;l<an;l++) C[u]+=A[i*an+l]*B[l*bn+j]; } return;}double*inv(double*a,intn)f",M[i*n+j]); printf("\n"); } printf("\n");}main()lf,Ys=%.5lf,Zs=%.5lf,t=%.5lf,w=%.5lf,k=%.5lf\n",Xs0,Ys0,Zs0,t,w,k);}printf("\n外方位元素為：\n");printf("Xs=%.5lf,Ys=%.5lf,Zs=%.5lf,t=%.5lf,w=%.5lf,k=%.5lf\n",Xs0,Ys0,Zs0,t,w,k);fclose(fp);}程序運(yùn)行的結(jié)果為： 五．實(shí)驗(yàn)總結(jié)：通過(guò)這次的實(shí)驗(yàn)我學(xué)到了很多的東西，通過(guò)編程加深了對(duì)攝影測(cè)量空間后方交會(huì)相關(guān)知識(shí)的理解。在老師的嚴(yán)格要求下，我翻閱了很多的C語(yǔ)言的書(shū)籍，看了