当前位置:首页 >> 电力/水利 >>

投影坐标转换


摘自武测毕业论文(赖增先)——我也是从网上弄来,转载时请务必保留此此 信息。

第二节 平面坐标基准转换
由于海上和陆地上在测量时,使用不同的坐标系和不同参考椭球,而且采用的 投影也不同,使得我们获得的数据不统一,必须进行坐标转换。

§3·2·1 欧拉角
设有两个空间直角坐标系,分别为 O-XYZ 和 O-X '

; Y ' Z ' ,为了便于讨论其相应 坐标轴间的变换,设其原点相同如图所示,选择 ε
x

、ε y 、ε z 为欧拉角,又称旋转参

数,经过三次旋转,使两个坐标系重合,既: (图见下页 A) 首先,绕 O Z ' 轴,将 O X ' 轴旋转到 OX
0

轴,所转的角为 ε z ;

其次,绕 OY 0 轴,将 O Z ' 轴旋转到 OZ 0 轴,所转的角为 ε y ; 最后,绕 OX 轴,将 O Z
0

轴旋转到 OZ 轴,所转的角为 ε x ; Z Z
0

Z'

X' X X
0

O

Y

0

Y Y'

图A 因此有 X Y Z = R 1 ( ε x )R 2 ( ε y )R 3 ( ε z ) X' Y' Z'

式中 R 1 ( ε x ) 2 ( ε y ) 3 ( ε z )为旋转矩阵,其表达式在 ε 、 ε y 、 ε z 很 、R 、R 小时可以最终表示为:

X Y Z =

1 -ε z

εz
1 - ε
x

εy

X' Y' Z' 公式 1

ε
1

x

εy

§3·2·2 不同三维空间直角坐标系的变换模型
GPS 测量的 WGS—84 属地心坐标系, 1980 年国家大地坐标系和 1954 年北京坐 而 标系属参心坐标系,他们所对应得空间直角坐标系是不同的,这里将讨论不同空间 直角坐标系的变换模型。 如图 B 两个空间直角坐标系分别为 O-XYZ 和 O ' -X ' Y ' Z ' , 其坐标系原点不同则 存在三个平移参数 ? X 0 、 ? Y 0 、 ? Z 0 ,他们表示 O ' - X ' Y ' Z ' 坐标系原点 O ' 相对于 O-XYZ 坐标系原点 O 在三个坐标轴上的分量;又当各坐标轴相互不平行时,既存在 三个旋转参数 ε
x

、ε y 、ε z 。 Z O' X'

Z'

Y' O Y

X 考虑到两个坐标系的平移和旋转以及尺度参数可得公式如下: X Y Z =(1+m) X' Y' Z'
?X0

1 -ε z

εz
1 - ε
x

εy

X' Y' Z'

ε
1

x

εy

+

?Y0

公式一

?Z 0

式中共有七个变换参数 ? X 0 、 ? Y 0 、 ? Z 0 、ε

x

、ε y 、ε z 、m,简称此公式为布

尔莎七参数变换公式,是坐标变换中一个非常重要的公式。七参数变换公式,除了 布尔莎公式外,还有莫洛琴斯基公式和范氏公式。这三种公式,它们之间的七个参 数相差很大,但各自构成完整的数学模型,参数间存在着明确的解析关系,可以相 互间转换。分别用它们来换算点的坐标时,其结果是完全相同的。因此,这三个公 式是等价的。我国的地心坐标变换参数地心二号是七个变换参数,既采用布尔莎公 式。 当公式一中 ε
x

= ε y = ε z =m=0,既称之为三参数公式。三参数公式表明两个空间

直角坐标系尺度一致,且两个坐标轴相互平行。我国地心坐标变换参数地心一号系 三个变换参数。同理在公式一中,略去某些参数,可分别得到四参数、五参数、六 参数等坐标变换参数。公式一中的变换参数,一般利用公共点上的两套空间直角坐 标系坐标值(X,Y,Z) i 和(X
'

,Y ' , Z ' ) i 即可采用最小二乘法解得。

