CN105373797A - 基于平均池化稀疏编码的煤岩识别方法 - Google Patents

Abstract

本发明公开了一种基于平均池化稀疏编码的煤岩识别方法，该方法从煤岩图像数据中学习煤岩的结构基元，所学习到结构基元捕捉到了煤岩图像本质的结构特征，因而具有很强的鉴别能力和对成像环境变化的鲁棒性，从而使得该方法具有很高的识别稳定性和识别正确率，能为自动化采掘、自动化放煤、自动化选矸等生产过程提供可靠的煤岩识别信息。

Description

基于平均池化稀疏编码的煤岩识别方法

技术领域

本发明涉及一种基于平均池化稀疏编码的煤岩识别方法，属于煤岩识别领域。

背景技术

煤岩识别即用一种方法自动识别出煤岩对象为煤或岩石。在煤炭生产过程中，煤岩识别技术可广泛应用于滚筒采煤、掘进、放顶煤开采、原煤选矸石等生产环节，对于减少采掘工作面作业人员或者实现无人化作业、减轻工人劳动强度、改善作业环境、实现煤矿安全高效生产具有重要意义。

有多种方法已应用于煤岩识别，如自然γ射线探测、雷达探测、应力截齿、红外探测、有功功率监测、震动检测、声音检测、粉尘检测、记忆截割等，但这些方法存在以下问题：①需要在现有设备上加装各类传感器获取信息，导致装置结构复杂，成本高。②采煤机滚筒、掘进机等设备在生产过程中受力复杂、振动剧烈、磨损严重、粉尘大，传感器部署比较困难，容易导致机械构件、传感器和电气线路受到损坏，装置可靠性差。③对于不同类型机械设备，传感器的最佳类型和信号拾取点的选择存在较大区别，需要进行个性化定制，系统的适应性差。

为解决上述问题，图像技术也越来越受到重视并研发了一些基于图像技术的煤岩识别方法，然而已有方法都是用人工主观设计的图像特征或者图像特征的组合来进行煤岩识别，人工设计的特征往往并不能准确抓住煤岩图像本质结构致使对因成像条件变化引起的图像数据变化不具有具鲁棒性，因而在识别稳定性和识别正确率上还有很大的不足。

需要一种解决或至少改善现有技术中固有的一个或多个问题的煤岩识别方法，以提高煤岩识别率和识别稳定性。

发明内容

因此，本发明的目的在于提供一种基于平均池化稀疏编码的煤岩识别方法，该方法从煤岩图像数据中学习煤岩的结构基元，所学习到结构基元捕捉到了煤岩图像本质的结构特征，因而具有很强的鉴别能力和对成像环境变化的鲁棒性，从而使得该方法具有很高的识别稳定性和识别正确率，能为自动化采掘、自动化放煤、自动化选矸等生产过程提供可靠的煤岩识别信息。

根据一种实施例形式，提供一种基于平均池化稀疏编码的煤岩识别方法，包括如下步骤：

A.采集一张已知煤(或岩)对象的图像；

B.从采集的图像中抽取N个图像块[x₁，x₂，...x_N]＝X∈R^p×N，p为图像块向量化后的维数；

C.用抽取的煤(或岩)图像块x_i通过解优化问题求出煤(或岩)图像的基元矩阵D＝[d₁，d₂，...d_K]∈R^p×K和N个图像块的稀疏系数矩阵U＝[u₁，u₂，...u_N]∈R^K×N，u_i＝[u_i1，u_i2，...u_iK]^T；

D.求煤(或岩)图像的基元响应分布z＝[z₁，z₂，...z_K]，第k个基元响应z_k为

$z_k=\frac{1}{T}{\mathrm{\Σ}}_{i=1}^N{u}_{ik}^2$ 其中， $T={\mathrm{\Σ}}_{k=1}^K{u}_{ik}^2$ 为归一化常数；

E.对于待识别的煤岩图像，用与步骤B相同的方法抽取N个图像块y_i，用步骤C中求得的D表达图像块y_i，通过解优化问题求出每一个图像块的系数 $v_i^{*}={\[v_{i1},v_{i2},...v_{iK}\]}^T;$

F.用与步骤D相同的方法求出待识别煤岩图像的基元响应分布

G.用度量与已知煤(或岩)对象的相似性，当距离大于设定的阈值时，判定为煤(或岩)，否则为岩(或煤)；

进一步特定的但非限制性地，步骤C的优化方法为：

C1.给D赋初值，设置迭代次数；

C2.固定D，用求出所有图像块的系数U；

C3.固定U，求 $D\∈\arg {\min }_D{\mathrm{\Σ}}_{i=1}^N\frac{1}{2}\left|\right|x_i-{\mathrm{Du}}_i|{|}_2^2;$

C4.C2和C3交替进行直到迭代结束。

附图说明

通过以下说明，附图实施例变得显而已见，其仅以结合附图描述的至少一种优选但非限制性实施例的示例方式给出。

图1是本发明所述煤岩识别方法的原理示意图。

具体实施方案

