本发明公开了一种基于压缩感知的坯布织物疵点检测方法，首先，将采集到的大幅面待测坯布织物图像分块，并统计各个图像块的灰度直方图Ht，并采用观测矩阵Φ获取其M×1维的压缩测量值yt；其次，以无疵点参照图像块Xref的正常灰度分布区间为基准，根据Ht的灰度分布信息，对yt进行量化，以消除光照渐变的影响，得到灰度量化后的压缩测量值yt_det；最后，判断待测的图像块Xt中是否含有疵点，并采用基追踪去噪算法对疵点信息进行恢复与定位。采用本发明的方法进行坯布织物疵点检测，可使大幅面坯布疵点在线自动检测所存储和传输的数据量压缩至原来的50％以上，疵点检出率达90％以上；可有效地提高系统运行的效率，降低系统的成本。
1.一种基于压缩感知的坯布织物疵点检测方法，其特征在于，具体过程为：步骤1：从采集的坯布织物图像中截取一幅无疵点图像，将其划分为互不重叠的p×q大小的图像块Xi，其中，i＝1,2,…n并统计其正常灰度分布区间，计算出训练图像块Xj灰度量化后的压缩测量值，随后计算训练图像块Xj灰度量化后的压缩测量值的均值ymean；
训练图像块Xj灰度量化后的压缩测量值的均值ymean的计算过程为：步骤1.1：随机地选定一个图像块，记作Xref，并统计Xref的灰度直方图Href；
步骤1.2：选取步骤1.1所得Href中小于等于总像素个数2％的灰度级像素个数作为阈值thresh，找出Href中灰度级像素个数大于thresh的所有像素点，并按其灰度级从小到大的顺序排列成一个集合，记为S>thresh，即S>thresh＝{j|Href(j)>thresh}；
步骤1.3：将步骤1.2所得的S>thresh中彼此相邻像素间灰度值两两相减，得集合SUB；将SUB中小于等于2对应的像素位置置为1，大于2的像素位置置为0；
步骤1.4：对步骤1.3中所得SUB中连续出现1的位置区域界定为S>thresh的子区间，若仅得到1个子区间，则将此子区间对应的像素灰度分布范围作为Xref的正常灰度分布区间；否则，选取对应灰度直方图中所含的灰度级像素点个数最多的子区间，将其对应的像素灰度分布范围作为参照图像块Xref的正常灰度分布区间，记为[g1,g2]；
步骤1.5：计算步骤1.2所得S>thresh中所有像素的灰度均值，记为μnum，即
其中，nj表示S>thresh中j灰度级像素的个数；
步骤1.6：找出步骤1.1所得Href中大于步骤1.5所得μnum的所有灰度级像素，并按灰度级从小到大的顺序排列成一个集合，记为 即 将所得 中彼此相邻像素间灰度值两两相减，得到集合SUMm；将SUMm中小于等于2对应的像素位置置为
步骤1.7：将步骤1.6所得SUMm中连续出现1的位置区域界定为 的子区间，若仅得到1个子区间，则将此区间对应的像素灰度分布范围作为Xref的主导像素的灰度分布区间，否则，选取对应灰度直方图中所含的灰度级像素点个数最多的子区间，将其对应的像素灰度分布范围作为Xref的主导像素的灰度分布区间，记为[gm1,gm2]；
步骤1.8：从Xi中，随机地选取C幅图像块Xj作为训练图像块，其中，i＝1,2,…n，j＝1,
2,...C，C的取值不宜低于10幅；分别统计训练图像块Xj的灰度直方图Hj；
步骤1.9：采用观测矩阵Φ，按照公式(1)计算训练图像块Xj的M×1维的压缩测量值yj：yj＝Φxj (1)
其中，xj是由二维图像信号Xj按列优先的顺序排列而成的N×1维的向量，yj为xj在压缩域中M×1维的压缩测量值；
步骤1.10：利用步骤1.1～1.4的方法计算选定的C幅训练图像块Xj的灰度分布区间[gj1,gj2]；以参照图像块Xref的灰度分布区间[g1,g2]为基准，根据 训练图像块的压缩测量值yj，按照公式(2)计算训练图像块Xj在[gj1,gj2]区间灰度量化后的压缩测量值yj_norm：其中，φn表示观测矩阵Φ的第n列；
步骤1.11：按照公式(3)计算训练图像块Xj灰度量化后的压缩测量值yj_norm的均值ymean：步骤2：采用极大似然估计法，计算训练图像块Xj灰度量化后的压缩测量值的正态分布参数μnorm和σnorm；
步骤3：将待测坯布织物图像划分成互不重叠的p×q大小的图像块Xt，其中，，统计其灰度直方图Ht，计算待测图像块Xt灰度量化后的压缩测量值yt_det；
步骤4：将步骤1、2和步骤3所得的ymean、μnorm、σnorm和yt_det的值带入下式，判断待测的图像块Xt中是否含有疵点，式中，ymean为训练图像块Xj灰度量化后的压缩测量值的均值；
μnorm和σnorm为训练图像块Xj灰度量化后的压缩测量值的正态分布参数；
yt_det为待测图像块Xt灰度量化后的压缩测量值；
c为一个正的常数，其取值与检测的精确度有关，根据正态分布数 理统计理论中的“3σ”原则，c的取值范围为[2,3]；
如果待测图像块Xt灰度量化后的压缩测量值yt_det与无疵点的训练图像块学习得到的参数ymean、μnorm和σnorm满足条件：则认为待测图像块Xt含有疵点，即置为1，否则，不含
步骤5：若待测的图像块Xt中无疵点，返回，处理其他待测图像块；若含有疵点，则进行疵点定位。
2.根据权利要求1所述的一种基于压缩感知的坯布织物疵点检测方法，其特征在于，步骤3中，待测图像块Xt灰度量化后的压缩测量值yt_det的计算过程为：步骤3.1：根据待测图像块Xt的灰度直方图Ht，参照步骤1.2～步骤1.4的方法计算待测图像块Xt的灰度分布区间，记为[gt1,gt2]；
步骤3.2：以参照图像块Xref的正常灰度分布区间[g1,g2]为基准，参照公式(1)计算待测图像块Xt的压缩测量值yt后，根据待测图像块Xt的压缩测量值yt，参照公式(2)计算Xt在[gt1,gt2]区间量化后的压缩测量值yt_det。
3.根据权利要求1所述的一种基于压缩感知的坯布织物疵点检测方法，其特征在于，步骤5中，疵点定位的过程为：步骤5.1：根据待测图像块Xt的灰度直方图Ht，利用步骤1.5～1.7的方法计算Xt主导像素的灰度分布区间，记为[gt_m1,gt_m2]；以参照图像块Xref主导像素的灰度分布区间[gm1,gm2]为基准，参照公式(1)计算待测图像块Xt的M×1维的压缩测量值yt后，根据待测图像块Xt的压缩测量值yt，参照公式(2)计算Xt在[gt_m1,gt_m2]区间灰度量化后的压缩测量值yt_mdet；
步骤5.2：将步骤5.1所得结果带入公式(4)中滤除yt_mdet中无疵点信息后，得待测图像块Xt中的疵点信息yt_def：yt_def＝yt_mdet-ymean (4)
步骤5.3：将步骤5.2所得yt_def和观测矩阵Φ带入公式(5)中，对获取的疵点信息yt_def进行恢复，得到N×1维的信号其中，λ是一个正的常数，取值范围为[2,4]，其取值需考虑信息的恢复精度与计算速度间的权衡；
步骤5.4：将步骤5.3所得的 按列优先的顺序排列成与原始图像块同维数的p×q维矩阵Xt_rec，利用公式(6)对Xt_rec进行3×3窗口的中值滤波，得到Xt_f：Xt_f(i,j)＝median{Xt_rec(i-3,j-3),…,Xt_rec(i,j),…Xt_rec(i+3,j+3)} (6)步骤5.5：利用公式(7)对步骤5.4所得Xt_f进行二值化处理，得疵点的位置及轮廓信息Xt_b，即完成疵点定位：
一种基于压缩感知的坯布织物疵点检测方法
技术领域
[0001] 本发明属于数字图像目标检测技术领域，具体涉及一种基于压缩感知的坯布织物疵点检测方法。
背景技术
[0002] 机器视觉和数字图像处理技术的飞速发展为基于机器视觉的自动化疵点检测奠定了重要的基础。尽管织物疵点检测的自动化研究已有三十余年的发展历史，然而，面对国际纺织市场竞争的日趋激烈，企业对坯布疵点检测也提出了新的需求：针对大幅面坯布的疵点检测，要求在检出坯布疵点的同时，能降低数据存储和传输量，提高系统运行的效率，降低系统的成本。而现有的技术方法难以满足这一新的需求。因此，提高大幅面坯布疵点的检出率和运行效率、降低系统的成本已成为坯布疵点自动检测系统中不可缺少的一个核心技术。
[0003] 压缩感知(Compressed Sensing，CS)理论是Donoho等人于2006年正式提出的一种新型的“全息”采样理论。该理论已证明：当信号是稀疏的或在某种变换下可以稀疏表示时，采用满足约束等距性的观测矩阵提取少量的的压缩测量值即可完美地保存原始信号的信息；在信号重构阶段，通过相应的重构算法即可从少量的压缩测量值中恢复出原始信号，这为压缩域的目标检测提供了理论依据。
发明内容
[0004] 本发明的目的是提供一种基于压缩感知的坯布织物疵点检测方法，能够满足大幅面坯布疵点的检出率，解决现有织物疵点检测技术中因大量的冗余信息导致系统运行效率低、成本高的问题。
[0005] 本发明所采用的技术方案是，一种基于压缩感知的坯布织物疵点检测方法，具体过程为：
[0006] 步骤1：从采集的坯布织物图像中截取一幅无疵点图像，将其划分为互不重叠的p×q大小的图像块Xi(i＝1,2,…n)，并统计其正常灰度分布区间，计算出训练图像块Xj灰度量化后的压缩测量值yj，随后计算其压缩测量值的均值ymean；
[0007] 步骤2：采用极大似然估计法，计算训练图像块Xj的压缩测量值的正态分布参数μnorm和σnorm；
[0008] 步骤3：将待测坯布织物图像划分成互不重叠的p×q大小的图像块Xt(t＝1,2….n)，统计其灰度直方图Ht，计算待测图像块Xt灰度量化后的压缩测量值yt_det；
[0009] 步骤4：将步骤1、2和步骤3所得的ymean、μnorm、σnorm和yt_det的值带入下式，判断待测的图像块Xt中是否含有疵点，
[0010]
[0011] 式中，ymean为计算训练图像块Xj灰度量化后的压缩测量值的均值；
[0012] μnorm和σnorm为训练图像块Xj灰度量化后的压缩测量值的正态分布参数；
[0013] yt_det为待测图像块Xt灰度量化后的压缩测量值；
[0014] c为一个正的常数，其取值与检测的精确度有关，根据正态分布数理统计理论中的“3σ”原则，c的取值范围为[2,3]；
[0015] 如果待测图像块Xt的压缩测量值yt_det与无疵点的训练图像块学习得到的参数ymean、μnorm和σnorm满足条件： 则认为待测图像块Xt含有疵点，即置为1，否则，不含疵点，置为0；
[0016] 步骤5：若待测的图像块Xt中无疵点，返回，处理其他待测图像块；若含有疵点，则进行疵点定位。
[0017] 本发明的特征还在于，
[0018] 步骤1中，训练图像块Xj灰度量化后的压缩测量值的均值ymean的计算过程为：
[0019] 步骤1.1：随机地选定一个图像块，记作Xref，并统计Xref的灰度直方图Href；
[0020] 步骤1.2：选取步骤1.1所得Href中小于等于总像素个数1％-2％的灰度级像素个数作为阈值thresh，找出Href中灰度级像素个数大于thresh的所有像素点，并按其灰度级从小到大的顺序排列成一个集合，记为S＞thresh，即S＞thresh＝{j|Href(j)＞thresh}；
[0021] 步骤1.3：将步骤1.2所得的S＞thresh中彼此相邻像素间灰度值两两相减，得集合SUB；将SUB中小于等于2对应的像素位置置为1，大于2的像素位置置为0；
[0022] 步骤1.4：对步骤1.3中所得SUB中连续出现1的位置区域界定为S＞thresh的子区间，若仅得到1个子区间，则将此子区间对应的像素灰度分布范围作为Xref的正常灰度分布区间；否则，选取对应灰度直方图中所含的灰度级像素点个数最多的子区间，将其对应的像素灰度分布范围作为参照图像块Xref的灰度分布区间，记为[g1,g2]；
[0023] 步骤1.5：计算步骤1.2所得S＞thresh中所有像素的灰度均值，记为
[0024] μnum，即
[0025] 其中，nj表示S＞thresh中j灰度级像素的个数；
[0026] 步骤1.6：找出步骤1.1所得Href中大于步骤1.5所得μnum的所有灰度级像素，并按灰度级从小到大的顺序排列成一个集合，记为 即 将所得中彼此相邻像素间灰度值两两相减，得到集合SUMm；将SUMm中小于等于2对应的像素位置置为1，大于2的像素位置置为0；
[0027] 步骤1.7：将步骤1.6所得SUMm中连续出现1的位置区域界定为 的子区间，若仅得到1个子区间，则将此区间对应的像素灰度分布范围作为Xref的主导像素的灰度分布区间，否则，选取对应灰度直方图中所含的灰度级像素点个数最多的子区间，将其对应的像素灰度分布范围作为Xref的主导像素的灰度分布区间，记为[gm1,gm2]；
[0028] 步骤1.8：从Xi(i＝1，2，…n)中，随机地选取C幅图像块Xj(j＝1,2,...C)作为训练图像块，其中，C的取值不宜低于10幅；分别统计训练图像块Xj的灰度直方图Hj；
[0029] 步骤1.9：采用观测矩阵Φ，按照公式(1)计算训练图像块Xj的M×1维的压缩测量值yj：
[0030] yj＝Φxj (1)
[0031] 其中，xj是由二维图像信号Xj按列优先的顺序排列而成的N×1维的向量，yj为xj在压缩域中M×1维的压缩测量值；
[0032] 步骤1.10：利用步骤1.1～1.4的方法计算选定的C幅训练图像块Xj的灰度分布区间[gj1,gj2]；以参照图像块Xref的灰度分布区间[g1,g2]为基准，根据训练图像块的压缩测量值yj，按照公式(2)计算训练图像块Xj在[gj1,gj2]区间量化后的压缩测量值yj_norm：
[0033]
[0034] 其中，φn表示观测矩阵Φ的第n列；
[0035] 步骤1.11：按照公式(3)计算训练图像块Xj灰度量化后的压缩测量值yj_norm的均值ymean：
[0036]
[0037] 步骤3中，待测图像块Xt灰度量化后的压缩测量值yt_det的计算过程为：
[0038] 步骤3.1：根据待测图像块Xt的灰度直方图Ht，参照步骤1.2～步骤1.4的方法计算待测图像块Xt的灰度分布区间，记为[gt1,gt2]；
[0039] 步骤3.2：以参照图像块Xref的正常灰度分布区间[g1,g2]为基准，根据待测图像块Xt的压缩测量值yt，参照公式(2)计算Xt在[gt1,gt2]区间量化后的压缩测量值yt_det。
[0040] 步骤5中，疵点定位的过程为：
[0041] 步骤5.1：根据待测图像块Xt的灰度直方图Ht，利用步骤1.5～1.7的方法计算Xt主导像素的灰度分布区间，记为[gt_m1,gt_m2]；以参照图像块Xref主导像素的灰度分布区间[gm1,gm2]为基准，计算Xt在[gt_m1,gt_m2]区间量化后的压缩测量值yt_mdet；
[0042] 步骤5.2：将步骤5.1所得结果带入公式(4)中滤除yt_mdet中无疵点信息后，得待测图像块Xt中的疵点信息yt_def：
[0043] yt_def＝yt_mdet-ymean (4)
[0044] 步骤5.3：将步骤5.2所得yt_def和观测矩阵Φ带入公式(5)中，对获取的疵点信息yt_def进行恢复，得到N×1维的信号
[0045]
[0046] 其中，λ是一个正的常数，取值范围为[2,4]，其取值需考虑信息的恢复精度与计算速度间的权衡；
[0047] 步骤5.4：将步骤5.3所得的 按列优先的顺序排列成与原始图像块同维数的p×q维矩阵Xt_rec，利用公式(6)对Xt_rec进行3×3窗口的中值滤波，得到Xt_f：
[0048] Xt_f(i,j)＝median{Xt_rec(i-3,j-3),…,Xt_rec(i,j),…Xt_rec(i+3,j+3)} (6)[0049] 步骤5.5：利用公式(7)对步骤5.4所得Xt_f进行二值化处理，得疵点的位置及轮廓信息Xt_b，即完成疵点定位：
[0050]
[0051] 本发明的有益效果是，
[0052] 1.本发明利用坯布疵点的空域稀疏性，根据CS理论对坯布图像信息进行压缩，能够将疵点检测系统的信息量降至原来的50％-70％之间；
[0053] 2.本发明根据压缩测量值判别疵点的有无，直接滤除压缩测量值中无疵点信息，仅对含有疵点的信息进行恢复与定位，在保证误检率低于10％的前提下，能够达到90％以上的疵点检出率；
[0054] 3.本发明利用坯布织物图像灰度分布信息，通过自适应量化方法，可有效地克服不可避免的光照渐变以及疵点区与背景间灰度分布交界处的模糊性带来的不利影响。
附图说明
[0055] 图1是利用本发明方法在大幅面坯布织物图像上截取的无疵点图像块及其量化结果图；
[0056] 图2是利用本发明方法在大幅面坯布织物图像上截取的疵点图像块及其量化结果图；
[0057] 图3是利用本发明方法在大幅面坯布织物图像上截取的无疵点图像块和疵点图像块分别对应的 值的分布；
[0058] 图4是利用本发明方法在大幅面坯布织物图像上截取的含有疵点信息的待测图像块中疵点的恢复及其定位结果图。
具体实施方式
[0059] 下面结合附图和具体实施方式对本发明进行详细说明。
[0060] 本发明方法的原理是：利用坯布中疵点具有稀疏性的特点，首先，将采集到的大幅面待测坯布织物图像划分为互不重叠的p×q大小的图像块Xt，分别统计各个图像块的灰度直方图Ht，并采用观测矩阵Φ获取其M×1维的压缩测量值yt；接着，以无疵点参照图像块Xref的正常灰度分布区间为基准，根据Ht的灰度分布信息，对yt进行量化，以消除光照渐变的影响，得到灰度量化后的压缩测量值yt_det；然后根据对无疵点训练图像块学习得到的参数，判断待测的图像块Xt中是否含有疵点，若无疵点，则返回，重新开始新的待测坯布织物图像的疵点检测，否则，以无疵点参照图像块Xref的主导像素灰度分布区间为基准，根据Ht的主导像素灰度分布信息，对yt进行量化，以消除疵点与无疵点间灰度分布交界处的模糊性，得到灰度量化后的压缩测量值yt_mdet；最后，从yt_mdet中提取疵点信息得到yt_def，根据yt_def和观测矩阵Φ，采用基追踪去噪算法对疵点信息进行恢复与定位。
[0061] 本发明的坯布织物疵点检测方法的实施可以分为三个阶段：参数学习阶段、疵点判别阶段和疵点定位阶段，具体按照以下步骤实施：
[0062] 根据压缩感知理论，选用高斯随机矩阵作为观测矩阵。设：Φ表示M×N维的高斯随机矩阵，其中，N＝p×q，M≥O(Kdeflog(N/Kdef))，Kdef是依据经验估计的待测织物图像块中最大的疵点所占的像素点个数； 服从均值为μnorm、方差为 的正态分布。
[0063] 第一阶段：参数学习阶段
[0064] 步骤1：从采集的坯布织物图像中截取一幅无疵点图像，并将其划分为互不重叠的p×q大小的图像块Xi(i＝1,2,…n)，为了保证图像块中疵点的稀疏性，其中p×q≥5×Kdef；
[0065] 步骤2：选定无疵点的参照图像块xref，计算其正常灰度分布区间[g1,g2]和主导像素的灰度分布区间[gm1,gm2]，为消除光照渐变的影响以及疵点区与背景间灰度分布交界处的模糊性提供参照，具体流程如下：
[0066] 1)从Xi(i＝1,2,…n)中，随机地选定一个无疵点图像块作为参照图像块，记作Xref，并统计Xref的灰度直方图Href；
[0067] 2)滤除灰度直方图Href中各灰度级像素个数占总像素个数1％-2％的像素点，然后，从中计算参照图像块Xref的正常灰度分布区间[g1,g2]，具体流程如下：
[0068] a.选取Href中小于等于总像素个数1％-2％的灰度级像素个数作为阈值thresh；
[0069] b.找出Href中灰度级像素个数大于thresh的所有像素点，并按其灰度级从小到大的顺序排列成一个集合，记为S＞thresh，即S＞thresh＝{j|Href(j)＞thresh}；
[0070] c.令S＞thresh中间灰度值两两相减(即后一个像素灰度值减前一个像素灰度值)，得到一个结果集合，记为SUB；
[0071] d.将SUB中小于等于2对应的像素位置置为1，大于2的像素位置置为0；
[0072] e.对连续出现1的位置区域界定为S＞thresh的子区间，若仅得到1个子区间，则将此子区间对应的像素灰度分布范围作为Xref的正常灰度分布区间，否则，选取对应灰度直方图中所含的灰度级像素点个数最多的子区间，将其对应的像素灰度分布范围作为Xref的正常灰度分布区间，记为[g1,g2]。
[0073] 3)计算参照图像块Xref主导像素的灰度分布区间[gm1,gm2]，具体流程如下：
[0074] a.计算S＞thresh中所有像素的灰度均值，记为μnum，即 其中，nj表示S＞thresh中j灰度级像素的个数；
[0075] b.找出Href中大于μnum的所有灰度级像素，并按灰度级从小到大的顺序排列成一个集合，记为 即
[0076] c.令 中彼此相邻像素间灰度值两两相减(即后一个像素灰度值减前一个像素灰度值)，得到一个结果集合，记为SUMm；
[0077] d.将SUMm中小于等于2对应的像素位置置为1，大于2的像素位置置为0；
[0078] e.对连续出现1的位置区域界定为 的子区间，若仅得到1个子区间，则将此区间对应的像素灰度分布范围作为Xref的主导像素的灰度分布区间，否则，选取对应灰度直方图中所含的灰度级像素点个数最多的子区间，将其对应的像素灰度分布范围作为Xref的主导像素的灰度分布区间，记为[gm1,gm2]。
[0079] 步骤3：选取无疵点的训练图像块Xj(j＝1,2,...C)，分别统计其灰度直方图Hj，并采用观测矩阵Φ获取其M×1维的压缩测量值yj，具体流程如下：
[0080] 1)从Xi(i＝1，2，…n)中，随机地选取C幅无疵点图像块Xj(j＝1,2,...C)作为训练图像块，其中，C的取值不宜低于10幅；
[0081] 2)分别统计训练图像块Xj的灰度直方图Hj；
[0082] 3)采用观测矩阵Φ，按照公式(1)计算训练图像块Xj的M×1维的压缩测量值yj：
[0083] yj＝Φxj (1)
[0084] 其中，xj是由二维图像信号Xj按列优先的顺序排列而成的N×1维的向量，yj为xj在压缩域中M×1维的压缩测量值；
[0085] 步骤4：分别计算训练图像块Xj(j＝1,2,...C)的正常灰度分布区间[gj1,gj2]，并以参照图像块Xref的正常灰度分布区间[g1,g2]为基准，计算训练图像块Xj灰度量化后的压缩测量值yj_norm，具体流程如下；
[0086] 1)参照步骤2中第2)步的方法依次计算选定的C幅训练图像块Xj的正常灰度分布区间[gj1,gj2]；
[0087] 2)以参照图像块Xref的正常灰度分布区间[g1,g2]为基准，根据训练图像块的压缩测量值yj，按照公式(2)计算训练图像块Xj在[gj1,gj2]区间量化后的压缩测量值yj_norm：
[0088]
[0089] 其中，φn表示观测矩阵Φ的第n列。
[0090] 步骤5：按照公式(3)计算训练图像块Xj灰度量化后的压缩测量值yj_norm的均值ymean：
[0091]
[0092] 步骤6：采用极大似然估计法，计算训练图像块Xj灰度量化后的压缩测量值的正态分布参数μnorm和σnorm，为后续判别有无疵点提供参照，其计算公式为：
[0093]
[0094]
[0095] 利用本发明方法在大幅面坯布织物图像上截取的无疵点图像块及其量化结果图如图1所示；
[0096] 第二阶段：疵点判别阶段
[0097] 步骤7：将采集待测坯布织物图像划分为互不重叠的图像块Xt，统计其灰度直方图Ht，并采用观测矩阵Φ获取其M×1维的压缩测量值yt，具体流程如下：
[0098] 1)将采集到的待测坯布织物图像划分为互不重叠的p×q大小的图像块Xt(t＝1,2,…n)，若不能等分，则不足补齐。为了保证图像块中疵点的稀疏性，其中p×q≥5×Kdef；
[0099] 2)分别统计各个图像块的灰度直方图Ht，并参照公式(1)计算待测图像块Xt的压缩测量值yt；
[0100] 步骤8：根据待测图像块Xt的灰度直方图Ht，计算待测图像块Xt的灰度分布区间[gt1,gt2]，并以参照图像块Xref的正常灰度分布区间[g1,g2]为基准，计算待测图像块Xt灰度量化后的压缩测量值yt_det，具体流程如下：
[0101] 1)根据待测图像块Xt的灰度直方图Ht，参照步骤2中第2)步的方法计算待测图像块Xt的灰度分布区间，记为[gt1,gt2]；
[0102] 2)以参照图像块Xref的正常灰度分布区间[g1,g2]为基准，根据待测图像块Xt的压缩测量值yt，参照公式(2)计算Xt在[gt1,gt2]区间量化后的压缩测量值yt_det；利用本发明方法在大幅面坯布织物图像上截取的疵点图像块及其量化结果图如图2所示；
[0103] 步骤9：根据步骤5和步骤6中对无疵点的训练图像块Xj(j＝1,2,...C)学习得到的参数ymean、μnorm和σnorm，按照公式(6)，判断待测的图像块Xt中是否含有疵点：
[0104]
[0105] 其中，c为一个正的常数，其取值与检测的精确度有关，根据数理统计理论中的“3σ”原则，c可以取3。在实际应用中，根据检测的精确度与速度间的权衡，c取[2,3]之间为宜；
[0106] 如果待测图像块Xt_def的压缩测量值yt_det与无疵点的训练图像块学习得到的参数ymean、μnorm和σnorm满足条件： 则认为待测图像块Xt_def含有疵点，即置为1，否则，不含疵点，置为0；利用本发明方法在大幅面坯布织物图像上截取的无疵点图像块和疵点图像块分别对应的 值的分布如图3所示；
[0107] 步骤10：若待测的图像块Xt中无疵点，返回，重新进入第二阶段处理其他待测图像块；若含有疵点，则进入第三阶段的处理；
[0108] 第三阶段：疵点定位阶段
[0109] 步骤11：根据待测图像块Xt的灰度直方图Ht，计算Xt主导像素的灰度分布区间[gt_m1,gt_m2]，并以参照图像块Xref主导像素的灰度分布区间[gm1,gm2]为基准，计算Xt灰度量化后的压缩测量值yt_mdet，具体流程如下：
[0110] 1)根据待测图像块Xt的灰度直方图Ht，参照步骤2中第3)步的方法计算Xt主导像素的灰度分布区间，记为[gt_m1,gt_m2]；
[0111] 2)以参照图像块Xref主导像素的灰度分布区间[gm1,gm2]为基准，根据待测图像块Xt的压缩测量值yt，参照公式(2)计算Xt在[gt_m1,gt_m2]区间量化后的压缩测量值yt_mdet；
[0112] 步骤12：按照公式(7)，滤除yt_mdet中无疵点信息：
[0113] yt_def＝yt_mdet-ymean (7)
[0114] 步骤13：根据提取的疵点信息yt_def和观测矩阵Φ，采用基追踪去噪算法，对疵点信号进行恢复，并通过重构、中值滤波以及二值化处理完成疵点定位，具体流程如下：
[0115] 1)按照公式(8)，对获取的疵点信息yt_def进行恢复，得到N×1维的信号[0116]
[0117] 其中，λ是一个正的常数，取值范围为[2,4]，其取值需考虑信息的恢复精度与计算速度间的权衡；
[0118] 2)将得到的 按列优先的顺序排列成与原始图像块同维数的p×q维矩阵Xt_rec；
[0119] 3)按照公式(9)，对Xt_rec进行3×3窗口的中值滤波，得到Xt_f：
[0120] Xt_f(i,j)＝median{Xt_rec(i-3,j-3),…,Xt_rec(i,j),…Xt_rec(i+3,j+3)} (9)[0121] 4)按公式(10)，对Xt_f进行二值化处理得到疵点的位置及轮廓信息Xt_b，即完成疵点定位：
[0122]
[0123] 利用本发明方法在大幅面坯布织物图像上截取的含有疵点信息的待测图像块中疵点的恢复及其定位结果图如图4所示。
