【正文】 for(i=0。 for(j=0。coutX[i][j]endl。couti+1。cout= X[i][j] 。 } else if(j==1) { coutY。couti+1。j3。i4。 } } cout已知地面四個點的坐標為：\n。 if(j==0) { coutxi+1= setw(8)x[i][j] 。j2。i4。 cout已知像點坐標為：\n。 double sumXYZ[3]={0}。 double result2[6]={0}。//定義常數(shù)項 const double pi=。//定義了一個系數(shù)陣 double AT[6][8]={0}。//用于存放像點估計值。//像片比例尺為1:40000 double R[3][3]={0}。//設定未知數(shù)(XS,YS,ZS,三個角度)初始值 double f=。 int i,j,num=0。int main(){ double x[4][2]={,}。templatetypename T1,typename T2void transpose (T1*mat1,T2*mat2,int a,int b)。includefstreaminclude iomanipconst int n=6?？傊兄x老師給了我這次鍛煉自己、提升自己的機會。在編程過程中為了確保計算的準確性，我還利用matlab同步編寫了一個相同功能的計算單張相片空間后方交會的迭代程序，并以此來驗算原c++程序的正確性。在本次實驗中，我遇到了很多困難，但是都一一得到了解決，在解決困難的過程中的編程能力得到了提高，對其所涉及到的知識的印象也得到了加深。11. 由誤差方程的計算式：以及單位權中誤差的計算式（式中：表示多余觀測數(shù)）以及平差中6個參數(shù)的協(xié)因數(shù)陣最終得到參數(shù)的中誤差為三、 算法流程否是否是輸入原始數(shù)據(jù)像點坐標計算，系統(tǒng)誤差改正確定外方位元素初始值組成旋轉(zhuǎn)矩陣R逐點組成誤差方程式并法化所有像點完否計算改正后的外方位元素計算中誤差，輸出成果，結(jié)束解法方程，求外方位元素改正數(shù)外方位元素改正數(shù)是否小于限差是結(jié)束并提示錯誤信息迭代次數(shù)小于n四、 源程序源程序代碼請見附頁五、 計算結(jié)果在經(jīng)過三次迭代之后得到的最終成果如表1所示表1，最終計算結(jié)果Xs（米）Ys（米）Zs（米）ψ（弧度）ω（弧度）κ（弧度）六、 結(jié)果分析由計算結(jié)果可知在拍攝照片瞬間，攝影中心在地面攝影測量坐標系中的坐標為（，）（單位：米）。10. 將求得的外方位元素改正數(shù)與規(guī)定的限差比較，若小于限差，則迭代結(jié)束。9. 用前次迭代取得的近似值，加本次迭代的改正數(shù)，計算外方位元素的新值。由常數(shù)項計算公式：得到常數(shù)項矩陣計算式為：7. 計算法方程的是系數(shù)矩陣和常數(shù)項，組成法方程式。在豎直攝影且地面控制點大體對稱分布的情況下，按如下方法確定初始值，即式中：為攝影比例尺分母；為控制點個數(shù) 4. 用三個角元素的初始值按下式計算各方向余弦值，組成旋轉(zhuǎn)矩陣矩陣中各元素的計算公式如下：5. 逐點計算像點坐標的近似值。2. 測量控制點的像點坐標并作系統(tǒng)誤差改正。二、 計算原理1. 獲取已知數(shù)據(jù)。單像空間后方交會姓名： 學號： 時間： Echo did this for you .2013/4/25目錄一、 作業(yè)任務 3 二、 計算原理 3 三、 算法流程 7 四、 源程序 8 五、 計算結(jié)果 8 六、 結(jié)果分析 8 七、 心得與體會 8 八、 附頁 8 1. c++程序 8 2. C++程序截圖 15 3. matlb程序 16 一、 作業(yè)任務已知條件：攝影機主距f=，x0=0，y0=0, 像片比例尺為1:40000，有四對點的像點坐標與相應的地面坐標如下表。 點號像點坐標地面坐標x(mm)y(mm)X(m)Y(m)Z(m)1234以單像空間后方交會方法，求解該像片的外方位元素。從航攝資料中查取平均航高與攝影機主距；獲取控制點的地面測量坐標并轉(zhuǎn)換為地面攝測坐標。3. 確定未知數(shù)的初始值。利用未知數(shù)的近似值和控制點的地面坐標，帶入以下共線方程式， 逐點計算像點坐標的近似值、6. 逐點計算誤差方程式的系數(shù)和常數(shù)項，組成誤差方程式。8. 解法方程，求得外方位元素的改正數(shù)。式中：代表迭代次數(shù)。否則用新的近似值重復4~9滿足要求為止。表2，精度評定結(jié)果　參數(shù)XsYsZsψωκ中誤差 七、 心得與體會通過本次實驗，我對單張像片的空間后方交會的計算原理及實現(xiàn)過程有了很深的認識，并在牢記各種計算公式的推導過程基礎上做到了熟練應用。這個程序本身還有一些不足，例如:在迭代過程中的矩陣的相乘、輸出等，沒有采用編寫函數(shù)的形式，而是逐一進行計算，這是程序整體看起來較臃腫。現(xiàn)在編寫完這個程序，心中很自豪，也有一些遺憾，因為我想用c來編寫，因為c在暑期實習用過之后，時隔近一年在沒有復習，不免把學會的知識又還給了老師，本想通過這次實驗對c進行復習，但是由于時間原因以及一些個人因素，沒能達到自己原本的目標。八、 附頁1. c++程