解析空中三角測量程序設計

解析空中三角測量程序設計

一、      目的要求

1、      通過程序的設計深入理解解析空中三角測量的整個過程
2、      提高應用程序設計語言解決問題的能力。

二、      資料及用具

一臺微機，一套調試程序所用的數據。

三、      實習內容及作業步驟

用VB（或者C語言）程序設計語言編寫單航帶區域網概算的程序，具有包含一些三方面的內容：
（一）      連續法像對的相對定向
  目標：求解各模型點在各模型像空間輔助坐標系中的坐標（Xi,Yi,Zi）
  步驟：
1、      用連續法相對定向求解相對定向元素（bx，by，bz，φ，ω，κ）。
2、      用前方交會公式計算各模型點坐標。
（二）      航帶網模型的建立
   目標：求出航帶內各模型點在航帶統一坐標系（第一個像片的像空間坐標系）中的坐標（Xi‘,Yi‘,Zi‘）。
   步驟：
1、      根據模型間的三個公共點求解個模型的縮放系數ki
2、      通過各模型縮放系數ki將模型內各模型點轉換到航帶統一坐標系。
（三）      模型的絕對定向
  目標：求出航帶內各模型點在大地坐標系中的坐標（Xti,Yti,Zti）。
  步驟：
1、      根據兩個已知控制點確定大地坐標系與地面攝測坐標系之間的轉換參數（a、b）。
2、      將所有控制點的大地坐標（Xt,Yt,Zt）通過轉換參數（a、b）轉換為地面攝測坐標系（Xtp,Ytp,Ztp）。
3、      根據已知控制點的地面攝測坐標（Xtp,Ytp,Ztp）求解絕對定向元素（ΔX，ΔY，ΔZ，Φ，Ω，Κ，λ）。
4、      根據絕對定向元素（ΔX，ΔY，ΔZ，Φ，Ω，Κ，λ）將各模型點在航帶統一坐標系中的坐標（Xi‘,Yi‘,Zi‘）轉換為地面攝測坐標系（Xtp,Ytp,Ztp）。
5、      將各模型點的地面攝測坐標（Xtp,Ytp,Ztp）通過轉換參數（a、b）轉換為所有控制點的大地坐標（Xt,Yt,Zt）。

四、      上交成果

調試完成的程序一份及計算結果。

附：調試數據：

1、      基本數據

攝影機主距：f＝153.033mm
            bx＝200mm

2、      像對坐標數據（單位：微米）

像對1：    701         702    
    1201      648790      735260     -230980      788550  
    1818      113860      595800     -771380      679190
    1202      241050      586260     -644640      661950
    1204      578030      223960     -327700      277590
    1203      256140      187820     -648360      260160
    1205      606230     -104550     -315660      -51440
    1206      278340     -565020     -660020     -487200
    1052      179220     -757800     -765430     -670400
     410      478510     -637940     -466930     -569840
     399      975930     -700850       16690     -658940
     192      917380      -16160       -4600       18540
     370      803150      758570      -76440      802960
    1118       94870      709670     -785850      795830
     194       35030      -34550     -877140       50930
-398       19790     -843710     -925660     -745370

像對2：     703         702    
     400     -568240     -631500      426790     -634450
    1207     -601170      642970      323920      637060
     399     -980300     -630600       16690     -658940
     192     -963100       51030       -4600       18540
     370     -989450      836700      -76410      802880
    1209       -8790      641530      921460      679700
     401      -29060     -915900      981740     -884160
-190        1460        5490      965780       38900

像對3：    703         704    
    1826      931930      729560      -26020      680090
     402      907100     -785750       20330     -836630
    1055      753660     -838360     -134160     -896670
     411      397770     -725220     -500070     -791380
    1211       49010       50160     -875490      -17320
     188      918060       48330      -10150       11950
    1225      890340      544420      -58540      499270
    1210      624000      444520     -317940      392160
    1209       -8690      641540     -945780      560530
     401      -29020     -915940     -930070    -1004070
-190        1460        5460     -921980      -63410

3、控制點數據

點號        X坐標          Y坐標            Z坐標
==================================================
 1201      24204.689      46604.652         46.251
 1205      24343.683      45111.263         48.198
 1206      24965.047      44309.253         49.340
 1209      22167.904      46432.515         46.457

4、檢查點數據

點號        X坐標          Y坐標          Z坐標
==================================================
 1118      25192.533      46608.059         48.936 
 1202      24941.046      46375.998         46.539
 1203      24945.705      45662.638         49.052
 1204      24369.486      45700.421         49.020
 1207      23218.292      46347.142         48.993