CN109117860A - 一种基于子空间投影和字典学习的图像分类方法 - Google Patents
一种基于子空间投影和字典学习的图像分类方法 Download PDFInfo
- Publication number
- CN109117860A CN109117860A CN201810677726.7A CN201810677726A CN109117860A CN 109117860 A CN109117860 A CN 109117860A CN 201810677726 A CN201810677726 A CN 201810677726A CN 109117860 A CN109117860 A CN 109117860A
- Authority
- CN
- China
- Prior art keywords
- sample
- dictionary
- test set
- matrix
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于子空间投影和字典学习的图像分类方法,首先通过带标签的训练集样本对判别字典进行初始化,然后使用判别字典预测测试样本的类别标签。选出伪标签具有高可信度的测试集样本,并用其与带标签的训练集样本学习一个低维子空间,在这个低维空间上更新判别字典。利用更新得到的判别字典重新对测试集样本进行分类,比较本轮迭代得到的伪标签和上一轮伪标签,将两轮得到的伪标签属性相同的样本称为稳定样本,如果在一轮迭代结束后稳定样本的数量超过了测试集样本数量的百分之八十,则迭代结束,输出该轮得到的伪标签作为分类的结果。与已有的领域自适应图像分类方法相比,本发明的算法能取得更高的分类准确率。
Description
技术领域
本发明涉及到模式识别中的领域自适应图像分类领域,尤其涉及一种基于子空间投影和字典学习的域自适应图像分类领域。
背景技术
在传统的模式识别和机器学习中,算法需要在一定的假设条件下才能有良好的表现,最常见的是:源域的训练集样本和目标域的测试集样本在分布上要保持一致。而在实际应用中,往往很难满足这个条件。例如:在图像分类领域,图像传感器的不同种类,拍摄角度的不同,不同的光线条件等,都会使训练集和测试集样本之间存在较大的分布差异。分布差异会导致训练好的分类器在实际处理测试集样本时表现不佳,大多数方法需要重新训练分类器,而训练分类器需要充足的带标签样本,在某些应用场景下,收集带标签的样本需要付出昂贵的代价甚至无法完成,领域自适应方法就致力于解决这个问题。
解决领域自适应问题的其中一种方法是基于特征的自适应,在图像分类领域,一幅图像会通过各种各样的特征提取方法将其从像素矩阵转换为特征向量,面对图像分布不同带来的差异,这类方法从特征向量着手,通过子空间投影,特征变换等方式来减小图像分布差异给分类带来的影响。较有代表性的是杨蒙(音)教授提出的FDDL字典学习方法和Boqing Gong等人提出的GFK方法。FDDL方法使用线性判别条件,对每一类的样本构造一个子字典,并通过整个字典将样本线性表示出来。GFK方法基于核方法,将数据嵌入到Grassmann流形中,取得了不错的效果。最近效果比较好的方法是S.Herath等人提出的额ILS方法,即学习一个不变的潜在子空间,将训练集样本和测试集样本投影到该子空间上,从而减小了分布的差异。
发明内容
本发明所要解决的技术问题针对背景技术中所涉及到的缺陷,提供一种基于子空间投影和字典学习的图像分类方法,解决了在不重新训练分类器的条件下,源域的训练集和目标域的测试集样本分布过大带来的分类器表现不佳的问题。
为了达到以上目的,本发明提出了一种基于子空间投影和字典学习的图像分类方法,包括以下步骤:
步骤1)使用带标签的训练集初始化判别字典并给测试集样本加上伪标签;
步骤2)选出具有高可信度的测试集样本;
步骤3)将高可信度的测试集样本和训练集样本分别投影到低维平面上,在低维平面上学习判别字典,同时联合学习投影矩阵;
步骤4)使用学习到的判别字典再次对测试集样本进行分类,即重新加上伪标签;
步骤5)计算稳定样本个数,若其数量大于等于测试集样本数量的80%则结束迭代过程,输出该轮伪标签作为结果;若小于80%则重复步骤2~5直到满足条件。
进一步的,所述步骤1)的详细步骤如下:
假设给定的训练集中包含C类样本,
即标签L={1,…,C},是训练集样本,是测试集样本,Ns,Nt分别是训练集和测试集样本的数量,n1,n2分别是训练集和测试集样本的特征维数;用 表示第k轮的判别字典,表示第j类的子字典,其中n为低维子空间的维数,K是整个字典的原子数,X=[x1,…,xN]∈RK×N是N个样本由判别字典D的稀疏表示;‖·‖F指矩阵的Frobenius范数,即指矩阵各元素平方和开平方;‖·‖2是矩阵的2范数,即指该矩阵的协方差矩阵最大特征值开平方;
步骤1.1)用带标签的训练集样本,通过线性判别字典学习的方法来初始化判别字典D(0),即通过最小化下式来使字典能够很好的使用原子线性表达出样本,并且能够通过稀疏表示系数X来对样本进行分类:
其中对第i个样本Xin和Xo1t分别为
步骤1.2)学习到字典D(1)后,给测试集样本加上伪标签,给第i个测试集样本加伪标签的方法如式(2):
其中是使用第j类子字典重建第i个测试集样本的稀疏系数。
进一步的,所述步骤2)的详细步骤如下:
由步骤1.2)给第i个测试集样本加上伪标签之后,通过计算每个训练集样本的置信度来衡量伪标签的可信度;
在第k轮迭代过程中,第i个样本被分类到第j类的置信度由式(3)定义:
即对第i个测试集样本,现得到长度为C的向量其值除了伪标签所属的项外都为零,即若共有C=5类,第i个样本伪标签为2,则其中为使用第k轮学习到的字典中的第j类子字典重建第i个测试集样本的重建误差,即在第k轮迭代学习到更新后的字典后,分别使用第j个子字典来重建第i个测试集样本,由式(4)给出:
其中是将测试集样本投影到低维子空间的投影矩阵;
为了得到伪标签可信度最大的Q%*Nt个测试集样本,以置信度为依据选择矩阵挑出所有类样本中可信度最大的前Q%个样本,W的每一个元素,即为零一变量,指出第i个样本在挑选第j类时是否被选中,若被选中则为1,否则为0,其中Q是一个固定常数,根据训练集大小而定,一般取50%;
挑选第j类样本的方法如式(5):
s.t.||Wj (k)||0=Q%*Nt,j=1,…,C (5)
其中是第j类样本的选择矩阵,它是一个对角矩阵,其对角线上的元素为选择矩阵第j列的元素;同样地,也是对角矩阵,其对角线上的元素是置信度矩阵C(k)的第j列;||Wj (k)||0是Wj (k)矩阵中非零元素的个数,即确保每一类中只有Q%被选出;
至此得到每一类的选择矩阵Wj (k),通过式(6)将选出的每一类高可信度样本加入字典的训练过程中,即更新用于下一轮训练的测试集样本:
进一步的,所述步骤3)的详细步骤如下:
步骤3.1)假设投影矩阵和n为投影后子空间的维数,则投影到低维子空间上后的训练集样本为测试集样本为
在这个低维平面上利用式(7)训练集样本和挑选之后的测试集样本学习字典,并联合学习投影矩阵Ps和Pt
且||xj||0≤T0,T0为稀疏系数
其中即将两个投影矩阵邻接在一起,同理
其次,即将Ys和Yt作为块矩阵的对角,其余的块矩阵为零,Trace指矩阵的迹,即矩阵主对角线上元素之和;
f1使用利用PsPt将分别将训练集样本和测试集样本投影到低维子空间中,并在低维子空间中联合学习判别字典D和‖·‖0指矩阵或向量的零范数,即向量中非零元素的个数,‖·‖F指矩阵的Frobenius范数,即指矩阵各元素平方和开平方;
f2是训练集样本和测试集样本的最大均值距离(MMD),衡量了训练集样本和测试集样本在分布上的差异,分布差异越小,f2的取值越接近0;其中矩阵由式(8)给出:
f3是将样本反投影的过程,期望在投影过程中损失尽可能少的样本信息;
步骤3.2)在求解式(7)的优化问题时引入矩阵 和B,
分别满足:
Ps=(YsAs)T(9)
Pt=(YtAt)T (10)
式(7)的优化方程可转换为:
且||xj||1≤T0,T0为稀疏系数
经过整理式(12)可得到:
需要求的变量有三个,在这里本发明采取固定其他先求其中一个的策略,即先固定求再迭代固定求当被固定时,它们可被看作常数,固式中只剩下为未知量;
首先将初始化,方法是:
先将K进行特征值分解得到:
K=VSVT (14)
其中,V是特征向量矩阵,S是一个对角矩阵,其对角线上的元素是矩阵K的特征值,本方法将As及At的列向量初始化为前n大的特征值所对应的特征向量,即对特征值从大到小排列之后,选取最大的前n个特征值对应的特征向量构成或这里的n是将为后平面的维数;
在得到的初始值后,式中还有未知量和由于
未知量可转换为D和我们利用多次迭代的方法求解,求解D和的具体步骤如下:
将D初始化为随机的单位向量,则式中只剩下为未知量,先将初始化为全O向量再采取逐列更新的方式来计算 为第i类样本系数表示的稀疏系数;在第h轮中,更新的公式为:
其中,是第h轮更新得到的稀疏系数是函数Q(x)对的偏导数,其中,Q(X)由下式给出:
由式(15)解得后,固定式中的同样利用逐列更新的方法来求解D,即每次更新字典D中的一列dj,更新dj的公式为:
其中,Z是经过投影的样本即
至此,求得了D(k)和利用下式求出
则固定和求解
由式(13)(14)可令:
可将式(13)转换为如下的形式:
{G}=argminGtrace(GTHG) (17)
其中G=[Gs,Gt],式(16)需要利用Stiefel流形来求解,可直接调用现成的软件包,求得G;
至此得到更新的公式为:
可由式(9)和(10)分别解出Ps和Pt,
至此,已由式(16)求得D(k),由式(9)和(10)分别解出Ps和Pt。
进一步的,所述步骤4)的详细步骤如下:
在得到本轮求得的字典D(k)以及投影矩阵后,用下式更新第i个测试集样本的伪标签:
进一步的,所述步骤5)的详细步骤如下:
将所有的测试集样本的伪标签都更新完后,将其与该样本上一轮的伪标签进行比较,若两轮得到的伪标签相同,则称该样本为稳定样本,计算稳定样本的数量;
若该数量大于测试集样本总数的Q%,则样本和字典都趋于稳定,结束迭代并将最后一轮测试集样本的伪标签输出,作为最后的结果;
若稳定样本数量小于测试集样本总数的Q%,则继续迭代直到不满足条件为止。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
1.能够通过联合学习子空间和字典的方法,显著的减小训练集和测试集样本之间的分布差异。
2.能够挖掘测试集样本中潜在的判别信息来加入训练分类器过程,即伪标签,并使其具有较高的可信度。
3.能够处理训练集样本和测试集样本异构的情况,即当训练集样本转化的特征向量与测试集样本转化的特征向量具有不同的维数时,多数算法无法处理,而本发明具有较好的应对方法。
附图说明
图1是本发明基于子空间投影和字典学习的图像分类的工作流程图。
具体实施方式
针对迁移学习中的领域自适应问题(源域为少量带标签的训练集图像,目标域为充足的不带标签的测试集图像,训练集和测试集图像之间存在较大的分布差异),我们的目标是通过算法,使得在训练集上训练的分类器能够在目标域的测试集样本上有良好的分类效果。
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明公开了一种基于子空间投影和字典学习的图像分类方法:
在使用训练集样本初始化字典后,使用字典给测试集样本加上伪标签,并挑选出伪标签具有高可信度的测试集样本加入迭代训练字典的过程。将带有标签的训练集样本和带有伪标签的测试集样本通过两个投影矩阵分别投影到低维子空间上,并在子空间上联合学习判别字典与投影矩阵,使用学习到的判别字典对所有的测试集样本进行重分类得到各自新的伪标签,则完成一轮迭代,每次迭代得到新的判别字典,新的投影矩阵和新的测试集样本伪标签。本发明将新伪标签与上轮伪标签相同的测试集样本称为稳定样本,每轮迭代结束后计算稳定样本的数量,如果稳定样本的数量小于测试集样本总数的百分之八十,则继续迭代过程,如果在一轮迭代之后稳定样本的数量占测试集样本总数的百分之八十及以上,则迭代结束,输出本轮得到的测试集样本伪标签作为最后的分类结果。
本发明所述的基于子空间投影和字典学习的图像分类方法,包括以下步骤:
(一)初始化判别字典并使用字典给测试集样本加上伪标签
假设给定的训练集中包含C类样本即标签L={1,…,C}, 是训练集样本,是测试集样本。Ns,Nt分别是训练集和测试集样本的数量。n1,n2分别是训练集和测试集样本的特征维数。用表示第k轮的判别字典,表示第j类的子字典,其中n为低维子空间的维数,K是整个字典的原子数。最后,X=[x1,…,xN]∈RK×N是N个样本由判别字典D的稀疏表示。‖·‖F指矩阵的Frobenius范数,即指矩阵各元素平方和开平方。‖·‖2是矩阵的2范数,即指该矩阵的协方差矩阵最大特征值开平方。
步骤1.1)用带标签的训练集样本,通过线性判别字典学习的方法来初始化判别字典D(0),即通过最小化下式来使字典能够很好的使用原子线性表达出样本,并且能够通过稀疏表示系数X来对样本进行分类:
其中对第i个样本Xin和X0ut分别为
可逐列对D进行更新得到D(1)
步骤1.2)学习到字典D(1)后,给测试集样本加上伪标签,给第i个测试集样本加伪标签的方法如式(2):
其中是使用第j类子字典重建第i个测试集样本的稀疏系数。
将所有的测试集样本加上标签之后得到即为第一轮的伪标签。
(二)选出高可信度的测试集样本
给第i个测试集样本加上伪标签之后,本发明通过计算每个训练集样本的置信度来衡量伪标签的可信度。
在第k轮迭代过程中,第i个样本被分类到第j类的置信度由式(3)定义:
即对第i个测试集样本,现得到长度为C的向量其值除了伪标签所属的项外都为零,即若共有C=5类,第i个样本伪标签为2,则
其中为使用第k轮学习到的字典中的第j类子字典重建第i个测试集样本的重建误差,即在第k轮迭代学习到更新后的字典后,分别使用其所有的子字典来重建第i个测试集样本。由式(4)给出:
其中是将测试集样本投影到低维子空间的投影矩阵
为了得到伪标签可信度最大的Q%*Nt个测试集样本,本发明以置信度为依据选择矩阵挑出所有类样本中可信度最大的前Q%个样本,W的每一个元素,即为零一变量,指出第i个样本在挑选第j类时是否被选中。若被选中则为1,否则为0,其中Q是一个固定常数,根据训练集大小而定,一般取50%。
挑选第j类样本的方法如式(5):
其中是第j类样本的选择矩阵,它是一个对角矩阵,其对角线上的元素为选择矩阵第j列的元素。同样地,也是对角矩阵,其对角线上的元素是置信度矩阵C(k)的第j列。||Wj (k)||0是Wj (k)矩阵中非零元素的个数,即确保每一类中只有Q%被选出。
至此得到每一类的选择矩阵Wj (k),通过式(6)将选出的每一类高可信度样本加入字典的训练过程中,即更新用于下一轮训练的测试集样本:
Yj (k+1)=[YtWj (k)],j=1,…,C (6)
(三)投影并联合学习投影矩阵和字典
假设投影矩阵和n为投影后子空间的维数,则投影到低维子空间上后的训练集样本为测试集样本为
在这个低维平面上利用式(7)训练集样本和挑选之后的测试集样本学习字典,并联合学习投影矩阵Ps和Pt
且||xj||0≤T0,T0为稀疏系数
其中即将两个投影矩阵邻接在一起,同理其次,即将Ys和Yt作为块矩阵的对角,其余的块矩阵为零。Trace指矩阵的迹,即矩阵主对角线上元素之和。
f1使用利用PsPt将分别将训练集样本和测试集样本投影到低维子空间中,并在低维子空间中联合学习判别字典D和指矩阵或向量的零范数,即向量中非零元素的个数,‖·‖F指矩阵的Frobenius范数,即指矩阵各元素平方和开平方。
f2是训练集样本和测试集样本的最大均值距离(MMD),衡量了训练集样本和测试集样本在分布上的差异,分布差异越小,f2的取值越接近0。其中矩阵由式(8)给出:
f3是将样本反投影的过程,期望在投影过程中损失尽可能少的样本信息
在求解式(7)的优化问题时引入矩阵 和B,分别满足:
Ps=(YsAs)T
(9)Pt=(YtAt)T(10)
于是式(7)的优化方程可转换为:
且||xj||1≤T0,T0为稀疏系数
经过整理式(12)可得到:
其中M为之前提到过的MMD矩阵。
需要求的变量有三个,在这里本发明采取固定其他先求其中一个的策略,即先固定求再迭代固定求当被固定时,它们可被看作常数,式中只剩下为未知量。
首先我们将初始化,方法是:
将K进行特征值分解得到:
K=VSVT(14)
其中,V是特征向量矩阵,S是一个对角矩阵,其对角线上的元素是矩阵K的特征值,本方法将As及At的列向量初始化为前n大的特征值所对应的特征向量,即对特征值从大到小排列之后,选取最大的前n个特征值对应的特征向量构成或这里的n是将为后平面的维数。
在得到的初始值后,式中还有未知量和由于
未知量可转换为D和我们利用多次迭代的方法求解,求解D和的具体步骤如下:
我们将D初始化为随机的单位向量,则式中只剩下为未知量,我们先将初始化为全O向量再采取逐列更新的方式来计算 为第i类样本系数表示的稀疏系数。
在第h轮中,更新的公式为:
其中,是第h轮更新得到的稀疏系数,是函数Q(x)对的偏导数,式中Q(x)由下式给出
由式(15)解得后,固定式中的同样利用逐列更新的方法来求解D,即每次更新字典D中的一列dj,更新dj的公式为:
其中,Z是经过投影的样本即
至此,我们已经将求得了D(k)和可利用下式求出
则我们可以固定和求解
由式(13)(14)可令:
可将式(13)转换为如下的形式:
{G}=argminGtrace(GTHG) (17)
其中G=[Gs,Gt],式(16)需要利用Stiefel流形来求解,可直接调用现成的软件包,求得G。
至此得到更新的公式为:
可由式(9)和(10)分别解出Ps和Pt。
至此,我们已由式(16)求得D(k),由式(9)和(10)分别解出Ps和Pt。
(四)更新测试集样本的伪标签
得到本轮求得的字典D(k)以及投影矩阵后,用下式更新第i个测试集样本的伪标签:
(五)判定稳定样本的数量
将所有的测试集样本的伪标签都更新完后,将其与该样本上一轮的伪标签进行比较,若两轮得到的伪标签相同,则称该样本为稳定样本,计算稳定样本的数量。
若该数量大于测试集样本总数的Q%,则样本和字典都趋于稳定,结束迭代并将最后一轮测试集样本的伪标签输出,作为最后的结果。
若稳定样本数量小于测试集样本总数的Q%,则继续迭代直到不满足条件为止。
将本发明所述的基于子空间投影和字典学习的图像分类方法在OFFICE和Caltech-256数据库上进行实验,并将实验结果与其他相关的域自适应方法进行对比分析。
OFFICE数据库包含三个域的图片,即有三组样本可作为训练集或测试集,分别为Amazon,Webcam,DSLR。Amazon(以下称A)数据集是从网络商店上下载来的一些图片,DSLR(以下称D)数据集是用SLR相机拍出的高分辨率照片,Webcam(以下称W)数据集是网络相机拍出来的低分辨率照片,Caltech-256(以下称C)数据集是加州理工大学收集的一些办公物品的照片。每张图片都提取它的SURF特征并使用K-means聚类,将其转换为一个800维的向量。我们把四个数据集分别当作源域的训练集和目标域的测试集,即可有种排列,我们选取了其中的一部分作为实验,实验的平均分类正确率如下:
表一:各方法在OFFICE和Caltech-256数据库上的分类结果
A→W | W→D | W→A | W→C | D→A | C→A | C→D | |
FDDL | 31.2 | 60.3 | 30.6 | 22.8 | 32.1 | 38.1 | 38.2 |
GFK | 35.7 | 71.2 | 35.5 | 29.3 | 36.2 | 40.4 | 41.1 |
ILS | 40.6 | 72.4 | 38.6 | 32.6 | 38.9 | 48.6 | 44.4 |
本发明 | 41.8 | 77.2 | 42.8 | 33.8 | 39.1 | 48.3 | 50.2 |
领域从表一中可以看出,本发明所提出的方法相对于其他方法,有较高的分类正确率。
本技术领域技术人员可以理解的是,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于子空间投影和字典学习的图像分类方法,其特征在于:包含以下步骤:
步骤1)使用带标签的训练集初始化判别字典并给测试集样本加上伪标签;
步骤2)选出具有高可信度的测试集样本;
步骤3)将高可信度的测试集样本和训练集样本分别投影到低维平面上,在低维平面上学习判别字典,同时联合学习投影矩阵;
步骤4)使用学习到的判别字典再次对测试集样本进行分类,即重新加上伪标签;
步骤5)计算稳定样本个数,若其数量大于等于测试集样本数量的80%则结束迭代过程,输出该轮伪标签作为结果;若小于80%则重复步骤2~5直到满足条件。
2.根据权利要求1所述的基于子空间投影和字典学习的图像分类方法,其特征在于:所述步骤1)的详细步骤如下:
假设给定的训练集中包含C类样本,
即标签L={1,…,C},是训练集样本,是测试集样本,Ns,Nt分别是训练集和测试集样本的数量,n1,n2分别是训练集和测试集样本的特征维数;
用表示第k轮的判别字典,表示第j类的子字典,其中n为低维子空间的维数,K是整个字典的原子数,X=[x1,…,xN]∈RK×N是N个样本由判别字典D的稀疏表示;||·||F指矩阵的Frobenius范数,即指矩阵各元素平方和开平方;||·||2是矩阵的2范数,即指该矩阵的协方差矩阵最大特征值开平方;
步骤1.1)用带标签的训练集样本,通过线性判别字典学习的方法来初始化判别字典D(0),即通过最小化下式来使字典能够很好的使用原子线性表达出样本,并且能够通过稀疏表示系数X来对样本进行分类:
其中对第i个样本Xin和Xout分别为
步骤1.2)学习到字典D(1)后,给测试集样本加上伪标签,给第i个测试集样本加伪标签的方法如式(2):
其中是使用第j类子字典重建第i个测试集样本的稀疏系数。
3.根据权利要求2所述的基于子空间投影和字典学习的图像分类方法,其特征在于:所述步骤2)的详细步骤如下:由步骤1.2)给第i个测试集样本加上伪标签之后,通过计算每个训练集样本的置信度来衡量伪标签的可信度;
在第k轮迭代过程中,第i个样本被分类到第j类的置信度由式(3)定义:
即对第i个测试集样本,现得到长度为C的向量其值除了伪标签所属的项外都为零,即若共有C=5类,第i个样本伪标签为2,则其中为使用第k轮学习到的字典中的第j类子字典重建第i个测试集样本的重建误差,即在第k轮迭代学习到更新后的字典后,分别使用第j个子字典来重建第i个测试集样本,由式(4)给出:
其中是将测试集样本投影到低维子空间的投影矩阵;
为了得到伪标签可信度最大的Q%*Nt个测试集样本,以置信度为依据选择矩阵挑出所有类样本中可信度最大的前Q%个样本,W的每一个元素,即为零一变量,指出第i个样本在挑选第j类时是否被选中,若被选中则为1,否则为0,其中Q是一个固定常数,根据训练集大小而定,一般取50%:
挑选第j类样本的方法如式(5):
其中是第j类样本的选择矩阵,它是一个对角矩阵,其对角线上的元素为选择矩阵第j列的元素;同样地,也是对角矩阵,其对角线上的元素是置信度矩阵C(k)的第j列;||Wj (k)||0是Wj (k)矩阵中非零元素的个数,即确保每一类中只有Q%被选出;
至此得到每一类的选择矩阵Wj (k),通过式(6)将选出的每一类高可信度样本加入字典的训练过程中,即更新用于下一轮训练的测试集样本:
Yj (k+1)=[YtWj (k)],j=1,…,C (6)。
4.根据权利要求1所述的基于子空间投影和字典学习的图像分类方法,其特征在于:所述步骤3)的详细步骤如下:
步骤3.1)假设投影矩阵和n为投影后子空间的维数,则投影到低维子空间上后的训练集样本为测试集样本为
在这个低维平面上利用式(7)训练集样本和挑选之后的测试集样本学习字典,并联合学习投影矩阵Ps和Pt
且||xj||0≤T0,T0为稀疏系数
其中即将两个投影矩阵邻接在一起,同理
其次,即将Ys和Yt作为块矩阵的对角,其余的块矩阵为零,Trace指矩阵的迹,即矩阵主对角线上元素之和;
f1使用利用PsPt将分别将训练集样本和测试集样本投影到低维子空间中,并在低维子空间中联合学习判别字典D和||·||0指矩阵或向量的零范数,即向量中非零元素的个数,||·||F指矩阵的Frobenius范数,即指矩阵各元素平方和开平方;
f2是训练集样本和测试集样本的最大均值距离(MMD),衡量了训练集样本和测试集样本在分布上的差异,分布差异越小,f2的取值越接近0;其中矩阵由式(8)给出:
f3是将样本反投影的过程,期望在投影过程中损失尽可能少的样本信息;
步骤3.2)在求解式(7)的优化问题时引入矩阵 和B,
分别满足:
Ps=(YsAs)T (9)
Pt=(YtAt)T (10)
式(7)的优化方程可转换为:
且||xj||1≤T0,T0为稀疏系数
经过整理式(12)可得到:
需要求的变量有三个,在这里本发明采取固定其他先求其中一个的策略,即先固定求再迭代固定求当被固定时,它们可被看作常数,固式中只剩下为未知量;
首先将初始化,方法是:
先将K进行特征值分解得到:
K=VSVT (14)
其中,V是特征向量矩阵,S是一个对角矩阵,其对角线上的元素是矩阵K的特征值,本方法将As及At的列向量初始化为前n大的特征值所对应的特征向量,即对特征值从大到小排列之后,选取最大的前n个特征值对应的特征向量构成或这里的n是将为后平面的维数;
在得到的初始值后,式中还有未知量和由于
未知量可转换为D和我们利用多次迭代的方法求解,求解D和的具体步骤如下:
将D初始化为随机的单位向量,则式中只剩下为未知量,先将初始化为全0向量再采取逐列更新的方式来计算 为第i类样本系数表示的稀疏系数;
在第h轮中,更新的公式为:
其中,是第h轮更新得到的稀疏系数,是函数Q(x)对的偏导数,
其中,Q(X)由下式给出:
由式(15)解得后,固定式中的同样利用逐列更新的方法来求解D,即每次更新字典D中的一列dj,更新dj的公式为:
其中,Z是经过投影的样本即
至此,求得了D(k)和利用下式求出
则固定和求解
由式(13)(14)可令:
可将式(13)转换为如下的形式:
{G}=arg minGtrace(GTHG) (17)
其中G=[Gs,Gt],式(16)需要利用Stiefel流形来求解,可直接调用现成的软件包,求得G;
至此得到更新的公式为:
可由式(9)和(10)分别解出Ps和Pt,
至此,已由式(16)求得D(k),由式(9)和(10)分别解出Ps和Pt。
5.根据权利要求1所述的基于子空间投影和字典学习的图像分类方法,其特征在于:所述步骤4)的详细步骤如下:
在得到本轮求得的字典D(k)以及投影矩阵Pt (k)后,用下式更新第i个测试集样本的伪标签:
6.根据权利要求1所述的基于子空间投影和字典学习的图像分类方法,其特征在于:所述步骤5)的详细步骤如下:
将所有的测试集样本的伪标签都更新完后,将其与该样本上一轮的伪标签进行比较,若两轮得到的伪标签相同,则称该样本为稳定样本,计算稳定样本的数量;
若该数量大于测试集样本总数的Q%,则样本和字典都趋于稳定,结束迭代并将最后一轮测试集样本的伪标签输出,作为最后的结果;
若稳定样本数量小于测试集样本总数的Q%,则继续迭代直到不满足条件为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810677726.7A CN109117860B (zh) | 2018-06-27 | 2018-06-27 | 一种基于子空间投影和字典学习的图像分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810677726.7A CN109117860B (zh) | 2018-06-27 | 2018-06-27 | 一种基于子空间投影和字典学习的图像分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109117860A true CN109117860A (zh) | 2019-01-01 |
CN109117860B CN109117860B (zh) | 2022-10-25 |
Family
ID=64822387
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810677726.7A Active CN109117860B (zh) | 2018-06-27 | 2018-06-27 | 一种基于子空间投影和字典学习的图像分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109117860B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110046669A (zh) * | 2019-04-22 | 2019-07-23 | 广东石油化工学院 | 基于素描图像的半耦合度量鉴别字典学习的行人检索方法 |
CN110443273A (zh) * | 2019-06-25 | 2019-11-12 | 武汉大学 | 一种用于自然图像跨类识别的对抗零样本学习方法 |
CN110728377A (zh) * | 2019-10-21 | 2020-01-24 | 山东大学 | 一种面向机电设备的智能故障诊断方法及系统 |
CN112183660A (zh) * | 2020-10-19 | 2021-01-05 | 浙江理工大学 | 一种基于类间模糊抑制的判别字典学习纺织图像分类方法 |
CN112183617A (zh) * | 2020-09-25 | 2021-01-05 | 电子科技大学 | 样本与类标签最大相关子空间rcs序列特征提取方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102609681A (zh) * | 2012-01-12 | 2012-07-25 | 北京大学 | 基于字典学习模型的人脸识别方法 |
CN107832786A (zh) * | 2017-10-31 | 2018-03-23 | 济南大学 | 一种基于字典学习的人脸识别分类方法 |
-
2018
- 2018-06-27 CN CN201810677726.7A patent/CN109117860B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102609681A (zh) * | 2012-01-12 | 2012-07-25 | 北京大学 | 基于字典学习模型的人脸识别方法 |
CN107832786A (zh) * | 2017-10-31 | 2018-03-23 | 济南大学 | 一种基于字典学习的人脸识别分类方法 |
Non-Patent Citations (3)
Title |
---|
MENG JIAN等: "Semi-supervised Bi-dictionary Learning Using Smooth Representation-Based Label Propagation", 《2015 INTERNATIONAL CONFERENCE ON CYBER-ENABLED DISTRIBUTED COMPUTING AND KNOWLEDGE DISCOVERY》 * |
孙志武等: "面向分类的增量字典学习算法", 《计算机工程》 * |
张志武: "基于二次学习的半监督字典学习软件缺陷预测", 《模式识别与人工智能》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110046669A (zh) * | 2019-04-22 | 2019-07-23 | 广东石油化工学院 | 基于素描图像的半耦合度量鉴别字典学习的行人检索方法 |
CN110443273A (zh) * | 2019-06-25 | 2019-11-12 | 武汉大学 | 一种用于自然图像跨类识别的对抗零样本学习方法 |
CN110443273B (zh) * | 2019-06-25 | 2022-06-14 | 武汉大学 | 一种用于自然图像跨类识别的对抗零样本学习方法 |
CN110728377A (zh) * | 2019-10-21 | 2020-01-24 | 山东大学 | 一种面向机电设备的智能故障诊断方法及系统 |
CN110728377B (zh) * | 2019-10-21 | 2020-06-09 | 山东大学 | 一种面向机电设备的智能故障诊断方法及系统 |
CN112183617A (zh) * | 2020-09-25 | 2021-01-05 | 电子科技大学 | 样本与类标签最大相关子空间rcs序列特征提取方法 |
CN112183617B (zh) * | 2020-09-25 | 2022-03-29 | 电子科技大学 | 样本与类标签最大相关子空间rcs序列特征提取方法 |
CN112183660A (zh) * | 2020-10-19 | 2021-01-05 | 浙江理工大学 | 一种基于类间模糊抑制的判别字典学习纺织图像分类方法 |
CN112183660B (zh) * | 2020-10-19 | 2024-02-06 | 浙江理工大学 | 一种基于类间模糊抑制的判别字典学习纺织图像分类方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109117860B (zh) | 2022-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109117860A (zh) | 一种基于子空间投影和字典学习的图像分类方法 | |
CN106326288B (zh) | 图像搜索方法及装置 | |
CN110348399B (zh) | 基于原型学习机制和多维残差网络的高光谱智能分类方法 | |
CN108510535A (zh) | 一种基于深度预测和增强子网络的高质量深度估计方法 | |
CN108319957A (zh) | 一种基于超点图的大规模点云语义分割方法 | |
CN105844635B (zh) | 一种基于结构字典的稀疏表示深度图像重建方法 | |
CN109190643A (zh) | 基于卷积神经网络中药识别方法及电子设备 | |
CN111625667A (zh) | 一种基于复杂背景图像的三维模型跨域检索方法及系统 | |
CN109376787B (zh) | 流形学习网络及基于其的计算机视觉图像集分类方法 | |
CN108509925B (zh) | 一种基于视觉词袋模型的行人重识别方法 | |
CN107862680B (zh) | 一种基于相关滤波器的目标跟踪优化方法 | |
CN109711401A (zh) | 一种基于Faster Rcnn的自然场景图像中的文本检测方法 | |
CN109902662A (zh) | 一种行人重识别方法、系统、装置和存储介质 | |
WO2022267388A1 (zh) | 红树林高光谱图像分类方法、装置、电子设备及存储介质 | |
CN112784782B (zh) | 一种基于多视角双注意网络的三维物体识别方法 | |
CN105654122B (zh) | 基于核函数匹配的空间金字塔物体识别方法 | |
CN106529586A (zh) | 基于补充文本特征的图像分类方法 | |
CN113095149A (zh) | 一种基于单张人脸图像的全头部纹理网络结构及生成方法 | |
CN109241813A (zh) | 用于非约束人脸识别的判别稀疏保持嵌入方法 | |
CN106250918A (zh) | 一种基于改进的推土距离的混合高斯模型匹配方法 | |
CN109857895B (zh) | 基于多环路视图卷积神经网络的立体视觉检索方法与系统 | |
CN111178229B (zh) | 一种基于深度学习的静脉显像方法和装置 | |
CN110543845A (zh) | 一种三维人脸的人脸级联回归模型训练方法及重建方法 | |
Filax et al. | Data for Image Recognition Tasks: An Efficient Tool for Fine-Grained Annotations. | |
CN109829377A (zh) | 一种基于深度余弦度量学习的行人重识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |