高斯投影正反算及換帶計算VB程序設計

1、高斯投影正反算及換帶計算vb程序設計 摘要 本設計主要闡述了高斯投影分帶以及高斯投影坐標正、反算的推導公式，從而根據(jù)公式來編寫基于vb語言基礎上的換帶及坐標轉(zhuǎn)換程序。作者系統(tǒng)介紹了測量中經(jīng)常使用的坐標系以及地圖投影的概念和高斯投影的具體含義，敘述了換帶和臨帶計算的原因以及它們在運算時的原理、過程，詳細敘述了在vb語言中實現(xiàn)的原理基礎以及代碼的編寫設計。 在設計中根據(jù)高斯的正反算公式寫出了基于vb語言的程序設計，其程序設計任務完成了由地理坐標向54平面坐標系和80平面坐標系轉(zhuǎn)換的功能，以及由54坐標系和80坐標系向地理坐標系轉(zhuǎn)換的功能，同時也有同一平面坐標系不同投影帶之間的換帶計算和同一平面坐標

2、系相同投影帶臨帶計算等相互轉(zhuǎn)換的功能。 關鍵詞：高斯投影、坐標正反算、換帶計算、臨帶換算、程序設計 5 程序設計 5.1界面設計 本程序要實現(xiàn)的功能是根據(jù)所選擇的橢球參數(shù)和指定的分帶情況，將已知地理坐標或高斯投影坐標經(jīng)正算和反算求得相應的高斯坐標和地理坐標，以及相應的換帶計算和臨帶計算。因此需要用一個框架控件來組織橢球參數(shù)、兩個框架分別組織分帶選擇和換算方式選擇，兩個框架組織地理坐標和高斯坐標，三個命令按鈕分別執(zhí)行投影計算、換帶和臨帶計算。程 1 序設計界面如圖5-19 圖5-1 高斯投影計算程序設計界面 命令按鈕屬性設置表如表5-1 表5-1 命令按鈕屬性設置表 對象 command1 co

3、mmand1 command2 command2 command3 command3 屬性 caption name caption name caption name 2 值 bl-xy cmdcalc 6-3 cmdchange 臨帶計算 cmdnear 選擇橢球框架內(nèi)控件的屬性值表5-2 表5-2 擇橢球框架內(nèi)控件的屬性值 對象 text1 text1 text2 text2 屬性 text name text name 值 6378245 txta 298.3 txtalfa 單選按鈕控件屬性設置表5-3 5-3 單選按鈕控件屬性設置表 對象 option1 option1 optio

4、n2 option2 option3 option3 option4 option4 屬性 caption name caption name caption name caption name 值 六度帶 opt6du 三度帶 opt3du 正算 opttran 反算 optcon 5.2程序代碼設計 在這里主要介紹高斯投影坐標轉(zhuǎn)換的正反算代碼設計，完整的代碼見附錄1所示。 5.2.1投影計算過程的正算子過程代碼設計 3 54系高斯投影正算子過程 public sub pro54() dim ll#, n#, a0#, a4#, a6#, a3#, a5#, cosb# cosb = cos

5、(b) ll = l - dotohu(l0) n = 6399698.902 - (21562.267 - (108.973 - 0.612 * cosb * cosb) * cosb * cosb) * cosb * cosb a0 = 32140.404 - (135.3302 - (0.7092 - 0.004 * cosb * cosb) * cosb * cosb) * cosb * cosb a4 = (0.25 + 0.00252 * cosb * cosb) * cosb * cosb - 0.04166 a6 = (0.166 * cosb * cosb - 0.084)

6、* cosb * cosb a3 = (0.3333333 + 0.001123 * cosb * cosb) * cosb * cosb - 0.1666667 a5 = 0.0083 - (0.1667 - (0.1968 + 0.004 * cosb * cosb) * cosb * cosb) * cosb * cosb x = 6367558.4969 * b - (a0 - (0.5 + (a4 + a6 * ll * ll) * ll * ll) * ll * ll * n) * sin(b) * cosb y = (1 + (a3 + a5 * ll * ll) * ll *

7、ll) * ll * n * cosb end sub 80系高斯投影正算子過程 public sub pro80() dim ll#, n#, a0#, a4#, a6#, a3#, a5#, cosb# cosb = cos(b) ll = l - dotohu(l0) n = 6399596.652 - (21565.045 - (108.996 - 0.603 * cosb * cosb) * cosb * cosb) * cosb * cosb a0 = 32144.5189 - (135.3646 - (0.7034 - 0.0041 * cosb * cosb) * cosb *

8、 4 cosb) * cosb * cosb a4 = (0.25 + 0.00253 * cosb * cosb) * cosb * cosb - 0.04167 a6 = (0.167 * cosb * cosb - 0.083) * cosb * cosb a3 = (0.3333333 + 0.001123 * cosb * cosb) * cosb * cosb - 0.1666667 a5 = 0.00878 - (0.1702 - 0.20382 * cosb * cosb) * cosb * cosb x = 6367452.1328 * b - (a0 - (0.5 + (a

9、4 + a6 * ll * ll) * ll * ll) * ll * ll * n) * sin(b) * cosb y = (1 + (a3 + a5 * ll * ll) * ll * ll) * ll * n * cosb end sub 5.2.2投影計算過程的反算子過程代碼設計 54系高斯投影反算子過程12 public sub conpro54() dim bf#, bet#, z#, nf#, b2#, b3#, b4#, b5#, cos2b#, cos2bf# bet = x / 6367558.4969 cos2b = cos(bet) * cos(bet) bf = b

10、et + (50221746 + (293622 + (2350 + 22 * cos2b) * cos2b) * cos2b) * 0.0000000001 * sin(bet) * cos(bet) cos2bf = cos(bf) * cos(bf) nf = 6399698.902 - (21562.267 - (108.973 - 0.612 * cos2bf) * cos2bf) * cos2bf z = y / (nf * cos(bf) b2 = (0.5 + 0.003369 * cos2bf) * sin(bf) * cos(bf) b3 = 0.333333 - (0.166667 - 0.001123 * cos2bf) * co