一种光学玻璃表面缺陷检测方法与流程

1.本发明涉及图像数据处理技术领域，具体涉及一种光学玻璃表面缺陷检测方法。


背景技术：

2.光学玻璃是光电技术产业的基础和重要组成部分。随着光学与电子信息科学、新材料科学的不断融合，作为光电子基础材料的光学玻璃在光传输、光储存和光电显示三大领域的应用更是突飞猛进，成为社会信息化尤其是光电信息技术发展的基础条件之一。由于光学玻璃生产工艺和环境的影响，会生产出有缺陷的光学玻璃，这严重影响了玻璃的结构性能，例如光学性质和高度光学均匀性。由此，光学玻璃缺陷检测在光学玻璃生产中是非常有必要的一步。
3.由于气泡缺陷多变的形态结构，导致气泡缺陷的边缘模糊较难分辨，采用传统的边缘检测方法并不能保证可以提取完全出气泡缺陷的边缘，导致无法识别出气泡缺陷的边缘，也就无法准确检测到玻璃图像中的是否存在气泡缺陷。另外，由于气泡表面的光照强度具有不均衡性，导致同一个气泡具有不同的亮度区域，而不同的亮度区域与正常玻璃区域的灰度值又比较接近，采用传统的阈值分割方法很难保证将气泡区域准确地识别出来。因此，现有方法无法准确地检测光学玻璃表面是否存在气泡缺陷。


技术实现要素：

