• <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棚膜拉伸強度估計
    久久久a久久爽久久v久久| 精品少妇黑人巨大在线播放 | 国产伦精品一区二区三区视频9| 欧美性感艳星| 又爽又黄无遮挡网站| 国产色婷婷99| 寂寞人妻少妇视频99o| 国产精品久久久久久精品电影| 少妇裸体淫交视频免费看高清| 一本久久精品| 乱系列少妇在线播放| 精品一区二区三区视频在线| 欧美激情久久久久久爽电影| 小蜜桃在线观看免费完整版高清| 日韩国内少妇激情av| 亚洲人成网站在线播放欧美日韩| 九九在线视频观看精品| av.在线天堂| 亚洲中文字幕日韩| 插逼视频在线观看| 亚洲激情五月婷婷啪啪| 国产精品综合久久久久久久免费| 免费看a级黄色片| 黄片无遮挡物在线观看| av国产免费在线观看| 国产精品综合久久久久久久免费| 99视频精品全部免费 在线| 欧美激情久久久久久爽电影| 成年版毛片免费区| 伦理电影大哥的女人| 中文资源天堂在线| 99热精品在线国产| 午夜视频国产福利| 日韩欧美在线乱码| 伦理电影大哥的女人| 国产亚洲精品久久久久久毛片| 久久久久久久久久成人| 大香蕉久久网| 五月伊人婷婷丁香| 女人十人毛片免费观看3o分钟| 日本黄色片子视频| 不卡视频在线观看欧美| www.色视频.com| 男女视频在线观看网站免费| 亚洲人成网站高清观看| 国产精品综合久久久久久久免费| 青春草亚洲视频在线观看| 人妻少妇偷人精品九色| 91精品一卡2卡3卡4卡| 大型黄色视频在线免费观看| 国产精品三级大全| 国产人妻一区二区三区在| 天堂网av新在线| 中文亚洲av片在线观看爽| 草草在线视频免费看| 日韩av在线大香蕉| 亚洲精品乱码久久久v下载方式| 大香蕉久久网| 亚洲一级一片aⅴ在线观看| 18禁在线无遮挡免费观看视频| 欧美极品一区二区三区四区| 美女脱内裤让男人舔精品视频 | 青春草亚洲视频在线观看| 韩国av在线不卡| 国产精品人妻久久久久久| 久久久午夜欧美精品| a级毛片a级免费在线| 国产成人aa在线观看| 99精品在免费线老司机午夜| 国产三级中文精品| 91久久精品国产一区二区三区| 一进一出抽搐动态| 中文资源天堂在线| 男女做爰动态图高潮gif福利片| 夜夜爽天天搞| 久久久精品大字幕| 1000部很黄的大片| 欧美性猛交╳xxx乱大交人| 日本撒尿小便嘘嘘汇集6| 欧美日韩国产亚洲二区| 国产精品av视频在线免费观看| 性插视频无遮挡在线免费观看| 中文字幕免费在线视频6| 老熟妇乱子伦视频在线观看| 99热只有精品国产| 午夜激情福利司机影院| 国产精品一区二区三区四区免费观看| 久久亚洲精品不卡| 99riav亚洲国产免费| 久久人妻av系列| 成年女人看的毛片在线观看| 国产一区亚洲一区在线观看| 久久久色成人| 日韩成人伦理影院| 婷婷色综合大香蕉| 18禁裸乳无遮挡免费网站照片| 久久久a久久爽久久v久久| 不卡一级毛片| 国产精品女同一区二区软件| 晚上一个人看的免费电影| 黄色配什么色好看| 亚洲av男天堂| 亚洲精品国产av成人精品| 国产色爽女视频免费观看| 国内精品美女久久久久久| 少妇的逼好多水| 午夜激情福利司机影院| 国产精品.久久久| 亚洲欧美成人综合另类久久久 | 99热这里只有精品一区| 亚洲第一区二区三区不卡| 变态另类丝袜制服| 女人被狂操c到高潮| 看免费成人av毛片| 日韩欧美三级三区| 国产麻豆成人av免费视频| 国产伦精品一区二区三区四那| 亚洲激情五月婷婷啪啪| 久久人人爽人人爽人人片va| 男人的好看免费观看在线视频| 国产精品永久免费网站| 日韩欧美在线乱码| avwww免费| 成人无遮挡网站| 国产在线精品亚洲第一网站| 国产av在哪里看| 亚洲av一区综合| 成人亚洲精品av一区二区| 一级黄色大片毛片| 成人特级av手机在线观看| 色视频www国产| 黄色配什么色好看| 久久精品影院6| 精品无人区乱码1区二区| 日韩国内少妇激情av| 搡女人真爽免费视频火全软件| 少妇猛男粗大的猛烈进出视频 | 有码 亚洲区| 免费在线观看成人毛片| 在线a可以看的网站| 欧美精品国产亚洲| 女人被狂操c到高潮| 久久久久久久久久久免费av| 爱豆传媒免费全集在线观看| 99热这里只有是精品在线观看| 国产成人freesex在线| 久久精品国产亚洲av香蕉五月| 禁无遮挡网站| 亚洲成av人片在线播放无| av国产免费在线观看| 午夜老司机福利剧场| 亚洲欧美日韩高清在线视频| 91av网一区二区| 国产精品久久久久久av不卡| a级一级毛片免费在线观看| 深夜a级毛片| 午夜福利在线观看吧| 黄色欧美视频在线观看| 欧美日本亚洲视频在线播放| 69av精品久久久久久| 中文字幕人妻熟人妻熟丝袜美| av在线蜜桃| 噜噜噜噜噜久久久久久91| 极品教师在线视频| 我的老师免费观看完整版| 中文字幕熟女人妻在线| 天堂影院成人在线观看| av在线观看视频网站免费| 毛片一级片免费看久久久久| 国产在线男女| 丰满的人妻完整版| 桃色一区二区三区在线观看| 永久网站在线| 国产黄片视频在线免费观看| 一夜夜www| 亚洲美女搞黄在线观看| 边亲边吃奶的免费视频| 欧美在线一区亚洲| 波多野结衣巨乳人妻| 欧美日韩精品成人综合77777| 最近手机中文字幕大全| 日韩在线高清观看一区二区三区| 99久久中文字幕三级久久日本| 亚洲欧美中文字幕日韩二区| 国产白丝娇喘喷水9色精品| 中文精品一卡2卡3卡4更新| 欧美色欧美亚洲另类二区| 国产精品国产高清国产av| 亚州av有码| 秋霞在线观看毛片| 网址你懂的国产日韩在线| АⅤ资源中文在线天堂| 日本三级黄在线观看| 一区二区三区高清视频在线| 搡老妇女老女人老熟妇| 一级毛片aaaaaa免费看小| 亚洲精品日韩av片在线观看| 国产精品久久久久久精品电影小说 | 嘟嘟电影网在线观看| 欧美区成人在线视频| 久久久精品94久久精品| 日本在线视频免费播放| 免费看a级黄色片| 能在线免费看毛片的网站| 国产成人午夜福利电影在线观看| 丰满的人妻完整版| 久久草成人影院| av.在线天堂| 免费看美女性在线毛片视频| 久久这里有精品视频免费| 18禁黄网站禁片免费观看直播| 性插视频无遮挡在线免费观看| 国产精品野战在线观看| 色噜噜av男人的天堂激情| 国产人妻一区二区三区在| 国产毛片a区久久久久| 国产av不卡久久| 免费看日本二区| 久久精品91蜜桃| 免费大片18禁| 成人欧美大片| 国内揄拍国产精品人妻在线| 熟女人妻精品中文字幕| 婷婷精品国产亚洲av| 色哟哟·www| 精品久久久久久成人av| 欧美精品一区二区大全| 欧美zozozo另类| 青春草视频在线免费观看| 欧美日韩精品成人综合77777| 麻豆成人午夜福利视频| 国产精品国产高清国产av| 精品久久国产蜜桃| 国产 一区 欧美 日韩| 2021天堂中文幕一二区在线观| 日韩欧美精品免费久久| 亚洲精华国产精华液的使用体验 | 国产av麻豆久久久久久久| 亚洲欧洲国产日韩| 久久人人爽人人片av| 日本成人三级电影网站| 欧美性感艳星| 午夜视频国产福利| 欧美精品国产亚洲| 国产又黄又爽又无遮挡在线| 热99在线观看视频| 婷婷色av中文字幕| 亚洲精品久久久久久婷婷小说 | 国产精品美女特级片免费视频播放器| 亚洲人成网站在线播放欧美日韩| 欧美在线一区亚洲| 久久久久久久久久成人| 淫秽高清视频在线观看| 日韩视频在线欧美| 日日撸夜夜添| 一卡2卡三卡四卡精品乱码亚洲| 天天躁夜夜躁狠狠久久av| 国产成人a∨麻豆精品| 人人妻人人看人人澡| 高清日韩中文字幕在线| 午夜精品国产一区二区电影 | 一本久久精品| h日本视频在线播放| 亚洲无线在线观看| 亚洲丝袜综合中文字幕| 日日撸夜夜添| 99久久人妻综合| 日产精品乱码卡一卡2卡三| 99riav亚洲国产免费| 国产探花极品一区二区| 精品久久国产蜜桃| 亚洲三级黄色毛片| 美女黄网站色视频| 国产在视频线在精品| 黄片无遮挡物在线观看| 国产中年淑女户外野战色| 成人特级黄色片久久久久久久| 校园人妻丝袜中文字幕| 尾随美女入室| 狂野欧美激情性xxxx在线观看| 99热这里只有是精品50| 五月伊人婷婷丁香| 成人二区视频| 国产三级在线视频| 美女 人体艺术 gogo| 亚洲第一区二区三区不卡| 少妇猛男粗大的猛烈进出视频 | 少妇猛男粗大的猛烈进出视频 | 美女内射精品一级片tv| 看免费成人av毛片| 激情 狠狠 欧美| 国产成人精品一,二区 | 亚洲人成网站在线观看播放| 美女国产视频在线观看| 特大巨黑吊av在线直播| 国产真实伦视频高清在线观看| 国产精品一区www在线观看| 在线观看免费视频日本深夜| 亚洲五月天丁香| 国产伦在线观看视频一区| 变态另类丝袜制服| 国产探花在线观看一区二区| 五月玫瑰六月丁香| 久久久久久久亚洲中文字幕| 久久久久久久久大av| 亚洲av不卡在线观看| 日产精品乱码卡一卡2卡三| 亚洲国产精品成人久久小说 | 在线免费观看不下载黄p国产| 日本-黄色视频高清免费观看| 丝袜美腿在线中文| 久久鲁丝午夜福利片| 我要看日韩黄色一级片| 欧美日本视频| 舔av片在线| 99精品在免费线老司机午夜| 久久人人爽人人爽人人片va| 99热精品在线国产| 五月玫瑰六月丁香| 精品久久久噜噜| 免费av毛片视频| 国产精品久久电影中文字幕| 91aial.com中文字幕在线观看| 国产美女午夜福利| 人妻夜夜爽99麻豆av| 亚洲va在线va天堂va国产| 网址你懂的国产日韩在线| 最近中文字幕高清免费大全6| 欧美日本视频| 久久综合国产亚洲精品| 伦精品一区二区三区| 激情 狠狠 欧美| 亚洲av中文字字幕乱码综合| 亚洲乱码一区二区免费版| 一进一出抽搐动态| 亚洲精品亚洲一区二区| 热99re8久久精品国产| 国产免费一级a男人的天堂| 最后的刺客免费高清国语| 熟妇人妻久久中文字幕3abv| 色综合亚洲欧美另类图片| 午夜福利在线观看吧| 国产精品无大码| 97热精品久久久久久| 麻豆成人av视频| 国产国拍精品亚洲av在线观看| 精品人妻一区二区三区麻豆| 最好的美女福利视频网| 男人的好看免费观看在线视频| 国产探花极品一区二区| 国产成人精品婷婷| 欧美成人精品欧美一级黄| 久久精品91蜜桃| 两个人视频免费观看高清| 99热这里只有是精品50| 草草在线视频免费看| 国产视频内射| 欧美不卡视频在线免费观看| av在线天堂中文字幕| 99热精品在线国产| 丰满的人妻完整版| 色视频www国产| 午夜福利高清视频| 一边摸一边抽搐一进一小说| 欧美+日韩+精品| 亚洲av男天堂| 国产日本99.免费观看| 成年女人看的毛片在线观看| а√天堂www在线а√下载| 亚洲av成人精品一区久久| 国产真实伦视频高清在线观看| 性色avwww在线观看| 欧美激情国产日韩精品一区| 丰满的人妻完整版| 免费在线观看成人毛片| 国产精品一区二区在线观看99 | 狂野欧美白嫩少妇大欣赏| 成年版毛片免费区| 欧美一区二区国产精品久久精品| 午夜福利成人在线免费观看| 国产乱人视频| 美女内射精品一级片tv| 91精品一卡2卡3卡4卡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 麻豆国产av国片精品| 三级经典国产精品| 我要看日韩黄色一级片| 国产成人精品一,二区 | 日日摸夜夜添夜夜添av毛片| 青青草视频在线视频观看| 欧美日韩国产亚洲二区| 成人午夜高清在线视频| 国产高清不卡午夜福利| 亚洲av免费高清在线观看| 亚洲中文字幕一区二区三区有码在线看| 毛片女人毛片| 精品人妻偷拍中文字幕| 99热这里只有是精品在线观看| 欧美日本亚洲视频在线播放| 女人被狂操c到高潮| 日本色播在线视频| 99热这里只有精品一区| 国产精品免费一区二区三区在线| 国产日韩欧美在线精品| 亚洲欧洲日产国产| 久久久久免费精品人妻一区二区| 亚洲七黄色美女视频| 伦理电影大哥的女人| 美女内射精品一级片tv| 国产一级毛片七仙女欲春2| 麻豆成人午夜福利视频| 99国产极品粉嫩在线观看| 欧美3d第一页| 少妇丰满av| 在线播放国产精品三级| 国产高清视频在线观看网站| 中文字幕av在线有码专区| av在线天堂中文字幕| 国产探花极品一区二区| 午夜福利在线观看吧| av.在线天堂| 亚洲国产精品成人久久小说 | 亚洲久久久久久中文字幕| 精品一区二区免费观看| 午夜精品一区二区三区免费看| 有码 亚洲区| 欧美日本视频| 免费人成视频x8x8入口观看| 成人特级黄色片久久久久久久| 好男人在线观看高清免费视频| 精品人妻偷拍中文字幕| 久久精品国产99精品国产亚洲性色| 少妇被粗大猛烈的视频| 国产av一区在线观看免费| 国产私拍福利视频在线观看| av专区在线播放| 国产精品爽爽va在线观看网站| 国产午夜精品论理片| 欧美激情在线99| 亚洲成人av在线免费| 日本欧美国产在线视频| 神马国产精品三级电影在线观看| 老司机福利观看| 久久99热这里只有精品18| 十八禁国产超污无遮挡网站| 97在线视频观看| 欧美+日韩+精品| 国产成人a∨麻豆精品| 成人特级黄色片久久久久久久| 国产人妻一区二区三区在| 蜜臀久久99精品久久宅男| 性欧美人与动物交配| 久久热精品热| 国产黄片视频在线免费观看| 99视频精品全部免费 在线| 毛片女人毛片| 2022亚洲国产成人精品| 国产精品久久电影中文字幕| 国产精品人妻久久久影院| 波多野结衣巨乳人妻| 午夜亚洲福利在线播放| 在线a可以看的网站| 精品一区二区三区视频在线| 村上凉子中文字幕在线| 国产乱人偷精品视频| 五月玫瑰六月丁香| 日日摸夜夜添夜夜添av毛片| 欧美成人a在线观看| 寂寞人妻少妇视频99o| 中文字幕精品亚洲无线码一区| 国产真实伦视频高清在线观看| 久久精品国产亚洲av天美| 亚洲国产欧美在线一区| 国产精品人妻久久久久久| av在线播放精品| 亚洲成人av在线免费| av女优亚洲男人天堂| 免费看a级黄色片| 蜜桃久久精品国产亚洲av| 国产精品国产三级国产av玫瑰| 亚洲av中文字字幕乱码综合| 一区二区三区四区激情视频 | 一进一出抽搐动态| 亚洲精品成人久久久久久| 级片在线观看| 亚洲国产欧美人成| 九草在线视频观看| 国产亚洲精品久久久com| 神马国产精品三级电影在线观看| 国产精品久久电影中文字幕| 一边亲一边摸免费视频| 天天一区二区日本电影三级| 免费观看的影片在线观看| 天天躁夜夜躁狠狠久久av| 中出人妻视频一区二区| 国产精品,欧美在线| 亚洲人成网站在线观看播放| 一个人观看的视频www高清免费观看| 日本-黄色视频高清免费观看| 麻豆乱淫一区二区| www.色视频.com| 免费一级毛片在线播放高清视频| 女人十人毛片免费观看3o分钟| 国产午夜福利久久久久久| 久久99热6这里只有精品| 成人国产麻豆网| 亚洲在线自拍视频| 免费观看人在逋| 插逼视频在线观看| 欧美色欧美亚洲另类二区| 久久精品国产亚洲av香蕉五月| 国产精品电影一区二区三区| 国产高清有码在线观看视频| 成年女人看的毛片在线观看| 久久精品久久久久久久性| 免费av毛片视频| 麻豆一二三区av精品| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩一区二区视频在线观看视频在线 | 一本精品99久久精品77| 噜噜噜噜噜久久久久久91| 亚洲经典国产精华液单| 亚洲人与动物交配视频| 国产精品av视频在线免费观看| 成人欧美大片| 亚洲欧美日韩卡通动漫| 99久久久亚洲精品蜜臀av| 又爽又黄无遮挡网站| 在线免费观看的www视频| 亚洲欧美成人精品一区二区| 国产精品久久视频播放| 神马国产精品三级电影在线观看| 国产av一区在线观看免费| a级毛片a级免费在线| 免费av不卡在线播放| 人妻少妇偷人精品九色| 日韩精品有码人妻一区| 欧美潮喷喷水| 国产成人a∨麻豆精品| 在线观看av片永久免费下载| 九九爱精品视频在线观看| 亚洲国产欧美人成| 国产精品一区二区在线观看99 | 久久久久久国产a免费观看| 伦理电影大哥的女人| 成人av在线播放网站| 免费观看a级毛片全部| 久久久久久久亚洲中文字幕| 日本av手机在线免费观看| 色哟哟哟哟哟哟| 晚上一个人看的免费电影| АⅤ资源中文在线天堂| 少妇高潮的动态图| 国产麻豆成人av免费视频| 日韩制服骚丝袜av| 18+在线观看网站| 好男人视频免费观看在线| 中文字幕av成人在线电影| 日本在线视频免费播放| 91久久精品国产一区二区三区| 国产高清不卡午夜福利| 夜夜夜夜夜久久久久| 免费电影在线观看免费观看| 特级一级黄色大片| 久久午夜亚洲精品久久| 久久精品国产自在天天线| 免费一级毛片在线播放高清视频| av卡一久久| 亚洲精品日韩在线中文字幕 | 狂野欧美激情性xxxx在线观看| 在线a可以看的网站| 青春草国产在线视频 | 女人被狂操c到高潮| 国产精品一区二区三区四区久久| 看免费成人av毛片| 国产美女午夜福利| 99热网站在线观看| 特级一级黄色大片| 成人av在线播放网站| 婷婷六月久久综合丁香| 日日摸夜夜添夜夜添av毛片| 99久国产av精品国产电影| 黄片wwwwww| 国产欧美日韩精品一区二区| 毛片一级片免费看久久久久| 国产av在哪里看| 国产日韩欧美在线精品| 丝袜美腿在线中文| 欧美又色又爽又黄视频| 国产色爽女视频免费观看| 日日干狠狠操夜夜爽| 久久精品人妻少妇| 精品人妻一区二区三区麻豆| 日韩成人av中文字幕在线观看| 亚洲成人久久性| 色尼玛亚洲综合影院| 欧美一区二区亚洲| 麻豆久久精品国产亚洲av| 欧洲精品卡2卡3卡4卡5卡区| 国产视频首页在线观看| 插逼视频在线观看| 亚洲欧美精品综合久久99| 麻豆国产97在线/欧美| 小说图片视频综合网站| 国产一级毛片在线| 精品久久久久久久久av| 精品不卡国产一区二区三区| 嘟嘟电影网在线观看| 国国产精品蜜臀av免费| 九九在线视频观看精品| 如何舔出高潮| 国产精品爽爽va在线观看网站| 国产极品精品免费视频能看的| 中文字幕制服av| 日韩强制内射视频|