基于圖像處理的蘋果大小分級(jí)研究

　　摘 要：為了實(shí)現(xiàn)蘋果的智能分揀，在GB/T 10651-2008《鮮蘋果》的標(biāo)準(zhǔn)下，采用圖像處理方法對(duì)分揀系統(tǒng)中的蘋果大小進(jìn)行分級(jí)。首先，針對(duì)同一種類的蘋果，對(duì)攝取到的蘋果表面圖像進(jìn)行圖像預(yù)處理，得到平滑且無噪聲影響的蘋果二值圖;然后，通過尋找每個(gè)蘋果5幅表面圖像中最小外接矩形長(zhǎng)寬比最接近1的蘋果圖像，依此作為獲取蘋果最大橫截面直徑的重要依據(jù);最后，通過尋找蘋果圖像質(zhì)心，將質(zhì)心與目標(biāo)像素點(diǎn)的最大距離作為半徑，用最小外接圓的方法得到蘋果的最大橫截面直徑。實(shí)驗(yàn)研究表明：選取的30個(gè)紅富士蘋果，最終得到系統(tǒng)分級(jí)正確率為93.3%，系統(tǒng)的平均分級(jí)時(shí)間在1.5 s以內(nèi)，滿足實(shí)際蘋果分級(jí)的需要。所提出的方法可為按照行業(yè)標(biāo)準(zhǔn)對(duì)蘋果進(jìn)行大小分級(jí)提供參考。

　　關(guān)鍵詞：圖像處理;最小外接矩形;最小外接圓;蘋果;分級(jí)

　　中國蘋果產(chǎn)銷量居世界首位[1]，傳統(tǒng)的蘋果分級(jí)由人工完成，長(zhǎng)期以來，其劣勢(shì)逐漸顯現(xiàn)[2]。因此，智能分揀是當(dāng)今蘋果分級(jí)的主要任務(wù)，而大小又是蘋果分級(jí)的重要參考指標(biāo)。國內(nèi)外專家做了許多基于圖像處理的蘋果大小分級(jí)研究，其中PAULUS等[3]提出以表面積、直徑和體積作為蘋果大小分級(jí)的指標(biāo)。應(yīng)義斌等[4]通過尋找柑橘的最小外接矩形的最大橫徑數(shù)來求取柑橘的大小，隨之也被應(yīng)用到了蘋果大小的檢測(cè)中。馮斌等[5]通過果梗側(cè)外形確定軸向，取垂直軸向的最大寬度作為水果大小。

　　陳艷軍等[6]通過比較3幅蘋果圖像的最大直徑來確定蘋果實(shí)際大小，還通過多元線性回歸原理找到3個(gè)最大直徑與蘋果實(shí)際大小的擬合關(guān)系。除此之外，蘋果大小檢測(cè)還有當(dāng)量直徑法、投影面積法和平均半徑法等[7]。而根據(jù)鮮蘋果分級(jí)標(biāo)準(zhǔn)[8]，需選取蘋果最大橫截面直徑作為測(cè)量依據(jù)。因此，本文選取垂直拍攝果徑面的圖像作為獲取最大橫截面直徑的參考，再利用最小外接圓得到蘋果的直徑，最后依據(jù)分級(jí)標(biāo)準(zhǔn)判斷大小等級(jí)。

　　1 蘋果分揀系統(tǒng)

　　蘋果分揀系統(tǒng)主要由水果輸送裝置、圖像采集系統(tǒng)、圖像處理系統(tǒng)和分揀系統(tǒng)組成[9]。其中水果輸送裝置選用市場(chǎng)上廣泛應(yīng)用的滾子式輸送結(jié)構(gòu)，如圖1所示。其主要工作原理是：整個(gè)裝置由一條循環(huán)輸送鏈帶動(dòng)，摩擦表面的滾子由于摩擦作用繞滾子軸轉(zhuǎn)動(dòng)并前進(jìn)，而相鄰的2個(gè)滾子之間有一“凹型槽”，蘋果正好處于“凹型槽”上，隨著滾子翻轉(zhuǎn)滾動(dòng)并向前輸送[10]，蘋果的所有表面都會(huì)被相機(jī)獲取，這也保證了蘋果采集的完整性。水果輸送過程中的實(shí)際狀態(tài)如圖2所示。將相機(jī)采集到的圖像傳送到處理系統(tǒng)進(jìn)行分析，獲得蘋果等級(jí)結(jié)果，然后控制模塊將指令傳遞到分級(jí)執(zhí)行機(jī)構(gòu)，使蘋果在對(duì)應(yīng)分類級(jí)別的位置落下，從而實(shí)現(xiàn)分級(jí)[11]。

