空间大地坐标系与平面直角坐标系转换公式

1、空间大地坐标系与平面直角坐 标系转换公式231坐标系的分类正如前面所提及的,所谓坐标系指的是描述空间位置的表达形式,即米用什么方法来表示空 间位置。人们为了描述空间位置，采用了多种方 法，从而也产生了不同的坐标系，如直角坐标系、 极坐标系等。在测量中常用的坐标系有以下几种：一、空间直角坐标系空间直角坐标系的坐标系原点位于参考椭球的中心，Z轴指向参考椭球的北极，X轴指向 起始子午面与赤道的交点，Y轴位于赤道面上且 按右手系与X轴呈90夹角。某点在空间中的 坐标可用该点在此坐标系的各个坐标轴上的投 影来表示。空间直角坐标系可用图 2-3来表示：图2-3空间直角坐标系、空间大地坐标系空间大地坐标系是

2、采用大地经、纬度和大地 高来描述空间位置的。纬度是空间的点与参考椭 球面的法线与赤道面的夹角；经度是空间中的点 与参考椭球的自转轴所在的面与参考椭球的起 始子午面的夹角；大地高是空间点沿参考椭球的 法线方向到参考椭球面的距离。空间大地坐标系 可用图2-4来表示：三、平面直角坐标系平面直角坐标系是利用投影变换，将空间坐 标空间直角坐标或空间大地坐标通过某种数学 变换映射到平面上，这种变换又称为投影变换。 投影变换的方法有很多，如横轴墨卡托投影、UTM 投影、兰勃特投影等。在我国采用的是高斯-克吕格 投影也称为高斯投影。UTM投影和高斯投影都是 横轴墨卡托投影的特例，只是投影的个别参数不 同而已。

3、高斯投影是一种横轴、椭圆柱面、等角投影。 从几何意义上讲，是一种横轴椭圆柱正切投影。 如图左侧所示，设想有一个椭圆柱面横套在椭球 外面，并与某一子午线相切（此子午线称为中央 子午线或轴子午线），椭球轴的中心轴CC通过 椭球中心而与地轴垂直。高斯投影满足以下两个条件：1、它是正形投影；2、中央子午线投影后应为x轴，且长度保持 不变。将中央子午线东西各一定经差（一般为6度 或3度）范围内的地区投影到椭圆柱面上， 再将 此柱面沿某一棱线展开，便构成了高斯平面直角 坐标系，如下图2 5右侧所示。图2-5高斯投影x方向指北，y方向指东。可见，高斯投影存在长度变形，为使其在测 图和用图时影响很小，应相隔一

4、定的地区，另立 中央子午线，采取分带投影的办法。我国国家测 量规定采用六度带和三度带两种分带方法。六度 带和三度带与中央子午线存在如下关系：L中 = 6N 3 ； L中 = 3n其中，N n分别为6度带和3度带的带号。 另外，为了避免y出现负号，规定y值认为 地加上500000m;又为了区别不同投影带，前面 还要冠以带号，如第20号六度带中， y=-200.25m，则成果表中写为y假定= 20499799.75m。x值在北半球总显正值，就无需 改变其观测值了。1、空间直角坐标系与空间大地坐标系间的转换图2 6表示了空间直角坐标系与空间大 地坐标系之间的关系。ZY在相同的基准下空间大地坐标系向空

5、间直 角坐标系的转换公式为：X (N H ) cos Bcos LY (N H )cos B sin L(2-1)Z N(1 e2) H sinB式中，N W，a为椭球的长半轴，N为椭球的卯酉圈曲率半径a=6378.137kmW . 1e2 sin2 Ba2b2e为椭球的第一偏心率,b为椭球的短半轴b=6356.7523141km在相同的基准下空间直角坐标系向空间大地坐标系的转换公式为2r丄 丄 ae sin BB arctg tg 1-Z W(2-2)YL arctg :RcosHNcosB式中R X2 Y2 Z22、空间坐标系与平面直角坐标系间的转换空间坐标系与平面直角坐标系间的转换采 用

6、的是投影变换的方法。在我国一般采用的是高 斯投影。因为高斯投影和UTM投影都是横轴墨卡 托的特例，因此，高斯投影和UTM投影都可以套 用横轴墨卡托投影的投影公式。横轴墨卡托投影的投影的正反算公式可参 见有关资料，它们的区别在于轴子午线投影到平 面上后，其长度的系数，对于高斯投影，系数为 1,对于UTM投影，其系数为0.9996。3、变动高程归化面的影响用户在建立地方独立坐标系时，有时变动高程归化面，这将产生一个新椭球，这就必须计算 新常数，新椭球常数按下列方法和步骤进行：1)新椭球是在国家坐标系的参考椭球上扩大形成的，它的扁率应与国家坐标系参考椭球 的扁率相等，即a a。2)计算该坐标系中央地