应该指出,当进行两种不同空间直角坐标系变换时,坐标变换的精度除取决于 坐标变换的数学模型和求解变换参数的公共点坐标精度外,还和公共点的多少、几 何形状结构有关。鉴于地面网可能存在一定的系统误差,且在不同区域并非完全一 样,所以采用分区变换参数,分区进行坐标转换,可以提高坐标变换精度。无论是 从我国的多普勒网还是 GPS 网,利用布尔莎公式求解和地面大地网间得变换参数, 分区变换均较明显地提高了坐标变换的精度。

§3·2·3 不同三维大地坐标系的变换模型
对于不同的三维大地坐标系的变换模型,除了上节的七个变换参数外,还应增 加两个变换参数, ,这就是两个大地坐标系所对应的地球椭球参数的不同。不同大地 坐标的变换公式,又称大地坐标微分公式或变换椭球微分公式。当包括旋转参数和 尺度参数时,称为广义大地坐标微分公式或广义变换椭球微分公式。 空间一点的空间直角坐标与大地坐标关系式是: X (N+H)cosBcosL Y Z = (N+H)cosBsinL [N (1-e 2 )+H]sinB 公式二

式中 N 为卯酉圈曲率半径。在这个公式中当已知 L,B,H 时,求 X,Y,Z 是非常简 单的,只要代入公式即可。当已知 X,Y,Z 时反求 L,B,H 则可以采用直接解或迭 代解法,解算时对公式做些变化即可。 由公式二最终我们可以得到不同三维大地坐标系的变换公式;

dL dB =

-

sin L ρ '' ( N + H ) cos B
sin B cos L '' ρ M +H

cos L ρ '' ( N + H ) cos B sin B sin L '' ρ N+H

0
cos B ρ '' M +H

?X0

?Y0

+

dH

cosBcosL

cosBsinL

sinB

?Z 0

N (1 ? e 2 ) + H tgB cos L N+H

N (1 ? e 2 ) + H tgB sin L N+H ( N + H ) ? Ne 2 sin 2 B cos L M +H Ne 2 sin B cos B cos L

-1

ε

x

-

( N + H ) ? Ne 2 sin 2 B sin L M +H Ne 2 sin B cos B sin L

0

εy

+

-

ρ"
0

ρ"
0

0

εz

-

N e 2 sin B cos Bρ " M +H

m+

N e 2 sin B cos Bρ " ( M + H )a

N+H-Ne 2 sin 2 B 0
M (2 ? e 2 sin 2 B) sin B cos Bρ " ( M + H )(1 ? f )
M (1 ? e 2 sin 2 B ) sin 2 B 1? f

-

N (1 ? e 2 sin 2 B ) a

da 公式三 df

式中 dL、dB 以弧度秒为单位,等式右端 L、B、H 均以换算前坐标值代入。公式 三也就是顾及七个参数和椭球大小变化的广义大地坐标微分公式或广义变换椭球微 分公式。略去旋转参数和尺度变化参数的影响,即为一般的大地坐标微分公式或椭

球微分公式。 利用一些公共点上两套大地坐标系坐标值,采用最小二乘法可解得变换参数。

§3·2·4 不同两维大地坐标系的变换模型
在三维不同大地坐标系的变换模型中, 当进行 WGS—84 和我国参心大地坐标系的 变换时,由于后者大地高的精度不高(一般在 3m 左右的误差) ,加之又难以确定其 方差和协方差,因此,也可以考虑选择二维大地坐标系的变换模型。 所谓二维大地坐标系,即当大地高 H 为零时的椭球面上的大地坐标系。其变换 模型,只要在公式三中,将 H=0 代入即可得到。 将此公式用于 GPS 网和地面网联合平差时,如果顾及地面网的系统性观测误差 对网的定向的影响时, 则可在椭球面上引入一个附加旋转参数 dA,以使两网更好的配 合。由于 dA 产生的对 dL、dB 的影响加于公式右端。 为了便于比较 GPS 网和地面网的大地坐标,若在将 GPS 网的 X、Y、Z 反算 L、B、 H 时,采用了地面网的椭球参数,即两网相应的椭球参数已化为一致,则公式中不再 有 da、df 项。

