徐艷華，周榮亞

（陜西鐵路工程職業(yè)技術(shù)學(xué)院，陜西渭南 714000）

在數(shù)字圖像處理領(lǐng)域，圖像邊緣[1]作為數(shù)字圖像最為重要的特征之一，可以方便地表示出一個圖像區(qū)域塊的開始與結(jié)束。從數(shù)學(xué)角度來看，圖像的邊界上存在大量梯度突變的特征點，而邊緣檢測[2]的主要功能則是查找出圖像邊緣灰度不連續(xù)的部分。目前，邊緣檢測已被廣泛應(yīng)用到圖像處理、圖像分析與應(yīng)用和機(jī)器視覺[3-4]等領(lǐng)域。

近年來，邊緣檢測已成為圖像處理研究領(lǐng)域的熱點，關(guān)于邊緣檢測的算法也較多。目前常見的邊緣檢測算法都是基于微分的思想，常見的微分算子有Robert、Soble、Log 和Canny[5]等。這些微分算子對計算機(jī)性能要求低，所需計算性能和存儲空間較小。但這些算法自身性能偏低，對復(fù)雜圖像或模糊圖像進(jìn)行檢索時表現(xiàn)較差，因此需對傳統(tǒng)算法進(jìn)行改進(jìn)。

小波分析是對微分算法進(jìn)行改進(jìn)的有效方法。其是一種局部算法，在對圖像的局部進(jìn)行分析時性能較好，同時可以在時域和頻域?qū)D像特征進(jìn)行處理。該文結(jié)合小波變換[6]與數(shù)學(xué)形態(tài)學(xué)[7]邊緣檢測的優(yōu)點，首先基于Haar 小波濾波器[8]對圖像進(jìn)行分解；然后利用Canny 檢測算子對圖像的高頻部分進(jìn)行邊緣檢測，而對于圖像的低頻子分量部分，則引入形態(tài)學(xué)邊緣檢測算法。

1 基于Haar小波變換的圖像邊緣檢測

1.1 一維非均勻Haar小波原理

作為Daubechies小波家族中最簡單的一種，Haar小波亦被稱為一階Daubechies 小波。在工程應(yīng)用中，傳統(tǒng)模式的均勻Haar 小波存在一些局限性，即對眾多工程中的邊界問題適用性較差。因此該文為提升檢測算法的精確性，采用非均勻Haar 小波[9]。下面將對非均勻Haar 小波原理進(jìn)行介紹：

非均勻Haar 小波定義的引入離不開子空間的非均勻劃分。令M=2J，伸縮因子定義為j=0,1,…,J，平移因子定義為k=0,1,…,m-1,m=2j。小波的數(shù)量由i=m+k+1 所決定，其中i的上限為2M，則定義Haar 尺度函數(shù)為：

其中，χ[a,b](x)表示[a,b]上的特征函數(shù)，對于非均勻子空間的劃分，令x∈[a,b]，規(guī)定第l個子區(qū)間的長度Δxl=xl-xl-l，其中l(wèi)的取值范圍為1,2,…,2M。假定Δxl+1=qΔxl，這里q＞1 為指定常數(shù)，所有子區(qū)間長度和為：

利用Δxl+1=qΔxl，類推迭代可得到非均勻子空間劃分的網(wǎng)格點為xl=a+lql-1/q2M-1。

以上述非均勻子空間的劃分為基礎(chǔ)，可引入一維非均勻Haar 小波的定義如下：

其中，hi(x)為第i個非均勻Haar 小波函數(shù)。當(dāng)i＞1 時，上式中ξ1(i)=x(2kμ),ξ2(i)=x[(2k+1)μ],ξ3(i)=x[2(k+1)μ]。

1.2 非均勻Haar小波算法在圖像分解中的優(yōu)勢

利用小波變換可以將基礎(chǔ)圖像進(jìn)行分解，該文選用Haar小波濾波器對圖像進(jìn)行分解。對圖像分解的本質(zhì)是對積分方程進(jìn)行求解，該文使用的Haar 小波變換在積分方程求解中具有較強的優(yōu)勢[10-11]。

在方程求解中基于中點配置法和非均勻Haar小波求解積分方程，其中將配置點選取如下：

假設(shè)求解的積分方程為：

由積分方程的相關(guān)數(shù)理知識可知，方程的精確解應(yīng)為：

