[发明专利]一种不确定环境下表面粗糙度在线测量方法

摘要

本发明提供一种不确定环境下表面粗糙度在线测量方法，针对实际环境中利用机器视觉进行表面粗糙度测量时环境光的影响，提出了通过建立环境光、表面特征量与表面粗糙度关系模型的方法来计算待测表面粗糙度值。利用本方法，可以显著提高在不同外界环境光条件下对表面粗糙度的测量精度。本发明通过下述技术方案予以实现：首先根据视觉图像利用算术平均偏差法和灰度共生矩阵法得到表面特征量：表面的统计特征和纹理特征，并用光学照度计实时监测环境光的变化情况，在此基础上建立环境光、表面特征量与表面粗糙度关系模型，最后利用模型确定的函数计算表面粗糙度，利用本方法显著提高了在线测量精度。

主权项

1、一种不确定环境下表面粗糙度在线测量方法，其特征在于，利用机器视觉技术，在线测量不确定环境下工件表面的粗糙度，具体步骤如下：(1)确定光下建模：A、取表面粗糙度值已知的第一工件置于生产现场，第一工件的表面粗糙度真实值为R_a；调节现场LED光源，使此时环境光亮度值为T_w，利用CCD摄像机获取工件表面的图像，并将获取的图像转换为灰度图；B、根据步骤A获取的灰度图，提取采样区域的灰度均值μ_G作为第一工件表面图像第一特征量T_G1，计算公式如下：TG1=μG=(Σi=1N1Σj=1N2gij)/N1N2]]>①其中g_ij为像素灰度值，i，j为自然数；N₁，N₂分别为采样区域行、列值；C、根据步骤A获取的灰度图，提取采样区域灰度共生矩阵的能量值E作为第一工件表面图像第二特征量T_G2，表征采样区域灰度值的空间相关性信息，计算公式如下：TG2=E=Σi=0L-1Σj=0L-1P‾2(i,j,d,θ)]]>②其中，P‾(i,j,d,θ)=P(i,j,d,θ)/Σi=0L-1Σj=0L-1P(i,j,d,θ)]]>③P(i，j，d，θ)是灰度图中具有同样亮度或接近亮度的象素间的位置分布特性的联合概率密度，计算公式如下：P(i，j，d，θ)＝ψ{[(x，y)，(x+D_x，y+D_y)]|f(x，y)＝i，f(x+D_x，y+D_y)＝j}  ④式中：ψ表示该集合中元素的数目；x，y为自然数，是N×N图像的像元坐标；i，j＝0，1，…，L-1，L为自然数，是灰度等级为L的灰度值；D_x，D_y为根据灰度矩阵生成步长d和方向θ得到的位置偏移量，其中d为自然数，θ可取0°，45°，90°，135°四个方向；D、更换第一工件为另一个与其结构及加工方式均相同，但表面粗糙度值不同的工件，重复步骤A～步骤C，并记录各个步骤得到的数据；E、重复步骤D至少N₁次，N₁为自然数，且N₁≥5；并根据每次得到的数据，利用数值分析中的拟合方法，建立确定光下的数学模型：R_a＝f(T_G1，T_G2)                            ⑤F、模型检验：取表面粗糙度值已知的第二工件置于生产现场，第二工件表面粗糙度真实值为R；利用CCD摄像机获取工件表面的图像，将获取的图像转换为灰度图，并根据步骤B和C计算第二工件表面图像第一特征量T_G1′和第二工件表面图像第二特征量T_G2′；G、将步骤F得到的第二工件表面图像第一特征量T_G1′和第二工件表面图像第二特征量T_G2′，代入步骤E确定的函数模型⑤，得到第二工件表面粗糙度计算值R′_a；H、对比第二工件表面粗糙度计算值R′_a和第二工件表面粗糙度真实值R，判断结果是否符合精度要求；当结果不符合精度要求时，增大步骤E中N₁的值，重复步骤E～步骤H，直到结果符合精度要求为止；当结果符合精度要求时，则可判断步骤E得到的数学模型⑤为确定光下的正确函数模型，进入下一步骤；(2)不确定光下建模：I、调节现场LED光源，改变步骤A中的环境光亮度值T_w，重复步骤A～步骤H至少N₂次，N₂为自然数，且N₂≥8；J、根据每次的环境光亮度值T_w以及由步骤H确定得到的确定光下的数学模型，再次利用数值分析中的拟合方法，建立不确定光下的函数模型：R_a＝f(T_W，T_G1，T_G2)                              ⑥K、模型检验：取表面粗糙度值已知的第三工件置于生产现场，第三工件表面粗糙度真实值为R₃；利用CCD摄像机获取工件表面的图像，将获取的图像转换为灰度图，并根据步骤B和C计算第三工件表面图像第一特征量T_G3和第三工件表面图像第二特征量T_G4；L、用照度计获取此时生产现场环境光亮度值，记为T_w3；M、将步骤F得到的第三工件表面图像第一特征量T_G3和第三工件表面图像第二特征量T_G4以及步骤G得到的环境光亮度值T_w3，代入步骤J确定的函数模型⑥，得到第三工件表面粗糙度计算值R′₃；N、对比第三工件表面粗糙度计算值R′₃和第三工件表面粗糙度真实值R₃，判断结果是否符合精度要求；当结果不符合精度要求时，增大步骤I中N₂的值，重复步骤I～步骤N，直到结果符合精度要求为止；当结果符合精度要求时，则可判断步骤J得到的数学模型⑥为不确定光下的正确函数模型，进入下一步骤；(3)在线应用：O、将待测工件置于生产现场，用CCD摄像机获取待测工件表面图像，并将获取的图像转换为灰度图；P、根据步骤B和C计算待测工件表面图像第一特征量T″_G1和待测工件表面图像第二特征量T″_G2；Q、用照度计获取生产现场环境光亮度值，记为T″_w；R、根据步骤J确定的函数模型⑥，代入步骤P得到的待测工件表面图像第一特征量T″_G1和待测工件表面图像第二特征量T″_G2、以及步骤Q得到的生产现场环境光亮度值T″_w，计算出待测工件表面粗糙度值R″_a。