0
  • 聊天消息
  • 系统消息
  • 评论与回复
登录后你可以
  • 下载海量资料
  • 学习在线课程
  • 观看技术视频
  • 写文章/发帖/加入社区
会员中心
创作中心

完善资料让更多小伙伴认识你,还能领取20积分哦,立即完善>

3天内不再提示

基于直线空间变换的线结构光光平面标定方法

3D视觉工坊 来源:机器人视觉 作者:机器人视觉 2022-11-14 11:59 次阅读
加入交流群
微信小助手二维码

扫码添加小助手

加入工程师交流群

线结构光非接触式测量技术因测量精度高得到了广泛的应用,而其光平面标定结果会直接影响系统的测量精度。为了得到较高的标定精度,很多专家学者对此进行了研究。现在主要的标定方法有基于平面靶标法、基于平面直线拟合及直线消影原理法,其中关于平面靶标法的研究较为广泛。

研究团队提出了一种基于直线空间变换的线结构光光平面标定方法。该方法根据平面单应性矩阵快速求解靶标面的二维直线方程,并根据靶标坐标系的特殊性构建三维Plücker矩阵形式直线方程,利用直线进行光平面拟合;充分利用整条激光线数据,不需要求解靶标面方程和投影面方程,另外采用直线进行拟合时,利用滑动平均法和正交回归拟合法,保证了前期数据点的抗噪性,增加了标定结果的准确性,从而所提方法满足实际线结构光的测量要求。

一、线结构光视觉测量模型

线结构光测量装置主要由激光投射器和相机组成,其几何结构如下图所示。OwXwYwZw为世界坐标系, OcXcYcZc为相机坐标系, ouv为图像坐标系,点 P为激光线与待测物体交线上任意一点,其在图像平面上对应的投影为点 p。

5614c00e-62a7-11ed-8abf-dac502259ad0.jpg

根据建立的相机成像透视模型,图像坐标系下的点 (u,v)经相机射影变换后与世界坐标系下的点 (xw,yw,zw)的对应关系为

5636b290-62a7-11ed-8abf-dac502259ad0.png

相机坐标系下,线结构光的光平面方程可表示为:

A⋅xc+B⋅yc+C⋅zc+D=0

二、基于直线空间变换的光平面标定方法

2.1 平面光条中心线提取

实际生产环境中,光条成像受外界光源、靶标面不均匀反射影响,存在很多噪声,它们会影响光条中心的提取。采用互相关算法,利用邻域光条纹的相关联性提取光条,即在场景本身深度与物体表面反射特性未发生突变的场合下,条纹灰度的变化是相对平稳的。因此,相邻条纹中心点间的坐标位置与灰度值均具有一定的相似性,利用这种相似性可以在很大的程度上降低噪声的干扰,从而能够提高中心条纹提取精度。以行为单位,使用灰度重心法提取光条中心。另外在光条作用于平面成像为一条直线的前提下,利用五点滑动平均(5MA)法对提取的光条中心点进行预处理,实现过程如下图所示。

5642a104-62a7-11ed-8abf-dac502259ad0.jpg

5MA法取邻近5点均值来消除光条点的随机误差,即由当前点 p=(xm,ym)与上下临近4个点的平均值确定 p'=(x'm,y'm),公式为:

565c919a-62a7-11ed-8abf-dac502259ad0.png

直线拟合常用方法有代数距离和欧氏距离2种。代数距离只考虑光条点集 y与观察值 y˙之间的距离约束,忽略了 x方向的误差,且无法表示竖直方向的直线,据文献,正交回归法具有较好的拟合效果与旋转不变性。为此采用点法式正交回归法拟合光条中心线,以点到直线的欧氏距离建立新约束函数 e:

567642ca-62a7-11ed-8abf-dac502259ad0.png

568bbdf8-62a7-11ed-8abf-dac502259ad0.jpg

2.2 空间直线方程求解

在二维射影变换中,平面直线方程 a⋅u+b⋅v+c=0的齐次表达式为:

56a9d2e8-62a7-11ed-8abf-dac502259ad0.png

满足中心投影映射的两个平面之间存在单应性矩阵 H变换关系,如图3所示,两平面中对应点 x'与 x满足的对应关系为:

