测绘基准与坐标系统转换

大地基准、坐标系统转换及软件系统开发报告主要内容大地基准简介独立坐标系的建立软件系统的简介坐标系统转换一、大地基准简介大地水准面：通过平均海水面的水准面。参考椭球面：最佳拟合于某一区域性大地水准面的旋转椭球面。正常椭球：参考椭球具有确定的质量和绕短轴恒定的角速度。平均地球椭球：其表面最佳拟合于全球大地水准面的旋转椭球。1.1地固坐标系与地球椭球赤道平面OPM大地经度L大地纬度BnLB起始子午面（首子午面）SZXY参心坐标系：以参考椭球为旋转椭球（实质是一个二维的曲面坐标系）地心坐标系：以地球质心为椭球中心二者均属于地固坐标系1954年北京坐标系：是前苏联1942年普尔科夫坐标系的延伸。它是依照1953年我国东北边境内三角点与前苏联天文大地网联测形成一等三角锁经局部平差后确定的、随后扩展、加密而遍及全国。（1）克拉索夫斯基椭球(a=6378245m,a=1:298.3)（2）多点定位：（3）参考定向时：令（4）大地点高程是以1956年青岛验潮站求出的黄海平均海水面为基准。（5）对应的参考椭球面与我国大地水准面存在着自西向东明显的系统性倾斜，在东部最大达65m，全国平均29m。（6）椭球定向不明确，短轴既不指向国际协议原点CIO，也不指向我国地极原点JYD1968。一般而言，对于精度要求较高的控制网不宜采用两个以上国家控制点的54系的坐标作为固定数据。54点控制点位移、局部地壳运动等各种情况均会使实际点位与当初所赋的坐标不一致。0zyx1.2我国的参心坐标系—1954北京坐标系（1）参考椭球：UGG1975(a=6378685m,a=1:291.6)（2）多点定位（3）定向明确，短轴指向我国地极原点JYD1968。（4）大地点高程是以1956年青岛验潮站求出的黄海平均海水面为基准。（5）整体平差，计算5万多点成果。1.3我国的参心坐标系—1980年国家大地坐标系1.4我国的地心坐标系—CGCS2000ChinaGeodeticCoordinateSystem2000，CGCS20001954北京坐标系和1980西安坐标系在我国的经济建设和国防建设中发挥了巨大作用。同时也应当看到，随着情况的变化和时间的推移，这两个以经典测量技术为基础的局部大地坐标系，目前已经不能适应科学技术特别是空间技术发展，不能适应我国经济建设和国防建设需要。我国大地坐标系的更新换代，是经济建设、国防建设、社会发展和科技发展的客观需要。以地球质量中心为原点的地心大地坐标系，是当今空间时代全球通用的基本大地坐标系。以空间技术为基础的地心大地坐标系，是我国新一代大地坐标系的适宜选择。地心大地坐标系可以满足大地测量、地球物理、天文、导航和航天应用以及经济、社会发展的广泛需求。多年以来，我国测绘、地震部门和科学院有关单位为建立我国新一代大地坐标系作了大量基础性工作，近年又先后建成全国GPS一、二级网，国家GPSA、B级网，中国地壳运动观测网络和许多地壳形变网，为地心大地坐标系的实现奠定了较好的基础。1.4我国的地心坐标系—CGCS2000大地坐标系的定义包括坐标系的原点、3个坐标轴的指向、尺度以及地球椭球的4个基本常数的定义。2000国家大地坐标系的原点为包括海洋和大气的整个地球的质量中心；2000国家大地坐标系的Z轴由原点指向历元2000.0的地球参考极的方向，该历元的指向由国际时间局给定的历元为1984.0的初始指向推算，定向的时间演化保证相对于地壳不产生残余的全球旋转，X轴由原点指向格林尼治参考子午线与地球赤道面(历元2000.0)的交点，Y轴与Z轴、X轴构成右手正交坐标系。采用广义相对论意义下的尺度。2000国家大地坐标系采用的地球椭球基本常数短半径b(m)6356752.31414极曲率半径c(m)6399593.62586第一偏心率e0.0818191910428第一偏心率平方e20.00669438002290第二偏心率e’0.0820944381519第二偏心率平方0.006739496775481/4子午圈的长度Q(m)10001965.7293椭球平均半径R1(m)6371008.77138相同表面积的球半径R2(m)6371007.18092相同体积的球半径R3(m)6371000.78997椭球的正常位U0(m2s-2)62636851.7149动力形状因子J20.001082629832258球谐系数J4-0.00000237091126球谐系数J60.00000000608347球谐系数J8-0.000000000014270.00344978650678赤道正常重力值γe（伽）9.7803253361两极正常重力值γp（伽）9.8321849379正常重力平均值γ（伽）9.7976432224纬度45度的正常重力值γ45°（伽）9.806197769522/mabGM2000国家大地坐标系采用的地球椭球基本常数1.5WGS84坐标系—WorldGeodeticSystem—1984CoordinateSystem一种国际上采用的地心坐标系。坐标原点为地球质心，其地心空间直角坐标系的Z轴指向BIH（国际时间）1984.O定义的协议地球极（CTP)方向，X轴指向BIH1984.0的零子午面和CTP赤道的交点，Y轴与Z轴、X轴垂直构成右手坐标系，称为1984年世界大地坐标系统。WGS-84采用的椭球是国际大地测量与地球物理联合会第17届大会大地测量常数推荐值，其四个基本参数长半径：a=6378137±2（m）；地球引力常数：GM=3986005×108m3s-2±0.6×108m3s-2；正常化二阶带谐系数：C20=-484.16685×10-6±1.3×10-9；J2=108263×10-8地球自转角速度：ω=7292115×10-11rads-1±0.150×10-11rads-11.6国家统一坐标系《国家三角测量和精密导线测量规范》规定：“所有大地测量的观测成果，都须归化到参考椭球面上，……所有国家大地点均按高斯正形投影计算其在6°带内的平面直角坐标(一般称为高斯一克吕格平面坐标)。在1：1万和更大比例尺测图的地区，还应加算其在3°带内的平面直角坐标”，根据这种规定所采用的高斯投影6°带或3°带坐标系统，我们通常称作国家统一坐标系统。工程控制网采用国家统一坐标系统的结果，常常改变控制网各条边的真实长度，引起长度的变形，这对于大比例尺地形测图和工程测量是十分不利的。为了有效地控制投影长度变形，就需要分析长度变形的来源和容许数值，国家统一坐标系统的适用程度和范围，研究独立坐标系统建立方法。二、独立坐标系的建立工程控制网是直接服务于城镇和工矿地区大比例尺测图和工程施工测量的，因此，由控制网所提供的距离应尽可能保持其真实性。这样，地面施测的距离可以直接绘图，图纸上量取的距离也可直接标设于实地。所以就控制网的实用性而言“长度综合变形愈小愈好”。对此，《城市测量规范》、《工程测量规范》均对控制网的长度综合变形的容许范围作了明确规定，一致确立了平面控制网的坐标系统应该保证长度综合变形不超过2．5cm／km(相对变形为1：40000)这一原则。这样的长度变形，与四等平面控制网边长的必需精度相适应，对于测量精度为1：5000--1：2000的施工放样，也能起到良好的控制作用。将实地测量的真实长度归化到国家统一的椭球面上时，应加如下改正数：式中：RA——长度所在方向的椭球曲率半径；Hm——长度所在高程面对于椭球面的高差；s——实地测量的水平距离。然后再将椭球面上的长度投影至高斯平面，加入如下改正数：式中；R—测区中点的平均曲率半径ym—距离的两端点横坐标的均值地面上的一段距离，经过2次改正计算，被改变了真实长度。这种高斯投影平面上的长度与地面长度之差，我们称之为长度综合变形，其计算公式为：sRHsAmSRySm222sRHSRyAmm2222.1长度变形的产生和容许数值为了方便计算，又不致损害必要精度，可以将椭球视为圆球，取圆球半径RA=R=6371km，又取不同投影面上的同一距离近似相等，即S=s，将上式写成相对变形的形式，则为式中，y表示测区中心的横坐标(自然值)，H表示测区平均高程，y与H均以km作单位。采用国家统一坐标系统所产生的长度综合变形，与测区所处投影带内的位置和测区平均高程（或生产水平面高程）有关。利用上述公式；可以便捷地计算出测区所用坐标系统的长度相对变形大小。5210)7.1500123.0(HyssRHSRyAmm2222.1长度变形的产生和容许数值将长度综合变形的容许数值1：40000代入上式，即可得到方程取测区中心的坐标为横轴，取测区平均高程H为纵轴，就可以画出相对变形恒为容许数值的两条曲线。这两条曲线就是适用于控制测量的投影带范围临界线，或者说两条曲线之间的区域就是适用于城镇测图和工程测量的投影带范围，利用右图可以直观形象地判断国家3°带统一坐标系是否适合于本测区的需要。如果根据本测区的平均高程和3°带的y坐标所确定的位置，处于两曲线以外的“不适用区”就应该考虑另行选择坐标系统。16.0)10(78.042yH2.2国家统一坐标系统的局限性为了保证投影的长度综合变形不超过2．5cm／km(1：40000)，必将极大降低3°带国家统一坐标系统的适用范围。从图可以看出，在3°投影带的大部分区域，长度综合变形均超出了《规范》中规定定的范围，均不适宜测图和工程测量的需求。所以如何选择能够抵偿长度综合变形的坐标系统，就是一个必须解决的问题。如果说高斯投影分带是限制长度变形的有效措施，那么正确地选择工程控制网的坐标系统，则是抵偿长度综合变形的有效途径。以椭球面作为计算表面(投影面)、按正形投影3°分带计算的平面坐标称为国家统一坐标系统，那么对投影面或投影带作出其它选择，计算出的坐标就不再属于国家统一坐标系统，我们称这种坐标系统为局部坐标系统，有时也称为独立坐标系统或实用坐标系统。下面讨论几种实用（局部）坐标系统可供选择的方案。2.3实用工程坐标系建立方法2.3.1抵偿坐标系—改变投影面不改变中央子午线选择“抵偿高程面”作为投影面，将距离由较高的高程面化算至较低的椭球面时，长度总是减小的；将椭球面上的距离化算至高斯平面时，长度总是增加的。所以两个投影过程对长度变形具有抵偿的性质。如果适当选择椭球的半径，使距离化算到这个椭球面上所减小的数值，恰好等于由这个椭球面化算至高斯平面所增加的数值，那么高斯平面上的距离同实地距离就一致了。这个适当半径的椭球面，就称为“抵偿高程面”。欲使长度综合变形得以抵偿，必须：代入数据得：H＝785y2（y以百公里，H以m作单位）例如，某地中心在高斯投影3°带的坐标y＝91km，该地平均高程为400m，按式算得：H＝785×0.91＝650m，H抵＝400m-650m＝-250mSRysRHmAm2223.2.1抵偿坐标系—改变投影面不改变中央子午线抵偿面位置确定后，就可以选择其中一个国家大地点作“原点”，保持它在3°带的国家统一坐标值(x0，y0)不变，而将其它大地控制点坐标(x，y)换算到抵偿高程面相应的坐标系中去。换算公式为：式中：R为该地平均纬度处的椭球平均曲率半径。经过上式换算的大地控制点坐标就可以作为控制测量的起算数据。抵偿高程面作为投影面30带高斯平面直角坐标系称为抵偿坐标系。RHyyyyRHxxxx抵抵抵抵)()(00保持国家统一的椭球面作投影面不变，选择“任意投影带”，按高斯投影计算平面直角坐标。不同投影带的出现，是因为选择了不同经度的中央子午线的缘故。合理选择中央子午线的位置，使长度投影到该投影带所产生的变形，恰好抵偿这一长度投影到椭球面所产生的变形，此时高斯投影平面上的长度仍和实地长度保持一致。我们称这种抵偿长度变形的投影带为“任意投影带”。为了确定任意投影带的中央子午线位置，需要引入经度差l”。取高斯投影坐标正算公式式中B、L——测区中