1]所述一維的低通濾波器h。和高通濾波器h1的序列長度均為奇數(shù)，并且h。序列和 hi序列均為中屯、對稱序列。
[0052] 所述一維的低通濾波器h。。、hw和高通濾波器h1。、hii的序列長度相等并且均為偶 數(shù)，hw序列為h。。序列的反轉(zhuǎn)序列，h1。序列通過對h。。序列進行反轉(zhuǎn)操作然后對偶數(shù)序號采 樣點的序列元素進行符號取反操作而得到，hii序列為h1。序列的反轉(zhuǎn)序列。
[0053] 在給定兩組序列長度相等的一維濾波器hf、h,和濾波對象X的前提下，所述雙重列 濾波包括W下步驟：
[0054] STC1.對X進行行對稱延拓得到X。,，如果X的總行數(shù)和總列數(shù)分別為RowsX和 ColsX，hf和h,的序列長度為len，那么X。的總行數(shù)和總列數(shù)分別為（RowsX+2Xlen)和 ColsX，并且X。,的第1，2，…，len行的數(shù)據(jù)內(nèi)容分別用X的第len，len-1，…，1行的數(shù)據(jù) 內(nèi)容進行填充，X。,的第len+l，len+2,…，len+RowsX行的數(shù)據(jù)內(nèi)容分別用X的第1，2,…， RowsX行的數(shù)據(jù)內(nèi)容進行填充，Xgx的第len+RowsX+1，len+RowsX+2，…，len+RowsX+len行 的數(shù)據(jù)內(nèi)容分別用X的第RowsX，RowsX-1，…，RowsX-(len-l)行的數(shù)據(jù)內(nèi)容進行填充；
[0055] STC2.抽取hf中序號為奇數(shù)的濾波系數(shù)構(gòu)成子濾波器hfO,抽取hf中序號為偶數(shù) 的濾波系數(shù)構(gòu)成子濾波器hfe，抽取h,中序號為奇數(shù)的濾波系數(shù)構(gòu)成子濾波器h,0,抽取h, 中序號為偶數(shù)的濾波系數(shù)構(gòu)成子濾波器h,e;
[0056] STC3.對步驟STC1所述的X。,進行修正，即舍棄X。的起始兩行和末尾兩行，從而 得到X'。,；
[0057] STC4.抽取步驟STC3所述X'中行號為奇數(shù)的數(shù)據(jù)構(gòu)成X'J)，抽取步驟STC3 所述X'。中行號為偶數(shù)的數(shù)據(jù)構(gòu)成X'。擊；
[0058] STC5.抽取步驟STC4所述X'。擊中行號為偶數(shù)的數(shù)據(jù)構(gòu)成X'JV，抽取步驟STC4 所述X'J)中行號為偶數(shù)的數(shù)據(jù)構(gòu)成X'wB，抽取步驟STC4所述X'山中行號為奇數(shù)的 數(shù)據(jù)構(gòu)成X' "C，抽取步驟STC4所述X'。,0中行號為奇數(shù)的數(shù)據(jù)構(gòu)成X'。,；
[0059] STC6.用步驟STC2所述的子濾波器hfO對步驟ST巧所述的X'。占進行逐列濾波， 從而得到ResB;
[0060] STC7.用步驟STC2所述的子濾波器hfe對步驟ST巧所述的X'。占進行逐列濾波， 從而得到ResD;
[0061] STC8.用步驟STC2所述的子濾波器h,o對步驟ST巧所述的X' 進行逐列濾波， 從而得到ResA;
[0062] STC9.用步驟STC2所述的子濾波器h,e對步驟ST巧所述的X' 進行逐列濾波， 從而得到ResC;
[0063] STC10.把步驟STC6所述的ResB和步驟STC7所述的ResD中行號相同且列號也相 同的元素相加得到ResBD，即ResBD=ResB+ResD;
[0064] STCll.把步驟STC8所述的ResA和步驟STC9所述的ResC中行號相同且列號也相 同的元素相加得到ResAC，即ResAC=ResA+ResC;
[00財 STC12.申請一塊總行數(shù)為（RowsAC+Rows抓）且總列數(shù)為ColsAC或Cols抓的數(shù) 據(jù)存儲區(qū)FinalRes，其中RowsAC表示步驟STC11所述ResAC的總行數(shù)，Rows抓表示步驟STC10所述Res抓的總行數(shù)，ColsAC表示步驟STC11所述ResAC的總列數(shù)，Cols抓表示步 驟STC10所述ResBD的總列數(shù)，并且滿足RowsAC=Rows抓和ColsAC=Cols抓的條件；
[0066] STC13.把一維濾波器hf和h,視為兩個向量，如果該兩個向量的內(nèi)積大于0,那么 步驟STC12所述FinalRes的第1，3,5,7，…，（RowsAC+Rows抓-1)行的數(shù)據(jù)內(nèi)容分別用步 驟STC10所述ResBD的第1，2, 3,4，…，Rows抓行的數(shù)據(jù)內(nèi)容進行填充，步驟STC12所述 FinalRes的第2,4,6,8,…，（RowsAC+RowsBD)行的數(shù)據(jù)內(nèi)容分別用步驟STC11所述ResAC 的第1，2,3,4,…，RowsAC行的數(shù)據(jù)內(nèi)容進行填充，如果該兩個向量的內(nèi)積小于或等于0,那 么步驟STC12所述FinalRes的第1，3,5,7，…，（RowsAC+Rows抓-1)行的數(shù)據(jù)內(nèi)容分別用 步驟STC11所述ResAC的第1，2,3,4,…，RowsAC行的數(shù)據(jù)內(nèi)容進行填充，步驟STC12所述 FinalRes的第2,4,6,8,…，（RowsAC+RowsBD)行的數(shù)據(jù)內(nèi)容分別用步驟STC10所述ResBD 的第1，2,3,4,…，Rows抓行的數(shù)據(jù)內(nèi)容進行填充，其中RowsAC表示步驟STCll所述ResAC 的總行數(shù)，Rows抓表示步驟STC10所述ResBD的總行數(shù)；
[0067] STC14.完成雙重列濾波操作，其最終結(jié)果輸出為步驟STC12所述的FinalRes。
[0068] 在給定兩組序列長度相等的一維濾波器hf、h,和濾波對象X的前提下，所述雙重行 濾波包括W下步驟：
[0069] STR1.對X進行列對稱延拓得到X。,，如果X的總行數(shù)和總列數(shù)分別為RowsX 和ColsX，hf和h,的序列長度為len，那么XK的總行數(shù)和總列數(shù)分別為RowsX和 (ColsX+2Xlen)，并且X。,的第1，2，…，len列的數(shù)據(jù)內(nèi)容分別用X的第len，len-1，…， 1列的數(shù)據(jù)內(nèi)容進行填充，X。,的第len+1，len+2，…，len+ColsX列的數(shù)據(jù)內(nèi)容分別用X的第1，2，…，ColsX列的數(shù)據(jù)內(nèi)容進行填充，XJ勺第len+ColsX+1，len+ColsX+2，…， len+ColsX+len列的數(shù)據(jù)內(nèi)容分別用X的第ColsX，ColsX-l，…，ColsX-(len-l)列的數(shù)據(jù) 內(nèi)容進行填充；
[0070]STR2.抽取hf中序號為奇數(shù)的濾波系數(shù)構(gòu)成子濾波器hfO,抽取hf中序號為偶數(shù) 的濾波系數(shù)構(gòu)成子濾波器hfe，抽取h,中序號為奇數(shù)的濾波系數(shù)構(gòu)成子濾波器h,0,抽取h, 中序號為偶數(shù)的濾波系數(shù)構(gòu)成子濾波器h,e;
[OCm]STR3.對步驟STRl所述的Xji行修正，即舍棄X。的起始兩列和末尾兩列，從而 得到X'。,；
[0072] STR4.抽取步驟STR3所述X'中列號為奇數(shù)的數(shù)據(jù)構(gòu)成X'J)，抽取步驟STR3 所述X'。中列號為偶數(shù)的數(shù)據(jù)構(gòu)成X'。擊；
[0073] ST貼.抽取步驟STR4所述X'。擊中列號為偶數(shù)的數(shù)據(jù)構(gòu)成X'JV，抽取步驟STR4 所述X'J)中列號為偶數(shù)的數(shù)據(jù)構(gòu)成X'。占，抽取步驟STR4所述X'山中列號為奇數(shù)的 數(shù)據(jù)構(gòu)成X' "C，抽取步驟STR4所述X'。,0中列號為奇數(shù)的數(shù)據(jù)構(gòu)成X'。,；
[0074] STR6.用步驟STR2所述的子濾波器hfO對步驟STR5所述的X'。占進行逐行濾波， 從而得到ResB;
[00巧]STR7.用步驟STR2所述的子濾波器hfe對步驟STR5所述的X'。占進行逐行濾波， 從而得到ResD;
[0076] STR8.用步驟STR2所述的子濾波器h,o對步驟STR5所述的X' 進行逐行濾波， 從而得到ResA;
[0077] STR9.用步驟STR2所述的子濾波器h,e對步驟STR5所述的X' 進行逐行濾波， 從而得到ResC;
[0078] STR10.把步驟STR6所述的ResB和步驟STR7所述的ResD中行號相同且列號也相 同的元素相加得到ResBD，即ResBD=ResB+ResD;
[0079] STRll.把步驟STR8所述的ResA和步驟STR9所述的ResC中行號相同且列號也相 同的元素相加得到ResAC，即ResAC=ResA+ResC;
[0080] STR12.申請一塊總行數(shù)為RowsAC或Rows抓且總列數(shù)為（ColsAC+Cols抓）的數(shù) 據(jù)存儲區(qū)FinalRes，其中RowsAC表示步驟STR11所述ResAC的總行數(shù)，Rows抓表示步驟 STR10所述Res抓的總行數(shù)，ColsAC表示步驟STR11所述ResAC的總列數(shù)，Cols抓表示步 驟STR10所述ResBD的總列數(shù)，并且滿足RowsAC=Rows抓和ColsAC=Cols抓的條件；
[0081]STR13.把一維濾波器hf和h,視為兩個向量，如果該兩個向量的內(nèi)積大于0,那么 步驟STR12所述FinalRes的第1，3,5,7，…，（ColsAC+Cols抓-1)列的數(shù)據(jù)內(nèi)容分別用步 驟STR10所述ResBD的第1，2, 3,4，…，Cols抓列的數(shù)據(jù)內(nèi)容進行填充，步驟STR12所述 FinalRes的第2,4,6,8,…，（ColsAC+ColsBD)列的數(shù)據(jù)內(nèi)容分別用步驟STR11所述ResAC 的第1，2,3,4,…，ColsAC列的數(shù)據(jù)內(nèi)容進行填充，如果該兩個向量的內(nèi)積小于或等于0,那 么步驟STR12所述FinalRes的第1，3,5,7，…，（ColsAC+Cols抓-1)列的數(shù)據(jù)內(nèi)容分別用 步驟STR11所述ResAC的第1，2,3,4,…，ColsAC列的數(shù)據(jù)內(nèi)容進行填充，步驟STR12所述 FinalRes的第2,4,6,8,…，（ColsAC+ColsBD)列的數(shù)據(jù)內(nèi)容分別用步驟STR10所述ResBD 的第1，2,3,4,…，Cols抓列的數(shù)據(jù)內(nèi)容進行填充，其中ColsAC表示步驟STRll所述ResAC 的總列數(shù)，Cols抓表示步驟STR10所述ResBD的總列數(shù)；
[0082] STR14.完成雙重行濾波操作，其最終結(jié)果輸出為步驟STR12所述的FinalRes。
[0083] 所述計算均方根小波系數(shù)子帶在服從廣義伽瑪分布條件下的指數(shù)形狀參數(shù)a、形 狀參數(shù)T和尺度參數(shù)A包括W下步驟：
[0084] STEP1.把T。=0. 05作為T參數(shù)的初始迭代值，進入W下步驟STEP2-STEP4的 循環(huán)迭代過程；
[0085] STEP2.把T作為自變量代入并計算
[0086]
[0087] 其中X康示均方根小波系數(shù)子帶中的第i個均方根小波系數(shù)，Num表示均方根小 波系數(shù)子帶中的均方根小波系數(shù)總數(shù)，In表示W(wǎng)自然常數(shù)e為底數(shù)的對數(shù)函數(shù)，即自然對 數(shù)；
[0088] STEP3.把T作為自變量代入并計算
[0089]
[0090] 其中Xi表示均方根小波系數(shù)子帶中的第i個均方根小波系數(shù)，Num表示均方根小 波系數(shù)子帶中的均方根小波系數(shù)總數(shù)，In表示W(wǎng)自然常數(shù)e為底數(shù)的對數(shù)函數(shù)，即自然對 數(shù)；
[00川 STEP4.計算
，若IA TI大于計算精度e，則通過公式T=T-I A T 來更新T，然后執(zhí)行步驟STEP2-STEP4的循環(huán)迭代過程，若IATI小于或等于計算精度e，則把當前T值作為最終T值并進入步驟STEP5 ;
[0092] STEP5.根據(jù)步驟STEP4獲得的最終T值，計算得到最終的a值和A值，最終的 a值，A值的計算公式分別為
[0093]
[0094]
[0095] 其中X康示均方根小波系數(shù)子帶中的第i個均方根小波系數(shù)，Num表示均方根小 波系數(shù)子帶中的均方根小波系數(shù)總數(shù)，In表示W(wǎng)自然常數(shù)e為底數(shù)的對數(shù)函數(shù)，即自然對 數(shù)。
[0096] 本發(fā)明所述6XNX3維的特征列向量的數(shù)據(jù)結(jié)構(gòu)采用[a1_1，T1_1，Ah，a1_2， T1-2，人 1-2>a1-3，T1-3，人 1-3，a1-4>T1-4> 人 1-4，a1-5，T1-5，人 1-5，a1-6，T1-6，人 1-6>a2-1，T2-1，、2-1'a2-2，T2-2，^2-2，a2-3，了2-3，^2-3，日 2-4，了2-4，^2-4，日 2-5，了3-5> ^2-5，日2-6，了2-6> 、2-6'???，aN-1，TN-1，^N-1，aN-2，TN-2，^N-2，aN-3，TN-3，^N-3，aN-4，TN-4，^N-4，aN-5，TN-5， 入,_5,aw_e，Tw_e，Aw_e]T的格式，特征列向量中的每一個a元素、T元素、A元素可W分別 表示為ai_。、Ti_。、。，其中下標i表示兩路小波變換的級數(shù)序號，i的取值為1，2,…，N， n表示均方根小波系數(shù)子帶均值和均方根小波系數(shù)子帶方差的乘積值由大到小順序排列的 排列序號，n的取值為1，2,…，6,ai_。、T。和A。分別表示煤巖樣本子圖在經(jīng)過第i級 兩路小波變換W后的6個均方根小波系數(shù)子帶中均方根小波系數(shù)子帶均值和均方根小波 系數(shù)子帶方差的乘積值由大到小順序排第n個的均方根小波系數(shù)子帶在服從廣義伽瑪分 布條件下的指數(shù)形狀參數(shù)a、形狀參數(shù)T和尺度參數(shù)A。
[0097] 所述比較特征列向量K,,與特征列向量、K.,.之間的相似度，從而判定待識別子圖 Uy所屬煤巖類別的過程包括W下步驟：
[009引S1.計算特征列向量K,,與特征列向量K,.、K,.之間的相似度，特征列向量K,,與特征 列向量之間的相似度計算公式，特征列向量K,，與特征列向量K.,.之間的相似度計算公式分 別為
[0099]
[0101] 其中*5>?於乂.)表示特征列向量^,與特征列向量^.,之間的相似腹訊*長，^；) 表示特征列向量K,，與特征列向量之間的相似度，特征列向量和K,中的下標k分別表 示煤炭樣本序號和巖石樣本序號，k的取值為1，2，…，m，N表示兩路小波變換的總級 數(shù)乂,、K,和^；中的每一個元素可W分別表示為uXp、ckp和rkp，p的取值為l，2，…， 化XNX：3)，uxp表示K,,的第P個元素，ckp表示K.,的第P個元素，rkp表示Ki.的第P個元素， 函數(shù)
在r(Z)函數(shù)表達式中，Z為自變量，t為積分變量，積分區(qū)間為[0， + -)，e表示自然常數(shù)，函蠻
，在1]^ (Z)函數(shù)表達式中，Z為自變量，r' (Z)為 r(Z)的導函數(shù)；
[0102] S2.進行相似度比較，如^
那么判定待識別 子圖11,屬于煤炭，如男
那么判定待識別子圖11,屬于巖 石。
【附圖說明】
[0103] 圖1是基于變換域廣義伽瑪分布的煤巖識別方法的基本流程圖；
[0104] 圖2是本發(fā)明所述N級兩路小波變換的基本流程圖；
[0105] 圖3是本發(fā)明所述雙重列濾波的基本流程圖；
[0106] 圖4是本發(fā)明所述雙重行濾波的基本流程圖；
[0107] 圖5是本發(fā)明所述計算均方根小波系數(shù)子帶在服從廣義伽瑪分布條件下的指數(shù) 形狀參數(shù)a、形狀參數(shù)T和尺度參數(shù)A的基本流程圖。
【具體實施方式】
[0108] 在對我國山西、河南等地主要煤種和巖種的圖像進行實驗分析的基礎上，本發(fā)明 提出了基于變換域廣義伽瑪分布的煤巖識別方法，該方法可W有效判別煤巖與巖石。
[0109] 下面結(jié)合附圖和【具體實施方式】對本發(fā)明作進一步的詳細描述。
[0110] 參照圖1，基于變換域廣