而所有的配置點均可由式（6）計算得到。將積分方程的非均勻Haar 小波解表示為級數(shù)的形式可得到：

式（9）中，ai式為小波系數(shù)向量a的對應(yīng)元素，將式（9）代入式（7）中可得到：

假設(shè)Gi為：

則由式（5）的非均勻Haar 小波定義可得：

將配置點xc代入矩陣方程：

通過式（14）和（15）即可求出向量a=(ai)。再將ai代入到式（9）中，即可得到方程的Haar 小波解。由計算結(jié)果可知，方程的精確解和數(shù)值解的誤差在10-5～10-4之間，由此可以說明非均勻Haar 小波在求解積分方程時誤差較小，表明Haar 小波算法在圖像分解過程中可達(dá)到較高精度。

1.3 基于Haar小波的圖像邊緣檢測算法

該文以Haar小波變換為基礎(chǔ)利用Haar小波濾波器將原圖像進(jìn)行分解。針對不同的子分量圖像，選擇不同的邊緣檢測算法，然后對不同子分量的結(jié)果進(jìn)行重構(gòu)以提升算法性能。最終提出的算法模型主要由三部分組成，分別是Haar小波分頻處理、Canny算子邊緣檢測和抗噪型形態(tài)學(xué)邊緣檢測，整體的算法流程如圖1所示。下面將對每個具體的組成部分進(jìn)行介紹：

圖1 該文算法流程

1）基于Haar 小波的分頻與重構(gòu)：該文提出算法在進(jìn)行原始輸入圖像頻率分解與最終的合成重構(gòu)時，采用的方法為一級Haar 小波分解方法，具體分解信息如圖2所示。由圖2可知，所提出算法將輸入的原始圖像分為4個子集合向量，用字母表示為A、V、H、D。

圖2 一級Haar小波分解圖

2）Canny 算子邊緣檢測：邊緣檢測操作簡單，但其需要處理的細(xì)節(jié)信息較多。通常基于圖像中形狀區(qū)域的梯度信息，找出灰度值變化較大的點，并對圖像的邊界區(qū)域進(jìn)行識別和提取，其具體實現(xiàn)過程如圖3 所示。將圖像作為二維數(shù)據(jù)，實際上就是對二維數(shù)據(jù)取梯度值即可得到所需的圖像數(shù)據(jù)矢量值。而圖像數(shù)據(jù)梯度值的本質(zhì)則是對二維變量的每個變量單獨求偏導(dǎo)，結(jié)果如式（16）所示：

圖3 邊緣檢測一般步驟

該文對原始圖像分解后得到的高頻子分量圖像部分使用Canny 算子作為該算法的邊緣檢測算子，文獻(xiàn)[12-16]中的結(jié)果均表示該算子的檢測效果最佳。基于Gauss 濾波器[12]的強大計算能力，可以準(zhǔn)確得出圖像的梯度信息。通過尋找基于局部梯度的最大強度邊界信息，Canny 算子能夠?qū)崿F(xiàn)邊緣信息的準(zhǔn)確檢測。針對噪聲問題，Canny 算子設(shè)置了雙閾值模式，其可以保證邊緣檢測時對圖像邊緣信息的準(zhǔn)確提取。該模式能夠顯著提高邊緣檢測的精確度，并可準(zhǔn)確提取出待測圖像的邊緣信息。閾值按照大小可以分為大閾值和小閾值，若梯度強度小于檢測大閾值的要求，則需要對大小邊緣點進(jìn)行擬合，該檢測點即確定為弱邊緣點[13]。正是因為該算子使用大小兩個閾值進(jìn)行邊緣點檢測，所以該算法的檢測性能相較于其他兩個算法更優(yōu)。

3）抗噪型形態(tài)學(xué)邊緣檢測：通過上文分析可知，雖然算法操作簡便，但檢測時的圖像噪聲對檢測結(jié)果影響較大。為解決這一問題，提升算法的抗噪能力，該文引入一種抗噪型形態(tài)學(xué)邊緣檢測算法來提升所提算法在噪聲干擾情況下的性能?；拘螒B(tài)學(xué)運算如開運算和腐蝕運算[14-16]可以較好地抑制圖像中的峰值噪聲；反之，基本形態(tài)學(xué)運算如閉運算和膨脹運算則可以較好地抑制圖像中的底谷噪聲。由此可見，基于形態(tài)學(xué)的理論對上述算子進(jìn)行處理，可大幅度增強算法的魯棒性，且圖像噪聲對結(jié)果的影響也會減小。加入抗噪因子后的算子分為如下幾種類型：

