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

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

Info

Publication number
CN112923887A
CN112923887A CN202110082467.5A CN202110082467A CN112923887A CN 112923887 A CN112923887 A CN 112923887A CN 202110082467 A CN202110082467 A CN 202110082467A CN 112923887 A CN112923887 A CN 112923887A
Authority
CN
China
Prior art keywords
format
control network
nonlinear
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.)
Granted
Application number
CN202110082467.5A
Other languages
English (en)
Other versions
CN112923887B (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 CN202110082467.5A priority Critical patent/CN112923887B/zh
Publication of CN112923887A publication Critical patent/CN112923887A/zh
Application granted granted Critical
Publication of CN112923887B publication Critical patent/CN112923887B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
    • G01B21/04Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness by measuring coordinates of points
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD

Landscapes

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

Abstract

本发明公开了任意平面测量控制网非线性最小二乘平差数据检验方法,提出了间接平差模型转化为条件平差模型的矩阵运算方法,采用了先间接平差、后条件平差的2套平差的数据处理方法,使得数据处理过程引入了较全面的质量检查机制,实现了软件的查错和质量检查功能。

Description

任意平面测量控制网非线性最小二乘平差数据检验方法
技术领域
本发明属于工程测量技术领域,具体涉及任意平面测量控制网非线性最小二乘平差数据检验方法。
背景技术
现有的测量平差与数据处理软件,都只能处理常规网形,比如导线网、三角网(锁)、大地四边形、中点多边形、角度和距离交会等等。对于按照最小二乘平差原理组建的任意控制网网形,由于概略坐标的解算问题没有彻底解决,导致后续的平差计算无法进行下去。为解决此问题,现有的测量平差与数据处理软件要求前期的外业测量不能按照最小二乘平差原理任意布设控制网的网形,控制网的网形只能选取常规网形,导致外业作业严重受限,不能根据实际情况灵活布网,工作效率大幅度下降。或者可以任意布网,但后续需要人工完成概略坐标的计算,无法实现全过程的信息化处理。实现高效、全面、彻底的工程测量控制网的质量检查,对确保成果的精度和可靠性至关重要,一直都是工程测量领域关注的重点。现有的测量平差与数据处理软件系统具有可靠性指标(多余观测分量)的计算,但还没能实现完整的闭合差检查。
发明内容
本发明目的在于提供任意平面测量控制网非线性最小二乘平差数据检验方法,用于解决上解决上述场景的问题,如现有的测量平差与数据处理软件,都只能处理常规网形,比如导线网、三角网(锁)、大地四边形、中点多边形、角度和距离交会等等。对于按照最小二乘平差原理组建的任意控制网网形,由于概略坐标的解算问题没有彻底解决,导致后续的平差计算无法进行下去。为解决此问题,现有的测量平差与数据处理软件要求前期的外业测量不能按照最小二乘平差原理任意布设控制网的网形,控制网的网形只能选取常规网形,导致外业作业严重受限,不能根据实际情况灵活布网,工作效率大幅度下降。或者可以任意布网,但后续需要人工完成概略坐标的计算,无法实现全过程的信息化处理。实现高效、全面、彻底的工程测量控制网的质量检查,对确保成果的精度和可靠性至关重要,一直都是工程测量领域关注的重点。现有的测量平差与数据处理软件系统具有可靠性指标(多余观测分量)的计算,但还没能实现完整的闭合差检查。
为实现上述目的,本发明的技术方案是:
任意平面测量控制网非线性最小二乘平差数据检验方法,包括以下步骤:
S1:提出快速序列二次规划算法,用于求解非线性最小二乘平差模型,使得控制点概略坐标解算与最小二乘平差同步进行,更新其快速步长调整算法,得到“原线性最小二乘平差模型-转换为非线性规划模型-再转换为序列二次规划模型-又转换为序列线性平差模型”这一组合的适于计算机编程计算的非线性组合算法;
S2:在步骤S1的基础上,提出间接平差模型转化为条件平差模型的矩阵运算方法,采用先间接平差、后条件平差的两套平差的数据处理方法,使得数据处理过程引入全面的质量检查机制;
S3:在步骤S2的基础上,推导基于坐标变换的快速相对点位误差椭圆在长半轴E和短半轴F、纵坐标X和横坐标Y、边长的纵向和横向方向上的算法;
S4:在步骤S3的基础上,测量控制网非线性最小二乘平差数据处理软件系统生成Autodesk AutoCAD格式脚本文件的控制网平面图;
步骤S1中,适于计算机编程计算的非线性组合算法数学模型如下:
原非线性最小二乘平差模型:
Figure BDA0002909835390000021
转换为非线性规划模型:
Figure BDA0002909835390000022
再转换为序列二次规划模型:
Figure BDA0002909835390000023
又转换为序列线性平差模型:
Figure BDA0002909835390000031
式中
Figure BDA0002909835390000032
为观测值的平差值向量,L为观测值向量,V为观测值的改正数向量,B为平差值方程中参数向量的系数矩阵,C为平差值方程的常数向量,
Figure BDA0002909835390000033
为参数的平差值向量,X0为参数的概略值向量,
Figure BDA0002909835390000034
为参数的改正数向量,P为观测值的权阵,Q为观测值的协因数阵,D为观测值的协方差阵,σ0为单位权中误差,Xk为非线性规划迭代算法第k次迭代的参数过程值,xk为非线性规划迭代算法第k次迭代的下降方向,λk为非线性规划迭代算法第k次迭代的下降方向步长,f为任意非线性函数,n为观测值个数,t为参数个数,r为多余观测数。
进一步的,步骤S2中,建立的间接平差的数学模型:
Figure BDA0002909835390000035
间接平差转换为条件平差的过程方法:
建立增广矩阵
Figure BDA0002909835390000036
对其进行行变换
Figure BDA0002909835390000037
使其转化为
Figure BDA0002909835390000038
可得
Figure BDA0002909835390000039
完成向条件平差模型的转换。
转换后的条件平差的数学模型:
Figure BDA0002909835390000041
式中x为参数的改正数向量,B为误差方程中x的系数矩阵,l为误差方程的常系数向量,C为增广矩阵,I为单位阵,R为行变换左乘矩阵,V为观测值的改正数向量,P为观测值的权阵,A为条件方程的系数矩阵,W为条件方程的闭合差列向量,D为观测值的协方差阵,Q为观测值的协因数阵,σ0为单位权中误差,n为观测值个数,t为参数个数,r为多余观测数。
进一步的,步骤S3中,控制网中任意一边坐标增量的矩阵表达为:
Figure BDA0002909835390000042
若把间接平差的参数向量的协因数分解表达为:
Figure BDA0002909835390000043
则控制网中任意一边坐标增量的协因数阵,既相对点位协因数阵为:
Figure BDA0002909835390000044
坐标系(E,F)旋转ψ角度后转换为坐标系(X,Y)的矩阵关系:
Figure BDA0002909835390000051
坐标向量
Figure BDA0002909835390000052
的协因数阵为:
Figure BDA0002909835390000053
针对平面测量控制网中的任意一边的相对误差椭圆,通过上述坐标系旋转变换方法,建立坐标系(E,F)、坐标系(ΔX,ΔY)、坐标系(横,纵)间任意2坐标系的转换关系,其中
Figure BDA0002909835390000054
转换后的协因数标量关系为:
Figure BDA0002909835390000055
Figure BDA0002909835390000056
式中x和y为平面直角坐标,Δx和Δy为坐标增量,Q为观测值的协因数阵,X和Y为平面直角坐标系的坐标轴方向,E和F分别为相对误差椭圆的长半轴和短半轴,
Figure BDA0002909835390000057
为坐标系旋转角,ψ为相对误差椭圆长半轴的坐标方位角,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 BDA0002909835390000081
转换为非线性规划模型:
Figure BDA0002909835390000091
再转换为序列二次规划模型:
Figure BDA0002909835390000092
又转换为序列线性平差模型:
Figure BDA0002909835390000093
式中
Figure BDA0002909835390000094
为观测值的平差值向量,L为观测值向量,V为观测值的改正数向量,B为平差值方程中参数向量的系数矩阵,C为平差值方程的常数向量,
Figure BDA0002909835390000095
为参数的平差值向量,X0为参数的概略值向量,
Figure BDA0002909835390000096
为参数的改正数向量,P为观测值的权阵,Q为观测值的协因数阵,D为观测值的协方差阵,σ0为单位权中误差,Xk为非线性规划迭代算法第k次迭代的参数过程值,xk为非线性规划迭代算法第k次迭代的下降方向,λk为非线性规划迭代算法第k次迭代的下降方向步长,f为任意非线性函数,n为观测值个数,t为参数个数,r为多余观测数。
进一步的,步骤S2中,建立的间接平差的数学模型:
Figure BDA0002909835390000097
间接平差转换为条件平差的过程方法:
建立增广矩阵
Figure BDA0002909835390000101
对其进行行变换
Figure BDA0002909835390000102
使其转化为
Figure BDA0002909835390000103
可得
Figure BDA0002909835390000104
完成向条件平差模型的转换。
转换后的条件平差的数学模型:
Figure BDA0002909835390000105
式中x为参数的改正数向量,B为误差方程中x的系数矩阵,l为误差方程的常系数向量,C为增广矩阵,I为单位阵,R为行变换左乘矩阵,V为观测值的改正数向量,P为观测值的权阵,A为条件方程的系数矩阵,W为条件方程的闭合差列向量,D为观测值的协方差阵,Q为观测值的协因数阵,σ0为单位权中误差,n为观测值个数,t为参数个数,r为多余观测数。
进一步的,步骤S3中,控制网中任意一边坐标增量的矩阵表达为:
Figure BDA0002909835390000106
若把间接平差的参数向量的协因数分解表达为:
Figure BDA0002909835390000111
则控制网中任意一边坐标增量的协因数阵,既相对点位协因数阵为:
Figure BDA0002909835390000112
坐标系(E,F)旋转ψ角度后转换为坐标系(X,Y)的矩阵关系:
Figure BDA0002909835390000113
坐标向量
Figure BDA0002909835390000114
的协因数阵为:
Figure BDA0002909835390000115
如图2所示,针对平面测量控制网中的任意一边的相对误差椭圆,通过上述坐标系旋转变换方法,建立坐标系(E,F)、坐标系(ΔX,ΔY)、坐标系(横,纵)间任意2坐标系的转换关系,其中
Figure BDA0002909835390000116
转换后的协因数标量关系为:
Figure BDA0002909835390000121
Figure BDA0002909835390000122
式中x和y为平面直角坐标,Δx和Δy为坐标增量,Q为观测值的协因数阵,X和Y为平面直角坐标系的坐标轴方向,E和F分别为相对误差椭圆的长半轴和短半轴,
Figure BDA0002909835390000123
为坐标系旋转角,ψ为相对误差椭圆长半轴的坐标方位角,S是任一边的边长,α为任一边的坐标方位角,σ0为单位权中误差,m为中误差,ρ″为弧度化秒常数206265,K为计算过程变量。
进一步的,步骤S4中测量控制网非线性最小二乘平差数据处理软件系统生成Autodesk AutoCAD格式脚本文件的控制网平面图中:
数据处理软件系统获取测量控制网非线性最小二乘平差数据,并识别其实时数据格式;
数据处理软件系统中预存有测量控制网非线性最小二乘平差数据的标准数据格式,数据处理软件系统将所述实时数据格式与所述标准数据格式进行格式匹配判断,若格式匹配判断的结果为所述实时数据格式与所述标准数据格式匹配,则数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据,若格式匹配判断的结果为所述实时数据格式与所述标准数据格式不匹配,则数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据;
数据处理软件系统以实际接收的测量控制网非线性最小二乘平差数据为基准,输出Autodesk AutoCAD格式脚本文件的控制网平面图。
进一步的,所述格式匹配判断具体为:
数据处理软件系统通过常用判断单元将所述实时数据格式与所述标准数据格式进行第一次格式匹配判断,若第一次格式匹配判断结果为所述实时数据格式与所述标准数据格式匹配,则以第一次格式匹配判断结果作为实际的格式匹配判断结果,数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据,若第一次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配,则数据处理软件系统通过备用判断单元将所述实时数据格式与所述标准数据格式进行第二次格式匹配判断,若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式匹配,则以第二次格式匹配判断结果作为实际的格式匹配判断结果,数据处理软件系统接收所述实时数据格式下的测量控制网非线性最小二乘平差数据,若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配,数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据。
进一步的,若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配,数据处理软件系统发出测量控制网非线性最小二乘平差数据的数据格式异常报警。
进一步的,若第二次格式匹配判断结果为所述实时数据格式与所述标准数据格式不匹配,数据处理软件系统拒收所述实时数据格式下的测量控制网非线性最小二乘平差数据,通过异常数据自动备份模块对数据处理软件系统拒收的测量控制网非线性最小二乘平差数据进行自动备份。
进一步的,步骤S4中测量控制网非线性最小二乘平差数据处理软件系统生成Autodesk AutoCAD格式脚本文件的控制网平面图具体为:
测量控制网非线性最小二乘平差数据处理软件系统将数据处理过程中点位及其精度数据,依据AutoCAD命令格式、控制网平面图的图示规则,批量生成控制网平面图的Autodesk AutoCAD格式脚本文件,然后在Autodesk AutoCAD软件系统里运行脚本方式,完成自动绘制控制网平面图,最后保存为DWG格式图形文件。
应用算例:
在图3的三角网中,A、B为已知点,CD为基线边,又已知BE的方位角αBE,观测了全部角度值列于角度观测值数据表,试求各待定点坐标的平差值和F点的点位中误差。
已知起算数据表
Figure BDA0002909835390000131
Figure BDA0002909835390000141
角度观测值数据表
Figure BDA0002909835390000142
解:
①适于计算机编程计算的“原线性最小二乘平差模型-转换为非线性规划模型-再转换为序列二次规划模型-又转换为序列线性平差模型”非线性组合算法,应用于任意平面网形的测量控制网非线性最小二乘平差与数据处理软件系统,输出成果表如下:
Figure BDA0002909835390000143
Figure BDA0002909835390000151
该表为:观测值、平差值、改正数、可靠性指标一览表
点位及其精度成果表
Figure BDA0002909835390000152
②采用先间接平差、后条件平差的2套平差的数据处理方法,数据处理过程具有较全面的质量检查机制,实现了软件的查错和质量检查功能,输出成果表如下:
Figure BDA0002909835390000161
③基于坐标变换的快速相对点位误差椭圆在长半轴E和短半轴F、纵坐标X和横坐标Y、边长的纵向和横向等方向上的算法,适应最新版的各类工程测量规范对精度评定的需求,输出成果表如下:
相对点位及其精度成果表1
Figure BDA0002909835390000162
Figure BDA0002909835390000171
相对点位及其精度成果表2
Figure BDA0002909835390000172
④测量平差与数据处理软件系统通过生成Autodesk AutoCAD格式脚本文件,实现了Autodesk AutoCAD格式测量控制网平面图的批量绘制,生成的控制网图(含点位、相对点位误差椭圆)如图4所示。
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。

Claims (6)

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

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110082467.5A CN112923887B (zh) 2021-01-21 2021-01-21 任意平面测量控制网非线性最小二乘平差数据检验方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110082467.5A CN112923887B (zh) 2021-01-21 2021-01-21 任意平面测量控制网非线性最小二乘平差数据检验方法

Publications (2)

Publication Number Publication Date
CN112923887A true CN112923887A (zh) 2021-06-08
CN112923887B CN112923887B (zh) 2021-08-20

Family

ID=76164061

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110082467.5A Active CN112923887B (zh) 2021-01-21 2021-01-21 任意平面测量控制网非线性最小二乘平差数据检验方法

Country Status (1)

Country Link
CN (1) CN112923887B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104048645A (zh) * 2014-06-25 2014-09-17 山东理工大学 线性拟合地面扫描点云整体定向方法
CN104933292A (zh) * 2015-05-08 2015-09-23 同济大学 一种结构化总体最小二乘条件平差方法
US20160275416A1 (en) * 2015-03-20 2016-09-22 Nec Laboratories America, Inc. Fast Distributed Nonnegative Matrix Factorization and Completion for Big Data Analytics
CN106980598A (zh) * 2017-04-11 2017-07-25 中铁第勘察设计院集团有限公司 自由测站边角交会网的三维平差严密方法
CN108279643A (zh) * 2018-01-29 2018-07-13 西南交通大学 一种基于测量点与自适应差分进化算法的工件姿态调整方法
CN111578881A (zh) * 2020-06-04 2020-08-25 中国电建集团华东勘测设计研究院有限公司 拱坝垂线监测数据的精确处理方法
CN112212843A (zh) * 2020-10-20 2021-01-12 中铁二院工程集团有限责任公司 一种铁路cpⅱ和cpⅲ控制网同步测量与联合平差的方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104048645A (zh) * 2014-06-25 2014-09-17 山东理工大学 线性拟合地面扫描点云整体定向方法
US20160275416A1 (en) * 2015-03-20 2016-09-22 Nec Laboratories America, Inc. Fast Distributed Nonnegative Matrix Factorization and Completion for Big Data Analytics
CN104933292A (zh) * 2015-05-08 2015-09-23 同济大学 一种结构化总体最小二乘条件平差方法
CN106980598A (zh) * 2017-04-11 2017-07-25 中铁第勘察设计院集团有限公司 自由测站边角交会网的三维平差严密方法
CN108279643A (zh) * 2018-01-29 2018-07-13 西南交通大学 一种基于测量点与自适应差分进化算法的工件姿态调整方法
CN111578881A (zh) * 2020-06-04 2020-08-25 中国电建集团华东勘测设计研究院有限公司 拱坝垂线监测数据的精确处理方法
CN112212843A (zh) * 2020-10-20 2021-01-12 中铁二院工程集团有限责任公司 一种铁路cpⅱ和cpⅲ控制网同步测量与联合平差的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
周拥军等: "附参数的条件平差与按行独立的加权总体最小二乘法估计的一致性研究", 《测绘学报》 *
游为等: "水准网闭合差自动解算的新方法", 《测绘工程》 *
赵一晗等: "控制网闭合环搜索算法的探讨", 《铁道勘察》 *

Also Published As

Publication number Publication date
CN112923887B (zh) 2021-08-20

Similar Documents

Publication Publication Date Title
CN104965483B (zh) 一种基于稳健设计的多轴数控机床加工精度保持性优化方法
CN111027010A (zh) 一种钢构件圆柱体拟合算法
CN109492347B (zh) 一种三要素模型描述空间目标光学散射特性的方法
Wang et al. A screw axis identification method for serial robot calibration based on the POE model
US20230384755A1 (en) Datum selection methods and systems for minimizing hole position errors in group hole machining of large components
CN104156519A (zh) 一种面向加工精度可靠度提升的多轴数控机床几何精度设计方法
CN112109084A (zh) 基于机器人关节角度补偿的末端位置补偿方法及其应用
Wang et al. 3D machining allowance analysis method for the large thin-walled aerospace component
CN112643658A (zh) 基于sir降维dh模型的串联机器人自适应误差建模的标定方法
CN112923887B (zh) 任意平面测量控制网非线性最小二乘平差数据检验方法
CN112906198B (zh) 任意平面测量控制网非线性最小二乘平差数据处理方法
CN112906197B (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
CN115752243A (zh) 一种测量数据融合方法
An et al. Geometric properties estimation from line point clouds using Gaussian-weighted discrete derivatives
Tian et al. An improved method for NURBS surface based on particle swarm optimization BP neural network
CN114332173A (zh) 一种基于点云数据的螺旋桨叶片曲面三维重建方法
CN108827300A (zh) 一种基于视觉的装备姿态位置测量方法及系统
Kuang A new approach to the optimal second-order design of geodetic networks
Dong et al. An on-machine precision measurement method for API threads
CN112613196A (zh) 一种五轴机器人姿态逆解方法及相关设备
Lu et al. Error analysis and accuracy optimization for the hand-eye calibration of line structured light
HRDINA et al. DUAL NUMBERS ARITHMETIC IN MULTI AXIS MACHINE ERROR MODELING.
CN115953477B (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