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

1、测绘数据处理测绘数据处理 测绘工程教研室测绘工程教研室土地科学技术学院土地科学技术学院一、自由网平差概述一、自由网平差概述 在控制网的经典间接平差中，必须具有在控制网的经典间接平差中，必须具有足够的起算数据。例如，在水准网的间接平差足够的起算数据。例如，在水准网的间接平差中，必须至少已知某一点的高程；在测角网的中，必须至少已知某一点的高程；在测角网的间接平差中，必须至少已知某一点的坐标、某间接平差中，必须至少已知某一点的坐标、某一条边的坐标方位角即某一条边的边长，等等一条边的坐标方位角即某一条边的边长，等等。下面将讨论无起算数据的平差方法，即自由。下面将讨论无起算数据的平差方法，即自由网平差。

2、网平差。当网中有足够的起算数据时，经典间接平差当网中有足够的起算数据时，经典间接平差的误差方程为的误差方程为 (1-7-1）系数矩阵系数矩阵B最大线性无关的行列向量的个最大线性无关的行列向量的个数，及数，及B矩阵的秩矩阵的秩R(B)等于未知参数等于未知参数 的个数的个数t.即即 (1-7-2)在最小二乘准则下，得其法方程为在最小二乘准则下，得其法方程为 (1-7-3)其中其中N= PB,W= 。此时，系数阵。此时，系数阵N为满秩方阵，即为满秩方阵，即det(N) ,N为非奇异阵，有唯一解，其解为为非奇异阵，有唯一解，其解为 (1-7-4)当平差网没有起算数据时，网中所有的点均为待定点。设未知当

3、平差网没有起算数据时，网中所有的点均为待定点。设未知参数的个数为参数的个数为u,误差方程为误差方程为 (1-7-5)组成的法方程为组成的法方程为 (1-7-6) 由于det(N)=0,故N为奇异阵，其凯利逆 不存在，此时如仍按经典平差公式接，将不可能得出唯一的解。 令B的列满秩数为 （B）,B的实际秩数为R(B)，d= (B)-R(B),d即为秩亏数。对于法方程系数矩阵N，必然也有d=(N)-R(N). 如果d=0,就是经典平差问题；当d 时，就是所谓的秩亏自由网平差问题。 在实用上，产生秩亏得主要原因是不设起算数据，而且选定网中高程、坐标等作为平差的未知参数，所以秩亏自由网平差也叫无固定数据

4、的的自由网平差，简称自由网平差。 具体问题中矩阵B或N的秩亏数d，虽然可以通过计算R(B)或R(N)求出，实际上并不需要这样做，究其原因可知，秩亏数d就是网中必要的起算数据个数。且有：就是网中必要的起算数据个数。且有： 二、秩亏自由网平差思路二、秩亏自由网平差思路 为了求得未知参数的唯一确定解，除了遵循最小二乘准则外为了求得未知参数的唯一确定解，除了遵循最小二乘准则外，还必须增加新的约束条件，从而达到求得唯一解的目的，还必须增加新的约束条件，从而达到求得唯一解的目的。由于约束条件不同，秩亏自由网平差可分为如下几种情。由于约束条件不同，秩亏自由网平差可分为如下几种情况：况：（1）、经典自由网平差

5、。它是在假设网中有）、经典自由网平差。它是在假设网中有d个必要起算数据个必要起算数据的条件下，求定未知参数的最佳估计。这种方法早就已为的条件下，求定未知参数的最佳估计。这种方法早就已为人们所熟知。不难理解，该法的平差结果未知参数人们所熟知。不难理解，该法的平差结果未知参数X的解的解及其协因数阵及其协因数阵 ）将随着假设的）将随着假设的d个必要起算数据的不同而个必要起算数据的不同而不同，即随着已知点位置的改变而改变。不同，即随着已知点位置的改变而改变。（2）、秩亏网平差。它是在最小二乘）、秩亏网平差。它是在最小二乘 和最小范数和最小范数的条件的条件 下求定未知参数的最佳估值。下求定未知参数的最佳

6、估值。（3）、加权秩亏网平差。它是在最小二乘）、加权秩亏网平差。它是在最小二乘 和加权最和加权最小范数的条件小范数的条件 下求定未知参数的最佳估值。式下求定未知参数的最佳估值。式中，中， 为表示未知参数稳定程度的先验权矩阵。为表示未知参数稳定程度的先验权矩阵。（4）、拟稳平差。若将平差网中的未知参数分为两类，即）、拟稳平差。若将平差网中的未知参数分为两类，即 (sd) (1-7-7) 式中，式中， 是非拟稳点的未知参数，是非拟稳点的未知参数， 是拟稳点的未知参数。这样是拟稳点的未知参数。这样拟稳平差是在拟稳平差是在 和和 求定未知参数的最佳估求定未知参数的最佳估值。值。 由上可知，三种秩亏自由

7、网平差均遵循由上可知，三种秩亏自由网平差均遵循 的原则，对的原则，对于同一平差问题，它们将有相同的法方程，三种自由网平差的解于同一平差问题，它们将有相同的法方程，三种自由网平差的解均能满足法方程式均能满足法方程式(1-7-6)，它它们都是这一相同法方程多组解，它它们都是这一相同法方程多组解中的一个特解。它们之间的不同只是由于各自对解向量中的一个特解。它们之间的不同只是由于各自对解向量x所加的所加的限制条件不同引起的，即由于各自所加的最小范数条件不同，因限制条件不同引起的，即由于各自所加的最小范数条件不同，因此得到了不同的解向量。此得到了不同的解向量。 由于秩亏网平差与拟稳平差都是加权秩亏网平差

8、的特殊情况，由于秩亏网平差与拟稳平差都是加权秩亏网平差的特殊情况，其区别仅在于各自选择了不同的先验权阵其区别仅在于各自选择了不同的先验权阵 。所以我们将先介。所以我们将先介绍加权秩亏网平差，然后再介绍秩亏网平差和拟稳平差。自由网绍加权秩亏网平差，然后再介绍秩亏网平差和拟稳平差。自由网平差的方法很多，本节我们只介绍附加条件法。平差的方法很多，本节我们只介绍附加条件法。 三、自由网平差基准三、自由网平差基准（一经典自由网平差基准（一经典自由网平差基准 实际上是固定某个实际的点的位置，固定一个实际的方位角，实际上是固定某个实际的点的位置，固定一个实际的方位角，一条实际的边长等来定义的。一条实际的边长

9、等来定义的。（二秩亏自由网平差基准（二秩亏自由网平差基准 与经典的自由网不同，秩亏自由网平差基准是通过对整个网点与经典的自由网不同，秩亏自由网平差基准是通过对整个网点的坐标或部分网点的坐标进行某种约束条件来定义的，这种的坐标或部分网点的坐标进行某种约束条件来定义的，这种约束实际上是固定某个虚拟点的位置，固定某个虚拟方向、虚拟约束实际上是固定某个虚拟点的位置，固定某个虚拟方向、虚拟距离等。距离等。 对三维网，这种约束将要求：对三维网，这种约束将要求：.平差前后网点重心不变即固定重心点坐标）；平差前后网点重心不变即固定重心点坐标）；.平差前后控制网相对其重心不绕平差前后控制网相对其重心不绕X、Y、

10、Z轴旋转即固定重心轴旋转即固定重心 与所有网点连线的平均方位角与天顶距）；与所有网点连线的平均方位角与天顶距）；.平差前后所有网点相对重心的平均距离不变即固定与重心的平差前后所有网点相对重心的平均距离不变即固定与重心的 平均距离）。平均距离）。对平面网，这种约束将要求：对平面网，这种约束将要求：.平差前后网点重心坐标不变；平差前后网点重心坐标不变；.平差前后各网点与重心连心的平均方位角不变；平差前后各网点与重心连心的平均方位角不变；.平差前后所有网点相对重心的平均距离不变即固定与重心的平差前后所有网点相对重心的平均距离不变即固定与重心的 平均距离）。平均距离）。对于一维的高程网，这种约束是使平

11、差前后网店的平均高程保持对于一维的高程网，这种约束是使平差前后网店的平均高程保持 不变。不变。这些约束条件我们称之为重心基准条件。这些约束条件我们称之为重心基准条件。（三加权秩亏自由网平差基准（三加权秩亏自由网平差基准和秩亏自由网平差基准类似，但应考虑各网点的权重，采用了带和秩亏自由网平差基准类似，但应考虑各网点的权重，采用了带权重心基准条件。权重心基准条件。（四拟稳平差基准（四拟稳平差基准也和秩亏自由网平差基准类似，但仅仅是采用所有拟稳点的重心也和秩亏自由网平差基准类似，但仅仅是采用所有拟稳点的重心基准条件。基准条件。 四、加权秩亏自由网平差四、加权秩亏自由网平差（一平差计算公式（一平差计算

12、公式 加权秩亏网平差就是求相容法方程：加权秩亏网平差就是求相容法方程： (1-7-8) 的加权最小范数解。的加权最小范数解。 附加条件法的基本思想是：由于网中没有起始数据，平差时多附加条件法的基本思想是：由于网中没有起始数据，平差时多选了选了d个未知参数。现在个未知参数。现在u个未知参数之间适当给定个未知参数之间适当给定d个附加条件个附加条件式，即在平差问题的函数模型中加入式，即在平差问题的函数模型中加入d个未知参数的限制条件个未知参数的限制条件方程，从而可以按附有限制条件的间接平差法求解。方程，从而可以按附有限制条件的间接平差法求解。 等价于约束条件等价于约束条件 的限制条件方程为的限制条件

13、方程为 (1-7-9)式中式中 BG=0 (1-7-10)故加权秩亏网平差函数模型为故加权秩亏网平差函数模型为 (1-7-11)此处的系数矩阵此处的系数矩阵B不是列满矩阵，而是列亏矩阵。不是列满矩阵，而是列亏矩阵。将式将式1-7-11组成法方程，得组成法方程，得 (1-7-12)式中式中 ， 因因N为降秩方阵，无正常逆，所以为降秩方阵，无正常逆，所以必须对法方程作适当变动。将式必须对法方程作适当变动。将式(1-7-12)中第二个方程左乘后中第二个方程左乘后再加到第一个方程上去，即得变形后的法方程为再加到第一个方程上去，即得变形后的法方程为 (1-7-13)式中，式中，解法方程，得解法方程，得

14、(1-7-14) (1-7-15)可以证明证明略），当G满足条件BG=0时，连系数向量K必等于零。故可简化为 (1-7-16)将代入式(1-7-11)，可求得V，再根据 即可求得个未知参数的平差值 需要说明的是，在实际计算时，附加阵G不仅要满足BG=0的条件，还要满足条件 (1-7-17)也即在实际计算前，尚需要把G阵进行标准化，是的满足式(1-7-17)所述条件。（二精度评定公式（二精度评定公式1、验后单位权方差、验后单位权方差 (1-7-18)式中，式中， 可以直接计算，也可用可以直接计算，也可用 计算得到计算得到2、协因数阵、协因数阵（1未知参数的协因数阵为未知参数的协因数阵为 (1-7

15、-19)当对当对G阵进行标准化后，由于阵进行标准化后，由于 ，故，故 式式(1-7-19)可进一步简化为可进一步简化为 (1-7-20) （2观测数据平差值的协因数阵为观测数据平差值的协因数阵为 (1-7-21) （三常见附加阵（三常见附加阵G采用附加条件法进行秩亏自由网平差计算非常容易，但必须预先采用附加条件法进行秩亏自由网平差计算非常容易，但必须预先构建附加阵构建附加阵G。下面是几种常见的附加阵。下面是几种常见的附加阵G。1、水准网、水准网设有设有u个点，那么个点，那么 (1-7-22) 2、平面测角网、平面测角网按角度平差，设有按角度平差，设有m个点，那么个点，那么 (1-7-23)3、

16、平面测角网、平面测角网按方向平差，设有按方向平差，设有m个点，那么个点，那么 （1-7-24）4、三维测角网、三维测角网 （1-7-25） 对于平面测边网、边角网和导线网，只要将式对于平面测边网、边角网和导线网，只要将式(1-7-23)或式或式(1-7-24)中的第四行划去，剩下的三行中的第四行划去，剩下的三行u列的阵，即分别为按列的阵，即分别为按角度平差时的附加阵。对于三维测边网，只要将式角度平差时的附加阵。对于三维测边网，只要将式1-7-25的第七行划去，剩下的的第七行划去，剩下的6三行三行u列的阵即为三维测边网平差时的列的阵即为三维测边网平差时的附加阵。附加阵。 很明显，上述的附加阵很明

17、显，上述的附加阵G均未标准化，即只是满足了均未标准化，即只是满足了BG=0，但尚未满足的条件。但尚未满足的条件。 阵标准化阵标准化1、用原始阵、用原始阵 和和 阵，求出相应的阵阵，求出相应的阵 ；2、设、设 中第中第i行主对角元素为行主对角元素为gii，把原始阵，把原始阵 相应相应的第的第i行数据均乘以行数据均乘以 即可得到标准化阵的相应数据；即可得到标准化阵的相应数据；3、原始阵中每一行数据均按、原始阵中每一行数据均按2所述做同样变换，即可得所述做同样变换，即可得到标准化阵。到标准化阵。例：有测边网如下图所示，各点的近似坐标列于下表。若又例：有测边网如下图所示，各点的近似坐标列于下表。若又已

18、知待定点的先验权阵为已知待定点的先验权阵为试求加权秩亏网平差的标准化阵。试求加权秩亏网平差的标准化阵。 表表1-7-1解：解： 计算秩亏网平差时的标准化阵，即计算满足计算秩亏网平差时的标准化阵，即计算满足BG=0和和 计算网的加权重心点坐标。计算网的加权重心点坐标。 点号 4 3 2 1（2 2计算以加权重心点坐标为坐标原点的各待定点的坐标值列于表计算以加权重心点坐标为坐标原点的各待定点的坐标值列于表1-7-2.1-7-2.表表1-7-21-7-2点号 /km4321(3)原始阵 确定。由于是测角网，根据式(1-7-23),即可得到测边网原始阵 （按角度平差）（4求解 （5标准化阵确定把原始阵

19、中的第一行、第二行、第三行分别乘以，即可求得标准化阵为 五、秩亏自由网平差五、秩亏自由网平差 秩亏网是在最小二乘和最小范数的条件下求定未知参数的最佳秩亏网是在最小二乘和最小范数的条件下求定未知参数的最佳估值。也可叙述为秩亏网平差是求相容法方程在最小范数估值。也可叙述为秩亏网平差是求相容法方程在最小范数 (1-7-26) 条件下的解。它是加权秩亏网平差条件下的解。它是加权秩亏网平差 时的特时的特例。因而，秩亏网平差的各种计算公式均可由加权秩亏网例。因而，秩亏网平差的各种计算公式均可由加权秩亏网平差时的计算公式列出。平差时的计算公式列出。（一平差计算公式（一平差计算公式等价于约束条件等价于约束条件

20、 的限制条件方程为的限制条件方程为 （1-7-27）式中 BG=0 （1-7-28）故秩亏网平差的函数模型为 （1-7-29）此处的系数矩阵B不是列满矩阵，而是列亏矩阵。将式（1-7-29组成法方程，得 (1-7-30)式中， 因N为降秩方阵，无正常逆，所以必须对法方程作适当变动。将式(1-7-30)中第二个方程左乘后再加到第一个方程上去，记得变形后的法方程为 (1-7-31) ，式中，解法方程，得 (1-7-32) (1-7-33)可以证明证明略），当G满足条件BG=0时，连系数向量K必等于零。故可简化为 (1-7-34)将代入式(1-7-29)，可求得V，再根据 即可求得个未知参数的平差值

21、。需要说明的是，在实际计算时，附加阵G不仅要满足BG=0的条件，还要满足条件 (1-7-35)也即在实际计算前，尚需要把G阵进行标准化，满足式(1-7-35)所述条件。（二精度评定公式（二精度评定公式（ ）1、验后单位权方差、验后单位权方差 (1-7-36)式中，式中， 可以直接计算，也可用可以直接计算，也可用 计算得到计算得到2、协因数阵、协因数阵（1未知参数的协因数阵为未知参数的协因数阵为 (1-7-37)当对当对G阵进行标准化后，由于阵进行标准化后，由于 ，故，故 式式(1-7-19)可进一步简化为可进一步简化为 (1-7-38) （2观测数据平差值的协因数阵为观测数据平差值的协因数阵为

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

23、P2P3P4（1计算网的重心点坐标 （2计算以加权重心点坐标为坐标原点的各待定点的坐标值 点号 /km P1 P2 P3 P4 (3)原始 阵确定。（4求解 （5标准化阵确定把原始阵中的第一行、第二行、第三行、第四行分别乘以 六、拟稳平差六、拟稳平差 拟稳平差是在最小二乘和最小范数局拟稳平差是在最小二乘和最小范数局部解向量的范数最小）部解向量的范数最小） 的条件下，的条件下，求定位置参数的最佳估值。也可叙述为：拟求定位置参数的最佳估值。也可叙述为：拟稳平差是相容法方程稳平差是相容法方程 ， （1-7-40） 在最小范数条件下的解。可见，它是加在最小范数条件下的解。可见，它是加权秩亏网平差取权秩

24、亏网平差取 时的特例。因而，拟稳平差的各种计算公式时的特例。因而，拟稳平差的各种计算公式也可由加权秩亏网平差的计算公式直接引出也可由加权秩亏网平差的计算公式直接引出。 （一平差计算公式 在式1-7-40第二式中，取 （1-7-41） 并设 （1-7-42）那么 （1-7-43）若令 （1-7-44）则可得拟稳平差的函数模型为 （1-7-45）此处的系数矩阵B不是列满矩阵，而是列亏矩阵。将式1-7-45组成法方程，得 （1-7-46）式中， ，。 因为降秩方阵，无正常逆，所以必须对法方程做适当变动。将1-7-46中的第二个方程左乘 后再加到第一个方程上去，既得变形后的法方程为 （ 1-7-47） 式中，。 解法方程，得