V對步驟（17)所述的低頻部分小波系 數(shù)子帶CoefL2-進(jìn)行雙重行濾波，從而得到低頻部分小波系數(shù)子帶CoefLLW，其中下標(biāo)i表 示兩路小波變換的級數(shù)序號；
[0168](23)利用步驟（19)所述的高頻部分小波系數(shù)子帶CoefHLw，分別構(gòu)造煤巖圖像在 經(jīng)過第i級兩路小波變換W后的第1個方向的第1路高頻部分小波系數(shù)子帶Co卸、煤巖圖 像在經(jīng)過第i級兩路小波變換W后的第1個方向的第2路高頻部分小波系數(shù)子帶Co或?,、 煤巖圖像在經(jīng)過第i級兩路小波變換W后的第6個方向的第1路高頻部分小波系數(shù)子帶 。€扭&、煤巖圖像在經(jīng)過第i級兩路小波變換W后的第6個方向的第2路高頻部分小波系數(shù) 子帶Cog枉。，它們中各元素的計算公式分別為
[017引上述各式中U和V分別表示Coe尤,、扣游。的行號和列號，Coe尤,(w)(v)、Co訴,(M)(v)、Co訴6(")(v)和仿與6(")(v)分別表示打，典,、C。游、C。訴。 和Coe把6中行號為U、列號為V的元素，CoefH^_i(2u-l) (2v-l)表示Coefa。-沖行號為化-1、 列號為2v-l的元素，CoefH^_i(2u) (2v)表示CoefH^_沖行號為2u、列號為2v的元素， CoefHL2-i(2u-l) (2v)表示CoefHL2-沖行號為化-1、列號為 2v的元素，CoefHL2-i(2u) (2v-l) 表示Coefaiw中行號為化、列號為2v-l的元素，u和V的最小取值為1，u的最大取值為 行數(shù)的一半，V的最大取值為Coef 列數(shù)的一半；
[0174] (24)利用步驟（20)所述的高頻部分小波系數(shù)子帶Coe片HW，分別構(gòu)造煤巖圖像在 經(jīng)過第i級兩路小波變換W后的第3個方向的第1路高頻部分小波系數(shù)子帶、煤巖 圖像在經(jīng)過第i級兩路小波變換W后的第3個方向的第2路高頻部分小波系數(shù)子帶Coe記，、 煤巖圖像在經(jīng)過第i級兩路小波變換W后的第4個方向的第1路高頻部分小波系數(shù)子帶 Co訴4、煤巖圖像在經(jīng)過第i級兩路小波變換W后的第4個方向的第2路高頻部分小波系數(shù) 子帶Coe充4，它們中各元素的計算公式分別為
[017引上述各式中U和V分別表示Coe記;、Coe記;、Coe記4、Coe記4的行號 和列號，Coe記sWW、Co說3(")(小Coe尤4(m)(v)和Coe/ii(M)(v)分別表示 Coe記3、Coe居、Co於4和Coe后中行號為U、列號為V的元素，CoefLH2-i腳-1)仍-1)表 示Coefm2_沖行號為化-1、列號為2v-l的元素，Coef\H2_i(2u) (2v)表示006片112_沖行號為 2u、列號為2v的元素，Coefu2_i(2u-l) (2v)表示Coefu2_i中行號為化-1、列號為2v的元素， Coefu2_i(2u) (2v-l)表示Coefm2_i中行號為化、列號為2v-l的元素，U和V的最小取值為 1，U的最大取值為006片112_1總行數(shù)的一半，V的最大取值為Coef\H2_^列數(shù)的一半；
[0180] (25)利用步驟（21)所述的高頻部分小波系數(shù)子帶CoefsHw，分別構(gòu)造煤巖圖像在 經(jīng)過第i級兩路小波變換W后的第2個方向的第1路高頻部分小波系數(shù)子帶Coe/i;、煤巖 圖像在經(jīng)過第i級兩路小波變換W后的第2個方向的第2路高頻部分小波系數(shù)子帶Co從?；、 煤巖圖像在經(jīng)過第i級兩路小波變換W后的第5個方向的第1路高頻部分小波系數(shù)子帶 Co蛛li、煤巖圖像在經(jīng)過第i級兩路小波變換W后的第5個方向的第2路高頻部分小波系 數(shù)子帶Coey占，它們中各元素的計算公式分別為
[01財上述各式中U和V分別表示Co說2、Coe后、Co與5、Coe店的行號 和列號，Coe技(M)(v)、Coe后(w)(v)、Coe技(w)(v)和Coe記(M)(v)分別表示 Coe記2、Coe托2、Co典5和Coe店中行號為U、列號為V的元素，CoefHH2-i腳-1)仍-1)表 示CoefHH2_i中行號為化-1、列號為2v-l的元素，CoefHH2_i(2u) (2v)表示(：06'胃_1中行號為 2u、列號為2v的元素，CoefHH2_i(2u-l) (2v)表示CoefHH2_i中行號為化-1、列號為2v的元素， C0efHH2_i(2u) (2v-l)表示CoefHH2_i中行號為化、列號為2v-l的元素，U和V的最小取值為 1，u的最大取值為006'胃_1總行數(shù)的一半，V的最大取值為Coef胃_1總列數(shù)的一半；
[0186] (26)釋放臨時數(shù)據(jù)存儲區(qū)TempCoef所占用的內(nèi)存空間，重新動態(tài)申請與步驟 (22)所述Co噸。_斯數(shù)據(jù)結(jié)構(gòu)完全一致的臨時數(shù)據(jù)存儲區(qū)TempCoef，并將其數(shù)據(jù)內(nèi)容初始 化為Coe片。的數(shù)據(jù)內(nèi)容；
[0187] (27)兩路小波變換的級數(shù)序號自增1，即i=i+1 ;
[018引 （28)轉(zhuǎn)到步驟（14)，進(jìn)行下一級兩路小波變換；
[0189] (29)釋放臨時數(shù)據(jù)存儲區(qū)TempCoef所占用的內(nèi)存空間，完成N級兩路小波變換。
[0190]所述一維的低通濾波器h。和高通濾波器h1的序列長度均為奇數(shù)，并且h。序列和 hi序列均為中屯、對稱序列。
[0191] 所述一維的低通濾波器h。。、hw和高通濾波器h1。、h。的序列長度相等并且均為偶 數(shù)，hw序列為h。。序列的反轉(zhuǎn)序列，h1。序列通過對h。。序列進(jìn)行反轉(zhuǎn)操作然后對偶數(shù)序號采 樣點的序列元素進(jìn)行符號取反操作而得到，hii序列為h1。序列的反轉(zhuǎn)序列。
[019引如圖3所示，在給定兩組序列長度相等的一維濾波器hf、h郝濾波對象X的前提 下，所述雙重列濾波包括W下步驟：
[0193] STC1.對X進(jìn)行行對稱延拓得到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)容進(jìn)行填充，X。,的第len+l，len+2,…，len+RowsX行的數(shù)據(jù)內(nèi)容分別用X的第1，2,…， RowsX行的數(shù)據(jù)內(nèi)容進(jìn)行填充，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)容進(jìn)行填充；
[0194] 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;
[0195] STC3.對步驟STC1所述的X。,進(jìn)行修正，即舍棄X的起始兩行和末尾兩行，從而 得到X'。,；
[0196] STC4.抽取步驟STC3所述X'中行號為奇數(shù)的數(shù)據(jù)構(gòu)成X'J)，抽取步驟STC3 所述X'。中行號為偶數(shù)的數(shù)據(jù)構(gòu)成X'。擊；
[0197] STC5.抽取步驟STC4所述X'。擊中行號為偶數(shù)的數(shù)據(jù)構(gòu)成X'JV，抽取步驟STC4 所述X'J)中行號為偶數(shù)的數(shù)據(jù)構(gòu)成X'抽取步驟STC4所述X'山中行號為奇數(shù)的 數(shù)據(jù)構(gòu)成X' "C，抽取步驟STC4所述X'。,0中行號為奇數(shù)的數(shù)據(jù)構(gòu)成X'。,；
[0198] STC6.用步驟STC2所述的子濾波器hfO對步驟ST巧所述的X'。占進(jìn)行逐列濾波， 從而得到ResB;
[0199] STC7.用步驟STC2所述的子濾波器hfe對步驟ST巧所述的X'。占進(jìn)行逐列濾波， 從而得到ResD;
[0200] STC8.用步驟STC2所述的子濾波器h,o對步驟ST巧所述的X' 進(jìn)行逐列濾波， 從而得到ResA;
[0201] STC9.用步驟STC2所述的子濾波器h,e對步驟ST巧所述的X'c,C進(jìn)行逐列濾波， 從而得到ResC;
[0202] STC10.把步驟STC6所述的ResB和步驟STC7所述的ResD中行號相同且列號也相 同的元素相加得到ResBD，即ResBD=ResB+ResD;
[0203] STCll.把步驟STC8所述的ResA和步驟STC9所述的ResC中行號相同且列號也相 同的元素相加得到ResAC，即ResAC=ResA+ResC;
[0204] STC12.申請一塊總行數(shù)為（RowsAC+RowsBD)且總列數(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抓的條件； [020引 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)容進(jìn)行填充，步驟STC12所述 FinalRes的第2,4,6,8,…，（RowsAC+RowsBD)行的數(shù)據(jù)內(nèi)容分別用步驟STC11所述ResAC 的第1，2,3,4,…，RowsAC行的數(shù)據(jù)內(nèi)容進(jìn)行填充，如果該兩個向量的內(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)容進(jìn)行填充，步驟STC12所述 FinalRes的第2,4,6,8,…，（RowsAC+RowsBD)行的數(shù)據(jù)內(nèi)容分別用步驟STC10所述ResBD 的第1，2,3,4,…，Rows抓行的數(shù)據(jù)內(nèi)容進(jìn)行填充，其中RowsAC表示步驟STCll所述ResAC 的總行數(shù)，Rows抓表示步驟STC10所述ResBD的總行數(shù)；
[0206] STC14.完成雙重列濾波操作，其最終結(jié)果輸出為步驟STC12所述的FinalRes。
[0207] 如圖4所示，在給定兩組序列長度相等的一維濾波器hf、h,和濾波對象X的前提 下，所述雙重行濾波包括W下步驟：
[020引 STR1.對X進(jìn)行列對稱延拓得到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)容進(jìn)行填充，X。,的第len+1，len+2，…，len+ColsX列的數(shù)據(jù)內(nèi)容分別用X的第1，2，…，ColsX列的數(shù)據(jù)內(nèi)容進(jìn)行填充，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)容進(jìn)行填充；
[0209] 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;
[0210] STR3.對步驟STR1所述的Xji行修正，即舍棄X。的起始兩列和末尾兩列，從而 得到X'。,；
[0211] STR4.抽取步驟STR3所述X'中列號為奇數(shù)的數(shù)據(jù)構(gòu)成X'J)，抽取步驟STR3 所述X'。中列號為偶數(shù)的數(shù)據(jù)構(gòu)成X'。擊；
[0212] 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'。,；
[0213] STR6.用步驟STR2所述的子濾波器hfO對步驟STR5所述的X'。占進(jìn)行逐行濾波， 從而得到ResB;
[0214] STR7.用步驟STR2所述的子濾波器hfe對步驟STR5所述的X'。占進(jìn)行逐行濾波， 從而得到ResD;
[0215] STR8.用步驟STR2所述的子濾波器h,o對步驟STR5所述的X' 進(jìn)行逐行濾波， 從而得到ResA;
[0216] STR9.用步驟STR2所述的子濾波器h,e對步驟STR5所述的X'c,C進(jìn)行逐行濾波， 從而得到ResC;
[0217] STR10.把步驟STR6所述的ResB和步驟STR7所述的ResD中行號相同且列號也相 同的元素相加得到ResBD，即ResBD=ResB+ResD;
[021引STR11.把步驟STR8所述的ResA和步驟STR9所述的ResC中行號相同且列號也相 同的元素相加得到ResAC，即ResAC=ResA+ResC;
[0219] 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抓的條件；
[0220] 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)容進(jìn)行填充，步驟STR12所述 FinalRes的第2,4,6,8,…，（ColsAC+ColsBD)列的數(shù)據(jù)內(nèi)容分別用步驟STR11所述ResAC 的第1，2,3,4,…，ColsAC列的數(shù)據(jù)內(nèi)容進(jìn)行填充，如果該兩個向量的內(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)容進(jìn)行填充，步驟STR12所述 FinalRes的第2,4,6,8,…，（ColsAC+ColsBD)列的數(shù)據(jù)內(nèi)容分別用步驟STR10所述ResBD 的第1，2,3,4,…，Cols抓列的數(shù)據(jù)內(nèi)容進(jìn)行填充，其中ColsAC表示步驟STR11所述ResAC 的總列數(shù)，Cols抓表示步驟STR10所述ResBD的總列數(shù)；
[0221] STR14.完成雙重行濾波操作，其最終結(jié)果輸出為步驟STR12所述的FinalRes。
[0222] 如圖5所示，所述計算均方根小波系數(shù)子帶在服從廣義伽瑪分布條件下的指數(shù)形 狀參數(shù)a、形狀參數(shù)T和尺度參數(shù)A包括W下步驟：
[0223] STEP1.把T。= 0. 05作為T參數(shù)的初始迭代值，進(jìn)入W下步驟STEP2-STEP4的 循環(huán)迭代過程；
[02