本发明属图像分析处理
技术领域：
：，涉及一种基于判别共享字典的机织物纹理重构方法，应用于机织物纹理表征和重构领域。
背景技术：
：：基于字典学习的稀疏表征模型是计算机视觉、模式识别和纹理重构领域的研究热点，这主要是由于该模型是通过自适应的学习字典而不是预定义字典(如dct字典)中部分列向量(原子)的线性组合来近似纹理图像，并且该模型在不同领域都取得了优异成果，包括压缩感知、图像去噪、图像修复以及织物纹理重构等，这表明基于稀疏编码的学习字典能自发地“学习”到数据本身的规律，与主成分分解及其衍生算法不同，稀疏表征模型并没有对基向量施加正交约束，因此该模型能更加灵活地适应真实世界中的数据分布。具体地，对于含有n个样本的数据集表示其中一个样本，稀疏表征模型假设y＝dx+e，其中是一个含有k个列向量(列向量又被称为原子)的超完备字典，e作为噪声项可提高模型的鲁棒性，在系数更新阶段，稀疏系数通过求解追踪问题获得：其中||x||0表示矩阵x中非零元素的个数。解决追踪问题的传统算法是贪婪算法(如匹配追踪(matchingpursuit,mp)、正交匹配追踪(orthogonalmatchingpursuit,omp))或者凸松弛算法(如基追踪(basispursuit,bp))。字典更新阶段，字典d可通过l2范数优化问题求解：常用的解决以上字典学习问题的算法为k-svd、mod、doublesparsity和在线字典学习(onlinedictionarylearning,odl)等等。织物纹理分析方面，文献1(unsupervisedfabricdefectsegmentationusinglocalpatchapproximation.thejournalofthetextileinstitute2016；107:800-809)首次采用字典学习方法表征织物纹理，并将其应用于织物疵点检测，取得了较高准确率，召回率也在可接受范围之内，表明模型具有优异的鲁棒性能，但其模型针对某一特定纹理的重复表征不能保证所学习的字典和稀疏系数的唯一性；文献2(基于子窗口字典学习的机织物纹理表征及应用[j].东华大学学报(自然科学版),2019,45(03):375-380)提出的基于子窗口字典学习模型将整幅图像划分为多个子窗口样本，并将子窗口样本展成列向量，所有的列向量联合组成灰度数据矩阵，并对字典个数和子窗口大小进行优化，结果表明，应用学习得到的字典能近似重构机织物纹理样本图像；文献3(fabrictexturerepresentationusingthestablelearneddiscretecosinetransformdictionary.textileresearchjournal2019；89:294-310)提出的稳定典学习模型能够在多次重复实验中学习到一致的字典和稀疏系数，另外，该文献还分析了织物组织结构参数对纹理表征结果的影响，但其模型只能应用于某一种机织物纹理的表征，未能实现多种类别的机织物纹理同时进行分析。因此，尽管上述已有文献所涉及的织物纹理分析方法有一定的适用性，但都只能针对某一类机织物纹理有效，不能对多种类别的机织物纹理同时分析处理，具有明显的局限和不足，因而在实用性和广泛性方面大打折扣。技术实现要素：本发明提供一种基于判别共享字典的机织物纹理重构方法，目的是解决现有技术中的机织物纹理分析方法只能针对某一类机织物纹理有效，不能对多种类别的机织物纹理同时分析处理的问题。本发明提供的一种基于判别共享字典的机织物纹理重构方法，能对多种类别的机织物纹理同时进行分析，分别提取他们的判别特征和共享特征，可以应用于一般情况下的纹理分析和纹理重构问题；特别是针对机织物纹理的重构，能够达到较高的图像相似度。为了达到上述目的，本发明采用的技术方案为：一种基于判别共享字典的机织物纹理重构方法，包括如下步骤：步骤1：对于含有c类样本的机织物纹理图像集，先对其中的第c(c＝1，2，…，c)类机织物纹理样本图像ic进行直方图均衡化增强纹理特征，再将样本图像ic划分成nc个w×w的图像块，并映射每个图像块至w2×1的一维列向量，将nc个一维列向量进行合并得到矩阵则包含c类样本的数据集其中和表示矩阵属于实数空间并指明了维度，w为划分得到的图像块的尺寸，c为机织物纹理图像集的样本种类数量，c为机织物纹理图像样本类别；步骤2：建立判别共享字典学习模型(dsdl)；所述判别共享字典学习模型包括fisher判别项和系数判别项其中，为完备字典，为完备的稀疏系数矩阵；步骤3：设置结构化字典d的原子个数，共享字典d0的原子个数，以及参数η1，η2，η3的值，对于包含c类机织物纹理样本的数据集y＝[y1，…，y0，…，yc]，首先使用在线字典学习算法(onlinedictionarylearning,odl)计算每一类样本的初始字典d1，…，d0，…，dc，通过随机选取包含c类样本的数据集y的列向量得到初始的共享字典d0，则初始的完备字典然后使用匹配追踪算法计算初始的稀疏系数矩阵在此基础上输入数据集y对dsdl模型进行优化求解，得到适应于该数据集的结构化字典d、共享字典d0和完备的稀疏系数矩阵步骤4：对新的机织物纹理测试图像ii，先按步骤1的方法对其进行预处理得到矩阵yi，再基于步骤3中得到的结构化字典d和共享字典d0利用匹配追踪算法计算其对应的稀疏系数步骤5：利用结构化字典d、共享字典d0和稀疏系数进行重构得到y′i；步骤6：提取矩阵y′i中的每一列重新组合成图像块，再将组合成的图像块拼接成重构图像i′i；所述步骤2中建立的判别共享字典学习模型为：其中dc为第c类机织物纹理样本图像对应的类别子字典且每一列都经过归一化处理，d＝[d1，…，d0，…dc]为结构化字典；d0为共享字典且每一列都经过归一化处理，为完备字典；为完备的稀疏系数矩阵；x为包含c类样本的数据集y基于结构化字典d的稀疏系数；xc为第c类样本yc基于结构化字典d的稀疏系数；为第c类样本基于类别子字典dj的稀疏系数；x0为包含c类样本的数据集y基于共享字典d0的稀疏系数；mc，m和m0分别为xc，x和x0的均值向量构成的均值矩阵；η1，η2和η3为预先定义好的参数，控制模型中各约束项所占的比重；表示矩阵a的f范数，即矩阵a中各元素的平方和；||a||0表示矩阵a的核范数，即矩阵a的奇异值的和；||a||1表示矩阵a的l1范数，即矩阵a中各元素绝对值之和。作为优选的技术方案：如上所述的一种基于判别共享字典的机织物纹理重构方法，步骤1中所述一维列向量由划分所得机织物图像块所有列与列之间首尾相连组成。如上所述的一种基于判别共享字典的机织物纹理重构方法，步骤1中所述将n0个一维列向量进行合并得到矩阵具体为令这nc个一维列向量分别为矩阵如上所述的一种基于判别共享字典的机织物纹理重构方法，由于机织物是由经纬纱线垂直交织而成，其纹理特征具有一定的正交特性，因此用于表征其判别特征的每一个类别子字典d0为正交字典，即d0满足：其中表示dc的第i个原子，表示dc的第j个原子。如上所述的一种基于判别共享字典的机织物纹理重构方法，在所述判别共享字典学习模型中机织物纹理被分为共享特征和判别特征，共享特征的线性子空间不与判别特征的线性子空间产生交叉。如上所述的一种基于判别共享字典的机织物纹理重构方法，||d0||x为对d0施加的低秩约束。如上所述的一种基于判别共享字典的机织物纹理重构方法，步骤3中对所构建的dsdl模型进行优化求解的方法为：针对dsdl模型中涉及到的三个变量d.d0和在固定任意两个变量的情况下对第三个变量进行优化，并交替迭代，当每次迭代过程中目标函数值都收敛，或者达到预设的最大迭代次数，则停止迭代。如上所述的一种基于判别共享字典的机织物纹理重构方法，所述机织物纹理样本图像和机织物纹理测试图像均为位深度为8位以上的灰度图像。有益效果：本发明的一种基于判别共享字典的机织物纹理重构方法，能对多种类别的机织物纹理同时进行分析，分别提取他们的判别特征和共享特征，可以应用于一般情况下的纹理分析和纹理重构问题；特别是针对机织物纹理的重构，能够达到较高的图像相似度。附图说明图1是本发明的流程图；图2是本发明所提出模型的示意图；图3是机织物纹理样本示意图；图4是完备字典可视化图；图5是完备的稀疏系数矩阵可视化图；图6是不同类别子字典对不同类别样本的重构图像；图7是基于不同类别子字典的峰值信噪比指标折线图。具体实施方式下面结合具体实施方式，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解，在阅读了本发明讲授的内容之后，本领域技术人员可以对本发明作各种改动或修改，这些等价形式同样落于本申请所附权利要求书所限定的范围。一种基于判别共享字典的机织物纹理重构方法，包括如下步骤：步骤1：对于含有c类样本的机织物纹理图像集，先对其中的的第c(c＝1，2，…，c)类机织物纹理样本图像ic(位深度为8位以上的灰度图像)进行直方图均衡化增强纹理特征，再将样本图像iv划分成nv个w×w的图像块，并映射每个图像块至w2×1一维列向量，将nv个一维列向量进行合并得到矩阵则包含c类样本的数据集其中和表示矩阵属于实数空间并指明了维度，w为划分得到的图像块的尺寸，c为机织物纹理图像集的样本种类数量，c为机织物纹理图像样本类别，所述一维列向量由划分所得机织物图像块所有列与列之间首尾相连组成，所述将nc个一维列向量进行合并得到矩阵具体为令这nc个一维列向量分别为矩阵步骤2：建立判别共享字典学习模型dsdl，所述判别共享字典学习模型包括fisher判别项和系数判别项判别共享字典学习模型为：其中dc为第c类机织物纹理样本图像对应的类别子字典且每一列都经过归一化处理，d＝[d1，…，d0，…dc]为结构化字典；d0为共享字典且每一列都经过归一化处理，为完备字典；为完备的稀疏系数矩阵；x为包含c类样本的数据集y基于结构化字典d的稀疏系数；xc为第c类样本yc基于结构化字典d的稀疏系数；为第c类样本基于类别子字典dj的稀疏系数；x0为包含c类样本的数据集y基于共享字典d0的稀疏系数；mc，m和m0分别为xc，x和x0的均值向量构成的均值矩阵；η1，η2和η3为预先定义好的参数，控制模型中各约束项所占的比重；表示矩阵a的f范数，即矩阵a中各元素的平方和；||a||x表示矩阵a的核范数，即矩阵a的奇异值的和；||a||1表示矩阵a的l1范数，即矩阵a中各元素绝对值之和；dc为正交字典，即dc满足：其中表示dc的第i个原子，表示dc的第j个原子；||d0||x为对d0施加的低秩约束；在上述判别共享字典学习模型中机织物纹理被分为共享特征和判别特征，共享特征的线性子空间不与判别特征的线性子空间产生交叉；步骤3：设置结构化字典d的原子个数，共享字典d0的原子个数，以及参数η1，η2，η3的值，对于包含c类机织物纹理样本的数据集y＝[y1，…，y0，…，yc]，首先使用在线字典学习算法(onlinedictionarylearning,odl)计算每一类样本的初始字典d1，…，d0，…，dc，通过随机选取包含c类样本的数据集y的列向量得到初始的共享字典d0，则初始的完备字典然后使用匹配追踪算法计算初始的稀疏系数矩阵在此基础上输入数据集y对dsdl模型进行优化求解，得到适应于该数据集的结构化字典d、共享字典d0和完备的稀疏系数矩阵如图2所示；其中，对dsdl模型进行优化求解的方法为：针对dsdl模型中涉及到的三个变量d、d0和在固定任意两个变量的情况下对第三个变量进行优化，并交替迭代，当每次迭代过程中目标函数值都收敛，或者达到预设的最大迭代次数，则停止迭代；步骤4：对新的机织物纹理测试图像ii(位深度为8位以上的灰度图像)，先按步骤1的方法对其进行预处理得到矩阵yi，再基于步骤3中得到的结构化字典d和共享字典d0利用匹配追踪算法计算其对应的稀疏系数步骤5：利用结构化字典d、共享字典d0和稀疏系数进行重构得到y′i；步骤6：提取矩阵y′i中的每一列重新组合成图像块，再将组合成的图像块拼接成重构图像i′i。下面结合附图，对本发明的实施例进行具体描述。实施例只用于对本发明进行进一步说明，不能理解为对本发明保护范围的限定，该领域的技术人员可根据上述发明的内容对本发明做出一些非本质的改进和调整；图1为本发明整体流程图。本发明的方法具体运行的硬件和编程语言并不受限制，可用任何语言编写完成，为此其他工作模式不再赘述；本实施例搭建于一台具有inteli7cpu和8g内存的基于ubuntu系统的计算机，并用matlab编程语言实现了基于判别共享字典的字典学习程序，实现了本发明的方法；采集平纹、斜纹、缎纹和蜂巢组织4种不同类别的位深度为8位的机织物纹理图像各一块，记为i0(c＝1，2，3，4)，如图3所示，其尺寸为512×512像素；分别对每一个机织物纹理图像作直方图均衡化处理，并将其划分成64×64的图像块，对每一个图像块按照所有列向量首尾相连的方式映射为一维列向量，则一共可得到64个长度为4096像素的一维列向量；将这64个一维列向量进行合并得到矩阵则包含所有4类图像的数据集建立判别共享字典学习模型dsdl：本实施例中c＝4，类别子字典d0的原子个数为50，共享字典d0的原子个数为10，参数η1＝0.001，η2＝0.1，η3＝0.01；对于包含4类机织物纹理样本的数据集y＝[y1，y2，y3，y4]，首先使用在线字典学习算法(onlinedictionarylearning,odl)计算每一类样本的初始字典d1，d2，d3d4，通过随机选取矩阵y的列向量得到初始的共享字典d0，则初始的完备字典然后使用匹配追踪算法计算初始的稀疏系数矩阵在此基础上针对dsdl模型中涉及到的三个变量d、d0和在固定任意两个变量的情况下对第三个变量优化，并交替迭代20次，得到的完备字典如图4所示，前四行为四个类别子字典d0，最后一行是共享字典d0，完备的稀疏系数矩阵如图5所示；对所有4种类别的机织物纹理图像，分别利用不同类别子字典、共享字典利用匹配追踪算法计算对应的稀疏系数；利用完备字典中的不同类别子字典、共享字典和对应的稀疏系数对这4种类别的机织物纹理样本图像分别进行重构，提取重构矩阵中的每一列重新组合成图像块，再将组合成的图像块拼接成重构图像，重构图像如图6所示，可以看出对应类别子字典只对该类机织物纹理图像有较好的重构效果而对其他类别机织物纹理图像重构效果较差，这说明该模型实现了对多种类别机织物纹理的同时表征，并且体现了该模型的判别性；计算不同重构组合下得到的重构图像与原图像之间的峰值信噪比(psnr)指标，结果如图7所示，图7定量地揭示了模型的判别性能。当前第1页12当前第1页12