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

    基于無(wú)人機(jī)遙感圖像紋理與植被指數(shù)的土壤含鹽量反演

    2023-08-22 06:56:22向友珍李汪洋安嘉琪陳俊英
    關(guān)鍵詞:鹽漬化含鹽量植被指數(shù)

    向友珍 李汪洋 臺(tái) 翔 安嘉琪 王 辛 陳俊英

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

    0 引言

    土壤鹽漬化是指在自然和人為作用下,鹽分在土壤表層不斷積累的現(xiàn)象,是造成土壤退化、農(nóng)業(yè)減產(chǎn)和生態(tài)環(huán)境惡化的重要因素[1]。河套灌區(qū)土壤鹽漬化問(wèn)題嚴(yán)重,對(duì)鹽漬化土壤進(jìn)行治理和開發(fā)利用是農(nóng)業(yè)可持續(xù)發(fā)展與生態(tài)文明建設(shè)高質(zhì)量發(fā)展的有力保障[2]。實(shí)時(shí)、高效地監(jiān)測(cè)土壤鹽漬化狀況,是治理土壤鹽漬化的重要前提。

    人工定點(diǎn)采樣法是土壤鹽漬化監(jiān)測(cè)的方式之一,該方法監(jiān)測(cè)周期長(zhǎng)、工作量大,且監(jiān)測(cè)的精度與采樣點(diǎn)的數(shù)量、布設(shè)方式有很大關(guān)系,很難反映鹽漬化在空間上的實(shí)際分布。隨著科學(xué)技術(shù)的不斷發(fā)展,遙感被廣泛應(yīng)用于鹽漬化研究和監(jiān)測(cè),為大面積時(shí)空分布動(dòng)態(tài)化的鹽漬化監(jiān)測(cè)提供了新的途徑。然而廣泛可用的衛(wèi)星圖像無(wú)法提供高空間分辨率和數(shù)據(jù)采集時(shí)間的靈活性,在農(nóng)業(yè)應(yīng)用方面有很大的局限性[3]。

    除衛(wèi)星遙感外,利用無(wú)人機(jī)平臺(tái)搭載的小型多光譜遙感設(shè)備對(duì)區(qū)域土壤鹽漬化監(jiān)測(cè)也是重要手段[4-5]。與衛(wèi)星圖像相比,基于無(wú)人機(jī)多光譜圖像的鹽漬化監(jiān)測(cè)具有空間分辨率高、光譜分辨率強(qiáng)、波段連續(xù)性強(qiáng)等優(yōu)點(diǎn),可以獲得多維、高精度的鹽漬化檢測(cè)信息,實(shí)現(xiàn)對(duì)土壤鹽漬化的動(dòng)態(tài)監(jiān)測(cè);同時(shí),無(wú)人機(jī)傳感器攜帶多光譜相機(jī)平臺(tái),具有實(shí)時(shí)、高分辨率、移動(dòng)性和靈活性的優(yōu)勢(shì),在高精度鹽漬化監(jiān)測(cè)方面具有優(yōu)勢(shì)[6]。劉旭輝等[7]利用無(wú)人機(jī)多光譜遙感影像建立的機(jī)器學(xué)習(xí)模型反演土壤含鹽量,得出不同季節(jié)的土壤含鹽量會(huì)有所不同。劉楠等[8]將局部區(qū)域的相關(guān)性分析和多維光譜變化特性采用二進(jìn)制編碼的形式來(lái)表征不同類型地物紋理特征,實(shí)現(xiàn)了區(qū)分和判別影像上不同紋理的目的,并證明該方法對(duì)提取地物紋理特性具有一定可行性。萬(wàn)亮等[9]利用無(wú)人機(jī)多光譜融合植被指數(shù)和紋理特征提高了水稻含水量預(yù)測(cè)結(jié)果;周聰?shù)萚10]通過(guò)實(shí)踐證明,高分辨率遙感圖像的紋理特征可以作為估算植被生長(zhǎng)參數(shù)的更有效指標(biāo)。陳鵬飛等[11]通過(guò)實(shí)踐證實(shí)了基于無(wú)人機(jī)多光譜圖像的剔除土壤背景和增加紋理特征處理可以提高棉花植株氮濃度的反演精度。但是基于無(wú)人機(jī)多光譜圖像剔除土壤背景和紋理特征是否可以提高土壤含鹽量精度方面研究較少。

    本文以河套灌區(qū)沙壕渠灌域4塊不同土壤鹽漬化程度的典型鹽漬化試驗(yàn)地為試驗(yàn)研究區(qū),以灌區(qū)典型地物為研究對(duì)象,在灌區(qū)作物生長(zhǎng)期內(nèi)采集試驗(yàn)區(qū)鹽漬土樣本并獲取土壤鹽漬化信息。同時(shí)利用無(wú)人機(jī)多光譜采集遙感數(shù)據(jù),對(duì)比研究基于鹽漬化土壤和耐鹽作物光譜的土壤含鹽量反演模型的性能,建立基于無(wú)人機(jī)多光譜遙感平臺(tái)的區(qū)域鹽漬化監(jiān)測(cè)體系,為區(qū)域鹽漬化研究提供可靠依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    本研究區(qū)域位于巴彥淖爾市沙壕渠灌域(40°52′~41°00′N,107°05′~107°10′E),隸屬于河套灌區(qū)解放閘灌域,如圖1所示。灌域因不合理的灌排方式,氣候、土質(zhì)和地貌等因素的綜合影響,土壤鹽漬化問(wèn)題突出。區(qū)內(nèi)種植作物以向日葵、玉米等耐鹽糧油作物為主。沙壕渠灌域土壤類型為粉壤土、砂壤土和壤土。沙壕渠灌域面積約52.4 km2,南北跨度約15.1 km,東西橫跨4.2 km,為典型的溫帶大陸性氣候,多年平均氣溫3.7~7.6℃,多年平均降水量210~290 mm,多年平均年蒸發(fā)量2 100~3 080 mm,多年平均日照時(shí)長(zhǎng)為3 000~3 200 h。

    1.2 試驗(yàn)地布設(shè)

    將試驗(yàn)地以土壤含鹽量為基礎(chǔ)從低到高依次分為1、2、3、4,每塊試驗(yàn)地面積為15.5 hm2左右。4塊試驗(yàn)地主要種植作物為向日葵和玉米。每塊試驗(yàn)地均勻布設(shè)30個(gè)土壤采樣點(diǎn)。采樣點(diǎn)布設(shè)見文獻(xiàn)[12]。

    1.3 光譜圖像采集和處理

    使用的無(wú)人機(jī)為深圳市大疆創(chuàng)新科技有限公司生產(chǎn)的M600型六旋翼無(wú)人機(jī),其最大上升速度5 m/s,最大下降速度3 m/s,最大飛行速度18 m/s,飛行承載質(zhì)量6 000 g,飛行高度2 500 m,單次飛行時(shí)間35~40 min。多光譜遙感相機(jī)采用美國(guó)Tetracam公司生產(chǎn)的6通道Micro-MCA多光譜相機(jī),包括藍(lán)光波段、綠光波段、紅光波段、紅邊波段、近紅外1波段、近紅外2波段共6個(gè)遙感波段,波長(zhǎng)分別為490、550、680、720、800、900 nm。試驗(yàn)時(shí)間為2022年7月16—20日。每次試驗(yàn)均在11:00—14:00進(jìn)行,試驗(yàn)日晴朗無(wú)風(fēng),以確保充分的輻射強(qiáng)度,盡量減小植被陰影對(duì)光譜的影響。根據(jù)提前規(guī)劃好的航線,設(shè)置無(wú)人機(jī)飛行高度120 m,對(duì)應(yīng)多光譜相機(jī)分辨率為6.5 cm,相機(jī)拍攝速率為18~19幅/min,每次試驗(yàn)均設(shè)有白板進(jìn)行圖像標(biāo)定。

    1.4 光譜指數(shù)計(jì)算

    光譜指數(shù)是綜合考慮地物的各波段光譜特征,對(duì)不同波段反射率進(jìn)行數(shù)學(xué)變換組合,以增強(qiáng)地物特定的信息[13]。分別選擇8種植被指數(shù)和10種鹽分指數(shù),其計(jì)算公式如表1所示。

    表1 光譜指數(shù)Tab.1 Spectral index

    1.5 灰度共生矩陣法

    灰度共生矩陣(GLCM)是一種基于統(tǒng)計(jì)數(shù)據(jù)的圖像紋理特征提取方法,由HARALICK等[26]提出。

    GLCM的元素是不同組合出現(xiàn)的頻數(shù),為了計(jì)算圖像紋理的特征參數(shù),對(duì)GLCM進(jìn)行歸一化。已知i、j分別為(x,y)、(x+Δx,y+Δy)的像素值,(i,j)的個(gè)數(shù)為V(i,j),將歸一化GLCM記為P,則

    (1)

    (2)

    式中Vi,j——矩陣元素N——矩陣元素個(gè)數(shù)

    基于歸一化灰度共生矩陣,計(jì)算圖像紋理特征參數(shù),本研究采用8個(gè)紋理特征參數(shù),其計(jì)算公式和特性如表2所示。

    表2 圖像紋理特征參數(shù)Tab.2 Texture feature parameters of images

    1.6 Otsu算法

    Otsu算法是一種計(jì)算簡(jiǎn)單、自適應(yīng)強(qiáng)且已得到最廣泛使用的圖像閾值自動(dòng)選取方法[29]。設(shè)某一灰度級(jí)對(duì)應(yīng)閾值為T,類間方差計(jì)算公式為

    (3)

    式中Wb——閾值T下背景占整幅圖像的比重

    Wf——閾值T下前景占整幅圖像的比重

    1.7 全子集篩選法

    全子集篩選法是利用全子集回歸分析,對(duì)自變量不同的組合,用最小二乘法進(jìn)行建模分析,篩選最優(yōu)的變量組合。選擇最優(yōu)模型的評(píng)價(jià)標(biāo)準(zhǔn)為:①似然函數(shù)最大化。②模型未知參數(shù)最小化[30]。本文通過(guò)R編程語(yǔ)言進(jìn)行全子集篩選,利用決定系數(shù)R2和貝葉斯信息準(zhǔn)則(Bayesian information criterion,BIC)來(lái)評(píng)價(jià)篩選結(jié)果,對(duì)比分析R2越大、BIC越小的篩選結(jié)果,為最優(yōu)變量組合。

    1.8 建模策略

    建模策略分別為:未剔除土壤背景的光譜指數(shù)(策略1)、剔除土壤背景后的光譜指數(shù) (策略2)、未剔除土壤背景的光譜指數(shù)+圖像紋理特征(策略3)、 剔除土壤背景的光譜指數(shù)+圖像紋理特征(策略4)。

    1.9 模型方法和評(píng)價(jià)指標(biāo)

    極限學(xué)習(xí)機(jī)(Extreme learning machine, ELM)是一種基于最小二乘學(xué)習(xí)算法的隱層前饋網(wǎng)絡(luò)[31]。ELM收斂速度比傳統(tǒng)算法快,因?yàn)樗鼰o(wú)需迭代即可學(xué)習(xí),同時(shí),隨機(jī)隱藏節(jié)點(diǎn)保證了全局逼近能力[32]。本文ELM模型采用R語(yǔ)言elmNNRcpp包構(gòu)建。支持向量機(jī)(Support vector machines,SVM)是用于監(jiān)督學(xué)習(xí)的強(qiáng)大計(jì)算工具[33]。本研究采用非線性支持向量機(jī)模型,利用R語(yǔ)言e1071包構(gòu)建。通過(guò)決定系數(shù)R2、均方根誤差(RMSE)、標(biāo)準(zhǔn)均方根誤差(NRMSE)指標(biāo)評(píng)價(jià)模型精度。R2越接近1,RMSE越接近0,說(shuō)明模型效果越好。

    本文利用R語(yǔ)言編程構(gòu)建SVM和ELM土壤含鹽量反演模型,其中實(shí)測(cè)土壤含鹽量為因變量,光譜變量為自變量。將深度0~20 cm平均土壤含鹽量按照比例2∶1隨機(jī)劃分建模集和驗(yàn)證集,通過(guò)調(diào)整參數(shù),獲得每個(gè)條件下最佳模型。建模流程如圖2所示。

    圖2 建模流程圖Fig.2 Model building flowchart

    2 結(jié)果與分析

    2.1 土壤含鹽量

    對(duì)試驗(yàn)地總計(jì)87個(gè)采樣點(diǎn)的土壤含鹽量進(jìn)行統(tǒng)計(jì),將計(jì)算表層和深度10~20 cm實(shí)測(cè)土壤含鹽量的平均值作為深度10~20 cm的平均土壤含鹽量。將深度0~20 cm平均土壤含鹽量按照比例2∶1隨機(jī)劃分建模集和驗(yàn)證集,結(jié)果如表3所示。將各總集、建模集、驗(yàn)證集樣本點(diǎn)土壤含鹽量劃分為4個(gè)等級(jí):非鹽土(D1,SSC(土壤含鹽量)小于等于0.2%)、輕度鹽漬化(D2,SSC為(0.2%、0.5%])、重度鹽漬化(D3,SSC為(0.5%、1%])和鹽土(D4,SSC大于1.0%)[34]。非鹽土、輕度鹽漬化、重度鹽漬化和鹽漬土占比分別為34.1%、55.7%、9.1%和1.1%。含鹽量變異系數(shù)均處于中等差異(變異系數(shù)CV反映樣點(diǎn)值的離散程度,CV<0.1為弱變異性;0.11為強(qiáng)變異性)[34]。

    表3 土壤含鹽量特征統(tǒng)計(jì)分析Tab.3 Statistical analysis of soil salt content characteristics

    對(duì)各鹽漬化等級(jí)下建模集、驗(yàn)證集和總集的含鹽量進(jìn)行統(tǒng)計(jì)分析,如圖3(圖中1.5IOR表示1.5倍的四分位距)所示。從圖3可以看出,D1、D2、D3等級(jí)的建模集、驗(yàn)證集和總集的含鹽量分布、值域和均值相近,確保建模集和驗(yàn)證集數(shù)據(jù)的代表性。

    圖3 土壤含鹽量特征統(tǒng)計(jì)Fig.3 Statistical map of soil salinity characteristics

    2.2 不同植被指數(shù)分類結(jié)果

    基于近紅外和可見光波段構(gòu)建4個(gè)植被指數(shù),利用對(duì)土壤背景敏感的4個(gè)植被指數(shù)(NDVI、MSAVI、DVI、CRSI)進(jìn)行圖像分類,驗(yàn)證不同指數(shù)的圖像分類精度。為了對(duì)比4個(gè)植被指數(shù)圖像分類的精度,從87個(gè)采樣點(diǎn)中隨機(jī)選取40個(gè)樣本,進(jìn)行如圖4所示,4個(gè)植被指數(shù)的精度存在顯著差異,其中NDVI的圖像分類精度最高,40個(gè)樣本的總體精度均大于93.4%,均值為96.9%;Kappa系數(shù)在0.84~0.99之間,均值為0.92。MSAVI的分類精度略低于NDVI,總體精度和Kappa系數(shù)均值分別為95.8%和0.91;DVI和CRSI 2個(gè)植被指數(shù)的分類效果相對(duì)較差,總體精度均值小于91%,Kappa系數(shù)均值均小于0.88?;谝陨戏治?NDVI圖像分類結(jié)果表現(xiàn)最優(yōu)。

    圖4 基于不同植被指數(shù)的圖像分類結(jié)果Fig.4 Evaluation based on image classification results of different vegetation indices

    分類結(jié)果評(píng)價(jià)。分別對(duì)40個(gè)樣本的多光譜圖像進(jìn)行2.1節(jié)的操作,獲得分類結(jié)果,并在ENVI 5.3中對(duì)分類結(jié)果進(jìn)行評(píng)價(jià),得到4個(gè)植被指數(shù)的總體精度和Kappa系數(shù)。

    2.3 圖像紋理計(jì)算

    灰度共生矩陣具有豐富的特征參數(shù),能從不同的角度對(duì)紋理進(jìn)行細(xì)致刻畫[35]。由鹽漬化土壤覆蓋植被冠層的光譜特征分析可知,植被冠層光譜中紅光波段與土壤含鹽量的相關(guān)性最高,對(duì)鹽漬化的變化有最顯著的響應(yīng),故選擇波段Band3計(jì)算多光譜圖像紋理。為了更好地體現(xiàn)植株間冠層圖像的紋理差異,紋理特征參數(shù)計(jì)算的窗口選擇尺度為11像素×11像素,尺寸為71.5 cm;每個(gè)特征參數(shù)都有4個(gè)不同方向的值,取其平均值作為方向無(wú)關(guān)的特征值。將采樣點(diǎn)的多光譜圖像輸入ENVI 5.3軟件,利用二階統(tǒng)計(jì)濾波工具計(jì)算2個(gè)波段灰度共生矩陣的8個(gè)特征參數(shù)。如圖5、6所示,分別為剔除土壤背景的植被冠層圖像(E1)和原始多光譜圖像(E2)處理的可見光示意圖,均值(MEA)、方差(VAR)、均勻性(HOM)、對(duì)比度(CON)、差異(DIS)、熵(ENT)、二階矩(SEC)和相關(guān)性(COR)特征參數(shù)灰度圖。從圖中可看出,灰度共生矩陣的特征參數(shù)表現(xiàn)出圖像豐富的紋理特征,剔除土壤背景前后的多光譜圖像紋理特征存在明顯差異,統(tǒng)計(jì)灰度圖的紋理特征參數(shù)數(shù)值,對(duì)紋理特征進(jìn)行進(jìn)一步分析。

    圖5 剔除土壤背景多光譜圖像灰度共生矩陣特征參數(shù)灰度圖Fig.5 Gray scale of characteristic parameters of co-occurrence matrix of multispectral image with soil background eliminated

    圖6 原始多光譜圖像灰度共生矩陣特征參數(shù)灰度圖Fig.6 Gray level co-occurrence matrix characteristic parameters gray level

    2.4 圖像紋理特征

    圖像紋理是用于識(shí)別圖像中感興趣的對(duì)象或區(qū)域的重要特征之一[36],為了對(duì)比剔除土壤背景以及不同鹽漬化等級(jí)下多光譜圖像紋理特征的差異,對(duì)E1和E2處理的多光譜圖像Band3的紋理特征參數(shù)進(jìn)行統(tǒng)計(jì),計(jì)算2個(gè)處理下3個(gè)鹽漬化等級(jí)樣本點(diǎn)多光譜圖像的紋理特征參數(shù)的平均值,如表4所示。

    表4 多光譜圖像紋理特征參數(shù)Tab.4 Multispectral image texture feature parameters

    分別對(duì)比2個(gè)處理下3個(gè)鹽漬化等級(jí)的紋理特征參數(shù)平均值。整體上,相比剔除土壤背景前的E2,剔除土壤背景后的E1處理的紋理特征值中,SEC、CON、MEA特征值減小,DIS、VAR、ENT、COR、HOM特征值增大。其中CON、ENT、COR 3個(gè)參數(shù)的變化最顯著,CON減小表示圖像紋理溝紋變淺,圖像清晰度下降;ENT變大,表示圖像紋理增多;COR增大,表示圖像紋理變粗。

    2.5 基于全子集回歸篩選光譜自變量

    通過(guò)R語(yǔ)言全子集回歸模型算法,分別對(duì)E1和E2處理的 a (植被指數(shù))和b (植被指數(shù)和圖像紋理特征參數(shù)) 數(shù)據(jù)集進(jìn)行全子集變量篩選,篩選評(píng)價(jià)指標(biāo)選擇R2和 BIC 兩個(gè)指標(biāo)。綜合衡量R2最大和 BIC最小的變量組合,得到每個(gè)數(shù)據(jù)集中最優(yōu)的變量組,結(jié)果如表5所示。

    表5 全子集篩選最優(yōu)變量組合統(tǒng)計(jì)Tab.5 Optimal variable combination statistics screened by full subset method

    對(duì)比E1處理下2個(gè)變量種類的篩選結(jié)果,變量組a篩選得到的敏感變量組含有NDVI、SRVI、NDSI 3個(gè)植被指數(shù);變量組b的敏感變量組中包含NDVI、

    SRVI、ARVI、TVI 4個(gè)植被指數(shù),篩選掉了紋理特征參數(shù)。對(duì)比E2處理下2個(gè)變量種類的篩選結(jié)果,變量組a篩選得到的敏感變量組包含DVI、EVI、SRVI 3個(gè)植被指數(shù);變量組b敏感變量組中包含RVI、NDSI兩個(gè)植被指數(shù)和ENT、COR、SEC 3個(gè)圖像紋理特征參數(shù)。

    2.6 基于圖像紋理特征和植被指數(shù)的SSC反演模型

    通過(guò)不同植被指數(shù)分類結(jié)果分析可以得出NDVI可以達(dá)到良好的圖像分類效果。因此,利用NDVI對(duì)多光譜圖像進(jìn)行分類,并對(duì)圖像分類得到的土壤像元進(jìn)行掩膜處理,剔除土壤像元,獲得純凈度更高的植被冠層圖像。分別將E1和E2輸入ENVI 5.3軟件提取6個(gè)波段的反射率。將多光譜可見光和近紅外1波段的反射率代入植被計(jì)算公式(表1),得到相應(yīng)的植被指數(shù),涉及的植被指數(shù)有NDVI、DVI、RVI、MSAVI、ARVI、EVI、CRSI。并基于全子集回歸變量模型篩選結(jié)果,構(gòu)建基于敏感變量組的SSC反演模型。

    基于E1和E2處理下a、b兩種變量組全子集法篩選的變量結(jié)果,構(gòu)建基于敏感變量的SVM和ELM土壤含鹽量反演模型,各模型建模集和驗(yàn)證集的R2和NRMSE如表6所示。

    表6 基于不同變量組的SVM、ELM反演模型精度Tab.6 SVM and ELM inversion model based on different variable groups

    基于全子集篩選的2種敏感變量組,利用SVM和ELM兩種機(jī)器學(xué)習(xí)方法,構(gòu)建E1和E2處理共8個(gè)土壤含鹽量反演模型,驗(yàn)證集實(shí)測(cè)值和模擬值如圖7所示。

    圖7 土壤含鹽量實(shí)測(cè)值和模擬值對(duì)比Fig.7 Comparison of measured and simulated soil salt contents

    E1、E2處理下2個(gè)變量組的SVM和ELM機(jī)器學(xué)習(xí)土壤含鹽量反演模型中,ELM算法的E1-a、E1-b、E2-a、E2-b模型的R2分別為0.602、0.657、0.584、0.681,均大于對(duì)應(yīng)的SVM模型。ELM模型的標(biāo)準(zhǔn)均方根誤差分別為0.181、0.156、0.201、0.191,均小于SVM模型驗(yàn)證集的標(biāo)準(zhǔn)均方根誤差。綜上,ELM模型的表現(xiàn)更好。

    3 討論

    表7 剔除土壤背景前后紋理特征參數(shù)平均值Tab.7 Average value of image texture characteristic parameters before and after removing soil background

    本研究?jī)H基于無(wú)人機(jī)遙感多光譜反射率、圖像紋理特征與植被指數(shù)的土壤含鹽量反演,對(duì)于作物的耕種方式、作物不同生育期受土壤鹽漬化的影響,灌水量對(duì)土壤含鹽量的影響尚未考慮。本研究建立了特定情形下的土壤含鹽量反演模型,模型結(jié)果較好,但模型是否具有廣泛的適用性,不僅需要更加深入的研究,還需要綜合考慮更多因素,進(jìn)一步提高反演模型精度。

    4 結(jié)論

    猜你喜歡
    鹽漬化含鹽量植被指數(shù)
    蔬菜大棚土壤鹽漬化成因及防治措施
    含鹽量及含水率對(duì)鹽漬土凍脹規(guī)律影響試驗(yàn)研究*
    土地質(zhì)量地球化學(xué)調(diào)查成果在判定土壤鹽漬化、沙化中的應(yīng)用
    黃河三角洲鹽漬土有機(jī)氮組成及氮有效性對(duì)土壤含鹽量的響應(yīng)*
    甘肅蘇干湖濕地土壤鹽漬化、地下水位埋深及其對(duì)生態(tài)環(huán)境的影響
    什么是水的含鹽量?
    AMSR_2微波植被指數(shù)在黃河流域的適用性對(duì)比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測(cè)精度比較研究
    瑪納斯河流域土壤鹽漬化時(shí)空動(dòng)態(tài)變化
    秦陵陪葬坑土遺址安全含鹽量探究
    男女无遮挡免费网站观看| 成年女人毛片免费观看观看9 | 欧美日韩成人在线一区二区| 91大片在线观看| 中文字幕精品免费在线观看视频| 国产亚洲精品第一综合不卡| 免费观看人在逋| 国产极品粉嫩免费观看在线| 成人影院久久| 动漫黄色视频在线观看| 午夜福利乱码中文字幕| 久久免费观看电影| 少妇的丰满在线观看| 午夜日韩欧美国产| 女人高潮潮喷娇喘18禁视频| 免费女性裸体啪啪无遮挡网站| 一区二区av电影网| 99九九在线精品视频| 中文字幕制服av| 精品乱码久久久久久99久播| 日本av免费视频播放| 中文字幕最新亚洲高清| 国产精品免费视频内射| 精品视频人人做人人爽| 一级,二级,三级黄色视频| 大片电影免费在线观看免费| 777久久人妻少妇嫩草av网站| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品人人爽人人爽视色| 日韩制服骚丝袜av| 我要看黄色一级片免费的| 97在线人人人人妻| 91成年电影在线观看| 99香蕉大伊视频| 青青草视频在线视频观看| 黑人巨大精品欧美一区二区mp4| 国产精品香港三级国产av潘金莲| 中亚洲国语对白在线视频| 久久久久国产一级毛片高清牌| 色视频在线一区二区三区| 人妻久久中文字幕网| 亚洲中文字幕日韩| 老熟妇仑乱视频hdxx| 免费少妇av软件| 人人妻人人澡人人爽人人夜夜| 国产黄频视频在线观看| 久久久久国产一级毛片高清牌| 久久青草综合色| 亚洲av片天天在线观看| 国产成人a∨麻豆精品| 色婷婷久久久亚洲欧美| 黑人巨大精品欧美一区二区mp4| 国产在线一区二区三区精| 一区二区日韩欧美中文字幕| 亚洲成av片中文字幕在线观看| 亚洲人成77777在线视频| 伊人亚洲综合成人网| 国产黄频视频在线观看| 免费不卡黄色视频| 午夜成年电影在线免费观看| 一级毛片女人18水好多| 国产免费一区二区三区四区乱码| 下体分泌物呈黄色| 少妇粗大呻吟视频| 天天躁狠狠躁夜夜躁狠狠躁| 黄片小视频在线播放| 国内毛片毛片毛片毛片毛片| 亚洲一区二区三区欧美精品| 国产区一区二久久| a级毛片黄视频| h视频一区二区三区| 欧美精品人与动牲交sv欧美| 国产一区二区 视频在线| 天堂俺去俺来也www色官网| 国产伦理片在线播放av一区| 成年人黄色毛片网站| 精品一区二区三区av网在线观看 | 国产精品一区二区精品视频观看| 日韩欧美一区视频在线观看| 免费看十八禁软件| 啦啦啦在线免费观看视频4| 久久久久久久国产电影| 少妇裸体淫交视频免费看高清 | 久久久久国产精品人妻一区二区| 亚洲精品美女久久久久99蜜臀| 欧美日韩视频精品一区| 国产亚洲精品一区二区www | 亚洲精品一区蜜桃| 国产精品av久久久久免费| 成人国产av品久久久| 精品国产乱子伦一区二区三区 | 国产精品 国内视频| 色综合欧美亚洲国产小说| 欧美精品亚洲一区二区| 日本精品一区二区三区蜜桃| 免费观看a级毛片全部| 啦啦啦在线免费观看视频4| 久久久久久久大尺度免费视频| 亚洲精品一区蜜桃| 久久av网站| 久久九九热精品免费| 精品国产乱子伦一区二区三区 | 777米奇影视久久| 9热在线视频观看99| 久久久精品免费免费高清| 视频在线观看一区二区三区| 久久青草综合色| 国产深夜福利视频在线观看| 久久免费观看电影| 中文字幕高清在线视频| 精品亚洲成a人片在线观看| 天堂中文最新版在线下载| av网站在线播放免费| 黑人巨大精品欧美一区二区mp4| 99国产精品免费福利视频| 久久精品久久久久久噜噜老黄| 99九九在线精品视频| 日本欧美视频一区| 欧美日韩视频精品一区| 99久久99久久久精品蜜桃| 日本91视频免费播放| 欧美精品人与动牲交sv欧美| 久久 成人 亚洲| 男女午夜视频在线观看| 亚洲精品日韩在线中文字幕| 中文字幕高清在线视频| 午夜福利一区二区在线看| 亚洲五月婷婷丁香| 美国免费a级毛片| 欧美在线黄色| 一区二区三区精品91| 窝窝影院91人妻| av在线app专区| 纵有疾风起免费观看全集完整版| 亚洲 欧美一区二区三区| 久久精品国产a三级三级三级| 免费女性裸体啪啪无遮挡网站| 法律面前人人平等表现在哪些方面 | 日本精品一区二区三区蜜桃| 首页视频小说图片口味搜索| 一边摸一边做爽爽视频免费| 日韩视频一区二区在线观看| www.精华液| 精品国产超薄肉色丝袜足j| 色94色欧美一区二区| 一本—道久久a久久精品蜜桃钙片| 成年人午夜在线观看视频| kizo精华| 人人妻人人澡人人看| 免费在线观看视频国产中文字幕亚洲 | 99久久国产精品久久久| 久久久久国产精品人妻一区二区| 黑丝袜美女国产一区| 亚洲精品国产一区二区精华液| 黄网站色视频无遮挡免费观看| 欧美黑人欧美精品刺激| 性少妇av在线| 亚洲精品成人av观看孕妇| 日本黄色日本黄色录像| 国产激情久久老熟女| 动漫黄色视频在线观看| 满18在线观看网站| 一区二区三区四区激情视频| av欧美777| 日本av免费视频播放| 一区二区三区乱码不卡18| 亚洲精品久久成人aⅴ小说| 男人操女人黄网站| 色老头精品视频在线观看| 免费在线观看完整版高清| 亚洲欧美成人综合另类久久久| 热99国产精品久久久久久7| 久久av网站| 色婷婷av一区二区三区视频| 一区二区三区激情视频| 久久人妻熟女aⅴ| 成人影院久久| 欧美av亚洲av综合av国产av| 嫁个100分男人电影在线观看| 国产成人a∨麻豆精品| 青春草视频在线免费观看| 欧美精品一区二区免费开放| 女性被躁到高潮视频| 亚洲av美国av| 777米奇影视久久| 妹子高潮喷水视频| cao死你这个sao货| 咕卡用的链子| 十八禁网站免费在线| 国产成人欧美在线观看 | 在线十欧美十亚洲十日本专区| 宅男免费午夜| e午夜精品久久久久久久| 丝袜脚勾引网站| avwww免费| 母亲3免费完整高清在线观看| 午夜免费鲁丝| 亚洲第一av免费看| 三级毛片av免费| 亚洲欧美精品综合一区二区三区| 丝袜美腿诱惑在线| 啦啦啦中文免费视频观看日本| 天天躁狠狠躁夜夜躁狠狠躁| 黄色片一级片一级黄色片| 9热在线视频观看99| 美女高潮到喷水免费观看| 久久久国产一区二区| 老司机靠b影院| 大型av网站在线播放| 精品国内亚洲2022精品成人 | 精品福利观看| 天天添夜夜摸| 久久久久久久精品精品| 国产欧美日韩一区二区三 | 免费在线观看影片大全网站| 久久精品成人免费网站| 欧美 亚洲 国产 日韩一| 美女高潮喷水抽搐中文字幕| 免费久久久久久久精品成人欧美视频| 国产精品免费视频内射| a级片在线免费高清观看视频| 波多野结衣一区麻豆| 五月天丁香电影| 日本精品一区二区三区蜜桃| 色综合欧美亚洲国产小说| 飞空精品影院首页| 汤姆久久久久久久影院中文字幕| 日韩免费高清中文字幕av| 人妻一区二区av| 亚洲免费av在线视频| 麻豆乱淫一区二区| 国产福利在线免费观看视频| 91九色精品人成在线观看| 欧美黄色片欧美黄色片| 欧美另类一区| 国产xxxxx性猛交| 欧美人与性动交α欧美精品济南到| 欧美激情久久久久久爽电影 | 久久这里只有精品19| 777久久人妻少妇嫩草av网站| 久久女婷五月综合色啪小说| av国产精品久久久久影院| 久热这里只有精品99| 一区二区av电影网| 欧美日韩亚洲高清精品| 嫁个100分男人电影在线观看| 亚洲自偷自拍图片 自拍| a级片在线免费高清观看视频| 爱豆传媒免费全集在线观看| 亚洲 国产 在线| 亚洲精品日韩在线中文字幕| 下体分泌物呈黄色| 51午夜福利影视在线观看| 国精品久久久久久国模美| 国产精品久久久人人做人人爽| 亚洲伊人久久精品综合| 蜜桃国产av成人99| 国产黄频视频在线观看| 一级黄色大片毛片| 国产精品久久久久久人妻精品电影 | 亚洲欧美清纯卡通| 久久精品人人爽人人爽视色| 一区二区三区四区激情视频| 欧美变态另类bdsm刘玥| 午夜激情久久久久久久| 在线观看人妻少妇| 色婷婷av一区二区三区视频| 丝瓜视频免费看黄片| 亚洲五月婷婷丁香| 丝袜脚勾引网站| 美女脱内裤让男人舔精品视频| 淫妇啪啪啪对白视频 | 亚洲第一欧美日韩一区二区三区 | 国产真人三级小视频在线观看| 久久免费观看电影| 久久99热这里只频精品6学生| 一个人免费看片子| 夜夜夜夜夜久久久久| 免费在线观看影片大全网站| 香蕉国产在线看| av免费在线观看网站| 99国产综合亚洲精品| 日日夜夜操网爽| 亚洲第一青青草原| www.熟女人妻精品国产| 少妇精品久久久久久久| 色婷婷av一区二区三区视频| 久久精品国产a三级三级三级| 桃红色精品国产亚洲av| 真人做人爱边吃奶动态| 亚洲中文字幕日韩| 欧美国产精品va在线观看不卡| 日韩 欧美 亚洲 中文字幕| 高清黄色对白视频在线免费看| 精品一区在线观看国产| 黄网站色视频无遮挡免费观看| 在线精品无人区一区二区三| 成人黄色视频免费在线看| 午夜福利视频精品| 亚洲国产看品久久| 亚洲av国产av综合av卡| av国产精品久久久久影院| 日本91视频免费播放| 亚洲一区二区三区欧美精品| 中文字幕av电影在线播放| 老司机午夜十八禁免费视频| 中文字幕精品免费在线观看视频| av视频免费观看在线观看| 最新在线观看一区二区三区| 国产日韩欧美亚洲二区| 两人在一起打扑克的视频| 纯流量卡能插随身wifi吗| 青草久久国产| 亚洲中文字幕日韩| 免费人妻精品一区二区三区视频| 一区二区三区激情视频| 国产色视频综合| 亚洲五月色婷婷综合| 亚洲一码二码三码区别大吗| 一边摸一边抽搐一进一出视频| 99国产综合亚洲精品| 精品高清国产在线一区| 国产男女超爽视频在线观看| 国产免费一区二区三区四区乱码| 少妇猛男粗大的猛烈进出视频| 久久亚洲国产成人精品v| 欧美国产精品一级二级三级| 国产精品久久久人人做人人爽| 欧美日韩亚洲国产一区二区在线观看 | 中国美女看黄片| 欧美激情久久久久久爽电影 | 国产精品免费大片| 精品人妻在线不人妻| 亚洲第一欧美日韩一区二区三区 | 宅男免费午夜| 国产精品一区二区免费欧美 | 亚洲av片天天在线观看| 成在线人永久免费视频| 午夜91福利影院| 人人妻人人澡人人爽人人夜夜| 日日夜夜操网爽| 十八禁网站免费在线| 久久人人爽av亚洲精品天堂| 欧美黑人欧美精品刺激| 久久99热这里只频精品6学生| 99精品欧美一区二区三区四区| 成人18禁高潮啪啪吃奶动态图| 热99久久久久精品小说推荐| 精品第一国产精品| www.999成人在线观看| 精品第一国产精品| 亚洲成人国产一区在线观看| 亚洲av电影在线观看一区二区三区| 日韩欧美一区二区三区在线观看 | 黄片大片在线免费观看| 欧美日韩黄片免| 香蕉国产在线看| 999久久久精品免费观看国产| 亚洲人成电影观看| 性色av一级| 成人手机av| 另类精品久久| 久久精品久久久久久噜噜老黄| 纯流量卡能插随身wifi吗| 99精品久久久久人妻精品| 老司机靠b影院| 午夜精品久久久久久毛片777| 久久影院123| 久久天堂一区二区三区四区| 久久久久久久久免费视频了| 在线天堂中文资源库| 窝窝影院91人妻| 国产一区二区激情短视频 | 免费女性裸体啪啪无遮挡网站| 日韩免费高清中文字幕av| 丝袜美腿诱惑在线| 久久久久网色| 麻豆av在线久日| 亚洲欧美日韩另类电影网站| 嫁个100分男人电影在线观看| 黄网站色视频无遮挡免费观看| 超色免费av| 亚洲成人免费电影在线观看| 国产一区二区三区av在线| 在线观看www视频免费| 两性夫妻黄色片| 美女视频免费永久观看网站| 午夜老司机福利片| 人人妻,人人澡人人爽秒播| e午夜精品久久久久久久| 在线观看人妻少妇| 美女扒开内裤让男人捅视频| 他把我摸到了高潮在线观看 | 亚洲五月色婷婷综合| 桃红色精品国产亚洲av| 婷婷成人精品国产| 在线观看免费午夜福利视频| 新久久久久国产一级毛片| 国产精品亚洲av一区麻豆| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧洲日产国产| 久久中文看片网| 久久人人爽av亚洲精品天堂| 五月开心婷婷网| 国产成人精品久久二区二区免费| 免费观看av网站的网址| 亚洲国产精品成人久久小说| 超色免费av| 国产一级毛片在线| 亚洲人成电影免费在线| 精品亚洲乱码少妇综合久久| 中文字幕色久视频| h视频一区二区三区| 欧美精品av麻豆av| 啦啦啦 在线观看视频| 热re99久久国产66热| 99精国产麻豆久久婷婷| 亚洲国产精品成人久久小说| 亚洲五月婷婷丁香| 自线自在国产av| 最新在线观看一区二区三区| 久久亚洲精品不卡| 9色porny在线观看| 午夜福利免费观看在线| 欧美人与性动交α欧美软件| 色精品久久人妻99蜜桃| 国产精品久久久久久精品古装| 午夜精品久久久久久毛片777| 国产亚洲午夜精品一区二区久久| 热99久久久久精品小说推荐| 免费在线观看日本一区| 俄罗斯特黄特色一大片| 黄色视频在线播放观看不卡| 黄片大片在线免费观看| 日韩制服丝袜自拍偷拍| 各种免费的搞黄视频| 午夜福利乱码中文字幕| 午夜免费鲁丝| 精品国产一区二区三区久久久樱花| 亚洲精品国产av成人精品| 高清黄色对白视频在线免费看| 久热爱精品视频在线9| 国产色视频综合| 天天躁夜夜躁狠狠躁躁| 大码成人一级视频| a级片在线免费高清观看视频| 精品国内亚洲2022精品成人 | 天天躁狠狠躁夜夜躁狠狠躁| 丰满人妻熟妇乱又伦精品不卡| 欧美黑人欧美精品刺激| 亚洲欧美一区二区三区久久| 啪啪无遮挡十八禁网站| 欧美亚洲日本最大视频资源| 亚洲av男天堂| 美女午夜性视频免费| 搡老乐熟女国产| 黄色视频不卡| 国产精品熟女久久久久浪| 久久久久国内视频| 伊人久久大香线蕉亚洲五| 亚洲精品乱久久久久久| 亚洲精品一二三| 免费一级毛片在线播放高清视频 | 午夜激情久久久久久久| 深夜精品福利| 人妻人人澡人人爽人人| 免费看十八禁软件| 久久人人97超碰香蕉20202| 老司机午夜福利在线观看视频 | 女性生殖器流出的白浆| www日本在线高清视频| 亚洲精品日韩在线中文字幕| 亚洲全国av大片| 亚洲成国产人片在线观看| 欧美黑人精品巨大| 国产av国产精品国产| 肉色欧美久久久久久久蜜桃| 亚洲欧美日韩另类电影网站| 国产成人av教育| 免费少妇av软件| 丰满迷人的少妇在线观看| 久久久国产精品麻豆| 在线 av 中文字幕| 成人黄色视频免费在线看| videos熟女内射| 久久久精品94久久精品| 国产在线免费精品| 女性被躁到高潮视频| 成人手机av| 日本五十路高清| 日本猛色少妇xxxxx猛交久久| 亚洲久久久国产精品| 天天躁日日躁夜夜躁夜夜| 日韩人妻精品一区2区三区| 999精品在线视频| 免费不卡黄色视频| 国产不卡av网站在线观看| 美女福利国产在线| 电影成人av| 国产欧美日韩一区二区三区在线| 人妻久久中文字幕网| 91精品三级在线观看| 亚洲精品国产av成人精品| 汤姆久久久久久久影院中文字幕| 国产在线一区二区三区精| 法律面前人人平等表现在哪些方面 | 婷婷丁香在线五月| 十八禁网站网址无遮挡| 高清欧美精品videossex| 一区二区三区精品91| 日本五十路高清| av一本久久久久| 国产一区二区激情短视频 | 一区在线观看完整版| 欧美精品高潮呻吟av久久| 久久午夜综合久久蜜桃| 美女脱内裤让男人舔精品视频| 亚洲五月色婷婷综合| 午夜老司机福利片| 国产精品国产av在线观看| 免费一级毛片在线播放高清视频 | 高清在线国产一区| 精品久久久久久久毛片微露脸 | 亚洲欧美一区二区三区黑人| 蜜桃国产av成人99| 91麻豆av在线| 亚洲国产日韩一区二区| 老司机亚洲免费影院| 成年人黄色毛片网站| 欧美av亚洲av综合av国产av| 色视频在线一区二区三区| 国产伦人伦偷精品视频| 欧美亚洲日本最大视频资源| 亚洲综合色网址| 十八禁网站免费在线| 老司机亚洲免费影院| 国产熟女午夜一区二区三区| 精品一区二区三卡| 啦啦啦视频在线资源免费观看| 大型av网站在线播放| 老熟妇仑乱视频hdxx| 日本猛色少妇xxxxx猛交久久| 大香蕉久久网| 一二三四在线观看免费中文在| 久久久精品国产亚洲av高清涩受| 久久久欧美国产精品| 99九九在线精品视频| 欧美日韩成人在线一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品久久午夜乱码| 人成视频在线观看免费观看| 在线亚洲精品国产二区图片欧美| 午夜成年电影在线免费观看| 国产男女超爽视频在线观看| 天堂俺去俺来也www色官网| 亚洲精品乱久久久久久| 国产在线视频一区二区| 下体分泌物呈黄色| 亚洲国产av新网站| 亚洲,欧美精品.| 无限看片的www在线观看| 欧美日韩视频精品一区| avwww免费| 国产成人精品久久二区二区免费| 最新的欧美精品一区二区| 亚洲成人免费av在线播放| 国产日韩一区二区三区精品不卡| 久久人妻熟女aⅴ| 12—13女人毛片做爰片一| 亚洲欧洲日产国产| 黑丝袜美女国产一区| 动漫黄色视频在线观看| 老司机福利观看| 欧美日韩中文字幕国产精品一区二区三区 | 在线观看舔阴道视频| 青草久久国产| 美女中出高潮动态图| 丰满迷人的少妇在线观看| 国产免费福利视频在线观看| 狠狠精品人妻久久久久久综合| 亚洲avbb在线观看| 精品久久久久久久毛片微露脸 | 老鸭窝网址在线观看| 91字幕亚洲| 久久久久精品国产欧美久久久 | 大片免费播放器 马上看| 桃红色精品国产亚洲av| 欧美97在线视频| 午夜日韩欧美国产| 两性夫妻黄色片| 一区福利在线观看| 9191精品国产免费久久| 国产黄频视频在线观看| 久久狼人影院| 亚洲欧美一区二区三区久久| 后天国语完整版免费观看| 人妻 亚洲 视频| 亚洲欧美清纯卡通| 操出白浆在线播放| 久久久欧美国产精品| 亚洲va日本ⅴa欧美va伊人久久 | 激情视频va一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 女性被躁到高潮视频| 一区二区av电影网| 99精品欧美一区二区三区四区| 老熟妇乱子伦视频在线观看 | 国产高清视频在线播放一区 | 视频区欧美日本亚洲| 精品国产超薄肉色丝袜足j| www.精华液| 中文字幕另类日韩欧美亚洲嫩草| 三上悠亚av全集在线观看| 亚洲精品国产色婷婷电影| 狠狠狠狠99中文字幕|