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

1、大地坐标与直角空间坐标转换计算公式一、参心大地坐标与参心空间直角坐标转换1名词解释：A:参心空间直角坐标系：a) 以参心0为坐标原点；b) Z轴与参考椭球的短轴(旋转轴)相重合；c) X轴与起始子午面和赤道的交线重合；d) Y轴在赤道面上与 X轴垂直，构成右手直角坐标系O-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 )* cosB* cosLY =(N +H )* cosB* sin L »Z =N * (1 _e2) +H* sin B”公式中，N为椭球面卯酉圈的曲率半径，e为椭球的第一偏心率，a、b椭球的长短半径,f椭球扁率，W为第一辅助系数a2 -b22* f -1e或e =afW = . (1 -e2*sin2BN aW西安80椭球参数：长半轴 a=6378140 土 5( m)短半轴 b=6356755.2882m扁率a =1/298.2573参心空间直角坐标转换参心大地坐标

3、YL =arcta nZ*(N H)XB = arctan(X2 Y2)* N *(1 -e2) HcosB二高斯投影及高斯直角坐标系1、高斯投影概述高斯-克吕格投影的条件：1.是正形投影；2.中央子午线不变形高斯投影的性质：1.投影后角度不变；2.长度比与点位有关，与方向无关；3.离中央子午线越远变形越大为控制投影后的长度变形，采用分带投影的方法。常用3度带或6度带分带，城市或工程控制网坐标可采用不按 3度带中央子午线的任意带。2、高斯投影正算公式：x =X Nsin BcosBl2 Nsin Bcos高斯投影反算公式： B(5 t2 9 2 4 4)l4224N5246sin Bcos5

4、B(61 -58t2 t4)l6720y = N cosBl 虫遇脚-t22)l36 cos5 B(5 -18t2 t414 2 -58t2 2)l5120cosBfy2占 5 28t2 24t4 6 2 8 2t2Nf丿B =Bf75-2 - 9 2t2NfNf61 90t245t：3601坐标转换简介坐标系统之间的坐标转换既包括不同的参心坐标之间的转换，或者不同的地心坐标系之间的转换，也包括参心坐标系与地心坐标系之间的转换以及相同坐标系的直角坐标与大地坐 标之间的坐标转换，还有大地坐标与高斯平面坐标之间的转换。在两个空间角直坐标系中， 假设其分别为 O-XYZ和O-XYZ，如果两个坐标系的

5、原点相同，通过三次旋转，就可以使两 个坐标系重合；如果两个直角坐标系的原点不在同一个位置，通过坐标轴的平移和旋转可以取得一致；如果两个坐标系的尺度也不尽一致，就需要再增加一个尺度变化参数；而对于大地坐标和高斯投影平面坐标之间的转换，则需要通过高斯投影正算和高斯投影反算，通过使用中央子午线的经度和不同的参考椭球以及不同的投影面的选择来实现坐标的转换。如何使用ArcGIS实现 WGS8经纬度坐标到 BJ54高斯投影坐标的转换？这是很多从事 GIS工作或者测绘工作者普遍遇到的问题。本文目的在于帮助用户解决这个问题。我们通常说的WGS-84坐标是指经纬度这种坐标表示方法，北京54坐标通常是指经过高 斯

6、投影的平面直角坐标这种坐标表示方法。为什么要进行坐标转换？我们先来看两组参数， 如表1所示：表1 BJ54与 WGS8基准参数参考椭球体长半轴短半轴扁率BJ54基准参数Krasovsky_194063782456356863.0188298.3WGS8基 准参数WGS 8463781376356752.3142298.257224很显然，WGS8与 BJ54是两种不同的大地基准面，不同的参考椭球体，因而两种地图下， 同一个点的坐标是不同的，无论是三度带六度带坐标还是经纬度坐标都是不同的。当要把GPS接收到的点（WGS8坐标系统的）叠加到 BJ54坐标系统的底图上，那就会发现这些GPS点不能准确

