王東霞, 許偉昶

(濟源職業(yè)技術(shù)學(xué)院,河南 濟源 459000)

0 引 言

圖像邊緣一般是指圖像中像素灰度有階躍變化或屋頂狀變化的那些像素的集合，可以粗略地分為階躍邊緣和屋頂邊緣，它存在于目標與背景、目標與目標、區(qū)域與區(qū)域、基元與基元之間[1]。邊緣檢測在圖像預(yù)處理中有非常大的作用，其在圖像處理及人臉識別等方面的應(yīng)用也是非常廣泛。目標對象的輪廓可以被它清晰的描繪出來，所含概的信息量非常大，如階躍性質(zhì)、形狀等，使讀者清晰明了。在邊緣檢測算法研究中，提出了比較多的檢測算法，如: Sobel算子﹑Canny算子[2]和梯度算子等，但這些算法通常較難得到良好的檢測結(jié)果，往往存在較差的抗干擾能力以及對方向性過于敏感等不足之處?；谛〔ㄗ儞Q邊緣檢測方法[3-4]具有良好的局域化特性，它是檢測邊緣很好的工具，但是該方法邊緣信息難以保留完全。數(shù)學(xué)形態(tài)學(xué)[5]是一種非線性的濾波方法，應(yīng)用數(shù)學(xué)形態(tài)學(xué)可以簡化圖像數(shù)據(jù)，保持他們基本的形狀特征。形態(tài)學(xué)在檢測圖像邊緣的連續(xù)性以及各個方向性方面都優(yōu)于其他的傳統(tǒng)方法[6-8]，由于單一的結(jié)構(gòu)元素只能檢測出與結(jié)構(gòu)元素同方向的邊緣信息，而對其它方向的邊緣并不敏感。不過，由于傳統(tǒng)的數(shù)學(xué)形態(tài)學(xué)檢測算法忽略了邊緣的方向性，只針對圖像邊緣的強度信息進行處理。這種處理方式具有一定的局限性，存在圖像檢測出的邊緣相對較寬以及分辨率低的不良現(xiàn)象。而且，使用單一的形態(tài)學(xué)梯度算法來對邊緣信息進行提取，雖然可以抑制噪聲，但也會存在對部分形態(tài)學(xué)梯度幅值與噪聲相近或低于噪聲的低強度邊緣造成損壞。

針對上述傳統(tǒng)圖像邊緣檢測算法的不足，文中提出了一種基于改進的數(shù)學(xué)形態(tài)邊緣檢測算法。實驗仿真表明，該算法在檢測圖像邊緣信息的完整性方面具有良好的效果，同時，能確保圖像邊緣的平滑性。

1 傳統(tǒng)形態(tài)邊緣檢測算子

數(shù)學(xué)形態(tài)學(xué)是一門建立在嚴格數(shù)學(xué)集合論基礎(chǔ)上的學(xué)科，它的基本思想是用具有一定形態(tài)的結(jié)構(gòu)元素去度量和提取圖像中的對應(yīng)形狀以達到對圖像分析和識別的目的[11]。數(shù)學(xué)形態(tài)學(xué)在圖像數(shù)據(jù)的簡化方面具有優(yōu)勢，同時能夠確保圖像數(shù)據(jù)的形狀特征不改變，并將不相關(guān)的結(jié)構(gòu)進行清理。數(shù)學(xué)形態(tài)學(xué)可以進行相關(guān)的運算，其主要的運算規(guī)則有：膨脹、腐蝕、開啟以及閉合。

腐蝕運算的主要作用是使物體沿著邊界逐漸減小的過程，膨脹則是合并的過程，它們之間的關(guān)系是互為對偶運算。

數(shù)學(xué)運算定義[9]如下：

(1) 腐蝕運算：

FΘB(x,y)=min{F(x-s,y-t)-B(s,t)}

(1)

(2) 膨脹運算：

F⊕B(x,y)=max{F(x-s,y-t)+B(s,t)}

(2)

式中：F(x,y)表示的是灰度圖像；B(s,t)表示的是結(jié)構(gòu)元素。

閉運算主要用來對物體內(nèi)部細小空洞進行填充，開運算的功能則是用來將細小的物體刪除，并將邊緣某些不必要的內(nèi)容清除。相應(yīng)的，閉運算和開運算也是互為對偶的運算。

數(shù)學(xué)公式定義如下：

(3) 開運算：

F°B=(FΘB)⊕B

(3)

(4) 閉運算：

F·B=(F⊕B)ΘB

(4)

