• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于機(jī)器視覺的馬鈴薯黑心病檢測機(jī)構(gòu)設(shè)計(jì)與試驗(yàn)

    2017-06-05 15:00:27彭彥昆魏文松鄭曉春王文秀
    農(nóng)業(yè)工程學(xué)報 2017年5期
    關(guān)鍵詞:心病塊莖光源

    田 芳,彭彥昆,魏文松,鄭曉春,王文秀

    基于機(jī)器視覺的馬鈴薯黑心病檢測機(jī)構(gòu)設(shè)計(jì)與試驗(yàn)

    田 芳,彭彥昆※,魏文松,鄭曉春,王文秀

    (中國農(nóng)業(yè)大學(xué)工學(xué)院,國家農(nóng)產(chǎn)品加工技術(shù)裝備研發(fā)分中心,北京 100083)

    馬鈴薯黑心病是一種難以用肉眼分辨的內(nèi)部病害。在生產(chǎn)加工中,為了準(zhǔn)確檢出黑心病個體,該研究基于光在馬鈴薯塊莖組織中的透射特性設(shè)計(jì)了一種馬鈴薯黑心病檢測機(jī)構(gòu)。機(jī)構(gòu)主要包括馬鈴薯塊莖承托槽、光源單元以及光源位置調(diào)節(jié)軌道等部分,光源單元包括2個LED面光源,形成多點(diǎn)透射的結(jié)構(gòu),光源波長為705 nm。文中取50個正常樣品和29個黑心樣品采集圖像分析,經(jīng)過預(yù)處理,分別提取馬鈴薯目標(biāo)圖像和馬鈴薯高灰度區(qū)域的灰度平均值,以二者的相對比值作為黑心病的判別參數(shù)。通過與單點(diǎn)透射模式所采集的正常樣品圖像的相對比值比較,由該機(jī)構(gòu)所采集的正常圖像的相對比值普遍偏小,表明所獲得圖像的灰度均勻性更好,同時得到光源較佳夾角為45°。在該機(jī)構(gòu)下獲得的正常樣品和黑心樣品的灰度平均值相對比值的分類閾值為0.220,得到2種樣品的判斷正確率分別為98%和96.6%,且黑心樣品的相對比值更大,表明使用該機(jī)構(gòu)可以檢測黑心馬鈴薯樣品。該研究為馬鈴薯黑心病檢測機(jī)構(gòu)的設(shè)計(jì)提供了參考。

    機(jī)器視覺;無損檢測;設(shè)計(jì);馬鈴薯塊莖;黑心病;檢測機(jī)構(gòu);承托槽;光源單元

    0 引 言

    馬鈴薯是世界四大主糧之一,營養(yǎng)豐富,素有“地下蘋果”、“第二面包”的美譽(yù)。近年來,馬鈴薯加工業(yè)在中國迅猛發(fā)展,正在向規(guī)模化,產(chǎn)業(yè)化以及高技術(shù)水平方向發(fā)展,但是馬鈴薯原薯品質(zhì)合格率低等因素限制了馬鈴薯加工業(yè)的快速發(fā)展[1]。在諸多種類的馬鈴薯病害中,黑心病是一種主要的內(nèi)部病害,會明顯降低馬鈴薯的食用、種用及加工用價值[2],嚴(yán)重影響原薯品質(zhì)合格率。在生產(chǎn)線上,由于馬鈴薯黑心病個體難以憑肉眼檢出,混入加工原料中會嚴(yán)重破壞整體的品質(zhì),因此開發(fā)馬鈴薯黑心病檢測機(jī)構(gòu)是有必要的,可以在二次加工前對原料馬鈴薯進(jìn)行逐一篩選,剔除黑心病個體,提高原料薯的合格率。

    國內(nèi)外針對馬鈴薯等果蔬農(nóng)產(chǎn)品所進(jìn)行的品質(zhì)無損檢測研究主要是使用基于可見/近紅外光譜技術(shù)[3-7]、高光譜圖像技術(shù)[8-11]或機(jī)器視覺技術(shù)[12-16]等設(shè)計(jì)的光學(xué)檢測平臺。根據(jù)光源的照射方式,光學(xué)檢測技術(shù)主要分為透射式檢測技術(shù)和反射式檢測技術(shù)[17-18]。López-Maestresalas等[19]基于反射式高光譜圖像技術(shù)設(shè)計(jì)了馬鈴薯黑心病檢測機(jī)構(gòu),分別分析了在400~1 000和1 000~2 500 nm波段范圍內(nèi)采集的高光譜反射圖像,然后建立了偏最小二乘判別分析模型,在靜態(tài)環(huán)境下可以較準(zhǔn)確地判斷出樣品中的黑心病個體。在國內(nèi),Zhou等[20]根據(jù)513~850 nm波段范圍內(nèi)馬鈴薯的透射式可見/近紅外光譜特征判斷其內(nèi)部黑心病,對黑心病的判別準(zhǔn)確率為97.11%。金瑞等[21]和李小昱等[22]分別應(yīng)用反射式高光譜圖像及光譜信息融合技術(shù)和透射式機(jī)器視覺技術(shù)進(jìn)行馬鈴薯內(nèi)外部多品質(zhì)指標(biāo)的檢測,對內(nèi)部黑心以及發(fā)芽、綠皮等特征的判別有較高的準(zhǔn)確率。上述針對馬鈴薯黑心病等品質(zhì)檢測所做的研究,主要是基于所搭建的靜態(tài)環(huán)境下的檢測機(jī)構(gòu),被測馬鈴薯樣品的物理狀態(tài)穩(wěn)定,干擾因素少,但是其檢測機(jī)構(gòu)還未與工業(yè)在線生產(chǎn)需求相匹配。另外,由于馬鈴薯形狀和尺寸不規(guī)則,在進(jìn)行反射式或者透射式光學(xué)檢測的時候,容易產(chǎn)生光照不均勻,漏光等現(xiàn)象,對品質(zhì)檢測造成較大的干擾。

    馬鈴薯黑心病主要表現(xiàn)為在薯塊中心的髓部組織細(xì)胞死亡,病變區(qū)域的塊莖組織極其硬,形狀呈不規(guī)則狀,但邊緣清晰,其光透過能力與馬鈴薯塊莖的正常組織不同。因此本文基于透射式機(jī)器視覺技術(shù)的原理設(shè)計(jì)了馬鈴薯黑心病檢測機(jī)構(gòu),通過分析馬鈴薯的尺寸和形態(tài)特征確定檢測機(jī)構(gòu)中光源單元的結(jié)構(gòu)及角度,并相應(yīng)確定馬鈴薯塊莖承托槽的結(jié)構(gòu)。然后分析比較單點(diǎn)透射光源下馬鈴薯塊莖樣品的透射圖像和使用本文所設(shè)計(jì)的檢測機(jī)構(gòu)采集的樣品透射圖像的灰度分布特征等參數(shù),調(diào)整檢測機(jī)構(gòu)到達(dá)最佳狀態(tài),為馬鈴薯黑心病在線檢測機(jī)構(gòu)的設(shè)計(jì)提供參考。

    1 馬鈴薯黑心病檢測機(jī)構(gòu)的結(jié)構(gòu)及原理

    本文所設(shè)計(jì)的馬鈴薯黑心病檢測機(jī)構(gòu)主要由馬鈴薯塊莖承托槽,光源單元,散熱片和底座四部分組成,結(jié)構(gòu)如圖1a所示。光源單元主要包括一對透射光源,其出射光會聚于承托槽的豎直中心上。承托槽上設(shè)計(jì)有可供光源單元位置調(diào)節(jié)的軌道(圖1b所示),用于適應(yīng)對不同尺寸馬鈴薯塊莖的檢測。

    圖1 檢測機(jī)構(gòu)的整體結(jié)構(gòu)示意圖Fig.1 Overall structure diagram of detection mechanism

    在黑心病檢測時,首先將馬鈴薯塊莖平穩(wěn)置于承托槽1上面,然后根據(jù)被測馬鈴薯的尺寸調(diào)節(jié)光源2在光源位置調(diào)節(jié)軌道5上的相對位置,使光源到樣品的垂直距離最近。由于所用光源為大功率LED光源,檢測時需要使用散熱片3配合散熱,保持光源發(fā)光能量的穩(wěn)定性。整個檢測機(jī)構(gòu)通過底座4進(jìn)行固定安置,經(jīng)過適當(dāng)調(diào)節(jié)可匹配在線系統(tǒng)的工作。

    2 關(guān)鍵部件的結(jié)構(gòu)設(shè)計(jì)及參數(shù)分析

    2.1 馬鈴薯塊莖尺寸及形狀特征分析

    由于栽培條件等因素的不同,即使為同一個品種的馬鈴薯,且生長在同一植株上,其形狀和尺寸也會有差異,只能根據(jù)統(tǒng)計(jì)學(xué)規(guī)律進(jìn)行歸納[23-26]。在工業(yè)生產(chǎn)中,為獲得高品質(zhì)的馬鈴薯加工產(chǎn)品,人們對馬鈴薯塊莖的大小和形狀均有嚴(yán)格的要求。根據(jù)農(nóng)業(yè)部標(biāo)準(zhǔn)NY/T 1605-2008《中華人民共和國農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)——加工用馬鈴薯 油炸》的規(guī)定,加工用馬鈴薯的薯形要求為近似圓形或長橢圓形,即側(cè)面形狀同樣為近似圓形或橢圓形;直徑要求為4.0~10.0 cm[27]。根據(jù)該要求,本文所設(shè)計(jì)機(jī)構(gòu)的檢測對象主要為側(cè)面形狀是近似圓形或橢圓形,直徑范圍在40.0~100.0 mm的馬鈴薯塊莖。圖2所示為馬鈴薯塊莖自然放置時側(cè)面尺寸測量示意圖,本文使用擬合高度與擬合寬度的比值(H/W)表示樣品側(cè)視形狀的圓形度。

    2.2 光源型號的選擇

    馬鈴薯黑心病是由于塊莖內(nèi)部中心組織病變導(dǎo)致的生理性病害,根據(jù)文獻(xiàn)[28]研究發(fā)現(xiàn),正常馬鈴薯組織和黑心病組織的光譜透射率不同,黑心病組織對光的吸收能力較正常組織更強(qiáng),尤其是在700 nm波長附近處,二者光吸收能力的差別最為明顯。因此,本文選擇波長為705 nm的大功率LED陣列作為檢測機(jī)構(gòu)的透射光源,所用光源為圓形面光源。

    圖2 馬鈴薯塊莖側(cè)面圖像的尺寸Fig.2 Potato tuber size measured in side image

    2.3 承托槽及光源單元結(jié)構(gòu)的設(shè)計(jì)

    本文首先針對側(cè)面形狀為近似圓形的馬鈴薯塊莖的檢測機(jī)構(gòu)進(jìn)行設(shè)計(jì)。為提高光源透射樣品的能量分布均勻性,選用2個LED面光源作為馬鈴薯黑心病檢測的透射光源,取馬鈴薯塊莖左下方和右下方兩部分區(qū)域?yàn)楣庠慈肷潼c(diǎn)。

    為使光源發(fā)光能量的利用率達(dá)到最大,且盡量減少光源漏光,將光源位置設(shè)置在馬鈴薯塊莖與承托槽相切切點(diǎn)的位置上,并假設(shè)不同尺寸的馬鈴薯塊莖的擬合圓共有同一條外切線,因此光源位置可以根據(jù)馬鈴薯尺寸的變化在該切線位置上進(jìn)行調(diào)節(jié)。如圖3所示,兩個光源的位置A1(A2)和B1(B2)由2個因素決定:1)為被測馬鈴薯塊莖樣品的側(cè)面擬合圓半徑;2)為與樣品擬合圓的切線相對應(yīng)的半徑法線O1A1(O2A2)和豎直中心線l的夾角α值的大小。在該透射式檢測機(jī)構(gòu)中,α的取值范圍為0~90°,α為0°時樣品擬合圓的切線與中心線l垂直,為90°時樣品擬合圓的切線與中心線l平行。

    在多光源的透射式光學(xué)檢測機(jī)構(gòu)中,光源與樣品的相對位置會影響光源能量在樣品中的分布,從而影響檢測結(jié)果的準(zhǔn)確性。馮尚坤等[29]在進(jìn)行基于透射式可見/近紅外光譜檢測蘋果內(nèi)可溶性固形物成分的研究時,分析比較了不同光源擺放形式下的檢測結(jié)果,發(fā)現(xiàn)光源與檢測探頭之間的夾角越小,且為多光源組合時的光照模式最佳。因此在本文中,α角度值越大,光源在樣品中的能量分布越均勻,越有利于獲得高質(zhì)量的馬鈴薯塊莖樣品透射圖像。

    同時,在相同α角度下,大尺寸馬鈴薯塊莖的中心比小尺寸馬鈴薯塊莖的中心在承托槽上的相對位置更高。根據(jù)三角函數(shù)公式,在圖3中,由于

    式中O1M為任一半徑大于最小半徑的馬鈴薯塊莖擬合圓圓心O1到切線交點(diǎn)M的距離,mm;O2M為最小半徑的馬鈴薯塊莖擬合圓圓心O2到切線交點(diǎn)M的距離,mm;O1O2為2個擬合圓的圓心距,mm;O1A1為圓心O1到切點(diǎn)A1的距離,mm;O2A2為圓心O2到切點(diǎn)A2的距離,mm;r1為任一尺寸的馬鈴薯塊莖的半徑,mm;r2為最小尺寸的馬鈴薯塊莖的半徑,mm;α為與樣品擬合圓的切線相對應(yīng)的半徑法線O1A1(O2A2)和豎直中心線l的夾角,(°)。其中r2=20.0 mm。

    圖3 光源位置設(shè)計(jì)Fig.3 Design of light source position

    由此得到擬合圓O1到擬合圓O2上最低點(diǎn)C2的距離

    式中O1C2為任一尺寸的馬鈴薯塊莖的中心到最小擬合圓尺寸的馬鈴薯塊莖底端的距離,mm;O2C2為最小半徑的馬鈴薯塊莖的中心到底端的距離,mm;O1O2為2個擬合圓的圓心距,mm。

    從式(2)中可以看出,O1C2由夾角α值和該馬鈴薯的擬合圓半徑r1決定。當(dāng)r1一定時,α越小,O1C2越小,馬鈴薯樣品的擬合圓心O1相對于最低點(diǎn)C2的高度差越小。同理,當(dāng)α角度越大,O1C2的值越大,此時不同尺寸的馬鈴薯塊莖樣品到機(jī)器視覺成像裝置的距離的差別更大,不利于光源能量的充分利用。因此,本文中α角度適當(dāng)減小有利于提高不同樣品檢測時的光源能量利用率。

    本文令夾角α的值為45°,此時,2個光源在樣品中的能量分布均勻性以及樣品擬合圓心相對最小尺寸樣品(r2=20.0 mm)底端的高度差在理論上達(dá)到了一個協(xié)調(diào)的狀態(tài),得到承托槽的結(jié)構(gòu)如圖4所示。圖4a中夾角β值的大小為45°,線段BC和EF分別為任意尺寸馬鈴薯擬合圓的外切線,圓弧︵CE設(shè)計(jì)為承托槽的底面,根據(jù)尺寸最小的馬鈴薯擬合圓得到,其半徑為20.0 mm;圓弧︵AB和︵FG根據(jù)尺寸最大的馬鈴薯擬合圓得到,其半徑為50.0 mm;三段圓弧分別與兩段切線在交點(diǎn)處相切。圖4b為根據(jù)該結(jié)構(gòu)所作的機(jī)械結(jié)構(gòu)示意圖,P為相應(yīng)的光源位置調(diào)節(jié)軌道。

    圖4 馬鈴薯塊莖承托槽設(shè)計(jì)Fig.4 Design of supporting groove for potato tuber

    在該結(jié)構(gòu)的承托槽上,可用于檢測擬合圓直徑范圍為40.0~100.0 mm的馬鈴薯塊莖樣品;兩個光源的位置可根據(jù)被測樣品的尺寸范圍分別在BC段和EF段進(jìn)行調(diào)節(jié),使得光源位于樣品與承托槽的切點(diǎn)上。根據(jù)公式(2),馬鈴薯塊莖擬合圓半徑為50.0 mm時,其圓心O′到承托槽最低點(diǎn)D的距離O′D=62.43 mm;當(dāng)擬合圓半徑為20.0 mm時,其圓心O到承托槽最低點(diǎn)C的距離OC= 20.0 mm。該結(jié)構(gòu)下馬鈴薯塊莖擬合圓圓心相對承托槽最低點(diǎn)的高度變化范圍為20.0~62.43 mm。

    另外,在承托槽上,為了提高光源能量的利用率,將面光源的出光孔直徑設(shè)計(jì)為15.0 mm。圖5所示為承托槽光源相對馬鈴薯塊莖擬合圓的位置示意。圖中N為光源A0B0的中點(diǎn),且為樣品擬合圓在承托槽上的切點(diǎn)。

    光源上半段A0N在水平線上的投影ED與擬合圓O0水平方向的半徑FO0的位置關(guān)系如公式(3)所示

    式中FO0為擬合圓O0在水平方向的半徑,mm;FE為光源上端點(diǎn)A0在FO0上的投影點(diǎn)到擬合圓O0左端點(diǎn)F的距離,mm;r表示擬合圓O0的半徑,mm;θ為與樣品擬合圓的切線相對應(yīng)的半徑法線O0N和豎直中心線l的夾角,(°)。其中θ=45°。

    圖5 承托槽光源相對馬鈴薯塊莖擬合圓的位置示意圖Fig.5 Schematic diagram of light sources’ position of supporting groove relevant to potato tuber fitting circle

    由此得到光源上端點(diǎn)A0與馬鈴薯擬合圓左端點(diǎn)F的相對位置如公式(4)所示

    當(dāng)FE>0時,樣品擬合圓半徑r>18.12 mm。本文中用于黑心病檢測的馬鈴薯樣品最小半徑為20.0 mm,滿足該區(qū)間要求。此時,光源入射到被測馬鈴薯樣品上的能量可以充分應(yīng)用于樣品透射,可減少光源漏光以及能量浪費(fèi)和干擾,保證了所采集的透射圖像的亮度分布均勻性。

    當(dāng)所檢測的馬鈴薯塊莖樣品側(cè)面擬合形狀為橢圓時,根據(jù)數(shù)理學(xué)統(tǒng)計(jì),幾種主要品種馬鈴薯塊莖的擬合高度H與擬合寬度W的比值范圍為0.6~1.0[30]。如圖6所示,分別為圓擬合和橢圓擬合馬鈴薯樣品,圓擬合樣品的半徑長度與橢圓擬合樣品的長半軸長度相同。從圖中看出,隨著橢圓擬合的馬鈴薯塊莖的短軸與長軸的比值變小,其中心到承托槽最低點(diǎn)的距離OminD2減小,此時擬合橢圓Omin與承托槽兩側(cè)切點(diǎn)N0的位置相對于擬合圓Omax的切點(diǎn)N的位置抬高。因此,將本文根據(jù)圓擬合樣品所設(shè)計(jì)的檢測機(jī)構(gòu)用于橢圓形馬鈴薯塊莖樣品的檢測時,樣品同樣可以獲得分布均勻的光源能量,且其中心位置相對于承托槽底面的高度變化范圍在圓擬合樣品中心的高度變化范圍內(nèi),可以保證圖像采集時的光源能量利用率的提高。

    圖6 圓擬合和橢圓擬合的馬鈴薯塊莖樣品參數(shù)比較Fig.6 Parameters comparison between fitting circle and fitting ellipse of potato tuber

    3 性能試驗(yàn)設(shè)計(jì)及分析

    3.1 樣品準(zhǔn)備

    試驗(yàn)所用馬鈴薯品種為克新四號,樣品表皮新鮮,表面無機(jī)械損傷,無病蟲害。根據(jù)馬鈴薯黑心病生成的原理對樣品進(jìn)行處理。馬鈴薯塊莖在高溫缺氧的環(huán)境下容易產(chǎn)生黑心病,在制作黑心病時,首先將樣品使用自封袋單個封裝,封裝時排掉袋中多余的空氣,然后放入恒溫箱中貯存,恒溫箱的溫度被設(shè)定為37 ℃,處理24 d后取出。

    取79個樣品(50個正常樣品,29個制作的黑心病樣品)進(jìn)行試驗(yàn)。使用精度為0.01 mm的游標(biāo)卡尺測量樣品的尺寸,統(tǒng)計(jì)結(jié)果如表1所示。從表中可以看出,所選樣品的尺寸分布范圍比較均衡,形態(tài)特征包括橢圓擬合體和圓擬合體。

    表1 馬鈴薯塊莖樣品尺寸統(tǒng)計(jì)結(jié)果Table1 Statistical result of potato tuber samples’ sizes

    3.2 試驗(yàn)儀器

    本文使用試制的檢測機(jī)構(gòu)進(jìn)行試驗(yàn)。所用光源為兩個中心波長是705 nm的LED圓形面光源,功率均為50 W,光源在承托槽上的出光位置可以沿調(diào)節(jié)軌道進(jìn)行調(diào)節(jié)。采用Mako G-125B黑白相機(jī)采集樣品的透射圖像,相機(jī)的分辨率為1 292×964像素,相對檢測機(jī)構(gòu)的安裝高度為200 mm。使用Matlab R2012a軟件進(jìn)行圖像處理和數(shù)據(jù)分析,從而評價所設(shè)計(jì)檢測機(jī)構(gòu)的工作性能。

    3.3 試驗(yàn)方法

    首先使用本文設(shè)計(jì)的檢測機(jī)構(gòu)進(jìn)行多點(diǎn)透射(multipointtransmission,MT)模式(圖7a所示)下的試驗(yàn),依次采集79個馬鈴薯樣品的透射圖像,試驗(yàn)中根據(jù)樣品尺寸調(diào)節(jié)檢測機(jī)構(gòu)上光源的位置。然后,為比較該機(jī)構(gòu)透射光源的能量在樣品中的分布均勻性,進(jìn)行單點(diǎn)透射模式(single-point transmission,ST)下的試驗(yàn),光源位置如圖7b所示,將兩個光源并排組合置放在水平面上,使其位于相機(jī)視野的豎直中心下方,調(diào)節(jié)相機(jī)和光源的相對高度與多點(diǎn)透射模式時一致。由于黑心馬鈴薯發(fā)生病變,內(nèi)部組織的密度不均勻,且不同樣品發(fā)生黑心病變的部位各有不同,對光源能量在組織內(nèi)部的正常傳播造成干擾,不利于光源透射均勻性的分析,因此只依次采集50個正常樣品的透射圖像用于分析。

    圖7 光源透射模式示意圖Fig.7 Schematic diagram of light source transmission mode

    對所采集的圖像分別在Matlab軟件中進(jìn)行去背景處理,提取出馬鈴薯目標(biāo)圖像。根據(jù)公式(5)計(jì)算馬鈴薯目標(biāo)圖像的灰度平均值。

    式中AVG表示灰度平均值;pixel(i, j)表示馬鈴薯目標(biāo)區(qū)域每個像素點(diǎn)的灰度值,其中i和j分別為圖像上橫向和縱向的像素點(diǎn)坐標(biāo);N表示圖像上馬鈴薯目標(biāo)區(qū)域的像素點(diǎn)總數(shù)。

    以所得灰度平均值A(chǔ)VG為閾值對每幅圖像進(jìn)行二值化處理,并以所得二值圖像為掩膜提取馬鈴薯圖像的高灰度區(qū)域。對所得的圖像按公式(5)分別計(jì)算目標(biāo)區(qū)域的灰度平均值。然后按照公式(6)計(jì)算灰度平均值的相對比值。

    式中RT表示灰度平均值的相對比值;AVG0和AVG1分別表示馬鈴薯原目標(biāo)圖像和二值化掩膜圖像的灰度平均值。

    3.4 結(jié)果及討論

    對所采集的79個馬鈴薯樣品(50個正常樣品,29個黑心樣品)圖像進(jìn)行處理之后,所得結(jié)果如圖8所示。從圖中可以看出,ST模式采集的正常樣品圖像灰度分布不均勻,中心較四周邊緣更亮;MT模式下的正常樣品圖像灰度分布相對更均勻,黑心樣品的圖像灰度分布不均勻,灰度值低的區(qū)域?qū)?yīng)于內(nèi)部發(fā)生黑心病變的區(qū)域。根據(jù)公式(5)分別計(jì)算所有樣品圖像中馬鈴薯原目標(biāo)圖像和二值化掩膜圖像的灰度平均值,并按照公式(6)計(jì)算灰度平均值的相對比值,圖9所示為樣品的灰度平均值的相對比值統(tǒng)計(jì)坐標(biāo)圖。

    圖8 兩種透射模式下的馬鈴薯樣品原圖像及二值化掩膜圖像Fig.8 Images of potato samples and after binarization and masking under two transmission modes

    圖9 樣品圖像的灰度平均值比例Fig.9 Grayscale distribution of sample images under two transmission modes

    根據(jù)圖9可以看出,三類樣品的灰度平均值相對比值RT存在明顯的分布區(qū)間,其中正常樣品在MT模式所采集的數(shù)據(jù),分布范圍為0.012~0.250;正常樣品在ST模式所采集的數(shù)據(jù),分布范圍為0.210~0.497;黑心樣品在MT模式所采集的數(shù)據(jù),分布范圍為0.204~0.433。取ST模式和MT模式下正常樣品RT值的重合點(diǎn)以及重合區(qū)間等分點(diǎn)為分類閾值,分析兩種模式下樣品的分類正確率以及分類總正確率,結(jié)果如圖10a所示。同樣,取MT模式下正常樣品和黑心樣品RT值的重合點(diǎn)及重合區(qū)間等分點(diǎn)為分類閾值,分析兩種樣品的分類正確率及總正確率,結(jié)果如圖10b所示。

    圖10 樣品的最佳分類閾值分析Fig.10 Analysis of samples’ optimal classification threshold

    從圖10a中可以看出,當(dāng)分類閾值為0.216時,正常樣品在ST和MT模式下的分類正確率均為98%,且總正確率達(dá)到最大,為98%,因此取0.216為最佳分類閾值。在圖10b中,當(dāng)分類閾值為0.220時,MT模式下正常樣品和黑心樣品的分類正確率分別為98%和96.6%,此時總正確率最大,為97.5%,因此取0.220為最佳分類閾值。

    以上分析結(jié)果表明,該檢測機(jī)構(gòu)能夠用于馬鈴薯黑心病的檢測,為在線加工中的原料薯品質(zhì)檢測裝置的設(shè)計(jì)提供了參考。檢測中由于馬鈴薯形狀不規(guī)則等因素所造成的漏光或圖像灰度不均勻等問題,可以通過提高圖像處理算法的處理精度進(jìn)行改善;同時增加預(yù)測參數(shù),進(jìn)一步提高黑心病判斷結(jié)果的可靠性。

    4 結(jié) 論

    本文為適應(yīng)工業(yè)生產(chǎn)中馬鈴薯在線品質(zhì)檢測需求,基于馬鈴薯塊莖的光透射特性設(shè)計(jì)了馬鈴薯黑心病的檢測機(jī)構(gòu)。

    1)通過分析馬鈴薯的外形和尺寸特征,將馬鈴薯塊莖側(cè)面形狀擬合為圓,設(shè)計(jì)了樣品擬合直徑范圍為40.0~100.0 mm時光源單元和承托槽的結(jié)構(gòu);檢測機(jī)構(gòu)上光源的位置可以根據(jù)被測馬鈴薯塊莖樣品的尺寸進(jìn)行調(diào)整。

    2)通過數(shù)學(xué)計(jì)算,確定了當(dāng)擬合圓切線的法線與豎直中心線的夾角為45°時,圓心相對承托槽底部的高度變化范圍為20~62.43 mm,光源能量可被更充分的應(yīng)用于樣品的透射,提高了光透射的均勻性,并且盡量減少了樣品邊緣處的光能量泄漏。

    3)針對一些側(cè)面擬合形狀為橢圓的馬鈴薯塊莖樣品,分析了其與承托槽的切點(diǎn)位置及其中心的高度變化范圍。通過和圓擬合樣品比較,確定了上述指標(biāo)均在圓擬合樣品的范圍內(nèi),應(yīng)用該檢測機(jī)構(gòu)可以采集橢圓形狀擬合的馬鈴薯樣品透射圖像。

    4)經(jīng)試驗(yàn)驗(yàn)證,在相同光源功率的條件下,使用該檢測機(jī)構(gòu)(多點(diǎn)透射模式)檢測正常樣品和單點(diǎn)透射模式下檢測正常樣品時灰度平均值相對比值的分類閾值為0.216,此時樣品的判斷正確率均為98%,小于該閾值的樣品數(shù)據(jù)為該檢測機(jī)構(gòu)獲得,表明該檢測機(jī)構(gòu)成像具有更好的灰度均勻性;另外,該檢測機(jī)構(gòu)檢測正常樣品和黑心樣品的分類閾值為0.220,此時正常樣品和黑心樣品的判斷正確率分別為98%和96.6%,大于該閾值的樣品為黑心樣品,表明該檢測機(jī)構(gòu)可用于檢測黑心馬鈴薯樣品。該機(jī)構(gòu)為馬鈴薯黑心病檢測機(jī)構(gòu)的設(shè)計(jì)提供了參考。

    [1] 孫東升. 中國馬鈴薯加工業(yè)現(xiàn)狀和發(fā)展趨勢[J]. 農(nóng)業(yè)展望,2009,5(11):20-22.

    [2] 孟令文. 馬鈴薯黑心病的發(fā)病條件及安全貯藏技術(shù)[J]. 園藝與種苗,1994(6):51-52.

    [3] Rady A M, Guyer D E. Rapid and/or nondestructive quality evaluation methods for potatoes: A review[J]. Computers and Electronics in Agriculture, 2015, 117(1): 31-48.

    [4] 郭志明,黃文倩,陳全勝,等. 蘋果腐心病的透射光譜在線檢測系統(tǒng)設(shè)計(jì)及試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報,2016,32(6):283-288.

    Guo Zhiming, Huang Wenqian, Chen Quansheng, et al. Design and test of on-line detection system for apple core rot disease based on transmitted spectrum[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(6): 283-288. (in Chinese with English abstract)

    [5] Vanoli M, Rizzolo A, Spinelli L, et al. Non-destructive detection of internal brown spot in potato tubers by time-resolved reflectance spectroscopy: preliminary results on a susceptible cultivar[C]//International Conference of Agricultural Engineering. CIGR-Ageng 2012, Valencia, Spain, 2012, Article No. P1370.

    [6] Alander J, Bochko V, Martinkauppi J B, et al. Optical sensoring of internal hollow heart related defects of potatoes[J]. IFAC Proceedings Volumes, 2013, 46(18): 24-28.

    [7] 王欣,謝錦春,韓東海,等. 水果內(nèi)部品質(zhì)在線近紅外分析儀的研制[J]. 現(xiàn)代科學(xué)儀器,2009,12(6):11-13.

    Wang Xin, Xie Jinchun, Han Donghai, et al. Development of near infrared spectroscopy on-line analyzer for fruit internal quality[J]. Modern Scientific Instruments, 2009, 12(6): 11-13. (in Chinese with English abstract)

    [8] Dacal-Nieto A, Formella A, Carrión P, et al. Non–destructive detection of hollow heart in potatoes using hyperspectral imaging[C]//International Conference on Computer Analysis of Images and Patterns. Springer Berlin Heidelberg, 2011: 180-187.

    [9] 胡鵬程,孫曄,吳海倫,等. 高光譜圖像對白蘿卜糠心的無損檢測[J]. 食品科學(xué),2015,36(12):171-176.

    Hu Pengcheng, Sun Ye, Wu Hailun, et al. Detecting hollowness of white radish based on hyperspectral imaging[J]. Food Science, 2015, 36(12): 171-176. (in Chinese with English abstract)

    [10] 張嬙,潘磊慶,吳林蔚,等. 利用反射和半透射高光譜圖像檢測水蜜桃早期冷害[J]. 食品科學(xué),2014,35(4):71-76.

    Zhang Qiang, Pan Leiqing, Wu Linwei, et al. Detecting chilling injury at early stage of ‘xiahui 5’ honey peach by hyperspectral reflectance and half-transmittance imaging[J]. Food Science, 2014, 35(4): 71-76. (in Chinese with English abstract)

    [11] Keresztes J C, Goodarzi M, Saeys W. Real-time pixel based early apple bruise detection using short wave infrared hyperspectral imaging in combination with calibration and glare correction techniques[J]. Food Control, 2016, 66(1): 215-226.

    [12] Noordam J C, Otten G W, Timmermans T J M, et al. High-speed potato grading and quality inspection based on a color vision system[C]//Electronic imaging. International Society for Optics and Photonics, 2000: 206-217.

    [13] 郝敏,麻碩士. 基于機(jī)器視覺的馬鈴薯單薯質(zhì)量檢測技術(shù)研究[J]. 農(nóng)機(jī)化研究,2009(9):61-63.

    Hao Min, Ma Shuoshi. Study on potato weight detection technology based on computer vision theory[J]. Journal of Agricultural Mechanization Research, 2009(9): 61-63. (in Chinese with English abstract)

    [14] Nandi C S, Tudu B, Koley C. An automated machine vision based system for fruit sorting and grading[C]//Sixth International Conference on Sensing Technology. IEEE, 2012: 195-200.

    [15] Cubero S, Aleixos N, Moltó E, et al. Advances in machine vision applications for automatic inspection and quality evaluation of fruits and vegetables[J]. Food and Bioprocess Technology, 2011, 4(4): 487-504.

    [16] Mizushima A, Lu R. An image segmentation method for apple sorting and grading using support vector machine and Otsu’s method[J]. Computers and Electronics in Agriculture, 2013, 94(1): 29-37.

    [17] 申曉彥,王鑒. 用于視覺檢測的光源照明系統(tǒng)分析[J]. 燈與照明,2009,33(3):7-9.

    Shen Xiaoyan, Wang Jian. The light source illumination system analysis used in visual inspection[J]. Light & Lighting, 2009, 33(3): 7-9. (in Chinese with English abstract)

    [18] 徐惠榮,應(yīng)義斌. 近紅外圖像和光譜分析技術(shù)在農(nóng)產(chǎn)品品質(zhì)無損檢測中的應(yīng)用與展望[J]. 浙江大學(xué)學(xué)報:農(nóng)業(yè)與生命科學(xué)版,2002,28(4):460-464.

    Xu Huirong, Ying Yibin. Application and prospect of near infrared imaging and spectroscopy analysis on quality inspecting of agricultural products[J]. Journal of Zhejiang University: Agric & Life Sci, 2002, 28(4): 460-464. (in Chinese with English abstract)

    [19] López-Maestresalas A, Keresztes J C, Goodarzi M, et al. Non-destructive detection of blackspot in potatoes by Vis-NIR and SWIR hyperspectral imaging[J]. Food Control, 2016, 70(1): 229-241.

    [20] Zhou Z, Zeng S, Li X, et al. Nondestructive detection of blackheart in potato by visible/near infrared transmittance spectroscopy[J]. Journal of Spectroscopy, 2015, 2015(1): 1-9.

    [21] 金瑞,李小昱,顏伊蕓,等. 基于高光譜圖像和光譜信息融合的馬鈴薯多指標(biāo)檢測方法[J]. 農(nóng)業(yè)工程學(xué)報,2015,31(16):258-263.

    Jin Rui, Li Xiaoyu, Yan Yiyun, et al. Detection method of multi-target recognition of potato based on fusion of hyperspectral imaging and spectral information[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(16): 258-263. (in Chinese with English abstract)

    [22] 李小昱,陶海龍,高海龍,等. 馬鈴薯缺陷透射和反射機(jī)器視覺檢測方法分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2014,45(5):191-196.

    Li Xiaoyu, Tao Hailong, Gao Hailong, et al. Comparison of transmission and reflection imaging technologies to detect potato defects based on machine vision technology[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(5): 191-196. (in Chinese with English abstract)

    [23] 石林榕,趙武云,吳建民,等. 切片造型技術(shù)在農(nóng)產(chǎn)品有限元分析中的應(yīng)用[J]. 中國農(nóng)機(jī)化學(xué)報,2013,34(3):95-98.

    Shi Linrong, Zhao Wuyun, Wu Jianming, et al. Application of slice modeling technology in finite element analysis of agricultural products[J]. Journal of Chinese Agricultural Mechanization, 2013, 34(3): 95-98. (in Chinese with English abstract)

    [24] 鄭冠楠,譚豫之,張俊雄,等. 基于計(jì)算機(jī)視覺的馬鈴薯自動檢測分級[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2009,40(4):166-168,156.

    Zheng Guannan, Tan Yuzhi, Zhang Junxiong, et al. Automatic detecting and grading method of potatoes with computer vision[J]. Transactions of the Chinese Society for Agricultural Machinery, 2009, 40(4): 166-168, 156. (in Chinese with English abstract)

    [25] Tabatabaeefar A. Size and shape of potato tubers[J]. International Agrophysics, 2002, 16(4): 301-306.

    [26] Tao Y, Morrow C T, Heinemann P H, et al. Fourier-based separation technique for shape grading of potatoes using machine vision[J]. Transactions of the ASAE, 1995, 38(3): 949-957.

    [27] NY/T 1605-2008. 加工用馬鈴薯油炸[S].

    [28] 周竹,李小昱,高海龍,等. 漫反射和透射光譜檢測馬鈴薯黑心病的比較[J]. 農(nóng)業(yè)工程學(xué)報,2012,28(11):237-242.

    Zhou Zhu, Li Xiaoyu, Gao Hailong, et al. Comparison of diffuse reflection and transmission mode of visible/near infrared spectroscopy for detecting black heart of potato[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(11): 237-242. (in Chinese with English abstract)

    [29] 馮尚坤,徐海菊. 基于可見/近紅外光譜無損檢測蘋果可溶性固形物的光照優(yōu)化[J]. 食品工業(yè)科技,2014,35(16):64-66.

    Feng Shangkun, Xu Haiju. Optimization of light distribution in non-destructive analysis of soluble solids content of apple base on visual/near-infrared spectroscopy[J]. Science and Technology of Food Industry, 2014, 35(16): 64-66. (in Chinese with English abstract)

    [30] 劉春香. 馬鈴薯塊莖外形與力學(xué)流變學(xué)性質(zhì)研究與應(yīng)用[D].哈爾濱:東北農(nóng)業(yè)大學(xué),2006.

    Liu Chunxiang. Studies and Application on the Shape and Mechanical Rheological Properties of Potato Tuber[D]. Harbin: Northeast Agricultural University, 2006. (in Chinese with English abstract)

    Design and experiment of detection mechanism for potato blackheart based on machine vision

    Tian Fang, Peng Yankun※, Wei Wensong, Zheng Xiaochun, Wang Wenxiu
    (National R&D Center for Agro-Processing Equipments, College of Engineering, China Agricultural University, Beijing 100083, China)

    Potato blackheart is a kind of disease that happens in the inner of potato tuber. The potato with blackheart is difficult to be distinguished from the surface with the naked eye, which affects the quality of potato processed products seriously. In order to detect the potato with blackheart on the production line, a detection mechanism for potato blackheart was designed based on the light transmission properties of potato tuber internal organization in this study. The detection mechanism was composed of the supporting groove, the light source unit, the cooling fin, and the base. Firstly, the characteristics of size and morphology of potato tuber were analyzed based on the statistics rule. The potato tuber studied in this paper was round or elliptic fitted from the side direction, and the fitting diameter range was 40.0-100.0 mm. Then the light source wavelength was chosen as 705 nm according to the result of the research about potato black heart’s high spectral transmission detection. Two area light sources were used in the detection mechanism, constituting multipoint transmission structure. The mechanism was designed against the fitted circle of potato tuber sample in the first place. Two common tangent lines for any size of potato tubers were defined, and the locations of 2 tangent points on each fitting circle were set as the light source locations. The value of the included angle formed by the radius passing the tangent point and the vertical centerline was a determinant to the light source location. A proper included angle could lead to uniform illumination and less light leaking to the potato sample, and improve the consistency of imaging in the camera. The angle value was chosen as 45° by the analysis in this study. After that, the structure of the supporting groove was obtained. The height range between the center of the fitted circle of potato tuber and the bottom of the supporting groove was 20.0-62.43 mm. It could be accepted as a certain range to ensure the consistency of imaging. At the same time, the energy of the area light source installed on the supporting groove was completely absorbed by the potato tuber for its horizontal projection was exactly in the range of potato tuber horizontal radius. When the mechanism was used in the detection of fitted ellipse of potato tuber, the light source position, the included angle and the size of the supporting groove all met the detection demands by the analysis. In order to test the performance of the detection mechanism, a batch of fresh potato samples were chosen and sealed in zip-lock bags, then stored in an incubator at the temperature of 37°. Twenty-nine black heart samples were obtained after 24 d.79 potato tuber samples, including 50 normal samples and 29 black heart samples were chosen to collect the transmission images. The transmission images of 50 normal potato samples under single-point transmission mode were collected at the same time. Then the samples were sealed in zip-lock bags and stored in an incubator at the temperature of 37°. They were taken out and the light transmission images were collected by the detection mechanism after 24 d. Twenty-nine blackheart samples were obtained through the verification. After the background elimination processing with MATLAB, the potato tuber target image and the high grayscale value area after binarization and masking were extracted. The value of their average grayscale value relative ratio was acquired to distinguish the uniformity of the light transmission and potato black heart. The discriminant accuracy of the normal and black heart samples were 98% and 96.6%, respectively. Results showed that, the mechanism introduced in the paper could be used to detect potato black heart. A reference is provided by the design of the detection mechanism in this study.

    machine vision; nondestructive examination; design; potato tuber; black heart; detection mechanism; supporting groove; light source unit

    10.11975/j.issn.1002-6819.2017.05.041

    S226.5

    A

    1002-6819(2017)-05-0287-08

    田 芳,彭彥昆,魏文松,鄭曉春,王文秀. 基于機(jī)器視覺的馬鈴薯黑心病檢測機(jī)構(gòu)設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(5):287-294.

    10.11975/j.issn.1002-6819.2017.05.041 http://www.tcsae.org

    Tian Fang, Peng Yankun, Wei Wensong, Zheng Xiaochun, Wang Wenxiu. Design and experiment of detection mechanism for potato blackheart based on machine vision[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(5): 287-294. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.05.041 http://www.tcsae.org

    2016-08-02

    2017-01-17

    國家重點(diǎn)研發(fā)計(jì)劃課題(2016YFD0400905)

    田 芳,女(土家族),湖北人,博士生,主要從事農(nóng)畜產(chǎn)品質(zhì)安全無損檢測研究。北京 中國農(nóng)業(yè)大學(xué)工學(xué)院,100083。

    Email:tianfang@cau.edu.cn

    ※通信作者:彭彥昆,男(漢族),山東人,教授,博士生導(dǎo)師,主要從事農(nóng)畜產(chǎn)品品質(zhì)安全無損檢測研究。北京 中國農(nóng)業(yè)大學(xué)工學(xué)院,100083。Email:ypeng@cau.edu.cn

    猜你喜歡
    心病塊莖光源
    短期窖藏下馬鈴薯塊莖品質(zhì)指標(biāo)的變化趨勢
    中老年保健(2021年8期)2021-08-24 06:24:32
    《光源與照明》征訂啟事
    光源與照明(2019年3期)2019-06-15 09:21:04
    關(guān)注孩子的“心病”
    中華家教(2018年8期)2018-09-25 03:23:00
    古人醫(yī)心病
    特別健康(2018年4期)2018-07-03 00:38:22
    塊根塊莖類植物細(xì)胞懸浮培養(yǎng)技術(shù)與應(yīng)用
    綠色光源
    白及須根與塊莖的多糖成分比較研究
    西紅柿的心病
    小說月刊(2015年9期)2015-04-23 08:48:17
    科技連載:LED 引領(lǐng)第三次照明革命——與傳統(tǒng)照明的對比(一)
    免费观看性生交大片5| 久久久久久久久久久免费av| 中文字幕人妻丝袜一区二区 | h视频一区二区三区| 亚洲四区av| 亚洲,一卡二卡三卡| 九色亚洲精品在线播放| 麻豆精品久久久久久蜜桃| 丝袜在线中文字幕| 国产午夜精品一二区理论片| 久久国产精品大桥未久av| 超色免费av| 成年女人毛片免费观看观看9 | 女人被躁到高潮嗷嗷叫费观| 亚洲国产毛片av蜜桃av| 午夜福利乱码中文字幕| 欧美日本中文国产一区发布| 亚洲国产精品一区二区三区在线| 一区在线观看完整版| 国产熟女午夜一区二区三区| 大片电影免费在线观看免费| 丰满迷人的少妇在线观看| 色视频在线一区二区三区| av女优亚洲男人天堂| 午夜影院在线不卡| 一级片免费观看大全| 亚洲av福利一区| av电影中文网址| 亚洲精品一区蜜桃| 久久人妻熟女aⅴ| 最近最新中文字幕免费大全7| 一二三四中文在线观看免费高清| 男人操女人黄网站| av国产精品久久久久影院| 激情视频va一区二区三区| 日韩在线高清观看一区二区三区| 亚洲精品日韩在线中文字幕| 色视频在线一区二区三区| 啦啦啦啦在线视频资源| 老汉色av国产亚洲站长工具| 视频区图区小说| 最近中文字幕高清免费大全6| 另类精品久久| 日韩熟女老妇一区二区性免费视频| 国产熟女午夜一区二区三区| 黄色毛片三级朝国网站| 高清在线视频一区二区三区| 久久久久国产网址| 国产精品 国内视频| 亚洲一码二码三码区别大吗| 老汉色∧v一级毛片| 欧美 日韩 精品 国产| 高清在线视频一区二区三区| 国产又色又爽无遮挡免| 黄色配什么色好看| 国产av精品麻豆| 亚洲一区二区三区欧美精品| 2022亚洲国产成人精品| 哪个播放器可以免费观看大片| 婷婷色综合大香蕉| 可以免费在线观看a视频的电影网站 | 午夜福利视频在线观看免费| 黄色一级大片看看| 桃花免费在线播放| 国产成人av激情在线播放| 欧美日韩精品成人综合77777| a级毛片在线看网站| 只有这里有精品99| 亚洲成国产人片在线观看| 热re99久久国产66热| www.自偷自拍.com| 十八禁网站网址无遮挡| 九草在线视频观看| 欧美精品一区二区大全| 婷婷色综合大香蕉| 亚洲 欧美一区二区三区| 激情视频va一区二区三区| 久久久欧美国产精品| 久久久亚洲精品成人影院| 成人国语在线视频| 中文字幕av电影在线播放| 一二三四在线观看免费中文在| 黑人巨大精品欧美一区二区蜜桃| 欧美黄色片欧美黄色片| 欧美精品一区二区大全| 国产激情久久老熟女| 亚洲国产最新在线播放| 亚洲精品中文字幕在线视频| 午夜免费观看性视频| 国产一区二区激情短视频 | 亚洲人成77777在线视频| 国产精品无大码| 国产人伦9x9x在线观看 | 亚洲国产精品一区二区三区在线| 26uuu在线亚洲综合色| 国产日韩一区二区三区精品不卡| 丝袜脚勾引网站| 亚洲精品国产av蜜桃| 国产色婷婷99| 精品国产乱码久久久久久男人| av片东京热男人的天堂| 精品国产国语对白av| 国产深夜福利视频在线观看| 国产精品.久久久| 如日韩欧美国产精品一区二区三区| www.精华液| 日本欧美国产在线视频| 一个人免费看片子| 国产精品一区二区在线观看99| av免费观看日本| 久久久久久久精品精品| 婷婷成人精品国产| 中文字幕人妻丝袜一区二区 | 中文字幕人妻熟女乱码| 丰满迷人的少妇在线观看| 高清欧美精品videossex| 亚洲国产毛片av蜜桃av| 免费高清在线观看日韩| 波多野结衣一区麻豆| 伦理电影免费视频| 老司机亚洲免费影院| 97在线人人人人妻| 岛国毛片在线播放| 1024视频免费在线观看| av免费在线看不卡| 久久青草综合色| 亚洲五月色婷婷综合| 亚洲国产毛片av蜜桃av| 看十八女毛片水多多多| 香蕉国产在线看| 国产精品一区二区在线不卡| 熟妇人妻不卡中文字幕| 成人毛片a级毛片在线播放| 久久ye,这里只有精品| 在线观看www视频免费| 日韩熟女老妇一区二区性免费视频| av免费在线看不卡| 性色av一级| 看十八女毛片水多多多| 欧美少妇被猛烈插入视频| 建设人人有责人人尽责人人享有的| 亚洲精品美女久久av网站| 亚洲国产av新网站| 亚洲欧美精品综合一区二区三区 | 又粗又硬又长又爽又黄的视频| 欧美精品亚洲一区二区| videosex国产| 激情五月婷婷亚洲| 亚洲伊人色综图| 久久99一区二区三区| 91久久精品国产一区二区三区| 免费黄网站久久成人精品| 岛国毛片在线播放| 国产精品 欧美亚洲| 国产精品久久久久久精品古装| 国产免费现黄频在线看| 七月丁香在线播放| 国产成人免费无遮挡视频| 亚洲国产精品一区三区| 精品人妻偷拍中文字幕| 美女大奶头黄色视频| 青草久久国产| 中国三级夫妇交换| 久久午夜福利片| 18禁动态无遮挡网站| 欧美+日韩+精品| 国产精品偷伦视频观看了| 久久青草综合色| 飞空精品影院首页| 久久婷婷青草| 日本av手机在线免费观看| av电影中文网址| 久久久久久人人人人人| 五月开心婷婷网| 免费在线观看黄色视频的| av天堂久久9| 最近2019中文字幕mv第一页| 免费看不卡的av| 亚洲一码二码三码区别大吗| 精品国产乱码久久久久久小说| 中国国产av一级| 国产精品国产三级国产专区5o| 精品人妻在线不人妻| 一级片免费观看大全| 爱豆传媒免费全集在线观看| 久久 成人 亚洲| 久热久热在线精品观看| 国产亚洲av片在线观看秒播厂| 巨乳人妻的诱惑在线观看| 午夜免费观看性视频| 午夜福利网站1000一区二区三区| 中国三级夫妇交换| 一级毛片我不卡| 亚洲国产精品国产精品| 久久精品人人爽人人爽视色| 成人18禁高潮啪啪吃奶动态图| 欧美成人午夜免费资源| 国产精品 国内视频| 亚洲 欧美一区二区三区| 欧美97在线视频| 69精品国产乱码久久久| 少妇人妻精品综合一区二区| 久久久国产精品麻豆| 亚洲欧洲日产国产| av在线老鸭窝| 欧美bdsm另类| 久久久精品区二区三区| 天堂中文最新版在线下载| 国产成人av激情在线播放| 欧美精品高潮呻吟av久久| 91aial.com中文字幕在线观看| 又黄又粗又硬又大视频| 黄色一级大片看看| 久久午夜综合久久蜜桃| 国产乱人偷精品视频| 天天躁夜夜躁狠狠久久av| 久久韩国三级中文字幕| 久久精品熟女亚洲av麻豆精品| 美女主播在线视频| 日韩一区二区三区影片| 边亲边吃奶的免费视频| 亚洲国产最新在线播放| 国产成人精品福利久久| 中国国产av一级| 久久久国产欧美日韩av| 人人妻人人澡人人看| 欧美精品亚洲一区二区| 亚洲欧美一区二区三区黑人 | 亚洲,欧美精品.| 色婷婷av一区二区三区视频| 精品少妇一区二区三区视频日本电影 | 欧美精品高潮呻吟av久久| 国产 一区精品| 99久国产av精品国产电影| 欧美+日韩+精品| 又粗又硬又长又爽又黄的视频| 国产欧美日韩一区二区三区在线| 91午夜精品亚洲一区二区三区| 毛片一级片免费看久久久久| 欧美人与善性xxx| 熟女少妇亚洲综合色aaa.| 熟妇人妻不卡中文字幕| 哪个播放器可以免费观看大片| 国产精品熟女久久久久浪| 老汉色av国产亚洲站长工具| 国产日韩一区二区三区精品不卡| 日韩欧美一区视频在线观看| 成人国语在线视频| 国产精品久久久久久久久免| 国产极品天堂在线| 黄色视频在线播放观看不卡| 男女免费视频国产| 春色校园在线视频观看| 人人妻人人澡人人爽人人夜夜| 制服人妻中文乱码| h视频一区二区三区| 亚洲伊人久久精品综合| 中文字幕亚洲精品专区| 在线观看免费日韩欧美大片| 欧美日韩av久久| 久久国产精品大桥未久av| 欧美日韩视频高清一区二区三区二| 免费少妇av软件| 久久精品国产自在天天线| 久久精品夜色国产| 国产成人精品久久二区二区91 | 中文乱码字字幕精品一区二区三区| 看免费成人av毛片| 欧美少妇被猛烈插入视频| 丝袜喷水一区| 一级毛片 在线播放| 国产精品久久久久久久久免| 亚洲美女黄色视频免费看| 色婷婷久久久亚洲欧美| 熟女电影av网| 好男人视频免费观看在线| 十八禁网站网址无遮挡| 国产精品女同一区二区软件| 国产精品不卡视频一区二区| 国产人伦9x9x在线观看 | 国产成人av激情在线播放| 亚洲精品久久成人aⅴ小说| 菩萨蛮人人尽说江南好唐韦庄| 性色av一级| a级毛片在线看网站| 2021少妇久久久久久久久久久| 久久久久久久久久久免费av| 777米奇影视久久| 日韩中文字幕欧美一区二区 | 久久久久久久国产电影| 国产一区二区三区综合在线观看| 美女主播在线视频| 久久人人爽av亚洲精品天堂| 蜜桃国产av成人99| 性高湖久久久久久久久免费观看| 亚洲综合色惰| 免费av中文字幕在线| 国产日韩欧美视频二区| 国产成人一区二区在线| 国产一区二区在线观看av| 午夜福利影视在线免费观看| 欧美av亚洲av综合av国产av | 欧美日本中文国产一区发布| 中文字幕人妻熟女乱码| 国产成人免费无遮挡视频| 极品少妇高潮喷水抽搐| 中文天堂在线官网| 两个人看的免费小视频| 亚洲美女搞黄在线观看| 永久免费av网站大全| 色视频在线一区二区三区| 考比视频在线观看| 国产精品欧美亚洲77777| 国产日韩欧美在线精品| 亚洲国产欧美在线一区| 亚洲国产看品久久| 免费不卡的大黄色大毛片视频在线观看| 韩国精品一区二区三区| 久久久精品国产亚洲av高清涩受| www.自偷自拍.com| 日韩一卡2卡3卡4卡2021年| 午夜精品国产一区二区电影| 两性夫妻黄色片| 韩国av在线不卡| 中文天堂在线官网| 最近2019中文字幕mv第一页| 1024视频免费在线观看| 欧美日韩av久久| 欧美日韩综合久久久久久| 少妇猛男粗大的猛烈进出视频| 十八禁网站网址无遮挡| 久久久久人妻精品一区果冻| 国产亚洲av片在线观看秒播厂| 日产精品乱码卡一卡2卡三| 日韩av免费高清视频| 久久女婷五月综合色啪小说| 18禁裸乳无遮挡动漫免费视频| 久久久久久久大尺度免费视频| 午夜福利网站1000一区二区三区| 精品少妇内射三级| 欧美97在线视频| 国产国语露脸激情在线看| 久久精品夜色国产| 人体艺术视频欧美日本| 亚洲av电影在线观看一区二区三区| 丁香六月天网| 69精品国产乱码久久久| 日韩中文字幕欧美一区二区 | 欧美成人午夜精品| 成人亚洲欧美一区二区av| 校园人妻丝袜中文字幕| 久久精品aⅴ一区二区三区四区 | 欧美国产精品va在线观看不卡| av有码第一页| 免费av中文字幕在线| videos熟女内射| 一区在线观看完整版| 亚洲四区av| 看免费成人av毛片| 激情五月婷婷亚洲| 免费黄网站久久成人精品| 亚洲精品自拍成人| 国产精品熟女久久久久浪| 国产白丝娇喘喷水9色精品| av天堂久久9| 国产日韩欧美亚洲二区| 日本色播在线视频| www日本在线高清视频| 熟女av电影| 日本色播在线视频| 中文乱码字字幕精品一区二区三区| 日本黄色日本黄色录像| 亚洲精品国产一区二区精华液| 在线亚洲精品国产二区图片欧美| 满18在线观看网站| 久久韩国三级中文字幕| 午夜福利,免费看| 久久影院123| 亚洲欧美成人精品一区二区| 黄色怎么调成土黄色| 亚洲精品美女久久久久99蜜臀 | 午夜激情久久久久久久| 考比视频在线观看| 亚洲情色 制服丝袜| 99国产精品免费福利视频| 久久99热这里只频精品6学生| 久久鲁丝午夜福利片| 久久青草综合色| 精品酒店卫生间| 国产1区2区3区精品| 亚洲成色77777| 午夜老司机福利剧场| 免费观看在线日韩| 日韩 亚洲 欧美在线| 天天躁夜夜躁狠狠久久av| 少妇 在线观看| 免费高清在线观看视频在线观看| 亚洲国产欧美日韩在线播放| 日韩伦理黄色片| 亚洲伊人色综图| 亚洲一级一片aⅴ在线观看| 国产日韩欧美在线精品| 91久久精品国产一区二区三区| 熟女av电影| 天美传媒精品一区二区| 亚洲av电影在线进入| 亚洲精品国产av成人精品| av.在线天堂| 亚洲四区av| 成人毛片a级毛片在线播放| 侵犯人妻中文字幕一二三四区| 一区二区三区精品91| 丝袜在线中文字幕| 国产成人91sexporn| 久久久a久久爽久久v久久| av网站在线播放免费| 少妇被粗大的猛进出69影院| 久久久久久久久免费视频了| 成人免费观看视频高清| 人妻人人澡人人爽人人| 美女大奶头黄色视频| 一边摸一边做爽爽视频免费| 18在线观看网站| 国产无遮挡羞羞视频在线观看| 五月伊人婷婷丁香| 精品人妻一区二区三区麻豆| 中文精品一卡2卡3卡4更新| 久久婷婷青草| 综合色丁香网| 90打野战视频偷拍视频| 女性被躁到高潮视频| 日韩av不卡免费在线播放| 亚洲精品在线美女| 日韩不卡一区二区三区视频在线| 午夜影院在线不卡| 久久ye,这里只有精品| 亚洲在久久综合| 国产免费视频播放在线视频| 中文字幕另类日韩欧美亚洲嫩草| 下体分泌物呈黄色| 午夜av观看不卡| 老熟女久久久| 成人国语在线视频| 午夜免费男女啪啪视频观看| 亚洲,一卡二卡三卡| 99国产综合亚洲精品| 免费日韩欧美在线观看| 国产黄色免费在线视频| 狠狠精品人妻久久久久久综合| 久久久久久伊人网av| 精品国产一区二区三区四区第35| 国产人伦9x9x在线观看 | 欧美国产精品va在线观看不卡| 日韩人妻精品一区2区三区| 天堂俺去俺来也www色官网| 婷婷色综合大香蕉| 欧美精品高潮呻吟av久久| 男女高潮啪啪啪动态图| 中文字幕亚洲精品专区| av免费在线看不卡| 亚洲人成网站在线观看播放| 啦啦啦视频在线资源免费观看| 叶爱在线成人免费视频播放| 国产黄色视频一区二区在线观看| 街头女战士在线观看网站| 三上悠亚av全集在线观看| 亚洲婷婷狠狠爱综合网| 国产成人精品久久二区二区91 | 高清在线视频一区二区三区| 夫妻午夜视频| 欧美激情极品国产一区二区三区| 超碰成人久久| 在线观看免费日韩欧美大片| 伊人久久国产一区二区| 国产精品三级大全| 欧美激情极品国产一区二区三区| 日本欧美国产在线视频| 中文欧美无线码| 久久精品国产a三级三级三级| 精品久久蜜臀av无| 人妻一区二区av| 国产免费一区二区三区四区乱码| 边亲边吃奶的免费视频| 飞空精品影院首页| 久久久久国产一级毛片高清牌| 欧美变态另类bdsm刘玥| 一区二区三区四区激情视频| 久久人妻熟女aⅴ| 欧美成人午夜精品| 亚洲国产av影院在线观看| 国产色婷婷99| 校园人妻丝袜中文字幕| 久久人妻熟女aⅴ| 麻豆精品久久久久久蜜桃| 亚洲美女视频黄频| 久热这里只有精品99| 一区二区三区四区激情视频| 满18在线观看网站| 香蕉精品网在线| 成年女人毛片免费观看观看9 | 如何舔出高潮| 国产男人的电影天堂91| 亚洲在久久综合| 又黄又粗又硬又大视频| 在线观看免费日韩欧美大片| 91成人精品电影| 久久午夜综合久久蜜桃| 99精国产麻豆久久婷婷| 下体分泌物呈黄色| 日韩伦理黄色片| 美女中出高潮动态图| 日本av手机在线免费观看| 免费观看在线日韩| 亚洲av综合色区一区| 亚洲精品乱久久久久久| 国产伦理片在线播放av一区| 国产在视频线精品| 中文欧美无线码| 一个人免费看片子| 天堂8中文在线网| 一级毛片黄色毛片免费观看视频| 日韩一区二区三区影片| 国产97色在线日韩免费| 久久av网站| 午夜老司机福利剧场| 国产精品国产三级专区第一集| 伊人久久国产一区二区| 在线观看国产h片| 最近中文字幕2019免费版| 一级a爱视频在线免费观看| 精品国产乱码久久久久久小说| 男女无遮挡免费网站观看| av在线播放精品| 天天操日日干夜夜撸| 三上悠亚av全集在线观看| 熟女少妇亚洲综合色aaa.| 大香蕉久久网| 曰老女人黄片| 日韩不卡一区二区三区视频在线| 80岁老熟妇乱子伦牲交| 国产免费一区二区三区四区乱码| 国产极品天堂在线| 久久久a久久爽久久v久久| 久久综合国产亚洲精品| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久免费高清国产稀缺| 日韩三级伦理在线观看| 国产在线视频一区二区| 9热在线视频观看99| 亚洲欧洲国产日韩| 亚洲精华国产精华液的使用体验| 国语对白做爰xxxⅹ性视频网站| 亚洲精品,欧美精品| 天天躁狠狠躁夜夜躁狠狠躁| 色视频在线一区二区三区| 最近最新中文字幕大全免费视频 | 青春草国产在线视频| 亚洲一级一片aⅴ在线观看| 欧美人与善性xxx| 老司机影院成人| 久久精品久久久久久噜噜老黄| 秋霞伦理黄片| 午夜91福利影院| 大香蕉久久成人网| 丝袜喷水一区| av国产精品久久久久影院| 美女大奶头黄色视频| 亚洲精品一区蜜桃| 久久综合国产亚洲精品| 大话2 男鬼变身卡| 午夜精品国产一区二区电影| 男女国产视频网站| 你懂的网址亚洲精品在线观看| 999精品在线视频| 一级片'在线观看视频| 欧美激情极品国产一区二区三区| 亚洲欧美精品综合一区二区三区 | 午夜精品国产一区二区电影| www.熟女人妻精品国产| 国产精品av久久久久免费| 亚洲精品中文字幕在线视频| 免费久久久久久久精品成人欧美视频| 国产成人精品在线电影| 亚洲欧美中文字幕日韩二区| 久久久亚洲精品成人影院| 国产一级毛片在线| 久久久久久伊人网av| 欧美精品国产亚洲| 国产淫语在线视频| 人人澡人人妻人| 国产av精品麻豆| 九色亚洲精品在线播放| 街头女战士在线观看网站| 亚洲精品国产av蜜桃| 中文字幕色久视频| 国产毛片在线视频| 精品久久久精品久久久| 亚洲精品自拍成人| 亚洲伊人色综图| 国产人伦9x9x在线观看 | 国产精品三级大全| 不卡视频在线观看欧美| 91在线精品国自产拍蜜月| 久久精品国产亚洲av高清一级| 韩国高清视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 午夜福利乱码中文字幕| 黄色怎么调成土黄色| 丝袜美腿诱惑在线| 国产精品一二三区在线看| 性少妇av在线| 免费看不卡的av| 中文字幕av电影在线播放| 午夜福利网站1000一区二区三区| 久久久久久人人人人人|