下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認(rèn)領(lǐng)
文檔簡介
本文格式為Word版,下載可任意編輯——攝影測量試驗報告(空間后方交會—前方交會)空間后方交會-空間前方交會程序編程試驗一.試驗?zāi)康囊?/p>
把握運用空間后方交會-空間前方交會求解地面點的空間位置。學(xué)會運用空間后方交會的原理,根據(jù)所給控制點的地面攝影測量坐標(biāo)系坐標(biāo)以及相應(yīng)的像平面坐標(biāo)系中的坐標(biāo),利用計算機編程語言實現(xiàn)空間后方交會的過程,完成所給像對中兩張像片各自的外方位元素的求解。然后根據(jù)空間后方交會所得的兩張像片的內(nèi)外方位元素,利用同名像點在左右像片上的坐標(biāo),求解其對應(yīng)的地面點在攝影測量坐標(biāo)系中的坐標(biāo),并完成精度評定過程,利用計算機編程語言實現(xiàn)此過程。
二.儀器用具
計算機、編程軟件(MATLAB)
三.試驗數(shù)據(jù)
試驗數(shù)據(jù)包含四個地面控制點(GCP)的地面攝影測量坐標(biāo)及在左右像片中的像平面坐標(biāo)。此四對坐標(biāo)運用最小二乘法求解左右像片的外方位元素,即完成了空間后方的過程。另外還給出了5對地面點在左右像片中的像平面坐標(biāo)和左右像片的內(nèi)方位元素。試驗數(shù)據(jù)如下:
左片點號xGCP1GCP2GCP3GCP41234516.01288.5613.36282.2451.75814.61849.8886.24348.135y79.96381.134-79.37-80.02781.555-0.231-0.792-1.346-79.962x-73.93-5.252y78.70678.184X5083.2055780.025210.8795909.264Y5852.0995906.3654258.4464314.283Z527.925571.549461.81455.484右片地面攝影測量坐標(biāo)-79.122-78.879-9.887-39.953-76.016-42.201-7.706-80.08978.4630.036-1.022-2.112-44.438-79.736內(nèi)方位元素:f=152.000mm,x0=0,y0=0
四.試驗框圖
輸入GCP的像點坐標(biāo)xy
確定初始值ψ=ω=κ=0,Xs,Ys,Zs
計算旋轉(zhuǎn)矩陣R計算像點在像空間坐標(biāo)系中的近似值(x),(y),并組成誤差方程的常數(shù)項l計算誤差方程的系數(shù)項組成系數(shù)矩陣A組成法方程式,計算系數(shù)A’A常數(shù)項A’L求解外方位元素計算ψ、ω、κ、Xs、Ys、Zs改正后的值
計算完畢
此過程完成空間后方交會求解像片的外方位元素,其中改正數(shù)小于限差(0.00003,相當(dāng)于0.1’的角度值)為止。在這個過程中采用迭代的方法,是外方位元素逐漸收斂于理論值,每次迭代所得的改正數(shù)都應(yīng)加到上一次的初始值之中。
改正數(shù)是否小大于于限差?小于
輸入所需計算點的像平面坐標(biāo)x1,y1;x2,y2
根據(jù)后方交會所得的旋轉(zhuǎn)矩陣R1,R2計算像點在左右像空間輔助坐標(biāo)系中的坐標(biāo)X1Y1Z1,X2Y2Z2
計算攝影基線的三個坐標(biāo)分量BxByBz
計算個點在左右像片中的的投影系數(shù)N1N2
計算地面所求點在地面攝影測量坐標(biāo)系中的坐標(biāo)XAYAZA
計算完畢,精度評定
在空間后方交會中運用的數(shù)學(xué)模型為共線方程
確定Xs,Ys,Zs的初始值時,對于左片可取地面左邊兩個GCP的坐標(biāo)的平均值作為左片Xs和Ys的初始值,取右邊兩個GCP的坐標(biāo)平均值作為右片Xs和Ys的初始值。Zs可取地面所有GCP的Z坐標(biāo)的平均值再加上航高??臻g前方交會的數(shù)學(xué)模型為:
五.試驗源代碼
functionMain_KJQHFJH()globalRg1g2mGacb1b2;m=10000;a=5;c=4;
feval(@shuru);%調(diào)用shuru()shurujcp()函數(shù)完成像點及feval(@shurujcp);ìP有關(guān)數(shù)據(jù)的輸入XYZ=feval(@MQZqianfangjh);%調(diào)用MQZqianfangjh()函數(shù)完成空間前方、%%%%%%單位權(quán)中誤差%%%%%后方交會計算解得外方位元素
globalV1V2;%由于以上三個函數(shù)定義在外部文件中故需VV=[];%用feval()完成調(diào)用過程fori=1:2*c
VV(i)=V1(i);VV(2*i+1)=V2(i);end
m0=sqrt(VV*(VV')/(2*c-6));
disp('單位權(quán)中誤差m0為正負(fù):');disp(m0);%計算單位權(quán)中誤差并將其輸出顯示
輸入GCP像點坐標(biāo)及地面攝影測量坐標(biāo)系坐標(biāo)的函數(shù)和輸入所求點像點坐標(biāo)函數(shù):
functionshurujcp()globalcm;
m=input('攝影比例尺:');%輸入GCP像點坐標(biāo)數(shù)據(jù)函數(shù)并分別將其c=input('GCP的總數(shù)=');%存入到不同的矩陣之中disp('GCP左片像框標(biāo)坐標(biāo):');globalg1;g1=zeros(c,2);i=1;
whilei<=c
m=input('x=');n=input('y=');
g1(i,1)=m;g1(i,2)=n;i=i+1;end
disp('GCP右片像框標(biāo)坐標(biāo):');globalg2;g2=zeros(c,2);i=1;
whilei<=c
m=input('x=');n=input('y=');
g2(i,1)=m;g2(i,2)=n;i=i+1;end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%functionshuru()globala;
a=input('計算總像對點數(shù)=');%完成想計算所需的像平面坐標(biāo)globalb1;%坐標(biāo)輸入,存入不同的矩陣中b1=zeros(a,2);
disp('左片像點坐標(biāo):')i=1;
whilei<=a
m=input('x=');n=input('y=');
b1(i,1)=m;b1(i,2)=n;i=i+1;end%%
globalb2;b2=zeros(a,2);
disp('右片像點坐標(biāo):')i=1;
whilei<=a
m=input('x=');n=input('y=');
b2(i,1)=m;b2(i,2)=n;i=i+1;end%%globalc;
c=input('GCP的總數(shù)=');
disp('GCP攝影測量系坐標(biāo):')globalG;G=zeros(3,c);i=1;
whilei<=c
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 醫(yī)用支撐床墊項目可行性實施報告
- 醫(yī)院用推車項目可行性實施報告
- 內(nèi)衣相關(guān)項目實施方案
- 冰箱貼裝飾磁鐵相關(guān)項目實施方案
- 西醫(yī)綜合分類模擬題284
- 臨床醫(yī)學(xué)綜合能力(西醫(yī))分類模擬題584
- 卷管機玻璃加工機械項目可行性實施報告
- 辦公室用座椅市場環(huán)境與對策分析
- 帶有內(nèi)窺鏡的便攜式挖耳勺商業(yè)發(fā)展計劃書
- 兒童帽市場環(huán)境與對策分析
- 四人的劇本殺
- 燃?xì)鈽I(yè)務(wù)代辦授權(quán)書模板
- 巖石的強度理論及破壞判據(jù)
- 現(xiàn)代設(shè)計史知識點(最終)要點
- 李志斌《五猖會》pptPPT通用課件
- 遵紀(jì)守法主題班會課件共16張PPT(共16張PPT)
- 徐州工程學(xué)院專升本(英語)科目考試題庫(含歷年重點題)
- 2022年北京市海淀初二英語期中試卷
- 多功能清障車工作裝置及液壓系統(tǒng)設(shè)計
- 春考語文必背
- 項目2 “博客”頁面制作
評論
0/150
提交評論