首页 / 专利库 / 电脑编程 / 图像处理 / 基于可控塔式分解和字典学习的煤岩识别方法

基于可控塔式分解和字典学习的岩识别方法

阅读:1058发布:2020-05-18

专利汇可以提供基于可控塔式分解和字典学习的岩识别方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了基于可控塔式分解和字典学习的 煤 岩识别方法。在样本训练阶段,采集多幅已知类别煤岩图像并进行预处理,对每幅预处理后的图像进行多级可控塔式分解得到多个方向子带,提取方向子带的相关参数并按一定顺序排列形成 特征向量 ,合并特征向量构成样本特征矩阵,在预设参数的条件下进行字典学习得到字典矩阵、权重矩阵和偏移行向量;在煤岩识别阶段,采集未知类别测试图像,与样本训练阶段的 图像处理 方式类似,对测试图像进行一系列处理,最终得到测试特征向量。根据测试特征向量、样本训练阶段的字典矩阵、权重矩阵、偏移行向量和字典学习的相关参数判断测试图像所属的煤岩类型。本发明可靠性好,准确识别率高,运行速度快。,下面是基于可控塔式分解和字典学习的岩识别方法专利的具体信息内容。

1.基于可控塔式分解和字典学习的岩识别方法,其特征在于,包括以下步骤:
QS1.在样本训练阶段,采集相同光照条件下的m幅煤炭样本图像和m幅岩石样本图像,截取不含非煤岩背景的子图像并对其进行灰度化处理,处理后的煤炭样本子图和岩石样本子图分别记为c1,c2,…,cm和r1,r2,…,rm;
QS2.采用面向可控塔式分解的3阶滤波器组分别对子图c1,c2,…,cm和r1,r2,…,rm进行N级可控塔式分解,每幅子图经过N级可控塔式分解以后得到4N个方向子带,子图经过第i级可控塔式分解以后的第d个方向子带记为 其中i为可控塔式分解的级数序号,i=1,
2,…,N,d为方向序号,d=0,1,2,3,分别对应0, 这4个方向;
QS3.分别由c1,c2,…,cm,r1,r2,…,rm的4N个方向子带提取得到2m个20N维特征向量,然后构建样本特征矩阵Y=[y1,y2,…,ym,ym+1,ym+2,…,y2m]∈R20N×2m及其对应的类别标签向量u=[1,1,…,1,2,2,…,2]∈R1×2m,其中Y的前m列y1,y2,…,ym分别为c1,c2,…,cm的特征向量,Y的后m列ym+1,ym+2,…,y2m分别为r1,r2,…,rm的特征向量,u是由2m个整数构成的行向量,在u中,1代表煤炭,2代表岩石,Y的第k列20N维特征向量yk所属的煤岩类别由u的第k个整数元素uk表示,k=1,2,…,2m;
QS4.设置参数η1,η2,η3和字典原子数2τ,其中η1>0,η2>0,η3>0,τ≤m,利用Y和u进行字典学习,得到字典矩阵D、权重矩阵W和偏移行向量b;
QS5.在煤岩识别阶段,获取相同光照条件下拍摄的未知类别煤岩图像,截取不含非煤岩背景的子图像并对其进行灰度化处理,得到未知类别子图tx;
QS6.对tx进行N级可控塔式分解,得到4N个方向子带;
QS7.由步骤QS6所述的4N个方向子带,提取tx的20N维特征向量yx;
QS8.利用yx,η2和步骤QS4中字典学习完成以后得到的D,W,b,执行煤岩类别的判定。
2.根据权利要求1所述的基于可控塔式分解和字典学习的煤岩识别方法,其特征在于,所述20N维特征向量的数据结构采用
的格式,T表示转置运算,20N维特征向量中的每一个α元素、β元素、γ元素、μ元素和δ元素可以分别表示为 和 i为可控塔式分解的级数序号,i=
1,2,…,N,d为方向序号,d=0,1,2,3,分别对应0, 这4个方向, 和 分别为的子带系数绝对值的均值和方差, 和 分别指 在满足子带系数服从概
密度函数为
的非对称广义高斯分布条件下的α,β和γ参数,其中随机变量x为子带系数,t为积分变量,[0,+∞)为积分区间,e为自然常数。
3.根据权利要求1所述的基于可控塔式分解和字典学习的煤岩识别方法,其特征在于,所述字典学习包括以下步骤:
QS31.构建初始字典矩阵D∈R20N×2τ,其中D的前τ列和后τ列分别从Y的前m列和后m列中随机抽取得到;
QS32.构建稀疏矩阵H∈R2τ×2m并用H=(DTD+η2I2τ)-1DTY初始化,其中I2τ为2τ阶单位矩阵,T表示转置运算,-1表示求逆运算,构建权重矩阵W∈R2τ×2和2维偏移行向量b并分别初始化为零矩阵和零向量;
QS33.把H的每一列视为1个2τ维数据样本,以u为与H相对应的类别标签向量,η1为正则化参数,通过学习线性支持向量机更新权重矩阵W和偏移行向量b;
QS34.对于H的前m列数据的更新,反复求问题描述为
2τ×1
的解h'∈R ,并用h'更新H的第k列hk,其中yk为Y的第k列,k=1,2,…,m,w1和w2分别为W的第1列和第2列,b1和b2分别为b的第1个元素和第2个元素,T表示转置运算;
QS35.对于H的后m列数据的更新,反复求问题描述为
的解h*∈R2τ×1,并用h*更新H的第k列hk,其中yk为Y的第k列,k=(m+1),(m+2),…,2m,w1和w2分别为W的第1列和第2列,b1和b2分别为b的第1个元素和第2个元素,T表示转置运算;
QS36.求问题描述为
* 20N×2τ *
的解D∈R ,其中||·||F为矩阵的Frobenius范数, 为D的第j列,j=1,2,…,2τ;
QS37.把H的每一列视为1个2τ维数据样本,以u为与H相对应的类别标签向量,η1为正则化参数,通过学习线性支持向量机更新权重矩阵W和偏移行向量b;
QS38.如果满足 其中||·||F为矩阵的Frobenius范数,那么用D*更新
D,接着执行步骤QS34-QS37;如果满足 那么用D*更新D;
QS39.完成字典学习,输出D,W和b。
4.根据权利要求1所述的基于可控塔式分解和字典学习的煤岩识别方法,其特征在于,所述执行煤岩类别的判定包括以下步骤:
QS401.由hx=(DTD+η2I2τ)-1DTyx计算得到稀疏向量hx,其中T表示转置运算,I2τ为2τ阶单位矩阵,-1表示求逆运算;
QS402.如果满足(w1Thx+b1)≥(w2Thx+b2),那么判定tx为煤炭;反之,则判定tx为岩石。

