一种润滑油磨粒质量检测方法

摘要

本发明涉及图像处理技术领域，具体涉及一种润滑油磨粒质量检测方法。方法包括：获取磨粒图像中的各磨粒区域的长宽比；根据各磨粒区域的轮廓像素点，得到弯曲度；根据各磨粒区域的最小拟合圆的面积，得到球状度；根据各磨粒区域的灰度共生矩阵，得到粗糙度；根据各磨粒区域中的像素点，得到空洞指标；根据K‑means算法以及各磨粒区域的长宽比、弯曲度、球状度、粗糙度和空洞指标，得到各磨粒类别对应的最优簇；根据各最优簇对应的各磨粒区域的面积，得到润滑油磨粒质量指标；若润滑油磨粒质量指标大于质量阈值，则判定磨粒的质量良好。本发明能够更加准确的判断润滑油中磨粒的质量。

权利要求

1.一种润滑油磨粒质量检测方法，其特征在于，该方法包括以下步骤：

获取待检测润滑油样本对应的磨粒图像，所述润滑油样本是从机器中取样得到的，所述磨粒图像包含多个磨粒区域；

获取各磨粒区域对应的长宽比；根据各磨粒区域对应的轮廓像素点，得到各磨粒区域对应的弯曲度，所述磨粒区域对应的轮廓像素点为磨粒区域最外层的像素点；根据各磨粒区域对应的最小拟合圆的面积，得到磨粒区域对应的球状度；根据各磨粒区域对应的灰度共生矩阵，得到各磨粒区域对应的粗糙度；根据各磨粒区域中的像素点，得到各磨粒区域对应的空洞指标；

所述空洞指标的获取方法为：对于任一磨粒区域：使用canny算子对该磨粒区域进行边缘检测，得到该磨粒区域中存在的各边缘；将边缘上的像素点记为边缘像素点；

对于该磨粒区域内的任一边缘：判断该边缘上的边缘像素点对应的八邻域内存在的边缘像素点数量是否大于等于2，若是，则将对应的边缘像素点标记为0，若否，则将对应的边缘像素点标记为1；计算该边缘上的各边缘像素点的标记值的累加和，若累加和为0，则判定该边缘为闭合边缘，若累加和不为0，则判定该边缘不是闭合边缘；

若该磨粒区域内存在闭合边缘，则将该磨粒区域对应的空洞指标赋值为1；若该磨粒区域内不存在闭合边缘，则将该磨粒区域对应的空洞指标赋值为0；

根据磨粒类别数量、K‑means算法以及各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，对各磨粒区域进行分类，得到各磨粒类别对应的最优簇；

所述得到各磨粒类别对应的最优簇的方法为：根据各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，将各磨粒区域映射到多维空间当中，得到各磨粒区域在多维空间中对应的磨粒点；

在多维空间中获取设定数量个初始质心；所述设定数量为磨粒类别数量；

根据各磨粒点对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，计算各磨粒点与各初始质心之间的特征距离；

基于特征距离，利用K‑means聚类算法对各磨粒点进行聚类，得到各磨粒类别对应的初始簇；

获取各初始簇对应的聚类中心的长宽比、弯曲度、球状度、粗糙度和空洞指标；

计算各初始簇中各磨粒点与对应的聚类中心之间的特征距离；计算各初始簇中各磨粒点与对应的聚类中心之间的特征距离的和，作为第一次聚类的簇内差异指标；

根据各初始簇对应的聚类中心的长宽比、弯曲度、球状度、粗糙度和空洞指标，得到第一次聚类的簇间差异指标；

获取正常滑动磨粒与球形磨粒对应的初始簇；根据正常滑动磨粒与球形磨粒对应的初始簇，得到正常滑动磨粒和球形磨粒对应的初始簇的轮廓评价指标和中心区域评价指标；

其中，所述得到正常滑动磨粒和球形磨粒对应的初始簇的轮廓评价指标和中心区域评价指标的方法为：对于正常滑动磨粒和球形磨粒中任一磨粒类别对应的初始簇：

对于该磨粒类别对应的初始簇对应的任一磨粒区域：获取该磨粒区域对应的最小拟合圆，以及最小拟合圆的圆心，将所述圆心记为该磨粒区域对应的轮廓中心；

根据该磨粒类别对应的初始簇对应的各磨粒区域的各轮廓像素点与对应磨粒区域的轮廓中心的欧氏距离，计算该磨粒类别对应的初始簇的轮廓评价指标；

对于该磨粒类别对应的初始簇对应的任一磨粒区域：以该磨粒区域的最小拟合圆的圆心为目标圆心，以最小拟合圆的半径的一半为目标半径作圆，得到目标圆形区域；计算该磨粒区域的目标圆形区域内各像素点的灰度值的平均值，作为该磨粒区域的中心灰度均值；

计算该磨粒类别对应的初始簇对应的各磨粒区域的中心灰度均值的平均值，作为该磨粒类别对应的初始簇的中心区域评价指标；

所述轮廓评价指标的计算公式为：

其中，F为该磨粒类别对应的初始簇的轮廓评价指标，U为该磨粒类别对应的初始簇对应的各磨粒区域的数量， 为该磨粒类别对应的初始簇对应的第u个磨粒区域的轮廓像素点的数量， 为该磨粒类别对应的初始簇对应的第u个磨粒区域的平滑指标， 为所述第u个磨粒区域对应的第i个轮廓像素点与所述第u个磨粒区域对应的轮廓中心的欧氏距离，为所述第u个磨粒区域对应的第i+1个轮廓像素点与所述第u个磨粒区域对应的轮廓中心的欧氏距离， 为取绝对值函数；

根据所述轮廓评价指标和中心区域评价指标，得到第一次聚类的聚类评价因子；所述正常滑动磨粒和球形磨粒为各磨粒类别中的两种磨粒类别；

所述聚类评价因子的计算公式为：

其中，为球形磨粒对应的初始簇的轮廓评价指标， 为球形磨粒对应的初始簇的中心区域评价指标， 为正常滑动磨粒对应的初始簇的轮廓评价指标， 为正常滑动磨粒对应的初始簇的中心区域评价指标；

根据第一次聚类的簇内差异指标、簇间差异指标以及聚类评价因子，得到第一次聚类的聚类效果指标；

