大地坐标直角空间坐标转换计算公式
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
大地坐标与直角空间坐标转换计算公式
一、参心大地坐标与参心空间直角坐标转换
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 :地面点沿椭球法线至椭球面的距离为大地高H ; e) 地面点的点位用(B ,L ,H )表示。
2 参心大地坐标转换为参心空间直角坐标:
⎪⎭
⎪
⎬⎫
+-=+=+=B H e N Z L B H N Y L B H N X sin *])1(*[sin *cos *)(cos *cos *)(2 公式中,N 为椭球面卯酉圈的曲率半径,e 为椭球的第一偏心率,a 、b 椭球的长短半径,f 椭球扁率,W 为第一辅助系数
a b a e 2
2-=
或 f f e 1*2-= W
a
N B W e =
-=22
sin *1(
西安80椭球参数:
长半轴a=6378140±5(m )
短半轴b=6356755.2882m 扁 率α=1/298.257
3 参心空间直角坐标转换参心大地坐标
[
]
N
B
Y X H H
e N Y X H N Z B X Y
L -+=
+-++==cos ))1(**)()
(*arctan()
arctan(2
22
2
2
二 高斯投影及高斯直角坐标系
1、高斯投影概述
高斯-克吕格投影的条件:1. 是正形投影;2. 中央子午线不变形
高斯投影的性质:1. 投影后角度不变;2. 长度比与点位有关,与方向无关; 3. 离中央子午线越远变形越大
为控制投影后的长度变形,采用分带投影的方法。常用3度带或6度带分带,城市或工程控制网坐标可采用不按3度带中央子午线的任意带。
2、高斯投影正算公式:
52224253
2236
425442232)5814185(cos 120
)1(cos 6
cos )5861(cos sin 720
495(cos sin 24cos sin 2l t t t B N l t B N
Bl N y l t t B B N l t B B N
Bl B N X x ηηηηη-++-++-+=+-+++-++
=)
3、高斯投影反算公式:
()()()⎥⎥⎦
⎤⎪⎪⎭⎫ ⎝⎛+++⎪⎪⎭
⎫
⎝⎛⎢⎣⎡-++-⎪⎪⎭⎫ ⎝⎛-=⎥⎥⎥⎥⎦
⎤⎪⎪⎭⎫ ⎝⎛+++++⎢⎢⎣⎡⎪⎪⎭
⎫
⎝⎛++-⎪⎪⎭⎫ ⎝⎛=
4422
22224222422
22
4590613601 935121128624285120
1 )21(6
11cos 1f f f f
f f f f f
f
f
f f
f f f f f f
f f f
f N y t t N y t t N y y M t B B N
y t t t N y t N y B l ηηηηη
1 坐标转换简介
坐标系统之间的坐标转换既包括不同的参心坐标之间的转换,或者不同的地心坐标系之间的转换,也包括参心坐标系与地心坐标系之间的转换以及相同坐标系的直角坐标与大地坐标之间的坐标转换,还有大地坐标与高斯平面坐标之间的转换。在两个空间角直坐标系中,假设其分别为O--XYZ 和O--XYZ ,如果两个坐标系的原点相同,通过三次旋转,就可以使两个坐标系重合;如果两个直角坐标系的原点不在同一个位置,通过坐标轴的平移和旋转可以取得一致;如果两个坐标系的尺度也不尽一致,就需要再增加一个尺度变化参数;而对于大地坐标和高斯投影平面坐标之间的转换,则需要通过高斯投影正算和高斯投影反算,通过使用中央子午线的经度和不同的参考椭球以及不同的投影面的选择来实现坐标的转换。 如何使用ArcGIS 实现WGS84经纬度坐标到BJ54高斯投影坐标的转换?这是很多从事GIS 工作或者测绘工作者普遍遇到的问题。本文目的在于帮助用户解决这个问题。 我们通常说的WGS-84坐标是指经纬度这种坐标表示方法,北京54坐标通常是指经过高斯投影的平面直角坐标这种坐标表示方法。为什么要进行坐标转换?我们先来看两组参数,如表1所示:
表1 BJ54与WGS84基准参数
很显然,WGS84与BJ54是两种不同的大地基准面,不同的参考椭球体,因而两种地图下,同一个点的坐标是不同的,无论是三度带六度带坐标还是经纬度坐标都是不同的。当要把GPS 接收到的点(WGS84坐标系统的)叠加到BJ54坐标系统的底图上,那就会发现这些GPS 点不能准确的在它该在的地方,即“与实际地点发生了偏移”。这就要求把这些GPS 点从WGS84的坐标系统转换成BJ54的坐标系统了。 有关WGS84与BJ54的坐标转换问题,实质是WGS-84椭球体到BJ54椭球体的转换问题。如果我们是需要把WGS84的经纬度坐标转换成BJ54的高斯投影坐标,那就还会涉及到投影变换问题。因此,这个转换过程,一般的GPS 数据处理软件都是采用下述步骤进行的: 1)(B ,L )84——(X ,Y ,Z )84,空间大地坐标到空间直角坐标的转换。
2)(X ,Y ,Z )84——(X ,Y ,Z )54,坐标基准的转换,即Datum 转换。通常有三种