式中：Θ用來描述的是腐蝕運算；⊕用來描述的是膨脹運算。形狀決定圖像中的對象及其特征。針對時域空間形狀的研究是數(shù)學(xué)形態(tài)學(xué)的主要方向，因此，在圖像處理方面應(yīng)用形態(tài)學(xué)是非常適合的。膨脹運算、腐蝕運算、開啟運算以及閉合運算都是基于集合的運算。其中，結(jié)構(gòu)元素用來調(diào)節(jié)圖像特征變換的幾何結(jié)構(gòu)，扮演著至關(guān)重要的角色。

在膨脹運算、腐蝕運算、開啟運算以及閉合運算的基礎(chǔ)上，能夠推演出灰度形態(tài)梯度邊緣檢測算子、灰度腐蝕邊緣檢測算子和灰度膨脹邊緣檢測算子三個常見算子。它們的數(shù)學(xué)定義如下：

灰度膨脹邊緣檢測算子

Cgd=(f⊕g)(x,y)-f(x,y)

(5)

灰度腐蝕邊緣檢測算子

Cge=f(x,y)-(fΘg)(x,y)

(6)

灰度形態(tài)梯度邊緣檢測算子

Cgsun=(f⊕g)(x,y)-(fΘg)(x,y)

(7)

通過分析，可以看出上面的三個算子都是屬于非線性差分算子，其本質(zhì)在某種程度上是對傳統(tǒng)線性差分算子的推廣應(yīng)用。式(5)用來描述形態(tài)膨脹，這種算子處理得到的圖像的邊緣比較模糊、產(chǎn)生的邊緣信號也比較弱。式(6)用來描述形態(tài)腐蝕，這種算子處理得到的圖像邊緣信號比較強，不過在噪聲方面有所增強。文獻(10)中修正了灰度形態(tài)學(xué)梯度邊緣檢測算子：

(8)

文獻中對邊緣檢測算子進行數(shù)學(xué)定義：

E(x,y)=Ed(x,y)+0.5Edec(x,y)

(9)

通過上面的公式能夠清晰地反映出文獻(10)在檢測結(jié)果的模糊性方面具有一定的優(yōu)勢，能夠降低模糊性，是因為該算子對原圖像的一些邊緣細節(jié)信息進行了迭加。因為光源照射等原因，圖像中具有更低的噪聲梯度幅值以及形態(tài)學(xué)梯度幅值。因此，在圖像處理中提取邊緣時，如果只應(yīng)用式(9)梯度閾值方法進行則往往容易對這類邊緣造成損壞。

2 形態(tài)邊緣檢測算子的改進

2.1 算子原理

針對傳統(tǒng)的形態(tài)邊緣檢測算子所存在的種種不足，文中提出了一種抗噪型檢測算子，是對傳統(tǒng)的檢測算子的改進。

通過對形態(tài)學(xué)邊緣檢測算子的研究可以看出，如果出現(xiàn)噪聲點的尺度小于檢測算子的結(jié)構(gòu)元尺度的情況，則在去除圖像噪聲方面形態(tài)邊緣檢測算子可以顯現(xiàn)出明顯的效果。雖然這種大尺度的結(jié)構(gòu)元素具有很強的噪聲去除能力，便于對原圖像中對象的整體輪廓進行定位，不過，也存在一定的不足，往往會將原圖像中某些細小的邊緣信息當作噪聲去處理。而如果出現(xiàn)檢測算子的結(jié)構(gòu)元尺度小于噪聲點的尺度的情況，則存在對圖像中的噪聲難以消除的現(xiàn)象，雖然其在圖像細節(jié)保持方面具有很好的效果。

針對上述存在的不足問題，文獻中給出了多尺度多結(jié)構(gòu)元邊緣檢測算子來解決該問題。其數(shù)學(xué)定義：

(10)

式中：EN(x,y)用來描述多尺度圖像邊緣檢測結(jié)果；EM(x,y)用來描述多結(jié)構(gòu)元圖像邊緣檢測結(jié)果；ai用來描述加權(quán)系數(shù)。不過，因為計算復(fù)雜度高等因素的影響，該算子一直難以在實際的研究中推廣應(yīng)用。

鑒于此，本文在該算子的基礎(chǔ)上進行了改進，即使用雙結(jié)構(gòu)不同尺度的兩個結(jié)構(gòu)元素來對圖像進行交替順序形態(tài)開閉濾波，對噪聲采用平滑圖像的方式去除處理，對單一結(jié)構(gòu)元素處理所存在的問題進行了根本性的解決。

改進的雙結(jié)構(gòu)元素多尺度形態(tài)學(xué)邊緣檢測算子的數(shù)學(xué)定義：

E(F)=(F°B)⊕A-(F·B)ΘA

(11)

式中，A的用鉆石型5×5結(jié)構(gòu)元素來描述，而B用十字型3×3結(jié)構(gòu)元素來描述，它們都是基于不同尺度的結(jié)構(gòu)素。

