基于椭球膨胀法建立独立坐标系在EXCEL VBA的实现.pdf

《基于椭球膨胀法建立独立坐标系在EXCEL VBA的实现.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第5期矿‘山测量No．52015年10月MINESURVEYING0et．2015doi：10．3969／j．issn．1001—358X．2015．05．011基于椭球膨胀法建立独立坐标系在EXCELVBA的实现术赵淑湘(甘肃林业职业技术学院，甘肃天水741020)摘要：文中介绍了基于椭球膨胀法建立独立坐标系的方法和步骤，并在EXCELVBA二次开发平台下编写了椭球长半径改正值计算和大地坐标转换的自定义函数，从而实现了自动化运算功能，不但操作简单，而且便于管理，在测绘生产中具有很好的实用性。关键词：椭球膨胀法；独立坐标系；EXCELVBA；实现中图分类号：P209文献标识码：B文章编号：

2、1001—358X【2015)05—0033—02在城市坐标系的选择与建设中，若采用国家统dL=0(4)一的3。带高斯正形投影，则规定长度变形值应不大：—(~／1兰(5)于2．5cm／km，否则要求采用独立坐标系¨。建立一esin丑)一独立坐标系的方法较多，其中椭球膨胀法是常用的N=—．=E／,=(6)一种，且该方法更有利于NGSS技术测量。基于1一esin此，作者在EXCELVBA二次开发平台下编写了椭球式(3)、式(6)、式(7)中，为子午圈曲率半径；N为膨胀法建立独立坐标系的自定义函数，具有操作简卯酉圈曲率半径。单、使用方便、易于数据管理等优点。从式(5)可以看出，椭球膨胀后不会引起

3、经度的1计算公式变化。椭球膨胀法的基本思路是膨胀后椭球与原椭球2在EXCELVBA中的实现的中心保持不变，方向保持不变，椭球扁率保持不变2．1代码设计‘(即d。=0)，仅对椭球长半径进行了改正，其改正量在Excel2003界面下，点击菜单栏的“工具的计算主要有以下3种方法乜]：(T)”一“宏(M)”一VisualBasic编辑器(V)Ah+(1)椭球长半径的改正值为投影面大地高H(dF11”，打开MicrosoftVisualBaic界面，并在该界面下=日)。点击菜单栏的“插入(I)”-+“模块(M)”，打开模块(2)投影面的大地高日等于卯酉曲率半径Ⅳ的代码窗口，编写如下代码：变化量(△Ⅳ

4、=)，椭球长半径改正值为OptionBase1设置数组默认下标下界为1d=,／1一e2sinB·H(1)函数tqpz()实现椭球长半径改正值的计算(3)以基点上参考椭球的平均曲率半径的变化PrivateFunctiontqpz(wd#，gc#，yc#，Optional值求出椭球长半径的改正值tqi％=1，Optionalj％=1)d：L垒．H(2)Dime2#，b拌，a#(4)，斜(4)，da#，af(1To2，101一eTo2)式(1)、式(2)中，曰为大地纬度；日为投影面大地高；a(1)=6378245：f(1)=298．3’北京54”e为原椭球的第一偏心率。a(2)=6378140：

5、f(2)=298．257西安80”由国家控制点的大地坐标系转换到独立坐标系a(3)=6378137：f(3)=298．257222101国中的计算公式为∞家2000”d=de2sinB—B(3)a(4)=6378137If(4)=298．257223563WGS84”}基金项目：甘肃林业职业技术学院科技计划资助fgsly一201504)af(1，1)=”da”：af(1，2)=”al”33第5期矿山测量Iftqi<>1Andtqi<>2Andtqi<>3And0．．5)tqi<>4Thentqi=4dB=e2Sin(wd)Cos(wd)／((nl+n)wd：jdzh(wd)Sqr(1一e2

6、Sin(wd)2))dab=a(tqi)一l／f(tqi)：ita(tqi)wj(1)=Round(jdzh(dB+wd，2，2)，9)e2=(a(tqi)2一b2)／a(tqi)2wj(2)=jdIfj：1Then椭球长半径的改正量为投影面大zbzh=wj地高HEndFunctionda=gc+yc2．2界面运行Elselfj=2Then投影面的大地高H等于卯自定义函数的使用方法同EXCEL内置函数。酉曲率半径N的变化量对于函数tqpz()的使用：首先在任意空白单元格连da=Sqr(1一e2$Sin(wd)2)(gc+续选定2行2列的区域，比如Dl：E2；然后在编辑栏yc)输入公式=tq

7、pz(A2，B2，c2，2，3)，同时按住[Ctr1]Else似基点上参考椭球的平均曲率半径的变+[Shift]+[Enter]键，则返回椭球长半径改正值和化值求得椭球长半径的改正值膨胀后椭球长半径的值；对于zbzh()的使用：在任意da=(1一e2Sin(wd)2)／Sqr(1一空白单元格连续选定1行2列的区域，比如D6：E6，e2)：I=(gc+yc)然后在编辑栏输入公式=zbzh(B6，C6，$E$2，2)，End