大地坐标与直角空间坐标转换计算公式
空间直角坐标系、大地坐标系、平面坐标系、高斯平面直角坐标系

空间直角坐标系、大地坐标系、平面坐标系、高斯平面直角坐标系本篇学习了空间直角坐标系、大地坐标系、平面坐标系、高斯平面直角坐标系。
这个个坐标系有时很容易弄混淆!(一)空间直角坐标系空间直角坐标系的坐标原点位于参考椭球的中心,Z轴指向参考椭球的北极,X轴指向起始子午面与赤道的交点,Y轴位于赤道面上切按右手系于X轴呈90度夹角,某点中的坐标可用该点在此坐标系的各个坐标轴上的投影来表示。
空间直角坐标系可用如下图所示:(二)大地坐标系大地坐标系是采用大地纬度、经度和大地高程来描述空间位置的。
纬度是空间的点与参考椭球面的法线与赤道面的夹角;经度是空间的点与参考椭球的自转轴所在的面与参考椭球的起始子午面的夹角;大地高程是空间的点沿着参考椭球的法线方向到参考椭球面的距离。
地面点的高程和国家高程基准(1)绝对高程。
地面点沿垂线方向至大地水准面的距离称为绝对高程或称海拔。
过去我国采用青岛验潮站(tide gauge station)1950~1956年观测成果求得的黄海平均海水面作为高程的零点,称为“1956年黄海高程系”(Huanghai height system 1956水准原点高程为72.289m)。
后经复查,发现该高程系的验潮资料时间过短,准确性较差,改用青岛验潮站1950~1979年的观测资料重新推算,并命名为“1985年国家高程基准”(Chinese height datum 1985)。
国家水准原点(leveling origin高程为72.260m)设于青岛市观象山附近,作为我国高程测量的依据。
它的高程值是以“1985年国家高程基准”所确定的平均海水面为零点测算而得。
在使用原“1956年黄海高程系”的高程成果时,应注意将其换算为新的高程基准系统。
(2)相对高程。
地面点沿铅垂线方向至任意假定的水准面的距离称为该点的相对高程,亦称假定高程。
在图l—5中,地面点A和B的相对高程分别为H'A 和H'B 。
大地坐标与直角空间坐标转换计算公式

