華梓錚, 華澤璽

（西南交通大學(xué) 電氣工程學(xué)院,四川 成都610031）

邊緣存在于圖像的不規(guī)則結(jié)構(gòu)和不平穩(wěn)情形中,包含了圖像的大部分信息.邊緣檢測(cè)是圖像處理和分析的關(guān)鍵步驟,對(duì)后續(xù)高層次的特征描述、匹配和識(shí)別等有著重大的影響,為圖像去噪等后續(xù)處理提供了依據(jù)［1］.邊緣檢測(cè)技術(shù)是 L.G.Roberts［2］最早開(kāi)始系統(tǒng)研究.如今,常用的邊緣檢測(cè)算法有Canny 算子［3-5］、Sobel 算子［6-7］、模板匹配法、變換域法［8-12］以及改進(jìn)的 Pal算法［13］等.由于邊緣檢測(cè)的固有復(fù)雜性,這些方法在抗噪性能和邊緣定位方面仍然不能完全適應(yīng)復(fù)雜的實(shí)際應(yīng)用需求.

為檢測(cè)到豐富的邊緣信息,本文提出基于非下采樣輪廓波變換（NSCT）［14］的自適應(yīng)閾值圖像邊緣檢測(cè)算法.NSCT是在輪廓波變換（CT）［15］的基礎(chǔ)上提出的,解決了傳統(tǒng)CT不具備平移不變性的缺點(diǎn).本文首先介紹了NSCT的構(gòu)造方法及其結(jié)構(gòu)特點(diǎn),根據(jù)NSCT能夠很好地提取圖像的高頻信息等特點(diǎn),將NSCT應(yīng)用于邊緣檢測(cè),即對(duì)利用NSCT提取的高頻圖像和低頻圖像分別處理得到高頻邊緣信息和低頻邊緣信息,利用低頻邊緣信息對(duì)高頻邊緣信息進(jìn)行修正,并去除修正后邊緣信息中孤立噪聲點(diǎn)以得到最終的邊緣圖像,并通過(guò)實(shí)驗(yàn)驗(yàn)證提出方法的有效性.

1 NSCT原理

CT主要是由方向?yàn)V波器組和拉普拉斯塔式濾波器組構(gòu)造而成,但是由于CT中的方向?yàn)V波器和拉普拉斯分解都有下采樣的操作,使得CT不具有平移不變性.因此文獻(xiàn)［14］給出了其下采樣的形式,即利用非下采樣分解和非下采樣濾波器組構(gòu)造出NSCT,如圖1所示.

文獻(xiàn)［14］利用àtrous算法實(shí)現(xiàn)并生成一種靈活的多尺度、多方向且平移不變的圖像分解算法,具有快速實(shí)現(xiàn)算法.算法的關(guān)鍵是不可分兩通道非下采樣濾波帶,其中非下采樣塔形分解結(jié)構(gòu)通過(guò)多級(jí)迭代的方式實(shí)現(xiàn),通過(guò)濾波器,圖像被分為1個(gè)二維低頻子帶和1個(gè)二維高頻子帶.要實(shí)現(xiàn)多級(jí)結(jié)構(gòu),對(duì)低頻子帶繼續(xù)迭代濾波即可.NSCT設(shè)計(jì)的非下采樣濾波帶放寬了約束條件,生成的濾波帶具有更好的頻域選擇性和規(guī)則性.相比CT可更好地進(jìn)行子帶方向分解.而相比小波分解,NSCT同樣是多尺度多分辨分解,但NSCT同時(shí)還是多方向分解,NSCT方向分解的靈活性可以使得圖像在分解過(guò)程中得到更多的細(xì)節(jié)信息.

圖1（a）顯示了 NSCT 分解流程圖,圖1（b）給出了分解后的系數(shù)矩陣分布示意圖,其中分解層數(shù)為3層,高頻層的方向數(shù)目分別為4和8的情形.

2 邊緣檢測(cè)算法