式中：A為小尺度的結(jié)構(gòu)元素，它的優(yōu)點是能夠?qū)D像的邊緣信息很好的保持住，但它也存在不足之處，即對噪聲去除處理欠缺；而B為大尺度的結(jié)構(gòu)元素，它的優(yōu)點是對噪聲去除處理效果好，但也存在不足之處就是會將一些邊緣信息去除。因此，在文中提出的改進算法中，充分的將A、B的優(yōu)點保存下來，而將其不足去除，這樣不但能有效地保持好邊緣信息，而且對噪聲的處理也非常有效果。

2.2 改進的邊緣檢測算法的描述

上面對改進的形態(tài)邊緣檢測算子的原理進行了說明。針對該改進的算法的執(zhí)行流程用圖1描述。

圖1 算法流程圖

本文改進的算法的執(zhí)行步驟如下：

Step1：先將圖像輸入，然后利用形態(tài)開閉平滑對圖像進行處理，目的是將圖像中的噪聲進行去除。

Step2：計算多尺度邊緣結(jié)果以及雙結(jié)構(gòu)元素邊緣。

Step3：在提取圖像的邊緣信息時采用不同方向的結(jié)構(gòu)元素。然后將各個方向的邊緣進行加權(quán)求和處理，最終就可以計算出圖像的邊緣信息。利用公式(12) 求得加權(quán)和結(jié)果：

(12)

式中，ai用來描述加權(quán)系數(shù)。

Step4：二值化處理得到的邊緣信息，最后得到檢測的邊緣。

3 實驗結(jié)果與分析

在Matlab7.0實驗環(huán)境下，根據(jù)上面所提出的方法，進行圖像邊緣檢測試驗，本文算法簡單、易于實現(xiàn)。給出了兩個實驗，① 將算法應(yīng)用在灰度圖像中，② 本文算法應(yīng)用在彩色圖像邊緣檢測中，其中兩幅圖像均為噪聲圖像。為了檢測算法的有效性，本文給出了幾種常見的邊緣檢測算法結(jié)果圖，見圖2及圖3。在本文仿真實驗中，選用兩種結(jié)構(gòu)元素類型為十字型結(jié)構(gòu)3×3和鉆石型5×5結(jié)構(gòu)型。

從實驗結(jié)果不難看出，在初始邊緣檢測結(jié)果(邊緣圖像)中，原始圖像模糊的部分邊緣像素的連續(xù)性較差。在更大尺度上檢測到的邊緣是對初始檢測結(jié)果的有效補充。這種補充是建立在利用圖像信息基礎(chǔ)上的，并且因為尺度較大，在補充邊緣信息的同時，有效地抑制了噪聲的干擾。

(a) 噪聲圖像

(b) 小波方法

(c) 傳統(tǒng)形態(tài)學(xué)檢測

(d) 文獻10檢測

(e) 本文算法

(a) 噪聲圖像

(b) 小波方法

(c) 傳統(tǒng)形態(tài)學(xué)檢測

(d) 文獻10檢測

(e) 本文算法

4 結(jié) 語

針對傳統(tǒng)邊緣檢測算法存在的邊緣分辨率較低、低強度邊緣保護能力較差等問題，提出了基于數(shù)學(xué)形態(tài)學(xué)彩色圖像邊緣檢測改進算法。該方法將利用數(shù)學(xué)形態(tài)四種運算，膨脹、腐蝕、開、閉等變換以及它們的組合，并根據(jù)不同的結(jié)構(gòu)元素的尺度大小和結(jié)構(gòu)元類型，給出了一種改進的形態(tài)學(xué)抗噪型邊緣檢測算子，有效地檢測出完整的圖像邊緣信息，并保持圖像邊緣的平滑性。實驗結(jié)果證明了，與傳統(tǒng)的邊緣檢測算法相比，本文提出的算法能有效提取準確的邊緣信息，而且又具有很強的抗噪性，是一種有效的邊緣檢測方法。更加充分地利用了圖像信息，提高了邊緣檢測的正確率。

[1] 孫江雄,谷達寶. 基于小波域多狀態(tài)隱馬爾科夫樹模型多尺度文本圖像分割[J]. 模式識別, 2004,37(7):1315-1324.

SUN J X,GU D B. A multiscale edge detection algorithm based on wavelet domain vector hidden Markov tree model[J].PatternRecognition,2004,37(7):1315-1324.

[2] 張 琳,薄 鵬. 小波域內(nèi)基于尺度乘法的邊緣檢測[J].模式識別, 2002,23(14):1771-1784.

ZHANG L, BAO P. Edge detection by scale multiplication in wavelet domain[J]. Pattern Recognition Letters, 2002,23(14):1771-1784.

