• <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作用的光譜研究
    久久毛片免费看一区二区三区| 2021少妇久久久久久久久久久| 亚洲成人手机| 久久人人97超碰香蕉20202| 日韩欧美一区视频在线观看| 黑丝袜美女国产一区| 亚洲一卡2卡3卡4卡5卡精品中文| svipshipincom国产片| 日韩欧美精品免费久久| 新久久久久国产一级毛片| 一本一本久久a久久精品综合妖精| 人人妻人人澡人人爽人人夜夜| 母亲3免费完整高清在线观看| 国产精品秋霞免费鲁丝片| 91老司机精品| 久久精品久久久久久久性| 赤兔流量卡办理| 成年美女黄网站色视频大全免费| 精品国产乱码久久久久久小说| 高清av免费在线| 久久精品国产综合久久久| 久久久久久人妻| 久久亚洲国产成人精品v| 一级片免费观看大全| 亚洲免费av在线视频| 午夜日韩欧美国产| 熟妇人妻不卡中文字幕| 亚洲成人av在线免费| 国产极品粉嫩免费观看在线| 丝袜美足系列| 三上悠亚av全集在线观看| 一本一本久久a久久精品综合妖精| 国产精品 国内视频| 午夜久久久在线观看| 国产精品熟女久久久久浪| 涩涩av久久男人的天堂| 亚洲欧美一区二区三区黑人| 一边亲一边摸免费视频| 侵犯人妻中文字幕一二三四区| 丝袜喷水一区| 久久99一区二区三区| 18禁裸乳无遮挡动漫免费视频| 国产麻豆69| 黑丝袜美女国产一区| 91aial.com中文字幕在线观看| xxx大片免费视频| 一本久久精品| 日韩电影二区| 日韩大码丰满熟妇| 久久99热这里只频精品6学生| 国产精品熟女久久久久浪| 精品亚洲成国产av| 亚洲一码二码三码区别大吗| www.精华液| 丝袜脚勾引网站| 美女午夜性视频免费| av.在线天堂| 久久久久视频综合| 久热爱精品视频在线9| 久久精品aⅴ一区二区三区四区| 精品国产乱码久久久久久男人| 中文字幕人妻熟女乱码| 丰满少妇做爰视频| 一级爰片在线观看| 一级片免费观看大全| 黑人猛操日本美女一级片| 欧美精品亚洲一区二区| 男女床上黄色一级片免费看| 中文字幕人妻熟女乱码| 无限看片的www在线观看| 91成人精品电影| 国产av精品麻豆| 亚洲av男天堂| 丁香六月欧美| 精品第一国产精品| 日韩av在线免费看完整版不卡| 免费高清在线观看日韩| 国产欧美亚洲国产| 亚洲精品美女久久久久99蜜臀 | 男女边吃奶边做爰视频| 免费高清在线观看日韩| 老汉色av国产亚洲站长工具| www.av在线官网国产| 蜜桃在线观看..| 亚洲精华国产精华液的使用体验| 欧美黄色片欧美黄色片| 人人妻人人澡人人爽人人夜夜| 成人免费观看视频高清| 伦理电影大哥的女人| 久久天堂一区二区三区四区| 亚洲国产毛片av蜜桃av| 91成人精品电影| 搡老岳熟女国产| 巨乳人妻的诱惑在线观看| 午夜福利,免费看| 亚洲欧美一区二区三区黑人| 最新的欧美精品一区二区| 七月丁香在线播放| 9色porny在线观看| 大片免费播放器 马上看| 久久精品亚洲熟妇少妇任你| 日日爽夜夜爽网站| 两个人看的免费小视频| 国产一区二区 视频在线| 在线观看免费高清a一片| 日韩成人av中文字幕在线观看| 亚洲综合色网址| 日韩成人av中文字幕在线观看| 日韩免费高清中文字幕av| 国产成人欧美| 一区二区日韩欧美中文字幕| 老司机影院毛片| 亚洲色图 男人天堂 中文字幕| 一本—道久久a久久精品蜜桃钙片| 少妇被粗大猛烈的视频| 一级片免费观看大全| 叶爱在线成人免费视频播放| 叶爱在线成人免费视频播放| 亚洲七黄色美女视频| 亚洲精品,欧美精品| 国产男人的电影天堂91| 欧美成人午夜精品| 国产熟女欧美一区二区| 亚洲精品国产一区二区精华液| 热re99久久国产66热| 新久久久久国产一级毛片| 最近最新中文字幕大全免费视频 | 亚洲色图综合在线观看| tube8黄色片| 老鸭窝网址在线观看| 伦理电影免费视频| 久久精品国产a三级三级三级| 久久精品人人爽人人爽视色| 午夜久久久在线观看| 亚洲熟女毛片儿| 亚洲精品美女久久久久99蜜臀 | 国产视频首页在线观看| 飞空精品影院首页| 精品国产露脸久久av麻豆| 午夜久久久在线观看| 亚洲欧美激情在线| videosex国产| 男女高潮啪啪啪动态图| 十八禁高潮呻吟视频| 纵有疾风起免费观看全集完整版| 亚洲精品国产色婷婷电影| 啦啦啦啦在线视频资源| 久久久久久久久久久免费av| 国产亚洲午夜精品一区二区久久| 国产日韩欧美视频二区| 人体艺术视频欧美日本| 一级毛片电影观看| 日日啪夜夜爽| 国产黄频视频在线观看| 一本大道久久a久久精品| 人妻 亚洲 视频| 午夜av观看不卡| 99久久99久久久精品蜜桃| 久久狼人影院| 一区在线观看完整版| 一区二区三区乱码不卡18| 欧美黑人精品巨大| 午夜老司机福利片| 丰满饥渴人妻一区二区三| 九草在线视频观看| 亚洲国产欧美在线一区| 欧美日韩亚洲高清精品| 国产精品久久久久久人妻精品电影 | 精品一区在线观看国产| 2021少妇久久久久久久久久久| 精品亚洲乱码少妇综合久久| 日韩精品免费视频一区二区三区| 日韩av在线免费看完整版不卡| 9色porny在线观看| 美女午夜性视频免费| 国产一区二区三区av在线| 女人爽到高潮嗷嗷叫在线视频| 无遮挡黄片免费观看| 老鸭窝网址在线观看| 日本爱情动作片www.在线观看| 欧美精品亚洲一区二区| 性高湖久久久久久久久免费观看| 午夜老司机福利片| 天堂俺去俺来也www色官网| 国产高清不卡午夜福利| 久久这里只有精品19| 午夜免费鲁丝| 免费日韩欧美在线观看| 在线 av 中文字幕| 一边摸一边抽搐一进一出视频| 丝袜美足系列| 久久精品久久久久久久性| 国产av一区二区精品久久| 国产精品一区二区在线观看99| 狠狠婷婷综合久久久久久88av| 黄片小视频在线播放| 狠狠婷婷综合久久久久久88av| 汤姆久久久久久久影院中文字幕| 一本色道久久久久久精品综合| svipshipincom国产片| 9色porny在线观看| 国产精品 欧美亚洲| 又粗又硬又长又爽又黄的视频| 久久99精品国语久久久| 色婷婷av一区二区三区视频| 色视频在线一区二区三区| 最近的中文字幕免费完整| 国产亚洲av高清不卡| 美女国产高潮福利片在线看| 国产精品国产三级国产专区5o| 中文字幕色久视频| 国产日韩一区二区三区精品不卡| netflix在线观看网站| 伊人久久大香线蕉亚洲五| 亚洲中文av在线| 欧美精品一区二区免费开放| 亚洲av福利一区| 亚洲成国产人片在线观看| 男女国产视频网站| 久久精品久久精品一区二区三区| 女人爽到高潮嗷嗷叫在线视频| netflix在线观看网站| 建设人人有责人人尽责人人享有的| 国产日韩欧美亚洲二区| 国产男女超爽视频在线观看| 亚洲欧美成人综合另类久久久| 免费久久久久久久精品成人欧美视频| 精品久久久精品久久久| 精品亚洲成a人片在线观看| 日韩一本色道免费dvd| 啦啦啦在线免费观看视频4| 丝袜在线中文字幕| 美女中出高潮动态图| 韩国精品一区二区三区| 亚洲av电影在线进入| 最近最新中文字幕免费大全7| 午夜久久久在线观看| 九色亚洲精品在线播放| 热99久久久久精品小说推荐| 午夜av观看不卡| 嫩草影院入口| 新久久久久国产一级毛片| 夫妻性生交免费视频一级片| 精品久久久久久电影网| 精品少妇内射三级| av在线app专区| 肉色欧美久久久久久久蜜桃| 色综合欧美亚洲国产小说| 一边摸一边做爽爽视频免费| 亚洲av日韩精品久久久久久密 | 国产欧美亚洲国产| 久久久精品94久久精品| 日韩 欧美 亚洲 中文字幕| 亚洲欧美成人综合另类久久久| 成年人午夜在线观看视频| 女性生殖器流出的白浆| 精品一区在线观看国产| 亚洲第一av免费看| 亚洲成国产人片在线观看| 18禁动态无遮挡网站| 色精品久久人妻99蜜桃| 久久久久网色| 精品免费久久久久久久清纯 | 丝袜在线中文字幕| av在线观看视频网站免费| 2021少妇久久久久久久久久久| 国产精品二区激情视频| av一本久久久久| av在线老鸭窝| 女人精品久久久久毛片| 亚洲国产精品国产精品| 视频区图区小说| 最近2019中文字幕mv第一页| 久久毛片免费看一区二区三区| 满18在线观看网站| 丰满乱子伦码专区| 欧美亚洲日本最大视频资源| 国产一区二区三区综合在线观看| 精品国产乱码久久久久久小说| 久久99热这里只频精品6学生| 日韩一区二区三区影片| 久久天堂一区二区三区四区| 少妇人妻 视频| 少妇 在线观看| 波野结衣二区三区在线| 一边亲一边摸免费视频| 免费看av在线观看网站| 欧美最新免费一区二区三区| www.熟女人妻精品国产| 国产亚洲午夜精品一区二区久久| 王馨瑶露胸无遮挡在线观看| 国产色婷婷99| 91国产中文字幕| 精品国产国语对白av| 久久精品aⅴ一区二区三区四区| 免费人妻精品一区二区三区视频| 最近中文字幕高清免费大全6| 久久久精品区二区三区| 亚洲成色77777| 国产亚洲欧美精品永久| 午夜精品国产一区二区电影| 日韩人妻精品一区2区三区| 国产av国产精品国产| 日本欧美视频一区| 午夜av观看不卡| 王馨瑶露胸无遮挡在线观看| 国产淫语在线视频| 欧美亚洲 丝袜 人妻 在线| 高清不卡的av网站| 成年av动漫网址| 99国产综合亚洲精品| 免费久久久久久久精品成人欧美视频| 最近2019中文字幕mv第一页| 少妇被粗大的猛进出69影院| 如日韩欧美国产精品一区二区三区| 久久 成人 亚洲| 丁香六月天网| av网站免费在线观看视频| 我的亚洲天堂| 国产xxxxx性猛交| 亚洲色图 男人天堂 中文字幕| 亚洲国产欧美网| 人妻 亚洲 视频| 男人操女人黄网站| 久久精品国产综合久久久| 天天影视国产精品| 免费女性裸体啪啪无遮挡网站| 欧美日韩视频精品一区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av电影在线观看一区二区三区| 国产成人a∨麻豆精品| 少妇 在线观看| 涩涩av久久男人的天堂| 国产一区二区三区综合在线观看| 制服丝袜香蕉在线| 成人国产麻豆网| 国产人伦9x9x在线观看| av不卡在线播放| 久久久久久久国产电影| 1024香蕉在线观看| 亚洲国产精品一区二区三区在线| 如日韩欧美国产精品一区二区三区| 国产97色在线日韩免费| 久久97久久精品| 麻豆精品久久久久久蜜桃| 波多野结衣av一区二区av| 国产欧美亚洲国产| av在线观看视频网站免费| 大片电影免费在线观看免费| 久久精品国产综合久久久| 晚上一个人看的免费电影| 国产精品久久久人人做人人爽| 夫妻性生交免费视频一级片| 久久免费观看电影| 成年女人毛片免费观看观看9 | 2021少妇久久久久久久久久久| 久久精品久久久久久久性| 一区二区日韩欧美中文字幕| 精品国产乱码久久久久久男人| 欧美日韩视频精品一区| 国产高清国产精品国产三级| 亚洲国产精品成人久久小说| 国产日韩欧美亚洲二区| 日韩熟女老妇一区二区性免费视频| 色精品久久人妻99蜜桃| 久久久国产一区二区| 亚洲综合精品二区| 熟妇人妻不卡中文字幕| 亚洲成av片中文字幕在线观看| 亚洲精品aⅴ在线观看| 桃花免费在线播放| netflix在线观看网站| 午夜精品国产一区二区电影| 91精品三级在线观看| 大香蕉久久成人网| 国产精品一区二区在线不卡| 国产黄色视频一区二区在线观看| 91成人精品电影| 国产av精品麻豆| 色婷婷av一区二区三区视频| 91精品伊人久久大香线蕉| 女的被弄到高潮叫床怎么办| 午夜福利,免费看| av在线观看视频网站免费| 男女国产视频网站| 亚洲欧美一区二区三区国产| 在线观看一区二区三区激情| 叶爱在线成人免费视频播放| 日韩一卡2卡3卡4卡2021年| 嫩草影院入口| 国产在线一区二区三区精| 人人妻人人澡人人爽人人夜夜| 亚洲色图 男人天堂 中文字幕| 久久久精品94久久精品| 欧美精品亚洲一区二区| 久久精品aⅴ一区二区三区四区| 成年美女黄网站色视频大全免费| 亚洲欧美成人综合另类久久久| 国产成人av激情在线播放| 亚洲av综合色区一区| 欧美黑人欧美精品刺激| 91国产中文字幕| 色视频在线一区二区三区| 美女视频免费永久观看网站| 国产在线视频一区二区| 最黄视频免费看| 黄色毛片三级朝国网站| 熟妇人妻不卡中文字幕| 午夜福利免费观看在线| 日本wwww免费看| 久久精品国产综合久久久| 涩涩av久久男人的天堂| 国产av码专区亚洲av| 各种免费的搞黄视频| 99久国产av精品国产电影| 亚洲国产看品久久| 久久免费观看电影| 国产精品久久久久久精品古装| 精品人妻在线不人妻| 亚洲精品视频女| 国产精品女同一区二区软件| 超碰97精品在线观看| 国产欧美亚洲国产| 精品国产乱码久久久久久小说| 亚洲色图综合在线观看| 日韩制服丝袜自拍偷拍| 国产av精品麻豆| 亚洲国产欧美在线一区| 欧美日韩av久久| 亚洲第一青青草原| 一区在线观看完整版| 国产成人a∨麻豆精品| 涩涩av久久男人的天堂| 国产精品秋霞免费鲁丝片| 久久 成人 亚洲| 国产亚洲精品第一综合不卡| 久久99一区二区三区| 国产成人啪精品午夜网站| 国产乱来视频区| 亚洲精品美女久久av网站| 日韩制服骚丝袜av| 欧美 亚洲 国产 日韩一| 欧美黑人精品巨大| 亚洲精品,欧美精品| 卡戴珊不雅视频在线播放| 久久久国产欧美日韩av| 母亲3免费完整高清在线观看| 国产激情久久老熟女| 日韩一区二区三区影片| 亚洲国产日韩一区二区| 国产成人免费观看mmmm| 9热在线视频观看99| 亚洲,欧美,日韩| 少妇人妻 视频| 在现免费观看毛片| 亚洲人成网站在线观看播放| 中文欧美无线码| netflix在线观看网站| 成年动漫av网址| 国产一区二区三区av在线| 好男人视频免费观看在线| 美女中出高潮动态图| 街头女战士在线观看网站| 午夜日韩欧美国产| 精品一区二区三卡| 香蕉国产在线看| 精品久久久精品久久久| 久久久精品免费免费高清| 婷婷色综合大香蕉| 2021少妇久久久久久久久久久| 久久精品熟女亚洲av麻豆精品| 亚洲精品视频女| 国产成人a∨麻豆精品| 欧美97在线视频| 国产精品一国产av| 天天躁日日躁夜夜躁夜夜| 亚洲人成网站在线观看播放| 在线观看免费午夜福利视频| 国产精品嫩草影院av在线观看| 少妇被粗大的猛进出69影院| 国产欧美日韩综合在线一区二区| 精品卡一卡二卡四卡免费| 亚洲精品一二三| 中文精品一卡2卡3卡4更新| 免费看av在线观看网站| 精品国产露脸久久av麻豆| 精品一区在线观看国产| 成年av动漫网址| 中文天堂在线官网| 国产av国产精品国产| 日韩电影二区| av福利片在线| 亚洲欧美一区二区三区国产| 一级片免费观看大全| 老司机在亚洲福利影院| 亚洲av中文av极速乱| 老司机在亚洲福利影院| 免费观看av网站的网址| 日韩欧美精品免费久久| 欧美老熟妇乱子伦牲交| 日韩一卡2卡3卡4卡2021年| 一级黄片播放器| 精品国产一区二区三区久久久樱花| 成人三级做爰电影| 久久久久久久国产电影| 天天操日日干夜夜撸| av国产精品久久久久影院| 美女中出高潮动态图| 免费在线观看黄色视频的| 国产精品一区二区精品视频观看| 日本黄色日本黄色录像| 国产亚洲欧美精品永久| 亚洲美女搞黄在线观看| 麻豆av在线久日| 午夜91福利影院| 欧美人与善性xxx| av网站在线播放免费| 国产深夜福利视频在线观看| 午夜日本视频在线| 中文字幕精品免费在线观看视频| 国产不卡av网站在线观看| 成人国产av品久久久| 黑人猛操日本美女一级片| 亚洲国产精品国产精品| 国产精品成人在线| 热99国产精品久久久久久7| 欧美激情极品国产一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 妹子高潮喷水视频| 永久免费av网站大全| 午夜日韩欧美国产| 男女床上黄色一级片免费看| 国产男女内射视频| 男女午夜视频在线观看| 国产熟女午夜一区二区三区| 国产精品免费视频内射| 51午夜福利影视在线观看| 午夜福利视频精品| 高清视频免费观看一区二区| 国产精品免费大片| 亚洲成人免费av在线播放| 99香蕉大伊视频| 秋霞在线观看毛片| 纵有疾风起免费观看全集完整版| 国精品久久久久久国模美| 亚洲精华国产精华液的使用体验| 亚洲伊人久久精品综合| 老汉色∧v一级毛片| 亚洲av中文av极速乱| 国产精品熟女久久久久浪| 夫妻性生交免费视频一级片| xxxhd国产人妻xxx| 国产又色又爽无遮挡免| 大片电影免费在线观看免费| 自线自在国产av| 国产精品一国产av| 一级黄片播放器| 男的添女的下面高潮视频| 又黄又粗又硬又大视频| 亚洲久久久国产精品| 久久ye,这里只有精品| 精品视频人人做人人爽| 亚洲一级一片aⅴ在线观看| 久久久久久人妻| 69精品国产乱码久久久| 韩国高清视频一区二区三区| 日韩av不卡免费在线播放| 五月天丁香电影| 黄频高清免费视频| 丰满少妇做爰视频| 国产精品一区二区在线观看99| 国产日韩一区二区三区精品不卡| 一个人免费看片子| 操美女的视频在线观看| 男人舔女人的私密视频| 别揉我奶头~嗯~啊~动态视频 | 国产熟女欧美一区二区| 午夜免费鲁丝| 男女床上黄色一级片免费看| 国产精品一区二区在线不卡| 午夜免费观看性视频| 国产男女内射视频| 美女午夜性视频免费| 久久久久国产一级毛片高清牌| 黑人巨大精品欧美一区二区蜜桃| 街头女战士在线观看网站| 亚洲精品成人av观看孕妇| 日韩av在线免费看完整版不卡| 亚洲自偷自拍图片 自拍| 免费高清在线观看日韩| 午夜日韩欧美国产| 午夜福利免费观看在线| 中文字幕高清在线视频| 成人亚洲精品一区在线观看| 少妇猛男粗大的猛烈进出视频| 午夜福利,免费看| 久久热在线av| 亚洲,欧美,日韩| 色视频在线一区二区三区| 99re6热这里在线精品视频| 我的亚洲天堂| 国产一区二区三区综合在线观看| 制服人妻中文乱码| 免费在线观看完整版高清| 精品第一国产精品| 国产精品二区激情视频| 日本av手机在线免费观看| 男女床上黄色一级片免费看| 天美传媒精品一区二区| 欧美日本中文国产一区发布| 日日爽夜夜爽网站| 日韩熟女老妇一区二区性免费视频| 色视频在线一区二区三区| 日本黄色日本黄色录像| 国产av国产精品国产|