高斯投影正反算及换带计算VB程序设计.doc

《高斯投影正反算及换带计算VB程序设计.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、摘要本设计主要阐述了高斯投影分带以及高斯投影坐标正、反算的推导公式，从而根据公式来编写基于VB语言基础上的换带及坐标转换程序。作者系统介绍了测量中经常使用的坐标系以及地图投影的概念和高斯投影的具体含义，叙述了换带和临带计算的原因以及它们在运算时的原理、过程，详细叙述了在VB语言中实现的原理基础以及代码的编写设计。在设计中根据高斯的正反算公式写出了基于VB语言的程序设计，其程序设计任务完成了由地理坐标向54平面坐标系和80平面坐标系转换的功能，以及由54坐标系和80坐标系向地理坐标系转换的功能，同时也有同一平面坐标系不同投影带之

2、间的换带计算和同一平面坐标系相同投影带临带计算等相互转换的功能。关键词：高斯投影、坐标正反算、换带计算、临带换算、程序设计5程序设计5.1界面设计本程序要实现的功能是根据所选择的椭球参数和指定的分带情况，将已知地理坐标或高斯投影坐标经正算和反算求得相应的高斯坐标和地理坐标，以及相应的换带计算和临带计算。因此需要用一个框架控件来组织椭球参数、两个框架分别组织分带选择和换算方式选择，两个框架组织地理坐标和高斯坐标，三28个命令按钮分别执行投影计算、换带和临带计算。程序设计界面如图5-1[9]图5-1高斯投影计算程序设计界面命令按钮

3、属性设置表如表5-1表5-1命令按钮属性设置表对象属性值Command1CaptionBL->xyCommand1NamecmdCalcCommand2Caption6->3Command2NamecmdChangeCommand3Caption临带计算Command3NamecmdNear28选择椭球框架内控件的属性值表5-2表5-2择椭球框架内控件的属性值对象属性值Text1Text6378245Text1NametxtaText2Text298.3Text2Nametxtalfa单选按钮控件属性设置表5-35-3单选按钮

4、控件属性设置表对象属性值Option1Caption六度带Option1Nameopt6duOption2Caption三度带Option2NameOpt3duOption3Caption正算Option3NameoptTranOption4Caption反算Option4NameoptCon5.2程序代码设计在这里主要介绍高斯投影坐标转换的正反算代码设计，完整的代码见附录1所示。5.2.1投影计算过程的正算子过程代码设计28①54系高斯投影正算子过程PublicSubPro54()Dimll#,N#,a0#,a4#,a6#,

5、a3#,a5#,cosB#cosB=Cos(B)ll=L-DoToHu(L0)N=6399698.902-(21562.267-(108.973-0.612*cosB*cosB)*cosB*cosB)*cosB*cosBa0=32140.404-(135.3302-(0.7092-0.004*cosB*cosB)*cosB*cosB)*cosB*cosBa4=(0.25+0.00252*cosB*cosB)*cosB*cosB-0.04166a6=(0.166*cosB*cosB-0.084)*cosB*cosBa3=(0.3

6、333333+0.001123*cosB*cosB)*cosB*cosB-0.1666667a5=0.0083-(0.1667-(0.1968+0.004*cosB*cosB)*cosB*cosB)*cosB*cosBX=6367558.4969*B-(a0-(0.5+(a4+a6*ll*ll)*ll*ll)*ll*ll*N)*Sin(B)*cosBY=(1+(a3+a5*ll*ll)*ll*ll)*ll*N*cosBEndSub②80系高斯投影正算子过程PublicSubPro80()Dimll#,N#,a0#,a4#,a6

7、#,a3#,a5#,cosB#cosB=Cos(B)ll=L-DoToHu(L0)N=6399596.652-(21565.045-(108.996-0.603*cosB*cosB)*cosB*cosB)*cosB*cosBa0=32144.5189-(135.3646-(0.7034-0.0041*cosB*cosB)*cosB*28cosB)*cosB*cosBa4=(0.25+0.00253*cosB*cosB)*cosB*cosB-0.04167a6=(0.167*cosB*cosB-0.083)*cosB*cosBa

8、3=(0.3333333+0.001123*cosB*cosB)*cosB*cosB-0.1666667a5=0.00878-(0.1702-0.20382*cosB*cosB)*cosB*cosBX=6367452.1328*B-(a0-(0.5+(a4+a6*ll*ll)