基于同心圆标志的空间平面夹角测量 下载: 1007次
1 引言
角度是生产制造等各领域中的重要物理参量,其测量技术广泛应用于精密制造系统中的高精度测角[1-2]、三维轮廓检测系统中的旋转角测量[3]、木材节点几何重构时纤维角的测量[4],以及卷烟包装盒外形检测[5]等。因此,物体间夹角的实时准确测量具有重要的应用价值。
测角法分为接触式和非接触式两类,接触式测角法由于测量速度快、操作繁琐等原因限制了其应用。早期的非接触式测量是基于自准直仪技术的测角法。自准直仪是一种光学仪器,通过比较目镜十字线位置与经物体反射后的十字线图像间的关系来得出角度值,该方法对被测物表面有镜面反光特性要求,因而应用场合受限[6]。后来,为测量空间大尺寸异面角度,研究人员提出了基于线结构光传播特性的视觉测量方法,该方法多用于大型平台的生产制造,但该方法需要附加结构光投射设备[7]。此外还有拍摄带有标记的被测物图像的视觉测角法,如:在机器人手臂末端贴上靶标,基于单目视觉求解机器人的协同位姿,该方法精度较高,但计算过程复杂[8];在平面粘贴带有圆心的圆形标志,测量平面间的夹角,但这种方法依赖圆形标志的真实圆心[9]。
为实现夹角的实时准确测量,摆脱以往方法过分依赖标志自身属性的特点,并降低图像处理的复杂度,本文提出一种基于同心圆标志的空间平面夹角测量方法。首先通过拍照获取两平面上各组同心圆标志图像,采用阈值分割法结合椭圆检测技术提取椭圆像素点,以拟合曲线方程;然后依据射影变换原理,得到同心圆标志所在平面上的消隐线方程,以此找出两平面法线的消隐点,再结合已知的相机内部参数,最终可计算出两平面间的夹角值。
2 空间平面夹角测量原理
2.1 同心圆标志的物像关系
假设空间中有一正圆曲线,拍照的相机一般采用针孔成像模型[10],成像模型可描述如下:各坐标系中点以相应齐次坐标表示,即以(
式中:
由于拍摄角度等原因,空间正圆曲线经成像模型映射到二维图像平面中的像会变形为椭圆[11-12],产生离心偏差[13],椭圆中心不再是正圆圆心像点[14]。两组同心圆标志的像是4条椭圆曲线,其曲线一般形式可假设为
2.2 消隐线计算原理
相机成像方式是一种射影变换,广泛应用于机器视觉测量中[15]。交比是射影变换中的重要不变量[16],假设有一条直线
式中:
图 1. 物像关系示意图。(a)空间直线与标志交点;(b)直线像与标志像交点
Fig. 1. Diagrams of relationship between object and image. (a) Intersection points of space straight line and marks; (b) intersection points of straight line image and mark images
当空间同心圆
式中:
根据射影变换的交比不变性可知,
由射影几何可知,空间中方向相同的直线有相同的无穷远点,同一平面的多个无穷远点可组成无穷远直线,其像即为消隐线,平面上有且只有一条消隐线[9];在像平面上取
与标记出圆心的圆形标志相比,本文采用的同心圆标志易于图像处理,只需识别出两条成像后的同心圆曲线即可,无需识别圆心;与圆形标志依赖真实圆心像点和曲线之间的关系得到消隐线不同,本文用每条直线与内外圆4个交点间的交比关系求得
2.3 夹角测量的几何原理
已知一平面上的消隐线为
式中:
已知两方向不同的直线的消隐点坐标和相机内部参数
式中:
当
式中:
当
式中:
将(5)式代入(8)式,所得公式对于任意常数变量
由式(9)可解出(
2.4 测量方案及步骤
本文方法无需知道标志的内外圆直径或其比例,仅需标定出相机内部参数即可,整个步骤流程如
拍取图像时需事先在待测两平面上分别贴上同心圆标志,无需标记出圆心;用已知内部参数的相机,在保证可拍摄到标志的前提下,任意拍摄一张图像;再对图像进行阈值分割,提取椭圆曲线的像素点,用最小二乘法拟合曲线方程得到参数
本文方法具有稳健性好、图像处理简单的特点,且整个过程多为线性运算,计算过程较为简洁,对拍照位置无特殊要求,灵活性高,因此可以实现在线实时测量。
3 两平面夹角测量方法实验
本文以MATLAB为平台进行实验,实验分为仿真实验和真实图像实验。
3.1 仿真实验
进行仿真实验时,相机内部参数矩阵
利用基于法向量空间圆曲线方程生成两组同心圆,通过法向量之间的关系可以得到两组同心圆所在平面的夹角,作为理论值。由于内外圆之间的关系对测量精度有一定的影响[17],因此内外半径应保持适当距离。随机生成相机的旋转矩阵
准备工作完成后,拟合出4条椭圆曲线的方程
第一次设置理论值为钝角,
旋转矩阵
图 3. 第一次仿真图。(a)待测物仿真图;(b)待测物上标志对应的像;(c)部分消隐点
Fig. 3. First simulation diagrams.(a) Simulated objects to be measured; (b) images of marks on object to be measured; (c) partial vanishing points
表 1. 第一次仿真实验数据
Table 1. Data of first simulation experiment
|
图 4. 第二次仿真图。(a)待测物仿真图;(b)待测物上标志对应的像;(c)部分消隐点
Fig. 4. First simulation diagrams. (a) Simulated objects to be measured; (b) Images of marks on object to be measured; (c) partial vanishing points
表 2. 第二次仿真实验数据
Table 2. Data of second simulation experiment
|
Note: 10 and 7 is the radius of excircle and inner circle, respectively. 在此次仿真中,拟合消隐线过程选取的直线条数
3.2 真实图像实验
同心圆标志用A4纸制作,虽不需要知道内外圆半径值或其比例关系,但根据前文所述,为了保证精度,适当设置外圆半径约为1 cm,内圆半径约为0.5 cm。实验过程中将标志粘在待测平面上,如
图 5. 实际拍摄图片。(a)(b)待测平面图像;(c)(d)待测平面上标志的像
Fig. 5. Pictures taken practically. (a)(b) Plane images to be measured; (c)(d)images of marks on plane to be measured
实验共测量了6个不同的角度,保证了空间平面夹角的一般性;两平面间的夹角先使用300 mm数字式测量投影仪进行测量,此仪器角度测量分辨率为0.01°,测量结果可作为角度标准值进行对比。然后使用内部参数已经标定出的通用摄像工具拍照,最后用本文方法计算夹角值,结果如
表 3. 实验结果
Table 3. Experimental results
|
从实验测量结果可以看出,测量绝对误差小于1°,相对误差小于1%,精度满足实际应用需求。本文方法多为线性运算,计算速度快,可满足生产制造等领域中夹角实时测量需求。经过分析,测量过程误差产生的原因主要为:标志所在平面不平整,拍照工具本身畸变和拍照获取图像的质量等,且本文采用的拟合方法也会带来一定的误差。
总的来说,与接触式测量方法(如以三坐标测量机为代表的空间平面测角法,其测量原理是通过测得平面上多个点来拟合平面方程,以此得到平面法线,进而得到两平面间的夹角值,测量精度较高,但由于是接触式测量,对环境和操作要求较高)相比,本文方法成本低,操作简便,算法简单,精度可以接受,可用于导轨垂直度检测和机械臂空间精准定位等场合。
4 结论
采用简单易识别的同心圆作为识别标志,图像处理过程简便。根据射影变换原理,提出基于同心圆标志的空间平面夹角测量方法。无需知道内外圆半径或其比例关系,只需标定出相机内部参数矩阵,对拍照位置无要求,灵活性高。实验结果表明,本文方法简单可行,测量绝对误差小于1°,相对误差不超过1%。总体满足生产制造等领域中夹角测量的要求。
本文研究还存在一些不足之处,如拍摄图片中标志曲线缺损,待测面之间的夹角过小,待测夹角表面反光等复杂情况下的测量问题,同心圆标志内外圆半径比例关系对测量结果的影响等。这些问题的研究将在后续工作中展开,以进一步提高本文方法的测量精度,并使其更好地适应复杂的检测条件。
[3] 王小娟, 刘丙才. 基于投影光刀的关节面轮廓旋转检测[J]. 激光与光电子学进展, 2018, 55(9): 091205.
[5] 张勍, 夏营威, 王澍, 等. 基于计算机视觉的卷烟包装方正度测量方法[J]. 烟草科技, 2017, 50(11): 87-92.
[7] 胡文川, 裘祖荣, 张国雄. 基于线结构光的大尺寸空间角度检测算法设计与实现[J]. 光电子·激光, 2013, 24(2): 329-335.
[8] 吕耀宇, 顾营迎, 高瞻宇, 等. 空间协同位姿单目视觉测量系统设计与实验[J]. 激光与光电子学进展, 2017, 54(12): 121505.
[9] 雷经发, 王德麾, 袁中凡. 基于机器视觉的平面夹角测量方法[J]. 华南理工大学学报(自然科学版), 2011, 39(8): 54-59.
[10] 赵小川, 何灏, 唐弘毅. MATLAB计算机视觉实战[M]. 北京: 清华大学出版社, 2018: 173- 187.
Zhao XC, HeH, Tang HY. Computer vision in action with MATLAB[M]. Beijing: Tsinghua University Press, 2018: 173- 187.
[12] 杨忞, 达飞鹏. 基于椭圆外切矩形性质的圆形标志点检测[J]. 光学学报, 2018, 38(12): 1215007.
[13] 安晶晶, 高楠, 陈超, 等. 圆环标识的离心偏差校正方法与评价[J]. 激光与光电子学进展, 2017, 54(4): 041203.
[14] 陈天飞, 赵吉宾, 王银灵, 等. 基于射影变换圆阵靶标中心像点的计算[J]. 仪器仪表学报, 2015, 36(4): 895-902.
[15] 张广军. 机器视觉[M]. 北京: 科学出版社, 2005: 14- 27.
Zhang GJ. Machine vision[M]. Beijing: Science Press, 2005: 14- 27.
[17] 陈新禹, 马孜, 胡英, 等. 视觉测量中圆形标记点的高精度定位[J]. 光电子·激光, 2013, 24(8): 1524-1528.
Article Outline
雷经发, 陈志强, 张淼, 李永玲, 孙虹. 基于同心圆标志的空间平面夹角测量[J]. 激光与光电子学进展, 2019, 56(15): 151504. Jingfa Lei, Zhiqiang Chen, Miao Zhang, Yongling Li, Hong Sun. Measurement of Space Plane Angle Based on Concentric Circle Mark[J]. Laser & Optoelectronics Progress, 2019, 56(15): 151504.