高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

成形砂轮激光修整的多轮廓图像合成检测方法

刘宣呈 陈根余 操坤 曹明月 梅枫

引用本文:
Citation:

成形砂轮激光修整的多轮廓图像合成检测方法

    通讯作者: 陈根余, hdgychen@163.com
  • 中图分类号: TP391;TG665;TG74+3

Multi-profile images synthesis detection method for laser dressing of forming grinding wheels

    Corresponding author: CHEN Genyu, hdgychen@163.com ;
  • CLC number: TP391;TG665;TG74+3

  • 摘要: 为了解决复杂轮廓成形砂轮激光修整过程中在线检测难、轮廓检测精度低、速度慢等问题, 搭建了一套基于砂轮激光修整机床的图像采集系统, 提出了一种多轮廓图像合成检测方法, 定义了轮廓波动带的概念。首先利用引导滤波和灰度变换预处理砂轮轮廓图像, 增强轮廓特征; 再使用全局阈值分割和反二值化得到二值图像, 将二值图像序列合成为带有像素频率信息的轮廓频率图像, 通过阈值分割提取出轮廓波动带; 最后使用随机抽样一致算法结合最小二乘法对波动带最外侧的边缘点集进行拟合, 完成直线和圆弧特征的检测。结果表明, 该检测方法与单一截面视觉检测法相比, 检测效率接近, 但误差范围缩小了66.67%;与磨削复映检测法相比, 检测结果存在一个固定偏差, 平均误差在7 μm左右, 但检测效率提高了91.67%。该方法能高精度、高效地检测成形砂轮轮廓, 为成形砂轮轮廓检测方法提供了新思路。
  • 图 1  视觉装夹运动部分

    Figure 1.  Visual clamping motion section

    图 2  检测系统整体架构

    Figure 2.  Overall architecture of the detection system

    图 3  轮廓检测方法

    Figure 3.  Profile detection method

    图 4  检测系统关键算法

    Figure 4.  Key algorithms for detection systems

    图 5  砂轮图像直方图

    Figure 5.  Histogram of grinding wheel images

    图 6  分段线性变换曲线

    Figure 6.  Segmented linear transformation curve

    图 7  图像预处理效果

    Figure 7.  Images pre-processing effect

    图 8  轮廓合成算法流程图

    Figure 8.  Profile synthesis algorithm flow chart

    图 9  阈值分割和反二值化效果

    Figure 9.  Effects of thresholding segmentation and inverse binarization

    图 10  轮廓图像中杂质干扰(放大)

    Figure 10.  Impurity interference in profile images (zoom)

    图 11  两种轮廓图像合成场景

    Figure 11.  Two kinds of profile images synthesis scenes

    图 12  轮廓图像合成算法效果

    Figure 12.  Effect of profile images synthesis algorithm

    图 13  频率图像局部放大

    Figure 13.  Frequency image partial zoom

    图 14  轮廓波动带提取效果

    Figure 14.  Profile fluctuation with extraction effect

    图 15  拟合效果对比

    Figure 15.  Comparison of fitting effects

    图 16  凹V型轮廓成形砂轮实物及尺寸情况

    Figure 16.  Concave V-slot forming grinding wheel and dimensional situation

    图 17  磨削复映装置

    Figure 17.  Grinding remapping device

    图 18  石墨板

    Figure 18.  Graphite sheet

    图 19  在线检测平台

    Figure 19.  Online detection platform

    图 20  影像测量软件主界面

    Figure 20.  Image measurement software main interface

    图 21  3种测量方式检测过程

    Figure 21.  Three measurement methods detection process

    表 1  V型槽角度D检测结果

    Table 1.  Detection results of V-slot angle D

    the number of V-slot the multi-profile images synthesis detection method/mm the single-section visual detection method/mm the grinding remapping detection method/mm
    1 44.213 44.156 44.281
    2 44.087 43.915 44.042
    3 44.245 44.323 44.203
    4 44.076 44.031 44.109
    5 44.070 44.264 44.115
    6 44.208 44.344 44.293
    7 44.168 44.199 44.156
    8 44.092 43.984 44.067
    下载: 导出CSV

    表 2  V型槽深度H检测结果

    Table 2.  Detection results of V-slot depth H

    the number of V-slot the multi-profile images synthesis detection method/mm the single-section visual detection method/mm the grinding remapping detection method/mm
    1 0.5251 0.5274 0.5183
    2 0.5294 0.5396 0.5206
    3 0.5342 0.5362 0.5299
    4 0.5297 0.5373 0.5222
    5 0.5300 0.5450 0.5255
    6 0.525 0.5338 0.5177
    7 0.5229 0.5287 0.5197
    8 0.5281 0.5471 0.5207
    下载: 导出CSV

    表 3  槽底圆弧半径R检测结果

    Table 3.  Detection results of V-slot bottom arc radius R

    the number of V-slot the multi-profile images synthesis detection method/mm the single-section visual detection method/mm the grinding remapping detection method/mm
    1 0.1664 0.1611 0.1749
    2 0.1624 0.1492 0.1686
    3 0.1586 0.161 0.1672
    4 0.1621 0.1566 0.1729
    5 0.1597 0.1584 0.1671
    6 0.1638 0.1571 0.1717
    7 0.1662 0.1623 0.1706
    8 0.1642 0.1528 0.1712
    下载: 导出CSV

    表 4  检测效率

    Table 4.  Detection efficiency

    the multi-profile images synthesis detection method the single-section visual detection method the grinding remapping detection method
    time spent/s 30 20 360
    下载: 导出CSV

    表 5  检测误差

    Table 5.  Detection errors

    detection features detection methods the number of V-slot
    1 2 3 4 5 6 7 8
    angel/(°) the single-section visual detection method -0.265 -0.157 0.16 0.122 0.169 0.071 0.043 -0.203
    the multi-profile images synthesis detection method -0.082 0.065 0.091 -0.033 -0.105 -0.085 0.012 0.025
    depth/mm the single-section visual detection method 0.0091 0.019 0.0063 0.0151 0.0195 0.0161 0.009 0.0264
    the multi-profile images synthesis detection method 0.0068 0.0088 0.0043 0.0075 0.0045 0.0073 0.0032 0.0074
    radius/mm the single-section visual detection method -0.0138 -0.0194 -0.0062 -0.0163 -0.0087 -0.0146 -0.0083 -0.0184
    the multi-profile images synthesis detection method -0.0085 -0.0062 -0.0086 -0.0108 -0.0074 -0.0079 -0.0044 -0.007
    下载: 导出CSV

    表 6  检测误差的统计指标

    Table 6.  Statistical indicators of the detection errors

    detection features statistical indicators the multi-profile images synthesis detection method the single- section visual detection method
    mean value/mm -0.0076 -0.0132
    radius range/mm 0.0064 0.0132
    variance 3.0925×10-6 2.13×10-5
    mean value/mm 0.0062 0.01506
    depth range/mm 0.0056 0.0201
    variance 3.36937×10-6 3.92×10-5
    下载: 导出CSV
  • [1] 王龙, 汪刘应, 刘顾, 等. 成型砂轮磨削齿轮的运动学模型构建与齿轮精度分析[J]. 制造技术与机床, 2019(12): 112-116.

    WANG L, WANG L Y, LIU G, et al. Kinematic model construction and gear precision analysis of gear machined by form grinding[J]. Manufacturing Technology & Machine Tool, 2019(12): 112-116(in Chinese). 
    [2] 王昊. 金刚石成形砂轮视觉检测精度分析与视觉系统改进研究[D]. 长沙: 湖南大学, 2020: 1-55.

    WANG H. Precision analysis of diamond wheel visual detection and improvement research of visual system[D]. Changsha: Hunan University, 2020: 1-55(in Chinese).
    [3] 邓辉, 朱鹏程, 王林青, 等. 砂轮表面测量方法研究进展[J]. 机械科学与技术, 2020, 39(4): 599-607.

    DENG H, ZHU P Ch, WANG L Q, et al. Research progress in measurement methods of grinding wheel surface[J]. Mechanical Science and Technology for Aerospace Engineering, 2020, 39(4): 599-607(in Chinese). 
    [4] 彭檐波, 陈根余, 周聪, 等. 基于图像处理的成形砂轮激光切向整形系统算法与实现[J]. 应用激光, 2018, 38(3): 454-460.

    PENG Y B, CHEN G Y, ZHOU C, et al. Algorithm and implementation of laser tangential shaping grinding wheel system based on image processing[J]. Applied Laser, 2018, 38(3): 454-460(in Chin-ese). 
    [5] 彭檐波. 基于机器视觉的激光切向整形V形凹面砂轮测控系统设计[D]. 长沙: 湖南大学, 2018: 1-67.

    PENG Y B. Design of measuring and controlling system for laser tangential shaping of V-shaped concave grinding wheels based on machine vision[D]. Changsha: Hunan University, 2018: 1-67(in Chin-ese).
    [6] 类志杰. 基于机器视觉的凸圆弧砂轮边缘检测及定位技术研究[D]. 长沙: 湖南大学, 2019: 1-66.

    LEI Zh J. Research on edge detection and positioning technology of convex arc grinding wheel based on machine vision[D]. Changsha: Hunan University, 2019: 1-66(in Chinese).
    [7] 牛牧, 许黎明, 赵达, 等. 基于工件轮廓图像的砂轮磨损在线检测方法[J]. 上海交通大学学报, 2021, 55(3): 221-228.

    NIU M, XU L M, ZHAO D, et al. Workpiece-contour-image based wheel wear online detection methodology[J]. Journal of Shanghai Jiao Tong University, 2021, 55(3): 221-228(in Chinese). 
    [8] 牛牧. 曲线磨削砂轮修整方法和修整工艺的研究[D]. 上海: 上海交通大学, 2020: 1-96.

    NIU M. Research on dressing method and dressing process of profile grinding wheel[D]. Shanghai: Shanghai Jiao Tong University, 2020: 1-96(in Chinese).
    [9] 母德强, 郑金涛, 司苏美, 等. 基于激光三角法的砂轮磨损量检测[J]. 长春工业大学学报, 2021, 42(3): 200-204.

    MU D Q, ZHENG J T, SI S M, et al. Research on wear measurement of grinding wheel based on laser triangulation[J]. Journal of Changchun University of Technology, 2021, 42(3): 200-204 (in Chin-ese). 
    [10] 郑金涛. 磨床砂轮磨损量与表面轮廓检测技术研究[D]. 长春: 长春工业大学, 2021: 1-47.

    ZHENG J T. Research on measurement technology of grinding wheel wear and surface profile for grinding machine[D]. Changchun: Changchun University of Technology, 2021: 1-47(in Chinese).
    [11] 周炼, 安晨辉, 侯晶, 等. 圆弧金刚石砂轮三维几何形貌的在位检测和误差评价[J]. 光学精密工程, 2017, 25(12): 3079-3088.

    ZHOU L, AN Ch H, HOU J, et al. 3D geometric topographic mea-surement and error evaluation of arc diamond wheel[J]. Optics and Precision Engineering, 2017, 25(12): 3079-3088 (in Chinese). 
    [12] 师超钰, 朱建辉, 孙冠男, 等. 成形磨砂轮圆弧廓形关键参数在位检测方法及试验研究[J]. 中国测试, 2021, 47(3): 36-42.

    SHI Ch Y, ZHU J H, SUN G N, et al. In-situ measuring method and experimental research on key circular-arc profile parameters of form grinding wheel[J]. China Measurement & Test, 2021, 47(3): 36-42(in Chinese). 
    [13] 邵升阳. 基于数字图像处理的砂轮磨损检测系统及方法研究[D]. 厦门: 厦门大学, 2020: 1-62.

    SHAO Sh Y. Research on grinding wheel wear detection system and method based on digital image processing[D]. Xiamen: Xiamen University, 2020: 1-62(in Chinese).
    [14]

    XU L M, XU K Z, CHAI Y D. Identification of grinding wheel wear signature by a wavelet packet decomposition method[J]. Journal of Shanghai Jiao Tong University (Science Edition), 2010, 15(3): 323-328.
    [15] 胡一星. 基于全局图像的轮廓曲线磨削加工误差动态检测与补偿[D]. 上海: 上海交通大学, 2018: 1-57.

    HU Y X. Global image based online error detection and compensation of profile grinding[D]. Shanghai: Shanghai Jiao Tong University, 2018: 1-57(in Chinese).
    [16] 胡一星, 范帆, 许黎明, 等. 砂轮廓形原位视觉检测方法与试验研究[J]. 机械制造, 2018, 56(5): 80-82.

    HU Y X, FAN F, XU L M, et al. In-situ visual detection method and experimental research of the abrasive profile[J]. Machinery, 2018, 56(5): 80-82 (in Chinese). 
    [17]

    LI T Y, QIU Z J, TANG J J. Research on measurement method of grinding wheel profile based on image mosaic[J]. Measurement Science and Technology, 2019, 31(3): 035402.
    [18] 李天一. 基于机器视觉的砂轮轮廓测量方法的研究[D]. 天津: 天津大学, 2019: 1-76.

    LI T Y. Research on grinding wheel profile measurement based on machine vision[D]. Tianjin: Tianjin University, 2019: 1-76(in Chinese).
    [19] 郭佩瑜, 张宝华. 基于引导滤波和模糊算法的红外背景抑制算法[J]. 激光技术, 2018, 42(6): 854-858.

    GUO P Y, ZHANG B H. Infrared background suppression algorithm based on guided filtering and fuzzy algorithm[J]. Laser Technology, 2018, 42(6): 854-858(in Chinese). 
    [20] 王浩, 张叶, 沈宏海, 等. 图像增强算法综述[J]. 中国光学, 2017, 10(4): 438-448.

    WANG H, ZHANG Y, SHEN H H, et al. Review of image enhancement algorithms[J]. Chinese Optics, 2017, 10(4): 438-448(in Chinese). 
    [21] 蔡旭明, 李笑, 刘玉县, 等. 基于灰度直方图的激光光斑中心定位算法[J]. 激光技术, 2023, 47(2): 273-279.

    CAI X M, LI X, LIU Y X, et al. Laser spot center location algorithm based on gray histogram[J]. Laser Technology, 2023, 47(2): 273-279 (in Chinese). 
    [22] 吴一全, 孟天亮, 吴诗婳. 图像阈值分割方法研究进展20年(1994—2014)[J]. 数据采集与处理, 2015, 30(1): 1-23.

    WU Y Q, MENG T L, WU Sh H. Research progress of image thresholding methods in recent 20 years(1994-2014)[J]. Journal of Data Acquisition and Processing, 2015, 30(1): 1-23(in Chinese). 
    [23] 张伟杰, 刘立君, 张红兴. 基于图像形态学激光模具裂纹修复技术研究[J]. 激光技术, 2016, 40(2): 189-194.

    ZHNAG W J, LIU L J, ZHANG H X. Repair techniques of dies with laser based on image morphological processing[J]. Laser Technology, 2016, 40(2): 189-194 (in Chinese). 
    [24] 袁清珂, 张振亚, 毕庆. 改进RANSAC算法在直线拟合中的应用[J]. 组合机床与自动化加工技术, 2015(1): 123-125.

    YUAN Q K, ZHANG Zh Y, BI Q. Linear fitting application based on the improved RANSAC algorithm[J]. Modular Machine Tool & Automatic Manufacturing Technique, 2015(1): 123-125(in Ch-inese). 
  • [1] 冯煦张瑞瑛周萍李松 . 大功率半导体线激光图像处理方法研究. 激光技术, 2010, 34(5): 624-627. doi: 10.3969/j.issn.1001-3806.2010.O5.013
    [2] 郭佩瑜张宝华 . 基于引导滤波和模糊算法的红外背景抑制算法. 激光技术, 2018, 42(6): 854-858. doi: 10.7510/jgjs.issn.1001-3806.2018.06.024
    [3] 陈勤霞武文成艾斯卡尔·艾木都拉 . 一种基于多尺度点状目标建模的检测算法. 激光技术, 2020, 44(4): 520-524. doi: 10.7510/jgjs.issn.1001-3806.2020.04.021
    [4] 蔡旭明李笑刘玉县何春华林俊杰 . 基于灰度直方图的激光光斑中心定位算法. 激光技术, 2023, 47(2): 273-279. doi: 10.7510/jgjs.issn.1001-3806.2023.02.018
    [5] 闵永智王红霞康飞党建武 . 基于图像式传感器的铁路轨距检测系统研究. 激光技术, 2015, 39(3): 344-348. doi: 10.7510/jgjs.issn.1001-3806.2015.03.014
    [6] 李文龙戈海龙任远成巍 . 图像处理技术在激光熔池温度检测的应用. 激光技术, 2018, 42(5): 599-604. doi: 10.7510/jgjs.issn.1001-3806.2018.05.004
    [7] 曹国强刘禹廷王琳霖 . 改进非下采样轮廓波在散斑条纹中的滤波处理. 激光技术, 2019, 43(2): 269-274. doi: 10.7510/jgjs.issn.1001-3806.2019.02.022
    [8] 吴家洲刘君施佳文张胜 . 激光焊缝图像分割与颜色识别方法研究. 激光技术, 2023, 47(5): 723-728. doi: 10.7510/jgjs.issn.1001-3806.2023.05.022
    [9] 未永飞杜正春姚振强 . 中值滤波在激光雷达点云数据预处理中的应用. 激光技术, 2009, 33(2): 213-216.
    [10] 崔治邓曙光肖卫初 . 利用HSSIM和残差比阈值的3维激光扫描图像去噪. 激光技术, 2015, 39(5): 669-673. doi: 10.7510/jgjs.issn.1001-3806.2015.05.018
    [11] 张海庄姚梅雷萍李鹏曾庆平 . 远场激光光斑图像处理方法研究. 激光技术, 2013, 37(4): 460-463. doi: 10.7510/jgjs.issn.1001-3806.2013.04.010
    [12] 顾国庆王开福燕新九 . 基于同态滤波的电子散斑干涉图像处理. 激光技术, 2010, 34(6): 750-752,797. doi: 10.3969/j.issn.1001-3806.2010.06.009
    [13] 张宝华刘鹤 . 基于区域定位与轮廓分割的红外目标检测. 激光技术, 2015, 39(6): 840-844. doi: 10.7510/jgjs.issn.1001-3806.2015.06.023
    [14] 刘逸飞苏亚姚晓天崔省伟杨丽君周聪聪何松 . OCT无创血糖检测图像处理最优化方法研究. 激光技术, 2023, 47(2): 178-184. doi: 10.7510/jgjs.issn.1001-3806.2023.02.004
    [15] 汤敏王惠南 . 激光扫描共聚焦显微镜图像的计算机处理. 激光技术, 2007, 31(5): 558-560.
    [16] 苏平牛燕雄李大乾牛海莎李易难张超 . 基于面阵CCD的激光告警系统的图像采集与处理. 激光技术, 2013, 37(3): 394-399. doi: 10.7510/jgjs.issn.1001-3806.2013.03.028
    [17] 孟宇帆张丽君何长涛肖婧阳宁静冯国英韩敬华 . 基于图像处理的激光清洗飞机蒙皮特性和机制研究. 激光技术, 2024, 48(3): 303-311. doi: 10.7510/jgjs.issn.1001-3806.2024.03.002
    [18] 张怡霄杜惊雷高福华姚军曾阳素郭永康 . 分数域啁啾滤波及其在数字图像处理中的应用. 激光技术, 2003, 27(1): 78-80.
    [19] 李昌海叶玉堂沈淦松徐伟叶涵姚景昭 . 基于图像轮廓分析的LCD线路缺陷检测. 激光技术, 2013, 37(2): 207-210. doi: 10.7510/jgjs.issn.1001-3806.2013.02.017
    [20] 李旭寒董安国封建湖 . 基于多级引导滤波器的图像区域融合算法. 激光技术, 2016, 40(5): 756-761. doi: 10.7510/jgjs.issn.1001-3806.2016.05.029
  • 加载中
