大地坐标与直角空间坐标转换计算公式

1、大地坐标与直角空间坐标转换计算公式一、参心大地坐标与参心空间直角坐标转换1 名词解释：A：参心空间直角坐标系：a) 以参心 0 为坐标原点；b) Z 轴与参考椭球的短轴（旋转轴）相重合；c) X 轴与起始子午面和赤道的交线重合；d) Y 轴在赤道面上与 X 轴垂直 . 构成右手直角坐标系 0-XYZ；e) 地面点 P 的点位用（）表示；B：参心大地坐标系：a) 以参考椭球的中心为坐标原点 . 椭球的短轴与参考椭球旋转轴重合；b)大地纬度B：以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B；c)大地经度L：以过地面点的椭球子午面与起始子午面之间的夹角为大地经度L；d)大地高 H：地面点沿椭球法

2、线至椭球面的距离为大地高H；e) 地面点的点位用（）表示。2 参心大地坐标转换为参心空间直角坐标：X( NH ) * cos B * cos LY(NH ) * cos B * sin LZ N * (1e2 )H * sin B公式中 .N 为椭球面卯酉圈的曲率半径.e 为椭球的第一偏心率.a 、 b 椭球的长短半径.f椭球扁率 .W 为第一辅助系数ea2b22 * f 1a或 ef22 BW(1e * sinaNW西安 80 椭球参数：长半轴 a=63781405（ m）短半轴 b=扁 率 =1/3 参心空间直角坐标转换参心大地坐标L arctan(Y )XB arctan(Z* (N H

3、)( X 2)Y 2 ) * N * (1 e2 ) HX 2Y 2NHcosB二 高斯投影及高斯直角坐标系1、高斯投影概述高斯 - 克吕格投影的条件：1.是正形投影； 2.中央子午线不变形高斯投影的性质：1.投影后角度不变；2.长度比与点位有关. 与方向无关；3. 离中央子午线越远变形越大为控制投影后的长度变形. 采用分带投影的方法。常用 3 度带或 6 度带分带 . 城市或工程控制网坐标可采用不按3 度带中央子午线的任意带。2、高斯投影正算公式 ：xXN sin B cosBl 2N sin B cos3 B(5 t 29 24 4） l 4224N sin B cos5 B(6158t2

4、t 4 )l 6720N cos3 B(1yN cosBlt 22 )l 36Ncos5 B(5 18t 2t 4 14258t 2 2 )l 51203、高斯投影反算公式：2l1y11 (1 2t 2ff2 )ycos B fN f6N f41528t f224t 4f 6 f28f2 t 2fy120N f2B B ft fy122922yy15 3t ffft f2M fN f12N f416190t 2f45t 4fy360N f1 坐标转换简介坐标系统之间的坐标转换既包括不同的参心坐标之间的转换. 或者不同的地心坐标系之间的转换 . 也包括参心坐标系与地心坐标系之间的转换以及相同坐标

5、系的直角坐标与大地坐标之间的坐标转换. 还有大地坐标与高斯平面坐标之间的转换。在两个空间角直坐标系中.假设其分别为O-XYZ 和 O-XYZ. 如果两个坐标系的原点相同. 通过三次旋转. 就可以使两个坐标系重合；如果两个直角坐标系的原点不在同一个位置. 通过坐标轴的平移和旋转可以取得一致；如果两个坐标系的尺度也不尽一致. 就需要再增加一个尺度变化参数；而对于大地坐标和高斯投影平面坐标之间的转换. 则需要通过高斯投影正算和高斯投影反算. 通过使用中央子午线的经度和不同的参考椭球以及不同的投影面的选择来实现坐标的转换。如何使用 ArcGIS 实现 WGS84经纬度坐标到BJ54 高斯投影坐标的转换

6、这是很多从事GIS工作或者测绘工作者普遍遇到的问题。本文目的在于帮助用户解决这个问题。我们通常说的WGS-84坐标是指经纬度这种坐标表示方法. 北京 54 坐标通常是指经过高斯投影的平面直角坐标这种坐标表示方法。为什么要进行坐标转换我们先来看两组参数. 如表1所示：表 1 BJ54 与 WGS84基准参数参考椭球体长半轴短半轴扁率BJ54 基准参数Krasovsky_19406378245WGS84基准参数WGS 846378137很显然 .WGS84与 BJ54 是两种不同的大地基准面. 不同的参考椭球体. 因而两种地图下 . 同一个点的坐标是不同的 . 无论是三度带六度带坐标还是经纬度坐标

7、都是不同的。当要把GPS接收到的点（ WGS84坐标系统的）叠加到BJ54 坐标系统的底图上. 那就会发现这些GPS点不能准确的在它该在的地方. 即“与实际地点发生了偏移”。这就要求把这些GPS点从 WGS84的坐标系统转换成BJ54 的坐标系统了。有关 WGS84与 BJ54 的坐标转换问题 . 实质是 WGS-84椭球体到 BJ54 椭球体的转换问题。如果我们是需要把 WGS84的经纬度坐标转换成 BJ54 的高斯投影坐标 . 那就还会涉及到投影变换问题。因此 . 这个转换过程 . 一般的 GPS数据处理软件都是采用下述步骤进行的：1）（） 84（） 84. 空间大地坐标到空间直角坐标的转

8、换。2）（） 84（） 54. 坐标基准的转换. 即 Datum 转换。通常有三种转换方法：七参数、简化三参数、 Molodensky 。3）（） 54（） 54. 空间直角坐标到空间大地坐标的转换。4）（） 54（） 54.高斯投影正算。从以上步骤不难看出 . 转换的关键是第二步 . 转换的参数。 鉴于我国曾使用不同的坐标基准（BJ54、State80 、Correct54 ）. 各地的重力值又有很大差异 . 所以很难确定一套适合全国且精度较好的转换参数。 在 WGS-84坐标和北京 54 坐标之间是不存在一套转换参数可以全国通用的 . 在每个地方会不一样。必须了解 . 在不同的椭球之间的转

9、换是不严密的。那么 . 两个椭球间的坐标转换应该是怎样的呢一般而言比较严密的是用七参数法. 即 3 个平移因子（X平移 .Y 平移 .Z 平移）.3 个旋转因子（ X 旋转 .Y 旋转 .Z 旋转） . 一个比例因子（也叫尺度变化K）。国内参数来源的途径不多 . 一般当地测绘部门会有。通行的做法是： 在工作区内找三个以上的已知点. 利用已知点的 BJ54 坐标和所测 WGS84坐标 . 通过一定的数学模型 . 求解七参数。 若多选几个已知点 . 通过平差的方法可以获得较好的精度。 如果区域范围不大 . 最远点间的距离不大于 30Km（经验值） .这可以用三参数 . 即只考虑 3 个平移因子（

10、X 平移 .Y 平移 .Z 平移） . 而将旋转因子及比例因子（ X 旋转 .Y 旋转 .Z 旋转 . 尺度变化 K）都视为 0. 所以三参数只是七参数的一种特例。北京54 和西安 80 也是两种不同的大地基准面. 不同的参考椭球体. 他们之间的转换也是同理。在ArcGIS 中提供了三参数、七参数转换法。而在同一个椭球里的转换都是严密的. 在同一个椭球的不同坐标系中转换需要用到四参数转换. 举个例子 . 在深圳既有北京54 坐标又有深圳坐标. 在这两种坐标之间转换就用到四参数. 计算四参数需要两个已知点2 ArcGIS 坐标转换例子应注意问题使用 ArcGIS 如何实现WGS84经纬度坐标到B

11、J54 高斯投影坐标的转换呢在ArcGIS 中 .这个坐标转换步骤简化了. 用户只需要两个步骤就能够直接从最初的WGS84经纬度坐标转换到 BJ54 高斯投影坐标。这就是ArcGIS 的强大之处。接下来 . 我们做一个例子。 假设我们已经知道了7 参数 . 应该如何操作呢在具体的操作前.请大家一定注意以下三点：WGS84的经纬度坐标值是用度来表示. 而不能是度分秒表示七参数的平移因子单位是米. 旋转因子单位是秒. 比例因子单位是百万。在 ArcGIS 中 .7 参数法的名字是 Coordinate_Frame 方法。有人在用 ArcGIS 进行不同椭球体间的坐标转换时. 转换出来的结果不对.

12、然后就写文章说变形如何如何. 很可能是由于他们没有注意上面这三个关键的问题造成的。转换步骤a、定义 7 参数的地理转换（Create Custom Geographic Transformation）在 Arctool中打开 Create Custom Geographic Transformation工具 . 如图 1 所示：在弹出的窗口中 . 输入一个转换的名字 . 如 wgs84ToBJ54。在定义地理转换方法下面 . 在 Method 中选择合适的转换方法如 COORDINATE_FRAME然.后输入平移参数、旋转角度和比例因子.如图 2所示：b、投影变换打开工具箱下的Projecti

13、ons and TransformationsFeatureProject.在弹出的窗口中输入要转换的数据以及Output Coordinate System.然后输入第一步自定义的地理坐标系如 wgs84ToBJ54. 开始投影变换 . 如图 3 所示：点击“确定”. 完成坐标转换。3 结束语我国现已启用新的坐标系统2000 国家大地坐标系.2000 国家大地坐标系与现行国家大地坐标系转换、衔接的过渡期仍需一段较长时期. 在实际工作、工程中还遇到不同坐标系之间转换 . 本文针对在生产中从事测绘工作遇到的坐标转换问题提供解决方法和经验. 希望对同行有所参考。【参考文献】1孔祥元、郭际明、刘宗泉

14、. 大地测量学基础. 武汉大学出版社. 第一版 .2001 年 9 月2 李征航、黄劲松 . GPS测量与数据处理 . 武汉大学出版社 . 第一版 .2005 年 3 月3 MAPGIS 使用教程大地坐标空间直角坐标转换(2009-10-22 21:12:41)转载标签：杂谈程序计算大地坐标与空间直角坐标转换#include #include #define PIdouble a,b,c,B,L,N,e,X,Y,Z,W,H;int choice;double B1=,B2=;double delta=;int main()printf(please insert long r an);scan

15、f(%f,&a);printf(please insert short r bn );scanf(%f,&b);e=sqrt(a*a-b*b)/a;c=a*a/b;printf(1kongzhitodadin2daditokongzhi0exit);scanf(%d,&choice);while(choice!=0) if(choice=2)printf(jingduL);scanf(%f,&L);printf(weiduB);scanf(%f,&B);printf(gaoduH);scanf(%f,&H);W=sqrt(1-e*e*sin(B)*sin(B);N=a/W;X=(N+H)*cos(B)*cos(L);Y=(N+H)*cos(B)*sin(L);Z=(N*(1-e*e)+H)*sin(B);printf(X=%f,Y=%f,Z=%f,X,Y,Z);if(choice=1)printf(zuobiao:nX=);scanf(%f,&X);printf(zuobiao:nY=);scanf(%f,&Y);printf(zuobiao:nZ=);scanf(%f,&Z);L=atan(Y/X);