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

    基于高光譜數(shù)據(jù)研究應(yīng)用近紅外相機(jī)加裝濾光片實(shí)現(xiàn)玉米葉片水分測(cè)量的關(guān)鍵參數(shù)

    2021-10-17 01:39:38李朋成趙龍蓮李軍會(huì)
    光譜學(xué)與光譜分析 2021年10期
    關(guān)鍵詞:濾光片波長(zhǎng)光譜

    李朋成, 劉 含, 趙龍蓮, 李軍會(huì)*

    1. 中國農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院, 北京 100083 2. 中國農(nóng)業(yè)大學(xué), 教育部現(xiàn)代精細(xì)農(nóng)業(yè)系統(tǒng)集成研究重點(diǎn)實(shí)驗(yàn), 北京 100083

    引 言

    植物葉片是光合作用的重要器官, 葉片含水量的變化會(huì)影響到葉綠素活體的改變, 水分虧缺會(huì)直接影響作物的生長(zhǎng)[1]。 實(shí)現(xiàn)活體作物葉片不同區(qū)域以及近地冠層的水分含量檢測(cè), 可為研究作物生長(zhǎng)中的水分脅迫響應(yīng)機(jī)制, 以及旱情監(jiān)測(cè)、 精準(zhǔn)灌溉等提供有效技術(shù)手段[2]。 當(dāng)前, 植物葉片水分多采用烘干法等化學(xué)測(cè)量方法, 這些方法費(fèi)時(shí)費(fèi)力且不可重復(fù), 而光譜技術(shù)具有精確、 無損、 快速等特點(diǎn)[3]。 近紅外高光譜成像技術(shù)是新一代光電檢測(cè)技術(shù), 同時(shí)結(jié)合了近紅外和高光譜成像的優(yōu)點(diǎn), 既包含了樣品的豐富化學(xué)特征信息, 又可得到像素級(jí)圖譜信息[4], 但儀器具有大體積、 大重量、 低光通量的特性, 需要外部光源, 難以在田間實(shí)現(xiàn)活體葉片檢測(cè)。

    目前, 近紅外相機(jī)具有小體積、 小重量, 可調(diào)節(jié)的光圈和大口徑的鏡頭使得不需要輔助光源也能實(shí)現(xiàn)在線測(cè)量[5]。 由于近紅外相機(jī)不具有波長(zhǎng)分辨能力, 一般多應(yīng)用于影像分析[6], 在近紅外相機(jī)鏡頭前加裝不同中心波長(zhǎng)的帶通濾光片, 使選擇的波段與待測(cè)物質(zhì)成分具有強(qiáng)相關(guān)性, 便可設(shè)計(jì)出有波長(zhǎng)分辨能力的多光譜儀器[7-8]; 如高建樹[9]等基于水和冰在近紅外光譜下的不同反射率, 選擇1 160和1 260 nm兩個(gè)波段的濾光片定量識(shí)別機(jī)翼積冰厚度; Ishii[10]等觀測(cè)脂肪對(duì)近紅外光的吸收強(qiáng)度, 發(fā)現(xiàn)在1 200 nm波段量化動(dòng)脈粥樣硬化的潛力, 開發(fā)出血管的光譜成像系統(tǒng); 在植物葉片水分、 氮素等成分檢測(cè)方面的研究鮮有文獻(xiàn)報(bào)道。

    本研究基于活體玉米葉片的近紅外高光譜數(shù)據(jù), 依據(jù)濾光片透光分布函數(shù)得到不同帶寬的仿真數(shù)據(jù), 應(yīng)用插值方法得到固定帶寬下的中心波長(zhǎng)偏移仿真數(shù)據(jù), 對(duì)近紅外相機(jī)加裝濾光片實(shí)現(xiàn)活體玉米葉片水分測(cè)量的關(guān)鍵參數(shù)進(jìn)行研究, 為多光譜成像系統(tǒng)檢測(cè)裝置提供技術(shù)參數(shù)支撐和環(huán)境工作范圍建議。 具體研究?jī)?nèi)容包括: (1)特征波長(zhǎng)位置及個(gè)數(shù)篩選; (2)不同帶寬數(shù)據(jù)對(duì)建模結(jié)果的影響; (3)中心波長(zhǎng)偏移對(duì)模型預(yù)測(cè)性能的影響。

    1 實(shí)驗(yàn)部分

    1.1 儀器

    Gaia Sorter近紅外高光譜儀器(四川雙利合譜科技有限公司)。 其系統(tǒng)配置為: 4個(gè)溴鎢燈組成的均勻光源; 56 mm定焦近紅外相機(jī)鏡頭; N17E光譜儀, 其光譜范圍為900~1 700 nm; AVT(320×256)探測(cè)頭; 計(jì)算機(jī)箱等。 其CCD相機(jī)工作電壓: DC12V; 采集軟件: Spec View(SpecView Ltd., Uckfield, UK)。 采集模式為線掃描, 高光譜儀器工作機(jī)制: 實(shí)驗(yàn)樣品放置在電子控制移動(dòng)平臺(tái)上, 通過推掃法獲得待測(cè)樣品的高光譜立方體信息。 電子天平, 精度0.000 1 g。 電熱恒溫干燥箱, 精度0.1 ℃。

    1.2 樣品

    實(shí)驗(yàn)樣品為不同生育期(苗期、 拔節(jié)期、 抽穗期、 灌漿期、 成熟期)的鮮活玉米葉片共60個(gè), 采集于中國農(nóng)業(yè)大學(xué)上莊實(shí)驗(yàn)站。 運(yùn)行SpecVIEW軟件, 系統(tǒng)預(yù)熱30 min后進(jìn)行聚焦等一系列操作, 最終確定電動(dòng)平臺(tái)移動(dòng)速度為0.54 cm·s-1, 曝光時(shí)間為0.09 s。 采集時(shí)將樣品放置在推掃式傳送帶上, 每次放置一個(gè)葉片, 共采集得到60幅原始高光譜圖像。

    采用烘干法測(cè)量樣品含水量。 將樣品掃描后測(cè)其鮮重, 然后在105 ℃恒溫箱中烘干4 h。 為避免過熱快速吸收水蒸氣, 造成測(cè)量精度的誤差, 干燥冷卻至室溫平衡3 h后測(cè)量干重。 鮮活玉米葉片含水量在60%~90%之間。

    1.3 方法

    1.3.1 高光譜數(shù)據(jù)校正與預(yù)處理

    為消除CCD相機(jī)的暗電流以及不同波長(zhǎng)下光源強(qiáng)度分布不均勻造成的噪聲干擾, 對(duì)采集到的原始高光譜圖像進(jìn)行校正。 相機(jī)鏡頭完全被相機(jī)鏡頭蓋蓋住采集的圖像為暗背景數(shù)據(jù), 高光譜儀的相機(jī)對(duì)準(zhǔn)標(biāo)準(zhǔn)白板采集的圖像為標(biāo)準(zhǔn)白板數(shù)據(jù)。 校正公式如(1)所示

    (1)

    式(1)中,Ai為校正后的光譜吸光度數(shù)據(jù);B為暗背景數(shù)據(jù), 理論0%反射率;W為標(biāo)準(zhǔn)白板數(shù)據(jù): 理論99%反射率。

    使用ENVI5.1采集不同葉肉的ROI(region of interest)并計(jì)算求得近紅外高光譜樣品數(shù)(60)×波長(zhǎng)數(shù)(256)平均吸光度的光譜數(shù)據(jù)矩陣。

    為消除光譜的基線漂移、 噪聲、 樣本大小以及環(huán)境因素引起的譜線漂移, 常常在特征分析前采用預(yù)處理來濾除干擾信號(hào)[11]。 采用標(biāo)準(zhǔn)正態(tài)變換(standard normal variate, SNV)對(duì)高光譜數(shù)據(jù)矩陣進(jìn)行預(yù)處理, SNV使用平均值和標(biāo)準(zhǔn)偏差作為校正參數(shù), 可在很大程度上消除樣品表面散射影響。

    1.3.2 特征波長(zhǎng)篩選

    選擇出與水分信息相關(guān)性高且個(gè)數(shù)少的特征波長(zhǎng), 是應(yīng)用近紅外相機(jī)加裝濾光片實(shí)現(xiàn)水分測(cè)量的關(guān)鍵, 可在滿足檢測(cè)準(zhǔn)確度基礎(chǔ)上簡(jiǎn)化裝置結(jié)構(gòu)。 本文基于高分辨率的近紅外高光譜數(shù)據(jù)結(jié)合逐步回歸方法進(jìn)行水分檢測(cè)所需特征波長(zhǎng)的篩選。 逐步回歸是一種較優(yōu)的選元方法可用來消除相關(guān)性低和共線性度高的波長(zhǎng), 每引入一個(gè)波長(zhǎng), 逐個(gè)對(duì)已選入的波長(zhǎng)進(jìn)行偏F檢驗(yàn), 直到新引入的波長(zhǎng)不再顯著[12]。 利用SPSS 24軟件的逐步回歸方法實(shí)現(xiàn)特征波長(zhǎng)篩選。

    1.3.3 不同帶寬數(shù)據(jù)的仿真獲取

    選定特征波長(zhǎng)即濾光片中心波長(zhǎng)后, 波長(zhǎng)帶寬寬度是另一個(gè)實(shí)現(xiàn)水分測(cè)量的關(guān)鍵參數(shù), 帶寬過寬雖然可以提高光通量但分辨率會(huì)降低, 在近紅外波段體現(xiàn)出不同基團(tuán)信息的耦合重疊, 直接影響建模效果; 目前市場(chǎng)上的濾光片帶寬參數(shù)一般為10~300 nm。 本研究結(jié)合不同型號(hào)濾光片實(shí)際帶寬參數(shù), 基于鮮活玉米葉片的近紅外高光譜數(shù)據(jù), 仿真得到不同帶寬的數(shù)據(jù), 通過建模效果對(duì)比得出合適的帶寬范圍, 為濾光片選型提供依據(jù)。

    依據(jù)濾光片的透光近似分布函數(shù), 代入中心波長(zhǎng)和帶寬參數(shù), 可得到相應(yīng)波長(zhǎng)和帶寬下的透光分布曲線, 然后, 將對(duì)應(yīng)波長(zhǎng)下的高光譜數(shù)據(jù)與分布曲線作卷積運(yùn)算, 即可得到不同中心波長(zhǎng)和帶寬下符合濾光片實(shí)際透光特性的仿真數(shù)據(jù)。 其中, 帶通濾光片透光近似分布函數(shù)見式(2)[13]; 卷積計(jì)算見式(3)。

    (2)

    (3)

    式(3)中,Si為透過濾光片的吸光度仿真數(shù)據(jù);Ai為濾光片帶寬范圍內(nèi)對(duì)應(yīng)波長(zhǎng)下的高光譜吸光度數(shù)據(jù)。

    1.3.4 中心波長(zhǎng)偏移數(shù)據(jù)的仿真獲取

    加裝的帶通濾光片采用技術(shù)成熟的薄膜干涉濾光片產(chǎn)品, 該濾光片膜層材料的熱膨脹系數(shù)(coefficient of thermal expansion, CTE)和折射率(refractive index)會(huì)隨著工作環(huán)境溫度變化, 導(dǎo)致中心波長(zhǎng)發(fā)生偏移[14], 一般情況下, 溫度影響造成的中心波長(zhǎng)偏移小于10-3nm·℃-1[15]。 工作環(huán)境溫度變化過大, 會(huì)導(dǎo)致中心波長(zhǎng)偏移過大, 使固定中心波長(zhǎng)下所建模型的預(yù)測(cè)性能大大降低。 為推測(cè)加裝濾光片所建模型受環(huán)境溫度的影響, 應(yīng)用固定中心波長(zhǎng)和帶寬所建模型, 預(yù)測(cè)分析不同偏移量數(shù)據(jù)的誤差情況, 得出多光譜成像檢測(cè)裝置的環(huán)境工作溫度范圍。

    符合濾光片實(shí)際偏移特性的仿真數(shù)據(jù)獲取過程如下:

    (1)基于鮮活玉米葉片高波長(zhǎng)分辨率的近紅外高光譜數(shù)據(jù), 通過線性插值計(jì)算得到更高波長(zhǎng)分辨數(shù)據(jù), 插值方法見式(4)。

    (2)固定帶寬下, 將中心波長(zhǎng)偏移后的波長(zhǎng)代入式(2), 得到相應(yīng)波長(zhǎng)和帶寬下的透光分布曲線。

    (3)將相應(yīng)波長(zhǎng)對(duì)應(yīng)的高光譜數(shù)據(jù)與分布曲線應(yīng)用式(3)作卷積運(yùn)算, 得到符合濾光片實(shí)際偏移特性的仿真數(shù)據(jù)。

    (4)

    對(duì)于(xi,xi+1)范圍內(nèi)任意波長(zhǎng)x, 吸光度y可由式(4)計(jì)算得到。 其中, (xi,yi)表示波長(zhǎng)和吸光度, (xi+1,yi+1)為后一個(gè)波長(zhǎng)所對(duì)應(yīng)的波長(zhǎng)和吸光度。 利用MATLAB 2018a中“spline”函數(shù)實(shí)現(xiàn)插值計(jì)算。

    1.3.5 模型的評(píng)價(jià)指標(biāo)

    使用決定系數(shù)(determination coefficient,R2)和均方根誤差(root mean squared error, RMSE)[16]對(duì)模型性能進(jìn)行評(píng)價(jià), RMSE單位為百分點(diǎn)。

    2 結(jié)果與討論

    2.1 特征波長(zhǎng)篩選

    應(yīng)用1.3.2中介紹的逐步回歸方法對(duì)玉米鮮活葉片近紅外高光譜數(shù)據(jù)進(jìn)行特征波長(zhǎng)篩選, 對(duì)p≤0.05的F值引入相關(guān)波長(zhǎng),F(xiàn)值(p≥0.1)時(shí)剔除該波長(zhǎng)。 篩選出的特征波長(zhǎng)及其對(duì)應(yīng)回歸模型的R2和RMSE結(jié)果的情況見表1。

    表1 逐步回歸模型結(jié)果

    由表1可以看出, 當(dāng)選擇2個(gè)或3個(gè)特征波長(zhǎng)時(shí)所建模型已能較好實(shí)現(xiàn)水分測(cè)量。 考慮到多光譜檢測(cè)裝置的結(jié)構(gòu)和成本, 選擇1 151和1 403 nm作為與水分密切相關(guān)的2個(gè)特征波長(zhǎng)。

    選用市場(chǎng)上現(xiàn)有濾光片其中心波長(zhǎng)為1 150和1 400 nm, 為驗(yàn)證是否合適, 將60份不同生育期的樣品以2∶1比例隨機(jī)劃分為建模集和預(yù)測(cè)集, 兩個(gè)模型的R2和RMSE情況見表2。

    表2 基于特征波長(zhǎng)的建模和預(yù)測(cè)結(jié)果

    表2結(jié)果表明兩組相鄰波長(zhǎng)建模和驗(yàn)證結(jié)果的R2和RMSE沒有明顯差異。 由近紅外基團(tuán)吸收特征可知[16], 1 150 nm附近主要包含CH基團(tuán)倍頻吸收, 而1 400 nm附近主要包含OH和CH基團(tuán)的倍頻及合頻吸收。 選擇這兩個(gè)特征波長(zhǎng)建立回歸模型時(shí), CH基團(tuán)的吸收信息可實(shí)現(xiàn)相互抵消, 主要顯現(xiàn)水分的特征吸收, 應(yīng)用數(shù)學(xué)方法選擇的特征波長(zhǎng)具有基團(tuán)吸收特征的物理解釋。 因此, 本文所選擇特征波長(zhǎng)合理可行。

    2.2 不同帶寬數(shù)據(jù)對(duì)建模結(jié)果的影響

    將2.1中篩選的特征波長(zhǎng)(λc=1 150和1 400 nm)和不同帶寬(FWHM=10, 25, …, 300 nm)參數(shù)代入式(2)(透光近似分布函數(shù)), 得到相應(yīng)波長(zhǎng)和帶寬下的透光率的分布曲線; 并由式(3), 將對(duì)應(yīng)波長(zhǎng)下的高光譜數(shù)據(jù)與分布曲線作卷積運(yùn)算, 得到不同中心波長(zhǎng)和帶寬下符合濾光片實(shí)際透光特性的仿真數(shù)據(jù)。 應(yīng)用不同參數(shù)下的仿真數(shù)據(jù), 建立多元線性回歸模型(multivarate linear regression, MLR), 不同模型結(jié)果情況見表3。

    表3 不同帶寬建模結(jié)果

    由表3結(jié)果可以看出, 帶寬小于100 nm的四個(gè)模型, 其R2和RMSE沒有明顯差異, 帶寬大于100 nm模型的決定系數(shù)和準(zhǔn)確度都有明顯降低。 帶寬小于100 nm時(shí), 所建模型預(yù)測(cè)水分含量的均方根誤差可在1.5個(gè)百分點(diǎn)以內(nèi), 滿足一般農(nóng)業(yè)分析要求。 為進(jìn)一步驗(yàn)證模型預(yù)測(cè)水分的可靠性, 以及所選特征波長(zhǎng)和帶寬范圍可行, 將25 nm帶寬的仿真數(shù)據(jù)集劃分為建模集與預(yù)測(cè)集, 劃分方式同2.1節(jié)一致, 其中, 建模集的R2和RMSE分別為0.968和1.186; 預(yù)測(cè)集的R2和RMSE分別為0.960和1.298, 圖1為建模集和預(yù)測(cè)集的水分參考值與預(yù)測(cè)值散點(diǎn)圖。 圖1結(jié)果表明, 選擇1 150和1 400 nm作為特征波長(zhǎng), 在帶寬為25 nm下所建水分含量模型, 可實(shí)現(xiàn)玉米葉片水分的準(zhǔn)確檢測(cè)。

    圖1 25 nm帶寬下建模集(a)和預(yù)測(cè)集(b)的水分參考值與預(yù)測(cè)值散點(diǎn)圖

    2.3 中心波長(zhǎng)偏移對(duì)模型預(yù)測(cè)性能的影響

    應(yīng)用式(4), 得到波長(zhǎng)分辨率為0.01~0.5 nm的數(shù)據(jù); 將偏移后的中心波長(zhǎng)和帶寬(25 nm)代入式(2), 其中, 偏移量依次設(shè)置為(0.01, 0.02, …, 0.5 nm), 得到中心波長(zhǎng)偏移后的透光分布曲線; 將相應(yīng)波長(zhǎng)的高光譜插值數(shù)據(jù)應(yīng)用式(3)作卷積運(yùn)算, 得到符合濾光片實(shí)際偏移特性的仿真數(shù)據(jù)。 應(yīng)用表3中, 中心波長(zhǎng)(λc=1 150和1 400 nm)和帶寬(25 nm)的模型, 預(yù)測(cè)不同偏移量數(shù)據(jù)的結(jié)果情況見表4。

    由表4中數(shù)據(jù)可以看出, 隨著中心波長(zhǎng)偏移逐漸增大, 模型預(yù)測(cè)誤差也隨之增大, 表明波長(zhǎng)準(zhǔn)確性對(duì)近紅外分析非常關(guān)鍵, 波長(zhǎng)偏移過大時(shí)模型預(yù)測(cè)能力顯著下降。 但濾光片中心波長(zhǎng)偏移受溫度影響較小, 而且玉米一般是夏季室外作物, 通常在室外環(huán)境使用干涉濾光片的溫度在0和+50 ℃之間, 由1.3.4節(jié)中介紹的干涉濾光片中心波長(zhǎng)偏移與溫度的定量關(guān)系可得, 外界環(huán)境溫度在50 ℃范圍內(nèi)變化, 等價(jià)于中心波長(zhǎng)漂移量為0.05 nm, 而在此漂移量下的預(yù)測(cè)誤差僅高出不漂移下預(yù)測(cè)誤差3%左右, 可忽略, 詳細(xì)結(jié)果情況見表4。

    表4 中心波長(zhǎng)(1 150和1 400 nm)和帶寬(25 nm)模型預(yù)測(cè)不同偏移量數(shù)據(jù)的結(jié)果

    3 結(jié) 論

    將小體積、 小重量、 高光通量近紅外相機(jī)加裝濾光片, 形成多光譜成像檢測(cè)裝置提供技術(shù)參數(shù)支撐和環(huán)境工作范圍, 從而實(shí)現(xiàn)田間活體葉片水分含量成像檢測(cè); 通過應(yīng)用鮮活玉米葉片的近紅外高光譜數(shù)據(jù)進(jìn)行仿真研究, 結(jié)果表明: 濾光片特征波長(zhǎng)分別為1 150和1 400 nm, 帶寬小于100 nm時(shí), 可實(shí)現(xiàn)準(zhǔn)確檢測(cè), 環(huán)境溫度在50 ℃范圍內(nèi)對(duì)檢測(cè)結(jié)果影響不大。 基于本文研究結(jié)果, 實(shí)際裝置的搭建工作已開展, 裝置的實(shí)現(xiàn)可為現(xiàn)代農(nóng)業(yè)作物生理、 生產(chǎn)研究提供新型、 有效手段。

    猜你喜歡
    濾光片波長(zhǎng)光譜
    HPLC-PDA雙波長(zhǎng)法同時(shí)測(cè)定四季草片中沒食子酸和槲皮苷的含量
    基于三維Saab變換的高光譜圖像壓縮方法
    雙波長(zhǎng)激光治療慢性牙周炎的療效觀察
    基于等效層理論的薄膜濾光片中心波長(zhǎng)消偏振膜系設(shè)計(jì)
    可調(diào)諧濾光片的透射譜研究
    日本研發(fā)出可完全覆蓋可見光波長(zhǎng)的LED光源
    中國照明(2016年4期)2016-05-17 06:16:15
    星載近紅外高光譜CO2遙感進(jìn)展
    近地層紫外光通信“日盲”濾光片技術(shù)研究
    激光與紅外(2015年8期)2015-03-23 06:07:08
    便攜式多用途光波波長(zhǎng)測(cè)量?jī)x
    苦味酸與牛血清蛋白相互作用的光譜研究
    色吧在线观看| 久久这里只有精品19| 国产精品.久久久| 狠狠精品人妻久久久久久综合| 视频中文字幕在线观看| av福利片在线| 欧美日韩综合久久久久久| 九九爱精品视频在线观看| 精品国产国语对白av| 亚洲欧洲日产国产| 男女啪啪激烈高潮av片| 欧美老熟妇乱子伦牲交| 香蕉精品网在线| 9色porny在线观看| 视频区图区小说| 亚洲人成77777在线视频| 欧美精品高潮呻吟av久久| 99久国产av精品国产电影| 免费久久久久久久精品成人欧美视频 | 亚洲精品中文字幕在线视频| 校园人妻丝袜中文字幕| 亚洲欧美成人综合另类久久久| 国产有黄有色有爽视频| 久久狼人影院| 亚洲欧美日韩另类电影网站| 久久精品久久精品一区二区三区| 亚洲国产av影院在线观看| 午夜福利网站1000一区二区三区| 18禁国产床啪视频网站| 婷婷色av中文字幕| 丰满少妇做爰视频| √禁漫天堂资源中文www| 制服诱惑二区| 在线天堂最新版资源| 成年动漫av网址| 少妇熟女欧美另类| 亚洲欧美日韩卡通动漫| 黑人高潮一二区| 成人毛片a级毛片在线播放| 久久久久久久久久成人| 国产日韩欧美亚洲二区| 欧美精品亚洲一区二区| 搡老乐熟女国产| 亚洲av成人精品一二三区| 高清黄色对白视频在线免费看| 人妻一区二区av| 亚洲精品一区蜜桃| 少妇人妻精品综合一区二区| av卡一久久| 久久人人爽人人爽人人片va| 如日韩欧美国产精品一区二区三区| 夫妻性生交免费视频一级片| 考比视频在线观看| 狠狠精品人妻久久久久久综合| 老司机亚洲免费影院| 天堂8中文在线网| 一级毛片 在线播放| 久久99一区二区三区| 捣出白浆h1v1| 国产色婷婷99| 99热这里只有是精品在线观看| 亚洲成人手机| 国产精品免费大片| 制服诱惑二区| 永久网站在线| 99国产精品免费福利视频| 免费少妇av软件| 性色avwww在线观看| 丰满饥渴人妻一区二区三| 久久久久久久国产电影| 蜜桃在线观看..| 乱人伦中国视频| 久久亚洲国产成人精品v| 少妇高潮的动态图| 高清黄色对白视频在线免费看| 成人手机av| 国产男人的电影天堂91| 人妻人人澡人人爽人人| a级毛色黄片| 日韩制服丝袜自拍偷拍| 蜜桃在线观看..| 久久国产精品大桥未久av| 男的添女的下面高潮视频| 国产成人精品久久久久久| 制服诱惑二区| 免费观看av网站的网址| 欧美激情国产日韩精品一区| 日韩电影二区| 老女人水多毛片| 久久久精品区二区三区| 丝瓜视频免费看黄片| 亚洲人成网站在线观看播放| 亚洲人成网站在线观看播放| 乱码一卡2卡4卡精品| 如何舔出高潮| 51国产日韩欧美| 亚洲欧美清纯卡通| 视频中文字幕在线观看| 国产男女内射视频| 亚洲国产av影院在线观看| 在线观看免费日韩欧美大片| 国产福利在线免费观看视频| 色94色欧美一区二区| 人妻一区二区av| 成年动漫av网址| 涩涩av久久男人的天堂| 如何舔出高潮| 国产乱来视频区| 建设人人有责人人尽责人人享有的| 国语对白做爰xxxⅹ性视频网站| 国产精品一区二区在线不卡| 精品久久久久久电影网| 国产精品一国产av| 一级,二级,三级黄色视频| 90打野战视频偷拍视频| 欧美老熟妇乱子伦牲交| 久久精品熟女亚洲av麻豆精品| 九九爱精品视频在线观看| 亚洲成人一二三区av| 99热全是精品| 午夜激情av网站| 99久国产av精品国产电影| 久久精品国产综合久久久 | 欧美97在线视频| 伊人久久国产一区二区| 亚洲av综合色区一区| 在线免费观看不下载黄p国产| 国产在视频线精品| 丝袜美足系列| 一区二区三区四区激情视频| 激情五月婷婷亚洲| 国产亚洲欧美精品永久| 日韩大片免费观看网站| 日韩制服骚丝袜av| av福利片在线| 成人黄色视频免费在线看| 久久99精品国语久久久| 婷婷色麻豆天堂久久| 国产成人精品无人区| 欧美日韩综合久久久久久| 女人久久www免费人成看片| 亚洲成av片中文字幕在线观看 | 午夜91福利影院| 少妇猛男粗大的猛烈进出视频| 狂野欧美激情性xxxx在线观看| 国产激情久久老熟女| 一级毛片我不卡| 成年人午夜在线观看视频| 2021少妇久久久久久久久久久| 欧美日韩国产mv在线观看视频| 午夜福利视频精品| 成人18禁高潮啪啪吃奶动态图| 亚洲伊人久久精品综合| 精品一区二区免费观看| 搡女人真爽免费视频火全软件| 欧美日韩精品成人综合77777| 午夜av观看不卡| 日本av免费视频播放| av播播在线观看一区| 一区二区三区精品91| 一区二区三区精品91| 免费人妻精品一区二区三区视频| 亚洲av中文av极速乱| 精品一品国产午夜福利视频| tube8黄色片| 在线观看一区二区三区激情| 少妇精品久久久久久久| 国产又爽黄色视频| 一级爰片在线观看| 国产一区有黄有色的免费视频| 91成人精品电影| av在线app专区| 丝袜喷水一区| av在线老鸭窝| 久久国产精品男人的天堂亚洲 | 欧美精品一区二区免费开放| 在线观看免费日韩欧美大片| 亚洲精品,欧美精品| 国产精品一国产av| 日韩中字成人| 亚洲国产av影院在线观看| 亚洲av男天堂| 又粗又硬又长又爽又黄的视频| 中文精品一卡2卡3卡4更新| 一边亲一边摸免费视频| 国产视频首页在线观看| 欧美性感艳星| xxxhd国产人妻xxx| av国产精品久久久久影院| 高清不卡的av网站| 咕卡用的链子| 99久国产av精品国产电影| 国产 精品1| 少妇人妻精品综合一区二区| 热re99久久精品国产66热6| 亚洲精品中文字幕在线视频| 欧美日韩亚洲高清精品| 高清毛片免费看| 97在线人人人人妻| 精品人妻一区二区三区麻豆| 老女人水多毛片| 国产色爽女视频免费观看| 久久精品aⅴ一区二区三区四区 | 老司机亚洲免费影院| 亚洲欧美一区二区三区国产| 欧美成人午夜精品| 欧美丝袜亚洲另类| 热re99久久精品国产66热6| 999精品在线视频| 男人舔女人的私密视频| 内地一区二区视频在线| 亚洲伊人色综图| 又黄又爽又刺激的免费视频.| 日韩电影二区| 免费在线观看完整版高清| 欧美激情极品国产一区二区三区 | 大码成人一级视频| 国产精品三级大全| 亚洲精品久久午夜乱码| 天堂俺去俺来也www色官网| 少妇的逼水好多| 国产精品一区二区在线观看99| av免费在线看不卡| 国精品久久久久久国模美| 女的被弄到高潮叫床怎么办| 18禁国产床啪视频网站| 日本av免费视频播放| 久久午夜综合久久蜜桃| 日韩 亚洲 欧美在线| kizo精华| 日韩人妻精品一区2区三区| 亚洲精品久久久久久婷婷小说| 亚洲美女视频黄频| 多毛熟女@视频| 男男h啪啪无遮挡| 美女视频免费永久观看网站| 国产精品国产av在线观看| 亚洲天堂av无毛| 久久久久久久精品精品| av卡一久久| 国产成人欧美| 日本黄色日本黄色录像| 视频区图区小说| 成人无遮挡网站| 内地一区二区视频在线| 欧美最新免费一区二区三区| 国产精品国产三级国产av玫瑰| 在线观看免费高清a一片| 久热久热在线精品观看| 人人妻人人爽人人添夜夜欢视频| 中文欧美无线码| 伦理电影免费视频| 中文乱码字字幕精品一区二区三区| 日韩制服丝袜自拍偷拍| 国产熟女欧美一区二区| 26uuu在线亚洲综合色| 国产免费一区二区三区四区乱码| 久久久久久久国产电影| 国产熟女欧美一区二区| 2018国产大陆天天弄谢| 中文字幕av电影在线播放| 国产精品一二三区在线看| 欧美+日韩+精品| 国产 一区精品| 午夜91福利影院| 久久狼人影院| 国产免费一级a男人的天堂| 自拍欧美九色日韩亚洲蝌蚪91| 日韩一区二区三区影片| 久久精品人人爽人人爽视色| 一本色道久久久久久精品综合| 色婷婷久久久亚洲欧美| 久久久a久久爽久久v久久| 日韩一本色道免费dvd| 国产综合精华液| 观看美女的网站| 欧美成人午夜免费资源| 精品国产一区二区三区久久久樱花| 日日撸夜夜添| 久久久久精品性色| 美女中出高潮动态图| 人体艺术视频欧美日本| 国产精品三级大全| 亚洲欧美日韩另类电影网站| 日本黄色日本黄色录像| 日韩不卡一区二区三区视频在线| av在线播放精品| 伊人久久国产一区二区| 天堂中文最新版在线下载| 狂野欧美激情性xxxx在线观看| 国产爽快片一区二区三区| 欧美xxⅹ黑人| 国产亚洲午夜精品一区二区久久| 最后的刺客免费高清国语| 成年动漫av网址| 18在线观看网站| 成人漫画全彩无遮挡| 久久久久久久国产电影| 久久鲁丝午夜福利片| 少妇被粗大猛烈的视频| 精品久久久精品久久久| 考比视频在线观看| 日韩伦理黄色片| 色婷婷久久久亚洲欧美| 国产免费又黄又爽又色| 亚洲欧美精品自产自拍| 九色成人免费人妻av| 9热在线视频观看99| 午夜福利乱码中文字幕| 久久99蜜桃精品久久| 青春草亚洲视频在线观看| 最近2019中文字幕mv第一页| 欧美激情 高清一区二区三区| 久久女婷五月综合色啪小说| 桃花免费在线播放| 国产精品嫩草影院av在线观看| 国产在线一区二区三区精| 看免费成人av毛片| 日本av免费视频播放| 搡老乐熟女国产| 九色亚洲精品在线播放| 免费看光身美女| 亚洲精品乱久久久久久| 少妇熟女欧美另类| 亚洲三级黄色毛片| 九九在线视频观看精品| 中文字幕人妻熟女乱码| 视频中文字幕在线观看| 欧美3d第一页| 欧美最新免费一区二区三区| 免费久久久久久久精品成人欧美视频 | 亚洲欧洲国产日韩| 91精品伊人久久大香线蕉| 最近2019中文字幕mv第一页| 伦精品一区二区三区| 美女主播在线视频| 国产av国产精品国产| 中文字幕亚洲精品专区| 看免费成人av毛片| 99香蕉大伊视频| 看免费av毛片| 黄色 视频免费看| 国产探花极品一区二区| 少妇高潮的动态图| 成人手机av| 成人18禁高潮啪啪吃奶动态图| 亚洲精品成人av观看孕妇| 亚洲av综合色区一区| 日韩欧美一区视频在线观看| 久久精品国产亚洲av涩爱| 国产福利在线免费观看视频| 啦啦啦在线观看免费高清www| 又黄又粗又硬又大视频| 久久精品熟女亚洲av麻豆精品| 女人久久www免费人成看片| 亚洲熟女精品中文字幕| 哪个播放器可以免费观看大片| 国产成人av激情在线播放| 亚洲精品中文字幕在线视频| 久久精品熟女亚洲av麻豆精品| 久久久国产精品麻豆| 一区在线观看完整版| 欧美日韩视频高清一区二区三区二| av女优亚洲男人天堂| 久久精品国产a三级三级三级| 久久亚洲国产成人精品v| 国产精品久久久久久久电影| 亚洲av福利一区| 女性被躁到高潮视频| 97在线人人人人妻| av.在线天堂| 精品一区在线观看国产| 又黄又粗又硬又大视频| 日日啪夜夜爽| 日韩,欧美,国产一区二区三区| 午夜av观看不卡| 黄色怎么调成土黄色| 欧美3d第一页| 午夜福利,免费看| 亚洲精品久久久久久婷婷小说| 丰满乱子伦码专区| 亚洲人与动物交配视频| 两个人免费观看高清视频| a级毛片在线看网站| 欧美另类一区| 一级黄片播放器| 国产成人午夜福利电影在线观看| 欧美精品人与动牲交sv欧美| 国产男女内射视频| a 毛片基地| 欧美变态另类bdsm刘玥| 亚洲精品日韩在线中文字幕| 欧美少妇被猛烈插入视频| 97精品久久久久久久久久精品| 少妇的逼水好多| 蜜臀久久99精品久久宅男| 精品一区在线观看国产| 人妻系列 视频| 国产成人一区二区在线| 亚洲精品乱码久久久久久按摩| 丰满乱子伦码专区| 婷婷色av中文字幕| 看免费成人av毛片| 日韩不卡一区二区三区视频在线| 国产精品久久久久久av不卡| 亚洲伊人色综图| 国产片特级美女逼逼视频| 99久久人妻综合| 国产精品久久久久久久电影| 啦啦啦中文免费视频观看日本| 久久国产精品男人的天堂亚洲 | 国产av国产精品国产| 午夜精品国产一区二区电影| 久久精品人人爽人人爽视色| 免费播放大片免费观看视频在线观看| 美女福利国产在线| 老熟女久久久| 国产精品一区二区在线观看99| 日韩制服骚丝袜av| 国产精品久久久久成人av| 国产色婷婷99| 男女边摸边吃奶| 亚洲精品国产色婷婷电影| 国产精品一二三区在线看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产深夜福利视频在线观看| 免费高清在线观看日韩| 亚洲精品色激情综合| 国产精品一区www在线观看| 18+在线观看网站| av免费观看日本| 亚洲国产最新在线播放| 免费少妇av软件| 亚洲精品自拍成人| 久久青草综合色| 国产日韩欧美亚洲二区| 成人毛片60女人毛片免费| 青春草亚洲视频在线观看| 久久久久久人人人人人| 99国产综合亚洲精品| 国产精品免费大片| a级毛色黄片| 成人免费观看视频高清| 国产一区亚洲一区在线观看| 久久久久国产精品人妻一区二区| 人妻一区二区av| 亚洲美女视频黄频| 老司机亚洲免费影院| 中国美白少妇内射xxxbb| 中文字幕人妻熟女乱码| 亚洲欧美清纯卡通| 国产午夜精品一二区理论片| 一区在线观看完整版| 性色avwww在线观看| 80岁老熟妇乱子伦牲交| 欧美人与性动交α欧美精品济南到 | 中国美白少妇内射xxxbb| 国产国语露脸激情在线看| 免费日韩欧美在线观看| 尾随美女入室| 日本vs欧美在线观看视频| 97人妻天天添夜夜摸| 在线精品无人区一区二区三| 久久久久久伊人网av| 在现免费观看毛片| 亚洲欧美成人精品一区二区| 日韩不卡一区二区三区视频在线| 精品少妇黑人巨大在线播放| 久久精品久久久久久噜噜老黄| 国产在视频线精品| 赤兔流量卡办理| 一区二区av电影网| 国产 精品1| 精品酒店卫生间| 日韩av免费高清视频| 久热久热在线精品观看| 午夜福利在线观看免费完整高清在| 丰满迷人的少妇在线观看| 日日啪夜夜爽| 婷婷色麻豆天堂久久| 大片免费播放器 马上看| 青春草视频在线免费观看| 如日韩欧美国产精品一区二区三区| 精品久久久精品久久久| 国产不卡av网站在线观看| 极品人妻少妇av视频| 少妇 在线观看| 久久久久久久久久人人人人人人| 激情五月婷婷亚洲| 国产在线一区二区三区精| 亚洲精品乱码久久久久久按摩| 热re99久久国产66热| 欧美变态另类bdsm刘玥| 精品一区二区免费观看| 99热全是精品| 国产有黄有色有爽视频| 国产综合精华液| 日韩视频在线欧美| 国产一区二区在线观看日韩| 亚洲欧美中文字幕日韩二区| av有码第一页| 色哟哟·www| 韩国精品一区二区三区 | 精品99又大又爽又粗少妇毛片| 最黄视频免费看| 女的被弄到高潮叫床怎么办| 国产成人精品无人区| 91精品伊人久久大香线蕉| 另类亚洲欧美激情| 成人国语在线视频| 妹子高潮喷水视频| 高清黄色对白视频在线免费看| 黄片播放在线免费| 国产一区有黄有色的免费视频| 久久久久视频综合| 国产精品久久久av美女十八| 各种免费的搞黄视频| 国产一区有黄有色的免费视频| 日韩三级伦理在线观看| 国产色婷婷99| 亚洲国产av影院在线观看| 综合色丁香网| 在线免费观看不下载黄p国产| 在线观看免费日韩欧美大片| 中文字幕精品免费在线观看视频 | 久久久久久伊人网av| 久久久精品区二区三区| 精品国产乱码久久久久久小说| 九九爱精品视频在线观看| 草草在线视频免费看| 少妇精品久久久久久久| 国产国拍精品亚洲av在线观看| 人妻 亚洲 视频| 欧美另类一区| 蜜桃在线观看..| 黑人欧美特级aaaaaa片| 国产一区二区三区综合在线观看 | 中国美白少妇内射xxxbb| 成人18禁高潮啪啪吃奶动态图| 国产午夜精品一二区理论片| 久久热在线av| 亚洲色图综合在线观看| av在线app专区| 又大又黄又爽视频免费| 亚洲美女搞黄在线观看| 国产黄频视频在线观看| a 毛片基地| 国产成人精品久久久久久| 丝袜美足系列| 国产欧美亚洲国产| 国产亚洲最大av| 国产成人精品一,二区| 欧美人与性动交α欧美精品济南到 | 婷婷色av中文字幕| 欧美亚洲 丝袜 人妻 在线| 熟女人妻精品中文字幕| 男女免费视频国产| 午夜福利影视在线免费观看| 亚洲av电影在线进入| 卡戴珊不雅视频在线播放| 久久精品国产亚洲av涩爱| 美女xxoo啪啪120秒动态图| 亚洲人成网站在线观看播放| 久久久国产一区二区| 一级爰片在线观看| 我的女老师完整版在线观看| 亚洲美女搞黄在线观看| 国产欧美日韩综合在线一区二区| 蜜桃国产av成人99| 两性夫妻黄色片 | 国产高清三级在线| 观看av在线不卡| 亚洲欧美日韩另类电影网站| 全区人妻精品视频| 国产av国产精品国产| 国产午夜精品一二区理论片| 七月丁香在线播放| 晚上一个人看的免费电影| 久久久亚洲精品成人影院| 夜夜骑夜夜射夜夜干| 免费大片黄手机在线观看| 欧美另类一区| 天天影视国产精品| 精品一区二区三区四区五区乱码 | 黄色怎么调成土黄色| 日韩一区二区视频免费看| 久久久久精品久久久久真实原创| 最近的中文字幕免费完整| 欧美日韩一区二区视频在线观看视频在线| 午夜免费鲁丝| 97超碰精品成人国产| xxx大片免费视频| 久久久久久人妻| 国产精品无大码| 秋霞在线观看毛片| 交换朋友夫妻互换小说| 亚洲精品国产av成人精品| av有码第一页| 精品少妇久久久久久888优播| 日韩欧美一区视频在线观看| 99久久中文字幕三级久久日本| 女的被弄到高潮叫床怎么办| 亚洲精品久久久久久婷婷小说| 色婷婷久久久亚洲欧美| 国产男女超爽视频在线观看| 亚洲 欧美一区二区三区| 国产又色又爽无遮挡免| 久久久国产一区二区| 尾随美女入室| 卡戴珊不雅视频在线播放| 永久网站在线| 99热这里只有是精品在线观看| 人人妻人人澡人人看| 伊人久久国产一区二区| 国产成人精品婷婷|