GPS控制网数据预处理无约束平差及其成果转换的探讨.docVIP

GPS控制网数据预处理无约束平差及其成果转换的探讨 　　中图分类号：P228.4 　　摘要：本文从无约束平差的数学模型、高斯投影转换、确定地方椭球、转换到地方独立坐标系及精度评定几个方面探讨了无约束平差及其成果转换，具有较高的价值和系统性，供参考。 　　关键词：无约束平差；GPS；坐标系 　　 　　GPS基线向量观测值是WGS-84坐标系中的三维空间直角坐标差，通过对GPS控制网进行同步环、异步环闭合差及重复观测基线较差的比较，检验闭合差及较差是否均在规定限范围之内这几步之后，三维或二维GPS网无约束平差所得的单位权、长度、方位的精度真反映了GPS观测所具有的随机误差统计结果，因而可以衡量CPS测量成果的精度。在GPS网维或二维无约束平差后，如果平差后单位权中误差，基线坐标分量的改正数在毫米或最大在2-3cm以内，一般可认为GPS网无明显粗差;所以，GPS网的无约束平差，相当于统测量中对观测数据的检核。 　　对不含地面数据的GPS网在WGS-84坐标系中进行三维平差有如下作用： 　　1.检查GPS基线向量是否有粗差和明显的系统误差，并考察GPS网的内部精度和 　　基线向量的观测精度。 　　2.提供GPS网中各点相对于某点的大地高差数据，以便进一步求定各点的大地高 　　用于进行高程转换。 　　3.提供网中各点WGS-84坐标，用于和地面己知坐标数据进行坐标转换。 　　1无约束平差的数学模型 　　GPS网三维无约束平差数学模型： 　　 　　 　　 　　如果把公式(1)式写为矩阵形式，则 　　 (2) 　　其中，系数矩阵A为 　　 　　 　　 　　因为GPS网是三维控制网，所以，在对GPS网进行平差时，应有三个位置基准，一个尺度基准和三个方位基准。而GPS基线向量包含有WGS一84坐标系中的尺度和方位信息，只缺位置基准，所以一般可以取某点单点定位的三维坐标或者某个国家点的坐标作为位置基准。 　　若选择第K点为基准点，则基准条件为 　　 　　 　　 　　若取基线解算软件给出的观测值的方差阵，则，综合公式(2)和(3)可得到待估参数的计算公式为， 　　(4) 　　而，所以的协因数阵为 　　(5) 　　也可以按秩亏平差法求解，即以网的重心为位置基准，仍以GPS基线向量观测值的尺度和方位为尺度基准和方位基准，将位置基准方程写为 　　（6） 　　其中G应满足 　　 (7) 　　通常取为 　　 　　 　　 　　按秩亏平差原理可得 　　 (8) 　　因为GPS单点定位的精度不高。因此，经GPS网平差后，各点在WGS-84坐标系中的坐标值精度较低。但他们相对于网的位置基准的精度很高。从这个意义上说，GPS网平差后各点的坐标并不属于真正的WGS一84坐标系，而是属于一种独立坐标系。 　　如果需要得到GPS点的大地坐标可按下式转换得到 　　 　　 　　 　　其中大地纬度B由第二式按迭代法求得迭代时B的初值。 　　2高斯投影转换 　　GPS测量得到的是空间三维直角坐标，需经过坐标变换、高斯投影后刁‘能得到所需的参考椭球面上的高斯平面直角坐标。经过了高斯投影，它必然产生两种变形，即高程归算变形和高斯投影变形。将GPS地面观测的长度归算到参考椭球面上产生的变形称高程归算变形，可按式(10)计算： 　　(10) 　　式中观测边的平均大地高。 　　―长度所在方向参考椭球面法截弧的曲率半径; 　　D―实测水平距离 　　再将参考椭球面上的长度投影到高斯平面上产生的变形称高斯投影变形，可按式(11)计算 　　(11) 　　式中―测线在高斯平面上离中央子午线垂距的平均值 　　一测线两端平均纬度处参考椭球面的平均曲率半径 　　S―参考椭球面上的长度 　　这样，地面上一段距离经过两次改正，被改变了真实长度，这种高斯投影面上的长度与地面真实长度之差称为长度综合变形，可按式(12)计算： 　　(12) 　　为实际计算方便，又不至损害必要的精度，取，取不同投影面上的同一距离近似相等，即，写成相对变形的形式为： 　　 （13） 　　公式(13)表明：采用国家统一坐标系统所产生的长度综合变形与测区所处投影带内的位置和测区平均高程有关，利用式(3.13)可以方便地计算出己知测区内长度相对变形的大小。 　　1.不同投影带的出现是因为选择了不同经度的中央子午线的缘故，如果我们合理选择中央子午线的位置，使长度投影到该投影带所产生的变形恰好抵偿这一长度投影到椭球面上产生的变形，则高斯平面上的长度也能够和实际长度保持一致，避免长度变形。我们称这种能够抵偿长度变形的投影带为“抵偿投影带”。令为零，可以计算出新投影带中央子午线离开测区中央的距离： 　　 （14）