测绘数据处理自由网平差课件

测绘数据处理

2023/3/141一、自由网平差概述在控制网的经典间接平差中，必须具有足够的起算数据。例如，在水准网的间接平差中，必须至少已知某一点的高程；在测角网的间接平差中，必须至少已知某一点的坐标、某一条边的坐标方位角即某一条边的边长，等等。下面将讨论无起算数据的平差方法，即自由网平差。当网中有足够的起算数据时，经典间接平差的误差方程为

(1-7-1）系数矩阵B最大线性无关的行（列）向量的个数，及B矩阵的秩R(B)等于未知参数的个数t.即

(1-7-2)2023/3/142在最小二乘准则下，得其法方程为

(1-7-3)其中N=PB,W=。此时，系数阵N为满秩方阵，即det(N),N为非奇异阵，有唯一解，其解为

(1-7-4)当平差网没有起算数据时，网中所有的点均为待定点。设未知参数的个数为u,误差方程为

(1-7-5)组成的法方程为2023/3/143d就是网中必要的起算数据个数。且有：

二、秩亏自由网平差思路为了求得未知参数的唯一确定解，除了遵循最小二乘准则外，还必须增加新的约束条件，从而达到求得唯一解的目的。由于约束条件不同，秩亏自由网平差可分为如下几种情况：（1）、经典自由网平差。它是在假设网中有d个必要起算数据的条件下，求定未知参数的最佳估计。这种方法早就已为人们所熟知。不难理解，该法的平差结果（未知参数X的解及其协因数阵）将随着假设的d个必要起算数据的不同而不同，即随着已知点位置的改变而改变。2023/3/145（2）、秩亏网平差。它是在最小二乘和最小范数的条件下求定未知参数的最佳估值。（3）、加权秩亏网平差。它是在最小二乘和加权最小范数的条件下求定未知参数的最佳估值。式中，为表示未知参数稳定程度的先验权矩阵。（4）、拟稳平差。若将平差网中的未知参数分为两类，即

(s>d)(1-7-7)

式中，是非拟稳点的未知参数，是拟稳点的未知参数。这样拟稳平差是在和求定未知参数的最佳估值。2023/3/146由上可知，三种秩亏自由网平差均遵循的原则，对于同一平差问题，它们将有相同的法方程，三种自由网平差的解均能满足法方程式(1-7-6)，它它们都是这一相同法方程多组解中的一个特解。它们之间的不同只是由于各自对解向量x所加的限制条件不同引起的，即由于各自所加的最小范数条件不同，因此得到了不同的解向量。由于秩亏网平差与拟稳平差都是加权秩亏网平差的特殊情况，其区别仅在于各自选择了不同的先验权阵。所以我们将先介绍加权秩亏网平差，然后再介绍秩亏网平差和拟稳平差。自由网平差的方法很多，本节我们只介绍附加条件法。2023/3/147对三维网，这种约束将要求：.平差前后网点重心不变（即固定重心点坐标）；.平差前后控制网相对其重心不绕X、Y、Z轴旋转（即固定重心与所有网点连线的平均方位角与天顶距）；.平差前后所有网点相对重心的平均距离不变（即固定与重心的平均距离）。对平面网，这种约束将要求：.平差前后网点重心坐标不变；.平差前后各网点与重心连心的平均方位角不变；.平差前后所有网点相对重心的平均距离不变（即固定与重心的平均距离）。对于一维的高程网，这种约束是使平差前后网店的平均高程保持不变。这些约束条件我们称之为重心基准条件。2023/3/149（三）加权秩亏自由网平差基准和秩亏自由网平差基准类似，但应考虑各网点的权重，采用了带权重心基准条件。（四）拟稳平差基准也和秩亏自由网平差基准类似，但仅仅是采用所有拟稳点的重心基准条件。

四、加权秩亏自由网平差（一）平差计算公式加权秩亏网平差就是求相容法方程：

(1-7-8)的加权最小范数解。

附加条件法的基本思想是：由于网中没有起始数据，平差时多选了d个未知参数。现在u个未知参数之间适当给定d个附加条件2023/3/1410式，即在平差问题的函数模型中加入d个未知参数的限制条件方程，从而可以按附有限制条件的间接平差法求解。等价于约束条件的限制条件方程为

(1-7-9)式中BG=0(1-7-10)故加权秩亏网平差函数模型为

(1-7-11)2023/3/1411

(1-7-15)可以证明（证明略），当G满足条件BG=0时，连系数向量K必等于零。故可简化为

(1-7-16)将代入式(1-7-11)，可求得V，再根据即可求得个未知参数的平差值需要说明的是，在实际计算时，附加阵G不仅要满足BG=0的条件，还要满足条件

(1-7-17)也即在实际计算前，尚需要把G阵进行标准化，是的满足式(1-7-17)所述条件。2023/3/1413（二）精度评定公式1、验后单位权方差

(1-7-18)式中，可以直接计算，也可用计算得到2、协因数阵（1）未知参数的协因数阵为

(1-7-19)当对G阵进行标准化后，由于，故式(1-7-19)可进一步简化为(1-7-20)

（2）观测数据平差值的协因数阵为

(1-7-21)2023/3/1414

（三）常见附加阵G采用附加条件法进行秩亏自由网平差计算非常容易，但必须预先构建附加阵G。下面是几种常见的附加阵G。1、水准网设有u个点，则

(1-7-22)

2、平面测角网按角度平差，设有m个点，则

(1-7-23)2023/3/1415对于平面测边网、边角网和导线网，只要将式(1-7-23)或式(1-7-24)中的第四行划去，剩下的三行u列的阵，即分别为按角度平差时的附加阵。对于三维测边网，只要将式（1-7-25）的第七行划去，剩下的6三行u列的阵即为三维测边网平差时的附加阵。很明显，上述的附加阵G均未标准化，即只是满足了BG=0，但尚未满足的条件。2023/3/1417

阵标准化1、用原始阵和阵，求出相应的阵；2、设中第i行主对角元素为gii，把原始阵相应的第i行数据均乘以即可得到标准化阵的相应数据；3、原始阵中每一行数据均按（2）所述做同样变换，即可得到标准化阵。例：有测边网如下图所示，各点的近似坐标列于下表。若又已知待定点的先验权阵为2023/3/1418试求加权秩亏网平差的标准化阵。表1-7-1解：计算秩亏网平差时的标准化阵，即计算满足BG=0和计算网的加权重心点坐标。

点号43212023/3/1419（4）求解

（5）标准化阵确定把原始阵中的第一行、第二行、第三行分别乘以，即可求得标准化阵为

2023/3/1421

五、秩亏自由网平差秩亏网是在最小二乘和最小范数的条件下求定未知参数的最佳估值。也可叙述为秩亏网平差是求相容法方程在最小范数

(1-7-26)条件下的解。它是加权秩亏网平差时的特例。因此，秩亏网平差的各种计算公式均可由加权秩亏网平差时的计算公式列出。（一）平差计算公式等价于约束条件的限制条件方程为2023/3/1422（1-7-27）式中

BG=0（1-7-28）故秩亏网平差的函数模型为

（1-7-29）此处的系数矩阵B不是列满矩阵，而是列亏矩阵。将式（1-7-29）组成法方程，得

(1-7-30)式中，2023/3/1423（二）精度评定公式（）1、验后单位权方差

(1-7-36)式中，可以直接计算，也可用计算得到2、协因数阵（1）未知参数的协因数阵为

(1-7-37)当对G阵进行标准化后，由于，故式(1-7-19)可进一步简化为(1-7-38)

（2）观测数据平差值的协因数阵为

(1-7-39)2023/3/1425

（三）阵标准化1、用原始阵和阵，求出相应的阵；2、设中第i行主对角元素为gii，把原始阵相应的第i行数据均乘以即可得到标准化阵的相应数据；3、原始阵中每一行数据均按（2）所述做同样变换，即可得到标准化阵。例：测角网如下图，全网4个待定点的坐标值列于下表.当采用秩亏网测角平差时，求其标准化阵。2023/3/1426(3)原始阵确定。（4）求解

（5）标准化阵确定把原始阵中的第一行、第二行、第三行、第四行分别乘以

2023/3/1429六、拟稳平差拟稳平差是在最小二乘和最小范数（局部解向量的范数最小）的条件下，求定位置参数的最佳估值。也可叙述为：拟稳平差是相容法方程，（1-7-40）

在最小范数条件下的解。可见，它是加权秩亏网平差取时的特例。因此，拟稳平差的各种计算公式也可由加权秩亏网平差的计算公式直接引出。2023/3/1430

（一）平差计算公式

在式（1-7-40）第二式中，取（1-7-41）并设（1-7-42）则（1-7-43）若令（1-7-44）则可得拟稳平差的函数模型为（1-7-45）此处的系数矩阵B不是列满矩阵，而是列亏矩阵。将式（1-7-45）组成法方程，得

2023/3/1431

（1-7-46）式中，，。因为降秩方阵，无正常逆，所以必须对法方程做适当变动。将（1-7-46）中的第二个方程左乘后再加到第一个方程上去，既得变形后的法方程为（1-7-47）式中，，。解法方程，得（1-7-48）（1-7-49）可以证明,当满足条件时，联系数向量K必等于零。故式可简化为（1-7-50）代入式（1-7-45），可求得V，即可求得各未知参数的平差值。2023/3/1432在实际计算时，还要对附加阵提出进一步的要求（1-7-51）在实际计算前，也需要把阵进行标准化，使得满足式（1-7-51）所述条件。（二）精度评定公式

1.验后单位权方差（1-7-52）式中，可以直接计算，也可以计算得到

2.协因数阵（1）未知参数的协因数阵为（1-7-53）当对阵进行标准化后，由于，可证明,故式（1-7-53）可进一步简化为（1-7-54）2023/3/1433（2）观测数据平差值的协因数阵为