CN110728286A - 一种基于火花图像的砂带磨削材料去除率识别方法 - Google Patents

Abstract

本发明公开了一种基于火花图像的砂带磨削材料去除率识别方法，步骤包括：1)使用CCD工业相机采集磨削过程中产生的磨削火花图像；2)对采集的磨削火花图像进行预处理；3)对磨削火花图像进行分割；4)提取磨削火花的特征量化值，包括分别提取磨削火花的亮度特征、颜色特征、面积特征以及轮廓特征；5)基于支持向量回归算法建立材料去除率识别模型。本发明的方法，有效避开砂带磨削参数的复杂耦合关系，并可以有效考虑砂带磨损状态的影响，为砂带磨削控制提供了一种新的思路与方法。

Description

一种基于火花图像的砂带磨削材料去除率识别方法

技术领域

本发明属于砂带磨削精度控制技术领域，涉及一种基于火花图像的砂带磨削材料去除率识别方法。

背景技术

砂带磨削因其加工效率高，成本低而广泛应用于工业领域，砂带磨削材料去除率的控制始终是该领域关注的问题。

传统理论的材料去除率预测模型表达式是：

其中，r表示材料去除率，C_A表示修正系数，K_A表示阻力系数，K_t表示砂带耐用程度，V_s表示砂带线速度，V_w表示工件进给速度，L_ω表示加工区域的宽度，F_A表示法向磨削压力。该材料去除率预测模型从砂带磨削机理出发，通过检测加工过程中的法向磨削力来识别材料去除率，其存在的不足是，一方面，将加工过程中磨削工具磨损简化为常量K_t，忽视了磨削工具磨损状态随时间变化的因素；另一方面，在模型中磨削参数之间呈高度耦合性，难以梳理清楚各参数之间的耦合关系。由于上述的不足，导致砂带磨削控制时应用上述材料去除率预测模型时带来较大的误差，直接影响控制精度。

发明内容

本发明的目的是提供一种基于火花图像的砂带磨削材料去除率识别方法，解决了现有技术采用材料去除率预测模型，结果误差较大，直接影响控制精度，无法准确描述磨削工具磨损状态的问题。

本发明所采用的技术方案是，一种基于火花图像的砂带磨削材料去除率识别方法，按照以下步骤实施：

步骤1：使用CCD工业相机采集磨削过程中产生的磨削火花图像，

该磨削火花图像是RGB模式的彩色图像类型，称为RGB火花图像；

步骤2：对采集的磨削火花图像进行预处理，

对上述的磨削火花图像进行预处理，以过滤掉干扰因素；

步骤3：对磨削火花图像进行分割；

步骤4：提取磨削火花的特征量化值，

磨削火花的特征量化值，包括分别提取磨削火花的亮度特征、颜色特征、面积特征以及轮廓特征；

步骤5：基于支持向量回归算法建立材料去除率识别模型，

5.1)根据步骤4提取出的砂带的磨削火花图像中四种特征的量化数值，建立对应的材料去除率构造数据集，见下式(13)：

D＝{c_1i,c_2i,c_3i,c_4i,M_i}i＝1,2,...,L (13)

其中，i表示样本数，L表示样本的总数量，c₁为亮度特征，c₂表示颜色特征，c₃表示面积特征，c₄表示边缘特征，M表示对应的材料去除率值，对该数据集的数值按组进行归一化处理，具体操作如下式(14)：

其中，d^*表示归一化后的样本数据，d为原始样本数据，d_min为此组样本数据的最小值，d_max为此组数据最大值；对归一化后的数据集重新进行分类，其中80％作为训练集，另外20％作为测试集；

5.2)基于支持向量回归算法训练基于磨削火花图像特征的材料去除率识别模型，确定得到回归模型，最终得到材料去除率识别模型。

本发明的有益效果是，该材料去除率识别方法，根据获取到的不同材料去除率下的砂带磨削火花图像，首先进行预处理以滤去背景中的砂带机等干扰因素，然后提取出亮度、颜色、面积、轮廓四种特征的量化数值，并构造数据集，利用SVR算法对各个材料去除率所对应的四种特征量化值的数据进行训练，基于训练所得的模型进行材料去除率的识别。本发明通过砂带磨削过程中产生的火花图像特征建立的材料去除率识别模型，一方面该方法可以有效避开复杂的磨削参数耦合关系；另一方面通过火花图像特征能够有效反映磨削工具不同磨损状态下材料去除率的差异，能够准确描述磨削工具磨损状态的实际。

附图说明

图1为本发明方法根据砂带磨削火花图像建立材料去除率识别模型的整体流程图；

