最近一直有人向我提問很多二值圖像分析相關(guān)的問題，特別選擇了兩個(gè)典型的輪廓分析問題。進(jìn)行分析與編碼實(shí)現(xiàn)與演示，廢話不多說，先看第一個(gè)問題。

問題一

描述如下：

想找到工具盤中間缺少的幾個(gè)點(diǎn)，統(tǒng)計(jì)出可以看到的工件數(shù)目

仔細(xì)分析圖像發(fā)現(xiàn)，中間都毫無另外的有個(gè)白色很亮的圓圈，這個(gè)給了我兩個(gè)思路

可以通過霍夫變換檢測(cè)圓來提取到

可以通過二值圖像分析來提取 + 輪廓分析來提取到這些點(diǎn)

得到這些輪廓點(diǎn)之后通過分析整個(gè)輪廓區(qū)域得到傾斜角度，進(jìn)行糾偏，然后通過X與Y投影進(jìn)行分割，得到每個(gè)零件的中心位置坐標(biāo)，根據(jù)每一行的間隔設(shè)置閾值，從而實(shí)現(xiàn)缺少部分部分的標(biāo)出與件數(shù)統(tǒng)計(jì)，確定了這樣的思路以后，我就開始了寫代碼。代碼實(shí)現(xiàn)是基于輪廓分析的思路，因?yàn)檫@個(gè)方法，用的閾值比較少，有利于算法穩(wěn)定性檢測(cè)。演示各部輸出。二值化處理之后(形態(tài)學(xué)處理)：

輪廓發(fā)現(xiàn)與校正角度之后

投影分析與統(tǒng)計(jì)結(jié)果如下：

此外基于霍夫也是可以嘗試的，霍夫的二值化效果也比較好，顯示如下：

感興趣的同學(xué)可以自己繼續(xù)嘗試下去。

問題二

描述如下：

如何統(tǒng)計(jì)下圖中的對(duì)象個(gè)數(shù)，原圖如下

看到這個(gè)圖像之后，個(gè)人覺得解決十分簡(jiǎn)單，基于最外層輪廓發(fā)現(xiàn)即可，無需樹形結(jié)構(gòu)與層次分析，集合圖像形態(tài)學(xué)分析或者距離變換就可以得到，最終代碼的運(yùn)行結(jié)果如下：

代碼

問題1的代碼如下(已經(jīng)添加各步驟注釋了)：

src=cv.imread("D:/images/zsxq/zsxq_01.jpg")
cv.imshow("input",src)

#二值化處理
gray=cv.cvtColor(src,cv.COLOR_BGR2GRAY)
se=cv.getStructuringElement(cv.MORPH_RECT,(3,3))
dst=cv.morphologyEx(gray,cv.MORPH_GRADIENT,se)
ret,binary=cv.threshold(dst,0,255,cv.THRESH_OTSU|cv.THRESH_BINARY)

#形態(tài)學(xué)處理
se=cv.getStructuringElement(cv.MORPH_ELLIPSE,(5,5))
binary=cv.morphologyEx(binary,cv.MORPH_OPEN,se)
se=cv.getStructuringElement(cv.MORPH_ELLIPSE,(10,10))
binary=cv.morphologyEx(binary,cv.MORPH_CLOSE,se)
cv.imshow("binary",binary)

#輪廓分析
contours,hireachy=cv.findContours(binary,cv.RETR_EXTERNAL,cv.CHAIN_APPROX_SIMPLE)
result=np.zeros_like(binary)
total=0
forcntinrange(len(contours)):
area=cv.contourArea(contours[cnt])
ifarea

	

	問題2的完整代碼如下：

	
importcv2ascv
importnumpyasnp

src=cv.imread("D:/images/zsxq/zsxq_02.jpg")
cv.imshow("input",src)
src=cv.GaussianBlur(src,(3,3),0)
gray=cv.cvtColor(src,cv.COLOR_BGR2GRAY)
edge=cv.Canny(src,50,100)

se=cv.getStructuringElement(cv.MORPH_ELLIPSE,(10,10))
binary=cv.morphologyEx(edge,cv.MORPH_CLOSE,se)
contours,hireachy=cv.findContours(binary,cv.RETR_EXTERNAL,cv.CHAIN_APPROX_SIMPLE)
count=0
forcntinrange(len(contours)):
area=cv.contourArea(contours[cnt])
ifarea

	

	

	審核編輯：湯梓紅