為了提取圖像中豐富的邊緣信息,本文利用NSCT能夠很好的逼近圖像中的方向這一優(yōu)良特性,以及邊緣特征主要體現(xiàn)在高頻成分中這一原理,設(shè)計(jì)一種基于NSCT的自適應(yīng)閾值圖像邊緣檢測(cè)算法.該算法首先對(duì)圖像進(jìn)行NSCT正變換,分別對(duì)高頻系數(shù)矩陣和低頻系數(shù)矩陣進(jìn)行處理:先去掉高頻系數(shù)矩陣中較小的系數(shù)值,并對(duì)處理的系數(shù)矩陣進(jìn)行NSCT反變換以得到高頻圖像,然后對(duì)高頻圖像處理得到高頻邊緣圖像；對(duì)低頻系數(shù)矩陣直接處理得到低頻邊緣圖像；而后利用低頻邊緣圖像對(duì)高頻邊緣圖像進(jìn)行補(bǔ)償,并針對(duì)補(bǔ)償后的邊緣圖像去除孤立噪聲點(diǎn)以得到最終的邊緣圖像.算法框圖如圖2所示.具體算法步驟為:

1）若輸入圖像為彩色圖像,則先將圖像轉(zhuǎn)換到Y(jié)UV空間,取其Y通道信息I；而若輸入圖像為灰度圖像,則直接令圖像其為I；

2）采用NSCT對(duì)輸入圖像I進(jìn)行多尺度多方向分解,即

其中,N+表示非下采樣輪廓波正變換,Qij代表圖像I的非下采樣輪廓波系數(shù)矩陣的集合,i表示尺度（i＝1,2,…,k）,j表示方向（j ＝1,2,…l,）；

3）利用自適應(yīng)閾值得到高頻圖像A.

（a）首先將其中低頻系數(shù)置為0,而后根據(jù)自適應(yīng)閾值去掉系數(shù)矩陣集合中的第2尺度層到第k尺度層系數(shù)中較小的系數(shù)值,僅保留較大的系數(shù)值,即:對(duì)于任意方向的第2尺度層到第k尺度層系數(shù)取值為

其中,ε為系數(shù)矩陣 Qij中以（r,c）點(diǎn)為中心的11*11的局部塊ψij的平均值,即為

這樣保留的均是較大的系數(shù),可以準(zhǔn)確地提取圖像的高頻信息以供提取豐富的邊緣信息.令H（Q）表示經(jīng)過(guò)上述處理后得到的系數(shù)矩陣集合.

（b）再對(duì)處理后的系數(shù)矩陣集合H（Q）進(jìn)行非下采樣輪廓波反變換得到高頻圖像,即

其中,N-表示非下采樣輪廓波反變換；

4）高頻圖像反映圖像的輪廓信息,在邊緣處取值較大,故對(duì)高頻圖像A采用閾值方法進(jìn)行處理,得到初始的高頻邊緣圖像B1.

其中,δ為自適應(yīng)閾值,具體計(jì)算方法為

其中,T為在高頻圖像A的平均值,λ為一乘子系數(shù),取值大于1.

5）為充分利用圖像低頻信息,補(bǔ)償高頻信息處理后出現(xiàn)的斷點(diǎn),NSCT內(nèi)的每個(gè)系數(shù)矩陣均與圖像大小一致,故利用canny算子直接對(duì)低頻系數(shù)矩陣Q11進(jìn)行處理得到低頻邊緣圖像B2,將低頻邊緣圖像B2作為對(duì)高頻邊緣圖像B1的修正和補(bǔ)償信息.而后對(duì)補(bǔ)償?shù)倪吘増D像進(jìn)行形態(tài)學(xué)處理以去除些微的孤立噪聲點(diǎn),得到最終的邊緣圖像B.

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

實(shí)驗(yàn)中,NSCT的分解層數(shù)為3層,高頻層上的方向數(shù)目分別為4和8,采用文獻(xiàn)［16］中的評(píng)價(jià)指標(biāo)評(píng)價(jià)檢測(cè)算法的性能,包括:CEN表示邊緣圖像中連續(xù)邊緣像素總數(shù),TEN表示邊緣圖像中總的像素?cái)?shù),R表示CEN和TEN的比率,定義為R＝CEN/TEN.R越大,說(shuō)明邊緣圖像的邊緣連續(xù)性越好.

首先將標(biāo)準(zhǔn)測(cè)試圖像lena作為邊緣檢測(cè)圖像,如圖3（a）所示,并將本文算法檢測(cè)結(jié)果與傳統(tǒng)的sobel算子、canny算子、log算子和prewitt算子結(jié)果做對(duì)比,分別如圖3（b） ～（f）所示.

