平面四参数法坐标转换在EXCEL VBA中的实现.pdf

《平面四参数法坐标转换在EXCEL VBA中的实现.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第1期矿山测量No．12014年2月MINESURVEYINGFeb．2014doi：10．3969／j．issn．1001—358X．2014．01．27平面四参数法坐标转换在EXCELVBA中的实现赵淑湘(甘肃林业职业技术学院，甘肃天水741020)摘要：介绍了平面四参数法坐标转换的方法与步骤，并基于EXCELVBA平台编写了参数求解和坐标转换的自定义函数。文中基于ExcelVBA平台编写的自定义函数简单实用，可供同行参考。关键词：平面四参数法；坐标转换；EXECLVBA实现中图分类号：P209文献标识码：B文章编号：1001—358

2、X(2014)01—008l一04在测绘工程中，往往会遇到不同坐标系统之间10的转换问题，常见的有WGS84坐标与国家坐标之Ol△间、不同国家坐标系之间以及国家坐标与其它地方1O△l，坐标之间的转换等。严格的转换方法为经典的七参Ol1122nn～mCOS0[数法，一般使用Bursa、Mondensky等数学模型，但这●●●●●●些模型对公共点的坐标精度、高程系统公共点的图10形强度以及数量均有严格要求；而事实上，在工程中01的现有已知资料或多或少地出现某些缺陷，而使用(2)平面四参数法相对于经典七参数法而言，条件较低，1122～nnVyz

3、···如令：V=[]比较容易实现。l0一ylI微软公司推出的MicrosoftOficeExcel是一款测～01绘技术人员比较熟悉且经常使用的自动化办公软ll10—Y2件，其数值运算和数据管理功能均很强大。因此，作2B=01Y2L=者基于ExcelVba平台编写了平面四参数法坐标转2，_，换的自定义函数。较目前常用的坐标转换软件而●●●●●●‰～言，操作简单，易于掌握，便于使用，还可以很好的进10一y，行数据管理，在测绘工程中具有很好的实用性。01Y墨一△1平面四参数法坐标转换的方法△yM=二维坐标转换是通过平移、旋转和尺度变化四，nco

4、sot个参数进行的，其转换公式为：，nSln+[]则误差方程式可表示为：=一L(3、若观测值均为等权观测，即P=E，按最小二乘△原理，可推导出转换参数解算的公式为：『10△l，一一=()曰L(4)【01Y，imCOS单位权方差为：。=．4-(V)／(2N一4)(5)mSln式中：AX、AY为平移参数；m和Ot分别为尺度2在EXCEL中的实现和旋转参数。2．1代码实现求解转换参数的误差方程为：在Excel2003界面下，点击菜单栏的“工具(T)”-÷“宏(M)”一VisualBasic编辑器(V)Alt+81第1期矿山测量2014年2月F1

5、1”，打开MicrosoftVisualBaie界面，并在该界面下sc(1，2)=x2一xl点击菜单栏的“插入(I)”一“模块(M)”，打开模块sc(2，2)=y2一y1代码窗口，编写参数求解(cs())和坐标转换(zhzhsc(3，2)=0())两个自定义函数的如下代码：sc(4，2)=1#OptionBase1ElsePublicConstPI=3．14159265358979ReDimb(n，4)：ReDimbt(4，n)：ReDim1，cs()为坐标转换参数求解的自定义函数(n)Funetioncsqj(zhqAsRange，zh

6、hAsRange)As’组成B数组Varianti=0Dimb()，bt()，l()，cs，bbtn，zhqi，zhhi，sck=1(4，2)ForEachzhqiInzhqDimxl!，y11，x21，y21，i％，j％，k％i=i+1sc(1，1)=”△X(米)”：se(2，1)=”△YIf(i+2)Mod2<>0Then(米)”b(k，i+2)=zhqi．Value8(3(3，1)=”旋转角a(秒)”：se(4，1)=”尺b(k4-1，4)=zhqi．Value度m”Else统计公共点的个数b(k，i+2)=一zhqi．Valuei

7、=0b(k+1，3)=zhqi．ValueForEachzhqiInzhqEndIfi=i+1IfiMod2=0ThenNextzhqik=k2n=ii=0’当n=2时，只有一个已知公共点EndIfIfn=2ThenNextzhqii=0Fori=】TOnForEachzhqiInzhqForj=1To4i=i+1Ifj：1AndiMod2<>0ThenIfi=lThenb(i，j)=1xl：zhqi．ValueElseIfJ=1AndiMod2=0ThenElseb(i，j)=0yl=zhqi．ValueElseIf=2AndiMod2

8、<>0ThenEndIfb(i，j=0NextzhqiElseIf·=2AndiMod2=0Theni=0b(i，j)=1ForEachzhhiInzhhEndIfi=i+1NextJIfi=