图1是本发明所述煤岩识别方法的原理示意图，主要包含3层：图像层，编码层及池化层，图像层给编码层提供输入，本实施例用从灰度图像中抽取图像块作为编码层的输入，也可以从图像中抽取图像特征集合如sift特征作为输入；编码层计算每个图像块用从图像数据中学习到的K个基元表达时的表达系数，若采用的计算方法使得表达系数中的非零元素很少，则称为稀疏编码，本实施例用l₁-norm优化计算表达系数，因而所得的系数是稀疏的；池化层计算所有表达系数的统计特性进而获得输入图像的特征表达，本实施例用了均值统计特性，具体的实施步骤如下：

A.从煤岩识别任务的现场如采煤工作面采集包含煤和岩的图像，从中截取只包含煤或者岩的图像区域，然后归一化到合适的大小如32*32像素大小作为样本图像；

B.从样本图像中抽取N个图像块[x₁，x₂，...x_N]＝X∈R^p×N，p为图像块向量化后的维数，如图像块大小取6*6像素，以步长为2个像素在样本图像中滑动采样图像块，对每个图像块向量进行标准化处理：用去掉图像亮度的均值，以消除亮度变化的影响，用对图像块向量进行归一化处理，其中，1_p表示p维全1向量，η为常数值；

C.用抽取的煤(或岩)图像块x_i通过解优化问题求出煤(或岩)图像的基元矩阵D＝[d₁，d₂，...d_K]∈R^p×K和N个图像块的稀疏系数矩阵U＝[u₁，u₂，...u_N]∈R^K×N，u_i＝[u_i1，u_i2，...u_iK]^T；由于采用的是l₁-norm优化，求出的u_i非常稀疏，即非零元素很少。

求解方法可采用交替最小化D和U的方法，即按下列步骤进行处理：

C1.给D赋初值，设置迭代次数；

C2.固定D，用求出所有图像块的稀疏系数U；

C3.固定U，求 $D\∈\arg {\min }_D{\mathrm{\Σ}}_{i=1}^N\frac{1}{2}\left|\right|x_i-{\mathrm{Du}}_i|{|}_2^2;$

C4.C2和C3交替进行直到迭代结束。

步骤C2的优化可用近似梯度算法，对每一个图像块x_i对应的系数u_i采用下列步骤优化：

1.给系数u赋初值，设置迭代次数；

2.在每一次迭代：

u←u+ξD^T(x-Du)，ξ为迭代步长；

$u\[k\]=\left\{\begin{array}{cc}u\[k\]-\mathrm{\λ}& ifu\[k\]\≥\mathrm{\λ}\\ u\[k\]+\mathrm{\λ}& ifu\[k\]\≤-\mathrm{\λ}\\ 0& otherwise\end{array}\right.,$ λ为给定的参数，k为基元元素索引；

3.重复2直到迭代结束。

步骤C3的优化可采用块坐标下降算法，用下列步骤优化：

1.B←XU^T，C←UU^T

2.Fork＝1，2，...，K

$d_k\←\frac{1}{C\[k,k\]}(b_k-{\mathrm{Dc}}_k)+d_k,$

$d_k\←\frac{1}{max\left(\right|\left|d_k\right|{|}_2),1)}{d}_k$

3.重复2直到收敛。

D.求煤(或岩)图像的基元响应分布z＝[z₁，z₂，...z_K]，将每个基元对所有图像块的响应求平均，即为该基元的响应，其中，为归一化常数；

E.对于待识别的煤岩图像，用与步骤B相同的方法抽取N个图像块y_i，用步骤C中求得的D表达图像块y_i，通过解优化问题求出每一个图像块的系数 $v_i^{*}={\[v_{i1},v_{i2},...v_{iK}\]}^T;$

F.用与步骤D相同的方法求出待识别煤岩图像的基元响应分布

G.用度量与已知煤(或岩)对象的相似性，当距离大于设定的阈值时，判定为煤(或岩)，否则为岩(或煤)。

Claims (2)

1.一种基于均值池化稀疏编码的煤岩识别方法，其特征在于包括以下步骤：

A.采集一张已知煤(或岩)对象的图像；

B.从采集的图像中抽取N个图像块[x₁，x₂，...x_N]＝X∈R^p×N，p为图像块向量化后的维数；

C.用抽取的煤(或岩)图像块x_i通过解优化问题求出煤(或岩)图像的基元矩阵D＝[d₁，d₂，...d_K]∈R^p×K和N个图像块的稀疏系数矩阵U＝[u₁，u₂，...u_N]∈R^K×N，u_i＝[u_i1，u_i2，...u_iK]^T；

D.求煤(或岩)图像的基元响应分布z＝[z₁，z₂，...z_K]，第k个基元响应z_k为

其中，为归一化常数；

E.对于待识别的煤岩图像，用与步骤B相同的方法抽取N个图像块y_i，用步骤C中求得的D表达图像块y_i，通过解优化问题求出每一个图像块的系数

F.用与步骤D相同的方法求出待识别煤岩图像的基元响应分布

G.用度量与已知煤(或岩)对象的相似性，当距离大于设定的阈值时，判定为煤(或岩)，否则为岩(或煤)。

2.根据权利要求1所述方法，其特征在于步骤C的优化方法为：

C1.给D赋初值，设定迭代次数；

C2.固定D，用求出所有图像块的稀疏系数U；

C3.固定U，求

C4.C2和C3交替进行直到迭代结束。