一种基于T-MAC的摄影测量光束法平差优化方法与流程

一种基于t-mac的摄影测量光束法平差优化方法
技术领域
1.本发明涉及一种摄影测量光束法平差方法，特别涉及一种基于t-mac的摄影测量光束法平差优化方法，属于数字化测量技术领域。


背景技术：

2.随着数字测量技术的发展，航空制造业对飞机部件在装配过程中测量的精度、范围和效率提出了更高的要求，因此，在不接触待测部件的情况下，满足此类需求的工业摄影测量技术引起了极大的关注。摄影测量作为一种数字化测量方法，具有高精度、高效率、便携等优点，已在航空工业中广泛应用。随着制造业的自动化发展趋势，摄影测量的研究重点已逐渐向自动化、柔性和拓展性上转移。传统的测量过程中，相机位姿的求解依赖于编码点的匹配，在图像的两两拼接过程中会形成累计误差，降低测量结果精度。同时，标记点的粘贴和清除工作易对被测物表面造成损坏。
3.因此，发明替代编码点拼接的高精度相机定向和测量平差方法，对提高摄影测量精度及自动化程度、保护被测物体有着重大意义。


技术实现要素：

4.本发明的目的在于提供一种基于t-mac的摄影测量光束法平差优化方法，以解决上述背景技术中提出的问题。
5.为实现上述目的，本发明提供如下技术方案：一种基于t-mac的摄影测量光束法平差优化方法，包括以下步骤：
6.s1：选取待测部件的cad模型设计坐标系作为基坐标系，建立各测量坐标系与基坐标系之间的转换关系，标定摄影测量像空间坐标系，将测量结果统一至基坐标系下；
7.s2：控制机器人带动相机，围绕待测部件多角度拍摄，提取相片中像点坐标信息，列出共线方程及对应的误差方程；
8.s3：记录不同摄站下t-mac的位姿参数，通过坐标系标定结果求解基坐标系下的相机外参，建立约束方程；
9.s4：利用光束法平差整体求解空间点三维坐标及改正后的相机外参。
10.作为本发明的一种优选技术方案，所述步骤s1中建立各测量坐标系与基坐标系之间的转换关系，标定摄影测量像空间坐标系包括以下步骤：
11.a1：将t-mac与相机固联在移动机器人末端，则像空间坐标系与t-mac坐标系之间的转换关系是固定的；
12.a2：利用激光跟踪仪测量3个及以上基准点，记基准点在基坐标系下的坐标为p
ibase
，在测量坐标系下坐标为p
imeasurement
，激光跟踪仪测量坐标系与基坐标系之间的转换关系可由转换关系求得；
13.a3：t-mac在测量坐标系下的位姿m
measurement
可根据转换关系计算；
14.a4：记t-mac在基坐标系下位姿参数向量为根据关系式可以求得t-mac在基坐标系下位姿参数；
15.a5：计算摄影测量像空间坐标系相对于基坐标系的转换矩阵
16.作为本发明的一种优选技术方案，所述步骤s2中提取相片中像点坐标信息，列出共线方程及对应的误差方程包括以下步骤：
17.b1：控制机器人带动相机，围绕待测部件多角度拍摄，并提取相片中像点坐标信息；
18.b2：根据共线方程列出误差方程组
[0019][0020]
b3：根据图像处理获取物方点pj在j张相片上同名像点坐标，列出误差方程组：
[0021][0022]
作为本发明的一种优选技术方案，所述步骤s3中记录不同摄站下t-mac的位姿参数，通过坐标系标定结果求解基坐标系下的相机外参，建立约束方程包括以下步骤：
[0023]
c1：记录不同摄站下t-mac的位姿参数，根据s1中坐标系标定结果转换得到基坐标系下的相机外参；
[0024]
c2：根据t-mac测得j张相片的外方位元素，建立约束方程组
[0025][0026]
作为本发明的一种优选技术方案，所述步骤s4中利用光束法平差整体求解空间点三维坐标及改正后的相机外参包括以下步骤：
[0027]
d1：综合s2和s3中误差方程组，将总的误差方程写为v＝ax-l；
[0028]
d2：整体求解待测点的三维空间坐标，其最小二乘解为
[0029][0030]
d3：将解得的各参数改正数值，补偿至误差方程组中，迭代求解直到所有外方位元素和待测点三维坐标的改正数均小于阈值。
[0031]
与现有技术相比，本发明的有益效果是：本发明一种基于t-mac的摄影测量光束法平差优化方法，采用移动机器人搭载工业摄影测量相机对待测部件进行拍摄，t-mac激光跟踪系统追踪相机的位置姿态，替代编码点拼接的方式完成相机定向，将相机外方位元素与待测点三维坐标共同作为平差参数建立约束，改进光束法平差模型，实现待测点的高精度、快速测量。并使用激光投点器在部件表面投射激光点阵作为识别对象，可彻底避免因标记点的粘贴和清除工作对待测部件造成损坏，实现无接触式测量。
附图说明
[0032]
图1本发明的系统示意图；
[0033]
图2本发明的整体流程框图。
具体实施方式
[0034]
对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0035]
请参阅图1-2，本发明提供了一种基于t-mac的摄影测量光束法平差优化方法，步骤1：选取待测部件的cad模型设计坐标系作为基坐标系，建立各测量坐标系与基坐标系之间的转换关系，标定摄影测量像空间坐标系，将测量结果统一至基坐标系下。
[0036]
通过3个及以上基准点建立激光跟踪仪测量坐标系与基坐标系之间的转换关系设基准点在基坐标系下的坐标为p
ibase
，在测量坐标系下坐标为p
imeasurement
，则存在转换关系
[0037][0038]
用m
measurement
表示t-mac在测量坐标系下的位姿，根据转换关系可得其在基坐标系下位姿为
[0039][0040]
像空间坐标系与t-mac坐标系之间的转换关系固定，像空间坐标系相对于基坐标系的转换矩阵为
[0041][0042]
步骤2：控制机器人带动相机，围绕待测部件多角度拍摄，提取相片中像点坐标信息，列出共线方程及对应的误差方程。
[0043]
假设相机无畸变，三维空间内物方点pi处汇集了来自j个摄站的j条光线，列出j个
共线方程：
[0044][0045]
内方位元素(x0,y0,f)认为是已标定过的已知值，将共线方程线性化后得到像点坐标误差方程式为：
[0046][0047]
上式中，v
x
、vy为像点坐标残差，为用近似参数计算的像点坐标。
[0048]
令x1＝[δx
s δy
s δz
s δψ δω δκ]
t
，x2＝[δx δy δz]
t
，系数阵像点坐标改正数方程式的一般式为：
[0049]
v＝a1x1 a2x
2-l
[0050]
根据图像处理获取物方点pj在j张相片上同名像点坐标，列误差方程组为：
[0051][0052]
式中：x1、x2、
……
、xj分别代表j张相片外方位元素的改正数向量，x
j 1
表示待测点三维坐标的改正数向量；a
11
、a
12
、
……
、a
j2
分别为误差方程的系数阵；l1、l2、
……
、lj为常数向量。
[0053]
步骤3：记录不同摄站下t-mac的位姿参数，通过坐标系标定结果求解基坐标系下的相机外参，建立约束方程。
[0054]
设由t-mac测得的第一张相片的外方位元素，经坐标系转换后在基坐标系下的值外方位元素初始估值为则有真值
[0055]
[0056]
其相应的误差方程为
[0057]vj 1
＝x
1-l
j 1
[0058]
以此类推，由t-mac测得j张相片的外方位元素，并转换至基坐标系下，可得误差方程组：
[0059][0060]
步骤4：利用光束法平差整体求解空间点三维坐标及改正后的相机外参。
[0061]
令系数阵(a为2j
×
(j 1)矩阵)。
[0062]
综合步骤3和步骤4中误差方程组，可将总的误差方程写为
[0063]
v＝ax-l
[0064]
把待测点的像坐标及由t-mac测得的相片外方位元素共同视为观测量，整体求解待测点的三维空间坐标，只要有足够的观测值，其最小二乘解为
[0065][0066]
式中，p为和观测值的权，实际计算中常取单位阵；σ
x
为x的协方差阵。
[0067]
将解得的各参数改正数值，补偿至误差方程组中，迭代求解直到所有外方位元素和待测点三维坐标的改正数均小于阈值。
[0068]
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。