Python圖像銳化與邊緣檢測之Sobel與Laplacian算子詳解
一.Sobel算子
Sobel算子是一種用于邊緣檢測的離散微分算子,它結(jié)合了高斯平滑和微分求導(dǎo)。該算子用于計算圖像明暗程度近似值,根據(jù)圖像邊緣旁邊明暗程度把該區(qū)域內(nèi)超過某個數(shù)的特定點記為邊緣。Sobel算子在Prewitt算子的基礎(chǔ)上增加了權(quán)重的概念,認(rèn)為相鄰點的距離遠近對當(dāng)前像素點的影響是不同的,距離越近的像素點對應(yīng)當(dāng)前像素的影響越大,從而實現(xiàn)圖像銳化并突出邊緣輪廓[1-4]。
Sobel算子的邊緣定位更準(zhǔn)確,常用于噪聲較多、灰度漸變的圖像。其算法模板如公式(1)所示,其中dx表示水平方向,dy表示垂直方向[3]。
其像素計算公式如下:
Sobel算子像素的最終計算公式如下:
Sobel算子根據(jù)像素點上下、左右鄰點灰度加權(quán)差,在邊緣處達到極值這一現(xiàn)象檢測邊緣。對噪聲具有平滑作用,提供較為精確的邊緣方向信息。因為Sobel算子結(jié)合了高斯平滑和微分求導(dǎo)(分化),因此結(jié)果會具有更多的抗噪性,當(dāng)對精度要求不是很高時,Sobel算子是一種較為常用的邊緣檢測方法。
Python和OpenCV將Sobel算子封裝在Sobel()函數(shù)中,其函數(shù)原型如下所示:
dst = Sobel(src, ddepth, dx, dy[, dst[, ksize[, scale[, delta[, borderType]]]]])
- – src表示輸入圖像
- – dst表示輸出的邊緣圖,其大小和通道數(shù)與輸入圖像相同
- – ddepth表示目標(biāo)圖像所需的深度,針對不同的輸入圖像,輸出目標(biāo)圖像有不同的深度
- – dx表示x方向上的差分階數(shù),取值1或 0
- – dy表示y方向上的差分階數(shù),取值1或0
- – ksize表示Sobel算子的大小,其值必須是正數(shù)和奇數(shù)
- – scale表示縮放導(dǎo)數(shù)的比例常數(shù),默認(rèn)情況下沒有伸縮系數(shù)
- – delta表示將結(jié)果存入目標(biāo)圖像之前,添加到結(jié)果中的可選增量值
- – borderType表示邊框模式,更多詳細(xì)信息查閱BorderTypes
注意,在進行Sobel算子處理之后,還需要調(diào)用convertScaleAbs()函數(shù)計算絕對值,并將圖像轉(zhuǎn)換為8位圖進行顯示。其算法原型如下:
dst = convertScaleAbs(src[, dst[, alpha[, beta]]])
- – src表示原數(shù)組
- – dst表示輸出數(shù)組,深度為8位
- – alpha表示比例因子
- – beta表示原數(shù)組元素按比例縮放后添加的值
Sobel算子的實現(xiàn)代碼如下所示。
# -*- coding: utf-8 -*- # By:Eastmount import cv2 import numpy as np import matplotlib.pyplot as plt #讀取圖像 img = cv2.imread('luo.png') lenna_img = cv2.cvtColor(img,cv2.COLOR_BGR2RGB) #灰度化處理圖像 grayImage = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY) #Sobel算子 x = cv2.Sobel(grayImage, cv2.CV_16S, 1, 0) #對x求一階導(dǎo) y = cv2.Sobel(grayImage, cv2.CV_16S, 0, 1) #對y求一階導(dǎo) absX = cv2.convertScaleAbs(x) absY = cv2.convertScaleAbs(y) Sobel = cv2.addWeighted(absX, 0.5, absY, 0.5, 0) #用來正常顯示中文標(biāo)簽 plt.rcParams['font.sans-serif']=['SimHei'] #顯示圖形 titles = ['原始圖像', 'Sobel算子'] images = [lenna_img, Sobel] for i in range(2): plt.subplot(1,2,i+1), plt.imshow(images[i], 'gray') plt.title(titles[i]) plt.xticks([]),plt.yticks([]) plt.show()
其運行結(jié)果如圖1所示:
二.Laplacian算子
拉普拉斯(Laplacian)算子是n維歐幾里德空間中的一個二階微分算子,常用于圖像增強領(lǐng)域和邊緣提取。它通過灰度差分計算鄰域內(nèi)的像素,基本流程是:
- 判斷圖像中心像素灰度值與它周圍其他像素的灰度值;
- 如果中心像素的灰度更高,則提升中心像素的灰度;
- 反之降低中心像素的灰度,從而實現(xiàn)圖像銳化操作。
在算法實現(xiàn)過程中,Laplacian算子通過對鄰域中心像素的四方向或八方向求梯度,再將梯度相加起來判斷中心像素灰度與鄰域內(nèi)其他像素灰度的關(guān)系,最后通過梯度運算的結(jié)果對像素灰度進行調(diào)整[2]。
一個連續(xù)的二元函數(shù)f(x,y),其拉普拉斯運算定義為:
Laplacian算子分為四鄰域和八鄰域,四鄰域是對鄰域中心像素的四方向求梯度,八鄰域是對八方向求梯度。其中,四鄰域模板如公式(5)所示:
其像素的計算公式可以簡化為:
通過模板可以發(fā)現(xiàn),當(dāng)鄰域內(nèi)像素灰度相同時,模板的卷積運算結(jié)果為0;當(dāng)中心像素灰度高于鄰域內(nèi)其他像素的平均灰度時,模板的卷積運算結(jié)果為正數(shù);當(dāng)中心像素的灰度低于鄰域內(nèi)其他像素的平均灰度時,模板的卷積為負(fù)數(shù)。對卷積運算的結(jié)果用適當(dāng)?shù)乃ト跻蜃犹幚聿⒓釉谠行南袼厣?,就可以實現(xiàn)圖像的銳化處理。
Laplacian算子的八鄰域模板如下:
其像素的計算公式可以簡化為:
Python和OpenCV將Laplacian算子封裝在Laplacian()函數(shù)中,其函數(shù)原型如下所示:
dst = Laplacian(src, ddepth[, dst[, ksize[, scale[, delta[, borderType]]]]])
- – src表示輸入圖像
- – dst表示輸出的邊緣圖,其大小和通道數(shù)與輸入圖像相同
- – ddepth表示目標(biāo)圖像所需的深度
- – ksize表示用于計算二階導(dǎo)數(shù)的濾波器的孔徑大小,其值必須是正數(shù)和奇數(shù),且默認(rèn)值為1,更多詳細(xì)信息查閱getDerivKernels
- – scale表示計算拉普拉斯算子值的可選比例因子。默認(rèn)值為1,更多詳細(xì)信息查閱getDerivKernels
- – delta表示將結(jié)果存入目標(biāo)圖像之前,添加到結(jié)果中的可選增量值,默認(rèn)值為0
- – borderType表示邊框模式,更多詳細(xì)信息查閱BorderTypes
注意,Laplacian算子其實主要是利用Sobel算子的運算,通過加上Sobel算子運算出的圖像x方向和y方向上的導(dǎo)數(shù),得到輸入圖像的圖像銳化結(jié)果。
同時,在進行Laplacian算子處理之后,還需要調(diào)用convertScaleAbs()函數(shù)計算絕對值,并將圖像轉(zhuǎn)換為8位圖進行顯示。其算法原型如下:
dst = convertScaleAbs(src[, dst[, alpha[, beta]]])
- – src表示原數(shù)組
- – dst表示輸出數(shù)組,深度為8位
- – alpha表示比例因子
- – beta表示原數(shù)組元素按比例縮放后添加的值
當(dāng)ksize=1時,Laplacian()函數(shù)采用3×3的孔徑(四鄰域模板)進行變換處理。下面的代碼是采用ksize=3的Laplacian算子進行圖像銳化處理,其代碼如下:
其運行結(jié)果如圖2所示:
# -*- coding: utf-8 -*- # By:Eastmount import cv2 import numpy as np import matplotlib.pyplot as plt #讀取圖像 img = cv2.imread('luo.png') lenna_img = cv2.cvtColor(img,cv2.COLOR_BGR2RGB) #灰度化處理圖像 grayImage = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY) #拉普拉斯算法 dst = cv2.Laplacian(grayImage, cv2.CV_16S, ksize = 3) Laplacian = cv2.convertScaleAbs(dst) #用來正常顯示中文標(biāo)簽 plt.rcParams['font.sans-serif']=['SimHei'] #顯示圖形 titles = ['原始圖像', 'Laplacian算子'] images = [lenna_img, Laplacian] for i in range(2): plt.subplot(1,2,i+1), plt.imshow(images[i], 'gray') plt.title(titles[i]) plt.xticks([]),plt.yticks([]) plt.show()
邊緣檢測算法主要是基于圖像強度的一階和二階導(dǎo)數(shù),但導(dǎo)數(shù)通常對噪聲很敏感,因此需要采用濾波器來過濾噪聲,并調(diào)用圖像增強或閾值化算法進行處理,最后再進行邊緣檢測。下面是采用高斯濾波去噪和閾值化處理之后,再進行邊緣檢測的過程,并對比了四種常見的邊緣提取算法。
# -*- coding: utf-8 -*- # By:Eastmount import cv2 import numpy as np import matplotlib.pyplot as plt #讀取圖像 img = cv2.imread('luo.png') lenna_img = cv2.cvtColor(img, cv2.COLOR_BGR2RGB) #灰度化處理圖像 grayImage = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY) #高斯濾波 gaussianBlur = cv2.GaussianBlur(grayImage, (3,3), 0) #閾值處理 ret, binary = cv2.threshold(gaussianBlur, 127, 255, cv2.THRESH_BINARY) #Roberts算子 kernelx = np.array([[-1,0],[0,1]], dtype=int) kernely = np.array([[0,-1],[1,0]], dtype=int) x = cv2.filter2D(binary, cv2.CV_16S, kernelx) y = cv2.filter2D(binary, cv2.CV_16S, kernely) absX = cv2.convertScaleAbs(x) absY = cv2.convertScaleAbs(y) Roberts = cv2.addWeighted(absX, 0.5, absY, 0.5, 0) #Prewitt算子 kernelx = np.array([[1,1,1],[0,0,0],[-1,-1,-1]], dtype=int) kernely = np.array([[-1,0,1],[-1,0,1],[-1,0,1]], dtype=int) x = cv2.filter2D(binary, cv2.CV_16S, kernelx) y = cv2.filter2D(binary, cv2.CV_16S, kernely) absX = cv2.convertScaleAbs(x) absY = cv2.convertScaleAbs(y) Prewitt = cv2.addWeighted(absX,0.5,absY,0.5,0) #Sobel算子 x = cv2.Sobel(binary, cv2.CV_16S, 1, 0) y = cv2.Sobel(binary, cv2.CV_16S, 0, 1) absX = cv2.convertScaleAbs(x) absY = cv2.convertScaleAbs(y) Sobel = cv2.addWeighted(absX, 0.5, absY, 0.5, 0) #拉普拉斯算法 dst = cv2.Laplacian(binary, cv2.CV_16S, ksize = 3) Laplacian = cv2.convertScaleAbs(dst) #效果圖 titles = ['Source Image', 'Binary Image', 'Roberts Image', 'Prewitt Image','Sobel Image', 'Laplacian Image'] images = [lenna_img, binary, Roberts, Prewitt, Sobel, Laplacian] for i in np.arange(6): plt.subplot(2,3,i+1),plt.imshow(images[i],'gray') plt.title(titles[i]) plt.xticks([]),plt.yticks([]) plt.show()
輸出結(jié)果如圖3所示。其中,Laplacian算子對噪聲比較敏感,由于其算法可能會出現(xiàn)雙像素邊界,常用來判斷邊緣像素位于圖像的明區(qū)或暗區(qū),很少用于邊緣檢測;Robert算子對陡峭的低噪聲圖像效果較好,尤其是邊緣正負(fù)45度較多的圖像,但定位準(zhǔn)確率較差;Prewitt算子對灰度漸變的圖像邊緣提取效果較好,而沒有考慮相鄰點的距離遠近對當(dāng)前像素點的影響;Sobel算子考慮了綜合因素,對噪聲較多的圖像處理效果更好。
三.總結(jié)
本文主要介紹圖像銳化和邊緣檢測知識,詳細(xì)講解了Sobel算子和Laplacian算子,并通過小珞珞圖像進行邊緣輪廓提取。圖像銳化和邊緣提取技術(shù)可以消除圖像中的噪聲,提取圖像信息中用來表征圖像的一些變量,為圖像識別提供基礎(chǔ)。
到此這篇關(guān)于Python圖像銳化與邊緣檢測之Sobel與Laplacian算子詳解的文章就介紹到這了,更多相關(guān)Python Sobel Laplacian內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!
相關(guān)文章
VSCode下配置python調(diào)試運行環(huán)境的方法
這篇文章主要介紹了VSCode下配置python調(diào)試運行環(huán)境的方法,需要的朋友可以參考下2018-04-04python 與GO中操作slice,list的方式實例代碼
這篇文章主要介紹了python 與GO中操作slice,list的方式實例代碼的相關(guān)資料,需要的朋友可以參考下2017-03-03利用Python創(chuàng)建API服務(wù)器并處理RESTful請求
在軟件開發(fā)實踐中,構(gòu)建API服務(wù)器是一項基礎(chǔ)且重要的任務(wù),本文將介紹如何使用Python中的Flask框架創(chuàng)建一個API服務(wù)器,并展示如何處理不同的RESTful請求方法,感興趣的小伙伴可以了解下2024-02-02python高效過濾出文件夾下指定文件名結(jié)尾的文件實例
今天小編就為大家分享一篇python高效過濾出文件夾下指定文件名結(jié)尾的文件實例,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧2018-10-10