劉繼承 張 琳 董青松 石祥華
(1.東北石油大學(xué),黑龍江 大慶 163318;2.大慶油田天然氣分公司,黑龍江 大慶 163000)
邊緣檢測技術(shù)在圖像分析與識別領(lǐng)域至關(guān)重要,尋找一種最優(yōu)邊緣檢測算法是近年來此領(lǐng)域的熱門課題之一。邊緣檢測是通過檢測圖像區(qū)域間的變化實現(xiàn)的,這些變化有灰度、顏色及紋理特征等。通常是將彩色圖像轉(zhuǎn)換為灰度圖像,再利用傳統(tǒng)的邊緣檢測算子(如一階微分算子——Prewitt算子、Roberts算子和Sobel算子,及二階微分算子——Laplace算子等),計算梯度的幅值和方向。這些算子易于實現(xiàn),而且實時性好,但是在抗噪和單像素邊緣定位方面效果不太好,尤其是對噪聲的干擾非常敏感。上述算法中,基于多尺度小波分析的圖像邊緣檢測算法在抑制噪聲和提取邊緣定位方面表現(xiàn)良好,受到了國內(nèi)外學(xué)者的廣泛重視[1]。1983年有學(xué)者提出尺度空間思想,對邊緣檢測中的多尺度和多分辨進(jìn)行了深入研究。1992年又有學(xué)者提出小波變換多尺度邊緣檢測方法,同時將小波邊緣檢測方法與LOG算子和Canny最優(yōu)檢測算子在小波意義下統(tǒng)一,更明確地表達(dá)了多尺度思想在邊緣檢測中的重要意義[1]。
多尺度小波變換圖像邊緣檢測算法存在的問題是大尺度下圖像邊緣細(xì)節(jié)信息會有丟失,導(dǎo)致檢測到邊緣點定位不準(zhǔn)確而偏離實際邊緣;小尺度下圖像邊緣點又易受噪聲影響導(dǎo)致檢測到偽邊緣[2]。選取合適的尺度并采用合理的邊緣細(xì)化算法是平衡多尺度下邊緣檢測優(yōu)/缺點的主要思路。筆者在多尺度小波變換得到的邊緣灰度圖像中選擇最優(yōu)尺度灰度圖像,再選擇合適的灰度門限,將圖像轉(zhuǎn)換為二值圖像后,用邊緣細(xì)化算法將邊緣準(zhǔn)確定位到單像素寬度,以期較傳統(tǒng)邊緣檢測算子在噪聲容限和單像素邊緣定位方面有所提高。
對于一幅普通圖像而言,邊緣就是圖像顏色或灰度發(fā)生突變之處,邊緣檢測的目的就是檢測出圖像中突變發(fā)生的位置。傳統(tǒng)的邊緣檢測算子對圖像進(jìn)行邊緣檢測的基本原理就是求微分運(yùn)算,找到圖像中的突變點作為圖像邊緣,原理簡單、操作容易,但是存在的最大問題就是對噪聲敏感。
較好地解決噪聲敏感問題的方法之一是對圖像進(jìn)行多尺度邊緣檢測,這需要首先對圖像進(jìn)行多尺度表達(dá)。對圖像進(jìn)行小波變換是獲得圖像多尺度表達(dá)的一種方法,小波變換能夠把圖像分解成多種尺度成分,并對大小不同的尺度成分采用相應(yīng)的時域或空域取樣步長,從而不斷地聚焦到圖像的任意微小細(xì)節(jié)、間斷點、奇異點和邊緣[3]。小波變換所獨(dú)有的這種多尺度特性恰好可以用于多尺度下圖像的邊緣檢測[4]。人在觀察目標(biāo)時,隨著分辨率的增加,越走近目標(biāo)觀察獲得的信息就越豐富;反之,獲得的信息越少[5]。在噪聲抑制與圖像細(xì)節(jié)表現(xiàn)之間尋求平衡是多尺度分析的主要任務(wù)。由于小波變換在各尺度上都提供了圖像的邊緣信息,所以稱為多尺度邊緣,小波變換是提取多尺度邊緣的基礎(chǔ)。
小波變換是將一個函數(shù)表示為一個尺度成分和這個尺度下的一個小的波動成分的一種運(yùn)算。
對每個像素點(n,m)計算相角Af(2j,n,m)的正切值:
確定閾值T>0,對于Mf(2j,n,m),如果Mf(2j,n,m)≥T,Mf(2j,n,m)取得局部最大值,即此時的(n,m)為模極大值點。因為梯度的局部極大值對應(yīng)著圖像f(n,m)的銳變處,而且圖像的邊緣也處于圖像的銳變處,所以圖像f(n,m)的梯度局部極大值點就對應(yīng)著圖像的邊緣點[7]。
局部梯度幅值最大就是圖像的局部高頻信息,圖像的噪聲也是局部高頻信息,經(jīng)過小波變換后也可能產(chǎn)生邊緣,這種邊緣為偽邊緣。圖像進(jìn)行小波變換后無論在哪個尺度上,尖銳邊緣都有很大的信號值;相反,噪聲點的信號值會隨著尺度的增大而衰減[8]。
從圖1所示的實驗結(jié)果可以看出:隨著尺度的增大,圖像的細(xì)節(jié)逐漸減少消失,邊緣逐漸平滑。這正是多尺度小波變換對高頻信息進(jìn)行平滑濾波的體現(xiàn)。但同時也存在圖像細(xì)節(jié)丟失的缺陷,直觀感覺就是圖像變模糊。
圖1 實驗結(jié)果
經(jīng)小波邊緣提取得到的圖像邊緣比較粗,難以達(dá)到單像素的精度,需要對圖像的邊緣進(jìn)行細(xì)化。圖像邊緣細(xì)化是圖像處理中的基本技術(shù),它要求完整地保存圖像的拓?fù)浣Y(jié)構(gòu),以便于代替原始圖像進(jìn)行識別和處理。將一個圖像的主要邊緣清晰、完整地提取出來,得到細(xì)化的、完整的邊緣,將為圖像檢索、目標(biāo)分割及識別等后續(xù)處理帶來極大的便利[9]。
此處的邊緣細(xì)化算法采用改進(jìn)的Zhang-suen算法,其基本思想是在細(xì)化過程中不斷移動3×3的模板,使它與圖像中的各點重合。在掃描圖像的過程中,定義待刪除的像素為P1,P1周圍的近鄰像素位置關(guān)系如圖2所示。
定義n(P1)和s(P1)兩個參數(shù),其中n(P1)是與P1相鄰的非零像素個數(shù),s(P1)是沿著P2→P3→P4→P5→P6→P7→P8→P9→P2的順序由0過渡到1的總次數(shù)。對于考察點P1=1,如果同時滿足2≤n(P1)≤6、s(P1)=1、P2P4P8=0且P2P6P8=0這4個條件,則可將P1刪除。細(xì)化過程是將滿足條件的像素點予以刪除。反復(fù)迭代直到再也沒有像素點滿足上述條件為止,此時完成檢驗,得到了細(xì)化后的圖像邊緣。這樣,所有檢測到的邊緣可以定位到單像素寬度,完整地保留了圖像的邊緣信息[10]。
圖3給出了不存在噪聲時,傳統(tǒng)邊緣檢測算子檢測的結(jié)果。
圖3 無噪聲圖像傳統(tǒng)邊緣檢測結(jié)果
由圖3可以看出,傳統(tǒng)的邊緣檢測算子在沒有噪聲的情況下可以得到較好的邊緣檢測效果。實際的數(shù)字圖像總是存在一定功率的噪聲,如椒鹽噪聲及高斯噪聲等,其中椒鹽噪聲可以通過中值濾波較好地去除,高斯噪聲是圖像噪聲容限方面要考慮的主要噪聲。圖4給出了存在高斯噪聲(μ=0,σ=0.01)時傳統(tǒng)的邊緣檢測算子的檢測的結(jié)果。
圖4 傳統(tǒng)方法檢測的含高斯噪聲圖像結(jié)果
可以看出,在有高斯噪聲存在的情況下,傳統(tǒng)的邊緣檢測方法檢測的圖像邊緣出現(xiàn)缺失,偽邊緣被檢測出大范圍出現(xiàn),對邊緣檢測結(jié)果造成了極大破壞。采用基于二進(jìn)小波變換的多尺度邊緣細(xì)化檢測方法檢測得到的邊緣如圖5所示。
從視覺直觀感受來看,加入高斯噪聲的傳統(tǒng)邊緣檢測方法的檢測結(jié)果受噪聲干擾很大;而采用筆者提出的方法后,圖像邊緣的連續(xù)性和準(zhǔn)確性都比傳統(tǒng)方法的檢測結(jié)果有較大改善。
圖5 檢測結(jié)果
在此,定量地來探討圖像邊緣檢測結(jié)果的比較。邊緣檢測過程會產(chǎn)生3個主要類型的誤差,即邊緣點丟失、邊緣點定位偏離以及將噪聲波動誤認(rèn)為是邊緣點。以Canny最優(yōu)邊緣檢測準(zhǔn)則為參考依據(jù),Canny最佳邊緣檢測三準(zhǔn)則具體如下:
a. 最優(yōu)檢測。對于真實存在的邊緣不漏檢,當(dāng)然也不會把非邊緣點檢出,使得輸出信噪比最大。
b. 最優(yōu)檢測精度。所得邊緣點的位置與實際邊緣點的位置最近。
c. 檢測點與邊緣點一一對應(yīng)。每個實際存在的邊緣點與檢測到的邊緣點一一對應(yīng)。
信噪比SNR越大,能夠準(zhǔn)確檢測到的真實邊角和較少含有偽邊界的可能性也就越大。檢測精度L為檢測到的邊界與真實邊界之間倒數(shù)的數(shù)學(xué)期望,檢測精度越高,測量誤差越小。偽邊界平均距離M為隨機(jī)噪聲與檢測函數(shù)卷積之后偽邊界出現(xiàn)的平均距離,偽邊界平均距離越長測量結(jié)果中出現(xiàn)偽邊界的個數(shù)就越少[12]。
以主要性能指標(biāo)信噪比SNR作為邊緣檢測評價標(biāo)準(zhǔn),將沒有噪聲存在時的Canny邊緣檢測算子檢測的邊緣作為實際邊緣,與存在(μ=0,σ=0.01)高斯噪聲時傳統(tǒng)的最優(yōu)邊緣檢測算法Canny邊緣檢測和筆者所提方法進(jìn)行比較,Canny邊緣檢測的SNR=1.005 0,小波多尺度邊緣檢測的SNR=1.534 8。
筆者研究了基于二進(jìn)小波變換的圖像多尺度邊緣檢測及其邊緣細(xì)化算法,在與傳統(tǒng)邊緣檢測方法進(jìn)行比較的基礎(chǔ)上,對傳統(tǒng)方法在無噪聲和有噪聲存在時邊緣檢測的結(jié)果進(jìn)行比較,并就筆者所提方法與傳統(tǒng)邊緣檢測方法的標(biāo)信噪比進(jìn)行了定量分析。實驗結(jié)果表明:基于二進(jìn)小波變換的圖像邊緣提取細(xì)化算法較傳統(tǒng)邊緣檢測方法在噪聲容限和邊緣細(xì)化到單像素寬度方面有明顯改善。與傳統(tǒng)邊緣算法相比,小波變換邊緣檢測細(xì)化算法具有定位精確和計算量小的優(yōu)點,對噪聲的敏感程度也遠(yuǎn)低于傳統(tǒng)邊緣算法。