7、区的新椭球平均曲 率半径和新椭球长半轴。新椭球平均曲率半径为： 2 1 2RRmHmMNHmT3)；HmHm(2.10)V W W1 e SIn Bm7式中Hm该地区平均大地咼;Bm该地区的平均纬度。新椭球的长半轴按下式计算:aR1e2 sin2 Bm、1e2(2.11)将新的椭球参数代入，就可以进行投影的正 反计算了。二、坐标系统的转换方法不同坐标系统的转换本质上是不同基准间的转换，不同基准间的转换方法有很多，其中最为常用的有布尔沙模型，又称为七参数转换法 七参数转换法是：设两空间直角坐标系间有七个转换参数：3个平移参数X y z、3个旋转参数x y z和1 个尺度参数k。比如，由空间直角坐

8、标系 A转换 到空间直角坐标系B可采用下面的公式：XxXYy(1 k)YZBzZ234 GPS测量中常用的坐标系统 一、世界大地坐标系WGS-84WGS-84坐标系是目前GPS所采用的坐标系 统，GPS所发布的星历参数和历书参数等都是基 于此坐标系统的。WGS-84坐标系统的全称是 World Geodical System-84 （世界大地坐标系-84）,它是一个地 心地固坐标系统。 WGS-84坐标系统由美国国防 部制图局建立，于1987年取代了当时GPS所采 用的坐标系统WGS-72坐标系统而成为现在GPS 所使用的坐标系统。WGS-84坐标系的坐标原点位于地球的质 心，Z轴指向BIH

9、1984.0定义的协议地球极方 向，X轴指向BIH 1984.0的启始子午面和赤道 的交点，Y轴与X轴和Z轴构成右手系。WGS-8 4系所采用椭球参数为见表2.1。二、1954年北京坐标系1954年北京坐标系是我国目前广泛采用的 大地测量坐标系。该坐标系源自于原苏联采用过 的1942年普尔科夫坐标系。该坐标系米用的参 考椭球是克拉索夫斯基椭球。该椭球的参数见表 2.1。遗憾的是该椭球并未依据当时我国的天文 观测资料进行重新定位，而是由前苏联西伯利亚 地区的一等锁经我国的东北地区传算过来的，该坐标系的高程异常是以前苏联1955年大地水准 面重新平差的结果为起算值，按我国天文水准路 线推算出来的，

10、而高程又是以1956年青岛验潮 站的黄海平均海水面为基准。由于当时条件的限制1954年北京坐标系存 在着很多缺点主要表现在以下几个方面：1.克拉索夫斯基椭球参数同现代精确的椭 球参数的差异较大，并且不包含表示地球物理特 性的参数，因而给理论和实际工作带来了许多不便。2椭球定向不十分明确，椭球的短半轴既不 指向国际通用的 CIO极，也不指向目前我国使 用的JYD极。参考椭球面与我国大地水准面呈 西高东低的系统性倾斜，东部高程异常达60余米，最大达67米。3.该坐标系统的大地点坐标是经过局部分 区平差得到的。因此全国的天文大地控制点实际 上不能形成一个整体，区与区之间有较大的隙 距,如在有的接合部

11、中同一点在不同区的坐标值 相差1-2米，不同分区的尺度差异也很大，而 且坐标传递是从东北到西北和西南，后一区是以 前一区的最弱部作为坐标起算点，因而一等锁具 有明显的坐标积累误差。三、1980年西安大地坐标系1978年我国决定重新对全国天文大地网施 行整体平差，并且建立新的国家大地坐标系统。 整体平差在新大地坐标系统中进行，这个坐标系 统就是1980年西安大地坐标系统。1980年西 安大地坐标系统所采用的地球椭球参数的四个 几何和物理参数采用了 IAG 1975年的推荐值，见表2.1中的西安80椭球的短轴平行于地球的自转轴（由地球质 心指向1968.0 JYD地极原点方向），起始子午 面平行于

12、格林尼治平均天文子午面，椭球面同似 大地水准面在我国境内符合最好，高程系统以 1956年黄海平均海水面为高程起算基准。四、几种常用的坐标系统的几何和物理参数下表列出了几种常用的坐标系统的几何和 物理参数，用户需要时可以查阅：表2.1 GPS测量中常用的坐标系统的几何和 物理参数坐标系WGS-84北京54西安80a(m)637813763782456378140C2.o(或 f )-484.16685X 10-6(或1/298.25722(或1/298.3)(或1/298.257)3563)J21.08263X 10-33(rad s-1)7.292115 X/八_5107.2921155X 1

13、0-GM(m s-2)3.986005 X 10143.986005X 10142.4 GPS高程系统在测量中常用的高程系统有大地高系统、正 咼系统和正常咼系统。2.4.1大地高系统大地高系统是以参考椭球面为基准面的高 程系统，某点的大地高是该点到通过该点的参考 椭球的法线与参考椭球面的交点间的距离。 大地 高也称为椭球高。大地高一般用符号 H表示。大地高是一个纯几何量，不具有物理意义，同一个点在不同的基准下具有不同的大地高通常，GPS接收机单点定位得到的高程为WGS-84F的大地高。2.4.2正高系统正高系统是以大地水准面为基准面的高程 系统，某点的正高是该点到通过该点的铅垂线与 大地水准面的交点之间的距离。正高用符号Hg表示。243正常咼正常咼系统是以似大地水准面为基准的咼 程系统，某点的正常高是该点到通过该点的铅垂 线