56bd51f6-62a7-11ed-8abf-dac502259ad0.png

利用图像上光条中心点与靶标面上点之间的单应性关系和直线方程的齐次表达式,可推导得到两平面上光条中心线单应性变换关系为:

56ce8c82-62a7-11ed-8abf-dac502259ad0.png

如图4所示,在建立的空间直线变换模型中,靶标面上所有点具有特殊性,其在靶标空间坐标系中满足zw=0,为此可将原先的二维直线方程上的点全部拓展至三维。选取光条上任意两点 A=(xa,ya,za,1)T,B=(xb,yb,zb,1)T,其中za=zb=0,将空间直线的方程组解析式转换成Plücker矩阵形式,便于矩阵运算,直线的Plücker表达形式 Lw不受取点大小影响。Lw的表达式为:

56dc7bc6-62a7-11ed-8abf-dac502259ad0.png

靶标与相机坐标系的转换关系可用相机成像模型(8)式求解得到:

56ed6f6c-62a7-11ed-8abf-dac502259ad0.png

图像上圆形标记点与靶标坐标系下的圆形标记点之间存在N点透视关系PnP(perspective N-point),因此可建立两者之间的重投影误差函数,使用Levenberg-Marquardt算法优化重投影误差函数,获得非线性最优解 R和 T。根据坐标系转换,得到了两坐标系下对应点的关系,表达式为:

571eb018-62a7-11ed-8abf-dac502259ad0.png

进一步可推出空间直线变换对应关系为:

5733a3ba-62a7-11ed-8abf-dac502259ad0.png

5746140a-62a7-11ed-8abf-dac502259ad0.jpg

2.3 光平面方程计算

激光线 Lc为光平面 π与靶标平面交线,其与光平面共面,根据线面约束方程,可知:

Lc⋅π=0

多次移动靶标平面( n≥2),得到多组光平面上的直线 Lci,建立线面方程组:

LT⋅π=0,L=(Lc1,Lc2,…,Lcn)

由于参与拟合光条直线方程的数据存在提取误差,可增加激光线数量,将多个线面约束方程求解问题转换成超静定齐次方程组的最小二乘解问题,利用SVD分解 L(L=UDVT),其中 U为酉矩阵, D为奇异值沿对角线从大到小构成的对角矩阵, V为酉矩阵。D中最小奇异值在 VT对应的特征向量为平面方程的系数 π,从而完成光平面方程参数求解。

三、试验结果

57759acc-62a7-11ed-8abf-dac502259ad0.jpg

标准台阶块实验。(a)台阶块及特征点;(b)台阶点云;(c)截面形状

57a03106-62a7-11ed-8abf-dac502259ad0.jpg

圆柱实验。(a)圆柱样件;(b)拟合圆

作者:翟鹏, 崔海华, 胡广露, 张益华, 靳宇婷, 黄怡.

审核编辑:郭婷


声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉
  • 激光
    +关注

    关注

    21

    文章

    3578

    浏览量

    69090
  • 相机
    +关注

    关注

    5

    文章

    1535

    浏览量

    55410

原文标题:视觉 | 面向线结构光测量的直线空间变换光平面标定方法

文章出处:【微信号:3D视觉工坊,微信公众号:3D视觉工坊】欢迎添加关注!文章转载请注明出处。

收藏 人收藏
加入交流群
微信小助手二维码

扫码添加小助手