⼤地坐标与直⾓空间坐标转换计算公式⼤地坐标与直⾓空间坐标转换计算公式⼀、参⼼⼤地坐标与参⼼空间直⾓坐标转换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* cosLY =(N +H )* cosB* sin L ?Z =[N * (I _e2) +H]* sin B”公式中,N为椭球⾯卯⾣圈的曲率半径,e为椭球的第⼀偏⼼率,a、b椭球的长短半径,f椭球扁率,W为第⼀辅助系数a2 -b22* f -1e 或e =a fW = . (1 -g*sin2BN aW西安80椭球参数:长半轴 a=6378140⼟ 5( m)短半轴 b=6356755.2882m扁率a =1/298.2573参⼼空间直⾓坐标转换参⼼⼤地坐标Z* (N + H) (X2 Y2)* N* (1 -e2) HX2 Y2cosB⼆⾼斯投影及⾼斯直⾓坐标系1、⾼斯投影概述⾼斯-克吕格投影的条件:1.是正形投影;2.中央⼦午线不变形⾼斯投影的性质: 1.投影后⾓度不变; 2.长度⽐与点位有关,与⽅向⽆关;3.离中央⼦午线越远变形越⼤为控制投影后的长度变形,采⽤分带投影的⽅法。
常⽤3度带或6度带分带,城市或⼯程控制⽹坐标可采⽤不按 3度带中央⼦午线的任意带。
空间大地坐标系及平面直角坐标系转换公式

§2.3.1 坐标系的分类正如前面所提及的,所谓坐标系指的是描述空间位置的表达形式,即采用什么方法来表示空间位置。
人们为了描述空间位置,采用了多种方法,从而也产生了不同的坐标系,如直角坐标系、极坐标系等。
在测量中常用的坐标系有以下几种:一、空间直角坐标系空间直角坐标系的坐标系原点位于参考椭球的中心,Z 轴指向参考椭球的北极,X 轴指向起始子午面与赤道的交点,Y 轴位于赤道面上且按右手系与X 轴呈90°夹角。
某点在空间中的坐标可用该点在此坐标系的各个坐标轴上的投影来表示。
空间直角坐标系可用图2-3来表示:图2-3 空间直角坐标系二、空间大地坐标系空间大地坐标系是采用大地经、纬度和大地高来描述空间位置的。
纬度是空间的点与参考椭球面的法线与赤道面的夹角;经度是空间中的点与参考椭球的自转轴所在的面与参考椭球的起始子午面的夹角;大地高是空间点沿参考椭球的法线方向到参考椭球面的距离。
空间大地坐标系可用图2-4来表示:图2-4空间大地坐标系三、平面直角坐标系平面直角坐标系是利用投影变换,将空间坐标空间直角坐标或空间大地坐标通过某种数学变换映射到平面上,这种变换又称为投影变换。
投影变换的方法有很多,如横轴墨卡托投影、UTM 投影、兰勃特投影等。
在我XX 用的是高斯-克吕格投影也称为高斯投影。
UTM 投影和高斯投影都是横轴墨卡托投影的特例,只是投影的个别参数不同而已。
高斯投影是一种横轴、椭圆柱面、等角投影。
从几何意义上讲,是一种横轴椭圆柱正切投影。
如图左侧所示,设想有一个椭圆柱面横套在椭球外面,并与某一子午线相切〔此子午线称为中央子午线或轴子午线〕,椭球轴的中心轴CC ’通过椭球中心而与地轴垂直。
高斯投影满足以下两个条件:1、 它是正形投影;2、 中央子午线投影后应为x 轴,且长度保持不变。
将中央子午线东西各一定经差〔一般为6度或3度〕X 围内的地区投影到椭圆柱面上,再将此柱面沿某一棱线展开,便构成了高斯平面直角坐标系,如以下图2-5右侧所示。
直角坐标系和大地坐标系转换

直角坐标系和大地坐标系的转换
在地理信息系统和测量领域中,直角坐标系和大地坐标系是两种常用的坐标系统。
直角坐标系是平面直角坐标系,由水平的x轴和垂直的y轴构成,可以用来表示平面上的点的位置,通常以米为单位。
而大地坐标系则是一种用来描述地球上点的位置的坐标系统,通常是经度(Longitude)和纬度(Latitude)的组合。
直角坐标系到大地坐标系的转换
直角坐标系到大地坐标系的转换涉及到高等数学的知识,主要是利用球面三角学的相关技巧。
在进行转换之前,需要知道点在直角坐标系中的坐标值,以及直角坐标系的原点。
然后,可以通过一系列的数学运算,将点的直角坐标值转换为大地坐标系中的经度和纬度。
大地坐标系到直角坐标系的转换
大地坐标系到直角坐标系的转换相对直接一些。
给定一个点的经度和纬度,我们可以利用地球的半径及球面三角学的相关公式,将该点的经度和纬度转换为直角坐标系中的坐标值。
这种转换可以帮助我们将地球表面上的点的位置转换为平面直角坐标系中的表示,便于进行地理信息系统中的测量和计算。
应用
直角坐标系和大地坐标系的转换在地理信息系统、地图制作、导航系统等领域都有着重要的应用。
通过这种转换,我们可以方便地将地球上的点的位置在不同坐标系统之间进行转换,从而实现不同系统之间的数据交换和信息共享。
总的来说,直角坐标系和大地坐标系的转换是地理信息系统和测量领域中的重要技术,对于地球表面上点的位置的表示和计算具有重要意义,能够为人类的地理信息分析和决策提供便利。
常用坐标系之间的关系与转换

7.5 常用坐标系之间的关系与转换一、大地坐标系和空间大地直角坐标系及其关系大地坐标系用大地纬度企丈地经度L 和丈地髙H 来表示点的位置°这种坐标系是经 典大地测量甬:両用座标紊7屜据地图投影的理论,大地坐标系可以通过一定的投影转 化为投影平面上的直角坐标系,为地形测图和工程测量提供控制基础。
同时,这种坐标系 还是研究地球形状和大小的 种有用坐标系°所以大地坐标系在大地测量中始终有着重要 的作用.空间大地直角坐标系是-种以地球质心为原点购亘墮®坐标系,一般用X 、化Z 表 示点BSSTSTT 逐碇SS 範菇飞両H 绕禎扭转冻其轨道平面随时通过 地球质心。
对它们的跟踪观测也以地球质心为坐标原点,所以空间大地直角坐标系是卫星 大地测量中一种常用的基本坐标系。
现今,利用卫星大地测量的手段*可以迅速地测定点的空间大地直角坐拯,广泛应用于导航定位等空间技术。
同时经过数学变换,还可求岀点 的大地坐标I 用以加强和扩展地面大地网,进行岛屿和洲际联测,使传统的大地测量方法 发生了深刻的变化,所以空间大地宜角坐标系对现今大地测量的发展’具有重要的意义。
、大地坐标系和空间大地直角坐标系的转换如图7- 23所示’尸点的位置用空间 大地直角坐标〔X, Y, Z)表示,其相应 的大地坐标为(E, L)a 将该图与图?一5上式表明了 2种基本坐标系之间的关系。
加以比较可见,图7-5中的子午椭圆平面 相当于图7-23中的OJVP 平面.其中 PPz=Z.相当于图7-5中的j7;OP 3相当 丫于图7-5中的仏两平面的经度乙可视为相同,等于"叽 于是可以直接写岀X=jrcQsi f Y=jrsinL, Z=y将式(7-21).式(7-20)分别代入上式, 井考虑式(7-26)得X=Ncos^cosZr ”Y =NcQsBsinL > (7—78)Z=N (1—护〉sin^ ;BB 7-231.由大地坐标求空间大地直角坐标当已知椭球面上任一点P 的大地坐标(B, L)时,可以按式(7-78)直接求该点的 空间大地直角坐标(X, Y, Z)。
高斯平面直角坐标与大地坐标的相互转换——高斯投影的正算与反算.

