張冬妍 張瑞 韓睿 曹軍

(東北林業(yè)大學(xué)，哈爾濱，150040)

榛子形似板栗，果仁口感香酥、營(yíng)養(yǎng)價(jià)值高，被贊譽(yù)為“堅(jiān)果之王”。作為林業(yè)副產(chǎn)品中的優(yōu)質(zhì)資源，在烘焙食用、油料加工等方面的需求量極大[1-3]。受生長(zhǎng)環(huán)境、保存條件等影響，榛子仁中常有干癟、蟲(chóng)蛀、發(fā)霉、腐爛等缺陷。以往對(duì)其分類識(shí)別，主要依靠人工分揀，但人工分揀效率低，且評(píng)判標(biāo)準(zhǔn)不一，很容易出現(xiàn)檢測(cè)結(jié)果的差異。且榛子仁個(gè)體間形狀大致相近，肉眼可見(jiàn)差別不大；由于產(chǎn)地、種類及采摘期不同等原因，對(duì)于榛子仁的識(shí)別、分類難度會(huì)進(jìn)一步加大。采用圖像處理技術(shù)[4-7]，可提高榛子仁缺陷識(shí)別效率和準(zhǔn)確率，減少人力資源的浪費(fèi)。

由于榛子仁表面種皮粗糙，圖像邊界模糊，采用傳統(tǒng)的閾值分割難以清晰地將樣本與背景分開(kāi)。而模糊C均值聚類(FCM)算法選取的聚類值具有非概率特性，聚類結(jié)果靈活[8]，可以把目標(biāo)與背景間邊界模糊不清的榛子仁個(gè)體，完整地從整體圖像中提取出來(lái)。但傳統(tǒng)的模糊C均值聚類算法對(duì)噪聲不敏感。本研究以榛子仁為檢測(cè)樣本，采用模糊C均值聚類算法進(jìn)行圖像分割；利用飛蛾撲火(MFO)算法[9-11]改進(jìn)其目標(biāo)函數(shù)；利用函數(shù)對(duì)個(gè)體樣本邊緣提取[12-14]，標(biāo)記邊緣拐點(diǎn)位置，計(jì)算拐點(diǎn)個(gè)數(shù)；對(duì)邊緣圖像進(jìn)行霍夫(Hough)變換的橢圓曲線擬合，標(biāo)記并輸出飽滿籽粒個(gè)數(shù)；依據(jù)試驗(yàn)數(shù)據(jù)，分析應(yīng)用改進(jìn)的模糊C均值聚類算法和霍夫變換對(duì)榛子仁缺陷檢測(cè)的效果。旨在為準(zhǔn)確有效地對(duì)榛子仁中的缺陷籽粒進(jìn)行識(shí)別檢測(cè)、提高榛子仁加工過(guò)程中的分揀效率提供參考。

1 材料與方法

1.1 試驗(yàn)材料

榛子仁：試驗(yàn)樣本選自大興安嶺地區(qū)的野生榛子。根據(jù)樣本實(shí)際情況，結(jié)合表1的定義標(biāo)準(zhǔn)[15-16]，將其劃分為飽滿、干癟、霉斑、蟲(chóng)蛀、腐爛5個(gè)不同種類，除飽滿粒外，其他種類均為缺陷。從已有1 000顆榛子仁中，分類選取350顆(包括飽滿粒100顆、干癟粒100顆、霉斑粒50顆、蟲(chóng)蛀粒50顆、腐爛粒50顆)，作為試驗(yàn)樣本進(jìn)行圖像采集。

表1 堅(jiān)果類食品術(shù)語(yǔ)和定義

1.2 模糊C均值聚類算法

模糊C均值聚類(FCM)算法[17]由Bezdek在1993年首次提出，主要是通過(guò)隸屬度優(yōu)化目標(biāo)函數(shù)，定量表示樣本間關(guān)系，盡可能減小類內(nèi)差別，增大類間差別，以便自動(dòng)對(duì)樣本數(shù)據(jù)進(jìn)行分類。先將該算法的聚類類別定義為C，當(dāng)中包括N個(gè)L維向量數(shù)據(jù)，用x(k)表示，表示數(shù)據(jù)元素xj被劃分為第i類的概率，即元素隸屬度記作uij，可以通過(guò)式(1)求最小值。

