攝影測量后方交會VC實現(xiàn)代碼及實習(xí)報告

1、遙感信息工程學(xué)院 09031班 2009302590123 吳煜暉攝影測量課間實習(xí)單片空間后方交會班級：09031姓名：吳煜暉學(xué)號：20093025901232011-10-8一、實習(xí)原理單像空間后方交會的基本思想是：以單幅影像為基礎(chǔ)，從該影像所覆蓋地面范圍內(nèi)若干控制點的已知地面坐標(biāo)和相應(yīng)點的像坐標(biāo)測值出發(fā)，根據(jù)共線條件方程，解求該影像在航空攝影時刻的外方位元素XS，YS，ZS，。由于空間后方交會所采用的數(shù)學(xué)模型共線方程是非線性函數(shù)，為了便于外方位元素的解求，須首先對共線方程進行線性化，然后進行后方交會，最后在精度評定。二、實習(xí)過程1、實習(xí)所用數(shù)據(jù)本次實習(xí)數(shù)據(jù)采用課本P.44 27題所給數(shù)據(jù)。

2、如下圖：2、程序流程圖及界面設(shè)計本程序程序框圖如下：輸入原始數(shù)據(jù)歸算像點坐標(biāo)x,y計算和確定初值XS0，YS0，ZS0，0，0，0組成旋轉(zhuǎn)矩陣計算(x)，(y)和lx，ly逐點組成誤差方程式并法化所有點完否？解法方程，求未知數(shù)改正數(shù)計算改正后的外方位元素未知數(shù)改正數(shù)限差否？整理并輸出計算結(jié)果正常結(jié)束迭代次數(shù)小于限差否？輸出中間結(jié)果和出錯信息非正常結(jié)束本人采用Visual C+6.0編此程序，利用MFC來設(shè)計程序。主程序頁面設(shè)計如下：子窗口（即進行計算后所得結(jié)果頁面）設(shè)計如下：3、程序代碼本程序代碼較多，在此講部分重要代碼列出，其余代碼參見程序源代碼。對話框類頭文件內(nèi)聲明的類成員及函數(shù)（后來增加

3、的）：public:void houfangjiaohui();double fi,w,k; /影像外方位角double Xs0,Ys0,Zs0; /后方交會所求解double a3,b3,c3;double mx6,m0;void inv(double *a,int n);/*正定矩陣求逆*/void transpose(double *m1, double *m2, int m, int n); /矩陣轉(zhuǎn)置void mult(double *m1, double *m2, double *result, int i_1, int j_12, int j_2);/矩陣相乘部分核心代碼，在此講

4、后方交會計算代碼給出：void CFinalworkDlg:houfangjiaohui()double t_fk;fi=w=k=0.0;int i; /中間變量double t_x4,t_y4,t_X4,t_Y4,t_Z4;/將單位換成米t_x0=m_x1/1000.0; t_x1=m_x2/1000.0; t_x2=m_x3/1000.0; t_x3=m_x4/1000.0;t_y0=m_y1/1000.0; t_y1=m_y2/1000.0; t_y2=m_y3/1000.0; t_y3=m_y4/1000.0;t_X0=m_X1; t_X1=m_X2; t_X2=m_X3; t_X3=

5、m_X4;t_Y0=m_Y1; t_Y1=m_Y2; t_Y2=m_Y3; t_Y3=m_Y4;t_Z0=m_Z1; t_Z1=m_Z2; t_Z2=m_Z3; t_Z3=m_Z4;Xs0=Ys0=Zs0=0.0;for(i=0;i4;i+)Xs0+=t_Xi;Ys0+=t_Yi;/確定未知數(shù)初始值Xs0/=4;Ys0/=4;t_fk=m_fk/1000.0;Zs0=t_fk*m_mk; double x00=m_x00;double y00=m_y00;double A8*6;double AT6*8;double ATA6*6;double L8;double ATL6*1;double

6、Xo4,Yo4,Zo4,Xom,Yom,Zom;double Delta6;while(1)a0=cos(fi)*cos(k)-sin(fi)*sin(w)*sin(k);a1=-cos(fi)*sin(k)-sin(fi)*sin(w)*cos(k);a2=-sin(fi)*cos(w);b0=cos(w)*sin(k);b1=cos(w)*cos(k);b2=-sin(w);c0=sin(fi)*cos(k)+cos(fi)*sin(w)*sin(k);c1=-sin(fi)*sin(k)+cos(fi)*sin(w)*cos(k);c2=cos(fi)*cos(w);for(i=0;i0

7、.00001|fabs(Delta1)0.00001|fabs(Delta2)0.00001|fabs(Delta3)0.00001|fabs(Delta4)0.00001|fabs(Delta5)0.00001)=0)double V8; double s=0.0; mult(A,Delta,V,8,6,1); for(i=0;i8;i+) Vi-=Li; s+=Vi*Vi; m0=sqrt(s/2.0);/單位權(quán)中誤差 for(i=0;i6;i+)/精度評定 mxi=sqrt(ATA7*i)*m0; break;三、實習(xí)結(jié)果程序主界面如下：在該界面中，我們可以手動輸入各項數(shù)據(jù)，也可以清空所有數(shù)據(jù)以及將所有空恢復(fù)默認數(shù)