從圖3可以看出,sobel算子和prewitt算子的檢測(cè)結(jié)果丟失了大量的邊緣信息,此處為在MATLAB庫(kù)中sobel算子和prewitt算子默認(rèn)閾值的檢測(cè)結(jié)果.利用canny算子的邊緣檢測(cè)效果較好,但邊緣過(guò)于細(xì)膩,一些非邊緣信息也被檢測(cè)出來(lái),導(dǎo)致圖像主要輪廓辨別不清楚,且圖像中的部分棱角邊緣被平滑,導(dǎo)致圖像中的棱角信息缺失.利用log算子的邊緣檢測(cè)效果相對(duì)來(lái)說(shuō)更好一些,但是存在不少的噪聲點(diǎn),并且丟失了部分邊緣信息,邊緣部分的連貫性不足.而本文所提出的方法與傳統(tǒng)的邊緣檢測(cè)方法相比魯棒性更好,邊緣定位更準(zhǔn)確,且能比傳統(tǒng)的邊緣檢測(cè)算子獲得更多的邊緣輪廓信息.同時(shí),表1對(duì)比給出了不同檢測(cè)算法對(duì)Lena圖像檢測(cè)結(jié)果的定量評(píng)價(jià)指標(biāo).從表1中可以看出,本文方法不僅檢測(cè)出的邊緣圖像中總的像素?cái)?shù)最多,說(shuō)明本文方法不僅能夠檢測(cè)出強(qiáng)邊緣,而且對(duì)于較弱的邊緣信息也能夠檢測(cè)出來(lái).同時(shí)本文方法的邊緣連續(xù)比率最高,說(shuō)明本文方法在檢測(cè)出的邊緣的連續(xù)性最好.

表1 Lena圖像的不同檢測(cè)方法性能比較Table 1 Performance comparisons of different edge detection approaches for image“Lena”

為測(cè)試算法的通用性,下面選擇flower和bird 2幅彩色圖像作為測(cè)試圖像,分別如圖4（a）和圖5（a）所示.將本文算法檢測(cè)結(jié)果與傳統(tǒng)的sobel算子、canny算子、log算子和prewitt算子結(jié)果做對(duì)比,分別如圖4（b） ～（f）及圖5（b） ～（f）所示.從圖4及圖5可以看出,sobel算子和prewitt算子的檢測(cè)結(jié)果丟失了較多的邊緣信息,邊緣不夠連續(xù),出現(xiàn)了較多的斷裂的小邊緣情形.利用canny算子的邊緣檢測(cè)效果較sobel算子和prewitt算子較好,但檢測(cè)出的邊緣圖像的主要輪廓辨別不清楚,不能描述圖像中的主要邊緣信息,且圖像中的部分棱角邊緣被平滑,導(dǎo)致圖像中的棱角信息缺失.利用log算子的邊緣檢測(cè)效果與canny算子的檢測(cè)效果類(lèi)似,但是存在不少的噪聲點(diǎn),并且丟失了部分邊緣信息,邊緣部分的連貫性不足.從圖4（e）及圖5（e）可以看出,本文所提出的方法與傳統(tǒng)的邊緣檢測(cè)方法相比魯棒性更好,邊緣定位更準(zhǔn)確,且能比傳統(tǒng)的邊緣檢測(cè)算子獲得更多的邊緣輪廓信息,同時(shí)較好的保持了邊緣部分的連貫性,使得邊緣圖像能夠更好的反應(yīng)圖像的輪廓信息.

表2給出了不同檢測(cè)算法對(duì)2副彩色圖像檢測(cè)結(jié)果的定量評(píng)價(jià)指標(biāo).從表2可以看出,本文方法不僅檢測(cè)出的邊緣圖像中總的像素?cái)?shù)最多,說(shuō)明本文方法不僅能夠檢測(cè)出強(qiáng)邊緣,而且對(duì)于較弱的邊緣信息也能夠檢測(cè)出來(lái).雖然針對(duì)flower圖像的邊緣檢測(cè)結(jié)果中,canny算子的比率R要略高于本文算法的比率R,但是本文算法能夠檢測(cè)出更多有效邊緣信息.在其余的邊緣檢測(cè)結(jié)果中,本文方法的邊緣連續(xù)比率最高,說(shuō)明本文方法在檢測(cè)出的邊緣的連續(xù)性最好.

