何飛佳， 李慶武,2， 韓 輝， 張建清， 譚顯江

(1.河海大學(xué) 物聯(lián)網(wǎng)工程學(xué)院，江蘇 常州 213022；2.常州市傳感網(wǎng)與環(huán)境感知重點(diǎn)實(shí)驗(yàn)室，江蘇 常州 213022;3.長(zhǎng)江地球物理探測(cè)(武漢)有限公司，湖北 武漢 430010)

?

基于同態(tài)濾波與Curvelet變換的鉆孔圖像自適應(yīng)增強(qiáng)*

何飛佳1， 李慶武1,2， 韓 輝1， 張建清3， 譚顯江3

(1.河海大學(xué) 物聯(lián)網(wǎng)工程學(xué)院，江蘇 常州 213022；2.常州市傳感網(wǎng)與環(huán)境感知重點(diǎn)實(shí)驗(yàn)室，江蘇 常州 213022;3.長(zhǎng)江地球物理探測(cè)(武漢)有限公司，湖北 武漢 430010)

針對(duì)巖石數(shù)字鉆孔圖像存在的光照不均、圖像中巖石表面邊緣細(xì)節(jié)模糊等情況，提出了一種鉆孔圖像自適應(yīng)增強(qiáng)算法。對(duì)原圖進(jìn)行同態(tài)濾波；使用Curvelet變換分解原圖與濾波后的圖像，對(duì)兩者的低頻子帶使用系數(shù)直方圖匹配算法，將前者與后者的直方圖進(jìn)行匹配，改善光照不均的影響；對(duì)原圖的高頻子帶使用自適應(yīng)的閾值進(jìn)行濾波，同時(shí)利用自適應(yīng)增強(qiáng)函數(shù)進(jìn)行增強(qiáng)；使用Curvelet反變換重構(gòu)得到增強(qiáng)后的圖像。實(shí)驗(yàn)結(jié)果表明：算法可以有效地改善鉆孔圖像光照不均的問(wèn)題，增強(qiáng)圖像中物體的邊緣信息，在主觀視覺(jué)效果和圖像客觀評(píng)價(jià)指標(biāo)上相對(duì)于其他算法均有一定優(yōu)勢(shì)。

鉆孔圖像； 圖像增強(qiáng)； Curvelet變換； 系數(shù)直方圖匹配； 自適應(yīng)增強(qiáng)

0 引 言

鉆孔電視技術(shù)被廣泛應(yīng)用于工程地質(zhì)、巖土工程、土木工程等多個(gè)領(lǐng)域[1]。在實(shí)際應(yīng)用中，由于拍攝條件限制以及巖石表面凹凸不平等原因，巖石鉆孔圖像通常存在光照不均、圖像中物體邊緣細(xì)節(jié)模糊等情況，圖像中光照強(qiáng)的高光區(qū)細(xì)節(jié)信息失真，而光照弱的暗光區(qū)細(xì)節(jié)信息較難顯現(xiàn)。

為了消除光照不均的影響，同時(shí)提高圖像局部細(xì)節(jié)信息，需要對(duì)圖像進(jìn)行增強(qiáng)處理。圖像增強(qiáng)方法主要分為基于空域和基于頻域兩大類?？沼蛟鰪?qiáng)算法雖能提高圖像整體對(duì)比度，但會(huì)降低圖像高光區(qū)和暗光區(qū)的灰度值區(qū)分度，造成兩類區(qū)域細(xì)節(jié)信息丟失。頻域圖像增強(qiáng)方法的主要思想是將圖像分解為低頻子帶與高頻子帶，通過(guò)對(duì)含有大量圖像基本信息的低頻子帶增強(qiáng)以提高圖像對(duì)比度，同時(shí)抑制高頻子帶中含有的噪聲并且增強(qiáng)圖像邊緣細(xì)節(jié)信息。Curvelet變換屬于頻域算法[2～7]，由Candes E J和Donoh D等人在小波變換的基礎(chǔ)上提出的一種新的多尺度幾何分析方法，具有各向異性的特點(diǎn)，對(duì)圖像幾何特征有更強(qiáng)的表達(dá)能力，能夠很好地提取圖像邊緣特征。