图(21) / 表(6)
计量
  • 文章访问数:  1950
  • HTML全文浏览量:  1033
  • PDF下载量:  13
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-06-05
  • 录用日期:  2023-07-06
  • 刊出日期:  2024-05-25

成形砂轮激光修整的多轮廓图像合成检测方法

    通讯作者: 陈根余, hdgychen@163.com
  • 1. 湖南大学 机械与运载工程学院, 长沙 410082, 中国
  • 2. 湖南大学 激光研究所, 长沙 410082, 中国

摘要: 为了解决复杂轮廓成形砂轮激光修整过程中在线检测难、轮廓检测精度低、速度慢等问题, 搭建了一套基于砂轮激光修整机床的图像采集系统, 提出了一种多轮廓图像合成检测方法, 定义了轮廓波动带的概念。首先利用引导滤波和灰度变换预处理砂轮轮廓图像, 增强轮廓特征; 再使用全局阈值分割和反二值化得到二值图像, 将二值图像序列合成为带有像素频率信息的轮廓频率图像, 通过阈值分割提取出轮廓波动带; 最后使用随机抽样一致算法结合最小二乘法对波动带最外侧的边缘点集进行拟合, 完成直线和圆弧特征的检测。结果表明, 该检测方法与单一截面视觉检测法相比, 检测效率接近, 但误差范围缩小了66.67%;与磨削复映检测法相比, 检测结果存在一个固定偏差, 平均误差在7 μm左右, 但检测效率提高了91.67%。该方法能高精度、高效地检测成形砂轮轮廓, 为成形砂轮轮廓检测方法提供了新思路。