图2为本发明方法根据砂带磨削火花图像四种特征提取量化的流程图；

图3为本发明方法基于磨削火花图像特征建立材料去除率识别模型的具体过程流程图。

具体实施方式

下面结合附图和具体实施方式对本发明作进一步详细说明。

参照图1，本发明基于火花图像的砂带磨削材料去除率识别方法，按照以下步骤实施：

步骤1：使用CCD工业相机采集磨削过程中产生的磨削火花图像；

采集到的磨削火花图像是RGB模式的彩色图像类型，称为RGB火花图像，即“RGB火花图像”表示相机采集的未经任何处理的图像；

步骤2：对采集的磨削火花图像进行预处理，

由于CCD工业相机所拍摄到的磨削火花图像包含砂带机等干扰因素，会对特征提取产生影响，因此需要对上述的磨削火花图像进行预处理，以过滤掉干扰因素，具体过程如下：

对RGB火花图像进行灰度值化处理，如下式(2)：

Gray(i,j)＝0.11·R(i,j)+0.59·G(i,j)+0.3·B(i,j) (2)

其中，Gray(i,j)表示转换得到的灰度图像，R(i,j)表示RGB火花图像中像素点(i,j)的红色分量，同样，G(i,j)表示绿色分量，B(i,j)表示蓝色分量；

步骤3：对磨削火花图像进行分割，

对采集到的RGB火花图像进行逐像素操作，判断与该像素对应的步骤2所得的灰度图像的像素灰度值与设定阈值的大小，若小于阈值，则将RGB火花图像中该像素值设为0；否则，不进行任何操作；

完成对RGB火花图像中所有像素的上述操作，获得单一背景的RGB火花图像；

该“单一背景的RGB火花图像”表示经过预处理的，将火花从背景中分离出来的、只包含火花的RGB图像。

步骤4：提取磨削火花的特征量化值，

参照图2，磨削火花的特征量化值，主要包括提取磨削火花的亮度特征、颜色特征、面积特征以及轮廓特征，具体过程是，

4.1)提取磨削火花的亮度特征，采取如下方式：

将步骤3预处理后得到的单一背景的RGB火花图像转换到HSV模型下，表达式如下：

v＝max (5)

其中，max表示某像素RGB分量中的最大值，min表示某像素RGB分量中的最小值，h表示HSV模型下的色调，s表示HSV模型下的饱和度，v表示HSV模型下的亮度；

统计HSV模型中亮度分量的像素个数，即表示该材料去除率下磨削火花图像的亮度特征。

4.2)提取磨削火花的颜色特征，采取如下方式：

其中，p_ij表示单一背景的RGB火花图像第j个像素的第i个颜色分量，N表示单一背景的RGB火花图像的像素总数，μ_i表示单一背景的RGB火花图像中三个分量的平均值，作为该材料去除率下磨削火花图像的颜色特征值。

4.3)提取磨削火花的面积特征，采取如下方式：

先对单一背景的RGB火花图像采用最大类间方差法进行阈值分割，再根据确定的阈值将单一背景的RGB火花图像进行二值化，得到像素值为0和1的二值图像，统计该二值图像中像素值为1的像素个数，即表征该材料去除率下磨削火花图像的面积特征。

4.4)提取磨削火花的轮廓特征，采取如下方式：

首先，对单一背景的RGB火花图像进行高斯滤波降噪处理，如下式：

其中，H(x,y,σ)是高斯滤波函数，σ为平滑系数，(x,y)表示像素的坐标，即行号和列号；式(7)中的(x,y)表示输出图像像素坐标，(x₀,y₀)表示输入(即单一背景的RGB火花图像)像素坐标，f(x,y)为单一背景的RGB火花图像描述函数，然后，对滤波降噪后的图像P(x,y)采用Sobel算子求取梯度幅值与方向角，以获得火花的边缘图像G(x,y)，如下式：

其中，G_x(x,y)表示像素(x,y)的水平梯度值；G_y(x,y)表示像素(x,y)的垂直梯度值；G(x,y)表示像素(x,y)的梯度幅值，也可以表示输出的边缘图像；T(x,y)表示像素(x y)的方向角，

对边缘图像G(x,y)进行二值化，对梯度幅值大于设定阈值的像素值设为1，梯度幅值小于设定阈值的像素值设为0；统计二值化后的边缘图像中像素值为1的像素点数量，即表征该材料去除率下磨削火花图像的轮廓特征。

至此，提取出了砂带的磨削火花图像中四种特征的量化数值，接下来就可以建立材料去除率预测模型。

步骤5：基于支持向量回归算法建立材料去除率识别模型，