同態(tài)濾波對(duì)改善光照不均的圖像對(duì)比度具有良好的效果，但同時(shí)也會(huì)模糊圖像紋理細(xì)節(jié)。經(jīng)Curvelet分解后，圖像概貌主要反映在Curvelet低頻系數(shù)上，圖像紋理細(xì)節(jié)主要反映在Curvelet高頻系數(shù)上。鑒于此，本文提出了一種基于同態(tài)濾波和Curvelet變換系數(shù)直方圖匹配的鉆孔圖像自適應(yīng)增強(qiáng)方法，以同態(tài)濾波處理后圖像的系數(shù)直方圖為基準(zhǔn)，將原圖分解所得的低頻子帶系數(shù)直方圖匹配于前者，改善原圖中光照不均問(wèn)題,同時(shí)保留了紋理細(xì)節(jié)；對(duì)原圖分解所得的高頻子帶使用自適應(yīng)增強(qiáng)算法，增強(qiáng)圖像邊緣紋理信息，最終實(shí)現(xiàn)了巖石鉆孔圖像的自適應(yīng)增強(qiáng)。

1 Curvelet變換

圖1 Curvelet變換頻帶劃分示意圖

本文采用基于“Wrapping”頻域特殊采樣的卷繞規(guī)則方法，即,形式的快速離散Curvelet變換算法[9]，其步驟如下：

1)將圖像看作笛卡爾坐標(biāo)系下的二維函數(shù)f，對(duì)f進(jìn)行二維快速傅里葉變換(FFT)計(jì)算得到；

2 基于同態(tài)濾波和Curvelet變換系數(shù)直方圖匹配增強(qiáng)

2.1 同態(tài)濾波

當(dāng)用f(x,y)二維函數(shù)形式表示空間域圖像時(shí)，可以由2個(gè)分量表征：入射光源總量i(x,y)和物體反射光的總量r(x,y)，即

f(x,y)=i(x,y)·r(x,y)

(1)

考慮到所獲得的鉆孔圖像可看作外部光源的入射光與巖石表面的反射光兩者的乘積，由此引入同態(tài)濾波來(lái)改善圖像中光照不均的問(wèn)題，通過(guò)適當(dāng)降低入射量函數(shù)i(x,y)的影響，同時(shí)增強(qiáng)反射率函數(shù)r(x,y)的頻譜成分達(dá)到增強(qiáng)巖石鉆孔圖像的目的。具體流程如圖2所示。

圖2 同態(tài)濾波流程

其中，影響高頻與低頻分量的高通濾波器設(shè)計(jì)如下

(2)

式中rH,rL分別為高通截止頻率和低通截止頻率；d(u,v)=[(u-u0)2+(v-v0)2]1/2為頻率(u,v)到濾波中心(u0,v0)的距離；C和D共同決定了高斯濾波器的頻率寬度。

2.2 低頻子帶系數(shù)直方圖匹配

設(shè)c為歸一化至區(qū)間[0,1]的原圖Curvelet低頻子帶C0中的系數(shù)，d為歸一化至區(qū)間[0,1]的同態(tài)濾波后圖像的Curvelet低頻子帶D0中的系數(shù)，ck和dq分別為系數(shù)c和系數(shù)d經(jīng)量化級(jí)數(shù)為L(zhǎng)的均勻量化后的系數(shù)。pc(ci)為ck中第i個(gè)系數(shù)個(gè)數(shù)ni占ck中系數(shù)總數(shù)n的比例，pd(di)為期望輸出的系數(shù)直方圖所具有的概率密度函數(shù)。對(duì)ck和dq均進(jìn)行一次均衡化處理，得到以下2個(gè)累積分布函數(shù)，即

(3)

(4)

根據(jù)直方圖匹配的原理[11]，有G(dq)=T(ck)，因此,dq必須滿足條件

dq=G-1(sk)=G-1[T(ck)]

(5)

即將原圖對(duì)應(yīng)低頻子帶中的每個(gè)均衡后的系數(shù)值sk映射為直方圖匹配后的低頻子帶中的相應(yīng)dq值。經(jīng)系數(shù)反歸一化后，最終得到直方圖匹配后的低頻子帶結(jié)果。

2.3 高頻子帶系數(shù)的自適應(yīng)增強(qiáng)

經(jīng)Curvelet分解后，高頻子帶部分包含大量圖像邊緣信息與噪聲，如果直接對(duì)其進(jìn)行增強(qiáng)，那么噪聲也隨之增強(qiáng)。因此，這里需要設(shè)定一個(gè)閾值 ，將大于閾值的高頻子帶系數(shù)視為邊緣細(xì)節(jié)信號(hào)，并進(jìn)行增強(qiáng)處理；將小于閾值的高頻子帶系數(shù)視為噪聲，并進(jìn)行衰減處理，即可以實(shí)現(xiàn)對(duì)數(shù)字鉆孔圖像邊緣細(xì)節(jié)信息進(jìn)行增強(qiáng)，提高清晰度，同時(shí)具有降低噪聲的作用。

