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

    夏玉米葉面積指數(shù)遙感反演研究

    2016-10-27 02:04:11李彥榮杜靈通
    農(nóng)業(yè)機械學(xué)報 2016年9期
    關(guān)鍵詞:全生育期植被指數(shù)夏玉米

    劉 珺 龐 鑫 李彥榮 杜靈通

    (1.太原理工大學(xué)安全技術(shù)與工程博士后科研流動站, 太原 030024;2.山西省交通科學(xué)研究院博士后科研工作站, 太原 030006; 3.太原理工大學(xué)礦業(yè)工程學(xué)院, 太原 030024;4.寧夏大學(xué)西北退化生態(tài)系統(tǒng)恢復(fù)與重建教育部重點實驗室, 銀川 750021)

    ?

    夏玉米葉面積指數(shù)遙感反演研究

    劉珺1,2龐鑫3李彥榮3杜靈通4

    (1.太原理工大學(xué)安全技術(shù)與工程博士后科研流動站, 太原 030024;2.山西省交通科學(xué)研究院博士后科研工作站, 太原 030006; 3.太原理工大學(xué)礦業(yè)工程學(xué)院, 太原 030024;4.寧夏大學(xué)西北退化生態(tài)系統(tǒng)恢復(fù)與重建教育部重點實驗室, 銀川 750021)

    利用LAI-2000的觀測數(shù)據(jù)與基于HJ衛(wèi)星遙感數(shù)據(jù)生成的植被指數(shù),結(jié)合3種常用的回歸模型,構(gòu)造了夏玉米分別按全生育期、不同生育階段和閾值分段的葉面積指數(shù)(Leaf area index, LAI)反演模型;獲取了3種模式下LAI的最優(yōu)反演模型;在驗證和評價各模型可靠性之后,生成了夏玉米在營養(yǎng)生長期、抽雄期和生殖生長期的LAI分布圖;并將基于HJ影像反演得到的LAIHJ與MODIS LAI產(chǎn)品(MOD15A2)LAIM進行了對比。研究發(fā)現(xiàn),與各種通用植被指數(shù)相比,針對HJ CCD數(shù)據(jù)構(gòu)造的環(huán)境植被指數(shù)HJVI與LAI的相關(guān)性在3種反演模式中均為最佳。HJVI與全生育期LAI的相關(guān)性達到0.875,在不同生育階段與LAI的相關(guān)性也高于其他植被指數(shù)(營養(yǎng)生長期線性模型最佳,決定系數(shù)為0.769;抽雄期對數(shù)模型最佳,決定系數(shù)為0.783;生殖生長期指數(shù)模型最佳,決定系數(shù)為0.703)。普適性植被指數(shù)中,OSAVI適用于夏玉米生長前中期的LAI反演,NDVI適用于夏玉米生長后期的LAI反演。在夏玉米全生育期內(nèi),各植被指數(shù)與LAI的相關(guān)性整體較高,但最優(yōu)回歸模型出現(xiàn)在按不同生育階段反演的模式中。LAI小于3時EVI為精度最佳指數(shù)(決定系數(shù)為0.358),LAI不小于3時OSAVI為精度最佳指數(shù)(決定系數(shù)為0.515)。在夏玉米3個生育階段,LAIM與LAIHJ的相關(guān)性分別達到0.732、0.761、0.661。HJ遙感數(shù)據(jù)具有較強的LAI反演能力,其高時間和高空間分辨率的特征可以使其代替?zhèn)鹘y(tǒng)的中分辨率遙感數(shù)據(jù)而成為農(nóng)業(yè)遙感研究的重要數(shù)據(jù)源。

    夏玉米; 葉面積指數(shù); 反演; 植被指數(shù); HJ衛(wèi)星

    引言

    葉面積指數(shù)(Leaf area index, LAI)被定義為地面面積上綠色葉子的單面面積之和,是描述植被生長狀況的重要參數(shù)之一[1],在植被碳、水循環(huán)中起著重要作用[2]。作物是生態(tài)系統(tǒng)的重要組成部分,LAI是反映作物產(chǎn)量、光合作用能力、生長健康狀況的重要指標[3]。

    傳統(tǒng)的葉面積指數(shù)一般采用直接測量法,直接測量法包括收獲測量法、落葉收集法和異速生長測量法,但由于該方法耗時、耗力、工作量大,且對農(nóng)作物具有一定的破壞性,因此其應(yīng)用受到一定的限制。

    利用遙感技術(shù)反演LAI已成為近年來植被遙感領(lǐng)域的研究熱點。利用遙感反演LAI的方法主要為經(jīng)驗?zāi)P头ê臀锢砟P头?。物理模型法包括幾何光學(xué)模型法、輻射傳輸模型法和混合模型法。它是目前反演全球1 km尺度LAI產(chǎn)品的主流算法。雖然物理模型法具有較高的LAI反演精度,但其受傳感器覆蓋范圍和獲取時相的制約較大[4],且物理模型法需要的輸入?yún)?shù)不但難以獲取,而且本身帶有一定誤差,同時其反演精度受模型模擬精度、查表精度、實測數(shù)據(jù)采集精度等因素的影響[5]。而經(jīng)驗?zāi)P头ㄊ抢脤崪yLAI與植被指數(shù)VI(Vegetation index)建立可靠的回歸模型,從而間接求得大范圍內(nèi)LAI的方法。由于植被指數(shù)本身具有克服干擾因素的特點,因此經(jīng)驗?zāi)P头ê喕斯趯觾?nèi)光子的復(fù)雜傳輸過程,減少了參數(shù)的不確定性,提高了LAI估算精度,是廣域范圍內(nèi)低成本估算LAI的有效方法,已成為進行小區(qū)域LAI反演的常用方法[6]。

    由于遙感數(shù)據(jù)捕獲的空間變異特性有利于調(diào)整區(qū)域內(nèi)作物和土壤特性,因此高時間、高空間分辨率的遙感數(shù)據(jù)在經(jīng)驗?zāi)P椭芯哂胁豢商娲奈恢肹7]。地物在可見光的反射可用來控制土壤背景的干擾,而在近紅外波段的反射則可探測大范圍內(nèi)的動態(tài)波動[8]。因此,由可見光和近紅外波段組合而成的植被指數(shù)被廣泛地應(yīng)用于LAI反演研究中。其中歸一化植被指數(shù)(Normalized difference vegetation index, NDVI)是目前應(yīng)用最廣的植被指數(shù),但是它受土壤背景影響較大,且容易達到飽和;為了彌補這個缺點,提出了土壤調(diào)節(jié)植被指數(shù)(Soil adjusted vegetation index, SAVI)、優(yōu)化的土壤調(diào)節(jié)植被指數(shù)(Optimized soil adjusted vegetation index, OSAVI)和修正的土壤調(diào)節(jié)植被指數(shù)(Modified soil adjusted vegetation index, MSAVI)[9-11];而為了提高LAI較大時的敏感度、減少大氣干擾,增強型植被指數(shù)(Enhanced vegetation index, EVI)也被較多地應(yīng)用于LAI反演[12-13]。在LAI反演研究中,學(xué)者們還應(yīng)用了其他植被指數(shù)。如譚昌偉等[14]通過測定不同時期的夏玉米冠層光譜反射率與LAI,探索了比值植被指數(shù)(Ratio vegetation index, RVI)等10種植被指數(shù)與LAI的相關(guān)性,利用R810/R560與LAI之間的極其顯著相關(guān)性,高精度地反演了夏玉米的葉面積指數(shù)。趙虎等[15]利用NDVI、EVI、RVI、SAVI等反演了冬小麥葉面積指數(shù),證明了葉面積指數(shù)和不同植被指數(shù)之間有較強的相關(guān)關(guān)系。YANG等[16]利用NDVI、EVI反演夏玉米LAI,通過比較發(fā)現(xiàn)EVI的反演精度高于NDVI。賈玉秋等[17]在分析GF-1和Landsat-8譜段特征的基礎(chǔ)上,利用紅光、藍光和綠光的不同組合構(gòu)成的歸一化植被指數(shù)(NDVI、GNDVI、BNDVI、GRNDVI、GBNDVI、BRNDVI、PNDVI)分別反演了玉米的葉面積指數(shù),說明采用高空間分辨率的數(shù)據(jù)可以提高LAI估算精度。趙娟等[18]采用分生育期反演的辦法,分別基于寬波段植被指數(shù)(OSAVI、NDVI、EVI、TVI、NDGI、RVI)和窄波段植被指數(shù)(NDVI705、mRVI705、mNDVI705、VOGI)反演了不同生育期的冬小麥LAI,有效改善了冬小麥LAI的反演精度。蘇偉等[19-20]用激光雷達數(shù)據(jù)精確研究了農(nóng)作物冠層參數(shù)與葉面積指數(shù)之間的關(guān)系??梢钥闯?,已有研究集中于2方面:利用窄波段植被指數(shù)探求高精度反演LAI的模型;通過細化和改進反演方法,利用寬波段植被指數(shù)估算LAI。這些研究都存在一些問題:與寬波段植被指數(shù)比較,構(gòu)建窄波段植被指數(shù)對數(shù)據(jù)來源要求較高,不易獲得,且數(shù)據(jù)處理工作繁重;大多數(shù)研究使用的是已經(jīng)過驗證的TM、SPOT、MODIS等國外的遙感數(shù)據(jù),較少使用國產(chǎn)數(shù)據(jù)。

    HJ衛(wèi)星是我國自主研發(fā)的資源環(huán)境衛(wèi)星,已被應(yīng)用于各個領(lǐng)域,其中農(nóng)業(yè)領(lǐng)域的研究也不少,但是針對LAI反演的精細研究尚不多見。為了驗證HJ衛(wèi)星在LAI估算方面的能力,同時得到更高的夏玉米LAI反演精度,本文以太原市為例,利用我國自主研發(fā)的HJ衛(wèi)星CCD遙感數(shù)據(jù)生成的寬波段植被指數(shù),分別從全生育期、不同生育階段和LAI閾值分段3個角度反演夏玉米LAI。通過探索不同模式下各植被指數(shù)與實測夏玉米LAI之間的相關(guān)關(guān)系,建立回歸模型,進一步探尋在遙感數(shù)據(jù)空間分辨率一定的前提下,有效提高夏玉米葉面積指數(shù)反演精度的途徑和方式,并驗證我國HJ衛(wèi)星在精細農(nóng)業(yè)方面的應(yīng)用能力。

    1 研究區(qū)概況

    圖1 太原市HJ假彩色(432)合成圖及采樣點Fig.1 HJ-1 composited image for band 4,3,2 of Taiyuan city and sample points

    太原地區(qū)(111°30′E~113°09′ E, 37° 27′N~38°25′ N,圖1)下轄太原市市轄區(qū)、古交市、清徐縣、婁煩縣、陽曲縣等地區(qū)。區(qū)域內(nèi)糧食作物面積為7.33萬hm2,其中夏玉米播種面積約5萬hm2,占全省夏玉米面積的68%。區(qū)域內(nèi)種植的夏玉米主要為春播中晚熟型,其生育期為4月中旬—9月下旬。具體為:4月中旬播種,5月初出苗,5月中旬—6月中旬為夏玉米拔節(jié)期,之后隨著溫度逐漸升高,夏玉米生長加速,到7月中旬夏玉米進入抽雄期,8月中旬—9月上旬為灌漿期,9月中下旬時,夏玉米逐漸成熟,到9月底、10月初時收獲。全生育期103~133 d,為一年一熟的純作夏玉米。

    2 數(shù)據(jù)與方法

    2.1數(shù)據(jù)獲取

    2.1.1夏玉米LAI觀測

    利用美國LI-COR公司生產(chǎn)的LAI-2000型植被冠層儀進行野外夏玉米LAI數(shù)據(jù)采集和測量。為了全面獲取研究區(qū)夏玉米的葉面積指數(shù),根據(jù)夏玉米種植面積和密度,在研究區(qū)所轄行政區(qū)(太原市市轄區(qū)、陽曲縣、婁煩縣、古交市和清徐縣)內(nèi)共選取了18個夏玉米地塊的162個點進行葉面積指數(shù)測量(圖1),樣點LAI最大值為4.20,最小值為0.90,均值為2.45,標準差為0.67。樣方的選擇是根據(jù)夏玉米長勢和植株密度確定的,參考HJ數(shù)據(jù)分辨率,將地塊大小設(shè)置為30 m×30 m,并利用Mobil Mapper 6 GPS觀測和記錄了定位點信息。測量日期為2015年7月中旬—9月下旬,共進行了3次。

    2.1.2植被指數(shù)的選擇

    選用NDVI、EVI、差值植被指數(shù)(Difference vegetation index, DVI)、RVI、OSAVI、3波段植被指數(shù)(Triangular vegetation index, TVI)6種應(yīng)用廣泛的典型植被指數(shù)[14-21]和專門針對HJ CCD數(shù)據(jù)構(gòu)造的環(huán)境植被指數(shù)(Huanjing vegetation index, HJVI)[22]估算夏玉米LAI。各植被指數(shù)計算公式及來源如表1所示。

    表1 各植被指數(shù)計算公式Tab.1 Overview of spectral vegetation indices

    注:ρNIR、ρR、ρG、ρB分別為植被在近紅外、紅光、綠光、藍光波段的反射率。

    2.1.3HJ CCD遙感數(shù)據(jù)

    為了將遙感影像與地面實測數(shù)據(jù)進行匹配,從中國資源衛(wèi)星應(yīng)用中心(http:∥www.cresda.com/n16/index.html)獲取了與地面LAI觀測日期同步的HJ-1 CCD影像,分別得到7月29日、8月22日、9月19日太原市HJ影像。利用采集的地面控制點對影像進行了幾何校正、輻射校正、大氣校正、研究區(qū)裁剪等預(yù)處理,其中幾何校正誤差控制在0.5個像元內(nèi),輻射校正根據(jù)資源衛(wèi)星應(yīng)用中心提供的定標參數(shù)進行,大氣校正采用6S模型進行,利用研究區(qū)矢量圖對大氣校正后的影像進行裁剪,最終得到研究區(qū)各時相反射率影像。

    根據(jù)表1中的公式,生成各生育期的NDVI、EVI、DVI、RVI、OSAVI、TVI影像。為了減小因幾何校正造成的定位誤差,夏玉米樣點植被指數(shù)取GPS采樣點周圍3像元×3像元的平均值。

    2.1.4MODIS LAI產(chǎn)品

    從美國地質(zhì)勘探局(http:∥lpdaac.usgs.gov/)下載研究區(qū)MODIS葉面積指數(shù)和光合有效輻射產(chǎn)品(MOD15A2),軌道號為h26v05,產(chǎn)品空間分辨率為1 km,時間分辨率為8 d。利用ENVI軟件對下載的數(shù)據(jù)進行了轉(zhuǎn)投影、裁剪等處理,得到研究區(qū)葉面積指數(shù)影像。

    2.2LAI估算模型的構(gòu)建與比較

    將實測LAI按照2∶1的比例分開,選取其中的98個數(shù)據(jù)用于模型反演。通過獲取相應(yīng)點的HJ植被指數(shù),構(gòu)建夏玉米LAI反演模型。分別按照全生育期、不同生育階段和LAI閾值分段3種方法進行建模。

    2.2.1全生育期建模

    在多種的回歸模型中,當樣本量較大時,線性模型被認為是最具說服力的回歸模型,因此,研究構(gòu)建了全生育期內(nèi)VI與LAI的線性回歸模型。將98個樣點的NDVI、EVI、DVI、RVI、OSAVI、TVI與夏玉米全生育期的LAI值一一對應(yīng),進行回歸分析,得到整個生育期內(nèi)各植被指數(shù)與LAI的線性相關(guān)性及其回歸模型(表2),并以決定系數(shù)R2為指標評價各模型的可靠性。

    從表2中可以看出,各線性模型中,HJVI與LAI的相關(guān)性最高;其余植被指數(shù)與LAI的相關(guān)性由大到小分別為RVI、EVI、NDVI、OSAVI、DVI、TVI。

    2.2.2不同生育階段建模

    雖然全生育期反演得到了較好的LAI反演精度,但細節(jié)性不足。其反演結(jié)果具有較好的整體性,但不足以為監(jiān)測作物長勢、估算作物產(chǎn)量提供詳細信息。為了提高夏玉米LAI反演精度,使研究結(jié)果可以得到更好地推廣應(yīng)用,按照不同生育階段對夏玉米LAI進行反演。

    表2 全生育期各植被指數(shù)與葉面積指數(shù)線性相關(guān)性Tab.2 Correlativity between VIs and LAI in whole growth period

    作物的生育過程可以分為營養(yǎng)生長期、生殖生長期和抽雄期3個階段。營養(yǎng)生長期是綠色植物根、莖、葉等營養(yǎng)器官的生長階段。生殖生長期是植物花、果實、種子等生殖器官形成的時期。對夏玉米而言,抽雄期為營養(yǎng)生長和生殖生長的轉(zhuǎn)折點,也是作物葉面積指數(shù)達到峰值的時期。

    根據(jù)營養(yǎng)生長期、抽雄期和生殖生長期3個階段對夏玉米LAI進行反演建模。在實地調(diào)查研究區(qū)夏玉米生長情況的基礎(chǔ)上,3個生育期分別選擇7月29日、8月22日和9月14日。分別將各生育期植被指數(shù)與對應(yīng)點實測LAI值進行回歸分析,得到不同生育階段各植被指數(shù)與LAI的線性、對數(shù)、指數(shù)3種回歸模型,如圖2所示。

    圖2 不同生育階段各植被指數(shù)與葉面積指數(shù)3種模型相關(guān)性Fig.2 Correlativity between VIs and LAI in different growth stages of summer maize

    從圖中可以看出,3個生育期內(nèi),相關(guān)性最高的回歸模型分別為HJVI與LAI在營養(yǎng)生長期的線性模型(R2=0.769)、抽雄期的對數(shù)模型(R2=0.783)、生殖生長期的指數(shù)模型(R2=0.703)。除HJVI之外,營養(yǎng)生長期內(nèi),普適性植被指數(shù)中,OSAVI與LAI的線性模型相關(guān)性最高,R2=0.724,其次為NDVI與LAI的線性模型、EVI與LAI的指數(shù)模型、TVI與LAI的線性模型、RVI與LAI的線性模型、DVI與LAI的線性模型,其中DVI與LAI的指數(shù)關(guān)系最差(R2=0.486)。抽雄期內(nèi),OSAVI與LAI的相關(guān)性最高,且二者的相關(guān)性明顯優(yōu)于其他植被指數(shù)與LAI的相關(guān)性。在OSAVI與LAI的3種模型中,指數(shù)模型相關(guān)性最佳,R2為0.749,其余依次為NDVI與LAI的指數(shù)模型、DVI與LAI的對數(shù)模型、RVI與LAI的對數(shù)模型、EVI與LAI的的指數(shù)模型、TVI與LAI的線性模型。其中DVI、RVI與LAI的相關(guān)性差異不明顯,TVI與LAI的指數(shù)關(guān)系最差,R2=0.352。生殖生長期內(nèi),NDVI與LAI的線性模型相關(guān)性最佳,R2=0.700,其他回歸模型依次為RVI與LAI的線性模型(R2=0.605)、OSAVI的線性模型(R2=0.558)、TVI與LAI的線性模型(R2=0.543)、DVI與LAI的對數(shù)模型(R2=0.495)、EVI的線性模型(R2=0.493)。生殖生長期內(nèi),EVI與LAI的對數(shù)關(guān)系最差,R2=0.454。

    2.2.3按葉面積指數(shù)閾值分段建模

    以LAI等于3為分界值,按照LAI小于3和LAI不小于3分段對夏玉米葉面積指數(shù)反演進行建模,選擇線性模型表征二者的相關(guān)性,得到各植被指數(shù)與LAI的相關(guān)性(表3)。

    表3 按LAI閾值分段時植被指數(shù)與葉面積指數(shù)的 線性相關(guān)性Tab.3 Linear correlativity between VIs and LAI according to LAI values

    從表中可以看出,在LAI小于3和LAI不小于3兩種情況下,HJVI與LAI的線性回歸模型的相關(guān)性遠高于其他植被指數(shù)。除HJVI外,LAI小于3時,EVI與LAI的線性相關(guān)性最佳,其他植被指數(shù)與LAI的線性關(guān)系從大到小依次為NDVI、RVI、OSAVI、DVI、TVI。LAI不小于3時,除HJVI外,OSAVI與LAI線性相關(guān)最強,其余由大到小分別為DVI、TVI、EVI和NDVI,RVI與LAI相關(guān)性最弱,R2僅為0.096,說明RVI不適合LAI較大時的反演研究。

    3 討論

    3.1LAI建模方式比較

    3種反演模式中,HJVI具有其他植被指數(shù)不可比擬的優(yōu)勢,這與HJVI的構(gòu)造有關(guān),表明針對特定遙感數(shù)據(jù)構(gòu)建的植被指數(shù)可以充分發(fā)揮該數(shù)據(jù)的特點。

    各普適性植被指數(shù)與夏玉米LAI在全生育期內(nèi)的線性相關(guān)性整體較高,R2最小為0.450。這也是大多數(shù)研究植被LAI反演的學(xué)者們普遍采用全生育期反演的LAI原因。

    分生育階段反演的模式下:在營養(yǎng)生長期和抽雄期,OSAVI與LAI的相關(guān)性最強,并且在這2個時期內(nèi),OSAVI與LAI的線性相關(guān)性高于其他各植被指數(shù)與LAI在全生育期內(nèi)的相關(guān)性,說明OSAVI適合用于夏玉米生長初期和中期的LAI反演。緊隨其后的是NDVI,NDVI與LAI在營養(yǎng)生長期的線性相關(guān)性也大于二者在全生育期內(nèi)的相關(guān)程度,而在抽雄期則相反。這與NDVI在作物生長旺盛時期容易飽和溢出而造成其值與LAI值對應(yīng)不準確有關(guān)。而在生殖生長期內(nèi)葉片開始逐步失水、發(fā)黃,NDVI值不易造成飽和。因此,在此期間內(nèi)NDVI與LAI的相關(guān)性顯著高于其他植被指數(shù)與LAI的相關(guān)程度,同時也高于NDVI與LAI在全生育期內(nèi)的相關(guān)性。綜上所述,NDVI更適用于作物生長前期和后期的LAI反演。

    LAI為3是植被葉面積指數(shù)高低的分界點,LAI小于3時對應(yīng)的各植被指數(shù)較小,LAI不小于3時對應(yīng)的各植被指數(shù)相應(yīng)也較大。LAI小于3時,EVI與LAI的相關(guān)性較強,這是由于EVI克服了土壤背景和大氣干擾的影響,其值更真實地反映了植被的生長狀況。此外,NDVI與LAI的相關(guān)性也證實了該植被指數(shù)在LAI較小時的反演能力。LAI不小于3時,OSAVI與LAI的相關(guān)性高于二者在LAI小于3時的相關(guān)性,而且優(yōu)于其他植被指數(shù)與LAI的相關(guān)性;其次,由于植被在近紅外波段和紅波段的反射率差異較大,使得DVI、TVI與LAI的對應(yīng)關(guān)系也較好。除OSAVI、DVI和TVI之外,其余植被指數(shù)與LAI不小于3時的相關(guān)性均小于其與LAI小于3時的相關(guān)性。

    綜合比較3種反演方式,不難發(fā)現(xiàn),按不同生育階段反演的LAI最優(yōu)精度高于其他2種方式的最優(yōu)反演精度;全生育期的各植被指數(shù)對LAI的反演精度總體較高,但最佳精度不及分生育階段的值;EVI、NDVI、RVI適合于LAI小于3時的反演,OSAVI、DVI、TVI則適合于LAI不小于3時的反演;3種模式下,HJVI均具有最優(yōu)的LAI反演精度,且明顯優(yōu)于其他普適性植被指數(shù)。

    3.2模型驗證

    為了檢驗LAI各反演模型的可靠性,研究利用各生育階段HJVI與LAI的反演模型得到LAI的最佳反演值,并將實測LAI作為真值,利用建模樣本以外的其余64個樣點LAI實測值與基于HJVI反演的LAI值進行回歸分析,以決定系數(shù)R2和均方根誤差(RMSE)為指標對其進行精度評價,并利用Z-檢驗對2組數(shù)據(jù)進行統(tǒng)計分析。得到各生育期LAI反演值(LAIR)與真實值(LAIT)之間的相關(guān)關(guān)系(圖3)。

    從圖3中可以看出,利用HJVI反演的營養(yǎng)生長期、抽雄期和生殖生長期的LAIR與LAIT之間的線性決定系數(shù)R2分別為0.695、0.681、0.591,RMSE分別為0.271、0.217、0.213(圖3a、3b、3c);利用HJVI反演的全生育期的LAIR與LAI真值之間的線性決定系數(shù)R2為0.837,RMSE為0.409(圖3d)。

    圖3 不同生育期真實葉面積指數(shù)與基于HJVI反演的葉面積指數(shù)相關(guān)性分析Fig.3 Correlation of measured LAI and retrieval vaule based on HJVI at different stages

    綜合分析圖3可以得出,雖然各生育期相關(guān)性不同,但其RMSE最大值為0.271,說明基于HJVI的夏玉米不同生育期LAI反演模型能真實反映夏玉米的長勢及其變化情況,其反演結(jié)果較好。

    3.3LAI成圖

    在驗證各反演模型可靠性的基礎(chǔ)上,利用最佳反演模型分別生成太原市夏玉米的LAI空間分布圖(圖4)。

    圖4 研究區(qū)夏玉米不同生育期葉面積指數(shù)空間分布Fig.4 LAI spatial distribution maps of summer maize at different stages in Taiyuan

    圖5 不同生育期LAIHJ與LAIM相關(guān)性分析Fig.5 Correlation of LAIHJ and LAIM

    從圖4中可以看出,營養(yǎng)生長階段夏玉米LAI取值在0.5~3.0之間,且夏玉米LAI值大多集中于0.5~2.0之間,均值為1.90。LAI介于1.5~3.0之間的較大值主要集中在太原市東北方向,這是由于該區(qū)夏玉米播種時間較其他區(qū)域早,因此其植株生長較其他區(qū)域茂盛。在抽雄期,夏玉米LAI值主要集中于1.5~3.5之間,均值為2.74。LAI大于3的值位于太原市南部的清徐縣,研究區(qū)西部的婁煩縣和古交市夏玉米LAI總體較小。夏玉米進入成熟期后,玉米果實逐步形成,LAI有所降低,主要集中在1.0~2.5之間,均值為1.79。其中LAI較大值主要集中于南部的清徐縣。

    3.4MODIS LAI產(chǎn)品驗證

    為了進一步驗證反演模型的準確性,研究將反演得到的夏玉米3個典型生長期LAI影像重采樣到1 km,并隨機選取了143個點,將基于HJ遙感影像的夏玉米營養(yǎng)生長期、抽雄期和生殖生長期的葉面積指數(shù)(LAIHJ)與對應(yīng)時相的MODIS葉面積指數(shù)產(chǎn)品(MOD15A2,LAIM)進行回歸分析,得到夏玉米營養(yǎng)生長期、抽雄期和生殖生長期LAIHJ與LAIM的相關(guān)關(guān)系(圖5)。

    而從圖5中可以看出,夏玉米3個典型生長期LAIM與LAIHJ呈現(xiàn)出較好的一致性,其中抽雄期二者的相關(guān)性最好,R2達到0.761(圖5b);營養(yǎng)生長期次之,R2為0.732(圖5a);生殖生長期二者相關(guān)性也較好,R2為0.661(圖5c)。對143個樣點的均值進行統(tǒng)計發(fā)現(xiàn),營養(yǎng)生長期LAIHJ均值為2.02,LAIM均值為1.75;抽雄期LAIHJ均值為3.06,LAIM均值為2.75;生殖生長期LAIHJ均值為2.79,LAIM均值為2.63。3個夏玉米生長期LAIM均不同程度地低于LAIHJ,這與MODIS空間分辨率低、像元混合嚴重的情況有關(guān)。這種高估或低估的現(xiàn)象也被其他學(xué)者證實過[3,29]。

    4 結(jié)論

    (1)以夏玉米為研究對象,利用實測LAI數(shù)據(jù)和基于HJ衛(wèi)星生成的HJVI、NDVI、EVI、DVI、RVI、OSAVI、TVI數(shù)據(jù),構(gòu)造了全生育期、不同生育階段、LAI閾值分段3種模式下的夏玉米LAI反演模型。3種反演模式下,HJVI均具有最優(yōu)的LAI估算精度。

    (2)3種夏玉米LAI反演模式中,全生育期內(nèi)各植被指數(shù)與LAI的相關(guān)性整體較高,但其細節(jié)性不強,無法精確反映作物在不同生育階段的生長狀況。而最優(yōu)反演精度出現(xiàn)在按不同生育期反演的模式中。說明分不同生育階段反演能有效提高LAI的估算精度,有利于更準確地監(jiān)測作物生長狀況,估算其產(chǎn)量。

    (3)分生育期反演模式中各普適性植被指數(shù)與LAI的相關(guān)性差異較大。營養(yǎng)生長期,OSAVI與LAI線性相關(guān)性最佳;抽雄期,OSAVI與LAI指數(shù)模型相關(guān)性最佳;生殖生長期,NDVI與LAI線性相關(guān)性最佳。說明在夏玉米生長前期OSAVI具有較強的LAI反演能力。

    (4)按LAI閾值分段反演的精度是3種反演模式中最差的,各普適性植被指數(shù)中反演精度最高為0.515(LAI不小于3,OSAVI與LAI),最低精度僅為0.096(LAI不小于3,RVI與LAI)。這種反演模式不利于利用LAI研究作物長勢和產(chǎn)量。

    (5)LAIM與LAIHJ雖然有差異,但在對精度要求不高的情況下,MODIS LAI產(chǎn)品亦可反映作物的生長狀況。

    (6)HJ遙感數(shù)據(jù)具有較強的LAI反演能力,其高時間和高空間分辨率的特征可以使其代替?zhèn)鹘y(tǒng)的中分辨率遙感數(shù)據(jù)而成為農(nóng)業(yè)遙感研究的重要數(shù)據(jù)源。

    1MYNENI R B, HOFFMAN S, KNYAZIKHIN Y, et al. Global products of vegetation leaf area and fraction absorbed PAR from year one of MODIS data[J]. Remote Sensing of Environment, 2002, 83(1): 214-231.

    2劉洋,劉榮高,陳鏡明,等.葉面積指數(shù)遙感反演研究進展與展望[J].地球信息科學(xué)學(xué)報, 2013,15(5): 734-743.

    LIU Yang, LIU Ronggao,CHEN Jingming, et al. Current status and perspectives of leaf area index retrieval from optical remore sensing data[J].Journalof Geo-Information Science, 2013,15(5): 734-743. (in Chinese)

    3LE Maire G, MARSDEN C, VERHOEF W, et al. Leaf area index estimation with MODIS reflectance time series and model inversion during full rotations of eucalyptus plantations[J]. Remote Sensing of Environment,2011,115(2):586-599.

    4LIANG S,ZHAO X,LIU S, et al. A long-term global land surface satellite (GLASS) data-set for environmental studies[J].International Journal of Digital Earth,2013,6(1):5-33.

    5劉玉琴,孟慶巖,沙晉明,等.基于HJ-1A CCD1數(shù)據(jù)的臺灣相思樹葉面積指數(shù)反演[J].地球信息科學(xué)學(xué)報,2014,16(4):645-652.

    LIU Yuqin,MENG Qingyan,SHA Jinming,et al. Comparison of different methods for retrieving acacia rachii leaf area index based on HJ-1A CCD1 imagery[J]. Journal of Geo-Information Science,2014,16(4):645-652. (in Chinese)

    6趙靜,李靜,柳欽火,等.聯(lián)合HJ-1 /CCD 和Landsat8 /OLI數(shù)據(jù)反演黑河中游葉面積指數(shù)[J].遙感學(xué)報,2015,19(5) :733-749.

    ZHAO Jing, LI Jing, LIU Qinhuo, et al. Leaf area index inversion combining with HJ-1/CCD and Landsat 8 /oli data in the middle reach of the Heihe river basin[J]. Journal of Remote Sensing, 2015,19(5):733-749. (in Chinese)

    7GUéRIF M, DUKE C L. Adjustment procedures of a crop model to the site specific characteristics of soil and crop using remote sensing data assimilation[J]. Agriculture, Ecosystems & Environment, 2000,81(1):57-69.

    8PINTY B, LAVERGNE T, WIDLOWSKI J L,et al. On the need to observe vegetation canopies in the near-infrared to estimate visible light absorption[J]. Remote Sensing of Environment, 2009,113(1):10-23.

    9HUETE A R. A soil-adjusted vegetation index(SAVI)[J]. Remote Sensing of Environment, 1988,25(3):295-309.

    10RONDEAUX G, STEVEN M, BARET F. Optimization of soil-adjusted vegetation indices[J]. Remote Sensing of Environment, 1996,55(2):95-107.

    11QI J, CHEHBOUNI A, HUETE A R, et al. A modified soil adjusted vegetation index [J]. Remote Sensing of Environment,1994,48(2): 119-126.

    12HUETE A, DIDAN K, MIURA T, et al. Overview of the radiometric and biophysical performance of the MODIS vegetation indices[J]. Remote Sensing of Environment, 2002,83(1): 195-213.

    13JIANG Z, HUETE A, DIDAN K, et al. Development of a two-band enhanced vegetation index without a blue band[J]. Remote Sensing of Environment, 2008,112(10):3833-3845.

    14譚昌偉,黃義德,黃文江,等.夏玉米葉面積指數(shù)的高光譜遙感植被指數(shù)法研究[J].安徽農(nóng)業(yè)大學(xué)學(xué)報,2004,31(4):392-397.

    TAN Changwei,HUANG Yide, HUANG Wenjiang, et al.Study on colony leaf area index of summer maize by remote sensing vegetation indexes method[J]. Journal of Anhui Agricultural University, 2004,31(4):392-397. (in Chinese)

    15趙虎,裴志遠,馬尚杰,等. 利用HJ-1-A/B CCD2數(shù)據(jù)反演冬小麥葉面積指數(shù)[J].農(nóng)業(yè)工程學(xué)報,2012,28(10): 172-176.

    ZHAO Hu, PEI Zhiyuan, MA Shangjie, et al. Retrieving leaf area index of winter wheat using HJ-1-A/B CCD2 data[J].Transactions of the CSAE, 2012,28(10): 172-176. (in Chinese)

    16YANG Fei, SUN Jiulin, FANG Honliang,et al. Comparison of different methods for corn LAI estimation over northeastern China[J].International Journal of Applied Earth Observation and Geoinformation, 2012,18:462-471.

    17賈玉秋,李冰,程永政,等. 基于GF-1 與 Landsat-8多光譜遙感影像的玉米LAI反演比較[J].農(nóng)業(yè)工程學(xué)報,2015, 31(9): 173-179.

    JIA Yuqiu, LI Bing, CHENG Yongzheng, et al.Comparison between GF-1 image and Landsat-8 images in monitoring maize lai[J].Transactions of the CSAE, 2015, 31(9): 173-179. (in Chinese)

    18趙娟,黃文江,張耀鴻,等. 冬小麥不同生育時期葉面積指數(shù)反演方法[J].光譜學(xué)與光譜分析,2013,33(9): 2546-2552.

    ZHAO Juan, HUANG Wenjiang, ZHANG Yaohong, et al. Inversion of leaf area index different growth stages in winter wheat[J].Spectroscopy and Spectral Analysis, 2013,33(9): 2546-2552. (in Chinese)

    19蘇偉,郭皓,趙冬玲,等.基于優(yōu)化PROSAIL葉傾角分布函數(shù)的玉米LAI反演方法[J].農(nóng)業(yè)機械學(xué)報,2016,47(3):234-241.

    SU Wei,GUO Hao,ZHAO Dongling,et al.Leaf area index retrivel for maize canopy using optimized leaf angle disribution function of PROSAIL model[J].Transactions of the Chinese Society for Agricultural Machinery,2016,47(3):234-241. (in Chinese)

    20蘇偉, 展郡鴿, 張明政,等. 基于機載LiDAR數(shù)據(jù)的農(nóng)作物葉面積指數(shù)估算方法研究[J]. 農(nóng)業(yè)機械學(xué)報, 2016,47 (3):272-277.

    SU Wei, ZHAN Junge, ZHANG Mingzheng, et al. Crop leaf area index estimation based on airborne LiDAR data[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(3):272-277. (in Chinese)

    21ZHENG G, MOSKAL L M. Retrieving leaf area index (LAI) using remote sensing: theories, methods and sensors[J]. Sensors, 2009, 9(4): 2719-2745.

    22BARET F, GUYOT G. Potentials and limits of vegetation indices for LAI and APAR assessment[J].Remote Sensing of Environment, 1991, 35(2): 161-173.

    23張瀛,孟慶巖,武佳麗,等.基于環(huán)境星CCD數(shù)據(jù)的環(huán)境植被指數(shù)及葉面積指數(shù)反演研究[J].光譜學(xué)與光譜分析,2011,31(10):2789-2793.

    ZHANG Ying, MENG Qingyan,WU Jiali,et al. Study of environmental vegetation index based on environment satellite CCD data and LAI inersion[J]. Spectroscopy and Spectral Analysis,2011,31(10):2789-2793. (in Chinese)

    24ROUSE J W, HAAS R H, SCHELL J A, et al. Monitoring vegetation systems in the great plains with erts[J]. NASA Special Publication, 1974,351:309.

    25LIU H Q, HUETE A. A feedback based modification of the NDVI to minimize canopy background and atmospheric noise[J]. IEEE Transactions on Geoscience and Remote Sensing, 1995(2): 457-465.

    26JORDAN C F. Derivation of leaf-area index from quality of light on the forest floor[J]. Ecology, 1969,50(4): 663-666.

    27PEARSON R L, MILLER L D. Remote mapping of standing crop biomass for estimation of the productivity of the shortgrass prairie, pawnee national grasslands, colorado[C]∥Proceedings of the 8th International Symposium on Remote Sensing of Environment, 1972:1357-1381.

    28BROGE N H, LEBLANC E. Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density [J].Remote Sensing of Environment, 2001,76(2): 156-172.

    29FENSHOLT Rasmus, SANDHOLT Inge, RASMUSSEN M S, et al. Evaluation of MODIS LAI,fAPAR and the reation between fAPAR and NDVI in a semi-arid environment using insitume asurements[J]. Remote Sensing of Environment, 2004,91(3-4): 490-507.

    Inversion Study on Leaf Area Index of Summer Maize Using Remote Sensing

    Liu Jun1,2Pang Xin3Li Yanrong3Du Lingtong4

    (1.PostdoctoralScientificResearchStationofSafetyTechnologyandEngineering,TaiyuanUniversityofTechnology,Taiyuan030024,China2.PostdoctoralResearchWorkstationofShanxiTransportationResearchInstitute,Taiyuan030006,China3.CollegeofMiningEngineering,TaiyuanUniversityofTechnology,Taiyuan030024,China4.KeyLaboratoryforRestorationandReconstructionofDegradedEcosysteminNorth-WesternChina,MinistryofEducation,NingxiaUniversity,Yinchuan750021,China)

    The observation data of LAI-2000 and vegetation index was generated by satellite remote sensing data of HJ, combining three kinds of commonly used regression model. LAI (Leaf area index) inversion model was constructed according to growth period, growth stage and threshold boundaries of summer maize, respectively. The optimal LAI inversion model was acquired based on the above three modes. The summer corn LAI scatter grams of the vegetative growth period, the tasseling stage as well as the reproductive stage were generated after verification and evaluation of the model reliability. The productions of MODIS LAI (MOD15A2) were verified by LAIHJbased on the inversion of model HJ image. According to the survey, except HJVI, during the whole growth period of summer maize, a linear model of RVI with LAI was regarded as the best fitting model (R2= 0.662); during the vegetative growth period, a linear model of OSAVI with LAI was regarded as the best fitting model (R2=0.724); at the tasseling stage, index model of OSAVI with LAI was regarded as the best fitting model(R2=0.749); at the reproductive stage, a linear model of NDVI with LAI was regarded as the best fitting model(R2=0.700). The correlation of HJVI and LAI at the growth period achieved to 0.875, and the correlation at different growth stages with LAI is higher than the other vegetation indexes (during the vegetative period,R2=0.769; at the tasseling stage,R2=0.783; at the reproductive stage,R2=0.703). EVI is the best index when LAI is less than 3 (R2=0.358), while OSAVI is the best when LAI is more than 3(R2=0.515). During the three reproductive periods,R2of LAIMand LAIHJis 0.732, 0.761 and 0.661. Conclusions were drawn: the inversion method of LAI at different stages is optimal. HJVI shows obvious advantage for LAI inversion ability. The production of MODIS LAI could be used for crop monitoring in special situation. The study not only broadens the mode of inversion LAI using vegetation index, but also confirms the importance of HJ data in agricultural field.

    summer maize; leaf area index; inversion; vegetation index; HJ satellite

    10.6041/j.issn.1000-1298.2016.09.042

    2015-11-16

    2016-03-25

    國家自然科學(xué)基金項目(41201438)、山西省青年研究基金項目(2014021032-1)和太原理工大學(xué)?;痦椖?2013Z016)

    劉珺(1981— ),女,講師,博士,主要從事資源環(huán)境遙感研究,E-mail: 8886355@163.com

    TP79

    A

    1000-1298(2016)09-0309-09

    猜你喜歡
    全生育期植被指數(shù)夏玉米
    小麥收割之后 如何種植夏玉米才能高產(chǎn)
    AMSR_2微波植被指數(shù)在黃河流域的適用性對比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測精度比較研究
    夏玉米高產(chǎn)的關(guān)鍵栽培技術(shù)措施
    湖南省2016年審定通過的水稻新品種(下)
    中國稻米(2017年2期)2017-04-28 08:00:06
    2016年靖遠縣春玉米全生育期農(nóng)業(yè)氣象條件分析
    湖南省2015年審定通過的水稻新品種(4)
    中國稻米(2016年2期)2016-06-29 09:53:29
    主要植被指數(shù)在生態(tài)環(huán)評中的作用
    西藏科技(2015年1期)2015-09-26 12:09:29
    基于MODIS數(shù)據(jù)的植被指數(shù)與植被覆蓋度關(guān)系研究
    黃驊市2013年度冬小麥全生育期農(nóng)業(yè)氣象條件分析
    日韩大尺度精品在线看网址| 女同久久另类99精品国产91| 成人亚洲欧美一区二区av| 国产高清视频在线播放一区| 五月玫瑰六月丁香| 欧美zozozo另类| 少妇人妻精品综合一区二区 | 久久99热6这里只有精品| 国产又黄又爽又无遮挡在线| h日本视频在线播放| 成人特级黄色片久久久久久久| 婷婷亚洲欧美| 热99re8久久精品国产| 91精品国产九色| 女人被狂操c到高潮| 精品午夜福利视频在线观看一区| 又爽又黄无遮挡网站| 国产视频内射| 日日干狠狠操夜夜爽| 亚洲美女黄片视频| 最新中文字幕久久久久| 国产私拍福利视频在线观看| 精品久久久久久久久亚洲| 97碰自拍视频| 深夜a级毛片| 日韩精品有码人妻一区| 91麻豆精品激情在线观看国产| 国产精品一区二区三区四区久久| 国产一区二区亚洲精品在线观看| 国国产精品蜜臀av免费| 最近视频中文字幕2019在线8| av在线播放精品| 亚洲精品国产av成人精品 | 国产伦在线观看视频一区| 日韩一本色道免费dvd| 国产黄a三级三级三级人| 亚洲av免费高清在线观看| 天堂影院成人在线观看| 成年女人看的毛片在线观看| 天堂动漫精品| 国产真实乱freesex| 国产探花在线观看一区二区| 一区二区三区高清视频在线| 国产一区二区亚洲精品在线观看| aaaaa片日本免费| 精品久久久久久久人妻蜜臀av| 亚洲成人久久性| 精品乱码久久久久久99久播| 在现免费观看毛片| 亚洲在线观看片| 三级男女做爰猛烈吃奶摸视频| 国产精品久久久久久久电影| 色噜噜av男人的天堂激情| 国产精品国产高清国产av| 日韩欧美免费精品| 亚洲自偷自拍三级| 久久婷婷人人爽人人干人人爱| 一区福利在线观看| 亚洲国产精品sss在线观看| 深夜a级毛片| 免费电影在线观看免费观看| 有码 亚洲区| 久久久久久国产a免费观看| 亚洲av一区综合| 亚洲精品影视一区二区三区av| 亚洲最大成人手机在线| 欧美xxxx性猛交bbbb| 99久国产av精品| 精品99又大又爽又粗少妇毛片| 亚洲在线观看片| 免费观看在线日韩| 色哟哟·www| 乱系列少妇在线播放| 国产精品电影一区二区三区| 波多野结衣高清无吗| 久久久久久大精品| 99热6这里只有精品| 久久久久久久久中文| 国产三级中文精品| 日韩成人av中文字幕在线观看 | 嫩草影视91久久| 久久久久久久午夜电影| 女人被狂操c到高潮| 日韩,欧美,国产一区二区三区 | 久久精品人妻少妇| 我要搜黄色片| av黄色大香蕉| 日本熟妇午夜| 国产成人91sexporn| 在线观看午夜福利视频| 国产欧美日韩精品一区二区| 自拍偷自拍亚洲精品老妇| 99九九线精品视频在线观看视频| 色av中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 老熟妇乱子伦视频在线观看| 国产私拍福利视频在线观看| 狂野欧美白嫩少妇大欣赏| 成人精品一区二区免费| 亚洲成人中文字幕在线播放| 六月丁香七月| 日韩大尺度精品在线看网址| 午夜福利18| 校园人妻丝袜中文字幕| 内地一区二区视频在线| 不卡一级毛片| 免费大片18禁| 久久精品国产鲁丝片午夜精品| 国产老妇女一区| 亚洲av免费在线观看| 禁无遮挡网站| av在线亚洲专区| 美女免费视频网站| 18禁裸乳无遮挡免费网站照片| 亚洲在线自拍视频| 亚洲乱码一区二区免费版| 日本爱情动作片www.在线观看 | 天堂√8在线中文| 一进一出抽搐动态| 成人高潮视频无遮挡免费网站| 最新中文字幕久久久久| 简卡轻食公司| 特级一级黄色大片| 国产精品一及| 欧美色视频一区免费| 欧美高清成人免费视频www| 无遮挡黄片免费观看| 亚洲av熟女| 久久久久精品国产欧美久久久| 精品熟女少妇av免费看| 国产精品三级大全| 国产视频内射| 白带黄色成豆腐渣| 久久人人爽人人爽人人片va| 性欧美人与动物交配| 白带黄色成豆腐渣| 亚洲av成人av| 免费看日本二区| 人人妻人人澡人人爽人人夜夜 | 在线播放国产精品三级| 日本免费一区二区三区高清不卡| 黄色视频,在线免费观看| 欧美另类亚洲清纯唯美| 亚洲成人久久性| 99九九线精品视频在线观看视频| 久久精品国产鲁丝片午夜精品| 色哟哟哟哟哟哟| 我要看日韩黄色一级片| 欧美另类亚洲清纯唯美| 日韩一本色道免费dvd| 搡老妇女老女人老熟妇| 精品一区二区三区av网在线观看| 亚洲真实伦在线观看| 日日摸夜夜添夜夜爱| 精品久久久久久久久久免费视频| 91av网一区二区| 国产成人影院久久av| 啦啦啦韩国在线观看视频| 日产精品乱码卡一卡2卡三| av视频在线观看入口| 午夜久久久久精精品| 国产精品综合久久久久久久免费| or卡值多少钱| 桃色一区二区三区在线观看| 三级男女做爰猛烈吃奶摸视频| 天堂动漫精品| 欧美另类亚洲清纯唯美| 亚洲精品亚洲一区二区| 欧美xxxx黑人xx丫x性爽| 日日摸夜夜添夜夜添小说| 欧美日韩在线观看h| 18禁在线无遮挡免费观看视频 | 黄色视频,在线免费观看| 亚洲不卡免费看| 日韩欧美 国产精品| 日本爱情动作片www.在线观看 | 欧美性猛交╳xxx乱大交人| 午夜精品国产一区二区电影 | 欧美日韩在线观看h| 最后的刺客免费高清国语| 欧美潮喷喷水| 久久久久久伊人网av| 三级毛片av免费| 又黄又爽又免费观看的视频| 国产精品一区二区三区四区久久| 亚洲人成网站在线播| 国产精品一区二区三区四区免费观看 | 久久九九热精品免费| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av成人精品一区久久| 日韩大尺度精品在线看网址| 亚洲中文字幕一区二区三区有码在线看| 国内揄拍国产精品人妻在线| 日韩成人av中文字幕在线观看 | av天堂中文字幕网| 嫩草影视91久久| 国产精品一区二区三区四区久久| 18禁裸乳无遮挡免费网站照片| 特大巨黑吊av在线直播| 内射极品少妇av片p| 又粗又爽又猛毛片免费看| 天堂动漫精品| 俺也久久电影网| 亚洲一区高清亚洲精品| 午夜视频国产福利| 1000部很黄的大片| 亚洲国产欧洲综合997久久,| 亚洲av电影不卡..在线观看| 此物有八面人人有两片| 春色校园在线视频观看| 国产午夜精品论理片| 亚洲一级一片aⅴ在线观看| 秋霞在线观看毛片| 精品一区二区三区人妻视频| 夜夜爽天天搞| 国产免费男女视频| 久久久久久大精品| 中文字幕久久专区| www.色视频.com| 久久精品夜色国产| 国产精品一区二区性色av| 狠狠狠狠99中文字幕| 国产成人影院久久av| 国产色爽女视频免费观看| 女同久久另类99精品国产91| 精品一区二区免费观看| 免费搜索国产男女视频| 人妻少妇偷人精品九色| 99九九线精品视频在线观看视频| 午夜精品在线福利| 天堂网av新在线| 午夜福利在线观看免费完整高清在 | 精品久久久久久久久久久久久| 日韩精品中文字幕看吧| 亚洲成人精品中文字幕电影| 国产一区二区激情短视频| 国产淫片久久久久久久久| 亚洲av成人精品一区久久| av在线蜜桃| 99热只有精品国产| 午夜免费激情av| 精品久久久久久久久久久久久| 国产精品久久久久久久电影| 一级毛片电影观看 | 欧美成人a在线观看| 天堂av国产一区二区熟女人妻| 久久亚洲精品不卡| 干丝袜人妻中文字幕| 12—13女人毛片做爰片一| 69av精品久久久久久| 久久久欧美国产精品| 日韩欧美三级三区| 久久精品国产亚洲av天美| 成人av一区二区三区在线看| 久久久精品94久久精品| 又爽又黄无遮挡网站| 国产在线男女| 99热这里只有是精品50| 成年版毛片免费区| 高清日韩中文字幕在线| 高清毛片免费观看视频网站| 晚上一个人看的免费电影| 白带黄色成豆腐渣| 国产精品一及| 亚洲国产色片| 在线国产一区二区在线| 男人舔奶头视频| 黑人高潮一二区| 简卡轻食公司| 真人做人爱边吃奶动态| 日日啪夜夜撸| 综合色丁香网| 1024手机看黄色片| 国产精品野战在线观看| 国产精品一区二区三区四区免费观看 | 99热只有精品国产| 晚上一个人看的免费电影| 男女啪啪激烈高潮av片| 国产欧美日韩一区二区精品| 三级经典国产精品| 免费看美女性在线毛片视频| 亚洲人成网站在线观看播放| 变态另类丝袜制服| 亚洲在线观看片| 久久久成人免费电影| 啦啦啦啦在线视频资源| 国产乱人视频| 在线免费十八禁| 中文亚洲av片在线观看爽| 99国产精品一区二区蜜桃av| 成人高潮视频无遮挡免费网站| 亚洲最大成人av| 欧美一级a爱片免费观看看| 热99在线观看视频| 亚洲美女搞黄在线观看 | 男女那种视频在线观看| 久久久成人免费电影| 特大巨黑吊av在线直播| 3wmmmm亚洲av在线观看| 国产精品伦人一区二区| 久久国产乱子免费精品| 久久99热6这里只有精品| 在线观看av片永久免费下载| 日本欧美国产在线视频| 寂寞人妻少妇视频99o| 欧美在线一区亚洲| 亚洲国产精品合色在线| 嫩草影院新地址| 18+在线观看网站| 国产亚洲91精品色在线| 亚洲av成人精品一区久久| 日韩欧美一区二区三区在线观看| 精品久久久久久成人av| 久久亚洲精品不卡| 日本在线视频免费播放| www.色视频.com| 中文字幕免费在线视频6| 2021天堂中文幕一二区在线观| 国产黄色视频一区二区在线观看 | 国产精品永久免费网站| avwww免费| 禁无遮挡网站| 伊人久久精品亚洲午夜| 国产高清有码在线观看视频| 男女那种视频在线观看| 深夜a级毛片| 女人十人毛片免费观看3o分钟| 99在线人妻在线中文字幕| 精品久久久久久久久久久久久| 日韩人妻高清精品专区| 亚洲国产精品成人综合色| 国产综合懂色| 无遮挡黄片免费观看| 最后的刺客免费高清国语| 真实男女啪啪啪动态图| 国产高清视频在线播放一区| av在线天堂中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 女同久久另类99精品国产91| 亚洲熟妇熟女久久| 国产视频内射| 黄色日韩在线| 99热只有精品国产| 国产视频内射| 在线a可以看的网站| 久久精品综合一区二区三区| 精品午夜福利视频在线观看一区| 成人亚洲欧美一区二区av| 一级毛片电影观看 | 少妇人妻一区二区三区视频| 女人被狂操c到高潮| 自拍偷自拍亚洲精品老妇| 久久精品国产亚洲网站| 国产男靠女视频免费网站| 最近2019中文字幕mv第一页| 国产三级中文精品| 色视频www国产| av天堂在线播放| 热99re8久久精品国产| 亚洲激情五月婷婷啪啪| 国产精品美女特级片免费视频播放器| 国产爱豆传媒在线观看| 日韩 亚洲 欧美在线| 日韩人妻高清精品专区| 麻豆精品久久久久久蜜桃| 欧美日韩在线观看h| 丝袜美腿在线中文| 观看免费一级毛片| 又粗又爽又猛毛片免费看| 亚洲av熟女| 亚洲av成人精品一区久久| 日本与韩国留学比较| 精品久久久久久久末码| 美女被艹到高潮喷水动态| 婷婷色综合大香蕉| 给我免费播放毛片高清在线观看| 亚洲av中文字字幕乱码综合| 欧美bdsm另类| 别揉我奶头~嗯~啊~动态视频| 尾随美女入室| 亚洲av免费在线观看| 一级av片app| 成人av一区二区三区在线看| 久久韩国三级中文字幕| 日韩欧美在线乱码| 亚州av有码| 精品99又大又爽又粗少妇毛片| 亚洲av中文av极速乱| 久久久久九九精品影院| 免费看光身美女| 欧美一区二区国产精品久久精品| 国产91av在线免费观看| 此物有八面人人有两片| 国产精华一区二区三区| 欧美绝顶高潮抽搐喷水| 国产av不卡久久| 女的被弄到高潮叫床怎么办| 亚洲激情五月婷婷啪啪| 日韩亚洲欧美综合| 国内精品美女久久久久久| 亚洲国产色片| 自拍偷自拍亚洲精品老妇| 亚洲精品456在线播放app| 精华霜和精华液先用哪个| 亚洲人成网站在线播| 国产蜜桃级精品一区二区三区| 午夜久久久久精精品| 国产成年人精品一区二区| 欧美最黄视频在线播放免费| 亚洲精品国产成人久久av| 亚洲无线观看免费| 久久6这里有精品| .国产精品久久| 一级毛片我不卡| 三级国产精品欧美在线观看| 99热全是精品| 久久精品久久久久久噜噜老黄 | 蜜桃亚洲精品一区二区三区| 精品国产三级普通话版| 精品一区二区三区av网在线观看| 国产一区二区在线av高清观看| 18+在线观看网站| 男女之事视频高清在线观看| 看十八女毛片水多多多| 日本欧美国产在线视频| 久久久午夜欧美精品| 男人狂女人下面高潮的视频| 色综合亚洲欧美另类图片| 成人美女网站在线观看视频| 噜噜噜噜噜久久久久久91| 搡老熟女国产l中国老女人| 国产高清有码在线观看视频| 熟女电影av网| 成人性生交大片免费视频hd| av在线观看视频网站免费| 特大巨黑吊av在线直播| 亚洲一级一片aⅴ在线观看| 成人欧美大片| 国产精品女同一区二区软件| 人人妻人人澡人人爽人人夜夜 | 女同久久另类99精品国产91| 别揉我奶头~嗯~啊~动态视频| 午夜福利在线在线| 人妻久久中文字幕网| 久久婷婷人人爽人人干人人爱| 国产在线精品亚洲第一网站| 色综合站精品国产| 不卡视频在线观看欧美| 国产精品久久视频播放| 麻豆久久精品国产亚洲av| 99热全是精品| 久久亚洲精品不卡| 中文字幕久久专区| 听说在线观看完整版免费高清| 日韩国内少妇激情av| 一区二区三区免费毛片| 日本免费一区二区三区高清不卡| 91精品国产九色| 亚洲综合色惰| 成人精品一区二区免费| 69av精品久久久久久| 亚洲电影在线观看av| 热99re8久久精品国产| 亚洲欧美日韩无卡精品| 一进一出抽搐动态| 看非洲黑人一级黄片| 精品人妻一区二区三区麻豆 | 精品久久久久久久久久久久久| 久久人人爽人人爽人人片va| 欧美+日韩+精品| 欧美成人一区二区免费高清观看| 国产女主播在线喷水免费视频网站 | 看免费成人av毛片| 九色成人免费人妻av| 亚洲在线自拍视频| 一级a爱片免费观看的视频| 午夜影院日韩av| 亚洲国产精品成人久久小说 | 久久久久国产网址| 日日撸夜夜添| 亚洲美女搞黄在线观看 | 校园春色视频在线观看| 精华霜和精华液先用哪个| 亚洲天堂国产精品一区在线| 国产精品爽爽va在线观看网站| 大香蕉久久网| 国产真实乱freesex| 大型黄色视频在线免费观看| 1024手机看黄色片| 国产高清视频在线播放一区| 国产午夜精品久久久久久一区二区三区 | 日韩强制内射视频| 老司机福利观看| 国产真实伦视频高清在线观看| 久久久久性生活片| 亚洲精品亚洲一区二区| 超碰av人人做人人爽久久| 亚洲欧美中文字幕日韩二区| 你懂的网址亚洲精品在线观看 | 亚洲欧美日韩东京热| 色综合站精品国产| 欧美潮喷喷水| 哪里可以看免费的av片| 精品久久久噜噜| 天堂动漫精品| 天堂网av新在线| 国产欧美日韩精品亚洲av| 亚洲在线观看片| 国产成年人精品一区二区| 国产色婷婷99| 亚洲熟妇熟女久久| 美女被艹到高潮喷水动态| 久久久久久久午夜电影| 午夜视频国产福利| a级毛片a级免费在线| 中国美白少妇内射xxxbb| 偷拍熟女少妇极品色| 久久久久免费精品人妻一区二区| 中文字幕熟女人妻在线| 中国美白少妇内射xxxbb| 免费大片18禁| 亚洲国产精品合色在线| 午夜福利在线观看吧| 国产免费一级a男人的天堂| 久久久久性生活片| eeuss影院久久| 久久精品国产亚洲av天美| 99久久精品热视频| 看黄色毛片网站| 夜夜爽天天搞| 亚洲精品日韩在线中文字幕 | av天堂中文字幕网| 成人亚洲精品av一区二区| 美女免费视频网站| 亚洲av电影不卡..在线观看| 搡女人真爽免费视频火全软件 | 你懂的网址亚洲精品在线观看 | 草草在线视频免费看| 精品人妻一区二区三区麻豆 | 久久久国产成人免费| 国产乱人视频| 国产亚洲欧美98| 3wmmmm亚洲av在线观看| 久久精品国产亚洲av天美| 亚洲一级一片aⅴ在线观看| 啦啦啦观看免费观看视频高清| 久久久精品94久久精品| 国产伦在线观看视频一区| 欧美激情在线99| 三级男女做爰猛烈吃奶摸视频| 国产亚洲精品久久久com| 在线观看av片永久免费下载| 麻豆av噜噜一区二区三区| 男女下面进入的视频免费午夜| 久久久久久久久久久丰满| 久久久久久久久久黄片| 啦啦啦韩国在线观看视频| 欧美性感艳星| 三级经典国产精品| 亚洲最大成人手机在线| 老熟妇仑乱视频hdxx| 国产亚洲欧美98| 亚洲欧美成人综合另类久久久 | 俺也久久电影网| 91久久精品电影网| aaaaa片日本免费| www日本黄色视频网| 久久久久久久久久久丰满| 欧美性猛交╳xxx乱大交人| 亚洲第一区二区三区不卡| 亚洲av不卡在线观看| 亚洲内射少妇av| .国产精品久久| 桃色一区二区三区在线观看| 亚洲国产欧洲综合997久久,| 欧美不卡视频在线免费观看| 精品无人区乱码1区二区| 精品乱码久久久久久99久播| 成人高潮视频无遮挡免费网站| 日韩 亚洲 欧美在线| 高清日韩中文字幕在线| 欧美绝顶高潮抽搐喷水| 乱系列少妇在线播放| 久久久久久伊人网av| 内射极品少妇av片p| 午夜爱爱视频在线播放| 久久久久久伊人网av| 两性午夜刺激爽爽歪歪视频在线观看| 成人av在线播放网站| 亚洲美女视频黄频| 欧美日韩国产亚洲二区| 亚洲国产精品成人久久小说 | 国产精华一区二区三区| 黄片wwwwww| 国产亚洲精品综合一区在线观看| 亚洲四区av| 热99在线观看视频| 欧美日韩国产亚洲二区| 成人性生交大片免费视频hd| 99热网站在线观看| 国产精品99久久久久久久久| 精品一区二区三区视频在线| 丰满的人妻完整版| 亚洲,欧美,日韩| 国产美女午夜福利| 亚洲av.av天堂| 国产探花极品一区二区| 少妇人妻精品综合一区二区 | 少妇裸体淫交视频免费看高清| 91久久精品电影网| 老师上课跳d突然被开到最大视频| 狠狠狠狠99中文字幕| 日日啪夜夜撸| 联通29元200g的流量卡| 18禁裸乳无遮挡免费网站照片|