一種基于圖像處理的胡蘿卜表面缺陷檢測(cè)方法
【專(zhuān)利摘要】本發(fā)明公開(kāi)了一種基于圖像處理的胡蘿卜表面缺陷檢測(cè)方法，包括圖像預(yù)處理和影響胡蘿卜外觀品質(zhì)的須根、彎曲、開(kāi)裂和斷裂4種缺陷的檢測(cè)方法。首先對(duì)胡蘿卜圖像進(jìn)行預(yù)處理，然后分別對(duì)須根、彎曲、開(kāi)裂和斷裂進(jìn)行檢測(cè)。須根檢測(cè)通過(guò)采用凹點(diǎn)檢測(cè)法對(duì)胡蘿卜輪廓上的點(diǎn)進(jìn)行逐點(diǎn)判斷；彎曲檢測(cè)通過(guò)獲取胡蘿卜區(qū)域的外接凸多邊形，并計(jì)算外接凸多邊形與胡蘿卜區(qū)域面積之比進(jìn)行；開(kāi)裂檢測(cè)通過(guò)在R分量圖像上進(jìn)行Canny邊緣檢測(cè)和Hough變換進(jìn)行。在對(duì)胡蘿卜朝向判斷的基礎(chǔ)上分別對(duì)上下部斷裂進(jìn)行檢測(cè)，上部斷裂檢測(cè)通過(guò)計(jì)算所截取的上下橫徑比得到，下部斷裂通過(guò)計(jì)算斷面的彎曲程度進(jìn)行。通過(guò)采用圖像處理技術(shù)對(duì)胡蘿卜表面缺陷進(jìn)行自動(dòng)檢測(cè)，可大大提高胡蘿卜外觀品質(zhì)檢測(cè)效率和準(zhǔn)確率，節(jié)省大量人力成本。
【專(zhuān)利說(shuō)明】
-種基于圖像處理的胡蘿卜表面缺陷檢測(cè)方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及一種在農(nóng)業(yè)生產(chǎn)加工領(lǐng)域內(nèi)使用的農(nóng)產(chǎn)品外觀品質(zhì)的檢測(cè)方法，具體 的說(shuō)涉及一種胡蘿卜須根、彎曲、開(kāi)裂和斷裂等表面缺陷的量化檢測(cè)方法。
【背景技術(shù)】
[0002] 胡蘿卜的分級(jí)銷(xiāo)售有助于提高胡蘿卜的市場(chǎng)競(jìng)爭(zhēng)力以及增加經(jīng)濟(jì)效益，目前一些 胡蘿卜生產(chǎn)加工企業(yè)主要依靠手工方式進(jìn)行檢測(cè)和分級(jí)，采用手工方式進(jìn)行檢測(cè)和分級(jí)的 方法增加了勞動(dòng)力成本，不能保證分選結(jié)果的準(zhǔn)確率與一致性，且生產(chǎn)效率低，隨著勞動(dòng)力 成本的提高，傳統(tǒng)的手工方式無(wú)法滿足當(dāng)前胡蘿卜加工企業(yè)大規(guī)模生產(chǎn)的需求。
[0003] 計(jì)算機(jī)視覺(jué)提供了一種高效、低成本、高準(zhǔn)確率的果蔬自動(dòng)檢測(cè)和分級(jí)方法，采用 圖像處理技術(shù)檢測(cè)胡蘿卜外觀品質(zhì)，不僅可以減少勞動(dòng)力成本，而且能夠獲得更高精度的 結(jié)果，使得分級(jí)具有更高的一致性。雖然目前已經(jīng)有大量基于圖像處理的果蔬外觀品質(zhì)檢 測(cè)方法，但主要集中在球型和類(lèi)球型果蔬上，如蘋(píng)果、柑橘、馬鈴薯、梨等。由于不同果蔬外 觀表現(xiàn)不同，影響其外觀品質(zhì)的性狀也會(huì)不同，同時(shí)由于研究方法固有的專(zhuān)用性和局限性， 在胡蘿卜等長(zhǎng)型果蔬上的應(yīng)用受到限制，因此需要針對(duì)胡蘿卜外觀品質(zhì)特征，研究適合于 其外觀性狀的檢測(cè)方法。
[0004] 目前國(guó)內(nèi)還沒(méi)有相對(duì)成熟的基于圖像處理的胡蘿卜分級(jí)檢測(cè)方法，相關(guān)研究是一 個(gè)具有挑戰(zhàn)性的課題。根據(jù)中華人民共和國(guó)國(guó)內(nèi)貿(mào)易標(biāo)準(zhǔn)（SB/T10450-2007)胡蘿卜銷(xiāo)售 等級(jí)要求（中華人民共和國(guó)商務(wù)部2007. 12. 28發(fā)布)標(biāo)準(zhǔn)，影響胡蘿卜銷(xiāo)售等級(jí)的關(guān)鍵指標(biāo) 有開(kāi)裂、彎曲、須根等，但其衡量準(zhǔn)則比較模糊，為胡蘿卜在線檢測(cè)帶來(lái)了較大的困難，對(duì)胡 蘿卜進(jìn)行準(zhǔn)確檢測(cè)需要對(duì)這些關(guān)鍵指標(biāo)進(jìn)行量化。因此采用圖像處理方法對(duì)影響胡蘿卜外 觀品質(zhì)的關(guān)鍵指標(biāo)進(jìn)行量化檢測(cè)是非常有意義的。 基于此，本案提出一種基于圖像處理的胡蘿卜表面缺陷檢測(cè)方法。