判断第一次聚类的聚类效果指标是否大于效果阈值，若是，则将第一次聚类得到各磨粒类别对应的初始簇记为最优簇；若否，则从多维空间中重新选取设定数量个初始质心，并对多维空间中的各磨粒点进行第二次聚类，得到各磨粒类别对应的第二初始簇；计算第二次聚类的聚类效果指标，判断第二次聚类的聚类效果指标是否大于效果阈值，若是，则将各磨粒类别对应的第二初始簇记为最优簇，若否，则继续重新选取设定数量个初始质心，对多维空间中各磨粒点进行重新聚类，直到聚类效果指标第一次大于效果阈值为止；将最后一次聚类得到的各磨粒类别对应的簇记为最优簇；

所述聚类效果指标的计算公式如下：

其中，P为聚类效果指标，p1为簇内差异指标，p2为簇间差异指标，p3为聚类评价因子，exp( )为以e为底的指数函数；

根据各磨粒类别对应的最优簇对应的各磨粒区域的面积，得到润滑油磨粒质量指标；

若润滑油磨粒质量指标大于质量阈值，则判定待检测润滑油样本中磨粒的质量良好。

2.根据权利要求1所述的一种润滑油磨粒质量检测方法，其特征在于，所述得到各磨粒区域对应的弯曲度，包括：对于任一磨粒区域：

对于该磨粒区域对应的轮廓上的任一轮廓像素点：获取该磨粒区域对应的轮廓上与该轮廓像素点的直线距离最远的轮廓像素点，将所述与该轮廓像素点距离最远的轮廓像素点记为该轮廓像素点对应的目标像素点；利用弗洛伊德算法获取该轮廓像素点与对应的目标像素点之间的内距离；所述内距离为该轮廓像素点与对应的目标像素点在该磨粒区域内部的最短路径的长度；

计算该磨粒区域对应的轮廓上的各轮廓像素点与对应的目标像素点的直线距离与对应的内距离的比值的和，作为该磨粒区域对应的弯曲度。

3.根据权利要求1所述的一种润滑油磨粒质量检测方法，其特征在于，所述得到磨粒区域对应的球状度，包括：对于任一磨粒区域：计算该磨粒区域的面积与对应的最小拟合圆的面积的比值，作为该磨粒区域对应的球状度。

4.根据权利要求1所述的一种润滑油磨粒质量检测方法，其特征在于，所述得到各磨粒区域对应的粗糙度，包括：对于任一磨粒区域：

根据该磨粒区域对应的灰度共生矩阵，计算该磨粒区域对应的灰度共生矩阵对应的熵和逆方差；

根据所述熵和逆方差，计算该磨粒区域对应的粗糙度；所述计算该磨粒区域对应的粗糙度的公式如下：其中，R为该磨粒区域对应的粗糙度，Ent为该磨粒区域对应的灰度共生矩阵的熵，IDM为该磨粒区域对应的灰度共生矩阵的逆方差，exp( )为以e为底的指数函数。

5.根据权利要求1所述的一种润滑油磨粒质量检测方法，其特征在于，计算各磨粒点与各初始质心之间的特征距离的公式如下：其中， 为第j个磨粒点与第q个初始质心之间的特征距离， 为第j个磨粒点对应的长宽比， 为第j个磨粒点对应的弯曲度， 为第j个磨粒点对应的球状度， 为第j个磨粒点对应的粗糙度， 为第j个磨粒点对应的空洞指标， 为第q个初始质心对应的长宽比，为第q个初始质心对应的弯曲度， 为第q个初始质心对应的球状度， 为第q个初始质心对应的粗糙度， 为第q个初始质心对应的空洞指标。

6.根据权利要求1所述的一种润滑油磨粒质量检测方法，其特征在于，所述簇间差异指标的计算公式如下：其中， 为第n1个初始簇对应的聚类中心与第n2个初始簇对应的聚类中心之间的特征距离，K为磨粒类别的数量。

7.根据权利要求1所述的一种润滑油磨粒质量检测方法，其特征在于，得到润滑油磨粒质量指标的公式如下：其中，V为润滑油磨粒质量指标， 为第 个磨粒类别对润滑油质量的影响程度， 为第 个磨粒类别对应的最优簇对应的所有磨粒区域的面积之和，exp( )为以e为底的指数函数。

说明书

一种润滑油磨粒质量检测方法

技术领域

[0001] 本发明涉及图像处理技术领域，具体涉及一种润滑油磨粒质量检测方法。

背景技术

[0002] 在机械系统中，例如齿轮系统等重要部件都是浸入在润滑油系统中的，润滑油可以缓解机械设备各元器件之间的磨损程度以及因摩擦所引起的不必要的能源损耗，增加机器的使用寿命。它能够缓冲器件之间的相互摩擦，减少磨损的发生，它的重要性就如同是“血液”之于人体，是不可替代的。考虑到设备元器件之间的互相摩擦产生的磨损颗粒会悬浮在润滑油中，它会混淆在润滑油中随润滑油在润滑系统中不断地循环，会影响到润滑油的润滑效果。因此，需定期对润滑油的磨粒质量进行检测来判断润滑油的健康状态。

[0003] 现有的往往通过传感器对润滑油中的磨粒进行检测，但是传感器成本较高，并且只能给出油液中磨粒浓度的信息，仅能区别铁磁性颗粒和非铁磁性颗粒，无法确定磨粒的具体类型；由于不同种类的磨粒所产生的影响是不同的，因此通过传感器难以准确的对润滑油中磨粒的质量进行检测。

发明内容

[0004] 为了解决上述技术问题，本发明的目的在于提供一种润滑油磨粒质量检测方法，所采用的技术方案具体如下：

[0005] 本发明提供了一种润滑油磨粒质量检测方法包括以下步骤：

[0006] 获取待检测润滑油样本对应的磨粒图像，所述润滑油样本是从机器中取样得到的，所述磨粒图像包含多个磨粒区域；

[0007] 获取各磨粒区域对应的长宽比；根据各磨粒区域对应的轮廓像素点，得到各磨粒区域对应的弯曲度，所述磨粒区域对应的轮廓像素点为磨粒区域最外层的像素点；根据各磨粒区域对应的最小拟合圆的面积，得到磨粒区域对应的球状度；根据各磨粒区域对应的灰度共生矩阵，得到各磨粒区域对应的粗糙度；根据各磨粒区域中的像素点，得到各磨粒区域对应的空洞指标；

[0008] 根据磨粒类别数量、K‑means算法以及各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，对各磨粒区域进行分类，得到各磨粒类别对应的最优簇；

[0009] 根据各磨粒类别对应的最优簇对应的各磨粒区域的面积，得到润滑油磨粒质量指标；若润滑油磨粒质量指标大于质量阈值，则判定待检测润滑油样本中磨粒的质量良好。

