田國富，宋江波

　　(沈陽工業大學 機械工程學院，遼寧 沈陽 110870)

      摘要： 針對傳統的邊緣檢測算子對噪聲敏感問題，提出了一種基于多尺度多結構元素形態學邊緣檢測算子。該算子在傳統形態學邊緣檢測的基礎上，通過改變結構元素的形狀和尺度，在很好抑制噪聲的前提下，實現圖像精細邊緣提取。MATLAB仿真結果表明，該邊緣檢測算子具有更強的去噪能力，定位準確，保留了更多的圖像細節，優于傳統的邊緣檢測算子。

　　關鍵詞：數學形態學；多尺度多結構元素；邊緣檢測

0引言

　　圖像最基本的特征是邊緣，它包含著大量的圖像信息，因此對圖像邊緣進行完整、精準的檢測在圖像處理過程中有著非常重要的意義。圖像邊緣檢測就是定位和提取邊緣信息，為圖像分析、目標識別和機器視覺等做前期準備［1］。由于噪聲對圖像邊緣檢測帶來的不利影響，在邊緣檢測過程中會丟失大量的圖像細節信息，因此圖像處理的基本問題就是在盡可能降低噪聲干擾的同時提高檢測精度。

　　傳統的邊緣檢測算子起步比較早，如Roberts算子、Sobel算子、Canny算子等一階微分算子，以及Laplacian算子、LOG（高斯—拉普拉斯）算子等二階微分算子。這些算法大多是在以要處理像素為中心的鄰域中進行灰度分析，實現對圖像邊緣的提?。?］，計算簡單，實時性好，有良好的邊緣檢測性。但由于微分運算的性質，對噪聲比較敏感，抗噪聲性能差，提取的邊緣不夠精細，因此在實際應用中受到限制［3］。近年來，很多學者提出了基于數學形態學的邊緣檢測算法，在此基礎上，本文針對單一尺度的結構元素容易導致在對目標圖像進行邊緣提取時定位不夠準確、邊緣信息丟失、去噪效果達不到預期要求等問題，提出了一種基于多尺度多結構元素的形態學邊緣檢測算法。

1數學形態學基礎

　　數學形態學是建立在數學理論之上的新興學科，理論基礎是集合論，基本運算是像素間的邏輯關系而非代數關系。其基本思想是：用一個具有一定形狀的結構元素去度量和提取圖像中的對應形狀以達到對圖像進行分析和識別的目的［4］。

　　數學形態學可以去除目標圖像中不相關的形態結構屬性，而保留本質的結構屬性，達到簡化圖像數據的目的。目前，數學形態學已廣泛應用于機器視覺、顯微圖像分析、醫學處理等領域 ［5］。

　　1.1數學形態學的基本運算

　　數學形態學基本運算主要有膨脹、腐蝕、開、閉運算，開運算和閉運算都是由膨脹和腐蝕運算組合而成。

　　設b(m,n)是結構元素，f(x,y)為輸入的灰度圖像，b的定義域為Db，f的定義域為Df，則結構元素b對灰度圖像f的形態學運算如下：

　?。?）形態學膨脹運算

　　圖像f被b膨脹，記作f⊕b，定義為：

　　f⊕b=max{f(x－m,y－n)+b(m,n)

　　|(x－m,y－n)∈Df;(m,n)∈Db} (1)

　?。?）形態學腐蝕運算

　　圖像f被b腐蝕，記作fΘb，定義為：

　　fΘb=min{f(x+m,y+n)-b(m,n)

　　|(x+m,y+n)∈Df;(m,n)∈Db}(2)

　?。?）形態學開運算

　　圖像f被b腐蝕后再用b來膨脹腐蝕結果，記作fb，定義為：

　　fb=f(x,y)Θb(m,n)⊕b(m,n)(3)

　　(4)形態學閉運算

　　圖像f被b膨脹后再用b來腐蝕膨脹結果，記作f·b，定義為：

　　f·b=f(x,y)⊕b(m,n)Θb(m,n)(4)

　　1.2形態學濾波

　　數字圖像的噪聲主要來源于圖像的獲取和傳輸過程，對圖像邊緣檢測產生不利的影響。因此在對圖像進行邊緣檢測之前必須對其進行平滑處理，去除噪聲的影響［6］。常用的圖像去噪方法有中值濾波、小波去噪、數學形態學濾波等。