【發(fā)明內(nèi)容】

[0005] 本發(fā)明所要解決的技術(shù)問(wèn)題就是提供一種基于計(jì)算機(jī)圖像處理技術(shù)對(duì)胡蘿卜須 根、彎曲、開(kāi)裂和斷裂進(jìn)行量化檢測(cè)的方法。
[0006] 本發(fā)明采用如下技術(shù)方案： 胡蘿卜圖像預(yù)處理方法，該方法步驟如下： 從胡蘿卜源圖像中提取R、G、B三個(gè)分量圖像，采用公式（1)對(duì)R、G、B三個(gè)分量圖像進(jìn) 行加權(quán)求和得到灰度圖像Gray (其中R、G、B分別為RGB圖像的三個(gè)分量圖像)，將復(fù)雜背 景轉(zhuǎn)換為黑色，然后采用全局閾值對(duì)灰度圖像進(jìn)行圖像分割，得到二值圖像BW。 = -3G - 0.SB (1〕
[0007] 基于圖像處理的胡蘿卜須根檢測(cè)方法包括如下步驟： (1) 在二值圖像BW上進(jìn)行輪廓提取，獲取輪廓上點(diǎn)的集合Pl ; (2) 對(duì)Pl中各點(diǎn)逐點(diǎn)判斷是否為凹點(diǎn)，得到可能凹點(diǎn)的集合P2 ;判斷凹點(diǎn)采用的方法 是：對(duì)當(dāng)前檢測(cè)點(diǎn)，取與其等距的前后兩點(diǎn)及其中點(diǎn)，然后判斷中點(diǎn)是否在胡蘿卜內(nèi)部，若 不在胡蘿卜區(qū)域內(nèi)，則此點(diǎn)可能為凹點(diǎn)； (3) 對(duì)P2中的各點(diǎn)分別計(jì)算凹陷程度：分別連接當(dāng)前點(diǎn)與其前驅(qū)點(diǎn)和后繼點(diǎn)，計(jì)算兩 條線段之間的夾角Θ，若Θ小于設(shè)定的閾值，則認(rèn)為當(dāng)前點(diǎn)為須根所在位置，夾角Θ的計(jì) 算方法如公式（2)所示：
夾角Θ越小，說(shuō)明此點(diǎn)凹陷程度越大； (4) 當(dāng)夾角Θ小于一定的閾值，則該點(diǎn)為須根可能所在位置，得到的點(diǎn)的集合記為P3; (5) 對(duì)P3中的鄰接點(diǎn)進(jìn)行合并，對(duì)鄰接點(diǎn)保留凹陷程度最大的點(diǎn)，所檢測(cè)到的點(diǎn)的個(gè) 數(shù)即為須根個(gè)數(shù)。
[0008] 基于圖像處理的胡蘿卜彎曲檢測(cè)方法，該方法包括以下步驟： (1) 在二值化圖像BW上提取胡蘿卜輪廓； (2) 根據(jù)胡蘿卜輪廓得到其擬合外接凸多邊形，計(jì)算外接凸多邊形面積計(jì)為Area2 ; (3) 在二值圖像BW上，計(jì)算胡蘿卜區(qū)域面積Areal ; (4) 計(jì)算形狀規(guī)則度CurveRate=Areal/Area2 e [0. 0, 1. 0];此值越大，形狀越規(guī)則； (5) 當(dāng)形狀規(guī)則度大于一定的閾值，判斷該胡蘿卜嚴(yán)重彎曲，為次品； 本發(fā)明中，形狀規(guī)則度CurveRate閾值設(shè)定為0. 9,當(dāng)CurveRate小于0. 9,判斷該胡蘿 卜為彎曲。
[0009] 基于圖像處理的胡蘿卜開(kāi)裂檢測(cè)方法，該方法包括以下步驟： (1) 獲取胡蘿卜圖像的R分量圖像，使用canny算子對(duì)R分量圖像進(jìn)行邊緣檢測(cè)，得到 邊緣圖像C ; (2) 對(duì)二值圖像BW進(jìn)行腐蝕操作去除邊界部分； (3) 將邊緣圖像C與二值圖像BW進(jìn)行與運(yùn)算，只保留內(nèi)部曲線； (4) 采用霍夫變換檢測(cè)曲線，將長(zhǎng)度大于一定閾值的曲線判斷為開(kāi)裂； 基于圖像處理的胡蘿卜朝向檢測(cè)方法，該方法步驟如下： (1) 在二值圖像BW上，分別從距離胡蘿卜兩端長(zhǎng)度的1/10處，取兩段橫徑的長(zhǎng)度，分別 記為Wl和W2 ; (2) 計(jì)算兩個(gè)橫徑的比值wRate=Wl/W2 ; (3) 根據(jù)上下橫徑比wRate來(lái)判斷胡蘿卜的朝向，若wRate > 1且胡蘿卜垂直放置，則 胡蘿卜朝上，否則朝下；若《Rate彡1且胡蘿卜水平放置，則胡蘿卜朝左，否則朝右。
[0010] 基于圖像處理的胡蘿卜上部斷裂檢測(cè)方法，該方法包括如下步驟： (1) 根據(jù)胡蘿卜朝向檢測(cè)中獲得的Wl和W2值，取較小值記為WMin，較大值記為WMax ; (2) 計(jì)算較小值與較大值之比：mRate=wMin/wMax e [0, 1]; (3) 根據(jù)mRate的值來(lái)判斷胡蘿卜是否上部斷裂，此值若大于給定的閾值，則判斷為上 部斷裂； 本發(fā)明中，mRate閾值設(shè)定為0. 8,若mRate的值大于0. 8,則判斷該胡蘿卜為上部斷裂。
[0011] 基于圖像處理的胡蘿卜下部斷裂檢測(cè)方法，該方法包括以下步驟： (1)在二值圖像BW上提取胡蘿卜輪廓，提取輪廓上的點(diǎn)并保存到數(shù)組Pl中； (2) 根據(jù)胡蘿卜朝向判斷結(jié)果，在胡蘿卜上部1/5區(qū)域內(nèi)使用局部極值法判斷胡蘿卜 斷面區(qū)域的兩個(gè)端點(diǎn)，具體做法是在區(qū)域內(nèi)計(jì)算位于胡蘿卜兩側(cè)的邊界點(diǎn)，記為A和B ; (3) 在數(shù)組Pl內(nèi)獲取A、B之間的所有的點(diǎn)，即斷面區(qū)域內(nèi)所有的點(diǎn)，并保存到數(shù)組P2 中； (4) 計(jì)算斷面的彎曲程度：依次計(jì)算數(shù)組P2中的各點(diǎn)到線段AB的距離，并計(jì)算其均值 得到斷面的彎曲程度cRate (式3);
其中：N表示斷面上點(diǎn)的個(gè)數(shù)，di農(nóng)不弟i個(gè)點(diǎn)到線段AB的距離； (5) 根據(jù)cRate的值來(lái)判斷胡蘿卜是否為下部斷裂，此值若小于給定的閾值，則判斷胡 蘿卜為下部斷裂。
[0012] 本發(fā)明的有益效果是： 本發(fā)明所公開(kāi)的一種胡蘿卜須根、彎曲、開(kāi)裂和斷裂的檢測(cè)方法，通過(guò)攝像頭實(shí)時(shí)獲取 胡蘿卜圖像，基于圖像處理技術(shù)對(duì)獲取的圖像進(jìn)行缺陷測(cè)量和檢測(cè)，從而克服人工檢測(cè)帶 來(lái)的主觀性，使量化檢測(cè)更加客觀、科學(xué)。應(yīng)用于農(nóng)業(yè)生產(chǎn)、品質(zhì)分級(jí)與商貿(mào)流通，可大大提 高生產(chǎn)效率和降低勞動(dòng)力成本。
【附圖說(shuō)明】
[0013] 圖1是本發(fā)明的缺陷檢測(cè)整體流程圖。
[0014] 圖2是本發(fā)明的須根檢測(cè)流程圖。
[0015] 圖3是本發(fā)明的凹點(diǎn)檢測(cè)原理示意圖。
[0016] 圖4是本發(fā)明的上部斷裂檢測(cè)示意圖。
[0017] 圖5是本發(fā)明的下部斷裂檢測(cè)原理示意圖。
【具體實(shí)施方式】
[0018] 為了使本發(fā)明的目的、技術(shù)方案及優(yōu)點(diǎn)更加清楚明白，以下結(jié)合附圖及實(shí)施例，對(duì) 本發(fā)明進(jìn)行進(jìn)一步詳細(xì)說(shuō)明。應(yīng)當(dāng)理解，此處所描述的具體實(shí)施例僅用以解釋本發(fā)明，并不 用于限定本發(fā)明。
[0019] 實(shí)施例，請(qǐng)參照?qǐng)D1-圖5 :圖1示意性地給出了本發(fā)明所公開(kāi)的胡蘿卜表面缺陷 檢測(cè)整體流程；圖2示意性的示出了本發(fā)明所公開(kāi)的一種胡蘿卜須根圖像量化檢測(cè)方法流 程圖；圖3示意性的示出了本發(fā)明所公開(kāi)的胡蘿卜須根檢測(cè)方法中凹點(diǎn)檢測(cè)原理；圖4示 意性地給出了本發(fā)明所公開(kāi)的胡蘿卜上部斷裂檢測(cè)示意圖；圖5示意性地給出了本發(fā)明所 公開(kāi)的胡蘿卜下部斷裂檢測(cè)原理圖。
[0020] 首先采用CCD攝像頭拍攝胡蘿卜圖像，圖像分辨率為320X240且胡蘿卜為垂直放 置。接著依據(jù)圖1所示流程依次對(duì)所采集的圖像進(jìn)行如下處理： 首先對(duì)圖像進(jìn)行預(yù)處理，從所采集的圖像中提取R、G、B三個(gè)分量圖像，采用 2. 5R-2G-0. 5B得到灰度圖像(其中R、G、B分別為RGB圖像的三個(gè)分量圖像)，然后采用全局 閾值(此處為100)對(duì)灰度圖像進(jìn)行二值化，得到二值圖像BW。接下來(lái)根據(jù)前面的技術(shù)方案 分別對(duì)須根、彎曲、開(kāi)裂和斷裂進(jìn)行檢測(cè)。
[0021 ] 依照?qǐng)D2所示流程，對(duì)胡蘿卜須根依據(jù)以下步驟進(jìn)行檢測(cè)： (1) 在二值圖像BW上進(jìn)行輪廓提取，得到輪廓上的點(diǎn)的集合Pl ; (2) 逐點(diǎn)判斷Pl中各點(diǎn)是否為凹點(diǎn)，得到可能凹點(diǎn)的集合P2。設(shè)當(dāng)前檢測(cè)點(diǎn)為A，取其 等距(距離為20像素）的前驅(qū)點(diǎn)和后繼點(diǎn)分別為B、C，如圖3所示。然后判斷點(diǎn)B和C的中 點(diǎn)M是否在胡蘿卜內(nèi)部，若點(diǎn)M不在胡蘿卜區(qū)域內(nèi)，此點(diǎn)為可能凹點(diǎn)； (3) 對(duì)P2中的各點(diǎn)分別計(jì)算凹陷程度：計(jì)算直線AB與AC所成夾角Θ的大小，夾角Θ 越小，說(shuō)明此點(diǎn)凹陷程度越大。若Θ小于設(shè)定的閾值，則認(rèn)為點(diǎn)A為須根所在位置； (4) 由于B和C與A等距的，為減少計(jì)算量，在實(shí)現(xiàn)過(guò)程中可直接通過(guò)BC的長(zhǎng)度a來(lái)判 斷點(diǎn)A的凹陷程度，a的值越小，該點(diǎn)的凹陷程度越大。若a的值小于設(shè)定的閾值T，則認(rèn)為 A點(diǎn)為須根所在位置，得到的點(diǎn)的集合記為P3 ; (5) 對(duì)P3中的鄰接點(diǎn)進(jìn)行合并，所檢測(cè)到的點(diǎn)的個(gè)數(shù)即為須根個(gè)數(shù)； 本案中，步驟（4)中閾值T設(shè)置為38,若T小于38,則A點(diǎn)為須根所在位置。
[0022] 胡蘿卜彎曲的檢測(cè)步驟如下： (1) 在二值化圖像BW上提取胡蘿卜輪廓； (2) 計(jì)算外接凸多邊形面積Area2和胡蘿卜區(qū)域面積Areal ; (3) 計(jì)算形狀規(guī)則度CurveRate=Areal/Area2 e [0. 0, 1. 0];此值越大，形狀越規(guī)則；若 形狀規(guī)則度CurveRate小于0. 9,則判斷該胡蘿卜為彎曲。
[0023] 胡蘿卜開(kāi)裂檢測(cè)方法的步驟如下： (1) 獲取胡蘿卜圖像的R分量圖像，使用canny算子對(duì)R分量圖像進(jìn)行邊緣檢測(cè)，得到 邊緣圖像C ; (2) 對(duì)二值圖像BW進(jìn)行腐蝕操作去除邊界部分； (3) 將邊緣圖像C與二值圖像BW進(jìn)行與運(yùn)算，只保留內(nèi)部曲線； (4) 采用霍夫變換檢測(cè)曲線，將長(zhǎng)度大于一定閾值的曲線判斷為開(kāi)裂；若曲線長(zhǎng)度大于 10,則判定胡蘿卜有開(kāi)裂。
[0024] 在進(jìn)行斷裂檢測(cè)前，首先要判斷胡蘿卜的朝向，朝向判斷的步驟如下： (1) 在二值圖像BW上，分別從距離胡蘿卜兩端長(zhǎng)度的1/10處，取兩段橫徑的長(zhǎng)度，分別 記為Wl和W2,如圖4所示； (2) 計(jì)算兩個(gè)橫徑的比值wRate=Wl/W2 ; (3) 根據(jù)上下橫徑比wRate的值來(lái)判斷胡蘿卜的朝向，若wRate > 1，則胡蘿卜大頭朝 上，否則朝下； 胡蘿卜上部斷裂檢測(cè)步驟如下： (1) 根據(jù)胡蘿卜朝向檢測(cè)中所獲取的Wl和W2值，較小值記為wMin，較大值記為wMax ; (2) 計(jì)算 wMin 與 wMax 的比值：mRate=wMin/wMax e [0, 1]; (3) 根據(jù)mRate的值來(lái)判斷胡蘿卜是否為上部斷裂，若值大于0. 8,則判斷胡蘿卜為上 部斷裂。
[0025] 基于圖像處理的胡蘿卜下部斷裂檢測(cè)方法，該方法包括以下步驟： (1) 在二值圖像BW上提取胡蘿卜輪廓，并提取輪廓上的點(diǎn)保存到數(shù)組Pl中； (2) 根據(jù)胡蘿卜朝向判斷結(jié)果，參照?qǐng)D5,在胡蘿卜上部1/5區(qū)域內(nèi)使用局部極值法計(jì) 算胡蘿卜斷面區(qū)域的兩個(gè)端點(diǎn)，記為A和B ;具體做法是在上部1/5區(qū)域內(nèi)，分別計(jì)算最左 側(cè)的點(diǎn)A和最右側(cè)的點(diǎn)B ; (3) 獲取A、B之間所有的點(diǎn)，并保存到數(shù)組P2中； (4) 計(jì)算斷面的彎曲程度：依次計(jì)算數(shù)組P2中的各點(diǎn)到線段AB的距離，并計(jì)算其均值 得到斷面的彎曲程度cRate (式3); (5) 根據(jù)cRate的值來(lái)判斷胡蘿卜是否為下部斷裂，若其值小于12,則判斷胡蘿卜為下 部斷裂。
【主權(quán)項(xiàng)】
1. 一種基于圖像處理的胡蘿卜表面缺陷檢測(cè)方法，其特征在于：該方法包括圖像預(yù)處 理方法和影響胡蘿卜外觀品質(zhì)的須根、彎曲、開(kāi)裂和斷裂等表面缺陷檢測(cè)方法。2. 根據(jù)權(quán)利要求1所述的圖像預(yù)處理方法，其特征在于：對(duì)胡蘿卜圖像的R、G、B三個(gè) 分量圖像進(jìn)行加權(quán)求和得到灰度圖像，然后對(duì)灰度圖像采用全局閾值進(jìn)行圖像分割得到二 值圖像。3. 根據(jù)權(quán)利要求1所述的胡蘿卜須根檢測(cè)，其特征在于：提取胡蘿卜輪廓上的點(diǎn)序列， 根據(jù)凹點(diǎn)判斷原理逐點(diǎn)判斷輪廓上的點(diǎn)是否為凹點(diǎn)，若為凹點(diǎn)則進(jìn)一步計(jì)算凹陷程度，當(dāng) 某點(diǎn)凹陷程度大于給定的閾值，則判定該點(diǎn)為須根所在位置；然后對(duì)所檢測(cè)到的鄰接點(diǎn)進(jìn) 行合并，得到的點(diǎn)的個(gè)數(shù)即為須根個(gè)數(shù)。4. 根據(jù)權(quán)利要求1所述的胡蘿卜彎曲檢測(cè)方法，其特征在于：根據(jù)胡蘿卜輪廓得到其 擬合外接凸多邊形并計(jì)算其面積，將形狀規(guī)則度定義為胡蘿卜區(qū)域面積與外接凸多邊形面 積比，然后根據(jù)形狀規(guī)則度的值判斷胡蘿卜彎曲程度。5. 根據(jù)權(quán)利要求1所述的胡蘿卜開(kāi)裂檢測(cè)方法，其特征在于：首先使用canny算子對(duì)R 分量圖像進(jìn)行邊緣檢測(cè)，得到邊緣圖像；接著對(duì)二值圖像進(jìn)行腐蝕操作去除邊界部分，然后 將邊緣圖像與二值圖像進(jìn)行"與"運(yùn)算去除邊界，只保留內(nèi)部曲線；最后采用霍夫變換檢測(cè) 曲線，將長(zhǎng)度大于給定閾值的曲線判斷為開(kāi)裂。6. 根據(jù)權(quán)利要求1所述的胡蘿卜斷裂檢測(cè)方法，其特征在于：在對(duì)胡蘿卜朝向判斷的 基礎(chǔ)上分別對(duì)上部和下部斷裂進(jìn)行檢測(cè)。7. 根據(jù)權(quán)利要求6所述的胡蘿卜朝向判斷方法，其特征在于：分別從胡蘿卜兩端取兩 個(gè)橫截面，計(jì)算兩個(gè)橫徑的長(zhǎng)度，通過(guò)比較兩個(gè)橫徑長(zhǎng)度的大小來(lái)判斷胡蘿卜的朝向。8. 根據(jù)權(quán)利要求6所述的胡蘿卜上部斷裂檢測(cè)方法，其特征在于：采用權(quán)利要求7所 述方法從與兩端等距的位置計(jì)算兩個(gè)橫徑的長(zhǎng)度，并計(jì)算較小值與較大值之間的比值，并 根據(jù)設(shè)定的閾值來(lái)判斷是否為上部斷裂。9. 根據(jù)權(quán)利要求6所述的胡蘿卜下部斷裂檢測(cè)方法，其特征在于：從胡蘿卜上部獲取 斷面區(qū)域，然后計(jì)算斷面的彎曲程度，并根據(jù)設(shè)定的閾值來(lái)判斷是否為下部斷裂。
【文檔編號(hào)】G01N21/952GK105891231SQ201510037746
【公開(kāi)日】2016年8月24日
【申請(qǐng)日】2015年1月26日
【發(fā)明人】鄧立苗, 韓仲志, 扈志強(qiáng), 耿琪超
【申請(qǐng)人】青島農(nóng)業(yè)大學(xué), 青島有田農(nóng)業(yè)發(fā)展有限公司, 青島大谷農(nóng)業(yè)信息有限公司