[0010] 优选的，所述得到各磨粒区域对应的弯曲度，包括：

[0011] 对于任一磨粒区域：

[0012] 对于该磨粒区域对应的轮廓上的任一轮廓像素点：获取该磨粒区域对应的轮廓上与该轮廓像素点的直线距离最远的轮廓像素点，将所述与该轮廓像素点距离最远的轮廓像素点记为该轮廓像素点对应的目标像素点；利用弗洛伊德算法获取该轮廓像素点与对应的目标像素点之间的内距离；所述内距离为该轮廓像素点与对应的目标像素点在该磨粒区域内部的最短路径的长度；

[0013] 计算该磨粒区域对应的轮廓上的各轮廓像素点与对应的目标像素点的直线距离与对应的内距离的比值的和，作为该磨粒区域对应的弯曲度。

[0014] 优选的，所述得到磨粒区域对应的球状度，包括：

[0015] 对于任一磨粒区域：计算该磨粒区域的面积与对应的最小拟合圆的面积的比值，作为该磨粒区域对应的球状度。

[0016] 优选的，所述得到各磨粒区域对应的粗糙度，包括：

[0017] 对于任一磨粒区域：

[0018] 根据该磨粒区域对应的灰度共生矩阵，计算该磨粒区域对应的灰度共生矩阵对应的熵和逆方差；

[0019] 根据所述熵和逆方差，计算该磨粒区域对应的粗糙度；所述计算该磨粒区域对应的粗糙度的公式如下：

[0020]

[0021] 其中，R为该磨粒区域对应的粗糙度，Ent为该磨粒区域对应的灰度共生矩阵的熵，IDM为该磨粒区域对应的灰度共生矩阵的逆方差，exp( )为以e为底的指数函数。

[0022] 优选的，所述得到各磨粒区域对应的空洞指标，包括：

[0023] 对于任一磨粒区域：

[0024] 使用canny算子对该磨粒区域进行边缘检测，得到该磨粒区域中存在的各边缘；将边缘上的像素点记为边缘像素点；

[0025] 对于该磨粒区域内的任一边缘：判断该边缘上的边缘像素点对应的八邻域内存在的边缘像素点数量是否大于等于2，若是，则将对应的边缘像素点标记为0，若否，则将对应的边缘像素点标记为1；计算该边缘上的各边缘像素点的标记值的累加和，若累加和为0，则判定该边缘为闭合边缘，若累加和不为0，则判定该边缘不是闭合边缘；

[0026] 若该磨粒区域内存在闭合边缘，则将该磨粒区域对应的空洞指标赋值为1；若该磨粒区域内不存在闭合边缘，则将该磨粒区域对应的空洞指标赋值为0。

[0027] 优选的，根据磨粒类别数量、K‑means算法以及各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，对各磨粒区域进行分类，得到各磨粒类别对应的最优簇，包括：

[0028] 根据各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，将各磨粒区域映射到多维空间当中，得到各磨粒区域在多维空间中对应的磨粒点；

[0029] 在多维空间中获取设定数量个初始质心；所述设定数量为磨粒类别数量；

[0030] 根据各磨粒点对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，计算各磨粒点与各初始质心之间的特征距离；

[0031] 基于各磨粒点与各初始质心之间的特征距离，利用K‑means聚类算法对各磨粒点进行聚类，得到各磨粒类别对应的初始簇；

[0032] 获取各初始簇对应的聚类中心的长宽比、弯曲度、球状度、粗糙度和空洞指标；

[0033] 计算各初始簇中各磨粒点与对应的聚类中心之间的特征距离；计算各初始簇中各磨粒点与对应的聚类中心之间的特征距离的和，作为第一次聚类的簇内差异指标；

[0034] 根据各初始簇对应的聚类中心的长宽比、弯曲度、球状度、粗糙度和空洞指标，得到第一次聚类的簇间差异指标；

[0035] 获取正常滑动磨粒与球形磨粒对应的初始簇；根据正常滑动磨粒与球形磨粒对应的初始簇，得到正常滑动磨粒和球形磨粒对应的初始簇的轮廓评价指标和中心区域评价指标；根据所述轮廓评价指标和中心区域评价指标，得到第一次聚类的聚类评价因子；所述正常滑动磨粒和球形磨粒为各磨粒类别中的两种磨粒类别；

[0036] 根据第一次聚类的簇内差异指标、簇间差异指标以及聚类评价因子，得到第一次聚类的聚类效果指标；

[0037] 判断第一次聚类的聚类效果指标是否大于效果阈值，若是，则将第一次聚类得到各磨粒类别对应的初始簇记为最优簇；若否，则从多维空间中重新选取设定数量个初始质心，并对多维空间中的各磨粒点进行第二次聚类，得到各磨粒类别对应的第二初始簇；计算第二次聚类的聚类效果指标，判断第二次聚类的聚类效果指标是否大于效果阈值，若是，则将各磨粒类别对应的第二初始簇记为最优簇，若否，则继续重新选取设定数量个初始质心，对多维空间中各磨粒点进行重新聚类，直到聚类效果指标第一次大于效果阈值为止；将最后一次聚类得到的各磨粒类别对应的簇记为最优簇；

[0038] 所述聚类效果指标的计算公式如下：

[0039]

[0040] 其中，P为聚类效果指标，p1为簇内差异指标，p2为簇间差异指标，p3为聚类评价因子，exp( )为以e为底的指数函数。

[0041] 优选的，计算各磨粒点与各初始质心之间的特征距离的公式如下：

[0042]

[0043] 其中， 为第j个磨粒点与第q个初始质心之间的特征距离， 为第j个磨粒点对应的长宽比， 为第j个磨粒点对应的弯曲度， 为第j个磨粒点对应的球状度， 为第j个磨粒点对应的粗糙度， 为第j个磨粒点对应的空洞指标， 为第q个初始质心对应的长宽比， 为第q个初始质心对应的弯曲度， 为第q个初始质心对应的球状度， 为第q个初始质心对应的粗糙度， 为第q个初始质心对应的空洞指标。

[0044] 优选的，所述簇间差异指标的计算公式如下：

[0045]

[0046] 其中， 为第n1个初始簇对应的聚类中心与第n2个初始簇对应的聚类中心之间的特征距离，K为磨粒类别的数量；

[0047] 所述聚类评价因子的计算公式如下：

[0048]

