基于EXCEL VBA平台抵偿坐标计算的实现.pdf

《基于EXCEL VBA平台抵偿坐标计算的实现.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第5期矿山测量NO．52014年10月MINESURVEYING0ct．2O14doi：10．3969／j．issn．1001—358X．2014．05．08{计算机应用}基于EXCELVBA平台抵偿坐标计算的实现赵淑湘(甘肃林业职业技术学院，甘肃天水741020)摘要：文中主要介绍了抵偿坐标的计算步骤，并在EXCELVBA平台下，通过开发自定义函数，实现了抵偿坐标计算的自动化、高效化和准确性，对城市坐标系的选择与建立具有一定的使用价值。关键词：EXCELVBA；抵偿坐标；实现中图分类号：P209文献标识标：B文章编号：1001—358X(2014)05—0

2、025—03在城市坐标系的选择与建设中，若采用国家统Xol(6)一的3。带高斯正形投影，则规定长度变形值应不大Y=Y+g(Y—ro)Jo于2．5cm／km，否则需要建立新的坐标系，该坐标系式中：X、卜国家统一坐标系的坐标；除长度变形值符合要求外，还应与国家统一坐标系X。、l，c——抵偿坐标系的坐标；之间能够相互转换。城市坐标系的建立方法较多，X。、Iio——缩放原点。可以选在测区中心的整其中抵偿坐标是最常用的一种，其原理是选择一个数值或控制点。和椭球体面相平行的抵偿高程归化面，然后仍然按(7)抵偿坐标系化算国家统一坐标系统一的投影带进行长度改化和方向改化，使

3、测区中X=Xf，c一Xo)1(7)心的长度变形接近于零，边缘地区的长度变形不大Y：】，c+g(ro—ro、)J‘于2．5cm／km，以达到满足城市工程建设的需要。(8)抵偿后变形值计算与可行性分析1抵偿坐标的计算步骤一+爰㈩(1)计算测区中心附近起算的横坐标2在EXCELVBA中的实现ym=1．8563×AAcosq~(1)式中：Ym——测区中心附近的横向坐标(km)；2．1代码编写△A——测区中心附近与中央子午线的经差，在Excel2003界面下，点击菜单栏的“工具(分)；(T)”一“宏(M)”一“VisualBasic编辑器(V)”，或直——测区中心附近

4、的纬度。接按组合键Alt+F11”，打开MicrosoftVisualBaic界(2)变形计算面，并在该界面下点击菜单栏的“插入(I)”一“模块(M)”，在模块代码窗口中，编写用于角度单位转换、一一—++爰‘2)抵偿坐标计算和国家与抵偿坐标之间相互转换的自(3)测区中心相对于抵偿高程归化面的高程定义函数jdzh()、dezb()和gjydczh()，其代码如下：y2函数jdzh()根据输入的角度单位，转换为需要H。薏×1000(m)(3)输出的角度单位，各参数含义为：(4)抵偿高程归化面相对于参考椭球面的高程jd一需要转换的角度值H。=(日+h)一日(4)，

5、srdw一可选项，省略值为0，表示输入的角度单式中：(+h)——测区平均高程面。位：0一度．分秒，1一度，其它值一弧度(5)缩放系数§cdw一可选项，省略值为0，表示输出的角度单Ho位：0一弧度，1一度，其它值一度·分·秒·g：(5)PublicFunctionjdzh(jd样，Optionalsrdw％=0，(6)国家统一坐标系化算抵偿坐标系Optionalsedw％=0)25第5期矿山测量2014午l0月Constpi#=3．14159265358979a(2)=6378137：f(2)=298．257222101国家Dimd％，，m#，nl％2000f

6、h=Sgn(jd)：jd=Abs(jd)dcxx(0，0)=”国家坐标系统：”SelectCasesrdwdcxx(1，0)=”投影方式：”Case0输入单位为度．分秒dcxx(2，0)=”测区中心附近纬度B：”d=Int(jd)dexx(3，0)=”测区中心附近经度L：”f=Int((jd—d)$100)dcxx(4，0)=”测区平均高程H(m)：”m=((jd—d)100一f)100dcxx(5，0)=”测区地球平均曲率半径Rm：”jdzh=(d+f／60+m／3600)dcxx(6，0)=”测区边缘起算的横坐标YmCase1输入单位为度(km)：”jd

7、zh=jddcxx(7，0)=”长度变形值6(cm／km)：”CaseElse输人单位为弧度dcxx(8，0)=”测区中心相对于抵偿高程归化jdzh=jd$180／pi面的高程Hc(m)：”EndSelectdcxx(9，0)=”抵偿高程归化面相对于参考椭Ifscdw=0Then输出单位为弧度球面的高程H0(m)：”jdzh=jdzhpi／180fhdcxx(10，0)=”抵偿后的长度变形值8(cm／ElseIfscdw=1Then输出单位为度km)：”jdzh=jdzh}fhdexx(11，0)=”缩放系数q：”Else输出单位为度。分秒”dh=Int((

8、jdzh(jd，0，1)一1．5)／3)+1带d=I