GPS测量常用坐标系统及相互转换

1、GPS 测量常用坐标系统及坐标转换摘要：本文GPS测量常用坐标系统，以及GPS静态、动态测量中坐标变换的参数和方法。关键词：GPS；坐标系统；坐标转换GPS （Global Positioning System）即全球定位系统，是由美国建立的一个卫星导航定 位系统。它具有全球性、全天候、连续性和实时性的精密三维导航与定位功能，现已广泛用 于大地测量、工程测量、航空摄影测量以及地形测量等各个方面。相对于常规测量来说， GPS 测量具有测量精度高、测站间无需通视、观测时间短、仪器操作简便、全天候作业、可提供 三维坐标等特点。大大地提高了测量效率和精度。但是由于坐标系统的不同，面临着大量的 坐标转换

2、问题。对GPS技术的推广使用造成了一定的障碍。本文就GPS测量常用坐标系统及 坐标转换的原理和方法，根据作者的理解介绍如下。一、GPS测量常用坐标系统及投影 一个完整的坐标系统是由坐标系和基准两方面要素所构成的。坐标系指的是描述空间位置的 表达形式，而基准指的是为描述空间位置而定义的一系列点、线、面。在大地测量中的基准 一般是指为确定点在空间中的位置，而采用的地球椭球或参考椭球的几何参数和物理参数， 及其在空间的定位、定向方式，以及在描述空间位置时所采用的单位长度的定义。大地基准 面是利用特定椭球体对特定地区地球表面的逼近，每个国家或地区均有各自的大地基准面， 因此相对同一地理位置，不同的大地

3、基准面，它们的经纬度坐标是有差异的。基准面是在椭 球体基础上建立的，但椭球体不能代表基准面，同样的椭球体能定义不同的基准面。1、坐标系统的分类1.1、空间直角坐标系空间直角坐标系的坐标系原点位于参考椭球的中心，Z轴指向参考椭球的北极，X轴指 向起始子午面与赤道的交点，Y轴位于赤道面上，且按右手系与X轴呈90夹角。某点在空 间中的坐标可用该点在此坐标系的各个坐标轴上的投影来表示。1.2、空间大地坐标系空间大地坐标系是采用大地经度（L）、大地纬度（B）和大地高（H）来描述空间位置 的。纬度是空间的点与参考椭球面的法线与赤道面的夹角，经度是空间中的点与参考椭球的 自转轴所在的面与参考椭球的起始子午面

4、的夹角，大地高是空间点沿参考椭球的法线方向到 参考椭球面的距离。1.3、平面直角坐标系 平面直角坐标系是利用投影变换，将空间坐标（空间直角坐标或空间大地坐标）通过某 种数学变换映射到平面上，这种变换又称为投影变换。投影变换的方法有很多，女口UTM 投影，在我国采用的是高斯-克吕格投影，也称为高斯投影。2、高斯-克吕格投影高斯-克吕格（Gauss-Kruger）投影简称“高斯投影”，又名等角横切椭圆柱投影”，设 想用一个椭圆柱横切于椭球面上投影带的中央子午线，按上述投影条件，将中央子午线两侧 一定经差范围内的椭球面正形投影于椭圆柱面。将椭圆柱面沿过南北极的母线剪开展平，即 为高斯投影平面。取中央

5、子午线与赤道交点的投影为原点，中央子午线的投影为纵坐标x 轴，赤道的投影为横坐标y轴，构成高斯克吕格平面直角坐标系。高斯-克吕格投影在长度 和面积上变形很小，中央经线无变形，自中央经线向投影带边缘，变形逐渐增加，变形最大 之处在投影带内赤道的两端。为了限制高斯投影中长度变形，按一定经差将地球椭球面划分 成若干投影带。通常按经差 6 度或 3 度分为六度带或三度带。六度带自 0 度子午线起每隔经 差 6 度自西向东分带，带号依次编为第 1、260 带。三度带是在六度带的基础上分成的， 它的中央子午线与六度带的中央子午线和分带子午线重合，即自 1.5 度子午线起每隔经差3 度自西向东分带，带号依次

