测量平差在测绘学科中的应用

测量平差在测绘学科中的应用测量平差与其他学科一样，是由于生产的需要而产生的，并在生产实践的过程中，随着科学技术的进步而发展。近代测量平差的内容非常丰富，其主要特点是，观测值概念广义化了，从处理随机独立的观测数据，展到可以处理随机相关的数据；扩展了经典测量平差的数学模型，从满秩平差问题，发展到降秩平差问题；从仅处理随机变量，发展到一并处理随机过程；从侧重于平差函数模型的研究，发展到也重视随机模型的研究；从不顾及模型误差，发展到顾及模型误差，针对最小二乘估计的局限性，提出了有偏估计和稳健估计。测量平差的基本任务是处理一系列带有偶然误差的观测值，求出未知量的最可靠值(平差值)，并评定测量成果的精度。测量平差中经典的估计准则是高斯创立的最小二乘估计准则。测量平差在进行数据处理时建立的函数模型一般都是确定的函数关系，即各种观测量之间都有明确的函数关系，例如：边长、角度与坐标之间的函数关系；水准网平差中的高程与高差之间的函数关系；GPS数据处理中的GPS卫星的伪距以及已知的卫星位置与接收机所在点的三个坐标之间，载波相位观测量以及已知的卫星位置与接收机所在点的三个坐标之间都是确定的函数关系；大地高、正常高与高程异常之间的函数关系式；卫星受摄动的轨道与六个轨道根数之间等等。1测量平差在变形监测中的应用在测量工作的实践和科学研究的活动中，变形观测占有重要的位置，而平差对于变形监测中的数据处理有着十分重要的作用。在工程建筑物的兴建中，从工程施工开始到竣工，以及建成后整个工程的运营期间都要不断的对工程建筑物进行监测，以便掌握工程建筑物变形的情况，及时发现问题，保证工程建筑物的安全，不论绝对网还是相对网，在观测期间网点位置均不能认为是没有变动的，即网中任意一点的稳定性必须进行检验。所谓对给定的控制网考察其可监测性，就是要预期该网可能监测到的最小变形量及方向。假定各观测点第一期真值为X1，第二期各点真值为X2，两期观测期间发生的位移量真值d̅，则：X2=X1+d̅(a)第一期的自由网平差的误差方程及基准条件方程:{V1=A1X1−l1GTX1=0误差方程解为：{1X=Ñ−1A1TP1l11XQ=Ñ−1−GGT第二期的自由网平差时，其误差方程与原来自由网平差第二期的误差方程相同，但其平差基准发生变化，是仍采用第一期的平差基准，以便保证基准一致性。现在将第二期自由网平差的误差方程及第一期的基准条件组合：{V2=A2X2−l2GTX1=0(b)将式(a)带入式(b)的第二式得：GTX2=GTd̅(c)式(c)就是用第二期近似高程或坐标值的改正值及真位移量表示的第一期基准。将(b)的第一式与(c)联合组成误差方程组得：{V2′=A2X2′−l2GTX2′=GTd̅(d)将(d)的第一式在最小二乘条件VTPV=min下求解得到法方程组{NX2′=A2TP2l2GTX2′=GTd̅(e)其中N不存在逆矩阵，将（e）第一式两边乘以G并加到第二式可以得到(N+GGT)X2′=ATP2l2+GGTd̅由于(N+GGT)可逆，解之。并代入以下式子{2X=Ñ−1A2TP2l22XQ=Ñ−1−GGT就有X2′=X2−Ñ−1GGTX2+Ñ−1GGTd̅(f)其中Ň=N+GGT，因为真位移量可以近似表示为：d̅=X2−X1(g)将(g)代入(f)就可以得出第二期观测数据在第一期基准下平差后的近似高程或坐标值的改正数，其值为X2′=X2—Ñ−1GGTX1，则同一基准下的位移量计算值为d̅‘=X2′−X1=d−Ñ−1GGTd̅，其中d̅‘=X2−X1。d是两期观测资料分别平差时的各点位移量，d̅‘就是基准一致性前提下，推导出的两期观测平差后各点的位移向量。在多期观测数据中如何合理地判断点的稳定性和计算位移量，这值得讨论。以往对多期观测数据的处理都是认为稳定点在不同观测期间将不发生变化，即网型不变，这只是一种理想化状态，但是实际中网型可能发生变化。如某期观测时部分稳定点被破坏，或者是对被破坏点重新埋设，此时网型都发生变化。平差时的基准也随之发生变化，已不是原来的基准。1.1监测网稳定性分析对于以上问题的解决我们可以设计如下一个观测网型。如在图1中，共有n个点，若作了m期观测，现在欲判断第i~j两期的发生位移点及位移量的大小。其中第j期观测时t号点被破坏，与t号点相关的几个观测量没有观测，此时网型发生变化。这就形成两期观测的基准不一致。同时对每个点的稳定程度也是未知的，即各个点稳定的权未知。图一水准网网型监测网稳定性分析思路为：(1)对多期观测数据作自由网整体平差，将各点在各期间视为互不相同的点，各观测周期数据看成相互独立。(2)计算各期的X和XQ。(3)对第i期的平差资料进行相似变换，解决网型不一致的情况。(4)计算位移量dij和协因数阵Qij。(5)计算i~j期间的合理参考系，并对位移量dij和协因数阵Qij再作相似变换，解决计算基准与实际基准不相符的情况。(6)再用平均间隙法作稳定点的判断。1.2监测网稳定性分析的基本理论（1）各期观测数据的自由网整体平差的误差方程组假定对图1的网形作了m期的沉降观测，在m期观测中网型可能发生变化，但观测精度相同。假定每期观测量为r个。11112212mmmnXAlAlXVAlX1122mmPPPP，1122mmAAAA其中，1X=(12...iinixxx)T是表示第i期点的近似高程改正值。il=(12...iimilll)T是第i期的观测数据，A11A22…Amm是各期的系数阵，权阵Pi只是第i期各观测值的权。（2）组成法方程并求解:将误差方程在最小二乘VTPV=min的条件下求解得出法方程:NX＝W(h)其中N＝ATPA，W＝ATPL1122mmNNNN，Nii=AiiTPiAii因自由网平差时N阵是奇异矩阵，不存在凯利逆，由法方程(h)求出的解不唯一。其解为:1XNW(i)协因数阵为1111221XmmNNQN，为求出方程(i)的唯一解，再给定一个最小范数条件minTXX，则可求出X的唯一解。1()TTTXAPAGGAPl，1()TTTXQAPAGGGG，其中，G满足条件0AG，TGGI；那么第i期的解为1()iTTTiiiiiiiiXQAPAGGGG1()TTTiiiiiiiiiiiiXAPAGGAPl。（3）对第i期的平差资料进行相似变换：因为前面己假设在第j期观测时第t号点被破坏。现在要将第i期平差值转换到第j期基准下，根据相似变换公式121(())TTXIHDHDX，2111(())(()TTTTXXQIHDHDQIHDHD有：iiTiiXSXSS后前后前XX，QQ，其中1TSIHHwHHw，11(111)TnHn，1101nnw。因为第j观测时t号点被破坏，欲将第i期观测的平差资料变换到第j期的基准下，这相似变换公式中稳定点权阵w应取j期的基准，其中第t项就应为0。（4）计算第i、j观测期间的位移量ijd及协因数阵ijQ。第i、j观测期间的位移量ijd及协因数阵ijQ计算公式如下：ijijjdXXdXXQQQ后i后，，在上式中，ijd表示的是第i、j两观测期间共同存在点的位移量，dQ是位移量ijd的协因数阵。（5）我们求尽量与实际相符合的参考系：与实际相符合的参考系可以表示为0,TTTDXDHw这是一个未知基准，TD参考系的系数阵，TH参考系中各点看成等权时的系数阵，w是参考系的权阵，其作用是对参考系中各点在平差中赋予不同的权重，为对角阵在这里是一个待求量。对于这里的高程网取为：1(1)1(1111)1TnHn，1()(1)()()iidkcwkcdkcdk，其中一般取为单位权方差，c为某一合适的常数。那么1(1)(1)(1)()TTdkIHHwkHHwKdk在上面的几个式中，k表示的是迭代次数。当(1)()iidkdk时停止迭代。一般取一个适当小的数。此时就求出了参考系各点的权阵w。计算出合乎实际的参考系权阵后，就可以对由自由网平差计算的位移和协因数阵作相似变换。变换公式如下：1(1)(1)TTSIHHwkHHwKDdSdDTddQSQSijdXXQQQ（6）利用平均间隙法判断各点的稳定性。利用相似变换后的位移量及协因数阵就可以判断点的稳定性。平均间隙法流程如下：图二平均间隙法计算流程对所有的参考点，都进行同样的分解，计算所有的iw、()iR，然后在所有的iw中取最大的一个，它相应的点为不稳定点，利用该点对应的()iR作图形一致性检验，若通过则终止，否则重复上述过程。2测量平差在GPS中的应用我们以GPS高程拟合的精度分析为例来谈其在GPS中的应用。2.1GPS高程方法在测量中常用的高程系统有以参考椭球面为基准面的大地高系统，一般用符号H表示；以大地水准面为基准面的正高系统，用符号gH表示；以似大地水准面为基准的正常高系统，用符号rH表示。高程系统间的相互关系如图所示：图三高程系统间的关系大地水准面到参考椭球面的距离，称为大地水准面差距，记为gh。大地高与正高之间的关系可以表示为：ggHHh似大地水准面到参考椭球面的距离，称为高程异常，记为。大地高与正常高之间的关系可以表示为：rHH由于采用GPS观测所得到的是点在WGS一84坐标系中的大地高，为了确定出正高或正常高，需要有大地水准面差距或高程异常数据。而我国常用的正常高(rH)则须有一定精度的高程异常值，才能保证由大地高求得。2.2高程拟合法高程拟合就是利用在范围不大的区域中，高程异常具有一定的几何相关性这一原理，由已知点的gH、rH在一定的数学模型和统计准则下求出未知点的高程异常，从而求出待定点的正常高。若要用零次多项式进行高程拟合时，要确定1个参数，因此，需要l个以上的已知点；若要采用一次多项式进行高程拟合，要确定3个参数，需要3个以丘的已知点；若要采用二次多项式进行高程拟合，要确定6个参数，则需要6个以上的己知点。将高程异常表示为下面多项式的形式：零次多项式：0a一次多项式：012aadBadL二次多项式：22012345aadBadLadBadLadBdL其中：0dBBB0dLLL01BBn01LLn（n为GPS的点数）利用公共点上GPS测定的大地高和水准测量测定的正常高计算出该点上的高程异常，存在一个这样的公共点，就可以依据上式列出一个方程：22012345iiiiiiiaadBadLadBadLadBdL若共存在m个这样的公共点，则可列出m个方程：221011213141511aadBadLadBadLadBdL222012223242522aadBadLadBadLadBdL……22012345mmmmmmmaadBadLadBadLadBdL即有：VAxL其中：22111111222222222111mmmmmmdBdLdBdLdBdLdBdLdBdLdBdLAdBdLdBdLdBdL012345Txaaaaaa12TmV通过最小二乘法可以求解出多项式的系数：1()()TTxAPAAPL其中：P为权阵，它可以根据水准高程和GPS所测得的大地高的精度来加以确定。按上述方法便可以确定计算点的高程，其精度主要决定于GPS测量的精度。这一方法的优点是概念明了，计算简单，精度高。不过，为描述大地水准面的细节，它需要布设均匀的、密度充分的GPS观测点，并且在这些点上，需要同时具