[3] 付永慶,王詠勝.一種基于數(shù)學(xué)形態(tài)學(xué)的灰度圖像邊緣檢測算法[J].哈爾濱工程大學(xué)學(xué)報,2005,26(5):685-687.

FU Yong-qing,WANG Yong-sheng.An algorithm for edge detection of gray-scale image based on mathematical morphology[J]. Journal of Harbin Engineering University,2005,26(5):685-687.

[4] 楊 暉,張繼武.數(shù)學(xué)形態(tài)學(xué)在圖像邊緣檢測中的應(yīng)用研究[J].遼寧大學(xué)學(xué)報,2005,32(1):50-53.

YANG Hui,ZHANG Ji-wu. Research on Application of Mathematical Morphology in Edge Detection of Image[J]. Journal of Liaoning University,2005,32(1):50-53.

[5] 劉 清，林士勝.基于數(shù)學(xué)形態(tài)學(xué)的圖像邊緣近側(cè)算法[J].華南理工大學(xué)學(xué)報，2008，36(9):113-116.

LIU Qing,LIN Shi-sheng. Image Edge Processing based on Mathematical-morphology[J]. Journal of South China University of Technology，2008，36(9):113-116.

[6] 林世毅，蘇廣川，陳 東，等.基于小波變換和數(shù)學(xué)形態(tài)學(xué)的邊緣檢測法[J].儀器儀表學(xué)報，2004，25(4)：685-687.

LIN Shi-yi,SU Guang-chuan,CHEN Dong,etal. Edge Detection Based on Wavelet Transform and Morphology[J]. Journal of Scientific Instrument，2004，25(4)：685-687.

[7] 崔 屹.圖像處理與分析——數(shù)學(xué)形態(tài)學(xué)方法及應(yīng)用[M].背景：科學(xué)技術(shù)出版社，2000:56-201.

[8] Canny J E. A computional approach to edge detection[J].IEEE Trans on PAM 1,1985,8(6):679-698.

[9] 羅圣敏. 基于數(shù)學(xué)形態(tài)學(xué)的自適應(yīng)邊緣檢測新算法[J]. 計算機仿真，2009，29(4):997-1000.

LUO Sheng-min. Adaptive algorithm of edge detection based on mathematical morphology[J]. Computer Simulation，2009,29(4):997-1000.

[10] 王萬良,周 明,徐新黎.基于改進粒子群算法的離子膜車間調(diào)度問題研究[J].控制與決策,2010,25(7):1021-1025.

WANG Wan-liang,ZHOU Ming, XU Xin-li. Based on improved particle swarm caustic Shop Scheduling Problem[J]. Control and Decision,2010,25(7):1021-1025.

[11] 何怡剛,祝文姬,周炎濤.基于粒子群算法的模擬電路故障診斷方法[J].電工技術(shù)學(xué)報,2010,25(6):163-171.

HE Yi-gang,ZHU Wen-ji,ZHOU Yan-tao. Analog circuit fault diagnosis method based on particle swarm optimization algorithm[J]. Journal of Electrician Technique,2010,25(6):163-171.

[12] 鄒 恩,陳建國,李祥飛.一種改進的變尺度混沌優(yōu)化方法及其仿真研究[J].系統(tǒng)仿真學(xué)報,2006,18(9):2426-2430.

ZOU En, CHEN Jian-guo, LI Xiang-fei. Research on an improved mutative scale chaos optimization method and its simulation[J]. Journal of system simulation,2006,18(9):2426-2430.

[13] 滕 皓,邵闊義,曹愛增,等. 量子遺傳算法的變尺度混沌優(yōu)化策略研究[J]. 計算機應(yīng)用研究,2009,26(2):543-546.

TENG Hao, SHAO Kuo-yi, CAO Ai-zeng,etal.Study on mutative scale chaos optimization strategy of quantum genetic algorithm[J]. The research and application of computer,2009,26(2):543-546.

[14] 滕 皓,曹愛增,楊炳儒.一種改進變尺度混沌優(yōu)化的模糊量子遺傳算法[J].計算機工程,2010,36(7):175-177.

TENG Hao, CAO Ai-zeng, YANG Bing-ru. An improved mutative scale chaos optimization fuzzy quantum genetic algorithm[J]. Computer Engineering,2010,36(7):175-177.

[15] 閆海霞，趙曉暉. 基于書序形態(tài)學(xué)的邊緣檢測方法[J]. 計算機應(yīng)用研究，2008,25(11):3497-3499.

YAN Hai-xia, ZHAO Xiao-hui. Detection method of book order morphology based edge[J]. The research and application of computer，2008,25(11):3497-3499.