§3·2·5 不同二维高斯投影平面坐标系的坐标转换
由高斯投影正算公式可得: ?x ?x dx dL ?L ?B ?y ?y dy = dB ?L ?B 式中等号右端偏导数由高斯投影正算公式得: ?x =NsinBcosBl ?L ?x 1 =M[1+ (1-2sin 2 B)l 2 ] ?B 2 ?y 1 =NcosB[1+ (1-2sin 2 B) l 2 ] ?L 2 ?y =-MsinBl ?B

公式四

公式五

上式中,l=L-L 0 , L 0 为中央子午线得大地经度,公式五中 dB、dL 见公式三。对 于不同二维高斯投影平面坐标系坐标差的模型,可以由公式三和公式四给出。

§3·2·6

同一参考系统下的高斯直角坐标、

大地坐标、空间直角坐标之间的相互转换

当运用了§3·2·1~§3·2·5 后,我们就可以将一个系统的坐标,转化到另 一个系统的对应结果。可以完成对应之间的坐标转化,但是,如果高斯直角坐标、 大地坐标、空间直角坐标之间的相互转换就需要用本节的内容。 3·2·6·1 高斯直角坐标同大地坐标之间的转化 完成两者之间的转化,要用高斯正、反算方法。这里并不详细介绍,只给出其数 学模型: 1 1 N*t(5- t 2 +9 η 2 +4 η 4 )cos 4 B*l 4 正算公式:x=X+ N*t*cos 2 B*l 2 + 2 24 1 + N*t(61-58 t 2 + t 4 +270 η 2 -330 η 2 t 2 )cos 6 * l 6 720 1 1 y=N*cosB*l+ N(1-t 2 + η 2 )cos 3 B* l 3 + N(5-18 t 2 6 120 + t 4 +14 η 反算公式:B=B f tf 2M f N f
2

-58 η -y
2

2

t 2 )cos 5 B*l 5

+24

tf 24M f N
3 f

(5+3t 2f + η 2 -9t 2f η 2 )y 4 f f

-

tf 720M f N
5 f

(61+90t 2f +45t 4f )y 6

l=

1 1 y(1+2t 2f + η 2 )y 3 f 3 N f cos B f 6 N f cos B f
+ 1 (5+28t 2f +24t 4f +6 η 2 +8t 2f η 2 ) y 5 f f 120 N cos B f
5 f

通过高斯正、反算,可以将高斯直角坐标化算为大地坐标,或将大地坐标化算为 高斯直角坐标。 3·2·6·2 大地坐标同空间直角坐标的化算 空间一点的空间直角坐标与大地坐标关系式是: X (N+H)cosBcosL Y Z = (N+H)cosBsinL [N (1-e 2 )+H]sinB

式中 N 为卯酉圈曲率半径。在这个公式中当已知 L,B,H 时,求 X,Y,Z 是非常简 单的,只要待入公式即可。当已知 X,Y,Z 时反求 L,B,H 则可以采用直接解或迭 代解法,解算时对公式做些变化即可。 3·2·6·3 高斯直角坐标与空间直角坐标的相互转换

关于这两者的转换关系并没有直接给出,但是,我们可以利用两者同大地坐标 的转换关系,先把其中的一种化算成大地坐标,然后再由大地坐标转换成另一种坐 标。

