54和80坐标系浅谈及转换工具

故事结束

发送原创文章至admin@caup.net，可申请注册邀请码

您需要登录才可以下载或查看，没有帐号？注册方式

x

文章来自：测绘之家
软件来自：网上搜的

coord坐标转换软件.zip (570.82 KB, 下载次数: 73)

一、北京54坐标系

1、定义

北京54坐标系为参心大地坐标系，大地上的一点可用经度L54、纬度M54和大地高H54定位，它是以克拉索夫斯基椭球为基础，经局部平差后产生的坐标系，其坐标详细定义可参见参考文献[朱华统1990]。

2、历史

新中国成立以后，我国大地测量进入了全面发展时期，在全国范围内开展了正规的，全面的大地测量和测图工作，迫切需要建立一个参心大地坐标系。由于当时的“一边倒”政治趋向，故我国采用了前苏联的克拉索夫斯基椭球参数，并与前苏联1942年坐标系进行联测，通过计算建立了我国大地坐标系，定名为1954年北京坐标系。因此，1954年北京坐标系可以认为是前苏联1942年坐标系的延伸。它的原点不在北京而是在前苏联的普尔科沃。它是将我国一等锁与原苏联远东一等锁相连接，然后以连接处呼玛、吉拉宁、东宁基线网扩大边端点的原苏联1942年普尔科沃坐标系的坐标为起算数据，平差我国东北及东部区一等锁，这样传算过来的坐标系就定名为1954年北京坐标系。

3、特点

a.属参心大地坐标系；

b.采用克拉索夫斯基椭球的两个几何参数；

c.大地原点在原苏联的普尔科沃；

d.采用多点定位法进行椭球定位；

e.高程基准为1956年青岛验潮站求出的黄海平均海水面；

f.高程异常以原苏联1955年大地水准面重新平差结果为起算数据。按我国天文水准路线推算而得。椭球坐标参数：长半轴a=6378245m；短半轴=6356863.0188m；扁率α=1/298.3。

4、缺点

ａ.椭球参数有较大误差。克拉索夫斯基椭球差数与现代精确的椭球参数相比，长半轴约大109m。

b.参考椭球面与我国大地水准面存在着自西向东明显的系统性的倾斜，在东部地区大地水准面差距最大达+60m。这使得大比例尺地图反映地面的精度受到影响，同时也对观测量元素的归算提出了严格的要求。

c.几何大地测量和物理大地测量应用的参考面不统一。我国在处理重力数据时采用赫尔默特1900～1909年正常重力公式，与这个公式相应的赫尔默特扁球不是旋转椭球，它与克拉索夫斯基椭球是不一致的，这给实际工作带来了麻烦。

d.定向不明确。椭球短半轴的指向既不是国际普遍采用的国际协议（原点）CIO（ConventionalInternationalOrigin），也不是我国地极原点JYD1968.0；起始大地子午面也不是国际时间局BIH（BureauInternationaldeIHeure）所定义的格林尼治平均天文台子午面，从而给坐标换算带来一些不便和误差。

为此，我国在1978年在西安召开了“全国天文大地网整体平差会议”，提出了建立属于我国自己的大地坐标系，即后来的1980西安坐标系。但时至今日，北京54坐标系仍然是在我国使用最为广泛的坐标系。

二、西安80坐标系

1978年4月在西安召开全国天文大地网平差会议，确定重新定位，建立我国新的坐标系。为此有了1980年国家大地坐标系。1980年国家大地坐标系采用地球椭球基本参数为1975年国际大地测量与地球物理联合会第十六届大会推荐的数据。该坐标系的大地原点设在我国中部的陕西省泾阳县永乐镇，位于西安市西北方向约60公里，故称1980年西安坐标系，又简称西安大地原点。基准面采用青岛大港验潮站1952－1979年确定的黄海平均海水面（即1985国家高程基准）。西安80坐标系是为了进行全国天文大地网整体平差而建立的。