　　數學形態學中，由于開運算可以去除比結構元素更小的明亮細節，閉運算可以去除比結構元素更小的暗細節，所以它們經常組合在一起用來平滑圖像并去除噪聲［7］。但由于閉運算的擴展性和開運算的反擴展性使得兩者單獨使用時并不能獲得良好的濾波效果，因此本文采用兩者的均值來實現噪聲圖像的平滑處理。假設經過去噪之后的圖像用Filter來表示，則公式表示如下：

　　Filter1={［(f⊕b)Θb］Θb}⊕b=(f·b) b(5)

　　Filter2={［(fΘb)⊕b］⊕b}Θb=(f b)·b(6)

　　Filter=12(Filter1+Filter2)

　　=1/2［(f·b) b］+［(f b)·b］(7)

　　2數學形態學邊緣檢測算法及其改進

　　2.1數學形態學邊緣檢測算法

　　對上文介紹的形態學的基本運算進行組合，便可得到基本的形態學邊緣檢測算子的數學表達式，設E為邊緣圖像。

　　若采用膨脹運算，則邊緣檢測算子為：

　　E1(x,y)=f(x,y)⊕b(m,n)－f(x,y)(8)

　　若采用腐蝕運算，則邊緣檢測算子為：

　　E2(x,y)=f(x,y)－f(x,y)Θb(m,n)(9)

　　若采用膨脹腐蝕復合運算，則邊緣檢測算子為：

　　E3(x,y)=f(x,y)⊕b(m,n)－f(x,y)Θb(m,n)(10)

　　由于形態邊緣檢測算子都屬于非線性差分算子，因此對噪聲都比較敏感，而且檢測出的邊緣與結構元素b緊密相關。利用數學形態學腐蝕和開運算組成的邊緣檢測器可以抑制圖像中的峰值噪聲，而利用數學形態學膨脹和閉運算可以抑制圖像中的低谷噪聲。因此對上述形態學邊緣檢測算子改進后得到抗噪型邊緣檢測算子。

　　(1)抗噪膨脹型：

　　E4(x,y)=［f(x,y)b(m,n)］⊕b(m,n)－

　　［f(x,y)b(m,n)］·b(m,n)(11)

　　(2)抗噪腐蝕型：

　　E5(x,y)=［f(x,y)·b(m,n)］b(m,n)－

　?。踗(x,y)·b(m,n)］Θb(m,n)(12)

　　(3)抗噪膨脹腐蝕型：

　　E6(x,y)=［f(x,y)b(m,n)］⊕b(m,n)－

　?。踗(x,y)·b(m,n)］Θb(m,n)(13)

　　2.2多結構元素邊緣檢測算法

　　上述式(11)～(13)介紹的改進的數學形態學邊緣檢測算子可以很好地去除噪聲和保持邊緣細節，很大程度上減少了對噪聲的敏感程度，但各算子均使用單一結構元素，因此檢測到的邊緣比較弱、不連續，而且很難檢測出復雜形狀的邊緣。文獻［8］針對以上問題提出了基于多結構元素的邊緣檢測算法，公式如下：

　　E7(x,y)=［f(x,y)Θb1(m,n)⊕b2(m,n)］⊕

　　b3(m,n)－［f(x,y)⊕b1(m,n)Θb2(m,n)］Θb3(m,n)(14)

　　結構元素b1、b2、b3均采用單一尺度3×3的方形結構元素,分別對原圖像和噪聲圖像采用文獻［8］算法進行邊緣檢測。仿真試驗結果如圖1所示。

　　上圖中，(a)為原圖像，(b)、(c)、(d)分別為采用文獻［8］的檢測算法對原圖像、密度0.02的噪聲圖像、密度為0.05的噪聲圖像提取邊緣的結果。對于文獻［8］提出的多結構元素的邊緣檢測算法，b1、b2、b3均采用單一尺度的結構元素。由仿真結果可以看出，圖像沒有添加椒鹽噪聲時，檢測到的圖像邊緣清晰，保留了圖像的邊緣細節，當添加噪聲之后，雖檢測出的圖像邊緣較平滑，但圖像中混有噪聲點，并且通過仿真試驗可知隨著噪聲密度的增大，邊緣圖像中的噪聲點將越來越多，當噪聲點多到一定程度之后，就會與圖像邊緣粘黏在一起，形成偽邊緣，導致不能準確地檢測圖像的輪廓邊緣。

　　2.3多尺度多結構元素邊緣檢測算法

　　本文結合式(13)提出另一種抗噪膨脹腐蝕型邊緣檢測算子：圖2各算法邊緣檢測結果對比

　　E8(x,y)=［f(x,y)·b(m,n)］⊕b(m,n)－

　?。踗(x,y)b(m,n)］Θb(m,n)(15)

