《坐标转换软件(毕业论文).doc》由会员分享,可在线阅读,更多相关《坐标转换软件(毕业论文).doc(56页珍藏版)》请在三一办公上搜索。
1、摘 要坐标转换已经不是一个新的课题了,随着与人们生活密切相关的测绘事业的迅速发展,全球一体化的形成,越来越多的要求全球测绘资料形成统一规范,尤其是坐标系统的统一。由于各测量单位工作目的不同,所选择的椭球参考系也会有所不同,出现了许多不同形式的坐标系,例如WGS-84坐标系、国家80坐标系、北京54坐标系、独立地方坐标及各种城建坐标。在同一坐标系下坐标的表示方式又有空间直角坐标、大地坐标、平面坐标。根据不同的测绘需求,需要将不同的坐标系下的坐标进行相互转换,在这些坐标转换的过程中既会运用到同一坐标系下的坐标转换模型,又会用到不同参考系下各坐标系间的坐标转换模型。本文研究了各坐标系之间坐标转换的方
2、法,包括同一坐标系下不同坐标(大地坐标,空间直角坐标,平面坐标)之间的相互转换和不同参考系下各坐标系(WGS84坐标系,北京54坐标系,国家80坐标系)之间的相互转换。深入研究了GPS的WGS84坐标系及中国常用的1954年北京坐标系和1980年国家坐标系三种坐标系之间的大地坐标、空间直角坐标及平面坐标的转换模型,并根据不同的转换模型研究了多种转换参数的计算方法。基于这些理论并结合相关的转换模型及算法编程实现了坐标转换系统,用实例对该坐标转换系统进行了测试,分析了实测数据和坐标转换系统一致计算数据之间差异原因,并提出了适当的避免办法。测试结果与预期,符合需求分析的要求。最后总结了该系统的特点和
3、不足,并提出了下一步研究的方向。关键词:地球椭球体,坐标转换模型,布尔沙-沃尔夫(Bursa-Wolf)模型,莫洛金斯基模型AbstractCoordinate conversion is not a new subject.With the rapid development of surveying and the formation of global integration,Increasing demand for the formation of a standardized global mapping data,Especially in the respect of coor
4、dinate system.Since the different purpose of the different units of measurement,The Choice of ellipsoidal reference system will be different too.There were many different forms of coordinate system.For example,WGS84 coordinate system,Xian80 coordinate system,Beijing 54 coordinate system,local coordi
5、nates etc.There are also have three forms to representation coordinates in the same coordinate system,like rectangular space coordinates,geodetic coordinates,plane coordinates.Depending on thedifferent mapping needs,there needs to coordinate conversionThen coordinate system transformation method is
6、researched.Especial transformation method of different coordinate systems,Geocentric Cartesian coordinates of the same coordinate system and it Geodetic System,Geocentric Cartesian coordinates of the different earth models,Geodetic System and map are researched.The simple models and the computing me
7、thod of transformation Parameter with Geoentric Cartesian coordinates and Geodetic system among WGS-84,XiAn80 and BeiJing54 coordinate systems are studied.Also the transformation algorithm betweenIt studied the reliability of transformational data in different coordinate system and same coordinate s
8、ystem.The system test results Consistent with the expected,it also meet the requirements of demand analysis.Finally,this paper sum up the characteristics and the lack of this system,also propose further research directions.Keywords:Earth ellipsoid,Coordinate transformation model,Bursa-Wolf model, Mo
9、lodensky model 目 录摘 要1Abstract21绪论11.1 研究的背景及意义11.2 坐标系间转化的研究现状21.3研究的发展方向32地球椭球与坐标系之基本理论42.1基本概念42.2三种常用坐标系122.3基准162.4常用的坐标系统173 坐标转换的理论与方法243.1坐标系偏移量及转换方法253.2坐标转换模型283.3大地坐标与空间直角坐标之间的转换343.4大地坐标与平面坐标之间的转换364坐标转换系统程序设计及其实现374.1 Visual Basic简介374.2 程序设计思路384.3本章小结395 检测实例与数据分析445.1 参数求解445.2 精度评定4
10、56 结论及展望45参考文献49致谢51附录52 1绪论1.1研究的背景及意义坐标转换已经不是一个新的课题了,随着与人们生活密切相关的测绘事业的迅速发展,全球一体化的形成,越来越多的全球测绘资料形成规范统一,尤其是坐标系统的统一方面。原始的大地测量主要是通过光学仪器进行,这样不可避免的会受到近地面大气和地球曲率的影响,在通视条件下也受到很大的限制,从而对全球测绘资料的一体化产生了一定的约束。另外,由于每一个国家的大地坐标系的建立和发展具有它特殊的历史特性,使得目前仅常用的大地坐标系就超过150个。即使是在同一个国家,在不同的历史时期由于经济的发展变化或习惯的改变也会采用不同的坐标系统。例如:在
11、我国建国之初,为了让尽快搞好基础建设,我国采用了克拉索椭球与我国的实际情况相结合产生了北京54坐标系;而随着测绘事业的发展北京54坐标系的缺陷也随之被表露的越来越明显,因此我国开始研究并使用80坐标系。在实际生活中,由于国家建设的急需在一些地方还来不及布设国家级的大地控制网,而是先建立局部的独立坐标系,然后再将这些独立坐标系转换到国家统一的大地控制网中,而这些坐标系的变换都离不开坐标转换。在国际上,随着第一代卫星导航系统NNSS于1964年在美国研制成功,测绘资料的全球一体化不再是一个梦想。经美国国防部批准,第二代卫星导航系统的开发研究工作于1972年开始,即现在所说的GPS。这套卫星导航系统
12、满足全球范围、全天候、连续实时以及三维导航和定位的要求。正是由于这套卫星导航系统的这些特点,它就很快被广大测绘工作者所接受。但由于坐标系统的不同,对GPS技术的推广使用造成一定的障碍。这样坐标转换的问题再一次被提到了很重要的位置。未来描述卫星运动,处理观测数据和表示测站位置,需要建立与之相应的坐标系统,及协议天球坐标系和协议地球坐标系。其中协议地球坐标系采用的是1984年大地坐标系(Word Geodetic System 1984WGS)该坐标系统的重要参数为:长半轴a=6378137;扁率f=1:298.257223563。而我国采用的坐标系并不是WGS84坐标系,而是BJ54坐标系和西安
13、80坐标系以及地方独立坐标系,这些坐标系是与前苏联的 1942年普尔科沃坐标系相关联的,这些坐标系的主要参数都为:长半轴a=6378245;扁率f=1:298.3。这就使得同一点的坐标值因所在坐标系不同而不同,使得测绘资料在我国的适用范围受到很大的限制,并且对GPS系统在我国的广泛使用造成了一定的约束性,对测绘事业在我国的发展极为不利。为了解决这个问题,测绘工作者做了大量的研究工作,并且已经有许多使用的商品化坐标转换软件投入市场,其中以七参数法用的最为普遍。但是由于此方法中参数较多,对于一部分精度要求不高的GPS用户来说不容易获得足够的已知点进行参数和坐标求解。而且国内现在所使用的坐标转换系统
14、大多是国外软件,出了成本较高以外,也不利于广泛使用,为了方便我国的测绘工作者并解决好这些问题,提出了坐标转换这个课题。1.2 坐标系间转化的研究现状目前,GPS定位技术已经被生活中各个领域广泛使用,GPS观测值是基于以地球的质心为原点的空间直角坐标系,对于采用北京54坐标、西安80坐标或者其他地方独立坐标而言,就需要解决这些不同椭球参考系下的坐标转换问题。经典的七参数法是现在坐标转换中常用的坐标转换方法,其数学模型一般为Mondensky模型、Bursa模型等,这些转换模型已经广泛应用在各种工程测绘中,但在高程系统、坐标精度、公共点的图形数量以及强度上均有缺陷,还需针对不同的转换要求采用不同精
15、度的参数经行坐标转换运算。1.3 研究的发展方向关于坐标系的相关课题研究将随着各方面研究的投入和加大而逐渐深入,在我国国防建设、国民经济等方面实际生产中发挥巨大的作用。(1)大地基准现代化是当前大地测量的根本任务,涉及到大地测量中当前所有学科的基础理论及其应用。我国已经建设成新一代的地心坐标系统和相应的大地基准控制网,并使得大地基准与多个学科交叉使用和发展,大地基准现代化的研究与发展将决定我国测绘学科的综合发展水平,也直接影响我国空间技术的发展水平和国防建设的发展水平。(2)随着大地测量技术的融合,对于基准中的几何量和物理量能够通过不同的模型有机的结合起来,而不是目前分类描述:比如利用全球参考
16、框架建立了高精度的物理高程基础。在将来全球所中地面观测体系的支撑下,将会实现真正意义的全球统一物理几何基准。(3)将我国不同部门、不同时期施测的多个平面和高程分离的大地控制网通过空间大地测量和数据处理技术,科学的整合为全国统一的整体的国家三位控制网,将原来大地测量中所采用分离的几何与物理参数,进行科学的统一的整合;将我国非地心大地坐标框架整体的科学的转换为地心大地坐标框架;提高我国大地坐标框架的地心坐标精度和我国重力基本点的精度;提高还俩昂数据的处理技术,研究现代数据处理技术,提高成果的精度和可靠性。2地球椭球与坐标系之基本理论2.1基本概念讨论与地球相关的坐标转换问题是离不开地球,人们对地球
17、形状和大小的认识,经过了一个相当长的历史过程。地球的自然表面的形状非常复杂的,有高山和深谷,崎岖不平。从整个地球外表来看,由于地球的体积较大,地球表面的高低与整个地球比较起来,显得非常小。同时地球表面上海洋的面积约占地球总面积的3/4,而且比较平坦,所以讨论地球形状时并不是指地球自然表面的形状,而是指有假想的大地水准面所包围的大地球体的形状。所谓大地水准面,它是一个假想的、与平均海面相吻合的水准面。如果将它向大陆延伸,并且使这一延伸面始终保持在任何地方都与该地的铅垂线正交,则这一连续的、无痕迹的、无棱角的闭合水准面,叫做大地水准面,被大地水准面所围成的几何体,是理想的地球形状,叫做大地球体。2
18、.1.1地球椭球无论是参考椭球抑或是总地球椭球,从几何上来说都是旋转椭球,亦即由长半轴为a、短半轴为b的椭圆绕其短轴而成的几何形体,如图2.1所示: 图 2.1 地球椭球椭球体与地球形体非常接近,椭球面是一个现状规则的数学表面。在其上可以做严密的计算,而且推算出的元素同大地水准面上的相应元素非常接近。(1) 总地球椭球与参考椭球总地球椭球:在全球范围内与大地体最为密合的椭球。理论上总地球椭球应该只有一个,但是因为各国地理位置不同、即采用不同的资料,会得到不同的参数,也就会得到不同的总地球椭球和地心坐标系。 参考椭球:为了大地测量的实际需要,各个国家和地区只根据局部的天文、大地和重力测量资料,研
19、究局部大地水准面的情况,确定一个与总椭球接近的椭球,以表示地球的大小,作为处理大地测量成果的依据。这样的椭球只能较好地接近局部地区的大地水准面,不能反应大地体的情况。(2) 椭球的定位和定向椭球的定位定向的实现包括以下几个方面的内容:确定椭球的现状和大小及其物理特性,及四个基本参数;确定椭球的中心位置,即参考椭球的定位;确定以椭球中心为原点的空间直角坐标系坐标轴的指向,及参考椭球的定向;确定大地原点。椭球定位:就是确定椭球中心的位置。分为两类:局部定位和地心定位。局部定位:要求在一定范围内与欧球面与大地水准面最佳密合对椭球中心位置无特殊要求;地心定位:要求在全球范围内椭球面和大地水准面最佳符合
20、,同时要求椭球中心与地心重合或最为接近。无论局部还是地心定位,都不需满足两个平行条件;椭球短半轴平行于地球自转轴;大地起始子午面平行于起始子午面。定位又分为一点定位和多点定位。椭球定向:实质是指椭球旋转轴的指向。2.1.2坐标系坐标系是定义坐标如何实现的一套理论方法,包括定义原点、基本平面和坐标轴的指向,同时还包括基本的数学和物理模型。坐标系根据原点位置的不同可分为参心坐标系、地心坐标系、站心坐标系。这三种坐标系都与地球体固连在一起,与地球同步运动,因而都是地固坐标系。另外,原点在地心的坐标系称为地心地固坐标系。与地固坐标系相对应的是与地球自转无关的天球坐标系或惯性坐标系。坐标系从其表现形式上
21、可以分为空间直角坐标系、空间大地坐标系、站心直角坐标系、极坐标系和曲面坐标系等。从维数上可分为二维坐标系、三维坐标系、多维坐标系等。坐标系分类如图2.2所示: 图2.2 地球坐标系统分类及相互关系2.1.3参考系与参考框架(1)参考系:参考系是在一定观测时间内有特定类型观测量推到该参考系中点位坐标的理论、方法及采用模型和常数的总成。一个参考系包括:一组模型和常数;一套理论和数据处理方法。坐标参考系统,分为天球坐标系与地球坐标系。天球坐标系用于研究天体和人造卫星的定位与运动。地球坐标系用于研究地球上的物体的定位与运动,是以旋转椭球为参照体建立起来的坐标系统,分为大地坐标系(图2.3)和空间直角坐
22、标系(图2.4)两种形式。 图2.3 大地坐标系 图2.4 空间直角坐标系(2)参考框架:参考框架是坐标参考系的实现。参考系定义明确且紧密,抽象难把握。需要通过一些具体直观的点来描述或反映某一特定的坐标参考系,这些满足特定坐标参考系的点就是人们通常所说的坐标参考框架。一般而言,只要涉及与空间位置有关的问题,就会涉及参考系;而涉及参考系就会涉及参考框架。考虑坐标变化时,还需要一个时间历元,故时间尺度也是坐标参考框架的一部分。参考框架是通过大地测量手段确定的固定在地面上的控制网所构建的,分为坐标参考框架、高程参考框架、重力参考框架。2.1.4参心坐标系1980年国家地理坐标系、1954年北京坐标系
23、均是属于参心坐标系。建立参心坐标系,包括如下几个方面内容了:(l)确定某一地球椭球体形状和大小。(2)确定该地球椭球体中心位置,简称定位。(3)确定该地球椭球体中心为原点的空间直角坐标系坐标轴的指向,简称定向。(4)确定大地原点。各个国家或地区,为了处理大地测量成果,计算点位坐标,测绘地图和进行工程建设,需要建立一个适合本国的地理坐标系。早期建立的大地坐标系都是利用天文观测、天文大地水准面测量和重力大地水准面高度差测量的方法,设定地面坐标的原点(即大地原点),建立天文大地坐标网,然后通过相对地面坐标原点及天文大地坐标网点进行弧度测量而建立的局部坐标系,它的地球椭球体的定位和定向是依据地面参考点
24、-大地原点来实现的,即相对定位,它使得在一定范围内地球椭球体表面与大地水准面有最佳的符合。当前,世界上100多个国家和地区己建立200多个参心坐标系。这些坐标系选用的参考椭球体参数不同,椭球体的实际定位和定向不同。坐标系原点不同。它们只能满足各自的需要。不是全球统一的坐标系,故参心坐标系又称为局部坐标系。参心坐标系与地心坐标系的关系如图2.5所示: 图2.5: 参心坐标系与地心坐标系的关系2.1.5地心坐标系 以地球质心为坐标系原点的坐标系,称为地心坐标系。地心坐标系包括地心大地坐标系和地心直角坐标系。建立一个地心坐标系的三个条件:(l)确定地球椭球体。这个地球椭球体的大小a和形状f要同大地球
25、体最佳吻合。(2)地心的定位和定向。坐标系原点建立于地球质心,首子午面与国际时间局(BIH)平均零子午面重合,Z轴同国际协议地极(CIP)的极轴相重合。(3)尺度。采用标准的国际米作为测量长度的尺度。一、地心大地坐标系大地坐标系又称地理坐标系。地理坐标系指的是:由赤道和格林经线为基准圈的球面坐标系。地球椭球体表面上任意一点的地理坐标,可以用地理纬度、地理经度和大地高H来表示。如图2.6所示: 图2.6:地心坐标系地面上某点的地理经度是这样确定的:它是格林子午线和该点子午线之间所截的赤道短弧或此短弧所对的球心角或极角;它以格林子午线为基准,向东或向西,由0到180计量,算至该点子午线作为该点的地
26、理经度。地理经度简称经度。有格林子午线向东计量的叫东经,缀以“E”,向西计量的叫西经,缀以“W”。例如:我国北京的地理经度L=11628E。地面上某点的地理纬度是这样确定的:它是该点的椭圆子午线的法线与赤道面的交角,作为该点的地理纬度,它以赤道面为基准,向北(极)或向南(极),由0到90计量。地理纬度简称为纬度。由赤道面向北计量的叫北纬,缀以“N”;向南计量的叫南纬,缀以“S。例如:我国北京的地理纬度B=3954N。大地高H是大地点沿椭球法线方向至椭球体面的距离。从椭球体面起算,向外为正,向内为负。由于地球模型不同,世界上有过许多种地心坐标系。如:WGS60、WGS66、WGS72、WGS84
27、、PE90等。二、地心空间直角坐标系地心空间直角坐标系定义:原点O与地球质心重合,z轴指向地球北极,X轴指向格林尼治子午面与地球赤道的交点E,Y轴垂直于XOZ平面并与XZ轴构成右手坐标系,任意一点的位置可用(X、Y、Z)坐标系来表示,它们是与地心大地坐标系相对应的。2000中国大地坐标系是中国为迎接21世纪经济的持续发展,为信息化发展提供一个基础地理平台,更科学的描述动态的地球,特别是随着全球定位系统等空间大地测量技术的不断发展和完善,而推出的地心三维坐标系。在下面将具体介绍2000中国大地坐标系。2.1.6国际地球参考系与国际地球参考框架(1)国际地球参考系:国际地球参考系是由国际地球自转服
28、务局定义的一种协议地球坐标参考系统。它根据一定要求,在地壳表面建立地面观测台站进行空间大地测量,并根据协议地球残系的定义,采用一组国际推荐的模型和常数系统,对观测数据进行处理,解算出各观测台站在某一历元的台站坐标和速度场。它的定义:原点位于地球质心,地球质心是包括海洋和大气的整个地球的质量中心;长度单位为m(SI),是在广义相对论框架下的定义;坐标轴的初始定向与国际时间局BIH1984.0历元的定义一致;定向的时间演化相对于地壳不产生残余的全球性旋转,即要满足无净旋转条件。ITRS是目前国际上最精确、最稳定的全球性地心坐标系。(2)国际地球参考框架: 国际地球参考系的实现方式是国际地球参考框架
29、,即它是通过一组固定于地球表面而且只作线性运动的大地点的坐标及坐标变化速率组成,这种站点装备有不同的空间大地测量系统并由IERS中心局的地球参考框架部负责建立和维护。目前的ITRF已有ITRF88、ITRF89、ITRF90、ITRF91、ITRF92、ITRF93、ITRF94、ITRF96、ITRF97、ITRF2000、ITRF2005。常用的有 ITRF94、ITRF97、ITRF2000、ITRF2005。2.2三种常用坐标2.2.1空间坐标过空间点O作三条相互垂直的具有相同单位长度的轴,它们都以O作为原点。这三条数轴分别称为X轴、Y轴、Z轴,这三条轴统称为坐标轴。各轴之间的顺序要符
30、合右手法则,即用手握住Z轴,让右手的四指从X轴的正向以90度转向Y轴的正向,这时右手大拇指所指的方向就是Z轴的方向。这样的三个坐标轴构成了空间直角坐标系。三条坐标轴中的任意两条都可以确定一个平面,称为坐标面。它们是:由X轴及Z轴确定的OXZ平面,由X轴及Y轴确定的OXY平面,由Y轴及Z轴确定的OYZ平面。空间被这三个相互垂直的坐标面分成了八个部分,每个部分称为一个卦限。每一卦限位于Z,X,Y轴的正半轴,从第一挂限开始,在OXY平面上方的卦限,按顺时针方向依次称为第一、四、三、二卦限,下方的卦限依次为第五、八、七、六卦限。设点M为空间的某一定点,过点M分别作为三个平面分别垂直于x、y、z轴,依次
31、交x、y、z轴于P、Q、R,设点P、Q、R在x、y、z轴上的坐标读数分别为x、y、z,那么就得到与点M唯一对应确定的有序实数组(x,y,z),该有序实数组(x,y,z)称为点M的坐标,记作M(x,y,z),这样M点的空间坐标就确定了,其中x、y、z分别为点M的横坐标值、纵坐标值、竖坐标值。而在测量学中,空间直角坐标系的定义是这样的:以椭球中心作为空间直角坐标系的坐标原点,以赤道WAE和起始子午面NGS的交线为X轴,以椭球向北为正的短轴为Z轴,而以赤道面上与X轴正交的方向为Y轴,构成以右手定则为准则的空间坐标系O-XYZ。地面点P曲的点位用OP曲在这三个坐标轴上的投影X、Y、Z表示,(x,y,z
32、)也就是OP曲的空间直角坐标。它们是与大地坐标系一一对应的。空间直角坐标系的坐标原点与参考椭球的中心重合,Z轴正向指向参考椭球的北极,X轴正向指向起始子午面的与赤道的交点,Y轴按右手系与X轴呈90度夹角且位于赤道面上。某点在空间中的坐标可用该点在空间坐标系的各个坐标轴上的投影来表示,如下图2.7所示: 图2.7: 空间直角坐标系2.2.2平面坐标平面直角坐标系是利用投影,将空间坐标通过某种数学变换映射到平面上,这种变换成为投影变换。投影变换的方法有很多,如Lambuda投影,UTM投影等,在我国一般采用的是高斯克吕格投影,也称为高斯投影。地形测图以及许多的测量定位应用在现实中是我们常见的平面直
33、角坐标。对于一个国家或较大区域,应按照一定的数学法则将参考椭球面上的各点的大地经纬度投影为平面上相对应点的平面直角坐标。由于地球椭球面是不可展曲面,所以无论采用什么样的投影都会产生一定的变形。投影变形一般分为长度变形、角度变形和面积变形这三种。根据控制测量的任务和目的应当采用等角投影。在采用正形投影时,还要求长度和面积变形不大,并且能用简单的公式来计算这些变形而带来的改正数。为了解决这些矛盾,测量上往往是将一个大的区域按照一定规律分成若干个小的区域。每个区域单独投影,并组成自身的直角坐标系。然后,再将这些带用简单的数学方法联系起来,从而组成统一系统。目前测量上广泛采用的是高斯投影,它是一种正形
34、投影,它的特点是:没有角度变形,在不同点的长度比随点位而异,但在同一点上个方向的长度比相同。高斯克吕格正形投影又称横轴椭圆柱投影,即椭圆柱内面横套在地球椭球的外表面,椭圆柱的中心通过椭球的中心,并在某一中央子午线上相切,该中央子午线就是高斯平面直角坐标系的X轴,X轴没有长度变形,赤道在椭圆柱上的投影是高斯平面直角坐标系的Y轴,把椭球柱展开,就得到以(X,Y)为坐标的高斯平面直角坐标系。 2.2.3大地坐标大地坐标系以大地经度L、大地纬度B、大地高H来表示空间某一点的位置。如图2.6,O表示椭球中心,NGS为起始大地子午面,WAE为赤道面,地面点P曲的法线P曲Kp交椭球面于P点,NGS为P地的子
35、午面。则地面点P曲的大地坐标定义为:大地纬度是B-P曲Kp,与赤道面WAE的夹角,从赤道面开始起算,以正为北(0度90度),称为北纬,以南为负(0度-90度),称为南纬。大地经度是L-P地的子午面NGS与起始子午面NGS所构成的二面角,以东为正(0度180度),称为东经,以西为负(0度-180度),称为西经;大地高H-P地沿法线方向到椭球体的P地P,从椭球面为起算标准,以外为正,以内为负。大地高H与水准测量中的正常高HN或正高HN有以下关系H=HN+N=HV+式中是高程异常,N是大地水准面差距。若将地面点P地沿垂线投影到水准面上,设投影点为P,再将P沿法线投影到椭球面上得到Po如图2-3,则沿
36、垂线的距离P地P就是HN,而PPo就是大地水准差距N。2.3 基准所为基准就是指为了描述空间中的某个位置而定义的面、点和线,在大地测量中,基准指描述地球形状的参考椭球的参数,如参考椭球的长半轴,短半轴,以及参考椭球空间中的定位和定向,还有在描述这些位置时单位长度的定义。不同的坐标系统会使用的基准也不同。在大地测量中,根据参考椭球所选原点的位置不同,可以非为地心坐标系和参心坐标系。基准面是利用某一特定椭球体对特定地区地球表面的无线逼近,因此每个国家或地区都有各自不同的基准面,如我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国常用的两个大地基准面。我国参照前苏联从1953年起采用克拉索
37、夫斯基椭球体建立了我国的北京54坐标系,1978年采用国际大地测量协会推荐的1975地球椭球体建立了我国的新的大地坐标系西安80坐标系。目前大地测量基本上仍以北京54坐标系作为参考,西安80坐标和北京54之间的转换可查阅国家测绘局公布的对照表。WGS1984基准面采用WGS84椭球体,它是地心坐标系,即以地球的质心作为椭球体原点,目前GPS的测量数据大多以WGS1984为基准。2.4 常用的坐标系统2.4.1 WGS84坐标系统GPS的测量成果,于1987年1月10日开始采用WGS-84世界大地坐标系以取代1987年以前所采用的WGS-72。WGS-84坐标系是由美国国防部制图局依据卫星定位测
38、量成果而建立的一种协议地球坐标系 。它是GPS卫星广播星历和精密星历的参考系。WGS-84坐标系理论上是一个以地球质心为坐标原点的地心坐标系,其坐标系的定向与BIH1984.0所定义的方向一致,其Z轴指向此BIH系统所定义的协议地极的方向,X轴指向国际时间句(BIH)1984.0的零度子午面与赤道的交点。定义如下:原点是指地球的质量中心;Z轴是指向国际时间局(BIH)1984.0定义的协议地球极方向;X轴是指向BIH1984.0的子午面和赤道面的交点;Y轴是一个在符合右手准则从地心与地球赤道的交点。WGS84椭球体的基本参数为:长半径 a= 6378137.000m+2m短半径 b= 6356
39、755.2881575287m扁率 = 1:298.257第一偏心率平方= 0.006 694 384 999 588椭球偏扁率 f=1/298.257223563地球自转角速度 =7.292115 * 10-5 rad.s-1 地球引力常数GM=398600.5km3.s-2WGS84坐标采用WGS84椭球,其常数为国际大地测量与地球物理学联合会(IUGG)第17届大会的推荐值。2.4.2 2000大地坐标系统2000中国大地坐标系是中国为迎接21世纪经济的持续发展,为信息化发展提供一个基础地理平台,更科学的描述动态的地球,特别是随着全球定位系统等空间大地测量技术的不断发展和完善,而推出的地
40、心三维坐标系。世界各国都在不断更新和完善各自的大地坐标系统及其框架。地心坐标系统及其框架正在逐渐取代传统的非地心大地坐标系统及其框架。空间大地测量技术的发展,对大地基准的现代化提出了新的要求,也提供了实现的可能。中国启用地心3维坐标系统为国家大地坐标系统,对我国国土范围内所考虑对象的空间位置(不论该对象是处于静态还是动态)都需要一个全国统一的、协调一致的大地坐标系统和大地坐标框架。3维坐标系统的原点可以是参心的,也可以是地心的。若采用以地球质心为大地坐标系的原点,即采用地心坐标系的优势是明显的。因为这种坐标系统是阐明地球上各种地理和物理现象,特别是空间物体运动的本始参照系。但长期以来由于人类不
41、能精确确定地心的位置,因而较少使用。但目前利用空间技术等手段已可在厘米级确定它的位置。因此采用地心坐标系在当今既有它的必要性,也有可能性。鉴于中国经济、社会和科学技术的发展需求和可能,中国政府决定采用3维大地坐标系统,既从2008年7月1日起正式启用中国大地坐标系(CGCS2000)作为国家法定的坐标系,作为我国新一代的平面基准。我国采用的3维地心坐标系统CGCS2000 的定义和迄今为止比较符合客观实际的国际通用地面参考系统(ITRS)的定义在原则上保持一致。一个大地坐标系由原点、坐标轴指向和尺度定义。CGCS2000的定义满足下列条件:(1)它是地心坐标系,其原点在包括海洋和大气的整个地球
42、的质量中心;(2)它的长度单位为米(SI),这一尺度与地心局部框架的地心坐标时一致;(3)它的初始定向由1984.0时国际时间局(BIH)定向给定;(4)它的定向的时间的演化由整个地球的水平构造运动无净旋转的条件保证。(5)与此定义相应,存在一个直角坐标系XYZ,其原点在地球质心,Z轴指向参考极方向;X轴为参考子午面与过原点且同Z轴正交的赤道面的交线;Y轴与Z、X轴构成右手直角坐标系。如图2.8所示: 图2.8 CGCS2000CGCS2000参考椭球采用如下4个定义常数:(1)地球赤道半径:a = 6378 137m(2)地球(包括大气)的地心引力常数:GM = 3.986 004 418
43、(3)地球的动力形状因子:= 1.082 629 832 258 (4)地球旋转速度: = 7.292 115 CGCS2000参考椭球为一旋转椭球,其几何中心与坐标系的原点重合,旋转轴与坐标系的Z轴重合。参考椭球面是大地坐标(经纬度和高程)的参考面 。另一方面,CGCS2000参考椭球为一正常椭球,即椭球面为位函数U=const.的等位面。正常椭球面是地球正常重力场(正常重力矢量)的参考面。所以,CGCS2000参考椭球又可以称为CGCS2000正常椭球。应当指出,定义常数 与常数a,GM,f,在理论上是等价的,其中f为椭球的几何扁率;而且在这里,当f采用1:298.257 101时,这两组
44、常数至给定的有效位是一致的。2000中国大地坐标框架建设主要策略:第一层次为CGCS2000连续运行GPS网(即连续运行CORS基站),我国维持CGCS2000主要依靠连续运行GPS观测站,它们是CGCS2000的骨架,其坐标精度为毫米级,速度精度为1mm/a.第二层次为“2000国家GPS大地控制网”包括中国全部领土和领海的高精度GPS网点,即全国GPS一、二级网、国家GPS A、B级网、地壳运动监测网和地壳运动观测网络工程网,共约2500多点,其三维地心坐标精度约为3cm。第三层次为全国天文大地控制网(约有50千点),它是CGCS2000的加密框架。它由全国天文大地网与2000GPS控制网
45、联合平差后的网点坐标体现,三维点位误差约为0.3m,大地高误差不超过0.5m。2.4.3 西安80坐标系1980年西安坐标系是经全国天文大地网整体平差而建立的另一个参心坐标系,采用的参考椭球比较合适,坐标轴的指向明确,参考椭球面与我国的大地水准面吻合的较好,椭球定位比较符合中国的实际,大地原点位于陕西泾阳县永乐镇。在西安市以北60KM,可简称西安原点。大地经纬度的概略值是:Lo=108度55分,B=34度32分。大地点高程是以1956年青岛验潮站求出的黄海平均海水面为基准。其基本参数为:长半轴 a = 6378140m短半轴 b = 6356755.2881575287 m扁率 = 1:298
46、.257第一偏心率平方= 0.006694384999588地心引力常数GM=3.986005*1014m3s-21980年西安坐标系建立以后,进行了全国天文大地网整体平差。通过平差消除了分区局部平差和逐级控制造成的影响,提高了平差结果的精度,平差后获得约50千点的坐标,构成了1980 年西安坐标系的基本参考框架。但全国天文大地网整体平差时,采用了大地原点固定的、以天文方位角与起始边作控制的单点自由网平差法,天文方位角与起始边的系统误差无条件的带入整网平差结果中,整体平差后的全国天文大地网存在较大的系统性扭差 1980年西安坐标系仍然存在以下问题:一是只能提供2维坐标,不能提供高精度3维坐标;
47、二是采用了国际大地测量协会(IAG)1975年推荐的椭球,该椭球与IERS推荐的椭球相比,长半轴大了3m这可能引起约量级的的长度误差;三是椭球短轴指向JYD1968.0极原点,与国际上通用的椭球短轴指向不一致;四是椭球定位没有顾及到占中国全部国土面积近三分之一的海域范围。尽管1980年西安坐标系比1954年北京坐标系有所改善,但并没有发生实质性的变化。2.4.4 北京54坐标系1954年北京坐标系是通过我国东北呼玛、吉拉林、东宁3个基线网与前苏联远东大地控制网相连接,将前苏联1942年普尔科沃坐标系延伸至我国的一个坐标系。于是1954年北京坐标系实际上是前苏联1942年普尔科沃坐标系在我国的扩展和延伸。其基本参数为:长半轴 a = 6378245m短半轴 b = 6356863.0188m扁率 =1:298.3第一偏心率平方= 0.00669342162297第二偏心率平方e2= 0.00673852541468 1954年北京坐标系建立以后,我国天文大地网采取边布设边平差的方式,获取了40千点的坐标,从而构成了1954年北京坐标系的基本参考框架。1954年北京坐标系存在以下问题:(1)采用的