表2 2幅彩色圖像的不同檢測(cè)方法性能比較Table 2 Performance comparisons of different edge detection approaches for two color images

本文利用NSCT分別提取圖像的高頻和低頻信息,對(duì)高頻信息利用自適應(yīng)閾值方法處理來(lái)獲得高頻邊緣圖像,同時(shí)也對(duì)低頻信息進(jìn)行處理得到低頻邊緣圖像,然后利用低頻邊緣對(duì)高頻邊緣進(jìn)行補(bǔ)償,得到最終的邊緣檢測(cè)圖像.在提取高頻信息時(shí)采用自適應(yīng)閾值方法,考慮了不同區(qū)域的能量不同這一重要特點(diǎn),便于更好的提取高頻邊緣特征.實(shí)驗(yàn)結(jié)果表明本文算法不僅能夠較好的提取灰度圖像的邊緣,同樣也適用于彩色圖像,較傳統(tǒng)的邊緣檢測(cè)算法能夠更好地定位邊緣,且提取到更為完整的圖像輪廓信息.

［1］鄭秀清,何坤,張健.基于Prewitt算子的TV圖像去噪算法［J］.四川師范大學(xué)學(xué)報(bào):自然科學(xué)版,2013,36（3）:456-462.

［2］ Roberts L G.Machine perception of three-dimension solids［C］//Optimal and Electro-Optimal Information Processing.MA:MIT Press,1965:99-197.

［3］ Canny J.A computational approach to edge detection［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence,1986,PAMI-8（6）:679-698.

［4］王娜,李霞.一種新的改進(jìn)Canny邊緣檢測(cè)算法［J］.深圳大學(xué)學(xué)報(bào):理工版,2005,22（2）:149-153.

［5］雒濤,鄭喜鳳,丁鐵夫.改進(jìn)的自適應(yīng)閾值Canny邊緣檢測(cè)［J］.光電工程,2009,36（11）:106-112.

［6］何春華,張雪飛,胡迎春.基于改進(jìn)Sobel算子的邊緣檢測(cè)算法的研究［J］.光學(xué)技術(shù),2012,38（3）:323-327.

［7］孫濤,王婧,王波.改進(jìn)的Sobel算子彩色圖像邊緣檢測(cè)［J］.計(jì)算機(jī)技術(shù)與應(yīng)用,2013,39（2）:128-133.

［8］張悅庭,孟曉峰,尹忠科,等.基于Contourlet模極大值的圖像邊緣檢測(cè)［J］.鐵道學(xué)報(bào),2008,30（5）:41-45.

［9］胡曉輝,張曉穎,陳俊蓮.一種融合小波變換和數(shù)學(xué)形態(tài)學(xué)的圖像邊緣檢測(cè)算法［J］.鐵道學(xué)報(bào),2011,33（3）:45-48.

［10］李哲濤,李仁發(fā),謝丼雄.基于全向小波的圖像邊緣檢測(cè)算法［J］.電子學(xué)報(bào),2012,40（12）:2451-2455.

［11］劉帥奇,胡紹海,肖揚(yáng),等.基于局部混合濾波的SAR圖像邊緣檢測(cè)［J］.電子與信息學(xué)報(bào),2013,35（5）:1120-1127.

［12］文武,苗放.基于小波的多尺度邊緣檢測(cè)在遙感圖像處理中的應(yīng)用［J］.微電子學(xué)與計(jì)算機(jī),2013,30（2）:18-20.

［13］張世華,宋振明.圖像模糊邊緣檢測(cè)算法的Delphi實(shí)現(xiàn)［J］.四川師范大學(xué)學(xué)報(bào):自然科學(xué)版,2003,36（5）:534-537.

［14］ Cunha A L,Zhou J,Do M N.The nonsubsampled contourlet transform:theory,design,and applications ［J］.IEEE Transactions on Image Processing,2006,15（10）:3089-3101.

［15］ Do M N,Vetterli M.The contourlet transform:an efficient directional multiresolution image representation ［J］.IEEE Transactions on Image Processing,2005,14（12）:2091-2106.

［16］張光玉,龔光珍,朱維樂(lè).基于克隆算法的彩色圖像邊緣檢測(cè)新算法［J］.電子學(xué)報(bào),2006,34（4）:702-707.