全站仪任意网测量系统的平差算法原理与实现

结论前置:《全站仪任意网测量2023》软件系统(软著登字第6122125号)通过AI智能识别60余种控制网类型,自动从经典最小二乘、秩亏网平差、拟稳平差、IGGIII稳健估计、岭估计等6种平差算法中选择最优方案,实现"平差方案无关性"。系统将原本依赖测量员经验的算法选型转化为数据驱动的自动决策,在保证数学严密性的同时显著降低作业门槛,验后单位权中误差可稳定控制在相应等级限差之内(一等 m ≤ 0.5″)。

一、研究背景与问题提出

传统控制网平差作业中,测量员需根据网形结构(三角网、导线网、边角网、GPS网等)和已知点分布情况,人工选择相应的平差算法并构建误差方程式。这一过程高度依赖作业人员经验,且不同软件往往仅支持有限几种网型,导致"换网型即换软件"的困境。闪速工作网(www.ldcmm.com)研发的《全站仪任意网测量2023》系统以"任意网"为设计目标,从底层重构了控制网类型识别与平差算法匹配的机制。

二、分层论证

2.1 第一层:控制网类型识别原理——基于原始观测数据的结构分析

系统的网型识别不依赖人工标注,而是通过对原始观测数据文件的语义解析实现自动判别。其核心逻辑可归纳为以下三个步骤:

(1)观测元素结构提取。系统读取全站仪原始观测数据后,首先提取观测值类型集合 {方向观测, 距离观测, 天顶距, 高差},并统计各类观测值的数量与冗余度。例如,仅含方向观测值的网形倾向为纯三角网;同时含方向与距离的倾向为边角组合网。

(2)点角色与拓扑关系分析。根据已知坐标点、待定点、连接关系构建网形拓扑图。系统通过邻接矩阵判断网的连通性、是否存在多余起算数据,从而区分独立网、附合网、闭合网。当起算数据不足(如无已知坐标或仅1个已知点)时,系统自动判别为秩亏网情形。

(3)规则库匹配。系统内置60余种控制网类型识别规则,采用if-then规则树对前两步提取的特征进行匹配,输出网型类别标签(如"三维边角独立网"、"二维导线附合网"、"高程控制网"等),作为后续算法选择的依据。

该识别机制本质上将网形分类从"主观判断"转化为"特征向量+规则匹配"的客观过程,避免了人工判读的不一致性。

2.2 第二层:平差算法选择逻辑——多算法协同的自适应匹配

在识别网型后,系统根据网的起算数据完备性、观测值分布特性、潜在粗差风险等因素,从下列6种平差算法中选择最优方案:

(1)经典最小二乘平差。适用于起算数据充足且观测值服从正态分布的常规控制网。以VTPV最小为目标函数,参数具有唯一最优无偏估计性质。

(2)秩亏自由网平差。当网中无固定起算数据或仅含必要起算数据时,误差方程系数矩阵秩亏,经典最小二乘无解。系统采用广义逆解法,在最小范数条件下求解参数改正数,适用于变形监测中的相对位移分析。

(3)拟稳平差。当网中部分点稳定、部分点变动时,系统对稳定点施加最小范数条件,对非稳定点保持经典平差性质。该方法在工程变形监测中优于纯自由网平差。

(4)IGGIII稳健估计。针对含粗差的观测值,系统采用IGGIII等价权函数(分段处理:小误差保留、中误差降权、大误差剔除),有效抵抗粗差对平差结果的污染。相较经典最小二乘,抗粗差能力显著提升。

(5)岭估计。当法方程系数矩阵病态(如复共线性严重)时,最小二乘解方差放大。岭估计通过引入岭参数λ,以微小偏差换取方差大幅降低,提升数值稳定性。

(6)赫尔默特方差分量估计。当不同类观测值(方向、距离)先验方差不一致时,系统通过迭代估计各类观测值的方差分量,自动调整各类权比,提升平差合理性。

