CN104376313B - 用图像局部曲线方向分布识别煤岩的方法 - Google Patents
用图像局部曲线方向分布识别煤岩的方法 Download PDFInfo
- Publication number
- CN104376313B CN104376313B CN201410746274.5A CN201410746274A CN104376313B CN 104376313 B CN104376313 B CN 104376313B CN 201410746274 A CN201410746274 A CN 201410746274A CN 104376313 B CN104376313 B CN 104376313B
- Authority
- CN
- China
- Prior art keywords
- coal
- image
- calculated
- curve
- pixel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/42—Global feature extraction by analysis of the whole pattern, e.g. using frequency domain transformations or autocorrelation
- G06V10/422—Global feature extraction by analysis of the whole pattern, e.g. using frequency domain transformations or autocorrelation for representing the structure of the pattern or shape of an object therefor
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种用图像局部曲线方向分布识别煤岩的方法,识别过程分2步:A.计算已知煤岩对象的图像Ic、Ir和待识别的图像Io的局部曲线方向分布特征Hc、Hr和Ho;B.根据Ho与Hc、Hr的相似程度来确定待识别的图像Io的对象是煤或者岩石。本发明的方法提高了煤岩识别稳定性和识别正确率,为自动化采掘、自动化放煤、自动化选矸等生产过程提供可靠的煤岩识别信息,解决现有的基于图像技术的煤岩识别方法识别稳定性和识别正确率较低的问题。
Description
技术领域
本发明涉及一种用图像局部曲线方向分布识别煤岩的方法,属于图像识别技术领域。
背景技术
煤岩识别即用一种方法自动识别出煤岩对象为煤或岩石。在煤炭生产过程中,煤岩识别技术可广泛应用于滚筒采煤、掘进、放顶煤开采、原煤选矸石等生产环节,对于减少采掘工作面作业人员、减轻工人劳动强度、改善作业环境、实现煤矿安全高效生产具有重要意义。
已有多种煤岩识别方法,如自然γ射线探测法、雷达探测法、应力截齿法、红外探测法、有功功率监测法、震动检测法、声音检测法、粉尘检测法、记忆截割法等,但这些方法存在以下问题:①需要在现有设备上加装各类传感器获取信息,导致装置结构复杂,成本高。②采煤机滚筒、掘进机等设备在生产过程中受力复杂、振动剧烈、磨损严重、粉尘大,传感器部署比较困难,容易导致机械构件、传感器和电气线路受到损坏,装置可靠性差。③对于不同类型机械设备,传感器的最佳类型和信号拾取点的选择存在较大区别,需要进行个性化定制,系统的适应性差。
通过对煤、岩石块状样本的观察,发现煤和岩石的表面纹理在局部呈现曲线结构而其方向有很大的差异,当用成像技术获取煤岩信息时,这种差异性就隐含在图像数据里,因此,提取出这种特征可用来区分煤岩。已有的基于图像技术的煤岩识别方法在识别稳定性和识别正确率上还有很大的不足。
因此,需要一种解决或至少改善现有技术中固有的一个或多个问题的煤岩识别方法。
发明内容
本发明的目的在于提供一种用图像局部曲线方向分布识别煤岩的方法,该方法可以提高煤岩识别稳定性和识别正确率,为自动化采掘、自动化放煤、自动化选矸等生产过程提供可靠的煤岩识别信息,解决现有的基于图像技术的煤岩识别方法识别稳定性和识别正确率较低的问题。
根据一种实施例形式,提供一种用图像局部曲线方向分布识别煤岩的方法,给定一幅或者多幅已知煤岩对象图像Ic、Ir和一幅或者多幅待识别煤岩对象的图像Io,识别过程包括以下步骤:
A.用相同的计算方法分别计算所述图像Ic、Ir和Io的局部曲线方向分布Hc、Hr和Ho,所述计算方法具体包括如下步骤:
A1.计算描述像素点(x,y)邻域ρ内图像结构的矩阵T;
A2.计算像素点(x,y)邻域ρ内的梯度平均方向φ;
A3.以像素点(x,y)为中心,建立局部曲线度量坐标系vw,v沿着像素邻域ρ的梯度平均方向φ,w沿垂直于φ方向的切线方向,所述度量坐标系下经过点(v=p,w=0)的曲线可由描述,其中,k为曲线曲率;
A4.将所述度量坐标系下的曲线变换到坐标系pw’,变换关系为:
A5.在pw′坐标系下计算描述像素点(x,y)邻域ρ内图像结构的矩阵T′,并计算所述矩阵T′所描述的结构与度量坐标系下的曲线最相似时的特征值λ1,λ2;
A6.计算像素点(x,y)ρ邻域内结构为所述度量坐标系下的曲线时的测量置信度C;
A7.重复步骤A1~A6,遍历图像中的所有像素点;
A8.计算图像局部曲线方向分布H:
其中,hφ为由步骤A2所计算的角度值为φ时的概率,N2=N×N为图像的大小,M为当C>Cth时的像素个数,Cth为测量置信度阈值,mφ为当C>Cth时角度值为φ的像素个数;
B.计算Ho分别与Hc、Hr的相似性,并根据Ho与Hc、Hr的相似性关系,判别所述待识别图像Io的对象是煤或者岩石,具体包括以下步骤:
B1.计算Ho分别与Hc、Hr的相似性do-c和do-r
B2.根据以下判别规则进行识别,ε为给定阈值:
当do-c-do-r>ε时,为煤;
当do-c-do-r<-ε时,为岩;
当-ε≤do-c-do-r≤ε时,无法判断。
在进一步特定的但非限制性的形式中,所述步骤A1中,所述矩阵T为:
其中,Ix,Iy分别表示图像I中像素点(x,y)沿X方向和Y方向的微分,F表示在像素点(x,y)邻域ρ内的平均操作。
在进一步特定的但非限制性的形式中,所述步骤A2中,梯度平均方向φ为:
当Txx<Tyy,Txy<0,
当Txx<Tyy,Txy>0,
在进一步特定的但非限制性的形式中,所述的步骤A5中,所述矩阵T′为:
其中,Ip,Iw′表示沿p,w’方向的微分,Iv,Iw表示沿v,w方向的微分;
在进一步特定的但非限制性的形式中,所述的步骤A5中,所述矩阵T′所描述的结构与度量坐标系下的曲线最相似时的特征值λ1,λ2为:
该特征值是在k=k0时取得。
附图说明
通过以下说明,附图实施例变得显而易见,其仅以结合附图描述的至少一种优选但非限制性实施例的示例方式给出。
图1是本发明所述煤岩识别方法的流程图;
图2是本发明中xy坐标和vw坐标变换关系示意图。
具体实施方式
纹理是煤岩图像最主要的特征之一,煤岩图像纹理的局部方向是纹理特征的主要元素,通过观察和大量的实验研究发现,将煤岩图像的局部纹理用曲线结构来表征,曲线的方向分布能够有效地识别煤岩对象。
因此,本发明用煤岩图像的局部曲线方向分布特征来识别煤岩对象。图1是本发明方法的识别流程,具体步骤如下:
A.计算已知煤岩对象的图像Ic、Ir和待识别的图像Io的局部曲线方向分布Hc、Hr和Ho,为了进一步提高精度,可以分别计算多幅已知对象的煤图像和岩石图像的局部曲线方向分布,然后取均值作为最终的分布特征;
B.根据Ho与Hc、Hr的相似程度来确定待识别的图像Io的对象是煤或者岩石。
图像Ic、Ir和Io的局部曲线方向分布由以下步骤A1~A8分别计算,每次计算只需更换不同的图像即可:
A1.计算描述像素点(x,y)邻域ρ内图像结构的矩阵T
其中,Ix,Iy分别表示图像I中像素点(x,y)沿X方向和Y方向的微分,F表示在像素点(x,y)邻域ρ内的平均操作(本发明中ρ=5)。F表示的平均操作可用与高斯函数的卷积来实现;Ix和Iy分别用和计算,为高斯函数分别在X方向和Y方向的微分,σ表示其尺度大小(本发明中σ=1);
A2.计算像素点(x,y)邻域ρ内的梯度平均方向φ,本发明中φ∈[0,π];
当Txx<Tyy,Txy<0,
当Txx<Tyy,Txy>0,
A3.以像素点(x,y)为中心建立局部曲线度量坐标系vw,v沿着像素点(x,y)邻域ρ内的梯度平均方向φ,w沿垂直于φ方向的切线方向,如图2所示。该度量坐标系下经过点(v=p,w=0)的曲线可由描述,其中,k为曲线曲率;
A4.将度量坐标系下的曲线变换到坐标系pw’,变换关系为
A5.在pw’坐标系下计算描述像素点(x,y)邻域ρ内图像结构的矩阵T’
其中:Ip,Iw′表示沿p,w’方向的微分,经坐标反变换即w’=w,得到用vw坐标系表示的矩阵T’,其中Iv,Iw表示沿v,w方向的微分, 则矩阵T’的特征值分别为:
λ′1+λ′2表征了像素点(x,y)邻域ρ内总微分能量,λ′1表征了平均梯度方向即v方向上的微分能量,λ′2则表征了垂直于平均梯度方向的即w方向的微分能量,最小化λ′2即取曲率k=k0时,像素点(x,y)邻域ρ内图像结构最有可能是曲率为k0的曲线,此时矩阵T’特征值为:
计算时将vw坐标系变换到xy坐标系,两者的变换关系为下式,如图2所示。
v=ξcosφ+ηsinφ
w=-ξsinφ+ηcosφ
(ξ,η)为像素点在xy坐标系下的坐标。
A6.计算像素点(x,y)ρ邻域内结构为所述度量坐标系下的曲线时的测量置信度C;
C∈[0,1],C越大表明测量的结果越可信;
A7.重复A1~A6,遍历图像中的所有像素点;
A8.计算图像局部曲线方向分布H
其中,hφ为由步骤A2所计算的角度值为φ时的概率,N2=N×N为图像的大小,M为当C>Cth时的像素个数,Cth为测量置信度阈值,mφ为当C>Cth时角度值为φ的像素个数。Ic、Ir和待识别图像Io的局部曲线方向分布Hc、Hr和Ho均由步骤A1A8分别计算,每次计算时,只要更换不同的图像数据即可;
判别图像Io的对象是煤或者岩的步骤如下:
B1.计算Ho分别与Hc、Hr的相似性do-c和do-r
B2.根据以下判别规则进行识别,为增加鲁棒性,设定阈值参数ε,
当do-c-do-r>ε时,为煤;
当do-c-do-r<-ε时,为岩;
当-ε≤do-c-do-r≤ε时,无法判断。
Claims (1)
1.一种用图像局部曲线方向分布识别煤岩的方法,给定一幅或者多幅已知煤岩对象的图像Ic、Ir和一幅或者多幅待识别煤岩对象的图像Io,其特征在于所述的方法包括以下步骤:
A.用相同的计算方法分别计算所述图像Ic、Ir和Io的局部曲线方向分布Hc、Hr和Ho,所述计算方法具体包括如下步骤:
A1.计算描述像素点(x,y)邻域ρ内图像结构的矩阵T:
其中,Ix,Iy分别表示图像I中像素点(x,y)沿X方向和Y方向的微分,F表示在像素点(x,y)邻域ρ内的平均操作;
A2.计算像素点(x,y)邻域ρ内的梯度平均方向φ:
当Txx<Tyy,Txy<0,
当Txx<Tyy,Txy>0,
A3.以像素点(x,y)为中心,建立局部曲线度量坐标系vw,v沿着像素邻域ρ的梯度平均方向φ,w沿垂直于φ方向的切线方向,所述度量坐标系下经过点(v=p,w=0)的曲线可由描述,其中,k为曲线曲率;
A4.将所述度量坐标系下的曲线变换到坐标系pw’,变换关系为:
w’=w
A5.在pw’坐标系下计算描述像素点(x,y)邻域ρ内图像结构的矩阵T′:
其中,Iv,Iw表示沿v,w方向的微分,Ip,Iw’表示沿p,w’方向的微分;
则矩阵T′所描述的vw坐标系下像素点(x,y)邻域ρ内的图像曲线结构曲率k=k0时的特征值λ1、λ2为:
A6.计算像素点(x,y)ρ邻域内曲线结构曲率为k时的测量置信度C:
A7.重复步骤A1~A6,遍历图像中的所有像素点;
A8.计算图像局部曲线方向分布H:
其中,hφ为由步骤A2所计算的角度值为φ时的概率,N2=N×N为图像的大小,M为当C>Cth时的像素个数,Cth为测量置信度阈值,mφ为当C>Cth时角度值为φ的像素个数;
B.计算Ho分别与Hc、Hr的相似性,并根据Ho与Hc、Hr的相似性关系,判别所述待识别图像Io的对象是煤或者岩石,具体包括以下步骤:
B1.计算Ho分别与Hc、Hr的相似性do-c和do-r:
B2.根据以下判别规则进行识别,ε为给定阈值:
当do-c-do-c>ε时,为煤;
当do-c-do-r<-ε时,为岩;
当-ε≤do-c-do-c≤ε时,无法判断。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410746274.5A CN104376313B (zh) | 2014-12-10 | 2014-12-10 | 用图像局部曲线方向分布识别煤岩的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410746274.5A CN104376313B (zh) | 2014-12-10 | 2014-12-10 | 用图像局部曲线方向分布识别煤岩的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104376313A CN104376313A (zh) | 2015-02-25 |
CN104376313B true CN104376313B (zh) | 2017-05-17 |
Family
ID=52555211
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410746274.5A Expired - Fee Related CN104376313B (zh) | 2014-12-10 | 2014-12-10 | 用图像局部曲线方向分布识别煤岩的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104376313B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6666521B1 (en) * | 1999-05-11 | 2003-12-23 | American Mining Electronics, Inc. | System for controlling cutting horizons for continuous type mining machines |
CN103927553A (zh) * | 2014-05-05 | 2014-07-16 | 中国矿业大学(北京) | 基于多尺度微纹理和对比度联合分布的煤岩识别方法 |
CN103942576A (zh) * | 2014-04-09 | 2014-07-23 | 中国矿业大学(北京) | 一种用空域多尺度随机特征识别煤岩的方法 |
-
2014
- 2014-12-10 CN CN201410746274.5A patent/CN104376313B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6666521B1 (en) * | 1999-05-11 | 2003-12-23 | American Mining Electronics, Inc. | System for controlling cutting horizons for continuous type mining machines |
CN103942576A (zh) * | 2014-04-09 | 2014-07-23 | 中国矿业大学(北京) | 一种用空域多尺度随机特征识别煤岩的方法 |
CN103927553A (zh) * | 2014-05-05 | 2014-07-16 | 中国矿业大学(北京) | 基于多尺度微纹理和对比度联合分布的煤岩识别方法 |
Non-Patent Citations (3)
Title |
---|
New Equipment of Distinguishing Rock from Coal Based on Statistical Analysis of Fast Fourier Transform;GU Tao,LI Xu;《Global Congress on Intelligent System》;20090821;全文 * |
基于支持向量机的煤岩图像特征抽取与分类识别;孙继平等;《煤炭学报》;20130930;第38卷;全文 * |
综放采煤自动化煤岩识别传感器的研究;张国军等;《传感器与微系统》;20110125;第30卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104376313A (zh) | 2015-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Automated estimation of reinforced precast concrete rebar positions using colored laser scan data | |
CN103927514B (zh) | 一种基于随机局部图像特征的煤岩识别方法 | |
CN102509087B (zh) | 一种基于图像灰度共生矩阵的煤岩识别方法 | |
CN103927528B (zh) | 基于紧邻域像素灰度联合分布特征的煤岩识别方法 | |
CN102768022B (zh) | 采用数码照相技术的隧道围岩变形检测方法 | |
CN103927553B (zh) | 基于多尺度微纹理和对比度联合分布的煤岩识别方法 | |
CN104134074B (zh) | 一种基于激光扫描的煤岩识别方法 | |
CN104751192A (zh) | 基于图像块共生特征的煤岩识别方法 | |
AU2020102336A4 (en) | Method, system and device for identifying water-filled goaf area in coal mine | |
CN102913285B (zh) | 矿井顶板危险报警的方法 | |
Houshmand et al. | Rock type classification based on petrophysical, geochemical, and core imaging data using machine and deep learning techniques | |
CN103942576B (zh) | 一种用空域多尺度随机特征识别煤岩的方法 | |
Walton et al. | An approach for automated lithological classification of point clouds | |
CN104751193B (zh) | 基于距离约束相似性的煤岩识别方法 | |
Sundaresan et al. | Robustness of change detection algorithms in the presence of registration errors | |
CN104376313B (zh) | 用图像局部曲线方向分布识别煤岩的方法 | |
CN104778461A (zh) | 基于相似测度学习的煤岩识别方法 | |
CN109543268A (zh) | 基于kriging模型的TBM推进力主要影响因素的识别方法 | |
CN104463098A (zh) | 用图像的结构张量方向直方图特征识别煤岩 | |
Jones et al. | A method for determining field accuracy of mobile scanning devices for geomechanics applications | |
CN105350963B (zh) | 一种基于相关性度量学习的煤岩识别方法 | |
Tantum et al. | Comparison of algorithms for land mine detection and discrimination using ground penetrating radar | |
Bahrampour et al. | Ground characterization and roof mapping: Online sensor signal-based change detection | |
Lee et al. | Peak detection with pile-up rejection using multiple-template cross-correlation for MWD (measurement while drilling) | |
Zhu et al. | Research on feature vector construction and classification recognition of typical mine microseismic waveforms |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170517 Termination date: 20191210 |
|
CF01 | Termination of patent right due to non-payment of annual fee |