预览加载中,请您耐心等待几秒...
1/6
2/6
3/6
4/6
5/6
6/6

在线预览结束,喜欢就下载吧,查找使用更方便

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

GPS经纬度与84坐标系统 摘要:文介绍了GPS测量中常用的几种坐标系统,论述了用电子表格EXCEL中的公式编 辑器,进行高斯投影正反算的公式编辑与计算,坐标转换的公式编辑与计算方法。在测量的 数据处理中,简单、方便。在GPS测量中通常采用两类坐标系统,一类是在空间固定的坐 标系统,另一类是与地球体相固联的坐标系统,也称固定坐标系统。如:WGS-84世界大地 坐标系和1980年西安大地坐标系。在实际使用中需要根据坐标系统间的转换参数进行坐标 系统的变换,来求出所使用的坐标系统的坐标。这样更有利于表达地面控制点的位置和处理 GPS观测成果,因此在GPS测量中得到了广泛的应用。 1坐标系统的介绍 1.1WGS-84坐标系统 WGS—84坐标系是目前GPS所采用的坐标系统,是由美国国防部制图局建 立,于1987年取代了当时GPS所采用的坐标系统(WGS—72坐标系统)而成为 GPS目前所使用的坐标系统。 WGS—84坐标系的坐标原点位于地球的质心,Z轴指向BIHl984.0定义的协 议地球极方向,X轴指向BIHl984.0的起始子午面和赤道的交点,Y轴与X轴和 Z轴构成右手系。WGS—84系所采用椭球参数为:a=6378138m;f=1/ 298.257223563。 1.21954年北京坐标系 1954年北京坐标系是我国目前广泛采用的大地测量坐标系。该坐标系源自 于原苏联采用过的1942年普尔科夫坐标系。建国前,我国没有统一的大地坐标 系统,建国初期,在苏联专家的建议下,我国根据当时的具体情况,建立起了全 国统一的1954年北京坐标系。该坐标采用的参考椭球是克拉索夫斯基椭球,该 椭球的参数为:a=6378245m;f=1/298.3。 该椭球并未依据当时我国的天文观测资料进行重新定位。而是由前苏联西伯 利亚地区的一等锁,经我国的东北地区传算过来的,该坐标的高程异常是以前苏 联1955年大地水准面重新平差的结果为起算值,按我国天文水准路线推算出来 的,而高程又是以1956年青岛验潮站的黄海平均海水面为基准。 1.31980年西安坐标系 1980年西安坐标系采用了全面描述椭球性质的四个基本参数a、GM、J2、 ω。四个参数的数值采用的是1975年国际大地测量与地球物理联合会16届大会 的推荐值:a=6378140m;GM=3986005x108m3/s2;J2=1082.63x10-6; ω=7292115X10-11rad/s. 1980年西安坐标系的原点位于我国的中部,陕西西安市的附近。椭球的短 轴平行于由地球质心指向我国地极原点JYD1968。0的方向,起始大地子午面平 行与我国起始天文子午面。大地点的高程是1956年青岛验潮站的黄海平均海水 面为基准。 2坐标系统的转换 一般情况下,我们使用的是1954年北京坐标系或1980年西安坐标系,而 GPS测定的坐标是WGS-84坐标系坐标,需要进行坐标系转换。对于非测量专 业的工作人员来说,虽然GPS定位操作非常容易,但坐标转换则难以掌握, EXCEL是比较普及的电子表格软件,能够处理较复杂的数学运算,用它的公式 编辑功能,进行GPS坐标转换,会非常轻松自如。要进行坐标系转换,离不开 高斯投影换算,下面分别介绍用EXCEL进行换算的方法和GPS坐标转换方法。 2.1用EXCEL进行高斯投影换算 从经纬度BL换算到高斯平面直角坐标XY(高斯投影正算),或从XY换算 成BL(高斯投影反算),一般需要专用计算机软件完成,在目前流行的换算软件 中,大都需要一个点一个点地进行,不能成批量地完成,给实际工作中带来了许 多不便。但是,通过实验发现,用EXCEL可以很直观、方便地完成坐标换算工 作,只需要在EXCEL的相应单元格中输入相应的公式即可。下面以54坐标系 为例,介绍具体的计算方法。 完成经纬度BL到平面直角坐标XY的换算。在EXCEL中,选择输入公式 的起始单元格,例如:第2行第1列(A2格)为起始单元格,各单元格的格式 如下: 单元格;单元格内容;说明 A2;输入中央子午线,以度.分秒形式输入,如115度30分则输入115.30;起算 数据L0 B2; =INT(A2)+(INT(A2*100)-INT(A2)*100)/60+(A2*10000-INT(A2*100)*100)/3600; 把L0化成度 C2;以度小数形式输入纬度值,如38°14′20″则输入38.1420;起算数据B D2;以度小数形式输入经度值;起算数据L E2; =INT(C2)+(INT(C2*100)-INT(C2)*100)/60+(C2*10000-INT(C2*100)*100