4.为了解决上述现有光学玻璃表面缺陷检测不准确的技术问题，本发明的目的在于提供一种光学玻璃表面缺陷检测方法。
5.本发明提供了一种光学玻璃表面缺陷检测方法，包括以下步骤：获取待检测玻璃的表面图像，进而获取图像预处理后的待检测玻璃的灰度图像，将所述待检测玻璃的灰度图像等分为n个图像块，从而得到n个图像块对应的纹理特征向量；根据n个图像块对应的纹理特征向量，确定任意两个图像块之间的纹理特征距离，进而确定纹理特征距离的最优分割阈值；根据任意两个图像块之间的纹理特征距离以及纹理特征距离的最优分割阈值，确定各个正常图像块和各个第一疑似气泡缺陷图像块；根据各个第一疑似气泡缺陷图像块对应的多个纹理特征距离和纹理特征距离的最优分割阈值，确定各个第一疑似气泡缺陷图像块的第一缺陷程度指标；根据各个正常图像块对应的多个纹理特征距离、纹理特征向量，各个第一疑似气泡缺陷图像块对应的多个纹理特征距离、纹理特征向量以及纹理特征距离的最优分割阈值，确定第一缺陷程度指标阈值；根据各个第一疑似气泡缺陷图像块的第一缺陷程度指标和第一缺陷程度指标阈值，判断各个第一疑似气泡缺陷图像块中是否存在气泡缺陷。
6.进一步的，还包括：
获取各个第一疑似气泡缺陷图像块中存在气泡缺陷的第一疑似气泡缺陷图像块，并将存在气泡缺陷的各个第一疑似气泡缺陷图像块作为第二疑似气泡缺陷图像块；根据各个第二疑似气泡缺陷图像块内每个像素点的灰度值，确定每个像素点对应的梯度幅值，进而确定梯度幅值的最优分割阈值；根据每个像素点对应的梯度幅值以及梯度幅值的最优分割阈值，确定各目标像素点，所述目标像素点为梯度幅值大于最优分割阈值的像素点，进而对各目标像素点进行聚类操作，从而得到满足预设聚类条件的各个聚类簇，所述预设聚类条件为聚类簇内目标像素点的数目大于预设目标像素点数目；根据各个聚类簇内各目标像素点的坐标位置，确定各个聚类簇的最小外接矩形，进而根据各个聚类簇的最小外接矩形内每个目标像素点的灰度值，确定各个聚类簇的最小外接矩形的灰度丰富指标；获取灰度丰富指标阈值，根据各个聚类簇的最小外接矩形的灰度丰富指标和灰度丰富指标阈值，判断各个聚类簇的最小外接矩形中是否存在气泡缺陷。
7.进一步的，确定各个聚类簇的最小外接矩形的灰度丰富指标的步骤包括：根据各个聚类簇的最小外接矩形内每个目标像素点的灰度值，确定各个聚类簇的最小外接矩形内灰度类别的数目，进而确定各个灰度类别对应的像素灰度均值；根据各个聚类簇的最小外接矩形内每个目标像素点的灰度值，确定各个灰度类别的目标像素点在对应最小外接矩形中出现的频率；根据各个聚类簇的最小外接矩形内灰度类别的数目、各个灰度类别对应的像素灰度均值、各个灰度类别的像素点在对应最小外接矩形中出现的频率以及各个灰度类别的像素灰度，确定各个聚类簇的最小外接矩形的灰度丰富指标。
8.进一步的，确定各个聚类簇的最小外接矩形的灰度丰富指标的计算公式为：其中，为各个聚类簇的最小外接矩形的灰度丰富指标，b为各个聚类簇的最小外接矩形内灰度类别的数目，为第i个灰度类别的像素点在对应最小外接矩形中出现的频率，为第i个灰度类别的像素灰度，为各个灰度类别对应的像素灰度均值。
9.进一步的，确定每个像素点对应的梯度幅值的步骤包括：根据各个第二疑似气泡缺陷图像块内每个像素点的灰度值，确定每个像素点对应的水平方向的梯度和竖直方向的梯度；根据每个像素点对应的水平方向的梯度和竖直方向的梯度，确定每个像素点对应的梯度幅值。
10.进一步的，所述纹理特征向量包含4个纹理特征，分别为图像块对应的灰度共生矩阵的能量、熵值、对比度以及逆差矩，确定任意两个图像块之间的纹理特征距离的计算公式为：
其中，为第w个图像块与第v个图像块之间的纹理特征距离，为第w个图像块对应纹理特征向量内的第i个纹理特征，为第v个图像块对应纹理特征向量内的第i个纹理特征。
11.进一步的，确定各个正常图像块和各个第一疑似气泡缺陷图像块的步骤包括：统计每个图像块对应的多个纹理特征距离中大于最优分割阈值的纹理特征距离的数目，并统计每个图像块对应的多个纹理特征距离中不大于最优分割阈值的纹理特征距离的数目；若某个图像块对应的大于最优分割阈值的纹理特征距离的数目大于该图像块对应的不大于最优分割阈值的纹理特征距离的数目，则判定该图像块为第一疑似气泡缺陷图像块，否则，判定该图像块为正常图像块。
12.进一步的，确定各个第一疑似气泡缺陷图像块的第一缺陷程度指标的步骤包括：从每个第一疑似气泡缺陷图像块对应的多个纹理特征距离中，筛选出大于最优分割阈值的各纹理特征距离；根据每个第一疑似气泡缺陷图像块对应的大于纹理特征距离的最优分割阈值的各纹理特征距离，确定每个第一疑似气泡缺陷图像块对应的预设数目个较小的大于纹理特征距离的最优分割阈值的纹理特征距离；计算每个第一疑似气泡缺陷图像块对应的预设数目个较小的大于纹理特征距离的最优分割阈值的纹理特征距离的均值，将该均值作为对应第一疑似气泡缺陷图像块的第一缺陷程度指标。
13.进一步的，确定第一缺陷程度指标阈值的计算公式为：其中，t为第一缺陷程度指标阈值，为纹理特征距离的最优分割阈值，和为调节参数，为第一疑似气泡缺陷图像块的数目，为正常图像块的数目，为任意两个正常图像块对应的纹理特征距离的方差，为任意两个第一疑似气泡缺陷图像块对应的纹理特征距离的方差，为各个正常图像块对应的纹理特征向量均值与各个第一疑似气泡缺陷图像块对应的纹理特征向量均值之间的纹理特征距离。
14.进一步的，判断各个第一疑似气泡缺陷图像块中是否存在气泡缺陷的步骤包括：若任意一个第一疑似气泡缺陷图像块的第一缺陷程度指标大于第一缺陷程度指标阈值，则判定该第一疑似气泡缺陷图像块中存在气泡缺陷，否则，判定该第一疑似气泡缺陷图像块中不存在气泡缺陷。
15.本发明具有如下有益效果：本发明对待检测玻璃的表面图像进行图像预处理操作，图像预处理操作有利于后续得到更加准确的图像纹理特征信息，将待检测玻璃的灰度图像等分为n个图像块，是为了克服由于玻璃气泡自身的特点所造成的无法准确提取气泡缺陷边缘的问题，通过分析每个图像块的纹理特征信息，有助于提高了光学玻璃缺陷检测的准确性；通过计算任意两个图
像块之间的纹理特征距离，确定纹理特征距离的最优分割阈值，以最优分割阈值为基础，确定各个正常图像块和各个第一疑似气泡缺陷图像块，对n个图像块进行分析，将图像块划分为两个类别，分别为正常图像块和第一疑似气泡缺陷图像块，气泡缺陷很有可能存在于第一疑似气泡缺陷图像块内，其缩小气泡缺陷检测的范围，有助提高了光学玻璃缺陷检测的效率；通过计算各个第一疑似气泡缺陷图像块的第一缺陷程度指标和第一缺陷程度指标阈值，判断各个第一疑似气泡缺陷图像块中是否存在气泡缺陷，实现玻璃气泡缺陷检测，本发明利用各个正常图像块和各个第一疑似气泡缺陷图像块相关数据得到第一缺陷程度指标阈值，基于自适应方式确定第一缺陷程度指标阈值，有效提高了判断结果的准确性，进而提高了光学玻璃表面缺陷检测的准确性。
附图说明
16.为了更清楚地说明本发明实施例或现有技术中的技术方案和优点，下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它附图。
17.图1为本发明一种光学玻璃表面缺陷检测方法的流程图；图2为本发明实施例中判断各个聚类簇的最小外接矩形中是否存在气泡缺陷的流程图。
具体实施方式
18.为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效，以下结合附图及较佳实施例，对依据本发明提出的技术方案的具体实施方式、结构、特征及其功效，详细说明如下。在下述说明中，不同的“一个实施例”或“另一个实施例”指的不一定是同一个实施例。此外，一个或多个实施例中的特定特征、结构或特点可由任何合适形式组合。
19.除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。
20.本实施例提供了一种光学玻璃表面缺陷检测方法，如图1所示，该方法包括以下步骤：（1）获取待检测玻璃的表面图像，进而获取图像预处理后的待检测玻璃的灰度图像，将所述待检测玻璃的灰度图像等分为n个图像块，从而得到n个图像块对应的纹理特征向量，其步骤包括：（1-1）获取待检测玻璃的表面图像，进而获取图像预处理后的待检测玻璃的灰度图像。
21.本实施例检测玻璃板表面的气泡缺陷，首先利用相机拍摄待检测的玻璃板，得到待检测玻璃的表面图像，待检测玻璃的表面图像为的可见光图像，也就是rgb彩色图像（red blue green）。为了便于后续确定待检测玻璃的表面图像的图像特征信息，对待检测玻璃的表面图像进行图像预处理操作。具体为：将待检测玻璃的表面图像转换成灰度图像，得到待检测玻璃的灰度图像，然后采用高斯滤波对待检测玻璃的灰度图像进行噪声过滤处理，得到过滤处理后的待检测玻璃的灰度图像，之后采用直方图均衡化技术对过滤处
理后的待检测玻璃的灰度图像进行图像增强处理，最终得到图像预处理操作后的待检测玻璃的灰度图像。图像灰度化过程、高斯滤波过程以及直方图均衡化过程均为现有技术，不在本发明保护范围内，此处不再进行详细阐述。
22.（1-2）将所述待检测玻璃的灰度图像等分为n个图像块，从而得到n个图像块对应的纹理特征向量。
23.需要说明的是，由于玻璃的气泡缺陷的形态结构是多变的，多变的形态结构导致很难确定气泡缺陷边缘，而采用传统的边缘检测方法不能保证能够将气泡缺陷边缘准确地提取出来，不能识别出气泡缺陷边缘也就无法检测到玻璃图像的气泡缺陷。考虑到待检测玻璃的表面图像中出现气泡缺陷区域的纹理特征偏离正常区域的纹理特征，由此本实施例将通过分析不同图像块的纹理特征来判断待检测玻璃中是否存在气泡缺陷。
24.为了便于后续准确地确定气泡缺陷的具体位置，将待检测玻璃的灰度图像等分为n个图像块，可以将n设置为，每个图像块的尺寸是相等的。然后，采用灰度共生矩阵提取n个图像块的纹理特征，计算出每个图像块对应的灰度共生矩阵的能量asm、熵值ent、对比度con和逆差矩idm，将该4个数值作为图像块的纹理特征，并将该4个数值构成纹理特征向量，记为，，此时，本实施例得到了n个图像块对应的纹理特征向量。
25.（2）根据n个图像块对应的纹理特征向量，确定任意两个图像块之间的纹理特征距离，进而确定纹理特征距离的最优分割阈值，其步骤包括：（2-1）根据n个图像块对应的纹理特征向量，确定任意两个图像块之间的纹理特征距离。
26.在本实施例中，由于包含气泡缺陷的图像块与不包含缺陷的图像块具有不同的纹理特征，所以本实施例对不同的图像块对应的纹理特征向量进行异常检测，可以实现气泡缺陷检测。纹理特征向量包含4个纹理特征，分别为图像块对应的灰度共生矩阵的能量、熵值、对比度以及逆差矩。对于每个图像块而言，计算每个图像块对应的纹理特征向量与其余的n-1个图像块对应的纹理特征向量之间的纹理特征距离，也就是计算任意两个图像块之间的纹理特征距离，得到每个图像块对应的多个纹理特征距离，其计算公式为：其中，为第w个图像块与第v个图像块之间的纹理特征距离，为第w个图像块对应纹理特征向量内的第i个纹理特征，为第v个图像块对应纹理特征向量内的第i个纹理特征。
27.需要说明的是，任意两个图像块之间的纹理特征距离较大时，该两个图像块中一个是气泡缺陷图像块，另一个图像块是正常图像块；任意两个图像块之间的纹理特征距离较小时，该两个图像块有可能是同类型的图像块，也就是该两个图像块可能都为气泡缺陷图像块，也可能都为正常图像块。
28.（2-2）根据任意两个图像块之间的纹理特征距离，确定纹理特征距离的最优分割
阈值。
29.若待检测玻璃中存在气泡缺陷，可能不止出现一个气泡缺陷，即多个图像块内存在气泡缺陷。基于上述分析，本实施例采用otsu大津法（最大类间方差法）对每个图像块对应的多个纹理特征距离进行分割，得到纹理特征距离的最优分割阈值。otsu大津法确定最优分割阈值的过程为现有技术，不在本发明保护范围内，此处不再进行详细阐述。
30.（3）根据任意两个图像块之间的纹理特征距离以及纹理特征距离的最优分割阈值，确定各个正常图像块和各个第一疑似气泡缺陷图像块，进而根据各个第一疑似气泡缺陷图像块对应的多个纹理特征距离和最优分割阈值，确定各个第一疑似气泡缺陷图像块的第一缺陷程度指标，其步骤包括：（3-1）根据任意两个图像块之间的纹理特征距离以及纹理特征距离的最优分割阈值，确定各个正常图像块和各个第一疑似气泡缺陷图像块，其步骤包括：（3-1-1）统计每个图像块对应的多个纹理特征距离中大于最优分割阈值的纹理特征距离的数目，并统计每个图像块对应的多个纹理特征距离中不大于最优分割阈值的纹理特征距离的数目。
31.需要说明的是，为了避免待检测玻璃的表面图像中可能不止出现一个气泡缺陷的干扰，同时为了排除选择的图像块为正常图像块的干扰，本实施例通过统计每个图像块对应的多个纹理特征距离与最优分割阈值的关系，来区分n个图像块中的正常图像块和第一疑似气泡缺陷图像块。另外，将每个图像块对应的多个纹理特征距离中大于最优分割阈值的纹理特征距离的数目记为，并将每个图像块对应的多个纹理特征距离中不大于最优分割阈值的纹理特征距离的数目记为。
32.（3-1-2）若某个图像块对应的大于最优分割阈值的纹理特征距离的数目大于该图像块对应的不大于最优分割阈值的纹理特征距离的数目，则判定该图像块为第一疑似气泡缺陷图像块，否则，判定该图像块为正常图像块。
33.在本实施例中，若某个图像块对应的大于最优分割阈值的纹理特征距离的数目大于该图像块对应的不大于最优分割阈值的纹理特征距离的数目，判定该图像块为第一疑似气泡缺陷图像块，否则，判定该图像块为正常图像块。基于此，将n个图像块划分为个正常图像块和个第一疑似气泡缺陷图像块。
34.需要说明的是，一般情况下，玻璃图像中出现气泡缺陷的图像块数目要小于没有气泡缺陷正常的图像块数目。由此可知，当某个图像块为疑似气泡缺陷图像块时，计算该图像块与正常图像块之间的纹理特征距离的次数较多，但该图像块对应的纹理特征向量与正常图像对应的纹理特征向量的相似性较低，纹理特征距离较大，那么该图像块对应的大于最优分割阈值的纹理特征距离的数目就会较多；当某个图像块为正常图像块时，计算该图像块与正常图像块之间的纹理特征距离的次数较多，但该图像块对应的纹理特征向量与正常图像对应的纹理特征向量的相似性较高，纹理特征距离较小，那么该图像块对应的不大于最优分割阈值的纹理特征距离的数目较多。
35.（3-2）根据各个第一疑似气泡缺陷图像块对应的多个纹理特征距离和纹理特征距离的最优分割阈值，确定各个第一疑似气泡缺陷图像块的第一缺陷程度指标，其步骤包括：
（3-2-1）从每个第一疑似气泡缺陷图像块对应的多个纹理特征距离中，筛选出大于最优分割阈值的各纹理特征距离。
36.使每个第一疑似气泡缺陷图像块对应的各个纹理特征距离与纹理特征距离的最优分割阈值作比较，筛选出大于纹理特征距离的最优分割阈值的各纹理特征距离，从而得到每个第一疑似气泡缺陷图像块对应的大于纹理特征距离的最优分割阈值的各纹理特征距离。
37.（3-2-2）根据每个第一疑似气泡缺陷图像块对应的大于纹理特征距离的最优分割阈值的各纹理特征距离，确定每个第一疑似气泡缺陷图像块对应的预设数目个较小的大于纹理特征距离的最优分割阈值的纹理特征距离。
38.在本实施例中，对步骤（3-2-1）得到的每个第一疑似气泡缺陷图像块对应的大于纹理特征距离的最优分割阈值的各纹理特征距离按照从小到大的顺序进行排序，得到排序后的大于纹理特征距离的最优分割阈值的各纹理特征距离，选取预设数目个较小的大于纹理特征距离的最优分割阈值的纹理特征距离，本实施例中将预设数目设置。
39.（3-2-3）计算每个第一疑似气泡缺陷图像块对应的预设数目个较小的大于纹理特征距离的最优分割阈值的纹理特征距离的均值，并将该均值作为对应第一疑似气泡缺陷图像块的第一缺陷程度指标。
40.需要说明的是，个较小的大于纹理特征距离的最优分割阈值的纹理特征距离的均值能够表征对应第一疑似气泡缺陷图像块的纹理特征，选取预设数目较小的大于最优分割阈值的纹理特征距离是为了提高气泡缺检测的准确性，当第一疑似气泡缺陷图像块对应的较小的纹理特征距离都大于后续确定的第一缺陷程度指标阈值，那么该第一疑似气泡缺陷图像块内存在气泡缺陷的可能性就会较大。计算均值的过程为现有技术，不在本发明保护范围内，此处不再进行详细阐述。
41.（4）根据各个正常图像块对应的多个纹理特征距离、纹理特征向量，各个第一疑似气泡缺陷图像块对应的多个纹理特征距离、纹理特征向量以及纹理特征距离的最优分割阈值，确定第一缺陷程度指标阈值，进而根据各个第一疑似气泡缺陷图像块的第一缺陷程度指标和第一缺陷程度指标阈值，判断各个第一疑似气泡缺陷图像块中是否存在气泡缺陷。
42.（4-1）根据各个正常图像块对应的多个纹理特征距离、纹理特征向量，各个第一疑似气泡缺陷图像块对应的多个纹理特征距离、纹理特征向量以及纹理特征距离的最优分割阈值，确定第一缺陷程度指标阈值。
43.通常情况下，经过专业人员的经验来确定最终阈值或经过多次实验采用试凑法选择效果最好的阈值作为最终阈值，该方法所确定的阈值的精准性较差，很容易出现所确定的阈值不合理的情况。为了避免该方法所带来的弊端，本实施例通过自适应的方法确定第一缺陷程度指标阈值，其计算公式为：其中，t为第一缺陷程度指标阈值，为纹理特征距离的最优分割阈值，和为调节参数，根据经验值，，为第一疑似气泡缺陷图像块的数目，为正
常图像块的数目，为任意两个正常图像块对应的纹理特征距离的方差，为任意两个第一疑似气泡缺陷图像块对应的纹理特征距离的方差，为各个正常图像块对应的纹理特征向量均值与各个第一疑似气泡缺陷图像块对应的纹理特征向量均值之间的纹理特征距离，所述纹理特征向量均值是指各个图像块对应的4个纹理特征相加，再对各个图像块对应的相加后的4个纹理特征求均值，将各个图像块对应的每个纹理特征的均值作为各个图像块对应的纹理特征向量均值。
44.，为修正参量，若正常图像块的数目与第一疑似气泡缺陷图像块的数目的比值越小，说明将最优分割阈值作为第一缺陷程度指标阈值越合理，减少修正参量，反之，则增大修正参量，增大修正参量是为了降低图像块被认为是包含气泡缺陷的图像块的可能性，减少正常图像块和疑似气泡缺陷图像块被误分类的概率；若任意两个正常图像块和任意两个第一疑似气泡缺陷图像块对应的纹理特征距离的方差之和较小，说明各个正常图像块和各个第一疑似气泡缺陷图像块对应的纹理特征距离的稳定性都比较高，此时被误分类的概率会比较小，也就是将最优分割阈值作第一缺陷程度指标阈值进行分类时被误分类的可能性比较小；若各个正常图像块对应的纹理特征向量均值和各个第一疑似气泡缺陷图像块对应的纹理特征向量均值之间的纹理特征距离越小，说明正常图像块与第一疑似气泡缺陷图像块的差距越小，越难区分正常图像块与第一疑似气泡缺陷图像块，此时需要增大修正参量，增大第一缺陷程度指标阈值t，减少被误分类的概率。
45.（4-2）根据各个第一疑似气泡缺陷图像块的第一缺陷程度指标和第一缺陷程度指标阈值，判断各个第一疑似气泡缺陷图像块中是否存在气泡缺陷。
46.若任意一个第一疑似气泡缺陷图像块的第一缺陷程度指标大于第一缺陷程度指标阈值，则判定该第一疑似气泡缺陷图像块中存在气泡缺陷，否则，判定该第一疑似气泡缺陷图像块中不存在气泡缺陷。也就是当各个第一疑似气泡缺陷图像块的第一缺陷程度指标均不大于第一缺陷程度指标阈值时，说明待检测玻璃的表面图像中不存在气泡缺陷，否则，说明待检测玻璃中存在气泡缺陷。
47.当待检测玻璃中存在气泡缺陷，受到光照强度不均衡性的影响，同一个气泡表面会产生不同的亮度区域，玻璃的亮度区域与正常区域的灰度值比较接近，采用传统的边缘检测方法和阈值分割方法就不能保证将气泡缺陷全部检测出来。为了提高光学玻璃表面缺陷检测的准确性，对存在气泡缺陷的各第一疑似气泡缺陷图像块再次进行气泡缺陷检测，判断各个聚类簇的最小外接矩形中是否存在气泡缺陷的流程图如图2所示，其步骤包括：（4-2-1）获取各个第一疑似气泡缺陷图像块中存在气泡缺陷的第一疑似气泡缺陷图像块，并将存在气泡缺陷的各个第一疑似气泡缺陷图像块作为第二疑似气泡缺陷图像块。
48.在本实施例中，根据各个第一疑似气泡缺陷图像块的第一缺陷程度指标与第一缺陷程度指标阈值的比较结果，得到大于第一缺陷程度指标阈值的各第一疑似气泡缺陷图像
块，将大于第一缺陷程度指标阈值的第一疑似气泡缺陷图像块判定为存在气泡缺陷，从而得到各个第一疑似气泡缺陷图像块中存在气泡缺陷的第一疑似气泡缺陷图像块，并将存在气泡缺陷的各个第一疑似气泡缺陷图像块作为第二疑似气泡缺陷图像块。
49.（4-2-2）根据各个第二疑似气泡缺陷图像块内每个像素点的灰度值，确定每个像素点对应的梯度幅值，进而确定梯度幅值的最优分割阈值，其步骤包括：（4-2-2-1）根据各个第二疑似气泡缺陷图像块内每个像素点的灰度值，确定每个像素点对应的水平方向的梯度和竖直方向的梯度。
50.在本实施例中，通过各个第二疑似气泡缺陷图像块内每个像素点的灰度值，采用sobel算子，计算每个像素点对应的水平方向的梯度和竖直方向的梯度，并将每个像素点的水平方向的梯度记为，将竖直方向的梯度记为。sobel算子计算像素点梯度的过程为现有技术，不在本发明保护范围内，此处不再进行详细阐述。
51.（4-2-2-2）根据每个像素点对应的水平方向的梯度和竖直方向的梯度，确定每个像素点对应的梯度幅值。
52.通过步骤（4-2-2-1）得到的每个像素点对应的水平方向的梯度和竖直方向的梯度，计算每个像素点对应的梯度幅值，其计算公式为：，g为每个像素点对应的梯度幅值。
53.（4-2-2-3）根据每个像素点对应的梯度幅值，确定梯度幅值的最优分割阈值。
54.在本实施例中，为了便于后续确定目标像素点，基于每个像素点对应的梯度幅值，采用otsu大津法对每个像素点对应的梯度幅值进行分割处理，从而得到梯度幅值的最优分割阈值。otsu大津法对每个像素点对应的梯度幅值进行分割处理的过程为现有技术，不在本发明保护范围内，此处不再进行详细阐述。
55.（4-2-3）根据每个像素点对应的梯度幅值以及梯度幅值的最优分割阈值，确定各目标像素点，所述目标像素点为梯度幅值大于最优分割阈值的像素点，进而对各目标像素点进行聚类操作，从而得到满足预设聚类条件的各个聚类簇，所述预设聚类条件为聚类簇内目标像素点的数目大于预设目标像素点数目。
56.在本实施例中，对每个像素点对应的梯度幅值与梯度幅值的最优分割阈值作比较，得到梯度幅值大于最优分割阈值的各个像素点，将梯度幅值大于最优分割阈值的像素点作为目标像素点，从而得到各目标像素点。
57.为了判断各目标像素点是否具有空间位置上的相关性，也就是判断各目标像素点中是否是孤立存在的噪声点，本实施例基于各目标像素点的位置，采用dbscan聚类算法，dbscan聚类算法的邻域半径设置为3，对各目标像素点进行聚类操作，从而得到各个聚类簇。根据各个聚类簇内目标像素点的个数和预设目标像素点数目，将目标像素点的个数大于预设目标像素点数目的聚类簇剔除掉，得到满足预设聚类条件的各个聚类簇，所述预设聚类条件为聚类簇内目标像素点的数目大于预设数目，本实施例将预设数目设置为4。dbscan聚类算法实现聚类的过程为现有技术，不在本发明保护范围内，此处不再进行详细阐述。
58.（4-2-4）根据各个聚类簇内各目标像素点的坐标位置，确定各个聚类簇的最小外
接矩形，进而根据各个聚类簇的最小外接矩形内每个目标像素点的灰度值，确定各个聚类簇的最小外接矩形的灰度丰富指标，其步骤包括：首先，确定各个聚类簇的最小外接矩形，通过各个聚类簇中目标像素点的最大行坐标、最小行坐标、最大列坐标以及最小列坐标，构建各个聚类簇的最小外接矩形。构建最小外接矩形的过程为现有技术，不在本发明保护范围内，此处不再进行详细阐述。
59.根据各个聚类簇的最小外接矩形内每个目标像素点的灰度值，确定各个聚类簇的最小外接矩形的灰度丰富指标，其步骤包括：（4-2-4-1）根据各个聚类簇的最小外接矩形内每个目标像素点的灰度值，确定各个聚类簇的最小外接矩形内灰度类别的数目，进而确定各个灰度类别对应的像素灰度均值。
60.灰度类别是指目标像素点的不同灰度值，例如，某个聚类簇的最小外接矩形内每个目标像素点的灰度值为[3，5，6，6，8，8]，该聚类簇的最小外接矩形内灰度类别的数目为4个，分别为[3，5，6，8]，4个灰度类别对应的像素灰度均值的计算公式为(3+5+6+8)
÷
4=5.5，所以4个灰度类别对应的像素灰度均值为5.5。
[0061]
（4-2-4-2）根据各个聚类簇的最小外接矩形内每个目标像素点的灰度值，确定各个灰度类别的目标像素点在对应最小外接矩形中出现的频率。
[0062]
在本实施例中，统计各个灰度类别的像素点在各个目标像素点中的数目，将各个灰度类别的目标像素点数目除以目标像素点的总数目，将各个灰度类别的目标像素点数目比值作为各个灰度类别的目标像素点在对应最小外接矩形中出现的频率，得到各个灰度类别的像素点在对应最小外接矩形中出现的频率。
[0063]
（4-2-4-3）根据各个聚类簇的最小外接矩形内灰度类别的数目、各个灰度类别对应的像素灰度均值、各个灰度类别的像素点在对应最小外接矩形中出现的频率以及各个灰度类别的像素灰度，确定各个聚类簇的最小外接矩形的灰度丰富指标。
[0064]
灰度丰富指标是通过聚类簇内目标像素点的灰度级别数目、各个灰度类别的像素均值以及不同灰度类别的像素点在全部目标像素点中的占比，来评估聚类簇内目标像素点的灰度值的复杂情况，计算各个聚类簇的最小外接矩形的灰度丰富指标的计算公式为：其中，为各个聚类簇的最小外接矩形的灰度丰富指标，b为各个聚类簇的最小外接矩形内灰度类别的数目，为第i个灰度类别的像素点在对应最小外接矩形中出现的频率，为第i个灰度类别的像素灰度，为各个灰度类别对应的像素灰度均值。
[0065]
需要说明的是，某个聚类簇的最小外接矩形的灰度丰富指标越大，说明该聚类簇存在气泡缺陷的可能性越大，相反，则说明该聚类簇存在气泡缺陷的可能性越小。
[0066]
（4-2-5）获取灰度丰富指标阈值，根据各个聚类簇的最小外接矩形的灰度丰富指标和灰度丰富指标阈值，判断各个聚类簇的最小外接矩形中是否存在气泡缺陷。
[0067]
获取灰度丰富指标阈值，获取步骤（3-1-2）得到的u个正常图像块，将u设置为数值10，根据u个正常图像块内每个像素点的灰度值，参考步骤（4-2-2）至步骤（4-2-4）确定各个聚类簇的最小外接矩形的灰度丰富指标的过程，得到u个正常图像块对应的灰度丰富指标，
计算u个正常图像块对应的灰度丰富指标的平均值，记为，将1.5倍的u个正常图像块对应的灰度丰富指标的平均值作为灰度丰富指标阈值。通过待检测玻璃的灰度图像的各个正常区域的数据信息为基础，计算灰度丰富指标阈值，使所得到的灰度丰富指标阈值更加准确，提高了玻璃气泡缺陷检测的准确性。
[0068]
根据各个聚类簇的最小外接矩形的灰度丰富指标和灰度丰富指标阈值，判断各个聚类簇的最小外接矩形中是否存在气泡缺陷。将各个聚类簇的最小外接矩形的灰度丰富指标与灰度丰富指标阈值作比较，若某个聚类簇的最小外接矩形的灰度丰富指标大于灰度丰富指标阈值，则判定该聚类簇的最小外接矩形为气泡缺陷区域，否则，判定该聚类簇的最小外接矩形不为气泡缺陷区域。若各个聚类簇的最小外接矩形的灰度丰富指标均不大于灰度丰富指标阈值，则说明待检测玻璃不存在气泡缺陷。至此，本实施例完成了对待检测玻璃的气泡缺陷检测。
[0069]
以上所述实施例仅用以说明本技术的技术方案，而非对其限制；尽管参照前述实施例对本技术进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本技术各实施例技术方案的精神和范围，均应包含在本技术的保护范围之内。