開運算型：

閉運算型：

開閉運算型：

文中在對低頻子分量圖像部分進(jìn)行邊緣檢測時，使用的并非3 種檢測算子，而是加入了抗噪因子的抗噪算子，從而使該邊緣檢測算子抗圖像噪聲干擾的能力更強，其檢測結(jié)果表明進(jìn)一步優(yōu)化了圖像的邊界區(qū)域，且能夠獲得更為良好的視覺效果。

2 實驗分析

2.1 實驗條件

為驗證該文提出算法的有效性，基于Matlab R2016a軟件平臺對算法進(jìn)行仿真實驗，測試采用Lena圖。具體仿真的軟硬件環(huán)境參數(shù)設(shè)置，如表1所示。

表1 實驗軟硬件環(huán)境配置

2.2 實驗測試與結(jié)果分析

為測試該文提出的邊緣檢測算法的性能，首先采用所提算法對標(biāo)準(zhǔn)Lena 圖像進(jìn)行邊緣檢測。為方便對比，文中給出了傳統(tǒng)邊緣檢測算子Robert、Sobel、Log 和Canny 檢測的邊緣圖像，具體的實驗結(jié)果主觀示意圖如圖4 所示。

圖4 邊緣檢測算法主觀示意圖

通過觀察圖4 可知，Robert 和Sobel 邊緣檢測算法檢測的圖像邊緣不連續(xù)且細(xì)節(jié)少；Log 邊緣檢測算法檢測的圖像噪聲較多；Canny 邊緣檢測算法得到的圖像邊緣平滑，但細(xì)節(jié)信息模糊；只有該文提出的邊緣檢測算法檢測的邊緣連續(xù)且清晰，而圖像數(shù)據(jù)也得到了較好的保留。在經(jīng)典的Lena 中，可以看到人物的臉部輪廓、頭頂?shù)拿弊舆吘壱琅f清晰。

上述實驗也稱為主觀評價法，可以通過人眼直觀地對圖片質(zhì)量進(jìn)行評價。圖像質(zhì)量的客觀評價法，是指使用一些統(tǒng)計學(xué)指標(biāo)對圖像的數(shù)據(jù)進(jìn)行客觀計算以比較算法性能。該文使用到的客觀評價指標(biāo)有均方誤差、峰值信噪比以及平均梯度，使用這3個客觀指標(biāo)對文中算法和對比算法進(jìn)行計算，其結(jié)果如表2 所示。

表2 邊緣檢測算法客觀指標(biāo)

由表2 可知，該文算法在客觀評價的3 個指標(biāo)中均領(lǐng)先于其他3 種算法。綜上所述，通過主觀評測對比與客觀指標(biāo)評測對比可以發(fā)現(xiàn)文中提出的算法具有更強的抗噪聲干擾能力，在性能表現(xiàn)上優(yōu)于傳統(tǒng)邊緣檢測算法。

3 結(jié)束語

該文提出的邊緣檢測算法基于Haar 小波變換將原始圖像進(jìn)行分頻處理。使用非均勻Haar 小波變換，其積分方程的解相較傳統(tǒng)小波變換精度更高，因此提出算法具有更優(yōu)的檢測性能。與其他圖像邊緣檢測的算法相比較而言，文中提出算法的優(yōu)勢在于顯著提升了抗圖像噪聲干擾的能力，且檢測出的結(jié)果邊緣清晰，圖像邊界被進(jìn)一步凸顯出來，主觀視覺效果更優(yōu)。同時，對于背景對比度較低的圖像類型，該算法依舊能對邊緣進(jìn)行高質(zhì)量的檢測與提取?；贚ena 圖像對該文算法的性能進(jìn)行驗證，在主觀性驗證中表現(xiàn)良好且圖像邊緣細(xì)節(jié)較為清晰。在客觀性對比評價中指標(biāo)領(lǐng)先于對比算法，說明文中算法對圖像的檢測性能較好，圖像的細(xì)節(jié)表現(xiàn)力和抗噪能力更強。