攝影測(cè)量學(xué)單像空間后方交會(huì)編程實(shí)習(xí)報(bào)告-文庫吧
2025-07-20 00:19 本頁面
【正文】 ni1nYtiφ0=ω0=0式中：m為攝影比例尺分母,n為控制點(diǎn)個(gè)數(shù)；4. 計(jì)算旋轉(zhuǎn)矩陣R。利用角元素的近似值按公式計(jì)算方向余弦值a1,a2,a3,b1,b2,b3,c1,c2,c3，組成R陣。5. 逐點(diǎn)計(jì)算像點(diǎn)坐標(biāo)的近似值。利用未知數(shù)的近似值按共線條件方程計(jì)算控制點(diǎn)像點(diǎn)坐標(biāo)的近似值(x),(y)。6. 逐點(diǎn)計(jì)算誤差方程式的系數(shù)和常數(shù)項(xiàng)，組成誤差方程。7. 計(jì)算法方程的系數(shù)矩陣ATA與常數(shù)項(xiàng)ATL，組成法方程。8. 解求外方位元素。根據(jù)法方程，解求外方位元素的改正數(shù)，并與相應(yīng)的近似值求和，得到外方位元素新的近似值。9. 檢查計(jì)算是否收斂。將所求得的外方位元素的改正數(shù)與規(guī)定的限差比較，通常對(duì)?,ω,κ的改正數(shù)△φ,△ω,△κ給予限差，當(dāng)改正數(shù)小于限差時(shí)，迭代結(jié)束。否則用新的近似值重復(fù)（4）——（8）步驟計(jì)算，直到滿足要求為止。10. 空間后方交會(huì)的精度估計(jì)：按上述方法所求得的影像外方位元素的精度可以通過法方程式中未知數(shù)的系數(shù)矩陣的逆陣（ATA）1來解求，此時(shí)視像點(diǎn)坐標(biāo)為等精度不相關(guān)觀測(cè)值。因?yàn)?ATA）1中第i個(gè)主對(duì)角線上的元素Qii就是法方程式中第i個(gè)未知數(shù)的權(quán)倒數(shù)，若單位權(quán)中誤差為m0，則第i個(gè)未知數(shù)的中誤差為：mi=Qiim0當(dāng)參加空間后方交會(huì)的控制點(diǎn)有n個(gè)時(shí)，則單位權(quán)中誤差可按下式計(jì)算：m0=177。vv2n6流程圖如下:六、 程序代碼include include //矩陣運(yùn)算頭文件，來自網(wǎng)絡(luò)，已上傳void main(){ double Xs,Ys,Zs,q,w,k。 //外方位元素 double x0,y0,f。 //內(nèi)方位元素 double x[4],y[4]。 //影像坐標(biāo) double X0[4],Y0[4],Z0[4]。 //地面坐標(biāo) double m。 //比例尺 double a1,a2,a3,b1,b2,b3,c1,c2,c3。 //旋轉(zhuǎn)矩陣 int n=0。 //迭代次數(shù) //地面坐標(biāo) X0[0]=。 X0[1]=。 X0[2]=。 X0[3]=。 Y0[0]=。 Y0[1]=。 Y0[2]=。 Y0[3]=。 Z0[0]=。 Z0[1]=。 Z0[2]=。 Z0[3]=。 //影像坐標(biāo) x[0]=。 x[1]=。 x[2]=。 x[3]=。 y[0]=。 y[1]=。 y[2]=。 y[3]=。 //確定
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號(hào)-1