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

    基于嶺回歸的土壤全氮含量反演模型

    2022-11-26 06:57:52徐藝邵光成丁鳴鳴章二子何菁
    排灌機械工程學報 2022年11期
    關鍵詞:模型

    徐藝,邵光成*,丁鳴鳴,章二子,何菁

    (1. 河海大學農業(yè)工程與科學學院,江蘇 南京 210098; 2. 南京市水務局,江蘇 南京 210036; 3. 江寧淳化街道水務站,江蘇 南京 211122; 4. 南京市水利建筑工程檢測中心有限公司,江蘇 南京 210036)

    表層土壤全氮含量(total nitrogen,TN)為氮含量的總和,包括有機氮、硝態(tài)氮、亞硝態(tài)氮以及銨態(tài)氮,還包括部分聯氮、偶氮和疊氮等含氮化合物.傳統(tǒng)測定全氮含量的凱氏定氮法雖然精度高,然而步驟復雜、成本高,不適用于田間實際應用[1],而遙感監(jiān)測具有以下優(yōu)點:一方面,遙感不需要直接接觸觀測對象,從而避免對土壤表面作物的傷害.另一方面,可以開展長時間連續(xù)監(jiān)測,而且它的數據可以很容易地集成到地理信息系統(tǒng)中,以便后續(xù)進行分析[2].目前,常用的遙感技術包括衛(wèi)星遙感、無人機遙感、載人機遙感和地基遙感[3-4].相對于其他遙感平臺,無人機遙感具有成本低[5]、能獲取高空間分辨率影像、適合田塊尺度觀測等優(yōu)勢[6],在農情監(jiān)測中占據了重要位置[7].目前常采用的多光譜反演方法有偏最小二乘法、LASSO縮減系數法、主成分回歸等.其中,偏最小二乘回歸將關聯分析、主成分分析與多元線性回歸分析相結合[8],使用來自響應變量的信息以及預測因子進行特征轉換[9],曾用于熱帶森林的光譜和化學分析[10],并作為處理大型高光譜數據集的常用方法[11],但在采用最小二乘法求解回歸系數時發(fā)現光譜數據的自相關程度高,存在嚴重的多重共線性問題,反演結果不可靠[12].楊福芹等[13]基于多重共線性,篩選出對冬小麥氮營養(yǎng)指數相關性較高的圖像指數,再利用偏最小二乘法構建反演模型,建模集均方根誤差RMSE可達0.085 8,驗證集RMSE達0.187 1,預測精度較高;LAURIN等[14]采用偏最小二乘回歸(PLSR),利用高光譜數據與植被指數,對非洲熱帶雨林生物量進行建模,發(fā)現改進后模型精度(決定系數R2=0.70)優(yōu)于不考慮多重共線性的回歸模型(R2=0.64).

    高光譜數據通常由100多個帶寬為10 nm或更小的連續(xù)波段組成,波段之間相關性高,存在多重共線性和“維數災難”的問題,對反演模型的可靠性造成嚴重影響[3].為彌補高光譜數據中的多重共線性問題,改善預測模型精度,文中基于無人機多光譜高精度影像,獲取其光譜反射率數據,分析農田土壤表層全氮含量實測值與光譜反射率的多元線性回歸中的多重共線性問題,構建基于嶺回歸的無人機遙感影像反演土壤全氮含量預測模型,旨在探索一種兼顧反演精度與光譜數據多重共線性問題的方法,以便為無人機遙感土壤氮素營養(yǎng)診斷提供理論依據.

    1 材料與方法

    1.1 研究區(qū)概況

    試驗于2021年4月29日在位于南京市江寧區(qū)淳化街道某小型灌區(qū)插秧前農田進行.該農田位于119°4′0″~119°4′6″E,31°54′11″~34°54′17″N,農田總面積約為6 000 m2,無植被覆蓋,有利于觀測表層土壤光譜反射率.

    1.2 數據采集與處理

    1.2.1 數據采集

    試驗采用棋盤式布點法,選定10塊條田,每塊再選取6個取樣點,每個點位在表土層深度0~30 cm進行取樣,并使用紫外分光度計法測定土壤樣本的全氮含量,測點分布如圖1所示.

    圖1 采樣點分布圖Fig.1 Distribution of measurement sampling points

    農田遙感影像采用大疆P4 Mulitispectral無人機獲取.該型號無人機共搭載6個1/2.9 in CMOS影像傳感器,其中1個彩色傳感器用于常規(guī)可見光(RGB)成像,5個單色傳感器用于多光譜成像.單色傳感器前濾光片可通過波段:藍光波段(B),450±16 nm;綠光波段(G),560±16 nm;紅光波段(R),650±16 nm;紅邊波段(RE),730±16 nm;近紅外波段(NIR),840±26 nm.航拍在4月29日正午11:00—12:00進行,飛行高度50 m,飛行速度5 m/s,規(guī)劃航點299個,航向重疊率90%,旁向重疊率75%,主航線角度262°,主航線9條,云臺俯仰角-90°,共拍攝圖片1 794張,采用ArcGIS 10.7對圖片進行拼接.該航高下遙感影像的地面分辨率為2.65 cm/pixel,滿足精度要求.

    1.2.2 遙感影像拼接

    精靈4無人機懸停拍攝的圖片為帶有坐標的TIF格式,將其導入ArcGIS后,指定坐標系為WGS1984,利用ArcToolbox中的“鑲嵌”工具進行拼接,最終得到5張不同波段的農田影像.

    1.2.3 光譜反射率計算

    由無人機拍攝的原始影像數據點為像元值(digital number).為了得到對應波段的地表反射率,后續(xù)處理參照《P4 MultiSpectral 圖像處理指南》[15]進行.

    以求解藍光波段反射率RB為例(按文獻[15]格式)為

    (1)

    式中:Bluecamera為圖像信號值;pCamBlue為相機參數,查遙感圖像EXIF中XMP-drone-dji項可得1.355 955;BlueLS為藍光光強傳感器信號值,pLSBlue為藍光傳感器校準參數,兩者乘積可直接查XMP-drone-dji下的Irradiance得到,BlueLS×pLSBlue=12 629.98.

    Bluecamera可計算為

    (2)

    式中:IBlue和IBlackLevel分別為歸一化到值域[0,1]上的像素值和黑電平的值,其中IBlue計算方式參考式(3),IBlackLevel可以在遙感圖像信息中的EXIF-IFD0-BlackLevel得到,文中IBlackLevel=4 096/65 635=0.062 4;Bluegain為相機曝光時的增益參數;Blueetime為曝光時間.

    查圖像信息XMP-drone-dji與XMP-Camera項,可得Bluegain=1,Blueetime=0.364.

    (3)

    式中:DNBlue為藍光波段像元值.

    經過上述計算,可得到土壤樣本5個波段的光譜反射率R曲線,如圖2所示,圖中λ為波長.

    圖2 光譜反射率曲線Fig.2 Spectral reflectance curve

    2 結果與分析

    2.1 多元線性回歸建模

    假定以式(4)形式的回歸模型,反演表土全氮含量,即

    TNtopsoil=β1RB+β2RG+β3RR+

    β4RRE+β5RNIR+ε,

    (4)

    式中:TNtopsoil為表土全氮含量,g/kg;RB,RG,RR,RRE,RNIR分別為藍光波段、綠光波段、紅光波段、紅邊波段、近紅外波段的地表光譜反射率;β1—β5為回歸系數;ε為殘差.

    以五波段反射率為自變量,土壤全氮含量為因變量,50組有效數據作為樣本,利用SPSS 21.0軟件,采用最小二乘回歸求解系數β1—β5,結果為

    TNtopsoil=-8.437RB+18.310 4RG-28.724RR-
    18.569 3RRE+33.995RNIR.

    (5)

    該模型回歸系數R2達到了0.504,然而在顯著性t檢驗中,回歸系數均大于0.05,說明多元線性回歸系數不具有統(tǒng)計學意義.

    2.2 回歸系數不合理原因

    一般多元線性回歸模型可概化為

    y=Xβ+ε,

    (6)

    式中:X={x1,x2,…,xn}為解釋變量矩陣;y為被解釋變量矩陣;β為待估計的回歸系數;ε為殘差,滿足數學期望E(ε)=0.

    依據最小二乘法,回歸系數的參數估計值可表示為

    (7)

    多元線性回歸理論假定,參與回歸的各自變量間線性無關.若該假設不滿足,則會導致模型對誤差ε極敏感,回歸系數不可靠.然而,當自變量之間存在嚴重多重共線性問題時,XTX的行列式值接近于0,導致回歸系數估計值的解非常不穩(wěn)定.

    由實測反射率數據可得

    (8)

    使用最小二乘法計算回歸系數時,需要使用矩陣XTX的逆,而|XTX|=1.164 6×10-7接近于0,將導致計算出的回歸系數過大,且當解釋變量發(fā)生微小擾動時,回歸系數波動劇烈,甚至改變符號.為衡量自變量之間的多重共線性程度,MARQUARDT[16]于1970年提出方差膨脹因子(variance inflation factor, VIF),計算公式為

    (9)

    式中:VIFi為對應自變量Xi的方差膨脹系數,i=1,2,…,6.

    方差膨脹因子根據式(9)計算可得

    (10)

    由式(10)可知,與波段反射率RB,RG,RR,RRE,RNIR對應的方差膨脹因子分別為17.804 2,254.710 1,1 047.716 9,986.509 9,103.774 1.根據HOERL等[17]的研究,當VIF>10時,說明自變量之間存在多重共線性問題;若VIF>100,說明多重共線性現象嚴重.本次所測量的光譜數據中,有4個波段的VIF在100以上,說明光譜反射率數據存在較嚴重的多重共線性問題.

    2.3 嶺回歸建模

    嶺回歸是一種在自變量高度相關的情況下估計多元回歸模型系數的方法.該理論由HOERL等[18]首次提出.當線性回歸模型具有一些高度相關的獨立變量時,嶺回歸為最小二乘估計不精確的一種可靠的解決方案,它通過人為引入懲罰項kIp,回避了XTX行列式接近0的問題,即

    (11)

    k取值在0到1.當k=0時,嶺回歸退化為最小二乘估計.由于引入人為誤差項,該估計為有偏估計.HOERL等[18]證明,存在k值使得嶺回歸參數估計值的均方誤差小于最小二乘估計,并提出嶺跡法以確定合適的嶺回歸系數;嶺跡為所有標準化回歸系數與k的曲線圖,k∈[0,1].此外,HOERL等[18]提出了4個選擇最佳k值應滿足的條件:① 嶺跡線基本穩(wěn)定;② 回歸系數沒有不合理的數值大小;③ 回歸系數不再正負波動,符號變得合理;④ 殘差平方和相較多元線性回歸沒有顯著增加.

    圖3為嶺跡圖及各方差膨脹因子隨k的變化圖,圖中VIF1,VIF2,VIF3,VIF4,VIF5分別表示回歸系數β1—β5對應的方差膨脹因子.其中圖3a為嶺跡圖,βs1—βs5為標準化回歸系數,是將自變量矩陣與因變量進行z-score標準化[19]后回歸得到的系數,分別與回歸系數β1—β5對應.由嶺跡圖可知:① 當k=0時,回歸系數βs2,βs3,βs5的絕對值較穩(wěn)定時偏大,說明βs2,βs3,βs5被嚴重高估;隨著k值增大,βs1—βs5的絕對值都逐漸減小,且趨于穩(wěn)定.② 隨著k值增大,βs2由正轉負,最后達到穩(wěn)定.③ 在k≥0.025之后,標準化回歸系數沒有不合理的值.

    由圖3b—3f可知,5個方差膨脹因子均隨k增加而迅速減小,當k=0.04時,VIF均小于10,此時可認為多重共線性對預測模型的精度影響較小.

    圖4為R2、均方根誤差與P值隨k值的變化情況.嶺回歸與其他加入人為懲罰項的回歸方法類似,模型回歸擬合精度會因懲罰因子的增大而迅速降低,因此k不宜設置過大.圖4c中P1—P5分別對應回歸系數β1—β5的P值,用于評估不同嶺回歸系數下的顯著性水平[20].當P≤0.05時,可認為通過回歸系數的顯著性檢驗;當k≥0.12時,P1—P5均小于0.05,此時所有的回歸系數都通過顯著性檢驗.

    綜上所述,當k=0.12時,回歸決定系數R2從0.504降至0.408,且此時回歸系數趨于穩(wěn)定,數值合理,各波段方差膨脹因子均小于10,且P值小于0.05,表明所有回歸系數均通過顯著性檢驗.因此,取k=0.12作為最佳嶺回歸參數,此時表土全氮含量估算值可由五波段反射率表示為

    TNtopsoil=1.931 6RB-5.308 4RG-10.538 9RR-
    0.324 4RRE+13.746 1RNIR.

    (12)

    圖3 嶺跡圖和方差膨脹因子隨k值的變化圖Fig.3 Ridge trace and variation of variance inflation factors with k

    圖4 R2、均方根誤差與P值隨k值的變化Fig.4 Variation of R2, root mean square error and variation of P with k

    2.4 預測模型評價

    圖5為擬合結果,圖中TNp,TNt分別為土全氮的預測值、實測值.圖5a—5b是多元線性回歸與嶺回歸的反演結果.

    圖5 擬合結果Fig.5 Fitting results

    多元線性回歸反演全氮含量的R2為0.504,均方根誤差RMSE為0.472,均優(yōu)于k=0.12時嶺回歸預測值;然而由于多重共線性的影響,線性回歸在模型預測時表現較差.圖5c—5d是利用反演模型對驗證集數據進行預測的結果,其中驗證集數據由10個條田中各自隨機選取1個樣本點構成.可以看出,多元線性回歸預測模型的R2為0.645,RMSE為0.820,對部分點位全氮含量嚴重高估,建模效果雖好,但驗證效果差.而嶺回歸預測模型RMSE下降不明顯,有較好的預測效果.

    通過式(12)并結合農田光譜反射率數據圖,可以得到如圖6所示的完整表土全氮含量反演結果.利用表土全氮含量分布圖,可以快速獲取土壤全氮含量信息,對無人機遙感土壤氮素營養(yǎng)診斷和精準施肥具有重要意義.

    圖6 農田表土全氮含量反演結果Fig.6 Inversion results of total nitrogen content in farmland topsoil

    3 討 論

    反演模型中,近紅外波段反射率與表土全氮含量的相關度最高.這是因為近紅外波段光譜的信息來源于分子振動的倍頻和合頻,常用于含C-H,N-H,O-H等基團的有機物分析[21],而土壤中氮素絕大部分為有機結合態(tài),與有機質關聯密切[22],因此高分辨率的近紅外波段包含了土壤全氮含量的敏感譜區(qū).利用近紅外全譜波段對土壤全氮含量反演可以達到較高的精度.如李頡[23]利用12 500~3 600 cm-1的近紅外光譜數據建立的土壤全氮含量偏最小二乘回歸模型,決定系數可達到89.63%.利用多光譜數據進行全氮含量反演時,由于土壤全氮含量的取值區(qū)間較狹窄,方差較小,不利于穩(wěn)定預測模型的建立,但通過相關分析,篩選與全氮含量相關性強、顯著性高的敏感波段進行建模,仍能達到較高精度[24].

    文中研究嘗試了利用可見光波段(RGB)與紅邊(RE)、近紅外(NIR)波段光譜反射率數據,在不同嶺回歸系數下對土壤全氮含量進行預測,其中多元線性回歸模型是k=0時的特例.由于上述五波段數據獲取便利,雖然模型精度有所下降,但有助于大面積田塊全氮含量的快速、定性診斷.綜合反演精度與回歸系數顯著性兩方面考慮,選取k=0.12時的嶺回歸模型作為反演結果.

    4 結 論

    基于低空無人機搭載多光譜傳感器,通過獲取表土全氮含量和光譜反射率,研究分析了多元線性回歸在表土全氮含量光譜反演問題上的不足,揭示了光譜反射率數據特有的多重共線性問題及其應用限制.通過嶺回歸方法,以損失一定回歸精度為代價,得到一組對多重共線性不敏感的且穩(wěn)定的回歸系數.

    研究結果表明,通過表土全氮含量與光譜反射率多元線性回歸,發(fā)現波段RG,RR,RRE,RNIR對應的方差膨脹因子均大于100,自變量之間存在嚴重的多重共線性問題.而基于嶺回歸建模,當k=0.12時,回歸R2雖然有所降低,從0.504降至0.408,然而回歸系數趨于穩(wěn)定,數值合理,方差膨脹因子均小于10,且回歸系數之間差異具有統(tǒng)計學意義,說明基于光譜反射率信息反演土壤全氮含量的效果較好.但在實踐中獲取無人機影像進行應用時,還需要綜合考慮無人機不同飛行高度而引起的不同分辨率對反演結果的影響.

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務本地化模型
    適用于BDS-3 PPP的隨機模型
    提煉模型 突破難點
    函數模型及應用
    p150Glued在帕金森病模型中的表達及分布
    函數模型及應用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權M-估計的漸近分布
    3D打印中的模型分割與打包
    天堂√8在线中文| 波多野结衣高清无吗| 国产视频内射| 国产爱豆传媒在线观看| 白带黄色成豆腐渣| 免费看光身美女| 国产精品一区二区三区四区免费观看 | 亚洲人成网站在线播| 国产单亲对白刺激| 一级av片app| 女的被弄到高潮叫床怎么办 | 精品国内亚洲2022精品成人| 干丝袜人妻中文字幕| 桃色一区二区三区在线观看| 国产日本99.免费观看| 最新在线观看一区二区三区| 女人被狂操c到高潮| 一级av片app| 久久精品国产亚洲av天美| 国产精品久久电影中文字幕| 久久久久久久久久久丰满 | 免费观看在线日韩| 久久久久久九九精品二区国产| 国产蜜桃级精品一区二区三区| av视频在线观看入口| 国产中年淑女户外野战色| 日日摸夜夜添夜夜添小说| av.在线天堂| 高清日韩中文字幕在线| 男插女下体视频免费在线播放| 日本三级黄在线观看| 看免费成人av毛片| 午夜a级毛片| 嫩草影院新地址| 最新在线观看一区二区三区| 日韩,欧美,国产一区二区三区 | 直男gayav资源| 色综合色国产| 色综合亚洲欧美另类图片| 岛国在线免费视频观看| 免费黄网站久久成人精品| 成人高潮视频无遮挡免费网站| 精品久久久久久,| 我要看日韩黄色一级片| 一级av片app| 国产免费男女视频| 婷婷色综合大香蕉| 丝袜美腿在线中文| 黄色丝袜av网址大全| h日本视频在线播放| 男人舔女人下体高潮全视频| 欧美激情在线99| 美女高潮的动态| 在线播放无遮挡| avwww免费| 国产男人的电影天堂91| 国内精品宾馆在线| 亚洲在线自拍视频| 精品一区二区三区人妻视频| 少妇的逼好多水| 亚洲国产色片| 免费一级毛片在线播放高清视频| 国产三级中文精品| 深爱激情五月婷婷| 亚洲av不卡在线观看| 欧美日韩瑟瑟在线播放| 久久亚洲精品不卡| 久久久久国内视频| 日韩大尺度精品在线看网址| 97碰自拍视频| 一级毛片久久久久久久久女| 午夜福利成人在线免费观看| 波野结衣二区三区在线| 深夜精品福利| 18禁黄网站禁片午夜丰满| 国产av麻豆久久久久久久| 国产乱人视频| 偷拍熟女少妇极品色| 干丝袜人妻中文字幕| 欧美绝顶高潮抽搐喷水| 婷婷六月久久综合丁香| 99久久精品国产国产毛片| 亚洲三级黄色毛片| 神马国产精品三级电影在线观看| 观看美女的网站| 如何舔出高潮| 亚洲精品一卡2卡三卡4卡5卡| 99久久精品热视频| 国产精华一区二区三区| 国产伦精品一区二区三区四那| 日韩人妻高清精品专区| 国内毛片毛片毛片毛片毛片| 99热这里只有精品一区| 日本-黄色视频高清免费观看| 波多野结衣高清作品| 熟女人妻精品中文字幕| 欧美人与善性xxx| 18禁黄网站禁片免费观看直播| 女人十人毛片免费观看3o分钟| 一卡2卡三卡四卡精品乱码亚洲| 色综合色国产| 亚洲18禁久久av| 日韩亚洲欧美综合| 看片在线看免费视频| 99国产极品粉嫩在线观看| 一级毛片久久久久久久久女| 黄色丝袜av网址大全| 永久网站在线| 嫩草影院精品99| 性欧美人与动物交配| 黄色一级大片看看| 精品午夜福利在线看| 成人性生交大片免费视频hd| 国产伦在线观看视频一区| 18禁黄网站禁片免费观看直播| 亚洲五月天丁香| 中文在线观看免费www的网站| 日本黄色视频三级网站网址| 综合色av麻豆| 成人二区视频| 国产精品亚洲一级av第二区| 琪琪午夜伦伦电影理论片6080| 一进一出抽搐动态| 老司机福利观看| 日日夜夜操网爽| 日日摸夜夜添夜夜添av毛片 | 俄罗斯特黄特色一大片| 美女cb高潮喷水在线观看| 亚洲一区高清亚洲精品| 中文字幕久久专区| 变态另类成人亚洲欧美熟女| 亚洲国产精品sss在线观看| 精品久久国产蜜桃| 给我免费播放毛片高清在线观看| 免费黄网站久久成人精品| 成人永久免费在线观看视频| 色哟哟哟哟哟哟| 亚洲一级一片aⅴ在线观看| 少妇猛男粗大的猛烈进出视频 | 国产乱人伦免费视频| 久久久久久伊人网av| 麻豆成人午夜福利视频| 天天一区二区日本电影三级| 国产成年人精品一区二区| 人妻少妇偷人精品九色| 老熟妇乱子伦视频在线观看| 国产精品,欧美在线| 国产淫片久久久久久久久| 国产精品不卡视频一区二区| 久久99热这里只有精品18| 久久国内精品自在自线图片| 国产精品亚洲美女久久久| 国产伦精品一区二区三区视频9| 给我免费播放毛片高清在线观看| 欧美+日韩+精品| 亚洲美女视频黄频| 国产成人aa在线观看| 国产精品亚洲美女久久久| 久久久成人免费电影| 99在线视频只有这里精品首页| 国产精品人妻久久久久久| 久久精品久久久久久噜噜老黄 | 亚洲精品亚洲一区二区| 99久久无色码亚洲精品果冻| 一本精品99久久精品77| 日韩欧美国产一区二区入口| 欧美成人一区二区免费高清观看| 在线国产一区二区在线| 国内揄拍国产精品人妻在线| av女优亚洲男人天堂| 91狼人影院| 一区福利在线观看| 国产淫片久久久久久久久| 久久久久久久久久成人| 国产成年人精品一区二区| 在线免费观看不下载黄p国产 | 一边摸一边抽搐一进一小说| 色5月婷婷丁香| 热99在线观看视频| 可以在线观看的亚洲视频| 久久精品国产亚洲av涩爱 | 老熟妇仑乱视频hdxx| 久久久精品欧美日韩精品| 欧美3d第一页| 麻豆久久精品国产亚洲av| 99九九线精品视频在线观看视频| 九色国产91popny在线| 中文字幕av成人在线电影| 99在线人妻在线中文字幕| aaaaa片日本免费| 亚洲av二区三区四区| 久久中文看片网| 午夜精品久久久久久毛片777| 免费看光身美女| 五月玫瑰六月丁香| 啦啦啦啦在线视频资源| 成年女人看的毛片在线观看| 99九九线精品视频在线观看视频| 精品人妻一区二区三区麻豆 | 99久久九九国产精品国产免费| 美女高潮喷水抽搐中文字幕| 99热这里只有是精品在线观看| 色综合色国产| 国产精品爽爽va在线观看网站| 中文字幕免费在线视频6| 女同久久另类99精品国产91| 深爱激情五月婷婷| 可以在线观看毛片的网站| 精华霜和精华液先用哪个| 给我免费播放毛片高清在线观看| 国产精品一区二区三区四区免费观看 | 国产av不卡久久| 最新在线观看一区二区三区| 又紧又爽又黄一区二区| 男女之事视频高清在线观看| 精品久久久久久久人妻蜜臀av| av国产免费在线观看| 国产欧美日韩一区二区精品| 免费无遮挡裸体视频| 女生性感内裤真人,穿戴方法视频| 亚洲av电影不卡..在线观看| 2021天堂中文幕一二区在线观| 春色校园在线视频观看| 欧美+日韩+精品| 日本精品一区二区三区蜜桃| 免费黄网站久久成人精品| 成人国产一区最新在线观看| 春色校园在线视频观看| 国产精品国产高清国产av| 精品久久国产蜜桃| 日本一本二区三区精品| 国产主播在线观看一区二区| 在线看三级毛片| 亚洲人成网站在线播放欧美日韩| 岛国在线免费视频观看| av在线亚洲专区| 欧美黑人欧美精品刺激| 亚洲国产精品久久男人天堂| 久久久午夜欧美精品| a在线观看视频网站| 国产亚洲精品av在线| 99久久中文字幕三级久久日本| 亚洲一区高清亚洲精品| 国产三级中文精品| 麻豆国产av国片精品| 人妻制服诱惑在线中文字幕| 床上黄色一级片| 直男gayav资源| 亚洲美女黄片视频| bbb黄色大片| 中文亚洲av片在线观看爽| 免费看美女性在线毛片视频| 窝窝影院91人妻| 久久婷婷人人爽人人干人人爱| 国产中年淑女户外野战色| 日本黄色视频三级网站网址| 天天躁日日操中文字幕| 不卡视频在线观看欧美| or卡值多少钱| 日韩大尺度精品在线看网址| 色吧在线观看| 日本在线视频免费播放| av在线观看视频网站免费| 一本一本综合久久| .国产精品久久| 18禁在线播放成人免费| 成人精品一区二区免费| 亚洲精品粉嫩美女一区| 午夜精品一区二区三区免费看| 欧美日韩精品成人综合77777| 精品人妻熟女av久视频| 亚洲一区二区三区色噜噜| 欧美成人a在线观看| 亚洲国产日韩欧美精品在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲三级黄色毛片| 亚洲自偷自拍三级| 亚洲图色成人| 在线观看一区二区三区| 欧美丝袜亚洲另类 | 国产亚洲欧美98| 一区二区三区高清视频在线| 嫁个100分男人电影在线观看| 看免费成人av毛片| 亚洲av第一区精品v没综合| a级一级毛片免费在线观看| 国产91精品成人一区二区三区| 男人狂女人下面高潮的视频| 久久人人精品亚洲av| 亚洲四区av| 国产黄片美女视频| 欧美在线一区亚洲| 黄色一级大片看看| 最近中文字幕高清免费大全6 | 国产男人的电影天堂91| 99热这里只有是精品50| 99久久中文字幕三级久久日本| 中文字幕免费在线视频6| 亚洲乱码一区二区免费版| 亚洲成人久久爱视频| 国产午夜福利久久久久久| 色精品久久人妻99蜜桃| 波野结衣二区三区在线| 国产精品久久久久久av不卡| 少妇熟女aⅴ在线视频| 久久久精品欧美日韩精品| 日本一本二区三区精品| 99久久中文字幕三级久久日本| 人妻久久中文字幕网| 国产三级中文精品| 一级黄片播放器| 动漫黄色视频在线观看| 一边摸一边抽搐一进一小说| 波野结衣二区三区在线| 啪啪无遮挡十八禁网站| 少妇高潮的动态图| 女人十人毛片免费观看3o分钟| 国产成人一区二区在线| 欧美性猛交╳xxx乱大交人| 免费人成在线观看视频色| 99热只有精品国产| 99热这里只有是精品50| 亚洲国产欧洲综合997久久,| 人妻丰满熟妇av一区二区三区| 99热网站在线观看| 中文字幕免费在线视频6| 亚洲欧美日韩无卡精品| 极品教师在线视频| www日本黄色视频网| 99热这里只有是精品在线观看| 成人国产一区最新在线观看| 一级黄片播放器| 日韩一本色道免费dvd| 在现免费观看毛片| av在线蜜桃| 久久久久九九精品影院| 中文亚洲av片在线观看爽| 亚洲五月天丁香| 国产精品国产高清国产av| 精品免费久久久久久久清纯| 亚洲欧美日韩东京热| 欧美一区二区国产精品久久精品| 天堂av国产一区二区熟女人妻| 美女高潮喷水抽搐中文字幕| 国产在线精品亚洲第一网站| 欧美精品国产亚洲| 午夜精品一区二区三区免费看| 91麻豆精品激情在线观看国产| 成熟少妇高潮喷水视频| 热99在线观看视频| 天堂影院成人在线观看| 亚州av有码| 真实男女啪啪啪动态图| 十八禁网站免费在线| 欧美一级a爱片免费观看看| 美女黄网站色视频| 又黄又爽又刺激的免费视频.| 久久久久久久久大av| 尤物成人国产欧美一区二区三区| 97超视频在线观看视频| 成人美女网站在线观看视频| 天堂√8在线中文| 69人妻影院| 欧美高清性xxxxhd video| 欧美日本亚洲视频在线播放| 又爽又黄无遮挡网站| 午夜福利在线观看免费完整高清在 | 中文字幕av成人在线电影| 最近在线观看免费完整版| av.在线天堂| 99久久精品国产国产毛片| 欧美三级亚洲精品| 99久久精品国产国产毛片| 在线观看免费视频日本深夜| 日本欧美国产在线视频| 1000部很黄的大片| 午夜激情福利司机影院| 色在线成人网| 国产精品伦人一区二区| 精品国产三级普通话版| 午夜福利成人在线免费观看| 久久精品久久久久久噜噜老黄 | 久久精品影院6| 国产探花在线观看一区二区| 91久久精品国产一区二区三区| 成人一区二区视频在线观看| 亚洲内射少妇av| 看片在线看免费视频| 五月玫瑰六月丁香| 成人综合一区亚洲| 欧美性猛交黑人性爽| 99精品在免费线老司机午夜| 国产精品爽爽va在线观看网站| 人妻丰满熟妇av一区二区三区| 国产av不卡久久| 国产精品野战在线观看| 亚洲狠狠婷婷综合久久图片| 国产伦在线观看视频一区| 一夜夜www| 又黄又爽又免费观看的视频| 一边摸一边抽搐一进一小说| 国产熟女欧美一区二区| 日本欧美国产在线视频| 欧美+亚洲+日韩+国产| 亚洲欧美激情综合另类| 五月伊人婷婷丁香| 级片在线观看| 亚洲国产欧洲综合997久久,| 午夜免费男女啪啪视频观看 | aaaaa片日本免费| 午夜福利18| 91麻豆av在线| 久久婷婷人人爽人人干人人爱| 亚洲天堂国产精品一区在线| 欧美精品国产亚洲| 中文字幕高清在线视频| 日韩欧美免费精品| 伦理电影大哥的女人| 欧美成人a在线观看| 在线天堂最新版资源| 国内精品宾馆在线| 日日啪夜夜撸| 欧美成人免费av一区二区三区| 国产主播在线观看一区二区| 精品一区二区三区av网在线观看| 久久精品综合一区二区三区| 内射极品少妇av片p| 亚洲国产精品成人综合色| 国产高清不卡午夜福利| 亚洲av免费高清在线观看| 久久久久久九九精品二区国产| 国产成人av教育| 亚洲成人精品中文字幕电影| 欧美国产日韩亚洲一区| 亚洲国产欧洲综合997久久,| 国产精品不卡视频一区二区| 丝袜美腿在线中文| 亚洲中文字幕一区二区三区有码在线看| 午夜爱爱视频在线播放| 亚洲欧美日韩无卡精品| 国产精品久久电影中文字幕| videossex国产| 国产白丝娇喘喷水9色精品| 亚洲精品久久国产高清桃花| 欧美最新免费一区二区三区| 少妇裸体淫交视频免费看高清| 蜜桃久久精品国产亚洲av| 91午夜精品亚洲一区二区三区 | 欧美+日韩+精品| 日日撸夜夜添| 老司机福利观看| 日韩精品青青久久久久久| 国产精品嫩草影院av在线观看 | 九九久久精品国产亚洲av麻豆| 老司机午夜福利在线观看视频| 久久久久久九九精品二区国产| 性色avwww在线观看| 白带黄色成豆腐渣| 免费观看人在逋| 成人精品一区二区免费| 波多野结衣高清作品| 国产 一区 欧美 日韩| 久久精品久久久久久噜噜老黄 | 国产精品99久久久久久久久| 综合色av麻豆| 天天一区二区日本电影三级| 久久久久九九精品影院| 男女啪啪激烈高潮av片| 成人亚洲精品av一区二区| 高清日韩中文字幕在线| 精品久久久久久久人妻蜜臀av| 一级毛片久久久久久久久女| 日本成人三级电影网站| 国产真实乱freesex| 亚洲精品一区av在线观看| 成人亚洲精品av一区二区| 国产91精品成人一区二区三区| 国产欧美日韩精品一区二区| 欧美日韩乱码在线| 全区人妻精品视频| 人妻少妇偷人精品九色| 亚洲精品成人久久久久久| av黄色大香蕉| 露出奶头的视频| 成人高潮视频无遮挡免费网站| 亚洲精品久久国产高清桃花| 精品人妻偷拍中文字幕| 欧美日韩亚洲国产一区二区在线观看| 国产精品,欧美在线| 午夜免费男女啪啪视频观看 | 亚洲男人的天堂狠狠| 国产私拍福利视频在线观看| 午夜视频国产福利| 国产黄色小视频在线观看| 亚洲专区国产一区二区| 久久久国产成人免费| 精品人妻视频免费看| 黄色丝袜av网址大全| 亚洲内射少妇av| 日本熟妇午夜| 久久久久国内视频| 久久久午夜欧美精品| 精品午夜福利视频在线观看一区| 在线免费观看的www视频| 久久久久久久久中文| 久久6这里有精品| 亚洲一区二区三区色噜噜| 欧美xxxx性猛交bbbb| 色噜噜av男人的天堂激情| 国产在视频线在精品| 国产一区二区亚洲精品在线观看| 欧美绝顶高潮抽搐喷水| 午夜福利视频1000在线观看| 婷婷精品国产亚洲av在线| 精品午夜福利视频在线观看一区| 欧美精品国产亚洲| 国产午夜精品论理片| 一本久久中文字幕| 男人舔奶头视频| 99热这里只有精品一区| 国产精品一区二区三区四区久久| 内射极品少妇av片p| 久久午夜福利片| 国产爱豆传媒在线观看| 亚洲av电影不卡..在线观看| 波多野结衣高清作品| 久久精品国产亚洲av涩爱 | 搡老妇女老女人老熟妇| 日本成人三级电影网站| 国产视频内射| 日本 av在线| 免费人成在线观看视频色| 深夜精品福利| 免费人成在线观看视频色| 波多野结衣高清作品| 久久精品国产鲁丝片午夜精品 | 精品日产1卡2卡| 最好的美女福利视频网| 午夜精品久久久久久毛片777| 在线观看66精品国产| 高清毛片免费观看视频网站| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲无线在线观看| 国产精品一区二区免费欧美| 久久久久久久久久黄片| 九色国产91popny在线| .国产精品久久| 亚洲国产精品成人综合色| 亚洲国产日韩欧美精品在线观看| 床上黄色一级片| 99精品在免费线老司机午夜| 亚洲熟妇中文字幕五十中出| 又黄又爽又刺激的免费视频.| 99久久中文字幕三级久久日本| 久久久国产成人精品二区| 欧美又色又爽又黄视频| 国产精华一区二区三区| 国产高潮美女av| 一级黄片播放器| 日本黄色视频三级网站网址| av女优亚洲男人天堂| 欧美日本亚洲视频在线播放| 老司机午夜福利在线观看视频| 一区二区三区高清视频在线| 精品欧美国产一区二区三| 身体一侧抽搐| 美女 人体艺术 gogo| 日本黄色视频三级网站网址| 亚洲国产色片| 大型黄色视频在线免费观看| 国产成人av教育| 99国产极品粉嫩在线观看| 免费看日本二区| 午夜a级毛片| 又爽又黄a免费视频| 亚洲熟妇熟女久久| 午夜日韩欧美国产| 欧美日韩乱码在线| 亚洲经典国产精华液单| 久久久久九九精品影院| 色综合婷婷激情| 偷拍熟女少妇极品色| 伦理电影大哥的女人| 天天躁日日操中文字幕| 简卡轻食公司| 成人国产一区最新在线观看| 国产69精品久久久久777片| 国产真实伦视频高清在线观看 | 亚洲成人中文字幕在线播放| 直男gayav资源| 亚洲自偷自拍三级| 人人妻人人看人人澡| 欧美日韩乱码在线| 黄色日韩在线| 成人性生交大片免费视频hd| 国产精品无大码| 国内久久婷婷六月综合欲色啪| 国产一区二区在线观看日韩| 亚洲欧美激情综合另类| 午夜福利18| 亚洲人成网站在线播| 91麻豆精品激情在线观看国产| 国产高清三级在线| 国产欧美日韩精品亚洲av| 国产成人一区二区在线| 日本 欧美在线| 色哟哟哟哟哟哟| 精品无人区乱码1区二区| eeuss影院久久| 亚洲av成人精品一区久久| 哪里可以看免费的av片| 赤兔流量卡办理| 搡老妇女老女人老熟妇| 亚洲经典国产精华液单| 禁无遮挡网站|