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

    尺度效應對多光譜遙感監(jiān)測土壤水分影響研究

    2023-03-04 06:28:08羅亞東郭宇宏孫智鵬金煜龍陳俊英余衛(wèi)華
    節(jié)水灌溉 2023年2期
    關鍵詞:特征模型

    羅亞東,許 齊,郭宇宏,2,孫智鵬,金煜龍,陳俊英,2,余衛(wèi)華,2

    (1.西北農林科技大學水利與建筑工程學院,陜西 楊凌 712100;2.西北農林科技大學 旱區(qū)農業(yè)水土工程教育部重點實驗室,陜西 楊凌 712100)

    0 引言

    土壤含水率是分析冬小麥生育情況的主要依據,對于理解和解決農業(yè)生產、生態(tài)規(guī)劃以及水資源管理中的科學與實際問題至關重要[1]。傳統(tǒng)的土壤含水率監(jiān)測方法操作復雜、且數據滯后。遙感技術因其快速、簡便以及無損的特點,已廣泛應用于農田農情的監(jiān)測[2-7]?;谶b感的墑情監(jiān)測是農村水利科技工作的重點領域[8]。但是在無人機多光譜遙感實踐中,存在尺度效應影響監(jiān)測精度的問題。

    近年來,遙感監(jiān)測已廣泛應用于土壤含水率的研究,楊珺博[9]等人利用無人機遙感監(jiān)測數據反演土壤含水率,發(fā)現(xiàn)反演精度可達80%以上,且反演效果最好的土壤深度為10~20 cm。Cheng Minghan[10]等利用無人機多模態(tài)數據對玉米田間含水率進行4種模型建立與反演,得出多光譜數據在不同模型下均提供最精確的土壤水分估算。

    尺度效應是指在無人機遙感觀測中,隨著遙感分辨率的變化,遙感反演得到的數據與實際表現(xiàn)不一致的現(xiàn)象,作為遙感科學的核心問題,近年來受到了學者的廣泛關注[11,12]。歐盛華[13]等利用機載航空WIDAS 數據,發(fā)現(xiàn)尺度效應是空間異質性和反演函數非線性的綜合體現(xiàn),空間異質性和反演函數非線性程度的增加均會造成尺度效應的增大。徐凱健[14]通過對多光譜影像的森林樹種識別及其空間尺度響應的研究。結果表明,影像空間分辨率對區(qū)域樹種識別結果具有顯著影響。Zhang[15]等利用偏振遙感器,根據地物目標的空間異質性和成像分辨率帶來的限制,發(fā)現(xiàn)了尺度效應對地物偏振圖像的影響。Tang[16]等對SEBAL 模型的空間尺度效應進行探究,得出高低分辨率對觀測熱通量的影響??梢?,尺度效應對于無人機遙感觀測有一定影響,但目前尺度效應對于無人機遙感監(jiān)測土壤水分的研究還鮮有報道,在實踐中缺乏最佳的監(jiān)測尺度。在提取作物冠層信息過程中,土壤、陰影背景均會對冠層信息準確性造成影響[17,18]。楊帥[19]等基于RGRI指數法獲取的玉米冠層溫度與實測冠層溫度的相關性,得出不剔除土壤背景效果最差。王佳兒[20]等通過剔除土壤背景反演玉米田間含水率,發(fā)現(xiàn)數據的相關性和回歸模型的精度始終比未剔除土壤背景數據的相關性和模型精度差。紋理特征在冠層信息提取中中是一種相當重要的特征,它反映了圖像或物體本身的屬性[21]。王思宇[22]利用紋理特征為基礎,基于多光譜影像并構建玉米冠層FAPAR 估算模型,發(fā)現(xiàn)其精度較高。目前仍缺乏通過去土壤去陰影掩膜方法提取作物冠層信息。

    本文通過無人機搭載六波段多光譜相機在不同高度下獲取冬小麥紋理特征,剔除土壤背景、陰影背景,在此基礎上研究土壤含水率和植被紋理特征之間的關系。探究在不同分辨率下無人機遙感數據與土壤含水率的相關性并分析尺度效應對于監(jiān)測數據的影響。最后得出尺度效應對無人機多光譜遙感監(jiān)測土壤含水率的影響,獲取無人機監(jiān)測土壤水分的最佳飛行高度。

    1 材料與方法

    1.1 研究區(qū)域概況

    研究區(qū)域位于陜西省楊陵區(qū)陜西楊凌西北農林科技大學旱區(qū)農業(yè)水土工程教育部重點實驗室的灌溉試驗站(34°17'N,108°04'E,海拔521 m)。該地區(qū)為典型的溫帶大陸性季風季候,夏季炎熱,冬季寒冷,蒸發(fā)強度較大。年平均氣溫在12 ℃左右,多年平均降水量為630 mm,年平均蒸發(fā)量為884 mm,土壤干容重為1.38 g/cm3。

    1.2 試驗材料與設計

    本試驗冬小麥選用小偃22 號,于2020年10月20日進行人工播種,共播種7.68 kg種子,并且施用20 kg的復合肥料作為底肥。同時在其生長過程中嚴格控制水分條件,試驗田灌水上限分別采用田間持水量的50%、65%、80%和95%,設置4 個水分處理,3 個重復,共12 個小區(qū),以獲得不同田持下的土壤含水率數據,使數據更有代表性。每個小區(qū)的面積為4 m×4 m=16 m2,每個小區(qū)間隔150 cm,相鄰兩個小區(qū)之間設置保護行。小區(qū)內種植的冬小麥行距為25 cm,每一行播種40 g種子。為減少外界環(huán)境的干擾,每個小區(qū)都配有擋雨棚,圖1為試驗小區(qū)布置。

    圖1 實驗小區(qū)布置Fig.1 Layout of experimental plot

    1.3 數據采集

    1.3.1 土壤水分獲取

    分別在3 個生育期在在每個實驗小區(qū)取3 個采樣點A、B、C,在采樣點處用土鉆采集10、20、30、40、60 cm 深度的土樣,每個土樣30 g 左右,土樣采出后立即裝入鋁盒中。將采集完成的鋁盒放入干燥箱,在105 ℃干燥后稱重并計算土壤含水率。烘干法測定土壤質量含水率的計算式為:

    式中:W1為鋁盒加濕土壤質量;W2為鋁盒加干土質量;W3為空鋁盒質量。

    1.3.2 無人機平臺與冠層數據獲取

    在試驗田內選取晴朗天氣下的土壤樣本,以配合無人機多光譜遙感圖像的獲取。采集時間為冬小麥分蘗期(2020年12月24日)、拔節(jié)期(2021年3月23日)、抽穗期(2021年4月19日),每日11∶00-14∶00。

    本試驗采用由深圳市大疆創(chuàng)新科技有限公司生產的經緯Matrice 600 六旋翼無人機,搭載Micro-MCA 多光譜相機進行多光譜遙感圖像采集。Micro-MCA 相機共有6個鏡頭,分別對應6 個波段,波長分別為900 nm(近紅外)、800 nm(近紅外)、720 nm(紅邊)、680 nm(紅光)、550 nm(綠光)、490 nm(藍光)。無人機飛行高度設置為19、37、55、74 和92 m,對應的分辨率分別為10、20、30、40 和50 mm,鏡頭垂直向下采集光譜信息。拍攝之前在實驗區(qū)附近放置白板和黑板,用于研究區(qū)內進行標定,從而獲取冬小麥光譜反射率。

    1.4 多光譜圖像處理

    使用多光譜相機完成圖像采集后,使用PixelWrench2[23]軟件對分別對分蘗期、拔節(jié)期、抽穗期無人機多光譜遙感原始圖像進行預處理,將預處理后遙感圖像導入Pix4Dmapper 軟件,進行幾何校正、輻射定標、試驗區(qū)拼接處理,得到試驗區(qū)多光譜遙感拼接圖像。

    將拼接處理后的圖像導入ENVI5.3(64 bit)軟件中,剔除土壤和陰影背景。使用ENVI5.3(64 bit)軟件中Co-occurrence measure 工具,獲取多光譜圖像6 個波段基于二階概率統(tǒng)計濾波的均值、方差、協(xié)同性、對比度、相異性、信息熵、二階矩、相關性等八種紋理特征。即從每個多光譜圖像中提取48組紋理特征。

    1.5 灰色關聯(lián)法優(yōu)選紋理特征

    不同冠層光譜反射率和紋理特征之間差異較大,存在較強的共線性問題[23]。從理論上來[24]說,預測指標與參考變量之間的相關性越強,其預測結果越可靠[25]。采用全變量建模時,會導致模型預測不準確。

    基于灰色系統(tǒng)理論的灰色關聯(lián)分析是一種利用比較序列和參考序列之間的相似度來判別各變量相關性強弱的方法,其模型準確度高、數據流失性小[26,27]。

    本文使用灰色關聯(lián)法研究步驟如下:

    構建參考序列:

    構建比較序列:

    計算關聯(lián)系數:

    計算關聯(lián)度:

    式中:y表示土壤含水率;xa、xb表示一個生育期下的紋理特征;ξi代表關聯(lián)系數;si代表關聯(lián)度;n表示3 個生育期下12塊小區(qū)各自對應的土壤含水率,n=36;k為6 個波段下的的紋理特征總數,k=48;n為樣本組數,n=36;i表示3 組掩膜處理,i=3。

    1.6 反演模型構建

    本研究每組樣本經灰色關聯(lián)法優(yōu)選后為8個輸入變量,采用多元線性回歸(MLR)、前饋神經網絡(BPNN)和隨機森林(RF)3種回歸方法在5 個高度和3 種掩膜方法下分別建立土壤水分反演模型,共45個土壤水分模型。將36組樣本按2∶1的比例隨機取樣作為建模集和驗證集。即選取24 個樣本用于建模,其余12個用于驗證。

    MLR 利用線性模型來擬合多個自變量和因變量的關系,從而確定多元線性回歸模型的參數,回歸至原假設方程中,通過回歸方程來預測因變量的趨勢[28-30]。

    BPNN 作為一種根據誤差逆?zhèn)鞑ニ惴ㄓ柧毜亩啻吻梆伨W絡,它通過對人類神經元的功能進行模擬,可以儲存及學習大量的輸入數據和輸出數據,且不需對變量的映射關系進行描述,利用輸入和輸出數據建模,其對非線性系統(tǒng)具有很強的模擬能力[31,32]。本文搭建4 層BP 神經網絡,采用兩個隱含層,傳遞函數分別為雙曲正切S型函數和線性函數,節(jié)點數分別為8和2。

    RF 是一種基于分類樹的算法,通過對大量分類樹的匯總提高了模型的預測精度,是取代神經網絡等傳統(tǒng)機器學習方法的新的模型。當林中的樹數變大時,林的泛化誤差會收斂到極限。可以體現(xiàn)變量間的交互作用,具有很高的預測準確率,對異常值和噪聲具有很好的容忍度,且不容易出現(xiàn)過擬合[33,34]。

    構建反演模型,比較3種掩膜方法下各反演模型的擬合優(yōu)度R2和均方根誤差RMSE,分析尺度效應對于無人機遙感監(jiān)測土壤水分的影響。R2越接近1,RMSE越接近0,則反演模型的精度和準確性越好[23]。

    2 結果與分析

    2.1 土壤含水率

    通過取土烘干法獲取的不同生育期下每個小區(qū)土壤含水率梯度值,結果如圖2所示。

    從圖2可以看出,不同生育期的土壤含水率有所不同,各小區(qū)分別在分蘗期、拔節(jié)期和抽穗期的土壤含水率依次降低,土壤含水率的覆蓋范圍較大且具有代表性。

    圖2 各生育期各小區(qū)土壤含水率Fig.2 SWC of each plot at each fertility stage

    2.2 灰色關聯(lián)法處理結果

    獲得3組不同掩膜方法下的紋理特征,每組數據統(tǒng)計各自48 個紋理特征。本文通過灰色關聯(lián)法結果顯示紋理特征與土壤含水率的相關度,進而對各分辨率與掩膜方法進行初步評價與分析。

    使用灰色關聯(lián)法計算紋理特征與土壤含水率相關度,結果如圖3所示。

    圖3 灰色關聯(lián)法處理結果Fig.3 Results of the GRA method

    無掩膜(NM)處理中,10 mm 和30 mm 下紋理特征與土壤含水率相關性最顯著,20 mm和40 mm下紋理特征與土壤含水率相關性總體最佳,關聯(lián)度分布在0.7~0.8。去土壤掩膜(SRM)處理中,紋理特征與土壤含水率的相關性變化不大。去土壤去陰影掩膜(SSRM)處理中,總體相關性較差,分辨率為10 mm時紋理特征與土壤含水率相關性最顯著,分辨率為40 mm 時次之。

    從圖3可以看出,在10 mm 下,紋理特征與土壤含水率呈現(xiàn)最高相關性,其次則為30 mm、20 mm、40 mm、50 mm。這是由于飛行高度提高,分辨率隨之降低,混合像元在總像元中的比例上升,對冬小麥紋理特征造成影響,紋理特征與土壤含水率的相關性下降。

    在掩膜分析結果中,無掩膜與去土壤掩膜的結果差異較小,且優(yōu)于去土壤去陰影掩膜。本文認為在進行去土壤掩膜處理時,混合像元容易分類,而陰影和冠層的差異較??;而在進行去土壤去陰影處理時,混合像元難以分類,導致過多地保留陰影或刪除必要的冠層像元導致紋理特征不能準確反映冠層情況。

    2.3 三種掩膜方法與不同分辨率相關性分析

    2.3.1 無掩膜處理下反演結果

    圖4可以看出,多元線性回歸模型在分辨率為10 mm 時擬合優(yōu)度(R2)最高,其次是40 mm。BP 神經網絡模型在分辨率為10 mm 時R2最高,30 mm 次之。隨機森林模型在分辨率為10 mm時R2較高,其次是30 mm。

    圖4 無掩膜處理反演結果Fig.4 Regression results without mask

    在無掩膜處理中,3 種回歸方法的最佳反演結果均在10 mm 時取得較好結果,擬合優(yōu)度在0.7 以上,其中BP 神經網絡的擬合優(yōu)度為0.86,RMSE為0.011。值得注意的是,BP 神經網絡和隨機森林兩種方法在30 mm 分辨率時R2超過了10 mm,因此在10 mm和30 mm分辨率下,冠層情況在紋理特征中得到準確表達。

    同時,3組回歸模型中擬合優(yōu)度總體上呈現(xiàn)出先降低再升高再降低的趨勢,在分辨率降低到30 mm 至40 mm 時R2出現(xiàn)升高趨勢。在此掩膜處理下,多元線性回歸模型與BP 神經網絡的擬合優(yōu)度低于隨機森林模型。

    2.3.2 去土壤掩膜處理下反演結果

    將篩選后去土壤掩膜處理的36 組濾波信息分別使用多元線性回歸、BP 神經網絡、隨機森林進行反演處理,結果如圖5所示。

    由圖5得出,多元線性回歸模型在分辨率為40 mm時R2最高,其次是30 mm。BP 神經網絡模型在分辨率為30 mm 時R2最高,其次是40 mm。隨機森林模型在分辨率為10 mm,40 mm次之。

    圖5 去土壤掩膜處理反演結果Fig.5 Regression results with soil background removed by masking

    在去土壤掩膜處理中,3種回歸方法的最佳反演結果不盡相同,總體分布顯示在該掩膜處理下,40 mm分辨率下冠層情況在紋理特征中得到準確表達。

    3組回歸模型中擬合優(yōu)度總體上呈現(xiàn)出先降低再升高再降低的趨勢,除10 mm 分辨率的結果較好之外,30 mm 與40 mm分辨率下的反演結果擬合優(yōu)度較高。在去土壤掩膜處理之下,3種模型的擬合優(yōu)度均在0.6~0.7之間。

    2.3.3 去土壤去陰影掩膜處理下反演結果

    將篩選后去土壤去陰影掩膜處理的36 組濾波信息分別使用多元線性回歸、BP 神經網絡、隨機森林進行反演處理,結果如圖6所示。

    圖6 去土壤去陰影掩膜處理反演結果Fig.6 Regression results with soil and shadow background removed by masking

    圖6可以看出,在去土壤去陰影掩膜處理中,多元線性回歸模型在分辨率為10 mm 時R2最高,其次是20 mm。BP 神經網絡模型在分辨率為10 mm時R2最高,其次是40 mm。隨機森林模型在分辨率為10 mm,40 mm次之。

    在此種掩膜處理中,3 種回歸方法的擬合結果顯示,在10 mm與40 mm分辨率下取得相對較高的反演精度,見圖7。

    圖7 各掩膜方法下擬合優(yōu)度對比Fig.7 Comparison of the goodness of fit under different masking method

    3組回歸模型中擬合優(yōu)度總體上呈現(xiàn)出逐漸降低趨勢,值得注意的是,在去土壤去陰影掩膜處理之下,除10 mm 的擬合優(yōu)度高于0.7,其余分辨率在各自模型下的擬合優(yōu)度均在0.5左右。

    2.4 討論與分析

    分析反演結果可以得出在分辨率為10 mm和40 mm時,模型反演效果較好。本文認為飛行高度較低時,遙感圖像分辨率較高,能較為精準地表現(xiàn)冠層狀況,同樣使得冠層情況也能在紋理特征中得到準確表達;分辨率為40 mm 時,反演效果僅次于最高分辨率10 mm,本文認為該現(xiàn)象由于是混合像元在總像元中占比降低,當分辨率由10 mm下降至30 mm時,混合像元在在總像元中的比重增加,在用閾值進行背景分割時,會保留部分混合像元中的背景,或將部分混合像元中的作物剔除,導致紋理特征不能準確表達冠層信息。當分辨率下降至40 mm 時,分辨率與葉面寬度相近,此時混合像元占比反而下降,從而獲得較為準確的反演精度。由于遙感圖像中存在土壤、陰影和田間雜草等背景,在用紋理特征進行建模時作為無關變量干擾了模型精度。導致回歸模型的魯棒性較差。當高度增加,分辨率降低至與作物葉面尺度相近時,混合像元中背景為優(yōu)勢類別的像元被完全剔除,冠層為優(yōu)勢類別的像元被保留,并且較少的背景像元與冠層像元混合后對紋理特征的影響降低,能一定程度上減小背景對回歸模型的干擾。

    最后,在各掩膜處理下,3種建模方法的反演結果不盡相同。在無掩膜處理下,多元線性回歸的擬合精度最差;在去土壤掩膜處理下,隨機森林的擬合精度最差;而在去土壤陰影掩膜處理中,BP 神經網絡的擬合精度最差。本文認為原因如下:多元線性回歸對數據有比較嚴格的要求如:①隨機誤差項是一個期望值為0 的隨機變量;②解釋變量的所有觀測值,隨機誤差項有相同的方差;③隨機誤差項彼此不相關等。在無掩膜處理的紋理特征中,干擾值很多隨機誤差不可控,所以這些條件在無掩膜的情況下很難滿足。隨機森林在回歸分析中不能連順序輸出,也不能預測超出訓練集范圍的數據,導致在出現(xiàn)特定噪聲的數據集中建模時容易出現(xiàn)過擬合,在去土壤掩膜處理的紋理特征數據中,出現(xiàn)了陰影的噪聲數據,對隨機森林回歸造成了干擾,并且在用支持向量機進行監(jiān)督分類時,部分土壤并沒有被掩膜處理,也對隨機森林回歸造成了干擾。而BP 神經網絡存在局部極小化的問題,采取局部搜索的方法的容易陷入局部最小值, 并且BP神經網絡的結構自由度很大,目前并沒有完整的理論指導,因此本文采用的隱含層結構可能較適合于有完整信息的數據集,而對噪聲的容忍能力較高,去土壤去陰影掩膜刪去過多有效信息,所以在去土壤去陰影掩膜處理的紋理特征中效果較差。

    3 結論

    (1)本文通過對不同分辨率下獲得的紋理特征進行土壤含水率的反演,發(fā)現(xiàn)在拍攝分辨率為10、40 mm 時,反演模型的擬合優(yōu)度較好。但分辨率在40 mm 時對不同反演方法表現(xiàn)更穩(wěn)定。在實際應用中,過高分辨率下工作效率低,不能大范圍進行遙感監(jiān)測,為工作帶來不便。故本文認為無人機遙感監(jiān)測土壤水分的最佳分辨率為40 mm。

    (2)30 mm-NM-RF 模型為改組數據反演冬小麥土壤含水率的最佳反演模型,但在30 mm 分辨率下其他模型反演效果不佳。

    (3)對于3 種反演模型,表現(xiàn)最好的是NM-RF 模型,擬合優(yōu)度較高且在各分辨率下表現(xiàn)穩(wěn)定。綜合各掩膜處理方法表現(xiàn)最好反演模型是BPNN 模型,在大部分情況下表現(xiàn)穩(wěn)定有較好的反演效果,為機器學習在遙感監(jiān)測的應用驗證了可行性。

    (4)無掩膜處理和去土壤掩膜處理的結果優(yōu)于去土壤去陰影掩膜,且去土壤去陰影掩膜處理下分辨率較低時,反演效果較差。

    猜你喜歡
    特征模型
    一半模型
    抓住特征巧觀察
    重要模型『一線三等角』
    新型冠狀病毒及其流行病學特征認識
    重尾非線性自回歸模型自加權M-估計的漸近分布
    如何表達“特征”
    不忠誠的四個特征
    當代陜西(2019年10期)2019-06-03 10:12:04
    抓住特征巧觀察
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉換方法初步研究
    精品第一国产精品| 肉色欧美久久久久久久蜜桃| 天天操日日干夜夜撸| 亚洲av国产av综合av卡| 大香蕉久久成人网| 精品久久蜜臀av无| 啦啦啦在线免费观看视频4| 侵犯人妻中文字幕一二三四区| 亚洲欧美色中文字幕在线| 人妻 亚洲 视频| 狠狠婷婷综合久久久久久88av| 久久久久国内视频| 国产片内射在线| 精品少妇内射三级| 激情视频va一区二区三区| 久热这里只有精品99| 肉色欧美久久久久久久蜜桃| 国产单亲对白刺激| 成年人午夜在线观看视频| 涩涩av久久男人的天堂| 肉色欧美久久久久久久蜜桃| 2018国产大陆天天弄谢| 变态另类成人亚洲欧美熟女 | 久久久久视频综合| 国产成人av激情在线播放| 日本一区二区免费在线视频| 久久精品国产亚洲av香蕉五月 | 黄网站色视频无遮挡免费观看| 丝袜在线中文字幕| 妹子高潮喷水视频| 国产成人啪精品午夜网站| 亚洲精品国产色婷婷电影| 又黄又粗又硬又大视频| 黄色 视频免费看| videos熟女内射| 另类亚洲欧美激情| 国产有黄有色有爽视频| 亚洲色图 男人天堂 中文字幕| 国产精品久久电影中文字幕 | 热99re8久久精品国产| 丰满迷人的少妇在线观看| 国产精品熟女久久久久浪| 精品国产亚洲在线| 日韩精品免费视频一区二区三区| 日韩三级视频一区二区三区| 国产亚洲精品久久久久5区| 亚洲国产欧美一区二区综合| 久久亚洲真实| 老司机午夜十八禁免费视频| 国产片内射在线| 9191精品国产免费久久| 又紧又爽又黄一区二区| 老司机靠b影院| 两个人免费观看高清视频| 国产精品av久久久久免费| 青草久久国产| 精品午夜福利视频在线观看一区 | 精品国产亚洲在线| 久久精品aⅴ一区二区三区四区| 一个人免费在线观看的高清视频| 久久精品成人免费网站| 亚洲人成伊人成综合网2020| 一级片'在线观看视频| 国产成人免费观看mmmm| 91麻豆av在线| 色视频在线一区二区三区| 国产亚洲精品一区二区www | 午夜两性在线视频| 亚洲伊人色综图| 不卡一级毛片| 欧美av亚洲av综合av国产av| 国产欧美日韩精品亚洲av| 久久亚洲真实| 国产区一区二久久| 色尼玛亚洲综合影院| 中文字幕高清在线视频| 男人舔女人的私密视频| 国产有黄有色有爽视频| √禁漫天堂资源中文www| 国产精品久久久久成人av| 天天躁狠狠躁夜夜躁狠狠躁| 免费看十八禁软件| 黄色怎么调成土黄色| 飞空精品影院首页| 国产有黄有色有爽视频| 欧美另类亚洲清纯唯美| 另类精品久久| 99re6热这里在线精品视频| 亚洲熟女精品中文字幕| 女人爽到高潮嗷嗷叫在线视频| 无遮挡黄片免费观看| 免费女性裸体啪啪无遮挡网站| 亚洲国产欧美在线一区| 黑人巨大精品欧美一区二区mp4| 一进一出好大好爽视频| 最黄视频免费看| 日日摸夜夜添夜夜添小说| 韩国精品一区二区三区| 国产成人精品久久二区二区91| 国产亚洲午夜精品一区二区久久| 亚洲国产欧美日韩在线播放| 欧美激情高清一区二区三区| 黑人猛操日本美女一级片| 一级a爱视频在线免费观看| 在线永久观看黄色视频| 免费在线观看黄色视频的| 久久青草综合色| 久久人人97超碰香蕉20202| 露出奶头的视频| 日本wwww免费看| 日韩免费高清中文字幕av| 丁香欧美五月| 亚洲精品自拍成人| 在线观看免费高清a一片| 乱人伦中国视频| 国内毛片毛片毛片毛片毛片| 亚洲精品国产色婷婷电影| 另类精品久久| 国产亚洲精品久久久久5区| 美女扒开内裤让男人捅视频| 麻豆乱淫一区二区| 国产一区二区三区在线臀色熟女 | 国产欧美日韩综合在线一区二区| 9191精品国产免费久久| 少妇被粗大的猛进出69影院| 国产色视频综合| 脱女人内裤的视频| 国产精品麻豆人妻色哟哟久久| 久久久久久久精品吃奶| 精品国内亚洲2022精品成人 | 男女无遮挡免费网站观看| 亚洲一区二区三区欧美精品| 成人精品一区二区免费| 日韩精品免费视频一区二区三区| 亚洲人成电影免费在线| 两性午夜刺激爽爽歪歪视频在线观看 | av电影中文网址| 午夜福利欧美成人| 亚洲一卡2卡3卡4卡5卡精品中文| 老司机亚洲免费影院| 久久精品国产亚洲av香蕉五月 | 看免费av毛片| 日本欧美视频一区| 如日韩欧美国产精品一区二区三区| 国产1区2区3区精品| 国产男靠女视频免费网站| 亚洲av日韩在线播放| 成人影院久久| 午夜免费成人在线视频| 国产黄色免费在线视频| 精品卡一卡二卡四卡免费| 中文字幕制服av| 午夜福利乱码中文字幕| 国产成人av教育| 黑人欧美特级aaaaaa片| 国产成人欧美| 亚洲av日韩在线播放| 久久影院123| 狠狠精品人妻久久久久久综合| 久久中文字幕一级| 久久久久久久久免费视频了| 极品少妇高潮喷水抽搐| 免费人妻精品一区二区三区视频| av国产精品久久久久影院| 亚洲成人免费av在线播放| 夜夜爽天天搞| 人人妻人人添人人爽欧美一区卜| 国产精品影院久久| 纵有疾风起免费观看全集完整版| 久久人妻av系列| 天天躁夜夜躁狠狠躁躁| 狠狠狠狠99中文字幕| 制服人妻中文乱码| 日韩免费av在线播放| 久久精品国产亚洲av高清一级| 国产高清国产精品国产三级| 伊人久久大香线蕉亚洲五| 亚洲av第一区精品v没综合| 久久ye,这里只有精品| 国产欧美日韩一区二区三区在线| 久久人妻熟女aⅴ| 国产高清国产精品国产三级| 大型黄色视频在线免费观看| 狠狠狠狠99中文字幕| 精品一区二区三区av网在线观看 | 亚洲五月婷婷丁香| 精品国产国语对白av| 午夜福利一区二区在线看| 国产熟女午夜一区二区三区| 一本久久精品| 蜜桃在线观看..| 99热国产这里只有精品6| www日本在线高清视频| 欧美激情 高清一区二区三区| 国产精品美女特级片免费视频播放器 | av有码第一页| 午夜激情久久久久久久| 亚洲专区国产一区二区| 男女床上黄色一级片免费看| 亚洲伊人久久精品综合| 少妇 在线观看| 亚洲av日韩精品久久久久久密| 成年女人毛片免费观看观看9 | 国产精品影院久久| 麻豆av在线久日| 日日摸夜夜添夜夜添小说| av欧美777| 久久久国产成人免费| 一本大道久久a久久精品| 亚洲精品一二三| 两人在一起打扑克的视频| 国产精品免费视频内射| 欧美激情极品国产一区二区三区| 精品久久蜜臀av无| 久久婷婷成人综合色麻豆| 亚洲av第一区精品v没综合| 高潮久久久久久久久久久不卡| 丝袜喷水一区| 制服诱惑二区| 五月天丁香电影| 天堂动漫精品| 老司机福利观看| 国产欧美日韩一区二区三区在线| cao死你这个sao货| 黄色毛片三级朝国网站| 人妻 亚洲 视频| 极品人妻少妇av视频| 亚洲精品成人av观看孕妇| 色94色欧美一区二区| 成人18禁在线播放| 老熟女久久久| 国产黄色免费在线视频| av电影中文网址| 99re6热这里在线精品视频| 777米奇影视久久| 丝袜在线中文字幕| 交换朋友夫妻互换小说| 窝窝影院91人妻| 国产淫语在线视频| bbb黄色大片| av不卡在线播放| 99在线人妻在线中文字幕 | 国产老妇伦熟女老妇高清| 老汉色av国产亚洲站长工具| 人妻久久中文字幕网| 久久久久久人人人人人| 久久久久久久久久久久大奶| 精品少妇内射三级| 欧美激情 高清一区二区三区| 久久久久视频综合| 在线亚洲精品国产二区图片欧美| 欧美精品一区二区大全| 美女福利国产在线| 久久久欧美国产精品| 夜夜夜夜夜久久久久| 天堂俺去俺来也www色官网| 91av网站免费观看| xxxhd国产人妻xxx| 精品视频人人做人人爽| 极品人妻少妇av视频| 国产无遮挡羞羞视频在线观看| 久久99一区二区三区| 中文字幕av电影在线播放| 80岁老熟妇乱子伦牲交| 99热网站在线观看| 999久久久精品免费观看国产| 亚洲欧美一区二区三区黑人| 国产av国产精品国产| 欧美午夜高清在线| 久久久久国产一级毛片高清牌| 亚洲专区中文字幕在线| 视频在线观看一区二区三区| 成人永久免费在线观看视频 | 亚洲自偷自拍图片 自拍| 成年人免费黄色播放视频| 黑人猛操日本美女一级片| 精品亚洲成国产av| av福利片在线| 少妇 在线观看| 天堂中文最新版在线下载| 日日夜夜操网爽| 久久天堂一区二区三区四区| 可以免费在线观看a视频的电影网站| 久久久久久久大尺度免费视频| 最新的欧美精品一区二区| 免费日韩欧美在线观看| 成人免费观看视频高清| 91大片在线观看| 国产淫语在线视频| 欧美黄色淫秽网站| 天天躁夜夜躁狠狠躁躁| 久久精品亚洲av国产电影网| 午夜激情av网站| 欧美日韩亚洲综合一区二区三区_| xxxhd国产人妻xxx| 91麻豆精品激情在线观看国产 | 精品国产一区二区久久| 美女视频免费永久观看网站| av超薄肉色丝袜交足视频| 国产成人精品久久二区二区免费| 99国产精品免费福利视频| 成人18禁在线播放| 免费少妇av软件| 人人妻人人澡人人看| 日本黄色视频三级网站网址 | 伦理电影免费视频| 国产又爽黄色视频| 最新在线观看一区二区三区| 日韩免费av在线播放| 香蕉国产在线看| 亚洲精品一卡2卡三卡4卡5卡| 日本a在线网址| 搡老熟女国产l中国老女人| 亚洲精品粉嫩美女一区| 天堂8中文在线网| 日韩视频在线欧美| 99热网站在线观看| 好男人电影高清在线观看| 国产色视频综合| 亚洲五月色婷婷综合| 天天操日日干夜夜撸| 亚洲自偷自拍图片 自拍| 黑丝袜美女国产一区| 91精品国产国语对白视频| 99国产精品免费福利视频| 国产日韩一区二区三区精品不卡| av在线播放免费不卡| 成人国产av品久久久| 精品熟女少妇八av免费久了| 日韩一区二区三区影片| 电影成人av| 99国产精品一区二区三区| 丝袜美腿诱惑在线| 欧美日韩国产mv在线观看视频| 国产精品.久久久| 淫妇啪啪啪对白视频| 国产免费av片在线观看野外av| 91麻豆av在线| 亚洲国产毛片av蜜桃av| e午夜精品久久久久久久| 999久久久国产精品视频| 欧美人与性动交α欧美精品济南到| 女同久久另类99精品国产91| 交换朋友夫妻互换小说| 久久久久网色| 黄频高清免费视频| 国产成+人综合+亚洲专区| 露出奶头的视频| 男男h啪啪无遮挡| 岛国毛片在线播放| 国产99久久九九免费精品| 99re6热这里在线精品视频| 国产一区二区 视频在线| 在线永久观看黄色视频| 成人18禁在线播放| 两个人免费观看高清视频| 国产成人啪精品午夜网站| 精品一品国产午夜福利视频| 精品欧美一区二区三区在线| 久久久国产欧美日韩av| 在线av久久热| 97在线人人人人妻| 免费不卡黄色视频| 精品亚洲乱码少妇综合久久| 精品一区二区三区四区五区乱码| 国产深夜福利视频在线观看| 久久久久久久久久久久大奶| 亚洲一区中文字幕在线| 天天操日日干夜夜撸| a级毛片黄视频| 69av精品久久久久久 | 久久午夜亚洲精品久久| 亚洲精品成人av观看孕妇| 亚洲av国产av综合av卡| 国产精品亚洲一级av第二区| 国产成人影院久久av| 女同久久另类99精品国产91| 亚洲人成电影观看| 动漫黄色视频在线观看| tocl精华| 黄色 视频免费看| 美女主播在线视频| 三上悠亚av全集在线观看| 中文字幕人妻熟女乱码| 纯流量卡能插随身wifi吗| 久久久国产成人免费| 久久久久久人人人人人| 女警被强在线播放| 国产精品影院久久| 久久久久久亚洲精品国产蜜桃av| 亚洲美女黄片视频| 国产单亲对白刺激| 国产又色又爽无遮挡免费看| 熟女少妇亚洲综合色aaa.| 国产精品香港三级国产av潘金莲| √禁漫天堂资源中文www| 另类亚洲欧美激情| 最近最新中文字幕大全电影3 | 黄频高清免费视频| 岛国在线观看网站| 人成视频在线观看免费观看| 老熟妇乱子伦视频在线观看| 久久亚洲真实| 一边摸一边做爽爽视频免费| 香蕉国产在线看| 久久久久久人人人人人| 三级毛片av免费| 精品人妻1区二区| 国产精品亚洲av一区麻豆| 搡老岳熟女国产| 美女扒开内裤让男人捅视频| 国产精品久久久人人做人人爽| 欧美精品av麻豆av| 日本精品一区二区三区蜜桃| 日本av手机在线免费观看| 国产在线视频一区二区| 亚洲色图综合在线观看| 国产免费现黄频在线看| 天天躁日日躁夜夜躁夜夜| 国产伦人伦偷精品视频| 国产精品国产高清国产av | 亚洲欧洲精品一区二区精品久久久| 欧美午夜高清在线| 亚洲欧美日韩高清在线视频 | 国产精品国产高清国产av | 两个人看的免费小视频| 啦啦啦 在线观看视频| 9191精品国产免费久久| 亚洲欧美一区二区三区黑人| 国产精品1区2区在线观看. | 一本—道久久a久久精品蜜桃钙片| 国产日韩欧美在线精品| 男女高潮啪啪啪动态图| 一本一本久久a久久精品综合妖精| 一级毛片女人18水好多| 亚洲人成77777在线视频| 无限看片的www在线观看| 欧美精品一区二区大全| 精品少妇久久久久久888优播| 高清视频免费观看一区二区| 大片免费播放器 马上看| 黄频高清免费视频| 精品少妇一区二区三区视频日本电影| 久9热在线精品视频| 亚洲专区中文字幕在线| 狂野欧美激情性xxxx| 久久亚洲精品不卡| 怎么达到女性高潮| 国产有黄有色有爽视频| 在线观看66精品国产| 男人操女人黄网站| 99国产极品粉嫩在线观看| 日韩中文字幕视频在线看片| 国产一卡二卡三卡精品| 日日爽夜夜爽网站| 久久国产精品大桥未久av| 欧美激情久久久久久爽电影 | 人人妻,人人澡人人爽秒播| 狂野欧美激情性xxxx| 国产精品久久久人人做人人爽| 色视频在线一区二区三区| 亚洲av电影在线进入| 老司机深夜福利视频在线观看| av超薄肉色丝袜交足视频| 国产精品久久久久久精品古装| 精品久久久久久电影网| 欧美黑人精品巨大| 国产精品电影一区二区三区 | 狂野欧美激情性xxxx| 国产免费av片在线观看野外av| 精品少妇久久久久久888优播| 黑丝袜美女国产一区| 9热在线视频观看99| 一区二区三区精品91| 纯流量卡能插随身wifi吗| 国产单亲对白刺激| 黄片小视频在线播放| 热re99久久国产66热| 国产福利在线免费观看视频| 不卡一级毛片| 一本—道久久a久久精品蜜桃钙片| av视频免费观看在线观看| 人人妻人人添人人爽欧美一区卜| 亚洲九九香蕉| 在线观看66精品国产| 99热国产这里只有精品6| 搡老岳熟女国产| 丰满迷人的少妇在线观看| 人人妻人人澡人人爽人人夜夜| 成年人午夜在线观看视频| 两性夫妻黄色片| 亚洲少妇的诱惑av| 正在播放国产对白刺激| 国产精品一区二区在线观看99| 美女国产高潮福利片在线看| 亚洲色图 男人天堂 中文字幕| videos熟女内射| 国产片内射在线| 日韩熟女老妇一区二区性免费视频| 热99国产精品久久久久久7| 桃红色精品国产亚洲av| 新久久久久国产一级毛片| 男女高潮啪啪啪动态图| 久久精品人人爽人人爽视色| 下体分泌物呈黄色| 欧美人与性动交α欧美精品济南到| 少妇精品久久久久久久| 久久久久久久久免费视频了| 中文字幕另类日韩欧美亚洲嫩草| 日韩中文字幕欧美一区二区| 亚洲国产av影院在线观看| 国产高清激情床上av| 日日爽夜夜爽网站| 交换朋友夫妻互换小说| 国产欧美日韩一区二区三区在线| 高清av免费在线| 国产欧美日韩精品亚洲av| 国产精品久久久久成人av| 91国产中文字幕| videos熟女内射| 国产aⅴ精品一区二区三区波| 一个人免费在线观看的高清视频| 在线观看66精品国产| 久久这里只有精品19| 国产av一区二区精品久久| 午夜福利视频精品| 久久亚洲真实| 黑人巨大精品欧美一区二区蜜桃| 欧美国产精品一级二级三级| 国产成人精品在线电影| 一二三四在线观看免费中文在| 人人妻人人澡人人看| 日韩欧美一区视频在线观看| 看免费av毛片| 天天操日日干夜夜撸| 午夜精品国产一区二区电影| 亚洲欧洲日产国产| 亚洲中文av在线| 午夜激情av网站| 99久久国产精品久久久| 他把我摸到了高潮在线观看 | 母亲3免费完整高清在线观看| 国产欧美日韩一区二区三区在线| 国产成人啪精品午夜网站| 国产精品.久久久| 国产精品免费视频内射| 性少妇av在线| 宅男免费午夜| 欧美人与性动交α欧美软件| 麻豆国产av国片精品| 国产不卡一卡二| 久久国产亚洲av麻豆专区| 亚洲精华国产精华精| 伊人久久大香线蕉亚洲五| www.熟女人妻精品国产| 国产精品av久久久久免费| 动漫黄色视频在线观看| 亚洲国产欧美在线一区| 大香蕉久久网| 侵犯人妻中文字幕一二三四区| 99久久精品国产亚洲精品| 亚洲视频免费观看视频| 黑人操中国人逼视频| 嫁个100分男人电影在线观看| 久久久国产欧美日韩av| 一级毛片女人18水好多| av网站免费在线观看视频| 亚洲av成人不卡在线观看播放网| 亚洲精品国产一区二区精华液| 视频区欧美日本亚洲| 麻豆乱淫一区二区| av又黄又爽大尺度在线免费看| 日韩三级视频一区二区三区| 国产精品1区2区在线观看. | 欧美大码av| 中文字幕另类日韩欧美亚洲嫩草| 黄色 视频免费看| 法律面前人人平等表现在哪些方面| 国产国语露脸激情在线看| 另类亚洲欧美激情| 亚洲专区字幕在线| 别揉我奶头~嗯~啊~动态视频| 久久精品熟女亚洲av麻豆精品| 亚洲中文日韩欧美视频| 桃红色精品国产亚洲av| 日韩视频在线欧美| 精品少妇内射三级| 国产精品成人在线| 精品视频人人做人人爽| 黑人操中国人逼视频| tocl精华| 国产av一区二区精品久久| 国产在视频线精品| 69精品国产乱码久久久| 他把我摸到了高潮在线观看 | 久久性视频一级片| 精品一区二区三卡| 99国产精品99久久久久| 欧美亚洲 丝袜 人妻 在线| 国产真人三级小视频在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲一区二区三区欧美精品| 国产精品麻豆人妻色哟哟久久| 黄色毛片三级朝国网站| 十八禁人妻一区二区| 日韩 欧美 亚洲 中文字幕| 亚洲午夜理论影院| 青草久久国产| 欧美黑人精品巨大| 日韩免费高清中文字幕av| 日韩有码中文字幕| 国产91精品成人一区二区三区 | 成年人黄色毛片网站| 日韩欧美一区视频在线观看| 日韩免费高清中文字幕av|