i=1、2、3、…、C，j=1、2、3、…、N。

(1)

據(jù)經(jīng)驗(yàn)，一般賦值m=2，uij需滿足：

∑iuij=1，?i，i=1、2、3、…、C。

(2)

0<∑iuij

(3)

各變量偏導(dǎo)可通過(guò)以上約束條件求得式(1)最小時(shí)相應(yīng)的隸屬度變量uij與聚類中心vi。

uij=1/∑k[d2(xj,vi)/d2(xj,vk)]2/(m-1)，

k=1、2、3、…、C。

(4)

j=1、2、3、…、N。

(5)

1.3 飛蛾撲火算法

2015年，Seyedali首次提出飛蛾撲火(MFO)概念[18]，該算法尋優(yōu)能力良好，參數(shù)少。設(shè)飛蛾位置為空間變量，飛蛾不斷改變位置矢量飛行于超維空間。飛蛾和火焰兩者區(qū)別，為算法迭代與更新方式間的差異。飛蛾不斷進(jìn)行空間搜索，將飛蛾當(dāng)前獲最優(yōu)位置看作火焰。在該算法中，將更新位置的主要機(jī)制設(shè)為對(duì)數(shù)螺旋，公式為：

Mi=Di·ebθ·cos(2πθ)+Fj。

(6)

式中：θ是為-1～1間的隨機(jī)數(shù)；b為常數(shù)；Mi為飛蛾i；Fj為火焰j；Di為兩者間的間距，Di=|Fj-Mi|。

飛蛾路徑如圖1所示。θ為收斂常數(shù)，定義飛蛾的下一次位置靠近火焰程度，θ=-1時(shí)飛蛾離火焰最近、θ=1時(shí)飛蛾離火焰距離最遠(yuǎn)。

圖1 參數(shù)θ的分布

式(6)僅定義飛蛾撲火的過(guò)程，易收斂過(guò)快，致局部最優(yōu)。為了避免此種情況，需計(jì)算飛蛾適應(yīng)度值，火焰以適應(yīng)度值為標(biāo)準(zhǔn)排序，飛蛾通過(guò)式(6)更新位置逐次迭代。通過(guò)這種方式，每只飛蛾與一簇火焰對(duì)應(yīng)，初始飛蛾靠向最優(yōu)火焰，末尾飛蛾靠向最差火焰。同時(shí)，對(duì)于火焰采用自適應(yīng)，減少機(jī)制，增大效率，尋找最優(yōu)解。當(dāng)?shù)磷詈笠淮危瑢⒈Ａ粑ㄒ换鹧孀鳛槿肿顑?yōu)解。公式為：

f=r[N-t(N-1)/T]。

(7)

式中：f為當(dāng)前迭代次數(shù)的火焰數(shù)量；r為向下取整操作的函數(shù)；t為迭代次數(shù)；T為最大迭代次數(shù)；N為初始火焰數(shù)。

1.4 飛蛾撲火算法改進(jìn)的FCM算法原理

傳統(tǒng)FCM算法是通過(guò)類梯度下降方法求取目標(biāo)函數(shù)J(U,V)的最優(yōu)解。由于隸屬度變量U為模糊值且不斷變化，通過(guò)類梯度下降法求目標(biāo)函數(shù)的最優(yōu)解，函數(shù)收斂較慢，且易于陷入局部極小值化。因此，結(jié)合圖像像素點(diǎn)位置坐標(biāo)，并引入飛蛾撲火算法，對(duì)其進(jìn)行改進(jìn)。首先通過(guò)腐蝕膨脹開(kāi)閉運(yùn)算將位置固定，并求得各樣本聚類中心，而后在進(jìn)行模糊均值聚類時(shí)，通過(guò)式(8)求得聚類距離(D)。

D=R2+G2+B2+X2+Y2。

(8)

式中：R、G、B為RGB色彩圖像中顏色通道色彩分量；X、Y為當(dāng)前像素點(diǎn)到聚類中心點(diǎn)間距離，避免了目標(biāo)函數(shù)陷入局部極小值化。