[0049] 其中，为球形磨粒对应的初始簇的轮廓评价指标， 为球形磨粒对应的初始簇的中心区域评价指标， 为正常滑动磨粒对应的初始簇的轮廓评价指标， 为正常滑动磨粒对应的初始簇的中心区域评价指标。

[0050] 优选的，根据正常滑动磨粒对应的初始簇与球形磨粒对应的初始簇，得到正常滑动磨粒和球形磨粒对应的初始簇的轮廓评价指标和中心区域评价指标，包括：

[0051] 对于正常滑动磨粒和球形磨粒中任一磨粒类别对应的初始簇：

[0052] 对于该磨粒类别对应的初始簇对应的任一磨粒区域：获取该磨粒区域对应的最小拟合圆，以及最小拟合圆的圆心，将所述圆心记为该磨粒区域对应的轮廓中心；

[0053] 根据该磨粒类别对应的初始簇对应的各磨粒区域的各轮廓像素点与对应磨粒区域的轮廓中心的欧氏距离，计算该磨粒类别对应的初始簇的轮廓评价指标；

[0054] 对于该磨粒类别对应的初始簇对应的任一磨粒区域：以该磨粒区域的最小拟合圆的圆心为目标圆心，以最小拟合圆的半径的一半为目标半径作圆，得到目标圆形区域；计算该磨粒区域的目标圆形区域内各像素点的灰度值的平均值，作为该磨粒区域的中心灰度均值；

[0055] 计算该磨粒类别对应的初始簇对应的各磨粒区域的中心灰度均值的平均值，作为该磨粒类别对应的初始簇的中心区域评价指标；

[0056] 所述轮廓评价指标的计算公式为：

[0057]

[0058]

[0059] 其中，F为该磨粒类别对应的初始簇的轮廓评价指标，U为该磨粒类别对应的初始簇对应的各磨粒区域的数量， 为该磨粒类别对应的初始簇对应的第u个磨粒区域的轮廓像素点的数量， 为该磨粒类别对应的初始簇对应的第u个磨粒区域的平滑指标， 为所述第u个磨粒区域对应的第i个轮廓像素点与所述第u个磨粒区域对应的轮廓中心的欧氏距离， 为所述第u个磨粒区域对应的第i+1个轮廓像素点与所述第u个磨粒区域对应的轮廓中心的欧氏距离， 为取绝对值函数。

[0060] 优选的，得到润滑油磨粒质量指标的公式如下：

[0061]

[0062] 其中，V为润滑油磨粒质量指标， 为第 个磨粒类别对润滑油质量的影响程度，为第 个磨粒类别对应的最优簇对应的所有磨粒区域的面积之和，exp( )为以e为底的指数函数。

[0063] 本发明具有如下有益效果：

[0064] 本发明提供了一种润滑油磨粒质量检测方法，用于对润滑油中的磨粒进行分析；本发明依据不同磨粒类别的特点，获取磨粒图像中各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度以及空洞指标；本发明根据各磨粒区域对应的特征能够更好的将不同磨粒类别的磨粒区域进行区分，便于更好的判断润滑油的质量。本发明基于磨粒类别数量、K‑means算法以及各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，对各磨粒区域进行分类，得到各磨粒类别对应的最优簇；本发明在获取最优簇的过程中，不断变换初始质心进行重新聚类，并分析每一次聚类的聚类效果，以得到聚类效果好的最优簇，避免了由于初始质心选取的随机性而造成的聚类效果不理想的情况的发生。本发明结合不同磨粒类别对润滑油的不同影响，得到润滑油磨粒质量指标；最后根据润滑油磨粒质量指标判断润滑油中磨粒的质量的好坏。本发明基于不同磨粒类别的特征，并利用变换初始质心的K‑means算法，更加准确对图像中的磨粒进行分类，使本发明能够更加准确的对润滑油中磨粒的质量进行判断。

附图说明

[0065] 为了更清楚地说明本发明实施例或现有技术中的技术方案和优点，下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它附图。

[0066] 图1为本发明提供的一种润滑油磨粒质量检测方法的流程图。

具体实施方式

[0067] 为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功能效果，以下结合附图及较佳实施例，对依据本发明提出的一种润滑油磨粒质量检测方法进行详细说明如下。

[0068] 除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。

[0069] 下面结合附图具体的说明本发明所提供的一种润滑油磨粒质量检测方法的具体方案。

[0070] 一种润滑油磨粒质量检测方法实施例：

[0071] 如图1所示，本实施例的一种润滑油磨粒质量检测方法包括以下步骤：

[0072] 步骤S1，获取待检测润滑油样本对应的磨粒图像，所述润滑油样本是从机器中取样得到的，所述磨粒图像包含多个磨粒区域。

[0073] 在机械系统中，例如齿轮系统等重要部件都是浸入在润滑油系统中的，润滑油可以缓解机械设备各元器件之间的磨损程度以及因摩擦所引起的不必要的能源损耗，增加机器的使用寿命。它能够缓冲器件之间的相互摩擦，减少磨损的发生。考虑到设备元器件之间的互相摩擦产生的磨损颗粒就悬浮在润滑油中，它会混淆在润滑油中随润滑油在润滑系统中不断地循环，会影响润滑油的健康状况。因此，需定期对润滑油的磨粒质量进行检测来判断润滑油的健康状态，同时磨粒的质量在一定程度上也可以反映出机器的健康情况；本实施例提供了一种润滑油磨粒质量检测方法用于对润滑油中的磨粒进行分析，判断润滑油中磨粒的种类；由于不同的磨粒类别对润滑油的影响是不同的，因此根据润滑油中磨粒的类别来对磨粒的质量进行判断，进而根据润滑油中磨粒的质量来反映润滑油的健康状态。

[0074] 首先对机器中的润滑油进行取样，得到待检测润滑油样本；本实施例通过对待检测润滑油样本中的磨粒质量进行检测，以反映机器中润滑油的状态；为了对待检测润滑油样本中的单个磨粒的特征进行分析，本实施例首先对待检测润滑油样本进行稀释，然后使用电子显微镜对稀释后的待检测润滑油样本进行拍摄，得到初始图像，所述初始图像为灰度图像；最后利用高斯除噪的方法对初始图像进行处理，去除初始图像中的噪声，并将去噪后的初始图像记为待检测润滑油样本图像，其中高斯除噪的方法为现有技术，就不再赘述。

[0075] 接下来，对待检测润滑油样本图像中的磨粒区域进行提取，一个磨粒区域对应一个磨粒；具体的：