7、的在它该在的地方，即“与实际地点发生了偏移”。这就要求把这些GPS点从WGS8的坐标系统转换成 BJ54的坐标系统了。有关WGS8与 BJ54的坐标转换问题，实质是 WGS-84椭球体到BJ54椭球体的转换问题。 如果我们是需要把 WGS8的经纬度坐标转换成 BJ54的高斯投影坐标，那就还会涉及到投影 变换问题。因此，这个转换过程，一般的GPS数据处理软件都是采用下述步骤进行的：1）（ B, L） 84（ X，Y，Z）84，空间大地坐标到空间直角坐标的转换。2）（ X，Y，Z）84（X，Y，Z）54,坐标基准的转换，即Datum转换。通常有三种转换方法：七参数、简化三参数、Molodensky

8、。3）（ X, Y, Z） 54 （B, L） 54,空间直角坐标到空间大地坐标的转换。4）（ B, L） 54（ X，y）54，高斯投影正算。从以上步骤不难看出，转换的关键是第二步，转换的参数。鉴于我国曾使用不同的坐标基准（BJ54、State80、Correct54 ），各地的重力值又有很大差异，所以很难确定一套适合 全国且精度较好的转换参数。在WGS-84坐标和北京54坐标之间是不存在一套转换参数可以 全国通用的，在每个地方会不一样。必须了解，在不同的椭球之间的转换是不严密的。那么，两个椭球间的坐标转换应该是怎样的呢？ 一般而言比较严密的是用七参数法，即3个平移因子（X平移，Y平移，Z平

9、移），3个旋转因子（X旋转，Y旋转，Z旋转），一个比例因子（也叫尺度变化K）。国内参数来源的途径不多，一般当地测绘部门会有。通行的做法是：在工作区内找三个以上的已知点， 利用已知点的BJ54坐标和所测 WGS8坐标，通过一定的数学模型，求解七参数。若多选几 个已知点，通过平差的方法可以获得较好的精度。如果区域范围不大，最远点间的距离不大于30Km（经验值），这可以用三参数，即只考虑3个平移因子（X平移，Y平移，Z平移），而将旋转因子及比例因子（X旋转，Y旋转，Z旋转，尺度变化 K）都视为0，所以三参数只 是七参数的一种特例。北京 54和西安80也是两种不同的大地基准面，不同的参考椭球体， 他们

10、之间的转换也是同理。在ArcGIS中提供了三参数、七参数转换法。而在同一个椭球里的转换都是严密的，在同一个椭球的不同坐标系中转换需要用到四参数转换，举个例子，在深圳既有北京54坐标又有深圳坐标，在这两种坐标之间转换就用到四参数，计算四参数需 要两个已知点2 ArcGIS坐标转换例子2.1应注意问题使用ArcGIS如何实现 WGS8经纬度坐标到BJ54高斯投影坐标的转换呢？在 ArcGIS 中, 这个坐标转换步骤简化了， 用户只需要两个步骤就能够直接从最初的WGS84经纬度坐标转换到BJ54高斯投影坐标。这就是 ArcGIS的强大之处。接下来，我们做一个例子。假设我们已经知道了 7参数，应该如何

11、操作呢？在具体的操 作前，请大家一定注意以下三点：WGS8的经纬度坐标值是用度来表示，而不能是度分秒表示七参数的平移因子单位是米，旋转因子单位是秒，比例因子单位是百万。在ArcGIS中，7参数法的名字是 Coordin ate Frame 方法。有人在用ArcGIS进行不同椭球体间的坐标转换时，转换出来的结果不对，然后就写文 章说变形如何如何，很可能是由于他们没有注意上面这三个关键的问题造成的。2.2转换步骤a、定义 7 参数的地理转换（ Create Custom Geographic Transformation）在 Arctool 中打开 Create Custom Geographic

12、 Transformation工具，如图 1 所示：分析工具 尊服务器工員 轴空间規计工具 知数据管理工眞关系类答分布式地理数据库(Geodatabase) 霉制團综合氐團层和表视團:II 表环子类型念字段属性域笫工作空间爲常规i； 打包投影和变换 蛰删格 二备要索芦创建空间参考“创建自定交地理坐标)变换*走文投影芦转换坐标记法71毎拓扑&数据库在弹出的窗口中，输入一个转换的名字，如wgs84ToBJ54。在定义地理转换方法下面，在Method中选择合适的转换方法如COORDINATE_FRAM然后输入平移参数、旋转角度和比例因子，如图2所示:b、投影变换打开工具箱下的 Project

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

