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

    艾比湖流域鹽漬土含水量光譜特征分析與建模

    2017-04-18 07:09:02張海威張飛李哲井云清
    中國(guó)水土保持科學(xué) 2017年1期
    關(guān)鍵詞:鹽漬方根微分

    張海威,張飛,3?,李哲,井云清

    (1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院,830046,烏魯木齊;2.新疆大學(xué)綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,830046,烏魯木齊; 3.新疆智慧城市與環(huán)境建模普通高校重點(diǎn)實(shí)驗(yàn)室,830046,烏魯木齊)

    艾比湖流域鹽漬土含水量光譜特征分析與建模

    張海威1,2,張飛1,2,3?,李哲1,2,井云清1,2

    (1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院,830046,烏魯木齊;2.新疆大學(xué)綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,830046,烏魯木齊; 3.新疆智慧城市與環(huán)境建模普通高校重點(diǎn)實(shí)驗(yàn)室,830046,烏魯木齊)

    鹽漬土表層鹽分累積與土壤含水量有著直接的關(guān)系。為了建立干旱區(qū)鹽漬土含水量高光譜遙感監(jiān)測(cè)模型,本文以艾比湖流域不同含水量的鹽漬土為研究對(duì)象,采用光譜反射率變換和多元統(tǒng)計(jì)分析(MSAM)方法,對(duì)土壤含水量的光譜特征進(jìn)行分析和建模。結(jié)果表明:隨著土壤含水量的增加,土壤反射率呈下降趨勢(shì);在一定范圍內(nèi),波長(zhǎng)越長(zhǎng),土壤光譜反射率與含水量的相關(guān)性越高,其中1 937 nm處的土壤光譜反射率與含水量具有最高的相關(guān)性(r=-0.636)。對(duì)土壤光譜反射率進(jìn)行8種光譜反射率變換后,在此基礎(chǔ)上利用多元統(tǒng)計(jì)方法分析鹽漬土的不同含水量與光譜之間的相關(guān)性,篩選敏感波段,建立關(guān)系模型。得出對(duì)數(shù)一階微分(Logarithm First Order Differential)在波長(zhǎng)2 024和2 357 nm建立的模型以及均方根一階微分(Root Mean Square First Order Differential)在波長(zhǎng)1 972和2 357 nm建立的模型最優(yōu),相關(guān)系數(shù)r分別為0.894和0.865。基于上述模型作者構(gòu)建了一種耦合模型,其相關(guān)性r=0.926比對(duì)數(shù)一階微分模型提高了0.032,比均方根一階模型提高了0.061;因此,所構(gòu)建的鹽漬土含水量估算模型是可行的,可以為遙感反演提供理論參考,對(duì)高光譜遙感反演具有一定意義。

    光譜特征;含水量;模型;艾比湖流域

    干旱區(qū)鹽漬土研究一直是土壤學(xué)科的熱門問題之一。土壤含水量不僅是微觀氣象學(xué)中一個(gè)重要的能量平衡參數(shù),而且還是旱情監(jiān)測(cè)的重要參數(shù)之一[1],因此土壤含水量對(duì)土壤鹽漬化的研究具有特殊意義。土壤含水量不僅與次生鹽漬化有關(guān),還與鹽漬土的鹽分運(yùn)移有著密切的關(guān)系,是影響綠洲生態(tài)環(huán)境穩(wěn)定性的重要因素之一[2]。國(guó)內(nèi)外一些學(xué)者建立光譜與土壤含水量模型,通常利用微波遙感對(duì)地表層濕度直接進(jìn)行數(shù)據(jù)收集,或利用熱量距平植被指數(shù)、作物缺水指數(shù)等方法,獲取地表能量與作物生長(zhǎng)狀況建立模型[35]。當(dāng)前分析土壤理化性質(zhì)與光譜反射率之間關(guān)系,主要使用主成分分析、逐步多元線性回歸分析、偏最小二乘法回歸分析及支持向量機(jī)和人工神經(jīng)網(wǎng)絡(luò)方法[67]。Lu Ning等[8]采用PHI高光譜圖像數(shù)據(jù),利用多元回歸分析與最小二乘法等方法對(duì)土壤鹽漬化信息進(jìn)行了可視化。孫小艷等[9]使用手持便攜式光譜儀對(duì)小麥的光譜反射率進(jìn)行研究。夏學(xué)齊等[10]對(duì)不同理化性質(zhì)的土壤進(jìn)行了光譜分析與建模。吳亞坤等[11]通過變換光譜反射率與土壤EC之間構(gòu)建回歸模型,分析變化光譜反射率與EC關(guān)系。以上研究均是在野外進(jìn)行的實(shí)驗(yàn)研究,另外也有很多學(xué)者開展了室內(nèi)光譜研究,如劉煥軍等[12]對(duì)黑土土壤水分進(jìn)行了光譜響應(yīng)特征分析與建模,何挺等[13]對(duì)土壤水分光譜特征進(jìn)行對(duì)比分析,翁永玲等[14]基于PLSR方法對(duì)青海茶卡-共和盆地的土壤鹽分進(jìn)行高光譜遙感反演。他們都是在室內(nèi)進(jìn)行土壤含水量研究,均取得一定的成果。目前,利用實(shí)測(cè)光譜與土壤理化性質(zhì)建立估算模型,已成為鹽漬土研究的主流方法。

    筆者以室內(nèi)實(shí)測(cè)光譜數(shù)據(jù)和鹽漬土含水量為數(shù)據(jù)源,提出基于不同模型下的耦合建模思想,彌補(bǔ)單一性建模方法的不足。由于針對(duì)新疆典型生態(tài)脆弱區(qū)土壤含水量的光譜特性及建模研究還比較缺乏;因此,選擇新疆艾比湖流域?yàn)槭痉秴^(qū),以鹽漬土為研究對(duì)象,選用多元統(tǒng)計(jì)分析等方法,建立鹽漬土不同含水量估算模型,為治理與監(jiān)測(cè)鹽漬土提供技術(shù)支持。

    1 研究區(qū)概況

    艾比湖流域地處歐亞大陸腹地,位于E 79°53'~85°02',N 43°38'~45°52',三面環(huán)山,遠(yuǎn)離海洋,位于托里縣和烏蘇南部。研究區(qū)屬于典型的中溫帶干旱大陸氣候,以干旱少雨、氣溫變化劇烈為特點(diǎn)。年平均氣溫6.6~7.8℃,年降水量為116.0~169.2 mm。近40年來,由于入湖水量的減少,艾比湖湖面常年萎縮,湖面曾縮至不足500 km2,水位下降近23 m,干涸湖地也早已淪為鹽漠,成為浮塵天氣發(fā)源地。植被也依次形成旱生、超旱生、沙生、鹽生、水生等多種類型,主要植物類型有梭梭(Haloxylon ammodendron)、胡楊(Populus euphratica)、檉柳(Tamar-rix chinensis)和蘆葦(Phragmites australis)等[1516]。

    2 材料與方法

    2.1 數(shù)據(jù)的采集與處理

    2015年10月,對(duì)艾比湖流域進(jìn)行綜合調(diào)查,整個(gè)流域共布設(shè)有53個(gè)樣點(diǎn),用GPS對(duì)其進(jìn)行坐標(biāo)定位同時(shí)采集距地面0~10 cm的土壤樣品,采集土壤時(shí)用鋁盒裝土樣,這樣可以防止水分的散失。對(duì)采集回來的土樣,立即使用便攜式地物光譜(Field Spec Hi res,ASD,美國(guó))的人工光源采集每個(gè)土樣光譜曲線,該儀器光譜范圍為350~2 500 nm,光譜分辨率3 nm。容器半徑為2.5 cm,垂直探頭,探頭距離為5 cm。測(cè)量步驟為:首先開機(jī)預(yù)熱2~3 min、去除暗電流、利用白板定標(biāo)獲取絕對(duì)反射率、容器土樣表面刮平、收集光譜曲線5條;每采集5個(gè)樣點(diǎn)光譜要利用白板定標(biāo);后利用光譜儀的后處理軟件ASD View Spec Pro對(duì)采集的光譜曲線分析,對(duì)每個(gè)點(diǎn)號(hào)的5條光譜數(shù)據(jù)在Process下的Statistics內(nèi)的Mean命令進(jìn)行平均值計(jì)算,最后光譜數(shù)據(jù)的歸一化;最后在Origin 8.0軟件中進(jìn)行微分變換與多元統(tǒng)計(jì)分析。

    將土樣放入烘箱,在105℃下烘干12 h,放于干燥器中冷卻,然后進(jìn)行干土樣稱量,根據(jù)土壤含水量測(cè)定方法[17]計(jì)算出土壤含水量

    式中:msc為土壤含水量,g/kg,m為鋁盒質(zhì)量,g;m1為土壤烘干前質(zhì)量,g;m2為土壤烘干后質(zhì)量,g。

    2.2 光譜數(shù)據(jù)分析

    為更好地分析與展示光譜數(shù)據(jù)與土壤含水量之間的相關(guān)性,對(duì)土壤光譜反射率分別進(jìn)行倒數(shù)、對(duì)數(shù)和均方根等數(shù)學(xué)變換,然后再進(jìn)行一階、二階微分變換消除干擾。一階、二階微分方程[18]分別為:

    式中:λi+1、λi-1和λi為波長(zhǎng);R'λi為波長(zhǎng)λi的一階微分光譜;R″λi為波長(zhǎng)λi的二階微分光譜。將變換后的反射率一階微分、反射率二階微分、倒數(shù)一階微分、對(duì)數(shù)一階微分、均方根一階微分與原始反射率作為鹽漬土的高光譜指數(shù)。

    2.3 相關(guān)性分析

    通過每個(gè)波段計(jì)算鹽漬土的高光譜指數(shù)與土壤含水量的相關(guān)系數(shù),篩選出鹽漬土含水量的敏感波段[19]即

    式中:rj為高光譜指數(shù)與土壤含水量的相關(guān)系數(shù);j為波段;Rij為第i個(gè)樣品第j波段的反射率值;個(gè)土樣在第j波_段反射率的平均值;Wi為第i個(gè)土樣的水分含量,g;為土壤樣本含水量的平均值,g;n為鹽漬土樣品個(gè)數(shù)。

    3 結(jié)果與分析

    3.1 鹽漬土含水量光譜特征分析

    通過室內(nèi)分析發(fā)現(xiàn)艾比湖流域鹽漬土含水量的范圍為2.6%~25.6%。把含水量的范圍分為5個(gè)層次:0~5%、5~10%、10~15%、15~20%、20~25%,分別從各個(gè)層次隨機(jī)選取一條來分析不同含水量的光譜特征。不同含水量下鹽漬土光譜特征曲線趨勢(shì)和態(tài)勢(shì)很相似,但反射率高低不同。如圖1所示,總體來看,光譜曲線從400到1350 nm呈平緩上升趨勢(shì),但不同含水量?jī)A斜度不同。其中含水量為10.2%的光譜曲線在380~600 nm范圍內(nèi),土壤光譜反射率最高。含水量為9.7%的光譜曲線在380~1 350 nm范圍內(nèi),不符合隨著含水量的增加而反射率降低的規(guī)律,這可能是由于土壤中的含鹽量過高所導(dǎo)致的此現(xiàn)象,在600~1450 nm附近反射率曲線斜率變緩,在1 470和1 950 nm左右各有一個(gè)水分反射谷,其中在1 950 nm附近的反射谷比較明顯。不同土壤含水量反射谷程度不同,反射率呈先升高后降低的趨勢(shì)。22.1%、10.2%和9.7%3條不同含水量的光譜曲線,在1 950 nm附近出現(xiàn)很明顯的反射谷,且這3條曲線反射率差異很小,而4.6%含水量的曲線差異較高。從整體看來隨著含水量的增加反射率降低,不同含水量增減速率不同,這與顧燕等[20]和劉曾媛等[21]的研究一致。

    圖1 不同含水量鹽漬土光譜反射率曲線Fig.1 Spectral reflectance curves of salinized soil with different levels of water content

    3.2 敏感波段的選取與分析

    在Origin 8.0軟件中進(jìn)行8種數(shù)學(xué)變換:一階微分、二階微分、倒數(shù)、倒數(shù)一階微分、對(duì)數(shù)、對(duì)數(shù)一階微分、均方根、均方根一階微分,之后分別與土壤含水量建立相關(guān)關(guān)系,可以看出在不同波段上變化差異比較大(圖2)。

    圖2 土壤含水量和光譜指數(shù)的相關(guān)系數(shù)(R為反射率)Fig.2 Correlation coefficients between soil water content and spectral indices(R represents reflectance)

    由圖2可得,對(duì)光譜反射率進(jìn)行微分變換后,其與鹽漬土含水量的相關(guān)性要明顯高于原始反射率,既有明顯正相關(guān),又有負(fù)相關(guān),且均通過顯著性檢驗(yàn)。這由于經(jīng)過數(shù)學(xué)變換后,放大某些原始光譜數(shù)據(jù)中細(xì)小的波段信息,所以微分變換后的波段與鹽漬土含水量相關(guān)性比較高。據(jù)此可以篩選出鹽漬土含水量的敏感波段,確定相關(guān)系數(shù)最高值對(duì)應(yīng)的波段(表1)??梢?出各種數(shù)學(xué)變換篩選出的敏感波段波長(zhǎng)均>1 000 nm;一階微分、倒數(shù)一階、對(duì)數(shù)一階和均方根一階都有相同的敏感波段2 357 nm,且此敏感波段的相關(guān)系數(shù)較高。

    3.3 建模與分析

    根據(jù)相關(guān)分析結(jié)果得到鹽漬土含水量的敏感波段,建立相應(yīng)的土壤含水量的估算模型。由于實(shí)驗(yàn)誤差或數(shù)據(jù)異常,53個(gè)樣本中只有44個(gè)樣點(diǎn)進(jìn)行了建模與驗(yàn)證,然后從44個(gè)樣本中隨機(jī)選取34個(gè)樣品構(gòu)建土壤含水量的估算模型,以敏感波段作為自變量,實(shí)測(cè)的土壤含水量作為因變量,并選擇相關(guān)系數(shù)r較大的敏感波段,進(jìn)行多元線性回歸分析,建立鹽漬土含水量估算模型(表2)。

    表1 光譜反射率及其變換形式與土壤含水量敏感波段Tab.1 Spectral reflectance and its transformations,and sensitive bands of soil water content

    表2 不同自變量回歸模型Tab.2 Prediction models based on different characteristic bands

    可見,土壤含水量與所選敏感波段呈良好的線性關(guān)系,原始反射率和經(jīng)各數(shù)學(xué)變換后構(gòu)建的模型均通過了0.05顯著性檢驗(yàn)。其中原始反射率的回歸模型回歸系數(shù)最低(r=0.618)。一階微分回歸模型回歸系數(shù)最高(r=0.858)。依據(jù)回歸系數(shù)和線性回歸方程顯著性檢驗(yàn)統(tǒng)計(jì)量F值最大、均方根誤差最小的原則,本文選取一階微分,對(duì)數(shù)一階,均方根一階進(jìn)行構(gòu)建模型。

    3.3.1 模型驗(yàn)證

    模型的驗(yàn)證主要考慮估算模型的穩(wěn)定性和效果。利用剩下的10個(gè)真實(shí)值作為驗(yàn)證數(shù)據(jù),對(duì)上述所選模型對(duì)其進(jìn)行驗(yàn)證,對(duì)反演模型的估算效果進(jìn)行檢驗(yàn)其中一階微分模型精度r=0.554,對(duì)數(shù)一階微分模型精度r=0.894,均方根一階微分精度r= 0.865。根據(jù)以上結(jié)果剔除一階微分模型,選取對(duì)數(shù)一階微分和均方根一階微分模型進(jìn)行鹽漬土含水量的反演。從模型的估算值和真實(shí)值的線性擬合(圖3)可以看出,對(duì)數(shù)一階模型樣點(diǎn)基本上聚集在對(duì)角線附近,樣本相關(guān)系數(shù)r=0.894。相比之下均方根一階樣點(diǎn)沒有對(duì)數(shù)一階樣點(diǎn)那樣緊湊,樣點(diǎn)相關(guān)系數(shù)r=0.865。除個(gè)別異常點(diǎn),估算值和實(shí)測(cè)值的線性擬合程度均較好。本文建立的鹽漬土含水量遙感監(jiān)測(cè)模型估算能力通過了檢驗(yàn),在一定程度上可以用來估算鹽漬土含水量。

    3.3.2 多元線性回歸與構(gòu)建的耦合估算模型

    與原始反射率相比,不同變換后的光譜波段與含水量相關(guān)性有著不同程度的提高。經(jīng)過以往研究模型的驗(yàn)證可知:構(gòu)建的模型比較單一性,都是一種數(shù)學(xué)變換下的模型構(gòu)建,有一定的不足;因此本文基于以上模型選取不同變換下的敏感波段構(gòu)建模型。為了保證模型數(shù)據(jù)的一致性,模型選取的樣本和上述模型一致。構(gòu)建模型為y=-6.147x1+53.535x2-181.344x3+0.044,式中:y為鹽漬土含水量,x1為倒數(shù)一階變換下2 032 nm處的光譜反射率,x2為對(duì)數(shù)一階變換2 357 nm處的光譜反射率,x3為均方根一階變換2 357 nm處的光譜反射率。通過10個(gè)樣點(diǎn)對(duì)模型進(jìn)行精度驗(yàn)證,可以得到,模型的精度有所提高,相關(guān)系數(shù)r=0.926,比對(duì)數(shù)一階微分模型相關(guān)性r提高了0.032,比均方根一階微分模型相關(guān)系數(shù)如提高了0.061(表3)。

    圖3 模型估算值與實(shí)測(cè)值線性擬合圖(N為樣本數(shù),下同)Fig.3 Linearfit of measured values and predicted values with the model(N:number of samples,the same below)

    表3 不同回歸模型對(duì)比Tab.3 Comparison of different regression models

    相比之前的模型,此模型基于對(duì)數(shù)一階微分模型與比均方根一階微分模型的相關(guān)性最優(yōu)波段,且經(jīng)過驗(yàn)證發(fā)現(xiàn)耦合模型對(duì)鹽漬土含水量的估算是可行的(圖4)。

    圖4 模型估算值與實(shí)測(cè)值線性擬合圖Fig.4 Linearfit of measured values and predicted values with the model

    4 結(jié)論與討論

    1)鹽漬土不同土含水量光譜曲線的形態(tài)與趨勢(shì)基本一致;在可見光波段范圍內(nèi),不同含水量的鹽漬土反射率呈規(guī)律性變化,即隨著土壤含水量的增加,土壤反射率呈下降走勢(shì)。

    2)對(duì)鹽漬土含水量與光譜建模與驗(yàn)證得:經(jīng)建模與驗(yàn)證后得出對(duì)數(shù)一階的模型最好(r=0.804)其次是均方根一階(r=0.798)。對(duì)數(shù)一階微分模型的估算值與實(shí)測(cè)值的相關(guān)系數(shù)最高(r=0.894),其次是均方根一階微分模型,相關(guān)系數(shù)為r=0.865,即鹽漬土含水量模型估算良好。

    3)基于均方根一階微分模型與對(duì)數(shù)一階微分模型建立耦合模型,其方程:y=-6.147x1+ 53.535x2-181.344x3+0.044,其模型效果良好(r= 0.926)。證明了選取不同變換后的敏感波段進(jìn)行建模,對(duì)鹽漬土含水量估算是可行的。

    筆者提出的基于均方根一階微分模型與對(duì)數(shù)一階微分模型下的耦合模型,對(duì)艾比湖流域鹽漬土的治理與監(jiān)測(cè)有著重要意義和作用;但本研究只是在室內(nèi)的一種嘗試,并沒有在室外進(jìn)行驗(yàn)證,因此耦合模型有待于進(jìn)一步檢驗(yàn)與改進(jìn)。

    [1] 李美婷,武紅旗,蔣平安,等.利用土壤的近紅外光譜特征測(cè)定土壤含水量[J].光譜學(xué)與光譜分析,2012, 32(8):2117. LI Meiting,WU Hongqi,JIANG Pingan,et al.Measuring soil water content by using near infrared spectral characteristics of soil[J].Spectroscopy and Spectral Analysis, 2012,32(8):2117.

    [2] 中國(guó)土壤學(xué)會(huì).中國(guó)土壤學(xué)在前進(jìn)[M].北京:中國(guó)農(nóng)業(yè)科技出版社,1995:9. Soil Science Society of China.Progress of soil science in China[M].Beijing:China's Agricultural Science and Technology Press,1995:9.

    [3] ROBINOVE C J,CHAVEZ P S,GEHRING D,et al. Arid land monitoring using Landsat albedo difference images[J].Remote Sensing of Environment,1981,11:133.

    [4] HENRICKSEN B L.Reflections on drought:Ethiopia 1983-1984[J].International Journal of Remote Sensing,1986,7(11):1447.

    [5] 劉志明,張柏,晏明,等.土壤水分與干旱遙感研究的進(jìn)展與趨勢(shì)[J].地球科學(xué)進(jìn)展,2003,18(4):576. LIU Zhiming,ZHANG Bai,YAN Ming,et al.Some research advances and trends on soil moisture and drought monitoring by remote sensing[J].Advance in Earth Sciences,2003,18(4):576.

    [6] LIU Huanjun,ZHANG Yuanzhi,ZHANG Bai.Novel hyper spectral reflectance models for estimating black-soil organic matter in Northeast China[J].Environmental Monitoring and Assessment,2009,154(1/2/3/4): 147.

    [7] LI Y,DEMETRIADES-SHAH T H,KANEMASU E T, et al.Use of second derivatives of canopy reflectance for monitoring prairie vegetation over different soil backgrounds[J].Remote Sensing of Environment,1993,44 (1):81.

    [8] LU Ning,ZHANG Zhi,GAO Yang.Recognition and mapping of soil salinization in arid environment with hyper spectral data[C].Proceedings.2005 IEEE International Geo science and Remote Sensing Symposium,2005. IGArSS'05.IEEE,2005,6:4520.

    [9] 孫小艷,常學(xué)禮,張寧,等.不同取樣單元對(duì)干旱區(qū)綠洲小麥地上生物量光譜估算模型的影響[J].中國(guó)沙漠,2012,32(2):568. SUN Xiaoyan,CHANG Xueli,ZHANG Ning,et al.Impact of different sampling units on ground spectral modelsfor estimating wheat aboveground biomass in oases of arid zone[J].Journal of Desert Research,2012,32(2):568.

    [10] 夏學(xué)齊,季峻峰,陳駿,等.土壤理化參數(shù)的反射光譜分析[J].地學(xué)前緣,2009,16(4):354. XIA Xueqi,JI Junfeng,CHEN Jun,et al.Analysis of soil physical and chemical properties by reflectance spectroscopy[J].Earth Science Frontiers,2009,16 (4):354.

    [11] 吳亞坤,楊勁松,李曉明.基于光譜指數(shù)與EM38的土壤鹽分空間變異性研究[J].光譜學(xué)與光譜分析, 2009(4):1023. WU Yakun,YANG Jinsong,LI Xiaoming.Study on spatial variability of soil salinity based on spectral indices and EM38 readings[J].Spectroscopy and Spectral Analysis,2009,29(4):1023.

    [12] 劉煥軍,張柏,宋開山,等.黑土土壤水分光譜響應(yīng)特征與模型[J].中國(guó)科學(xué)院研究生院學(xué)報(bào),2008,25 (4):503. LIU Huanjun,ZHANG Bai,SONG Kaishan,et al.Soil moisture spectral analysis and its spectral model[J]. Journal of the Graduate School of the Chinese Academy of Sciences,2008,25(4):503.

    [13] 何挺,王靜,程燁,等.土壤水分光譜特征研究[J].土壤學(xué)報(bào),2006,43(6):1027. HE Ting,WANG Jing,CHENG Ye,et al.Soil water spectrum characteristic research[J].Acta Pedologica Sinica,2006,43(6):1027.

    [14] 翁永玲,戚浩平,方洪賓,等.基于PLSR方法的青海茶卡-共和盆地土壤鹽分高光譜遙感反演[J].土壤學(xué)報(bào),2010,47(6):1255. WENG Yongling,QIE Haoping,FANG Hongbin,et al. PLSR based hyper spectral remote sensing retrieval of soil salinity of Chaka·Gonghe basin in Qinghai province[J].Acta Pedologica Sinica,2010,47(6):1255.

    [15] 姜紅濤,特依拜,克里木,等.艾比湖流域NDWI變化及其與降水、溫度的關(guān)系[J].中國(guó)沙漠,2014,34 (6):1678. JIANG Hongtao,TIYIP,KELIMU,et al.Responses of NDVI to the variation of precipitation and temperature the Ebinur Lake Basin[J].Journal of Desert Research, 2014,34(6):1678.

    [16] 張小龍.艾比湖流域氣候變化及其徑流響應(yīng)[J].鹽湖研究,2011,19(2):11. ZHANG Xiaolong Climatic variation and runoff response in Ebinur Lake Catchment[J].Journal of Salt Lake ResearcH,2011,19(2):11.

    [17] 王昌佐,王紀(jì)華,王錦地,等.裸土表層含水量高光譜遙感的最佳波段選擇[J].遙感信息,2003(4):33. WANG Changzuo,WANG Jihua,WANG Jindi,et al. The choice of best detecting band for hyper spectral remote sensing on surface water content of bare soil[J]. Remote Sensing Information,2003.(4):33.

    [18] 關(guān)紅,賈科利,張至椅,等.鹽漬化土壤光譜特征分析與建模[J].國(guó)土資源遙感,2015,27(2)100. GUAN Hong,JIA Keli,ZHANG Zhinan.Research on remote sensing monitoring model of soil salinization based on spectrum characteristic analysis[J].Remote Sensing for Land and resources,2015,27(2):100.

    [19] 趙振亮,特依拜,丁建麗,等.新疆典型綠洲土壤電導(dǎo)率和PH值的光譜響應(yīng)特征[J].中國(guó)沙漠,2013,33 (5):1413. ZHAO Zhenliang,TIYIP,DING Jianli,et al.Characteristics of spectral responding to soil electrical conductivity PH in the typical Oasis of Xinjiang[J].Journal of Desert Research,2013,33(5):1413.

    [20] 顧燕,張鷹,李歡.基于實(shí)測(cè)光譜的潮灘土壤含水量遙感反演模型研究[J].濕地科學(xué),2013,11(2): 167. GU Yan,ZHANG Ying,LI Huan,et al.Remote sensing retrieval model on soil moisture content of tidal flat based on measured spectra[J].Wetland Science, 2012,11(2):168.

    [21] 劉曾媛,貢璐,張雪妮,等.克里雅河流域于田綠洲土壤酶活性c理化因子相關(guān)性分析[J].中國(guó)土壤與肥料,2014,4(7):35. LIU Zeng yuan,GONG Lu,ZHANG Xueni,et al.Soil enzyme activities and their relationship with physicochemical factors in Yutian oasis of the Keriya river watershed[J].Soil and Fertilizer Sciences in China,2014, 4(7):35.

    Modeling and analysis of spectral characteristic of soil water content in the salinized soil of Ebinur Lake Watershed

    ZHANG Haiwei1,2,ZHANG Fei1,2,3,LI Zhe1,2,JING Yunqing1,2
    (1.College of Resources&Environmental Science,Xinjiang University,830046,Urumqi,China; 2.Key Laboratory of Oasis Ecology,830046,Urumqi,China;3.General Institutes of Higher Learning Key Laboratory of Smart City and Environmental Modeling,Xinjiang University,830046,Urumqi,China)

    [Background]With the development of spatial information science,the hyperspectral remote sensing becomes more and more important in nowadays.Studying on spectral characteristic of soil water content is an important work,for it is the base of monitoring remote sensing.This study aims to investigate the spectral characteristic of soil water content and the relationship between the hyperspectral data and the soil water content,to explore a rapid and accurate method for estimating soil water contents, and to establish the hyperspectral remote sensing monitoring model for saline soil water in the arid area of the Ebinur Lake Watershed.[Methods]This study took saline soil with different water contents in Ebinur Lake Watershed as the research object,used the spectral reflectance transformation and multivariate statistical analysis methods(MSAM)to analyze the spectral characteristic of soil water content,and built models.[Results]The result showed that with the increase of soil water content,thereflectance of soil declined.In a certain range,the longer the wavelength,the higher the correlation between the soil spectral reflectance and soil water content;the soil spectral reflectance at the wavelength of 1 937 nm(r=-0.636)had the highest correlation with water content.After 8 transformation of soil spectral reflectance,the correlation of logarithmic first order differential sensitive band of 2 357 nm was the best(r=-0.808 6).Subsequently,MSAM were applied to analyze the correlation between spectrum and salinized soil with different water contents,then the spectral sensitive band of the soil was screened,and the correlation model was established.The result indicated that the model set up by Logarithm First Order Differential at the wavelength of 2 024 nm and 2 357 nm and the Root Mean Square First Order Differential at the wavelength of 1 972 nm and 2 357 nm was the best,the correlation coefficientrwas 0.894 and 0.865 respectively.Based on the above established models,the authors constructed a new model,and the correlationrof which was 0.926,increased 0.032 against the Logarithm First Order Differential model,and increased 0.061 compared to the Root Mean Square First Order Differential model.[Conclusions]Therefore,the estimation model of soil water content is feasible.In addition,this study provides a new model for the indoor hyperspectral estimation of soil water content in Ebinur Lake Watershed,and it also provides a theoretical and technical reference for the hyperspectral quantitative estimation of soil water content,which has certain guiding significance for hyperspectral remote sensing.

    spectrum characteristic;water content;modeling;Ebinur Lake Watershed

    O433.4

    :A

    :2096-2673(2017)01-0008-07

    10.16843/j.sswc.2017.01.002

    2016- 07- 03

    2016- 11- 27

    項(xiàng)目名稱:國(guó)家自然科學(xué)基金“新疆聯(lián)合基金本地優(yōu)秀青年人才培養(yǎng)專項(xiàng)”(U1503302)

    張海威(1990—),男,碩士。主要研究方向:干旱區(qū)遙感應(yīng)用。E-mail:1522634426@qq.com

    ?通信作者簡(jiǎn)介:張飛(1980—),男,博士,副教授。主要研究方向:干旱區(qū)資源環(huán)境遙感。E-mail:zhangfei3s@163.com

    猜你喜歡
    鹽漬方根微分
    燕麥對(duì)濱海鹽漬土的適應(yīng)性及納鹽改土效果
    方根拓展探究
    鹽漬池深度對(duì)漬菜理化特性和微生物的影響
    擬微分算子在Hp(ω)上的有界性
    上下解反向的脈沖微分包含解的存在性
    均方根嵌入式容積粒子PHD 多目標(biāo)跟蹤方法
    借助微分探求連續(xù)函數(shù)的極值點(diǎn)
    揭開心算方根之謎
    談鹽漬土地區(qū)建筑的地基處理設(shè)計(jì)
    對(duì)不定積分湊微分解法的再認(rèn)識(shí)
    丝袜脚勾引网站| 狠狠婷婷综合久久久久久88av| 天堂中文最新版在线下载| 亚洲人成电影观看| 亚洲av男天堂| 一区二区三区四区激情视频| 亚洲熟女毛片儿| 日韩欧美国产一区二区入口| 高清视频免费观看一区二区| 91成年电影在线观看| 日韩欧美免费精品| 色精品久久人妻99蜜桃| 色婷婷av一区二区三区视频| 脱女人内裤的视频| 老司机在亚洲福利影院| 久久久久久久国产电影| 两个人免费观看高清视频| 捣出白浆h1v1| 99国产精品免费福利视频| 久久人妻福利社区极品人妻图片| 视频区图区小说| 亚洲精品一区蜜桃| 欧美黑人欧美精品刺激| 国产欧美日韩精品亚洲av| 中文字幕人妻丝袜制服| 精品国产一区二区三区久久久樱花| 国产精品久久久久久人妻精品电影 | 亚洲精品一二三| 大陆偷拍与自拍| 亚洲人成电影观看| 欧美国产精品一级二级三级| 欧美精品啪啪一区二区三区 | 亚洲精品粉嫩美女一区| 大香蕉久久成人网| 成人国语在线视频| 日日夜夜操网爽| 国产精品久久久久久精品电影小说| 老熟妇仑乱视频hdxx| 一区福利在线观看| 日韩视频一区二区在线观看| 亚洲色图 男人天堂 中文字幕| 国产男女内射视频| 9色porny在线观看| 黄色视频在线播放观看不卡| 午夜影院在线不卡| 99国产综合亚洲精品| 91国产中文字幕| 各种免费的搞黄视频| 久久99热这里只频精品6学生| 男人爽女人下面视频在线观看| 日韩 亚洲 欧美在线| 老熟女久久久| 国产精品九九99| 男女边摸边吃奶| av片东京热男人的天堂| 青春草亚洲视频在线观看| 一本大道久久a久久精品| 成人av一区二区三区在线看 | 欧美黄色淫秽网站| 国产免费视频播放在线视频| 亚洲欧美精品综合一区二区三区| av线在线观看网站| 亚洲精品国产精品久久久不卡| av超薄肉色丝袜交足视频| 夜夜夜夜夜久久久久| av一本久久久久| 在线观看舔阴道视频| 妹子高潮喷水视频| 欧美老熟妇乱子伦牲交| a级毛片在线看网站| 亚洲美女黄色视频免费看| 久久久国产一区二区| 精品卡一卡二卡四卡免费| 热re99久久精品国产66热6| 日本五十路高清| 中文字幕人妻丝袜一区二区| 香蕉国产在线看| 夫妻午夜视频| 女人爽到高潮嗷嗷叫在线视频| 在线亚洲精品国产二区图片欧美| 精品人妻熟女毛片av久久网站| 又大又爽又粗| 色综合欧美亚洲国产小说| av在线播放精品| 午夜两性在线视频| 国产免费现黄频在线看| 亚洲国产av新网站| 亚洲 欧美一区二区三区| 我要看黄色一级片免费的| 亚洲 国产 在线| 高清视频免费观看一区二区| 久久精品成人免费网站| 中文字幕人妻熟女乱码| 精品少妇久久久久久888优播| 欧美精品亚洲一区二区| 午夜影院在线不卡| 美女高潮到喷水免费观看| 欧美亚洲 丝袜 人妻 在线| 91麻豆av在线| 亚洲国产精品成人久久小说| 午夜免费观看性视频| 欧美激情高清一区二区三区| 首页视频小说图片口味搜索| 亚洲自偷自拍图片 自拍| 亚洲 国产 在线| 在线天堂中文资源库| 国产黄频视频在线观看| 两个人看的免费小视频| 亚洲欧美激情在线| 1024香蕉在线观看| 免费看十八禁软件| 久久影院123| 美女国产高潮福利片在线看| 精品亚洲成国产av| 欧美激情高清一区二区三区| 热99re8久久精品国产| 青春草亚洲视频在线观看| a级毛片黄视频| 午夜影院在线不卡| av网站免费在线观看视频| 91麻豆av在线| 在线 av 中文字幕| 亚洲av成人不卡在线观看播放网 | 99久久国产精品久久久| 中文字幕高清在线视频| 亚洲国产成人一精品久久久| 精品亚洲成a人片在线观看| 亚洲成人国产一区在线观看| 免费观看av网站的网址| 国产男人的电影天堂91| 曰老女人黄片| 考比视频在线观看| 欧美精品人与动牲交sv欧美| 日韩制服骚丝袜av| 亚洲专区字幕在线| 丝袜喷水一区| 精品国产乱码久久久久久男人| 精品国产一区二区久久| 女人爽到高潮嗷嗷叫在线视频| 国产在线一区二区三区精| av有码第一页| 夜夜夜夜夜久久久久| 亚洲国产精品一区三区| 狂野欧美激情性bbbbbb| 欧美黄色片欧美黄色片| 国产精品成人在线| 色视频在线一区二区三区| 国产精品久久久久久精品古装| 国产91精品成人一区二区三区 | 在线天堂中文资源库| 制服诱惑二区| 免费一级毛片在线播放高清视频 | 色老头精品视频在线观看| 亚洲欧美色中文字幕在线| 香蕉丝袜av| 成人国产av品久久久| 国产男人的电影天堂91| 一级a爱视频在线免费观看| 午夜免费观看性视频| 美女大奶头黄色视频| 午夜激情久久久久久久| 美女国产高潮福利片在线看| 欧美激情 高清一区二区三区| 亚洲五月婷婷丁香| 亚洲国产精品一区三区| 老司机在亚洲福利影院| 亚洲第一青青草原| av福利片在线| 超碰成人久久| 精品视频人人做人人爽| 久久久久国内视频| 欧美久久黑人一区二区| 亚洲专区中文字幕在线| 久久人妻熟女aⅴ| 久久久久国产一级毛片高清牌| 两性午夜刺激爽爽歪歪视频在线观看 | 少妇精品久久久久久久| 叶爱在线成人免费视频播放| 少妇的丰满在线观看| 欧美 亚洲 国产 日韩一| 日本五十路高清| 久久中文字幕一级| 狂野欧美激情性xxxx| 国产高清国产精品国产三级| 又紧又爽又黄一区二区| 国产欧美亚洲国产| 99热网站在线观看| 高清欧美精品videossex| 亚洲成人免费av在线播放| 免费女性裸体啪啪无遮挡网站| 亚洲成国产人片在线观看| 国产精品久久久久成人av| 国产日韩欧美在线精品| 在线看a的网站| 一级片'在线观看视频| 老司机在亚洲福利影院| 久久中文看片网| 欧美精品一区二区免费开放| 欧美另类一区| 啦啦啦啦在线视频资源| 亚洲九九香蕉| 正在播放国产对白刺激| 亚洲国产成人一精品久久久| 男女之事视频高清在线观看| 在线av久久热| 精品福利观看| 国产又色又爽无遮挡免| 久久精品成人免费网站| 99久久综合免费| 久久久国产欧美日韩av| 老鸭窝网址在线观看| 多毛熟女@视频| 成年人午夜在线观看视频| 国产主播在线观看一区二区| 成人黄色视频免费在线看| 黄色a级毛片大全视频| 成人国产一区最新在线观看| 午夜福利免费观看在线| 可以免费在线观看a视频的电影网站| av国产精品久久久久影院| 在线十欧美十亚洲十日本专区| 日本五十路高清| 亚洲,欧美精品.| 一区二区三区精品91| 精品少妇久久久久久888优播| 午夜精品国产一区二区电影| 亚洲av电影在线进入| 日本撒尿小便嘘嘘汇集6| 大片免费播放器 马上看| 亚洲成人免费av在线播放| 纵有疾风起免费观看全集完整版| 青青草视频在线视频观看| 777久久人妻少妇嫩草av网站| 1024香蕉在线观看| 极品人妻少妇av视频| 成人影院久久| 热re99久久精品国产66热6| 亚洲精品粉嫩美女一区| 爱豆传媒免费全集在线观看| 亚洲中文字幕日韩| 天天添夜夜摸| 国产片内射在线| 亚洲成人免费av在线播放| 国产一区二区三区av在线| 久久精品亚洲熟妇少妇任你| 国产精品一区二区在线观看99| 亚洲欧洲精品一区二区精品久久久| 精品亚洲成a人片在线观看| 成年人黄色毛片网站| 老司机亚洲免费影院| 日本黄色日本黄色录像| 青春草亚洲视频在线观看| 欧美av亚洲av综合av国产av| 麻豆国产av国片精品| a 毛片基地| 日韩熟女老妇一区二区性免费视频| av一本久久久久| 亚洲成人免费av在线播放| 男人操女人黄网站| 国产日韩欧美视频二区| 欧美精品一区二区大全| 在线观看免费高清a一片| 国产一级毛片在线| 叶爱在线成人免费视频播放| 高清在线国产一区| 嫩草影视91久久| 男女之事视频高清在线观看| 精品一区二区三区av网在线观看 | 午夜激情av网站| 欧美日韩亚洲高清精品| 久久人妻福利社区极品人妻图片| 18禁裸乳无遮挡动漫免费视频| 男男h啪啪无遮挡| 涩涩av久久男人的天堂| 大香蕉久久成人网| 午夜福利影视在线免费观看| 美女高潮喷水抽搐中文字幕| e午夜精品久久久久久久| 欧美日韩av久久| 中文字幕人妻丝袜一区二区| 两人在一起打扑克的视频| 老司机在亚洲福利影院| av有码第一页| 色综合欧美亚洲国产小说| 亚洲欧美清纯卡通| av视频免费观看在线观看| 男人操女人黄网站| 老司机在亚洲福利影院| 大香蕉久久网| 免费av中文字幕在线| 亚洲成人免费电影在线观看| 成在线人永久免费视频| 黑人猛操日本美女一级片| 老汉色∧v一级毛片| 中文字幕另类日韩欧美亚洲嫩草| 少妇的丰满在线观看| 伦理电影免费视频| 久久久久久久国产电影| 黑丝袜美女国产一区| 日本wwww免费看| 成年人午夜在线观看视频| 欧美黑人欧美精品刺激| 欧美日韩亚洲国产一区二区在线观看 | 12—13女人毛片做爰片一| 久久中文字幕一级| xxxhd国产人妻xxx| 美女中出高潮动态图| 国产欧美日韩一区二区精品| 我要看黄色一级片免费的| 蜜桃在线观看..| 成年人黄色毛片网站| 中文字幕精品免费在线观看视频| 亚洲第一青青草原| 亚洲第一欧美日韩一区二区三区 | 国产高清videossex| 男人舔女人的私密视频| 久久天躁狠狠躁夜夜2o2o| 窝窝影院91人妻| 婷婷色av中文字幕| 满18在线观看网站| av国产精品久久久久影院| 新久久久久国产一级毛片| 一级,二级,三级黄色视频| 中文字幕高清在线视频| 国产精品国产av在线观看| 国产极品粉嫩免费观看在线| 亚洲五月色婷婷综合| 久久国产精品大桥未久av| √禁漫天堂资源中文www| 天天躁夜夜躁狠狠躁躁| 麻豆av在线久日| netflix在线观看网站| 亚洲第一青青草原| 女警被强在线播放| 久久久久国产精品人妻一区二区| 叶爱在线成人免费视频播放| 亚洲精品久久午夜乱码| 免费少妇av软件| 欧美日韩一级在线毛片| 九色亚洲精品在线播放| 久久中文字幕一级| 人人澡人人妻人| 巨乳人妻的诱惑在线观看| 午夜精品久久久久久毛片777| 美女国产高潮福利片在线看| 免费高清在线观看日韩| 一级毛片女人18水好多| 午夜免费成人在线视频| 最新的欧美精品一区二区| netflix在线观看网站| 亚洲国产毛片av蜜桃av| 亚洲美女黄色视频免费看| 欧美 亚洲 国产 日韩一| 亚洲精品粉嫩美女一区| 精品少妇内射三级| 日韩三级视频一区二区三区| 操美女的视频在线观看| 亚洲精品中文字幕在线视频| 最新的欧美精品一区二区| 亚洲国产欧美网| 精品乱码久久久久久99久播| e午夜精品久久久久久久| 婷婷色av中文字幕| 最近最新免费中文字幕在线| 欧美性长视频在线观看| 90打野战视频偷拍视频| 国产精品欧美亚洲77777| 乱人伦中国视频| 91国产中文字幕| 人妻人人澡人人爽人人| 19禁男女啪啪无遮挡网站| 国产亚洲午夜精品一区二区久久| 国产在线观看jvid| 丝袜人妻中文字幕| 欧美精品啪啪一区二区三区 | 国产在线视频一区二区| av天堂久久9| 免费少妇av软件| 三级毛片av免费| 亚洲欧美激情在线| 俄罗斯特黄特色一大片| 欧美久久黑人一区二区| 欧美黑人精品巨大| 免费不卡黄色视频| 亚洲精品国产av蜜桃| 丝袜人妻中文字幕| 老司机亚洲免费影院| 丰满饥渴人妻一区二区三| 亚洲av电影在线进入| 精品国产国语对白av| 99国产综合亚洲精品| 中文字幕精品免费在线观看视频| 交换朋友夫妻互换小说| 亚洲专区国产一区二区| 国产黄频视频在线观看| 久久久久久人人人人人| e午夜精品久久久久久久| tube8黄色片| 色综合欧美亚洲国产小说| 青春草视频在线免费观看| 久久久久国产一级毛片高清牌| 精品一区二区三卡| 欧美精品av麻豆av| 十八禁网站网址无遮挡| 99久久99久久久精品蜜桃| 成在线人永久免费视频| 欧美激情久久久久久爽电影 | 婷婷丁香在线五月| 国产又色又爽无遮挡免| 脱女人内裤的视频| 国产av精品麻豆| 亚洲国产欧美网| 中文字幕另类日韩欧美亚洲嫩草| 色老头精品视频在线观看| a级毛片在线看网站| 亚洲熟女毛片儿| 国产欧美日韩一区二区精品| 久久人人97超碰香蕉20202| 一区二区av电影网| 无遮挡黄片免费观看| 777米奇影视久久| 亚洲精品久久成人aⅴ小说| 欧美少妇被猛烈插入视频| 亚洲成人国产一区在线观看| 啦啦啦在线免费观看视频4| 精品久久久久久久毛片微露脸 | 日本精品一区二区三区蜜桃| 91精品国产国语对白视频| 国产成人精品在线电影| www日本在线高清视频| 大型av网站在线播放| av免费在线观看网站| 一区二区三区精品91| 亚洲国产成人一精品久久久| 首页视频小说图片口味搜索| 少妇被粗大的猛进出69影院| 久久青草综合色| 色94色欧美一区二区| 中文字幕精品免费在线观看视频| 久久天躁狠狠躁夜夜2o2o| 999久久久精品免费观看国产| 中文字幕人妻熟女乱码| 丝袜喷水一区| 欧美精品av麻豆av| 亚洲精品国产精品久久久不卡| 欧美日韩国产mv在线观看视频| 国产精品av久久久久免费| 99国产精品一区二区三区| 一级,二级,三级黄色视频| 午夜两性在线视频| 久久久精品国产亚洲av高清涩受| 精品国产一区二区三区四区第35| av天堂久久9| 另类亚洲欧美激情| xxxhd国产人妻xxx| 精品国产乱码久久久久久男人| 曰老女人黄片| 青草久久国产| 欧美+亚洲+日韩+国产| 18禁裸乳无遮挡动漫免费视频| 亚洲国产看品久久| 两个人看的免费小视频| 久久久久精品人妻al黑| 一本大道久久a久久精品| 色精品久久人妻99蜜桃| 亚洲国产成人一精品久久久| 18禁裸乳无遮挡动漫免费视频| 亚洲成人免费av在线播放| 免费久久久久久久精品成人欧美视频| a级片在线免费高清观看视频| 窝窝影院91人妻| 国产亚洲精品久久久久5区| 777久久人妻少妇嫩草av网站| 中文欧美无线码| 乱人伦中国视频| 2018国产大陆天天弄谢| 免费在线观看完整版高清| 热re99久久国产66热| 国产成人精品无人区| 国产精品国产三级国产专区5o| 日韩有码中文字幕| 少妇精品久久久久久久| www.精华液| 亚洲国产精品999| 99国产精品99久久久久| 久久99一区二区三区| 亚洲中文字幕日韩| 色婷婷av一区二区三区视频| 日韩欧美一区二区三区在线观看 | 欧美日韩一级在线毛片| 国产一区有黄有色的免费视频| 中文字幕最新亚洲高清| 成人亚洲精品一区在线观看| 美女国产高潮福利片在线看| 国产熟女午夜一区二区三区| 少妇粗大呻吟视频| 免费观看a级毛片全部| 好男人电影高清在线观看| 午夜激情久久久久久久| 久久久久精品人妻al黑| 亚洲va日本ⅴa欧美va伊人久久 | 啦啦啦在线免费观看视频4| 性色av乱码一区二区三区2| 亚洲精品成人av观看孕妇| 最新的欧美精品一区二区| 性色av一级| 人人妻人人澡人人看| 成年动漫av网址| 亚洲色图综合在线观看| 满18在线观看网站| 12—13女人毛片做爰片一| 可以免费在线观看a视频的电影网站| a级片在线免费高清观看视频| 99久久99久久久精品蜜桃| 色精品久久人妻99蜜桃| 日韩制服丝袜自拍偷拍| √禁漫天堂资源中文www| 大香蕉久久网| 黄色视频在线播放观看不卡| 亚洲人成电影免费在线| 成年动漫av网址| 窝窝影院91人妻| 亚洲七黄色美女视频| 亚洲精品在线美女| 爱豆传媒免费全集在线观看| 99精品欧美一区二区三区四区| 国产亚洲精品一区二区www | 天天操日日干夜夜撸| 另类精品久久| 首页视频小说图片口味搜索| 高清视频免费观看一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美 亚洲 国产 日韩一| 久久精品亚洲av国产电影网| av在线老鸭窝| 精品一区二区三区av网在线观看 | 国产色视频综合| 一个人免费看片子| 色视频在线一区二区三区| 51午夜福利影视在线观看| 999精品在线视频| 啪啪无遮挡十八禁网站| 久久人人97超碰香蕉20202| 午夜成年电影在线免费观看| 色94色欧美一区二区| 97精品久久久久久久久久精品| 日本wwww免费看| 欧美午夜高清在线| 美国免费a级毛片| 一区二区三区精品91| 久久久国产欧美日韩av| 50天的宝宝边吃奶边哭怎么回事| 51午夜福利影视在线观看| 悠悠久久av| 大片电影免费在线观看免费| 亚洲精品粉嫩美女一区| av免费在线观看网站| 久久这里只有精品19| 1024香蕉在线观看| 又黄又粗又硬又大视频| 麻豆国产av国片精品| avwww免费| 丝袜美足系列| 一区二区三区精品91| 欧美精品高潮呻吟av久久| 50天的宝宝边吃奶边哭怎么回事| 国产福利在线免费观看视频| 纯流量卡能插随身wifi吗| 亚洲欧洲日产国产| 日韩视频一区二区在线观看| 啦啦啦啦在线视频资源| 人妻久久中文字幕网| 亚洲av日韩精品久久久久久密| 亚洲专区国产一区二区| 亚洲免费av在线视频| 男女床上黄色一级片免费看| 老司机影院成人| 视频区图区小说| 亚洲精品一二三| 不卡av一区二区三区| 日韩有码中文字幕| 国产不卡av网站在线观看| 欧美人与性动交α欧美软件| 最黄视频免费看| 在线观看免费视频网站a站| 亚洲av成人一区二区三| 青草久久国产| av欧美777| 精品卡一卡二卡四卡免费| 国产精品成人在线| 日韩电影二区| 一区二区三区四区激情视频| 亚洲精品自拍成人| 亚洲国产毛片av蜜桃av| 中文字幕av电影在线播放| 亚洲激情五月婷婷啪啪| av天堂在线播放| 久久久久久久精品精品| 精品福利永久在线观看| 韩国精品一区二区三区| 久久热在线av| 国产福利在线免费观看视频| 乱人伦中国视频| 高潮久久久久久久久久久不卡| 精品福利永久在线观看| videos熟女内射| 国产av国产精品国产| 亚洲精品久久久久久婷婷小说| 欧美精品高潮呻吟av久久| 少妇猛男粗大的猛烈进出视频| 免费人妻精品一区二区三区视频| 国产又爽黄色视频| 王馨瑶露胸无遮挡在线观看| 亚洲熟女精品中文字幕| 久久精品国产亚洲av香蕉五月 |