在Curvelet變換系數(shù)中大部分噪聲處于最細(xì)子帶。在最細(xì)子帶中采用經(jīng)典中值估計(jì)公式σ=median[abs(Ct)]/0.674 5對(duì)方差σ進(jìn)行估計(jì)。其中，Ct為Curvelet分解后的最細(xì)子帶系數(shù)。閾值設(shè)定為

(6)

式中σj,l為第j尺度，第l方向的噪聲方差，其值可由蒙特—卡洛估計(jì)法得到。

利用如下自適應(yīng)增強(qiáng)函數(shù)[12]進(jìn)行相應(yīng)的系數(shù)調(diào)整

(7)

2.4 算法實(shí)現(xiàn)步驟

數(shù)字鉆孔圖像自適應(yīng)增強(qiáng)算法步驟:

1)獲取鉆孔圖像I，并進(jìn)行同態(tài)濾波，得到圖像J。

2)對(duì)I進(jìn)行Curvelet變換，分解得到低頻子帶系數(shù)矩陣C0，以及其余高頻子帶系數(shù)Cj,l計(jì)算低頻子帶系數(shù)矩陣C0的系數(shù)直方圖HC0以及Cj,l的系數(shù)直方圖HCj,l。

3)對(duì)J進(jìn)行Curvelet變換，分解得到低頻子帶系數(shù)矩陣D0，以及其余高頻子帶系數(shù)矩陣Dj,l。然后計(jì)算D0的系數(shù)直方圖HD0以及Dj,l的系數(shù)直方圖HDj,l。

4)低頻子帶系數(shù)直方圖匹配：以圖像J為基準(zhǔn)，將步驟(2)中得到的C0的系數(shù)直方圖HC0匹配到步驟(3)中得到的D0的系數(shù)直方圖HD0上。

5)高頻子帶系數(shù)自適應(yīng)去噪增強(qiáng)：估計(jì)步驟(2)得到的高頻最細(xì)子帶系數(shù)矩陣C各方向的噪聲水平，按照式(6)設(shè)定閾值T，將Cj,l值小于閾值T的部分濾除，大于閾值T的部分利用式(7)對(duì)其增強(qiáng)得到系數(shù)矩陣j,l。

6) 對(duì)所有子帶系數(shù)進(jìn)行Curvelet逆變換,得到增強(qiáng)后的圖像。

3 實(shí)驗(yàn)結(jié)果與分析

為驗(yàn)證算法針對(duì)數(shù)字鉆孔圖像的有效性，利用本文提出的基于Curvelet變換的鉆孔圖像自適應(yīng)增強(qiáng)方法分別與雙向直方圖均衡[12]、同態(tài)濾波算法、非下采樣Contourlet變換(NSCT)域自適應(yīng)增強(qiáng)算法、Curvelet域增強(qiáng)算法等4種算法進(jìn)行比較。其中圖3(a),圖4(a)分別為某地區(qū)1號(hào)、2號(hào)場(chǎng)地的鉆孔圖像。下列鉆孔圖像的分辨率均為944 pixel×944 pixel。

實(shí)驗(yàn)環(huán)境：Intel Core i5—2450 CPU，8 GB內(nèi)存，64位win7操作系統(tǒng)， Matlab8.1.0(R2013a)。實(shí)驗(yàn)參數(shù)設(shè)置如下：同態(tài)濾波增強(qiáng)算法利用式(2)計(jì)算，取rH=4，rL=1.02，C=1，D=2。NSCT域自適應(yīng)增強(qiáng)算法利用文獻(xiàn)[4]給出的算法計(jì)算。Curvelet增強(qiáng)算法利用文獻(xiàn)[7]給出的算法計(jì)算，其中,低頻子帶增強(qiáng)模塊取a=2,b=4,高頻子帶增強(qiáng)模塊取b=0.13,c=10,d=0。本文算法中的同態(tài)濾波參數(shù)值設(shè)置同上。

圖3 1號(hào)場(chǎng)地鉆孔圖像及5種增強(qiáng)結(jié)果

圖4 2號(hào)場(chǎng)地鉆孔圖像及5種增強(qiáng)結(jié)果

圖3(a)、圖4(a)分別為1號(hào)、2號(hào)場(chǎng)地鉆孔圖像，在“凹槽”內(nèi)均存在光照不足的問(wèn)題，細(xì)節(jié)紋理不夠明顯且圖像整體光照度偏暗對(duì)比度低。圖3(b)，(d)，(e)和圖4(b)，(d)，(e)的對(duì)比度較原圖有所提高，但是“凹槽”內(nèi)紋理細(xì)節(jié)依舊無(wú)法辨識(shí)，例如圖3(b),(d),(e)左下角以及圖4左上角區(qū)域。圖3(c)、圖4(c)為經(jīng)過(guò)同態(tài)濾波處理后的效果圖，該方法雖改善了光照不均的問(wèn)題，然而圖像細(xì)部分清晰度沒(méi)有顯著提升。圖3(f)、圖4(f)為本文算法增強(qiáng)結(jié)果，增強(qiáng)后的圖像較為明亮，對(duì)比度較高，無(wú)論是外表面還是“凹槽”內(nèi)紋細(xì)節(jié)均清晰可見(jiàn)。

為進(jìn)一步說(shuō)明實(shí)驗(yàn)效果，對(duì)實(shí)驗(yàn)結(jié)果圖像進(jìn)行客觀評(píng)價(jià)。本文選取平均梯度、平均梯度指標(biāo)、EMEE、對(duì)比度作為圖像增強(qiáng)質(zhì)量的評(píng)價(jià)標(biāo)準(zhǔn)，另外選用運(yùn)行時(shí)間作為算法運(yùn)行效率的評(píng)價(jià)標(biāo)準(zhǔn)，分別從增強(qiáng)效果和處理時(shí)間2個(gè)維度比較各算法的性能。平均梯度的計(jì)算公式

(8)

式中Δxf(x,y)=f(x+1,y)-f(x,y);Δyf(x,y)=f(x,y+1)-f(x,y)；M和N分別為圖像的高度和寬度。

平均梯度指標(biāo)是用相鄰點(diǎn)的差分計(jì)算的梯度值，其值越大，圖像細(xì)節(jié)越多，清晰度越高。計(jì)算EMEE時(shí)，首先將圖像f(i,j)劃分為k1×k2個(gè)相同子塊。平均梯度指標(biāo)的計(jì)算公式

(9)

式中 fmax i,j和fmin i,j分別為子塊(i,j)中的最大值和最小值;一般取eps=0.000 1，以避免分母過(guò)小，子塊的大小取為8×8。

EMEE是基于熵的圖像對(duì)比度增強(qiáng)測(cè)量指標(biāo)，其值往往代表圖像一副圖像所包含的信息量，其值越大說(shuō)明圖像對(duì)比度越高，圖像清晰度越高，質(zhì)量越好。對(duì)比度計(jì)算公式為

(10)

對(duì)比度是衡量一幅圖像增強(qiáng)效果的重要指標(biāo)，對(duì)比度值越大，增強(qiáng)后的圖像視覺(jué)效果越好。

表1中的數(shù)據(jù)表明，由平均梯度和EMEE指標(biāo)來(lái)看，本文算法對(duì)圖3的增強(qiáng)效果優(yōu)于其他4種算法；對(duì)圖4的增強(qiáng)效果優(yōu)于同態(tài)濾波增強(qiáng)算法、NSCT增強(qiáng)算法、Curvelet增強(qiáng)算法。說(shuō)明本文算法對(duì)圖像中巖石表面邊緣細(xì)節(jié)的增強(qiáng)較明顯。從對(duì)比度指標(biāo)來(lái)看，本文算法的增強(qiáng)效果優(yōu)于同態(tài)濾波增強(qiáng)算法、Curvelet增強(qiáng)算法，略小于NSCT增強(qiáng)算法以及雙向直方圖增強(qiáng)算法。但是NSCT增強(qiáng)算法耗時(shí)較長(zhǎng)，不適用于對(duì)實(shí)時(shí)性要求較高的應(yīng)用場(chǎng)景，而雙向直方圖則存在過(guò)增強(qiáng)的問(wèn)題，原本“凹槽”內(nèi)的紋理細(xì)節(jié)信息丟失更加嚴(yán)重。因而，本文算法無(wú)論是從增強(qiáng)效果還是算法效率來(lái)看均存在一定優(yōu)勢(shì)。

表1 5種增強(qiáng)算法實(shí)驗(yàn)結(jié)果 s

4 結(jié) 論

提出了一種基于同態(tài)濾波方法和Curvelet變換系數(shù)直方圖匹配的鉆孔圖像自適應(yīng)增強(qiáng)算法。在Curvelet低頻子帶部分，將原圖系數(shù)直方圖匹配到經(jīng)同態(tài)濾波后的系數(shù)直方圖上，使前后兩者系數(shù)直方圖形狀近似，達(dá)到改善圖像光照不均的目的；Curvelet高頻子帶部分使用自適應(yīng)增強(qiáng)算法提高圖像邊緣紋理細(xì)節(jié)。實(shí)驗(yàn)結(jié)果表明：算法可以有效地增強(qiáng)鉆孔圖像中物體的邊緣細(xì)節(jié)，在主觀視覺(jué)效果和圖像客觀評(píng)價(jià)指標(biāo)上相對(duì)于其他算法均有一定優(yōu)勢(shì)。

