作業(yè)4空間后方交會(已修改)
2025-04-05 07:03 本頁面
【正文】 攝影測量作業(yè)報告之空間后方交會 陳聞亞 20080729 作業(yè)報告空間后方交會專 業(yè)： 測繪工程 班 級： 2008級(1)班 姓 名： 陳聞亞 指導教師： 陳強 2010 年 4 月 16 日 1 作業(yè)任務 32 作業(yè)思想 33 作業(yè)條件及數(shù)據(jù) 34 作業(yè)過程 35 源程序 46 計算結果 177心得體會與建議 17 1 作業(yè)任務計算近似垂直攝影情況下后方交會解。即利用攝影測量空間后方交會的方法，獲取相片的6個外方位元素。2作業(yè)思想利用攝影測量空間后方交會的方法求解。該方法的基本思想是利用至少三個一直地面控制點的坐標A（XA，YA，ZA）、B（XB，YB，ZB）C（XC，YC，ZC），與其影像上對應的三個像點的影像坐標a（xa，ya）、b（xb，yb）、c（xc，yc），根據(jù)共線方程，反求該相片的外方位元素XS、YS、ZS、φ、ω、κ。3作業(yè)條件及數(shù)據(jù)已知攝影機主距f=，四對點的像點坐標與相應的地面坐標列入下表：表1點號像點坐標地面坐標x（mm）y(mm)X(m)Y(m)Z(m)12.344作業(yè)過程 4．1 獲取已知數(shù)據(jù) 相片比例尺1/m=1：10000，內(nèi)方位元素f=，x0，y0；獲取控制點的地面測量坐標Xt、Yt、Zt。4．2 量測控制點的像點坐標： 本次作業(yè)中為已知。見表1。4．3 確定未知數(shù)的初始值： 在近似垂直攝影情況下，膠原素的初始值為0，即φ0 = ω0 = κ0=0；線元素中，ZS0=H=mf=，XS0、YS0的取值可用四個控制點坐標的平均值，即：XS0= =YS0= =4．4 計算旋轉矩陣R： 利用膠原素的近似值計算方向余弦值，組成R陣。4．5 逐點計算像點坐標的近似值： 利用未知數(shù)的近似值按共線方程式計算控制點像點坐標的近似值（x）（y）。4．6 組成誤差方程： 逐點計算誤差方程式的系數(shù)和常數(shù)項。4．7 組成法方程式： 計算法方程的系數(shù)矩陣ATA與常數(shù)項ATL。4．8 求解外方位元素： 根據(jù)法方程，由式X=（AtA）1 ATL解求外方位元素改正數(shù)，并與相應的近似值求和，得到外方位元素新的近似值。4．9 求解外方位元素： 將求得的外方位元素的改正數(shù)與規(guī)定的限差（）比較，小于限差則計算終止，知道滿足要求為止。5 源程序include include include const double PRECISION=1e5。typedef double DOUBLE[5]。int InputData(int amp。Num, DOUBLE *amp。Data,double amp。m,double amp。f)。int Resection(const int amp。Num,const DOUBLE *amp。Data,const double amp。m,const double amp。f)。int InverseMatrix(double *matrix,const int amp。row)。int main(int argc, char* argv[]){DOUBLE *Data=NULL。int Num。double f(0),m(0)。if(InputData(Num,Data,m,f)){ if (Data!=NULL) { delete []Data。 } return 1。}if(Resection(Num,Data,m,f)){ if (Data!=NULL) { delete []Data。 } return 1。}if (Data!=NULL){ delete []Data。}printf(解算完畢...\n)。do{ printf(計算結果保存于\\文件中\(zhòng)n 請選擇操作（輸入P打開結果數(shù)據(jù)，R打開原始數(shù)據(jù)，其它退出程序）：)。 fflush(stdin)。 //刷新輸入流 char order=getchar()。 if (39。P39。==order || 39。p39。==order) { system()。 } else if (39。R39。==order || 39。r39。==or
教學教案相關推薦
鄂ICP備17016276號-1