一种用于煤岩识别的b-cdtm纹理特征提取方法
【技术领域】
[0001] 本发明涉及一种用于煤岩识别的B-⑶TM纹理特征提取方法，属于煤岩识别技术 领域。
【背景技术】
[0002] 煤岩识别即用一种方法自动识别出煤岩对象为煤或岩。在煤炭生产过程中煤岩识 别技术可广泛应用于滚筒采煤、掘进、放顶煤开采、原煤选矸石等生产环节，对于减少采掘 工作面工作人员、减轻工人劳动强度、改善作业环境、实现煤矿安全高效生产具有重要的意 义
[0003] 目前，国内外煤岩识别主要有以下两类方法，第一种是煤层厚度的测量方法，通过 反馈的信息调整滚筒的高度，包括各种射线和电磁波，例如自然γ射线探测法、电磁波探 测法以及雷达探测法；另一种是煤和岩石界面的测量方法，通过采煤机自生的响应来调整 滚筒高度，包括应力截齿法、震动检测法、红外线检测法等。以上各种方法有各自的优缺点， 但都是采用传感器进行识别。一方面在不同的煤矿产区，需要选取不同的传感器，并需要安 装在不同的机器设备上；另一方面在煤炭开采环境中，传感器本身容易出现故障和失灵的 情况，无疑都加大了人力和物力的浪费。
[0004] 为解决上述问题，图像技术越来越受到重视并研发了一些图像技术的煤岩识 别方法，然而在已有的方法中，例如基于灰度共生矩阵（gray level co-occurrence matrix)方法，十字对角纹理矩阵（cross-diagonal texture matrix)方法，用来分析 图像纹理的特点.但是，此类方法都没有实际考虑煤岩图像识别的实时性，例如：煤矿 井下的采煤机需不间断的工作，因此对于用于图像技术进行煤岩识别需考虑识别的实 时性。对于煤岩图像识别中，分析图像纹理特点消耗了大量的时间，为了解决该问题， B-CDTM(binary-cross-diagonal texture matrix)纹理特征提取方法，通过实验比较，相 比灰度共生矩阵方法，十字对角纹理矩阵方法，可以很好的改善煤岩识别的实时性。
【发明内容】
[0005] 本发明的目的在于提供一种用于煤岩识别的B-⑶TM纹理特征提取方法，该纹理 特征提取方法很好的解决了的煤岩图像识别的实时性的要求，能够快速方便的提取出煤岩 图像的纹理特征，并且能够为煤岩识别提供准确的识别信息，为自动化采掘提供必要的信 息。用以解决现有的煤岩识别实时性不足，浪费人力物力的问题。
[0006] 为实现上述目的，本发明的方案是一种用于煤岩识别的B-CDTM纹理特征提取方 法，包括以下步骤：
[0007] (1)将煤岩样本灰度图像中的3X3的图像块的中心像素的八连通像素分为两组， 每组4个元素，即一组为正十字线上的四相邻像素和另一组对角线上的四相邻像素；
[0008] (2)根据所在像素值与中心像素值相比较，两组中每个元素分别取值0或1，即小 于中心像素的值则取〇,大于或者等于则取1 ;
[0009] (3)每组四个元素按照纹理单元值计算公式，获得中心像素点的十字纹理单元值 和对角纹理单元值，这些单元值的取值范围文〇--15之间的整数；
[0010] (4)计算图像中除边缘外的所有像素值的十字纹理单元值和对角纹理单元值，并 以十字纹理单元值为X坐边和以对角纹理单元值为Y坐标的矩阵的对应元素值累加，得到 一个B-CDTM矩阵，
[0011] (5)基于B-⑶TM矩阵，计算大量统计特征，最终将这些统计特征当做煤岩图像的 纹理特征，构造成特征向量。
[0012] 进一步地，根据本发明所述的煤岩图像的纹理特征的提取方法，步骤（3)中所述 的纹理单元值的具体计算公式如下：
[0013] 十字纹理单元值：.
[0014] 对角纹理单元值：
[0015] 其中五为3 x 3的图像块中正十字线上的四相邻像素组中元素的取值，即为〇或 1 ;其中五卩为3X3的图像块中对角线线上的四相邻像素组中元素的取值，即为0或1 ;
[0016] 进一步地，根据上述步骤中：纹理单元值的具体计算公式中的的起始 位置的不同，可得到四种十字对角纹理矩阵：T1、!'2、!'3、!' 4。
[0017] 进一步地，根据本发明所述的煤岩图像的纹理特征的提取方法，所述的步骤（3) 中的十字纹理单元值用NB CTU表示，对角纹理单元值用ΝΒ _表示。
[0018] 本发明达到的有益效果：本发明所述的一种用于煤岩识别的Β-⑶ΤΜ纹理特征提 取方法，本发明在原始CDTM纹理特征提取方法的基础上，将每组中每个像素的取值进行二 值化，即取〇或者1，对于纹理单元值的计算公式进行了改善，在保证原始CDTM纹理特征提 取方法的有效性的基础上，并且达到了良好的实时性。因此，考虑到煤矿井下工作环境的需 求，本发明更加适用于煤岩识别的纹理特征提取方法。
【附图说明】
[0019] 图1为本发明实施例提供的方法的步骤流程图；
[0020] 图2为⑶ΤΜ的方法的步骤流程图；
[0021] 图3为本发明所述的4种十字对角纹理矩阵的构造原理图；
[0022] 图4为本发明所述的方法与⑶ΤΜ方法的效果比对图；
[0023] 其中，图3中，说明了通过该变每组各个元素的位置，得到的所有可能的十字纹理 单元与对角纹理单元的值；图4,本发明所述的方法与CDTM方法的效果比对图：分别使用 Β-CDTM与CDTM提取320张64X64的煤岩样本灰度图像的纹理特征，并采用SRC (Sparse Representation-based Classification)进行分类，对纹理特征提取时间（s)，煤岩图像识 别率（％ )进行比对。
【具体实施方式】
[0024] 下面结合附图和【具体实施方式】对本发明作进一步的详细阐述；本发明的方案是一 种用于煤岩识别的B-CDTM纹理特征提取方法，包括以下步骤：
[0025] (1)如图1所示，将煤岩样本灰度图像中的3X3的图像块的中心像素的八连通像 素分为两组，每组4个元素，即一组为正十字线上的四相邻像素和另一组对角线上的四相 邻像素；
[0026] (2)如图1所示，根据所在像素值与中心像素值相比较，两组中每个元素分别取值 〇或1，即小于中心像素的值则取〇,大于或者等于则取1 ;
[0027] (3)如图1所示，每组四个元素按照纹理单元值计算公式，获得中心像素点的十字 纹理单元值和对角纹理单元值，这些单元值的取值范围文0--15之间的整数；
[0028] (4)如图1所示，计算图像中除边缘外的所有像素值的十字纹理单元值和对角纹 理单元值，并以十字纹理单元值为X坐边和以对角纹理单元值为Y坐标的矩阵的对应元素 值累加，得到一个B-⑶TM矩阵，
[0029] (5)基于B-⑶TM矩阵，计算大量统计特征，最终将这些统计特征当做煤岩图像的 纹理特征，构造成特征向量。
[0030] 进一步地，根据本发明所述的煤岩图像的纹理特征的提取方法，步骤（3)中所述 的纹理单元值的具体计算公式如下：
[0031] 十字纹理单元值：
[0032] 对角纹理单元值
[0033] 其中五e为3 X 3的图像块中正十字线上的四相邻像素组中元素的取值，即为0或 1 ;其中五卩为3X3的图像块中对角线线上的四相邻像素组中元素的取值，即为0或1 ;
[0034] 进一步地，如图3所示，根据上述步骤中：纹理单元值的具体计算公式中的 、五Λ的起始位置的不同，可得到四种十字对角纹理矩阵：T1、T2、T 3、T4。
[0035] 进一步地，根据本发明所述的煤岩图像的纹理特征的提取方法，所述的步骤（3) 中的十字纹理单元值用NB CTU表示，对角纹理单元值用ΝΒ _表示。
[0036] 进一步地，根据本发明所述的煤岩图像的纹理特征的提取方法，步骤（5)中所述 的纹理特征的具体计算过程如下：
[0037] a.将十字对角纹理矩阵Τ进行归一化处理，计算公式 其 中，Ta,b,。矩阵（i，j)对应的值为T(i，j);
[0051] 〇·将得到的特征值组成特征向量 q = {F1，F2，F3，F4，F5，F6，F7，F8，F9，F10，F11， F12, F13};
[0052] p.将得到的特征向量归一化，得到归一化的特征向量Q ;
[0053] q.将分别得到每张煤岩图像的十字对角纹理矩阵T1、T2、T 3、T4的特征向量Q \ Q2、 Q3、Q4。
【主权项】
1. 一种用于煤岩识别的B-CDTM纹理特征提取方法，其特征在于，所述的纹理特征提取 方法包括如下步骤： (1) 将煤岩样本灰度图像中的3X3的图像块的中心像素的八连通像素分为两组，每组 4个元素，即一组为正十字线上的四相邻像素和另一组对角线上的四相邻像素； (2) 根据所在像素值与中心像素值相比较，两组中每个元素分别取值O或1，即小于中 心像素的值则取〇,大于或者等于则取1 ; (3) 每组四个元素按照纹理单元值计算公式，获得中心像素点的十字纹理单元值和对 角纹理单元值，这些单元值的取值范围为0--15之间的整数； (4) 计算图像中除边缘外的所有像素值的十字纹理单元值和对角纹理单元值，并以十 字纹理单元值为X坐边和以对角纹理单元值为Y坐标的矩阵的对应元素值累加，得到一个 B-CDTM 矩阵， (5) 基于B-CDTM矩阵，计算大量统计特征，最终将这些统计特征当做煤岩图像的纹理 特征，构造成特征向量。2. 根据权利1所述的一种用于煤岩识别的B-CDTM纹理特征提取方法，其特征在于所述 步骤（3)中的四个元素的纹理单元值的计算公式是：其中五为3X3的图像块中正十字线上的四相邻像素组中元素的取值，即为0或1 ; 其中五^为3X3的图像块中对角线线上的四相邻像素组中元素的取值，即为0或1。3. 根据权利1所述的一种用于煤岩识别的B-CDTM纹理特征提取方法，其特征在与所述 的步骤（3)中的十字纹理单元值用Nbctu表示，对角纹理单元值用N bdtu表示。
【专利摘要】本发明提供了一种用于煤岩识别的B-CDTM纹理特征提取方法，步骤一，将煤岩样本灰度图像的3×3的图像块的中心像素的八连通像素分为两组，每组4个元素；步骤二，根据所在像素值与中心像素值相比较，两组中每个元素分别取值0或1；步骤三，获得中心像素点的十字纹理单元值和对角纹理单元值；步骤四，以十字纹理单元值为X坐边和以对角纹理单元值为Y坐标的矩阵的对应元素值累加，得到一个B-CDTM矩阵；步骤五，基于B-CDTM矩阵，计算大量统计特征，构造成特征向量。本发明通过改善原始CDTM的构造方式，使得提取煤岩图像的纹理特征的速度更加快捷，使得煤岩图像识别的实时性得到很好的改善。
【IPC分类】G06K9/46
【公开号】CN105354577
【申请号】CN201510695903
【发明人】孙继平, 杨坤
【申请人】中国矿业大学(北京)
【公开日】2016年2月24日
【申请日】2015年10月26日