　　對其改進得到多結構元素的邊緣檢測算子：

　　E9(x,y)=［f(x,y)⊕b1(m,n)Θ

　　b2(m,m)］⊕b3(m,n)－［f(x,y)Θ

　　b1(m,n)⊕b2(m,n)］Θb3(m,m)(16)

　　由于用較小尺寸的結構元素對圖像進行處理，能較好地保護圖像的細節，也可以濾除一定尺寸的噪聲；反之，用較大尺寸的結構元素對圖像進行處理時，能濾除較大尺寸的噪聲，但也會損失更多的圖像細節。針對不同尺度結構元素的優缺點，本文結合式(14)、(16)對文獻［8］的多結構元素邊緣檢測算子進行改進，提出先用較大結構元素b1和b2對圖像進行開、閉運算處理，然后用小尺度結構元素b3對圖像做膨脹和腐蝕操作。取結構元素b1、b2為不同尺度平坦的圓形結構元素，半徑分別為2和3，b3為R=1的菱形結構元素，得到一種多尺度多結構元素形態學邊緣檢測算子為：

　　E(x,y)=1/2{［f(x,y)Θb1(m,n)⊕b2(m,n)］

　　⊕b3(m,n)－［f(x,y)⊕b1(m,n)Θb2(m,n)］

　　Θb3(m,n)+［f(x,y)⊕b1(m,n)Θb2(m,n)］

　　⊕b3(m,n)－［f(x,y)Θb1(m,n)⊕b2(m,n)］

　　Θb3(m,n)}(17)

3試驗分析

　　為驗證本文提出的基于多尺度多結構元素算法的有效性，將其與Canny邊緣檢測算法、Sobel邊緣檢測算法和一般抗噪型邊緣檢測算法(式（13）)進行對比，用MATLAB對Lenna圖像和添加噪聲的Lenna圖像、Rice圖像進行仿真。

　　本文圖像邊緣檢測算法步驟如下：

　　(1)輸入Lenna、Rice灰度圖像，添加噪聲(密度為0.05的椒鹽噪聲)。

　　(2)用本文式(7)的方法對加入噪聲的圖像進行開—閉濾波運算，開、閉運算均采用3×3的方形結構元素。

　　(3)利用公式(17)對經過濾波處理后的圖像進行邊緣提取。其中，b1和b2分別是半徑為2和3的平坦的圓形結構元素，b3為R=1的菱形結構元素［9］。

　　實驗結果如圖2所示。

　　由以上檢測結果可以看出，對于Lenna圖像，無噪聲情況下，傳統的Sobel算子出現了很多的邊緣斷裂現象，Canny算子檢測出過多的邊緣細節，而一般的抗噪型邊緣檢測算子(式（13）)和本文的算法檢測的輪廓邊緣清晰，連貫性較好。對于添加噪聲之后的Lenna圖像，傳統的Sobel與Canny算子檢測結果出現很大的噪聲，不能清楚地觀察到圖像邊緣。一般的抗噪型邊緣檢測算子雖然抵抗噪聲能力比較強，但邊緣出現了很多缺口，平滑度達不到要求。對于Rice圖像，由于其邊緣輪廓比較簡單，無噪聲情況下，幾種檢測算子都顯示出不錯的邊緣檢測結果。但添加噪聲之后，Sobel算子已經完全不能識別圖像的邊緣，Canny算子不能準確地定位圖像輪廓邊緣，產生很多偽邊緣。一般的抗噪型算子邊緣不連貫，出現很多的缺口，平滑度也不夠好。本文提出的算法，不管是對原圖像還是添加噪聲的圖像，都顯示出很好的邊緣檢測性能，邊緣定位準確，輪廓清晰完整，邊緣細節更加豐富。分別對Lenna圖像和Rice圖像的檢測，也充分說明了本文檢測算法的普遍性和實用性。

4結束語

　　本文針對傳統的邊緣檢測算子抗噪能力差的特點，從數學形態學的基本思想出發，通過對單一尺寸或形狀結構元素的邊緣檢測算法分析與改進，提出了多尺度多結構元素邊緣提取算法，通過仿真結果驗證，本文算法可在抑制噪聲的同時，有效地提取圖像的邊緣，且提取的邊緣定位精度高，平滑度好，邊緣細節保持也比較完整，充分驗證了本文算法的優越性。

參考文獻

　?。?］ 范立南,韓曉微,王忠石,等. 基于多結構元的噪聲污染灰度圖像邊緣檢測研究［J］. 武漢大學學報(工學版),2003,36(3):8690.