根据椭球定位的基本原理，在建立西安80坐标系时有以下先决条件：

（1）大地原点在我国中部，具体地点是陕西省泾阳县永乐镇；（2）西安80坐标系是参心坐标系，椭球短轴Z轴平行于地球质心指向地极原点方向，大地起始子午面平行于格林尼治平均天文台子午面；X轴在大地起始子午面内与Z轴垂直指向经度0方向；Y轴与Z、X轴成右手坐标系；

（3）椭球参数采用IUG1975年大会推荐的参数，因而可得西安80椭球两个最常用的几何参数为：长半轴a=6378140±5（m）短半轴b=6356755.2882（m）扁率α=1/298.257第一偏心率平方=0.00669438499959第二偏心率平方=0.00673950181947椭球定位时按我国范围内高程异常值平方和最小为原则求解参数。

（4）多点定位；

（5）基准面采用青岛大港验潮站1952－1979年确定的黄海平均海水面（即1985国家高程基准）。

三、北京54与西安80坐标系的转换

1、北京54与西安80坐标转换原理与方法

西安80坐标系与北京54坐标系其实是一种椭球参数的转换，这种转换在同一个椭球里的转换都是严密的，而在不同的椭球之间的转换是不严密的，因此不存在一套转换参数可以全国通用的，在每个地方会不一样，因为它们是两个不同的椭球基准。

那么，两个椭球间的坐标转换，一般而言比较严密的是用七参数布尔莎模型，即X平移，Y平移，Z平移，X旋转（WX），Y旋转（WY），Z旋转（WZ），尺度变化（DM）。要求得七参数就需要在一个地区需要3个以上的已知点。如果区域范围不大，最远点间的距离不大于30Km（经验值），这可以用三参数，即X平移，Y平移，Z平移，而将X旋转，Y旋转，Z旋转，尺度变化面DM视为0。

方法如下（MAPGIS平台中）：

第一步：向地方测绘局（或其它地方）找本区域三个公共点坐标对（即54坐标x，y，z和80坐标x，y，z）；

第二步：将三个点的坐标对全部转换以弧度为单位。（菜单：投影转换/输入单点投影转换，计算出这三个点的弧度值并记录下来）

第三步：求公共点求操作系数（菜单：投影转换/坐标系转换）。如果求出转换系数后，记录下来。

第四步：编辑坐标转换系数。（菜单：投影转换/编辑坐标转换系数。）最后进行投影变换，“当前投影”输入80坐标系参数，“目的投影”输入54坐标系参数。进行转换时系统会自动调用曾编辑过的坐标转换系数。

2、北京54坐标到西安80坐标转换小结：

a.北京54和西安80是两种不同的大地基准面，不同的参考椭球体，因而两种地图下，同一个点的坐标是不同的，无论是三度带六度带坐标还是经纬度坐标都是不同的。

b.数字化后的得到的坐标其实不是WGS84的经纬度坐标，因为54和80的转换参数至今没有公布，一般的软件中都没有54或80投影系的选项，往往会选择WGS84投影。

c.WGS84、北京54、西安80之间，没有现成的公式来完成转换。

d.对于54或80坐标，从经纬度到平面坐标（三度带或六度带）的相互转换可以借助软件完成。（54和80也有经纬度，只不过我们都用其投影的直角坐标值罢了，不能看到经纬度就以为是wgs84的）

e.54和80间的转换，必须借助现有的点和两种坐标，推算出变换参数，再对待转换坐标进行转换。（均靠软件实现）

f.在选择参考点时，注意不能选取河流、等高线、地名、高程点，公路尽量不选。这些在两幅地图上变化很大，不能用作参考。而应该选择固定物，如电站，桥梁等。

坐标转换问题（附）

工程施工过程中，常常会遇到不同坐标系统间，坐标转换的问题。目前国内常见的转换有以下几种：1，大地坐标（BLH）对平面直角坐标（XYZ）；2，北京54全国80及WGS84坐标系的相互转换；3，任意两空间坐标系的转换。其中第2类可归入第三类中。所谓坐标转换的过程就是转换参数的求解过程。常用的方法有三参数法、四参数法和七参数法。以下对上述三种情况作详细描述如下：