引入飛蛾撲火算法，隨機(jī)選取飛蛾初始位置并計(jì)算其適應(yīng)度值，通過(guò)式(6)迭代，直至達(dá)到最大迭代次數(shù)，保留最佳位置作為新的聚類中心對(duì)圖像數(shù)據(jù)元素進(jìn)行重新聚類劃分。

2 結(jié)果與分析

2.1 改進(jìn)FCM算法的圖像分割

樣本圖像膨脹、腐蝕處理及中心點(diǎn)繪制過(guò)程見(jiàn)圖2。預(yù)處理過(guò)程中，由于部分樣本發(fā)霉、蟲(chóng)蛀等缺陷，導(dǎo)致二值化圖片有白色空洞，為了將其填補(bǔ)從而獲得更加完整準(zhǔn)確的圖像，對(duì)樣本圖像進(jìn)行膨脹處理。同時(shí)，針對(duì)試驗(yàn)過(guò)程中，樣本個(gè)體間存在黏連現(xiàn)象，對(duì)樣本圖像進(jìn)行腐蝕處理。通過(guò)不斷的膨脹腐蝕開(kāi)閉運(yùn)算，直至將樣本個(gè)體彼此間清晰分開(kāi)，繼而通過(guò)搜索連通域，選取連通域中心點(diǎn)。

圖2 樣本預(yù)處理圖像

用FCM改進(jìn)前后的兩種方法對(duì)整體樣本圖像進(jìn)行分割后得到樣本個(gè)體圖像，分別從中隨機(jī)選取飽滿、干癟、霉斑、蟲(chóng)蛀、腐爛5種圖像(見(jiàn)圖3、圖4)。對(duì)比可見(jiàn)，傳統(tǒng)方法所得圖像，細(xì)節(jié)處分割效果不佳，且存在像素失真現(xiàn)象；對(duì)于存在缺陷的干癟、霉斑、蟲(chóng)蛀等籽粒，分割邊界不夠清晰。改進(jìn)后的FCM算法分割所得圖像，邊緣更加清晰流暢，更有利于試驗(yàn)的后期處理。

圖3 傳統(tǒng)FCM算法處理圖像

圖4 飛蛾撲火優(yōu)化的FCM算法處理圖像

2.2 樣本圖像的邊緣提取

在圖像處理過(guò)程中，通過(guò)邊緣拐點(diǎn)分布情況，能更準(zhǔn)確地了解圖像的形狀特征，解決實(shí)際問(wèn)題[19]。試驗(yàn)中，飽滿樣本的邊緣較為平滑，而干癟、霉斑、蟲(chóng)蛀、腐爛樣本，由于表面缺陷導(dǎo)致其邊緣曲線平滑度較差。因此，對(duì)個(gè)體樣本進(jìn)行邊緣提取，并以通過(guò)拐點(diǎn)個(gè)數(shù)的多少判斷樣本表面光整度，從而對(duì)其達(dá)到缺陷識(shí)別效果。與其他算子(如索貝爾(Sobel)算子)相比，康尼(Canny)算子[20-21]因其抗噪能力強(qiáng)，能檢測(cè)到易被忽略的弱邊緣。因此，選取Canny算子提取樣本邊緣，以便更好識(shí)別出樣本邊緣輪廓間的細(xì)微差異。圖5為通過(guò)Canny算子對(duì)5種樣本圖像進(jìn)行邊緣檢測(cè)所得圖像，圖6為樣本拐點(diǎn)標(biāo)記圖像。

圖5 采用Canny算子邊緣提取樣本圖像

圖6 拐點(diǎn)標(biāo)記樣本圖像