说明书全文

基于可控塔式分解和字典学习的岩识别方法

技术领域

[0001] 本发明涉及基于可控塔式分解和字典学习的煤岩识别方法,属于图像识别技术领域。

背景技术

[0002] 煤岩识别是指通过某种手段辨别煤炭和岩石的相关技术,它在采煤机滚筒高度调整、综采放顶煤过程控制、选煤厂矸石初步筛选等方面有着良好的应用前景。自20世纪中叶以来,中国、美国、南非、澳大利亚等世界主要产煤国在寻找或尝试煤岩识别技术上做了大量工作,涌现出γ射线探测法、雷达探测法、红外探测法、震动检测法、声音检测法等多种煤岩识别方法。然而复杂的地质条件、多样的采煤工艺、恶劣的开采环境等因素使得现有煤岩识别方法的通用性不强,成熟的煤岩识别方法一直没有问世。煤炭和岩石的视觉差异明显,而且同一矿区或矿井的煤炭种类和岩石种类相对单一,再加上图像识别技术日新月异,这就给从图像分析度发明煤岩识别方法带来了机遇。

发明内容

[0003] 为了克服现有煤岩识别方法存在的不足,本发明提出基于可控塔式分解和字典学习的煤岩识别方法,该方法具有实时性强、识别率高、稳健性好等优点,能够进一步推进煤炭开采与加工的少人化甚至无人化。
[0004] 本发明所述的煤岩识别方法采用如下技术方案实现,包括样本训练阶段和煤岩识别阶段,具体步骤如下:
[0005] RS1.在样本训练阶段,采集相同光照条件下的m幅煤炭样本图像和m幅岩石样本图像,截取不含非煤岩背景的子图像并对其进行灰度化处理,处理后的煤炭样本子图和岩石样本子图分别记为c1,c2,…,cm和r1,r2,…,rm;
[0006] RS2.采用面向可控塔式分解的3阶滤波器组分别对子图c1,c2,…,cm和r1,r2,…,rm进行N级可控塔式分解,每幅子图经过N级可控塔式分解以后得到4N个方向子带,子图经过第i级可控塔式分解以后的第d个方向子带记为 其中i为可控塔式分解的级数序号,i=1,2,…,N,d为方向序号,d=0,1,2,3,分别对应0, 这4个方向;
[0007] RS3.分别由步骤RS2所述的所有子图的4N个方向子带提取得到2m个20N维特征向20N×2m
量,然后构建样本特征矩阵Y=[y1,y2,…,ym,ym+1,ym+2,…,y2m]∈R 及其对应的类别标签向量u=[1,1,…,1,2,2,…,2]∈R1×2m,其中Y的前m列y1,y2,…,ym分别为c1,c2,…,cm的特征向量,Y的后m列ym+1,ym+2,…,y2m分别为r1,r2,…,rm的特征向量,u是由2m个整数构成的行向量,在u中,1代表煤炭,2代表岩石,Y的第k列20N维特征向量yk所属的煤岩类别由u的第k个整数元素uk表示,k=1,2,…,2m;
[0008] RS4.设置参数η1,η2,η3和字典原子数2τ,其中η1>0,η2>0,η3>0,τ≤m,利用Y和u进行字典学习,得到字典矩阵D、权重矩阵W和偏移行向量b;
[0009] RS5.在煤岩识别阶段,获取相同光照条件下拍摄的未知类别煤岩图像,截取不含非煤岩背景的子图像并对其进行灰度化处理,得到未知类别子图tx;
[0010] RS6.与上述步骤RS2类似,对tx进行N级可控塔式分解,得到4N个方向子带;
[0011] RS7.由步骤RS6所述的4N个方向子带,提取tx的20N维特征向量yx;
[0012] RS8.利用yx,η2和步骤RS4中字典学习完成以后得到的D,W,b,执行煤岩类别的判定。
[0013] 步骤 RS 3和步 骤R S7所 述的 20 N维特 征向 量的数 据结 构采 用的格式,T表示转置运算,20N维特征向量中
的每一个α元素、β元素、γ元素、μ元素和δ元素可以分别表示为 和 i为可控塔式分解的级数序号,i=1,2,…,N,d为方向序号,d=0,1,2,3,分别对应0, 这
4个方向, 和 分别为 的子带系数绝对值的均值和方差, 和 分别指
在满足子带系数服从概率密度函数为
[0014]
[0015] 的非对称广义高斯分布条件下的α,β和γ参数,其中随机变量x为子带系数,t为积分变量,[0,+∞)为积分区间,e为自然常数。
[0016] 步骤RS4所述的字典学习包括以下步骤:
[0017] RS31.构建初始字典矩阵D∈R20N×2τ,其中D的前τ列和后τ列分别从Y的前m列和后m列中随机抽取得到;
[0018] RS32.构建稀疏矩阵H∈R2τ×2m并用H=(DTD+η2I2τ)-1DTY初始化,其中I2τ为2τ阶单位矩阵,T表示转置运算,-1表示求逆运算,构建权重矩阵W∈R2τ×2和2维偏移行向量b并分别初始化为零矩阵和零向量;
[0019] RS33.把H的每一列视为1个2τ维数据样本,以u为与H相对应的类别标签向量,η1为正则化参数,通过学习线性支持向量机更新权重矩阵W和偏移行向量b;
[0020] RS34.对于H的前m列数据的更新,反复求问题描述为
[0021]
[0022] 的解h'∈R2τ×1,并用h'更新H的第k列hk,其中yk为Y的第k列,k=1,2,…,m,w1和w2T分别为W的第1列和第2列,b1和b2分别为b的第1个元素和第2个元素,表示转置运算;
[0023] RS35.对于H的后m列数据的更新,反复求问题描述为
[0024]
[0025] 的解h*∈R2τ×1,并用h*更新H的第k列hk,其中yk为Y的第k列,k=(m+1),(m+2),…,2m,w1和w2分别为W的第1列和第2列,b1和b2分别为b的第1个元素和第2个元素,T表示转置运算;
[0026] RS36.求问题描述为
[0027]
[0028] 的解D*∈R20N×2τ,其中||·||F为矩阵的Frobenius范数, 为D*的第j列,j=1,2,…,2τ;
[0029] RS37.把H的每一列视为1个2τ维数据样本,以u为与H相对应的类别标签向量,η1为正则化参数,通过学习线性支持向量机更新权重矩阵W和偏移行向量b;
[0030] RS38.如果满足 其中||·||F为矩阵的Frobenius范数,那么用D*更新D,接着执行步骤RS34-RS37;如果满足 那么用D*更新D;
[0031] RS39.完成字典学习,输出D,W和b。
[0032] 步骤RS8所述的执行煤岩类别的判定,包括以下步骤:
[0033] RS401.由hx=(DTD+η2I2τ)-1DTyx计算得到稀疏向量hx,其中T表示转置运算,I2τ为2τ阶单位矩阵,-1表示求逆运算;
[0034] RS402.如果满足(w1Thx+b1)≥(w2Thx+b2),那么判定tx为煤炭;反之,则判定tx为岩石。附图说明
[0035] 图1是基于可控塔式分解和字典学习的煤岩识别方法的基本流程图
[0036] 图2是本发明所述字典学习的基本流程图。