[0076] 将待检测润滑油样本图像输入到训练好的语义分割网络中，得到仅包含磨粒区域的图像，记为磨粒图像，由于待检测润滑油样本中可能会存在许多磨粒，所以磨粒图像中包含多个磨粒区域，一个磨粒区域为对应磨粒在图像中所占的区域。所述语义分割网络为ResNet神经网络；该语义分割网络的训练集为大量的带有磨粒的润滑油图像（预处理后的图像），将润滑油图像中磨粒区域标记为1，其他区域标记为0，作为对应的标签；然后利用所述训练集和对应的标签对语义分割网络进行训练，得到训练好的语义分割网络，其训练过程中的损失函数为交叉熵损失函数。语义分割网络以及对应的训练过程为现有技术，就不再赘述。

[0077] 步骤S2，获取各磨粒区域对应的长宽比；根据各磨粒区域对应的轮廓像素点，得到各磨粒区域对应的弯曲度，所述磨粒区域对应的轮廓像素点为磨粒区域最外层的像素点；根据各磨粒区域对应的最小拟合圆的面积，得到磨粒区域对应的球状度；根据各磨粒区域对应的灰度共生矩阵，得到各磨粒区域对应的粗糙度；根据各磨粒区域中的像素点，得到各磨粒区域对应的空洞指标。

[0078] 根据步骤S1得到了磨粒图像，磨粒图像中包含了待检测润滑油样本中存在的各种磨粒对应的磨粒区域；本实施例中磨粒的类别包括七类，即正常滑动磨粒、严重滑动磨粒、切削磨粒、剥块磨粒、层状磨粒、球形磨粒以及氧化物磨粒；每种类别的磨粒都有其独特的特征可以将其进行区分，接下来本实施例分别对磨粒图像中的各磨粒区域的特征进行分析，便于后续对磨粒图像中的各磨粒区域进行分类，具体的：

[0079] 第一，获取各磨粒区域对应的长宽比。

[0080] 对于切割产生的磨粒（即切削磨粒），其形状是长向弯曲的；因此本实施例通过计算磨粒区域的长宽比来区分切削磨粒。对于任一磨粒区域：该磨粒区域对应的长宽比为该磨粒区域对应的最小外接矩形的长和宽的比值。

[0081] 至此，能够得到各磨粒区域对应的长宽比。

[0082] 第二，获取各磨粒区域对应的弯曲度。

[0083] 考虑到严重滑动产生的磨粒（即严重滑动磨粒）是细长形的，因此严重滑动磨粒和切削磨粒的长宽比都会比较大，所以根据长宽比无法将严重滑动磨粒和切削磨粒进行区分；但是切削磨粒的形状更散乱，更加弯曲，因此本实施例对各磨粒区域的轮廓进行分析，得到各磨粒区域对应的弯曲度，所述磨粒区域的轮廓为磨粒区域对应在最外层的像素点组成轮廓（即磨粒区域对应的边界像素点组成的轮廓），将轮廓上的像素点记为轮廓像素点。

[0084] 对于任一磨粒区域：

[0085] 对于该磨粒区域对应的轮廓上的任一轮廓像素点：获取该磨粒区域对应的轮廓上与该轮廓像素点的直线距离最远的轮廓像素点，将所述与该轮廓像素点距离最远的轮廓像素点记为该轮廓像素点对应的目标像素点；为了反映该磨粒区域的弯曲度，本实施例利用弗洛伊德算法获取该轮廓像素点与对应的目标像素点之间的内距离，即利用弗洛伊德算法获取该轮廓像素点与对应的目标像素点在磨粒区域内部的最短路径，然后计算所述最短路径的长度，作为该轮廓像素点与对应的目标像素点的内距离。所述弗洛伊德算法和计算内距离的方法为现有技术，不再赘述；所述直线距离为根据两点坐标得到的，例如欧氏距离。

[0086] 至此，能够得到该磨粒区域对应的各轮廓像素点与对应的目标像素点的内距离。考虑到如果该磨粒区域不弯曲，则各轮廓像素点与对应的目标像素点的直线距离和内距离应该是相近的，因此根据该磨粒区域对应的各轮廓像素点与对应的目标像素点的直线距离和内距离，计算该磨粒区域对应的弯曲度，即：

[0087]

[0088] 其中，W为该磨粒区域对应的弯曲度，m为该磨粒区域对应的轮廓像素点的数量，为该磨粒区域对应的第i个轮廓像素点与对应的目标像素点的直线距离， 为该磨粒区域对应的第i个轮廓像素点与对应的目标像素点的内距离。当各轮廓像素点与对应的目标像素点的内距离和直线距离的比值越大时，则对应的弯曲度越大，说明该磨粒越弯曲。

[0089] 至此，根据上述过程能够得到各磨粒区域对应的弯曲度。

[0090] 第三，获取各磨粒区域对应的球状度。

[0091] 对于球形磨粒，其形状为球形，因此对应的磨粒区域应近似圆形；本实施例获取各磨粒区域对应的面积（所述磨粒区域对应的面积为磨粒区域包含的属于该磨粒区域的像素点数量），然后获取各磨粒区域对应的最小拟合圆的面积；球形磨粒对应的磨粒区域的面积与对应的最小拟合圆的面积应该是相近的。对于任一磨粒区域：计算该磨粒区域对应的面积与对应的最小拟合圆的面积的比值，将所述比值作为该磨粒区域对应的球状度；当球状度越大时，说明该磨粒区域越接近圆形，即对应的磨粒越接近球形，越可能是球形磨粒。

[0092] 至此，能够得到各磨粒区域对应的球状度。

[0093] 第四，获取各磨粒区域对应的粗糙度。

[0094] 对于剥块磨粒和氧化物磨粒，剥块磨粒表面较为光滑，而氧化物磨粒表面非常粗糙，因此本实施例通过对各磨粒区域的纹理进行分析，判断磨粒区域的光滑情况；首先获取各磨粒区域对应的灰度共生矩阵，灰度共生矩阵可以分析表面的光滑度和粗糙度，便于区分剥块磨粒和氧化物磨粒；然后根据各磨粒区域对应的灰度共生矩阵，计算得到各磨粒区域对应的灰度共生矩阵对应的熵和逆方差（逆方差又称反差分矩阵），其中熵值可以反映纹理的杂乱程度，即粗糙程度，熵值越大，说明图像中的纹理越杂乱，图像越粗糙，熵值越小说明图像越光滑；而逆方差可以反映纹理的清晰度和规则程度，逆方差越大，说明图像中的纹理越清晰，规律性越强，易于描述，即图像越粗糙，逆方差越小，说明图像中的纹理不易于描述，即图像越光滑。本实施例中获取灰度共生矩阵，根据灰度共生矩阵计算对应的熵和逆方差的过程为现有技术，不再赘述。

