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

    高光譜圖像技術(shù)檢測(cè)梨表面農(nóng)藥殘留試驗(yàn)研究

    2011-12-06 02:52:02索少增劉翠玲吳靜珠陳興海孫曉榮吳勝男
    關(guān)鍵詞:毒死殘留量光譜

    索少增, 劉翠玲, 吳靜珠, 陳興海, 孫曉榮, 吳勝男

    (1.北京工商大學(xué)計(jì)算機(jī)與信息工程學(xué)院,北京 100048;2.北京卓立漢光儀器有限公司,北京 101149)

    高光譜圖像技術(shù)檢測(cè)梨表面農(nóng)藥殘留試驗(yàn)研究

    索少增1, 劉翠玲1, 吳靜珠1, 陳興海2, 孫曉榮1, 吳勝男1

    (1.北京工商大學(xué)計(jì)算機(jī)與信息工程學(xué)院,北京 100048;2.北京卓立漢光儀器有限公司,北京 101149)

    以滴有不同濃度毒死蜱和炔螨特農(nóng)藥的水晶皇冠梨為研究對(duì)象,探討高光譜圖像技術(shù)結(jié)合人工神經(jīng)網(wǎng)絡(luò)方法檢測(cè)水果表面農(nóng)藥殘留量的可行性.分別配制不同濃度分布的毒死蜱水溶液和炔螨特水溶液樣本各20個(gè),按100 μL和150 μL取農(nóng)藥溶液滴在梨表面,在835.467 8~1 648.356 8 nm范圍采集高光譜圖像,提取感興趣區(qū)域數(shù)據(jù)得到4組樣本數(shù)據(jù),每組20個(gè),每組隨機(jī)抽取5個(gè)樣本作為預(yù)測(cè)集,利用BP人工神經(jīng)網(wǎng)絡(luò)對(duì)每組數(shù)據(jù)分別建立數(shù)學(xué)模型.滴有150 μL和100 μL農(nóng)藥溶液區(qū)域的建模結(jié)果為:兩種農(nóng)藥的殘留樣本相關(guān)系數(shù)分別都大于0.99和0.95;RMSEC和RMSEP的最大值分別為0.634 9,1.323 9和1.742 5,3.441 7.結(jié)果表明:150 μL農(nóng)藥樣本區(qū)優(yōu)于100 μL農(nóng)藥樣本區(qū)的建模結(jié)果,高光譜圖像技術(shù)結(jié)合人工神經(jīng)網(wǎng)絡(luò)法檢測(cè)梨表面農(nóng)藥殘留量是可行的,為水果表面農(nóng)藥殘留量檢測(cè)提供了新方法.

    高光譜圖像;BP神經(jīng)網(wǎng)絡(luò);農(nóng)藥殘留;毒死蜱;炔螨特

    自1980年A.F.H.Goetz等人在對(duì)機(jī)載成像光譜儀改進(jìn)的遙感技術(shù)重要改革中提出高光譜圖像的概念以來,高光譜圖像技術(shù)已經(jīng)發(fā)展為覆蓋了上百條光譜通道、像素點(diǎn)攜帶波譜信息量豐富的高分辨率檢測(cè)技術(shù)[1].高光譜圖像系統(tǒng)包括數(shù)字?jǐn)z像機(jī)和光譜儀,圖像技術(shù)和光譜技術(shù)的結(jié)合使得高光譜圖像技術(shù)可以同時(shí)分析目標(biāo)樣本的物理和化學(xué)特征,因此,在食品工業(yè)中經(jīng)常被用來進(jìn)行樣本的識(shí)別、分類和瑕疵檢測(cè)[2].

    近年來,高光譜圖像技術(shù)被廣泛地用于對(duì)精度要求高的農(nóng)產(chǎn)品品質(zhì)檢測(cè)和食品安全方面檢測(cè)的實(shí)驗(yàn)室研究.洪添勝等[3]探討了高光譜圖像技術(shù)在雪花梨品質(zhì)無損檢測(cè)應(yīng)用中的可行性,得到了肯定的結(jié)果.陳全勝等[4]利用高光譜圖像技術(shù)對(duì)茶葉質(zhì)量等級(jí)進(jìn)行判別,使用反向傳播神經(jīng)網(wǎng)絡(luò)建立的模型獲得了較高的識(shí)別率.薛利紅等[5]分別使用PLS和PCR方法建立了可見近紅外高光譜圖像技術(shù)檢測(cè)菠菜葉片硝酸鹽含量的模型,通過相關(guān)系數(shù)和均方根誤差分析證明了此方法的可行性.鄒偉等[6]使用高光譜圖像技術(shù)進(jìn)行油菜籽品種鑒別,得到了分類結(jié)果.薛龍等[7]針對(duì)水果表面農(nóng)藥殘留,以滴有較高濃度的臍橙為研究對(duì)象,利用光譜范圍425~725 nm的高光譜圖像系統(tǒng)進(jìn)行檢測(cè),發(fā)現(xiàn)對(duì)較高濃度的農(nóng)藥殘留檢測(cè)效果較好.

    本文使用短波近紅外高光譜法結(jié)合BP人工神經(jīng)網(wǎng)絡(luò)檢測(cè)水晶皇冠梨表面微量農(nóng)藥殘留量,旨在探討高光譜圖像技術(shù)檢測(cè)微量農(nóng)藥殘留量的可行性,并為發(fā)現(xiàn)微量農(nóng)藥殘留量無損、快速、無接觸和無污染的安全檢測(cè)新方法[8-9]作出貢獻(xiàn).

    1 高光譜圖像技術(shù)

    高光譜圖像是一種新興的集數(shù)字?jǐn)z像機(jī)和光譜攝制儀于一體的無接觸式檢測(cè)技術(shù),高光譜圖像系統(tǒng)采集的數(shù)據(jù)是由圖像平面坐標(biāo)軸和波譜空間坐標(biāo)軸組成的三維的數(shù)據(jù)立方體空間,其中,圖像坐標(biāo)中的每個(gè)像素點(diǎn)都由數(shù)百張甚至更多不同波段信息的灰度圖像堆疊在一起而形成完整的波譜信息[10].

    高光譜圖像系統(tǒng)用二維陣列傳感器以一條窄帶掃描傳送帶上樣本平面圖像的同時(shí)獲得波譜信息,從而構(gòu)成三維數(shù)據(jù)空間體[11].其中x軸和y軸確定平面可視圖像的坐標(biāo),縱軸λ為波長信息的坐標(biāo)軸.兼圖像處理和光譜分析的優(yōu)點(diǎn),利用每個(gè)像素點(diǎn)都包含系統(tǒng)所覆蓋整個(gè)波段光譜數(shù)據(jù)的原理,通過在灰度圖像上選擇任何感興趣區(qū)域(region of interest,ROI)可提取波譜信息分析樣本,高光譜圖像原理如圖1.克服了傳統(tǒng)NIR檢測(cè)方法不能精確檢測(cè)到樣本ROI或?qū)颖緬呙璨痪鶆虻娜秉c(diǎn).

    圖1 高光譜圖像原理Fig.1 Principle of hyperspectral imaging

    2 材料與方法

    2.1 試驗(yàn)試劑

    農(nóng)藥炔螨特標(biāo)準(zhǔn)物質(zhì)(標(biāo)準(zhǔn)值1.00 mg/mL,以乙腈為溶劑)購自中國計(jì)量科學(xué)研究院,根據(jù)中華人民共和國農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn):NY1500.41.3~1500.41.6—2009,NY1500.50 ~ 1500.92—2009 規(guī)定,炔螨特最大殘留量為5 mg/kg.使用標(biāo)準(zhǔn)物質(zhì)在濃度范圍2~20 mg/kg配制20個(gè)樣本.

    農(nóng)藥毒死蜱標(biāo)準(zhǔn)物質(zhì)(標(biāo)準(zhǔn)值1.01 mg/mL,以甲醇為溶劑)同樣購自中國計(jì)量科學(xué)研究院,根據(jù)中華人民共和國農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn):NY1500.41.3~1500.41.6—2009,NY1500.50 ~ 1500.92—2009 規(guī)定,毒死蜱最大殘留量為1 mg/kg.使用標(biāo)準(zhǔn)物質(zhì)在濃度范圍0.03~6 mg/kg配制20個(gè)樣本.

    分別用移液槍取每種濃度的農(nóng)藥樣本150 μL和100 μL滴在購買的無公害水晶皇冠梨表面,放置8小時(shí)后進(jìn)行數(shù)據(jù)采集.在每個(gè)水晶皇冠梨樣本表面都滴有150 μL和100 μL的農(nóng)業(yè)標(biāo)準(zhǔn)規(guī)定的最大殘留量農(nóng)藥樣本與其他濃度農(nóng)藥樣品進(jìn)行比對(duì).

    2.2 儀器及參數(shù)設(shè)置

    高光譜圖像數(shù)據(jù)采集使用北京卓立漢光儀器有限公司HyperSIS-NIR近紅外增強(qiáng)型高光譜成像系統(tǒng).該系統(tǒng)主要由高光譜成像儀、CCD相機(jī)、光源、暗箱、計(jì)算機(jī)組成,如圖 2.光譜掃描范圍:835.467 8~1 648.356 8 nm,光譜分辨率為6.3 nm,掃描速度:60 images/s,焦平面陣列為320×256,曝光時(shí)間0.03 s,傳送帶的傳輸速度為:2.35 mm/s.

    圖2 北京卓立漢光儀器有限公司HyperSIS-NIR高光譜成像系統(tǒng)Fig.2 HyperSIS-NIR hyperspectral imaging system of Zolix Instruments CO.LTD

    2.3 高光譜圖像采集

    高光譜圖像數(shù)據(jù)的采集使用spectraSENS高光譜軟件.高光譜圖像數(shù)據(jù)分析軟件使用ENVI(environment for visualizing images)4.7,MATLAB 7.7 和OPUS 6.5.

    在采集高光譜圖像數(shù)據(jù)之前,為了克服光強(qiáng)分布弱的波段存在的圖像噪聲和暗電流的影響,首先掃描標(biāo)準(zhǔn)白板采集反射率為100%的全白標(biāo)定圖像Dwhite,然后蓋上攝像頭蓋采集反射率為0的全黑標(biāo)定圖像Ddark,進(jìn)行黑白標(biāo)定后再進(jìn)行樣本光譜圖像Dsample的采集,由此可根據(jù)公式1得到黑白標(biāo)定后的相對(duì)樣本光譜圖像的感興趣像素區(qū)域c或波段i處的反射率R:

    3 結(jié)果與分析

    3.1 感興趣區(qū)域波譜信息的提取

    采集高光譜圖像數(shù)據(jù)之后,根據(jù)特征波段灰度圖中明顯可見的滴有農(nóng)藥樣本的區(qū)域,以面積為80個(gè)像素的橢圓形提取感興趣區(qū)域,如圖3.另外,在沒有滴農(nóng)藥的皇冠梨表面區(qū)域,可以用同樣的方法提取農(nóng)藥含量為0 mg/kg的光譜曲線.

    圖3 感興趣區(qū)域的光譜信息提取Fig.3 Drawing spectrum data from RIO of hyperspectral imaging

    3.2 建立模型

    提取到感興趣區(qū)域的波譜數(shù)據(jù)后,對(duì)數(shù)據(jù)進(jìn)行21點(diǎn)平滑和矢量歸一化處理,處理后的譜圖如圖4.根據(jù)OPUS6.5軟件分析得出的波段貢獻(xiàn)率,毒死蜱農(nóng)殘樣本特征波段選擇為1 323.201 2~1 406.084 0.炔螨特農(nóng)殘樣本特征波段選取1 313.637 8~1 453.901 0和1 479.403 4~1 552.722 8兩個(gè)波段.

    圖4 預(yù)處理后4組樣本的波譜圖Fig.4 Spectra of four sample groups after Pre-Processing

    3.2.1 毒死蜱樣本的BP神經(jīng)網(wǎng)絡(luò)模型

    把每組校正集中的16個(gè)樣本做訓(xùn)練樣本,使用MATLAB7.7神經(jīng)網(wǎng)絡(luò)工具箱中Powell-Beale共軛梯度反向傳播算法訓(xùn)練函數(shù)(traincgb)分別為150 μL和100 μL兩組毒死蜱樣本的ROI波譜信息建立3層BP神經(jīng)網(wǎng)絡(luò)模型,以神經(jīng)網(wǎng)絡(luò)輸出和實(shí)際化學(xué)值的平均絕對(duì)誤差(mean squared error,MSE)為其性能函數(shù),經(jīng)多次調(diào)試,根據(jù)最優(yōu)結(jié)果將隱含層選為16個(gè)神經(jīng)元,將訓(xùn)練目標(biāo)誤差設(shè)為0.01,訓(xùn)練1 000次,訓(xùn)練好網(wǎng)絡(luò)后對(duì)每組的5個(gè)預(yù)測(cè)樣本進(jìn)行預(yù)測(cè),模型的運(yùn)行結(jié)果如下.150 μL組:R=0.992 4,RMSEC=0.268 3,RMSEP=1.742 5.100 μL 組:R=0.953 0,RMSEC=0.652 9,RMSEP=3.441 7.模型的擬合結(jié)果如圖5和圖6,預(yù)測(cè)結(jié)果如表1和表2.

    圖5 毒死蜱150 μL感興趣區(qū)域組樣本擬合結(jié)果Fig.5 Correlation plots between prediction and actual values of 150 μL ROI for Chlorpyrifos solution

    圖6 毒死蜱100 μL感興趣區(qū)域組樣本擬合結(jié)果Fig.6 Correlation plots between prediction and actual values of 100 μL ROI for Chlorpyrifos solution

    表1 毒死蜱150 μL感興趣區(qū)域BP神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)結(jié)果Tab.1 Predicting result of 150 μL ROI for Chlorpyrifossolution samples

    表2 毒死蜱100 μL感興趣區(qū)域BP神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)結(jié)果Tab.2 Predicting result of 100 μL ROI for Chlorpyrifos solution samples

    3.2.2 炔螨特樣本的BP神經(jīng)網(wǎng)絡(luò)模型

    用兩組校正集中的16個(gè)樣本做訓(xùn)練樣本,利用MATLAB7.7神經(jīng)網(wǎng)絡(luò)工具箱中量化共軛梯度(scaled conjugate gradient,SCG)反向傳播算法訓(xùn)練函數(shù)(trainscg訓(xùn)練函數(shù))分別為 150 μL 和 100 μL兩組炔螨特樣本的ROI波譜信息建立三層BP神經(jīng)網(wǎng)絡(luò)模型,以神經(jīng)網(wǎng)絡(luò)輸出和實(shí)際化學(xué)值的平均絕對(duì)誤差(mean squared error,MSE)為性能函數(shù),經(jīng)多次調(diào)試,根據(jù)最優(yōu)結(jié)果將隱含層選為16個(gè)神經(jīng)元,將訓(xùn)練目標(biāo)誤差設(shè)為0.01,訓(xùn)練1 000次,訓(xùn)練好網(wǎng)絡(luò)后對(duì)每組的5個(gè)預(yù)測(cè)樣本進(jìn)行預(yù)測(cè),模型的運(yùn)行結(jié)果如下.150 μL 組:R=0.993 5,RMSEC=0.634 9,RMSEP=1.593 2.100 μL 組:R=0.972 4,RMSEC=1.323 9,RMSEP=2.722 7.模型的擬合結(jié)果如圖7和圖8,預(yù)測(cè)結(jié)果如表3和表4.

    圖7 炔螨特150 μL感興趣區(qū)域組樣本擬合結(jié)果Fig.7 Correlation plots between prediction and actual values of 150 μL ROI for Propargite solution

    表3 炔螨特150 μL感興趣區(qū)域BP神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)結(jié)果Tab.3 Predicting result of 150 μL ROI forPropargite solution samples

    圖8 炔螨特100 μL感興趣區(qū)域組樣本擬合結(jié)果Fig.8 Correlation plots between prediction and actual values of 100 μL ROI for Propargite solution

    表4 炔螨特100 μL感興趣區(qū)域BP神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)結(jié)果Tab.4 Predicting result of 100 μL ROI for Propargite solution samples

    4 結(jié)論

    農(nóng)藥殘留檢測(cè)問題是食品安全領(lǐng)域的重要內(nèi)容.由于農(nóng)產(chǎn)品農(nóng)藥殘留為微量檢測(cè),對(duì)儀器精度和檢測(cè)方法的選擇有較高的要求.本文對(duì)高光譜圖像技術(shù)檢測(cè)水晶皇冠梨表面農(nóng)藥殘留量進(jìn)行了研究,充分利用人工神經(jīng)網(wǎng)絡(luò)強(qiáng)大的非線性映射和泛化能力建立復(fù)雜定量分析模型,對(duì)分析結(jié)果總結(jié)如下.

    1)特征波長的選取在通過主成分分析按貢獻(xiàn)確定的波段中仍需從多波段中篩選擇農(nóng)藥物質(zhì)的特征波段,選擇過程比較繁瑣,這就需要建立所有農(nóng)藥品種的譜庫,確定特征波段,由于農(nóng)藥品種極多,此任務(wù)需耗費(fèi)大量的時(shí)間和經(jīng)費(fèi).

    2)根據(jù)建模結(jié)果,尤其對(duì)預(yù)測(cè)集的預(yù)測(cè)結(jié)果(RMSEP)來看,滴有較多農(nóng)藥溶液的水晶皇冠梨的建模效果較好,此方法所建立模型的適用性仍然在一定程度上受到農(nóng)藥殘留量的制約,因此,如何優(yōu)化模型提高預(yù)測(cè)能力仍需要進(jìn)一步研究.

    3)在使用BP神經(jīng)網(wǎng)絡(luò)建模過程中,確定隱含層神經(jīng)元的數(shù)目是難點(diǎn)之一,隱含層神經(jīng)元數(shù)目太多,會(huì)出現(xiàn)過擬合,導(dǎo)致泛化能力很弱;隱含層數(shù)目太少,擬合不充分,模型質(zhì)量差.因此,需要多次調(diào)試確定隱含層神經(jīng)元數(shù).

    結(jié)果表明,高光譜圖像技術(shù)結(jié)合BP人工神經(jīng)網(wǎng)絡(luò)可用于農(nóng)藥殘留量的檢測(cè),無需對(duì)樣本進(jìn)行任何預(yù)處理即可進(jìn)行數(shù)據(jù)采集,高光譜圖像技術(shù)是農(nóng)藥殘留量無損、高精度和快速的檢測(cè)方法.

    [1]Antonio P,Jon A B,Joseph W B,et al.Recent advances in techniques for hyperspectral image processing[J].Remote Sensing of Environment,2009(113):110-122.

    [2]Mohammed K,Gamal E,Sun D W,et al.Application of NIR hyperspectral imaging for discrimination of lamb muscles[J].Journal of Food Engineering,2011(104):332-340.

    [3]洪添勝,喬軍,Wang Ning,等.基于高光譜圖像技術(shù)的雪花梨品質(zhì)無損檢測(cè)[J].農(nóng)業(yè)工程學(xué)報(bào),2007,23(2):151-155.

    [4]陳全勝,趙杰文,蔡健榮,等.利用高光譜圖像技術(shù)評(píng)判茶葉的質(zhì)量等級(jí)[J].光學(xué)學(xué)報(bào),2008,28(4):669-674.

    [5]薛利紅,楊林章.基于可見近紅外高光譜的菠菜硝酸鹽快速無損測(cè)定研究[J].光譜學(xué)與光譜分析,2009,29(4):926-930.

    [6]鄒偉,方慧,劉飛,等.基于高光譜圖像技術(shù)的油菜籽品種鑒別方法研究[J].浙江大學(xué)學(xué)報(bào):農(nóng)業(yè)與生命科學(xué)版,2011,37(2):175-180.

    [7]薛龍,黎靜,劉木華.基于高光譜圖像技術(shù)的水果表面農(nóng)藥殘留檢測(cè)試驗(yàn)研究[J].光學(xué)學(xué)報(bào),2008,28(12):2277-2280.

    [8]劉翠玲,隋淑霞,孫曉榮,等.近紅外光譜技術(shù)用于菠菜中毒死蜱殘留的定量分析研究[J].食品科學(xué),2008,29(7):356-358.

    [9]劉翠玲,隋淑霞,吳靜珠,等.近紅外光譜技術(shù)檢測(cè)溶液中毒死蜱含量試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2009,40(1):129-131.

    [10]Govindarajan K N,Lauren M G,Jeyamkondan S,et al.Visible/near-infrared hyperspectral imaging for beef tenderness prediction[J].Computers and Electronics in Agriculture,2008(64):225-233.

    [11]Millar S J,Whitworth M B,Chau A,et al.Mapping food composition using NIR hyperspectral imaging[J].New Food,ISSUE,2008(3):34-39.

    (責(zé)任編輯:王 寬)

    Detecting Pesticide Residue on Crystal Crown Pear Surface by Hyperspectral Imaging Technology Combined with Artificial Neural Network

    SUO Shao-zeng1,LIU Cui-ling1,WU Jing-zhu1,CHEN Xing-hai2,SUN Xiao-rong1,WU Sheng-nan1
    (1.School of Computer Science and Information Engineering,Beijing Technology and Business University,
    Beijing 100048,China;2.Zolix Instruments Co.Ltd.,Beijing 101149,China)

    The hyperspectral imaging technology combined with BP neural network was used to detect pesticide residues on surface of fruit,crystal crown pears,such as chlorpyrifos and propargite pesticide.Twenty samples were prepared for both chlorpyrifos and propargite with different concentrations.Dropping the 100 μL and 150 μL pesticide solution droplets in pear surface,the hyperspectral images were scanned in 835.467 8~1 648.356 8 nm band range.Twenty regions of interest(ROI)were available for four sample groups,from which five samples were randomly selected as predictor set.The results of models built by BP neural network for 100 μL and 150 μL pesticide solution region showed that correlation coefficient of the samples was greater than 0.99 and 0.95 respectively.RMSEC and RMSEP maximum value was respectively 0.6349,1.3239 and 1.7425,3.4417.The results indicated that 150 μL pesticide solution region was better than 100 μL in the modeling.Hyperspectral image technology combined with BP neural network method is feasible for detection of pesticide residues on pear surface.which provide a new method for detection of pesticide residues on fruit surface.

    hyperspectral imaging;BP neural network;pesticide residue;Chlorpyrifos;propargite

    TS207.3

    A

    1671-1513(2011)06-0073-05

    2011-10-30

    北京市自然科學(xué)基金資助項(xiàng)目(4073031);北京市優(yōu)秀人才資助項(xiàng)目(20081D0500300130).

    索少增,男,碩士研究生,研究方向?yàn)橹悄軝z測(cè)與數(shù)據(jù)融合理論及方法、食品安全檢測(cè)技術(shù);

    劉翠玲,女,教授,博士,主要從事智能測(cè)量技術(shù)與數(shù)據(jù)處理、系統(tǒng)建模與仿真、智能控制等方面的研究.通訊作者.

    猜你喜歡
    毒死殘留量光譜
    基于三維Saab變換的高光譜圖像壓縮方法
    氣相色譜串聯(lián)質(zhì)譜法測(cè)定茶葉中戊唑醇的殘留量
    HPLC-MS/MS法檢測(cè)花生中二嗪磷的殘留量
    黨參中二氧化硫殘留量的測(cè)定
    中成藥(2018年1期)2018-02-02 07:20:31
    夾竹桃:酷愛“模仿秀”的毒死人
    綠色中國(2017年19期)2017-05-25 13:25:14
    農(nóng)藥毒死牲畜 為何諸多說法
    星載近紅外高光譜CO2遙感進(jìn)展
    浙貝母中有機(jī)氯農(nóng)藥殘留量和二氧化硫殘留量分析
    苦味酸與牛血清蛋白相互作用的光譜研究
    鋱(Ⅲ)與PvdA作用的光譜研究
    久久久久久久国产电影| 亚洲综合色惰| 只有这里有精品99| 国产精品偷伦视频观看了| 日本av免费视频播放| 亚洲色图 男人天堂 中文字幕 | 91久久精品国产一区二区三区| 国产视频首页在线观看| 免费大片18禁| 在线观看免费高清a一片| 春色校园在线视频观看| 欧美xxxx性猛交bbbb| 天堂俺去俺来也www色官网| 美女中出高潮动态图| 免费av不卡在线播放| 欧美日韩视频精品一区| 欧美性感艳星| 国产精品国产三级国产专区5o| 日韩一区二区视频免费看| 你懂的网址亚洲精品在线观看| 精品国产一区二区久久| 在线观看美女被高潮喷水网站| 亚洲五月色婷婷综合| 亚洲第一av免费看| 成年美女黄网站色视频大全免费| 搡女人真爽免费视频火全软件| 亚洲色图综合在线观看| 国产精品久久久久久久久免| 亚洲精品美女久久久久99蜜臀 | 天天躁夜夜躁狠狠躁躁| 黄色一级大片看看| 精品国产乱码久久久久久小说| 校园人妻丝袜中文字幕| 国产亚洲一区二区精品| 午夜福利,免费看| 边亲边吃奶的免费视频| 日韩av免费高清视频| 国产亚洲精品久久久com| 国产乱人偷精品视频| 国产片特级美女逼逼视频| 久久精品国产a三级三级三级| 黄色配什么色好看| 26uuu在线亚洲综合色| 免费观看性生交大片5| 国产无遮挡羞羞视频在线观看| 亚洲精品一二三| 亚洲av国产av综合av卡| 女人精品久久久久毛片| 最后的刺客免费高清国语| 99久久人妻综合| 国产白丝娇喘喷水9色精品| 国产 一区精品| 多毛熟女@视频| 免费观看无遮挡的男女| 看免费av毛片| 亚洲色图 男人天堂 中文字幕 | 婷婷色综合大香蕉| 美女主播在线视频| 制服诱惑二区| 大香蕉久久成人网| 国产成人免费观看mmmm| 国国产精品蜜臀av免费| 免费观看无遮挡的男女| 侵犯人妻中文字幕一二三四区| 啦啦啦啦在线视频资源| 高清毛片免费看| 久热这里只有精品99| 久久久久国产网址| 亚洲第一区二区三区不卡| 成人黄色视频免费在线看| 国产亚洲欧美精品永久| 久久久久久久精品精品| 超碰97精品在线观看| 亚洲成av片中文字幕在线观看 | 99热网站在线观看| 婷婷色综合www| 在线观看一区二区三区激情| 天天操日日干夜夜撸| 一级a做视频免费观看| a 毛片基地| 夜夜爽夜夜爽视频| 日韩一本色道免费dvd| 一级毛片电影观看| 成年人午夜在线观看视频| 久久狼人影院| 日韩av在线免费看完整版不卡| 国产亚洲最大av| 国产精品国产三级国产专区5o| 欧美+日韩+精品| 五月伊人婷婷丁香| 大片免费播放器 马上看| 狠狠精品人妻久久久久久综合| 免费大片黄手机在线观看| 国产国语露脸激情在线看| 毛片一级片免费看久久久久| 97在线人人人人妻| 亚洲av日韩在线播放| 多毛熟女@视频| 一级黄片播放器| 免费在线观看完整版高清| 欧美日韩视频精品一区| 精品少妇内射三级| 日本91视频免费播放| 肉色欧美久久久久久久蜜桃| 亚洲国产精品专区欧美| 亚洲综合精品二区| 久久狼人影院| 大香蕉久久成人网| 免费观看无遮挡的男女| 老女人水多毛片| 成人亚洲精品一区在线观看| 精品国产一区二区久久| 一本大道久久a久久精品| 精品视频人人做人人爽| 欧美日韩成人在线一区二区| 人人妻人人添人人爽欧美一区卜| 成年av动漫网址| 国产在线免费精品| 草草在线视频免费看| 99九九在线精品视频| av在线老鸭窝| 国产精品 国内视频| 黑人猛操日本美女一级片| 黄网站色视频无遮挡免费观看| 国产成人一区二区在线| 亚洲av成人精品一二三区| 午夜福利影视在线免费观看| 久热这里只有精品99| 只有这里有精品99| 99久久综合免费| 乱码一卡2卡4卡精品| 免费观看av网站的网址| 黑人高潮一二区| 最近2019中文字幕mv第一页| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜av观看不卡| 亚洲,一卡二卡三卡| 97在线人人人人妻| 国产精品 国内视频| 观看美女的网站| 乱人伦中国视频| 免费观看av网站的网址| 丰满乱子伦码专区| 飞空精品影院首页| 日本免费在线观看一区| 欧美精品av麻豆av| 精品99又大又爽又粗少妇毛片| 熟女电影av网| 黑人猛操日本美女一级片| 999精品在线视频| 久久女婷五月综合色啪小说| 亚洲国产欧美在线一区| 亚洲五月色婷婷综合| 亚洲国产日韩一区二区| 亚洲精品色激情综合| 99热6这里只有精品| 国产探花极品一区二区| 草草在线视频免费看| 国产老妇伦熟女老妇高清| 国产不卡av网站在线观看| 美女主播在线视频| 老司机影院成人| 亚洲国产精品999| 亚洲欧美成人综合另类久久久| 亚洲四区av| 欧美老熟妇乱子伦牲交| 亚洲国产av影院在线观看| 老司机影院成人| 久久精品国产综合久久久 | 亚洲av电影在线观看一区二区三区| 亚洲精品中文字幕在线视频| 国产有黄有色有爽视频| 91久久精品国产一区二区三区| 丰满饥渴人妻一区二区三| 色94色欧美一区二区| 国产成人精品婷婷| 韩国高清视频一区二区三区| 99精国产麻豆久久婷婷| 青春草亚洲视频在线观看| 晚上一个人看的免费电影| 你懂的网址亚洲精品在线观看| 精品一区在线观看国产| 妹子高潮喷水视频| 全区人妻精品视频| 久久综合国产亚洲精品| 久久久国产欧美日韩av| 91国产中文字幕| 精品久久久精品久久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 熟女电影av网| 一本色道久久久久久精品综合| 校园人妻丝袜中文字幕| 99热这里只有是精品在线观看| 色94色欧美一区二区| 国产男女内射视频| 男人舔女人的私密视频| 母亲3免费完整高清在线观看 | 精品卡一卡二卡四卡免费| 成年动漫av网址| 欧美日韩精品成人综合77777| 午夜91福利影院| av视频免费观看在线观看| 国产一区二区在线观看日韩| 欧美国产精品va在线观看不卡| 97超碰精品成人国产| www.色视频.com| 久久久久国产精品人妻一区二区| 夫妻性生交免费视频一级片| 青春草亚洲视频在线观看| 精品亚洲乱码少妇综合久久| 日产精品乱码卡一卡2卡三| 少妇猛男粗大的猛烈进出视频| 国产精品久久久av美女十八| 国产 一区精品| 国产精品嫩草影院av在线观看| 国产成人一区二区在线| 啦啦啦在线观看免费高清www| 韩国高清视频一区二区三区| 免费观看在线日韩| 日日爽夜夜爽网站| 午夜视频国产福利| 亚洲综合色惰| 欧美亚洲日本最大视频资源| 制服丝袜香蕉在线| 9色porny在线观看| 午夜福利视频在线观看免费| 五月伊人婷婷丁香| 亚洲精品第二区| 少妇高潮的动态图| 在线观看国产h片| 高清在线视频一区二区三区| 婷婷色麻豆天堂久久| 好男人视频免费观看在线| 两个人免费观看高清视频| 国产极品粉嫩免费观看在线| 97精品久久久久久久久久精品| 老司机影院毛片| 一边亲一边摸免费视频| 国产深夜福利视频在线观看| 国产黄频视频在线观看| 女性被躁到高潮视频| 亚洲国产日韩一区二区| 狂野欧美激情性xxxx在线观看| 国产亚洲av片在线观看秒播厂| 亚洲av国产av综合av卡| 亚洲国产最新在线播放| 国产精品不卡视频一区二区| 精品午夜福利在线看| 欧美日韩亚洲高清精品| 日韩av在线免费看完整版不卡| 精品国产国语对白av| 久久av网站| xxx大片免费视频| 狠狠精品人妻久久久久久综合| 男女边摸边吃奶| 国产成人a∨麻豆精品| 午夜免费男女啪啪视频观看| 亚洲精品乱久久久久久| 丝瓜视频免费看黄片| 国产极品天堂在线| 欧美日韩av久久| 亚洲伊人久久精品综合| 亚洲精品,欧美精品| 国产在视频线精品| 五月开心婷婷网| 成人18禁高潮啪啪吃奶动态图| 在线观看一区二区三区激情| 日韩制服骚丝袜av| 蜜桃国产av成人99| 国产黄色免费在线视频| www日本在线高清视频| 亚洲经典国产精华液单| 日本爱情动作片www.在线观看| 免费人成在线观看视频色| 国产麻豆69| 久久99热6这里只有精品| 最近最新中文字幕大全免费视频 | 免费看不卡的av| 国产欧美另类精品又又久久亚洲欧美| 黑人巨大精品欧美一区二区蜜桃 | 乱码一卡2卡4卡精品| 欧美日韩精品成人综合77777| 侵犯人妻中文字幕一二三四区| 日本-黄色视频高清免费观看| 亚洲精品av麻豆狂野| 国产精品不卡视频一区二区| 国产男女内射视频| 日本欧美视频一区| 女人精品久久久久毛片| 韩国精品一区二区三区 | videos熟女内射| av福利片在线| 欧美丝袜亚洲另类| 欧美亚洲 丝袜 人妻 在线| 老司机影院成人| 大陆偷拍与自拍| 亚洲,一卡二卡三卡| 久久精品久久精品一区二区三区| 午夜视频国产福利| 一级a做视频免费观看| 免费黄频网站在线观看国产| 搡女人真爽免费视频火全软件| 亚洲精品中文字幕在线视频| 国产黄频视频在线观看| 国产乱来视频区| 在线看a的网站| 咕卡用的链子| 久久久国产精品麻豆| 人人妻人人添人人爽欧美一区卜| 国产伦理片在线播放av一区| 亚洲精品国产色婷婷电影| 久久狼人影院| 啦啦啦在线观看免费高清www| 新久久久久国产一级毛片| 最黄视频免费看| 国产精品国产三级国产专区5o| 亚洲av欧美aⅴ国产| 国产亚洲精品久久久com| 免费看光身美女| 久久午夜综合久久蜜桃| 精品人妻在线不人妻| 自拍欧美九色日韩亚洲蝌蚪91| 韩国精品一区二区三区 | 亚洲第一av免费看| 爱豆传媒免费全集在线观看| 一个人免费看片子| 成人国产av品久久久| av电影中文网址| 成人亚洲精品一区在线观看| 男男h啪啪无遮挡| 十八禁网站网址无遮挡| 日本爱情动作片www.在线观看| 亚洲,欧美精品.| 免费av中文字幕在线| 免费少妇av软件| 色吧在线观看| 精品亚洲成国产av| 一本色道久久久久久精品综合| 97在线人人人人妻| 色吧在线观看| 亚洲欧美色中文字幕在线| av免费在线看不卡| 在线观看人妻少妇| 亚洲国产av新网站| 极品少妇高潮喷水抽搐| 99久久中文字幕三级久久日本| 亚洲图色成人| 日韩精品有码人妻一区| 内地一区二区视频在线| 精品久久久久久电影网| 97在线人人人人妻| 成人国语在线视频| 亚洲第一av免费看| 国产激情久久老熟女| 国产xxxxx性猛交| 精品久久蜜臀av无| 亚洲三级黄色毛片| 黑人猛操日本美女一级片| 午夜免费观看性视频| 国产男女内射视频| 国产午夜精品一二区理论片| 精品第一国产精品| 免费看光身美女| av在线播放精品| 九草在线视频观看| 嫩草影院入口| 国产精品久久久久久久电影| 国内精品宾馆在线| 亚洲精品aⅴ在线观看| 亚洲精品日本国产第一区| 满18在线观看网站| 最近2019中文字幕mv第一页| 亚洲精品久久久久久婷婷小说| av女优亚洲男人天堂| 性高湖久久久久久久久免费观看| 欧美精品国产亚洲| 精品99又大又爽又粗少妇毛片| 日韩中字成人| 亚洲美女搞黄在线观看| 美女国产视频在线观看| 一级爰片在线观看| 精品久久国产蜜桃| 欧美老熟妇乱子伦牲交| 久久久久久久精品精品| 亚洲成色77777| 久久精品国产综合久久久 | 秋霞在线观看毛片| 亚洲色图 男人天堂 中文字幕 | 亚洲一区二区三区欧美精品| 欧美丝袜亚洲另类| 又大又黄又爽视频免费| 在线观看人妻少妇| 有码 亚洲区| 两个人免费观看高清视频| 丰满乱子伦码专区| 亚洲久久久国产精品| 免费观看性生交大片5| 亚洲中文av在线| 国产免费一区二区三区四区乱码| 欧美xxⅹ黑人| 久久久久国产精品人妻一区二区| 中文欧美无线码| 如何舔出高潮| 精品一区二区免费观看| 丝袜喷水一区| 日本猛色少妇xxxxx猛交久久| 80岁老熟妇乱子伦牲交| 黄片播放在线免费| 伦理电影大哥的女人| 王馨瑶露胸无遮挡在线观看| 亚洲欧美精品自产自拍| 欧美xxxx性猛交bbbb| 欧美亚洲 丝袜 人妻 在线| 超碰97精品在线观看| 观看av在线不卡| 中文字幕制服av| 搡女人真爽免费视频火全软件| 新久久久久国产一级毛片| 国产精品人妻久久久影院| 日韩精品免费视频一区二区三区 | 亚洲国产毛片av蜜桃av| √禁漫天堂资源中文www| kizo精华| 免费观看无遮挡的男女| 中文字幕人妻丝袜制服| 久久久久国产精品人妻一区二区| 成人黄色视频免费在线看| 久久人人97超碰香蕉20202| 亚洲av男天堂| 精品视频人人做人人爽| 两个人看的免费小视频| 国产在线免费精品| 春色校园在线视频观看| 1024视频免费在线观看| 飞空精品影院首页| 一级毛片 在线播放| 精品国产一区二区三区四区第35| 搡老乐熟女国产| 亚洲精品美女久久久久99蜜臀 | 成年美女黄网站色视频大全免费| 国产黄色视频一区二区在线观看| 成人黄色视频免费在线看| 午夜精品国产一区二区电影| 色婷婷av一区二区三区视频| 亚洲三级黄色毛片| 免费大片18禁| 久久久久久久久久人人人人人人| 亚洲av免费高清在线观看| 母亲3免费完整高清在线观看 | 精品亚洲乱码少妇综合久久| 国产亚洲欧美精品永久| 亚洲色图 男人天堂 中文字幕 | 亚洲精品自拍成人| 成人二区视频| 极品人妻少妇av视频| tube8黄色片| 一级毛片电影观看| 国产欧美亚洲国产| 一本大道久久a久久精品| 母亲3免费完整高清在线观看 | 精品人妻一区二区三区麻豆| 又粗又硬又长又爽又黄的视频| 97超碰精品成人国产| 亚洲精品av麻豆狂野| 国产熟女欧美一区二区| 久久精品夜色国产| 精品国产乱码久久久久久小说| 老女人水多毛片| 欧美国产精品一级二级三级| 久久鲁丝午夜福利片| 午夜91福利影院| 天天操日日干夜夜撸| 天堂俺去俺来也www色官网| 精品国产乱码久久久久久小说| 91成人精品电影| 国国产精品蜜臀av免费| 欧美激情极品国产一区二区三区 | 母亲3免费完整高清在线观看 | 国产色爽女视频免费观看| 亚洲欧洲日产国产| 日韩,欧美,国产一区二区三区| av卡一久久| 国产女主播在线喷水免费视频网站| 天堂俺去俺来也www色官网| 人人妻人人澡人人看| 国产日韩欧美在线精品| 久久热在线av| 国产男女超爽视频在线观看| 水蜜桃什么品种好| 免费在线观看黄色视频的| 国产日韩欧美在线精品| 人人澡人人妻人| 波野结衣二区三区在线| 欧美xxxx性猛交bbbb| 国产一区二区在线观看日韩| 国产xxxxx性猛交| 美女xxoo啪啪120秒动态图| 日韩熟女老妇一区二区性免费视频| 99久国产av精品国产电影| 一区二区三区精品91| 韩国精品一区二区三区 | 亚洲精华国产精华液的使用体验| 国产成人a∨麻豆精品| 永久免费av网站大全| 亚洲成人一二三区av| 99久久中文字幕三级久久日本| 久久毛片免费看一区二区三区| 国产成人免费无遮挡视频| 国产免费一级a男人的天堂| 成人二区视频| 五月玫瑰六月丁香| 十八禁网站网址无遮挡| 日本欧美国产在线视频| 91aial.com中文字幕在线观看| 国产亚洲精品第一综合不卡 | 国产精品.久久久| 亚洲av免费高清在线观看| 91国产中文字幕| 91成人精品电影| 青春草亚洲视频在线观看| 久久久久久久久久人人人人人人| 亚洲欧美日韩另类电影网站| 精品久久久久久电影网| 18+在线观看网站| 亚洲高清免费不卡视频| 亚洲激情五月婷婷啪啪| 亚洲伊人色综图| 男女国产视频网站| 老司机影院毛片| 一级黄片播放器| 国产成人91sexporn| 18+在线观看网站| videossex国产| 2022亚洲国产成人精品| 国产亚洲一区二区精品| 国产极品粉嫩免费观看在线| 亚洲人与动物交配视频| 1024视频免费在线观看| a 毛片基地| 人成视频在线观看免费观看| 下体分泌物呈黄色| 亚洲精品自拍成人| a级毛色黄片| 国产一区二区三区av在线| 婷婷成人精品国产| av福利片在线| 日日啪夜夜爽| 亚洲一级一片aⅴ在线观看| 亚洲欧美色中文字幕在线| 国产亚洲欧美精品永久| 久久这里有精品视频免费| 国产精品三级大全| 亚洲欧美成人精品一区二区| 视频中文字幕在线观看| 女人久久www免费人成看片| 99久久综合免费| 国产熟女欧美一区二区| 女性被躁到高潮视频| 久久免费观看电影| 99热这里只有是精品在线观看| 制服诱惑二区| 亚洲美女搞黄在线观看| 中国三级夫妇交换| 亚洲精品久久成人aⅴ小说| 久久ye,这里只有精品| 国产探花极品一区二区| 婷婷色综合www| 亚洲天堂av无毛| 在线观看免费视频网站a站| 91精品国产国语对白视频| 交换朋友夫妻互换小说| 亚洲av在线观看美女高潮| 久久午夜综合久久蜜桃| av.在线天堂| 欧美成人精品欧美一级黄| 久久久久视频综合| 欧美激情极品国产一区二区三区 | 男人操女人黄网站| 日韩精品免费视频一区二区三区 | 久久精品国产亚洲av涩爱| 欧美国产精品va在线观看不卡| 日韩伦理黄色片| 国产成人午夜福利电影在线观看| 最近手机中文字幕大全| 波多野结衣一区麻豆| 秋霞伦理黄片| 国产免费视频播放在线视频| 国产极品天堂在线| 亚洲第一av免费看| 久久国产亚洲av麻豆专区| 精品人妻偷拍中文字幕| 精品亚洲成a人片在线观看| 黄色怎么调成土黄色| 伊人亚洲综合成人网| 精品人妻一区二区三区麻豆| 日韩成人av中文字幕在线观看| 亚洲av欧美aⅴ国产| 极品少妇高潮喷水抽搐| 亚洲精品久久午夜乱码| 男女边摸边吃奶| 少妇精品久久久久久久| 天美传媒精品一区二区| 一级片免费观看大全| 亚洲少妇的诱惑av| 日本爱情动作片www.在线观看| 日韩av不卡免费在线播放| 丝袜喷水一区| 91精品国产国语对白视频| 蜜臀久久99精品久久宅男| 国产日韩欧美在线精品| 99久久综合免费| 亚洲熟女精品中文字幕| 少妇的逼水好多| 在现免费观看毛片| 菩萨蛮人人尽说江南好唐韦庄|