14、祥元、郭际明、刘宗泉大地测量学基础武汉大学出版社，第一版，2001年9月2 李征航、黄劲松GPS测量与数据处理武汉大学出版社，第一版，2005年3月3 MAPGIS使用教程大地坐标空间直角坐标转换(2009-10-22 21:12:41)转载标签：杂谈程序计算大地坐标与空间直角坐标转换#include <stdio.h>#include <math.h>#define PI 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(&q

15、uot;please insert long r an");scanf("%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("jin

16、gduL"); 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)

17、;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);/ double PP2=Z;B1=atan(Z/sqrt(X*X+Y*Y);delta=PI/(180*60*60*1000.0);while(B1-B2)>=delta

18、) B2=B1;B1=atan(Z+N*e*e*sin(B1)/sqrt(X*X+Y*Y);B=B1;H=Z/sin(B)-N*(1-e*e);printf("L=%f,B=%f,H=%f");else printf("enter error,please enter again"); printf("1kongzhitodadin2daditokongzhi0exit"); scanf("%d",&choice);return 0;高斯投影正、反算/ 高斯投影正、反算/6 度带宽 54 年北京坐标系/ 高

19、斯投影由经纬度 (Unit:DD) 反算大地坐标 (含带号， Unit:Metres) void GaussProjCal(double longitude, double latitude, double *X, doubl e *Y)int ProjNo=0; int ZoneWide; / 带宽double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yv al;double a,f, e2,ee, NN, T,C,A, M, iPI;iPI ZoneWide = 6;/6 度带宽a=6378245.0; f=1.0/2

20、98.3; /54 年北京坐标系参数/a=6378140.0; f=1/298.257; /80 年西安坐标系参数 ProjNo = (int)(longitude / ZoneWide) ;longitude0 = ProjNo * ZoneWide + ZoneWide / 2;longitude0 = longitude0 * iPI ; latitude0=0;longitude1 = longitude * iPI ;/ 经度转换为弧度latitude1 = latitude * iPI ; / 纬度转换为弧度 e2=2*f-f*f;ee=e2*(1.0-e2); NN=a/sqrt

21、(1.0-e2*sin(latitude1)*sin(latitude1);T=tan(latitude1)*tan(latitude1); C=ee*cos(latitude1)*cos(latitude1);A=(longitude1-longitude0)*cos(latitude1); M=a*(1-e2/4-3*e2*e2/64-5*e2*e2*e2/256)*latitude1-(3*e2/8+3*e2*e2/3 2+45*e2*e2*e2/1024)*sin(2*latitude1)+(15*e2*e2/256+45*e2*e2*e2/1024)*sin(4*latitude1)

22、-(35*e2*e2*e2/30 72)*sin(6*latitude1);xval = NN*(A+(1-T+C)*A*A*A/6+(5-18*T+T*T+72*C-58*ee)*A*A*A*A*A/120); yval = M+NN*tan(latitude1)*(A*A/2+(5-T+9*C+4*C*C)*A*A*A*A/24+(61-58*T+T*T+600*C-330*ee)*A*A*A*A*A*A/720);X0 = 1000000L*(ProjNo+1)+500000L;Y0 = 0;xval = xval+X0; yval = yval+Y0;*X = xval;*Y = yv

23、al;/ 高斯投影由大地坐标 (Unit:Metres) 反算经纬度 (Unit:DD) void GaussProjInvCal(double X, double Y, double *longitude, double *l atitude)int ProjNo; int ZoneWide; / 带宽 double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval; double e1,e2,f,a, ee, NN, T,C, M, D,R,u,fai, iPI;iPI a = 6378245.0;f = 1.0/29

24、8.3;/54 年北京坐标系参数/a=6378140.0; f=1/298.257; /80 年西安坐标系参数 ZoneWide = 6;/6 度带宽ProjNo = (int)(X/1000000L) ;/ 查找带号longitude0 = (ProjNo-1) * ZoneWide + ZoneWide / 2; longitude0 = longitude0 * iPI ;/ 中央经线X0 = ProjNo*1000000L+500000L;Y0 = 0;xval = X-X0;yval = Y-Y0;/ 带内大地坐标e2 = 2*f-f*f;e1 = (1.0-sqrt(1-e2)/(1.0+sqrt(1-e2);ee = e2/(1-e2);M = yval;u = M/(a*(1-e2/4-3*e2*e2/64-5*e2*e2*e2/256);fai = u+(3*e1/2-27*e1*e1*e1/3