大地坐标与直角空间坐标转换计算公式.docx
《大地坐标与直角空间坐标转换计算公式.docx》由会员分享，可在线阅读，更多相关《大地坐标与直角空间坐标转换计算公式.docx（15页珍藏版）》请在冰豆网上搜索。
大地坐标与直角空间坐标转换计算公式
大地坐标与直角空间坐标转换计算公式
一、参心大地坐标与参心空间直角坐标转换
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：
地面点沿椭球法线至椭球面的距离为大地高H;
e）地面点的点位用（B，L，H）表示。
2参心大地坐标转换为参心空间直角坐标:
X=（N+H）*cosB*cosL
Y=（N+H）*cosB*sinL»
Z=[N*（1_e2）+H]*sinB”
公式中，N为椭球面卯酉圈的曲率半径，e为椭球的第一偏心率，a、b椭球的长短半
径，f椭球扁率，W为第一辅助系数
、a2-b22*f-1
e或e=
af
W=（1-e%in2B
西安80椭球参数:
长半轴a=6378140土5（m）
短半轴b=6356755.2882m
扁率a=1/298.257
3参心空间直角坐标转换参心大地坐标
Y
=arctan（）
Z*（NH）
X
=arctan（）
（X2Y2）*N*（1-e2）Hi
X2y2_n
cosB
二高斯投影及高斯直角坐标系
1、高斯投影概述
高斯-克吕格投影的条件：
1.是正形投影；2.中央子午线不变形
高斯投影的性质：
1.投影后角度不变；2.长度比与点位有关，与方向无关；
3.离中央子午线越远变形越大
为控制投影后的长度变形，采用分带投影的方法。
常用3度带或6度带分带，城市或工
程控制网坐标可采用不按3度带中央子午线的任意带。
2、高斯投影正算公式：
x二XNsinBcosBl2—sinBcos3B（5-t29244）/
224
—sinBcos5B（61-58t2t4）l6
720
y=NcosBINcos3B（1-t22）l3
6
—cos5B（5-18t2t4142-58t22）l5
120
3、高斯投影反算公式:
1坐标转换简介
坐标系统之间的坐标转换既包括不同的参心坐标之间的转换，或者不同的地心坐标系之
间的转换，也包括参心坐标系与地心坐标系之间的转换以及相同坐标系的直角坐标与大地坐标之间的坐标转换，还有大地坐标与高斯平面坐标之间的转换。
在两个空间角直坐标系中，假设其分别为O--XYZ和O--XYZ，如果两个坐标系的原点相同，通过三次旋转，就可以使两个坐标系重合；如果两个直角坐标系的原点不在同一个位置，通过坐标轴的平移和旋转可以
取得一致；如果两个坐标系的尺度也不尽一致，就需要再增加一个尺度变化参数；而对于大
地坐标和高斯投影平面坐标之间的转换，则需要通过高斯投影正算和高斯投影反算，通过使
用中央子午线的经度和不同的参考椭球以及不同的投影面的选择来实现坐标的转换。
如何使用ArcGIS实现WGS8经纬度坐标到BJ54高斯投影坐标的转换？
这是很多从事GIS工作或者测绘工作者普遍遇到的问题。
本文目的在于帮助用户解决这个问题。
我们通常说的WGS-84坐标是指经纬度这种坐标表示方法，北京54坐标通常是指经过高斯投影的平面直角坐标这种坐标表示方法。
为什么要进行坐标转换？
我们先来看两组参数，如表1所示：
表1BJ54与WGS8基准参数
参考椭球体
长半轴
短半轴
扁率
BJ54基准参数
Krasovsky_1940
6378245
6356863.0188
298.3
WGS84基准参数
WGS84
6378137
6356752.3142
298.257224
很显然，WGS8与BJ54是两种不同的大地基准面，不同的参考椭球体，因而两种地图下，同一个点的坐标是不同的，无论是三度带六度带坐标还是经纬度坐标都是不同的。
当要把
GPS接收到的点（WGS8坐标系统的）叠加到BJ54坐标系统的底图上，那就会发现这些GPS
点不能准确的在它该在的地方，即“与实际地点发生了偏移”。
这就要求把这些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转换。
通常有三种
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平移），
3个旋转因子（X旋转，Y旋转，Z旋转），一个比例因子（也叫尺度变化K）。
国内参数来
源的途径不多，一般当地测绘部门会有。
通行的做法是：
在工作区内找三个以上的已知点，利用已知点的BJ54坐标和所测WGS8坐标，通过一定的数学模型，求解七参数。
若多选几个已知点，通过平差的方法可以获得较好的精度。
如果区域范围不大，最远点间的距离不大
于30Km（经验值），这可以用三参数，即只考虑3个平移因子（X平移，Y平移，Z平移），
而将旋转因子及比例因子（X旋转，Y旋转，Z旋转，尺度变化K）都视为0，所以三参数只是七参数的一种特例。
北京54和西安80也是两种不同的大地基准面，不同的参考椭球体，
他们之间的转换也是同理。
在ArcGIS中提供了三参数、七参数转换法。
而在同一个椭球里
的转换都是严密的，在同一个椭球的不同坐标系中转换需要用到四参数转换，举个例子，在
深圳既有北京54坐标又有深圳坐标，在这两种坐标之间转换就用到四参数，计算四参数需要两个已知点
2ArcGIS坐标转换例子
2.1应注意问题
使用ArcGIS如何实现WGS8经纬度坐标到BJ54高斯投影坐标的转换呢？
在ArcGIS中,这个坐标转换步骤简化了，用户只需要两个步骤就能够直接从最初的WGS84经纬度坐标转换
到BJ54高斯投影坐标。
这就是ArcGIS的强大之处。
接下来，我们做一个例子。
假设我们已经知道了7参数，应该如何操作呢？
在具体的操作前，请大家一定注意以下三点：
WGS8的经纬度坐标值是用度来表示，而不能是度分秒表示
七参数的平移因子单位是米，旋转因子单位是秒，比例因子单位是百万。
在ArcGIS中，7参数法的名字是CoordinateFrame方法。
有人在用ArcGIS进行不同椭球体间的坐标转换时，转换出来的结果不对，然后就写文章说变形如何如何，很可能是由于他们没有注意上面这三个关键的问题造成的。
2.2转换步骤
a、定义7参数的地理转换（CreateCustomGeographicTransformation）
在Arctool中打开CreateCustomGeographicTransformation工具，如图1所示：
奇分析工具
尊服务黠工貝
谕空间统计工具三辱数据管理工眞田关系类
分布式地理数据库（GeodaUbase）
备制團综合
&・團层和表视團
H令图表
II子类型
闪岭字段
F詹属性域
笫霍工作空何
舟令常规
■i；电打包
驚・投影和变换
Gu竜删格
［工抵要素
芦创建空间参老
"创建自定交地理（坐标）变换
“定文投誌
芦转鹽坐标记法
71拓扑
I直$数据库
在弹出的窗口中，输入一个转换的名字，如wgs84ToBJ54。
在定义地理转换方法下面，在
Method中选择合适的转换方法如COORDINATE_FRAM然后输入平移参数、旋转角度和比例
因子，如图2所示：
创建自定文临理｛坐标〕费换
地理坐标）喪换名称
WGS84TOBJ54
输入地理坐标系
GCS」脂「934甘
输出地理坐标系
GCS^eijing_1954呼
自疑文地理型标浹换
方法
匚OORBIiTATI_FRAlflE*
参数
名称值
x轴平移侏）EL
tSWOH）o
Z轴平移冰）0
X轴旅转砂〕0
Y轴能转裁〕0
7铀前转用卜、n
确定
取消
环境…
显示帮助>>
11
b、投影变换
打开工具箱下的ProjectionsandTransformations>Feature>Project，在弹出的窗口
中输入要转换的数据以及OutputCoordinateSystem，然后输入第一步自定义的地理坐标
系如wgs84ToBJ54,开始投影变换，如图3所示：
投影
输入数据集前要素类
C：
\arcgis\^rcTutorVlrsckingArk^Lyst\D^t&Mshp输入坐标系何选）
輸出数据集或要素类
d：
、我的\ArcGIS«£ault.gdb\14ke3_Project
输出坐标累
点击"确定"，完成坐标转换。
3结束语
我国现已启用新的坐标系统2000国家大地坐标系，2000国家大地坐标系与现行国家大
地坐标系转换、衔接的过渡期仍需一段较长时期，在实际工作、工程中还遇到不同坐标系之
间转换，本文针对在生产中从事测绘工作遇到的坐标转换问题提供解决方法和经验，希望对
同行有所参考。
【参考文献】
[1]孔祥元、郭际明、刘宗泉•《大地测量学基础》•武汉大学出版社，第一版，2001年9
月
[2]李征航、黄劲松•《GPS测量与数据处理》•武汉大学出版社，第一版，2005年3月
[3]MAPGIS使用教程
大地坐标空间直角坐标转换
（2009-10-2221:
12:
41）
转载▼
标签：
杂谈
程序计算大地坐标与空间直角坐标转换
#include
#include
#definePI3.1415926535897932384626433832795doublea,b,c,B,L,N,e,X,Y,Z,W,H;
intchoice;
doubleB1=0.0,B2=0.0;
doubledelta=0.0;
intmain（）
{
printf（"pleaseinsertlongra\n"）;
scanf（"%f",&a）;
printf（"pleaseinsertshortrb\n"）;
scanf（"%f",&b）;
e=sqrt（a*a-b*b）/a;
c=a*a/b;
printf（"1kongzhitodadi\n2daditokongzhi\0exit"）;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）;
//doublePP2=Z;
B1=atan（Z/sqrt（X*X+Y*Y））;
delta=PI/（180*60*60*1000.0）;
while（（B1-B2）>=delta）
{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"）;
}
elseprintf（"entererror,pleaseenteragain"）;printf（"1kongzhitodadi\n2daditokongzhi\0exit"）;scanf（"%d",&choice）;
}
return0;
}
高斯投影正、反算
//高斯投影正、反算
//////6度带宽54年北京坐标系
//高斯投影由经纬度（Unit:
DD）反算大地坐标（含带号，Unit:
Metres）voidGaussProjCal（doublelongitude,doublelatitude,double*X,double*Y）{
intProjNo=0;intZoneWide;////带宽doublelongitude1,latitude1,longitude0,latitude0,X0,Y0,xval,yval;doublea,f,e2,ee,NN,T,C,A,M,iPI;
iPI=0.0174532925199433;////3.1415926535898/180.0;ZoneWide=6;////6度带宽
a=6378245.0;f=1.0/298.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（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/32+45*e2*e2*e2/1024）*sin（2*latitude1）
+（15*e2*e2/256+45*e2*e2*e2/1024）*sin（4*latitude1）-（35*e2*e2*e2/3072）*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/12
0）;
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=yval;
}
//高斯投影由大地坐标（Unit:
Metres）反算经纬度（Unit:
DD）
voidGaussProjInvCal（doubleX,doubleY,double*longitude,double*latitude）
{
intProjNo;intZoneWide;////带宽
doublelongitude1,latitude1,longitude0,latitude0,X0,Y0,xval,yval;doublee1,e2,f,a,ee,NN,T,C,M,D,R,u,fai,iPI;
iPI=0.0174532925199433;////3.1415926535898/180.0;
a=6378245.0;f=1.0/298.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/32）*sin（2*u）+（21*e1*e1/16-55*e1*e1*e1*e1/32）*sin（4*u）
+（151*e1*e1*e1/96）*sin（6*u）+（1097*e1*e1*e1*e1/512）*sin（8*u）;
C=ee*cos（fai）*cos（fai）;
T=tan（fai）*tan（fai）;
NN=a/sqrt（1.0-e2*sin（fai）*sin（fai））;
R=a*（1-e2）/sqrt（（1-e2*sin（fai）*sin（fai））*（1-e2*sin（fai）*sin（fai））*（1-e2*sin（fai）*sin（fai）））;
D=xval/NN;
//计算经度（Longitude）纬度（Latitude）
longitude1=longitude0+（D-（1+2*T+C）*D*D*D/6+（5-2*C+28*T-3*C*C+8*ee+24*T*T）*D*D*D*D*D/120）/cos（fai）;
latitude1=fai-（NN*tan（fai）/R）*（D*D/2-（5+3*T+10*C-4*C*C-9*ee）*D*D*D*D/
24+（61+90*T+298*C+45*T*T-256*ee-3*C*C）*D*D*D*D*D*D/720）;
//转换为度DD
*longitude=longitude1/iPI;
*latitude=latitude1/iPI;
}
如果有需要程序的，可以直接跟我联系，呵呵
附：
高斯正反算参数
pi=0.0174532925探M）.0174532925199433//n
长半轴a=6378245.0;扁率f=1.0/298.3;//54年北京坐标系参数
长半轴a=6378140.0;扁率f=1/298.257;〃80年西安坐标系参数
长半轴a=6378137m扁率f=1：
298.257223563。
//WGS-84坐标系
升级会员 