6、编为三度带第 1、2120 带。高斯- 克吕格投影是按分带方法 各自进行投影，故各带坐标成独立系统。以中央经线投影为纵轴(X),赤道投影为横轴(y), 两轴交点即为各带的坐标原点。纵坐标以赤道为零起算，赤道以北为正，以南为负。我国位 于北半球，纵坐标均为正值。横坐标如以中央经线为零起算，中央经线以东为正，以西为负， 横坐标出现负值，使用不便，故规定将坐标纵轴西移500 公里当作起始轴，凡是带内的横坐 标值均加 500 公里。由于高斯-克吕格投影每一个投影带的坐标都是对本带坐标原点的相对 值，所以各带的坐标完全相同，为了区别某一坐标系统属于哪一带，在横轴坐标前加上带号， 如(4231898m,2

7、1655933m)，其中 21 即为带号。3、GPS 测量常用坐标系3.1、WGS-84WGS-84坐标系是目前GPS所采用的坐标系统，WGS-84坐标系的坐标原点位于地球的质 心， Z 轴指向 BIH1984.0 定义的协议地球极方向， X 轴指向 BIH1984.0 的启始子午面和赤道 的交点， Y 轴与 X 轴和 Z 轴构成右手系。WGS-84 系所采用椭球参数为：a = 6378137mf = 1/298.257223563C =484.16685 x 10-620w = 7.292115 x 10-5 rad - s-1GM = 398600.5km3 - s-23.2、1954 年

8、北京坐标系1954 年北京坐标系是我国目前广泛采用的大地测量坐标系。该坐标系采用的参考椭球 是克拉索夫斯基椭球，该椭球的参数为：a = 6378245mf =1/ 298.3由于当时条件的限制， 1954 年北京坐标系存在着很多缺点，主要表现在以下几个方面： 克拉索夫斯基椭球参数同现代精确的椭球参数的差异较大，并且不包含表示地球物理特性的 参数，因而给理论和实际工作带来了许多不便。 椭球定向不十分明确，椭球的短半轴既不指向国际通用的 CIO 极，也不指向目前我国使用的 JYD 极。参考椭球面与我国大地水准面呈西高东低的系统性倾斜，东部高程异常达60 余米， 最大达 67 米。该坐标系统的大地点

9、坐标是经过局部分区平差得到的，因此，全国的天文大地控制点实 际上不能形成一个整体，区与区之间有较大的隙距，如在有的接合部中，同一点在不同区的 坐标值相差 1-2 米，不同分区的尺度差异也很大，而且坐标传递是从东北到西北和西南，后 一区是以前一区的最弱部作为坐标起算点，因而一等锁具有明显的坐标积累误差。3.3、1980 年西安大地坐标系1978 年，我国决定重新对全国天文大地网施行整体平差，并且建立新的国家大地坐标系统 整体平差在新大地坐标系统中进行，这个坐标系统就是1980 年西安大地坐标系统。 1980 年 西安大地坐标系统所采用的地球椭球参数的四个几何和物理参数采用了 IAG 1975年的

10、推荐值，它们是a 二 6378140mGM - 3.986005 x 1014 m3 - s-2J 二 1.08263 x 10-32& = 7.292115 x 10-5 rad - s-1根据上面所给的参数，可算出1980 年西安大地坐标系所采用的参考椭球的扁率为：f = 1298.257椭球的短轴平行于地球的自转轴（由地球质心指向1968.0 JYD 地极原点方向），起始子午面 平行于格林尼治平均天文子午面，椭球面同似大地水准面在我国境内符合最好，高程系统以 1956 年黄海平均海水面为高程起算基准。二、坐标系统的转换方法GPS定位结果属于协议地球地心坐标系，即WGS-84坐标系，而我们