[0095] 本实施例根据各磨粒区域对应的灰度共生矩阵对应的熵和逆方差，计算各磨粒区域对应的粗糙度；对于任一磨粒区域，该磨粒区域对应的粗糙度的计算公式如下：

[0096]

[0097] 其中，R为该磨粒区域对应的粗糙度，Ent为该磨粒区域对应的灰度共生矩阵的熵，IDM为该磨粒区域对应的灰度共生矩阵的逆方差，exp( )为以e为底的指数函数。

[0098] 该磨粒区域对应的灰度共生矩阵的熵值越大，逆方差越大，说明该磨粒区域越粗糙，即对应的粗糙度越大；反之，该磨粒区域对应的灰度共生矩阵的熵值越小，逆方差越小，说明该磨粒区域越光滑，即对应的粗糙度越小。粗糙度越大说明磨粒区域为氧化物磨粒的可能性越大，粗糙度越小说明磨粒区域为剥块磨粒的可能性越大。

[0099] 至此，根据上述过程能够得到各磨粒区域对应的粗糙度。

[0100] 第五，获取各磨粒区域对应的空洞指标。

[0101] 对于层状磨粒，层状磨粒对应的磨粒区域内部会存在空洞，即内部会存在不属于该磨粒区域的区域，因此本实施例对各磨粒区域进行检测，判断磨粒区域中是否存在闭合边缘，若存在，则说明对应的磨粒区域内存在空洞，若不存在，则说明对应的磨粒区域内不存在空洞。接下来对各磨粒区域进行分析，得到各磨粒区域对应的空洞指标，具体的：

[0102] 对于任一磨粒区域：

[0103] 首先使用canny算子对该磨粒区域进行边缘检测，得到该磨粒区域中存在的各边缘，并将该磨粒区域中边缘上的像素点记为边缘像素点；考虑到闭合边缘上的每个边缘像素点的八邻域内至少存在两个边缘像素点，因此本实施例对该磨粒区域内的每个边缘分别进行分析，以判断该磨粒区域内的边缘是否为闭合边缘。

[0104] 对于该磨粒区域内的任一边缘：对该边缘上的边缘像素点进行遍历，判断该边缘上的边缘像素点对应的八邻域内存在的边缘像素点数量是否大于等于2，若是，则将对应的边缘像素点标记为0，若否，则将对应的边缘像素点标记为1；计算该边缘上的各边缘像素点的标记值的累加和，若累加和为0，则判定该边缘为闭合边缘，若累加和不为0，则判定该边缘不是闭合边缘。

[0105] 本实施例根据上述方法对该磨粒区域内的所有边缘都进行判断，若该磨粒区域内存在闭合边缘，则说明该磨粒区域内存在空洞，将该磨粒区域对应的空洞指标赋值为1；若该磨粒区域内不存在闭合边缘，则说明该磨粒区域内不存在空洞，将该磨粒区域对应的空洞指标赋值为0。

[0106] 至此，根据上述过程能够得到各磨粒区域对应的空洞指标。

[0107] 步骤S3，根据磨粒类别数量、K‑means算法以及各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，对各磨粒区域进行分类，得到各磨粒类别对应的最优簇。

[0108] 本实施例根据步骤S2得到了各磨粒区域对应的各特征值，即对应的长宽比、弯曲度、球状度、粗糙度和空洞指标（即各磨粒对应的各特征值）；接下来本实施例根据各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，对各磨粒区域进行分类，便于后续对润滑油的质量进行检测；具体的：

[0109] 首先根据各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，将各磨粒区域映射到多维空间当中，使各磨粒区域在多维空间中可以看作一个点，记为磨粒点；本实施例使用K‑means聚类算法对多维空间中的各磨粒点进行分类，其中聚类数量K为磨粒种类的数量，本实施例中磨粒一共有7种，因此K=7。

[0110] 首先在多维空间中获取七个初始质心，然后根据各磨粒点对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，计算各磨粒点与各初始质心之间的特征距离，即：

[0111]

[0112] 其中， 为第j个磨粒点与第q个初始质心之间的特征距离， 为第j个磨粒点对应的长宽比， 为第j个磨粒点对应的弯曲度， 为第j个磨粒点对应的球状度， 为第j个磨粒点对应的粗糙度， 为第j个磨粒点对应的空洞指标， 为第q个初始质心对应的长宽比， 为第q个初始质心对应的弯曲度， 为第q个初始质心对应的球状度， 为第q个初始质心对应的粗糙度， 为第q个初始质心对应的空洞指标。

[0113] 然后本实施例基于各磨粒点与各初始质心之间的特征距离，利用K‑means聚类算法对各磨粒点进行聚类（即对各磨粒区域进行聚类）；本实施例为了保证聚类得到的每一簇都能代表一种磨粒类别，因此在上述进行聚类时，在各磨粒点中加入已经分类好的样本磨粒对应的样本磨粒点，即每种类别选取10个样本磨粒，然后根据对应的特征值确定样本磨粒点（样本磨粒的具体数量可根据实际需要进行设置）；聚类完成后可以根据簇中包含的各样本磨粒点中占比最大磨粒类别，确定对应簇的磨粒类别（例如：对于任意一个簇，若该簇中包含的切削磨粒对应的样本磨粒点最多，则该簇的磨粒类别为切削磨粒）。

[0114] 至此，根据上述过程可以得到各磨粒类别对应的初始簇，一个初始簇中包含多个磨粒点，即包含多个磨粒区域。

[0115] 由于使用K‑means算法进行聚类时，其聚类效果与初始质心的选取有关，因此为了达到较好的聚类效果，本实施例通过对聚类完成后的每个簇进行分析，进而判断此时聚类的聚类效果，若聚类效果没有达到要求，则重新选取聚类中心，并重新进行聚类，直到聚类效果达到要求为止。

[0116] 本实施例在上述过程中已经完成了第一次聚类，得到了各磨粒类别对应的初始簇，接下来，本实施例分别对各磨粒类别对应的初始簇进行分析，具体的：

[0117] 第一，考虑到如果聚类效果好，那么每个初始簇内的磨粒点对应的各特征值应该相似，因此本实施例根据各初始簇内各磨粒点对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，计算此次聚类的簇内差异指标，即：

[0118]