1，大地坐标（BLH）对平面直角坐标（XYZ）

常规的转换应先确定转换参数，即椭球参数、分带标准（3度，6度）和中央子午线的经度。椭球参数就是指平面直角坐标系采用什么样的椭球基准，对应有不同的长短轴及扁率。一般的工程中3度带应用较为广泛。对于中央子午线的确定有两种方法，一是取平面直角坐标系中Y坐标的前两位*3，即可得到对应的中央子午线的经度。如x=3250212m，y=395121123m，则中央子午线的经度=39*3=117度。另一种方法是根据大地坐标经度，如果经度是在155.5~185.5度之间，那么对应的中央子午线的经度=（155.5+185.5）/2=117度，其他情况可以据此3度类推。

另外一些工程采用自身特殊的分带标准，则对应的参数确定不在上述之列。

确定参数之后，可以用软件进行转换，以下提供坐标转换的程序下载。

2，北京54全国80及WGS84坐标系的相互转换

这三个坐标系统是当前国内较为常用的，它们均采用不同的椭球基准。

其中北京54坐标系，属三心坐标系，大地原点在苏联的普而科沃，长轴6378245m，短轴6356863，扁率1/298.3；西安80坐标系，属三心坐标系，大地原点在陕西省径阳县永乐镇，长轴6378140m，短轴6356755，扁率1/298.25722101；WGS84坐标系，长轴6378137.000m，短轴6356752.314，扁率1/298.257223563。由于采用的椭球基准不一样，并且由于投影的局限性，使的全国各地并不存在一至的转换参数。对于这种转换由于量较大，有条件的话，一般都采用GPS联测已知点，应用GPS软件自动完成坐标的转换。当然若条件不许可，且有足够的重合点，也可以进行人工解算。详细方法见第三类。

3，任意两空间坐标系的转换

由于测量坐标系和施工坐标系采用不同的标准，要进行精确转换，必须知道至少3个重合点（即为在两坐标系中坐标均为已知的点。采用布尔莎模型进行求解。布尔莎公式:

对该公式进行变换等价得到：

解算这七个参数，至少要用到三个已知点（2个坐标系统的坐标都知道），采用间接平差模型进行解算：

其中：V为残差矩阵;

X为未知七参数;

A为系数矩阵;

解之：L为闭合差

解得七参数后，利用布尔莎公式就可以进行未知点的坐标转换了，每输入一组坐标值，就能求出它在新坐标系中的坐标。但是要想GPS观测成果用于工程或者测绘，还需要将地方直角坐标转换为大地坐标，最后还要转换为平面高斯坐标。

上述方法类同于我们的间接平差，解算起来较复杂，以下提供坐标转换程序，只需输入三个已知点的坐标即可求解出坐标转换的七个参数。如果已知点的数量较多，可以进行参数间的平差运算，则精度更高。

当已知点的数量只有两个时，我们可以采用简单变换法，此法较为方便易行，适于手算，只是精度受到一定的限制。

详细解算方程如下：

式中调x,y和x\'、y\'分别为新旧（或；旧新）网重合点的坐标，a、b、、k为变换参数，显然要解算出a、b、、k，必须至少有两个重合点，列出四个方程。

即可进行通常的参数平差，解求a、x、b、c、d各参数值。将之代人（3）式，可得各拟合点的残差（改正数）代人（2）式，可得待换点的坐标。

求出解算参数之后，可在Excel中，进行其余坐标的转换。

上次笔者用此法进行过80和54坐标的转换，由于当时没有多余的点可供验证和平差，所以转换精度不得而知，但转换之后各点的相对位置不变。估计，实际的转换误差应该是10m量级的。

还有一些情况是先将大地坐标转换为直角坐标，然后进行相关转换。