圖6中檢測(cè)到的拐點(diǎn)個(gè)數(shù)，飽滿粒為6、干癟粒為8、霉斑粒為32、蟲(chóng)蛀粒為24、腐爛粒為20。可見(jiàn)，飽滿粒、干癟粒，由于表面相對(duì)光滑，通過(guò)對(duì)其邊緣檢測(cè)，得到的拐點(diǎn)數(shù)在10個(gè)以下；而霉斑、蟲(chóng)蛀、腐爛粒，由于變質(zhì)、表面損壞等原因，拐點(diǎn)個(gè)數(shù)均大于10；因此，擬將邊緣圖像拐點(diǎn)數(shù)作為榛子仁缺陷檢測(cè)的區(qū)分指標(biāo)之一。為了使試驗(yàn)結(jié)果更具有普遍性、代表性，現(xiàn)從已采集的5種樣本圖像中選取清晰圖像進(jìn)行邊緣檢測(cè)處理，并求取拐點(diǎn)個(gè)數(shù)(見(jiàn)表2)。由表2可見(jiàn)：96.94%的飽滿粒檢測(cè)到的拐點(diǎn)數(shù)在10個(gè)以下，95.88%的干癟粒檢測(cè)到的拐點(diǎn)數(shù)在10個(gè)以下，其余3種缺陷(包括霉斑、蟲(chóng)蛀、腐爛粒)邊緣拐點(diǎn)個(gè)數(shù)10個(gè)以上所占比例均達(dá)到了95%以上。由此得出，在誤差允許范圍內(nèi)(小于0.05)，樣本圖像邊緣輪廓檢測(cè)到的拐點(diǎn)數(shù)大于10個(gè)，即可標(biāo)記為缺陷。

表2 數(shù)據(jù)集樣本邊緣輪廓檢測(cè)到的拐點(diǎn)個(gè)數(shù)及所占比例

2.3 樣本圖像邊緣的橢圓曲線擬合

飽滿樣本與其他幾種樣本相比，邊緣更加圓潤(rùn)，邊緣曲線接近橢圓，因此將提取到的邊緣圖像進(jìn)行橢圓曲線擬合。對(duì)于可以識(shí)別檢測(cè)到的橢圓圖像標(biāo)記為飽滿，否則標(biāo)記為缺陷?，F(xiàn)有的橢圓檢測(cè)方法中，最小二乘法[22]、依據(jù)Hough變換法[23-24]最為常用。但最小二乘法對(duì)于孤立點(diǎn)的判斷，易受到噪聲的影響，擬合誤差很大。因此在橢圓輪廓提取中，Hough變換法效果更好。

將圖6中樣本邊緣圖像通過(guò)Hough變換進(jìn)行橢圓曲線擬合，其中只有飽滿籽粒邊緣得到橢圓擬合曲線(見(jiàn)圖7)，其余籽粒均未檢測(cè)到橢圓擬合曲線，標(biāo)記為缺陷粒。

圖7 Hough變換邊緣橢圓擬合曲線圖像

由表3可見(jiàn)：飽滿樣本中有96.94%可以得到橢圓擬合曲線，干癟樣本中僅有3.09%可以得到橢圓擬合曲線，而霉斑、蟲(chóng)蛀、腐爛樣本均未輸出橢圓擬合曲線，未輸出擬合曲線樣本均被標(biāo)記為缺陷，在誤差允許范圍內(nèi)(小于0.05)，可以很好地實(shí)現(xiàn)對(duì)榛子仁缺陷準(zhǔn)確識(shí)別的目的。

表3 數(shù)據(jù)集樣本橢圓曲線擬合結(jié)果

3 結(jié)論

對(duì)榛子仁缺陷檢測(cè)提出了一種新的方法。通過(guò)結(jié)合像素點(diǎn)空間位置，并引入飛蛾撲火算法對(duì)FCM算法進(jìn)行改進(jìn)，從整體樣本圖像中獲得清晰且分割效果良好的樣本個(gè)體圖像。通過(guò)對(duì)樣本圖像形態(tài)學(xué)分析，利用拐點(diǎn)數(shù)的標(biāo)識(shí)計(jì)算，判斷榛子仁表面外輪廓光滑度及果仁飽滿程度，識(shí)別出樣本是否存在缺陷并分類標(biāo)記，從而實(shí)現(xiàn)榛子仁的缺陷檢測(cè)。結(jié)果表明，本研究所提的改進(jìn)方法對(duì)于榛子仁樣本缺陷識(shí)別準(zhǔn)確率達(dá)到96%以上。在未來(lái)的研究中，可推廣用于其他堅(jiān)果類的分類識(shí)別中，為林下產(chǎn)品加工行業(yè)的檢測(cè)與識(shí)別提供有力的技術(shù)支持。