[0119] 其中，p1为簇内差异指标，K为初始簇的数量（即磨粒类别的数量，所以K=7）， 为第n个磨粒类别对应的初始簇内磨粒点的数量， 为第n个磨粒类别对应的初始簇内第c个磨粒点与对应的聚类中心之间的特征距离；所述聚类中心对应的长宽比、弯曲度、球状度、粗糙度和空洞指标为对应初始簇内各磨粒点对应的长宽比的均值、弯曲度的均值、球状度的均值、粗糙度的均值以及空洞指标的均值；所述磨粒点与聚类中心之间的特征距离是根据磨粒点和聚类中心对应的长宽比、弯曲度、球状度、粗糙度和空洞指标得到的，与上述计算特征距离的方法相同，在此就不再阐述。

[0120] p1越大，说明此次聚类得到的各初始簇内各磨粒点之间的差异越大，相似度越低，即说明聚类效果越不好；反之，p1越小，说明此次聚类得到的各初始簇内各磨粒点之间的差异越小，相似度越高，即说明聚类效果越好。

[0121] 第二，如果聚类效果好，那么不同磨粒种类对应的初始簇内的磨粒点应越不相似，即不同磨粒类别对应的初始簇应越独立，因此本实施例获取每个初始簇对应的聚类中心的长宽比、弯曲度、球状度、粗糙度和空洞指标；然后计算任意一个初始簇对应的聚类中心与其他初始簇对应的聚类中心的特征距离，以得到此次聚类的簇间差异指标，即：

[0122]

[0123] 其中，p2为簇间差异指标， 为第n1个初始簇对应的聚类中心与第n2个初始簇对应的聚类中心之间的特征距离；上述公式反映了任意一个初始簇与其他初始簇之间的差异程度， 越大，则p2越大，说明此次聚类得到的各初始簇之间的差异越大，说明聚类效果越好；反之， 越小，则p2越小，说明此次聚类得到的各初始簇之间的差异越小，说明聚类效果越差。

[0124] 第三，考虑到正常滑动磨粒与球形磨粒仅用球状度进行区分很容易区分错误，因此本实施例将聚类结果中类别为正常滑动磨粒对应的初始簇与球形磨粒对应的初始簇提取出来，并对这两个初始簇进行分析，获取聚类评价因子。具体的：

[0125] 球形磨粒的边缘变化是平滑的，而正常滑动磨粒的边缘变化不是很平滑；并且由于球形磨粒较为立体，因此对应图像中心部分的灰度值较高（即对应的磨粒区域内中心区域的像素点的灰度值较大）；因此若球形磨粒对应的初始簇对应的各磨粒区域的轮廓越平滑、中心区域平均灰度值越大，说明聚类效果越好，除此之外，若正常滑动磨粒对应的初始簇对应的各磨粒区域的轮廓越不平滑、中心区域平均灰度值无明显变化，也说明聚类效果越好。所述初始簇对应的各磨粒区域即为初始簇中的各磨粒点对应的磨粒区域。

[0126] 本实施例首先对正常滑动磨粒和球形磨粒对应的初始簇对应的磨粒区域的轮廓的平滑度进行分析，得到正常滑动磨粒和球形磨粒对应的初始簇对应的磨粒区域的轮廓平滑度。对于正常滑动磨粒和球形磨粒中任一磨粒类别对应的初始簇：

[0127] 对于该磨粒类别对应的初始簇对应的任一磨粒区域：首先获取该磨粒区域对应的最小拟合圆，以及最小拟合圆的圆心，将所述圆心记为该磨粒区域对应的轮廓中心；然后计算该磨粒区域对应的各轮廓像素点与该磨粒区域对应的轮廓中心的欧氏距离。

[0128] 本实施例根据该磨粒类别对应的初始簇对应的各磨粒区域的各轮廓像素点与对应磨粒区域的轮廓中心的欧氏距离，计算该磨粒类别对应的初始簇的轮廓评价指标，即：

[0129]

[0130]

[0131] 其中，F为该磨粒类别对应的初始簇的轮廓评价指标，U为该磨粒类别对应的初始簇对应的各磨粒区域的数量（即对应的初始簇中磨粒点的数量）， 为该磨粒类别对应的初始簇对应的第u个磨粒区域的轮廓像素点的数量， 为该磨粒类别对应的初始簇对应的第u个磨粒区域的平滑指标， 为所述第u个磨粒区域对应的第i个轮廓像素点与所述第u个磨粒区域对应的轮廓中心的欧氏距离， 为所述第u个磨粒区域对应的第i+1个轮廓像素点与所述第u个磨粒区域对应的轮廓中心的欧氏距离， 为取绝对值函数。

[0132] 本实施例通过分析每个磨粒区域的轮廓像素点与对应的轮廓中心的欧氏距离的差异，来反映该磨粒类别对应的初始簇对应的各磨粒区域的轮廓是否平滑；若 越大，说明相邻轮廓像素点与轮廓中心的欧氏距离的差异越大，说明轮廓越不平滑，则 越大，反之，若 越小，说明相邻轮廓像素点与轮廓中心的欧氏距离的差异越小，说明轮廓越平滑，则越小。至此，根据上述过程能够得到正常滑动磨粒对应的初始簇的轮廓评价指标和球形磨粒对应的初始簇的轮廓评价指标。

[0133] 接着，对于正常滑动磨粒和球形磨粒中任一磨粒类别对应的初始簇：

[0134] 本实施例首先获取该磨粒类别对应的初始簇对应的各磨粒区域的中心区域，对于任一磨粒区域：以该磨粒区域的最小拟合圆的圆心为目标圆心，以最小拟合圆的半径的一半为目标半径作圆，得到一个目标圆形区域，将该磨粒区域中的目标圆形区域作为该磨粒区域的中心区域；计算该磨粒区域的中心区域内各像素点的灰度值的平均值，作为该磨粒区域的中心灰度均值。本实施例中各磨粒区域是步骤S1得到的磨粒图像中的区域，因此像素点的灰度值可直接从磨粒图像中得到。

[0135] 计算该磨粒类别对应的初始簇对应的各磨粒区域的中心灰度均值的平均值，作为该磨粒类别对应的初始簇的中心区域评价指标。

[0136] 至此本实施例能够得到正常滑动磨粒对应的初始簇的中心区域评价指标和球形磨粒对应的初始簇的中心区域评价指标。若球形磨粒对应的初始簇的中心区域评价指标越大，说明该初始簇对应的各磨粒区域的中心区域的平均灰度值越大，说明该初始簇的聚类效果越好；若正常滑动磨粒对应的初始簇的中心区域评价指标越大，说明该初始簇的聚类效果越不好。

