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

    基于SPA-MLR方法的土壤含水量光譜預(yù)測(cè)模型研究

    2018-07-04 08:02:16賈學(xué)勤馮美臣楊武德武改紅
    關(guān)鍵詞:微分反射率波段

    賈學(xué)勤,馮美臣,楊武德,王 超,孫 慧,武改紅,張 松

    (山西農(nóng)業(yè)大學(xué)農(nóng)學(xué)院, 山西 太谷 030801)

    土壤水分對(duì)植物的生長(zhǎng)發(fā)育具有重要的影響,與植物和外界環(huán)境間的物質(zhì)和能量交換密切相關(guān),也是植物生長(zhǎng)狀況監(jiān)測(cè)的重要指標(biāo)[1],因此土壤含水量的快速準(zhǔn)確測(cè)定對(duì)于農(nóng)業(yè)生產(chǎn)具有重要意義。傳統(tǒng)的測(cè)定土壤含水量的方法主要有烘干法、中子儀法以及時(shí)間域反射儀法(TDR)等,但耗時(shí)費(fèi)力,且測(cè)定范圍小,難以滿足現(xiàn)代農(nóng)業(yè)對(duì)大面積土壤含水量準(zhǔn)確及時(shí)測(cè)量的需求[2]。

    高光譜遙感監(jiān)測(cè)具有快速、準(zhǔn)確、無(wú)破壞的特點(diǎn),在監(jiān)測(cè)土壤含水量方面具有重要作用。隨著高光譜遙感技術(shù)的快速發(fā)展,光譜分辨率不斷提高,光譜信息極大豐富的同時(shí)也增大了背景和噪音的影響,進(jìn)而增加準(zhǔn)確提取特征光譜信息的難度,影響監(jiān)測(cè)模型的精度[3]。研究發(fā)現(xiàn)光譜數(shù)據(jù)變換在消除背景和噪音干擾,降低內(nèi)外部因子的影響,放大光譜吸收特征方面具有較好的效果,有利于光譜信息的挖掘和提取,進(jìn)而提高監(jiān)測(cè)模型的精度和穩(wěn)定性[4,5]。Kemper等[6]對(duì)土壤原始光譜反射率進(jìn)行一階微分處理,提高了土壤重金屬、粘粒含量等的監(jiān)測(cè)精度。何挺等[7]研究表明反射率對(duì)數(shù)一階微分變換,提高了土壤光譜與含水量的相關(guān)性。伊葉彪等[8]通過(guò)土壤光譜去包絡(luò)線后進(jìn)行一階微分,預(yù)測(cè)含水量R2達(dá)到了0.855。謝伯承等[9]利用反射率倒數(shù)的對(duì)數(shù)(A)、A值的一階微分和相對(duì)反射率反演土壤含水量,得出A的一階微分反演效果最佳,且土壤含水量與光譜反射率呈明顯非線性關(guān)系。姚艷敏等[10]對(duì)土壤光譜反射率進(jìn)行光譜變換,預(yù)測(cè)黑土土壤含水量,結(jié)果表明反射率進(jìn)行對(duì)數(shù)一階微分變換后建立的模型最好,R2達(dá)到了0.931。張俊華等[11]研究得出高光譜預(yù)測(cè)龜裂堿土土壤含水量最優(yōu)模型為使用反射率一階微分的冪函數(shù)模型,R2高達(dá)0.9447。目前,光譜變換已廣泛應(yīng)用于光譜數(shù)據(jù)處理中,但是不同研究中各種變換方法效果不一。本文以此為切入點(diǎn),在前人研究的基礎(chǔ)上,首先對(duì)原始光譜進(jìn)行11種常規(guī)變換,然后選擇連續(xù)投影算法(SPA)選擇特征波段,結(jié)合多元線性回歸(MLR)建立土壤含水量光譜監(jiān)測(cè)模型。旨在利用SPA-MLR實(shí)現(xiàn)不同光譜變換方法監(jiān)測(cè)土壤含水量,并分析比較不同變換處理后各SPA-MLR回歸模型的監(jiān)測(cè)精度,進(jìn)而選擇出土壤含水量高光譜監(jiān)測(cè)的最佳模型。

    1 材料與方法

    1.1 土壤樣品采集與處理

    本試驗(yàn)以山西省晉中市太谷縣山西農(nóng)業(yè)大學(xué)試驗(yàn)站的麥田土壤為研究對(duì)象,土壤質(zhì)地為黃土母質(zhì)發(fā)育而成的石灰性褐土,取表層土壤(0~20cm)為試驗(yàn)樣品。土壤樣品自然風(fēng)干后剔除雜質(zhì),再研磨過(guò)2mm孔徑篩后混勻。稱取72份,每份100g,放于金屬盒中,編號(hào)為1~72,然后進(jìn)行土壤含水量調(diào)配,具體方法如下:在36h內(nèi),每隔0.5h取一份稱量好的土樣,向土壤中緩慢注水,使土壤達(dá)到飽和狀態(tài),注水后室溫下使其自然蒸發(fā)。待最后一份土樣注水完成、表面自由水消失后對(duì)所有土樣進(jìn)行光譜和含水量測(cè)定,從而得到72份不同含水量的土壤樣本。其中土壤含水量用烘干法測(cè)得。

    1.2 土壤光譜采集及預(yù)處理

    本試驗(yàn)采用美國(guó)ASD公司的Field Spec Pro FR 型便攜式光譜儀測(cè)定土壤光譜,波段范圍為350~2500nm,在350~1000nm區(qū)間采樣間隔為1.4nm,在1000~2500nm區(qū)間采樣間隔為2nm。土壤光譜測(cè)定在暗室內(nèi)進(jìn)行,光源為50W鹵素?zé)?。光譜測(cè)量時(shí)光照方向與垂直方位夾角15°,光源距土壤表層40cm,采用10°視場(chǎng)角 ,探頭位于土壤表面上方20cm處。光譜采集前進(jìn)行標(biāo)準(zhǔn)白板校正,光譜采集時(shí)水平轉(zhuǎn)動(dòng)樣品5次,每次轉(zhuǎn)動(dòng)60°,每個(gè)方向采集10條光譜曲線,共獲得60個(gè)光譜數(shù)據(jù),剔除異常值后將采集的光譜數(shù)據(jù)平均作為該土樣反射率光譜。為消除噪聲影響,本研究采用Savitzky-Golay3點(diǎn)平滑對(duì)光譜進(jìn)行平滑去噪處理。由于350~400nm,2450~2500nm噪聲影響較大,本研究選擇400~2450nm波段范圍進(jìn)行研究。

    1.3 光譜變換算法

    本文在前人研究基礎(chǔ)上選取11種典型的數(shù)學(xué)變換算法對(duì)光譜數(shù)據(jù)進(jìn)行變換處理,如表1所示。

    1.4 高光譜模型建立及模型評(píng)價(jià)

    連續(xù)投影算法(SPA)可以有效降低光譜數(shù)據(jù)共線性,與偏最小二乘法(PLS)相比,可以最大程度地降低信息重疊度,有效地降低光譜模型復(fù)雜度[12]。且可以簡(jiǎn)化校正模型,縮短校正時(shí)間,已經(jīng)越來(lái)越多地應(yīng)用于特征波段的選取[13]。此外,還可以改善多元線性回歸模型(MLR)的建模條件。SPA-MLR回歸模型可以減小模型擬合復(fù)雜度,加快運(yùn)算速度[14]。因此本研究選用SPA提取特征波段,結(jié)合MLR建立土壤含水量高光譜預(yù)測(cè)模型。

    表1 11種光譜數(shù)據(jù)變換算法信息

    注:R為波段λ(nm)處的反射率值。

    Note:Ris the spectral reflectance value at the wavelengths ofλ(nm).

    為確保土壤含水量模型穩(wěn)定性和精度,將得到的土壤含水量樣本(72個(gè))隨機(jī)分為兩部分,選擇48個(gè)樣本用于構(gòu)建模型,剩余24個(gè)樣本用于模型驗(yàn)證。本研究采用相對(duì)預(yù)測(cè)偏差(RPD),決定系數(shù)(R2)以及均方根誤差(RMSE)作為模型預(yù)測(cè)精度的評(píng)價(jià)標(biāo)準(zhǔn)。其中,R2越接近于1,RMSE越小,表明模型具有高的預(yù)測(cè)精度和穩(wěn)定性。當(dāng)RPD大于2時(shí),模型預(yù)測(cè)效果準(zhǔn)確;當(dāng)RPD大于1.4且小于2,模型預(yù)測(cè)效果一般可以粗略估計(jì);當(dāng)RPD小于1.4,模型無(wú)法進(jìn)行準(zhǔn)確預(yù)測(cè)[15]。

    1.5 數(shù)據(jù)分析

    本研究使用ViewSpec Pro光譜分析軟件進(jìn)行預(yù)處理,Excel 2007和Unscramb 10.3軟件進(jìn)行光譜、變換處理,Matlab 7.0和SPSS軟件進(jìn)行特征波段提取以及模型建立,Origin 8軟件制圖。

    2 結(jié)果與分析

    2.1 不同含水量土壤光譜特征

    圖1為不同含水量土壤光譜反射率曲線圖,由圖可以看出在400~2450nm范圍內(nèi),隨著土壤含水量增加光譜反射率不斷降低。在可見(jiàn)光波段,光譜反射率隨著土壤含水量的變化幅度明顯小于近紅外波段。在可見(jiàn)光-近紅外波段(380~780nm),隨著波長(zhǎng)增加反射率有升高的趨勢(shì),但是升高斜率不同;在1400nm、1900nm附近有兩個(gè)明顯的水分吸收谷。

    圖1 不同含水量土壤的反射率曲線Fig.1 Spectral curves under of different soil moisture content

    2.2 基于SPA方法和不同光譜變化方法提取的敏感波段

    表2為光譜變換處理后基于SPA提取的土壤含水量最優(yōu)波段個(gè)數(shù)和具體波段。由表2可知,最優(yōu)波段個(gè)數(shù)都為3~5個(gè),但是使用不同變換算法處理后光譜最優(yōu)波段個(gè)數(shù)有一定差異,微分處理增加了最優(yōu)波段個(gè)數(shù)。對(duì)不同變換算法處理后提取的光譜特征波段進(jìn)行分析可知,提取的土壤含水量特征波段分布比較分散,特征波段在近紅外、短波紅外區(qū)域都有分布。不同變換算法處理后的特征波段有一定差異。

    表2 基于SPA和不同光譜變化方法的重要波段的提取

    2.3 基于SPA-MLR方法的土壤含水量光譜預(yù)測(cè)模型建立

    以原始光譜以及11種變換后基于SPA提取的特征波段為自變量,以土壤含水量為因變量,利用MLR建立的土壤含水量高光譜預(yù)測(cè)模型,如表3所示?;?1種變換處理建立的SPA-MLR模型監(jiān)測(cè)精度,除了T3比較低外,其余變換均優(yōu)于原始光譜所建的模型,表明原始光譜經(jīng)一定數(shù)學(xué)變換處理可以提高土壤含水量高光譜監(jiān)測(cè)精度。T6、T7、T8、T9、T10、T11的監(jiān)測(cè)模型精度均高于對(duì)應(yīng)的一階微分處理前T0、T1、T2、T3、T4、T5。表明一階微分變換能提高土壤含水量模型的精度。在所有數(shù)學(xué)變換中,T8的校正模型的監(jiān)測(cè)效果最好,其決定系數(shù)R2和RPD最大,分別達(dá)到了0.957和4.74,RMSE最小,為2.16。T11次之,其R2和RPD分別達(dá)到了0.944和4.15,RMSE為2.55。T3的校正模型監(jiān)測(cè)效果最差,R2和RPD都最低,分別為0.860和2.31,RMSE最大,為6.18。

    表3不同土壤含水量光譜預(yù)測(cè)模型

    Tab. 3 Predictive models of different soil moisture content

    變換算法Algorithms模型ModelsR2RMSE/%RPDT0y=417.660R401-415.185R1888+271.944R1459+19.3910.8873.632.79T1y=0.694R403+0.100R418-13.925R1708+6.812R1898+13.7530.9003.423.03T2y=312.330R1415-304.339R1887+13.2570.9183.103.37T3y=225.499R1465-178.004R1052-67.516R2257-27.9320.8606.182.31T4y=-332.373R1934-29.668R408+452.190R1007-62.9350.9152.484.36T5y=46.633R1898-79.946R1708+5.419R403+37.6670.9202.893.65T6y=108790.163R1548-92439.173R1690-10313.521R2311-0.6110.9422.614.02T7y=1720.981R1302-1389.581R828+253.218R2171-1976.432R1122-612.636R2031+1.2500.9142.793.78T8y=629.471R666+8158.075R1014-12587.075R1039-31671.078R1375+34808.678R1707-2.1520.9572.164.94T9y=14294.340R1384-3194.977R2322+8246.352R1681-1491.601R1125-12.9010.9436.114.74T10y=-123683.801R1439+11675.161R1648-18882.797R1043-9774.030R1822-6604.232R1143+19.1850.9353.943.82T11y=-9208.900R1522-4388.004R1796-2637.737R1039+3034.038R2202+137.782R434-2.1220.9442.554.15

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

    利用驗(yàn)證集樣本對(duì)土壤含水量預(yù)測(cè)模型分別進(jìn)行驗(yàn)證。原始光譜以及11種變換后基于SPA-MLR建立的土壤含水量監(jiān)測(cè)模型驗(yàn)證結(jié)果如表4所示。從表中可看出,所有數(shù)學(xué)變換處理后模型的驗(yàn)證效果均優(yōu)于原始光譜建立的模型,其中T8變換模型精度最高。為了更直觀地展示模型精確度,以T8為例制作預(yù)測(cè)值與實(shí)測(cè)值1∶1關(guān)系圖(圖2)。由圖可以看出,T8處理模型的擬合效果較好,其驗(yàn)證集R2最大為0.903,RMSE最小為3.41,RPD較大為2.95。

    表4 模型驗(yàn)證結(jié)果

    圖2 土壤含水量實(shí)測(cè)值與估測(cè)值1∶1圖Fig.2 The correlation between predicted values and measured values of soil water content

    3 討 論

    光譜數(shù)據(jù)變換可以有效地消除土壤光譜噪音和背景的影響,放大特征光譜信息。本研究選擇11種典型的變換算法對(duì)土壤光譜進(jìn)行變換,然后使用連續(xù)投影算法(SPA)提取特征波段,結(jié)合MLR建立土壤含水量光譜監(jiān)測(cè)模型,比較選擇土壤含水量光譜監(jiān)測(cè)最佳模型。有研究表明,特征波段選取個(gè)數(shù)為3~5個(gè)最佳,這樣不但可以準(zhǔn)確地表達(dá)出待測(cè)指標(biāo)的信息,而且可以有效避免過(guò)擬合現(xiàn)象的出現(xiàn)[16]。本研究利用SPA提取的特征波段個(gè)數(shù)均為3~5個(gè),所以不用再對(duì)特征波段進(jìn)一步篩選。

    從本研究建立的土壤含水量的校正模型及其驗(yàn)證結(jié)果可知,各變換處理后建立的模型都有較高的預(yù)測(cè)精度,其R2都高于0.7。其中除了T3外,各變換處理后建立的土壤含水量光譜預(yù)測(cè)模型都優(yōu)于原始光譜T0,表明數(shù)學(xué)變換能提高土壤含水量模型的預(yù)測(cè)精度。這主要是由于變換處理有效地減小了光譜信息中噪音和背景的影響,放大了特征光譜信息,有利于光譜信息的提取和挖掘[17]。目前,在作物冠層和土壤光譜研究中,已經(jīng)證明一階微分處理能有效提高模型的預(yù)測(cè)精度[18,19],在本研究中光譜一階微分變換建立的模型表現(xiàn)也比較好,僅次于反射率對(duì)數(shù)一階微分變換T8。由于土壤光譜受土壤質(zhì)地、有機(jī)質(zhì)含量、土壤顆粒大小等的影響較大,不同土壤的有機(jī)質(zhì)含量和質(zhì)地差別較大[20]。由于本試驗(yàn)是在單一土壤質(zhì)地和有機(jī)質(zhì)含量背景下完成的,因此反射率對(duì)數(shù)一階微分T8變換能否適用于其它背景條件下(不同質(zhì)地、有機(jī)質(zhì)含量等)的土壤光譜還有待進(jìn)一步的研究。

    4 結(jié) 論

    本文以人工模擬72個(gè)不同含水量的土壤樣本為研究對(duì)象,測(cè)定室內(nèi)光譜以及對(duì)應(yīng)的土壤含水量,基于11種常見(jiàn)的數(shù)學(xué)變換算法,建立土壤含水量的SPA-MLR光譜監(jiān)測(cè)模型,綜合比較分析模型校正集和驗(yàn)證集模型精度評(píng)價(jià)指標(biāo)(R2、RMSE、RPD)得出:各光譜變換建立的模型都有較好的預(yù)測(cè)精度,其中,反射率對(duì)數(shù)的一階微分T8建立的模型表現(xiàn)最好,校正集(R2=0.957,RMSE=2.16,RPD=4.74)和驗(yàn)證集(R2=0.903,RMSE=3.41,RPD=2.95)表現(xiàn)均最好。因此反射率對(duì)數(shù)一階微分變換處理后建立的SPA-MLR模型為土壤含水量高光譜監(jiān)測(cè)最佳模型。

    參考文獻(xiàn):

    [1] Yin Z, Lei T, Yan Q, et al. A near-infrared reflectance sensor for soil surface moisture measurement[J]. Computers and Electronics in Agriculture, 2013, 99: 101-107.

    [2] 孫越君, 鄭小坡, 秦其明, 等. 不同質(zhì)量含水量的土壤反射率光譜模擬模型[J]. 光譜學(xué)與光譜分析, 2015,35(8): 2236-2240.

    [3] Stenberg B, Rossel R A V, Mouazen A M,et al. Chapter five-visible and near infrared spectroscopy in soil science[J]. Advances in Agronomy, 2010,107: 163-215.

    [4] 王磊, 白由路, 盧艷麗, 等. 光譜數(shù)據(jù)變換對(duì)玉米氮素含量反演精度的影響[J]. 遙感技術(shù)與應(yīng)用, 2011,26(02): 220-225.

    [5] 喬星星, 馮美臣, 楊武德, 等. 變換光譜數(shù)據(jù)對(duì)土壤氮素PLSR模型的影響研究[J]. 地球信息科學(xué)學(xué)報(bào), 2016,18(8): 1123-1132.

    [6] Kemper T, Sommer S. Estimate of heavy metal contamination in soils after a mining accident using reflectance spectroscopy[J]. Environmental Science & Technology, 2002,36(12): 2742-2747.

    [7] 何挺, 王靜, 程燁, 等. 土壤水分光譜特征研究[J]. 土壤學(xué)報(bào), 2006,43(6): 1027-1032.

    [8] 尹業(yè)彪, 李霞, 趙釗, 等.沙質(zhì)土壤含水率高光譜預(yù)測(cè)模型建立及分析[J]. 遙感技術(shù)與應(yīng)用, 2011,26(3): 355-359.

    [9] 謝伯承. 基于高光譜遙感不同發(fā)生層土壤的光譜信息的提取研究[D]. 楊凌:西北農(nóng)林科技大學(xué),2004.

    [10] 姚艷敏, 魏娜, 唐鵬欽, 等. 黑土土壤水分高光譜特征及反演模型[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2011,27(8): 95-100.

    [11] 張俊華, 賈科利. 典型龜裂堿土土壤水分光譜特征及預(yù)測(cè)[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2015,26(3): 884-890.

    [12] Araújo M C U, Saldanha T C B, Galvao R K H, et al. The successive projections algorithm for variable selection in spectroscopic multicomponent analysis[J]. Chemometrics and Intelligent Laboratory Systems, 2001,57(2): 65-73.

    [13] 陳斌, 孟祥龍, 王豪. 連續(xù)投影算法在近紅外光譜校正模型優(yōu)化中的應(yīng)用[J]. 分析測(cè)試學(xué)報(bào), 2007,26(1): 66-69.

    [14] Galv?o R K H, Araújo M C U, Silva E C,et al. Cross-validation for the selection of spectral variables using the successive rojections algorithm[J]. Journal of the Brazilian Chemical Society, 2007,18(8): 1580-1584.

    [15] 劉秀英, 王力, 宋榮杰, 等. 黃綿土風(fēng)干過(guò)程中土壤含水率的光譜預(yù)測(cè)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2015,46(4): 266-272.

    [16] 高洪智, 盧啟鵬, 丁海泉, 等. 基于連續(xù)投影算法的土壤總氮近紅外特征波長(zhǎng)的選取[J]. 光譜學(xué)與光譜分析, 2009,29(11): 2951-2954.

    [17] 魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京: 中國(guó)農(nóng)業(yè)科技出版社,2000.

    [18] 王強(qiáng), 易秋香, 包安明, 等. 基于高光譜反射率的棉花冠層葉綠素密度估算[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2012,28(15): 125-132.

    [19] Brown D J, Shepherd K D, Walsh M G,et al. Global soil characterization with VNIR diffuse reflectance spectroscopy[J]. Geoderma, 2006,132(3): 273-290.

    [20] Muller E, Decamps H. Modeling soil moisture reflectance[J]. Remote Sensing of Environment, 2001,76(2): 173-180.

    猜你喜歡
    微分反射率波段
    春日暖陽(yáng)
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    具有顏色恒常性的光譜反射率重建
    擬微分算子在Hp(ω)上的有界性
    上下解反向的脈沖微分包含解的存在性
    化學(xué)腐蝕硅表面結(jié)構(gòu)反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    借助微分探求連續(xù)函數(shù)的極值點(diǎn)
    M87的多波段輻射過(guò)程及其能譜擬合
    日常維護(hù)對(duì)L 波段雷達(dá)的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    色吧在线观看| 听说在线观看完整版免费高清| 久久97久久精品| 99热这里只有精品一区| 亚州av有码| 亚洲精品国产成人久久av| 国产黄片视频在线免费观看| av一本久久久久| 日韩伦理黄色片| 日日摸夜夜添夜夜添av毛片| 秋霞在线观看毛片| 黄色怎么调成土黄色| 久久久精品94久久精品| 毛片女人毛片| 嘟嘟电影网在线观看| 国产黄色视频一区二区在线观看| 18禁在线无遮挡免费观看视频| 亚洲激情五月婷婷啪啪| 国产又色又爽无遮挡免| 真实男女啪啪啪动态图| 亚洲精品456在线播放app| 纵有疾风起免费观看全集完整版| 嘟嘟电影网在线观看| 久久6这里有精品| 少妇人妻一区二区三区视频| 大香蕉久久网| 久久国产乱子免费精品| 我要看日韩黄色一级片| 91久久精品国产一区二区成人| 久久久精品欧美日韩精品| 欧美日韩精品成人综合77777| 欧美xxⅹ黑人| 精品少妇黑人巨大在线播放| 赤兔流量卡办理| 秋霞伦理黄片| 亚洲精品乱码久久久v下载方式| 精品亚洲乱码少妇综合久久| 夜夜看夜夜爽夜夜摸| 精品国产三级普通话版| 王馨瑶露胸无遮挡在线观看| 国产大屁股一区二区在线视频| 老司机影院毛片| 看非洲黑人一级黄片| 国产高清国产精品国产三级 | 日韩在线高清观看一区二区三区| 丰满少妇做爰视频| av免费在线看不卡| 又大又黄又爽视频免费| 97人妻精品一区二区三区麻豆| 亚洲精品视频女| 我的女老师完整版在线观看| 日韩亚洲欧美综合| 小蜜桃在线观看免费完整版高清| 制服丝袜香蕉在线| 日韩欧美精品v在线| 国产 一区精品| 日韩国内少妇激情av| 直男gayav资源| 天堂网av新在线| 亚洲精品国产av蜜桃| 亚洲精品视频女| 免费大片18禁| 少妇被粗大猛烈的视频| 国产毛片在线视频| 狂野欧美白嫩少妇大欣赏| 国产精品国产三级国产av玫瑰| 国产 精品1| 免费大片18禁| 国产亚洲5aaaaa淫片| 久久女婷五月综合色啪小说 | 黄色视频在线播放观看不卡| 最后的刺客免费高清国语| 一级毛片久久久久久久久女| 国产精品国产三级国产av玫瑰| 国产成人精品福利久久| 中文字幕人妻熟人妻熟丝袜美| 噜噜噜噜噜久久久久久91| 哪个播放器可以免费观看大片| 天堂网av新在线| 久久精品国产亚洲av天美| 人人妻人人看人人澡| 肉色欧美久久久久久久蜜桃 | 在线观看免费高清a一片| 中文乱码字字幕精品一区二区三区| 蜜桃亚洲精品一区二区三区| 亚洲久久久久久中文字幕| 国产一级毛片在线| 六月丁香七月| 国产老妇女一区| videos熟女内射| 国产乱来视频区| 国产成人免费无遮挡视频| 亚洲最大成人av| 丝袜美腿在线中文| 精品人妻熟女av久视频| 成人亚洲精品av一区二区| 尤物成人国产欧美一区二区三区| 美女被艹到高潮喷水动态| 亚洲人成网站在线播| 欧美成人精品欧美一级黄| 在现免费观看毛片| 又大又黄又爽视频免费| 高清午夜精品一区二区三区| 国产精品麻豆人妻色哟哟久久| 国产精品久久久久久精品电影小说 | 国产精品秋霞免费鲁丝片| 国产精品人妻久久久影院| 美女xxoo啪啪120秒动态图| 大片免费播放器 马上看| 久久午夜福利片| 久久久久久伊人网av| 国产成人freesex在线| 久久精品国产a三级三级三级| 插阴视频在线观看视频| 精品久久久久久久久av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 美女视频免费永久观看网站| 午夜福利视频1000在线观看| 国产成人a∨麻豆精品| 婷婷色综合www| 国产成年人精品一区二区| av在线亚洲专区| 丝瓜视频免费看黄片| 少妇裸体淫交视频免费看高清| 亚洲精品日韩av片在线观看| 人人妻人人澡人人爽人人夜夜| 国产精品国产三级专区第一集| 久久精品综合一区二区三区| 午夜福利视频1000在线观看| 最近2019中文字幕mv第一页| 亚洲人与动物交配视频| 18禁在线无遮挡免费观看视频| 国产精品人妻久久久影院| 交换朋友夫妻互换小说| 欧美区成人在线视频| 你懂的网址亚洲精品在线观看| 亚洲精品一二三| 欧美最新免费一区二区三区| av在线观看视频网站免费| 国产av不卡久久| videossex国产| 午夜亚洲福利在线播放| 色吧在线观看| 国产女主播在线喷水免费视频网站| 亚洲在线观看片| 亚洲国产色片| 亚洲国产欧美在线一区| 精品午夜福利在线看| 免费观看无遮挡的男女| 中文字幕亚洲精品专区| 久久久久国产网址| 日韩不卡一区二区三区视频在线| 少妇的逼好多水| 久久久精品欧美日韩精品| 国产精品偷伦视频观看了| 色综合色国产| av在线app专区| 欧美3d第一页| 在线免费观看不下载黄p国产| 2021少妇久久久久久久久久久| 一级毛片我不卡| 成年免费大片在线观看| 国产精品蜜桃在线观看| 成人鲁丝片一二三区免费| 国产精品精品国产色婷婷| 美女主播在线视频| 亚洲精品视频女| av在线蜜桃| 亚洲av日韩在线播放| 26uuu在线亚洲综合色| 黄色怎么调成土黄色| 国产精品久久久久久精品电影小说 | 男女啪啪激烈高潮av片| 国产精品三级大全| 久久久精品欧美日韩精品| 欧美精品国产亚洲| 国产黄频视频在线观看| 高清视频免费观看一区二区| 另类亚洲欧美激情| 精品一区二区三卡| 久久精品熟女亚洲av麻豆精品| 国产熟女欧美一区二区| 午夜亚洲福利在线播放| 极品少妇高潮喷水抽搐| 1000部很黄的大片| 国产av不卡久久| 久久久久久久精品精品| 精品久久久久久久久av| 亚洲真实伦在线观看| 免费av不卡在线播放| av在线亚洲专区| 亚洲三级黄色毛片| 欧美激情久久久久久爽电影| 精华霜和精华液先用哪个| 最近中文字幕高清免费大全6| 下体分泌物呈黄色| 国产精品伦人一区二区| 高清毛片免费看| 男人爽女人下面视频在线观看| 亚洲国产最新在线播放| 男插女下体视频免费在线播放| 亚洲,一卡二卡三卡| 色婷婷久久久亚洲欧美| 亚洲丝袜综合中文字幕| 日韩av在线免费看完整版不卡| 久久久国产一区二区| 亚洲精品久久午夜乱码| 日韩欧美精品免费久久| 日本爱情动作片www.在线观看| 国产伦精品一区二区三区视频9| 熟女电影av网| 成人高潮视频无遮挡免费网站| 黄色一级大片看看| 国产午夜福利久久久久久| av播播在线观看一区| 3wmmmm亚洲av在线观看| 成人国产av品久久久| 色网站视频免费| 熟妇人妻不卡中文字幕| 国产爽快片一区二区三区| 十八禁网站网址无遮挡 | 高清毛片免费看| 欧美变态另类bdsm刘玥| 国产毛片在线视频| 人人妻人人看人人澡| 99热这里只有精品一区| 免费在线观看成人毛片| 建设人人有责人人尽责人人享有的 | 久久精品国产亚洲网站| 免费黄频网站在线观看国产| 精华霜和精华液先用哪个| 欧美少妇被猛烈插入视频| 午夜福利视频精品| 禁无遮挡网站| 亚洲欧美日韩卡通动漫| 中国国产av一级| 超碰av人人做人人爽久久| 亚洲aⅴ乱码一区二区在线播放| 久久精品久久精品一区二区三区| 卡戴珊不雅视频在线播放| 51国产日韩欧美| 欧美日本视频| 男女下面进入的视频免费午夜| 亚洲av.av天堂| xxx大片免费视频| 免费av毛片视频| 国产欧美另类精品又又久久亚洲欧美| 久久精品久久久久久噜噜老黄| 国产真实伦视频高清在线观看| av免费观看日本| 99热网站在线观看| 在线观看三级黄色| 日本爱情动作片www.在线观看| 国产精品99久久99久久久不卡 | 亚洲精品久久午夜乱码| 男人爽女人下面视频在线观看| 成年免费大片在线观看| 狠狠精品人妻久久久久久综合| 99热全是精品| 黄色视频在线播放观看不卡| 水蜜桃什么品种好| 久久国内精品自在自线图片| 最近2019中文字幕mv第一页| 夜夜爽夜夜爽视频| 色吧在线观看| av国产精品久久久久影院| 久热久热在线精品观看| 丝袜喷水一区| 一级爰片在线观看| 久久精品国产鲁丝片午夜精品| 成人美女网站在线观看视频| 欧美激情久久久久久爽电影| 亚洲色图av天堂| 国产乱人视频| 国产高清国产精品国产三级 | 亚洲aⅴ乱码一区二区在线播放| 久久国内精品自在自线图片| 下体分泌物呈黄色| 亚洲va在线va天堂va国产| 大香蕉97超碰在线| 99久国产av精品国产电影| 最近的中文字幕免费完整| 亚洲av不卡在线观看| 国产毛片在线视频| 我的女老师完整版在线观看| 国产精品99久久99久久久不卡 | 亚洲国产最新在线播放| 人妻 亚洲 视频| 久久99蜜桃精品久久| 日本欧美国产在线视频| 国产在线男女| 有码 亚洲区| 国产精品一二三区在线看| 女人久久www免费人成看片| 日本色播在线视频| 少妇猛男粗大的猛烈进出视频 | 亚洲精品国产av蜜桃| 内地一区二区视频在线| 如何舔出高潮| 亚洲国产精品国产精品| 蜜桃久久精品国产亚洲av| 国内精品宾馆在线| 欧美高清性xxxxhd video| 黄色配什么色好看| 亚洲精品久久久久久婷婷小说| 我的女老师完整版在线观看| 国产av国产精品国产| 深爱激情五月婷婷| 免费观看无遮挡的男女| 五月天丁香电影| 国产黄色视频一区二区在线观看| videossex国产| 久久久久网色| 蜜桃亚洲精品一区二区三区| 男人狂女人下面高潮的视频| 精品熟女少妇av免费看| 日本av手机在线免费观看| 99久久人妻综合| 久久热精品热| av在线播放精品| 一级爰片在线观看| 免费观看a级毛片全部| 国模一区二区三区四区视频| 成人二区视频| 国产伦精品一区二区三区四那| 五月开心婷婷网| 人人妻人人看人人澡| 亚洲欧美成人精品一区二区| 日韩,欧美,国产一区二区三区| 亚洲精品aⅴ在线观看| 深夜a级毛片| 我要看日韩黄色一级片| 亚洲精品中文字幕在线视频 | 插逼视频在线观看| 亚洲av电影在线观看一区二区三区 | 亚洲av成人精品一区久久| 亚洲精品日韩在线中文字幕| 欧美性感艳星| 国产精品久久久久久精品电影小说 | 人人妻人人看人人澡| 午夜精品一区二区三区免费看| 日韩在线高清观看一区二区三区| 国产有黄有色有爽视频| 国产精品无大码| 国语对白做爰xxxⅹ性视频网站| 男人添女人高潮全过程视频| 国产男女内射视频| 韩国av在线不卡| 国产成人freesex在线| 成人毛片60女人毛片免费| 精品午夜福利在线看| 免费高清在线观看视频在线观看| 亚洲欧美日韩无卡精品| 亚洲av.av天堂| 久久人人爽av亚洲精品天堂 | 制服丝袜香蕉在线| 啦啦啦中文免费视频观看日本| 亚洲欧美中文字幕日韩二区| 亚洲人成网站高清观看| 女人久久www免费人成看片| 国产精品国产av在线观看| 国产精品99久久久久久久久| 欧美变态另类bdsm刘玥| 色5月婷婷丁香| 麻豆成人av视频| av在线观看视频网站免费| 一区二区av电影网| 最近中文字幕高清免费大全6| 男人爽女人下面视频在线观看| 国产美女午夜福利| 日本欧美国产在线视频| 亚洲人与动物交配视频| 国产在线一区二区三区精| 日韩一区二区视频免费看| 国产永久视频网站| videossex国产| 高清在线视频一区二区三区| 国产精品.久久久| 日产精品乱码卡一卡2卡三| 日本午夜av视频| 亚洲欧洲日产国产| 欧美zozozo另类| 最近的中文字幕免费完整| 国产久久久一区二区三区| 麻豆成人av视频| 亚洲欧美成人综合另类久久久| 大片免费播放器 马上看| 色播亚洲综合网| 一本一本综合久久| 久久人人爽人人片av| 最近的中文字幕免费完整| 国产精品蜜桃在线观看| 国产毛片在线视频| 网址你懂的国产日韩在线| 只有这里有精品99| 国产淫片久久久久久久久| 亚洲av成人精品一区久久| 免费看a级黄色片| 男插女下体视频免费在线播放| 伊人久久精品亚洲午夜| 日本猛色少妇xxxxx猛交久久| 丰满少妇做爰视频| 精品99又大又爽又粗少妇毛片| 中文资源天堂在线| 日本-黄色视频高清免费观看| 在线观看av片永久免费下载| 亚洲av国产av综合av卡| 九九久久精品国产亚洲av麻豆| 夫妻午夜视频| 在线观看国产h片| 超碰97精品在线观看| 三级经典国产精品| 综合色av麻豆| 久久这里有精品视频免费| 免费观看av网站的网址| 91午夜精品亚洲一区二区三区| 亚洲国产欧美人成| 亚洲欧美日韩卡通动漫| 少妇 在线观看| 亚洲人成网站在线播| 麻豆成人av视频| 日韩欧美 国产精品| 美女主播在线视频| 国产黄片视频在线免费观看| 97在线人人人人妻| 精品久久久久久久久亚洲| 男男h啪啪无遮挡| 久久精品夜色国产| 久久久久久久久久人人人人人人| 欧美日韩视频高清一区二区三区二| 最近最新中文字幕大全电影3| 纵有疾风起免费观看全集完整版| 少妇裸体淫交视频免费看高清| 日韩欧美一区视频在线观看 | 日日摸夜夜添夜夜爱| 国产伦理片在线播放av一区| a级毛片免费高清观看在线播放| 国产成人a区在线观看| 亚洲国产av新网站| 69av精品久久久久久| 亚洲三级黄色毛片| 亚洲精品色激情综合| 久热这里只有精品99| 国产成年人精品一区二区| 亚洲国产精品国产精品| 色吧在线观看| 高清在线视频一区二区三区| 日韩中字成人| 国语对白做爰xxxⅹ性视频网站| 麻豆久久精品国产亚洲av| 一本久久精品| h日本视频在线播放| 免费观看性生交大片5| 欧美老熟妇乱子伦牲交| 激情五月婷婷亚洲| 久久精品熟女亚洲av麻豆精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产欧美亚洲国产| 18禁裸乳无遮挡免费网站照片| 交换朋友夫妻互换小说| 女人被狂操c到高潮| 亚洲欧美成人综合另类久久久| 男女那种视频在线观看| 搡老乐熟女国产| 国产一区二区三区综合在线观看 | 真实男女啪啪啪动态图| 亚洲精品国产色婷婷电影| 天天一区二区日本电影三级| 最近手机中文字幕大全| 国产黄色免费在线视频| 嘟嘟电影网在线观看| 少妇人妻久久综合中文| 欧美潮喷喷水| 免费观看性生交大片5| 插逼视频在线观看| 少妇 在线观看| 五月伊人婷婷丁香| 国产精品麻豆人妻色哟哟久久| a级毛片免费高清观看在线播放| 日韩在线高清观看一区二区三区| 亚洲内射少妇av| 男女无遮挡免费网站观看| 两个人的视频大全免费| 插阴视频在线观看视频| 国产伦精品一区二区三区四那| av在线观看视频网站免费| 一二三四中文在线观看免费高清| 晚上一个人看的免费电影| 三级男女做爰猛烈吃奶摸视频| 自拍欧美九色日韩亚洲蝌蚪91 | 久久精品国产鲁丝片午夜精品| 1000部很黄的大片| 午夜亚洲福利在线播放| 中文字幕人妻熟人妻熟丝袜美| 国产永久视频网站| 欧美性感艳星| 成年av动漫网址| 日日摸夜夜添夜夜添av毛片| av国产免费在线观看| 免费观看性生交大片5| 亚洲不卡免费看| 美女脱内裤让男人舔精品视频| 久久久国产一区二区| 女的被弄到高潮叫床怎么办| 老司机影院毛片| 亚洲精品影视一区二区三区av| 黄色一级大片看看| 国产免费视频播放在线视频| 纵有疾风起免费观看全集完整版| 一级毛片久久久久久久久女| 国产色爽女视频免费观看| 欧美三级亚洲精品| 国产男人的电影天堂91| 色婷婷久久久亚洲欧美| 久久久欧美国产精品| 2021天堂中文幕一二区在线观| 91精品伊人久久大香线蕉| 别揉我奶头 嗯啊视频| 天天躁夜夜躁狠狠久久av| 午夜亚洲福利在线播放| 亚洲天堂国产精品一区在线| 亚洲一区二区三区欧美精品 | 亚洲aⅴ乱码一区二区在线播放| 水蜜桃什么品种好| 看黄色毛片网站| 99热全是精品| 一级爰片在线观看| 在线观看av片永久免费下载| 午夜亚洲福利在线播放| 日韩电影二区| 看非洲黑人一级黄片| 丝袜脚勾引网站| 成人高潮视频无遮挡免费网站| 日韩伦理黄色片| 一区二区三区免费毛片| 免费少妇av软件| 国产黄频视频在线观看| 国产精品久久久久久av不卡| 欧美97在线视频| 精品人妻视频免费看| 精品久久久精品久久久| 亚洲精品一二三| 精品久久久久久久久av| 久热这里只有精品99| 国产成人91sexporn| 青春草国产在线视频| 人人妻人人澡人人爽人人夜夜| 国产一区亚洲一区在线观看| 一级黄片播放器| 免费av观看视频| 在线观看av片永久免费下载| 国精品久久久久久国模美| 日韩视频在线欧美| 亚洲精品日韩av片在线观看| 色播亚洲综合网| 午夜视频国产福利| 免费电影在线观看免费观看| 又黄又爽又刺激的免费视频.| 一区二区三区精品91| 在现免费观看毛片| 久久久久国产网址| 最近最新中文字幕免费大全7| 国产精品一二三区在线看| 亚洲电影在线观看av| 亚洲最大成人av| 欧美成人午夜免费资源| 色网站视频免费| 狂野欧美激情性xxxx在线观看| 啦啦啦啦在线视频资源| 精品熟女少妇av免费看| 久久6这里有精品| 亚洲美女搞黄在线观看| 两个人的视频大全免费| 亚洲av电影在线观看一区二区三区 | 99热6这里只有精品| 亚洲精品,欧美精品| 干丝袜人妻中文字幕| 日本免费在线观看一区| 亚洲av电影在线观看一区二区三区 | 亚洲av电影在线观看一区二区三区 | 黄色配什么色好看| 亚洲美女视频黄频| 日本爱情动作片www.在线观看| 欧美一级a爱片免费观看看| 少妇人妻 视频| 中文乱码字字幕精品一区二区三区| 一个人看视频在线观看www免费| 97超视频在线观看视频| 麻豆精品久久久久久蜜桃| 涩涩av久久男人的天堂| 男女啪啪激烈高潮av片| 69av精品久久久久久| 国产精品人妻久久久久久| 一级毛片电影观看| 欧美人与善性xxx| 一二三四中文在线观看免费高清| 精品酒店卫生间| 久久精品国产鲁丝片午夜精品| 日韩国内少妇激情av| 欧美xxxx黑人xx丫x性爽| 九九在线视频观看精品| 久久久色成人| 亚洲一区二区三区欧美精品 | 人人妻人人澡人人爽人人夜夜| 久久久久网色| www.色视频.com| 搡女人真爽免费视频火全软件| 涩涩av久久男人的天堂| 国产男女内射视频| 18+在线观看网站| 久久人人爽av亚洲精品天堂 | 国产一区二区三区综合在线观看 | 两个人的视频大全免费| 最近2019中文字幕mv第一页| 中文欧美无线码|