　　2 蘋果圖像預(yù)處理

　　2.1 圖像采集

　　根據(jù)蘋果在輸送中的實(shí)際運(yùn)動(dòng)狀態(tài)以及相機(jī)攝取到的多姿態(tài)蘋果圖像，在實(shí)際操作中發(fā)現(xiàn)，5幅圖像就能包含蘋果的所有表面，因此將獲取到的每個(gè)蘋果的5幅圖像作為研究對(duì)象，如圖3所示。

　　2.2 圖像二值化

　　將采集到的蘋果彩色圖像進(jìn)行二值化處理，得到像素值為0的黑色背景和像素值為1的目標(biāo)圖像[12]。具體方法為先將彩色圖像進(jìn)行灰度化處理，利用如下語句實(shí)現(xiàn)：

　　其次，利用Qtsu圖像分割將灰度圖轉(zhuǎn)換成二值圖，通過計(jì)算方差尋找最合適的灰度級(jí)作為閾值，使劃分區(qū)域的兩部分差別最大，Qtsu法的語句實(shí)現(xiàn)如下：

　　2.3 形態(tài)學(xué)去噪

　　圖像二值化后，目標(biāo)區(qū)域可能會(huì)有小的孔洞出現(xiàn)，因此用圖像形態(tài)學(xué)進(jìn)行去噪處理，填補(bǔ)小孔洞，細(xì)化邊緣[13]。使用圓形結(jié)構(gòu)進(jìn)行形態(tài)學(xué)閉操作處理，語句實(shí)現(xiàn)如下：

　　對(duì)于采集到果梗面或有缺陷的圖像，上述操作后還會(huì)有較大孔洞出現(xiàn)，因此使用imfill函數(shù)填充較大孔洞，語句實(shí)現(xiàn)如下。

　　圖像預(yù)處理結(jié)果如圖4所示。

　　3 最小外接圓分級(jí)

　　3.1 果徑面的確定

　　鮮蘋果分級(jí)標(biāo)準(zhǔn)規(guī)定，取蘋果最大橫截面直徑作為蘋果大小等級(jí)的判斷標(biāo)準(zhǔn)。通過觀察可以發(fā)現(xiàn)，由于蘋果是立體圖形，蘋果以正立形式(果梗上置)放置時(shí)，獲取到的蘋果主視圖圖像的最大直徑不一定是蘋果的最大橫截面直徑，而此時(shí)的俯視圖，也就是垂直于果徑拍攝的圖像，它的最大直徑可以作為大小等級(jí)劃分的重要依據(jù)。進(jìn)一步觀察發(fā)現(xiàn)，對(duì)于大多數(shù)蘋果，尤其是形狀不規(guī)則的蘋果而言，垂直于果徑拍攝的圖像比其他姿態(tài)時(shí)的圖像更加接近于圓形[14]。因此，可以利用圖像最小外接矩形長(zhǎng)寬比這一指標(biāo)來尋找最接近圓形的圖像[15]。而且長(zhǎng)寬比數(shù)值越接近于1，圖形越接近于圓，所以可用式(1)來求取圖像最小外接矩形長(zhǎng)寬比與1的差值，差值最小的圖像即為所確定的果徑面的圖像。

　　3.2 最小外接圓圓心確定

　　確定好果徑面后，就要對(duì)果徑面求取它的最大直徑，將其作為蘋果的實(shí)際大小，本文用最小外接圓法求取直徑。首先就要確定最小外接圓的圓心，將蘋果圖像的質(zhì)心作為最小外接圓的圓心，計(jì)算蘋果二值圖目標(biāo)區(qū)域的總像素?cái)?shù)，記作s，質(zhì)心坐標(biāo)可以通過求取二值圖目標(biāo)像素的平均像素點(diǎn)獲得[16]，如式(2)和式(3)所示。

　　3.3 最小外接圓半徑確定

　　確定好最小外接圓圓心后，通過尋找圖像坐標(biāo)與圓心坐標(biāo)的距離，距離值最大的即為最小外接圓半徑。首先，設(shè)初始半徑為2個(gè)像素點(diǎn)，再遍歷各個(gè)像素點(diǎn)，利用兩點(diǎn)間距離公式，如式(4)所示，計(jì)算此時(shí)像素點(diǎn)與圓心的距離，如果大于當(dāng)前半徑，則此距離作為新的半徑，最終迭代更新得到最小外接圓半徑。

　　最終在二值圖的基礎(chǔ)上，根據(jù)圓心和半徑畫出圖像的最小外接圓，如圖5所示。

　　3.4 圖像標(biāo)定