[0137] 本实施例根据上述过程得到了正常滑动磨粒对应的初始簇的轮廓评价指标和中心区域评价指标、球形磨粒对应的初始簇的轮廓评价指标和中心区域评价指标，二者均可以反映此次聚类的效果，因此根据得到的轮廓评价指标和中心区域评价指标，计算此次聚类的聚类评价因子，即：

[0138]

[0139] 其中，p3为聚类评价因子，为球形磨粒对应的初始簇的轮廓评价指标， 为球形磨粒对应的初始簇的中心区域评价指标， 为正常滑动磨粒对应的初始簇的轮廓评价指标， 为正常滑动磨粒对应的初始簇的中心区域评价指标。

[0140] 当 越小， 越大时，说明球形磨粒对应的初始簇对应的各磨粒区域的轮廓越光滑，中心区域的灰度值越大，说明该初始簇的聚类效果越好，此时p3越小；当 越大， 越小时，说明正常滑动磨粒对应的初始簇对应的各磨粒区域的轮廓越不光滑，中心区域的灰度值越小，说明该初始簇的聚类效果越好，此时p3越小。

[0141] 第四，本实施例根据上述过程得到了此次聚类的簇内差异指标、簇间差异指标以及聚类评价因子，这三个参数在不同层面反映了此次聚类的聚类效果；接下来本实施例将此次聚类的簇内差异指标、簇间差异指标以及聚类评价因子进行融合，得到聚类效果指标；本实施例利用聚类效果指标来综合反映聚类的效果。聚类效果指标的计算公式如下：

[0142]

[0143] 其中，P为聚类效果指标。

[0144] 当p1和p3越小，p2越大时，说明聚类效果越好，越能够将不同磨粒类别的磨粒区域区分开，此时P越大；反之，当p1和p3越大，p2越小时，说明聚类效果越差，越不能将不同磨粒类别的磨粒区域区分开，此时P越小。

[0145] 本实施例中此次聚类记为第一次聚类，本实施例判断第一次聚类的聚类效果指标是否大于效果阈值，若是，则判定第一次聚类的聚类效果已达到要求，则将第一次聚类得到各磨粒类别对应的初始簇记为最优簇；若否，则从多维空间中重新选取七个初始质心，并对多维空间中的各磨粒点进行重新聚类（即进行第二次聚类），再次得到各磨粒类别对应的第二初始簇，并计算第二次聚类的聚类效果指标，判断第二次聚类的聚类效果指标是否大于效果阈值，若是，则判定第二次聚类的聚类效果已达到要求，则将第二次聚类得到各磨粒类别对应的第二初始簇记为最优簇；若否，则继续重新选取七个初始质心，对多维空间中各磨粒点进行重新聚类，直到聚类效果指标第一次大于效果阈值为止，将最后一次聚类得到的各磨粒类别对应的簇记为最优簇。本实施例中所述效果阈值的大小由实施者根据经验进行设置。

[0146] 至此，本实施例得到了各磨粒类别对应的最优簇，一个最优簇中包含了属于同一磨粒类别的各磨粒点，即属于同一磨粒类别的各磨粒区域，实现了磨粒图像中各磨粒区域的分类。

[0147] 步骤S4，根据各磨粒类别对应的最优簇对应的各磨粒区域的面积，得到润滑油磨粒质量指标；若润滑油磨粒质量指标大于质量阈值，则判定待检测润滑油样本中磨粒的质量良好。

[0148] 本实施例根据步骤S3对磨粒图像中各磨粒区域进行了分类，得到了不同磨粒类别对应的最优簇。考虑到磨粒大小的不同对润滑油质量的影响也不同，此外，不同磨粒类别对润滑油的影响也是不同的，例如正常滑动磨粒对润滑油的影响较小，而其他磨粒类别对润滑油影响较大。

[0149] 本实施例根据不同磨粒类别对应的最优簇对应的各磨粒区域的面积以及不同磨粒类别对润滑油质量的影响程度，得到润滑油磨粒质量指标，所述润滑油磨粒质量指标可以反映待检测润滑油样本中磨粒的情况，进而可以反映待检测润滑油样本对应的润滑油的质量；计算润滑油磨粒质量指标的公式如下：

[0150]

[0151] 其中，V为润滑油磨粒质量指标， 为第 个磨粒类别对润滑油质量的影响程度，为第 个磨粒类别对应的最优簇对应的所有磨粒区域的面积之和（本实施例中磨粒区域的面积已经在步骤S2中获取过了，就不再赘述）。

[0152] 当V越大时，说明待检测润滑油样本中磨粒的质量越好，说明对应的机器中的润滑油的质量越好，也可以反映出机器的磨损较小；反之，当V越小时，说明待检测润滑油样本中磨粒的质量越差，说明对应的机器中润滑油的健康状态越差，也可以反映出机器的磨损较为严重。因此本实施例判断润滑油磨粒质量指标是否大于质量阈值，若大于，则判定待检测润滑油样本中磨粒的质量良好，即对应的机器中的润滑油的健康状态良好，若小于等于，则判定待检测润滑油样本中的磨粒质量较差，即对应的机器中的润滑油的健康状态差。

[0153] 本实施例中不同磨粒类别对润滑油质量的影响程度和所述质量阈值由实施者根据经验进行设置。

[0154] 本实施例提供了一种润滑油磨粒质量检测方法，用于对润滑油中的磨粒进行分析；本实施例依据不同磨粒类别的特点，获取磨粒图像中各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度以及空洞指标；本实施例根据各磨粒区域对应的特征能够更好的将不同磨粒类别的磨粒区域进行区分，便于更好的判断润滑油的质量。本实施例基于磨粒类别数量、K‑means算法以及各磨粒区域对应的长宽比、弯曲度、球状度、粗糙度和空洞指标，对各磨粒区域进行分类，得到各磨粒类别对应的最优簇；本实施例在获取最优簇的过程中，不断变换初始质心进行重新聚类，并分析每一次聚类的聚类效果，以得到聚类效果好的最优簇，避免了由于初始质心选取的随机性而造成的聚类效果不理想的情况。本实施例结合不同磨粒类别对润滑油的不同影响，得到润滑油磨粒质量指标；最后根据润滑油磨粒质量指标判断润滑油中磨粒的质量的好坏。本实施例基于不同磨粒类别的特征，并利用变换初始质心的K‑means算法，更加准确对图像中的磨粒进行分类，使本发明能够更加准确的对润滑油中磨粒的质量进行判断。

[0155] 需要说明的是：以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。