昆明冶金高等专科学校测绘学院 (4)计算公式
3 2 2 2 4 ( 5 3 t 9 t ) y f f f f 2M f N f 2 4M f N 3 f tf 2 4 6 (6 1 9 0t f 4 5t f ) y 7 2 0M f N 5 f 1 1 2 2 3 l y (1 2t f f ) y 3 N f co s B f 6 N f co s B f 1 2 5 (5 2 8t 2 t4 2 2 f 24 f 6 f 8 f t f )y 5 1 2 0N f co s B f B Bf tf y2 tf
式中:
2 e 2 cos2 B
t 2 tan2 B l (L L0) X为B对应子午线弧长 N为卯酉圈曲率半径 20626 5
昆明冶金高等专科学校测绘学院
2
高斯投影坐标反算公式
(1)高斯投影反算:
已知某点 x, y ,求该点 L, B ,即 x, y ( L, B) 的坐标变换。 (2)投影变换必须满足的条件
昆明冶金高等专科学校测绘学院
二、高斯投影坐标正反算得实用公式及算例
1 高斯投影坐标正算公式 (1)高斯投影正算: 已知某点的 L, B ,求该点的 x, y ,即 (2)投影变换必须满足的条件: 中央子午线投影后为直线; 中央子午线投影后长度不变; 投影具有正形性质,即正形投影条件。 (3)投影过程 在椭球面上有对称于中央子午线的两点 P1 和 P2 ,它们的大地坐标 分别为 ( L1 , B1 )或(l1 , B1)及 (L2 , B2)或(l2 , B2 ) 式中 l 为椭球面上点的经 度与中央子午线 ( L0 ) 的经度差:l L L0 ,点在中央子午线之东, l 为正,在西则为负,则投影后的平面坐标一定为P1 ( x1 , y1 ) 和 P2 ( x 2 , y 2 ) 。
空间直角坐标系与空间大地坐标系的相互转换及其C++源程序
空间直角坐标系与空间大地坐标系的相互转换1.空间直角坐标系/笛卡尔坐标系坐标轴相互正交的坐标系被称作笛卡尔坐标系。
三维笛卡尔坐标系也被称为空间直角坐标系。
在空间直角坐标系下,点的坐标可以用该点所对应的矢径在三个坐标轴上的投影长度来表示,只有确定了原地、三个坐标轴的指向和尺度,就定义了一个在三维空间描述点的位置的空间直角坐标系。
以椭球体中心O为原点,起始子午面与赤道面交线为X轴,在赤道面上与X轴正交的方向为Y轴,椭球体的旋转轴为Z轴构成右手坐标系O.XYZ,在该坐标系中,P点的位置用X,Y,Z表示。
在测量应用中,常将地球空间直角坐标系的坐标原点选在地球质心(地心坐标系)或参考椭球中心(参心坐标系),z轴指向地球北极,x轴指向起始子午面与地球赤道的交点,y轴垂直于XOZ面并构成右手坐标系。
空间直角坐标系2.空间大地坐标系由于空间直角坐标无法明确反映出点与地球之间的空间关系,为了解决这一问题,在测量中引入了大地基准,并据此定义了大地坐标系。
大地基准指的是用于定义地球参考椭球的一系列参数,包括如下常量:2.1椭球的大小和形状2.2椭球的短半轴的指向:通常与地球的平自转轴平息。
2.3椭球中心的位置:根据需要确定。
若为地心椭球,则其中心位于地球质心。
2.4本初子午线:通过固定平极和经度原点的天文子午线,通常为格林尼治子午线。
以大地基准为基础建立的坐标系被称为大地坐标系。
由于大地基准又以参考椭球为基准,因此,大地坐标系又被称为椭球坐标系。
大地坐标系是参心坐标系,其坐标原点位于参考椭球中心,以参考椭球面为基准面,用大地经度L、纬度B 和大地高H表示地面点位置。
过地面点P的子午面与起始子午面间的夹角叫P 点的大地经度。
由起始子午面起算,向东为正,叫东经(0°~180°),向西为负,叫西经(0°~-180°)。
过P点的椭球法线与赤道面的夹角叫P点的大地纬度。
由赤道面起算,向北为正,叫北纬(0°~90°),向南为负,叫南纬(0°~-90°)。
站心坐标系
X
Zi Z0
Y X 1 Zi
2.4.4 两个空间大地直角坐标系间的转换模型
略去尺度参数和旋转参数的乘积项,上式可进一 步简化为:
X i X 0 X i 0 Z Y X i X i Yi Y0 Yi Z 0 X Yi Yi Zi Z0 Zi Y X 0 Zi Zi
3、空间直角坐标系与站心地平直角坐标系的旋转矢量 之间的关系
若x、y和z为空间坐标系的旋转矢量, x、 y和z为站心坐标系的旋转矢量。顾及旋转矢量是 平移不变量,旋转关系与坐标矢量相同。
X sin B0 cosL0 Y sin B0 sin L0 Z cosB0
sin L0 c os L0
sin B0 cosL0 sin B0 sin L0
cosB0
sin L0 c os L0
0
cosB0 cosL0 x cosB0 sin L0 y
sin B0 z
2.4.3 站心地平坐标系及其应用
由上式得,空间直角坐标系到站心系的变换公式为:
x
X X P0 sin B0 cos L0
X 0 X i 0 Zi Yi X X i Y0 Yi Zi 0 X i Y Yi
Z0 Zi Yi X i 0 Y Zi
上式第二式常用于转换参数未知时,利用同 名点在两个坐标系中的坐标计算转换参数。
2.4.4 两个空间大地直角坐标系间的转换模型
(2). 绕站心系坐标轴的旋转向量有特殊意义
x
y
z 相当于平面控制网间的旋转角。
2.4.3 站心地平坐标系及其应用
(4). 计算卫星的高度角和方位角 卫星Q的方位角和高度角可用其站心坐标xQ、yQ计算。
AQ
tan
高斯正反算及空间直角坐标与大地地理坐标转换
高斯正反算及空间直角坐标与大地地理坐标转换一、实验目的与要求1.对以上理论内容的验证与应用。
2.通过学习掌握测绘软件开发过程与方法,初步具备测绘软件开发基本技能。
3.熟练掌握Visual C++编程环境的使用,了解其特点与程序开发过程,掌软件调试、测试的技术方法。
4.分析测绘程序设计技术课程中相关软件的结构和模块功能,掌握结构化程序设计方法和技术,掌握测绘数据处理问题的基本特点。
5.开发相关程序功能模块,独立完成相关问题概念结构分析、程序结构设计、模块设计、代码编写、调试、测试等工作。
二、实验安排1.实验时数12学时。
2.每实验小组可以由3~4人组成,或独立完成。
若由几个人完成程序设计,应进行合理的分工。
三、实验步骤和要点1.熟悉程序设计任务书的基本内容,调查了解软件需求状况,进行需求分析;2.进行总体设计。
根据所调查收集的资料和任务书的要求,对系统的硬件资源进行初步设计,提出硬件配置计划;进行软件总体设计,设计出软件程序功能的模块;3.根据总体设计的结果,进行详细设计,进行数据存储格式设计、算法等,写出逻辑代码;4.编写程序代码,调试运行;5.程序试运行。
最后同学们可根据自己的选题,写出软件开发设计书一份,打印程序代码和运行结果。
四实验原理高斯正反算:高斯正反算包括两部分内容:高斯正算和高斯反算。
简单的说就是大地地理坐标系坐标(B,L)与其对应的高斯平面直角坐标系坐标(x,y)之间的转换。
若已知大地地理坐标系坐标(B,L)解求对应的高斯平面直角坐标系坐标(x,y)称为高斯正算;反之,则为高斯反算。
空间直角坐标与大地地理坐标转换:地球表面可用一个椭球面表示。
设空间直角坐标系为OXYZ,当椭球的中心与空间直角坐标系原点重合,空间坐标系Z 轴与地球旋转重合(北极方向为正),X 轴正向经度为零时,就可以确定空间直角坐标系与大地地理坐标系的数学关系。
⎪⎩⎪⎨⎧+-=+=+=B H e N Z LB H N Y L B H N X sin ])1([sin cos )(cos cos )(2 式中 N 为卯酉圈曲率半径,B e a N 22sin 1-=; e 为椭球偏心率,222a b a e -=(a ,b 为椭球长半轴和短半轴)。
坐标转换之计算公式+7参+四参模型
坐标转换之计算公式一、参心大地坐标与参心空间直角坐标转换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 为第一辅助系数ab a e 22-= 或 f f e 1*2-= W a N BW e =-=22sin *1(3 参心空间直角坐标转换参心大地坐标[]N BY X H H e N Y X H N Z B XY L -+=+-++==cos ))1(**)()(*arctan()arctan(22222 二 高斯投影及高斯直角坐标系1、高斯投影概述高斯-克吕格投影的条件:1. 是正形投影;2. 中央子午线不变形高斯投影的性质:1. 投影后角度不变;2. 长度比与点位有关,与方向无关;3. 离中央子午线越远变形越大为控制投影后的长度变形,采用分带投影的方法。
常用3度带或6度带分带,城市或工程控制网坐标可采用不按3度带中央子午线的任意带。
2、高斯投影正算公式:522242532236425442232)5814185(cos 120)1(cos 6cos )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、高斯投影反算公式:()()()⎥⎥⎦⎤⎪⎪⎭⎫ ⎝⎛+++⎪⎪⎭⎫ ⎝⎛⎢⎣⎡-++-⎪⎪⎭⎫ ⎝⎛-=⎥⎥⎥⎥⎦⎤⎪⎪⎭⎫ ⎝⎛+++++⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛++-⎪⎪⎭⎫ ⎝⎛=442222224222422224590613601 9351211286242851201 )21(611cos 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、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
大地坐标与直角空间坐标转换计算公式Company Document number:WUUT-WUUY-WBBGB-BWYTT-1982GT大地坐标与直角空间坐标转换计算公式一、参心大地坐标与参心空间直角坐标转换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 参心大地坐标转换为参心空间直角坐标:公式中,N为椭球面卯酉圈的曲率半径,e为椭球的第一偏心率,a、b 椭球的长短半径,f 椭球扁率,W为第一辅助系数a ba e2 2-=或f fe1*2-=西安80椭球参数:长半轴a=6378140±5(m)短半轴b=扁率α=1/3 参心空间直角坐标转换参心大地坐标二高斯投影及高斯直角坐标系1、高斯投影概述高斯-克吕格投影的条件:1. 是正形投影;2. 中央子午线不变形高斯投影的性质:1. 投影后角度不变;2. 长度比与点位有关,与方向无关;3. 离中央子午线越远变形越大为控制投影后的长度变形,采用分带投影的方法。
常用3度带或6度带分带,城市或工程控制网坐标可采用不按3度带中央子午线的任意带。
2、高斯投影正算公式:3、高斯投影反算公式: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转换。
通常有三种转换方法:七参数、简化三参数、Molodensky。
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坐标和所测WGS84坐标,通过一定的数学模型,求解七参数。
若多选几个已知点,通过平差的方法可以获得较好的精度。
如果区域范围不大,最远点间的距离不大于30Km(经验值),这可以用三参数,即只考虑3个平移因子(X平移,Y平移,Z平移),而将旋转因子及比例因子(X旋转,Y旋转,Z旋转,尺度变化K)都视为0,所以三参数只是七参数的一种特例。
北京54和西安80也是两种不同的大地基准面,不同的参考椭球体,他们之间的转换也是同理。
在ArcGIS中提供了三参数、七参数转换法。
而在同一个椭球里的转换都是严密的,在同一个椭球的不同坐标系中转换需要用到四参数转换,举个例子,在深圳既有北京54坐标又有深圳坐标,在这两种坐标之间转换就用到四参数,计算四参数需要两个已知点2 ArcGIS坐标转换例子应注意问题使用ArcGIS如何实现WGS84经纬度坐标到BJ54高斯投影坐标的转换呢在ArcGIS中,这个坐标转换步骤简化了,用户只需要两个步骤就能够直接从最初的WGS84经纬度坐标转换到BJ54高斯投影坐标。
这就是ArcGIS的强大之处。
接下来,我们做一个例子。
假设我们已经知道了7参数,应该如何操作呢在具体的操作前,请大家一定注意以下三点:WGS84的经纬度坐标值是用度来表示,而不能是度分秒表示七参数的平移因子单位是米,旋转因子单位是秒,比例因子单位是百万。
在ArcGIS中,7参数法的名字是Coordinate_Frame 方法。
有人在用ArcGIS进行不同椭球体间的坐标转换时,转换出来的结果不对,然后就写文章说变形如何如何,很可能是由于他们没有注意上面这三个关键的问题造成的。
转换步骤a、定义7参数的地理转换(Create Custom Geographic Transformation)在Arctool中打开Create Custom Geographic Transformation工具,如图1所示:在弹出的窗口中,输入一个转换的名字,如wgs84ToBJ54。
在定义地理转换方法下面,在Method 中选择合适的转换方法如 COORDINATE_FRAME,然后输入平移参数、旋转角度和比例因子,如图2所示:b、投影变换打开工具箱下的Projections and Transformations>Feature>Project,在弹出的窗口中输入要转换的数据以及Output Coordinate System,然后输入第一步自定义的地理坐标系如wgs84ToBJ54,开始投影变换,如图3所示:点击“确定”,完成坐标转换。
3结束语我国现已启用新的坐标系统2000国家大地坐标系,2000国家大地坐标系与现行国家大地坐标系转换、衔接的过渡期仍需一段较长时期,在实际工作、工程中还遇到不同坐标系之间转换,本文针对在生产中从事测绘工作遇到的坐标转换问题提供解决方法和经验,希望对同行有所参考。
【参考文献】[1] 孔祥元、郭际明、刘宗泉.《大地测量学基础》.武汉大学出版社,第一版,2001年9月[2] 李征航、黄劲松.《GPS测量与数据处理》.武汉大学出版社,第一版,2005年3月[3] MAPGIS使用教程大地坐标空间直角坐标转换(2009-10-22 21:12:41)标签:程序计算大地坐标与空间直角坐标转换#include <>?printf("please insert long r a\n");scanf("%f",&a);printf("please insert short r b\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);//double PP2=Z;?B1=atan(Z/sqrt(X*X+Y*Y));?delta=PI/(180*60*60*;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");}else printf("enter error,please enter again"); printf("1kongzhitodadi\n2daditokongzhi\0exit"); scanf("%d",&choice);}return 0;}/高斯投影正、反算如果有需要程序的,可以直接跟我联系,呵呵附:高斯正反算参数。