　　系統(tǒng)通過最小外接圓法得到蘋果直徑的單位為像素，而蘋果實(shí)際直徑單位是毫米，因此需要進(jìn)行單位轉(zhuǎn)換。具體方法：選取一標(biāo)準(zhǔn)形狀的球體，利用圖像采集系統(tǒng)進(jìn)行球體的圖像采集，然后運(yùn)用計(jì)算機(jī)繪圖軟件，如Photoshop進(jìn)行球體直徑的測(cè)量，得到的球體直徑為d像素，而球體的實(shí)際直徑為D mm，則系統(tǒng)的標(biāo)定系數(shù)k如式(5)所示：

　　用標(biāo)定系數(shù)k與系統(tǒng)測(cè)量的蘋果直徑進(jìn)行轉(zhuǎn)換，就可以得到實(shí)際的蘋果直徑。

　　4 蘋果分級(jí)結(jié)果分析

　　4.1 單一蘋果果徑面的選取與測(cè)量

　　將一個(gè)蘋果樣本在光箱中進(jìn)行不同姿態(tài)拍攝，獲取5幅圖像(保證提取到所有表面圖像)，經(jīng)過圖像處理與分析后，選擇5幅圖像中最小外接矩形長(zhǎng)寬比最接近于1的圖像，也就是垂直于果徑拍攝的圖像作為果徑面，如表1所示。由表1可知5幅圖像中最小外接矩形長(zhǎng)寬比與1差值最小的是0.015，因此選擇圖像5作為果徑面。然后用最小外接圓方法測(cè)量圖像5的直徑大小，通過圖像標(biāo)定，最終系統(tǒng)測(cè)量蘋果樣本的大小為79.05 mm，而用游標(biāo)卡尺測(cè)量出實(shí)際大小為77.82 mm，兩尺寸相差不大，絕對(duì)誤差為1.23 mm。

　　4.2 蘋果大小分級(jí)結(jié)果

　　采用同樣的方法，隨機(jī)選取30個(gè)同一品種的紅富士蘋果作為研究對(duì)象，依次測(cè)量蘋果的直徑，參考《鮮蘋果》分級(jí)標(biāo)準(zhǔn)和實(shí)際選取的蘋果果徑大小，將各級(jí)之間的分割閾值確定為85，80，75 3個(gè)閾值，其中蘋果直徑大于85 mm的為一級(jí)果，80～85 mm的為二級(jí)果，75～<80 mm的為三級(jí)果，小于75 mm的為四級(jí)果。實(shí)驗(yàn)中，首先用游標(biāo)卡尺分別測(cè)量30個(gè)蘋果的實(shí)際果徑大小，并編號(hào)標(biāo)記;其次，按照編號(hào)依次獲取每個(gè)蘋果的5幅圖像;最后，系統(tǒng)確定果徑面并測(cè)量直徑大小，與蘋果實(shí)際直徑進(jìn)行對(duì)比，分出等級(jí)，如表2所示。

　　實(shí)驗(yàn)結(jié)果顯示，在對(duì)30個(gè)蘋果進(jìn)行大小分級(jí)后，有28個(gè)蘋果分級(jí)正確，2個(gè)蘋果分級(jí)錯(cuò)誤。因此，使用本方法最終得到的蘋果大小分級(jí)準(zhǔn)確率為93.3%。

　　5 結(jié) 語

　　本文對(duì)蘋果大小分級(jí)問題進(jìn)行了研究，首先對(duì)獲取到的蘋果表面圖像進(jìn)行預(yù)處理，得到二值化圖像，然后將垂直果徑的拍攝面圖像作為尋找蘋果最大橫截面直徑的重要依據(jù)，其中通過對(duì)相機(jī)獲取到的每個(gè)蘋果的5幅表面圖像進(jìn)行最小外接矩形長(zhǎng)寬比的比較，將比值最接近于1的圖像作為最終的果徑面。得到果徑面后，采用最小外接圓法得到蘋果的直徑，并將單位轉(zhuǎn)換成毫米，實(shí)驗(yàn)選取30個(gè)紅富士蘋果進(jìn)行實(shí)驗(yàn)驗(yàn)證，最終分級(jí)準(zhǔn)確率為93.3%，系統(tǒng)平均分級(jí)時(shí)間在1.5 s以內(nèi)，滿足實(shí)際蘋果分級(jí)的需要。但在蘋果每個(gè)表面都近似圓形的情況下，就會(huì)影響果徑面的選擇，進(jìn)而導(dǎo)致蘋果大小分級(jí)錯(cuò)誤。針對(duì)這種情況，如何選取果徑面是下一步研究的重點(diǎn)。

　　推薦閱讀：計(jì)算機(jī)圖形圖像處理技術(shù)的作用探究