CN112906198B - 任意平面测量控制网非线性最小二乘平差数据处理方法 - Google Patents

任意平面测量控制网非线性最小二乘平差数据处理方法 Download PDF

Info

Publication number
CN112906198B
CN112906198B CN202110084166.6A CN202110084166A CN112906198B CN 112906198 B CN112906198 B CN 112906198B CN 202110084166 A CN202110084166 A CN 202110084166A CN 112906198 B CN112906198 B CN 112906198B
Authority
CN
China
Prior art keywords
format
nonlinear
control network
adjustment
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110084166.6A
Other languages
English (en)
Other versions
CN112906198A (zh
Inventor
范文辉
范晓燕
游为
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202110084166.6A priority Critical patent/CN112906198B/zh
Publication of CN112906198A publication Critical patent/CN112906198A/zh
Application granted granted Critical
Publication of CN112906198B publication Critical patent/CN112906198B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Computing Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

本发明公开了任意平面测量控制网非线性最小二乘平差数据处理方法，提出了快速序列二次规划算法，用于求解非线性最小二乘平差模型，使得控制点概略坐标解算与最小二乘平差同步进行，最近更新了其快速步长调整算法，修正了算法的理论缺陷和漏洞，通过“原线性最小二乘平差模型‑转换为非线性规划模型‑再转换为序列二次规划模型‑又转换为序列线性平差模型”这一组合的适于计算机编程计算的非线性组合算法，从理论找出了适用于任意平面测量控制网的坐标自动解算的方法，解决了以往软件受网形限制不通用、需要手工干预计算控制点概略坐标的难题，实现了测量平差与数据处理软件系统通用性和全自动化。

Description

任意平面测量控制网非线性最小二乘平差数据处理方法
技术领域
本发明属于工程测量技术领域，具体涉及任意平面测量控制网非线性最小二乘平差数据处理方法。
背景技术
现有的测量平差与数据处理软件，都只能处理常规网形，比如导线网、三角网(锁)、大地四边形、中点多边形、角度和距离交会等等。对于按照最小二乘平差原理组建的任意控制网网形，由于概略坐标的解算问题没有彻底解决，导致后续的平差计算无法进行下去。为解决此问题，现有的测量平差与数据处理软件要求前期的外业测量不能按照最小二乘平差原理任意布设控制网的网形，控制网的网形只能选取常规网形，导致外业作业严重受限，不能根据实际情况灵活布网，工作效率大幅度下降。或者可以任意布网，但后续需要人工完成概略坐标的计算，无法实现全过程的信息化处理。
发明内容
本发明目的在于提供任意平面测量控制网非线性最小二乘平差数据处理方法，用于解决上解决上述场景的问题，如现有的测量平差与数据处理软件，都只能处理常规网形，比如导线网、三角网(锁)、大地四边形、中点多边形、角度和距离交会等等。对于按照最小二乘平差原理组建的任意控制网网形，由于概略坐标的解算问题没有彻底解决，导致后续的平差计算无法进行下去。为解决此问题，现有的测量平差与数据处理软件要求前期的外业测量不能按照最小二乘平差原理任意布设控制网的网形，控制网的网形只能选取常规网形，导致外业作业严重受限，不能根据实际情况灵活布网，工作效率大幅度下降。或者可以任意布网，但后续需要人工完成概略坐标的计算，无法实现全过程的信息化处理。
为实现上述目的，本发明的技术方案是：
任意平面测量控制网非线性最小二乘平差数据处理方法，包括以下步骤：
S1：提出快速序列二次规划算法，用于求解非线性最小二乘平差模型，使得控制点概略坐标解算与最小二乘平差同步进行，更新其快速步长调整算法，得到“原线性最小二乘平差模型-转换为非线性规划模型-再转换为序列二次规划模型-又转换为序列线性平差模型”这一组合的适于计算机编程计算的非线性组合算法；
S2：在步骤S1的基础上，提出间接平差模型转化为条件平差模型的矩阵运算方法，采用先间接平差、后条件平差的两套平差的数据处理方法，使得数据处理过程引入全面的质量检查机制；
S3：在步骤S2的基础上，推导基于坐标变换的快速相对点位误差椭圆在长半轴E和短半轴F、纵坐标X和横坐标Y、边长的纵向和横向方向上的算法；
S4：在步骤S3的基础上，测量控制网非线性最小二乘平差数据处理软件系统生成Autodesk AutoCAD格式脚本文件的控制网平面图；
步骤S1中，适于计算机编程计算的非线性组合算法数学模型如下：
原非线性最小二乘平差模型：
Figure BDA0002910310070000021
转换为非线性规划模型：
Figure BDA0002910310070000022
再转换为序列二次规划模型：
Figure BDA0002910310070000023
又转换为序列线性平差模型：
Figure BDA0002910310070000024
式中
Figure BDA0002910310070000031
为观测值的平差值向量，L为观测值向量，V为观测值的改正数向量，B为平差值方程中参数向量的系数矩阵，C为平差值方程的常数向量，
Figure BDA0002910310070000032
为参数的平差值向量，X0为参数的概略值向量，
Figure BDA0002910310070000033
为参数的改正数向量，P为观测值的权阵，Q为观测值的协因数阵，D为观测值的协方差阵，σ0为单位权中误差，Xk为非线性规划迭代算法第k次迭代的参数过程值，xk为非线性规划迭代算法第k次迭代的下降方向，λk为非线性规划迭代算法第k次迭代的下降方向步长，f为任意非线性函数，n为观测值个数，t为参数个数，r为多余观测数。
进一步的，步骤S2中，建立的间接平差的数学模型：
Figure BDA0002910310070000034
间接平差转换为条件平差的过程方法：
建立增广矩阵
Figure BDA0002910310070000035
对其进行行变换
Figure BDA0002910310070000036
使其转化为
Figure BDA0002910310070000037
可得
Figure BDA0002910310070000038
完成向条件平差模型的转换。
转换后的条件平差的数学模型：
Figure BDA0002910310070000039
式中x为参数的改正数向量，B为误差方程中x的系数矩阵，l为误差方程的常系数向量，C为增广矩阵，I为单位阵，R为行变换左乘矩阵，V为观测值的改正数向量，P为观测值的权阵，A为条件方程的系数矩阵，W为条件方程的闭合差列向量，D为观测值的协方差阵，Q为观测值的协因数阵，σ0为单位权中误差，n为观测值个数，t为参数个数，r为多余观测数。
进一步的，步骤S3中，控制网中任意一边坐标增量的矩阵表达为：
Figure BDA0002910310070000041
若把间接平差的参数向量的协因数分解表达为：
Figure BDA0002910310070000042
则控制网中任意一边坐标增量的协因数阵，既相对点位协因数阵为：
Figure BDA0002910310070000043
坐标系(E,F)旋转ψ角度后转换为坐标系(X,Y)的矩阵关系：
Figure BDA0002910310070000044
坐标向量
Figure BDA0002910310070000045
的协因数阵为：
Figure BDA0002910310070000051
针对平面测量控制网中的任意一边的相对误差椭圆，通过上述坐标系旋转变换方法，建立坐标系(E,F)、坐标系(ΔX,ΔY)、坐标系(横,纵)间任意2坐标系的转换关系，其中
Figure BDA0002910310070000052
转换后的协因数标量关系为：
Figure BDA0002910310070000053
Figure BDA0002910310070000054
式中x和y为平面直角坐标，Δx和Δy为坐标增量，Q为观测值的协因数阵，X和Y为平面直角坐标系的坐标轴方向，E和F分别为相对误差椭圆的长半轴和短半轴，
Figure BDA0002910310070000055
为坐标系旋转角，ψ为相对误差椭圆长半轴的坐标方位角，S是任一边的边长，α为任一边的坐标方位角，σ0为单位权中误差，m为中误差，ρ″为弧度化秒常数206265，K为计算过程变量。
进一步的，步骤S4中测量控制网非线性最小二乘平差数据处理软件系统生成Autodesk AutoCAD格式脚本文件的控制网平面图中：
数据处理软件系统获取测量控制网非线性最小二乘平差数据，并识别其实时数据格式；
数据处理软件系统中预存有测量控制网非线性最小二乘平差数据的标准数据格式，数据处理软件系统将所述实时数据格式与所述标准数据格式进行格式匹配判断，若格式匹配判断的结果为所述实时数据格式与所述标准数据格式匹配，则数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据，若格式匹配判断的结果为所述实时数据格式与所述标准数据格式不匹配，则数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据；
数据处理软件系统以实际接收的测量控制网非线性最小二乘平差数据为基准，输出Autodesk AutoCAD格式脚本文件的控制网平面图。
进一步的，所述格式匹配判断具体为：
数据处理软件系统通过常用判断单元将所述实时数据格式与所述标准数据格式进行第一次格式匹配判断，若第一次格式匹配判断结果为所述实时数据格式与所述标准数据格式匹配，则以第一次格式匹配判断结果作为实际的格式匹配判断结果，数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据，若第一次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，则数据处理软件系统通过备用判断单元将所述实时数据格式与所述标准数据格式进行第二次格式匹配判断，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式匹配，则以第二次格式匹配判断结果作为实际的格式匹配判断结果，数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据。
进一步的，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，数据处理软件系统发出测量控制网非线性最小二乘平差数据的数据格式异常报警。
进一步的，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据，通过异常数据自动备份模块对数据处理软件系统拒收的测量控制网非线性最小二乘平差数据进行自动备份。
进一步的，步骤S4中测量控制网非线性最小二乘平差数据处理软件系统生成Autodesk AutoCAD格式脚本文件的控制网平面图具体为：
测量控制网非线性最小二乘平差数据处理软件系统将数据处理过程中点位及其精度数据，依据AutoCAD命令格式、控制网面图的图示规则，批量生成控制网平面图的Autodesk AutoCAD格式脚本文件，然后在Autodesk AutoCAD软件系统里运行脚本方式，完成自动绘制控制网平面图，最后保存为DWG格式图形文件。
与现有技术相比，本发明所具有的有益效果为：
本方案的一个创新点在于，提出了快速序列二次规划算法，用于求解非线性最小二乘平差模型，使得控制点概略坐标解算与最小二乘平差同步进行，最近更新了其快速步长调整算法，修正了算法的理论缺陷和漏洞，通过“原线性最小二乘平差模型-转换为非线性规划模型-再转换为序列二次规划模型-又转换为序列线性平差模型”这一组合的适于计算机编程计算的非线性组合算法，从理论找出了适用于任意平面测量控制网的坐标自动解算的方法，解决了以往软件受网形限制不通用、需要手工干预计算控制点概略坐标的难题，实现了测量平差与数据处理软件系统通用性和全自动化。
本方案的一个创新点在于，提出了间接平差模型转化为条件平差模型的矩阵运算方法，采用了先间接平差、后条件平差的2套平差的数据处理方法，使得数据处理过程引入了较全面的质量检查机制，实现了软件的查错和质量检查功能。
本方案的一个创新点在于，推导了基于坐标变换的快速相对点位误差椭圆在长半轴E和短半轴F、纵坐标X和横坐标Y、边长的纵向和横向等方向上的算法，适应了最新版的各类工程测量规范对精度评定的需求。
本方案的一个创新点在于，测量平差与数据处理软件系统通过生成AutodeskAutoCAD格式脚本文件，实现了Autodesk AutoCAD格式测量控制网平面图的批量绘制。
附图说明
图1是本发明具体实施方式的实施例结构示意图。
图2是本发明具体实施方式的实施例的平面测量控制网中的任意一边的相对误差椭圆示意图。
图3是本发明具体实施方式的实施例的平面测量控制网中三角网示意图。
图4是本发明具体实施方式的实施例生成的控制网示意图。
具体实施方式
下面结合本发明的附图1-4，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1：
如图1所示，提出任意平面测量控制网非线性最小二乘平差数据处理方法，包括以下步骤：
S1：提出快速序列二次规划算法，用于求解非线性最小二乘平差模型，使得控制点概略坐标解算与最小二乘平差同步进行，更新其快速步长调整算法，得到“原线性最小二乘平差模型-转换为非线性规划模型-再转换为序列二次规划模型-又转换为序列线性平差模型”这一组合的适于计算机编程计算的非线性组合算法；
S2：在步骤S1的基础上，提出间接平差模型转化为条件平差模型的矩阵运算方法，采用先间接平差、后条件平差的两套平差的数据处理方法，使得数据处理过程引入全面的质量检查机制；
S3：在步骤S2的基础上，推导基于坐标变换的快速相对点位误差椭圆在长半轴E和短半轴F、纵坐标X和横坐标Y、边长的纵向和横向方向上的算法；
S4：在步骤S3的基础上，测量控制网非线性最小二乘平差数据处理软件系统生成Autodesk AutoCAD格式脚本文件的控制网平面图；
步骤S1中，适于计算机编程计算的非线性组合算法数学模型如下：
原非线性最小二乘平差模型：
Figure BDA0002910310070000081
转换为非线性规划模型：
Figure BDA0002910310070000082
再转换为序列二次规划模型：
Figure BDA0002910310070000091
又转换为序列线性平差模型：
Figure BDA0002910310070000092
式中
Figure BDA0002910310070000093
为观测值的平差值向量，L为观测值向量，V为观测值的改正数向量，B为平差值方程中参数向量的系数矩阵，C为平差值方程的常数向量，
Figure BDA0002910310070000094
为参数的平差值向量，X0为参数的概略值向量，
Figure BDA0002910310070000095
为参数的改正数向量，P为观测值的权阵，Q为观测值的协因数阵，D为观测值的协方差阵，σ0为单位权中误差，Xk为非线性规划迭代算法第k次迭代的参数过程值，xk为非线性规划迭代算法第k次迭代的下降方向，λk为非线性规划迭代算法第k次迭代的下降方向步长，f为任意非线性函数，n为观测值个数，t为参数个数，r为多余观测数。
进一步的，步骤S2中，建立的间接平差的数学模型：
Figure BDA0002910310070000096
间接平差转换为条件平差的过程方法：
建立增广矩阵
Figure BDA0002910310070000097
对其进行行变换
Figure BDA0002910310070000098
使其转化为
Figure BDA0002910310070000101
可得
Figure BDA0002910310070000102
完成向条件平差模型的转换。
转换后的条件平差的数学模型：
Figure BDA0002910310070000103
式中x为参数的改正数向量，B为误差方程中x的系数矩阵，l为误差方程的常系数向量，C为增广矩阵，I为单位阵，R为行变换左乘矩阵，V为观测值的改正数向量，P为观测值的权阵，A为条件方程的系数矩阵，W为条件方程的闭合差列向量，D为观测值的协方差阵，Q为观测值的协因数阵，σ0为单位权中误差，n为观测值个数，t为参数个数，r为多余观测数。
进一步的，步骤S3中，控制网中任意一边坐标增量的矩阵表达为：
Figure BDA0002910310070000104
若把间接平差的参数向量的协因数分解表达为：
Figure BDA0002910310070000105
则控制网中任意一边坐标增量的协因数阵，既相对点位协因数阵为：
Figure BDA0002910310070000111
坐标系(E,F)旋转ψ角度后转换为坐标系(X,Y)的矩阵关系：
Figure BDA0002910310070000112
坐标向量
Figure BDA0002910310070000113
的协因数阵为：
Figure BDA0002910310070000114
如图2所示，针对平面测量控制网中的任意一边的相对误差椭圆，通过上述坐标系旋转变换方法，建立坐标系(E,F)、坐标系(ΔX,ΔY)、坐标系(横,纵)间任意2坐标系的转换关系，其中
Figure BDA0002910310070000115
转换后的协因数标量关系为：
Figure BDA0002910310070000116
Figure BDA0002910310070000117
式中x和y为平面直角坐标，Δx和Δy为坐标增量，Q为观测值的协因数阵，X和Y为平面直角坐标系的坐标轴方向，E和F分别为相对误差椭圆的长半轴和短半轴，
Figure BDA0002910310070000121
为坐标系旋转角，ψ为相对误差椭圆长半轴的坐标方位角，S是任一边的边长，α为任一边的坐标方位角，σ0为单位权中误差，m为中误差，ρ″为弧度化秒常数206265，K为计算过程变量。
进一步的，步骤S4中测量控制网非线性最小二乘平差数据处理软件系统生成Autodesk AutoCAD格式脚本文件的控制网平面图中：
数据处理软件系统获取测量控制网非线性最小二乘平差数据，并识别其实时数据格式；
数据处理软件系统中预存有测量控制网非线性最小二乘平差数据的标准数据格式，数据处理软件系统将所述实时数据格式与所述标准数据格式进行格式匹配判断，若格式匹配判断的结果为所述实时数据格式与所述标准数据格式匹配，则数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据，若格式匹配判断的结果为所述实时数据格式与所述标准数据格式不匹配，则数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据；
数据处理软件系统以实际接收的测量控制网非线性最小二乘平差数据为基准，输出Autodesk AutoCAD格式脚本文件的控制网平面图。即测量控制网非线性最小二乘平差数据处理软件系统以实际观测数据为基准，输出Autodesk AutoCAD格式脚本文件的控制网平面图。
进一步的，所述格式匹配判断具体为：
数据处理软件系统通过常用判断单元将所述实时数据格式与所述标准数据格式进行第一次格式匹配判断，若第一次格式匹配判断结果为所述实时数据格式与所述标准数据格式匹配，则以第一次格式匹配判断结果作为实际的格式匹配判断结果，数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据，若第一次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，则数据处理软件系统通过备用判断单元将所述实时数据格式与所述标准数据格式进行第二次格式匹配判断，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式匹配，则以第二次格式匹配判断结果作为实际的格式匹配判断结果，数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据。
进一步的，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，数据处理软件系统发出测量控制网非线性最小二乘平差数据的数据格式异常报警。
进一步的，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据，通过异常数据自动备份模块对数据处理软件系统拒收的测量控制网非线性最小二乘平差数据进行自动备份。
进一步的，步骤S4中测量控制网非线性最小二乘平差数据处理软件系统生成Autodesk AutoCAD格式脚本文件的控制网平面图具体为：
测量控制网非线性最小二乘平差数据处理软件系统将数据处理过程中点位及其精度数据，依据AutoCAD命令格式、控制网面图的图示规则，批量生成控制网平面图的Autodesk AutoCAD格式脚本文件，然后在Autodesk AutoCAD软件系统里运行脚本方式，完成自动绘制控制网平面图，最后保存为DWG格式图形文件。
应用算例：
在图3的三角网中，A、B为已知点，CD为基线边，又已知BE的方位角αBE，观测了全部角度值列于角度观测值数据表，试求各待定点坐标的平差值和F点的点位中误差。
已知起算数据表
Figure BDA0002910310070000131
角度观测值数据表
Figure BDA0002910310070000132
Figure BDA0002910310070000141
解：
①适于计算机编程计算的“原线性最小二乘平差模型-转换为非线性规划模型-再转换为序列二次规划模型-又转换为序列线性平差模型”非线性组合算法，应用于任意平面网形的测量控制网非线性最小二乘平差与数据处理软件系统，输出成果表如下：
Figure BDA0002910310070000142
Figure BDA0002910310070000151
该表为：观测值、平差值、改正数、可靠性指标一览表
点位及其精度成果表
Figure BDA0002910310070000152
②采用先间接平差、后条件平差的2套平差的数据处理方法，数据处理过程具有较全面的质量检查机制，实现了软件的查错和质量检查功能，输出成果表如下：
Figure BDA0002910310070000153
③基于坐标变换的快速相对点位误差椭圆在长半轴E和短半轴F、纵坐标X和横坐标Y、边长的纵向和横向等方向上的算法，适应最新版的各类工程测量规范对精度评定的需求，输出成果表如下：
相对点位及其精度成果表1
Figure BDA0002910310070000161
相对点位及其精度成果表2
Figure BDA0002910310070000162
Figure BDA0002910310070000171
④测量平差与数据处理软件系统通过生成Autodesk AutoCAD格式脚本文件，实现了Autodesk AutoCAD格式测量控制网平面图的批量绘制，生成的控制网图(含点位、相对点位误差椭圆)如图4所示。
以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作用未超出本发明技术方案的范围时，均属于本发明的保护范围。

Claims (8)

1.任意平面测量控制网非线性最小二乘平差数据处理方法，其特征在于，包括以下步骤：
S1：提出快速序列二次规划算法，用于求解非线性最小二乘平差模型，使得控制点概略坐标解算与最小二乘平差同步进行，更新其快速步长调整算法，得到“原线性最小二乘平差模型-转换为非线性规划模型-再转换为序列二次规划模型-又转换为序列线性平差模型”这一组合的适于计算机编程计算的非线性组合算法；
S2：在步骤S1的基础上，提出间接平差模型转化为条件平差模型的矩阵运算方法，采用先间接平差、后条件平差的两套平差的数据处理方法，使得数据处理过程引入全面的质量检查机制；
S3：在步骤S2的基础上，推导基于坐标变换的快速相对点位误差椭圆在长半轴E和短半轴F、纵坐标X和横坐标Y、边长的纵向和横向方向上的算法；
S4：在步骤S3的基础上，测量控制网非线性最小二乘平差数据处理软件系统生成Autodesk AutoCAD格式脚本文件的控制网平面图；
步骤S1中，适于计算机编程计算的非线性组合算法数学模型如下：
原非线性最小二乘平差模型：
Figure FDA0003206884900000011
转换为非线性规划模型：
Figure FDA0003206884900000012
再转换为序列二次规划模型：
Figure FDA0003206884900000013
又转换为序列线性平差模型：
Figure FDA0003206884900000021
式中
Figure FDA0003206884900000022
为观测值的平差值向量，L为观测值向量，V为观测值的改正数向量，B为平差值方程中参数向量的系数矩阵，C为平差值方程的常数向量，
Figure FDA0003206884900000023
为参数的平差值向量，X0为参数的概略值向量，
Figure FDA0003206884900000024
为参数的改正数向量，P为观测值的权阵，Q为观测值的协因数阵，D为观测值的协方差阵，σ0为单位权中误差，Xk为非线性规划迭代算法第k次迭代的参数过程值，xk为非线性规划迭代算法第k次迭代的下降方向，λk为非线性规划迭代算法第k次迭代的下降方向步长，f为任意非线性函数，n为观测值个数，t为参数个数，r为多余观测数。
2.如权利要求1所述的任意平面测量控制网非线性最小二乘平差数据处理方法，其特征在于，步骤S2中，建立的间接平差的数学模型：
Figure FDA0003206884900000025
间接平差转换为条件平差的过程方法：
建立增广矩阵
Figure FDA0003206884900000026
对其进行行变换
Figure FDA0003206884900000027
使其转化为
Figure FDA0003206884900000028
可得
Figure FDA0003206884900000029
完成向条件平差模型的转换；
转换后的条件平差的数学模型：
Figure FDA0003206884900000031
式中x为参数的改正数向量，B为误差方程中x的系数矩阵，l为误差方程的常系数向量，C为增广矩阵，I为单位阵，R为行变换左乘矩阵，V为观测值的改正数向量，P为观测值的权阵，A为条件方程的系数矩阵，W为条件方程的闭合差列向量，D为观测值的协方差阵，Q为观测值的协因数阵，σ0为单位权中误差，n为观测值个数，t为参数个数，r为多余观测数。
3.如权利要求1所述的任意平面测量控制网非线性最小二乘平差数据处理方法，其特征在于，步骤S3中，控制网中任意一边坐标增量的矩阵表达为：
Figure FDA0003206884900000032
若把间接平差的参数向量的协因数分解表达为：
Figure FDA0003206884900000033
则控制网中任意一边坐标增量的协因数阵，既相对点位协因数阵为：
Figure FDA0003206884900000041
坐标系(E,F)旋转ψ角度后转换为坐标系(X,Y)的矩阵关系：
Figure FDA0003206884900000042
坐标向量
Figure FDA0003206884900000043
的协因数阵为：
Figure FDA0003206884900000044
针对平面测量控制网中的任意一边的相对误差椭圆，通过上述坐标系旋转变换方法，建立坐标系(E,F)、坐标系(ΔX,ΔY)、坐标系(横,纵)间任意2坐标系的转换关系，其中
Figure FDA0003206884900000045
转换后的协因数标量关系为：
Figure FDA0003206884900000046
Figure FDA0003206884900000047
式中x和y为平面直角坐标，Δx和Δy为坐标增量，Q为观测值的协因数阵，X和Y为平面直角坐标系的坐标轴方向，E和F分别为相对误差椭圆的长半轴和短半轴，
Figure FDA0003206884900000051
为坐标系旋转角，ψ为相对误差椭圆长半轴的坐标方位角，S是任一边的边长，α为任一边的坐标方位角，σ0为单位权中误差，m为中误差，ρ″为弧度化秒常数206265，K为计算过程变量。
4.如权利要求1所述的任意平面测量控制网非线性最小二乘平差数据处理方法，其特征在于，步骤S4中测量控制网非线性最小二乘平差数据处理软件系统生成AutodeskAutoCAD格式脚本文件的控制网平面图中：
数据处理软件系统获取测量控制网非线性最小二乘平差数据，并识别其实时数据格式；
数据处理软件系统中预存有测量控制网非线性最小二乘平差数据的标准数据格式，数据处理软件系统将所述实时数据格式与所述标准数据格式进行格式匹配判断，若格式匹配判断的结果为所述实时数据格式与所述标准数据格式匹配，则数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据，若格式匹配判断的结果为所述实时数据格式与所述标准数据格式不匹配，则数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据；
数据处理软件系统以实际接收的测量控制网非线性最小二乘平差数据为基准，输出Autodesk AutoCAD格式脚本文件的控制网平面图。
5.如权利要求4所述的任意平面测量控制网非线性最小二乘平差数据处理方法，其特征在于，所述格式匹配判断具体为：
数据处理软件系统通过常用判断单元将所述实时数据格式与所述标准数据格式进行第一次格式匹配判断，若第一次格式匹配判断结果为所述实时数据格式与所述标准数据格式匹配，则以第一次格式匹配判断结果作为实际的格式匹配判断结果，数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据，若第一次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，则数据处理软件系统通过备用判断单元将所述实时数据格式与所述标准数据格式进行第二次格式匹配判断，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式匹配，则以第二次格式匹配判断结果作为实际的格式匹配判断结果，数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据。
6.如权利要求5所述的任意平面测量控制网非线性最小二乘平差数据处理方法，其特征在于，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，数据处理软件系统发出测量控制网非线性最小二乘平差数据的数据格式异常报警。
7.如权利要求5所述的任意平面测量控制网非线性最小二乘平差数据处理方法，其特征在于，若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配，数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据，通过异常数据自动备份模块对数据处理软件系统拒收的测量控制网非线性最小二乘平差数据进行自动备份。
8.如权利要求1所述的任意平面测量控制网非线性最小二乘平差数据处理方法，其特征在于，步骤S4中测量控制网非线性最小二乘平差数据处理软件系统生成AutodeskAutoCAD格式脚本文件的控制网平面图具体为：
测量控制网非线性最小二乘平差数据处理软件系统将数据处理过程中点位及其精度数据，依据AutoCAD命令格式、控制网面图的图示规则，批量生成控制网平面图的AutodeskAutoCAD格式脚本文件，然后在Autodesk AutoCAD软件系统里运行脚本方式，完成自动绘制控制网平面图，最后保存为DWG格式图形文件。
CN202110084166.6A 2021-01-21 2021-01-21 任意平面测量控制网非线性最小二乘平差数据处理方法 Active CN112906198B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110084166.6A CN112906198B (zh) 2021-01-21 2021-01-21 任意平面测量控制网非线性最小二乘平差数据处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110084166.6A CN112906198B (zh) 2021-01-21 2021-01-21 任意平面测量控制网非线性最小二乘平差数据处理方法

Publications (2)

Publication Number Publication Date
CN112906198A CN112906198A (zh) 2021-06-04
CN112906198B true CN112906198B (zh) 2021-10-22

Family

ID=76118285

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110084166.6A Active CN112906198B (zh) 2021-01-21 2021-01-21 任意平面测量控制网非线性最小二乘平差数据处理方法

Country Status (1)

Country Link
CN (1) CN112906198B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104091064A (zh) * 2014-07-02 2014-10-08 北京航空航天大学 基于优化解空间搜索法的PS-DInSAR地表形变测量参数估计方法
CN110793508A (zh) * 2019-11-11 2020-02-14 中国人民解放军战略支援部队信息工程大学 三维控制网中控制点的观测数据处理方法与装置
CN111238454A (zh) * 2020-01-17 2020-06-05 中国人民解放军战略支援部队信息工程大学 一种空间三维坐标测量仪系统及其三维控制网测量方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2003289649A1 (en) * 2003-03-19 2004-10-11 Institute Of Acoutics, Chinese Academy Of Sciences Method and system for measuring flow layer velocities using correlation velocity measuring sonar
CN102607506B (zh) * 2012-03-01 2013-10-30 中国人民解放军空军工程设计研究局 高填方机场边坡单台全站仪的自由设站变形监测方法
US9122261B2 (en) * 2012-04-24 2015-09-01 Honeywell International Inc. Apparatus and method for real-time sequential quadratic programming in industrial process control systems
CN106354008B (zh) * 2016-08-10 2019-05-03 中铁第一勘察设计院集团有限公司 Cpiii控制网平差方法
CN107314763B (zh) * 2017-07-18 2019-05-07 上海海洋大学 一种基于制约函数非线性估计的卫星影像区域网平差方法
CN109752958B (zh) * 2019-02-20 2021-12-03 成都市勘察测绘研究院 一种地面控制网更新维护的自适应网型约束平差方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104091064A (zh) * 2014-07-02 2014-10-08 北京航空航天大学 基于优化解空间搜索法的PS-DInSAR地表形变测量参数估计方法
CN110793508A (zh) * 2019-11-11 2020-02-14 中国人民解放军战略支援部队信息工程大学 三维控制网中控制点的观测数据处理方法与装置
CN111238454A (zh) * 2020-01-17 2020-06-05 中国人民解放军战略支援部队信息工程大学 一种空间三维坐标测量仪系统及其三维控制网测量方法

Also Published As

Publication number Publication date
CN112906198A (zh) 2021-06-04

Similar Documents

Publication Publication Date Title
CN104965483B (zh) 一种基于稳健设计的多轴数控机床加工精度保持性优化方法
CN108873809B (zh) 一种螺旋锥齿轮的高阶齿面误差修正方法
CN103995496A (zh) 一种基于数字化测量的飞机部件高精度匹配零件加工方法
US20230384755A1 (en) Datum selection methods and systems for minimizing hole position errors in group hole machining of large components
CN110161965B (zh) 一种大型航天机匣斜孔的在机测量方法
CN112109084A (zh) 基于机器人关节角度补偿的末端位置补偿方法及其应用
Wang et al. 3D machining allowance analysis method for the large thin-walled aerospace component
CN110716194A (zh) 同时参考球与面的多波束三维激光成像系统的检校方法
CN110686655B (zh) 三维曲面工件的定位方法、装置及系统
CN112906198B (zh) 任意平面测量控制网非线性最小二乘平差数据处理方法
CN107066726B (zh) 一种数控机床旋转轴垂直度误差建模方法
CN111891176A (zh) 基于轨道几何状态测量仪的轨道测量数据整体平差方法
CN112923887B (zh) 任意平面测量控制网非线性最小二乘平差数据检验方法
CN112906197B (zh) 任意平面测量控制网基于坐标变换的数据处理方法
TW202302301A (zh) 機械手臂與輪廓感測器座標系相對關係之自動校正方法與系統
Xu et al. A new welding path planning method based on point cloud and deep learning
Pechenin et al. Development of a method of ICP algorithm accuracy improvement during shaped profiles and surfaces control
CN108638060B (zh) 多自由度机器人参数标定中冗余参数分析剔除方法
CN112643658A (zh) 基于sir降维dh模型的串联机器人自适应误差建模的标定方法
CN115795220A (zh) 一种数控加工中心工件坐标系计算方法及系统
Yang et al. Gap volume prediction for aircraft wing assembly
Kuang Second-order design: Shooting for maximum reliability
CN107367235B (zh) 一种红外面阵扫描系统的轴系误差标定方法
CN115953477B (zh) 坐标数据处理方法和装置
CN112308891B (zh) 一种近净成形类零件多公差约束的配准计算方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant