1、課程設計資料袋理學院學院（系、部）2013-2014學年第2學期課程名稱數字圖像處理指導教師職稱講師學生姓名專業班級信1f1102學號題目經典邊緣檢測算子比較成績起止日期2014年6月16日2014年6月20日序號材料名稱資料數量備注1課程設計任務書12課程設計說明書13源程序（電子文檔）145課程設計任務書2013-2014學年第2學期理學院學院（系、部）信息與計算科學專業1102班級課程名稱：數字圖像處理設計題目：經典邊緣檢測算子比較完成期限：自2014年6月16日至2014年6月20日共1周內容及任務將幾種經典邊緣算子進行比較，編寫出相應的程序。進度安排起止日期工作內容2014.6.16

2、2014.6.16選擇題目，查閱資料。2014.6.172014.6.19利用matlab軟件編寫相應的程序，并嘗試改進已有的算法。2014.6.202014.6.20整理相關的結果，書寫課程設計報告。主要參考資料1RafaelC.Gonzalez、RichardE.Woods.岡薩雷斯.數字圖像處理（第二版）電子工業出版社2陳傳波、金先級.數字圖彳t處理.北京：機械工業出版社，2004.73KennethR.Castleman者；朱志剛、石te機等譯.數子圖彳裒處理.北東：電子工業出版社，2002.24陳天華.數字圖彳t處理.北京：清華大學出版社，2007.65曹茂永.數字圖彳t處理.北京：

3、北京大學出版社，2007.96龔聲蓉、劉純平、王強等編著.數字圖像處理與分析.北京：清華大學出版社指導教師（簽字）：年月日系（教研室）主任（簽字）：數值分析課程設計說明書經典邊緣檢測算子比較起止日期：2014年6月16日至2014年6月20日學生姓名班級信息與計算科學1102班學號成績指導教師（簽字）理學院2014年6月20日摘要5關鍵字5引言1 .各種經典邊緣檢測算子原理簡介1.1 Roberts（羅伯特）邊緣檢測算子.1.2 Sobel（索貝爾）邊緣檢測算子1.3 Prewitt（普瑞維特）邊緣檢測算子.1.4 Laplacian（拉普拉斯）邊緣檢測算子1.5 Marr-Hildreth（

4、馬爾）邊緣檢測算子.1.6 canny（凱尼）邊緣檢測算子2 .各種經典邊緣檢測算子性能比較2.1 MATLAB程序仿真2.2 實驗效果比較3 .結論參考文獻：錯誤！未定義書簽。5677891012121516.1617附錄.經典邊緣檢測算子比較摘要：圖像邊緣檢測技術是圖像分割、目標識別、區域形態提取等圖像分析領域中十分重要的基礎。本文簡要介紹各種經典圖像邊緣檢測算子的基本原理，用Matlab仿真實驗結果表明各種算子的特點及對噪聲的敏感度，為學習和尋找更好的邊緣檢測方法提供參考價值。關鍵字：圖像處理；邊緣檢測；算子；比較引言圖像的邊緣時圖像最基本的特征之一。所謂邊緣（或邊沿）是指周圍像素灰度有

5、階躍性變化或“屋頂”變化的那些像素的集合。邊緣廣泛存在于物體與背景之間、物體與物體之間、基元與基元之間，因此它是圖像分割依賴的重要特征。圖像邊緣對圖像識別和計算機分析十分有用，邊緣能勾劃出目標物體，使觀察者一目了然；邊緣蘊含了豐富的內在信息（如方向、階躍性質、形狀等）。從本質上說，圖像邊緣是圖像局部特性不連續性（灰度突變、顏色突變、紋理結構突變等）的反應，它標志著一個區域的終結和另一個區域的開始。邊緣檢測技術是所有基于邊界分割的圖像分析方法的第一步，首先檢測出圖像局部特性的不連續性，再將它們連成邊界，這些邊界把圖像分成不同的區域，檢測出邊緣的圖像就可以進行特征提取和形狀分析。為了得到較好的邊緣

6、效果，現在已經有了很多的邊緣檢測算法以及一些邊緣檢測算子的改進算法。但各算子有自己的優缺點和適用領域。本文著重對一些經典邊緣檢測算子進行理論分析、實際驗證并對各自性能特點做出比較和評價，以便實際應用中更好地發揮其長處，為新方法的研究提供衡量尺度和改進依據。1 .各種經典邊緣檢測算子原理簡介圖像的邊緣對人的視覺具有重要的意義，一般而言，當人們看一個有邊緣的物體時，首先感覺到的便是邊緣?；叶然蚪Y構等信息的突變處稱為邊緣。邊緣是一個區域的結束，也是另一個區域的開始，利用該特征可以分割圖像。需要指出的是，檢測出的邊緣并不等同于實際目標的真實邊緣。由于圖像數據時二維的，而實際物體是三維的，從三維到二維的

7、投影必然會造成信息的丟失，再加上成像過程中的光照不均和噪聲等因素的影響，使得有邊緣的地方不一定能被檢測出來，而檢測出的邊緣也不一定代表實際邊緣。圖像的邊緣有方向和幅度兩個屬性，沿邊緣方向像素變化平緩，垂直于邊緣方向像素變化劇烈。邊緣上的這種變化可以用微分算子檢測出來，通常用一階或兩階導數來檢測邊緣，如下圖所以。不同的是一階導數認為最大值對應邊緣位置，而二階導數則以過零點對應邊緣位置?；谝浑A導數的邊緣檢測算子包括Roberts算子、Sobel算子、Prewit算子等，在算法實現過程中，通過2m2(Roberts算子)或者3父3模板作為核與圖像中的每個像素點做卷積和運算，然后選取合適的閾值以提取

8、邊緣。拉普拉斯邊緣檢測算子是基于二階導數的邊緣檢測算子，該算子對噪聲敏感。一種改進方式是先對圖像進行平滑處理，然后再應用二階導數的邊緣檢測算子，其代表是LOG算子。前邊介紹的邊緣檢測算子法是基于微分方法的，其依據是圖像的邊緣對應一階導數的極大值點和二階導數的過零點。Canny算子是另外一類邊緣檢測算子，它不是通過微分算子檢測邊緣，而是在滿足一定約束條件下推導出的邊緣檢測最優化算子。1.1 Roberts(羅伯特)邊緣檢測算子景物的邊緣總是以圖像中強度的突變形式出現的，所以景物邊緣包含著大量的信息。由于景物的邊緣具有十分復雜的形態，因此，最常用的邊緣檢測方法是所謂的“梯度檢測法”。設f(x,y)