§3·2·7 不同投影之间的转换
陆地地形图采用高斯投影系统,而海图采用的是墨卡托投影系统,这使得两种 地形图不能拼接。在此以由高斯投影系统转化到墨卡托投影系统,给出转化的思路 和过程: 首先利用布尔莎七参数变换公式 布尔莎七参数变换公式(公式一)将一参考系内的高斯坐标转换为另 一参考系内的高斯坐标;再利用高斯反算公式,将转换后的高斯直角坐标换算为大 高斯反算公式, 高斯反算公式 地坐标;然后在同一参考系内进行墨卡托投影 墨卡托投影将大地坐标化算为该平面内的平面坐 墨卡托投影 标。这里考虑的是高斯坐标同墨卡托平面坐标不在同一参考椭球上,如果在同一参 考椭球内,则第一步化算可以省略。


相关文章:
坐标系投影方式的选择及坐标转换
坐标系投影方式的选择及坐标转换_天文/地理_自然科学_专业资料。坐标系投影方式的选择及坐标转换 [摘要]通过对几种常用投影方式的分析对比,详细剖述了海外项目投影...
Arcgis地图投影和坐标转换方法
Arcgis 地图投影坐标转换方法 1、动态投影(ArcMap) 所谓动态投影指, ArcMap 中的 Data 的空间参考或是说坐标系统是默认为第一加 载到当前工作区的那个文件的...
利用ArcGIS进行地图投影和坐标转换的方法
利用ArcGIS进行地图投影坐标转换的方法_计算机软件及应用_IT/计算机_专业资料。利用 ArcGIS 进行地图投影坐标转换的方法 3、投影变换(ArcToolBox) 上面说了这么多...
投影坐标转换
投影坐标转换_天文/地理_自然科学_专业资料。投影坐标转换如何在 ArcView(ArcGIS)中将从 GPS 上获取的经纬度坐标进行投影转化?很简单,这个工作 流程是: 1.将 GPS...
利用ArcGIS进行地图投影和坐标转换的方法
利用ArcGIS 进行地图投影坐标转换的方法 1,动态投影(ArcMap) 所谓动态投影指, ArcMap 中的 Data 的空间参考或是说坐标系统是默认为第一加载到当前 工作区的那个...
坐标转换与投影浅析
坐标转换投影浅析_天文/地理_自然科学_专业资料。坐标转换投影浅析 摘要: 本文主要介绍不同坐标系下测绘成果进行高程和平面坐标系之间转换 的基本方法、地图投影...
ARCGIS中坐标转换及地理坐标、投影坐标定义
ARCGIS 中坐标转换及地理坐标、投影坐标定义 1、动态投影(ArcMap) 所谓动态投影指, ArcMap 中的 Data 的空间参考或是说坐标系统是默认为 第一加载到当前工作区的...
GIS系统中的常用坐标转换和地图投影功能实现
GIS 系统中的常用坐标转换和地图投影功能实现及其应用 摘要:本文论述了现实世界中的经纬度坐标和 GIS 系统的图形、设备坐标之间的转换,结 合在通信设备网管系统中的...
ArcGIS中坐标转换及地理坐标、投影坐标的定义
ARCGIS 中坐标转换及地理坐标、投影坐标的定义 1.ARCGIS 中坐标转换及地理坐标、投影坐标的定义 1.1 动态投影(ArcMap) 所谓动态投影指,ArcMap 中的 Data 的空间...
坐标转换与投影浅析
龙源期刊网 http://www.qikan.com.cn 坐标转换投影浅析 作者:张劲乾 陈志华 来源:《城市建设理论研究》2012 年第 35 期 摘要:本文主要介绍不同坐标系下测绘...
更多相关标签:
arcgis投影坐标转换 | 经纬高斯投影坐标转换 | 投影坐标转换地理坐标 | arcgis投影坐标系转换 | 高斯投影坐标转换软件 | 投影坐标系转换 | 高斯投影坐标转换 | 投影坐标转换成经纬度 |