[1] 王川嬰,Law K Tim.鉆孔攝像技術(shù)的發(fā)展與現(xiàn)狀[C]∥全國(guó)地面巖石工程學(xué)術(shù)會(huì)議暨中南地區(qū)巖石力學(xué)與工程學(xué)術(shù)會(huì)議,2005.

[2] Demirel H,Ozcinar C,Anbarjafari G.Satellite image contrast enhancement using discrete wavelet transform and singular value decomposition[J].IEEE Geoscience & Remote Sensing Letters,2010,7(2):333-337.

[3] Demirel H,Anbarjafari G.Image resolution enhancement by using discrete and stationary wavelet decomposition[J].IEEE Transactions on Image Processing,2011,20(5):1458-1460.

[4] 周 妍,李慶武,霍冠英.基于非下采樣Contourlet變換系數(shù)直方圖匹配的自適應(yīng)圖像增強(qiáng)[J].光學(xué)精密工程,2014,22(8):2214-2222.

[5] 沙宇恒,劉 芳,焦李成.基于非下采樣Contourlet變換的SAR圖像增強(qiáng)[J].電子與信息學(xué)報(bào),2009,31(7):1716-1721.

[6] 陳學(xué)堂,張 元,楊坤平,等.基于第二代曲波變換算法的檢測(cè)圖像增強(qiáng)[J].傳感器與微系統(tǒng),2008,27(12):8-10.

[7] 盛惠興,孟凡玲,李慶武,等.Curvelet變換域側(cè)掃聲納圖像增強(qiáng)算法[J].海洋測(cè)繪,2012,32(1):8-10.

[8] Candes E J,Donoho D L.New tight frames of curvelets and optimal representations of objects with piecewise C2 singularities[J].Communications on Pure and Applied Mathematics,2004,57(2):219-266.

[9] Candes E J,Demanet L,Donoho D L,et al.Fast discrete curvelet transforms[J].Multiscale Modeling & Simulation,2006,5(3):861-899.

[10] 岡薩雷斯.數(shù)字圖像處理[M].北京：電子工業(yè)出版社,2007.

[11] Soyel H,Mcowan P W.Automatic image enhancement using intrinsic geometrical information[J].Electronics Letters,2012,48(15):917-919.

[12] 陳 錢(qián),柏連發(fā),張保民.紅外圖像直方圖雙向均衡技術(shù)研究[J].紅外與毫米波學(xué)報(bào),2003,22(6):428-430.

Borehole image adaptive enhancement based on homomorphic filtering and Curvelet transform*

HE Fei-jia1, LI Qing-wu1,2, HAN Hui1, ZHANG Jian-qing3, TAN Xian-jiang3

(1.College of IOT Engineering,HoHai University,Changzhou 213022,China;2.Changzhou Key Laboratory of Sensor Networks and Environmental Sensing,Changzhou 213022,China;3.Changjiang Geophysical Exploration (Wuhan) Co Ltd,Wuhan 430010,China)

To solve problems such as uneven illumination and fuzzy details of rock surface edge are shown in digital borehole images,an adaptive image enhancement algorithm is proposed.Firstly,homomorphic filtering is used to original image.Then,the original image and the filtered image are decomposed by Curvelet transform while the coefficient histogram matching algorithm is used for the low frequency subband of the two.The histogram of the former is matched with the latter to improve effect of illumination unevenness.Next,the high-frequency subband is filtered by an adaptive threshold and is enhanced with an adaptive enhancement function.Finally,the enhanced image is reconstructed by Curvelet inverse transform.The experimental results show that the illumination unevenness is improved effectively and the edge information is enhanced obviously.Certain advantages compared with other algorithms are shown in both subjective and objective evaluation indices.

borehole image; image enhancement; Curvelet transform; coefficient histogram matching; adaptive enhancement

10.13873/J.1000—9787(2017)08—0145—04

2016—12—09

國(guó)家自然科學(xué)基金資助項(xiàng)目 (41306089); 江蘇省重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目(BE2016056); 常州市科技支撐計(jì)劃資助項(xiàng)目(CE20150068)

TP 391.41

A

1000—9787(2017)08—0145—04

何飛佳(1992-)，男，碩士研究生，研究方向?yàn)閿?shù)字圖像處理。

李慶武(1964-)，男，通訊作者，教授，博士生導(dǎo)師，研究方向?yàn)橹悄芨兄c圖像處理，E—mail:li_qingwu@163.com。