空間后方交會程序空間后方交會程序空間后方交會程序xxx公司空間后方交會程序文件編號：文件日期：修訂次數：第1.0次更改批準審核制定方案設計，管理制度實驗目的：掌握攝影測量空間后方交會的原理，利用計算機編程語言實現空間后方交會外方位元素的解算。儀器用具及已知數據文件：計算機windowsxp系統，編程軟件（VISUALC++），地面控制點在攝影測量坐標系中的坐標及其像點坐標文件。實驗內容：單張影像的空間后方交會：利用已知地面控制點數據及相應像點坐標根據共線方程反求影像的外方位元素。數學模型：共線條件方程式：求解過程：（1）獲取已知數據。從航攝資料中查取平均航高與攝影機主距；獲取控制點的地面測量坐標并轉換為地面攝影測量坐標。（2）量測控制點的像點坐標并做系統改正。（3）確定未知數的初始值。在豎直攝影且地面控制點大致分布均勻的情況下，按如下方法確定初始值，即：,,φ=ω=κ=0式中；m為攝影比例尺分母；n為控制點個數。（4）用三個角元素的初始值，計算個方向余弦，組成旋轉矩陣R。（5）逐點計算像點坐標的近似值。利用未知數的近似值和控制點的地面坐標代入共線方程式，逐點計算像點坐標的近似值（x）、（y）。（6）逐點計算誤差方程式的系數和常數項，組成誤差方程式。（7）計算法方程的系數矩陣和常數項，組成法方程式。（8）解法方程，求得外方位元素的改正數,，，dφ,dω,dκ。（9）用前次迭代取得的近似值，加本次迭代的改正數，計算外方位元素的新值。（10）將求得的外方位元素改正數與規定的限差比較，若小于限差則迭代結束。否則用新的近似值重復（4）～（9），直到滿足要求為止。四．實驗結果：程序的源代碼如下所示：#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);}程序運行的結果為： 五．實驗總結：通過這次的實驗我學到了很多的東西，通過編程加深了對攝影測量空間后方交會相關知識的理解。在老師的嚴格要求下，我翻閱了很多的C語言的書籍，看了