加入工程师交流群

    评论

    相关推荐
    热点推荐

    IMU+多相机高速联合自动标定方案

    “IMU+多相机高速联合自动标定仪”,20 秒完成全流程,一举填补国产高端标定装备市场断层。 机械臂全程抓取位姿,20秒完成多相机+IMU联合标定,人工零干预,产线节拍提升30-10
    发表于 10-23 14:04

    直线电机平台的优势解析

    与专业实力 作为专注非标定直线电机平台及精密传动控制解决方案达13年的厂家,积累了丰富的行业洞见与项目经验。服务超过500家来自不同领域的客户,使其深刻理解多样化的应用需求。汉诺的核心竞争力在于其强大的非标定制化设计
    的头像 发表于 08-04 18:19 602次阅读

    ZG02-D/SG系列 在线动态拟合全长直线度测量仪

    设备主要用于金属棒材、管材、陶瓷管材等产品需要在线快速测量直径和直线度检测,设备采用单台采集模块,占用空间小,适用于安装空间相对较小的工况,可实时显示产品通过后的整体直线度数据,搭配测
    发表于 08-04 15:09

    亳州ZG02-D/SY系列 离线扫描直线度测量仪

    、每米等多种直线度误差和外径。 1、基本原理 两组测头可组成一个坐标系,通过每个测头测量的被测物中心位置可计算出被测物中心在坐标系中的坐标位置 2、数据计算方法 测量仪通过标定建立空间
    发表于 07-31 15:33

    日本直线电机地铁系统的发展与改进

    地铁特别适用于地形较复杂或线网较为成熟的大城市。对于重庆这种地理条件特殊、地下空间受限、人口又比较密集的城市,尤其适合修建直线电机地铁。如果新建线与普通地铁线有多个换乘站,选择直线电机
    发表于 06-12 13:54

    JCMsuite应用:一维周期线

    和P偏振平面波照亮。JCMsuite计算近场分布。下图显示了当波长为193nm时,平面波从衬底侧垂直入射到结构内的近场强度。 S偏振光照明的近场强度 P偏振光照明的近场强度 后处理傅里叶
    发表于 06-10 08:48

    PanDao应用:输入中频空间公差,最小中频波长

    半径为6400公里的行星上的中频空间波。 条纹的非直边是由未优化设置的CNC抛光工艺中的中频空间误差引起的 预抛光光学元件上的同心圆状中频空间误差,由未优化的点接触式CNC研磨工
    发表于 06-04 08:46

    JCMsuite应用:孤立线

    斜入射S和P偏振平面波照亮。JCMsuite计算近场分布。下图显示了当波长为193nm时,平面波从衬底侧垂直入射到结构内的近场强度 S偏振光照明的近场强度 P偏振光照明的近场强度 后处理傅里叶
    发表于 05-30 08:48

    一文告诉你为什么不要随便在高速线旁边铺铜!

    传输线结构:带状线、微带线和共面波导。带状线是嵌入在两个参考平面之间的信号
    发表于 04-07 10:52

    永磁直线同步电机建模、特性分析及推力控制

    旋转运动和直线运动是两种最基本的运动形式。1831年世界上出现了第 一台旋转电机,旋转电机传动方式是通过中间转换装置将旋转运动变换直线 运动。1845英国人 Wheatstone研制出了世界上
    发表于 03-31 15:45

    一种分段气隙的CLLC变换平面变压器设计

    、应用场景等方法变换器进行了研究,由于天然的ZVS和ZCS备受行业和学者的青睐。就磁性器件方面,传统的磁性器件设计方法不能再满足对功率密度和性能的需求,适用于高频的平面变压器随着出现
    发表于 03-27 13:57

    测出直线度数据后 如何评估直线度误差的大小、分布和趋势?

    测物体上直线度误差的最大程度。 二、评估直线度误差的分布 直线度误差的分布描述了误差值在被测物体上的空间分布情况。评估直线度误差的分布时,可
    发表于 02-05 16:35

    全长直线度和米直线度如何测量?

    的高平行度和干涉原理,可以非常精确地测量直线度。适用于高精度要求的直线度检测,如机床导轨、光学平面的检测。 三坐标测量机(CMM):可以进行三维空间
    发表于 01-16 14:19

    JCMsuite应用:倾斜平面波传播透过光阑的传输

    这个例子计算入射平面波(在倾斜的入射角度)通过一个独立光阑的传播: 光阑几何结构 下图显示了计算得到的近场的矢量图。 近场的矢量图 几何图形在xy横截方向上具有透明的边界条件。因此,对于周期
    发表于 12-27 08:51

    直线模组安装调试与要点

    直线模组的安装位置需根据实际应用要求和机械装置的结构来选择,直线模组可以被安装在机械装置的各个部位中,并且可以与其他部件配合使用,而直线模组的安装调试
    的头像 发表于 12-10 17:52 1141次阅读
    <b class='flag-5'>直线</b>模组安装调试与要点