7 a/ i7 g5 F( t( V3 q
在航测中可能经常会遇到不知道如何选择正确的坐标系和坐标系之间的转换,现在我们针对于航测坐标系做详细的讲解。 ! M; r0 J: p; v2 H
首先简单介绍一下航测中地理坐标系、投影坐标系以及地图投影的概念: 地理坐标系:为球面坐标。 参考平面地是椭球面,坐标单位:经纬度;投影坐标系:为平面坐标。参考平面地是水平面,坐标单位:米、千米等;地理坐标转换到投影坐标的过程可理解为投影。(投影:将不规则的地球曲面转换为平面)。针对上面三个问题,将一一介绍。
8 L) q% o0 Y) X1 h. @" K 1、地理坐标系
0 O& I: H+ p8 [ 1.1 地球的三级逼近* ?1 L0 C8 z$ y
1.1.1大地水准面2 L# k1 g. G W0 `5 }
地球的自然表面有高山也有洼地,是崎岖不平的,我们要使用数学法则来描述他,就必须找到一个相对规则的数学面。
5 U5 G: [# u# j 大地水准面是地球表面的第一级逼近。假设当海水处于完全静止的平衡状态时,从海平面延伸到所有大陆下部,而与地球重力方向处处正交的一个连续、闭合的曲面,这就是大地水准面。
. p6 k9 ]. \( {$ I/ [ 1.1.2地球椭球体
5 d( u! a. G; ~ g' k4 i 大地水准面可以近似成一个规则成椭球体,但并不是完全规则,其形状接近一个扁率极小的椭圆绕短轴旋转所形成的规则椭球体,这个椭球体称为地球椭球体。它是地球的第二级逼近。
2 U; r4 n4 T8 j @5 o& d 下面列举了一些常见椭球体的参数。我国1952年以前采用海福特椭球体,从1953年起采用克拉索夫斯基椭球体。 1978年我国决定采用新椭球体GRS(1975),并以此建立了我国新的、独立的大地坐标系,对应ArcGIS里面的Xian_1980椭球体。从1980年开始采用新椭球体GRS(1980),这个椭球体参数与ArcGIS中的CGCS2000椭球体相同。 p* e9 ~, H! K9 Z) I$ F
1.1.3大地基准面1 Y+ M8 F6 k& w, V# u& i
确定了一个规则的椭球表面以后,我们会发现还有一个问题,参考椭球体是对地球的抽象,因此其并不能去地球表面完全重合,在设置参考椭球体的时候必然会出现有的地方贴近的好(参考椭球体与地球表面位置接近),有地地方贴近的不好的问题,因此这里还需要一个大地基准面来控制参考椭球和地球的相对位置。 这是地球表面的第三级逼近。有以下两类基准面:
4 n) q$ M$ K' \( T" T 地心基准面:由卫星数据得到,使用地球的质心作为原点,使用最广泛的是 WGS 1984。 I: x) J! ]3 p. I: a2 Z% M: W
区域基准面:特定区域内与地球表面吻合,大地原点是参考椭球与大地水准面相切的点,例如Beijing54、Xian80。我们通常称谓的Beijing54、Xian80坐标系实际上指的是我国的两个大地基准面。
( W' E% c; W B& k 我们通常说的参心大地坐标系和地心大地坐标系的区别就在于此。 - S U" F& g8 M! u: n, ?& D$ a
参心大地坐标系:指经过定位与定向后,地球椭球的中心不与地球质心重合而是接近地球质心。区域性大地坐标系。是我国基本测图和常规大地测量的基础。如Beijing54、Xian80。 1 h8 r) f& Q3 W
地心大地坐标系:指经过定位与定向后,地球椭球的中心与地球质心重合。如CGCS2000、WGS84。 3 W3 V5 }8 D; I
1.2地理坐标* F/ s! A7 C9 ~9 P. k8 k L9 C
地理坐标,就是用经线(子午线)、纬线、经度、纬度表示地面点位的球面坐标。 ; X* Q$ g# d; R4 w
一般地理坐标可分为三种,天文经纬度,大地经纬度,地心经纬度。通常地图上使用的经纬度都为大地经纬度,所以这里我介绍一下大地经纬度,其他两种要想了解的话可以百度一下,其实区别不大。 # X/ G' C/ A" q M) T
大地经纬度:
6 n2 ~. t6 q- d 大地经度:参考椭球面上某点的大地子午面与本初子午面间的两面角。东正西负;
( D+ r2 f! ^+ e" t' U 大地纬度 :参考椭球面上某点的法线与赤道平面的夹角。北正南负; 7 {, Z6 a8 F$ K; i1 ^" |
大地高: 指某点沿法线方向到参考椭球面的距离。 9 ?; }: ]0 Z( d4 Z. W' ^
看到这里,地理坐标系的思路基本明确的了吧!只需要参考椭球体参数以及大地基准面就可以确定地理坐标系。下面是Arcgis中对北京1954坐标系的说明。
' `3 o9 f* H4 x$ @, V 主要就是以下几个参数:
4 r# B8 k V' L7 i Prime Meridian(起始经度)Datum(大地基准面); " d; k" b; p5 p
D_Beijing_1954 Spheroid(参考椭球体); + n0 q- v& p/ j9 G/ R: `
Krasovsky_1940 (克拉索夫斯基椭球体)。 5 k: C+ }+ ]/ Q8 k1 L
2、投影坐标系
: r' Y) \$ R8 z) u# K6 v8 w 我们在选择坐标系的时候经常会发现以下情况: ( L# G' p2 a1 E* [0 P ]7 S
这一大堆1954坐标系究竟是什么鬼,beijing1954不是地理坐标系吗?为什么投影坐标系里也有?相信懵逼的不止我一个···
; W& P6 ~/ L7 k* c7 n+ `* K. U 首先,投影坐标系的生成是以地理坐标系为基准的,所以每个投影坐标系前面都会挂有地理坐标系。而地理坐标系后面的一串乱七八糟的,则是投影参数!
) t4 v6 j& ~9 Y8 B4 K. a+ n( `5 E GK_Zone与GK_CM的区别* y/ k- K. z! m
例如: 6 Z4 e- [. X. z K9 k! W# `
Beijing 1954 3 Degree GK CM 117E表示三度分带法的北京54坐标系,中央经线在东117度的分带坐标,横坐标前不加带号; + g7 y- [, q6 K3 E
Beijing 1954 3 Degree GK Zone 39表示三度分带法的北京54坐标系,中央经线在东117度的分带坐标,横坐标前加带号; ! g' }+ H' \ j \4 x
下面是Arcgis中对北京1954坐标系的说明:
5 @7 \6 ~% G& }& `9 p9 ~2 G' D4 D 可以看出横坐标前加带号在False—Easting比,不加带号在False—Easting多了两位带号39。 0 u' `2 R9 B5 ~# y1 ?- r7 N
注释:GK 是高斯克吕格,CM 是CentralMeridian 中央子午线,Zone是分带号,N是表示不显示带号。
" T5 q' A$ H; k8 Q: [ 2.1投影
! q: m' n" e/ b H 在地球椭球面和平面之间建立点与点之间函数关系的数学方法,称为地图投影。
/ P" X3 H( u: v1 q) | 地球椭球表面是一种不可能展开的曲面,要把这样一个曲面表现到平面上,就会发生裂隙或褶皱。在投影面上,可运用经纬线的“拉伸”或“压缩”(通过数学手段)来加以避免,以便形成一幅完整的地图。但不可避免会产生变形。 ; S) @8 q, _; h; k N1 M
地图投影的变形通常有:长度变形、面积变形和角度变形。在实际应用中,根据使用地图的目的,限定某种变形。
. F" p! G9 |7 i$ ^4 q/ b$ w 根据不同的需要,我们会选择不同的投影组合!
2 O7 i5 A8 N3 _5 r6 f 按变形性质分类: 8 b" M' y2 f2 g0 G$ @4 S5 v
等角投影:角度变形为零(Mercator) ;
' a: T0 @5 x9 s2 ]; P 等积投影:面积变形为零(Albers) ; * [3 e2 v0 f, }' T$ G) M
任意投影:长度、角度和面积都存在变形 其中,各种变形相互联系相互影响:等积与等角互斥,等积投影角度变形大,等角投影面积变形大。 4 U; O' ], W0 z" s C E
从投影面类型划分:
% _, I4 o8 y; d" u 横圆柱投影:投影面为横圆柱 ; * D' r7 [- i8 ?
圆锥投影:投影面为圆锥 ; # i A1 o9 s0 P4 Q
方位投影:投影面为平面 ;
4 Y, k c8 q. H9 ^0 m) |' C2 u/ U 从投影面与地球位置关系划分为: 2 w$ X# Q) ?! D% O! n
正轴投影:投影面中心轴与地轴相互重合 ;
; I0 ]' Z4 q, x8 e 斜轴投影:投影面中心轴与地轴斜向相交 ;
6 T) F, D. S* @0 |( b8 m% T! R8 h) H 横轴投影:投影面中心轴与地轴相互垂直 ; " E% d( G# P% ]4 m+ ]/ a
相切投影:投影面与椭球体相切 相割投影:投影面与椭球体相割。 , b! {. m( N- p
2.2我国常用投影5 u, q& f5 Q# l" @& M
2.2.1高斯-克吕格投影: m9 t% K* c! y# x$ F0 a$ P! n
我国基本比例尺地形图(1:100万、1:50万、1:25万、1:10万、1:5万、1:2.5万、1:1万、1:5000)除1:100万以外均采用高斯-克吕格Gauss-Kruger投影(横轴等角切圆柱投影)为地理基础。
, b! e0 z* h% p5 R7 ]3 p 高斯克吕格投影的特点:
) ]% ], Z8 w2 I% n( h 横轴等角切圆柱投影:
: j1 b. l; R( Y' l# P. f6 \ – 离开中央子午线越远,变形越大 $ A- p, C8 Q, J3 Z) K6 w
– 赤道是直线,离开赤道的纬线是弧线,凸向赤道
, l( H j7 r9 P: Y& \# K – 没有角度变形 – 长度和面积变形很小
1 v: g @- G1 H$ d 北京54和西安80投影坐标系的投影方式高斯投影特点:
$ c0 w9 G% i' v9 d3 ^- p – 中央子午线长度变形比为1 / q: p( E E2 H/ e* S' R' ~
– 在同一条经线上,长度变形随纬度的降低而增大,在赤道处为最大
! x# e# b) \7 ?& I X – 在同一条纬线上,长度变形随经差的增加而增大,且增大速度较快 9 o v6 J5 h! g4 [6 k
我们经常会听到6°分带,3°分带的说法。其实并不是所有投影都有分带,从下面一张图就可以看出,分带是高斯克吕格投影自带的。 " H7 u! d) z( V* t1 X+ f9 |
高斯-克吕格投影分带规定:该投影是国家基本比例尺地形图的数学基础,为控制变形,采用分带投影的方法,在比例尺1:2.5万—1:50万图上采用6°分带,对比例尺为1:1万及大于1:1万的图采用3°分带。 1 S3 N+ c5 }! T {7 z8 Z5 M
6°分带法:从格林威治零度经线起,每6°分为一个投影带,全球共分为60个投影带,东半球从东经0°—6°为第一带,中央经线为3°,依此类推,投影带号为1—30。其投影代号n和中央经线经度L0的计算公式为:L0=(6n—3)°;西半球投影带从180°回算到0°,编号为31—60,投影代号n和中央经线经度L0的计算公式为L0=360—(6n—3)°。
9 Q- m+ ^+ a/ j. a& X 3°分带法:从东经1°30′起,每3°为一带,将全球划分为120个投影带,东经1°30′—4°30′,…178°30′—西经178°30′,…1°30′—东经1°30′。
6 Z0 a6 }: i6 F* u1 j- M) X 东半球有60个投影带,编号1—60,各带中央经线计算公式:L0=3°n,中央经线为3°、6°…180°。西半球有60个投影带,编号1—60,各带中央经线计算公式:L0=360°—3°n,中央经线为西经177°、…3°、0°。
# t+ \- C3 o. E 为了便于地形图的测量作业,在高斯-克吕格投影带内布置了平面直角坐标系统,具体方法是,规定中央经线为X轴,赤道为Y轴,中央经线与赤道交点为坐标原点,x值在北半球为正,南半球为负,y值在中央经线以东为正,中央经线以西为负。由于我国疆域均在北半球,x值均为正值,为了避免y值出现负值,规定各投影带的坐标纵轴均西移500km,中央经线上原横坐标值由0变为500km。为了方便带间点位的区分,可以在每个点位横坐标y值的百千米位数前加上所在带号。 ) w" I% v( k0 j& S- b: [* W6 I
2.2.2其他投影
# T/ Q0 [) i$ w% B" o' C" x 1:100万地形图采用兰伯特Lambert投影(正轴等角割圆锥投影),其分幅原则与国际地理学会规定的全球统一使用的国际百万分之一地图投影保持一致。海上小于50万的地形图多用墨卡托Mercator投影(正轴等角圆柱投影)。 6 o( O+ J; M4 r. N* {
3、不同坐标之间的转换
: D& t) ]8 ^9 t- J9 j 坐标转换在测绘工作中经常会用到,特别是在处理原始数据的时候。在这里以COORD坐标转换软件为例,介绍一下两个方面:
8 C. t3 C5 E) T5 \, \ 1、如何进行两个坐标椭球系之间的数据转换,也就是求取七参数; - f/ t0 J3 _$ S- k- A) W+ g
2、经纬度转平面坐标的操作方法。
: x1 u+ V# o8 Q. M 下面演示一下该软件进行坐标转换的方法: 0 k) p q, F/ r3 y2 B
3.1、坐标系之间的转换 4 l$ {( B' L8 K, o! r* f# ^& ?6 l7 e
通过三个以上已知点计算七参数时的一些操作。
W: ], F: r1 V& `8 j$ U/ ~2 W& \: ]; b 图3.1:软件界面0 R2 g6 p# ^0 L# ]" b
! k# `$ H% S- h. c- f" w! q# M! U 3.1.1 参数的分类
: i; W: }/ y$ ?" k- _8 e7 X 3.1.1.1 三参数 3 t, a4 ] n* h' r
即 X 平移, Y 平移, Z 平移只需一个已知点即可。适用于小范围内使用,不同椭球间可以互转,转换后的坐标系方向与源坐标系方向一致,只有当源坐标系与当前坐标系的方向一致或在精度范围内才可使用。 9 k+ X. u% T) n2 _
3.1.1.2 四参数
3 c5 w* o+ B' `9 b! u 即 X 平移, Y 平移,坐标旋转,投影比例,需两个已知点,在高程精度不高的情况下使用,因为在四参数中,没有高程改正参数,在实际的测量中,会加上三参数中的高程改正参数或利用高程拟后来得到较高精度的高程。 4 q5 d- K4 k' o+ l
四参数属于同一椭球下的转换,当源坐标系与当前坐标系不一致时,如将WGS84经纬度转换到北京54坐标系,一般会利用以下的转换原理。
/ V, ^8 {3 c$ v0 ]! @4 u- N8 ]' Q 由于四参数中,没有高程改正,如果需要高程精度较高的用户,需再计算高程改正参数,高程改参数根据所使用已知点的数量又分为多种参数:
7 X6 S ^2 I e% M* s! j 1~2个高程已知点时,即Z平移,在使用中即为三参数中的Z平移参数; 1 f- A# O$ J+ y. U( A& h& ], F3 e W
3~4个高程已知点时,采用高程拟合中的平面拟合参数;
7 ]& y- o- \% n 6~7个高程已知点时,采用高程拟合中的曲面拟合参数。 ' u' p' J! [9 I1 K' H
在一般使用过程中,高程拟合参数可根据输入已知点的个数自动或手动选择计算各参数。 6 M* B$ n! N+ C$ m
3.1.1.3 七参数
" i* R5 |2 j! B2 c, d' w" x. H X 平移, Y 平移, Z 平移, X 旋转(WX), Y 旋转(WY), Z 旋转(WZ),尺度变化(DM )。需要三个点以上。
) Y" K. y* `* ~/ w, g 七参数属于不同椭球下的转换,适用于大范围,一般RTK的使用中,在做完控制静态测量后,可直接使用静态平差结果里的数据进行参数的计算,部分软件会旋转角度,比如南方RTK所使用的“工程之星”中,对于七参数就要求旋转角度不能大于10秒,否则只能用四参数+高程拟合参数。
* [2 o3 n4 o0 i# r9 I4 X }& r9 K9 P; M 3.2、COORD实例
' d* w H, o0 |7 Z 在这里,我们以国家80坐标系和北京54坐标系进行参数的计算,它们是属于不同椭球系间的转换,在不同的椭球之间的转换都是不严密,因此不存在一套转换参数可以全国通用的,在每个地方会不一样,因为它们是两个不同的椭球基准。那么,两个椭球间的坐标转换,一般而言比较严密的是用七参数布尔莎模型,即 X 平移, Y 平移, Z 平移, X 旋转(WX), Y 旋转(WY), Z 旋转(WZ),尺度变化(DM )。要求得七参数就需要在一个地区需要 3 个以上的已知点。 & C/ u' F& B# c B3 t
3.2.1 坐标系共点坐标数据 ! }9 H5 Q1 A; \/ o, o
表3.1:示例数据
9 K. y9 \( ~2 f" ], y1 P1 X+ i- F % |+ h) ~5 u6 O. D: V
3.2.2 投影参数
. `: b O, U& d6 n: | 投影方式:高斯-克吕格3度带;中央子午线:117度。 7 A6 O' l& _/ v: W. j, w" ^# J! n! g
3.2.3 根据七参数
* Q2 [6 z+ a9 s6 V! T3 W 2.3.1求北京54 坐标系--> 国家80坐标系的七参数 / \1 z2 r/ K, x7 v1 U5 F
根据两个已有坐标(例3个)求七参数;在操作之前必须先将中央子午线必须先输好。 6 l% b$ V. o2 E! [
图3.2:投影参数设置 8 T2 l. ~" T- }1 v+ t
输入两个坐标系三个或三个以上相对应的坐标,输入完后点击“计算”可以看到其中的点位精度,当点击“确定”时,软件将自动将七参数设置好。 6 M6 X' B- S; p( [6 k B4 s
图3.3:已知数据的输入 & \6 p( F# `2 c2 q6 ^8 S' ]! {; e
注意:请注意椭球的选择顺序,如果我们要计算北京54->国家80,在此软件中,这里的源坐标应该是国家80坐标系,反之亦反。 4 Y0 d0 D+ N( w/ \; m1 n
接下来就是选中“七参数转换”,设定转换坐标系。可以查看输出的值,这个值一定就是我们所需的,我们可以返回我的刚计算的七参数,看是否正确。 5 q8 H+ a% l A
图3.4:单点转换1 t' J7 G! Q+ Z t. b' `
& ?, ~' Q4 {. j4 a# ?+ q 检查数据,发现正确无误,可以使用;在计算过程中,我们可以随时保存我们的设置信息,下次使用这个转换时不用再输入七参数;也可以先记下七参数,用时再输入设置。 . \0 ]$ g# W; y0 o; M
这里我们只用到了单点转换,当数据比较多时,我们可以使用“文件转换”功能先将我们需转换的数据输入到记事本中,并给点点号,确认XYZ的位置,在这一项里,我们可以按刚输入记事本的的格式,在此定义好。此功能比较简单,可以自己看看。
5 A% i, S p0 \, k: K5 \* H$ q 3.2.4、国家80 --> 北京54
; `- u' _; w7 r2 i 操作步骤和(北京54 --> 国家80)差不多一致,只是在计算“七参数时”顺序相反,当我们要将国家80为北京54时,输入的源坐标应该是54坐标系,所以在这里,大家一定要注意。 8 s l5 N& e- m6 q
图3.5:已知数据的输入
" I8 l" I' Q' W% w$ K; n+ |, n 图3.6 单点计算
|9 d X. J$ h, E5 t 3.3、经纬度转成平面坐标 2 Z, G3 X( h# M% ?' |9 ^
3.3.1、使用方法 % b) C0 M' }% N1 T" `9 ^
首先打开软件,“设置”——“地形投影”——设置是属于3度带还是6度带,输入**经度,操作过程如下图:6 V3 N2 N- ?! v- W, F
8 ?% a/ | K" m
图3.6 软件界面 6 L0 y, W! ]& E7 G' Y- T* h' U
图3.7 打开地图投影
2 J) L: Y1 \# j3 O( Q6 @ 图3.8 选择相应的六度带或3度带,输入**经度 " D" p3 a/ d% ^+ g! q& C8 q( D
3.3.2、单点转换的方法 ' B; g. C5 B: |! S3 ^+ \
选择单点转换,设置好坐标类型,以及转换到的坐标系,然后输入经纬度,点击转换坐标,就可以得到转换后的坐标。
1 [6 @" j9 J# X 图3.9 单点转换操作
+ O# o- E! @3 Y 3.3.3、批量转换的方法
, { P! y# e" u( U. A) c" g8 r 选择文本转换,设置好坐标类型,以及转换到的坐标系,然后加载坐标文件,点击转换坐标,就可以得到批量转换后的坐标文件。
2 T/ X# ]: Y! E" i9 I6 B G 以上就是航测工作中经常用到的坐标系和坐标系转换问题以及使用COORD坐标转换软件的操作方法,希望对大家有所帮助。 g2 q! C `. t! l# Y% w$ p
; L1 M7 r& ?" g, C( ]3 A2 k, Z) y+ b
8 b( x9 v* J4 q/ I9 c! f. U
3 F9 R, d( S* @
: X9 Z I1 j5 T2 m0 y0 d |