高斯投影正反算及換帶計算VB程序設(shè)計

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

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

4、n2 option2 option3 option3 option4 option4 屬性 caption name caption name caption name caption name 值 六度帶 opt6du 三度帶 opt3du 正算 opttran 反算 optcon 5.2程序代碼設(shè)計 在這里主要介紹高斯投影坐標(biāo)轉(zhuǎn)換的正反算代碼設(shè)計，完整的代碼見附錄1所示。 5.2.1投影計算過程的正算子過程代碼設(shè)計 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投影計算過程的反算子過程代碼設(shè)計 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