WGS-84坐标与国家或地方坐标的转换在excel中的实现教学内容.doc

《WGS-84坐标与国家或地方坐标的转换在excel中的实现教学内容.doc》由会员分享，可在线阅读，更多相关《WGS-84坐标与国家或地方坐标的转换在excel中的实现教学内容.doc（30页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、Good is good, but better carries it.精益求精，善益求善。WGS-84坐标与国家或地方坐标的转换在excel中的实现-WGS-84坐标与国家或地方坐标的转换在excel中的实现摘要随着GPS技术的发展，精度的提高，其以全天候，高精度及操作简单的特点被越来越广泛的运用。GPS平差后结果为大地坐标，而工程中我们常用的为国家坐标系或地方独立坐标系，所以需要进行坐标转换。本文简要介绍了WGS-84坐标系和西安80坐标系及北京54坐标系等常用坐标系。通过空间直角坐标系和大地坐标系间的关系公式，用EXCEL表实现了两者的相互转换。接下来又介绍了两个不同空间直角坐标系的关系

2、，转换原理及模型，使用的是七参数布尔莎方法。用EXCEL表实现了不同空间直角坐标系间的互相转换，也就实现了WGS-84同其他坐标系间的互相转换，因为还涉及到换带计算，文中有添加了间接换带计算方法。关键词：EXCEL表，WGS-84坐标系，坐标转换ABSTRACTAlongwiththedevelopmentofGPStechnology,higherandhigherprecision,itismoreandmorewidelyappliedforitssimpleoperation,accuracyandall-weather.TheGPSresultsafteradjustmentbelo

3、ngtoWGS-84.Butcoordinateusedinengineeringisusuallyournationalcoordinatesystemorlocalindependentcoordinatesystem,soweneedcoordinatetransformation.ThispaperbrieflyintroducestheWGS-84coordinatesystemandxian80coordinatesystemandBeijing54coordinatesystemandcoordinatesystem.Throughtherelationalformulabetw

4、eentSpaceright-anglecoordinateandGeodeticcoordinates,weuseEXCELtabletoprocesscoordinateconversion.Thenintroducestwodifferentspaces,therelationship,theprincipleandthemodelparameterswhicharesevenBooleanSallymethod.EXCELachieveconversionbetweendeferentSpacecartesiancoordinatesystem,alsocanachieveWGS-84

5、withothercoordinate.Theroomisinvolvedinthecalculation,soaddedindirectlychangewithcalculationmethod.Keywords:EXCELtable,WGS-84coordinatesystem,coordinatetransformationWGS-84坐标与国家或地方坐标的相互转换在Excel中的实现1973年12月美国国防部批准海陆空三军联合研制了一种新的军用卫星导航系统即GPS卫星全球定位系统。随着GPS技术的发展，精度的提高，其以全天候，高精度及操作简单的特点被越来越广泛的运用。它能方便快捷地测定

6、出点位坐标，无论是操作上还是精度上，比全站仪等其他常规测量设备有明显的优越性。随着我国各地GPS差分台站的不断建立以及美国SA政策的取消，使得单机定位的精度大大提高，有的已经达到了亚米级精度，能够满足国土资源调查、土地利用更新、遥感监测、海域使用权清查等工作的应用。由此可以看出GPS定位技术有很大的应用潜力，近年来GPS接收机的小型化、小功耗以及卫星星座的不断改善与发展，技术也越来越成熟，实践证明，在缩短工期、降低成本和设计的灵活性方面，GPS测量较常规测量更为优越。GPS平差后的解属WGS-84坐标系，而一些工程网要求使用原坐标系统，以充分利用原有的地形图、工程图、地下管道施工图等资料，所以

7、进行坐标转换是必要的。如果原采用的是我国的8O坐标系或54坐标系(整体平差值)，我们需要把WGS-84坐标转化为国家坐标。很多情况下为了减少投影变形等原因，采用的是地方独立坐标系，所以我们也要掌握WGS-84坐标向地方独立坐标转换的方法。还有很多时候我们为了计算等需要还可能反过来，把我们现有的国家坐标或地方坐标转换为WGS-84坐标。坐标转换的方法有很多，其运算非常复杂，语言也很难理解，经常需要运用软件转换而且很多软件需要一个一个进行转换，但用EXCEL可以进行批量转换，极大的节省了时间。EXCEL表是大家熟知的一种比较普及的电子表格软件，其操作简单，可视化强，功能强大，不只能做简单的加减乘除