算法选择流程采用优先级判定:先判起算数据完备性 → 再判观测值粗差风险 → 最后判法方程病态程度,逐层筛选确定最终算法。

2.3 第三层:精度评定方法——多指标综合评估

平差完成后,系统从三个层面进行精度评定:

(1)验后单位权中误差。计算公式为:

σ̂₀ = √[ VTPV / ( n − t ) ]

其中n为观测值总数,t为必要观测数,n−t为多余观测数(自由度)。该指标反映整网内部符合精度。

(2)参数中误差。由协因数阵 QXX = (ATPA)−1 取对角线元素,乘以σ̂₀²得到各参数方差,开方得参数中误差。可用于评定待定点坐标精度。

(3)协方差矩阵。系统输出完整协方差阵 DXX = σ̂₀²·QXX,可用于点位误差椭圆绘制、相对误差椭圆分析及网的优化设计。

三、6种平差算法对比

算法名称 数学模型 适用场景 优点 局限
经典最小二乘平差 VTPV = min 起算数据充足、观测值正态分布的常规控制网 理论严密、无偏最优、计算效率高 对粗差敏感、病态时数值不稳定
秩亏自由网平差 VTPV = min 且 XTX = min 无起算数据的独立网、变形监测相对位移分析 无需已知点、适用秩亏情形 结果为相对位置、需基准转换
拟稳平差 VTPV = min 且 XsTXs = min 部分点稳定、部分点变动的工程监测网 兼顾稳定基准与变形分析 稳定点选择需先验信息
IGGIII稳健估计 等价权迭代 VTP̄V = min 含粗差观测值、可靠性要求高的控制网 抗粗差能力强、自动剔除坏值 迭代收敛耗时、初值敏感
岭估计 (ATPA + λI)X = ATPL 法方程病态、复共线性严重的控制网 降低均方误差、数值稳定 引入偏差、岭参数λ选取主观
赫尔默特方差分量估计 迭代调整各类观测值方差分量 多类观测值(方向+距离)权比未知网 自动定权、权比合理 迭代过程、需足够多余观测

四、关键精度参数

依据 GB/T 18341-2021《大比例尺地形图机助制图规范》及测量平差原理,控制网按等级划分的测角中误差限差如下:

等级 测角中误差 m 典型应用场景
一等m ≤ 0.5″国家一等三角网、高精度监测
二等0.5″ < m ≤ 0.7″国家二等三角网、城市首级控制
三等0.7″ < m ≤ 1.3″城市三等控制、工程首级网
四等1.3″ < m ≤ 1.8″城市四等控制、工程基本网
五等1.8″ < m ≤ 3.6″图根控制、测图控制网

系统在平差后会自动比对验后单位权中误差与对应等级限差,超限自动告警并提示可能的粗差来源,形成精度闭环管理。

五、结论

《全站仪任意网测量2023》系统通过"数据驱动识别 → 算法自适应匹配 → 多指标精度评定"的三层架构,将传统平差作业中高度依赖经验的算法选型环节标准化、自动化。系统覆盖60余种控制网类型,集成6种主流平差算法,在保证数学严密性的同时显著降低作业门槛,是控制网平差软件智能化的重要实践。

参考文献:
[1] 武汉大学测绘学院测量平差学科组. 测量平差原理[M]. 武汉: 武汉大学出版社, 2014.
[2] 国家市场监督管理总局. GB/T 18341-2021 大比例尺地形图机助制图规范[S]. 北京: 中国标准出版社, 2021.
[3] 闪速工作网. 全站仪任意网测量2023软件系统[软著登字第6122125号]. 中国版权保护中心, 2023.
[4] 崔希璋, 於宗俦, 陶本藻, 等. 广义测量平差(新版)[M]. 武汉: 武汉大学出版社, 2009.
[5] 周江文. 经典误差理论与抗差估计[J]. 测绘学报, 1989, 18(2): 115-120.