基于激光散斑图像多特征参数的表面粗糙度建模研究
0 引 言
表面粗糙度能对工件的配合性质、耐磨性、接触性能、疲劳强度等产生较大影响,最终影响到机器和零件的工作性能和寿命,对准确评估工件的表面粗糙度具有重要意义[1]。表面粗糙度测量的经典工具是触针式轮廓仪,它不仅可以测量表面粗糙度值,还可以记录表面轮廓。但是,触针式轮廓仪存在可能会划伤被测表面,测量过程耗时较长等缺点[2]。与接触式测量方法相比,非接触式光学测量法具有高效率、低损耗、耐高温等特点。非接触式光学测量法主要包括电镜法[3-4]、聚焦法[5-6]、干涉法[7-8]、散射法[9-10]、散斑法[11-12]等。Persson通过建立散斑对比度与均方根粗糙度(Rq)之间的关系,实现了对Rq≤0.2 μm的研磨工件表面粗糙度的测量,但是该方法的测量范围较小[13]。Patel等利用灰度共生矩阵提取了工件表面图像的12个特征参数,用机器学习算法实现了工件加工类型的识别,但是建立的模型存在特征冗余问题[14],Patel等还发现散斑二值图像白黑像素比与平磨工件的轮廓算术平均偏差(Ra)之间存在线性相关性,但适用性的广度仍需增强[15]。Chen等通过拟合散斑图像Tamura纹理粗糙度特征和Ra的函数关系,建立了基于最小二乘回归的表面粗糙度测量模型[16]。Goh等提取了铣削表面图像的多个特征,发现平均像素强度与表面算术平均高度(Sa)的相关性最高[17]。Haridas等通过建立散斑图像分形参数与表面粗糙度的关系,实现了对抛光表面的粗糙度测量[18]。Shao等提取了散斑图像的多个特征参数,建立了基于散斑图像特征参数的多元回归模型,实现了对1 μm≤Sa≤2 μm的冷轧带钢表面粗糙度的测量,与单一特征参数的模型相比,建立的多参数模型性能更优[19]。这些研究中,不同的方法具有不同的测量范围、灵敏度和稳定性。基于单一特征参数的表面粗糙度测量方法提取的图像信息较少,单一特征受环境影响时稳定性变差,随着表面粗糙度测量场景的复杂多变和精度要求的不断提高,该方法适应性在逐渐下降[19]。基于多特征参数的表面粗糙度测量方法提取的原始特征参数不一定都与表面粗糙度有良好的相关性,特征参数之间还可能存在特征冗余问题,不相关特征和冗余特征会耗费特征提取时间,增加模型的计算成本和复杂度,部分不相关特征和冗余特征还会降低模型的准确性和稳定性。因此,研究基于多特征参数的表面粗糙度测量方法在建模过程中如何剔除不相关特征和冗余特征,并筛选出一组特征,建立可实现不同加工类型表面粗糙度测量的模型具有重要的意义。
文中搭建了激光散斑图像采集系统,采集了平磨、卧铣、立铣、研磨标准试件的激光散斑图像,从试件的散斑图像中提取了多个特征参数。通过引入斯皮尔曼相关系数并制定简约规则,以及改进序列后向选择算法,解决了特征相关和特征冗余问题,并利用筛选出的这组特征,建立了可识别加工类型和测量不同加工工艺表面粗糙度的模型。
1 建模原理和过程
粗糙表面被激光照射时,反射波会在观察平面上形成亮斑和暗斑,这些随机分布的亮斑和暗斑被称为激光散斑。在统计学方法中,用光强的概率密度函数描述激光散斑[20]。描述如下:
式中:I(r)表示点r处散斑场的光强。
公式(1)中含有待定参数。因此,通过光强的概率密度函数直接建立表面粗糙度测量模型具有较大难度。此外,激光散斑携带大量物体表面的信息,表面粗糙度是重要的表面信息之一,故文中从采集的激光散斑图像中提取特征参数,建立了基于支持向量机的表面粗糙度测量模型。从激光散斑图像中提取的原始特征参数表示为Tn ={t1, t2, …, tn},其中n表示提取的原始特征参数个数。整个建模过程的示意图如图1所示。
1.1 散斑特征提取
激光散斑图像的一阶统计量是对整个图像中像素值的统计描述,通常用来表征图像的全局特征,不涉及图像的局部结构或纹理等信息。激光散斑图像的二阶统计量考虑了像素值之间的空间关系,可以用来表征图像的局部特征和纹理特征。文中用灰度共生矩阵法和灰度差分统计法提取了散斑图像的8个纹理特征;提取二值图像前景像素占比和灰度图像的灰度均值、灰度标准差、灰度均方根作为一阶统计特征。提取的特征参数如下:
1)用灰度共生矩阵法提取了激光散斑图像的5个纹理特征。它们的计算公式为:
式中:E为能量;S为熵;I为惯性矩;L为相关性;H为逆差矩;Gg为灰度级数;
2)用灰度差分统计法提取了激光散斑图像的3个纹理特征。它们的计算公式为:
式中:Mmean为平均值;Acon为对比度;Bent为熵;
3)二值图像的前景像素占比定义为二值图像中像素值为1的像素数与总像素数的比值用
式中:h(i, j)表示二值图像的像素值;Nx和Ny表示图像水平方向和竖直方向的像素数。
4)散斑灰度图像的灰度均值κ、灰度标准差σ、灰度均方根υ。它们的计算公式为:
式中:Ig(i, j)表示灰度图像各像素点的灰度值。
1.2 特征参数预筛选
提取的每个原始特征参数都提供了散斑图像的信息,但是,并非每一个原始特征参数都与表面粗糙度参数有良好的相关性。特征参数越多,提取特征参数的过程越复杂,模型数据维数也越高。斯皮尔曼相关系数有不依赖样本分布、鲁棒性较强等优点,可以衡量两个变量之间的单调相关性。因此,文中引入斯皮尔曼相关系数,制定简约规则对提取的散斑特征参数进行预筛选。斯皮尔曼相关系数计算公式如下:
式中:ρxy表示变量x和y的相关性;R(xi)和R(xj)表示xi和yi的位次;
1.3 基于支持向量机的模型建立
文中建立了基于支持向量机的表面粗糙度测量模型。支持向量机(SVM)是一种适合小样本问题的分类和回归算法。其原理如下。
给定N个样本数据集合U={(
式中:
式中:K(
将SVC推广到回归问题可以得到支持向量回归(SVR)。设在高维空间的线性回归函数为:
式中:f(
式中:ξi和ξi为松弛变量;ε为不敏感函数参数。引入拉格朗日乘子αi和αi,将公式(18)转化为对偶问题:
选用高斯径向基核函数,最终的SVR函数公式:
式中:g为核函数参数。模型有两个非常重要的超参数:正则化系数C与核函数参数g。
1.4 改进的序列后向选择算法优化模型
序列后向选择算法采用启发式搜索策略,能够有效地剔除冗余特征。其思想是:以特征全集为起点,每次剔除一个特征,使得剔除特征后的性能评价函数值达到最优。该算法对多个模型进行优化时,从各模型剔除的冗余特征不一定相同,这可能会导致各模型最终保留的的特征子集不同。文中针对该问题提出了改进的序列后向选择算法,步骤如下:
1)确定各模型的性能评价函数,把它们记为:{a1, a2, ··· , ak},其中k为模型个数;
2)确定各模型的性能评价函数之和,把它记为:γ = a1+a2+ , ··· , +ak;
3)设置各模型的性能评价函数的阈值,把它们记为:{η1, η2, ···, ηk};
4)获取各模型的特征全集,每个模型的输入由其特征全集组成,用特征全集训练各模型,计算γ的值;
5)减去一个特征,用剩余的特征训练各模型,计算γ的值。若γ更优,并且各模型的性能评价函数值在设置的阈值范围内,则保留特征,反之则剔除该特征;
6)重复步骤5),直到某个模型的性能评价函数值超出阈值或γ达到最优。
改进的序列后向选择算法能同时对多个模型进行优化,各模型每次剔除相同的特征。在解决特征冗余问题的同时,实现了各模型保留相同的特征子集。
1.5 模型评价标准
采用平均绝对百分比误差(MAPE)和均方根误差(RMSE)作为模型的性能评价指标,公式如下:
式中:yi为标准值;
2 实验设计
2.1 表面粗糙度试件选择
该实验中标准试件的材料和表面粗糙度参数如表1所示。
表 1. 材料和表面粗糙度参数
Table 1. Material and surface roughness parameters
|
2.2 实验装置设计
激光散斑图像采集装置示意图如图2所示,该装置主要包括半导体激光器(635 nm、1 mW)、CCD相机(190万 pixel)、远心镜头、激光器专用恒流电源、平移台。
采集的激光散斑图像与激光的入射角度、散射表面相对于激光的位置、激光功率及激光光斑的大小都有关系。该实验半导体激光器功率固定为1 mW,入射角度固定为与粗糙度试件表面法向呈10°,这确保了散射表面相对于激光始终处于相同的角度和位置,激光光斑大小始终不变。CCD相机和远心镜头用于记录散斑图案,带有镜头的CCD相机位于粗糙度试件表面的法线方向,与试件表面的距离固定为110 mm,这确保了观测平面在整个实验过程中都是一致的,并且相机接收到的镜面反射光较小,有效避免了相机信号饱和。实验中通过控制二维平移台移动试件,采集试件表面的散斑图像,每次曝光时都确保二维平移台和试件是静止的,整个过程中无需移动激光器,也不会影响其相对于散射表面的方向和位置。
2.3 散斑图像采集和数据集划分
实验选用平磨、卧铣、立铣、研磨四种不同加工工艺的粗糙度标准试件,共计15个标准试件。对每一个标准试件采集46张散斑图像,共计采集690张散斑图像。将采集的散斑图像划分为训练集和测试集,具体的样本划分如表2所示。
表 2. 数据集样本划分
Table 2. Dataset sample division
|
3 结果与分析
3.1 激光散斑图像预处理
激光照射粗糙度标准试件,相机采集到的图像中除了激光散斑图案外还有背景区域。为消除边缘效应,保留像素大小为500 pixel×500 pixel的区域作为激光散斑图案的有效部分。然后用高斯滤波过滤图像中的噪声,预处理后的部分激光散斑图像见图3。
图 3. 预处理后的激光散斑图像(平磨试件)。(a) R a =0.1 μm; (b) R a =0.2 μm; (c) R a =0.4 μm; (d) R a =0.8 μm
Fig. 3. Preprocessed laser speckle images (Plane grinding specimens). (a) Ra =0.1 μm; (b) Ra =0.2 μm; (c) Ra =0.4 μm; (d) Ra =0.8 μm
3.2 特征参数预筛选
从预处理后的激光散斑图像中提取了27个原始特征参数。灰度共生矩阵法提取了0°、45°、90°、135°方向的能量、熵、惯性矩、相关性和逆差矩特征,但是这样得到的特征参数过于繁多,取4个方向的平均值作为灰度共生矩阵的特征参数值。然后,引入斯皮尔曼相关系数,设置门限值Q为0.8,如果特征参数与各加工类型的表面粗糙度参数Ra之间的斯皮尔曼相关系数绝对值均高于门限值,则保留该特征参数,否则剔除该特征参数。各特征参数与表面粗糙度参数Ra之间的斯皮尔曼相关系数绝对值如表3所示。
表 3. 特征参数与表面粗糙度参数之间的斯皮尔曼相关系数绝对值
Table 3. The absolute value of Spearman's correlation coefficient between feature parameters and surface roughness parameter
|
从表3可以看出,L与平磨试件表面粗糙度参数的相关性低于门限值,Mmean、Acon、D与卧铣和立铣试件表面粗糙度参数的相关性低于门限值,因此剔除L、Mmean、Acon、D,4个特征参数。经过预筛选后,筛选出的8个特征参数均与各加工类型表面粗糙度参数Ra有强相关性,记为:TV ={E, S, I, H, Bent, κ, σ, υ}。
3.3 表面粗糙度测量模型及优化
经过特征参数的预筛选后,分别对平磨、卧铣、立铣、研磨试件表面建立基于多特征参数TV ={E, S, I, H, Bent, κ, σ, υ}的SVR模型,以实现四种加工类型表面粗糙度的测量,建模过程中不敏感函数参数取值为ε = 0.01。SVR模型的超参数和输入特征均会影响模型的性能。
优先考虑超参数对模型性能的影响,SVR模型的正则化系数C和核函数参数g会对模型的性能产生重要影响。文中利用5折交叉验证和网格搜索法对模型的超参数C和g进行寻优,各SVR模型最优的C和g如表4所示。将各加工类型的测试集数据分别输入到相应的SVR模型中,模型对测试集样本表面粗糙度预测的RMSE和MAPE如表5所示。
表 4. 8输入特征SVR模型的最优C和g
Table 4. Optimal C and g of the SVR model with 8 input feature parameters
|
表 5. 8输入特征SVR模型的预测误差
Table 5. Prediction error of the SVR model with 8 input feature parameters
|
从表5可知,表面粗糙度测量模型的输入特征为TV ={E, S, I, H, Bent, κ, σ, υ}时,在最优超参数下,表面粗糙度测量模型对平磨、卧铣、立铣、研磨测试集样本预测的平均绝对百分比误差分别为4.77%、3.72%、8.16%和3.88%。
此外,输入特征也会影响表面粗糙度测量模型的性能,采用文中改进的序列后向选择算法来解决输入特征的冗余问题。选择5折交叉验证的最小均方误差函数为各SVR模型的性能评价函数,各SVR模型性能评价函数的阈值上限设置为输入特征全集时性能评价函数值的1.1倍。该算法依次从各SVR模型中剔除了冗余特征H、κ和σ。此时,表面粗糙度测量模型的输入特征数为5,记为TW ={E, S, I, Bent, υ},各SVR模型最优的
表 6. 5输入特征SVR模型的最优C和g
Table 6. Optimal C and g of the SVR model with 5 input feature parameters
|
将各加工类型的测试集数据分别输入到相应的SVR模型中,模型对测试集样本表面粗糙度预测的RMSE和MAPE如表7所示。
表 7. 5输入特征SVR模型的预测误差
Table 7. Prediction error of the SVR model with 5 input feature parameters
|
文中改进的序列后向选择算法把各SVR模型的输入特征数从8降低到了5,剔除冗余特征的同时,各SVR模型保留了相同的特征子集。此外,从表5和表7可知,该方法还将平磨、卧铣、立铣、研磨SVR模型的MAPE分别降低了1.22%、0.62%、4.99%、1.61%,模型的RMSE也更低,具有更高的准确性和稳定性。
测试集样本表面粗糙度标准值和预测值如图4所示,预测结果的相对误差如图5所示。从图5可知,平磨0.1~0.8 μm粗糙度等级的样本,表面粗糙度预测的最大相对误差为9.85%;卧铣和立铣0.4~3.2 μm粗糙度等级的样本,表面粗糙度预测的最大相对误差为7.10%和7.95%;研磨0.025~0.1 μm粗糙度等级的样本,表面粗糙度预测的最大相对误差为9.21%,对各加工类型的表面粗糙度的最大相对误差均小于10%。
图 4. 测试集表面粗糙度标准值和预测值。(a)平磨;(b)卧铣;(c)立铣;(d)研磨
Fig. 4. Standard and predicted values of surface roughness for the test set. (a) Plane grinding; (b) Horizontal milling; (c) Vertical milling; (d) Grinding polishing
图 5. 表面粗糙度预测的相对误差图。(a)平磨;(b)卧铣;(c)立铣;(d)研磨
Fig. 5. Relative error graph of surface roughness prediction. (a) Plane grinding; (b) Horizontal milling; (c) Vertical milling; (d) Grinding polishing
前述已用筛选出的特征参数TW ={E, S, I, Bent, υ}建立了测量各加工类型表面粗糙度的模型,由于各模型的支持向量回归函数不同,因此需要确定加工类型后,才能代入相应的SVR模型进行表面粗糙度测量。基于特征参数TW ={E, S, I, Bent, υ}建立了SVM分类器和K最近邻(K-Nearest Neighbor, KNN)分类器。分类器对测试集样本加工类型的预测结果如表8和表9所示,表格的行表示实际类别,列表示预测类别。
表 8. SVM分类器预测结果的混淆矩阵
Table 8. The confusion matrix of the prediction results by the SVM classifier
|
表 9. KNN分类器预测结果的混淆矩阵
Table 9. The confusion matrix of the prediction results by the KNN classifier
|
从表8和表9可知,基于特征参数TW ={E, S, I, Bent, υ}建立的SVM分类器和KNN分类器对测试集样本加工类型的识别率分别为100%和97.78%。由于需要根据识别的加工类型代入相应的SVR模型实现表面粗糙度的测量,因此对加工类型识别率要求较高,通过对比后,文中选择建立的SVM分类器作为加工类型识别的模型。
4 结 论
文中提出了一种解决表面粗糙度建模过程中特征相关和特征冗余问题的方法。通过引入斯皮尔曼相关系数,制定简约规则,可以筛选出与各加工类型表面粗糙度参数均强相关的特征参数,改进的序列后向选择算法,解决了模型的特征冗余问题。实验结果表明,文中提出的方法不仅有效剔除了不相关特征和冗余特征,还可以筛选出一组特征TW ={E, S, I, Bent, υ},建立适用于测量各加工类型表面粗糙度的模型。建立的模型对四种加工类型的识别率为100%,对平磨、卧铣、立铣和研磨表面粗糙度预测的MAPE分别为3.55%、3.10%、3.17%、2.27%,对表面粗糙度有较高的测量精度,与剔除冗余特征前相比,平磨、卧铣、立铣和研磨SVR模型的MAPE分别降低了1.22%、0.62%、4.99%、1.61%,模型的准确性和稳定性更优。
[1] 陈苏婷, 张勇, 胡海锋. 基于激光散斑分形维数的表面粗糙度测量方法[J]. 中国激光, 2015, 42(4): 0408002.
[2] 刘颖, 朗 治 国, 唐 文 彦. 表面粗糙度光切显微镜测量系统的研制[J]. 红外与激光工程, 2012, 41(3): 775-779.
Liu Ying, Lang Zhiguo, Tang Wenyan. Development of measurement system about light-section microscope for surface roughness[J]. Infrared and Laser Engineering, 2012, 41(3): 775-779.
[12] 蒋磊, 刘恒彪, 李同保. 超辐射发光二极管的散斑自相关法表面粗糙度测量研究[J]. 红外与激光工程, 2019, 48(7): 0717003.
[15] Patel D R, Kiran M B. Non-contact surface roughness measurement using laser speckle technique[J]. IOP Conference Series: Materials Science and Engineering, 2020, 895(1): 012007.
[18] Haridas A, Crivoi A, Prabhathan P, et al. Fractal speckle image analysis f surface acterization of aerospace structures[C]International Conference on Optical & Photonics Engineering, 2017, 10449: 329336.
[20] Goodman J W. Speckle Phenomena in Optics: They Applications[M]. US: Roberts Company Publishers, 2014: 911.
Article Outline
吴鹏飞, 邓植中, 雷思琛, 谭振坤, 王姣. 基于激光散斑图像多特征参数的表面粗糙度建模研究[J]. 红外与激光工程, 2023, 52(12): 20230348. Pengfei Wu, Zhizhong Deng, Sichen Lei, Zhenkun Tan, Jiao Wang. Research on surface roughness modeling based on multiple feature parameters of laser speckle image[J]. Infrared and Laser Engineering, 2023, 52(12): 20230348.