11、使用的测量成果是 国家参心坐标系，需要进行坐标系转换。在相同的基准下，空间直角坐标系与空间大地坐标系间的转换，在GPS实际测量中意义 不大。空间坐标系与平面直角坐标系间的转换（高斯正反算），只要正确设置参考椭球的长 半轴和扁率、中央子午线的度数及变形系数（高斯投影中央子午线变形系数为1）、北（0） 东（500000）偏移值，一般GPS测量仪器自带软件都可自行完成。这里不再赘述，现仅就 GPS 测量中常见坐标转换问题介绍如下。1、静态测量中的坐标转换常规方法是 GPS 待测点与已知地方坐标控制网基线联测（基线越短越精确），通过网 平差求解GPS待测点的地方坐标。确定转换方程的关键是根据三个以上已

12、知参考点（两类坐 标系的坐标值都精确确定）用最小二乘法求解转换参数，常用布尔莎（Bursa）7参数法进行。即三个平移参数 X、 Y、 Z,三个旋转参数，比例因子m。XYZ布尔莎模型公式可表示为：XXAXXBAAY=Y+AY+ KYBAAZZAZZBAAm其中一 0-ZAYAX 一AK 二Z0-XYAAA-YX0ZAAA2、动态测量中的坐标转换2.1、已知点有地方坐标但无WGS84坐标已知控制点只有地方坐标（北京54 或西安 80），可采用的唯一方法是联测控制点， 取 得已知点WGS84坐标，求解转换参数。这种方法至少需要二个控制点（例如点0和点A）, 设基准站于点O，点A作为方位点。先将两点联

13、测得到的WGS84经纬度（B ,L ）、（B ,L ）投影为WGS84平面坐标O O 84 A A 84（X , Y ）、（X , Y ）；再按照静态测量方法分别求出的WGS84平面坐标（X , Y ）、O O 84 A A 84 O O 84（X , Y ）和已知地方坐标（X ,Y ）、仗A A 84O O LA YA ）l下的基线边长（X , Y ）和已知地方坐标（X ,Y ）、仗A A 84O O LA YA ）l下的基线边长S84、SL和方位角a 84a：L边长：S八AX 2 + AY 284 V=2 + AY 2a方位角：a 84 二 arctan=arctanAXXAXOAXXAX

14、OAYYYAYYY84A84O84LALOL其中：L求出WGS84坐标向地方坐标转换的平移参数AX、旋转参数a和比例因子m：AY平移参数：AYX _X _O一OY Y O 84O L旋转参数：a 二a -a比例因子：求出测将测量点（例如点B）WGS84经纬度匕，Lb 24投影为WGS84平面坐标（XB，J I求出测量点相对于基准点在WGS84平面上的坐标增量AX B 量点相对于基准点在WGS84平面上的坐标增量BAYB 84AXBAYBL，最后求出测量点在当地坐标系中的坐标。AX 一X 一X 一B=AXBAYBL，最后求出测量点在当地坐标系中的坐标。AX 一X 一X 一B=B一OAY D YD

15、 Y- n 坐标增量：B 84B 84O 84AXAY/、cos a sin a=v + m一 sin a cos aXBXOAX+AXBXO+AXBYY- AYAYY- AYBBL坐标：84LLLAXBAYB 842.2、已知点既有地方坐标又有WGS84坐标已知控制点既有地方坐标又有WGS84坐标，可用2.1所述方法任取两点求转换参数，然 后用转换参数求出任一测量点在地方坐标系中的坐标。2.3 只有一个已知点或无已知点只有一个已知点，用已知点的地方坐标与WGS84平面坐标之差作为平移参数（即三参 数）,然后用平移参数求出任一测量点在地方坐标系中的坐标。当工作区范围不大时（W30Km） 精度也能够满足要求。没有已知点的时候，可采取直接投影的方式，采用WGS84平面坐标，也可用北京54或 西安 80 的椭球参数，通过高斯正算得出该点平面坐标，这种坐标其实是假定坐标。此时和 只有一个已知点的情况一样。如果需要，此方法测得的准备可提供坐标平移、旋转，转换成统一北京54或西安80 地方坐标。3、GPS 高程转换GPS 高程是以 WGS84 椭球为基准的大地高，实际应用的地面高程是以似大地水准面起算 的正常高，两者存在高程异常。用水准测量的方法联测若干GPS点的正常高，求出高程异 常E