《2022年WGS-坐标与国家或地方坐标的转换在excel中的实现 .pdf》由会员分享，可在线阅读，更多相关《2022年WGS-坐标与国家或地方坐标的转换在excel中的实现 .pdf（23页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、WGS-84 坐标与国家或地方坐标的转换在excel 中的实现摘 要随着 GPS技术的发展，精度的提高，其以全天候，高精度及操作简单的特点被越来越广泛的运用。GPS平差后结果为大地坐标，而工程中我们常用的为国家坐标系或地方独立坐标系，所以需要进行坐标转换。本文简要介绍了WGS-84 坐标系和西安80 坐标系及北京54坐标系等常用坐标系。通过空间直角坐标系和大地坐标系间的关系公式，用 EXCEL 表实现了两者的相互转换。接下来又介绍了两个不同空间直角坐标系的关系，转换原理及模型，使用的是七参数布尔莎方法。用 EXCEL 表实现了不同空间直角坐标系间的互相转换，也就实现了WGS-84同其他坐标系间
2、的互相转换，因为还涉及到换带计算，文中有添加了间接换带计算方法。关键词：EXCEL 表，WGS-84坐标系，坐标转换ABSTRACT Along with the development of GPStechnology,higher and higher precision,it is more and morewidely applied for its simple operation,accuracy andall-weather.The GPS resultsafter adjustmentbelong to WGS-84.But coordinate used in enginee
3、ring is usuallyour national coordinate system or local independentcoordinate system,so weneed coordinatetransformation.This paper briefly introduces the WGS-84coordinate system andxi an 80 coordinate system and Beijing54 coordinatesystem and coordinate system.Throughthe relational formula between t
4、Space right-anglecoordinateand Geodetic coordinates,we use EXCEL table toprocess coordinateconversion.Then introduces twodifferent spaces,the relationship,theprinciple andthe model parameters which are seven Boolean Sally 名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 23 页 -method.EXCELachieve conversion between de
5、ferent Space cartesian coordinate system,also canachieve WGS-84 with other coordinate.The room isinvolved in thecalculation,so added indirectlychange with calculation method.Keywords:EXCEL table,WGS-84 coordinate system,coordinate transformation WGS-84坐标与国家或地方坐标的相互转换在Excel 中的实现1973年 12 月美国国防部批准海陆空三军
6、联合研制了一种新的军用卫星导航系统即GPS卫星全球定位系统。随着 GPS技术的发展，精度的提高，其以全天候，高精度及操作简单的特点被越来越广泛的运用。它能方便快捷地测定出点位坐标，无论是操作上还是精度上，比全站仪等其他常规测量设备有明显的优越性。随着我国各地GPS差分台站的不断建立以及美国 S A政策的取消，使得单机定位的精度大大提高，有的已经达到了亚米级精度，能够满足国土资源调查、土地利用更新、遥感监测、海域使用权清查等工作的应用。由此可以看出GPS定位技术有很大的应用潜力，近年来GPS接收机的小型化、小功耗以及卫星星座的不断改善与发展，技术也越来越成熟，实践证明，在缩短工期、降低成本和设计
7、的灵活性方面，GPS测量较常规测量更为优越。GPS平差后的解属WGS-84坐标系，而一些工程网要求使用原坐标系统，以充分利用原有的地形图、工程图、地下管道施工图等资料，所以进行坐标转换是必要的。如果原采用的是我国的8O坐标系或 54 坐标系(整体平差值)，我们需要把WGS-84坐标转化为国家坐标。很多情况下为了减少投影变形等原因，采用的是地方独立坐标系，所以我们也要掌握WGS-84坐标向地方独立坐标转换的方法。还有很多时候我们为了计算等需要还可能反过来，把我们现有的国家坐标或地方坐标转换为WGS-84坐标。坐标转换的方法有很多，其运算非常复杂，语言也很难理解，经常需要运用软件转换而且很多软件需
8、要一个一个进行转换，但用EXCEL可以进行批量转换，极大的节省了时间。名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页，共 23 页 -EXCEL 表是大家熟知的一种比较普及的电子表格软件，其操作简单，可视化强，功能强大，不只能做简单的加减乘除运算而且能够处理复杂的数学运算，迭加，矩阵都可以在EXCEL表中实现，用它来进行坐标转换轻松自如。本文就是运用EXCEL表进行相关计算以及坐标转换的实现。常用坐标系进行坐标转换，首先要适当了解各个坐标系统。下面简要介绍了几种常用坐标系，以便更好的理解坐标系间的关系及转换。坐标系GPS卫星定位测量是用三维地心坐标系WGS-84坐标系为依据来测定和表
9、示总的空间位置，它既可以用地心空间坐标系（X,Y,Z）表示，也可以用托球大地坐标系（B,L,H）表示。WGS-84 坐标系属于协议地球地心坐标系，最初是由美国国防部根据TRANSIT导航卫星系统的多普勒观测数据所建立的，1987 年 1 月作为 GPS卫星所发布的广播星历的坐标参照基准。该坐标系的原点是地球的质心，Z 轴指向 BIH1984.0 定义的协议地球极CTP方向，X轴指向 BIH1984.0 零度子午面和CTP对应的赤道的交点，Y轴和 Z、X轴构成右手坐标系，它的椭球如图1 所示。名师资料总结-精品资料欢迎下载-名师精心整理-第 3 页，共 23 页 -椭球的 4 个基本参数是：长半
10、轴 a=6378137m 地球引力常数（含大气层）GM=3986005 108m3s-2 正常化二阶带球谐系数2.0=-484.1668510-6 地球自转角速度=7292115 10-11rad/s 根据以上4 个参数可以进一步求得：地球扁率 =0.00335281066474 第一偏心率平方e2=0.0066943799013 第二偏心率平方e2=0.00673949674227赤道正常重力e=9.7803267714m/s2 名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页，共 23 页 -极正常重力p=9.8321863685m/s2 北京 54 坐标系北京 54 坐标系是一个
11、参心大地坐标系，原点是前苏联的普尔科沃，采用了克拉索夫斯基椭球参数，其中高程异常是以前苏联1955 年大地水准面差距重新平差的结果为依据，按我国的天文水准路线换算过来的。中华人民共和国建立后，在全国范围内开展了正规的、全面的大地测量和测图工作，由此迫切建立了北京54 坐标系，其与前苏联1942 年坐标系进行联测，通过计算建立了我国大地坐标系，因此北京 54 坐标系可以认为是前苏联1942 年坐标系的延伸。西安 80 坐标系西安 80 坐标系也是一个参心大地坐标系，原点为我国陕西省泾阳县永乐镇，采用国际大地测量和地球物理联合会1975 年推荐的四个地球椭球基本参数。其椭球短轴平行于地球质心指向我
12、国地极原点JYD1968.0 方向，大地起始子午面平行于格林尼治平均天文台的子午面。西安 80 坐标系与北京54 坐标系相比，椭球参数更加精确，椭球面也与我国大地水准面更加吻合，是在北京54 坐标系的基础上建立起来的。地方独立坐标系在一些城市测量和工程测量中，若直接采用国家坐标系可能会由于远离中央子午线，造成变形严重，误差太大，还有些特殊的测量如桥梁大坝等工程使用国家坐标系也很不方便，因此，常常会建立适合本地区的独立坐标系。空间直角坐标系和大地坐标系的互相转换坐标转换中不可避免的会用到空间直角坐标和大地坐标之间的转换。比如有时我们GPS平差后的结果为大地坐标系，与其他坐标转换时需要先转换为空间
13、直角坐标。当检核时，我名师资料总结-精品资料欢迎下载-名师精心整理-第 5 页，共 23 页 -们还需要反过来把空间直角坐标转换为大地坐标。所以先介绍一下两者之间的关系以及转换的实现。由图 2 我们可以推出空间直角坐标与椭球大地坐标之间的公式，进而实现两者的互相转换。以下是公式：表 1 椭球参数克拉索夫斯基椭球体1975 年国际椭球体WGS-84 椭球体a 6378245.0000000000 6378140.0000000000 6378137.0000000000 名师资料总结-精品资料欢迎下载-名师精心整理-第 6 页，共 23 页 -b 6356863.0187730473 63567
14、55.2881575287 6356752.3142 c 6399698.9017827110 6399596.6519880105 6399593.6258 1/298.3 1/298.257 1/298.257223563 e20.006693421622966 0.006694384999588 0.0066943799013 e20.006738525414683 0.006739501819473 0.00673949674227 每个参考椭球的形状和大小由五个基本参数来决定，长半轴a，短半轴 b，椭圆的扁率，椭圆的第一偏心率e,第二偏心率e。但只需知道其中两个参数就（至少知道一个长
15、度元素）可以算出其他元素。我们以 WGS-84 椭球为例，进行空间直角坐标同椭球大地坐标的相互转换。大地坐标系向空间直角坐标系的转换表 2 大地坐标向空间直角坐标转换的准备数据表 2 中 B,L,H 是以度分秒的格式输入的，如38 度 48 分 7.38137 秒写成 38.480738137 B,L,H列分别为把B,L,H 转成度的形式。表格 H3中输入公式：=INT(B3)+(INT(B3*100)-INT(B3)*100)/60+MOD(B3,INT(B3*100)/100)*1000/3600 表格 I3 中输入公式：名师资料总结-精品资料欢迎下载-名师精心整理-第 7 页，共 23
16、页 -=INT(C3)+(INT(C3*100)-INT(C3)*100)/60+MOD(C3,INT(C3*100)/100)*1000/3600 表格 J3 中输入公式：=INT(D3)+(INT(D3*100)-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
17、+J3)*COS(RADIANS(H3)*SIN(RADIANS(I3)计算出 Y 表格 F3 中输入公式=(M3*(1-L3)+J3)*SIN(RADIANS(H3)计算出 Z 隐藏计算过程，使界面更直观简洁.表 3 大地坐标向空间直角坐标转换大地坐标向空间直角坐标转换完毕。空间直角坐标系向大地坐标系的转换名师资料总结-精品资料欢迎下载-名师精心整理-第 8 页，共 23 页 -步骤一：求出B的值因此式中k=（3）ti为前一次迭代值，第一次迭代令ti=t0步骤二：求出L 步骤三：求出H 可以用 EXCEL 表完成相关计算及转换。一般迭代三四次就能算出B的值，在这里我们可以多迭代几次。以WGS
18、-84坐标系椭球为例。先输入已知量X,Y,Z,c,e2,e 2,直接拖动鼠标即可。表 4 空间直角坐标转换成大地坐标已知量名师资料总结-精品资料欢迎下载-名师精心整理-第 9 页，共 23 页 -根据式 3，计算出t0,P,k,并拖到鼠标填充表格。表 5 空间直角坐标转换成大地坐标的已知计算量G列输入公式H列输入公式I 列输入公式表 6 迭加量名师资料总结-精品资料欢迎下载-名师精心整理-第 10 页，共 23 页 -注：J2 中输入=H2，即给 t1赋值为 t0。J3 中输入公式=IF(ABS(J2-K2)0.00000001,K2)。IF(ABS(J2-K2)0,DEGREES(ATAN(
19、C3/B3),DEGREES(3.1415926+ATAN(C3/B3)为了界面美观简洁我们可以把不需要的列隐藏，即隐藏计算过程，最后界面如下：表 8 空间直角坐标向大地坐标转换方法2 我们可以根据已知数据进行精度评定。表 9 空间直角坐标和大地坐标精度评定J 列输入公式=(E3B3)/B3 名师资料总结-精品资料欢迎下载-名师精心整理-第 13 页，共 23 页 -K列输入公式=(F3-C3)/C3 L 列输入公式=(G3-D3)/D3 就此我们完成了大地坐标和空间直角坐标的转换。坐标系和其他坐标系的转换七参数布尔莎方法进行坐标转换GPS的测量结果和我国的西安80 或北京 54 系坐标相差很
20、多，有时甚至100 多米，由此可见我们必须进行坐标转换才行，以空间直角坐标转换为例。进行空间直角两坐标的转换，除对坐标原点实施三个平移参数外，当坐标轴间互不平行时还存在三个旋转角度参数，以及两个坐标系尺度不一样的一个尺度变换参数。我们选用布尔莎七参数公式进行推导，如果忽略旋转参数或尺度参数等，可以推出五参数，四参数，三参数转换公式。图 3 为布尔莎模型由布尔莎模型我们可以推出任意点在两坐标间的转换公式，如公式3 所示名师资料总结-精品资料欢迎下载-名师精心整理-第 14 页，共 23 页 -其中为平移参数，dK为尺度变换参数，为旋转后的坐标根据公式3 我们可以完成任意两个空间直角坐标系统相互转
21、换。当已知多个公共点（两个以上）时按最小二乘法求解转换参数。（4）式中 i=1,2,3,，N,若设,Y=,，（5）名师资料总结-精品资料欢迎下载-名师精心整理-第 15 页，共 23 页 -则公式 3 变为误差方程设观测值等全观测，则PLx=E，则法方程进而求出转换参数单位权方差协因数阵若不是等权观测，则单位权方差协因数阵为了方便计算及理解，在这里设为等权观测，然后把这些公式输入到相应的EXCEL表中，进行任意坐标转换的实现。已知三对WGS-84和北京 54 的对应空间直角坐标坐标（-2142730.263,4492272.368,3975216.033）（-2142730.16,449239
22、8.27,3975268.731）；（-2124722.344,4525873.151，3946900.096）（-2124722.276242,4525999.450226,4525873.151）；（-2161391.350,4508812.888）（-2161390.903938,4508938.774386,3946534.491616）对（-2138714.416,4541513.853,3921454.813）进行坐标转换。具体过程如下：名师资料总结-精品资料欢迎下载-名师精心整理-第 16 页，共 23 页 -输入已知坐标,如表 10 所示表 10 已知坐标在表格中输入公式求得L
23、X表格中的公式不用一一输入，只需拖拽即可。然后根据前文中的公式在表格中输入B如表 11 所示。名师资料总结-精品资料欢迎下载-名师精心整理-第 17 页，共 23 页 -表 11 B 矩阵注：B中需要输入Xi 等坐标的地方，需要输入公式，如K4中输入公式“=B2”。表 12 所求参数注：根据公式6，选中 P2:P8,输入公式名师资料总结-精品资料欢迎下载-名师精心整理-第 18 页，共 23 页 -=MMULT(MMULT(MINVERSE(MMULT(TRANSPOSE(H2:N10),H2:N10),TRANSPOSE(H2:N10),G2:G10)，按住 CTRL+SHIFT+ENTER
24、 后，计算出各个参数。在 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 精度评定名师资料总结-精品资料欢迎下载-名师精心整理-第 19 页，共 23 页 -精度评定步骤：V列按照公式5 输入，选中 V2:V8,输入=MMULT(H2:N10,P2:P8)-G2:G10，按 CTRL+SHIFT+ENTER；W列输入已知坐标对数，X列按照公式7
25、输入，选中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,此时只有五个参数需要求，可简化计算。换带转换计算以上方法还有一个漏洞，就是已知坐标必须是在同一带中，所
26、以当 WGS-84 坐标向国家，地方坐标转换时，如果已知坐标不在同一带，还需要考虑换带问题，把坐标都转到同一带下，一般选用有已知坐标比较多的带中。换带有两种方法，直接法和间接法，本文详细介绍一下应用高斯投影正反算公式间接进行换带计算的间接法。名师资料总结-精品资料欢迎下载-名师精心整理-第 20 页，共 23 页 -应用高斯投影正反算公式间接进行换带计算的间接法这种方法的实质是把椭球大地坐标作为过渡坐标。过程如下（x,y）1(B,l)(x,y)2具体通过一个实例说明一下。比如在中央子午线=123。的 1 带中，有一点的平面直角坐标 x1=5728374.726m,y1=+210198.193m
27、,现要求计算出在=129。的第 2 带的平面直角坐标。先由 x1,y1 换算 B1,L1。B1=51 38 43.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表可视化强，操作简单，是进行坐标转换的很好选择。进行不同空间
28、直角坐标系的转换时运用的布尔莎模型即公式3 有一个缺点：一般不同空间直角坐标系旋转角度很小，即设 sin =，cos=1，如果角度很大误差也会相应变大。我们平常运用的北京54 坐标或西安80 坐标一般是空间直角坐标进行高斯投影后的坐标，高程也不是单纯的Z 值。所以进行完WGS-84空间直角坐标与国家空间直角坐标转换后，要进名师资料总结-精品资料欢迎下载-名师精心整理-第 21 页，共 23 页 -行高斯投影。方法为国家空间直角坐标大地坐标按照国家椭球参数转换成大地坐标再进行高斯投影正算，求出平面坐标。这中间还涉及误差改正等问题，这里不再详细论述。致 谢在本论文的写作过程中，我的导师陈艳红老师倾
29、注了大量的心血，从选题到开题报告，从写作提纲，到一遍又一遍地指出每稿中的具体问题，严格把关，循循善诱，在此我表示衷心感谢。同时我还要感谢在我学习期间给我极大关心和支持的各位老师以及关心我的同学和朋友。写作毕业论文是一次再系统学习的过程，毕业论文的完成，同样也意味着新的学习生活的开始。参考文献1 王虎文.WGS 84 坐标系与任意坐标系坐标转换M.2004：05-08.2 张勤，王立，张永志.全球定位系统（GPS）测量原理与应用Z.长安大学地质工程与测绘工程学院讲义.3 孔祥云，郭际明.控制测量学 M.武汉大学出版社，1996 年 10 月第 1 版.4 张勤，李家权,等.GPS控制原理及应用M
30、.北京：科学出版社,2005 年 7 月第 1版.5 韩雪培，廖邦固.地方坐标系与国家坐标系转换方法探讨N.测绘通报，2004-10-20(10).6 廖超明,姜卫平,覃允森.一种有效的WGS-84坐标系与地方坐标系转换方法N.测绘通报，2008-03-18(3).名师资料总结-精品资料欢迎下载-名师精心整理-第 22 页，共 23 页 -7 刘曦光.EXCEL在 GPS坐标转换中的应用A.8 胡毓钜，龚剑文，黄伟，等地图投影M 北京：测绘出版社，1992 9 于崇文，等数学地质的方法及应用M 北京：冶金工10 汪云甲，等矿产资源评价及其应用研究M 徐州：中国矿业大学出版社，199811 刘承祚，等数学地质基本方法及应用M 北京：地质出版社，1981 名师资料总结-精品资料欢迎下载-名师精心整理-第 23 页，共 23 页 -