一种基于乘积量化学习策略的缺陷检测算法的制作方法

1.本发明涉及数据处理与模式识别技术领域，尤其涉及一种基于乘积量化学习策略的缺陷检测算法，主要用于图像的检测与识别，特别是一种根据各自在图像信息中所反映的不同特征，把图像中不同的缺陷种类进行区别和定位的方法。它利用计算机对图像进行定量分析，把图像或图像中的每个像元或区域划归为若干个类别中的某一种，以代替人的视觉判读。实现了图像的高准确率识别。


背景技术：

2.缺陷检测被广泛使用于布匹瑕疵检测、工件表面质量检测、航空航天领域等。传统的算法对规则缺陷以及场景比较简单的场合，能够很好工作，但是对特征不明显的、形状多样、场景比较混乱的场合，则不再适用。近年来，基于机器学习的识别算法越来越成熟，许多公司开始尝试把深度学习算法应用到工业场合中。
3.当今社会，随着计算机技术，人工智能等科学技术的出现和发展，以及研究的深入，出现了基于机器视觉技术的表面缺陷检测技术。这种技术的出现，大大提高了生产作业的效率，避免了因作业条件，主观判断等影响检测结果的准确性，实现能更好更精确地进行表面缺陷检测，更加快速的识别产品表面瑕疵缺陷。产品表面缺陷检测属于机器视觉技术的一种，就是利用计算机视觉模拟人类视觉的功能，从具体的实物进行图象的采集处理、计算、最终进行实际检测、控制和应用。产品的表面缺陷检测是机器视觉检测的一个重要部分，其检测的准确程度直接会影响产品最终的质量优劣。由于使用人工检测的方法早已不能满足生产和现代工艺生产制造的需求，而利用机器视觉检测很好地克服了这一点，表面缺陷检测系统的广泛应用促进了企业工厂产品高质量的生产与制造业智能自动化的发展。
4.而哈希(hashing)和量化(quantization)是当前两种基于图像内容对图像进行检索的重要算法，这两种算法都是基于近似搜索理念提出的。哈希的方法是通过哈希函数把向量x变换成二值码(海明码)b(例如0101110101011100)，然后将距离dist(x_1，x_2)近似成二值码距离(海明距离)dist(b_1，b_2)。二值码距离可以通过popcount快速计算。yunchao gong等人提出的迭代量化的方法(iterative quantization，itq)的出发点是把二值编码当成原向量的近似，利用欧氏距离旋转不变性的性质，建立了最小化二值编码重建旋转原向量误差的目标函数，寻找最优的旋转变换和二值编码。尽管直观看上去重建向量的方法比保相似、保距离或者保序的方法简单，近似得更强，但itq实际效果还是很不错，原因是保相似、保距离或者保序需要建立二元或者多元关系，计算复杂度很大，从而需要各种近似，使得最后的效果不如预期。
5.近似最近邻搜索以其在大型数据集上的高效检索性能成为近年来研究的热点。人工神经网络搜索的目的是查找那些与查询数据的欧几里得距离是基础数据集中最小的实例。在神经网络研究中计算查询向量之间的欧几里得距离。而基础数据集中的所有向量都涉及到大量的计算，这对于大尺度和高维情况来说是不可行的。为了消除这些计算，已经提出并开发了许多数据结构和算法，其中一类是基于提升检索结构性能的方法，主要方法大
多基于树形结构。另一类主要基于对数据本身的处理，包括哈希算法、矢量量化方法等。
6.量化在近似最近邻搜索研究工作中发挥着重要作用，它通过采用数据表示策略来解决问题。首先采用无监督算法来实现数据聚类任务，使原始数据能够被标记。然后，由聚类中心表示或重建每个集群中的实例。因此，查询数据与基础数据之间的距离计算可以转换为查询数据和基础数据所属的类中心之间的距离计算。换句话说，查询和数据库向量之间的距离由查询数据向量和基础数据向量所属的聚类中心之间的距离近似。
7.乘积量化把原来的向量空间分解为若干个低维向量空间的笛卡尔积，并对分解得到的低维向量空间分别做量化。这样每个向量就能由多个低维空间的量化编码组合表示。自量化算法提出以来，开发了许多扩展算法来提高其搜索性能。笛卡尔k均值扩展了乘积量化算法，并对经典k均值目标函数的线性映射矩阵施加了柱正交约束，使优化过程比传统k均值更有效率。与笛卡尔k均值仅针对子码本进行优化不同，最佳乘积量化在子码本和空间分解方面进行优化，同时找到最优的空间分解方案和旋转矩阵。通过平衡协方差矩阵的特征值，最佳乘积量化可以得到最优空间分解，但强烈的多模态分布可能无法从这种分解中受益。最优笛卡尔k均值不同于以前的编码方案，在该方案中只选择了相应的子码本中的一个子码字，最优笛卡尔k均值使用多个子码字对数据点的子向量进行编码，性能为在近似最近邻搜索研究实验中，也提出了相同的算法。
8.上述算法在产品量化算法的研究上都有了显著的改进，但所有这些算法仍然属于无监督学习框架，这可能会明确限制这些算法的性能。


