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

    高光譜成像技術(shù)檢測油茶果成熟度

    2022-09-01 08:02:16胡逸磊姜洪喆周宏平許林云
    食品科學(xué) 2022年16期
    關(guān)鍵詞:成熟度油茶正確率

    胡逸磊,姜洪喆,周宏平,許林云,鞠 皓,王 影

    (南京林業(yè)大學(xué)機(jī)械電子工程學(xué)院,江蘇 南京 210037)

    油茶果采摘期受環(huán)境因素影響較大,不同地區(qū)的天氣、氣候、積溫、土壤肥力等條件不同,造成油茶果的成熟期有早有晚,因此油茶果的采摘時間應(yīng)根據(jù)其具體成熟度而定。油茶果采收前的最后一個月是油脂積累的高峰期,這段時間油茶果含油量增幅最為明顯,同時伴隨著內(nèi)部營養(yǎng)物質(zhì)的相互轉(zhuǎn)化。然而在此時期油茶果的外部形態(tài)特征如形狀大小、鮮果質(zhì)量、顏色等趨于穩(wěn)定,無明顯變化,從而給茶農(nóng)判斷油茶果的成熟度和最佳采摘期帶來一定的困難。故迫切需要提出一種快速準(zhǔn)確地檢測油茶果成熟度的方法,為油茶果的精準(zhǔn)采收作業(yè)提供幫助。

    高光譜成像技術(shù)是一種集圖像和光譜于一體的技術(shù),圖像中的每個像素點都包含特定位置的光譜信息,可實現(xiàn)被測物各組分分布情況的可視化。近年來不少研究人員利用高光譜成像技術(shù)在檢測果品的成熟度及相關(guān)成熟度參數(shù)方面取得了重大進(jìn)展。Pu Hongbin等比較了兩臺成像光譜儀(光譜集I:600~1 000 nm和光譜集II:1 000~2 500 nm)檢測荔枝成熟度的能力,建立的偏最小二乘判別分析(partial least squaresdiscrimination analysis,PLS-DA)在光譜集I上對成熟和未熟荔枝的分類正確率比在光譜集II上高6.25%。Zhang Chu等采集了成熟、半熟和未熟草莓在兩組不同波段下的高光譜圖像,融合草莓的光譜特征和紋理特征建立的支持向量機(jī)(support vector machine,SVM)模型發(fā)現(xiàn):在441.1~1 013.97 nm,SVM模型識別正確率為95%,在941.46~1 578.13 nm,SVM模型識別正確率只有70.83%。Rosli等將柑橘樣品劃分成6 個區(qū)域,分別采集每個區(qū)域的光譜反射率,發(fā)現(xiàn)盡管生長期和成熟期的柑橘果皮顏色相差不大,利用粒子群優(yōu)化(particle swarm optimization,PSO)算法建立基于光譜反射率的分類模型能篩選出未成熟的柑橘。周宏平等通過比較不同成像波段(400~1 000 nm及900~1 700 nm)得到油茶籽光譜反射率數(shù)據(jù)建立的偏最小二乘回歸(partial least squares regression,PLSR)模型性能的優(yōu)劣,得出900~1 700 nm是油茶籽含油率最佳檢測波段的結(jié)論。這些研究結(jié)果為使用光譜分析技術(shù)檢測油茶果成熟度奠定了良好的理論基礎(chǔ)。

    本研究通過高光譜成像技術(shù)獲取油茶果的高光譜數(shù)據(jù)并結(jié)合化學(xué)計量學(xué)方法獲取油茶果的成熟度參數(shù)作為參照,使用分類算法完成對油茶果成熟度的判別,并驗證分類模型的正確率,旨在為油茶果最佳采摘期的判斷提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 材料

    油茶果樣品來自江蘇南京市江寧區(qū)金航油茶專業(yè)合作社的油茶林,采集的油茶果均屬于霜降籽品種(在霜降節(jié)氣前后成熟),具體采樣要求為:采集位于近似同一海拔高度的茶樹上向陽面的油茶果,且油茶果的大小相似,以最大程度保證每一批次的樣品處于同一成熟度水平。

    于2020年10月13號、10月20號、10月27號、11月3號、11月10號分別采集油茶鮮果各200 個,當(dāng)天帶回進(jìn)行油茶果高光譜圖像的采集,并于11月14號采集果殼開裂的油茶果20 個作為對照組進(jìn)行含油率的測量(11月14號茶樹上幾乎沒有果皮完整的油茶果,由于過熟導(dǎo)致果殼開裂,整個茶園的油茶果采收工作基本完畢,因此特意留出5 棵茶樹用于實驗樣本的采集)。

    1.2 儀器與設(shè)備

    南京林業(yè)大學(xué)生物質(zhì)包裝無損檢測實驗室搭建的高光譜成像無損檢測平臺如圖1所示,包括成像光譜儀(GaiaField-V10E-AZ4型,400~1 000 nm)、探測器(sCMOS相機(jī))、白色食品級傳送帶(HSIACSD800)、一套由12 只50 W的鹵素?zé)艉吐瓷漶讽斀M成的照明系統(tǒng)以及一臺計算機(jī),其中成像光譜儀的光譜分辨率為5 nm,被測物品置于傳送帶上的載物臺由步進(jìn)電機(jī)驅(qū)動,暗箱用于屏蔽外界雜散光對數(shù)據(jù)采集的干擾;111-101v-10G型游標(biāo)卡尺 桂林廣陸數(shù)字測控股份有限公司;BSM-220.4型分析天平 上海卓精電子科技有限公司;DHG-9101-2SA型干燥箱 常州朗越儀器制造有限公司;NAI-ZFCDY-6Z型脂肪測定儀 上海那艾精密儀器有限公司。

    圖1 高光譜成像系統(tǒng)Fig. 1 Schematic of hyperspectral imaging system

    1.3 方法

    1.3.1 高光譜圖像的采集和校正

    高光譜圖像數(shù)據(jù)獲取基于計算機(jī)上的Specview軟件,將高光譜儀器預(yù)熱30 min后進(jìn)行油茶果圖像采集。為避免獲取的圖像失真,經(jīng)過多次預(yù)實驗確定最佳的數(shù)據(jù)采集參數(shù)如下:光譜儀I的曝光物距為300 mm,曝光時間為1.2 ms,電控位移臺掃描速率為0.601 4 nm/s,掃描線實際長度為200 mm,圖像分辨率為800像素×666像素。每次采集2 個油茶果,采集完成后統(tǒng)一裁剪成230像素×290像素的單個油茶果高光譜圖像。

    由于高光譜圖像采集過程中存在暗電流的影響,而且不同波段下成像系統(tǒng)光源的強(qiáng)度分布也不均勻,從而導(dǎo)致獲取的高光譜圖像中含有較大的噪聲。因此要對其進(jìn)行黑白校正以消除暗電流的影響,校正方法如式(1)所示:

    式中:為校正后的漫反射光譜圖像數(shù)據(jù);為樣本原始的漫反射光譜圖像數(shù)據(jù);為暗圖像數(shù)據(jù);為白板的漫反射圖像數(shù)據(jù)。

    1.3.2 理化指標(biāo)測定

    分別測定油茶果的果高、果徑、鮮果質(zhì)量、鮮籽質(zhì)量以及烘干后的油茶果殼質(zhì)量,得到油茶果殼的含水率和出籽率計算如式(2)、(3)所示:

    式中:為油茶果殼含水率/%;為出籽率/%;為果殼烘干后質(zhì)量/g;為油茶鮮果質(zhì)量/g;為油茶鮮籽質(zhì)量/g。

    油茶果的含油率按照GB 5009.6—2016《食品中脂肪的測定》的方法測定,并按照黃佳聰?shù)鹊姆椒ㄓ嬎愫吐?,如式?)所示:

    式中:為含油率/%;為出油量/g;為茶仁質(zhì)量/g。

    1.3.3 曲率校正方法

    油茶果表面曲率較大,使得油茶果高光譜圖像每個像素點的光譜強(qiáng)度存在較大差異,因此必須校正反射強(qiáng)度的變化。Gowen等發(fā)現(xiàn),均值歸一化比最大值法和中值法更適合校正球形水果表面的光譜強(qiáng)度變化。本研究使用均值歸一化方法對高光譜圖像的像素點光譜反射率進(jìn)行逐一校正,并利用變異系數(shù)評價光譜校正的質(zhì)量,計算如式(5)、(6)所示:

    式中:A為高光譜圖像中單個像素點的光譜向量;為光譜向量中元素個數(shù)(即波段數(shù)量);A為單個像素點校正后的光譜向量;SD為光譜矩陣的標(biāo)準(zhǔn)差;MN為光譜矩陣的平均值;CV為變異系數(shù)。

    1.3.4 變量選擇方法

    高光譜數(shù)據(jù)量大且數(shù)據(jù)之間的共線性嚴(yán)重,影響模型的運(yùn)算速度。因此采用以下3 種方法提取有效信息變量進(jìn)行對比,從而得到最優(yōu)的變量選擇方法。

    連續(xù)投影算法(succesive projections algorithm,SPA)是一種前向選擇算法,通過在光譜中尋找最低限度冗余光譜信息變量集,使得變量之間的共線性最小化。該方法需預(yù)先設(shè)置選擇的變量數(shù)范圍,最終選擇的變量數(shù)在該范圍內(nèi)的均方根誤差(root mean square error,RMSE)最小。本研究中選擇的最佳變量數(shù)范圍為5~30,選擇出28 個特征波長。

    競爭性自適應(yīng)重加權(quán)算法(competitive adaptive reweighted sampling,CARS)是一種以回歸系數(shù)(regression coefficient,RC)作為變量重要性評價指標(biāo)的變量選擇方法。該方法利用自適應(yīng)重加權(quán)采樣技術(shù)結(jié)合指數(shù)衰減函數(shù)優(yōu)選出PLSR模型中RC絕對值大的變量點,并將交叉驗證選出個PLSR子集模型中交叉驗證均方根誤差(root mean square error of cross-validation,RMSECV)最小的子集定義為最優(yōu)變量子集。本研究中將蒙特卡洛采樣次數(shù)設(shè)置為5 000,每次運(yùn)行程序選擇的訓(xùn)練集和測試集的樣本比例為2∶1,篩選出28 個特征波長。

    遺傳算法(genetic algorithm,GA)是模擬達(dá)爾文生物進(jìn)化論的自然選擇和遺傳學(xué)機(jī)理的生物進(jìn)化過程的計算模型,是一種通過模擬自然進(jìn)化過程搜索最優(yōu)解的方法。本研究設(shè)置進(jìn)化代數(shù)為200,算法運(yùn)行次數(shù)50,種群大小為64,初始時平均5 個波長構(gòu)成一個染色體,染色體個數(shù)為30,變異概率為1%,篩選出26 個特征波長。

    1.3.5 模型建立與評價

    采用SVM、-最近鄰(-nearest neighbor,KNN)算法、隨機(jī)森林(random forest,RF)算法、PLS-DA四種分類方法建立油茶果成熟度分類模型。其中SVM是一種非線性模型,能夠有效避免樣本空間的維數(shù)災(zāi)難;KNN屬于簡單的線性分類方法,易于實現(xiàn);RF是一種集成算法,抗過擬合能力強(qiáng);PLS-DA也是一種線性分類方法,魯棒性較好。最后通過圖像分類的方法,以整幅油茶果的高光譜圖像為輸入量,搭建基于遷移學(xué)習(xí)的卷積神經(jīng)網(wǎng)絡(luò)模型,并討論模型性能。

    KNN算法的核心思想是用距離最近的個樣本數(shù)據(jù)的類別代表目標(biāo)數(shù)據(jù)的類別,其中值是超參數(shù),會很大程度上影響分類的效果。將值設(shè)置為[3,40]以2為步長遞增,以平均正確率為依據(jù),使用10折交叉驗證法進(jìn)行值尋優(yōu)。

    RF算法是以Bagging為框架的基于CART決策樹的集成算法,該算法使用隨機(jī)抽樣形成的多個分類器達(dá)到更低的泛化誤差。對RF算法分類效果影響最大的超參數(shù)主要包括Bagging框架參數(shù)N estimators和CART決策樹參數(shù)Max depth,其中N estimators為學(xué)習(xí)器的最大迭代次數(shù),Max depth為決策樹的最大深度。本研究中將N estimators的值設(shè)置為[10,200]以10為步長,將Max depth的值設(shè)置為[1,15]以1為步長進(jìn)行網(wǎng)格搜索。

    SVM算法中使用高斯核函數(shù),其中參數(shù)和對模型的性能影響很大,是懲罰系數(shù),表示對誤差的容忍度,是高斯核函數(shù)中的一個參數(shù),反映數(shù)據(jù)映射到新的特征空間后的分布狀況。本研究將設(shè)置為以2為底,在[-10,10]內(nèi)以指數(shù)0.1為步長遞增,將設(shè)置為以2為底,在[-10,5]內(nèi)以指數(shù)0.375為步長遞增進(jìn)行網(wǎng)格搜索。

    PLS-DA是根據(jù)樣品的光譜矩陣和類別矩陣分解出的載荷矩陣和得分矩陣計算測試樣本的類別信息矩陣,然后根據(jù)類別信息矩陣與類別標(biāo)簽的接近程度確定樣品所屬類別。該模型中的潛變量數(shù)N components為超參數(shù),該值過大,容易發(fā)生過擬合,該值過小,則模型的精度不足。本研究將N components設(shè)置在[1,20]之間進(jìn)行網(wǎng)格搜索。

    本實驗采用分類正確率作為模型評價指標(biāo),計算如式(7)所示。正確率越高,說明模型的預(yù)測性能越好。使用的數(shù)據(jù)處理軟件包括ENVI5.1、Matlab2014a、Python3.8。

    式中:Acc為正確率/%;TP為正確識別為正樣本的數(shù)量;TN為正確識別為負(fù)樣本的數(shù)量;FP為錯誤識別為正樣本的數(shù)量;FN為錯誤識別為負(fù)樣本的數(shù)量。

    2 結(jié)果與分析

    2.1 油茶果成熟期理化參數(shù)

    從每一批次中取10 個樣品進(jìn)行油茶果果高、果徑、鮮果質(zhì)量、鮮籽質(zhì)量、出籽率、含油率、果殼含水率等理化參數(shù)的統(tǒng)計,取每一批次各理化參數(shù)的平均值作為該批次樣品的理化值,計算每個理化參數(shù)在整個采摘期內(nèi)的平均值和標(biāo)準(zhǔn)差,結(jié)果如表1所示。

    可以發(fā)現(xiàn)油茶果的果高、果徑在整個采摘期內(nèi)出現(xiàn)些許的遞增,但遞增幅度很小,因此標(biāo)準(zhǔn)差較低,只有0.43 mm左右;鮮果質(zhì)量范圍在25~30 g之間,鮮籽質(zhì)量范圍在10~15 g之間;油茶果的果殼含水率在采摘期內(nèi)有下降的趨勢,但下降的幅度較??;油茶果出籽率在35%~40%之間波動,油茶果含油率隨著采摘時間的推遲而遞增,增幅達(dá)到59%,其標(biāo)準(zhǔn)差為5.55%,比其他理化參數(shù)的標(biāo)準(zhǔn)差高很多,可見若提前采摘油茶果將造成茶油產(chǎn)量的嚴(yán)重?fù)p失。并且從采樣時間上看,11月10號油茶林中部分油茶果的果殼開裂,油茶果進(jìn)入全面采收階段,至11月14號采收工作基本完畢,預(yù)留的5 棵茶樹上的大部分油茶果的果殼開裂,有的油茶果已脫落在地,符合油茶果過熟期特征,將此時采集的油茶果作為對照組,發(fā)現(xiàn)其含油率基本保持穩(wěn)定,說明11月10號左右油茶果進(jìn)入完全成熟期,此時采摘可使茶油產(chǎn)量最大化,并且陸續(xù)有油茶果果殼開裂并脫落導(dǎo)致難以收集,無法再推遲油茶林的整體采收時間。

    表1 油茶果理化參數(shù)隨采摘時間的變化Table 1 Changes in physicochemical parameters of C. oleifera with picking date

    因此本研究采用含油率作為油茶果成熟度的衡量指標(biāo),并結(jié)合表2給出其成熟度的定量判別標(biāo)準(zhǔn):采收前30 d左右,且含油率為(22.00±1.00)%的樣品作為成熟度I;油茶果采收前23 d左右,且含油率為(24.00±1.00)%的樣品作為成熟度II;茶果采收前16 d左右,且含油率為(27.00±2.00)%的樣品作為成熟度III;茶果采收前9 d左右,且含油率為(31.50±2.50)%的樣品作為成熟度IV;油茶果采收期間,且含油率為(35.00±1.00)%的樣品作為成熟度V。成熟度等級I~V的油茶果樣品的成熟度依次遞增,成熟度V的油茶果為完熟期樣品。本實驗不包含過熟期樣品,因為過熟期的油茶果果殼開裂容易脫落,憑借人眼很容易判別,無需借助高光譜設(shè)備。將以上5 種成熟度油茶果樣品的高光譜數(shù)據(jù)作為數(shù)據(jù)集,用作成熟度的分類。

    表2 油茶果含油率的變化Table 2 Changes in oil content in C. oleifera with picking date

    2.2 基于像素點的油茶果高光譜圖像曲率校正

    油茶果果面曲率較大造成高光譜成像后每個像素點的光譜強(qiáng)度差異較大,并且從油茶果圖像的中心區(qū)域往果子邊緣處的光譜強(qiáng)度逐漸減小,導(dǎo)致邊緣處的光譜反射率比中心處低,直接使用這種狀態(tài)下數(shù)據(jù)進(jìn)行建模分類的誤差較大,現(xiàn)使用均值歸一化方法進(jìn)行校正。歸一化方法在光譜曲線預(yù)處理中有重要的應(yīng)用,可以在一定程度上減小水果形狀對采集光譜的影響。

    圖2 校正前(a)、后(b)油茶果半徑方向像素點光譜曲線Fig. 2 Pixel spectral curves of C. oleifera in the radius direction before (a) and after (b) correction

    圖3 774 nm波段校正前(a)、后(b)油茶果像素點反射率分布Fig. 3 Reflectivity distribution of the pixel points of C. oleifera before (a) and after (b) correction at 774 nm

    圖2a是校正前從油茶果中心沿半徑方向至邊緣處89 個像素點的光譜曲線圖,可以看出由于油茶果曲率的影響,不同位置的光譜反射強(qiáng)度差異很大,光譜CV為0.244 0;圖2b是相同位置的像素點校正后的光譜曲線圖,可以發(fā)現(xiàn)每個像素點之間的光譜強(qiáng)度差異變小很多,CV為0.075 3。

    774 nm通道油茶果樣品的反射率與背景的反射率差異最大,因此該通道能清晰地顯示油茶果樣品的灰度圖。圖3是在774 nm波段校正前后油茶果反射率的偽彩色分布,可以直觀地看到校正后油茶果像素點的反射率分布比校正前更均勻,說明均值歸一化方法的有效性。

    2.3 不同成熟度油茶果光譜特征分析

    在對1 000 份油茶果樣本的高光譜圖像進(jìn)行曲率校正之后,使用掩膜的方法計算油茶果的光譜值:以反射率0.2作為分割原始高光譜圖像中油茶果樣本與背景的閾值,將背景區(qū)域的像素值置為0,樣品區(qū)域的像素值置為1生成掩膜圖像,將掩膜圖像與油茶果圖像做內(nèi)積運(yùn)算得到去除背景信息的油茶果高光譜圖像,最后將分割出來的油茶果樣品每個波長通道所有像素點反射率的平均值作為該樣品的最終光譜值,從而提取出1 000 份樣品的平均光譜曲線。

    由于總樣本中可能存在奇異樣本影響最終分類效果,使用馬氏距離法分別對不同成熟度等級的樣本進(jìn)行檢測并予以剔除,最終將剩余939 個樣品的光譜數(shù)據(jù)用于建立分類模型。

    圖4 939 個樣本的原始(a)及其曲率校正后(b)的光譜曲線Fig. 4 Spectral curves of 939 samples before (a) and after (b)curvature correction

    939 個樣本曲率校正前后的光譜曲線圖見圖4,校正后的數(shù)據(jù)分布更加集中,數(shù)據(jù)特征更加明顯。670 nm處的吸收峰主要與油茶果殼中的花青素和葉綠素有關(guān),970 nm處的吸收峰與果殼中的水分有關(guān)。將5 個不同成熟度等級的油茶果光譜數(shù)據(jù)分別取平均值,得到油茶果的成熟度曲線如圖5所示??梢园l(fā)現(xiàn),不同成熟度油茶果的光譜強(qiáng)度差異主要體現(xiàn)在500~630 nm和720~970 nm之間。在500~630 nm之間,成熟度I和成熟度II樣品的反射強(qiáng)度差別較小,隨著采樣時間的推遲,成熟度III和成熟度IV樣品的反射強(qiáng)度較前者有所降低,成熟度V樣品的反射強(qiáng)度降至最低,該波段范圍內(nèi)呈現(xiàn)出油茶果的光譜反射率隨成熟度的增加而降低的規(guī)律。在720~970 nm之間,成熟度II~V樣品的反射強(qiáng)度隨成熟度的增加而降低,而從成熟度I~V樣品的反射強(qiáng)度經(jīng)歷兩次先增高后降低的過程,可能是與果殼中內(nèi)源激素含量的動態(tài)變化有關(guān)。在5 個不同的采收時間點,成熟度曲線存在較明顯的區(qū)別,為成熟度分類提供了理論依據(jù)。

    圖5 成熟度曲線Fig. 5 Spectral curves with different maturities

    2.4 基于全波長的油茶果成熟度分類

    將5 種不同成熟度等級的油茶果樣品賦予類別標(biāo)簽1~5,分別采用KNN、RF、SVM、PLS-DA建立基于全波長光譜信息的油茶果成熟度分類模型,各模型參數(shù)設(shè)置如表3所示。在對各個模型進(jìn)行參數(shù)尋優(yōu)后,獲得了各模型在最優(yōu)參數(shù)下的分類結(jié)果如表3和圖6所示??梢园l(fā)現(xiàn),SVM相比于其他3 個模型的分類效果最好:訓(xùn)練集的分類正確率為98%,測試集的平均分類正確率為97%,對成熟度I樣品的分類正確率最高;其次是RF模型,對油茶果成熟度的平均分類正確率為86%,對成熟度I和成熟度IV樣品的分類正確率較高;PLS-DA模型的分類效果一般,測試集的正確率為76%;KNN模型在訓(xùn)練集和測試集上的分類正確率是所有模型中最低。

    表3 不同分類模型的分類結(jié)果的比較Table 3 Comparison of results of classification of different classification models

    圖6 不同算法對5 種成熟度樣品的分類正確率Fig. 6 Classification accuracy of different algorithms for samples at five maturity stages

    2.5 基于特征波長的油茶果成熟度分類

    通過降維的方式選取特征波長建立分類模型,在最大程度保留有效波段信息的基礎(chǔ)上,通過減少特征維度提高計算速度。由于在4 種分類模型中,SVM分類器對油茶果光譜信息的分類效果最好,為便于比較不同降維方法的優(yōu)劣,這里統(tǒng)一使用SVM建立分類模型。

    CARS和GA選出的特征變量數(shù)目較多,而且CARS屬于快速粗選算法,因此本研究嘗試在CARS和GA降維的基礎(chǔ)上,再使用SPA算法進(jìn)行二次降維,進(jìn)一步消除干擾波段,最大程度保留有效波段。表4是不同降維方法選擇特征變量建立的SVM分類模型結(jié)果??梢园l(fā)現(xiàn),使用CARS降維后的28 個特征波長的光譜數(shù)據(jù)建立的SVM模型效果最好:訓(xùn)練集平均正確率為91%,測試集平均正確率為82%,對成熟度IV樣品的識別率最高為91%,對成熟度II樣品的識別率較低為73%。經(jīng)SPA和GA降維后建立的分類模型的效果相差不大:訓(xùn)練集上對油茶果的分類正確率為86%左右,測試集上的識別率為80%左右。主成分分析(principal component analysis,PCA)-SVM模型的效果最差:訓(xùn)練集上的分類正確率為64%,測試集上的分類正確率只有55%,可見PCA方法將全波段光譜數(shù)據(jù)轉(zhuǎn)化為7 個PC的過程中即使保留了99%的信息量,結(jié)合SVM建立的分類模型的性能還是較差。CARS-SPASVM和GA-SPA-SVM都是經(jīng)過二次降維之后建立的分類模型,兩組數(shù)據(jù)集下的SVM模型分類效果相當(dāng):測試集的正確率為75%左右,并且兩個模型均對成熟度V樣品的分類正確率最高,對成熟度II和成熟度III樣品的分類正確率較低,究其原因,成熟度V樣品的光譜強(qiáng)度相對于其他4 個成熟度等級的樣品最低,而成熟度II和成熟度III樣品的光譜強(qiáng)度相差不大,在某些波段甚至重疊,導(dǎo)致這兩類樣品的誤判數(shù)量較多??偟膩碚f,無論通過哪種方法尋找特征變量建立分類模型,在測試集上對成熟度I和成熟度V樣品的識別率最高,因為這兩類樣品的采樣時間間隔最長,光譜特征差異最大,而對成熟度II、成熟度III和成熟度IV樣品的分類正確率相對低一些,因為這3 類樣品采樣時間較為集中,光譜特征的差異性相對較小。

    表4 不同降維方法的SVM分類模型結(jié)果Table 4 Results of SVM classification models with different dimensionreduction methods

    2.6 基于圖像特征和光譜特征融合的油茶果成熟度分類

    采用顏色直方圖表征油茶果高光譜圖像的顏色特征。由于常用的顏色空間是RGB和HSV空間,分別取高光譜圖像的第64(641.8 nm、紅)、46(552.3 nm、綠)、27(459.8 nm、藍(lán))3 個通道生成對應(yīng)油茶果樣品的偽彩色RGB圖像。同時采用灰度共生矩陣表征油茶果高光譜圖像的紋理特征,為了能更直觀地以灰度共生矩陣描述紋理狀況,從矩陣導(dǎo)出能量、對比度、相關(guān)度、熵和逆差距5 個參數(shù)。

    為驗證油茶果高光譜圖像的顏色特征和紋理特征是否能作為成熟度分類的表征因子,以及在油茶果的光譜特征數(shù)據(jù)中融入圖像特征數(shù)據(jù)是否能起到優(yōu)化分類模型、提高模型正確率的效果,分別建立基于油茶果的顏色特征、紋理特征、融合圖像特征和光譜特征的SVM分類模型,同時約定光譜特征采用經(jīng)CARS降維后的數(shù)據(jù)集,因為各種變量選擇方法中,經(jīng)CARS選擇的特征波段建立的SVM模型效果最好。由于顏色特征中灰度級出現(xiàn)頻次的量級差異較大,紋理特征中各個紋理參數(shù)的量級也不同,因此在建模之前先對數(shù)據(jù)進(jìn)行歸一化處理。

    表5是不同輸入特征下的SVM模型分類結(jié)果,可以發(fā)現(xiàn),使用單一的顏色特征建立的油茶果成熟度SVM模型測試集正確率為79%,相比使用單一的光譜特征建模的正確率降低了3.7%,并且該模型對成熟度IV和成熟度V樣品的識別率較高,分別為92%和95%,對成熟度II和成熟度III樣品的識別率較低。使用單一的紋理特征建立的油茶果成熟度SVM分類模型正確率不足50%,說明油茶果的紋理特征參數(shù)在油茶果成熟期間沒有發(fā)生較明顯的變化,并且油茶果的表皮光滑,紋理特征并不明顯,不能較好地作為油茶果成熟度的表征因子。進(jìn)一步地,融合油茶果高光譜圖像中的顏色特征和紋理特征建立的SVM模型正確率沒有使用單一的顏色特征建立的模型正確率高,原因是數(shù)據(jù)集中的紋理參數(shù)拉低了模型性能。融合圖像中的光譜特征和顏色特征建立的SVM模型正確率達(dá)到93%,說明顏色特征對使用光譜特征建立的分類模型起到優(yōu)化增益作用,原因是在光譜數(shù)據(jù)中加入顏色數(shù)據(jù),使得不同油茶果樣品的特征得到更為充分詳盡的描述,也使得不同成熟度等級油茶果的特征得到差異化表達(dá),該模型對成熟度I、IV和V樣品的識別率較高,對成熟度II和成熟度III的油茶果樣品識別率稍低。融合圖像特征(顏色及紋理)和光譜特征建立的油茶果SVM模型的正確率為88%,比使用單一的光譜特征建立的模型正確率提高了7.3%,比使用光譜特征和顏色特征融合建立的SVM模型正確率降低了5.4%,說明油茶果的光譜信息與紋理信息相比占絕對主導(dǎo),加入紋理數(shù)據(jù)后反而會導(dǎo)致信息冗余,使得模型性能有所降低。

    表5 不同輸入特征下的SVM分類結(jié)果Table 5 Results of SVM classification with different input characteristics

    總體上說,融合顏色特征和光譜特征建立的油茶果成熟度SVM分類模型的效果最優(yōu),并且模型對成熟度I、IV、V油茶果樣品的識別率普遍高于對成熟度II、III樣品的識別率,原因是成熟度I、IV、V油茶果樣品之間的特征參數(shù)差異性較大,易于分辨,而成熟度II、III樣品之間的特征差異相對較小。

    3 結(jié) 論

    對利用高光譜成像技術(shù)進(jìn)行油茶果成熟度檢測的可行性進(jìn)行探究,采用實驗室內(nèi)靜態(tài)測量的方式觀察和分析油茶果的光學(xué)特性,在此基礎(chǔ)上分析利用油茶果的光學(xué)特征判斷油茶果的成熟度的準(zhǔn)確性。目前結(jié)合田間動態(tài)測量討論自然環(huán)境下光譜檢測油茶果成熟度的可行性的田間試驗正處于方案論證階段,將在今年油茶果成熟期(10—11月)進(jìn)行,為達(dá)到油茶先熟先采、后熟后采、隨熟隨采的智能化采摘新階段奠定理論基礎(chǔ)。主要結(jié)論如下:1)分別使用KNN、RF、SVM、PLS-DA建立基于全波段光譜數(shù)據(jù)的油茶果成熟度判別模型,發(fā)現(xiàn)SVM模型的分類正確率最高:訓(xùn)練集正確率為98%,測試集正確率為97%。2)分別使用SPA、CARS、GA、CARS-SPA、GA-SPA對全光譜進(jìn)行降維,建立基于特征變量的SVM模型,發(fā)現(xiàn)經(jīng)過CARS方法選擇的特征波段建立的分類模型正確率最高:訓(xùn)練集分類正確率為91%,測試集的正確率為82%。3)提取油茶果高光譜圖像中的顏色特征和紋理特征,分別建立SVM分類模型后發(fā)現(xiàn),融合顏色特征和光譜特征建立的SVM模型的正確率高于使用單一的光譜特征(經(jīng)CARS降維)建立的模型正確率:訓(xùn)練集分類正確率為95%,測試集正確率為93%。

    猜你喜歡
    成熟度油茶正確率
    產(chǎn)品制造成熟度在型號批生產(chǎn)風(fēng)險管理中的應(yīng)用
    門診分診服務(wù)態(tài)度與正確率對護(hù)患關(guān)系的影響
    整機(jī)產(chǎn)品成熟度模型研究與建立
    不同成熟度野生水茄果實的種子萌發(fā)差異研究
    種子(2018年9期)2018-10-15 03:14:04
    剛好夠吃6天的香蕉
    油茶芽苗嫁接育苗技術(shù)
    高產(chǎn)油茶的整形修剪技術(shù)
    永勝油茶的深情厚意
    中國西部(2017年4期)2017-04-26 03:49:43
    “四化”舉措做大做強(qiáng)油茶產(chǎn)業(yè)
    生意
    eeuss影院久久| 最近中文字幕高清免费大全6| 亚洲18禁久久av| 亚洲精品一区av在线观看| 欧美日韩综合久久久久久| 国产成人freesex在线 | 女同久久另类99精品国产91| 欧美人与善性xxx| 欧美高清成人免费视频www| 国产在线精品亚洲第一网站| 女的被弄到高潮叫床怎么办| 欧美丝袜亚洲另类| 亚洲欧美中文字幕日韩二区| 亚洲最大成人手机在线| 99久国产av精品国产电影| 在线免费十八禁| 男人的好看免费观看在线视频| 欧美xxxx黑人xx丫x性爽| 免费人成视频x8x8入口观看| 亚洲五月天丁香| 亚洲五月天丁香| 日韩高清综合在线| 免费观看的影片在线观看| 免费观看的影片在线观看| 男人的好看免费观看在线视频| 国产精华一区二区三区| 在现免费观看毛片| 免费观看的影片在线观看| 在现免费观看毛片| 国产精品嫩草影院av在线观看| 99riav亚洲国产免费| 日日干狠狠操夜夜爽| 国产一区二区三区在线臀色熟女| 一级毛片久久久久久久久女| av视频在线观看入口| 国产精品久久久久久av不卡| 亚洲精品久久国产高清桃花| 亚洲精品国产av成人精品 | 老熟妇乱子伦视频在线观看| 亚洲av二区三区四区| 国产视频一区二区在线看| 三级国产精品欧美在线观看| 亚洲激情五月婷婷啪啪| 国产一区二区三区av在线 | 久久精品国产亚洲av香蕉五月| 日日撸夜夜添| 亚洲无线在线观看| 色哟哟·www| 卡戴珊不雅视频在线播放| 51国产日韩欧美| 欧美高清性xxxxhd video| 美女被艹到高潮喷水动态| 在线播放无遮挡| 日韩成人伦理影院| 亚洲最大成人中文| 大香蕉久久网| 九九在线视频观看精品| 成人毛片a级毛片在线播放| av在线观看视频网站免费| 搡老岳熟女国产| 高清日韩中文字幕在线| 精品乱码久久久久久99久播| 国产一区二区激情短视频| 色哟哟·www| 一区二区三区四区激情视频 | 男女视频在线观看网站免费| 丝袜美腿在线中文| 中国美女看黄片| 国产精品一区www在线观看| 麻豆av噜噜一区二区三区| 成人亚洲精品av一区二区| 亚洲图色成人| 国产激情偷乱视频一区二区| 日本熟妇午夜| 高清午夜精品一区二区三区 | 无遮挡黄片免费观看| 99久久精品热视频| 黑人高潮一二区| 黄色视频,在线免费观看| 嫩草影院入口| 国产精品一区www在线观看| 久久草成人影院| 色尼玛亚洲综合影院| 亚洲在线观看片| 国产午夜精品论理片| 变态另类成人亚洲欧美熟女| 国内少妇人妻偷人精品xxx网站| 亚洲成人精品中文字幕电影| 国产精品av视频在线免费观看| 精品午夜福利在线看| 欧美日韩乱码在线| videossex国产| 熟女人妻精品中文字幕| 蜜桃亚洲精品一区二区三区| 又爽又黄a免费视频| 大又大粗又爽又黄少妇毛片口| 午夜免费男女啪啪视频观看 | 欧美性感艳星| 亚洲va在线va天堂va国产| 在线免费十八禁| 如何舔出高潮| 黄色欧美视频在线观看| 中文资源天堂在线| 色综合色国产| 亚洲中文字幕一区二区三区有码在线看| 最好的美女福利视频网| avwww免费| 亚洲精品色激情综合| 国国产精品蜜臀av免费| 欧美日韩乱码在线| 美女内射精品一级片tv| 日韩亚洲欧美综合| 如何舔出高潮| 成人欧美大片| 成人高潮视频无遮挡免费网站| 免费av不卡在线播放| 身体一侧抽搐| 午夜福利18| 丰满的人妻完整版| АⅤ资源中文在线天堂| or卡值多少钱| 午夜免费激情av| 欧美xxxx性猛交bbbb| 欧美高清性xxxxhd video| 日日摸夜夜添夜夜添av毛片| 亚洲中文字幕一区二区三区有码在线看| 国产高清视频在线播放一区| 亚洲av中文字字幕乱码综合| 久久久久久久久久成人| 国产真实伦视频高清在线观看| 中文字幕熟女人妻在线| 99热全是精品| 国产精品人妻久久久久久| 一区二区三区高清视频在线| 亚洲av免费在线观看| 久久久久九九精品影院| 亚洲av不卡在线观看| 人妻制服诱惑在线中文字幕| 亚洲成人精品中文字幕电影| 深爱激情五月婷婷| 亚洲欧美日韩高清专用| 99久国产av精品国产电影| 天堂影院成人在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产av在哪里看| 少妇裸体淫交视频免费看高清| 欧美色欧美亚洲另类二区| 99久久精品一区二区三区| 三级国产精品欧美在线观看| 九色成人免费人妻av| 乱系列少妇在线播放| 啦啦啦观看免费观看视频高清| 久久久精品大字幕| 国产熟女欧美一区二区| 免费电影在线观看免费观看| 高清毛片免费观看视频网站| 精品久久久久久久久久久久久| 在线观看66精品国产| av免费在线看不卡| 亚洲欧美清纯卡通| 精品国内亚洲2022精品成人| 91久久精品电影网| 欧美日韩在线观看h| 亚洲av成人av| 久久久精品欧美日韩精品| 国产成人a∨麻豆精品| 成人亚洲精品av一区二区| 国产伦在线观看视频一区| 日韩精品青青久久久久久| 久久人人爽人人爽人人片va| 乱系列少妇在线播放| 午夜视频国产福利| 波多野结衣高清无吗| 久久久精品94久久精品| 亚洲精品国产av成人精品 | 99视频精品全部免费 在线| 亚洲av中文av极速乱| 97超视频在线观看视频| 国产一区二区三区av在线 | 99久久精品一区二区三区| 国语自产精品视频在线第100页| 日韩人妻高清精品专区| 日日撸夜夜添| 嫩草影院新地址| a级毛色黄片| 大又大粗又爽又黄少妇毛片口| 老熟妇乱子伦视频在线观看| 神马国产精品三级电影在线观看| 十八禁网站免费在线| 亚洲欧美日韩卡通动漫| 1024手机看黄色片| 国产av一区在线观看免费| 国产精品乱码一区二三区的特点| 男女啪啪激烈高潮av片| 黑人高潮一二区| 国产成人a∨麻豆精品| 一级av片app| 亚洲国产日韩欧美精品在线观看| 麻豆久久精品国产亚洲av| 久久这里只有精品中国| 最近视频中文字幕2019在线8| 日本一二三区视频观看| 内地一区二区视频在线| 一个人看视频在线观看www免费| 日本免费一区二区三区高清不卡| 一区福利在线观看| 欧美日韩乱码在线| 亚洲人与动物交配视频| 天堂av国产一区二区熟女人妻| 人妻夜夜爽99麻豆av| 波多野结衣高清作品| 国产成人影院久久av| 内地一区二区视频在线| 99久久九九国产精品国产免费| 深夜a级毛片| 噜噜噜噜噜久久久久久91| 久99久视频精品免费| 99热这里只有是精品50| 亚洲性久久影院| 最近最新中文字幕大全电影3| 亚洲欧美日韩卡通动漫| 最好的美女福利视频网| 天堂影院成人在线观看| 国产欧美日韩一区二区精品| 搡老岳熟女国产| 国产伦精品一区二区三区视频9| 欧美性感艳星| 国产中年淑女户外野战色| 尾随美女入室| 成人毛片a级毛片在线播放| 黑人高潮一二区| 午夜精品在线福利| 成人美女网站在线观看视频| 国产视频一区二区在线看| 日日撸夜夜添| 亚洲欧美日韩卡通动漫| 欧美中文日本在线观看视频| 日日摸夜夜添夜夜添av毛片| 波多野结衣高清无吗| 亚洲精品国产成人久久av| 久久精品夜色国产| 国产亚洲精品综合一区在线观看| 亚洲一区二区三区色噜噜| 一进一出好大好爽视频| 69av精品久久久久久| 欧美另类亚洲清纯唯美| 不卡视频在线观看欧美| 欧美日韩在线观看h| 日日干狠狠操夜夜爽| 成人av在线播放网站| a级一级毛片免费在线观看| 国产高清不卡午夜福利| 久久久成人免费电影| 欧美成人a在线观看| 亚洲av不卡在线观看| 看十八女毛片水多多多| 成人国产麻豆网| 午夜福利在线观看吧| 亚洲国产精品国产精品| 成年版毛片免费区| eeuss影院久久| 校园人妻丝袜中文字幕| 国产日本99.免费观看| 91狼人影院| 亚洲国产精品sss在线观看| 嫩草影院新地址| 国产v大片淫在线免费观看| 亚洲欧美日韩东京热| 色噜噜av男人的天堂激情| 欧美成人a在线观看| 亚洲人成网站在线播| 精品福利观看| 淫妇啪啪啪对白视频| 亚洲av二区三区四区| 久久精品夜夜夜夜夜久久蜜豆| 观看美女的网站| 国产v大片淫在线免费观看| 白带黄色成豆腐渣| av免费在线看不卡| 99riav亚洲国产免费| 亚洲人成网站高清观看| 国产高潮美女av| 精品乱码久久久久久99久播| 久久午夜亚洲精品久久| 啦啦啦韩国在线观看视频| 一进一出抽搐gif免费好疼| 免费不卡的大黄色大毛片视频在线观看 | 97人妻精品一区二区三区麻豆| 观看美女的网站| 又爽又黄a免费视频| 国产午夜福利久久久久久| 日韩强制内射视频| 极品教师在线视频| 人人妻人人澡人人爽人人夜夜 | 国产一区二区三区av在线 | 高清毛片免费观看视频网站| 男女那种视频在线观看| 中文亚洲av片在线观看爽| 自拍偷自拍亚洲精品老妇| а√天堂www在线а√下载| 别揉我奶头 嗯啊视频| 一个人看的www免费观看视频| 日本 av在线| 最近手机中文字幕大全| 啦啦啦观看免费观看视频高清| 91午夜精品亚洲一区二区三区| 深夜精品福利| 非洲黑人性xxxx精品又粗又长| 国产伦在线观看视频一区| 黄片wwwwww| 一边摸一边抽搐一进一小说| 中文字幕久久专区| 人人妻,人人澡人人爽秒播| 观看美女的网站| 国产aⅴ精品一区二区三区波| 麻豆久久精品国产亚洲av| 女生性感内裤真人,穿戴方法视频| 热99re8久久精品国产| 成人永久免费在线观看视频| 久久精品久久久久久噜噜老黄 | 亚洲成人av在线免费| 成人国产麻豆网| 亚洲四区av| 亚洲人成网站在线观看播放| 亚洲精品一区av在线观看| 国产精品久久久久久av不卡| 真实男女啪啪啪动态图| 亚洲成人久久性| 男女之事视频高清在线观看| 欧美另类亚洲清纯唯美| 午夜福利在线观看免费完整高清在 | 免费观看在线日韩| 国产黄a三级三级三级人| 欧美色欧美亚洲另类二区| 亚洲,欧美,日韩| 国产精品av视频在线免费观看| 欧美在线一区亚洲| 国产精品一区二区免费欧美| 亚洲精品在线观看二区| 亚洲婷婷狠狠爱综合网| 大又大粗又爽又黄少妇毛片口| 国产精品无大码| 国产成人影院久久av| 看片在线看免费视频| 18禁黄网站禁片免费观看直播| 国产成人影院久久av| 亚洲国产日韩欧美精品在线观看| 最新中文字幕久久久久| 精品一区二区三区视频在线| 少妇的逼水好多| 美女大奶头视频| 禁无遮挡网站| 日韩在线高清观看一区二区三区| 亚洲精品日韩av片在线观看| 亚洲精品久久国产高清桃花| 国产伦精品一区二区三区四那| av国产免费在线观看| 99riav亚洲国产免费| 日日摸夜夜添夜夜添小说| 亚洲成人精品中文字幕电影| 国产成人a∨麻豆精品| 真人做人爱边吃奶动态| 婷婷色综合大香蕉| 男插女下体视频免费在线播放| 久久久欧美国产精品| 91麻豆精品激情在线观看国产| 欧美xxxx性猛交bbbb| 成人高潮视频无遮挡免费网站| 免费不卡的大黄色大毛片视频在线观看 | 男女之事视频高清在线观看| 十八禁国产超污无遮挡网站| 欧美日本亚洲视频在线播放| 久久午夜福利片| 亚洲成av人片在线播放无| 国产一区二区三区av在线 | 亚洲国产欧美人成| 特大巨黑吊av在线直播| 欧美日韩综合久久久久久| 色综合站精品国产| 亚洲国产精品久久男人天堂| 久久热精品热| 少妇猛男粗大的猛烈进出视频 | 国模一区二区三区四区视频| 99热精品在线国产| 亚洲精品日韩av片在线观看| 久久精品夜色国产| 夜夜看夜夜爽夜夜摸| 天美传媒精品一区二区| 直男gayav资源| 91午夜精品亚洲一区二区三区| 岛国在线免费视频观看| 内地一区二区视频在线| 男女啪啪激烈高潮av片| 18禁在线无遮挡免费观看视频 | 亚洲熟妇熟女久久| 国产单亲对白刺激| 亚洲欧美成人精品一区二区| 在线观看av片永久免费下载| 一个人观看的视频www高清免费观看| 99热这里只有是精品50| 国产精品久久久久久亚洲av鲁大| 日韩av在线大香蕉| 国产探花在线观看一区二区| 波多野结衣高清无吗| 国模一区二区三区四区视频| 久久久欧美国产精品| 久久久久久久久大av| 一进一出抽搐动态| 大香蕉久久网| 九九热线精品视视频播放| 在线播放国产精品三级| 国产三级中文精品| 少妇猛男粗大的猛烈进出视频 | 最近手机中文字幕大全| 成人永久免费在线观看视频| 深夜精品福利| 成熟少妇高潮喷水视频| 久久欧美精品欧美久久欧美| 高清毛片免费观看视频网站| 特级一级黄色大片| 久久人人爽人人爽人人片va| 国产精品永久免费网站| 97碰自拍视频| 午夜福利在线观看吧| 亚洲成av人片在线播放无| 黑人高潮一二区| 精品99又大又爽又粗少妇毛片| 久久欧美精品欧美久久欧美| 亚洲精品久久国产高清桃花| 欧美+日韩+精品| 国产欧美日韩精品亚洲av| 亚洲不卡免费看| av在线观看视频网站免费| 在现免费观看毛片| 日韩精品中文字幕看吧| 欧美性猛交黑人性爽| 精品久久久噜噜| 亚洲精品一区av在线观看| 最近在线观看免费完整版| 免费无遮挡裸体视频| 久久久久久久午夜电影| 国产色爽女视频免费观看| 精华霜和精华液先用哪个| 在线免费十八禁| 精品熟女少妇av免费看| 69人妻影院| 人人妻人人看人人澡| 国产男人的电影天堂91| 欧美另类亚洲清纯唯美| 日产精品乱码卡一卡2卡三| 午夜精品在线福利| 人人妻人人澡人人爽人人夜夜 | 最近中文字幕高清免费大全6| 最近手机中文字幕大全| 搡老熟女国产l中国老女人| 日本精品一区二区三区蜜桃| 欧美极品一区二区三区四区| 一进一出抽搐gif免费好疼| 性色avwww在线观看| 国产成人freesex在线 | 亚洲人成网站高清观看| 日本色播在线视频| 亚洲精品国产av成人精品 | 久久6这里有精品| 亚洲精华国产精华液的使用体验 | 国国产精品蜜臀av免费| 尤物成人国产欧美一区二区三区| 国产午夜精品久久久久久一区二区三区 | 成人漫画全彩无遮挡| 美女被艹到高潮喷水动态| 国产精品人妻久久久影院| 99久久成人亚洲精品观看| 午夜福利18| 热99re8久久精品国产| 亚洲一区二区三区色噜噜| 男插女下体视频免费在线播放| 少妇高潮的动态图| 免费大片18禁| 国产精品电影一区二区三区| 国产精品av视频在线免费观看| 成年免费大片在线观看| 人妻夜夜爽99麻豆av| 高清毛片免费观看视频网站| 又黄又爽又免费观看的视频| 乱人视频在线观看| 日韩一区二区视频免费看| 在线免费观看的www视频| 性欧美人与动物交配| 99国产极品粉嫩在线观看| 一进一出抽搐动态| 悠悠久久av| 美女cb高潮喷水在线观看| 国产精品国产高清国产av| 免费观看人在逋| 亚洲欧美成人精品一区二区| 日韩三级伦理在线观看| 老师上课跳d突然被开到最大视频| 精品国产三级普通话版| 蜜桃亚洲精品一区二区三区| 午夜福利在线在线| 久久久国产成人免费| 欧美日韩一区二区视频在线观看视频在线 | 岛国在线免费视频观看| 波多野结衣巨乳人妻| 亚洲熟妇熟女久久| 最近最新中文字幕大全电影3| 久久欧美精品欧美久久欧美| 国产淫片久久久久久久久| 99视频精品全部免费 在线| 性插视频无遮挡在线免费观看| 不卡视频在线观看欧美| 美女xxoo啪啪120秒动态图| 国产欧美日韩精品一区二区| 深夜精品福利| 久久久a久久爽久久v久久| 久久中文看片网| 久久人妻av系列| 亚洲自偷自拍三级| 色尼玛亚洲综合影院| 日韩精品有码人妻一区| 一卡2卡三卡四卡精品乱码亚洲| 91av网一区二区| 欧美丝袜亚洲另类| 久久久久久九九精品二区国产| 在现免费观看毛片| 国产一级毛片七仙女欲春2| 免费观看的影片在线观看| 亚洲中文字幕一区二区三区有码在线看| 久久韩国三级中文字幕| 三级经典国产精品| 亚洲婷婷狠狠爱综合网| 人人妻人人澡人人爽人人夜夜 | 久久久久九九精品影院| 国产精品人妻久久久久久| 亚洲精品成人久久久久久| 干丝袜人妻中文字幕| 中文字幕av成人在线电影| 亚洲欧美日韩卡通动漫| 91狼人影院| 我要搜黄色片| 伦精品一区二区三区| 国产精品女同一区二区软件| 亚洲国产日韩欧美精品在线观看| 听说在线观看完整版免费高清| 午夜精品一区二区三区免费看| 天堂动漫精品| 精品久久久噜噜| 少妇的逼好多水| 午夜精品一区二区三区免费看| 99热只有精品国产| 国产综合懂色| 长腿黑丝高跟| 啦啦啦观看免费观看视频高清| 十八禁网站免费在线| 亚洲在线自拍视频| 国产成人福利小说| a级一级毛片免费在线观看| a级毛片a级免费在线| 中文资源天堂在线| 在线观看美女被高潮喷水网站| 国产伦在线观看视频一区| 亚洲七黄色美女视频| 黄色配什么色好看| 卡戴珊不雅视频在线播放| 女同久久另类99精品国产91| 春色校园在线视频观看| 成人毛片a级毛片在线播放| 国产中年淑女户外野战色| 蜜臀久久99精品久久宅男| 日本黄色片子视频| 国产三级在线视频| 国产精品久久久久久久久免| 亚洲欧美日韩无卡精品| 如何舔出高潮| 午夜a级毛片| 直男gayav资源| 国产精品,欧美在线| 亚洲在线自拍视频| 久久久久久国产a免费观看| 中文资源天堂在线| 99国产精品一区二区蜜桃av| 免费av不卡在线播放| 免费观看的影片在线观看| 九九在线视频观看精品| 男人和女人高潮做爰伦理| 91午夜精品亚洲一区二区三区| 97超级碰碰碰精品色视频在线观看| 国产真实伦视频高清在线观看| 中文字幕免费在线视频6| 91麻豆精品激情在线观看国产| 啦啦啦观看免费观看视频高清| a级毛片免费高清观看在线播放| 18+在线观看网站| 天美传媒精品一区二区| 网址你懂的国产日韩在线| 亚洲精品国产av成人精品 | 亚洲欧美清纯卡通| 国内少妇人妻偷人精品xxx网站| 97超级碰碰碰精品色视频在线观看| 美女大奶头视频| 高清毛片免费看| 免费高清视频大片| 91狼人影院| 亚洲高清免费不卡视频| 国产亚洲精品久久久久久毛片| 午夜精品在线福利| av中文乱码字幕在线| 特大巨黑吊av在线直播| 国产一区二区在线观看日韩| 中文字幕人妻熟人妻熟丝袜美| 国产av不卡久久| 亚洲性久久影院| 国产精品久久久久久av不卡| 国产午夜福利久久久久久| 三级毛片av免费|