《大地坐标及直角空间坐标转换计算公式》由会员分享，可在线阅读，更多相关《大地坐标及直角空间坐标转换计算公式（12页珍藏版）》请在金锄头文库上搜索。
1、大地坐标与直角空间坐标转换计算公式一、参心大地坐标与参心空间直角坐标转换1 名词解释：A:参心空间直角坐标系：a) 以参心 0 为坐标原点；b) Z 轴与参考椭球的短轴(旋转轴)相重合；c) X 轴与起始子午面和赤道的交线重合；d) Y 轴在赤道面上与 X 轴垂直.构成右手直角坐标系 0-XYZ；e) 地面点P的点位用(X.Y.Z)表示；B:参心大地坐标系：a) 以参考椭球的中心为坐标原点.椭球的短轴与参考椭球旋转轴重合；b) 大地纬度B:以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B；c) 大地经度L:以过地面点的椭球子午面与起始子午面之间的夹角为大地经度L;d) 大地高H:地面点沿椭球
2、法线至椭球面的距离为大地高H;e) 地面点的点位用(B.L.H)表示。2 参心大地坐标转换为参心空间直角坐标:X = (N + H) * cos B * cos LY = (N + H) * cos B *sin L Z = N *(1 - e 2) + H *sin B _公式中.N为椭球面卯酉圈的曲率半径.e为椭球的第一偏心率.a、b椭球的长短半径.f椭球扁率.W为第一辅助系数W=或XX80 椭球参数:长半轴 a=63781405（m）短半轴 b=6356755.2882m扁率 a=1/298.2573 参心空间直角坐标转换参心大地坐标B = arctan(Z *( N + H)v(X2
3、+ Y2)* In*(1 -e2)+ H：X 2 + Y 2cos B-N二 高斯投影及高斯直角坐标系1、高斯投影概述高斯-克吕格投影的条件：1. 是正形投影；2. 中央子午线不变形高斯投影的性质：1. 投影后角度不变；2. 长度比与点位有关.与方向无关； 3. 离中央子午线越远变形越大为控制投影后的长度变形.采用分带投影的方法。常用 3 度带或 6 度带分带.城市或工程 控制网坐标可采用不按3 度带中央子午线的任意带。2、高斯投影正算公式：NNx = X +sin B cos Bl2 +sin B cos3 B(5 -12 + 9 2 + 4q 4) 14224N+sin Bcos5 B(6
4、1 -58t2 +14)16720Ny = N cos Bl +cos3 B(1 -12 +q 2)136N+cos5 B(5 - 18t 2 +14 + 14n 2 - 58t 2n 2)151203、高斯投影反算公式：l =cos B( 1(、2y1 一一 (1 + 2t 2 +n 2)yIN f丿6 f fNf丿1+T20+ 28t 2 + 24t 4 + 6q 2 +212 彳If f f f f I N ItB = B - f f2M+3t 2 +n 2 -刘 f 2 怜 I2+占 C+90t 2+45t f INf1坐标转换简介坐标系统之间的坐标转换既包括不同的参心坐标之间的转换.
5、或者不同的地心坐标系之 间的转换.也包括参心坐标系与地心坐标系之间的转换以及相同坐标系的直角坐标与大地坐 标之间的坐标转换.还有大地坐标与高斯平面坐标之间的转换。在两个空间角直坐标系中.假 设其分别为O-XYZ和O-XYZ.如果两个坐标系的原点相同通过三次旋转就可以使两个坐 标系重合；如果两个直角坐标系的原点不在同一个位置.通过坐标轴的平移和旋转可以取得 一致；如果两个坐标系的尺度也不尽一致就需要再增加一个尺度变化参数；而对于大地坐 标和高斯投影平面坐标之间的转换则需要通过高斯投影正算和高斯投影反第通过使用中央 子午线的经度和不同的参考椭球以及不同的投影面的选择来实现坐标的转换。如何使用Arc
6、GIS实现WGS84经纬度坐标到BJ54高斯投影坐标的转换？这是很多从事 GIS工作或者测绘工作者普遍遇到的问题。本文目的在于帮助用户解决这个问题。我们通常说的WGS-84坐标是指经纬度这种坐标表示方法.54坐标通常是指经过高斯投 影的平面直角坐标这种坐标表示方法。为什么要进行坐标转换？我们先来看两组参数如表1 所示：表1 BJ54与WGS84基准参数参考椭球体长半轴短半轴扁率BJ54基准参数Krasovsky_194063782456356863.0188298.3WGS84基准参数WGS 8463781376356752.3142298.257224很显然.WGS84与BJ54是两种不同的
7、大地基准面.不同的参考椭球体.因而两种地图下.同一个 点的坐标是不同的无论是三度带六度带坐标还是经纬度坐标都是不同的。当要把GPS接收 到的点（WGS84坐标系统的）叠加到BJ54坐标系统的底图上.那就会发现这些GPS 点不能 准确的在它该在的地方.即“与实际地点发生了偏移”。这就要求把这些GPS 点从WGS84 的坐标系统转换成BJ54的坐标系统了。有关WGS84与BJ54的坐标转换问题.实质是WGS-84椭球体到BJ54椭球体的转换问题。 如果我们是需要把WGS84的经纬度坐标转换成BJ54的高斯投影坐标那就还会涉及到投影 变换问题。因此这个转换过程一般的GPS数据处理软件都是采用下述步骤
8、进行的：1）（B.L）84（X.Y.Z） 84.空间大地坐标到空间直角坐标的转换。2）（X.Y.Z） 84（X.Y.Z） 54.坐标基准的转换.即Datum转换。通常有三种转换方法：七参数、简化三参数、Molodensky。3）（X.Y.Z） 54（B.L）54.空间直角坐标到空间大地坐标的转换。4）（B.L）54（x.y）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 和XX80也是两种不同的大地基准面.不同的参考椭球体.他们之间的转换也是同理。在 ArcGIS中提供了三参数、七参数转换法。而在同一个椭球里的转换都是严密的.在同一个椭 球的不同坐标系中转换需要用到四参数转换举个例子在XX既有54坐标又有XX坐标在这 两种坐标之间转换就用到四参数.计算四参数需要两个已知点2 ArcGIS坐标转换例子2.1 应注意问题使用ArcGIS如何实现WGS84经纬度坐标到BJ54高斯投影坐标的转换呢？在ArcGIS 中这个坐标转换步骤简化了
11、用户只需要两个步骤就能够直接从最初的WGS84经纬度坐标 转换到BJ54高斯投影坐标。这就是ArcGIS的强大之处。接下来.我们做一个例子。假设我们已经知道了7参数.应该如何操作呢？在具体的操作 前.请大家一定注意以下三点：WGS84的经纬度坐标值是用度来表示而不能是度分秒表示 七参数的平移因子单位是米旋转因子单位是秒比例因子单位是百万。在ArcGIS中.7参数法的名字是Coordinate_Frame方法。有人在用ArcGIS进行不同椭球体间的坐标转换时.转换出来的结果不对.然后就写文章 说变形如何如何.很可能是由于他们没有注意上面这三个关键的问题造成的。2.2 转换步骤a、定义 7 参数的
12、地理转换（Create Custom Geographic Transformation）在 Arctool 中打开 Create Custom Geographic Transformation工具.如图 1 所示：-X野分析工具辱服劈器工具:辱空间统计工員奇數据管理工員匡爲关聚棊霉分布式地理 （Geo datsbase） 国毎Z制團综合E抵股层和表视燮匡;乘阴表E龜子类型匡醪字段E &属性域旺曙工作空间匡窃第规E赶打包E：程投覚和变换5騒栅格用霉要素:.创建空间参考创連自定文地理冬抵）愛换 仁走叟投影J.转换坐标记（去亩备拓扌卜E曜v数据库在弹出的窗口中.输入一个转换的名字.如wgs84T
13、oBJ54。在定义地理转换方法下面.在 Method中选择合适的转换方法如COORDINATE_FRAME.然后输入平移参数、旋转角度和 比例因子.如图2所示：刨礎自定更地理（坐标滾换地理理标）变换名称WG584TOBJ54磁地辔坐标罢辅出地理坐标莱c?CSi.Eeijing_1954 亘罡女地理理标变换方號COORDINATE rRAWE鬥)|?)|?1 诽(:俳诧徒尬 移移移转转往 平平平施症備 諦厠轴轴轴轴幼 X YZ X Yn.确定；| | 取消 | | 环境 | 显示帮助AAb、投影变换打开工具箱下的Projections and TransformationsFeatureProj
14、ect在弹出的窗口中输入要 转换的数据以及Output Coordinate System.然后输入第一步自定义的地理坐标系如 wgs84ToBJ54.开始投影变换.如图3所示：点击“确定” 完成坐标转换。3结束语我国现已启用新的坐标系统2000国家大地坐标系.2000国家大地坐标系与现行国家大地 坐标系转换、衔接的过渡期仍需一段较长时期在实际工作、工程中还遇到不同坐标系之间 转换本文针对在生产中从事测绘工作遇到的坐标转换问题提供解决方法和经验希望对同行 有所参考。【参考文献】1孔祥元、郭际明、X宗泉大地测量学基础.XX大学第一版.2001年9 月2 李征航、黄劲松.GPS测量与数据处理.XX
15、大学.第一版.2005年3月3 MAPGIS使用教程大地坐标空间直角坐标转换(2009-10-22 21:12:41) 标签：杂谈程序计算大地坐标与空间直角坐标转换#include #include #define PI 3.1415926535897932384626433832795 double a,b,c,B,L,N,e,X,Y,Z,W,H;int choice;double B1=0.0,B2=0.0;double delta=0.0;int main()printf(please insert long r an);scanf(%f,&a);printf(please insert short r bn );scanf(%f,&b)
