《最终给第九章GPS数据处理.ppt》由会员分享，可在线阅读，更多相关《最终给第九章GPS数据处理.ppt（115页珍藏版）》请在三一办公上搜索。
1、9.1 概述9.2 GPS基线向量的解算9.3 GPS定位成果的坐标转换(P158)9.4 基线向量网平差9.5 GPS高程,第九章 GPS定位测量数据处理（DATA PROCESSING FOR GPS）,基本处理过程：（1）数据传输（2）数据分流 从原始数据中，剔除无效观测值和冗余信息，形成各种数据文件，如星历文件、载波相位和伪距观测文件、测站信息文件。（3）GPS数据的预处理 对数据进行平滑滤波检验，剔除粗差；统一数据文件格式，将不同类型接收机的数据记录格式统一为标准化的文件格式，探测周跳，修复观测值。（4）基线向量解算 一般先采用三差模型法对基线向量进行预求解，然后再采用双差模型对基线
2、向量进行精确求解。对于20 km以下的基线，常采用所谓的固定双差解，即整周未知数取为整数后的基线平差解。而对于30 km以上的基线，一般采用浮点双差解，即整周未知数不取整，以实数作为整周值，所以也可称为实数解。,（5）GPS网平差与坐标转换概述 由同步观测和异步观测的基线向量互相联结构成GPS网，称为GPS基线向量网。由于存在观测 误差，网中由不同时段观测的基线向量组成的闭合图形存在不符值(闭合差)。因此，应在WGS-84坐标系统下，以GPS基线向量及其相应的方差阵作为观测信息，对GPS网进行平差计算，消除不符值，获得网中点的平差后的三维坐标、基线边长的平差值、基线向量观测值改正数及其对观测值
3、、点位坐标的精度评定。,9.2 GPS基线向量的解算,9.2.2 法方程的组成及解算,9.2.3 精度评定,9.2.4 基线向量解算结果分析1.观测值残差分析2.基线长度的精度3.双差固定解与双差实数解,9.2.5 GPS基线向量解算示例,9.3 GPS定位成果的坐标转换(P158),WGS-84大地坐标系是GPS卫星定位系统采用的大地坐标系，因而，所用利用GPS接收机进行测量计算的成果均属于WGS-84坐标系；而实用的测量成果往往是属于某一国家坐标系或地方坐标系（或叫局部的、参考坐标系）；参考坐标系与WGS-84坐标系之间一般存在着平移和旋转的关系。,9.3 GPS定位成果的坐标转换（续）,
4、WGS-84坐标转换为国家或地方坐标系的方法（8454）1.利用已知重合点的三维直角坐标进行坐标转换；2.利用已知重合点的三维大地坐标进行坐标转换；3.利用已知重合点的二维高斯平面坐标进行坐标转换；4.利用已知重合的二维大地坐标进行坐标转换。,本节仅讨论1、3两种方法，即不同空间直角坐标系或不同平面直角坐标系统之间的转换。1.利用已知重合点的三维直角坐标将GPS点的WGS-84坐标转换为国家坐标系中的坐标。用七参数法实现坐标转换；局部地区应用坐标差求解转换参数的方法（p27）；在GPS网的约束平差中实现坐标转换。2.利用已知重合点的二维高斯平面坐标将GPS点的WGS-84转换为国家坐标系中的坐
5、标。将 GPS点的大地坐标(B、L)按 WGS-84参考椭球和高斯正形投影公式换算为高斯平面坐标(x、y)；利用重合点（至少两个）的两套平面坐标值按平面坐标系统之间的转换方法将GPS点的高斯平面坐标转换为国家坐标系高斯平面坐标,9.3 GPS定位成果的坐标转换（续）,预备知识,常规测量,平面坐标系,高程坐标系,GPS测量,三维地心坐标系,地心空间直角坐标系、椭球大地坐标系,9.2 国家坐标系与地方独立坐标系,9.2.1 旋转椭球与参心坐标系 水准面：在地球重力场中，当水处于静止时的表面必定与重力方向(即铅垂线方向)处处正交。我们称这个与铅垂线正交的静止水平面为水准面。大地水准面：假设海水面处于
6、静止平衡状况，并将它一直沿伸到地球陆地内部形成一个闭合的水准面，用来表示地球的形状，我们将这个水准面称为大地水准面。,大地水准面是对地球的物理逼近，它可以较真实地反映地球的形状，但是地壳内部物质密度分布的不均匀，造成地面各点重力大小和方向不同，因此，与铅垂线处处正交的大地水准面是起伏不平的，因而它也很难以用简单的数学模型描述。要用它作为各种地面测量数据的计算基准面比较困难，必须寻找一个简单的适合测量计算的基准面。,大地水准面相当接近于一个规则的具有微小扁率的数学曲面旋转椭球。旋转椭球可用两个几何参数确定，即为椭球的长半径a和扁率f。这两个参数解决了椭球的形状和大小。为了将地面测量数据归算到椭球
7、面上，仅仅知道它的形状和大小是不够的，还必须确定它与大地水准面的相关位置，也就是所谓的椭球定位和定向。另外，为了从几何特性和物理特性两个方面来研究全球的形状，则还要使椭球与全球大地水准面结合最为密切。,为了研究局部球面的形状，且使地面测量数据归算至椭球的各项改正数最小，各个国家和地区分别选择和某一局部区域的大地水准面最为密合的椭球建立坐标系。这样选定和建立的椭球称为参考椭球，对应的坐标系称为参心坐标系。显然，该坐标系的中心一般和地球质心不一致，所以参心坐标系又称为非地心坐标系、局部坐标系或相对坐标系，由于参心坐标系处理局部区域数据带来的变形较小，所以，参心坐标系至今对大地测量仍有重要作用。同样
8、，参心坐标系可分为参心空间直角坐标系和参心大地坐标系。,参心空间直角坐标系 参心空间直角坐标系是：（1）以参心O为坐标原点；（2）Z轴与参考椭球的短轴(旋转轴)相重合；（3）X轴与起始子午面和赤道的交线重合；（4）Y轴在赤道面上与X轴垂直，构成右手直角坐标系 O-XYZ。地面点P的点位用(X，Y，Z)表示。,参心大地坐标系 参心大地坐标系是以参考椭球的中心为坐标原点，椭球的短轴与参考椭球旋转轴重合。大地纬度B以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B；大地经度L以过地面点的椭球子午面与起始子午面之间的夹角为大地经度L；大地高H地面点沿椭球法线至椭球面的距离为大地高H(如图)。地面点的点
9、位用(B，L，H)表示。,B=19.4705047353，L=110.0143500924，H=49.727,X=-2054198.29421，Y=5934981.51372，Z=2192959.41759,WGS-84椭球：,在同一参心坐标系中，地面点的参心空间直角坐标与相应的参心大地坐标之间存在如下的互换关系。（1）参心大地坐标转换为参心空间直角坐标：式中，N为椭球面卯酉圈的曲率半径，e为椭球的第一偏心率：,（2）空间直角坐标转换为大地坐标：,确定参考椭球是建立参心坐标系的主要依据。通常包括确定参考椭球的形状和大小，确定它的空间位置(参考椭球的定位与定向)，以及确定大地原点T的大地纬度BT
10、、大地经度LT及它至一相邻点的大地方位角AT。参考椭球的定位和定向是通过确定大地原点的大地经纬度、大地高和大地方位角来实现的，参考椭球一般采用“双平行”定向条件，即要求椭球的短轴与地球某一历元的自转轴平行，起始大地子午面与起始天文子午面平行。,9.2.2 P54北京和C80西安国家坐标系,目前我国常用的两个国家坐标系1954年北京坐标系和1980年西安坐标系，均是参心坐标系。,9.2.2.1 1954年北京坐标系 1954北京坐系采用了前苏联的克拉索夫斯基椭球体，其椭球参数是：长半轴a为9 378 245m，扁率f为1/298.3，其原点为原苏联的普尔科沃。1954年北京坐标系虽然是苏联194
11、2年坐标系的延伸，但也还不能说它们完全相同。因为该椭球的高程异常是以苏联1955年大地水准面重新平差结果为起算数据，按我国天文水准路线推算而得。而高程又是以1959年青岛验潮站的黄海平均海水面为基准。,1954年北京坐标系建立之后，在这个系统上，多年来，我国用该坐标系统完成了大量的测绘工作，获得了许多的测绘成果，在国家经济建设和国防建设的各个领域中发挥了巨大作用。但是，随着科学技术的发展，这个坐标系的先天弱点也显得越来越突出，难以适应现代科学研究、经济建设和国防尖端技术的需要，它的缺点主要表现在：(1)克拉索夫斯基椭球参数同现代精确的椭球参数相比，误差较大，长半径约大105109m，这不仅对研
12、究地球几何形状有影响，特别是该椭球参数只有两个几何参数，不包含表示物理特性的参数，不能满足现今理论研究和实际工作的需要，对于发展空间技术也带来诸多不便。,(2)椭球定向不明确，即不指向国际通用的CIO极，也不指向目前我国使用的JYD极，椭球定位实际上采用了前苏联的普尔科沃定位，该定位椭球面与我国的大地水准面呈系统性倾斜。东部高程异常达90余米。而我国东部地势平坦、经济发达，要求椭球面与大地水准面有较好的密合，但实际情况与此相反。(3)该坐标系统的大地点坐标是经局部平差逐次得到的，全国天文大地控制点坐标值实际上连不成一个统一的整体。不同区域的接合部之间存在较大隙距，同一点在不同区的坐标值相差12
13、m，不同区域的尺度差异也很大。而且坐标传递是从东北至西北西南，前一区的最弱点即为后一区的坐标起算点，因而坐标积累误差明显，这对于发展我国空间技术、国防建设和国家大规模经济建设不利，因此有必要建立新的大地坐标系统。,9.2.2.2 1980年西安坐标系 1978年，我国决定建立新的国家大地坐标系统，并且在新的大地坐标系统中进行全国天文大地网的整体平差，这个坐标系统定名为1980年西安大地坐标系统。1980年西安坐标系的大地原点设在我国的中部，处于陕西泾阳永乐镇，椭球参数采用1975年国际大地测量与地球物理联合会推荐值，它们为：椭球长半径 a=9378140m；重力场二阶带球谐系数 J2=1.08
14、29310-3；地心引力常数 GM=3.9890051014m3/s;地球自转角速度=7.29211510-5rad/s。因而可得80椭球两个最常用几何参数为:a=9378140m；f=1/298.257。,椭球定位按我国范围高程异常值平方和最小为原则求解参数。椭球的短轴平行于由地球质心指向1998.0地极原点(JYD)的方向，起始大地子午面平行于格林尼治天文台子午面。长度基准与国际统一长度基准一致。高程基准以青岛验潮站1959年黄海平均海水面为高程起算基准，水准原点高出黄海平均海水面72.289 m。1980年西安大地坐标系建立后，利用该坐标进行了全国天文大地网平差，提供全国统一的、精度较高
15、的1980年国家大地点坐标，据分析，它完全可以满足1/5000测图的需要。,9.2.2.3 新1954年北京坐标系 由于1980年西安坐标系与1954年北京坐标系的椭球参数和定位均不同，因而大地控制点在两坐标系中的坐标存在较大差异，最大的达100m以上，这将引起成果换算的不便和地形图图廓和方格线位置的变化，且已有的测绘成果大部分是1954年北京坐标系下的。所以，作为过渡，产生了所谓的新1954年北京坐标系。新1954年北京坐标系是通过将1980年西安坐标系的三个定位参数平移至克拉索夫斯基椭球中心，长半径与扁率仍取克拉索夫斯基椭球几何参数。而定位与1980年大地坐标系相同(即大地原点相同)，定向
16、也与1980椭球相同。因此，新1954年北京坐标系的精度和1980年坐标系精度相同，而坐标值与旧1954年北京坐标系的坐标接近。,9.2.3 地方独立坐标系 在我国许多城市测量与工程测量中，若直接采用国家坐标系，则可能会由于远离中央子午线，或由于测区平均高程较大，而导致长度投影变形较大，难以满足工程上或实用上的精度要求。另一方面，对于一些特殊的测量，如大桥施工测量、水利水坝测量、滑坡变形监测等，采用国家坐标系在实用中也会很不方便。因此，基于限制变形，以及方便实用、科学的目的，在许多城市和工程测量中，常常会建立适合本地区的地方独立坐标系。,建立地方独立坐标系，实际上就是通过一些元素的确定来决定地
17、方参考椭球与投影面。地方参考椭球一般选择与当地平均高程相对应的参考椭球，该椭球的中心、轴向和扁率与国家参考椭球相同，其椭球半径1增大为：1=+1,1=Hm+0式中：Hm为当地平均海拔高程，0为该地区的平均高程异常。而地方投影面的确定中，选取过测区中心的经线或某个起算点的经线作为独立中央子午线。以某个特定方便使用的点和方位为地方独立坐标系的起算原点和方位，并选取当地平均高程面Hm为投影面。,9.2.4 高斯平面直角坐标系和UTM,大地测量建立的大地坐标的重要作用之一是为测图服务，传统地图均为平面图，作为测图控制的大地点的坐标也必须是平面坐标。因此，需要将椭球面上各点的大地坐标，按照一定的数学规律
18、投影到平面上成为平面直角坐标。,由于地球椭球面是不可展的曲面，无论采用什么数学规律投影都会产生变形。因此，只能按照满足某种特定需要与用途，对一些变形加以限制，使其减小到适当程度，甚至为零。按变形性质，我们可以将投影分为等角投影、等面积投影、等距离投影以及任意投影。,等角投影也叫正形投影、相似投影。即该投影在小区域范围内使平面图形与椭球面上的图形保持相似。目前世界各国采用最广泛的高斯投影和墨卡托投影(UTM)均是正形投影。,高斯投影和墨卡托(UTM)投影具有如下特征：,(1)椭球面上任一角度，投影到平面上后保持不变。,(3)高斯投影的中央子午线长度比 m0=1，而UTM投影的m0=0.9999。
19、,(2)中央子午线投影为纵坐标轴。,3度带是在9度带的基础上分带，其带号n与相应的中央子午线经度L0的关系为：,高斯投影和UTM投影是正形投影，因此，其角度没有变形，而长度除中央子午线外均存在变形，距中央子午线越远，长度变形越大。为了限制长度变形，根据国际测量协会规定，将全球按一定经差分成若干带。,我国采用6度带或3度带。6度带是自零度子午线起每隔经度差9度自西向东分带，带号n与相应的中央子午线经度L0的关系为,L0=6n-3；n=1/6(L0+3),L0=3n；n=1/3*L0=2n-1,高斯克吕格投影,在我国范围内:3度带号从25-45，6度带号从13-23。,9.3 GPS定位测量中的坐
20、标转换,不同的测量成果均对应于各自的坐标系。GPS定位结果属于协议地球地心坐标系，即WGS-84坐标系，且通常以空间直角坐标(,)s或以椭球大地坐标(B，L，H)s的形式给出。而实用的常规地面测量成果或是属于国家的参心大地坐标系，或是属于地方独立坐标系。因此必须实现GPS成果的坐标系的转换。,另外，GPS相对定位所求得的GPS基线向量通常是以WGS-84坐标差的形式表示，对于这种特殊的坐标表示形式，应考虑其相应的转换模型。为了与传统测量成果一致，常将GPS成果投影到平面，形成GPS二维坐标系成果，因此还应考虑二维坐标转换。,9.3.1 空间直角坐标系与椭球大地坐标系的关系,正解,反解,9.3.
21、2 三维坐标转换模型,设有两三维空间直角坐标系OT-XTYTZT和OS-XSYSZS有图9-所示的关系。1、其坐标系原点不一致，存在三个平移量X0,Y0,Z0；2、且通常各坐标轴之间相互不平行，对应的坐标轴之间存在三个微小的旋转角x,y,z；3、两个坐标系的尺度也不一致，设OT-XTYTZT 的尺度为1，而OS-XSYSZS 的尺度为1+m。坐标转换模型有多种，应用最广的是7参数转换模型，主要有布尔沙模型和莫洛金斯基模型。,9.3.2.1 布尔沙(Bursa)模型 由下图知，任意点Pi在两坐标系中的坐标之间有如下关系：,布尔沙(Bursa)模型,在两个坐标系中，如果认为，受尺度和旋转影响的只是
22、任意点Pi与某一参考点Pk的坐标差，则有关系：,9.3.2.2 莫洛金斯基(MOOAHCK)模型,或者,莫洛金斯基模型,7参数坐标转换模型，除了布尔莎模型、莫洛金斯基模型外，还有维斯模型、范氏模型、武测模型，这些模型在表达形式上虽不尽相同，但参数间存在着明确的解析关系，可以相互进行转换，用它们分别换算其他点的坐标时，结果完全相同。因此，这几种转换模型是等价的。,GPS相对定位的基线向量通常是以三维坐标差的形式表示，按照上述布尔莎模型(或其他模型)列出两个点的坐标转换方程，将两式相减，就得到两点间的三维坐标差的转换模型为：,9.3.3 三维坐标差转换模型,记xij=xj-xi,yij=yj-yi
23、,zij=zj-zi，上式也可写为：,由于坐标差与平移参数无关，所以，三维坐标差中仅包含三个旋转参数和一个尺度参数，且由以上任何一个7参数转换模型得到的坐标转换模型完全相同。,无论利用哪种坐标转换模型，均必须已知相应的转换参数。如果不知道两个坐标系间的转换参数，则需根据两坐标系在共同点(公共点)的坐标(X,Y,Z)iS和(X,Y,Z)iT(i=1,2,3)，代入转换模型反求两个坐标系间的转换参数，然后利用所求得的转换参数再回代到模型中对另一部分点进行坐标转换。一般常用联合平差的方法求转换参数。可以将布尔莎转换模型写成如下形式：,9.3.4 联合平差确定转换参数,或者 Xis=CiT+XiT式中
24、 Xis=x y zTis，XiT=X Y ZTiT T=x0 y0 z0mxyzT是要确定的7个转换参数组成的向量。要确定7个参数，至少需要同时知道三个公共点在两坐标系的坐标值，利用最小二乘法对参数T进行求解。由于这两个坐标系分别属于卫星网的地心坐标系(S)和地面网的参心坐标系(T)，所以将这个求解过程称为卫星网与地面网的联合平差。,9.4 GPS网的三维平差（基线向量网平差）,GPS网是由GPS相对定位求得的基线向量构成的空间基线向量网，并在GPS网平差时，将这些基线向量及其协方差作为网平差的基本观测量。GPS空间三维基线向量网平差常采用以下几种平差类型:(1)三维无约束平差：GPS基线向
25、量本身已隐含尺度基准和方位基准，因此在三维平差中可只选某一点的固定坐标进行网平差，即所谓的无约束平差。三维无约束平差是GPS网平差中不可缺少的步骤，它可发现基线向量中存在的粗差、系统误差。通过检验发现基线向量随机模型误差，可客观评价GPS网本身的内符合精度。,(2)GPS网三维约束平差：以国家大地坐标系或地方坐标系的某些点的固定坐标、固定边长和方位为网的基准，将其作为平差中的约束条件，并在平差中考虑GPS网与地面网之间的转换系数。因此，这种形式的平差是在地面参考坐标系中进行的，故称为GPS三维约束平差。该平差后获得网的坐标已是国家大地坐标系或地方坐标系的坐标，因而约束平差是目前GPS网成果转换
26、行之有效的方法。,(3)三维联合平差：平差中除了GPS基线向量观测值和地面基准约束数据外，还包含了地面常规网观测值，如边长、方向、天文方位角、天顶距、水准高差乃至天文经纬度，将这些数据一并进行平差，也就是GPS网和地面观测数据的联合平差，其平差后网中点的坐标仍属地面坐标系框架下的。,GPS三维平差的主要流程图,GPS三维平差的主要流程图(续）,GPS网三维平差中，首先应进行三维无约束平差，平差后通过观测值改正数检验发现基线向量中是否存在粗差，并剔除含有粗差的基线向量，再重新进行平差，直至确定网中没有粗差后，应对单位权方差因子进行2检验，判断平差的基线向量随机模型是否存在误差，并对随机模型进行改
27、正，以提供较为合适的平差随机模型。再对GPS网进行约束平差或联合平差后，还应对平差中加入的转换参数进行显著性检验，对于不显著的参数应剔除，以免破坏平差方程的性态。,9.4 GPS基线向量网的二维平差,GPS基线向量网二维平差应在某一参考椭球面上，或是在某一投影平面坐标系上进行。因此，平差前，首先须将GPS三维基线向量观测值及其协方差阵转换投影至二维平差计算面。也就是从三维基线向量中提取二维信息，在平差计算面上构成一个二维GPS基线向量网。,GPS基线向量网二维平差也可分为无约束平差、约束平差和联合平差三类，平差原理及方法均与三维平差相同。由二维约束平差和联合平差获得的GPS平面成果，就是国家坐
28、标系下或地方坐标系下具有传统意义的控制成果。在平差中的约束条件往往是由地面网与GPS网重合的已知点坐标，这些作为基准的已知点的精度或它们之间的兼容性是必须保证的。否则由于基准本身误差太大互不兼容，将会导致平差后的GPS网产生严重变形，精度大大降低。因此，在平差结束前，应通过检验发现并淘汰精度低且不兼容的地面网已知点，再重新平差。,9.4.1 GPS基线向量网的二维投影变换 三维基线向量转换成二维基线向量，应避免地面网不准确的大地高引起尺度误差和网变形，保证GPS网转换后整体及相对几何关系不变。因此，可采用在一点上实行位置强制约束，在一条基线的空间方向上实行定向约束的三维转换方法，也可以是在一点
29、上实行位置强制约束，在一条基线的参考椭球面投影的法截弧和大地线方向上实行定向约束的准三维转换方法。转换后的GPS网与地面网在一个基准点上和一条基线上的方向完全重合一致，两网之间只存在尺度比差和残余定向差。,设地面控制位置基准点在国家大地坐标系中的大地坐标为(B0T、L0T、H0T)(H0T=h0+0)，由大地坐标与空间三维直角坐标关系式可得该点在国家空间直角坐标系下的坐标(X0T，Y0T，Z0T)。假定网中基准点的坐标为(B0，L0，H0)，不同于基准点的其他点为(Bi，Li，Hi)。而该基准点在GPS网的三维直角坐标为(X0S，Y0S,Z0S)，由此可求得GPS网平移至地面测量控制网基点的平
30、移参数为：X=X0T-X0S，Y=X0T-Y0S，Z=Z0T-Z0S,9.4.1.1 GPS基线向量网至地面网的平移,于是，可将GPS网中其他各点坐标经下式平移到国家大地坐标系中：XTi=XSi+X YTi=YSi+Y ZTi=ZSi+Z 利用大地坐标与空间直角坐标的反算公式可得各点在国家大地坐标系中的大地坐标Bi、Li、Hi。,9.4.1.2 三维GPS网至国家大地坐标系的二维投影变换 由平移变换已将GPS网与地面控制网在基准点上(见图9-，1号点为基准点)实现了重合，为使GPS网与地面控制网在方位上重合一致，可利用椭球大地测量学中的赫里斯托夫第一类微分公式。该公式给出当基准数据发生变化(起
31、始点的大地坐标dB、dL，方位dA、长度dS)时，相应的其他大地点坐标的变化值，实现两网在同一椭球面上的符合。,图9-7 GPS网的转换过程,上式即为赫里斯托夫第一微分式。,式中,在该公式中，对于已经平移变换、两网在基准上重合的网，有dB0=0，dL0=0。同时由于在进行三维至二维的投影变换时，往往难以准确确定两网的尺度差异，因此可将此变换留待约束(联合)平差时考虑，因而此时可设dS0=0，那么，此处赫里斯托夫第一类微分式就化简为：,式中,根据两网起始坐标方位角之差dA=A0T-A0S。由上式就可得GPS网各点在国家大地坐标系内与地面网起始基准点一致、起始方位一致(见图9-7中12方位)的坐标
32、。B1=B1+dB1 L1=L1+dL1 为了平差计算及科研生产实用方便，二维平差通常是在平面上进行。可利用高斯投影正算式将GPS各点由参心椭球坐标投影到高斯平面坐标系。,9.4.1.3 三维基线向量协方差阵至二维高斯平面协方差阵 除了将三维GPS网投影到二维平面上外，还应把相应的协方差阵变换到二维高斯平面。由空间直角坐标与椭球大地坐标的关系可得任一条基线的大地坐标差关系为：,式中：ai,bi,ci(i=B,L,H)为一阶偏导系数。,可将上式展开成三维级数式，并通过对其求偏导得到空间直角坐标差与大地坐标差之间的全微分式：,由于偏导系数矩阵是可逆的，于是可惟一地得到大地坐标差关于直角坐标差的微分
33、关系式：,按协方差传播律，可得到大地坐标差与直角坐标差之间的协方差转换公式：,而由高斯正算公式可得平面直角坐标差与椭球坐标的全微分式：,由此，可得三维空间基线向量到二维高斯平面的协方差阵:,9.4.2 GPS基线向量网的二维平差 由上述转换方法可将GPS基线向量及其协方差阵转换到二维国家大地平面坐标系下，网的二维平差即可在该平面上进行。,9.4.2.1 GPS二维约束平差的观测方程和约束条件方程 设二维基线向量观测值为Xij=(xij,yij)T，而待定坐标改正数dX=(dxi,dyi)T和尺度差参数m以及残余定向差参数d为平差未知数，则GPS基线向量的观测误差方程为：,当网中有已知点的坐标约
34、束时，则GPS网中与该点重合的点的基线向量的坐标改正数为零，即：当网中有边长约束时，则边长约束条件方程为：,这里的Sij即为GPS网的尺度标准。,当网中有已知方位角约束时，则其约束条件方程为：,9.4.2.2 GPS基线向量网与地面网的二维联合平差 GPS基线向量网与地面网的二维联合平差是在上述的GPS基线观测方程及坐标、边长、方位的约束条件方程的基础上再加上地面网的观测方向和观测边长的观测方程。方向观测值误差方程：,边长观测值误差方程：,GPS基线向量网的二维平差方法和平差过程均与三维网相同。,9.5 GPS高程,GPS相对定位高程方面的相对精度一般可达(23)-9；在绝对精度方面，实验表明
35、，对于10km以下的基线边长，可达几个厘米，如果在观测和计算时采用一些消除误差的措施，其精度将优于1cm。,高程系统简介：正高(HG)地面点沿通过该点的铅垂线到达大地水准面的距离；基准面大地水准面 正常高(Hr)地面点沿通过该点的铅垂线到似大地水准面的距离；基准面似大地水准面 大地高(HD)地面点沿通过该点椭球面法线到椭球面的距离；基准面参考椭球面,大地高与正常高的关系：即：大地高=正常高+高程异常,9.5.1 GPS水准 众所周知，实际应用中的地面点高程是以似大地水准面为起算面的正常高，而GPS高程是以WGS-84椭球面为基准的大地高。由前所述已知两者之间相差高程异常，显然，如果知道了各GP
36、S点的高程异常，则可由各GPS点的大地高求得各点的正常高。,我国似大地水准面主要是采用天文重力方法测定的，其精度为1m左右，因此很难直接由GPS大地高求得正常高。目前在小区域范围内，常采用GPS水准的方法较为精确地计算GPS点的正常高。所谓GPS水准就是在小区域范围的GPS网中，用水准测量的方法联测网中若干GPS点的正常高(这些联测点称为公共点)，那么根据各GPS点的大地高就可求得各公共点上的高程异常。然后由公共点的平面坐标和高程异常采用数值拟合计算方法，拟合出区域的似大地水准面，即可求出各点高程异常值，并由此求出各GPS点的正常高。,目前，国内外GPS水准主要是采用纯几何的曲面拟合法，即根据
37、区域内若干公共点上的高程异常值，构造某种曲面逼近似大地水准面，随着所构造的曲面不同，计算方法也不一样。其中，主要的方法有：平面拟合法、曲面拟合法、多面函数拟合法、样条函数法等。,9.5.3 GPS高程精度 影响GPS高程精度的主要有GPS大地高的精度、公共点几何水准的精度、GPS高程拟合的模型及方法、公共点的密度与分布等几个因素。具有高精度的GPS大地高是获得高精度GPS正常高的重要基础之一，因此必须采取措施以获得高精度的大地高，其中包括改善GPS星历的精度，提高GPS基线解算起算点坐标的精度，减弱对流层、电离层、多路径误差的影响等。,几何水准测量必须认真组织施测，保证提供具有足以满足精度要求
38、的相应等级的水准测量高程值。应根据不同测区，选用合适的拟合模型，以便使计算既准确又简便。均匀合理且足够地布设公共点，点位的分布和密度影响着GPS高程的精度。对于高差大于100m的测区，应加地形改正。对于大区域范围，可采用重力场模型加GPS水准的方法，拟合时对于不同趋势的区域，采用分区平差方法。理论分析和实践检验表明，在平原地区的局部GPS网，GPS水准可代替四等水准测量。在山区只要加地形改正，一般也可达到四等水准的精度。,介绍几种常用的GPS水准高程计算方法,1.绘等值线图法：,2.解析内插法：,3.曲面拟合,其步骤如下：,4.移动曲面拟合：,作业：教材第1、2、4、7、8、9、10、11、12题。,5.地形改正方法6多项式拟合精度,9.5.2 GPS重力高程,9.5.3 GPS三角高程,9.5.4 应用实例,布设原则,9.5.5 提高GPS水准精度的措施,1.提高大地高（差）测定的精度2.提高联测几何水准的精度3.提高转换参数的精度4.提高拟合计算的精度,9.6 精密基线解算软件简介,
