发明名称 基于紧邻域像素灰度联合分布特征的煤岩识别方法
摘要 本发明公开了一种基于紧邻域像素灰度联合分布特征的煤岩识别方法,该方法用紧邻域像素灰度联合分布特征来描述煤、岩图像特征信息,用学向量量化方法抽取煤岩特征信息空间关键字向量,用关键字向量直方图表征煤岩在不同形态下的模式;在识别时,待识别图像用与训练图像同样的方法抽取图像特征信息和建立直方图,然后和训练阶段学到的模式进行比较,用G统计来度量,用最近邻准则来识别。该方法用了煤、岩在不同照度、不同视点下的图像作为训练样本,受照度和成像视点变化影响小,识别率高,稳定性好。
申请公布号 CN103927528B 申请公布日期 2017.02.01
申请号 CN201410184275.5 申请日期 2014.05.05
申请人 中国矿业大学(北京) 发明人 伍云霞;刘毅
分类号 G06K9/00(2006.01)I;G06K9/46(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 代理人
主权项 一种基于紧邻域像素灰度联合分布特征的煤岩识别方法,其特征在于,包括以下步骤:A.将采集的煤、岩样本图像进行预处理并分成训练集SET_A和测试集SET_B将采集的每一张煤、岩样本图像,在图像的中心截取像素大小为N的子图像,并将每张子图像的灰度归一化为零均值和单位方差,选择处理过的一部分图像作为训练集SET_A,另一部分作为测试集SET_B;B.提取训练集SET_A和测试集SET_B中图像的特征信息向量分别提取训练集SET_A和测试集SET_B中每张煤图像I<sub>c</sub>和每张岩图像I<sub>r</sub>的特征信息向量{x<sub>i</sub>}<sub>c</sub>和{x<sub>i</sub>}<sub>r</sub>,所述{x<sub>i</sub>}<sub>c</sub>和{x<sub>i</sub>}<sub>r</sub>中的特征信息向量x<sub>i</sub>∈R<sup>n</sup>为该图像中像素<img file="FSB0000158187680000011.GIF" wi="555" he="63" />的紧邻域内的n个像素灰度的联合分布,即x<sub>i</sub>={I<sub>1</sub>,...,I<sub>n‑1</sub>,I<sub>i</sub>},I<sub>i</sub>为像素(i,j)的灰度值,I<sub>1</sub>,...,I<sub>n‑1</sub>为像素(i,j)邻域像素的灰度值,将x<sub>i</sub>归一化;C.从图像特征信息向量空间中提取K个关键字向量{w<sub>1</sub>,w<sub>2</sub>,...,w<sub>K</sub>},w<sub>s</sub>∈R<sup>n</sup>利用训练集中的特征信息向量{{x<sub>i</sub>}<sub>c</sub>,{x<sub>i</sub>}<sub>r</sub>}<sub>SET_A</sub>,首先用OLVQ1向量量化训练算法对关键字向量位置进行粗定位,再用LVQ2训练算法对关键字向量位置进行细调节,用测试集SET_B特征信息向量{{x<sub>i</sub>}<sub>c</sub>,{x<sub>i</sub>}<sub>r</sub>}<sub>SET_B</sub>对LVQ2的训练结果进行精度测试,直至所提取的K个关键字向量{w<sub>1</sub>,w<sub>2</sub>,...,w<sub>K</sub>}满足识别精度要求为止;D.将煤、岩模式用关键字向量直方图表征从训练集SET_A中分别选取煤、岩图像各M张,将选定的每一张图像的特征信息向量x<sub>i</sub>用与其最邻近的关键字向量标注,计算每个关键字向量出现的次数,作相对该图像总标注像素的归一化直方图,煤图像的归一化直方图用H<sub>c</sub>表示,岩图像的归一化直方图用H<sub>r</sub>表示,煤模式为所选定的M张煤图像的归一化直方图H<sub>c</sub>的集合,即{H<sub>c</sub>}<sub>M</sub>,岩模式为所选定的M张岩图像的归一化直方图H<sub>r</sub>的集合,即{H<sub>r</sub>}<sub>M</sub>;E.识别待识别图像给定待识别图像,在图像的中心截取像素大小为N的子图像,并将子图像的灰度归一化为零均值和单位方差,记处理过后的子图像为I<sub>x</sub>,用与步骤B中相同的方法提取I<sub>x</sub>的特征信息向量{x<sub>i</sub>}<sub>x</sub>,用与D中相同的方法对x<sub>i</sub>进行标注,计算I<sub>x</sub>关键字归一化直方图H<sub>x</sub>,用对数似然率即G统计度量H<sub>x</sub>与煤和岩石模式距离,计算公式如下:<maths num="0001"><math><![CDATA[<mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>H</mi><mi>x</mi></msub><mo>,</mo><mi>H</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>H</mi><mrow><mi>x</mi><mi>k</mi></mrow></msub><mi>log</mi><mi> </mi><msub><mi>H</mi><mi>k</mi></msub><mo>,</mo><mi>H</mi><mo>&Element;</mo><mo>{</mo><msub><mrow><mo>{</mo><msub><mi>H</mi><mi>c</mi></msub><mo>}</mo></mrow><mi>M</mi></msub><mo>,</mo><msub><mrow><mo>{</mo><msub><mi>H</mi><mi>r</mi></msub><mo>}</mo></mrow><mi>M</mi></msub><mo>}</mo></mrow>]]></math><img file="FSB0000158187680000012.GIF" wi="1089" he="155" /></maths>其中,H<sub>xk</sub>为待识别图像第k个关键字向量的概率,H<sub>k</sub>为煤模式或者岩模式第k个关键字向量的概率,使G(H<sub>x</sub>,H)最大的模式即为待识别图像模式。
地址 100083 北京市海淀区学院路丁11号中国矿业大学(北京)