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

    基于無人機多光譜與地面高光譜遙感的土壤主要養(yǎng)分含量估測

    2022-02-15 14:47:16楊棟膗李亞強刀劍王建雄
    江蘇農(nóng)業(yè)科學 2022年2期
    關鍵詞:高光譜土壤養(yǎng)分支持向量機

    楊棟膗 李亞強 刀劍 王建雄

    摘要:耕地土壤作為農(nóng)業(yè)資源的重要組成部分,是我國農(nóng)業(yè)生產(chǎn)的基本資源與保障條件。山原紅壤作為云南典型的土壤類型,在云南的農(nóng)業(yè)生產(chǎn)與經(jīng)濟發(fā)展中發(fā)揮著重要作用。通過多光譜與高光譜遙感結(jié)合技術,對山原紅壤主要養(yǎng)分含量的高光譜特性進行了研究分析,研究了土壤主要養(yǎng)分含量的多高光譜特性、篩選出了各養(yǎng)分含量的特征波段。而后,通過多光譜波段變換、植被指數(shù)合成、構(gòu)建高光譜土壤指數(shù),與山原紅壤主要養(yǎng)分含量進行了相關分析。通過分析得出,多光譜擬合模型效果整體不佳,僅鉀元素與合成光譜存在相關性,擬合模型R2大于0.5,有良好的擬合意義。高光譜指數(shù)相關性及擬合模型效果較多光譜都有所提升。其中,支持向量機擬合模型效果最佳,R2均在0.75以上,適合用于土壤養(yǎng)分含量估測。

    關鍵詞:土壤養(yǎng)分;多光譜;植被指數(shù);高光譜;支持向量機

    中圖分類號:P237;S127? 文獻標志碼: A

    文章編號:1002-1302(2022)02-0178-08

    收稿日期:2021-04-16

    基金項目:云南省教育廳科學研究基金(編號:2019Y0080、2020Y0177、2020Y186)。

    作者簡介:楊棟淏(1997—),男,云南曲靖人,碩士研究生,主要從事資源與環(huán)境遙感專業(yè)研究。E-mail:843765886@qq.com。

    通信作者:王建雄,博士,教授,碩士生導師,主要研究方向為資源與環(huán)境遙感。E-mail:jianxiongw@126.com。

    云南省地處高原山區(qū),耕地資源有限,且山原紅壤作為云南典型的土壤類型,研究山原紅壤的土壤肥力對云南的農(nóng)業(yè)生產(chǎn)具有重要的作用。近年來,遙感技術快速發(fā)展,對土壤調(diào)查、土地利用變化及土壤肥力監(jiān)測起到了獨特作用[1]。目前,針對土壤多光譜方面的研究大多針對土壤含鹽量、含水量與重金屬的反演模型建立,少有針對土壤主要養(yǎng)分含量的分析以及估測研究[2-4],并且國內(nèi)的研究多以土壤有機質(zhì)的反演估算為主[5-9],而且取得了較高的反演精度。例如,南鋒等以土壤有機質(zhì)含量為研究對象,研究分析后發(fā)現(xiàn),基于顯著波段構(gòu)建的偏最小二乘回歸(PLSR)模型預測精度高于基于全波段的PLSR模型[10]。在高光譜方面,欒福明等使用全波段地物光譜儀針對速效鉀養(yǎng)分含量進行了模型估測研究,發(fā)現(xiàn)小波分析處理能顯著提高模型精度[11]。蔣璐璐等分別采用偏最小二乘支持向量機和偏最小二乘-人工神經(jīng)網(wǎng)絡2種方法進行土壤磷、鉀元素含量預測建模,結(jié)果表明,中紅外波段對磷和鉀的預測效果更好[12]。盡管已有不少學者針對土壤養(yǎng)分含量的定量估測進行了研究,但不同土壤類型反射光譜存在差異,使得各土壤類型的養(yǎng)分估測模型不具備普遍適宜性。而目前針對山原紅壤光譜反射特性的研究不多,且針對其主要養(yǎng)分含量的估測更為少見。因此,本研究結(jié)合上述學者的研究方法,分析山原紅壤速效氮、磷、鉀含量的光譜反射特性,以期為今后云南的農(nóng)業(yè)生產(chǎn)提供幫助。

    1 數(shù)據(jù)獲取與預處理

    1.1 數(shù)據(jù)獲取

    1.1.1 土壤采集 綜合地理位置、無人機影像采集可行性等多方面因素,本研究以昆明市云南農(nóng)業(yè)大學后山試驗田為研究對象。試驗田為山原紅壤地塊,位于25°8′N、102°45′E,地勢走向為北高南低,海拔為1 959.17 m和1 923.05 m。試驗地塊為 2.33 hm2 的長方形區(qū)域。經(jīng)土壤區(qū)分后,選取符合山原紅壤特征的區(qū)域地塊進行密集采樣,選取20個采樣點采集土壤樣本,采樣時間為2019年11月16日。土樣采集分布示意圖見圖1。土壤樣本經(jīng)風干、研磨并通過 2 mm 孔篩,采用四分法分成2份,分別供化學分析和光譜測試用。

    1.1.2 多光譜數(shù)據(jù)采集 使用掛載派諾特鏡頭的

    大疆M100無人機飛行器,進行灰板校正后在晴朗少云天氣對實驗地塊進行多光譜影像采集。飛行高度為60 m,飛行航向重疊率為70%,旁向重疊率為75%,飛行速度在12 m/s以下,拍照間隔時間為 2 s。手持無人機多光譜相機對白板上放置好的20個土樣進行多光譜影像的室內(nèi)采集。

    1.1.3 高光譜數(shù)據(jù)采集 利用FieldSpec4 NG全波段地物光譜儀對所有土樣進行采集,光譜儀波長范圍350~2 500、350~1 000 nm波段樣本間隔 1.4 nm,1 001~2 500 nm波段樣本間隔1.1 nm,數(shù)據(jù)輸出通道數(shù)2 151個。測試前進行白板矯正,每個土樣樣本均采集10條光譜反射曲線,對各樣本的光譜反射值取平均值得到各土樣實際光譜反射值。

    1.2 數(shù)據(jù)預處理

    1.2.1 影像拼接 對無人機采集的照片進行剔除和篩選后利用Pix4D mapper軟件對照片進行拼接獲取實驗地的各波段影像、數(shù)字正射影像(DOM)、數(shù)字表面模型(DSM)、試驗區(qū)四波段多光譜影像。

    1.2.2 土壤養(yǎng)分測定 通過化學和物理方法結(jié)合對土樣的有機質(zhì)、速效氮磷鉀、全氮共5種土壤主要養(yǎng)分的含量檢測。

    1.2.3 光譜反射曲線優(yōu)化處理 由于地物光譜儀在低波段與高波段采用了不同傳感器,在光譜反射曲線中會存在斷點,因而用Viewspec Pro軟件對原始光譜反射率進行斷點修正,消除斷點陡變。同時使用Origin軟件中的Savitizky-Golay平滑處理,消除光譜曲線中的噪聲影響,處理后的光譜曲線仍然保持原土壤光譜曲線特征。

    2 結(jié)果與分析

    2.1 山原紅壤多光譜反演研究

    2.1.1 單波段土壤光譜反演研究

    2.1.1.1 土壤主要養(yǎng)分含量敏感波段分析 利用多光譜四波段土壤反射率與測定的土壤養(yǎng)分含量進行正態(tài)分布檢驗后,分別對土壤主要養(yǎng)分與各波段光譜進行Pearson相關分析,結(jié)果見表1。

    由表1可知,除速效磷和有機質(zhì)外,其余養(yǎng)分含量與光譜之間存在顯著或極顯著相關性。因此,對速效氮、磷和全氮3個指標進行回歸分析,對主要養(yǎng)分含量進行多光譜模型估測。

    2.1.1.2 回歸模型分析 分別使用單變量和雙變量回歸模型對土壤養(yǎng)分進行分析。經(jīng)試驗檢驗,建立的單變量回歸模型效果差,改進使用雙相關變量波段組合的方式對山原紅壤養(yǎng)分反射率進行研究。通過研究分析發(fā)現(xiàn),僅有速效鉀的養(yǎng)分含量與近紅外和綠波存在相關性,因此選取兩波段為自變量,速效鉀含量為因變量建立回歸模型。分別構(gòu)建土壤養(yǎng)分含量的線性回歸和非線性回歸方程(表2)。

    由表2可知,通過雙變量波段組合能夠有效提高土壤養(yǎng)分速效鉀擬合模型的R2,使得土壤反射率能夠更好地反映土壤養(yǎng)分。但整體的R2仍然較低,無法選定速效鉀的土壤反射率為回歸模型。

    2.1.2 合成指數(shù)土壤光譜反演研究 研究引入9種植被指數(shù):歸一化植被指數(shù)(NDVI)、比值植被指數(shù)(RVI)、差值植被指數(shù)(DVI)、紅邊歸一化植被指數(shù)(NDVI705)、紅綠比值指數(shù)(RG)、綠通道植被指數(shù)(GNDVI)、葉面葉綠素指數(shù)(LCI)、歸一化差異紅邊指數(shù)(NDRE)、土壤調(diào)節(jié)植被指數(shù)(OSAVI)對山原紅壤的養(yǎng)分進行相關性分析。經(jīng)檢驗,僅有速效鉀元素含量與較多的土壤指數(shù)(NDVI、RVI、DVI、RG、OSAVI)存在顯著相關性。

    2.1.2.1 單變量回歸分析 選取速效鉀元素為因變量,與其具有較好相關性的5種植被指數(shù)作物自變量,利用SPSS軟件進行單變量模型估測。通過分析得出,NDVI、DVI和OSAVI這3種土壤指數(shù)擬合模型的R2均大于0.5,但整體模型的擬合誤差都比較大,須要對養(yǎng)分含量進行更優(yōu)的模型擬合。

    2.1.2.2 多元回歸分析 選取單變量中R2最高的3種植被指數(shù),將模型兩兩組合構(gòu)建模型表達式,最終擬合3類曲線回歸模型如表3所示。

    由表3可知,多元變量曲線回歸的擬合程度更高,且選入變量越多,擬合程度越高,選擇三元變量三次擬合模型能較好地反演鉀元素養(yǎng)分含量。

    2.1.3 二維土壤光譜反演研究

    進行一維分析后,再通過多光譜二維光譜合成進行土壤養(yǎng)分含量分析,分別對波段進行加減乘除和歸一化進行二維處理。經(jīng)上述分析可知,養(yǎng)分含量反演中,鉀元素與波段和各合成指數(shù)的相關性較好。因此,利用合成波段分別與鉀元素進行相關性分析,結(jié)果如圖2所示。通過熱力圖(圖2)可以看出,K元素在不同模型中,相關性較高的波段組合如表4所示。將以上相關性較好的組合提出,構(gòu)建新的模型,在構(gòu)建模型中乘法擬合模型與歸一化擬合的模型R2低于0.5,因此不使用這2種方法擬合模型,通過對最優(yōu)合成波段進行再次組合和篩選,最終擬合模型詳見表5。由表5分析可得,進行組合后的三元模型擬合程度更好,除法擬合模型中的組合2和組合10構(gòu)成的反演模型能更好地反演鉀元素養(yǎng)分含量。

    2.2 山原紅壤高光譜反演研究

    2.2.1 單相關分析的養(yǎng)分模型估測

    在上面的多光譜遙感中,利用合成植被指數(shù)和波段組合模型能較好地反演鉀元素含量。在構(gòu)建高光譜指數(shù)時,利用對原始光譜反射率(R)進行不同形式的變換,最終選取包含原始反射率、反射率倒數(shù)、反射率對數(shù)、反射率平方根、反射率倒數(shù)的對數(shù)、反射率對數(shù)的倒數(shù)以及以上指數(shù)的一階導數(shù)共12種變換光譜作為此次研究的輸入光譜。

    通過對土壤速效氮含量與光譜指數(shù)的相關分析得到,研究中構(gòu)建的反射率、平方根、倒數(shù)、對數(shù)等幾種光譜指數(shù)對土壤主要養(yǎng)分含量相關性均不明顯,而經(jīng)過一階導數(shù)處理后的光譜指數(shù)相關性顯著提升(表6)。通過相同方法對其他主要養(yǎng)分含量進行相關分析后發(fā)現(xiàn)均存在此類現(xiàn)象,因此僅挑選具有特征波段研究特性的土壤合成指數(shù)進行研究。

    提取土壤各主要養(yǎng)分含量的特征波段后,選取與各主要養(yǎng)分含量相關系數(shù)最優(yōu)的波段及合成土壤指數(shù)波段組合進行擬合分析。分別進行一元回歸分析和多元回歸分析,一元回歸分析中除速效磷元素的曲線擬合模型R2大于0.7外,其他養(yǎng)分含量擬合模型的R2均不高,一元回歸模型無法完全表述變量信息。通過改進后使用多元回歸分析,除全氮元素的擬合模型R2有所提升外,其余土壤養(yǎng)分含量擬合模型沒有改善,也無法較好地表達變量貢獻率。因此,結(jié)合曲線擬合和多元回歸結(jié)果,再次構(gòu)建非線性回歸模型對土壤養(yǎng)分含量進行估測。

    綜合一元回歸模型和多元回歸模型,選取一元回歸分析中擬合值最優(yōu)的兩模型構(gòu)建土壤二元回歸曲線回歸方程,計算出各養(yǎng)分含量擬合曲線方程如表7所示。

    經(jīng)二元曲線擬合后的土壤養(yǎng)分模型整體較多元線性逐步回歸分析的擬合模型R2均有所提升,但RMSE檢驗值也相對上升。相對于一元回歸分析的擬合極值,除速效氮養(yǎng)分含量擬合值有所上升外,其他養(yǎng)分含量擬合值并未上升。

    利用相關分析法進行分析,一定程度上可以提取土壤主要養(yǎng)分含量的的特征波段范圍,為回歸分析提供依據(jù)。但選取相關系數(shù)極值點的敏感波段進行擬合的方法難以體現(xiàn)波段間的相關影響,僅靠相關系數(shù)極值選定變量構(gòu)建模型進行養(yǎng)分回歸估測這一方法存在一定片面性,取得的建模效果并不理想。

    經(jīng)研究發(fā)現(xiàn),機器學習算法在土壤有機質(zhì)等主要養(yǎng)分含量的估測研究中,能夠較好地預測精度[15]。結(jié)合前人的研究[16-17],本研究選取了偏最小二乘(PLSR)回歸與支持向量機(SVM)2種方法對山原紅壤主要養(yǎng)分含量進行估測研究。

    通過以上相關性分析發(fā)現(xiàn),除一階微分變換能較好地提升光譜反射曲線與養(yǎng)分含量的相關性外,其余光譜變換均無很好效果。因此,在PLSR回歸和SVM回歸分析中,只使用原始反射率和平方根一階導光譜指數(shù)進行分析。

    2.2.2 PLSR回歸養(yǎng)分含量反演研究

    通過The Unscrambler X多變量數(shù)據(jù)分析軟件對目標數(shù)據(jù)進行處理。以20個山原紅壤土壤樣本的原始光譜率及平方根一階導光譜指數(shù)為自變量,速效氮磷鉀及其余養(yǎng)分含量為因變量。各主要養(yǎng)分含量回歸分析結(jié)果如圖3至圖7所示。由表8和圖3至圖7可知,利用PLSR回歸分析得到的結(jié)果中,速效氮磷鉀元素一階微分變換指數(shù)的R2較原始反射率有所提高,但有機質(zhì)和全氮一階微分變換指數(shù)的R2降低,并且PLSR擬合的模型整體效果提升很小,不適用于山原紅壤養(yǎng)分含量估測。

    2.2.3 SVM回歸養(yǎng)分含量反演研究

    與PLSR回歸相同,選擇20個山原紅壤土壤原始反射率和平方根一階微分作為自變量,各土壤養(yǎng)分含量作為因變量,分別對5種養(yǎng)分含量進行建模分析,各主要養(yǎng)分含量回歸分析結(jié)果如圖8至圖12所示。通過對山原紅壤速效氮磷鉀及其他主要養(yǎng)分含量的SVM回歸分析可以看出,各主要養(yǎng)分含量在與光譜反射率平方根一階微分的回歸模型中,均獲得了較好的擬合度。R2依舊呈現(xiàn)一階微分變換指數(shù)高于原始光譜反射率這一現(xiàn)象,且一階微分變換光譜指數(shù)的R2較原始光譜反射率有明顯提升,RMSE檢驗值也相對降低??傊?,使用SVM回歸分析方法進行擬合的模型,整體效果均有所提升,比較其他方法更適合用于進行養(yǎng)分含量的估測。

    3 結(jié)論與討論

    3.1 結(jié)論

    分別使用多光譜和高光譜結(jié)合對山原紅壤主要養(yǎng)分含量進行估測。在多光譜養(yǎng)分含量估測中,山原紅壤在紅外波段較為敏感,其土壤反射率值高于其他波段。利用合成植被指數(shù)對比,除比值植被指數(shù)與差值植被指數(shù)波段較為明顯外,其他植被指數(shù)波動范圍較小,通過原始反射率、合成植被指數(shù)與土壤各養(yǎng)分含量的相關分析表明,合成波段的相關性高于單波段。應用多種回歸分析方法對與土壤養(yǎng)分含量相關性較好的多光譜植被指數(shù)進行分析,整體R2系數(shù)較低,無法較好地擬合土壤指數(shù)與土壤養(yǎng)分的關系,造成這一結(jié)果的原因初步估測為分析土壤養(yǎng)分樣本量較少或多光譜探測波段分辨率不夠高導致。在多光譜養(yǎng)分估測中,利用多種方式擬合了鉀元素的合成光譜回歸模型,其R2大于0.5,有良好的擬合意義。

    由于多光譜進行養(yǎng)分估測效果不理想,再利用高光譜對其進行分析。對特征波段與山原紅壤進行多種回歸模型擬合,單相關分析擬合的模型效果不好,無法較好地反映土壤養(yǎng)分含量。而后使用機器學習法進行分析,分別使用最小二乘法回歸模型和支持向量機回歸模型對不同高光譜反射曲線進行擬合分析,擬合結(jié)果更好,使用支持向量機擬合模型的R2較高,RMSE均方根誤差相對較好,比較適用于山原紅壤養(yǎng)分含量的估測,實現(xiàn)了土壤主要養(yǎng)分含量的高光譜估測。

    3.2 討論

    基于高光譜多變量的土壤養(yǎng)分含量回歸分析由于變量眾多因而分析較為復雜,基于機器學習的回歸算法能夠較好地解決多變量擬合回歸的過擬合現(xiàn)象,但需要大量的樣本數(shù)據(jù)。本研究相對其他研究而言樣本變量相對較少,因而未采用獨立樣本檢驗的方法對模型精度進行驗證。同時,由于樣本變量較少,在進行BP(back propagation)神經(jīng)網(wǎng)絡與隨機森林回歸分析時,未能獲得較好的擬合效果。因此在今后的多變量估測研究中,為保證擬合模型的擬合精度與獨立樣本檢驗的實現(xiàn),應盡可能提升樣本數(shù)量,提升樣本檢測的可能性。

    參考文獻:

    [1]邱 壑. 耕地土壤有機質(zhì)與速效氮磷鉀含量高光譜遙感反演研究[D]. 福州:福建農(nóng)林大學,2017.

    [2]王騰軍,趙明海,楊 耘,等. 多光譜影像的陜西大西溝礦區(qū)土壤重金屬含量反演[J]. 光譜學與光譜分析,2019,39(12):3880-3887.

    [3]張智韜,魏廣飛,姚志華,等. 基于無人機多光譜遙感的土壤含鹽量反演模型研究[J]. 農(nóng)業(yè)機械學報,2019,50(12):151-160.

    [4]張智韜,譚丞軒,許崇豪,等. 基于無人機多光譜遙感的玉米根域土壤含水率研究[J]. 農(nóng)業(yè)機械學報,2019,50(7):246-257.

    [5]趙 瑞,崔希民,劉 超. GF-5高光譜遙感影像的土壤有機質(zhì)含量反演估算研究[J]. 中國環(huán)境科學,2020,40(8):3539-3545.

    [6]劉煥軍,張美薇,楊昊軒,等. 多光譜遙感結(jié)合隨機森林算法反演耕作土壤有機質(zhì)含量[J]. 農(nóng)業(yè)工程學報,2020,36(10):134-140.

    [7]姜雪芹,葉 勤,林 怡,等. 基于諧波分析和高光譜遙感的土壤含水量反演研究[J]. 光學學報,2017,37(10):292-302.

    [8]馬 馳. 基于Landsat 8遙感圖像的長春中北部地區(qū)土壤有機質(zhì)含量反演[J]. 江蘇農(nóng)業(yè)科學,2016,44(12):415-419.

    [9]楊建鋒,馬軍成,王令超. 基于地貌類型的土壤有機質(zhì)多光譜遙感反演[J]. 地理空間信息,2015,13(2):48-50.

    [10]南 鋒,朱洪芬,畢如田.黃土高原煤礦區(qū)復墾農(nóng)田土壤有機質(zhì)含量的高光譜預測[J]. 中國農(nóng)業(yè)科學,2016,49(11):2126-2135.

    [11]欒福明,熊黑鋼,王 芳,等. 基于小波分析的土壤速效K含量高光譜反演[J]. 干旱區(qū)地理,2015,38(2):320-326.

    [12]蔣璐璐,張 瑜,王艷艷,等. 基于光譜技術的土壤養(yǎng)分快速測試方法研究[J]. 浙江大學學報(農(nóng)業(yè)與生命科學版),2010,36(4):445-450.

    [13]石萌萌. 全譜段成像儀的地表場景模擬[D]. 西安:西安科技大學,2017.

    [14]秦倩如. 基于高光譜的陜西省土壤游離鐵與碳酸鈣含量預測研究[D]. 楊凌:西北農(nóng)林科技大學,2019.

    [15]葉元元. 多金屬礦區(qū)土壤重金屬的高光譜定量估算研究[D]. 徐州:中國礦業(yè)大學,2014.

    [16]紀文君,李 曦,李成學,等. 基于全譜數(shù)據(jù)挖掘技術的土壤有機質(zhì)高光譜預測建模研究[J]. 光譜學與光譜分析,2012,32(9):2393-2398.

    [17]沈蘭芝,高懋芳,閆敬文,等. 基于SVR和PLSR的土壤有機質(zhì)高光譜估測模型研究[J]. 中國農(nóng)業(yè)信息,2019,31(1):58-71.

    猜你喜歡
    高光譜土壤養(yǎng)分支持向量機
    基于實測光譜的紅花葉綠素含量高光譜估算模型研究
    不同施肥模式對油茶植株營養(yǎng)生長和土壤養(yǎng)分的影響
    廣靈縣平川區(qū)土壤養(yǎng)分變化及施肥建議
    動態(tài)場景中的視覺目標識別方法分析
    論提高裝備故障預測準確度的方法途徑
    價值工程(2016年32期)2016-12-20 20:36:43
    基于熵技術的公共事業(yè)費最優(yōu)組合預測
    價值工程(2016年29期)2016-11-14 00:13:35
    一種基于引導濾波和MNF的高光譜遙感圖像分類方法
    軟件導刊(2016年9期)2016-11-07 22:19:22
    基于支持向量機的金融數(shù)據(jù)分析研究
    稻蟹共作模式下稻蟹產(chǎn)出與土壤理化性質(zhì)的研究
    基于高光譜響應的農(nóng)田大棚EVA棚膜拉伸強度估計
    精品久久久久久久久av| 午夜亚洲福利在线播放| 91久久精品电影网| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 村上凉子中文字幕在线| 欧美最黄视频在线播放免费| 日韩欧美精品免费久久 | 国产69精品久久久久777片| 亚洲av第一区精品v没综合| 一级a爱片免费观看的视频| 一本综合久久免费| xxxwww97欧美| 色精品久久人妻99蜜桃| 黄片小视频在线播放| 久久热精品热| 亚洲一区二区三区不卡视频| 色av中文字幕| 精品熟女少妇八av免费久了| 久久久久久国产a免费观看| 久久久久亚洲av毛片大全| 成人一区二区视频在线观看| 一级黄片播放器| 久久午夜福利片| 亚洲av免费在线观看| 少妇丰满av| 国产老妇女一区| 少妇熟女aⅴ在线视频| 91九色精品人成在线观看| 国产高清三级在线| 老司机午夜十八禁免费视频| 久久久久国内视频| 国内毛片毛片毛片毛片毛片| 国产91精品成人一区二区三区| 国产真实伦视频高清在线观看 | 成人永久免费在线观看视频| 亚洲精品色激情综合| 两个人的视频大全免费| 亚洲最大成人中文| 亚洲美女搞黄在线观看 | 成年版毛片免费区| 国产精品一区二区免费欧美| 久久久久性生活片| 国产精华一区二区三区| 亚洲av日韩精品久久久久久密| 九九久久精品国产亚洲av麻豆| 国产精品久久久久久亚洲av鲁大| 99视频精品全部免费 在线| 欧美最新免费一区二区三区 | 国产三级在线视频| 日本黄色视频三级网站网址| 欧美区成人在线视频| 欧美+日韩+精品| 中文在线观看免费www的网站| 两人在一起打扑克的视频| 国产成人影院久久av| 一个人免费在线观看电影| 嫩草影院精品99| 宅男免费午夜| 亚洲av美国av| 我的老师免费观看完整版| 久久久久久久久久成人| 国产伦精品一区二区三区四那| 国产色婷婷99| 真实男女啪啪啪动态图| 一本精品99久久精品77| 99久久成人亚洲精品观看| 精品99又大又爽又粗少妇毛片 | 别揉我奶头~嗯~啊~动态视频| 禁无遮挡网站| 久久久色成人| 午夜精品久久久久久毛片777| 久久午夜亚洲精品久久| 亚洲国产精品成人综合色| a级毛片免费高清观看在线播放| 色哟哟·www| 丰满人妻一区二区三区视频av| 丰满的人妻完整版| 欧美在线黄色| 亚洲成人精品中文字幕电影| 91午夜精品亚洲一区二区三区 | 网址你懂的国产日韩在线| 免费av观看视频| 国产精品不卡视频一区二区 | 高清日韩中文字幕在线| 国产午夜精品久久久久久一区二区三区 | 国产蜜桃级精品一区二区三区| 此物有八面人人有两片| 中文亚洲av片在线观看爽| .国产精品久久| 国产av麻豆久久久久久久| 国产精品久久电影中文字幕| 国产单亲对白刺激| 亚洲av不卡在线观看| 如何舔出高潮| 亚洲一区高清亚洲精品| 18禁黄网站禁片午夜丰满| 美女cb高潮喷水在线观看| 噜噜噜噜噜久久久久久91| 欧美精品国产亚洲| 亚洲av.av天堂| 日本 av在线| 此物有八面人人有两片| 搡老岳熟女国产| 国产一区二区在线观看日韩| 91av网一区二区| 男人舔奶头视频| 久久精品人妻少妇| 麻豆久久精品国产亚洲av| 精品一区二区三区人妻视频| 欧美三级亚洲精品| 麻豆久久精品国产亚洲av| 999久久久精品免费观看国产| 1000部很黄的大片| 不卡一级毛片| 亚洲成人久久性| 精品不卡国产一区二区三区| 成人欧美大片| 国产精品久久电影中文字幕| 99热这里只有是精品50| av专区在线播放| 99热这里只有是精品50| 免费一级毛片在线播放高清视频| 精品无人区乱码1区二区| 久久99热6这里只有精品| 一级黄片播放器| 性欧美人与动物交配| 国产美女午夜福利| 国产淫片久久久久久久久 | 真实男女啪啪啪动态图| 亚洲三级黄色毛片| 午夜精品一区二区三区免费看| 精品久久久久久久人妻蜜臀av| 看片在线看免费视频| 色综合婷婷激情| 99久久精品一区二区三区| 中文在线观看免费www的网站| 亚洲欧美日韩高清在线视频| 真人做人爱边吃奶动态| 免费一级毛片在线播放高清视频| 小蜜桃在线观看免费完整版高清| 黄色配什么色好看| 成人高潮视频无遮挡免费网站| 亚洲欧美日韩东京热| 757午夜福利合集在线观看| 黄色一级大片看看| 在线十欧美十亚洲十日本专区| 亚洲va日本ⅴa欧美va伊人久久| 国内精品一区二区在线观看| 99国产精品一区二区三区| 伦理电影大哥的女人| 在线观看免费视频日本深夜| 国产极品精品免费视频能看的| 全区人妻精品视频| 757午夜福利合集在线观看| 黄色一级大片看看| 日韩中文字幕欧美一区二区| 国产免费一级a男人的天堂| 99国产极品粉嫩在线观看| 亚洲国产精品999在线| 一区福利在线观看| 亚洲国产精品sss在线观看| 精品一区二区免费观看| 欧美黄色片欧美黄色片| 有码 亚洲区| 欧美成人性av电影在线观看| 99国产精品一区二区蜜桃av| 午夜老司机福利剧场| 18禁黄网站禁片午夜丰满| 亚洲一区高清亚洲精品| 国产在线精品亚洲第一网站| 欧美成人a在线观看| 国产精品久久视频播放| 久久九九热精品免费| 非洲黑人性xxxx精品又粗又长| 精品日产1卡2卡| 国产成人a区在线观看| av专区在线播放| 精品久久久久久久人妻蜜臀av| 国产成+人综合+亚洲专区| 国产精品99久久久久久久久| 好看av亚洲va欧美ⅴa在| 欧美日韩瑟瑟在线播放| 国内精品一区二区在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 日韩国内少妇激情av| 日韩精品中文字幕看吧| 欧美不卡视频在线免费观看| 欧美三级亚洲精品| 国产精品综合久久久久久久免费| 欧美+亚洲+日韩+国产| 亚洲真实伦在线观看| 99国产极品粉嫩在线观看| 亚洲成人精品中文字幕电影| 精品国产三级普通话版| 男女做爰动态图高潮gif福利片| 久久精品国产亚洲av香蕉五月| 91麻豆精品激情在线观看国产| 亚洲精品色激情综合| 综合色av麻豆| 黄色视频,在线免费观看| 亚洲av熟女| 亚洲精品影视一区二区三区av| 色播亚洲综合网| 99riav亚洲国产免费| 日本一二三区视频观看| 国产欧美日韩一区二区精品| bbb黄色大片| 亚洲欧美激情综合另类| 老女人水多毛片| 国产av麻豆久久久久久久| 日韩欧美在线乱码| 国产熟女xx| 一个人看视频在线观看www免费| 国产免费一级a男人的天堂| 日本黄大片高清| av在线蜜桃| 好男人电影高清在线观看| 一进一出抽搐gif免费好疼| 国产男靠女视频免费网站| 国产不卡一卡二| 91麻豆精品激情在线观看国产| 在线免费观看的www视频| 88av欧美| 国产精品日韩av在线免费观看| 99久久久亚洲精品蜜臀av| 欧美日韩黄片免| 久久久久亚洲av毛片大全| 亚洲国产精品sss在线观看| 色吧在线观看| av视频在线观看入口| 久久久久九九精品影院| 麻豆成人午夜福利视频| 久久午夜亚洲精品久久| 国产美女午夜福利| 色综合欧美亚洲国产小说| 亚洲欧美日韩高清在线视频| 亚洲av中文字字幕乱码综合| 久久精品国产亚洲av涩爱 | 99久久精品一区二区三区| 在线观看舔阴道视频| 成人高潮视频无遮挡免费网站| 在线观看66精品国产| 日韩人妻高清精品专区| 美女大奶头视频| 波多野结衣高清作品| 一个人观看的视频www高清免费观看| 午夜福利视频1000在线观看| 亚洲人与动物交配视频| 91麻豆av在线| 午夜精品一区二区三区免费看| 欧美乱色亚洲激情| 日本免费一区二区三区高清不卡| 精华霜和精华液先用哪个| 免费高清视频大片| 欧美成狂野欧美在线观看| 国产aⅴ精品一区二区三区波| 亚洲在线观看片| 国产午夜精品论理片| 麻豆成人av在线观看| 国产一区二区在线观看日韩| 91麻豆精品激情在线观看国产| 国产精品爽爽va在线观看网站| 99riav亚洲国产免费| 97超级碰碰碰精品色视频在线观看| 中文字幕免费在线视频6| 搞女人的毛片| 熟女电影av网| 欧美又色又爽又黄视频| 一a级毛片在线观看| 国产成人福利小说| 亚洲欧美日韩卡通动漫| 女人被狂操c到高潮| 国产高清三级在线| 国产精品乱码一区二三区的特点| 在线观看免费视频日本深夜| 精品人妻1区二区| 色综合亚洲欧美另类图片| 最近中文字幕高清免费大全6 | 久久人人精品亚洲av| 久久99热这里只有精品18| 在线免费观看的www视频| 偷拍熟女少妇极品色| 麻豆成人午夜福利视频| 日韩免费av在线播放| 97人妻精品一区二区三区麻豆| 亚洲av免费在线观看| 一个人观看的视频www高清免费观看| 能在线免费观看的黄片| 人妻久久中文字幕网| 人人妻人人看人人澡| av福利片在线观看| 神马国产精品三级电影在线观看| 男女床上黄色一级片免费看| 岛国在线免费视频观看| 中文字幕精品亚洲无线码一区| 免费黄网站久久成人精品 | 91麻豆精品激情在线观看国产| 嫩草影院新地址| 国产成人aa在线观看| 亚洲精品乱码久久久v下载方式| 真人做人爱边吃奶动态| .国产精品久久| 国产精品影院久久| АⅤ资源中文在线天堂| 女人十人毛片免费观看3o分钟| 成人av在线播放网站| 国内久久婷婷六月综合欲色啪| 国产精品亚洲美女久久久| 999久久久精品免费观看国产| 女生性感内裤真人,穿戴方法视频| 亚洲av第一区精品v没综合| 国产成人av教育| 国产色爽女视频免费观看| 亚洲欧美日韩东京热| 婷婷精品国产亚洲av在线| 俄罗斯特黄特色一大片| 国产淫片久久久久久久久 | 国产成+人综合+亚洲专区| 三级国产精品欧美在线观看| 在线观看一区二区三区| 一个人免费在线观看的高清视频| 国产亚洲欧美98| 天堂av国产一区二区熟女人妻| 国内精品美女久久久久久| 精华霜和精华液先用哪个| 俄罗斯特黄特色一大片| bbb黄色大片| 真实男女啪啪啪动态图| 1000部很黄的大片| 亚洲精品色激情综合| 老鸭窝网址在线观看| 国产亚洲欧美98| 国产精品人妻久久久久久| 精品久久久久久久末码| 久久99热6这里只有精品| 少妇的逼水好多| 日本与韩国留学比较| or卡值多少钱| 国产亚洲欧美98| 搡老熟女国产l中国老女人| 精品人妻视频免费看| 亚洲性夜色夜夜综合| 18禁黄网站禁片午夜丰满| 国产欧美日韩精品一区二区| 国产视频一区二区在线看| 国产私拍福利视频在线观看| 高清毛片免费观看视频网站| 天堂网av新在线| 亚洲专区中文字幕在线| 亚洲av日韩精品久久久久久密| 日本一二三区视频观看| 真实男女啪啪啪动态图| 成人毛片a级毛片在线播放| 成人特级av手机在线观看| 免费在线观看亚洲国产| 亚洲中文日韩欧美视频| 身体一侧抽搐| 国产一区二区激情短视频| bbb黄色大片| 欧美日本视频| 久久国产乱子伦精品免费另类| 最近视频中文字幕2019在线8| 日韩亚洲欧美综合| 精品久久久久久久久久免费视频| 亚洲 欧美 日韩 在线 免费| 美女黄网站色视频| 在现免费观看毛片| 男女那种视频在线观看| 99视频精品全部免费 在线| 热99在线观看视频| 免费大片18禁| xxxwww97欧美| 亚洲中文日韩欧美视频| 中文字幕人成人乱码亚洲影| 久久6这里有精品| 久9热在线精品视频| 久久久久性生活片| 亚洲av电影不卡..在线观看| 美女被艹到高潮喷水动态| 亚洲成人久久性| 波多野结衣高清作品| 亚洲成av人片免费观看| 精品国内亚洲2022精品成人| av在线蜜桃| 有码 亚洲区| 国产成人av教育| 看黄色毛片网站| 午夜福利在线在线| 有码 亚洲区| 免费在线观看成人毛片| 欧美日韩瑟瑟在线播放| 国内揄拍国产精品人妻在线| 91久久精品国产一区二区成人| 一卡2卡三卡四卡精品乱码亚洲| 国内精品一区二区在线观看| 少妇丰满av| 亚洲欧美清纯卡通| 国产精品一区二区免费欧美| 又紧又爽又黄一区二区| 国产高清激情床上av| 中亚洲国语对白在线视频| 国产久久久一区二区三区| 亚洲最大成人中文| 黄色女人牲交| 一本一本综合久久| 啦啦啦韩国在线观看视频| 神马国产精品三级电影在线观看| 中文字幕熟女人妻在线| 亚洲欧美激情综合另类| 看片在线看免费视频| 神马国产精品三级电影在线观看| 亚洲av电影在线进入| 精品午夜福利视频在线观看一区| 免费看光身美女| 黄色女人牲交| 亚洲最大成人手机在线| 国产探花在线观看一区二区| 久久精品人妻少妇| 观看美女的网站| 久久国产精品人妻蜜桃| 成人av一区二区三区在线看| 亚洲av日韩精品久久久久久密| 国产麻豆成人av免费视频| 免费看光身美女| 国产蜜桃级精品一区二区三区| 桃红色精品国产亚洲av| 国产激情偷乱视频一区二区| 久久精品国产自在天天线| 日本 欧美在线| 亚洲欧美日韩高清在线视频| 波多野结衣高清无吗| 51午夜福利影视在线观看| 亚洲人成电影免费在线| 久久精品国产99精品国产亚洲性色| 亚洲最大成人中文| 精品一区二区免费观看| 一边摸一边抽搐一进一小说| 久久人妻av系列| 国产一区二区在线观看日韩| 国产成人影院久久av| 脱女人内裤的视频| 国产精品影院久久| 久久天躁狠狠躁夜夜2o2o| 久久国产精品人妻蜜桃| 亚洲 国产 在线| 亚洲欧美日韩卡通动漫| 国产大屁股一区二区在线视频| 999久久久精品免费观看国产| 综合色av麻豆| 黄色丝袜av网址大全| 国产一级毛片七仙女欲春2| 国产精品日韩av在线免费观看| 欧美黄色淫秽网站| 欧美色视频一区免费| 久久伊人香网站| 日韩欧美国产一区二区入口| 成人午夜高清在线视频| 国产成人a区在线观看| 最好的美女福利视频网| 夜夜爽天天搞| av在线老鸭窝| 人人妻人人澡欧美一区二区| 性插视频无遮挡在线免费观看| 真实男女啪啪啪动态图| 日韩欧美国产一区二区入口| 一区二区三区免费毛片| 99精品在免费线老司机午夜| 日韩中字成人| 久99久视频精品免费| 激情在线观看视频在线高清| 日韩欧美三级三区| av黄色大香蕉| 久久国产乱子伦精品免费另类| 国模一区二区三区四区视频| 亚洲国产日韩欧美精品在线观看| 久久婷婷人人爽人人干人人爱| 精品无人区乱码1区二区| 中文亚洲av片在线观看爽| 午夜福利免费观看在线| 午夜久久久久精精品| 亚洲成人精品中文字幕电影| 男人和女人高潮做爰伦理| 欧美不卡视频在线免费观看| 国产三级黄色录像| 老女人水多毛片| 欧美bdsm另类| 久久99热6这里只有精品| 国产极品精品免费视频能看的| 69av精品久久久久久| 99久久成人亚洲精品观看| 精品不卡国产一区二区三区| 欧美一区二区亚洲| 欧美中文日本在线观看视频| 欧美日本视频| 亚洲国产欧洲综合997久久,| 亚洲国产精品久久男人天堂| 亚洲精品亚洲一区二区| 欧美在线一区亚洲| 国产亚洲精品综合一区在线观看| 国产私拍福利视频在线观看| 亚洲av电影不卡..在线观看| 3wmmmm亚洲av在线观看| 欧美色欧美亚洲另类二区| 嫩草影视91久久| 欧美乱妇无乱码| 久久久久国产精品人妻aⅴ院| or卡值多少钱| 99精品在免费线老司机午夜| 久久久久久久久大av| 久久草成人影院| 91久久精品电影网| 中国美女看黄片| 内射极品少妇av片p| 国产一区二区激情短视频| 中文字幕免费在线视频6| av中文乱码字幕在线| 麻豆国产av国片精品| 女同久久另类99精品国产91| 天美传媒精品一区二区| 亚洲av成人精品一区久久| 2021天堂中文幕一二区在线观| 欧美色视频一区免费| 小说图片视频综合网站| 免费看美女性在线毛片视频| 国产在线精品亚洲第一网站| 欧美bdsm另类| 亚洲成a人片在线一区二区| 日韩免费av在线播放| 一本一本综合久久| 女同久久另类99精品国产91| 欧美zozozo另类| 欧美乱色亚洲激情| 嫩草影院精品99| 看免费av毛片| 色吧在线观看| 内地一区二区视频在线| 午夜福利18| 亚洲成a人片在线一区二区| 日韩欧美国产一区二区入口| 久9热在线精品视频| 亚洲电影在线观看av| 久久亚洲真实| 成年版毛片免费区| 精品免费久久久久久久清纯| 真人做人爱边吃奶动态| 久久精品综合一区二区三区| 国产aⅴ精品一区二区三区波| 搡女人真爽免费视频火全软件 | 国产精品电影一区二区三区| 波多野结衣高清作品| 国产综合懂色| 日韩欧美一区二区三区在线观看| 色av中文字幕| 深夜a级毛片| 国产视频内射| 国内少妇人妻偷人精品xxx网站| 色精品久久人妻99蜜桃| 国产精品国产高清国产av| 日本a在线网址| 欧美激情国产日韩精品一区| 国产精品,欧美在线| 国产大屁股一区二区在线视频| 九九久久精品国产亚洲av麻豆| 欧美极品一区二区三区四区| 18禁黄网站禁片免费观看直播| 性欧美人与动物交配| 日韩欧美精品v在线| 亚洲成a人片在线一区二区| 久久久久亚洲av毛片大全| 99热6这里只有精品| 国产三级在线视频| 一进一出抽搐动态| 国产爱豆传媒在线观看| 禁无遮挡网站| 动漫黄色视频在线观看| 午夜老司机福利剧场| 国产老妇女一区| 老鸭窝网址在线观看| 日韩精品中文字幕看吧| 69av精品久久久久久| 国产精品乱码一区二三区的特点| 天美传媒精品一区二区| 日韩国内少妇激情av| 欧美成人性av电影在线观看| 免费看光身美女| 欧美在线一区亚洲| 精品福利观看| 老鸭窝网址在线观看| 欧美日韩综合久久久久久 | 亚洲精品日韩av片在线观看| 精品日产1卡2卡| avwww免费| 国产探花在线观看一区二区| 天美传媒精品一区二区| 色5月婷婷丁香| 51午夜福利影视在线观看| 少妇高潮的动态图| 国产亚洲av嫩草精品影院| 一进一出抽搐gif免费好疼| 老司机午夜十八禁免费视频| 男女下面进入的视频免费午夜| 精品久久久久久久久久久久久| 免费av毛片视频| 窝窝影院91人妻| 性色av乱码一区二区三区2| 黄色一级大片看看| www.www免费av| 91久久精品电影网| 亚洲成a人片在线一区二区| 久久精品久久久久久噜噜老黄 | 精品免费久久久久久久清纯| 欧美高清成人免费视频www| 18禁在线播放成人免费| 欧美激情国产日韩精品一区| 亚洲av一区综合|