8、运算而且能够处理复杂的数学运算，迭加，矩阵都可以在EXCEL表中实现，用它来进行坐标转换轻松自如。本文就是运用EXCEL表进行相关计算以及坐标转换的实现。常用坐标系进行坐标转换，首先要适当了解各个坐标系统。下面简要介绍了几种常用坐标系，以便更好的理解坐标系间的关系及转换。坐标系GPS卫星定位测量是用三维地心坐标系WGS-84坐标系为依据来测定和表示总的空间位置，它既可以用地心空间坐标系（X,Y,Z）表示，也可以用托球大地坐标系（B,L,H）表示。WGS-84坐标系属于协议地球地心坐标系，最初是由美国国防部根据TRANSIT导航卫星系统的多普勒观测数据所建立的，1987年1月作为GPS卫星所发布

9、的广播星历的坐标参照基准。该坐标系的原点是地球的质心，Z轴指向BIH1984.0定义的协议地球极CTP方向，X轴指向BIH1984.0零度子午面和CTP对应的赤道的交点，Y轴和Z、X轴构成右手坐标系，它的椭球如图1所示。椭球的4个基本参数是：长半轴a=6378137m地球引力常数（含大气层）GM=3986005108m3s-2正常化二阶带球谐系数2.0=-484.1668510-6地球自转角速度=729211510-11rad/s根据以上4个参数可以进一步求得：地球扁率=0.00335281066474第一偏心率平方e2=0.0066943799013第二偏心率平方e2=0.006739496

10、74227赤道正常重力e=9.7803267714m/s2极正常重力p=9.8321863685m/s2北京54坐标系北京54坐标系是一个参心大地坐标系，原点是前苏联的普尔科沃，采用了克拉索夫斯基椭球参数，其中高程异常是以前苏联1955年大地水准面差距重新平差的结果为依据，按我国的天文水准路线换算过来的。中华人民共和国建立后，在全国范围内开展了正规的、全面的大地测量和测图工作，由此迫切建立了北京54坐标系，其与前苏联1942年坐标系进行联测，通过计算建立了我国大地坐标系，因此北京54坐标系可以认为是前苏联1942年坐标系的延伸。西安80坐标系西安80坐标系也是一个参心大地坐标系，原点为我国陕西

11、省泾阳县永乐镇，采用国际大地测量和地球物理联合会1975年推荐的四个地球椭球基本参数。其椭球短轴平行于地球质心指向我国地极原点JYD1968.0方向，大地起始子午面平行于格林尼治平均天文台的子午面。西安80坐标系与北京54坐标系相比，椭球参数更加精确，椭球面也与我国大地水准面更加吻合，是在北京54坐标系的基础上建立起来的。地方独立坐标系在一些城市测量和工程测量中，若直接采用国家坐标系可能会由于远离中央子午线，造成变形严重，误差太大，还有些特殊的测量如桥梁大坝等工程使用国家坐标系也很不方便，因此，常常会建立适合本地区的独立坐标系。空间直角坐标系和大地坐标系的互相转换坐标转换中不可避免的会用到空间

12、直角坐标和大地坐标之间的转换。比如有时我们GPS平差后的结果为大地坐标系，与其他坐标转换时需要先转换为空间直角坐标。当检核时，我们还需要反过来把空间直角坐标转换为大地坐标。所以先介绍一下两者之间的关系以及转换的实现。由图2我们可以推出空间直角坐标与椭球大地坐标之间的公式，进而实现两者的互相转换。以下是公式：表1椭球参数克拉索夫斯基椭球体1975年国际椭球体WGS-84椭球体a6378245.00000000006378140.00000000006378137.0000000000b6356863.01877304736356755.28815752876356752.3142c6399698

13、.90178271106399596.65198801056399593.62581/298.31/298.2571/298.257223563e20.0066934216229660.0066943849995880.0066943799013e20.0067385254146830.0067395018194730.00673949674227每个参考椭球的形状和大小由五个基本参数来决定，长半轴a，短半轴b，椭圆的扁率，椭圆的第一偏心率e,第二偏心率e。但只需知道其中两个参数就（至少知道一个长度元素）可以算出其他元素。我们以WGS-84椭球为例，进行空间直角坐标同椭球大地坐标的相互转换。大