5.1)参照图3，根据步骤4提取出的砂带的磨削火花图像中四种特征的量化数值，建立对应的材料去除率构造数据集，见下式(13)：

D＝{c_1i,c_2i,c_3i,c_4i,M_i}i＝1,2,...,L (13)

其中，i表示样本数，L表示样本的总数量，c₁为亮度特征，c₂表示颜色特征，c₃表示面积特征，c₄表示边缘特征，M表示对应的材料去除率值，对该数据集的数值按组进行归一化处理，具体操作如下式(14)：

其中，d^*表示归一化后的样本数据，d为原始样本数据，d_min为此组样本数据的最小值，d_max为此组数据最大值；对归一化后的数据集重新进行分类，其中80％作为训练集，另外20％作为测试集；

5.2)基于支持向量回归算法训练基于磨削火花图像特征的材料去除率识别模型，具体过程如下：

5.2.1)根据支持向量回归算法，使用核函数的多特征输入的支持向量回归模型表示为：

Y＝w·φ(X)+b (15)

其中，w与b为模型的待求参数，Y为输出向量，即材料去除率M，X为输入特征向量{c₁,c₂,c₃,c₄}，φ(X)指通过核函数将输入特征映射到高维空间的特征向量；

在此，使用径向基核函数(也称高斯核函数)：

其中，δ为核函数的超参数，δ优选取值为0.25，d_i表示归一化的样本数据，d_c为核函数的中心；

5.2.2)求解最优化问题，如下式(17)：

其中，C表示惩罚因子，惩罚因子C优选取值为1，ξ_i和ξ_i'分别表示第i个样本的上下松弛变量，L为样本的总数，ε表示模型的容差值，s.t.表示服从于后面表达式；

5.2.3)使用拉格朗日乘子和拉格朗日对偶性质，则将式(17)转换为：

其中，α_i、α_i'、μ、μ'、μ_i、μ_i'是拉格朗日系数，在计算过程中迭代求取，函数L_g表示回归估计函数；最后结合序列最小优化算法对式(18)进行求解，获得使(18)值最小的参数的组合，求得满足式(18)的参数的值即为回归模型的参数，确定得到回归模型，即基于火花图像特征的材料去除率识别模型。(此处回归模型就是材料去除率识别模型)

本发明的方法，首先提取砂带磨削过程产生的火花图像特征，如面积特征、亮度特征、颜色特征、轮廓特征，再通过这些特征建立材料去除率识别模型，直接判断出砂带磨削过程中的材料去除率，有效避免砂带磨削理论模型中各因素的耦合影响和砂带磨损的影响，具有识别准确、效率高的优势。

Claims (6)

1.一种基于火花图像的砂带磨削材料去除率识别方法，其特征在于，按照以下步骤实施：

步骤1：使用CCD工业相机采集磨削过程中产生的磨削火花图像，

该磨削火花图像是RGB模式的彩色图像类型，称为RGB火花图像；

步骤2：对采集的磨削火花图像进行预处理，

对上述的磨削火花图像进行预处理，以过滤掉干扰因素；

步骤3：对磨削火花图像进行分割；

步骤4：提取磨削火花的特征量化值，

磨削火花的特征量化值，包括分别提取磨削火花的亮度特征、颜色特征、面积特征以及轮廓特征；

步骤5：基于支持向量回归算法建立材料去除率识别模型，

5.1)根据步骤4提取出的砂带的磨削火花图像中四种特征的量化数值，建立对应的材料去除率构造数据集，见下式(13)：

D＝{c_1i,c_2i,c_3i,c_4i,M_i}i＝1,2,...,L (13)

其中，i表示样本数，L表示样本的总数量，c₁为亮度特征，c₂表示颜色特征，c₃表示面积特征，c₄表示边缘特征，M表示对应的材料去除率值，对该数据集的数值按组进行归一化处理，具体操作如下式(14)：

其中，d^*表示归一化后的样本数据，d为原始样本数据，d_min为此组样本数据的最小值，d_max为此组数据最大值；对归一化后的数据集重新进行分类，其中80％作为训练集，另外20％作为测试集；

5.2)基于支持向量回归算法训练基于磨削火花图像特征的材料去除率识别模型，确定得到回归模型，即得到材料去除率识别模型。

2.根据权利要求1所述的基于火花图像的砂带磨削材料去除率识别方法，其特征在于：所述步骤2中，预处理的具体过程如下，

对RGB火花图像进行灰度值化处理，表达式如下式(2)：

Gray(i,j)＝0.11·R(i,j)+0.59·G(i,j)+0.3·B(i,j) (2)