具体实施方式

[0037] 在对我国山西、河南等地主要煤种和岩种的图像进行实验分析的基础上,本发明提出了基于可控塔式分解和字典学习的煤岩识别方法,该方法可以有效判别煤炭和岩石。下面结合附图和具体实施方式对本发明作进一步的详细描述。
[0038] 参照图1,基于可控塔式分解和字典学习的煤岩识别方法的具体步骤如下:
[0039] SS1.在样本训练阶段,采集相同光照条件下的m幅煤炭样本图像和m幅岩石样本图像,截取不含非煤岩背景的子图像并对其进行灰度化处理,处理后的煤炭样本子图和岩石样本子图分别记为c1,c2,…,cm和r1,r2,…,rm;
[0040] SS2.采用面向可控塔式分解的3阶滤波器组分别对子图c1,c2,…,cm和r1,r2,…,rm进行N级可控塔式分解,每幅子图经过N级可控塔式分解以后得到4N个方向子带,子图经过第i级可控塔式分解以后的第d个方向子带记为 其中i为可控塔式分解的级数序号,i=1,2,…,N,d为方向序号,d=0,1,2,3,分别对应0, 这4个方向;
[0041] SS3.分别由步骤SS2所述的所有子图的4N个方向子带提取得到2m个20N维特征向20N×2m
量,然后构建样本特征矩阵Y=[y1,y2,…,ym,ym+1,ym+2,…,y2m]∈R 及其对应的类别标签向量u=[1,1,…,1,2,2,…,2]∈R1×2m,其中Y的前m列y1,y2,…,ym分别为c1,c2,…,cm的特征向量,Y的后m列ym+1,ym+2,…,y2m分别为r1,r2,…,rm的特征向量,u是由2m个整数构成的行向量,在u中,1代表煤炭,2代表岩石,Y的第k列20N维特征向量yk所属的煤岩类别由u的第k个整数元素uk表示,k=1,2,…,2m;
[0042] SS4.设置参数η1,η2,η3和字典原子数2τ,其中η1>0,η2>0,η3>0,τ≤m,利用Y和u进行字典学习,得到字典矩阵D、权重矩阵W和偏移行向量b;
[0043] SS5.在煤岩识别阶段,获取相同光照条件下拍摄的未知类别煤岩图像,截取不含非煤岩背景的子图像并对其进行灰度化处理,得到未知类别子图tx;
[0044] SS6.与上述步骤SS2类似,对tx进行N级可控塔式分解,得到4N个方向子带;
[0045] SS7.由步骤SS6所述的4N个方向子带,提取tx的20N维特征向量yx;
[0046] SS8.利用yx,η2和步骤SS4中字典学习完成以后得到的D,W,b,执行煤岩类别的判定。
[0047] 步骤 SS 3和步 骤S S7所 述的 20 N维特 征向 量的数 据结 构采 用的格式,T表示转置运算,20N维特征向量中
的每一个α元素、β元素、γ元素、μ元素和δ元素可以分别表示为 和 i为可控塔式分解的级数序号,i=1,2,…,N,d为方向序号,d=0,1,2,3,分别对应0, 这
4个方向, 和 分别为 的子带系数绝对值的均值和方差, 和 分别指
在满足子带系数服从概率密度函数为
[0048]
[0049] 的非对称广义高斯分布条件下的α,β和γ参数,其中随机变量x为子带系数,t为积分变量,[0,+∞)为积分区间,e为自然常数。
[0050] 参照图2,步骤SS4所述的字典学习包括以下步骤:
[0051] SS31.构建初始字典矩阵D∈R20N×2τ,其中D的前τ列和后τ列分别从Y的前m列和后m列中随机抽取得到;
[0052] SS32.构建稀疏矩阵H∈R2τ×2m并用H=(DTD+η2I2τ)-1DTY初始化,其中I2τ为2τ阶单位矩阵,T表示转置运算,-1表示求逆运算,构建权重矩阵W∈R2τ×2和2维偏移行向量b并分别初始化为零矩阵和零向量;
[0053] SS33.把H的每一列视为1个2τ维数据样本,以u为与H相对应的类别标签向量,η1为正则化参数,通过学习线性支持向量机更新权重矩阵W和偏移行向量b;
[0054] SS34.对于H的前m列数据的更新,反复求问题描述为
[0055]
[0056]
[0057] 的解h'∈R2τ×1,并用h'更新H的第k列hk,其中yk为Y的第k列,k=1,2,…,m,w1和w2分别为W的第1列和第2列,b1和b2分别为b的第1个元素和第2个元素,T表示转置运算;
[0058] SS35.对于H的后m列数据的更新,反复求问题描述为
[0059]
[0060] 的解h*∈R2τ×1,并用h*更新H的第k列hk,其中yk为Y的第k列,k=(m+1),(m+2),…,2m,w1和w2分别为W的第1列和第2列,b1和b2分别为b的第1个元素和第2个元素,T表示转置运算;
[0061] SS36.求问题描述为
[0062]
[0063] 的解D*∈R20N×2τ,其中||·||F为矩阵的Frobenius范数, 为D*的第j列,j=1,2,…,2τ;
[0064] SS37.把H的每一列视为1个2τ维数据样本,以u为与H相对应的类别标签向量,η1为正则化参数,通过学习线性支持向量机更新权重矩阵W和偏移行向量b;
[0065] SS38.如果满足 其中||·||F为矩阵的Frobenius范数,那么用D*更新D,接着执行步骤SS34-SS37;如果满足 那么用D*更新D;
[0066] SS39.完成字典学习,输出D,W和b。
[0067] 步骤SS8所述的执行煤岩类别的判定,包括以下步骤:
[0068] SS401.由hx=(DTD+η2I2τ)-1DTyx计算得到稀疏向量hx,其中T表示转置运算,I2τ为2τ-1阶单位矩阵,表示求逆运算;
[0069] SS402.如果满足(w1Thx+b1)≥(w2Thx+b2),那么判定tx为煤炭;反之,则判定tx为岩石。
[0070] 需要指出的是,以上所述实施实例用于进一步说明本发明,实施实例不应被视为限制本发明的范围。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