English Abstract

    • 随着成形砂轮大量应用在高精度脆硬材料零件加工中,比如半导体、陶瓷、复合材料等材料的加工,砂轮轮廓精度的要求日益严格。通过检测成形砂轮轮廓,能够获取成形砂轮的轮廓尺寸,检查其精度是否满足要求; 同时为实际磨削加工中,工件型面的误差分析和补偿提供依据[1]; 另外还可以了解砂轮的轮廓磨损情况,定位轮廓发生磨损的位置,确定轮廓待修整余量,指导砂轮修整[2]。作为成形砂轮生产和修整环节中的一项关键技术,研究成形砂轮轮廓检测方法有重要作用和意义。

      目前常用的砂轮轮廓检测方法有磨削复映检测法、接触式检测法、传统光学投影检测法和单截面视觉检测法等[3]。PENG等人[4-5]设计了一套激光整形测控系统,结合磨削复映和机器视觉技术检测砂轮的轮廓,并根据检测结果对砂轮补偿整形。LEI等人[6]搭建了一套砂轮轮廓视觉检测系统,实现轮廓尺寸的在机检测,同时分析修整偏差并改善修整工艺以提高修整质量。NIU等人[7-8]设计了轮廓在线检测系统实时采集被磨削工件的图像,提取并测量实际轮廓,同时与理论轮廓对比,判断砂轮磨损状态和计算磨损量。MU等人[9-10]使用基恩士LK-G30激光位移传感器搭建一套实验检测系统,可以非接触地检测砂轮的磨损量,整周期测量时,其线性度小于0.6%。ZHOU等人[11]使用基恩士LK-H020激光位移传感器螺旋扫描砂轮的外轮廓获取轮廓数据并进行滤波和插值得到圆弧轮廓的点云数据,实现凸圆弧金刚石砂轮的轮廓检测和评价。SHI等人[12]使用2维激光位移传感器采集砂轮轮廓数据,建立点云模型,获取砂轮廓形的圆弧半径和圆度误差,正确反映成形砂轮的磨削工件的精度。SHAO等人[13]搭建了一种砂轮形貌采集系统,提取高精度的磨损平面,评价砂轮的磨损状态和磨削性能。XU等人[14-16]在曲线磨床上设计了在位测量系统,通过电荷耦合器件(charge-coupled device,CCD)相机采集轮廓图像,使用轮廓拟合测量出轮廓尺寸数据,测量结果与磨削复映法结果接近。LI等人[17-18]搭建了砂轮轮廓测量系统,进一步研究了自动对焦和图像拼接等算法,并通过检测圆柱形和半球形轮廓的砂轮进行验证,可实现砂轮轮廓的高效测量以及轮廓参数的评价。

      现有的检测方法只能满足精度和效率中的单方面需求,并主要针对于凸圆弧砂轮,对于其它轮廓的成形砂轮检测研究相对较少,且只测量某一截面轮廓,无法反映砂轮的全局轮廓信息。

      为同时满足成形砂轮激光修整过程中高效、高精度的检测需求,本文作者对成形砂轮轮廓的检测技术开展研究。搭建基于砂轮激光修整机床的图像采集系统对砂轮轮廓进行旋转采样,获取砂轮圆周方向的全局轮廓图像;提出一种多轮廓图像合成检测方法,有效利用砂轮全圆周的轮廓信息, 进一步提出轮廓波动带概念,通过轮廓波动带对砂轮进行轮廓尺寸检测。为了验证该方法的检测效果,将该本文中的检测方法与磨削复印检测法和单截面视觉检测法进行了对比分析。

    • 检测系统主要由图像采集部分、砂轮激光修整机视觉装夹运动部分、数控系统通讯和工控机及软件算法构成。图像采集部分由工业相机、工业相机镜头和照明光源三大组件组成。工业相机选用的是MV-CH120-10TM,具有1200万像素,负责采集和传输图像;镜头选用的是WKLT-40-100双侧远心镜头,负责物体成像;光源选用的是MHCL-65H-65-B远心平行光源和配套控制器,负责提升图像质量。视觉装夹运动部分负责将图像采集部分的硬件安装到砂轮激光修整机上,同时实现对焦和增加测量范围的功能,如图 1所示。工控机负责图像采集控制、机床通讯、检测软件运行等。检测系统示意图整体架构设计如图 2所示。

      图  1  视觉装夹运动部分

      Figure 1.  Visual clamping motion section

      图  2  检测系统整体架构

      Figure 2.  Overall architecture of the detection system

    • 多轮廓合成轮廓检测方法的技术路线如图 3所示。待测砂轮安装在修整机主轴上,视觉工控机与修整机数控系统通讯软件读取主轴角度,主轴每转过一定角度,软件控制相机采集一张照片,直到砂轮旋转一周,得到砂轮完整的圆周方向的轮廓图像序列。对轮廓图像序列进行预处理,增强轮廓特征并生成二值图像。再将二值图像序列合成为带像素频率信息的轮廓频率图像,在轮廓频率图像中选取一定灰度区间中的像素,提取出轮廓波动带,最后基于轮廓波动带, 实现轮廓尺寸测量。

      图  3  轮廓检测方法

      Figure 3.  Profile detection method

    • 本文中的轮廓检测算法需要基于圆周方向多个截面的轮廓图像,多轮廓图像合成算法是其核心,合成后的频率图像是后续检测和试验的基础。检测算法整体设计如图 4所示。

      图  4  检测系统关键算法

      Figure 4.  Key algorithms for detection systems

    • 引导滤波是近年来流行的新的保边滤波算法,能在去除大量噪点的同时保留砂轮轮廓边缘细节,避免梯度反转现象,提高检测精度,它假设图像内存在局部的线性关系,认为图像是非解析的2维函数,但在局部可以用解析的线性函数来表达[19],如下式所示:

      $ q_i=a_k I_i+b_k, \left(i \in Q_k\right) $

      (1)

      式中: qi表示滤波输出图像q的值;Ii为引导图像I (可以是原始图像)的值;Qk为第k个局部窗口;akbk在第k个局部窗口的常量系数。

      对式(1)求导可得∇qi=akIi,如果引导图像中出现边缘,则输出图像q中也存在边缘,故可以实现边缘保留。求解式(1)中akbk的值,为了减少失真,使用均方误差代价函数M(ak, bk) 进行优化,如下式所示:

      $ M\left(a_k, b_k\right)=\sum\limits_{i \in Q_k}\left[\left(a_k I_i+b_k-p_i\right)^2+\varepsilon a_k^2\right] $

      (2)

      式中: ε为补偿系数;pi为原始图像p的值。

      转化为最优化问题,即使式(2)的值最小,使用最小二乘法可计算akbk的值分别为:

      $ a_k=\frac{\sum\limits_{i \in Q_k}\left(I_i \bar{p}+\mu_k \bar{p}\right)}{|Q|\left(\sigma_k^2+\varepsilon\right)} $

      (3)

      $ b_k=\bar{p}-a_k \mu_k $

      (4)

      式中: |Q|为窗口像素数;σkμk为函数标准差和均值;p为第k个局部窗口内像素灰度均值。

      为了提高图像中背景和砂轮投影的对比度,突出砂轮轮廓特征,方便后续处理,采用直方图均衡改善图像中灰度的分布[20-21],增强局部的对比度,而不影响整体的对比度。砂轮轮廓图像的直方图如图 5所示。根据砂轮轮廓图像特征和此处对比度增强的需求,本文中图像增强使用分段线性变换,其变换曲线如图 6所示,图中两个锚点的位置可根据实际图像进行微调。

      图  5  砂轮图像直方图

      Figure 5.  Histogram of grinding wheel images

      图  6  分段线性变换曲线

      Figure 6.  Segmented linear transformation curve

      轮廓图像经过上述图像滤波和图像增强的预处理后,其效果如图 7所示。

      图  7  图像预处理效果

      Figure 7.  Images pre-processing effect

    • 针对传统视觉检测中单一截面不能反映砂轮全域轮廓信息、检测精度差的问题,本文作者提出了一种基于圆周方向多截面轮廓图像合成的检测方法,算法的典型流程如图 8所示。

      图  8  轮廓合成算法流程图

      Figure 8.  Profile synthesis algorithm flow chart

    • 预处理后的轮廓图像背景和砂轮轮廓在图像中的灰度差异较大,因此可以通过阈值分割判定轮廓边缘[22]。本文中使用一个全局阈值将砂轮的轮廓图像划分为背景和砂轮轮廓两个特征,且预处理后的图像亮度分布均匀。在全局阈值分割过程中,将砂轮本体部分的像素灰度值置为1,将背景部分的像素灰度值置为0,即反二值化操作,其数学表达式如下:

      $ G_0(x, y)=\left\{\begin{array}{l} 1, \left(0 \leqslant f_0(x, y)<T_0\right) \\ 0, \left(f_0(x, y) \geqslant T_0\right) \end{array}\right. $

      (5)

      式中: f0(x, y)为输入图像(x, y)点的灰度值; G0(x, y)为输出图像(x, y)点的像素值; T0为全局分割阈值。

      按照上述方法进行阈值分割和反二值化的效果如图 9所示。

      图  9  阈值分割和反二值化效果

      Figure 9.  Effects of thresholding segmentation and inverse binarization

    • 本文中通过开运算形态学操作,减少轮廓图像中的干扰杂点[23],如吸附在砂轮表面和镜头表面的灰尘和毛絮,如图 10所示。

      图  10  轮廓图像中杂质干扰(放大)

      Figure 10.  Impurity interference in profile images (zoom)

      开运算由是图像膨胀和图像腐蚀两个基本操作的组成, 其数学表达式为:

      $ G_1(x, y)=\max \limits_{\left(x^{\prime}, y^{\prime}\right) \in e\left(x^{\prime}, y^{\prime}\right) \neq 0}\left[f_1\left(x+x^{\prime}, y+y^{\prime}\right)\right] $

      (6)

      $ G_2(x, y)=\min \limits_{\left(x^{\prime}, y^{\prime}\right) \in e\left(x^{\prime}, y^{\prime}\right) \neq 0}\left[f_1\left(x+x^{\prime}, y+y^{\prime}\right)\right] $

      (7)

      式中: e代表结构元; (x, y)为结构元锚点的位置; (x′, y′)为结构元内值为1的元素相对其锚点的位置; G1(x, y)、G2(x, y)为输出图像; f1(x, y)为输入二值图像。

    • 经过二值化和形态学处理后的图像中,像素值为1的位置对应砂轮实体,像素值为0的位置为背景。多轮廓图像合成算法即计算图像序列中图像像素位置砂轮实体出现的频次,得出轮廓频率图像。本文中检测系统设计的多轮廓图像合成算法有两种场景,如图 11所示。一种是砂轮旋转一圈采集一组图像序列然后进行静态合成,另一种是砂轮持续旋转实时对砂轮轮廓进行旋转采样并进行动态合成。

      图  11  两种轮廓图像合成场景

      Figure 11.  Two kinds of profile images synthesis scenes

      静态合成算法只需要对一组圆周方向采样的砂轮轮廓图像转为二值图像序列,然后对序列中的图像相同像素位置求和即可。静态合成算法的数学表达式为:

      $ F_1(x, y)=\sum\limits_{j \in N} h_j(x, y) $

      (8)

      式中: N为输入二值化和形态学处理后的图像序列; hj(x, y)为输入序列中第j张图像二值反转和形态学处理后(x, y)点的像素值; F1(x, y)为合成后的频率图像(x, y)点的值。

      动态合成算法是在实时图像序列中动态合成最新的n张图像。算法处理流程为:每采集一张图像先进行二值化和形态学处理,实时图像序列中第1张~第n张图像的合成与静态合成算法一致,实时序列中第n+1张图像的合成需要先将动态序列第n+1张图像累加到频率图像中,然后频率图像减去第1张图像,则得到第2张~第n+1张的合成频率图。以此类推,在动态累加最新图像第j张到频率图像时,要减去第j-n张图像。动态合成算法的数学表达式为:

      $ F_2(x, y)=\left\{\begin{array}{l} \sum\limits_{j=1}^n h_j(x, y), (1 \leqslant j \leqslant n) \\ F_2(x, y)+h_j(x, y)-h_{j-n}(x, y), \\ (j>n) \end{array}\right. $

      (9)

      式中: n为合成窗口大小; hj(x, y)为输入序列中第k张图像二值反转后(x, y)点的像素值; hj-n(x, y)为输入序列中第j-n张图像二值反转和形态学处理后(x, y)点的像素值; F2(x, y)为输出频率图像(x, y) 点的值。

      砂轮圆周方向采样的轮廓图像序列经过多轮廓图像合成算法处理后的效果如图 12所示。图 12a图是轮廓图像序列中的某一帧; 图 12b图是256张轮廓图像合成后的频率图像。

      图  12  轮廓图像合成算法效果

      Figure 12.  Effect of profile images synthesis algorithm

    • 多轮廓图像合成后生成的频率图像本质上是一张灰度图,将频率图像进行局部放大,如图 13所示。频率图像中像素的值表示轮廓图像序列中该像素为砂轮实体的次数。本文作者将像素值为0的像素定义为绝对背景,即序列中所有的轮廓图像上该像素均对应背景;将像素值为图像序列长度的像素定义为绝对砂轮实体,即序列中所有的轮廓图像上该像素均对应砂轮实体;将频率图像中背景和砂轮轮廓之间的灰度过渡像素定义为轮廓波动带,轮廓波动带反映了圆周方向砂轮轮廓的综合信息。

      图  13  频率图像局部放大

      Figure 13.  Frequency image partial zoom

      为了简化后续的检测算法流程,先将频率波动带从频率图像中提取出来。由轮廓波动带的定义和上述分析可知,频率图像中灰度在一定范围内的像素即为轮廓波动带,因此使用阈值分割即可实现轮廓波动带的提取。此处的阈值分割数学表达式为:

      $ v(x, y)=\left\{\begin{array}{l} F(x, y), \left(T_1<F(x, y)<T_2\right) \\ 0, \left(F(x, y) \leqslant T_1 \text { or } F(x, y) \geqslant T_2\right) \end{array}\right. $

      (10)

      式中: v(x, y)为轮廓波动带图像中(x, y)的值;F(x, y)为频率图像(x, y)点的值,静态合成时等于F1(x, y),动态合成时等于F2(x, y);T1为阈值区间下边界;T2为阈值区间上边界。

      经实验,若m为图像序列的长度大小,在区间上边界阈值为0.005 m~0.02 m且区间下边界阈值为0.99 m~0.999 m时,轮廓波动带的精度较高。轮廓波动带提取的效果如图 14所示。

      图  14  轮廓波动带提取效果

      Figure 14.  Profile fluctuation with extraction effect

    • 成形砂轮轮廓通常由直线段和圆弧段组成。通过对从频率图像中提取得到的轮廓波动带最外侧边缘拟合,可以测量砂轮的最外侧轮廓。相较与传统的单一图像检测方法,它能更真实地反映砂轮磨削工件复映的轮廓。但是,轮廓波动带边缘可能存在一些异常凸起和凹坑,轮廓边缘可能会出现部分随机离群噪点,这对直线和圆弧的拟合算法的鲁棒性提出一定的要求。因此,本文中使用随机抽样一致算法(random sample consensus,RANSAC)结合最小二乘法进行特征拟合,计算流程如下[24]:(a)随机从原始数据集S中取出s个元素,s个元素刚好能确定待求数学模型,求解待求模型的假设模型;(b)根据步骤(a)中得出的假设模型,计算数据集中所有点的残差E,将残差小于最低残差阈值t的点加入内点集Si;(c)如果内点集Si的大小大于设定阈值T,使用内点集Si进行最小二乘法拟合,并结束;(d)如果内点集Si的大小小于设定阈值T,则重复上述3个步骤;(e)达到最大循环次数N后,取N次循环中最大的内点集Si进行最小二乘法拟合。

      上述过程中sE参数由被拟合的特征类型决定,对于直线拟合,两点即可确定一条直线,故在RANSAC中s=2,残差E为数据集中点到假设直线的距离,数学表达式为:

      $ E=\frac{a x_i+b y_i+c}{\sqrt{a^2+b^2}} $

      (11)

      式中: a, bc为假设直线模型ax+by+c=0的参数; (xi, yi)为点集中第i个点的坐标。

      对于圆和圆弧拟合,3点即可确定一个圆,故在RANSAC中s=3,残差E为数据集中点到假设圆的距离,数学表达为:

      $ E=\sqrt{\left(x_i-x_0\right)^2+\left(y_i-y_0\right)^2}-R $

      (12)

      式中:x0, y0和R为假设圆模型(x-x0)2+(y-y0)2=R2的参数。

      最大循环次数N可由理论推导,假设原始点集中内点的占比为p,则ps表示取s个点都是内点的概率,1-ps表示取s个点至少有一个为外点的概率,即单次采样失败概率,(1-ps)M表示M次循环采样均失败的概率,1-(1-ps)M则表示M次循环采样至少一次成功的概率。因此,假设采样成功率为w,则最少循环次数为ln(1-w)/ln(1-ps)。若假设点集中内点占比为0.9,采样成功率要求为0.99,则拟合直线时最大循环次数N=3,拟合圆弧时最大循环次数N=4即可。

      内点集样本数达标阈值T和最低残差阈值t对算法的鲁棒性有很大影响,需要根据使用场景确定。经实验,在本文的场景中,内点集样本数达标阈值T为原始数据集S大小的0.9倍,最低残差阈值t=30时,有较好的拟合效果。如图 15所示,横纵坐标XY的数值表示砂轮V槽轮廓中的直线部分的一组点在2维平面的位置。对比直接最小二乘法,结合RANSAC的最小二乘法对离群点有更好的抗干扰性,符合本文作者预期设计。

      图  15  拟合效果对比

      Figure 15.  Comparison of fitting effects

    • 本实验中选用凹V型成形砂轮作为被测对象,如图 16a所示,其轮廓形状尺寸如图 16b所示。轮廓具有直线和圆弧特征,是最基本的成形砂轮轮廓,适合作为本文中的验证对象。

      图  16  凹V型轮廓成形砂轮实物及尺寸情况

      Figure 16.  Concave V-slot forming grinding wheel and dimensional situation

      离线检测时,需要先使用砂轮激光修整机床上的磨削复映装置夹持复映板材进行磨削,磨削复映机构如图 17所示。通过向负方向进给磨削复映装置的进给机构P轴,即可在复映板材上复映出砂轮的轮廓。

      图  17  磨削复映装置

      Figure 17.  Grinding remapping device

      为了减少磨削复映板材时砂轮自身的磨损,以及提高磨削效率,本实验中所用复映板材为细石墨板,实物图如图 18所示。同时石墨具有良好的耐热、耐压和润滑性,较适合干式磨削。本文中所选用的石墨板尺寸为30 mm×30 mm×3 mm,纯度大于99.5%,密度为1.82 g/cm3,肖氏硬度为35HS,平均颗粒直径10 μm。磨削石墨试片时,砂轮线速率为20 m/s,进给速率为0.1 mm/min,冷却方式为吹起冷却。

      图  18  石墨板

      Figure 18.  Graphite sheet

      对复映有砂轮轮廓的石墨板,采用超景深显微镜的平面测量功能进行尺寸检测,具有最大20倍的光学放大倍率,搭配正面照射和背面照射双光源,支持全自动对焦、平面测量、3维测量等功能,满足高精度离线检测石墨片需求。

      视觉单一截面检测法与本文作者提出的多轮廓图像合成检测法都是通过在线检测平台获取图像,如图 19所示。视觉单一截面检测对照实验时,使用已有的影像测量软件进行静态检测,其主界面如图 20所示。该软件已广泛应用于多款精密影像测量仪设备,支持图像实时测量、自动边缘拟合、几何标注等功能,具有良好的测量精度。

      图  19  在线检测平台

      Figure 19.  Online detection platform

      图  20  影像测量软件主界面

      Figure 20.  Image measurement software main interface

    • 实验中以凹V型槽成形砂轮作为实验对象,使用砂轮激光修整机床的纳秒脉冲激光修整其轮廓形状,修整完成后分别使用磨削复映检测法、视觉单一截面检测法和本文作者提出的多轮廓图像合成检测法对砂轮上多个V型槽的轮廓尺寸进行检测,对比几种方法的检测结果,验证本文中检测系统的精度。

    • 经过标定保证检测精度之后,按照上述的试验方案,使用3种方法对成形砂轮轮廓8个V型槽进行检测,检测过程如图 21所示。

      图  21  3种测量方式检测过程

      Figure 21.  Three measurement methods detection process

      分别测量了V型槽角度D、深度H和槽底圆弧半径R,测量结果分别如表 1表 2表 3所示。表中,单一截面检测法的测量结果为4个截面测量结果的平均值。检测效率上,计算完整检测单个V型槽的平均检测时间结果如表 4所示。

      表 1  V型槽角度D检测结果

      Table 1.  Detection results of V-slot angle D

      the number of V-slot the multi-profile images synthesis detection method/mm the single-section visual detection method/mm the grinding remapping detection method/mm
      1 44.213 44.156 44.281
      2 44.087 43.915 44.042
      3 44.245 44.323 44.203
      4 44.076 44.031 44.109
      5 44.070 44.264 44.115
      6 44.208 44.344 44.293
      7 44.168 44.199 44.156
      8 44.092 43.984 44.067

      表 2  V型槽深度H检测结果

      Table 2.  Detection results of V-slot depth H

      the number of V-slot the multi-profile images synthesis detection method/mm the single-section visual detection method/mm the grinding remapping detection method/mm
      1 0.5251 0.5274 0.5183
      2 0.5294 0.5396 0.5206
      3 0.5342 0.5362 0.5299
      4 0.5297 0.5373 0.5222
      5 0.5300 0.5450 0.5255
      6 0.525 0.5338 0.5177
      7 0.5229 0.5287 0.5197
      8 0.5281 0.5471 0.5207

      表 3  槽底圆弧半径R检测结果

      Table 3.  Detection results of V-slot bottom arc radius R

      the number of V-slot the multi-profile images synthesis detection method/mm the single-section visual detection method/mm the grinding remapping detection method/mm
      1 0.1664 0.1611 0.1749
      2 0.1624 0.1492 0.1686
      3 0.1586 0.161 0.1672
      4 0.1621 0.1566 0.1729
      5 0.1597 0.1584 0.1671
      6 0.1638 0.1571 0.1717
      7 0.1662 0.1623 0.1706
      8 0.1642 0.1528 0.1712

      表 4  检测效率

      Table 4.  Detection efficiency

      the multi-profile images synthesis detection method the single-section visual detection method the grinding remapping detection method
      time spent/s 30 20 360
    • 3种特征检测结果中,多轮廓合成检测法的检测结果与磨削复映检测法更为贴近,而单一截面检测法的结果存在较大的随机波动。在V型槽深度的检测中,多轮廓合成检测和单一截面检测的结果普遍偏小;在底部圆弧半径的检测中,多轮廓合成检测和单一截面检测的结果普遍偏大,这是直接测量和间接测量之间的误差。一种原因可能是对石墨板进行磨削时,受石墨板的弹性变形和磨削过程热膨胀等因素影响,石墨板上复映的轮廓与砂轮轮廓不一致,产生固定误差; 另一种推测是砂轮表面未清理干净,仍有细小灰尘毛絮等附着物,同时砂轮轮廓存在不可避免的散射和反射,导致了最终采集的轮廓图像与实际轮廓之间存在误差。本文作者以工业界普遍认可的磨削复映法检测结果为基准,分别计算了单一轮廓检测和多轮廓检测两种方法的检测误差值,结果如表 5所示。

      表 5  检测误差

      Table 5.  Detection errors

      detection features detection methods the number of V-slot
      1 2 3 4 5 6 7 8
      angel/(°) the single-section visual detection method -0.265 -0.157 0.16 0.122 0.169 0.071 0.043 -0.203
      the multi-profile images synthesis detection method -0.082 0.065 0.091 -0.033 -0.105 -0.085 0.012 0.025
      depth/mm the single-section visual detection method 0.0091 0.019 0.0063 0.0151 0.0195 0.0161 0.009 0.0264
      the multi-profile images synthesis detection method 0.0068 0.0088 0.0043 0.0075 0.0045 0.0073 0.0032 0.0074
      radius/mm the single-section visual detection method -0.0138 -0.0194 -0.0062 -0.0163 -0.0087 -0.0146 -0.0083 -0.0184
      the multi-profile images synthesis detection method -0.0085 -0.0062 -0.0086 -0.0108 -0.0074 -0.0079 -0.0044 -0.007

      观察V型槽角度的检测误差,几组检测样本中,两种检测方法都有较大的偏差,其中多轮廓合成检测法与磨削复映检测法的测量值偏差相对较小,最大正负误差分别为+0.045°和-0.085°,而单一截面检测法的测量值偏差较大,最大正负误差分别为+0.149°和-0.127°。V型槽角度的测量通常是先对V型槽两侧边进行直线拟合,再计算两直线夹角,当直线长度较短时,直线的拟合精度较差,会导致角度测量结果偏差较大。

      槽底深度的检测误差和底部圆弧半径检测误差都分别偏向一个方向,但多轮廓合成检测误差更小,且误差的变化范围小。而单一截面检测,由于是仅抽取砂轮某一个截面的轮廓进行检测,随机性很大,导致其误差更大,且误差变化范围也更大。两种方法检测误差的统计指标如表 6所示。

      表 6  检测误差的统计指标

      Table 6.  Statistical indicators of the detection errors

      detection features statistical indicators the multi-profile images synthesis detection method the single- section visual detection method
      mean value/mm -0.0076 -0.0132
      radius range/mm 0.0064 0.0132
      variance 3.0925×10-6 2.13×10-5
      mean value/mm 0.0062 0.01506
      depth range/mm 0.0056 0.0201
      variance 3.36937×10-6 3.92×10-5

      由于直接检测和间接检测有方法上的差异,检测结果存在固定偏差值,导致存在7 μm左右的平均误差,但测量误差的极差可以在一定程度上反映检测精度。多轮廓合成检测法的误差极差最大值为6.4 μm,单一截面检测法的误差极差最大值达到20.1 μm,多轮廓合成法误差范围约为单一截面检测的1/3,同时误差方差要小一个数量级。综上分析,本文作者提出的多轮廓合成检测法效率远超磨削复映法,与单一截面检测法效率接近;检测结果上,与磨削复映法虽然有所差异,但误差值相对稳定,可通过补偿进行改善,相较于单一截面检测法,本文检测方法的误差更小,且误差变化范围约为±3 μm,检测精度更高。

    • 介绍了基于砂轮激光修整机床搭建的图像采集系统的组成和工作原理,通过分析砂轮激光修整的特点,提出了一种多轮廓图像合成检测方法,定义了轮廓波动带的概念。使用合适的图像预处理算法,分析多轮廓图像合成算法的实现,并从合成结果中提取轮廓波动带,进一步研究基于轮廓波动带的轮廓测量。通过对凹V型轮廓成形砂轮进行检测试验,与磨削复映法和单一截面检测法进行对比,结果表明:多轮廓图像合成检测方法的检测结果与磨削复映法的检测结果接近,检测效率大幅提升,由于直接检测和间接检测方法上的差异,存在固定误差值,但误差的变化范围较小,约为±3 μm,相较于单一截面检测法,多轮廓图像合成检测方法的检测误差显著下降,检测效率接近,同时稳定性更佳,符合预期设计。

      随着未来硬件设备的更新和算法效率的提升,本方法可与生成式人工智能(artificial intelligence generated content,AIGC)技术融合。借助辅助生成方法,通过输入关键词或需求,将待检测轮廓图像的尺寸精度和轮廓信息全部识别和提取出来,甚至形成专门的凹V型轮廓砂轮的数据集进行增强学习,进一步提高检测效率和精度,降低检测成本,但是要注意评判最终结果的真实性。

参考文献 (24)

目录

    /

    返回文章
    返回