技术实现要素：

9.本部分的目的在于概述本发明的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本技术的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊，而这种简化或省略不能用于限制本发明的范围。
10.鉴于上述和/或现有的产品量化算法中存在的问题，提出了本发明。
11.因此，本发明其中的一个目的是提供一种基于乘积量化学习策略的缺陷检测算法，其能够有效地减少每个子空间的量化误差，提高识别性能。
12.为解决上述技术问题，本发明提供如下技术方案：一种基于乘积量化学习策略的缺陷检测算法，其包括，s1：基于局部线性嵌入的图像降维算法对高维图像进行降维操作，得降维后数据，将该数据作为输入样本，并给每个输入样本对应赋予一个其所属类别的类别标签；s2：根据输入样本的类标签，基于最优反向预测算法来设定目标函数；s3：特征空间最优分解，以得到半监督笛卡尔k均值模型；s4：采用拉普拉斯正则化的最优反向预测算法构建半监督的量化模型；s5：模型的优化以及基于图像的缺陷检测算法的实现。
13.作为本发明所述基于乘积量化学习策略的缺陷检测算法的一种优选方案，其中：所述最优反向预测算法的目标函数包含两项，分别为聚类中心矩阵与标签矩阵均未知的无监督聚类算法以及标签已知的监督学习约束项；
14.目标函数中的标签变量采用1
‑
k编码方案，可以通过优化聚类中心矩阵和未知标签矩阵变量之间的最小平方损失函数来求解最优反向预测算法，构建的目标函数如下，
[0015][0016]
其中和分别是训练实例矩阵和标记矩阵，是未标记的数据矩阵，是未知的标签矩阵，η2是交易参数，y
(l)
和b使用1
‑
k编码方案，p表示实例的维数，n
l
和n
u
分别是标记和未标记实例的数量，k是簇的数量。
[0017]
作为本发明所述基于乘积量化学习策略的缺陷检测算法法的一种优选方案，其中：所述特征空间的最优分解包括，
[0018]
s31：采用乘积量化算法，使得每个子空间中的码字均通过采用k均值聚类生成，该聚类相对于b和c迭代地优化平方失真误差；
[0019]
s32：笛卡尔k均值通过在公式(1.2)中的映射矩阵c的列上施加正交约束来获得ann搜索任务的空间分解方法；
[0020][0021]
s.t.b∈{0，1}
k
and||b||＝1
ꢀꢀ
(0.3)
[0022][0023]
s33：聚类中心的正交约束保证聚类中心表示为c≡rd，其中r是旋转矩阵，r
t
r＝rr
t
＝i，将公式(1.2)重新表示为公式(1.3)，最小化(1.3)关于r，d和b，所获得的旋转矩阵r和聚类中心d再结合公式(1.3)以获得失真误差；
[0024][0025][0026][0027]
式中，m是子空间的编号；
[0028]
s34：使用公式(1.1)在量化过程中代入公式(1.3)并给出半监督的笛卡尔k均值算法，给定标记数据集和未标记数据集其中p是实例的维数，k是量化中心的数量，n
l
和n
u
是标记和未标记数据实例的数量；
[0029]
s35：基于输入空间分解策略，半监督笛卡尔k均值(ssck)可以表示为
[0030][0031][0032][0033]
s36：利用希尔伯特
‑
施密特不相关性准则对特征空间进行分解。
[0034]
作为本发明所述基于乘积量化学习策略的缺陷检测算法的一种优选方案，其中：将拉普拉斯正则化项引入上述公式(1.4)，并得到以下公式：
[0035][0036]
[0037][0038]
其中，μ是平均值向量输入数据；y
(l)
和b都是量化标签且未知，l是拉普拉斯矩阵，l＝w
‑
d；
[0039]
w是相似矩阵，d
ii
＝∑
j
w
i，j
。
[0040]
作为本发明所述基于乘积量化学习策略的缺陷检测算法的一种优选方案，其中：优化公式(1.5)，且仅保留与y
(l)
相关的项：
[0041][0042]
进一步计算得：
[0043][0044][0045]
作为本发明所述基于乘积量化学习策略的缺陷检测算法的一种优选方案，其中：在获得量化y
(l)
的基础上，通过计算每个量化簇中所有标记数据x
(l)
的平均值来获得簇中心d；基于预测的聚类中心d，能够通过采用knn聚类算法获得未标记数据的标签矩阵b，如下：
[0046][0047]
作为本发明所述基于乘积量化学习策略的缺陷检测算法的一种优选方案，其中：根据标记数据x
(l)
，未标记数据x
(u)
及其量化标签y
(l)
和b，更新聚类中心d，如下：
[0048][0049][0050]
作为本发明所述基于乘积量化学习策略的缺陷检测算法的一种优选方案，其中：采用笛卡尔k
‑
means算法来更新r，基于：
[0051][0052]
得到：
[0053][0054]
作为本发明所述基于乘积量化学习策略的缺陷检测算法的一种优选方案，其中：在得到r，d，y
(l)
和b的基础上，能够更新变量μ，如下：
[0055]
μ
←
mean(x
‑
rdy)
[0056]
x＝[x
(l) ηx
(u)
] y＝[y
(l) ηb]。
[0057]
有益效果
[0058]
与现有技术相比，本发明所提供的一种基于乘积量化学习策略的缺陷检测算法具有如下有益效果：
[0059]
(1)将标记数据集成到量化步骤中，以提供标签信息，减少数据重建误差。
[0060]
(2)局部线性嵌入是建立在标记数据的基础上，并添加到最优反向预测函数，以获得半监督笛卡尔k均值。
[0061]
(3)给出了优化半监督笛卡尔k均值的算法，使函数能够得到最小值。
附图说明
[0062]
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。其中：
[0063]
图1为所述基于半监督学习的缺陷检测算法的流程图；
[0064]
图2、3和4为缺陷数据相关的系统界面和效果图。
具体实施方式
[0065]
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合说明书附图对本发明的具体实施方式做详细的说明。
[0066]
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是本发明还可以采用其他不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广，因此本发明不受下面公开的具体实施例的限制。
[0067]
其次，此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例，也不是单独的或选择性的与其他实施例互相排斥的实施例。
[0068]
如图1所示，为本发明的一个实施例，该实施例提供了一种基于乘积量化学习策略的缺陷检测算法，其是一种半监督笛卡尔k均值算法。本发明将量化步骤中传统的最小平方损失函数替换为最优反向预测损失函数，有标记的数据需要先通过拉普拉斯正则化，然后再用于模型训练，这样就有效地减少每个子空间的量化误差，提高检索性能。
[0069]
本发明主要采用三阶段方式，第一阶段：确定缺陷的位置和种类，采用最优反向预测算法作为检测手段；第二阶段：对第一阶段检测出的各个结构图像中的缺陷位置做定位检测，也是采用基于目标检测的方式；第三阶段：在第二阶段的基础上，对检测图像进行裁剪，将其送到分类模型进行分类，以确定最终的分类和识别结果。
[0070]
本发明的核心主要包括两部分，首先根据标记数据构造一个权值矩阵，主要特征为计算具有相同标签和不同标签样本之间的相似度值，给出两个具有相同标签的样本之间的相似度计算一个较大的值，而属于不同集群的数据则被赋予一个较小的值。然后在损失函数中加入拉普拉斯正则器，构造半监督笛卡尔k均值目标函数。
[0071]
具体的，一种基于乘积量化学习策略的缺陷检测算法，其主要包括如下步骤：
[0072]
s1：数据降维和特征的提取。基于局部线性嵌入的图像降维算法对高维图像进行降维操作，然后把这些降维后的数据作为输入样本，并给每个输入样本对应赋予一个其所属类别的类别标签，然后对输入样本(降维后的图像)进行特征提取。
[0073]
对于图像特征的提取，本专利基于hog(histogram of oriented gridients)和sift算法的特征融合提取图像特征。通过多种图像特征融合的方式对图像进行检索，对于
彩色图像，既要提取其颜色特征，也要提取其灰度图像特征，并分别建立彩色图像和灰度图像数据库，同时检索彩色图像库和灰度图像库，将检测的结果分开显示。
[0074]
s2：根据输入样本的类别标签，基于最优反向预测算法来设定目标函数。
[0075]
在现有技术中，xu等人提出了一种最优反向预测算法
[3]
，其目标函数包含两个项：一个是传统的k均值算法公式(一种无监督聚类算法)，其中聚类中心矩阵和标签矩阵都不知道；另一个是监督学习约束项，它与传统的k均值公式类似，但标签是已知的。目标函数中的标签变量采用1
‑
k编码方案。可以通过相对于聚类中心矩阵和未知标签矩阵变量迭代地优化最小平方损失函数来求解最优反向预测算法。
[0076]
目标函数中的标签变量采用1
‑
k编码方案，可以通过优化聚类中心矩阵和未知标签矩阵变量之间的最小平方损失函数，使得聚类中心和未知标签类之间的差异不断缩小，以此来求解最优反向预测算法，构建的目标函数如下：
[0077][0078]
其中种分别是训练实例矩阵和标记矩阵，是未标记的数据矩阵，是未知的标签矩阵，η2是交易参数，y
(l)
和b使用1
‑
k编码方案，p表示实例的维数，n
l
和n
u
分别是标记和未标记实例的数量，k是簇的数量。最优反向预测算法是通过最优反向预测的概念统一几种有监督和无监督的训练原则：预测来自目标标签的输入，优化模型参数和任何缺失标签。监督最小二乘法、主成分分析、k均值聚类和归一化图切割都可以表示为相同训练原则的实例。
[0079]
s3：特征空间最优分解，以得到半监督笛卡尔k均值模型。所述特征空间的最优分解具体包括如下步骤：
[0080]
s31：采用乘积量化算法，使得每个子空间中的码字均通过采用k均值聚类生成，该聚类相对于b和c迭代地优化平方失真误差。本专利以基本的乘积量化算法作为基础，对现有的各量化算法进行分析和改进提升。
[0081]
在计算笛卡尔k均值的过程中，乘积量化(pq)算法被用于ann研究任务。高维度输入数据空间被平均分解并表示为m个低维度子空间的笛卡尔积。每个子空间可以通过使用传统的k均值算法生成码本。因此，将为每个子向量生成k个子码字，通过这种方式，m个子向量将生成仅具有o(kp)存储的k
m
簇，而如果我们采用传统方法通过k
‑
对整个输入数据进行编码意味着需要存储o(k
m
p)，具有相同数量的集群。同时，计算复杂度从o(k
m
p)降低到o(kp)。在pq算法中，每个子空间中的码字是通过采用k均值聚类生成的，该聚类相对于b和c迭代地优化平方失真误差。然而，pq算法没有提出如何获得ann搜索任务的最优空间分解的方法。
[0082]
s32：笛卡尔k均值通过在公式(1.2)中的映射矩阵c的列上施加正交约束来获得ann搜索任务的空间分解方法，解决了上述问题，以隐含地调整实例的维度信息，这使得关于b的优化在正交笛卡尔k均值
[4]
中变得更易处理。
[0083][0084]
s.t.b∈{0，1}
k and||b||＝1
ꢀꢀ
(0.14)
[0085]
[0086]
s33：聚类中心的正交约束保证聚类中心表示为c≡rd，其中r是旋转矩阵，r
t
r＝rr
t
＝i，因此公式(1.2)可以重新表示为公式(1.3)，最小化公式(1.3)关于r、d和b，所获得的最佳旋转矩阵r和聚类中心d可以帮助公式(1.3)获得较低的失真误差。
[0087][0088][0089][0090]
式中，m是子空间的编号。
[0091]
s34：构建模型与定义。半监督笛卡尔k均值计算，经典k均值，主成分分析(pca)和归一化切割可以认为是最优反向预测算法的特例。基于这个概念，可以使用公式(1.1)在量化过程中代入公式(1.3)并给出半监督的笛卡尔k均值算法。给定标记数据集和未标记数据集其中p是实例的维数，k是量化中心的数量，n
l
和n
u
是标记和未标记数据实例的数量。
[0092]
s35：基于输入空间分解策略，半监督笛卡尔k均值(ssck)可以表示如下：
[0093][0094][0095][0096]
s36：利用希尔伯特
‑
施密特不相关性准则对特征空间进行分解。
[0097]
本专利将利用基于希尔伯特
‑
施密特不相关性准则(hilbert
‑
schmidt independence criterion，hsic准则构建优化模型来实现特征空间的分解。要达到的目的是：不同子空间之间的数据维数是不相关的，而同一个子空间内部的数据维数是严格相关的。方案是：给定归一化后的训练数据集其中n为样本个数，d
★
为样本的维数，将数据x基于维数进行二等分割得到两个子矩阵和这里求解一个线性映射矩阵p1，使得映射后的协方差矩阵呈现对角矩阵块形式，即x1和x2是不相关的：
[0098][0099]
为此就必须使得x1和x2具有最小的相关性。为此本专利参考了周志华老师文章中的算法
‑‑‑
基于希尔伯特
‑
施密特不相关性准则(hilbert
‑
schmidt independence criterion，hsic)来构建数学模型求解映射矩阵p1，使得子空间x1和x2具有最小的相关性。给出优化模型如下所示：
[0100][0101]
其中h＝[h
ij
]
n
×
n
和δ
ij
＝1如果i＝＝j，否则δ
ij
＝0。关于模型最优解的求解，文献
[15]
中的mddm算法为我们提供了很好的解决思路。在得到最优值后，将x1和x2进行重新映射得到两个不相关的子空间和然后在x1和x2的基础上按照上述方法继续进行分解，由此可以分解成为4个不相关的子空间，而子空间内部数据是
相关的，依次类推，我们就可以得到不相关性子空间的最优分解，在模型建模过程中为了有效的降低码本与数据之间误差，使用旋转矩阵对数据进行相应优化旋转也是所必须的。在图(11)中我们给出了基于希尔伯特
‑
施密特不相关性准则对特征空间分解成为8个子空间的算法流程图。图中为原始数据，可称之为第0层，其中d
★
为数据的维数，n为样本数量。在第一层分解中将x1分解成为两层和其所对应的原始数据的维数下标分别为种然后对和分别进行分解为与和与其所对应的原始数据维数下标依次为和第三层分解算法依此类推。
[0102]
s4：采用拉普拉斯正则化的最优反向预测算法来构建半监督的量化模型；
[0103]
以量化算法中的笛卡尔k均值为例，构建半监督笛卡尔k均值模型，这里将其称之为第i类模型，其算法可以表示成为如下数学模型的优化问题：
[0104][0105][0106][0107]
在上述公式中x
(l)
和x
(u)
分别为有标签数据集和无标签数据集，y
(l)
和b分别对应是有标签数据集x
(l)
和无标签数据集x
(u)
对应的码本标签矩阵，码本标签矩采用的编码形式是1
‑
of
‑
k形式，l为拉普拉斯矩阵，λ和η为权重因子。在上式中模型的优化问题需要对4个变量，即旋转矩阵r、码本d，无标签数据集的码本标签矩阵b以及有标签数据集的码本标签矩阵y
(l)
，分别做迭代优化。
[0108]
基于空间概念：
[0109][0110]
||
·
||
f
代表是frobenius范数，公式(1.4)给出的标记数据x
(l)
，y
(l*)
和未标记数据x
(u)
，这样就形成了半监督笛卡尔k均值算法，更具体地说，公式(1.4)只是量化问题，是笛卡尔k均值的一部分。
[0111]
但是，暂时还不能使用公式(1.4)量化标记和未标记的数据，因为此时不知道量化标签矩阵y
(l*)
。当在聚类任务中使用最佳反向预测时，公式(1.1)中的标记矩阵y
(l)
的符号表示聚类标签或分类标签，其可以从已知的标记数据获得。然而，本专利提出的半监督笛卡尔k均值公式(1.4)中使用的最优反向预测，其作用是量化或编码数据时，标签矩阵y
(l)
用于指示量化标签，但它是未知的。
[0112]
换句话说，量化标签与聚类标签不同，聚类标签可以直接从标记数据中获得，但量化标签不能，因此公式(1.4)不能用于本专利提出的半监督笛卡尔k均值。为了解决这个问题，本专利进一步将拉普拉斯正则化项引入上述半监督笛卡尔k均值模型(见公式1.4)中，并得到以下公式：
[0113][0114][0115][0116]
其中和μ是平均值向量输入数据。y
(l)
和b都是量化标签且未知，l是拉普拉斯矩阵，l＝w
‑
d。w是相似矩阵，d
ii
＝∑
j
w
i，j
。
[0117]
s5：模型的优化以及基于图像的缺陷检测算法的实现。
[0118]
我们使用监督方法构造相似性矩阵w，如果x
i
和x
j
具有相同的聚类标记，则条目w
ij
将被赋予大的权重，并且如果x
i
和x
j
属于不同的聚类，则将给出小的权重。优化公式(1.5)通常是一项难以处理的工作，因为y
(l)
和b是具有1
‑
k编码方案的离散矩阵，离散优化是np难问题并且在非模块问题中是不同的。
[0119]
这里提供两种方法来解决这个问题。第一个：采用了详尽搜索方法。我们迭代地彻底检查y
(l)
中的所有元素，找到元素使得目标值最小化，并将所有其他元素设置为0；第二个：首先放松将离散变量y
(l)
变为连续变量并求解y
(l)
，然后利用约束优化笛卡尔k均值目标函数得到最优的y
(l)
，采用阈值法得到离散变量。
[0120]
一、更新y
(l)
：为了使用穷举搜索方法针对y
(l)
优化公式(1.5)，我们可以重写公式(1.5)如下(仅保留与y
(l)
相关的那些项)：
[0121][0122]
为了得到y
(l)
最优值，经过进一步计算如下公式：
[0123][0124][0125]
二、更新b：在获得量化y
(l)
之后，本专利通过计算每个量化簇中所有标记数据x
(l)
的平均值来获得簇中心d.基于预测的聚类中心d，然后通过采用knn聚类算法获得未标记数据的标签矩阵b.更准确地说，通过解决以下问题可以获得关于b的优化，这样才能得到更加鲁棒的特征矩阵：
[0126][0127][0128]
三、更新d：根据标记数据x
(l)
，未标记数据x
(u)
及其量化标签y
(l)
和b，为了适应不断变化的数据和调整聚类中心，我们可以通过求解以下内容来更新聚类中心d：
[0129][0130]
[0131]
四、更新r：已经提出了许多算法来解决正交约束优化问题，在实验中，本专利仍然采用笛卡尔k
‑
means算法来解决它的高性能和高效率，通过更新r，能够获得更好的优化策略。
[0132]
基于：
[0133][0134]
于是：
[0135][0136]
五、更新μ：得到r，d，y
(l)
和b后，为了得到最优的参数，本专利使用以下更新变量μ：
[0137]
μ
←
mean(x
‑
rdy)
[0138]
x＝[x
(l) ηx
(u)
] y＝[y
(l) ηb]
[0139]
实验分析。为了验证本专利提出的半监督笛卡尔k均值的性能，在公共数据集上进行了一系列ann搜索实验，验证在3个数据集上的表现。本实验拟从召回率、检索精度这2个指标来评价算法的性能。
[0140]
本发明实施例中，提出了一种基于笛卡尔k均值和正交最优反向预测算法的半监督笛卡尔k均值算法。为了提高传统笛卡尔k均值算法的性能，我们利用样本的标签信息来构造相似矩阵——那些标签相同的样本被赋予较高的权重值，同时将较低的权重值赋给标签不同的样本。标签相同的样本将被分配到量化步骤中的相同或相邻类中聚集。然而，对标签矩阵的目标函数进行优化是一个离散问题，很难解决。于是我们分别优化了方程，并在公共数据集上验证了该方程的性能。实验结果表明，该算法比比较算法性能更优。我们在半监督笛卡尔k均值算法的研究中取得了较好的结果。如图2、3和4所示，给出了系统和算法发效果图。
[0141]
应说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。