1、常用坐标系与高程系简介 2009-09-27 10:06:45| 分类： GIS 技术 | 标签： |字号大中小 订阅 坐标系的概念1.坐标系的定义： 如果空间上任意一点 P 的位置，可以用一组基于某一时间系统时刻 t 的空间结构的数学描述来确定，则这个空间结构可以称为坐标系，数学描述称为 P 点在该坐标系中的坐标。牛顿运动学原理要求坐标系是惯性的，惯性是每个物体所固有的当没有外力作用时保持静止或匀速直线运动的属性，基于这个特性，惯性坐标系的定义需与时间无关，通常这样的坐标系需要三个属性来描述(这应该是三维空间的本性吧)，首先一个是原点(O)，就是坐标系的中心点，第二个是过原点的任意直线(这里
2、称为 Z 轴)，第三个是过原点且与 Z 轴不重合的任意直线(这里称为 X 轴) ，如果 X 轴与 Z 轴垂直，会带来较优美的数学描述，我们称这样的坐标系是笛卡尔坐标系。P 点的位置可以用 P 到原点的距离 r，OP 与 Z 轴的夹角，OP 与 X 轴的夹角来描述 (当然也可以有其它等价描述) ，可以证明这个描述确定的 P 点是唯一的。 2.GPS 领域常用坐标系模型： 在 GPS 测量中，最常用的坐标系模型是协议地球坐标系，该坐标系随同地球一起旋转，讨论随地球一起自转的目标位置，用这类坐标系方便；另外一类是协议天球坐标系，这个坐标系随同太阳系一同旋转，与地球自转无关，讨论卫星轨道运动时，用这类
3、坐标系方便。 天球坐标系的定义是这样的，原点是地球质心(O)， Z 轴指向地球自转轴( 天极，向北为正)， X 轴指向春分点，根据春分点的定义可以证明 X 轴与 Z 轴互相垂直，且 X 轴在赤道面上，同时为数学描述方便，引入与 XOZ 成右手旋转关系的 Y 轴。因为地球自转轴受其它天体影响(日、月)在空间产生进动，使得春分点变化 (章动和岁差)，导致用“瞬时天极”定义的坐标系不断旋转，而旋转的坐标系表现出非惯性的特性，不能直接应用牛顿定律。我们可以用某一历元时刻的天极和春分点(协议天极和协议春分点)定义一个三轴指向不变的天球坐标系，称为固定极天球坐标系。 地球坐标系的定义是这样的，原点为地球质
4、心(O) ，Z 轴为地球自转轴，X 轴指向地球上赤道的某一固定“刚性”点，所谓“刚性”是指其自转速度与地球一致，同时也为数学描述方便，引入与 XOZ 成右手旋转关系的 Y 轴。地球不是一个严格刚性的球体，Z 轴在地球上随时间而变，称为极移，同天球坐标系一样，需要指定一个固定极为 Z 轴，这样的地球坐标系称为固定极地球坐标系。可以证明当观察地球上的物体时，该坐标系是惯性的。如果一个坐标系 OXYZ，O 不是地球质心，Z 轴与地球自转轴平行，则这个坐标系具有与地球相同的自转角速度，我们也把此类坐标系称为地球坐标系。 3.协议坐标系统： 那么，什么是“协议”坐标系呢？通常，理论上坐标系由定义的坐标原
5、点和坐标轴指向来确定。坐标系一经定义，任意几何点都具有唯一一组在该坐标系内的坐标值，反之，一组该坐标系内的坐标值就唯一定义了一个几何点。实际应用中，在已知若干参考点的坐标值后，通过观测又可反过来定义该坐标系。可以将前一种方式称为坐标系的理论定义。而由一系列已知点所定义的坐标系称为协议坐标系，这些已知参考点构成所谓的坐标框架。在点位坐标值不存在误差的情况下，这两种方式对坐标系的定义是一致的。事实上点位的坐标值通常是通过一定的测量手段得到，它们总是有误差的，由它们定义的协议坐标系与原来的理论定义的坐标系会有所不同，凡依据这些点测定的其它点位坐标值均属于这一协议坐标系而不属于理论定义的坐标系。由坐标
6、框架定义的固定极天球坐标系和固定极地球坐标系，称为协议天球坐标系和协议地球坐标系。 一个完整的坐标系统，除了定义坐标系外，还需要定义基准，所谓基准就是在指定坐标系中的尺度单位、基本的点、线、面(如椭球面、水准面等 )，本专题讨论点 P 的坐标在不同坐标系统的转换，主要是在 WGS-84 坐标系统和中国国家地方坐标系统的转换，下一章先讨论 WGS-84 坐标系统和中国国家地方坐标系统的定义。 WGS-84 大地坐标系统的几何定义是：原点位于地球质心，Z 轴指向 BIH1984.0 定义的协议地球极(CTP)方向，X 轴指向 BIH1984.0 的零子午面和 CTP 赤道的交点。对应与WGS-84
7、 大地坐标系有一 WGS-84 椭球，WGS-84 椭球及有关常数采用国际大地测量 (IAG)和地球物理联合会(IUGG)第 17 届大会大地测量常数的推荐值，四个基本常数为：长半轴a、地心引力常数 GM、地球重力场模型系数 C2.0、地球自转角速度 ，其它的椭球常数可以根据以上四个常数计算得到，如偏心率、扁率等。常用坐标系1、北京 54 坐标系北京 54 坐标系为参心大地坐标系，大地上的一点可用经度 L54、纬度 M54 和大地高 H54 定位，它是以克拉索夫斯基椭球为基础，经局部平差后产生的坐标系。1954 年北京坐标系的历史：新中国成立以后，我国大地测量进入了全面发展时期，再全国范围内开
8、展了正规的，全面的大地测量和测图工作，迫切需要建立一个参心大地坐标系。由于当时的“一边倒”政治趋向，故我国采用了前苏联的克拉索夫斯基椭球参数，并与前苏联 1942 年坐标系进行联测，通过计算建立了我国大地坐标系，定名为 1954 年北京坐标系。因此，1954 年北京坐标系可以认为是前苏联 1942 年坐标系的延伸。它的原点不在北京而是在前苏联的普尔科沃。北京 54 坐标系，属三心坐标系，长轴 6378245m，短轴 6356863，扁率 1/298.3；2、西安 80 坐标系1978 年 4 月在西安召开全国天文大地网平差会议，确定重新定位，建立我国新的坐标系。为此有了1980 年国家大地坐标
9、系。1980 年国家大地坐标系采用地球椭球基本参数为 1975 年国际大地测量与地球物理联合会第十六届大会推荐的数据，即 IAG75 地球椭球体。该坐标系的大地原点设在我国中部的陕西省泾阳县永乐镇，位于西安市西北方向约 60 公里，故称 1980 年西安坐标系，又简称西安大地原点。基准面采用青岛大港验潮站 19521979 年确定的黄海平均海水面（即 1985 国家高程基准）。西安 80 坐标系，属三心坐标系，长轴 6378140m，短轴 6356755，扁率 1/298.257221013、2000 国家大地坐标系的定义国家大地坐标系的定义包括坐标系的原点、三个坐标轴的指向、尺度以及地球椭球
10、的 4 个基本参数的定义。 2000 国家大地坐标系的原点为包括海洋和大气的整个地球的质量中心；2000 国家大地坐标系的 Z轴由原点指向历元 2000.0 的地球参考极的方向，该历元的指向由国际时间局给定的历元为 1984.0 的初始指向推算，定向的时间演化保证相对于地壳不产生残余的全球旋转，X 轴由原点指向格林尼治参考子午线与地球赤道面（历元 2000.0）的交点，Y 轴与 Z 轴、X 轴构成右手正交坐标系。采用广义相对论意义下的尺度。2000 国家大地坐标系，长半轴 6378137m，扁率 f=1/298.257222101，地心引力常数GM3.9860044181014m3s-2，自转
11、角速度 7.292l1510-5rads-1。4、WGS84 坐标系WGS84 坐标系（WorldGeodeticSystem ）是一种国际上采用的地心坐标系。坐标原点为地球质心，其地心空间直角坐标系的 Z 轴指向国际时间局（BIH）1984.0 定义的协议地极（CTP）方向，X 轴指向BIH1984.0 的协议子午面和 CTP 赤道的交点，Y 轴与 Z 轴、X 轴垂直构成右手坐标系，称为 1984 年世界大地坐标系。这是一个国际协议地球参考系统（ITRS ），是目前国际上统一采用的大地坐标系。GPS 广播星历是以 WGS-84 坐标系为根据的。WGS84 坐标系，长轴 6378137.000
12、m，短轴 6356752.314，扁率 1/298.257223563。 由于采用的椭球基准不一样，并且由于投影的局限性，使的全国各地并不存在一至的转换参数。对于这种转换由于量较大，有条件的话，一般都采用 GPS 联测已知点，应用 GPS 软件自动完成坐标的转换。当然若条件不许可，且有足够的重合点，也可以进行人工解算。注释：WGS84 经纬度坐标与西安 80 的坐标转换由 gisnewcomer2002-08-0810:05因为 7 参数 3 参数无法获得且没有已知点，想直接把 WGS84 的经纬度当作西安 80 的经纬度，经高斯投影得到西安 80 平面坐标。这样的简化转换在精度要求不高时是否
13、实用？浙江的 WGS84 经纬度和西安 80 经纬度到底相差多大？另外，有人能提供浙江 40 带的 7 参数/3 参数或者近似参数吗？-由 LLM2002-08-0813:21这要看你的精度要求到底有多高,在 1:400 万上应该没问题.LLM-由 gisnewcomer2002-08-0814:44地图是 1:10000 的，精度要求不高，但总不能误差 50 米以上吧。gisnewcomer-由 LLM2002-08-0816:01能差几百米 LLM-由 gisnewcomer2002-08-0817:25没这么严重吧？这岂不是束手无策了？gisnewcomer-由 chiaojg2002-
14、08-0818:58必须有已知点,然后求参数,否则误差很大 chiaojg-由三木2002-08-1018:43我想问一下如果在 mapgis 中直接用投影转换行不行？三木-由 LLM2002-08-1211:24在 mapgis 中也需知道参数.LLM-由 gisnewcomer2002-08-1213:15一个已知点行不行？WGS84 经纬度坐标没有高程能不能转换为西安 80 平面坐标？gisnewcomer-由 LLM2002-08-1213:45一个已知点只能平移,小范围可以.LLM-由三木2002-08-1215:27在 1：10000 的图上那最少需要几个已知点才可以进行坐标转换三
15、木-由 gisnewcomer2002-08-1310:55WGS84 经纬度坐标没有高程能不能转换为西安 80 平面坐标？几个已知点才好求出 7 参数？gisnewcomer-由 LLM2002-08-147:28在 1：10000 的图上那最少要一个已知点才可以平移转换,不同图平量可能不同.WGS84 经纬度坐标没有高程能转换为西安 80 平面坐标,这只是忽略了高程变化对转换中平面位置的影响.求七个参数不是一种好方法.LLM-由 gisnewcomer2002-08-1410:53但在转换 BLHXYZ 的公式中 H 是必需的呀，没有 H 好像不行的。如果转换时假设 H=0 好像说不过去。
16、难道任何 WGS 经纬度坐标相同而高程不同的点转换到西安 80 平面坐标会相同么？为什么求 7 参数不好？这不是标准转换公式么？LLM，我的 QQ 是 15971629，想向你请教一下这些问题。gisnewcomer-由 LLM2002-08-1422:26可用平均高程代替。只用能掌握的有限点求七个参数不可能准确。用这些点直接变换会更简单。LLM常用高程系1、我国常用高程系“1956 年黄海高程系”，是在 1956 年确定的。它是根据青岛验潮站 1950 年到 1956 年的黄海验潮资料，求出该站验潮井里横按铜丝的高度为 3.61 米，所以就确定这个钢丝以下 3.61 米处为黄海平均海水面。从
17、这个平均海水面起，于 1956 年推算出青岛水准原点的高程为 72.289 米。国家 85 高程基准其实也是黄海高程基准，只不过老的叫“1956 年黄海高程系统”，新的叫“1985 国家高程基准”，新的比旧的低 0.029m我国于 1956 年规定以黄海(青岛) 的多年平均海平面作为统一基面，为中国第一个国家高程系统，从而结束了过去高程系统繁杂的局面。但由于计算这个基面所依据的青岛验潮站的资料系列（1950年1956 年）较短等原因，中国测绘主管部门决定重新计算黄海平均海面，以青岛验潮站 1952 年1979年的潮汐观测资料为计算依据，并用精密水准测量接测位于青岛的中华人民共和国水准原点，得出
18、 1985年国家高程基准高程和 1956 年黄海高程的关系为：1985 年国家高程基准高程=1956 年黄海高程-0.029m。1985 年国家高程基准已于 1987 年 5 月开始启用，1956 年黄海高程系同时废止。各高程系统之间的关系：56 黄海高程基准：+0.00085 高程基准（最新的黄海高程）：56 高程基准-0.029吴淞高程系统：56 高程基准+1.688珠江高程系统：56 高程基准-0.586我国目前通用的高程基准是：85 高程基准2、EGM96EGM96 (Earth Gravitational Model 1996) is a geopotential model of
19、the Earth consisting of spherical harmonic coefficients complete to degree and order 360. It is a composite solution, consisting of: (1) a combination solution to degree and order 70, (2) a block diagonal solution from degree 171 to 621, and (3) the quadrature solution at degree 360. This model is t
20、he result of a collaboration between the National Imagery and Mapping Agency (NIMA), the NASA Goddard Space Flight Center (GSFC), and Ohio State University.The joint project took advantage of new surface gravity data from many different regions of the globe, including data newly released from the NI
21、MA archives. Major terrestrial gravity acquisitions by NIMA since 1992 include airborne gravity surveys over Greenland and parts of the Arctic and the Antarctic, surveyed by the Naval Research Lab (NRL) and cooperative gravity collection projects, several which were undertaken with the University of
22、 Leeds. These collection efforts have improved the data holdings over many of the worlds land areas, including Africa, Canada, parts of South America and Africa, Southeast Asia, Eastern Europe, and the former Soviet Union. In addition, there have been major efforts to improve NIMAs existing 30 mean
23、anomaly database through contributions over various countries in Asia.NIMA also computed and made available 3030 mean altimeter derived gravity anomalies from the GEOSAT Geodetic Mission. EGM96 also included altimeter derived anomalies derived from ERS-1 by Kort & Matrikelstyrelsen (KMS), (National
24、Survey and Cadastre, Denmark) over portions of the Arctic, and the Antarctic, as well as the altimeter derived anomalies of Schoene 1996 over the Weddell Sea.PGM2000A is an EGM96 derivative model that incorporates normal equations for the dynamic ocean topography implied by the POCM4B ocean circulat
25、ion model.注释：1985 国家高程基准与全球似大地水准面之间的系统差及其分布规律一文中指出：利用分布全国大陆范围的 GPS 网 949 个点的 GPS/水准数据和地球重力场模型 EGM96、DQM99A，求出 1985 国家高程基准点与 WGS84 定义的似大地水准面 之间有 35.7cm 的垂直偏差，1985 国家高程基准面的系统差自东向西、自南向北明显增大。坐标系统间的转换在实际应用中需要将 GPS 观测成果点位的 WGS-84 坐标转换为地面网的坐标，首先要把点位的WGS-84 坐标转换成国家（或地区）的大地坐标，然后再把大地坐标转换成高斯平面直角坐标。1)WGS-84 坐标转
26、换成国家大地坐标（图 2-6）：这是把地心坐标转换为参心坐标，其实质就是确定转换参数。通常由三个平移参数，三个旋转参数和一个尺度变化因素 组成 7 个转换参数，其直角坐标转换公式为（23）式：式中：X0,Y0,Z0 为平移参数k 为尺度变化因子X,Y,Z 为旋转参数,其大地坐标转换公式简略写成（2-4）式：在 GPS 卫星定位网的测量中，一般只需要进行两个坐标系之间作基线向量的转换，其转换公式为（2-5）式：上式中不需要平移参数，只要三个旋转参数和一个尺度变化因子即可2)大地坐标（ B，L）转换为高斯平面直角坐标。这种转换按高斯投影正算公式（2-6）式进行：式中：X0 为过 P 点的平行圈所截的中央子午线距赤道弧长；NP 为 P 点的卯酉圈半径；l 为过 P 点经度与投影带中央子午线经度之差；B 为 P 点的纬度；t=tgB；2=(e)2cos2B 。
