空間后方交會編程實習報告.doc

《空間后方交會編程實習報告.doc》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。

1、空間后方交會編程實習報告一實習目的用程序設(shè)計語言（VisualC++或者C語言）編寫一個完整的單片空間后方交會程序，通過對提供的試驗數(shù)據(jù)進行計算，輸出像片的外方位元素并評定精度。本實驗的目的在于讓學生深入理解單片空間后方交會的原理，體會在有多余觀測情況下，用最小二乘平差方法編程實現(xiàn)解求影像外方位元素的過程。通過上機調(diào)試程序加強動手能力的培養(yǎng)，通過對實驗結(jié)果的分析，增強學生綜合運用所學知識解決實際問題的能力。二實習內(nèi)容利用一定數(shù)量的地面控制點，根據(jù)共線條件方程求解像片外方位元素。三實習數(shù)據(jù)已知航攝儀的內(nèi)方位元素：fk=153.24mm，x0=y

2、0=0.0mm，攝影比例尺為1:50000；4個地面控制點的地面坐標及其對應像點的像片坐標：點號像片坐標(mm)地面點坐標(m)xyXYZ1-86.15-68.9936589.4125273.322195.172-53.4082.2137631.0831324.51728.693-14.78-76.6339100.9724934.982386.50410.4664.4340426.5430319.81757.31四實習原理如果我們知道每幅影像的6個外方位元素，就能確定被攝物體與航攝影像的關(guān)系。因此，如何獲取影像的外方位元素，一直是攝影測量工作

3、者所探討的問題?？刹扇〉姆椒ㄓ校豪美走_、全球定位系統(tǒng)（GPS）、慣性導航系統(tǒng)（INS）以及星相攝影機來獲取影像的外方位元素；也可以利用影像覆蓋范圍內(nèi)一定數(shù)量的控制點的空間坐標與攝影坐標，根據(jù)共線條件方程，反求該影像的外方位元素，這種方法稱為單幅影像的空間后方交會。單像空間后方交會的基本思想是：以單幅影像為基礎(chǔ)，從該影像所覆蓋地面范圍內(nèi)若干控制點的已知地面坐標和相應點的像坐標量測值出發(fā)，根據(jù)共線條件方程，解求該影像在航空攝影時刻的外方位元素Xs，Ys，Zs，t，w，k。五實習流程（1）獲取已知數(shù)據(jù)。從攝影資料中查取影像比例尺1/m，平均攝影距

4、離（航空攝影的航高、內(nèi)方位元素x0，y0，f；獲取控制點的空間坐標Xt，Yt，Zt。（2）量測控制點的像點坐標并進行必要的影像坐標系統(tǒng)誤差改正，得到像點坐標。（3）確定未知數(shù)的初始值。單像空間后方交會必須給出待定參數(shù)的初始值，在豎直航空攝影且地面控制點大體對稱分布的情況下，可按如下方法確定初始值：Z0s=H=m*f+ΣZi/4；X0s=Σxi/n;Y0s=ΣYi/n;t=ω=κ=0;式中：m為攝影比例尺分母；(4)計算旋轉(zhuǎn)矩陣R。利用角元素的近似值按下式計算方向余弦值a1,a2,a3,b1,b2,b3,c1,c2,c3，組成R陣。（5）逐點計

5、算像點坐標的近似值。利用未知數(shù)的近似值按共線條件方程計算控制點像點坐標的近似值（x）、（y）。（6）按下式逐點計算誤差方程式的系數(shù)和常數(shù)項，組成誤差方程。(7)計算法方程的系數(shù)矩陣ATA與常數(shù)項ATL，組成法方程；（8）解求外方位元素。根據(jù)法方程，解求外方位元素的改正數(shù)，并與相應的近似值求和，得到外方位元素新的近似值。（9）檢查計算是否收斂。將所求得的外方位元素的改正數(shù)與規(guī)定的限差比較，通常對t、ω、κ、Xs、Ys、Zs的改正數(shù)Δt，Δω，Δκ，ΔXs，ΔYs，ΔZs給予限差，當改正數(shù)小于限差時，迭代結(jié)束。否則用新的近似值重復（4）——（8）

6、步驟計算，直到滿足要求為止。（10）空間后方交會的精度估計：按上述方法所求得的影像外方位元素的精度可以通過法方程式中未知數(shù)的系數(shù)矩陣的逆陣（ATA）-1來解求，此時視像點坐標為等精度不相關(guān)觀測值。因為ATA）-1中第i個主對角線上的元素Qii就是法方程式中第i個未知數(shù)的權(quán)倒數(shù)，若單位權(quán)中誤差為m0，則第i個未知數(shù)的中誤差為：mi=0當參加空間后方交會的控制點有n個時，則單位權(quán)中誤差可按下式計算：六主要代碼與詳解voidR(doublet,doublew,doublek,double*a,double*b,double*c){a[0]=cos(

7、t)*cos(k)-sin(t)*sin(w)*sin(k);a[1]=-cos(t)*sin(k)-sin(t)*sin(w)*cos(k);a[2]=-sin(t)*cos(w);b[0]=cos(w)*sin(k);b[1]=cos(w)*cos(k);b[2]=-sin(w);c[0]=sin(t)*cos(k)+cos(t)*sin(w)*sin(k);c[1]=-sin(t)*sin(k)+cos(t)*sin(w)*cos(k);c[2]=cos(t)*cos(w);}//子函數(shù)計算旋轉(zhuǎn)矩陣R。利用角元素的近似值按下式計算方向余

8、弦值a1,a2,a3,b1,b2,b3,c1,c2,c3，組成R陣。voidmain(){inti,m,num;doublet,w,k,Xs,Ys,Zs,f;//六