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

    基于偏最小二乘與隨機(jī)森林的土壤鹽含量反演研究

    2021-05-11 05:53:51肖志云徐新宇
    安徽農(nóng)業(yè)科學(xué) 2021年8期
    關(guān)鍵詞:高光譜

    肖志云 徐新宇

    摘要 針對(duì)土默川平原地區(qū)的土壤鹽分含量提出了偏最小二乘與隨機(jī)森林相結(jié)合(RF-PLSR、PLSR-RF)對(duì)土壤鹽分含量進(jìn)行預(yù)測(cè)的回歸反演模型。該研究共采集45份土壤樣本,隨機(jī)選取35份為建模集,10份為驗(yàn)證集。試驗(yàn)首先對(duì)采集到的高光譜土壤圖像進(jìn)行分割處理提取出土壤在400~1 000 nm的原始反射光譜,其次對(duì)原始反射光譜進(jìn)行4種光譜變換(一階微分、多元散射校正的一階微分、SG平滑去噪的一階微分、對(duì)數(shù)的一階微分),并與土壤的實(shí)測(cè)鹽分量進(jìn)行相關(guān)性分析(CA),利用相關(guān)系數(shù)選取敏感波段,最后建立偏最小二乘與隨機(jī)森林結(jié)合的回歸反演模型。結(jié)果表明,與偏最小二乘回歸、隨機(jī)森林回歸單獨(dú)建模相比,2種模型結(jié)合后的預(yù)測(cè)精度有明顯的改善。光譜經(jīng)過對(duì)數(shù)的一階微分變換建立的PLSR-RF反演模型更為明顯,其建模集決定系數(shù)Rc2為0.852,均方根誤差RMSEc為0.102 g/kg,相對(duì)分析誤差RPDc為2.600,驗(yàn)證集決定系數(shù)Rv2為0.941,均方根誤差RMSEv為0.049 g/kg,相對(duì)分析誤差RPDv為4.117。

    關(guān)鍵詞 高光譜;土壤鹽含量;光譜變換;偏最小二乘回歸;隨機(jī)森林回歸

    中圖分類號(hào) TP391.4;TP79文獻(xiàn)標(biāo)識(shí)碼 A

    文章編號(hào) 0517-6611(2021)08-0010-06

    doi:10.3969/j.issn.0517-6611.2021.08.004

    開放科學(xué)(資源服務(wù))標(biāo)識(shí)碼(OSID):

    Research on Inversion of Soil Salt Content Based on Partial Least Squares Combined with Random Forest

    XIAO Zhi-yun1,2,XU Xin-yu1,2 (1.College of Electric Power,Inner Mongolia University of Technology,Huhhot,Inner Mongolia 010080;2.Inner Mongolia Key Laboratory of Mechatronic Control,Huhhot,Inner Mongolia 010051)

    Abstract Aiming at the soil salt content in the Tumochuan Plain,a regression inversion model combining partial least squares and random forest (RF-PLSR,PLSR-RF) to predict soil salt content was proposed.A total of 45 soil samples were collected in the study,35 of which were randomly selected as the modeling set and 10 of which were randomly selected as the verification set. The experiment first performed segmentation processing on the collected hyperspectral image of the soil to extract the original reflection spectrum of the soil at 400-1 000 nm,and then performed 4 kinds of spectral transformations on the original reflection spectrum (first-order differential,first-order differential of multiple scattering correction,SG smoothing Denoising first-order differential and logarithmic first-order differential). And it performed correlation analysis (CA) with the measured salt content of the soil,utilized the correlation coefficient to select the sensitive band,and finally established a regression model combining partial least squares and random forest. Compared with partial least square regression and random forest regression,the prediction accuracy of the combination of the two models was significantly improved. The PLSR-RF inversion model that established by the first-order differential transformation of the spectrum was more obvious. Its modeling set determination coefficient Rc2 was 0.852,the root mean square error RMSEc was 0.102 g/kg,and the relative analysis error RPDc was 2.600. The set determination coefficient Rv2 was 0.941,the root mean square error RMSEv was 0.049 g/kg,and the relative analysis error RPDv was 4.117.

    Key words Hyperspectral; Soil salt content; Spectral transformation;Partial least squares regression; Random forest regression

    土壤鹽堿化是目前世界面臨的最主要的環(huán)境問題之一,直接影響著農(nóng)業(yè)的可持續(xù)發(fā)展,而土默川平原地區(qū)是內(nèi)蒙古主要的糧食生產(chǎn)基地之一,由于特定的水文地質(zhì)條件、不合理的耕作和灌溉系統(tǒng),該地區(qū)出現(xiàn)了大面積的鹽堿地[1]。土地鹽漬化問題變得越來越嚴(yán)重,這嚴(yán)重影響了該地區(qū)農(nóng)牧民的收入和農(nóng)業(yè)生產(chǎn)[2]。因此研究土壤鹽含量具有重要的意義。

    隨著光譜技術(shù)的發(fā)展,越來越多的學(xué)者利用光譜技術(shù)對(duì)土壤鹽含量進(jìn)行了一定的反演研究[3-8],F(xiàn)arifteh等[9]研究發(fā)現(xiàn)土壤鹽分含量與光譜反射率之間存在顯著相關(guān)性,表明光譜之間的相似性隨土壤中鹽分濃度的增加而降低。Srivastava等[10]證明了光譜反射率對(duì)鹽度變化非常敏感。張智韜等[11]引入敏感波段組、光譜指數(shù)組、全變量組作為模型輸入變量,并建立了4種反演模型,通過對(duì)比分析得出基于光譜指數(shù)組的隨機(jī)森林鹽分反演模型在12個(gè)模型中反演效果最佳。陳俊英等[12]利用無人機(jī)搭載六波段多光譜相機(jī)和熱紅外成像儀獲取大田葵花土壤的遙感數(shù)據(jù),并同步采集區(qū)域內(nèi)不同土壤深度處的鹽分?jǐn)?shù)據(jù);通過建立反演模型得出鹽分指數(shù)和光譜指數(shù)作為變量組構(gòu)建的模型效果優(yōu)于植被指數(shù)變量組,且建立的支持向量機(jī)、反向傳播神經(jīng)網(wǎng)絡(luò)、極限學(xué)習(xí)機(jī)等機(jī)器學(xué)習(xí)鹽分反演模型比傳統(tǒng)偏最小二乘方法較優(yōu)。馬利芳等[13]利用在新疆阜康市實(shí)測(cè)的VIS-NIR光譜通過相關(guān)性分析選取特征波段建立的RF模型具有很好的預(yù)測(cè)效果。張賢龍等[14]對(duì)原始光譜進(jìn)行15種光譜變換,并構(gòu)造光譜指數(shù)對(duì)土壤鹽分含量進(jìn)行反演,得出基于倒數(shù)的對(duì)數(shù)光譜變換構(gòu)建歸一化植被指數(shù)建立的土壤鹽分反演模型精度最高。馬馳[15]對(duì)采集到的HJ-1A高光譜影像數(shù)據(jù)進(jìn)行大氣校正,并對(duì)校正之后的光譜進(jìn)行多種數(shù)學(xué)變換,結(jié)果表明一階微分和倒數(shù)的一階微分可以提高土壤反射率與實(shí)際鹽分含量的相關(guān)系數(shù)。

    綜上所述,國內(nèi)外針對(duì)土壤鹽分含量進(jìn)行的高光譜反演已經(jīng)做了一定的研究,肯定了高光譜在預(yù)測(cè)土壤鹽分含量的可行性。大多數(shù)都是采用多種光譜變換之后建立的回歸預(yù)測(cè)模型,其模型反演精度低,而利用多種光譜變換建立的2種回歸模型相結(jié)合反演土壤鹽含量的研究很少,所以該研究提出了用2種算法相結(jié)合的預(yù)測(cè)模型,首先對(duì)采集到的高光譜土壤圖像樣本進(jìn)行處理,提取出土壤的原始光譜數(shù)據(jù),其次對(duì)其進(jìn)行多種光譜變換,通過與實(shí)測(cè)土壤鹽分含量進(jìn)行相關(guān)性分析,選出敏感波段,最后建立偏最小二乘回歸、隨機(jī)森林回歸及2種模型相結(jié)合的回歸反演模型,以期為指導(dǎo)農(nóng)田種植提供參考。

    1 材料與方法

    1.1 研究區(qū)概況 察哈爾右翼中旗位于內(nèi)蒙古自治區(qū)烏蘭察布市,土壤類型為砂質(zhì)土,農(nóng)田土壤適宜馬鈴薯、玉米等多種農(nóng)作物種植,耕種制度為一年一熟。該試驗(yàn)研究地區(qū)位于察哈爾右翼中旗馬鈴薯示范基地,是內(nèi)蒙古馬鈴薯高產(chǎn)高效理論應(yīng)用技術(shù)創(chuàng)新團(tuán)隊(duì)的實(shí)驗(yàn)基地。

    1.2 土壤樣本采集與測(cè)定 樣本的采集通過野外調(diào)查取樣,根據(jù)該研究區(qū)耕種面積,共采集土樣45份,采集土壤樣本時(shí)應(yīng)均勻布局采樣點(diǎn),采樣深度為0~20 cm,將采集到的樣本土壤裝入密封袋內(nèi)帶回實(shí)驗(yàn)室,土壤樣品過篩去除雜質(zhì),每個(gè)樣本分為2份,一份利用手持式高光譜相機(jī)Specim IQ采集獲得目標(biāo)高光譜圖像數(shù)據(jù)樣本庫,數(shù)據(jù)樣本庫的采集在實(shí)驗(yàn)室用室內(nèi)拍攝系統(tǒng)(圖1)拍攝獲得土壤的高光譜圖像,高光譜相機(jī)拍攝時(shí),樣本和白板一起拍攝,白板校正可以消除環(huán)境不匹配的問題;另一份利用HM-WSYP土壤鹽分速測(cè)儀對(duì)其進(jìn)行測(cè)量,HM-WSYP土壤鹽分速測(cè)儀測(cè)量精度為±2%,故每份樣本測(cè)量3次取其平均數(shù)作為該樣本的實(shí)測(cè)鹽含量。將45個(gè)樣本隨機(jī)分為兩組,選取35個(gè)樣本用于建立反演回歸模型,10個(gè)樣本用于模型精度的驗(yàn)證(表1)。

    1.3 樣本的光譜獲取

    土壤高光譜圖像的獲取采用高光譜相機(jī)Specim IQ獲取,Specim IQ相機(jī)的波長為400~1 000 nm,光譜分辨率為3 nm,生成的高光譜圖像是3D數(shù)據(jù)結(jié)構(gòu),相機(jī)始終捕獲2D圖像,其分辨率為512×512像素,在光譜維度上記錄的光譜波段數(shù)量為204。該研究根據(jù)采集到樣本圖像中每個(gè)像素點(diǎn)的光譜曲線特征進(jìn)行圖像的分割及光譜的提取。不同的物體由于組成它們的分子結(jié)構(gòu)不同,故它們對(duì)電磁波能量的反射、吸收、透射隨波長的不同而不同。由于白板、背景和土壤的光譜反射曲線存在明顯不同,根據(jù)每個(gè)像素點(diǎn)光譜曲線的不同進(jìn)行分類,提取出土壤所有像素點(diǎn)的光譜曲線,并對(duì)土壤中的所有像素點(diǎn)的光譜反射率進(jìn)行算術(shù)平均作為此樣本的實(shí)際光譜曲線。

    1.4 光譜預(yù)處理

    高光譜圖像的光譜域噪聲在采集和傳輸?shù)雀魈幚憝h(huán)節(jié)都有可能被引入并交織呈現(xiàn)在高光譜圖像中。常用的光譜預(yù)處理方法包括SG平滑去噪、多元散射校正、倒數(shù)、對(duì)數(shù)、微分等。其中SG平滑變換可以減少雜點(diǎn),有效去除由于儀器噪聲和隨機(jī)誤差等原因?qū)е碌母哳l噪聲;多元散射校正可以消除土壤顆粒的不均勻,增強(qiáng)光譜與數(shù)據(jù)之間的相關(guān)性;倒數(shù)變換有利于顯示隱藏信息,增強(qiáng)光譜細(xì)節(jié),提高分辨率;微分變換可以消除其他背景干擾,提供比原始光譜更清晰的光譜輪廓變化和更高的分辨率[16-18]。針對(duì)所采集的土壤高光譜圖像特性,該研究采用一階微分(1D)、多元散射校正的一階微分(MSC+1D)、SG平滑去噪的一階微分(SG+1D)、對(duì)數(shù)的一階微分(LD)對(duì)原始光譜(R)進(jìn)行處理。

    1.5 敏感波段的選擇

    由于高光譜圖像的光譜波段范圍廣、波段窄、波段數(shù)量多,導(dǎo)致相鄰波段之間相關(guān)性較大,高光譜圖像中會(huì)存在較高的信息冗余,使得預(yù)測(cè)精度受到影響[19]。為了提取敏感波段,該研究對(duì)土壤含鹽量與光譜反射率的4種變換形式進(jìn)行相關(guān)性分析,并對(duì)相關(guān)系數(shù)進(jìn)行α=0.01水平的顯著性檢驗(yàn),相關(guān)系數(shù)越大表示相關(guān)性越高,此波段就越敏感。根據(jù)這一特性,選取相關(guān)系數(shù)超過顯著性檢驗(yàn)的波段為敏感波段。

    1.6 偏最小二乘與隨機(jī)森林模型的建立

    偏最小二乘回歸(PLSR)提供一種多對(duì)多線性回歸建模的方法,尤其是當(dāng)觀測(cè)數(shù)據(jù)的樣本量較少且變量之間存在多個(gè)相關(guān)性時(shí),用偏最小二乘回歸建立模型可以具有傳統(tǒng)的經(jīng)典回歸分析等方法所沒有的優(yōu)點(diǎn)[20-21],可有效簡化數(shù)據(jù)結(jié)構(gòu),進(jìn)而解決多個(gè)自變量之間高度線性相關(guān)的問題。

    隨機(jī)森林(RF)算法結(jié)構(gòu)清晰、易于解釋、運(yùn)行效率高,對(duì)于數(shù)據(jù)要求低,且具有很好的抗噪聲能力,能夠處理高維度數(shù)據(jù),不用做特征選擇,訓(xùn)練速度快,泛化能力強(qiáng),比較容易實(shí)現(xiàn)并行計(jì)算,不易出現(xiàn)過擬合問題,對(duì)于不平衡的數(shù)據(jù)來說,其還可以平衡誤差[22]。在土壤鹽含量的預(yù)測(cè)中,光譜的采集會(huì)受多種不確定因素的影響,有著高度的隨機(jī)性和非線性,隨機(jī)森林作為非參數(shù)預(yù)測(cè)模型由于不對(duì)模型的形式做假設(shè),可以在較大的函數(shù)空間內(nèi)對(duì)函數(shù)進(jìn)行擬合,因而可以更好地?cái)M合實(shí)際的模型結(jié)構(gòu),提高預(yù)測(cè)的精度。

    從統(tǒng)計(jì)學(xué)習(xí)理論的角度分析,單個(gè)預(yù)測(cè)算法往往是基于某一類假設(shè)空間。將多個(gè)預(yù)測(cè)算法相結(jié)合,可以擴(kuò)大假設(shè)空間,從而避免單個(gè)預(yù)測(cè)算法泛化性能不佳的風(fēng)險(xiǎn),達(dá)到更高的預(yù)測(cè)精度[23]。該研究提出的偏最小二乘與隨機(jī)森林的新模型(RF-PLSR、PLSR-RF)分為3個(gè)階段: ①首先利用隨機(jī)森林(或偏最小二乘)建立初始的預(yù)測(cè)模型,輸入波段xi(n)(i表示樣本,n表示波段)則可以得到隨機(jī)森林(或偏最小二乘)在訓(xùn)練樣本上的輸出預(yù)測(cè)值y^i1,用實(shí)際值yi減去預(yù)測(cè)值y^i1得到訓(xùn)練殘差yi2(yi2=yi-y^i1),然后輸入的波段xi(n)與訓(xùn)練殘差yi2進(jìn)行組合形成新的數(shù)據(jù)集;②采用偏最小二乘(或隨機(jī)森林)算法,對(duì)新的數(shù)據(jù)集進(jìn)行訓(xùn)練,輸入波段xi(n)則可以得到預(yù)測(cè)殘差y^i2;③將第1階段的預(yù)測(cè)結(jié)果y^i1與第2階段的預(yù)測(cè)結(jié)果y^i2相加即形成最終預(yù)測(cè)結(jié)果y^i(y^i=y^i1+y^i2)。

    該研究的模型即兩階段的模型相加,得到最終的預(yù)測(cè)模型。圖2為算法的原理圖。

    1.7 模型的檢驗(yàn)

    對(duì)模型精度和質(zhì)量的分析,通過計(jì)算比較模型的決定系數(shù)(R2)、均方根誤差(RMSE)和相對(duì)分析誤差(RPD)來進(jìn)行評(píng)價(jià)。R2的取值在[0,1]區(qū)間內(nèi),R2越接近于1,表示模型的擬合效果越好;R2越接近于0,表示模型的擬合效果越差。其計(jì)算公式如下:

    R2=ni=1(yi-y^i)2/ni=1(yi-y)2(1)

    式中,yi為實(shí)際測(cè)量值;y^i為預(yù)測(cè)值;y為實(shí)際值的平均;n為樣本數(shù)。

    RMSE是用來判定模型的預(yù)測(cè)能力,RMSE越小,模型精度越高,預(yù)測(cè)能力越好,其計(jì)算公式如下:

    RMSE=ni=1(yi-y^i)2/n(2)

    式中,yi為實(shí)際測(cè)量值;y^i為預(yù)測(cè)值;n為樣本數(shù)。

    RPD應(yīng)用較廣,可以一定程度上減少不同研究中預(yù)測(cè)樣本屬性值范圍差異的影響,有利于與其他相關(guān)研究對(duì)比分析。RPD越大說明所建模型越可靠,能夠用于模型分析。其計(jì)算公式如下:

    PRD=11-R2 (3)

    式中,R2為決定系數(shù)。

    該研究建模集的決定系數(shù)用Rc2表示,均方根誤差用RMSEc表示,相對(duì)分析誤差用RPDc表示;驗(yàn)證集的決定系數(shù)用Rv2表示,均方根誤差用RMSEv表示,相對(duì)分析誤差用RPDv表示。當(dāng)決定系數(shù)R2越高、RMSE越小、RPD越大時(shí),模型的反演準(zhǔn)確率越高,可靠性越強(qiáng),反之越低。

    2 結(jié)果與分析

    2.1 光譜與土壤鹽含量相關(guān)性分析

    土壤鹽含量分別與原始光譜及其4種變換形式(1D、MSC+1D、SG+1D、LD)進(jìn)行相關(guān)性分析并進(jìn)行α=0.01水平的顯著性檢驗(yàn),結(jié)果如圖3所示。由圖3可知,光譜進(jìn)行微分變換可以使敏感波段變明顯,由于土壤鹽含量與原始光譜相關(guān)性較小,沒有波段通過α=0.01的顯著性檢驗(yàn),故不適合進(jìn)行波段提取及建模估算;而土壤含鹽量與經(jīng)過1D、MSC+1D、SG+1D和LD處理后的光譜的相關(guān)性明顯提高,經(jīng)過1D光譜變換后敏感波段主要集中在467~549、811~854 nm,經(jīng)過MSC+1D光譜變換后的敏感波段主要集中在472~549、810~863 nm,經(jīng)過SG+1D光譜變換后的敏感波段主要集中在472 ~ 549、835 ~839 nm,經(jīng)過LD光譜變換后的敏感波段主要集中在467 ~ 549、729 ~ 863 nm。

    安徽農(nóng)業(yè)科學(xué) 2021年

    2.2 偏最小二乘與隨機(jī)森林模型的驗(yàn)證

    由于光譜范圍的首尾部分信噪比通常較低,數(shù)據(jù)分析過程中要考慮去除首尾部分,即400~450和900~1 000 nm的光譜。為了用采集到的土壤高光譜圖像對(duì)土壤實(shí)際含鹽量進(jìn)行預(yù)測(cè),該研究中所有回歸模型的自變量為選取相關(guān)系數(shù)由高到低的前10個(gè)敏感波段,因變量為樣本土壤的鹽分實(shí)測(cè)含量。PLSR回歸模型、RF回歸模型、RF-PLSR回歸模型和PLSR-RF回歸模型的檢驗(yàn)結(jié)果如表2所示。通過表2可以發(fā)現(xiàn),偏最小二乘(PLSR)建模集的決定系數(shù)Rc2相對(duì)偏低,均方根誤差RMSEc相對(duì)偏大,相對(duì)分析誤差RPDc沒有超過1.4,證明模型的擬合度和可靠性較差,預(yù)測(cè)精度較低。建立的RF-PLSR回歸模型和PLSR-RF回歸模型與PLSR和RF回歸模型相比,發(fā)現(xiàn)2種模型結(jié)合后建模集的決定系數(shù)Rc2有所提高,均方根誤差RMSEc有所降低,相對(duì)分析誤差RPDc都超過2.0。其中最優(yōu)模型為光譜經(jīng)過對(duì)數(shù)的一階微分變換建立的PLSR-RF反演模型,其建模集決定系數(shù)Rc2為0.852,均方根

    誤差RMSEc為0.102 g/kg,相對(duì)分析誤差RPDc為2.600,驗(yàn)證集決定系數(shù)Rv2為0.941,均方根誤差RMSEv為0.049 g/kg,相對(duì)分析誤差RPDv為4.117。

    從鹽含量與PLSR、RF、RF-PLSR、PLSR-RF回歸模型預(yù)測(cè)值與實(shí)測(cè)值的比較(圖4~7)可以看出,PLSR和RF回歸預(yù)測(cè)模型的驗(yàn)證集中有些樣本偏離1∶1線較為嚴(yán)重,而RF-PLSR、PLSR-RF回歸預(yù)測(cè)模型驗(yàn)樣本基本在1∶1線附近。說明將2個(gè)預(yù)測(cè)算法相結(jié)合可以提高預(yù)測(cè)精度。

    從研究結(jié)果來看,原始光譜的反射率經(jīng)過光譜變換后與實(shí)測(cè)含鹽量的相關(guān)性有明顯的提高,可以更好地突出敏感波段,建立的4種回歸反演模型中新提出的RF-PLSR和PLSR-RF回歸模型相比PLSR和RF回歸模型決定系數(shù)提高、均方根誤差降低、相對(duì)分析誤差變高,證明模型的擬合效果提高,精度變高,可靠性變強(qiáng),可以更好地反演土壤含鹽量。而通過LD光譜變換所建立的PLSR-RF回歸模型擬合效果最好,且預(yù)測(cè)精度最高,其建模集決定系數(shù)Rc2為0.852,均方根誤差RMSEc為0.102 g/kg,相對(duì)分析誤差RPDc為2.600,驗(yàn)證集決定系數(shù)Rv2為0.941,均方根誤差RMSEv為0.049 g/kg,相對(duì)分析誤差RPDv為4.117。由此可見,PLSR-RF回歸模型為今后預(yù)測(cè)土鹽含量提供了一個(gè)新的思路。

    通過表2中各預(yù)測(cè)方法的結(jié)果可知,該研究提出的偏最小二乘與隨機(jī)森林結(jié)合的預(yù)測(cè)方法取得了較高的預(yù)測(cè)精度,原因是采用了殘差學(xué)習(xí),如果使用偏最小二乘與隨機(jī)森林單獨(dú)建模預(yù)測(cè)則可能會(huì)導(dǎo)致訓(xùn)練數(shù)據(jù)中某些相關(guān)信息缺失,從統(tǒng)計(jì)學(xué)習(xí)理論的角度分析,采用2種模型預(yù)測(cè)可以擴(kuò)大模型的假設(shè)空間,盡可能地使假設(shè)空間包含數(shù)據(jù)之間的真實(shí)關(guān)系,而通過殘差學(xué)習(xí)對(duì)各個(gè)模型進(jìn)行組合,可以擴(kuò)大預(yù)測(cè)模型的假設(shè)空間,進(jìn)而在更大的假設(shè)空間內(nèi)對(duì)數(shù)據(jù)之間的真實(shí)關(guān)系進(jìn)行搜索,提高了模型對(duì)數(shù)據(jù)之間真實(shí)關(guān)系的逼近能力,從而提高模型的預(yù)測(cè)精度。在預(yù)測(cè)土壤鹽含量的研究中,該研究首次提出的2種模型的結(jié)合,與現(xiàn)有的高光譜反演土壤鹽含量的研究相比,2種模型的結(jié)合可以提高反演土壤鹽含量的精確度,為指導(dǎo)農(nóng)田提供參考。

    3 結(jié)論

    該研究結(jié)果表明,土默川平原土壤高光譜反射率進(jìn)行光譜變換可以更好地突出敏感波段;而采用對(duì)殘差進(jìn)行學(xué)習(xí)的2種模型的結(jié)合可以提高泛化能力和反演精度,模型的可靠性也提高。該研究提出的2種模型相結(jié)合不僅為研究土壤鹽含量提供了參考價(jià)值,同時(shí)也為今后研究土壤其他成分含量提供了一個(gè)新的思路。

    參考文獻(xiàn)

    [1]劉全明,成秋明,王學(xué),等.河套灌區(qū)土壤鹽漬化微波雷達(dá)反演[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(16):109-114.

    [2]郝遠(yuǎn)遠(yuǎn),徐旭,任東陽,等.河套灌區(qū)土壤水鹽和作物生長的HYDRUS-EPIC模型分布式模擬[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(11):110-116,315.

    [3]吳亞坤,劉廣明,蘇里坦,等.多源數(shù)據(jù)的區(qū)域土壤鹽漬化精確評(píng)估[J].光譜學(xué)與光譜分析,2018,38(11):3528-3533.

    [4]王濤,喻彩麗,姚娜,等. MLR和PLSR的沙壤土鹽分含量光譜檢測(cè)對(duì)比研究[J].干旱區(qū)地理,2018,41(6):1295-1302.

    [5]張俊華,賈萍萍,孫媛,等.基于高光譜特征的鹽漬化土壤不同土層鹽分離子含量預(yù)測(cè)[J].農(nóng)業(yè)工程學(xué)報(bào),2019,35(12):106-115.

    [6]陶培峰,王建華,李志忠,等.基于高光譜的土壤養(yǎng)分含量反演模型研究[J].地質(zhì)與資源,2020,29(1):68-75,84.

    [7]王丹陽,陳紅艷,王桂峰,等.無人機(jī)多光譜反演黃河口重度鹽漬土鹽分的研究[J].中國農(nóng)業(yè)科學(xué),2019,52(10):1698-1709.

    [8]張雅莉,塔西甫拉提·特依拜,阿爾達(dá)克·克里木,等. 基于Landsat8 OLI影像光譜的土壤鹽分估算模型研究[J].國土資源遙感 2018,30(1):87-94.

    [9]FARIFTEH J,VAN DER MEER F,CARRANZA E J M. Similarity measures for spectral discrimination of salt-affected soils[J]. International journal of remote sensing,2007,28(23):5273-5293.

    [10]SRIVASTAVA R,SETHI M,YADAV R K,et al. Visible-near infrared reflectance spectroscopy for rapid characterization of salt-affected soil in the indo-gangetic plains of Haryana,India[J]. Journal of the Indian society of remote sensing,2017,45(2):307-315.

    [11]張智韜,魏廣飛,姚志華,等.基于無人機(jī)多光譜遙感的土壤含鹽量反演模型研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2019,50(12):151-160.

    [12]陳俊英,姚志華,張智韜,等.大田葵花土壤含鹽量無人機(jī)遙感反演研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào), 2020,51(7):178-191.

    [13]馬利芳,熊黑鋼,張芳.基于野外 VIS-NIR 光譜的土壤鹽分主要離子預(yù)測(cè)[J].土壤,2020,52(1):188-194.

    [14]張賢龍,張飛,張海威,等.基于光譜變換的高光譜指數(shù)土壤鹽分反演模型優(yōu)選[J].農(nóng)業(yè)工程學(xué)報(bào),2018,34(1):110-117.

    [15]馬馳.基于HJ-1A高光譜影像的土壤鹽堿化遙感研究[J].干旱區(qū)資源與環(huán)境,2014,28(2):180-184.

    [16]石樸杰,王世東,張合兵,等.基于高光譜的復(fù)墾農(nóng)田土壤有機(jī)質(zhì)含量估測(cè)[J].土壤,2018,50(3):558-565.

    [17]于雷,洪永勝,耿雷,等.基于偏最小二乘回歸的土壤有機(jī)質(zhì)含量高光譜估算[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(14):103-109.

    [18]褚小立,袁洪福,陸婉珍.近紅外分析中光譜預(yù)處理及波長選擇方法進(jìn)展與應(yīng)用[J].化學(xué)進(jìn)展,2004,16(4):528-542.

    [19]張?zhí)栧?,李映,姜曄?深度學(xué)習(xí)在高光譜圖像分類領(lǐng)域的研究現(xiàn)狀與展望[J].自動(dòng)化學(xué)報(bào),2018,44(6):961-977.

    [20]蔣燁林,王讓會(huì),李焱,等.艾比湖流域不同土地覆蓋類型土壤養(yǎng)分高光譜反演模型研究[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào),2016,24(11):1555-1564.

    [21]高惠璇.兩個(gè)多重相關(guān)變量組的統(tǒng)計(jì)分析(3)(偏最小二乘回歸與PLS過程)[J].數(shù)理統(tǒng)計(jì)與管理,2002,21(2):58-64.

    [22]賈文超,戚蘭蘭,施凡,等.采用隨機(jī)森林改進(jìn)算法的WebShell檢測(cè)方法[J].計(jì)算機(jī)應(yīng)用研究,2018,35(5):1558-1561.

    [23]ZHOU Z H. Ensemble methods:Foundations and algorithms[M].Boca Raton,USA:CRC Press,2012.

    猜你喜歡
    高光譜
    高光譜技術(shù)在茶葉品種檢測(cè)中的應(yīng)用
    基于BP神經(jīng)網(wǎng)絡(luò)高光譜土壤氮素信息檢測(cè)研究
    科技資訊(2017年11期)2017-06-09 12:32:12
    高光譜和LiDAR數(shù)據(jù)融合在樹種識(shí)別上的應(yīng)用
    綠色科技(2017年8期)2017-05-22 14:47:06
    PCA降維和決策樹在多光譜圖像中的分類研究
    高光譜遙感技術(shù)在水環(huán)境監(jiān)測(cè)中的應(yīng)用探討
    基于實(shí)測(cè)高光譜數(shù)據(jù)的太湖湖濱帶典型植被分
    遙感圖像分類方法的研究
    基于實(shí)測(cè)光譜的紅花葉綠素含量高光譜估算模型研究
    一種基于引導(dǎo)濾波和MNF的高光譜遙感圖像分類方法
    吐絲期玉米倒伏后地面高光譜特征參數(shù)分析
    精品一区二区三区四区五区乱码 | 老司机影院成人| 国产av精品麻豆| 欧美+日韩+精品| 夜夜爽夜夜爽视频| 国产亚洲精品第一综合不卡 | 91午夜精品亚洲一区二区三区| 伦理电影大哥的女人| 久久久国产精品麻豆| 欧美日本中文国产一区发布| 97超碰精品成人国产| 亚洲人成网站在线观看播放| 亚洲欧美一区二区三区黑人 | 青青草视频在线视频观看| 美女国产高潮福利片在线看| 成人综合一区亚洲| 国产熟女欧美一区二区| 亚洲国产av新网站| 黄片播放在线免费| 另类精品久久| 精品少妇黑人巨大在线播放| 久久久亚洲精品成人影院| 成人国语在线视频| 中文字幕精品免费在线观看视频 | 人人妻人人澡人人看| 国产亚洲最大av| 国产免费视频播放在线视频| 妹子高潮喷水视频| 另类亚洲欧美激情| 国产福利在线免费观看视频| 99热网站在线观看| 一区二区三区乱码不卡18| 飞空精品影院首页| 精品一品国产午夜福利视频| 在线天堂最新版资源| 久久精品国产自在天天线| av天堂久久9| 五月天丁香电影| 涩涩av久久男人的天堂| 国产高清不卡午夜福利| 国产一区二区在线观看av| 人人澡人人妻人| 97精品久久久久久久久久精品| 午夜福利视频在线观看免费| av.在线天堂| 日韩不卡一区二区三区视频在线| 男男h啪啪无遮挡| 欧美3d第一页| 日本黄色日本黄色录像| 亚洲天堂av无毛| 国产精品一区二区在线不卡| 国产无遮挡羞羞视频在线观看| 久热这里只有精品99| 中文字幕人妻丝袜制服| 国产精品人妻久久久影院| 亚洲情色 制服丝袜| 国产亚洲欧美精品永久| 国产综合精华液| 国产亚洲精品久久久com| 麻豆乱淫一区二区| 新久久久久国产一级毛片| 80岁老熟妇乱子伦牲交| 90打野战视频偷拍视频| 天天躁夜夜躁狠狠久久av| 搡女人真爽免费视频火全软件| 久久久国产一区二区| 欧美性感艳星| 女人被躁到高潮嗷嗷叫费观| 少妇高潮的动态图| 亚洲国产精品专区欧美| 婷婷色综合www| 新久久久久国产一级毛片| 欧美最新免费一区二区三区| 日日撸夜夜添| 最后的刺客免费高清国语| 欧美成人午夜免费资源| 亚洲中文av在线| 一区在线观看完整版| 久久影院123| 人妻少妇偷人精品九色| 热99久久久久精品小说推荐| 香蕉国产在线看| 成人影院久久| 国产无遮挡羞羞视频在线观看| 蜜桃在线观看..| 一级,二级,三级黄色视频| 2018国产大陆天天弄谢| 久久精品夜色国产| 人人妻人人爽人人添夜夜欢视频| 婷婷色综合大香蕉| 国精品久久久久久国模美| 亚洲国产欧美日韩在线播放| 人人妻人人澡人人爽人人夜夜| 99香蕉大伊视频| 看免费av毛片| 熟妇人妻不卡中文字幕| 香蕉丝袜av| 十分钟在线观看高清视频www| 成人无遮挡网站| 久久精品国产鲁丝片午夜精品| 久久久久久久久久久免费av| 久久这里只有精品19| 亚洲精品av麻豆狂野| 婷婷成人精品国产| 看十八女毛片水多多多| 亚洲一级一片aⅴ在线观看| 夜夜骑夜夜射夜夜干| 新久久久久国产一级毛片| 九九在线视频观看精品| 五月玫瑰六月丁香| 精品一区二区三区视频在线| 国产精品国产av在线观看| 成人毛片60女人毛片免费| 视频在线观看一区二区三区| 亚洲伊人久久精品综合| 多毛熟女@视频| 啦啦啦中文免费视频观看日本| 人人妻人人添人人爽欧美一区卜| 一级片'在线观看视频| 午夜av观看不卡| 久久婷婷青草| 久久久久国产精品人妻一区二区| 国产福利在线免费观看视频| 18+在线观看网站| 在线亚洲精品国产二区图片欧美| 韩国av在线不卡| 美女国产高潮福利片在线看| 国产精品久久久久久精品古装| 丝袜在线中文字幕| 天天影视国产精品| 国产一区二区激情短视频 | 99久久中文字幕三级久久日本| 午夜福利影视在线免费观看| 人妻 亚洲 视频| 在线看a的网站| 高清av免费在线| 少妇人妻精品综合一区二区| 国产在线视频一区二区| 中国三级夫妇交换| 国产乱来视频区| 日本黄色日本黄色录像| 国产精品人妻久久久影院| 久久久久精品性色| 亚洲精品一二三| 一区在线观看完整版| 国产精品久久久久久av不卡| 精品人妻在线不人妻| freevideosex欧美| 人人澡人人妻人| 秋霞在线观看毛片| 国产一区有黄有色的免费视频| 男人添女人高潮全过程视频| 久久精品人人爽人人爽视色| 欧美亚洲 丝袜 人妻 在线| 亚洲综合色网址| 精品国产国语对白av| 妹子高潮喷水视频| 两个人免费观看高清视频| 欧美日韩一区二区视频在线观看视频在线| 欧美精品一区二区免费开放| 少妇的逼好多水| 最后的刺客免费高清国语| 久久久久久人人人人人| 午夜免费观看性视频| 欧美日韩综合久久久久久| 国产精品成人在线| 亚洲一级一片aⅴ在线观看| 亚洲婷婷狠狠爱综合网| 成年av动漫网址| 亚洲精品美女久久av网站| www.色视频.com| 色婷婷久久久亚洲欧美| 国产免费福利视频在线观看| 黑人猛操日本美女一级片| 在线 av 中文字幕| 国产日韩欧美在线精品| 国产成人精品在线电影| 日本与韩国留学比较| 亚洲精品乱久久久久久| 亚洲综合精品二区| 乱人伦中国视频| 熟女av电影| 日本免费在线观看一区| 99久久中文字幕三级久久日本| 欧美少妇被猛烈插入视频| 国产亚洲午夜精品一区二区久久| 久久久久人妻精品一区果冻| 国产亚洲最大av| 欧美国产精品一级二级三级| 国产欧美日韩综合在线一区二区| 超色免费av| 亚洲精品乱码久久久久久按摩| 狠狠婷婷综合久久久久久88av| 女性生殖器流出的白浆| 亚洲中文av在线| av有码第一页| 国产免费现黄频在线看| 日韩熟女老妇一区二区性免费视频| 久久久久精品人妻al黑| 欧美日韩一区二区视频在线观看视频在线| 日本wwww免费看| 精品少妇内射三级| 久久精品熟女亚洲av麻豆精品| 成人手机av| 99九九在线精品视频| 午夜免费观看性视频| 日本色播在线视频| 欧美精品亚洲一区二区| 亚洲精品美女久久久久99蜜臀 | 黄色 视频免费看| 插逼视频在线观看| 天堂8中文在线网| 久久精品国产鲁丝片午夜精品| 欧美日韩国产mv在线观看视频| 亚洲国产成人一精品久久久| 欧美激情极品国产一区二区三区 | 看免费成人av毛片| 亚洲av日韩在线播放| 十分钟在线观看高清视频www| 国产成人av激情在线播放| 国产无遮挡羞羞视频在线观看| 在线天堂中文资源库| 精品国产一区二区久久| 婷婷色综合www| 青春草视频在线免费观看| 男女免费视频国产| 各种免费的搞黄视频| 国产1区2区3区精品| 在线观看www视频免费| 80岁老熟妇乱子伦牲交| 纵有疾风起免费观看全集完整版| 国产一区二区在线观看av| 黄色毛片三级朝国网站| 中文字幕亚洲精品专区| 欧美人与善性xxx| 亚洲色图 男人天堂 中文字幕 | 亚洲欧美一区二区三区国产| 两个人免费观看高清视频| 国产福利在线免费观看视频| 十分钟在线观看高清视频www| 午夜老司机福利剧场| 极品少妇高潮喷水抽搐| 成人午夜精彩视频在线观看| 一级毛片我不卡| 国产69精品久久久久777片| 成人综合一区亚洲| 99视频精品全部免费 在线| 亚洲美女视频黄频| 在线观看免费高清a一片| 人妻系列 视频| 午夜免费男女啪啪视频观看| 亚洲精品乱码久久久久久按摩| 国国产精品蜜臀av免费| 国产成人aa在线观看| av视频免费观看在线观看| 亚洲精品乱码久久久久久按摩| 欧美国产精品一级二级三级| 欧美精品亚洲一区二区| 国产成人午夜福利电影在线观看| 少妇人妻 视频| 亚洲欧美中文字幕日韩二区| 大陆偷拍与自拍| 丰满迷人的少妇在线观看| 亚洲国产精品国产精品| 美女主播在线视频| 久久精品久久久久久久性| 日本-黄色视频高清免费观看| 亚洲情色 制服丝袜| 欧美精品高潮呻吟av久久| 国产精品一国产av| 国产黄色视频一区二区在线观看| 精品国产露脸久久av麻豆| 国产综合精华液| 精品国产一区二区三区四区第35| 国产麻豆69| 国产乱来视频区| 精品熟女少妇av免费看| 91精品国产国语对白视频| 少妇的逼水好多| 桃花免费在线播放| 欧美成人午夜免费资源| 免费观看a级毛片全部| 国产精品熟女久久久久浪| 国产在线免费精品| 天美传媒精品一区二区| 18禁国产床啪视频网站| 黑人高潮一二区| 国产成人欧美| 亚洲精品自拍成人| 夜夜骑夜夜射夜夜干| 精品人妻熟女毛片av久久网站| 日韩av免费高清视频| 久久精品国产a三级三级三级| av在线观看视频网站免费| av不卡在线播放| 国产精品.久久久| 亚洲欧洲精品一区二区精品久久久 | 国产精品熟女久久久久浪| 侵犯人妻中文字幕一二三四区| 免费av不卡在线播放| 美女xxoo啪啪120秒动态图| 夜夜爽夜夜爽视频| 亚洲成国产人片在线观看| 久久久久久人人人人人| 国产精品一国产av| 久久久精品区二区三区| 亚洲一区二区三区欧美精品| 国产 精品1| 国产av精品麻豆| 亚洲 欧美一区二区三区| 久久婷婷青草| 国产成人精品婷婷| 国产日韩欧美在线精品| 国产精品国产三级国产av玫瑰| 女的被弄到高潮叫床怎么办| 亚洲色图 男人天堂 中文字幕 | 99精国产麻豆久久婷婷| 丝袜美足系列| 成人免费观看视频高清| 啦啦啦中文免费视频观看日本| 国产成人午夜福利电影在线观看| 国产精品欧美亚洲77777| 久久久a久久爽久久v久久| 日本午夜av视频| 国产精品一区www在线观看| 精品国产一区二区三区久久久樱花| 免费观看性生交大片5| 赤兔流量卡办理| av又黄又爽大尺度在线免费看| 久久精品国产a三级三级三级| 国产精品一国产av| 国产无遮挡羞羞视频在线观看| 一级片免费观看大全| 一区二区三区乱码不卡18| 亚洲中文av在线| 综合色丁香网| 18+在线观看网站| 九色成人免费人妻av| 国产亚洲精品第一综合不卡 | 在线观看一区二区三区激情| 亚洲一级一片aⅴ在线观看| 免费黄网站久久成人精品| 超色免费av| 日韩 亚洲 欧美在线| 在线天堂最新版资源| 国产成人a∨麻豆精品| 久久av网站| 亚洲精品国产色婷婷电影| 久久精品熟女亚洲av麻豆精品| 18禁裸乳无遮挡动漫免费视频| 精品99又大又爽又粗少妇毛片| 午夜av观看不卡| 人妻系列 视频| a级毛片在线看网站| 老司机亚洲免费影院| 国产男女超爽视频在线观看| 99九九在线精品视频| 日韩熟女老妇一区二区性免费视频| 日韩成人av中文字幕在线观看| 国产精品人妻久久久久久| 久久人人爽人人爽人人片va| 午夜影院在线不卡| 国产成人精品婷婷| 日韩精品免费视频一区二区三区 | 欧美亚洲日本最大视频资源| 男人操女人黄网站| 欧美 日韩 精品 国产| 制服诱惑二区| 热re99久久精品国产66热6| 国产黄色免费在线视频| 亚洲精品成人av观看孕妇| 精品人妻偷拍中文字幕| 免费黄色在线免费观看| 亚洲精品乱码久久久久久按摩| 亚洲欧美中文字幕日韩二区| 欧美日韩一区二区视频在线观看视频在线| 天堂8中文在线网| 涩涩av久久男人的天堂| 欧美亚洲 丝袜 人妻 在线| av卡一久久| 久久狼人影院| 亚洲av成人精品一二三区| 人妻少妇偷人精品九色| 性高湖久久久久久久久免费观看| 精品一区在线观看国产| 黄色视频在线播放观看不卡| 91成人精品电影| 国产成人精品福利久久| 老司机影院毛片| 日本与韩国留学比较| 久久午夜综合久久蜜桃| 乱码一卡2卡4卡精品| 亚洲av国产av综合av卡| 一本久久精品| www.av在线官网国产| 大码成人一级视频| 亚洲国产成人一精品久久久| freevideosex欧美| 欧美国产精品va在线观看不卡| 在线精品无人区一区二区三| 啦啦啦在线观看免费高清www| 欧美丝袜亚洲另类| av网站免费在线观看视频| 一区二区三区四区激情视频| 国产成人午夜福利电影在线观看| 99久久中文字幕三级久久日本| 久久久精品区二区三区| 黄片无遮挡物在线观看| 黄色 视频免费看| av福利片在线| av黄色大香蕉| 爱豆传媒免费全集在线观看| 成年女人在线观看亚洲视频| 久久久久人妻精品一区果冻| 国产在线视频一区二区| 18禁在线无遮挡免费观看视频| 国产高清不卡午夜福利| 热99国产精品久久久久久7| 欧美xxⅹ黑人| 欧美xxxx性猛交bbbb| freevideosex欧美| 飞空精品影院首页| 日本vs欧美在线观看视频| 精品亚洲成国产av| 国产精品免费大片| 国产精品人妻久久久久久| 9色porny在线观看| 女人久久www免费人成看片| 好男人视频免费观看在线| 精品国产露脸久久av麻豆| 亚洲av免费高清在线观看| 人成视频在线观看免费观看| 麻豆精品久久久久久蜜桃| 黄片无遮挡物在线观看| 九九爱精品视频在线观看| 少妇人妻 视频| videos熟女内射| 18禁裸乳无遮挡动漫免费视频| 日韩视频在线欧美| av一本久久久久| 深夜精品福利| 亚洲精品av麻豆狂野| 免费高清在线观看日韩| 大香蕉97超碰在线| 母亲3免费完整高清在线观看 | 亚洲国产欧美日韩在线播放| 日韩制服骚丝袜av| 国产精品麻豆人妻色哟哟久久| 自拍欧美九色日韩亚洲蝌蚪91| 午夜免费观看性视频| 成人18禁高潮啪啪吃奶动态图| 99国产精品免费福利视频| 久久韩国三级中文字幕| 国产av码专区亚洲av| 亚洲av免费高清在线观看| 天天躁夜夜躁狠狠久久av| 国产又爽黄色视频| 午夜av观看不卡| 黄色 视频免费看| 草草在线视频免费看| 又黄又爽又刺激的免费视频.| 免费久久久久久久精品成人欧美视频 | 国产成人精品一,二区| 久久久久视频综合| 国产成人午夜福利电影在线观看| 一区二区三区四区激情视频| 九九在线视频观看精品| 国产黄色视频一区二区在线观看| 男人操女人黄网站| 丁香六月天网| 国产精品免费大片| 在线精品无人区一区二区三| 精品国产国语对白av| 国产在线免费精品| 久久精品久久精品一区二区三区| 黄色一级大片看看| 亚洲综合精品二区| 久久免费观看电影| 亚洲人与动物交配视频| 99香蕉大伊视频| 婷婷成人精品国产| 91国产中文字幕| 免费看av在线观看网站| www.色视频.com| 黑人巨大精品欧美一区二区蜜桃 | 免费黄网站久久成人精品| 成人国语在线视频| 18+在线观看网站| 亚洲 欧美一区二区三区| 人体艺术视频欧美日本| 亚洲国产精品成人久久小说| 国产精品无大码| 免费黄网站久久成人精品| 黑人高潮一二区| 午夜免费男女啪啪视频观看| 亚洲,一卡二卡三卡| 日本免费在线观看一区| 欧美成人午夜精品| 精品国产一区二区久久| 久久婷婷青草| 91在线精品国自产拍蜜月| 插逼视频在线观看| 侵犯人妻中文字幕一二三四区| 日本午夜av视频| 在线观看人妻少妇| 国产高清国产精品国产三级| 日韩熟女老妇一区二区性免费视频| 日本欧美国产在线视频| 国产精品成人在线| 极品少妇高潮喷水抽搐| www.熟女人妻精品国产 | 春色校园在线视频观看| 丝袜人妻中文字幕| 久久精品国产亚洲av天美| 久久人人爽人人爽人人片va| 99国产精品免费福利视频| 性色avwww在线观看| a级毛片在线看网站| 18禁动态无遮挡网站| 日韩三级伦理在线观看| kizo精华| av卡一久久| 欧美精品人与动牲交sv欧美| 欧美亚洲日本最大视频资源| 免费不卡的大黄色大毛片视频在线观看| 国产又色又爽无遮挡免| 在线观看免费高清a一片| 午夜福利网站1000一区二区三区| 久久精品久久久久久久性| 亚洲av男天堂| 夜夜爽夜夜爽视频| 啦啦啦在线观看免费高清www| 国产日韩欧美在线精品| 欧美日韩精品成人综合77777| 国产精品国产三级国产av玫瑰| 在线观看人妻少妇| 亚洲国产欧美在线一区| 久久精品国产鲁丝片午夜精品| 精品一区二区免费观看| 国产亚洲av片在线观看秒播厂| 日韩制服骚丝袜av| 国产精品蜜桃在线观看| 亚洲美女黄色视频免费看| 最近2019中文字幕mv第一页| 99国产综合亚洲精品| 在线观看美女被高潮喷水网站| 中国美白少妇内射xxxbb| 久久久久久久精品精品| a级片在线免费高清观看视频| 搡老乐熟女国产| 免费黄频网站在线观看国产| 久久这里只有精品19| av线在线观看网站| 日韩精品有码人妻一区| 黑人巨大精品欧美一区二区蜜桃 | 国产成人免费无遮挡视频| 国产在线免费精品| 丝袜美足系列| 久久久久久久久久久免费av| 久久国产亚洲av麻豆专区| 亚洲久久久国产精品| 欧美性感艳星| 久久久久精品人妻al黑| av在线播放精品| 看十八女毛片水多多多| 性色avwww在线观看| 男女免费视频国产| 久久久久网色| 边亲边吃奶的免费视频| 又黄又粗又硬又大视频| 日本免费在线观看一区| 色婷婷久久久亚洲欧美| 在线看a的网站| 国产成人aa在线观看| 亚洲精品日韩在线中文字幕| 欧美少妇被猛烈插入视频| 国产精品久久久av美女十八| 国产亚洲欧美精品永久| 久久国产精品男人的天堂亚洲 | 久久久久久人妻| av片东京热男人的天堂| 国产精品久久久av美女十八| 亚洲欧洲日产国产| 亚洲国产精品999| 日日爽夜夜爽网站| 亚洲经典国产精华液单| 最后的刺客免费高清国语| 欧美xxⅹ黑人| 狂野欧美激情性bbbbbb| 国产精品国产三级专区第一集| 亚洲精品456在线播放app| 亚洲美女视频黄频| 亚洲国产毛片av蜜桃av| 有码 亚洲区| 亚洲,欧美精品.| 高清欧美精品videossex| 69精品国产乱码久久久| 人妻 亚洲 视频| 成人无遮挡网站| 如日韩欧美国产精品一区二区三区| 欧美97在线视频| 久久韩国三级中文字幕| 最新的欧美精品一区二区| 少妇精品久久久久久久| 国精品久久久久久国模美| 1024视频免费在线观看| 国产色爽女视频免费观看| 色视频在线一区二区三区| 精品福利永久在线观看| 97在线人人人人妻| 久久久久久久精品精品| 日韩人妻精品一区2区三区| 乱码一卡2卡4卡精品| 91在线精品国自产拍蜜月| 1024视频免费在线观看|