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

    基于EFAST方法的蘋(píng)果葉片葉綠素含量估算

    2017-05-17 11:29:03楊福芹沙從術(shù)馮海寬韓瑞芳
    河南農(nóng)業(yè)科學(xué) 2017年5期
    關(guān)鍵詞:實(shí)測(cè)值波段葉綠素

    楊福芹,沙從術(shù),馮海寬,韓瑞芳,徐 平

    (1.河南工程學(xué)院 土木工程學(xué)院,河南 鄭州 451191; 2.國(guó)家農(nóng)業(yè)信息化工程技術(shù)研究中心,北京 100097)

    基于EFAST方法的蘋(píng)果葉片葉綠素含量估算

    楊福芹1,2,沙從術(shù)1*,馮海寬2,韓瑞芳1,徐 平1

    (1.河南工程學(xué)院 土木工程學(xué)院,河南 鄭州 451191; 2.國(guó)家農(nóng)業(yè)信息化工程技術(shù)研究中心,北京 100097)

    為了快速、準(zhǔn)確地估算葉綠素含量,使用2012年和2013年在山東省肥城市潮泉鎮(zhèn)獲取的整個(gè)生育期蘋(píng)果葉片葉綠素含量和配套的光譜數(shù)據(jù),利用PROSPECT模型和EFAST方法探討了對(duì)葉綠素含量敏感的波段,然后采用經(jīng)驗(yàn)統(tǒng)計(jì)方法實(shí)現(xiàn)了單波段高光譜對(duì)蘋(píng)果葉片葉綠素含量的監(jiān)測(cè)。結(jié)果表明:以571 nm和697 nm波段光譜參數(shù)為自變量所建立的估測(cè)模型擬合精度較高,其決定系數(shù)(R2)分別為0.71和0.69,均方根誤差(RMSE)分別為1.14、1.17 mg/dm2,相對(duì)誤差(RE)分別為-1.07%和-1.01%。以PROSPECT模型和EFAST方法整合篩選的敏感波段建立的估算模型監(jiān)測(cè)葉綠素含量效果較好,為利用高光譜技術(shù)監(jiān)測(cè)蘋(píng)果長(zhǎng)勢(shì)提供了理論依據(jù)。

    蘋(píng)果葉片; 高光譜; 葉綠素含量; PROSPECT模型; EFAST方法; 隨機(jī)森林

    葉綠素含量是監(jiān)測(cè)果樹(shù)長(zhǎng)勢(shì)的主要生化參數(shù),其變化不僅可以評(píng)價(jià)果樹(shù)的光合作用,而且是評(píng)價(jià)果樹(shù)是否受病害脅迫的指示劑[1-2]。利用高光譜反射特征快速、無(wú)損和準(zhǔn)確監(jiān)測(cè)葉綠素含量已成為監(jiān)測(cè)果樹(shù)長(zhǎng)勢(shì)的重要研究?jī)?nèi)容[3-4]。房賢一等[5]分別利用單變量和多元逐步回歸方法估算了蘋(píng)果盛果期冠層葉綠素含量,結(jié)果表明,用多元逐步回歸方法建立的模型可以更好地監(jiān)測(cè)葉綠素含量。李敏夏等[6]研究發(fā)現(xiàn),篩選出的對(duì)蘋(píng)果葉綠素含量敏感的一階微分波段為694 nm,其決定系數(shù)為0.59。梁爽等[7]分析了一階微分、紅邊位置和蘋(píng)果葉面葉綠素指數(shù)與葉綠素含量之間的相關(guān)性,結(jié)果表明,以葉面葉綠素指數(shù)為變量的估測(cè)模型擬合精度最高。朱西存等[8]利用高光譜紅邊參數(shù)信息,建立了不同物候期蘋(píng)果葉片的SPAD值估算模型,結(jié)果表明,以紅邊位置參數(shù)與葉片SPAD值之間的相關(guān)性最為顯著。冀榮華等[9]研究發(fā)現(xiàn),華北地區(qū)蘋(píng)果葉片葉綠素a含量與反射光譜在515~590 nm和688~715 nm波段內(nèi)高度相關(guān)。鄧小蕾等[10]利用小波變換技術(shù)建立了基于反射光譜的蘋(píng)果葉片葉綠素含量估算模型。劉京等[11]利用支持向量機(jī)建立了蘋(píng)果葉片葉綠素含量高光譜遙感估算模型。韓兆迎等[12]利用連續(xù)統(tǒng)去除法對(duì)原始光譜進(jìn)行處理,采用逐步回歸法和主成分分析法提取主成分,構(gòu)建了基于支持向量機(jī)算法的蘋(píng)果樹(shù)冠SPAD值高光譜估算模型。這些研究利用高光譜手段監(jiān)測(cè)蘋(píng)果葉片葉綠素含量,取得了一定的研究成果,但這些方法篩選的敏感波段比較多,為了提高模型的預(yù)測(cè)能力和穩(wěn)健性,本研究嘗試?yán)肞ROSPECT模型和EFAST方法整合篩選敏感波段,進(jìn)而利用高光譜技術(shù)估算蘋(píng)果葉片葉綠素含量,以期為利用高光譜技術(shù)開(kāi)展果樹(shù)長(zhǎng)勢(shì)及營(yíng)養(yǎng)診斷提供新的思路和方法。

    1 材料和方法

    1.1 研究區(qū)概況和供試材料

    田間試驗(yàn)分別于2012年和2013年在山東省肥城市潮泉鎮(zhèn)下寨村的國(guó)家農(nóng)業(yè)信息化工程技術(shù)研究中心的2個(gè)示范基地的果園(116°50′22″E,36°14′01″N)進(jìn)行。主栽品種為富士和嘎啦,樹(shù)高約3.0 m,果樹(shù)樹(shù)干高約0.5 m,樹(shù)形為紡錘形。在樹(shù)冠東西南北4個(gè)方向的樹(shù)冠外圍中部當(dāng)年延長(zhǎng)枝中部葉(帶葉柄)中分別取8片葉,每株共32片。將樣品迅速放入自封袋中,并將其放在移動(dòng)保溫冰箱中,以減少因呼吸導(dǎo)致的葉片變化。

    1.2 蘋(píng)果葉片光譜測(cè)定

    蘋(píng)果葉片光譜測(cè)定具體方法參考文獻(xiàn)[13-14]。

    1.3 蘋(píng)果葉片葉綠素含量測(cè)定

    用直徑0.6 cm的打孔器從測(cè)量光譜的4片蘋(píng)果葉片切取32小片樣品,樣品質(zhì)量0.2 g左右,然后將樣品放入95%無(wú)水乙醇溶液中,并在黑暗的環(huán)境下靜置24 h,葉片變成白綠色,最后采用紫外分光光度計(jì)測(cè)定蘋(píng)果葉片的葉綠素含量。

    1.4 PROSPECT模型

    著名的PROSPECT模型用于模擬葉片的反射率和透射率[15-16],光譜范圍為400~2 500 nm,間隔1 nm,該模型主要包括6個(gè)輸入?yún)?shù):葉片結(jié)構(gòu)參數(shù)(N)、葉綠素含量(Cab)、等效水厚度(Cw)、類胡蘿卜素含量(Car)、葉黃素含量(Cbrown)和干物質(zhì)含量(Cm),輸出參數(shù)為葉片的反射率和透射率。PROSPECT模型6個(gè)輸入?yún)?shù)的取值范圍如表1所示。

    表1 PROSPECT模型輸入?yún)?shù)的取值范圍

    1.5 EFAST方法

    全局敏感性分析方法主要采用擴(kuò)展傅里葉幅度敏感性檢驗(yàn)方法(extended Fourier amplitude sensitivity,EFAST)[17],它是經(jīng)典傅里葉振幅靈敏度測(cè)試[18-19]方法的擴(kuò)展,由Jacquemoud等[20-21]首先應(yīng)用于遙感領(lǐng)域。EFAST方法考慮了模型參數(shù)的變化范圍對(duì)輸出參量一階敏感性指數(shù)和總體敏感性指數(shù)的影響。一階敏感性指數(shù)指的是模型輸入?yún)?shù)xi的方差占總體方差的百分比,而總體敏感性指數(shù)指的是輸入?yún)?shù)xi以及所有輸入?yún)?shù)xm之間的相互作用的方差占總體方差的百分比。

    1.6 統(tǒng)計(jì)分析

    選取決定系數(shù)(R2)、均方根誤差(RMSE)、相對(duì)誤差(RE)作為評(píng)價(jià)建模與驗(yàn)證精度高低的指標(biāo)。R2表示模擬值與實(shí)測(cè)值的擬合程度,R2越接近于1,表明該擬合曲線精度越高;RMSE和RE反映了模擬值與實(shí)測(cè)值的偏離程度,其值越小,模型精度越高。

    2 結(jié)果與分析

    2.1 利用EFAST方法篩選葉綠素含量敏感波段

    為了研究蘋(píng)果葉片反射光譜與葉綠素含量之間的關(guān)系,采用PROSPECT模型隨機(jī)模擬了2 000組葉片光譜反射率數(shù)據(jù),采用EFAST方法對(duì)PROSPECT模型各個(gè)參數(shù)在400~2 500 nm波段范圍的葉片反射光譜進(jìn)行敏感性分析,利用SimLab軟件實(shí)現(xiàn)定量一階敏感性分析,PROSPECT模型中各參數(shù)變化按照表1中的參數(shù)范圍,并在取值范圍內(nèi)認(rèn)定為均勻分布,模擬結(jié)果如圖1所示。從圖1可以看出,在400~760 nm可見(jiàn)光范圍內(nèi),葉綠素含量對(duì)葉片反射率的影響非常大,一階敏感指數(shù)大于0.4的波段范圍分別為417~497 nm和532~726 nm,其中,葉綠素含量的反射光譜一階敏感指數(shù)較大值分別出現(xiàn)在藍(lán)波段(440 nm)、紅波段(571 nm)和近紅外波段(697 nm)。在660 nm附近葉綠素對(duì)葉片反射率的影響在可見(jiàn)光區(qū)域最低。綜上,波段440、571、697 nm處的反射光譜參數(shù)對(duì)葉綠素含量最為敏感。

    圖1 葉綠素含量反射光譜的敏感指數(shù)

    2.2 估算模型的構(gòu)建

    利用2013年采集的180個(gè)蘋(píng)果葉片葉綠素含量及對(duì)應(yīng)的光譜數(shù)據(jù),采用經(jīng)驗(yàn)回歸方法建立了3個(gè)單波段反射光譜敏感波段的葉片葉綠素含量估算模型。結(jié)果表明,處于440 nm波段的光譜參數(shù)與

    葉綠素含量完全不相關(guān),這是由于受環(huán)境因素的影響以及吸收的強(qiáng)度不一樣。從表2可以看出,處于571 nm和697 nm波段處的光譜參數(shù)與葉綠素含量都達(dá)到了極顯著相關(guān)水平(P<0.01),其R2分別為0.71和0.69,RMSE分別為1.14、1.17 mg/dm2,RE分別為-1.07%和-1.01%,預(yù)測(cè)值與實(shí)測(cè)值具有較高的一致性,建模精度較為理想,表明采用571 nm和697 nm波段處的光譜參數(shù)進(jìn)行葉綠素含量估算具有較好的可靠性。

    2.3 估算模型的驗(yàn)證

    為了進(jìn)一步驗(yàn)證模型的可靠性,利用2012年采集的297個(gè)山東地區(qū)蘋(píng)果相關(guān)數(shù)據(jù)進(jìn)行可靠性檢驗(yàn),如表2和圖2所示,571 nm和697 nm單波段算法的R2分別為0.11和0.19,RMSE分別為1.33、1.55 mg/dm2,RE分別為5.92%和13.20%。571 nm波段處的光譜參數(shù)估算的葉綠素含量與實(shí)測(cè)值之間存在很好的關(guān)系,預(yù)測(cè)值和實(shí)測(cè)值均勻分布在1∶1線的兩側(cè),但決定系數(shù)為0.11,僅達(dá)到顯著相關(guān)水平(P<0.05)。697 nm波段處的預(yù)測(cè)值和實(shí)測(cè)值之間的關(guān)系較好,其決定系數(shù)為0.19,達(dá)到極顯著相關(guān)水平(P<0.01),預(yù)測(cè)值和實(shí)測(cè)值多數(shù)處于1∶1線之上,預(yù)測(cè)值存在高估現(xiàn)象,可能由于氣候原因,2012年采集蘋(píng)果數(shù)據(jù)生育期與2013年采集蘋(píng)果數(shù)據(jù)生育期不完全一致,導(dǎo)致實(shí)測(cè)值普遍小于預(yù)測(cè)值。

    表2 基于EFAST方法的葉綠素含量建模及驗(yàn)證

    注:x—反射光譜(nm);y—葉綠素含量(mg/dm2)。

    a.571 nm; b.697 nm

    3 結(jié)論

    本研究使用2012年和2013年地面實(shí)測(cè)的全生育期蘋(píng)果葉片葉綠素含量和相應(yīng)的葉片光譜數(shù)據(jù),首先采用PROSPECT模型隨機(jī)模擬了2 000組數(shù)據(jù),然后采用EFAST方法篩選出對(duì)葉片葉綠素含量的敏感波段分別為440、571、697 nm,根據(jù)選擇的敏感波段利用經(jīng)驗(yàn)統(tǒng)計(jì)方法建立了3個(gè)單波段線性回歸方程,估算葉片葉綠素含量。結(jié)果表明,571、697 nm所建立的模型可以提高蘋(píng)果葉片葉綠素含量的估算精度,增加模型的實(shí)用性,具有一定的應(yīng)用前景。該估算模型的R2分別為0.71和0.69,RMSE分別為1.14、1.17 mg/dm2,RE分別為-1.07%和-1.01%。驗(yàn)證結(jié)果同樣表明,葉片葉綠素含量的預(yù)測(cè)值與實(shí)測(cè)值具有較好的一致性,其RMSE分別為1.33、1.55 mg/dm2,驗(yàn)證模型的RE分別為5.92%和13.20%,驗(yàn)證模型都具有較高的精度與可靠性。本試驗(yàn)結(jié)果表明,用PROSPECT模型和EFAST方法整合篩選的敏感波段可以較好地估算蘋(píng)果葉片葉綠素含量,這為利用高光譜特征參數(shù)監(jiān)測(cè)蘋(píng)果長(zhǎng)勢(shì)情況提供了新的理論依據(jù)。

    [1] 徐新剛,趙春江,王紀(jì)華,等.新型光譜曲線特征參數(shù)與水稻葉綠素含量間的關(guān)系研究[J].光譜學(xué)與光譜分析,2011,31(1):188-191.

    [2] 牛魯燕,孫家波,劉延忠,等.基于成像高光譜的小麥葉片葉綠素含量估測(cè)模型研究[J].河南農(nóng)業(yè)科學(xué),2016,45(4):150-154.

    [3] 王秀珍,王人潮,黃敬峰.微分光譜遙感及其在水稻農(nóng)學(xué)參數(shù)測(cè)定上的應(yīng)用研究[J].農(nóng)業(yè)工程學(xué)報(bào),2002,18(1):9-13.

    [4] 劉偉東,項(xiàng)月琴.高光譜數(shù)據(jù)與水稻指數(shù)及葉綠素密度的相關(guān)分析[J].遙感學(xué)報(bào),2000,4(4):279-283.

    [5] 房賢一,朱西存,王凌,等.基于高光譜的蘋(píng)果盛果期冠層葉綠素含量監(jiān)測(cè)研究[J].中國(guó)農(nóng)業(yè)科學(xué),2013,46(16):3504-3513.

    [6] 李敏夏,張林森,李丙智,等.蘋(píng)果葉片高光譜特性與葉綠素含量和 SPAD 值的關(guān)系[J].西北林學(xué)院學(xué)報(bào),2010,25(2):35-39.

    [7] 梁爽,趙庚星,朱西存.蘋(píng)果樹(shù)葉片葉綠素含量高光譜估測(cè)模型研究[J].光譜學(xué)與光譜分析,2012,32(5):1367-1370.

    [8] 朱西存,趙庚星,姜遠(yuǎn)茂,等.基于高光譜紅邊參數(shù)的不同物候期蘋(píng)果葉片的SPAD值估測(cè)[J].紅外,2011,32(12):31-38.

    [9] 冀榮華,鄭立華,鄧小蕾,等.基于反射光譜的蘋(píng)果葉片葉綠素和含水率預(yù)測(cè)模型[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2014,45(8):269-275.

    [10] 鄧小蕾,李民贊,鄭立華,等.基于反射光譜預(yù)處理的蘋(píng)果葉片葉綠素含量預(yù)測(cè)[J].農(nóng)業(yè)工程學(xué)報(bào),2014,30(14):140-147.

    [11] 劉京,常慶瑞,劉淼,等.基于SVR算法的蘋(píng)果葉片葉綠素含量高光譜反演[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2016,47(8):260-265.

    [12] 韓兆迎,朱西存,王凌,等.基于連續(xù)統(tǒng)去除法的蘋(píng)果樹(shù)冠 SPAD 高光譜估測(cè)[J].激光與光電子學(xué)進(jìn)展,2016,53(2):214-223.

    [13] 馮海寬,楊福芹,李振海,等.最優(yōu)權(quán)重組合模型和高光譜估算蘋(píng)果葉片全磷含量[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(7):173-180.

    [14] 馮海寬,李振海,金秀良,等.基于EFAST和PLS的蘋(píng)果葉片等效水厚度高光譜估算[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(12):165-171.

    [15] Jacquemoud S,Ustin S L,Verdebout J,etal.Estimating leaf biochemistry using the PROSPECT leaf optical properties model[J].Remote Sensing of Environment,1996,56(3):194-202.

    [16] Ustin S L,Roberts D A,Pinzon J,etal.Estimating canopy water content of chaparral shrubs using optical methods[J].Remote Sensing of Environment,1998,65(3):280-291.

    [17] Cukier R I,Levine H B,Shuler K E.Nonlinear sensitivity analysis of multiparameter model systems[J].Journal of Computational Physics,1978,26(1):1-42.

    [18] Ceccato P,Flasse S,Tarantola S,etal.Detecting vegetation leaf water content using reflectance in the optical domain[J].Remote Sensing of Environment,2001,77(1):22-33.

    [19] Ceccato P,Gobron N,Flasse S,etal.Designing a spectral index to estimate vegetation water content from remote sensing data:Part 1:Theoretical approach[J].Remote Sensing of Environment,2002,82(2):188-197.

    [20] Jacquemoud S,Baret F.PROSPECT:A model of leaf optical properties spectra[J].Remote Sensing of Environment,1990,34(2):75-91.

    [21] Jacquemoud S,Baret F,Andrieu B,etal.Extraction of vegetation biophysical parameters by inversion of the PROSPECT+SAIL models on sugar beet canopy reflectance data.Application to TM and AVIRIS sensors[J].Remote Sensing of Environment,1995,52(3):163-172.

    Estimation of Chlorophyll Content in Apple Leaves Using Hyperspectral Data Based on EFAST Method

    YANG Fuqin1,2,SHA Congshu1*,FENG Haikuan2,HAN Ruifang1,XU Ping1

    (1.College of Civil Engineering,Henan Institute of Engineering,Zhengzhou 451191,China; 2.National Engineering Research Center for Information Technology in Agriculture,Beijing 100097,China)

    In order to estimate chlorophyll content quickly and exactly,chlorophyll content parameters and the concurrent spectral reflectance of apple leaves were acquired in Chaoquan town,Feicheng city,Shandong province,during 2012 and 2013 apple growth seasons.Sensitive wavebands to chlorophyll content were selected using PROSPECT model and EFAST(extended Flourier amplitude sensitivity test) method,and then estimation models of chlorophyll content were built using empirical statistical methods.Results showed that the fitting accuracy of the estimation models using wavebands of 571 nm and 697 nm was higher.The determination coefficients(R2) were 0.71 and 0.69,the root mean square errors(RMSE) were 1.14,1.17 mg/dm2,and the relative errors(RE) were -1.07% and -1.01%.The model established by using PROSPECT model and EFAST method can predict the apple leaf chlorophyll content better,providing a theoretical basis for monitoring apple growth conditions by using hyperspectral technology.

    apple leaves; hyperspectra; chlorophyll content; PROSPECT model; extended Flourier amplitude sensitivity test; random forest

    2016-11-25

    國(guó)家自然科學(xué)基金項(xiàng)目(41601346);北京市自然科學(xué)基金項(xiàng)目(4141001);國(guó)家高技術(shù)研究發(fā)展計(jì)劃863課題(2011AA100703)

    楊福芹(1979-),女,河南安陽(yáng)人,講師,博士,主要從事農(nóng)業(yè)定量遙感研究。E-mail:yangfuqin0202@163.com

    *通訊作者:沙從術(shù)(1964-),男,安徽霍邱人,教授,碩士,主要從事測(cè)繪工程方面的教學(xué)與研究。 E-mail:shacongshu123@126.com

    S661.1

    A

    1004-3268(2017)05-0157-04

    猜你喜歡
    實(shí)測(cè)值波段葉綠素
    春日暖陽(yáng)
    ±800kV直流輸電工程合成電場(chǎng)夏季實(shí)測(cè)值與預(yù)測(cè)值比對(duì)分析
    提取葉綠素
    常用高溫軸承鋼的高溫硬度實(shí)測(cè)值與計(jì)算值的對(duì)比分析
    哈爾濱軸承(2020年1期)2020-11-03 09:16:22
    市售純牛奶和巴氏殺菌乳營(yíng)養(yǎng)成分分析
    桃樹(shù)葉綠素含量與SPAD值呈極顯著正相關(guān)
    一種基于實(shí)測(cè)值理論計(jì)算的導(dǎo)航臺(tái)電磁干擾分析方法
    電子制作(2018年23期)2018-12-26 01:01:22
    葉綠素家族概述
    M87的多波段輻射過(guò)程及其能譜擬合
    日常維護(hù)對(duì)L 波段雷達(dá)的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    十八禁人妻一区二区| 一本精品99久久精品77| 一区福利在线观看| 少妇人妻精品综合一区二区 | 欧美成人性av电影在线观看| 无遮挡黄片免费观看| 可以在线观看毛片的网站| 狂野欧美白嫩少妇大欣赏| 丰满人妻一区二区三区视频av| 欧美区成人在线视频| 日本在线视频免费播放| 毛片一级片免费看久久久久 | 18美女黄网站色大片免费观看| 久久精品夜夜夜夜夜久久蜜豆| 在线天堂最新版资源| 亚洲av成人精品一区久久| 亚洲男人的天堂狠狠| 免费看美女性在线毛片视频| 男女做爰动态图高潮gif福利片| 国产精品98久久久久久宅男小说| 精品国内亚洲2022精品成人| xxxwww97欧美| 男人舔女人下体高潮全视频| 午夜精品在线福利| 日本 av在线| 韩国av一区二区三区四区| 国产中年淑女户外野战色| 久久亚洲精品不卡| 日韩精品中文字幕看吧| 1024手机看黄色片| 搡老岳熟女国产| 国产又黄又爽又无遮挡在线| 在线观看美女被高潮喷水网站 | netflix在线观看网站| 亚洲精品影视一区二区三区av| 我的老师免费观看完整版| 日韩欧美国产在线观看| 精品免费久久久久久久清纯| 亚洲美女黄片视频| 国产探花在线观看一区二区| 久9热在线精品视频| 亚洲美女视频黄频| 搡老熟女国产l中国老女人| 在线免费观看不下载黄p国产 | 亚洲欧美激情综合另类| 麻豆成人午夜福利视频| 亚洲专区中文字幕在线| 在线播放无遮挡| 久久99热6这里只有精品| 国产精品爽爽va在线观看网站| 一进一出抽搐gif免费好疼| 夜夜躁狠狠躁天天躁| 中国美女看黄片| 亚洲成人中文字幕在线播放| 亚洲精品在线美女| 精品一区二区三区av网在线观看| 桃红色精品国产亚洲av| 噜噜噜噜噜久久久久久91| 国产一区二区在线观看日韩| 欧美激情国产日韩精品一区| 在线播放国产精品三级| 欧美成人a在线观看| 欧美一区二区国产精品久久精品| 国产精品综合久久久久久久免费| 99国产精品一区二区三区| 久久人妻av系列| 男人狂女人下面高潮的视频| 少妇的逼好多水| 国内少妇人妻偷人精品xxx网站| 国产黄片美女视频| 赤兔流量卡办理| 精品99又大又爽又粗少妇毛片 | 哪里可以看免费的av片| 婷婷色综合大香蕉| 国产精品98久久久久久宅男小说| 99热精品在线国产| 男人狂女人下面高潮的视频| 久9热在线精品视频| 欧美成人一区二区免费高清观看| or卡值多少钱| 十八禁网站免费在线| 国产不卡一卡二| 不卡一级毛片| 国产免费男女视频| 日日夜夜操网爽| 蜜桃久久精品国产亚洲av| 51国产日韩欧美| 性色avwww在线观看| av在线天堂中文字幕| 久久久久久久久大av| 国产熟女xx| www.熟女人妻精品国产| 看免费av毛片| 超碰av人人做人人爽久久| eeuss影院久久| 亚洲成av人片在线播放无| 可以在线观看毛片的网站| 真实男女啪啪啪动态图| 国产精品亚洲av一区麻豆| 99视频精品全部免费 在线| 久久久久久久久久成人| 91在线观看av| 亚洲天堂国产精品一区在线| 最新在线观看一区二区三区| 欧美xxxx黑人xx丫x性爽| 日韩欧美免费精品| 精品久久久久久久久av| 色在线成人网| 亚洲成人免费电影在线观看| 国产精品不卡视频一区二区 | 人妻丰满熟妇av一区二区三区| 内地一区二区视频在线| 国产69精品久久久久777片| 男人的好看免费观看在线视频| 99国产精品一区二区蜜桃av| 亚洲av第一区精品v没综合| 在线观看美女被高潮喷水网站 | 亚洲国产色片| 日韩高清综合在线| 免费人成视频x8x8入口观看| 日韩欧美免费精品| 99在线人妻在线中文字幕| 日韩成人在线观看一区二区三区| 欧美一区二区精品小视频在线| 我的老师免费观看完整版| 熟女电影av网| 国内精品一区二区在线观看| 欧美激情国产日韩精品一区| 国产伦精品一区二区三区四那| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 97热精品久久久久久| 美女免费视频网站| 18+在线观看网站| 久久性视频一级片| 精品久久久久久久久亚洲 | 久久精品综合一区二区三区| 亚洲av五月六月丁香网| 国产成人aa在线观看| 露出奶头的视频| 国内少妇人妻偷人精品xxx网站| 亚洲av不卡在线观看| 中国国产av一级| 久久精品人妻少妇| 男女那种视频在线观看| 亚洲,一卡二卡三卡| 2022亚洲国产成人精品| 成人一区二区视频在线观看| 国产av国产精品国产| 只有这里有精品99| 午夜福利高清视频| 永久网站在线| 亚洲成色77777| 三级经典国产精品| 大陆偷拍与自拍| .国产精品久久| 欧美最新免费一区二区三区| 99精国产麻豆久久婷婷| 97超碰精品成人国产| 欧美日本视频| 波多野结衣巨乳人妻| 日韩视频在线欧美| 人人妻人人看人人澡| 青春草视频在线免费观看| 中文在线观看免费www的网站| 久久人人爽人人片av| 一本久久精品| 精品少妇久久久久久888优播| 有码 亚洲区| 欧美少妇被猛烈插入视频| 亚洲av福利一区| 国产成人免费观看mmmm| av免费观看日本| 亚洲欧美精品专区久久| 久久99精品国语久久久| av在线蜜桃| 寂寞人妻少妇视频99o| 国产高清不卡午夜福利| 黄片无遮挡物在线观看| 精品人妻偷拍中文字幕| 女人久久www免费人成看片| 久久精品熟女亚洲av麻豆精品| 日韩电影二区| 夜夜爽夜夜爽视频| 涩涩av久久男人的天堂| 国产精品麻豆人妻色哟哟久久| 搡老乐熟女国产| 日韩av免费高清视频| 亚洲精品国产色婷婷电影| 夫妻午夜视频| 亚洲av电影在线观看一区二区三区 | 亚洲色图av天堂| 不卡视频在线观看欧美| 久久鲁丝午夜福利片| 岛国毛片在线播放| 99re6热这里在线精品视频| 搞女人的毛片| 亚洲av中文字字幕乱码综合| 乱系列少妇在线播放| 国产一区二区亚洲精品在线观看| 国产精品一区二区在线观看99| 亚洲美女搞黄在线观看| 久久国内精品自在自线图片| 国产精品一及| 成年女人看的毛片在线观看| 久久久精品免费免费高清| 国产精品99久久99久久久不卡 | 在线观看人妻少妇| 日韩 亚洲 欧美在线| 毛片一级片免费看久久久久| 大香蕉97超碰在线| 麻豆精品久久久久久蜜桃| 又爽又黄无遮挡网站| 麻豆精品久久久久久蜜桃| 国产精品久久久久久久久免| 国产成人精品一,二区| 国产精品人妻久久久久久| 看免费成人av毛片| 色网站视频免费| 人妻夜夜爽99麻豆av| 国产黄频视频在线观看| 久久久久久伊人网av| 欧美日韩国产mv在线观看视频 | 五月天丁香电影| 女的被弄到高潮叫床怎么办| 成人亚洲精品av一区二区| 中文字幕人妻熟人妻熟丝袜美| 欧美潮喷喷水| 国产探花极品一区二区| 一个人看的www免费观看视频| 男的添女的下面高潮视频| 免费看a级黄色片| 美女主播在线视频| 一本色道久久久久久精品综合| 综合色av麻豆| 免费观看无遮挡的男女| 亚洲国产精品成人久久小说| 嘟嘟电影网在线观看| 久久久久久久亚洲中文字幕| 日韩一区二区视频免费看| 精品酒店卫生间| 一级毛片黄色毛片免费观看视频| 天天躁夜夜躁狠狠久久av| 26uuu在线亚洲综合色| 少妇人妻 视频| 国产成人一区二区在线| 99精国产麻豆久久婷婷| 免费大片黄手机在线观看| 国产精品国产三级专区第一集| 麻豆成人av视频| 亚洲色图av天堂| 交换朋友夫妻互换小说| 婷婷色综合大香蕉| 亚洲内射少妇av| 在线观看免费高清a一片| 国产中年淑女户外野战色| 老女人水多毛片| 精品少妇黑人巨大在线播放| 成人毛片a级毛片在线播放| 日本wwww免费看| 国产视频首页在线观看| 婷婷色综合www| 亚洲成人一二三区av| 男女无遮挡免费网站观看| 人人妻人人澡人人爽人人夜夜| 免费看a级黄色片| 亚洲av福利一区| 性插视频无遮挡在线免费观看| 丝袜脚勾引网站| 午夜精品国产一区二区电影 | 久久热精品热| 爱豆传媒免费全集在线观看| 69人妻影院| 黄色视频在线播放观看不卡| 新久久久久国产一级毛片| 亚洲,欧美,日韩| 国产乱人视频| 国产欧美另类精品又又久久亚洲欧美| 国产高清有码在线观看视频| 免费黄色在线免费观看| 国产成年人精品一区二区| 国产高清有码在线观看视频| 免费看av在线观看网站| 亚洲丝袜综合中文字幕| 成人亚洲欧美一区二区av| 国产精品无大码| 一级a做视频免费观看| 插阴视频在线观看视频| 在线亚洲精品国产二区图片欧美 | 国产日韩欧美在线精品| 日本熟妇午夜| 国产极品天堂在线| 大码成人一级视频| 午夜福利视频精品| 国产一级毛片在线| 国产又色又爽无遮挡免| 国内精品美女久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 国产精品女同一区二区软件| 水蜜桃什么品种好| 免费少妇av软件| 成人亚洲欧美一区二区av| 精品久久久久久久久亚洲| 18禁在线播放成人免费| 亚洲欧美一区二区三区国产| 国产一区有黄有色的免费视频| 精品少妇黑人巨大在线播放| 色综合色国产| 亚洲精品日韩在线中文字幕| 国国产精品蜜臀av免费| 男人和女人高潮做爰伦理| 精品人妻一区二区三区麻豆| 一个人看的www免费观看视频| 国产 一区精品| 精品一区在线观看国产| 国产精品国产三级国产专区5o| 一级毛片aaaaaa免费看小| 成人亚洲欧美一区二区av| 草草在线视频免费看| 国产综合精华液| 国产精品99久久99久久久不卡 | 久久久久精品久久久久真实原创| 国产精品女同一区二区软件| 少妇熟女欧美另类| a级一级毛片免费在线观看| 国产免费一级a男人的天堂| 精品午夜福利在线看| 亚洲综合精品二区| 国产综合懂色| 精品国产一区二区三区久久久樱花 | 亚洲国产精品成人久久小说| 国产日韩欧美亚洲二区| av天堂中文字幕网| 亚洲国产av新网站| 天天躁日日操中文字幕| 狠狠精品人妻久久久久久综合| 免费看光身美女| 亚洲欧美日韩无卡精品| 国产精品精品国产色婷婷| 日韩一区二区三区影片| 小蜜桃在线观看免费完整版高清| 精品国产露脸久久av麻豆| 少妇被粗大猛烈的视频| 午夜免费鲁丝| 又爽又黄无遮挡网站| 欧美人与善性xxx| 在线观看一区二区三区| 亚洲av一区综合| 亚洲精品456在线播放app| 在线精品无人区一区二区三 | 人人妻人人澡人人爽人人夜夜| 精品国产三级普通话版| 嫩草影院入口| 尤物成人国产欧美一区二区三区| 波野结衣二区三区在线| 七月丁香在线播放| 毛片女人毛片| 黄色欧美视频在线观看| 免费观看av网站的网址| 国产亚洲91精品色在线| 亚洲最大成人手机在线| 国产一区有黄有色的免费视频| 美女脱内裤让男人舔精品视频| 一区二区av电影网| 日韩欧美精品免费久久| 亚洲天堂国产精品一区在线| 国产男女超爽视频在线观看| 我的女老师完整版在线观看| 亚洲四区av| 欧美变态另类bdsm刘玥| 天天躁夜夜躁狠狠久久av| 在线观看国产h片| 久久精品国产亚洲网站| 99re6热这里在线精品视频| eeuss影院久久| 一区二区三区免费毛片| 一级爰片在线观看| 日韩一本色道免费dvd| 欧美一区二区亚洲| 国产v大片淫在线免费观看| 少妇高潮的动态图| 国产精品久久久久久精品电影| 一区二区三区乱码不卡18| 国产淫语在线视频| 91久久精品电影网| 久久久久久国产a免费观看| 亚洲精品成人av观看孕妇| 亚洲,欧美,日韩| 精华霜和精华液先用哪个| 亚洲aⅴ乱码一区二区在线播放| 秋霞伦理黄片| 亚洲熟女精品中文字幕| 国产精品女同一区二区软件| 成年免费大片在线观看| 久久鲁丝午夜福利片| www.色视频.com| 亚洲av免费高清在线观看| 美女被艹到高潮喷水动态| 亚洲精品亚洲一区二区| 欧美日韩视频高清一区二区三区二| 女人被狂操c到高潮| 久久久久久久久久久免费av| 午夜福利在线观看免费完整高清在| 在线观看三级黄色| 亚洲人与动物交配视频| 亚洲,欧美,日韩| 精品国产三级普通话版| 国产亚洲av嫩草精品影院| 欧美一级a爱片免费观看看| 午夜福利网站1000一区二区三区| 日韩强制内射视频| 人妻系列 视频| 伊人久久精品亚洲午夜| 亚洲欧洲国产日韩| 18禁在线无遮挡免费观看视频| 国产成人a区在线观看| 精品人妻偷拍中文字幕| 欧美日韩视频高清一区二区三区二| 免费看日本二区| 国产老妇女一区| 久久精品综合一区二区三区| 亚洲精品国产av蜜桃| 国产免费福利视频在线观看| 听说在线观看完整版免费高清| 别揉我奶头 嗯啊视频| 国产女主播在线喷水免费视频网站| 久久ye,这里只有精品| 亚洲图色成人| 久久精品人妻少妇| 日日摸夜夜添夜夜爱| 国产色婷婷99| 国产精品无大码| 亚洲欧美清纯卡通| 久久亚洲国产成人精品v| av国产免费在线观看| 禁无遮挡网站| 国产又色又爽无遮挡免| av播播在线观看一区| 激情五月婷婷亚洲| 亚洲精品中文字幕在线视频 | 国产精品一区二区三区四区免费观看| 国产精品一及| 一级黄片播放器| 最近中文字幕2019免费版| 狂野欧美激情性xxxx在线观看| 成年av动漫网址| 男的添女的下面高潮视频| 色播亚洲综合网| 亚洲av中文字字幕乱码综合| 99热国产这里只有精品6| 午夜福利在线观看免费完整高清在| 国产精品成人在线| 午夜福利视频精品| 日韩一区二区视频免费看| 久久精品国产自在天天线| 国产片特级美女逼逼视频| 精品少妇黑人巨大在线播放| 免费不卡的大黄色大毛片视频在线观看| 色婷婷久久久亚洲欧美| 97超碰精品成人国产| 亚洲av.av天堂| 国产淫片久久久久久久久| 少妇高潮的动态图| 久久精品综合一区二区三区| 成年av动漫网址| 男女啪啪激烈高潮av片| 午夜老司机福利剧场| 亚洲无线观看免费| 爱豆传媒免费全集在线观看| 免费看av在线观看网站| 插逼视频在线观看| 校园人妻丝袜中文字幕| 精品一区二区三卡| 日韩中字成人| 狂野欧美激情性bbbbbb| 黄色日韩在线| 舔av片在线| av黄色大香蕉| 欧美日韩视频精品一区| 日韩亚洲欧美综合| 黄色一级大片看看| 日本猛色少妇xxxxx猛交久久| 久久久久久久久久久丰满| 国产黄色免费在线视频| 免费观看的影片在线观看| 亚洲最大成人中文| 国国产精品蜜臀av免费| 中文字幕人妻熟人妻熟丝袜美| 国产午夜精品久久久久久一区二区三区| 18禁裸乳无遮挡免费网站照片| 熟女av电影| 亚州av有码| 中文字幕免费在线视频6| 91午夜精品亚洲一区二区三区| 一区二区三区免费毛片| 久久久久久久久久成人| 国产黄a三级三级三级人| 亚洲美女搞黄在线观看| 亚洲激情五月婷婷啪啪| 欧美成人一区二区免费高清观看| 中文资源天堂在线| av又黄又爽大尺度在线免费看| 女的被弄到高潮叫床怎么办| 日韩人妻高清精品专区| 美女cb高潮喷水在线观看| 日韩,欧美,国产一区二区三区| 欧美国产精品一级二级三级 | 一区二区av电影网| 热re99久久精品国产66热6| 身体一侧抽搐| 欧美成人a在线观看| 欧美+日韩+精品| 国产av不卡久久| 国产乱人偷精品视频| 欧美日韩视频精品一区| 国产美女午夜福利| 午夜福利在线观看免费完整高清在| 亚洲aⅴ乱码一区二区在线播放| 熟女av电影| 久久久久久久午夜电影| 男女那种视频在线观看| 草草在线视频免费看| 国产精品久久久久久精品古装| 中文天堂在线官网| 午夜福利视频精品| 人妻 亚洲 视频| 又大又黄又爽视频免费| 亚洲av免费高清在线观看| 精品酒店卫生间| 日韩成人av中文字幕在线观看| 国产精品人妻久久久久久| 亚洲国产精品专区欧美| 大香蕉97超碰在线| 日日摸夜夜添夜夜添av毛片| 精品国产乱码久久久久久小说| 好男人视频免费观看在线| 中文字幕久久专区| 欧美日韩在线观看h| 人妻少妇偷人精品九色| 亚洲av二区三区四区| 在线观看美女被高潮喷水网站| 国产又色又爽无遮挡免| 人人妻人人看人人澡| 一级毛片电影观看| 尾随美女入室| 亚洲精品色激情综合| 国产一区有黄有色的免费视频| 国产国拍精品亚洲av在线观看| 波多野结衣巨乳人妻| 精品久久国产蜜桃| 亚洲精华国产精华液的使用体验| 在线观看美女被高潮喷水网站| 国产欧美亚洲国产| 好男人视频免费观看在线| 在线天堂最新版资源| 在线 av 中文字幕| 日韩av不卡免费在线播放| 久久精品久久久久久久性| av国产免费在线观看| 久久精品国产鲁丝片午夜精品| 国产成年人精品一区二区| 亚洲三级黄色毛片| 日韩强制内射视频| 亚洲在久久综合| 亚洲精品影视一区二区三区av| 肉色欧美久久久久久久蜜桃 | 我的老师免费观看完整版| 狂野欧美激情性bbbbbb| 亚洲精品,欧美精品| 亚洲国产成人一精品久久久| 亚洲伊人久久精品综合| 国产免费又黄又爽又色| 男女国产视频网站| 日韩国内少妇激情av| 国模一区二区三区四区视频| 日本黄大片高清| 久久精品夜色国产| 深夜a级毛片| 又爽又黄a免费视频| 亚洲自偷自拍三级| 午夜免费男女啪啪视频观看| 18禁裸乳无遮挡免费网站照片| 欧美日韩在线观看h| 成人一区二区视频在线观看| 久久午夜福利片| 国产午夜精品久久久久久一区二区三区| 欧美一级a爱片免费观看看| 国产免费视频播放在线视频| 婷婷色麻豆天堂久久| 国产大屁股一区二区在线视频| 日本av手机在线免费观看| 22中文网久久字幕| 亚洲天堂国产精品一区在线| 嫩草影院新地址| 亚洲激情五月婷婷啪啪| 丝袜喷水一区| 亚洲国产精品成人综合色| 国产极品天堂在线| 天堂网av新在线| 69av精品久久久久久| freevideosex欧美| 亚洲在线观看片| 午夜福利视频精品| 日本与韩国留学比较| 国产日韩欧美在线精品| 91久久精品国产一区二区成人| 蜜臀久久99精品久久宅男| 成人综合一区亚洲| 好男人视频免费观看在线| 欧美日韩在线观看h| 久热这里只有精品99| 久久久久九九精品影院| 国产一区二区在线观看日韩| kizo精华| 亚洲综合精品二区| 日韩欧美精品v在线|