主权项 |
一种基于协同稀疏编码的煤岩识别方法,其特征在于包括以下步骤:A.各采集一张已知煤和岩对象的图像;B.从煤和岩图像中各抽取N个图像块[x<sub>1</sub>,x<sub>2</sub>,...x<sub>N</sub>]=X∈R<sup>p×N</sup>,p为图像块向量化后的维数;C.分别用抽取的煤和岩图像块x<sub>i</sub>通过解优化问题<img file="FSA0000123054760000011.GIF" wi="757" he="79" />求出煤和岩图像的基元矩阵D<sub>c</sub>和D<sub>r</sub>,D<sub>c or r</sub>=[d<sub>1</sub>,d<sub>2</sub>,...d<sub>K</sub>]∈R<sup>p×K</sup>;D.通过解优化问题<img file="FSA0000123054760000012.GIF" wi="841" he="76" />分别求出煤和岩图像块x<sub>i</sub>用D=[D<sub>c</sub> D<sub>r</sub>]表达的每一个图像块的系数u<sub>i</sub>=[u<sub>i1</sub>,u<sub>i2</sub>,...u<sub>i2K</sub>]<sup>T</sup>;E.分别求出煤和岩图像的基元响应分布z<sub>c</sub>和z<sub>r</sub>,其中z<sub>c or r</sub>=[z<sub>1</sub>,z<sub>2</sub>,...z<sub>2K</sub>],第k个基元响应<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>z</mi><mi>k</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>T</mi></mfrac><msubsup><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msubsup><mi>u</mi><mrow><mi>i</mi><mi>k</mi></mrow><mn>2</mn></msubsup><mo>,</mo></mrow>]]></math><img file="FSA0000123054760000013.GIF" wi="378" he="82" /></maths>其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><msubsup><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mn>2</mn><mi>K</mi></mrow></msubsup><msubsup><mi>u</mi><mrow><mi>i</mi><mi>k</mi></mrow><mn>2</mn></msubsup></mrow>]]></math><img file="FSA0000123054760000014.GIF" wi="282" he="71" /></maths>为归一化参数;F.对于待识别的煤岩图像,用与步骤B相同的方法抽取N个图像块y<sub>i</sub>,用步骤C中求得的D=[D<sub>c</sub> D<sub>r</sub>]表达图像块y<sub>i</sub>,通过解优化问题<img file="FSA0000123054760000015.GIF" wi="835" he="77" />求出每一个图像块的系数<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>v</mi><mi>i</mi><mo>*</mo></msubsup><mo>=</mo><msup><mrow><mo>[</mo><msub><mi>v</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>v</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>,</mo><mo>...</mo><msub><mi>v</mi><mrow><mi>i</mi><mn>2</mn><mi>K</mi></mrow></msub><mo>]</mo></mrow><mi>T</mi></msup><mo>;</mo></mrow>]]></math><img file="FSA0000123054760000016.GIF" wi="479" he="71" /></maths>G.用与步骤E相同的方法求出待识别煤岩图像的基元响应分布<img file="FSA0000123054760000017.GIF" wi="392" he="68" />H.用<img file="FSA0000123054760000018.GIF" wi="351" he="83" />度量与已知煤和岩对象的相似性,当小于1时为煤,否则为岩,其中,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>ϵ</mi><mi>c</mi></msub><mo>=</mo><msubsup><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mn>2</mn><mi>K</mi></mrow></msubsup><mi>m</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>z</mi><mrow><mi>c</mi><mi>k</mi></mrow></msub><mo>,</mo><msub><mover><mi>z</mi><mo>‾</mo></mover><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>ϵ</mi><mi>r</mi></msub><mo>=</mo><msubsup><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mn>2</mn><mi>K</mi></mrow></msubsup><mi>m</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>z</mi><mrow><mi>r</mi><mi>k</mi></mrow></msub><mo>,</mo><msub><mover><mi>z</mi><mo>‾</mo></mover><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FSA0000123054760000019.GIF" wi="1038" he="80" /></maths>δ为给定的参数。 |