14、地坐标系向空间直角坐标系的转换表2大地坐标向空间直角坐标转换的准备数据表2中B,L,H是以度分秒的格式输入的，如38度48分7.38137秒写成38.480738137B,L,H列分别为把B,L,H转成度的形式。表格H3中输入公式：=INT(B3)+(INT(B3*100)-INT(B3)*100)/60+MOD(B3,INT(B3*100)/100)*1000/3600表格I3中输入公式：=INT(C3)+(INT(C3*100)-INT(C3)*100)/60+MOD(C3,INT(C3*100)/100)*1000/3600表格J3中输入公式：=INT(D3)+(INT(D3*100)-

15、INT(D3)*100)/60+MOD(D3,INT(D3*100)/100)*1000/3600表格M2为N的计算：输入公式=K3/SQRT(1-L3*POWER(SIN(RADIANS(H3),2)根据公式1计算出X,Y,Z。表格D3中输入公式=(M3+J3)*COS(RADIANS(H3)*COS(RADIANS(I3)计算出X表格E3中输入公式=(M3+J3)*COS(RADIANS(H3)*SIN(RADIANS(I3)计算出Y表格F3中输入公式=(M3*(1-L3)+J3)*SIN(RADIANS(H3)计算出Z隐藏计算过程，使界面更直观简洁.表3大地坐标向空间直角坐标转换大地坐标

16、向空间直角坐标转换完毕。空间直角坐标系向大地坐标系的转换步骤一：求出B的值因此式中k=（3）ti为前一次迭代值，第一次迭代令ti=t0步骤二：求出L步骤三：求出H可以用EXCEL表完成相关计算及转换。一般迭代三四次就能算出B的值，在这里我们可以多迭代几次。以WGS-84坐标系椭球为例。先输入已知量X,Y,Z,c,e2,e2,直接拖动鼠标即可。表4空间直角坐标转换成大地坐标已知量根据式3，计算出t0,P,k,并拖到鼠标填充表格。表5空间直角坐标转换成大地坐标的已知计算量G列输入公式H列输入公式I列输入公式表6迭加量注：J2中输入=H2，即给t1赋值为t0。J3中输入公式=IF(ABS(J2-K2

17、)0.00000001,K2)。IF(ABS(J2-K2)0,DEGREES(ATAN(C3/B3),DEGREES(3.1415926+ATAN(C3/B3)为了界面美观简洁我们可以把不需要的列隐藏，即隐藏计算过程，最后界面如下：表8空间直角坐标向大地坐标转换方法2我们可以根据已知数据进行精度评定。表9空间直角坐标和大地坐标精度评定J列输入公式=(E3B3)/B3K列输入公式=(F3-C3)/C3L列输入公式=(G3-D3)/D3就此我们完成了大地坐标和空间直角坐标的转换。坐标系和其他坐标系的转换七参数布尔莎方法进行坐标转换GPS的测量结果和我国的西安80或北京54系坐标相差很多，有时甚至1

18、00多米，由此可见我们必须进行坐标转换才行，以空间直角坐标转换为例。进行空间直角两坐标的转换，除对坐标原点实施三个平移参数外，当坐标轴间互不平行时还存在三个旋转角度参数，以及两个坐标系尺度不一样的一个尺度变换参数。我们选用布尔莎七参数公式进行推导，如果忽略旋转参数或尺度参数等，可以推出五参数，四参数，三参数转换公式。图3为布尔莎模型由布尔莎模型我们可以推出任意点在两坐标间的转换公式，如公式3所示其中为平移参数，dK为尺度变换参数，为旋转后的坐标根据公式3我们可以完成任意两个空间直角坐标系统相互转换。当已知多个公共点（两个以上）时按最小二乘法求解转换参数。（4）式中i=1,2,3，N,若设,Y=

19、,，（5）则公式3变为误差方程设观测值等全观测，则PLx=E，则法方程进而求出转换参数单位权方差协因数阵若不是等权观测，则单位权方差协因数阵为了方便计算及理解，在这里设为等权观测，然后把这些公式输入到相应的EXCEL表中，进行任意坐标转换的实现。已知三对WGS-84和北京54的对应空间直角坐标坐标（-2142730.263,4492272.368,3975216.033）（-2142730.16,4492398.27,3975268.731）；（-2124722.344,4525873.151，3946900.096）（-2124722.276242,4525999.450226,452587

20、3.151）；（-2161391.350,4508812.888）（-2161390.903938,4508938.774386,3946534.491616）对（-2138714.416,4541513.853,3921454.813）进行坐标转换。具体过程如下：输入已知坐标,如表10所示表10已知坐标在表格中输入公式求得LX表格中的公式不用一一输入，只需拖拽即可。然后根据前文中的公式在表格中输入B如表11所示。表11B矩阵注：B中需要输入Xi等坐标的地方，需要输入公式，如K4中输入公式“=B2”。表12所求参数注：根据公式6，选中P2:P8,输入公式=MMULT(MMULT(MINVERS

21、E(MMULT(TRANSPOSE(H2:N10),H2:N10),TRANSPOSE(H2:N10),G2:G10)，按住CTRL+SHIFT+ENTER后，计算出各个参数。在Q2:S4中输入，Xi,Yi,Zi为需要转换的坐标，输入时也要按公式输入。选中U2:U4,输入公式=P2:P4+MMULT(T2:T4,P5)+MMULT(Q2:S4,P6:P8)+T2:T4，按shift+ctrl+enter,计算出最终坐标。接下来根据公式7计算精度，进行精度评定。表13精度评定精度评定步骤：V列按照公式5输入，选中V2:V8,输入=MMULT(H2:N10,P2:P8)-G2:G10，按CTRL+

22、SHIFT+ENTER；W列输入已知坐标对数，X列按照公式7输入，选中X2,输入公式=MMULT(TRANSPOSE(V2:V8),V2:V8)/(3*W2-7)转换完毕。隐藏计算过程最终界面如下：表14不同空间直角坐标的转换以上是WGS-84与北京54两个直角空间坐标系的转换。如果需要其他转换只需把已知坐标改成WGS-84坐标和西安80等你需要转换的坐标，就可以完成WGS-84和西安80和北京54等空间直角坐标的转换。进行WGS-84与北京54、西安80坐标系间的转换一般只涉及X,Y坐标，此时我们可以设Z=0,z=0,此时只有五个参数需要求，可简化计算。换带转换计算以上方法还有一个漏洞，就是

23、已知坐标必须是在同一带中，所以当WGS-84坐标向国家，地方坐标转换时，如果已知坐标不在同一带，还需要考虑换带问题，把坐标都转到同一带下，一般选用有已知坐标比较多的带中。换带有两种方法，直接法和间接法，本文详细介绍一下应用高斯投影正反算公式间接进行换带计算的间接法。应用高斯投影正反算公式间接进行换带计算的间接法这种方法的实质是把椭球大地坐标作为过渡坐标。过程如下（x,y）1(B,l)(x,y)2具体通过一个实例说明一下。比如在中央子午线=123。的1带中，有一点的平面直角坐标x1=5728374.726m,y1=+210198.193m,现要求计算出在=129。的第2带的平面直角坐标。先由x1

24、,y1换算B1,L1。B1=513843.9024,L1=1260213.1362，L=123+1260213.1362=2490213.1362。L2=2490213.1362-129=1200213.1362，然后由B,L2换算x2=5728164.378m,y2=-205079.963m。就此我们完成了此次坐标的相互转换。结论及建议根据上述论述，我们可以运用EXCEL表方便快捷的进行大地空间坐标与空间直角坐标的转换以及任意空间直角坐标系的转换。EXCEL表可视化强，操作简单，是进行坐标转换的很好选择。进行不同空间直角坐标系的转换时运用的布尔莎模型即公式3有一个缺点：一般不同空间直角坐标系

25、旋转角度很小，即设sin=，cos=1，如果角度很大误差也会相应变大。我们平常运用的北京54坐标或西安80坐标一般是空间直角坐标进行高斯投影后的坐标，高程也不是单纯的Z值。所以进行完WGS-84空间直角坐标与国家空间直角坐标转换后，要进行高斯投影。方法为国家空间直角坐标大地坐标按照国家椭球参数转换成大地坐标再进行高斯投影正算，求出平面坐标。这中间还涉及误差改正等问题，这里不再详细论述。致谢在本论文的写作过程中，我的导师陈艳红老师倾注了大量的心血，从选题到开题报告，从写作提纲，到一遍又一遍地指出每稿中的具体问题，严格把关，循循善诱，在此我表示衷心感谢。同时我还要感谢在我学习期间给我极大关心和支持

26、的各位老师以及关心我的同学和朋友。写作毕业论文是一次再系统学习的过程，毕业论文的完成，同样也意味着新的学习生活的开始。参考文献1王虎文.WGS84坐标系与任意坐标系坐标转换M.2004：05-08.2张勤，王立，张永志.全球定位系统（GPS）测量原理与应用Z.长安大学地质工程与测绘工程学院讲义.3孔祥云，郭际明.控制测量学M.武汉大学出版社，1996年10月第1版.4张勤，李家权,等.GPS控制原理及应用M.北京：科学出版社,2005年7月第1版.5韩雪培，廖邦固.地方坐标系与国家坐标系转换方法探讨N.测绘通报，2004-10-20(10).6廖超明,姜卫平,覃允森.一种有效的WGS-84坐标系与地方坐标系转换方法N.测绘通报，2008-03-18(3).7刘曦光.EXCEL在GPS坐标转换中的应用A.8胡毓钜，龚剑文，黄伟，等地图投影M北京：测绘出版社，19929于崇文，等数学地质的方法及应用M北京：冶金工10汪云甲，等矿产资源评价及其应用研究M徐州：中国矿业大学出版社，199811刘承祚，等数学地质基本方法及应用M北京：地质出版社，1981-