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

    基于高光譜的草坪草葉綠素含量模擬估算

    2020-08-08 07:39:50楊軍銀柳小妮王洪偉王彩玲潘冬榮
    光譜學(xué)與光譜分析 2020年8期
    關(guān)鍵詞:模型

    紀(jì) 童,王 波,楊軍銀,柳小妮*,王洪偉,王彩玲,潘冬榮,徐 君

    1. 甘肅農(nóng)業(yè)大學(xué)草業(yè)學(xué)院,甘肅 蘭州 730070 2. 草業(yè)生態(tài)系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室(甘肅農(nóng)業(yè)大學(xué)),甘肅 蘭州 730070 3. 中國(guó)人民武裝警察部隊(duì)工程大學(xué),陜西 西安 710086 4. 西安石油大學(xué)計(jì)算機(jī)學(xué)院,陜西 西安 710065 5. 甘肅省草原技術(shù)推廣總站,甘肅 蘭州 730000 6. 西安航空學(xué)院,陜西 西安 710077

    引 言

    “草坪色澤”是評(píng)價(jià)草坪優(yōu)劣的重要指標(biāo),草坪色澤的確定,主要有目測(cè)法和包括勒克斯計(jì)法、 分光光度法和植物效能分析儀-PEA在內(nèi)的實(shí)測(cè)法。目測(cè)法評(píng)定結(jié)果因人而異,受主觀影響較大。實(shí)測(cè)法雖然結(jié)果準(zhǔn)確,但方法繁瑣,在景觀尺度上進(jìn)行大范圍測(cè)量較為困難。如何對(duì)草坪草進(jìn)行快速、 無(wú)損的估算對(duì)草坪業(yè)的發(fā)展尤為重要。

    近年來(lái),隨著遙感技術(shù)的快速發(fā)展,高光譜技術(shù)已成為現(xiàn)代遙感技術(shù)的重要組成部分,如利用光譜儀進(jìn)行植物光譜識(shí)別[1],以及對(duì)植物含水量[2]、 葉綠素含量[3]、 土壤有機(jī)質(zhì)[4]、 生物量等的估算和反演研究,這些成果有效解讀了地物光譜特征規(guī)律,為遙感監(jiān)測(cè),生產(chǎn)生活提供了理論支撐與技術(shù)指導(dǎo)。

    植被指數(shù)一直是研究植被特征最簡(jiǎn)單、 最常用的方法之一。結(jié)合植被指數(shù)定量估測(cè)植被葉綠素含量成為當(dāng)今研究熱點(diǎn)之一。韓浩坤等[5]建立了糜子(PanicummiliaceumL.)冠層葉片葉綠素含量的高光譜反射率模型,使用植被指數(shù)RVI,PSNDb和GNDVI750較為準(zhǔn)確地預(yù)測(cè)了糜子冠層葉綠素含量。張瑋等[2]發(fā)現(xiàn)反射率一階微分值構(gòu)建的多元回歸方程以及修正的綠色歸一化植被指數(shù)(mGNDVI)對(duì)雷竹(Phyllostachyspraecox)葉片葉綠素的擬合效果較好。竇志國(guó)等[6]發(fā)現(xiàn)高光譜數(shù)據(jù)全波段反演能較好地估算蘆葦粉大尾蚜蟲(chóng)(Hyalopteruspruni)為害下蘆葦(Phragmitesaustralis)葉綠素相對(duì)含量。

    近年來(lái)單因素線性回歸模型與多元線性回歸模型被廣泛的應(yīng)用于光譜模型建立,如周敏姑等[7]利用一元與多元線性回歸建立了SPAD-SAVI模型,并作為無(wú)人機(jī)遙感快速、 無(wú)損監(jiān)測(cè)冬小麥葉綠素的技術(shù)手段;鄭一立等[8]采用二元回歸、 多元回歸、 偏最小二乘法和基于主成分分析的BP神經(jīng)網(wǎng)絡(luò)方法建立了竹(Bambusoideae)葉片氮含量估測(cè)模型,發(fā)現(xiàn)基于主成分分析的BP神經(jīng)網(wǎng)絡(luò)方法估測(cè)效果較好。

    目前高光譜技術(shù)測(cè)定葉綠素相對(duì)含量主要應(yīng)用于農(nóng)作物[9]、 森林和部分天然草地[10],尚未見(jiàn)針對(duì)草坪草的研究報(bào)道。本研究以3種常用的草坪草為研究對(duì)象,在草坪生長(zhǎng)旺季(4月—6月)[11],使用TYS-A3500葉綠素儀和SOC710VP光譜儀,采集不同草坪草冠層葉綠素含量和對(duì)應(yīng)的光譜數(shù)據(jù),對(duì)葉綠素含量(SPAD)進(jìn)行倒數(shù)(1/SPAD)和倒數(shù)的對(duì)數(shù)(log1/SPAD)變換[12],探尋SPAD及其相關(guān)變換值與高光譜波段擬合度較高的植被指數(shù),并構(gòu)建SPAD及其相關(guān)變換值的估算模型,旨在為使用高光譜進(jìn)行草坪質(zhì)量評(píng)價(jià)提供新的思路。

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

    1.1 材料

    試驗(yàn)材料為3個(gè)常見(jiàn)的草坪草品種:“紅象”高羊茅(Festucaarundinaceacv. Hongxiang),“百靈鳥(niǎo)”多年生黑麥草(Loliumperennecv. Bailingniao)和“肯塔基”草地早熟禾(Poapratensiscv. Kentucky)。

    1.2 方法

    試驗(yàn)于2017年4月在甘肅農(nóng)業(yè)大學(xué)培養(yǎng)室中進(jìn)行,分為建模組(n=50),檢驗(yàn)組(n=15)。選擇籽粒飽滿的草坪草種子栽種于口徑一致的花盆(20 cm×20 cm)中。土壤來(lái)自甘肅農(nóng)業(yè)大學(xué)牧草試驗(yàn)站,有機(jī)質(zhì)含量0.84%,pH 7.5,速效氮95.05 mg·kg-1,速效磷7.32 mg·kg-1,速效鉀182.8 mg·kg-1。每盆施用磷酸銨30 g·m-2后,約播種 100粒種子。培養(yǎng)期間,隔3 d稱(chēng)重澆一次水,保證水量一致,水分恒定。苗期施尿素一次,施量為純氮5 g·m-2。1個(gè)月取樣測(cè)定葉片葉綠素(SPAD)含量,同時(shí)測(cè)定對(duì)應(yīng)的光譜數(shù)據(jù)。

    1.2.1 光譜數(shù)據(jù)獲取與校正

    所用儀器:美國(guó)SOC公司SOC710VP 便攜式可見(jiàn)/紅外成像光譜儀。光譜范圍:400~1 000 nm;分辨率:4.687 5 nm;波段:128個(gè);Dynamic Range:12-bit;每行像素:696 。

    使用SOC710VP于暗室測(cè)定光譜, 測(cè)量的結(jié)果為像素值,通過(guò)SRAnalysis軟件進(jìn)行反射率標(biāo)準(zhǔn)化,以獲取原始光譜反射率。測(cè)定時(shí)采用儀器自帶鎢燈光源,鏡頭距材料上方15 cm目標(biāo)物周?chē)辛己玫耐ㄒ晽l件,高度角10°以上無(wú)近遮擋物,測(cè)量期間附近應(yīng)無(wú)運(yùn)動(dòng)物體(避免外界自然光的影響)。使用反射參考板時(shí),不應(yīng)該接觸光學(xué)面,參考板每次清潔后重新定標(biāo)。取草坪草冠層9點(diǎn)光譜數(shù)據(jù)的平均值作為原始光譜反射率。

    1.2.2 草坪草葉片葉綠素含量測(cè)定

    所用儀器:TYS-A手持葉綠素儀。測(cè)量最小間隔:小于2 s;測(cè)量范圍:0.0~99.9 SPAD;精度:±1.0 SPAD (室溫下,SPAD值介乎0~50);重復(fù)性:±0.3 SPAD單位以內(nèi)(SPAD值介乎0~50);重現(xiàn)性:±0.5 SPAD單位以內(nèi)(SPAD值介乎0~50)。測(cè)定時(shí)不必?fù)p傷植物葉片,可直接檢測(cè),工作效率更高。

    為減少試驗(yàn)誤差,每測(cè)一組光譜數(shù)據(jù)后,立即使用TYS-A手持葉綠素儀進(jìn)行草坪草葉綠素含量的測(cè)定。隨機(jī)選取葉片,在距葉尖3 cm處、 葉中、 距葉尾3 cm處各測(cè)一次葉綠素SPAD含量值,最后取平均值作為樣品SPAD值。

    1.2.3 數(shù)據(jù)處理與分析

    選擇10個(gè)植被指數(shù)[13-15](表1),采用Person相關(guān)分析方法,進(jìn)行10個(gè)植被指數(shù)和草坪草葉片的3個(gè)葉綠素指標(biāo):SPAD,1/SPAD,log(1/SPAD)的相關(guān)分析。

    表1 光譜指數(shù)Table 1 Spectral index

    10個(gè)植被指數(shù)中,GI,CI和RGI對(duì)植物自身葉綠素含量變化較敏感,EVI,RVI和VARI可以有效矯正大氣溶膠影響,NDVI670,NDVI705和MSR705對(duì)冠層結(jié)構(gòu)非常敏感,PSRI可用于植被健康的監(jiān)測(cè)與檢測(cè)。

    以葉綠素指標(biāo)為因變量,選擇相關(guān)分析中系數(shù)最大的光譜指數(shù)為自變量,運(yùn)用單因素回歸和多元逐步線性回歸法,建立草坪草葉片SPAD相關(guān)指標(biāo)值的反演模型。為確定最優(yōu)模型,用驗(yàn)證樣本對(duì)模型進(jìn)行檢驗(yàn),篩選最優(yōu)模型。

    2 結(jié)果與討論

    2.1 不同草坪草的光譜反射率(REF)

    圖1為3種草坪草在400~1 000 nm波段的光譜反射率(REF)。

    圖1 不同草坪草光譜反射率(REF)Fig.1 Spectral reflectance of different turf grasses(REF)

    由圖1可以看出,3種草坪草光譜曲線變化趨勢(shì)相似,波峰、 波谷等出現(xiàn)的位置基本一致。在可見(jiàn)光波段(400~760 nm),3種草坪草在400~530 nm波段內(nèi)反射率近乎一致且變化平穩(wěn)。在420~450 nm波段,高羊茅、 多年生黑麥草和草地早熟禾都出現(xiàn)了2個(gè)小的反射波峰,隨之又迅速降低出現(xiàn)第1個(gè)波谷。680 nm處出現(xiàn)了第2個(gè)波谷,之后在760 nm附近均形成植物所特有的“紅邊”,這種紅邊現(xiàn)象雖然整體趨勢(shì)相差不大,但是不同植物之間REF還是有區(qū)別,高羊茅、 草地早熟禾在765 nm處REF達(dá)到最大值,分別為0.92,0.91,多年生黑麥草在760 nm處REF最大(0.80)。

    在近紅光波段(780~1 000 nm)內(nèi),REF一直都保持較高水平,在895 nm處出現(xiàn)了一個(gè)微小的吸收谷,但整體趨勢(shì)是下降的,草地早熟禾的REF在730~1 000 nm波段明顯大于高羊茅和多年生黑麥草。

    2.2 葉綠素指標(biāo)與不同光譜指數(shù)的相關(guān)性

    經(jīng)Person相關(guān)性矩陣分析(圖2),植被指數(shù)PSRI,VARI,RGI,GI,ARVI,CI與葉綠素指標(biāo)SPAD極顯著(p<0.01)相關(guān),MSR705和NDVI670與SPAD顯著(p<0.05)相關(guān),而EVI和NDVI705與SPAD不相關(guān)(p>0.05)。其中VARI,RGI和PSRI與SPAD的相關(guān)系數(shù)分別為0.7**,-0.64**和-0.59**,絕對(duì)值均大于0.59,相關(guān)性較高。與葉綠素指標(biāo)1/SPAD極顯著相關(guān)的植被指數(shù)有PSRI,VARI,RGI,GI,ARVI和CI,而EVI,NDVI705,MSR705和NDVI670與1/SPAD的相關(guān)性不顯著。RGI(0.88**),PSRI(0.85**)和VARI(-0.79**)與1/SPAD相關(guān)系數(shù)的絕對(duì)值均大于0.79,說(shuō)明植被指數(shù)RGI,PSRI和VARI與1/SPAD的相關(guān)性較高。與葉綠素指標(biāo)log(1/SPAD)極顯著相關(guān)的植被指數(shù)有PSRI,VARI,RGI,GI,ARVI和NDVI670;顯著相關(guān)的植被指數(shù)是CI;而EVI,NDVI705和MSR705與log(1/SPAD)相關(guān)性不顯著。其中VARI(-0.8**),RGI(0.8**)和PSRI(0.75**)的相關(guān)系數(shù)的絕對(duì)值較高,均大于0.75。

    圖2 相關(guān)性矩陣注:圖形分為兩部分,上三角為相關(guān)系數(shù)與顯著性檢驗(yàn),下三角為散點(diǎn)圖與擬合曲線,圖中星號(hào)代表顯著性檢驗(yàn)P值,無(wú)*與*代表差異性不顯著,**與***代表差異性顯著,各變量的相關(guān)關(guān)系可在兩變量的行列交叉處找到Fig.2 Person correlation matrixNote: The graph is divided into two parts, the upper triangle is the correlation coefficient and significance test, the lower triangle is the scatter diagram and fitting curve, in which the asterisk represents the significance test P value, no * and * represent the insignificant difference, ** and *** represent the significant difference, and the correlation of each variable can be found at the row and column intersection of the two variables

    2.3 反演模型的建立

    運(yùn)用單因素回歸和多元逐步線性回歸建立的草坪草葉片SPAD相關(guān)指標(biāo)值的反演模型如表2所示。

    表2表明,利用VARI,RGI和PSRI這3個(gè)植被指數(shù)建立的葉綠素指標(biāo)反演模型的擬合系數(shù)R2均較高,其中植被指數(shù)RGI非線性擬合葉綠素SPAD和1/SPAD指標(biāo)的效果最好,R2分別為0.646和0.870;而利用VARI非線性擬合log1/SPAD指標(biāo)的效果也較好,R2=0.657。

    表2 單因素回歸模型Table 3 Single factor regression model

    圖3是植被指數(shù)與3個(gè)葉綠素指標(biāo)多元逐步線性回歸分析結(jié)果。即將10個(gè)植被指數(shù)分別代入多元逐步回歸中,得出的方程如下:

    圖3 葉綠素指標(biāo)的多元線性回歸模型注:虛線面代表的是預(yù)測(cè)值,圖中的點(diǎn)是實(shí)際值;虛線面到點(diǎn)的垂直距離代表殘差值Fig.3 SPAD multiple linear regression modelNote: The dot-line-plan represents predicted values, and the points in the graph represent actual values;The vertical distance from the dot-line-plane to the point represents the residual value

    ySPAD=47.235xVARI-23.75xMSR705+20.481

    (R2=0.532,RMSE=0.150)

    y1/SPAD=0.161xRGI+0.007xGI-0.054

    (R2=0.817,RMSE=0.023)

    ylog(1/SPAD)=0.895xRGI+0.298xMSR705-2.08

    (R2=0.680,RMSE=0.600)

    其中多元回歸方程:y1/SPAD=0.161xRGI+0.007xGI-0.054的擬合系數(shù)R2最大,RMSE最小。說(shuō)明利用植被指數(shù)RGI和GI,通過(guò)多元線性回歸,可以很好的擬合葉綠素指標(biāo)1/SPAD。

    2.4 模型精度檢驗(yàn)

    為了進(jìn)一步驗(yàn)證模型的準(zhǔn)確性和穩(wěn)定性,選擇R2大于0.7的擬合模型進(jìn)行精度檢驗(yàn),結(jié)果如圖4所示。

    由圖4可知,植被指數(shù)RGI多元線性回歸擬合1/SPAD的R2=0.928 9,而均方根誤差RMSE=0.002;VARI擬合1/SPAD的R2=0.782 3,RMSE=0.057 7,RGI擬合1/SPAD的R2=0.825 8,RMSE=0.03,PSRI擬合1/SPAD的R2=0.908 4,RMSE為0.01。以上反演模型,擬合系數(shù)都較大,且擬合誤差較小,反演精度均較高。即利用植被指數(shù)RGI,VARI和 PSRI可以很好的擬合葉綠素倒數(shù)(1/SPAD)。

    圖4 草坪草葉片1/SPAD實(shí)測(cè)值與預(yù)測(cè)值比較Fig.4 Comparison of measured and predicted SPAD values of Turfgrass blade

    3 結(jié) 論

    3.1 不同草坪草的光譜特征

    植被葉片和冠層的光譜差異與植被葉片顏色、 葉片中的葉綠素含量、 冠層結(jié)構(gòu)、 植物長(zhǎng)勢(shì)等緊密相關(guān)。

    530~650 nm波段是綠光波段,草坪草葉片均呈綠色,3種草坪草光譜曲線整體趨勢(shì)相差不大,但不同種間反射率還是有區(qū)別。在730~1 000 nm波段,草地早熟禾光譜特征更為明顯,這可能是因?yàn)椴莸卦缡旌倘~片狹長(zhǎng)且枝葉茂盛、 蓋度大,葉與葉之間透視光反射導(dǎo)致近紅外反射增強(qiáng)的結(jié)果。

    3.2 葉綠素指標(biāo)與植被指數(shù)相關(guān)性分析

    使用GNDVI能夠較好的預(yù)測(cè)草坪草葉綠素含量,3個(gè)葉綠素指標(biāo)SPAD,1/SPAD,log(1/SPAD)與選取的植被指數(shù)相關(guān)性較強(qiáng),其中RGI,PSRI和VARI與3個(gè)葉綠素指標(biāo)均呈極顯著相關(guān),相關(guān)系數(shù)絕對(duì)值均大于0.65。

    利用相關(guān)性最高植被指數(shù)RGI,PSRI和VARI,通過(guò)單因素回歸與多元線性回歸,建立草坪草葉綠素含量的估算模型,其中葉綠素1/SPAD與植被指數(shù)擬合程度最高,log(1/SPAD)次之,原值SPAD最低,與李哲等[12]的研究結(jié)果類(lèi)似。植被指數(shù)RGI在所有擬合模型中的R2(0.870)均高于其他植被指數(shù)。

    3.3 葉綠素估算模型

    綜合估算精度與預(yù)測(cè)效果,雖然多元逐步回歸模型y1/SPAD=0.161xRGI+0.007xGI-0.054的模擬精度略低于RGI(1/SPAD),但模型擬合檢驗(yàn)效果最好,預(yù)測(cè)值和實(shí)測(cè)值更接近,可據(jù)其模擬估算草坪草葉綠素含量。而多元線性模型的精度與預(yù)測(cè)效果較好,可能是因?yàn)榫C合考慮了多個(gè)植被指數(shù),提高了回歸模型預(yù)測(cè)效果。余蛟洋等[14]也發(fā)現(xiàn)了多元線性回歸模型反演獼猴桃(Actinidia Chinensis)葉片SPAD值單因素回歸模型效果更好,與本研究結(jié)果相似。

    以上試驗(yàn)結(jié)果為草坪草葉綠素的快速估算提供了理論基礎(chǔ),后期擬進(jìn)行大田草坪草葉綠素含量的反演研究,為草坪草質(zhì)量評(píng)估提供更便利的方案。

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    一级片免费观看大全| 国产高清videossex| 中文精品一卡2卡3卡4更新| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜老司机福利片| 日韩人妻精品一区2区三区| 免费高清在线观看日韩| 黄色一级大片看看| 啦啦啦 在线观看视频| 成年动漫av网址| 91成人精品电影| 在线观看一区二区三区激情| 波野结衣二区三区在线| 国产欧美日韩综合在线一区二区| 免费在线观看日本一区| 国产精品免费视频内射| 久久精品aⅴ一区二区三区四区| 少妇裸体淫交视频免费看高清 | 一级毛片黄色毛片免费观看视频| 午夜老司机福利片| 99九九在线精品视频| 精品福利观看| 精品国产一区二区三区四区第35| 亚洲av电影在线进入| 只有这里有精品99| 亚洲精品国产色婷婷电影| 夜夜骑夜夜射夜夜干| 国产深夜福利视频在线观看| 国产免费现黄频在线看| 日本一区二区免费在线视频| 美女福利国产在线| 久久精品久久久久久久性| 欧美黑人精品巨大| 亚洲av日韩在线播放| 国产麻豆69| 国产精品一区二区免费欧美 | 午夜福利乱码中文字幕| 99精国产麻豆久久婷婷| 国产成人91sexporn| 我的亚洲天堂| 中文字幕制服av| 亚洲国产日韩一区二区| 美女主播在线视频| 久久青草综合色| 欧美成人午夜精品| 国产伦理片在线播放av一区| 国产精品一国产av| 50天的宝宝边吃奶边哭怎么回事| 欧美 日韩 精品 国产| 亚洲,欧美,日韩| 久久免费观看电影| 国产精品欧美亚洲77777| 最近手机中文字幕大全| 国产精品麻豆人妻色哟哟久久| 成人亚洲精品一区在线观看| 18禁裸乳无遮挡动漫免费视频| 亚洲精品国产av成人精品| 十八禁高潮呻吟视频| 亚洲av电影在线进入| 亚洲国产av影院在线观看| 国产成人精品久久二区二区免费| 亚洲图色成人| 国产伦人伦偷精品视频| 亚洲国产欧美在线一区| 欧美日韩av久久| 男女无遮挡免费网站观看| 亚洲欧美精品自产自拍| 黑人猛操日本美女一级片| 国产精品一区二区精品视频观看| 一本一本久久a久久精品综合妖精| 欧美日韩成人在线一区二区| 一区二区三区四区激情视频| av在线app专区| 欧美日韩综合久久久久久| 亚洲伊人久久精品综合| 老司机影院成人| 欧美97在线视频| 久久99一区二区三区| 成人国语在线视频| 亚洲情色 制服丝袜| 亚洲色图综合在线观看| 亚洲av片天天在线观看| 欧美中文综合在线视频| 成人三级做爰电影| 国产精品一区二区免费欧美 | 夜夜骑夜夜射夜夜干| 一区二区三区激情视频| 无遮挡黄片免费观看| 国产亚洲欧美在线一区二区| 成人手机av| 亚洲国产中文字幕在线视频| netflix在线观看网站| 亚洲精品中文字幕在线视频| kizo精华| 最近中文字幕2019免费版| 免费观看av网站的网址| 日韩人妻精品一区2区三区| 别揉我奶头~嗯~啊~动态视频 | xxxhd国产人妻xxx| 男男h啪啪无遮挡| 国语对白做爰xxxⅹ性视频网站| av电影中文网址| 国产99久久九九免费精品| 一边摸一边做爽爽视频免费| 99re6热这里在线精品视频| 免费在线观看黄色视频的| 91精品伊人久久大香线蕉| 亚洲av成人不卡在线观看播放网 | 一区二区av电影网| 一区二区三区乱码不卡18| 国产在线观看jvid| 一级,二级,三级黄色视频| 一级毛片电影观看| 亚洲精品在线美女| 国产av精品麻豆| 91麻豆av在线| 美女视频免费永久观看网站| 99精国产麻豆久久婷婷| 制服人妻中文乱码| 久久精品国产a三级三级三级| 叶爱在线成人免费视频播放| 少妇的丰满在线观看| 9色porny在线观看| 在线观看人妻少妇| 性色av乱码一区二区三区2| 一本色道久久久久久精品综合| 免费日韩欧美在线观看| 国产又色又爽无遮挡免| 中文字幕人妻熟女乱码| 最近手机中文字幕大全| 人妻一区二区av| 18禁黄网站禁片午夜丰满| 欧美日韩国产mv在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久久久久久久大奶| 人人妻人人添人人爽欧美一区卜| 男人舔女人的私密视频| 最近最新中文字幕大全免费视频 | 啦啦啦 在线观看视频| 色婷婷av一区二区三区视频| 美女国产高潮福利片在线看| 国产一区二区三区综合在线观看| 国产男人的电影天堂91| 激情视频va一区二区三区| 国产欧美亚洲国产| 精品久久蜜臀av无| 亚洲精品国产区一区二| 国产伦人伦偷精品视频| 久久99热这里只频精品6学生| 亚洲国产精品国产精品| 亚洲国产精品一区二区三区在线| 一级毛片 在线播放| 人人妻人人添人人爽欧美一区卜| 少妇 在线观看| 十八禁高潮呻吟视频| 一二三四社区在线视频社区8| 伊人亚洲综合成人网| 国产精品亚洲av一区麻豆| 看免费av毛片| 一边摸一边做爽爽视频免费| 久久精品熟女亚洲av麻豆精品| 叶爱在线成人免费视频播放| 少妇人妻 视频| 啦啦啦中文免费视频观看日本| 国产精品国产三级国产专区5o| 国产淫语在线视频| 亚洲九九香蕉| 少妇人妻久久综合中文| 国产黄色视频一区二区在线观看| 国产在线视频一区二区| 大码成人一级视频| 欧美成人精品欧美一级黄| 少妇裸体淫交视频免费看高清 | 国产日韩一区二区三区精品不卡| 女性被躁到高潮视频| 每晚都被弄得嗷嗷叫到高潮| 午夜福利乱码中文字幕| 久久女婷五月综合色啪小说| 亚洲美女黄色视频免费看| 中文字幕最新亚洲高清| 一边摸一边抽搐一进一出视频| 成人国产一区最新在线观看 | 欧美日韩视频高清一区二区三区二| 久热爱精品视频在线9| 精品视频人人做人人爽| 大话2 男鬼变身卡| 少妇的丰满在线观看| 欧美 亚洲 国产 日韩一| 女人高潮潮喷娇喘18禁视频| 精品国产乱码久久久久久男人| 国产亚洲欧美精品永久| 黑人巨大精品欧美一区二区蜜桃| 欧美xxⅹ黑人| 天堂8中文在线网| 丝袜喷水一区| 欧美日韩福利视频一区二区| 久热爱精品视频在线9| 女人久久www免费人成看片| 久久久久国产一级毛片高清牌| 黄网站色视频无遮挡免费观看| 90打野战视频偷拍视频| 国产精品一区二区免费欧美 | 美女脱内裤让男人舔精品视频| 成人影院久久| 亚洲 国产 在线| 老鸭窝网址在线观看| 真人做人爱边吃奶动态| 欧美老熟妇乱子伦牲交| e午夜精品久久久久久久| 亚洲五月色婷婷综合| 大香蕉久久网| 久久影院123| 国产野战对白在线观看| 亚洲五月婷婷丁香| 高清欧美精品videossex| 亚洲,欧美,日韩| 国产精品人妻久久久影院| 激情视频va一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲,一卡二卡三卡| 亚洲第一av免费看| 久久久精品国产亚洲av高清涩受| 国产亚洲一区二区精品| 伊人亚洲综合成人网| 男男h啪啪无遮挡| 久久毛片免费看一区二区三区| 日本wwww免费看| 亚洲专区国产一区二区| 免费观看av网站的网址| 亚洲av片天天在线观看| 捣出白浆h1v1| 制服人妻中文乱码| 久久精品aⅴ一区二区三区四区| 日本wwww免费看| 一边摸一边做爽爽视频免费| 亚洲av欧美aⅴ国产| 国产免费视频播放在线视频| 尾随美女入室| 午夜久久久在线观看| 国产爽快片一区二区三区| 男的添女的下面高潮视频| 色播在线永久视频| 亚洲精品美女久久久久99蜜臀 | 亚洲国产日韩一区二区| 在线观看国产h片| 丁香六月天网| 操出白浆在线播放| 一本综合久久免费| 日韩视频在线欧美| 国产熟女欧美一区二区| 精品少妇久久久久久888优播| 久久国产精品男人的天堂亚洲| 在线观看免费午夜福利视频| 久久精品久久久久久久性| 国产福利在线免费观看视频| 在线观看免费高清a一片| 午夜免费男女啪啪视频观看| 国产男女内射视频| 秋霞在线观看毛片| 考比视频在线观看| 久久精品成人免费网站| 制服人妻中文乱码| 在线观看免费视频网站a站| 菩萨蛮人人尽说江南好唐韦庄| 建设人人有责人人尽责人人享有的| 中文乱码字字幕精品一区二区三区| 亚洲国产精品999| 亚洲精品乱久久久久久| 老熟女久久久| 日本午夜av视频| 日韩 亚洲 欧美在线| 99精国产麻豆久久婷婷| 丝袜美足系列| 日韩熟女老妇一区二区性免费视频| 国产片内射在线| 天天躁狠狠躁夜夜躁狠狠躁| 少妇人妻久久综合中文| 伊人久久大香线蕉亚洲五| 嫁个100分男人电影在线观看 | 久久精品亚洲av国产电影网| 手机成人av网站| 波野结衣二区三区在线| 精品高清国产在线一区| 女人久久www免费人成看片| 久久国产精品人妻蜜桃| 18禁国产床啪视频网站| 夫妻午夜视频| 一个人免费看片子| 免费不卡黄色视频| av天堂在线播放| 亚洲熟女精品中文字幕| 欧美日韩黄片免| 亚洲激情五月婷婷啪啪| 亚洲精品一卡2卡三卡4卡5卡 | 18在线观看网站| 久久久久网色| 大型av网站在线播放| 亚洲成人手机| 国产欧美日韩精品亚洲av| 成人影院久久| 国产91精品成人一区二区三区 | 各种免费的搞黄视频| 久久人人97超碰香蕉20202| 成年人黄色毛片网站| 久久国产精品影院| 亚洲国产欧美一区二区综合| 国产片特级美女逼逼视频| 国产熟女欧美一区二区| 国产在视频线精品| 人体艺术视频欧美日本| 咕卡用的链子| 久久人妻福利社区极品人妻图片 | 欧美精品人与动牲交sv欧美| 下体分泌物呈黄色| 老鸭窝网址在线观看| 久久精品亚洲熟妇少妇任你| 成人亚洲欧美一区二区av| 黄色怎么调成土黄色| 国产精品国产av在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产深夜福利视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕制服av| 日本av免费视频播放| 一级黄片播放器| 欧美精品人与动牲交sv欧美| 国产精品一区二区精品视频观看| 亚洲激情五月婷婷啪啪| 亚洲五月色婷婷综合| 老鸭窝网址在线观看| 久久久久网色| 天天躁夜夜躁狠狠久久av| 国产老妇伦熟女老妇高清| 叶爱在线成人免费视频播放| 又黄又粗又硬又大视频| 香蕉丝袜av| 午夜日韩欧美国产| 国产成人精品在线电影| 91九色精品人成在线观看| www.av在线官网国产| 看免费av毛片| 性高湖久久久久久久久免费观看| 国产成人精品久久二区二区免费| 色网站视频免费| 91精品国产国语对白视频| 亚洲久久久国产精品| 成人影院久久| 欧美黑人欧美精品刺激| 国产成人av激情在线播放| kizo精华| 国产一区亚洲一区在线观看| 亚洲精品国产av蜜桃| 婷婷色麻豆天堂久久| 一区二区三区乱码不卡18| 午夜免费男女啪啪视频观看| 伊人亚洲综合成人网| 激情五月婷婷亚洲| 亚洲少妇的诱惑av| 国产一区二区三区av在线| 欧美日韩精品网址| 成人亚洲欧美一区二区av| 欧美中文综合在线视频| www.自偷自拍.com| 1024视频免费在线观看| 国产在视频线精品| 日韩,欧美,国产一区二区三区| 国产男女内射视频| 亚洲国产精品国产精品| 日本av免费视频播放| 九色亚洲精品在线播放| 美女脱内裤让男人舔精品视频| 久久人人爽av亚洲精品天堂| 少妇精品久久久久久久| 天天操日日干夜夜撸| 91老司机精品| 巨乳人妻的诱惑在线观看| 国产国语露脸激情在线看| 国产精品免费视频内射| 性色av乱码一区二区三区2| 91精品国产国语对白视频| kizo精华| 大码成人一级视频| 亚洲精品国产av蜜桃| 国产亚洲欧美在线一区二区| 午夜久久久在线观看| 丝瓜视频免费看黄片| 国产精品国产三级专区第一集| 2021少妇久久久久久久久久久| 国产欧美日韩精品亚洲av| 夫妻午夜视频| 日本91视频免费播放| 在线观看免费日韩欧美大片| 成年人午夜在线观看视频| 性高湖久久久久久久久免费观看| 国产人伦9x9x在线观看| 国产精品久久久av美女十八| 少妇人妻 视频| 热re99久久精品国产66热6| 久久精品国产亚洲av高清一级| 热re99久久国产66热| 国产xxxxx性猛交| 性少妇av在线| 老汉色∧v一级毛片| 亚洲人成电影观看| 高清欧美精品videossex| 成人免费观看视频高清| 免费在线观看日本一区| 18禁国产床啪视频网站| 性少妇av在线| 老汉色∧v一级毛片| 最近中文字幕2019免费版| 两个人看的免费小视频| 免费少妇av软件| 午夜老司机福利片| 日韩一区二区三区影片| 亚洲激情五月婷婷啪啪| 丝袜脚勾引网站| 肉色欧美久久久久久久蜜桃| 黄片播放在线免费| 成人午夜精彩视频在线观看| 国产成人一区二区三区免费视频网站 | 中文字幕av电影在线播放| 欧美精品高潮呻吟av久久| 欧美日韩视频精品一区| 亚洲av片天天在线观看| av有码第一页| 中文精品一卡2卡3卡4更新| 国产一区有黄有色的免费视频| 波多野结衣一区麻豆| 啦啦啦啦在线视频资源| 久久久久久久久久久久大奶| 亚洲自偷自拍图片 自拍| 十八禁网站网址无遮挡| 国产精品一区二区在线不卡| 赤兔流量卡办理| www.精华液| 国产又爽黄色视频| 国产97色在线日韩免费| 国产激情久久老熟女| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲国产最新在线播放| 久久精品aⅴ一区二区三区四区| 高清黄色对白视频在线免费看| 日日爽夜夜爽网站| 欧美+亚洲+日韩+国产| 国产在线一区二区三区精| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品一国产av| 久久精品久久久久久噜噜老黄| 精品卡一卡二卡四卡免费| 亚洲国产成人一精品久久久| 超色免费av| av电影中文网址| 人人妻,人人澡人人爽秒播 | av又黄又爽大尺度在线免费看| 久久精品久久精品一区二区三区| 日韩精品免费视频一区二区三区| 欧美亚洲 丝袜 人妻 在线| 男人爽女人下面视频在线观看| 99精国产麻豆久久婷婷| 美女扒开内裤让男人捅视频| 伊人久久大香线蕉亚洲五| 1024视频免费在线观看| 国产在视频线精品| 成人亚洲欧美一区二区av| 亚洲,一卡二卡三卡| 啦啦啦啦在线视频资源| 日韩大码丰满熟妇| av片东京热男人的天堂| 日本av手机在线免费观看| 自线自在国产av| 欧美日韩一级在线毛片| 激情视频va一区二区三区| 日韩大码丰满熟妇| 一二三四在线观看免费中文在| 成人免费观看视频高清| 久久午夜综合久久蜜桃| 无遮挡黄片免费观看| av片东京热男人的天堂| 热re99久久精品国产66热6| 欧美在线一区亚洲| 久久影院123| 天天躁夜夜躁狠狠躁躁| 交换朋友夫妻互换小说| 91九色精品人成在线观看| 老司机影院毛片| 咕卡用的链子| 亚洲人成77777在线视频| 国产99久久九九免费精品| 天天躁夜夜躁狠狠躁躁| 九草在线视频观看| 青春草亚洲视频在线观看| h视频一区二区三区| 国产精品亚洲av一区麻豆| 亚洲欧美中文字幕日韩二区| 国产成人一区二区在线| 操出白浆在线播放| 曰老女人黄片| 亚洲国产精品国产精品| 国产三级黄色录像| 国产人伦9x9x在线观看| 国产亚洲av片在线观看秒播厂| 波多野结衣一区麻豆| a 毛片基地| 90打野战视频偷拍视频| 精品国产国语对白av| 国产精品香港三级国产av潘金莲 | 男男h啪啪无遮挡| 国产精品.久久久| 久久国产亚洲av麻豆专区| 国产精品久久久av美女十八| 成年女人毛片免费观看观看9 | 午夜福利视频精品| 美女午夜性视频免费| 黄色片一级片一级黄色片| www.自偷自拍.com| 人人妻人人澡人人看| 日韩一区二区三区影片| 这个男人来自地球电影免费观看| 侵犯人妻中文字幕一二三四区| 色精品久久人妻99蜜桃| 亚洲av片天天在线观看| 性高湖久久久久久久久免费观看| 国产1区2区3区精品| 亚洲成av片中文字幕在线观看| 少妇人妻久久综合中文| 一区二区三区乱码不卡18| 成年人免费黄色播放视频| 高清av免费在线| 丝袜喷水一区| 极品人妻少妇av视频| 乱人伦中国视频| 精品久久久久久电影网| 亚洲国产精品国产精品| 18禁裸乳无遮挡动漫免费视频| 悠悠久久av| 色网站视频免费| www日本在线高清视频| 国产成人精品无人区| 91九色精品人成在线观看| 欧美xxⅹ黑人| 免费在线观看视频国产中文字幕亚洲 | 婷婷色综合www| 在线精品无人区一区二区三| 亚洲情色 制服丝袜| 亚洲精品一二三| 国产片内射在线| 久久亚洲精品不卡| 18禁裸乳无遮挡动漫免费视频| 国产一区有黄有色的免费视频| 免费不卡黄色视频| 国产真人三级小视频在线观看| 女警被强在线播放| 亚洲av成人不卡在线观看播放网 | 中文欧美无线码| 少妇猛男粗大的猛烈进出视频| 精品一品国产午夜福利视频| 久久久久久亚洲精品国产蜜桃av| 美女视频免费永久观看网站| 一级,二级,三级黄色视频| 亚洲精品国产色婷婷电影| 亚洲国产精品999| 伦理电影免费视频| 欧美av亚洲av综合av国产av| 久久久久国产一级毛片高清牌| 亚洲国产毛片av蜜桃av| 国产精品免费大片| 黑丝袜美女国产一区| 中文字幕av电影在线播放| 国产成人精品在线电影| 人人澡人人妻人| 精品亚洲乱码少妇综合久久| 女警被强在线播放| 好男人电影高清在线观看| 中文字幕色久视频| 夜夜骑夜夜射夜夜干| 99精国产麻豆久久婷婷| 亚洲午夜精品一区,二区,三区| 美女脱内裤让男人舔精品视频| 国产亚洲欧美精品永久| 欧美97在线视频| 亚洲欧洲日产国产| 亚洲第一av免费看| 少妇人妻久久综合中文| 久久亚洲精品不卡| √禁漫天堂资源中文www| 国产精品久久久av美女十八| 久久这里只有精品19| 亚洲精品国产av蜜桃| 成年女人毛片免费观看观看9 | 精品欧美一区二区三区在线| 国产亚洲av片在线观看秒播厂| 亚洲天堂av无毛| 天堂俺去俺来也www色官网| 波多野结衣一区麻豆| 亚洲国产精品国产精品| 亚洲精品自拍成人| 国精品久久久久久国模美| 美女脱内裤让男人舔精品视频| 日日夜夜操网爽| 国产高清不卡午夜福利| 亚洲中文日韩欧美视频| 女性生殖器流出的白浆| 久久久久久久精品精品| 亚洲第一青青草原| 久久免费观看电影| 亚洲一区中文字幕在线| 国产精品成人在线| 欧美成人精品欧美一级黄| 18禁裸乳无遮挡动漫免费视频| 赤兔流量卡办理| 中文精品一卡2卡3卡4更新| 午夜影院在线不卡| 久久久精品国产亚洲av高清涩受| 久热爱精品视频在线9|