9、是圖像灰度分布函數；s(x,y)是圖像邊緣的梯度值；?(x,y)是梯度的方向。1則有s(x,y)=Uf(x+n,y)-f(x,y)+f(x,y+n)-f(x,y)了2(1)(n=1,2，)中(x,y)=tan,Uf(x,y+n)-f(x,y)/f(x+n,y)-f(x,y)D(2)式(1)與式(2)可以得到圖像在(x,y)點處的梯度大小和梯度方向。將式(1)改寫為：1g(x,y)=Jf(x,y)-Jf(x+1,y+1)f+Jf(x+1,y)-Jf(x,y+1)門2(3)g(x,y麻為Roberts邊緣檢測算子。式中對f(x,y)等的平方根運算使該處理類似于人類視覺系統的發生過程。事實上Robe

10、rts邊緣檢測算子是一種利用局部差分方法尋找邊緣的算子，Robert梯度算子所采用的是對角方向相鄰兩像素值之差，所以用差分代替一階偏導，算子形式可表不如下：xf(x,y)=f(x,y)-f(x-1,y-1)«*(4)yf(x,y)=f(x-1,y)-f(x,y-1)上述算子對應的兩個2M2模板如圖(A)所示。實際應用中，圖像中的每個像素點都用這兩個模板進行卷積運算，為避免出現負值，在邊緣檢測時常提取其絕對值。10-101-10(b)圖(A)Robert算子模板1.2 Sobel(索貝爾)邊緣檢測算子該算子是由兩個卷積核gi(x,y)與g2(x,y)對原圖像f(x,y)進行卷積運算而得

11、到的。其數學表達式為：一MNMNS(x,y)=MAXF£f(m,n)g1(im,jn),£工f(m,n)g2(im,jn)(5)_m4nTm=1nT實際上Sobel邊緣算子所采用的算法是先進行加權平均，然后進行微分運算，我們可以用差分代替一階偏導，算子的計算方法如下：4f(x,y)=f(x-1,尸1十fxy1)fXt*)fx(y1,+1)f區X,+f陽(-y1,yf(x,y)=f(x-1，廣1)牙X1y+)f*)fX(y1,+1)fX(y+1,f)Hj(寸1,(6)Sobel算子垂直方向和水平方向的模板如圖(B)所示，前者可以檢測出圖像中的水平方向的邊緣，后者則可以檢測圖像

12、中垂直方向的邊緣。實際應用中，圖像中的每一個像素點都用這兩個卷積核進行卷積運算，取其最大值作為輸出。運算結果是一幅體現邊緣幅度的圖像。-1-2-1000121-101-202-301(b)(a)圖(B)Sobel算子模板1.3Prewitt(普瑞維特)邊緣檢測算子Prewitt邊緣檢測算子就是一種利用局部差分平均方法尋找邊緣的算子，它體現了三對像素點像素值之差的平均概念，因為平均能減少或消除噪聲，為此我們可以先求平均，再求差分，即利用所謂的平均差分來求梯度。用差分代替一階偏導可得算子形式如下：Xf(x,y)-lf(x1,y1)f(x,y1)f(x-1,y1),1-If(x1,y-1)f(x,y

13、-1)f(x-1,y-1)lyf(x,y)=1f(xT,yT)f(x-1,y)f(x-1,y1)-f(x1,y-1)f(x1,y)f(x1,y1)1Prewitt邊緣檢測算子的兩個模板如圖(C)所示，它的使用方法同Sobel算子一樣，圖像中的每個點都用這兩個核進行卷積，取得最大值作為輸出。Prewitt算子也產生一幅邊緣圖像。-1-1-1000111(b)圖(C)Prewitt10-1(a)10-110-1算子模板1.4Laplacian(拉普拉斯)邊緣檢測算子對于階躍狀邊緣，其二階導數在邊緣點出現過零交叉，即邊緣點兩旁的二階導數取異號，據此可以通過二階導數來檢測邊緣點。拉普拉斯邊緣檢測算子正

14、是對二維函數進行二階導數運算的標量算子，它的定義是：二2二22f(x,y)=9f(x,y)十七f(x,y)(8):x二y用差分代替二階偏導時，與前述三個一階導數算子不同，拉普拉斯算子的形式可表示如下：2f(x,y)=f(x+1,y)+f(x1,y)+f(x,y+1)+f(x,y1)4f(x,y)/f(x,y)=f(x-1,y-1)+f(x,y-1)+f(x+1,y-1)+f(x-1,y)+f(x+1,y)f(x-1,y1)f(x,y1)f(x1,y1)-8f(x,y)(9)拉普拉斯邊緣檢測算子的模板如圖(D)所示，模板的基本特征是中心位置的系數為正，其余位置的系數為負，且模板的系數之和為零。它

15、的使用方法是用圖中的兩個點陣之一作為卷積核，與原圖像進行卷積運算即可。拉普拉斯算子又是一個線性的移不變算子，它的傳遞函數在頻域空間的原點為零，因此，一個經拉普拉斯濾波過的圖像具有零平均灰度。拉普拉斯檢測模板的特點是各向同性，對孤立點及線端的檢測效果好，但邊緣方向信息丟失，對噪聲敏感，整體檢測效果不如梯度算子。因此，它很少直接用于邊緣檢測。但注意到與Sobel算子相比，對圖像進行處理時，拉普拉斯算子能使噪聲成分得到加強，對噪聲更敏感。0-10-14-10-10-1-1-1-18-1-1-1-1(a)(b)圖(D)Laplacian算子模板1.5Marr-Hildreth(馬爾)邊緣檢測算子實際應

16、用中，由于噪聲的影響，對噪聲敏感的邊緣檢測點檢測算法(如拉普拉斯算子法)可能會把噪聲當邊緣點檢測出來，而真正的邊緣點會被噪聲淹沒而未檢測出。為此Marr和Hildreth提出了馬爾算子，因為是基于高斯算子和拉普拉斯算子的，所以也稱高斯-拉普拉斯(LaplacianofGaussian,LoG)邊緣檢測算子，簡稱LoG算子。該方法是先采用高斯算子對原圖像進行平滑又降低了噪聲，孤立的噪聲點和較小的結構組織將被濾除由于平滑會導致邊緣的延展，因此在邊緣檢測時僅考慮那些具有局部最大值的點為邊緣點，這一點可以用拉普拉斯算子將邊緣點轉換成零交叉點，然后通過零交叉點的檢測來實現邊緣檢測。所謂零交叉點就是：如果

17、一個像素處的值小于一%，而此像素8-連通的各個像素都是大于仇(也是一個正數)，那么這個像素就是零交叉點。這樣還能克服拉普拉斯算子對噪聲敏感的缺點，減少了噪聲的影響。二維高斯函數為22.,、，xy、h(x,y)=exp(-)(10)2二則連續函數f(x,y)的LoG邊緣檢測算子定義為G(x,y)-42lh(x,y)f(x,y)1=J2h(x,y)*f(x,y)(11)(12)=H(x,y)*f(x,y)_22H(x,y)=、2h(x,y)=一丁exp(仃H(x,y)是一個軸對稱函數，其橫截面如圖(E)所示。所以當邊緣模糊或噪聲較大時，利用H(x,y)檢測可有效地消除尺度遠小于高斯分布因其中r2=

18、x2+y2,。是標準差。算子由于它相當的平滑，能減少噪聲的影響,過零點能提供較可靠的邊緣位置。LoG算子在(該濾波器中的高斯函數部分對圖像具有平滑作用,子仃的噪聲信號。高斯函數在空域和頻域內都具有平滑作用。該濾波器采用拉普拉斯算子2可以減少計算量。馬爾算子用到的卷積模板一般比較大(典型半徑為8-32個像素)，不過這些模板可以分解為一維卷積來快速計算。常用的LoG算子是5M5模板，如圖(F)。與其他邊緣檢測算子一樣，LoG算子也是先對邊緣做出假設，然后再這個假設下尋找邊緣像素。但LoG算子對邊緣的假設條件最少，因此它的應用范圍更廣。另外，其他邊緣檢測算子檢測得到的邊緣時不連續的，不規則的，還需要

19、連接這些邊緣，而LoG算子的結果沒有這個缺點。對于LoG算子邊緣檢測的結果可以通過高斯函數標準偏差仃來進行調整。即仃值越大，噪聲濾波效果越好，但同時也丟失了重要的邊緣信息，影響了邊緣檢測的性能；仃值越小，又有可能平滑不完全而留有太多的噪聲。因此，在不知道物體尺度和位置的情況下，很難準確確定濾波器的仃值。一般來說，使用大仃值的濾波器產生魯棒邊緣，小的仃值的濾波器產生精確定位的邊緣，兩者結合，能夠檢測出圖像的最佳邊緣。數學上以證明，馬爾算子是按零交叉檢測階躍狀邊緣的最佳算子。但在實際圖像中要注意到，高斯濾波的零交叉點不一定全部是邊緣點，還需要進一步對其真偽進行檢驗。-2-4-4-4-2-4080-

20、4-48248-4-4080-4-2-4-4-4-2圖(F)LoG算子是5x5模板1.6canny(凱尼)邊緣檢測算子前面介紹的邊緣檢測算子都是局域窗口梯度算子，由于它們對噪聲敏感，所以在處理實際圖像時效果并不是十分理想。根據邊緣檢測的有效性和定位的可靠性，Canny研究了最優邊緣檢測器所需的特性，給出了評價邊緣檢測性能優劣的三個指標： 高的準確性，在檢測的結果里應盡量多的包含真正的邊緣，而盡量少的包含假邊緣。 高的精確度，檢測到的邊緣應該在真正的邊界上。 單像素寬，要有很高的選擇性，對每個邊緣有唯一的響應。針對這二個指標，Canny提出了用于邊緣檢測的一階微分濾波器h(x)的三個最優化標準則

21、，即最大信噪比準則、最優過零點定位準則和單邊緣響應準則。具體如下：(a)信噪比準則SNRwQG(-x)h(x)dxJ；.723dx(16)式中，G(x)為邊緣函數；h(x)為帶寬為w的低通濾波器的脈沖響應；仃是高斯噪聲的均方差。(b)定位精確度準則L為邊緣的定位精度，定義如下:w.G(-x)h(x)dx-w'2h(x)dx(17)r、I_'.式中，G(x)和h(x)為G(x)和h(x)的一階導數；L是對邊緣定位精確程度的度量，L越大定位精度越高。(c)單邊緣響應準則要保證對但邊緣只有一個響應，檢測算子的脈沖響應導數的零交叉點平均距離應該滿:2h(x)dx足：DZCa(f)=n(

22、18)J(x)dx、，、式中，h(x)是h(x)的二階導數；f是進行邊緣檢測后的圖像。這三個準則是對前述邊緣檢測指標的定量描述。抑制噪聲和邊緣精確定位是無法同時得到滿足的，即邊緣檢測算法通過圖像平滑算子去除噪聲，勢必增加邊緣定位的不確定性；反之，若提高邊緣檢測算子對邊緣的敏感性，同時也提高了對噪聲的敏感性。因此，在實際應用中只能寄希望于再抑制噪聲和提高邊緣定位精度之間實現一個合理的折衷。值得慶幸的是，有一個線性算子可以在抵抗噪聲與邊緣檢測之間獲得一個最佳的折衷，這個算子就是高斯函數的一階導數。高斯函數與原圖的卷積達到了抵抗噪聲的作用，而求導數，則是檢測景物邊緣的手段。對于階躍形的邊緣，Cann

23、y推導出的最優邊緣檢測器的形狀與高斯函數的一階導數類似，因此Canny邊緣檢測器就是高斯函數的一階導數構成的。我們知道，高斯函數是圓對稱的，因此，Canny算子在邊緣方向上是對稱的，在垂直于邊緣的方向上是反對稱的。這就意味著該算子對最急劇變化方向上的邊緣特別敏感，但在沿邊緣方向上是不敏感的。1x2y2設二維高斯函數為G(x,y)=2exp(A)(19)2；一2"其中，仃是高斯函數的分布參數，可用以控制對圖像的平滑程度。最優階躍邊緣檢測算子是以卷積VG*f(x,y)為基礎的，邊緣強度為VG*f(x,y),而邊緣方向為pivG*f(x：y)r從高斯函數的定義可知，該函數是無限拖尾的，在實

24、際應用中，一般情況下是將原始模板截斷到有限尺寸No實驗表明，當N=bJ2。十1時，能夠獲得較好的邊緣檢測結果。下面給出Canny算子的具體實現。利用高斯函數的可分性，將VG的兩個濾波卷積模板分解為兩個一維的行列濾波器：其中-：G(x,y)x迨(x,y)y=kxexp(=kyexp(一2)exp(2)=l(x)h2(y)(20)2二2)exp(-2二2)=A(y)h2(x)(21)h1(x)=.kxexp(h2(x)=、.kexp(2x2"2x2二22),%(y)一樂yexp(),h2(y)=,kexp(22"2。),2-可見,“回df0:x-:G(x,y)*f令A(i,j)

25、=E；E2,a(i,j)=arctan:yEy(i,j)Ex(i,j)(22)(23)%(x)=xh(勾，1h(子)yh(功為常數。然后把這兩個模板分別與f(x,y)進行卷積，得到則A(i,j)反映邊緣強度，根據Canny的定義,區域中的最大值。這樣，來確定該點是否為邊緣點。a(i,j)為垂直于邊緣的方向。中心邊緣點為算子Gn與圖像f(x,y)的卷積在邊緣梯度方向上的就可以在每一點的梯度方向上判斷此點強度是否為其領域的最大值當一個像素滿足以下三個條件時，則被認為是圖像的邊緣點。1)該點的邊緣強度大于沿該點梯度方向的兩個相鄰像素點的邊緣強度;2）與該點梯度方向上相鄰兩點方向差小于45°

26、3）以該點為中心的3父3領域中的邊緣強度極大值小于某個閾值。此外，如果1）和2）同時被滿足，那么在梯度方向上的相鄰像素就從候選邊緣點中取消,條件3）相當于區域梯度最大值組成的閾值圖像與邊緣點進行匹配，這一過程消除了許多虛假的邊緣點。Canny邊緣檢測算子步驟如下：Stepl:用高斯濾波器對圖像進行濾波消噪；Step2：用一階偏導的有限差分來計算梯度的幅值和方向；Step3:對梯度幅值進行非極大值抑制；Step4：用雙閾值算法檢測和連接邊緣。2.各種經典邊緣檢測算子性能比較由上述邊緣算子產生的邊緣圖像看來很相似。它們看起來像一個繪畫者從圖片中做出的線條畫。Roberts算子是2M2算子，對具有陡

27、峭的低噪聲圖像響應最好。Sobel算子，Prewitt算子，Kirsch算子都是3M3算子，對灰度漸變和噪聲較多的圖像處理得較好。使用兩個掩模板組成邊緣檢測器時，通常取較大的幅度作為輸出值。這使得它們對邊緣的走向有些敏感。取它們的平方和的開方可以獲得性能更一致的全方位響應。這與真實的梯度值更接近。值得注意的是3M3的Sobel算子和Prewitt邊緣算子可擴展成八個方向，并且可以像使用Kirsch算子一樣獲得邊緣方向圖。2.1 MATLAB程序仿真為了分析上述幾種算子的效果，我們用MATLAB對lenna原始圖像以及分別加入高斯噪聲和椒鹽噪聲后的lenna圖像進行檢測。仿真實驗結果如下（程序見

28、附錄）：Lenna原始圖像（一）各算子對lenna原始圖像的處理效果如下所示:prewitt算子檢測此已rts算子檢測號ub以算子檢測canny算子檢測log算子檢測laplacian算子檢測（二）對lenna原始圖像加入高斯噪聲:加入高斯噪聲的lenna圖像各算子對加入高斯噪聲的lenna圖像的處理效果如下所示:rubsri號算子檢測號口bsl算子檢測prewitt算子檢測enriY算子檢測lug#子檢測I叩lacia奇子檢測（三）對lenna原始圖像加入椒鹽噪聲:加入椒鹽噪聲的lenna圖像各算子對加入椒鹽噪聲的lenna圖像的處理效果如下所示:2.2 實驗效果比較由于Roberts算子是

29、利用圖像的兩個對角線的相鄰像素之差進行梯度幅值的檢測，所以求得的是在差分點處梯度幅值的近似值，并且檢測水平和垂直方向邊緣的性能好于斜線方向的邊緣，檢測精度比較高，但容易丟失一部分邊緣，同時由于沒經過圖像平滑計算，因此不能抑制噪聲，但該算子對具有陡峭的低噪聲圖像響應最好。Prewitt算子和Sobel算子都是對圖像進行差分和濾波運算，僅在平滑部分的權值選擇上有些差異，因此兩者均對噪聲具有一定的抑制能力，但這種抗噪能力是通過像素平均來實現的，所以圖像產生了一定的模糊，而且還會檢測出一些偽邊緣，所以檢測精度比較低，該算子比較適合用于圖像邊緣灰度值比較尖銳且圖像噪聲比較小的情況。LoG算子首先通過高斯

30、函數對圖像進行平滑處理，因此對噪聲的抑制作用比較明顯，但同時也可能將原有的邊緣也平滑了，造成某些邊緣無法檢測到，比外高斯分布因子仃的選擇對圖像邊緣檢測效果有較大的影響，。越大，檢測到的圖像細節越豐富，但抗噪能力下降，從而出現偽邊緣，反之則抗噪能力提高，但邊緣精度下降，易丟失許多真邊緣，因此，對于不同圖像應選擇不同參數。Canny算子也采用高斯函數對圖像進行平滑處理，因此具有較強的去噪能力，但同樣存在容易平滑掉一些邊緣信息，其后所采用的一階微分算子的方向性較LoG算子要好，因此邊緣定位精度較高。該算子與其它邊緣檢測算子的不同之處在于，它使用2種不同的閾值分別檢測強邊緣和弱邊緣，并且僅當弱邊緣相連

31、時才將弱邊緣包含在輸出圖像中，因此這種方法較其它方法而言不容易被噪聲“填充”更容易檢測出真正的弱邊緣。通過對lenna圖的仿真實驗結果可以看出，該算子在上述幾種邊緣檢測算子當中效果最好。邊緣定位準確，連續性較好，虛假邊緣少且邊緣均具有單像素寬度。表一各種經典邊緣檢測算子的優缺點比較算子優缺點比較Roberts對具有陡峭的低噪聲的圖像處理效果較好，但利用Roberts算子提取邊緣的結果是邊緣比較粗，因此邊緣定位不是很準確。Sobel對灰度漸變和噪聲較多的圖像處理效果比較好，Sobel算子對邊緣定位比較準確。Prewitt對灰度漸變和噪聲較多的圖像處理效果較好。Laplacian對圖像中的階躍性邊

32、緣點定位準確，對噪聲非常敏感，丟失一部分邊緣的方向信息，造成一些不連續的檢測邊緣。LoGLoG算子經常出現雙邊緣像素邊界，而且該檢測方法對噪聲比較敏感，所以很少用LoG算子檢測邊緣，而是用來判斷邊緣像素是位于圖像的明區還zt哈區。Canny此方法不容易受噪聲的干擾，能夠檢測到真正的弱邊緣。在edge函數中，最有效的辿緣檢測方法是Canny方法。該方法的優點在于使用兩種不同的閾值分別檢測強邊緣和弱邊緣，并且僅當弱邊緣與強邊緣相連時，才將弱邊緣包含在輸出圖像中。因此，這種方法不容易被噪聲“填充”，跟容易檢測出真正的弱邊緣。3.結論目前，邊緣檢測技術在很多領域中都得到了廣泛的應用。如對醫學圖像的邊緣

33、特征提取、車牌識別、人臉檢測技術等。本文討論和比較了幾種常用的經典邊緣檢測算子?？梢钥闯鏊鼈兌疾皇蔷哂薪^對優勢的方法，在圖像的抗噪聲、圖像邊緣的定位情況、圖像部分邊緣是否可以檢測出，以及算法運行的速度等方面，各自存在優缺點。所以對于圖像的邊緣檢測，要根據不同的系統，針對不同的環境條件和要求，選擇合適的算子。希望本文總結的各種算法的優缺點及適用范圍，能為今后人們使用這些算子以及研究更好的邊緣檢測算子方法提供參考價值。參考文獻:1RafaelC.Gonzalez、RichardE.Woods.岡薩雷斯.數字圖像處理（第二版）電子工業出版社2陳傳波、金先級.數字圖彳t處理.北京：機械工業出版社，20

34、04.73KennethR.Castleman著；朱志剛、石定機等譯.數字圖彳裒處理.北京：電子工業出版社,2002.24陳天華.數字圖彳t處理.北京：清華大學出版社，2007.65曹茂永.數字圖彳t處理.北京：北京大學出版社，2007.96龔聲蓉、劉純平、王強等編著.數字圖像處理與分析.北京：清華大學出版社7朱虹.數字圖像處理基礎.北京：科學出版社，20058郎銳.數字圖像處理學VisualC+實現.北京：北京希望電子出版社，2002.129賈永紅.計算機圖像處理與分析.武漢：武漢大學出版社，200110韓曉軍.數字圖像處理技術與應用.北京：電子工業出版社，200911王耀南，李樹濤，毛建旭

35、.計算機圖像處理與識別技術.北京：高等教育出版社，200112余成波.數字圖像處理及MATLAB；現.重慶：重慶大學出版社，200313飛思科技產品研發中心.MATLAB6.5輔助圖像處理.北京：電子工業出版社，200314何斌等.VisualC+數字圖像處理（第二版）.北京：人民郵電出版社，200215孫兆林.MATLAB6.X圖像處理.北京：清華大學出版社，200216（法）麥特爾（MaitreH.）等著；孫洪譯.現代數字圖像處理.北京：電子工業出版社，200617劉洪鳴.圖像邊緣檢測技術的研究與比較.武警杭州指揮學院.高校講壇第1期.201018 AHILK.JAIN著；韓博、徐楓譯.數

36、字圖像處理基礎.北京：清華大學出版社，2006.1119許錄平.數字圖彳t處理.北京：科學出版社，200720劉直芳、王運瓊、朱敏編著.數字圖像處理與分析.北京：清華大學出版社，2006.821藍章禮、李益才、李艾星編著.數字圖像處理與圖像通信.北京：清華大學出版社，2009.322胡學龍、許開宇編著.數字圖彳t處理.北京：電子工業出版社，2006.923傅德勝、壽亦禾編著.圖形圖像處理學.南京：東南大學出版社，2001.1224賴志國等編著.Matlab圖像處理與應用（第二版）.北京：國防工業出版社，2007.425蔣先剛著.數字圖像模式識別工程軟件設計.北京：中國水利水電出版社，20082

37、6姚敏等編著.數字圖彳t處理.北京：機械工業出版社，2006.127周道炳、朱衛綱.幾種邊緣檢測算子的評估.指揮技術學院學報.第11卷（第一期），2000.228李安安.幾種圖像邊緣檢測算法的比較與展望.大眾科技.第12期，2009.1229鄭靜、梁少華、王騰.基于MATLA的真的邊緣檢測算子研究.電腦知識與技術.第6卷（第五期）2010.230 TPoggio,HVoorhees.AYuille.ARegularizedSolutiontoEdgeDetectionR.TechRepMA,RepAIM-833,MITArtificialIntellLab.1985.31 XiaoyiJian

38、g,HorstBunke.EdgeDetectioninRangeImageBasedonScanLianApproximationJ.ComputerVisionandImageUnderstanding,1999,73（2）:183-199.32 MallatSZhong.CharacterizationofsignalfrommultiscaleedgeJ.IEEETransPAMI,1992;:710-732附錄附錄(1)各算子對lenna原始圖像處理的仿真程序:functionjingdianI=imread('lenna.jpg');I=rgb2gray(I);B1

39、=edge(I,'roberts');B2=edge(I,'sobel');B3=edge(I,'prewitt');B4=edge(I,'canny');B5=edge(I,'log');h=0-10;-14-1;0-10;B6=conv2(I,h);subplot(2,3,1);imshow(B1);title('roberts算子檢測,);subplot(2,3,2);imshow(B2);title('sobel算子檢測)；subplot(2,3,3);imshow(B3);title('prewitt算子檢測,);subplot(2,3,4);imshow(B4);title('c