其中，Gray(i,j)表示转换得到的灰度图像，R(i,j)表示RGB火花图像中像素点(i,j)的红色分量，同样，G(i,j)表示绿色分量，B(i,j)表示蓝色分量。

3.根据权利要求1所述的基于火花图像的砂带磨削材料去除率识别方法，其特征在于：所述步骤3中，对采集到的RGB火花图像进行逐像素操作，判断与该像素对应的步骤2所得的灰度图像的像素灰度值与设定阈值的大小，若小于阈值，则将RGB火花图像中该像素值设为0；否则，不进行任何操作；

完成对RGB火花图像中所有像素的上述操作，获得单一背景的RGB火花图像。

4.根据权利要求1所述的基于火花图像的砂带磨削材料去除率识别方法，其特征在于：所述步骤4中，具体过程是，

4.1)提取磨削火花的亮度特征，采取如下方式：

将步骤3预处理后得到的单一背景的RGB火花图像转换到HSV模型下，表达式如下：

v＝max (5)

其中，max表示某像素RGB分量中的最大值，min表示某像素RGB分量中的最小值，h表示HSV模型下的色调，s表示HSV模型下的饱和度，v表示HSV模型下的亮度；

统计HSV模型中亮度分量的像素个数，即表示该材料去除率下磨削火花图像的亮度特征；

4.2)提取磨削火花的颜色特征，采取如下方式：

其中，p_ij表示单一背景的RGB火花图像第j个像素的第i个颜色分量，N表示单一背景的RGB火花图像的像素总数，μ_i表示单一背景的RGB火花图像中三个分量的平均值，作为该材料去除率下磨削火花图像的颜色特征值；

4.3)提取磨削火花的面积特征，采取如下方式：

先对单一背景的RGB火花图像采用最大类间方差法进行阈值分割，再根据确定的阈值将单一背景的RGB火花图像进行二值化，得到像素值为0和1的二值图像，统计该二值图像中像素值为1的像素个数，即表征该材料去除率下磨削火花图像的面积特征；

4.4)提取磨削火花的轮廓特征，采取如下方式：

首先，对单一背景的RGB火花图像进行高斯滤波降噪处理，如下式：

其中，H(x,y,σ)是高斯滤波函数，σ为平滑系数，(x,y)表示像素的坐标，即行号和列号；式(7)中的(x,y)表示输出图像像素坐标，(x₀,y₀)表示输入像素坐标，f(x,y)为单一背景的RGB火花图像描述函数；

然后，对滤波降噪后的图像P(x,y)采用Sobel算子求取梯度幅值与方向角，以获得火花的边缘图像G(x,y)，如下式：

其中，G_x(x,y)表示像素(x,y)的水平梯度值；G_y(x,y)表示像素(x,y)的垂直梯度值；G(x,y)表示像素(x,y)的梯度幅值，也可以表示输出的边缘图像；T(x,y)表示像素(x y)的方向角，

对边缘图像G(x,y)进行二值化，对梯度幅值大于设定阈值的像素值设为1，梯度幅值小于设定阈值的像素值设为0；统计二值化后的边缘图像中像素值为1的像素点数量，即表征该材料去除率下磨削火花图像的轮廓特征。

5.根据权利要求1所述的基于火花图像的砂带磨削材料去除率识别方法，其特征在于：所述步骤5.2)中，具体过程如下：

5.2.1)根据支持向量回归算法，使用核函数的多特征输入的支持向量回归模型表示为：

Y＝w·φ(X)+b (15)

其中，w与b为模型的待求参数，Y为输出向量，即材料去除率M，X为输入特征向量{c₁,c₂,c₃,c₄}，φ(X)指通过核函数将输入特征映射到高维空间的特征向量；

使用径向基核函数：

其中，δ为核函数的超参数，d_i表示归一化的样本数据，d_c为核函数的中心；

5.2.2)求解最优化问题，如下式(17)：

其中，C表示惩罚因子，ξ_i和ξ_i'分别表示第i个样本的上下松弛变量，L为样本的总数，ε表示模型的容差值，s.t.表示服从于后面表达式；

5.2.3)使用拉格朗日乘子和拉格朗日对偶性质，则将式(17)转换为：

其中，α_i、α_i'、μ、μ'、μ_i、μ_i'是拉格朗日系数，在计算过程中迭代求取，函数L_g表示回归估计函数；最后结合序列最小优化算法对式(18)进行求解，获得使(18)值最小的参数的组合，求得满足式(18)的参数的值即为回归模型的参数，确定得到回归模型。

6.根据权利要求5所述的基于火花图像的砂带磨削材料去除率识别方法，其特征在于：所述核函数的超参数δ取值为0.25，惩罚因子C取值为1。

Images