基于衛(wèi)星遙感與同步地面實(shí)驗(yàn)數(shù)據(jù)的冬小麥長(zhǎng)勢(shì)遙感監(jiān)測(cè)

蔡慶空 李二俊 王果 張迪 陳超

摘? 要：農(nóng)作物長(zhǎng)勢(shì)的有效監(jiān)測(cè)對(duì)于精準(zhǔn)農(nóng)業(yè)建設(shè)意義重大。基于Landsat8衛(wèi)星影像和同步實(shí)驗(yàn)數(shù)據(jù)，采用雙變量相關(guān)性分析法篩選植被指數(shù)，構(gòu)建LAI反演模型并進(jìn)行精度評(píng)定。結(jié)果表明：DVI與LAI的相關(guān)性最高，其次為SAVI和EVI，NRI與LAI相關(guān)性最低。基于DVI建立的模型精度最高，模型R2為0.77。研究區(qū)冬小麥長(zhǎng)勢(shì)總體良好，其中長(zhǎng)勢(shì)一般的區(qū)域占44.37%，長(zhǎng)勢(shì)良好的區(qū)域主要分布在扶風(fēng)縣中南部，武功縣中部，長(zhǎng)勢(shì)過旺地區(qū)主要分布在扶風(fēng)縣西部地區(qū)，武功縣西北部地區(qū)。

關(guān)鍵詞：冬小麥；葉面積指數(shù)；植被指數(shù)；遙感反演；長(zhǎng)勢(shì)監(jiān)測(cè)

中圖分類號(hào)：TP391? ? 文獻(xiàn)標(biāo)識(shí)碼：A? ? ? 文章編號(hào)：2096-4706（2023）06-0098-05

Remote Sensing Growth Monitoring of Winter Wheat Based on Satellite Remote Sensing and Synchronous Ground Experiment Data

CAI Qingkong1， LI Erjun2， WANG Guo1， ZHANG Di1， CHEN Chao1

（1.College of Civil Engineering， Henan University of Engineering， Zhengzhou? 451191， China;

2.School of Humanities， Political Science and Law， Henan University of Engineering， Zhengzhou? 451191， China）

Abstract： The effective monitoring of crop growth is of great significance to the construction of precision agriculture. Based on Landsat8 satellite image and synchronous experimental data， bivariate correlation analysis method is used to screen vegetation index， build LAI inversion model and evaluate accuracy. The results show that DVI has the highest correlation with LAI， followed by SAVI and EVI， and NRI has the lowest correlation with LAI. The built model based on DVI has the highest accuracy， with R2 of 0.77. The winter wheat in the study area is growing well in general， of which 44.37% is in the area with general growth. The areas with good growth are mainly distributed in the middle and south central parts of Fufeng County， the middle part of Wugong County， and the areas with excessive growth are mainly distributed in the west part of Fufeng County and the northwest part of Wugong County.

Keywords： winter wheat; leaf area index; vegetation index; remote sensing inversion; growth monitoring

0? 引? 言

葉面積指數(shù)（Leaf Area Index， LAI）作為農(nóng)作物長(zhǎng)勢(shì)監(jiān)測(cè)的重要指標(biāo)，其數(shù)值與產(chǎn)量的高低密切相關(guān)，在農(nóng)作物長(zhǎng)勢(shì)監(jiān)測(cè)、產(chǎn)量估算等領(lǐng)域應(yīng)用廣泛[1]。傳統(tǒng)的LAI獲取方法費(fèi)時(shí)費(fèi)力，成本高，而且需要破壞性采樣。與傳統(tǒng)方法相比，遙感技術(shù)具有現(xiàn)勢(shì)性、便捷性、時(shí)效性、宏觀性等優(yōu)勢(shì)，能快速、準(zhǔn)確、動(dòng)態(tài)地獲取農(nóng)作物的長(zhǎng)勢(shì)信息，且精度有保證，因此在LAI估測(cè)中表現(xiàn)出強(qiáng)大的優(yōu)勢(shì)[2]。目前，常用的LAI反演方法主要有兩種：物理模型法和經(jīng)驗(yàn)?zāi)Ｐ头╗3]。物理模型法是在物理學(xué)理論的基礎(chǔ)上，模擬植被冠層光譜特征、大氣以及土壤等對(duì)電磁波的作用，來實(shí)現(xiàn)葉面積指數(shù)的反演[4]。苗乃哲等采取物理模型法對(duì)北京市通州和順義冬小麥葉面積指數(shù)進(jìn)行反演[5]。赫曉慧等利用PROSAIL模型和遺傳算法優(yōu)化后的BP神經(jīng)網(wǎng)絡(luò)模型，對(duì)大豆種群的LAI進(jìn)行了反演[6]。王梟軒等針對(duì)物理模型抗噪能力差且容易過擬合的問題，提出一種PROSAIL模型結(jié)合VMG多元回歸模型反演冬小麥葉面積指數(shù)方法，提高了LAI的反演精度[7]。包剛等采用徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)的物理模型方法對(duì)草地LAI進(jìn)行了高光譜反演[8]。物理模型法能較好的解釋事物的本質(zhì)，但物理模型參數(shù)較多，且部分參數(shù)獲取困難，應(yīng)用受到一定的限制[4]。經(jīng)驗(yàn)?zāi)Ｐ头ň哂斜憬菪浴r(shí)效性、現(xiàn)勢(shì)性、宏觀性等優(yōu)勢(shì)，在農(nóng)作物長(zhǎng)勢(shì)監(jiān)測(cè)中應(yīng)用較為廣泛。束美艷等分析了基于地面實(shí)驗(yàn)采集的冬小麥冠層光譜數(shù)據(jù)，構(gòu)建紅邊抗水植被指數(shù)，提高了作物L(fēng)AI的反演精度[9]。孫華林等在正常播期和晚播條件下選取最佳植被指數(shù)建立了LAI估算模型[10]。Moriondo等利用冬小麥的NDVI數(shù)據(jù)對(duì)意大利兩個(gè)省份的冬小麥進(jìn)行長(zhǎng)勢(shì)監(jiān)測(cè)，取得了較高的精度[11]。吾木提·艾山江等基于野外實(shí)測(cè)的光譜數(shù)據(jù)，利用植被指數(shù)和優(yōu)化光譜指數(shù)，建立了冬小麥LAI的回歸模型，對(duì)其進(jìn)行精度驗(yàn)證[12]。李軍玲等利用野外試驗(yàn)獲取的光譜數(shù)據(jù)與LAI數(shù)據(jù)，建立了不同生育時(shí)期冬小麥葉面積指數(shù)的遙感反演模型[13]。經(jīng)驗(yàn)?zāi)Ｐ途哂斜憬菪浴r(shí)效性、現(xiàn)勢(shì)性、宏觀性等優(yōu)勢(shì)在LAI的反演中得到了廣泛的應(yīng)用。

綜合以上研究，由于受實(shí)驗(yàn)條件的限制，以往LAI的研究大多基于地面實(shí)驗(yàn)采集的光譜數(shù)據(jù)，構(gòu)建LAI反演模型，而沒有同步獲取的遙感影像數(shù)據(jù)，或者是有衛(wèi)星影像但無同步地面實(shí)驗(yàn)數(shù)據(jù)，導(dǎo)致無法將地面實(shí)驗(yàn)建立的LAI模型推廣應(yīng)用至大區(qū)域范圍的遙感影像上。論文選取陜西省武功縣和扶風(fēng)縣為研究區(qū)，綜合利用星（Landsat8衛(wèi)星影像）-地（地面實(shí)驗(yàn)數(shù)據(jù)）同步實(shí)驗(yàn)數(shù)據(jù)，構(gòu)建冬小麥LAI反演模型，結(jié)合Landsat8衛(wèi)星影像將建立的LAI模型推廣應(yīng)用于遙感影像上，實(shí)現(xiàn)大區(qū)域范圍的冬小麥長(zhǎng)勢(shì)遙感監(jiān)測(cè)。

1? 研究區(qū)和數(shù)據(jù)源

1.1? 研究區(qū)概況

選取陜西省武功縣和扶風(fēng)縣為研究區(qū)域，如圖1所示，地理位置介于北緯東經(jīng)間。氣候類型屬溫帶大陸性季風(fēng)氣候，具有春季溫暖多風(fēng)，夏季炎熱多雨，秋季涼爽多連陰雨，冬季寒冷少雪等特征，年均氣溫12.9 ℃，年均降水量635.1 mm。區(qū)內(nèi)地勢(shì)總體較為平坦，最高點(diǎn)海拔約為540米，最低點(diǎn)海拔約為418米，土壤肥沃，土質(zhì)類型主要為重壤土和中壤土。境內(nèi)有渭河，漆水河、湋河、漠浴河流經(jīng)，降水是水資源的重要來源，占莊稼年水需求量的63.9%。主要種植農(nóng)作物是小麥、玉米、油菜和各種蔬菜。

1.2? 數(shù)據(jù)獲取及預(yù)處理

論文采用的數(shù)據(jù)主要包括Landsat8 OLI遙感影像數(shù)據(jù)和同步獲取的地面實(shí)驗(yàn)數(shù)據(jù)。Landsat8影像獲取時(shí)間為2014年5月11日（http：//www.gscloud.cn/），影像空間分辨率為30米，幅寬185千米，影像條帶號(hào)為127/36，云覆蓋度0.3%，天氣晴朗無云。遙感數(shù)據(jù)的預(yù)處理主要包括輻射定標(biāo)、大氣校正、幾何糾正以及影像裁剪。采用野外實(shí)驗(yàn)時(shí)采集的地面控制點(diǎn)對(duì)影像進(jìn)行幾何糾正，糾正誤差控制在1個(gè)像元之內(nèi)。同步獲取的地面實(shí)驗(yàn)數(shù)據(jù)主要包括：采樣點(diǎn)的地面坐標(biāo)和葉面積指數(shù)，葉面積指數(shù)采用LAI2000進(jìn)行采集，共采集有效實(shí)驗(yàn)樣點(diǎn)37個(gè)，其中隨機(jī)選取25個(gè)樣點(diǎn)作為訓(xùn)練樣本，其余12個(gè)樣本作為檢驗(yàn)樣本。

2? 冬小麥LAI模型構(gòu)建及長(zhǎng)勢(shì)監(jiān)測(cè)

2.1? 植被指數(shù)與LAI相關(guān)性分析

在借鑒前人研究的基礎(chǔ)上，選取8個(gè)與LAI密切相關(guān)的植被指數(shù)，分別是歸一化植被指數(shù)（NDVI）、差值植被指數(shù)（DVI）、增強(qiáng)型植被指數(shù)（EVI）、土壤調(diào)節(jié)植被指數(shù)（SAVI）、優(yōu)化的土壤調(diào)節(jié)植被指數(shù)（OSAVI）、比值植被指數(shù)（RVI）、綠色歸一化植被指數(shù)（GNDVI）和作物氮反應(yīng)指數(shù)（NRI）[5，14，15]。基于植被指數(shù)的計(jì)算公式，通過波段運(yùn)算，得到研究區(qū)8個(gè)植被指數(shù)影像。基于采集的樣本點(diǎn)坐標(biāo)，將樣本點(diǎn)分別導(dǎo)入到各植被指數(shù)影像上，提取樣本點(diǎn)位置的植被指數(shù)，并將提取的植被指數(shù)分為訓(xùn)練樣本和檢驗(yàn)樣本。將訓(xùn)練集樣本點(diǎn)植被指數(shù)與葉面積指數(shù)導(dǎo)入SPSS軟件，繪制散點(diǎn)圖，并進(jìn)行Pearson雙變量相關(guān)性分析，分析結(jié)果如表1所示。

由表1可知，8個(gè)植被指數(shù)與LAI相關(guān)性均達(dá)到了顯著性檢驗(yàn)水平，其中DVI、SAVI和EVI與LAI相關(guān)性最高，分別為0.79、0.76和0.76，其次為OSAVI、GNDVI、NDVI和RVI，NRI與LAI相關(guān)性最差，相關(guān)系數(shù)僅為0.53。因此，后續(xù)研究中采用NDVI、DVI、EVI、GNDVI、OSAVI、RVI和SAVI進(jìn)行進(jìn)一步分析。

2.2? 構(gòu)建LAI反演模型及精度評(píng)定

2.2.1? 建立LAI反演模型

基于訓(xùn)練集樣本，將LAI與植被指數(shù)在SPSS軟件中進(jìn)行回歸分析。對(duì)于每一個(gè)植被指數(shù)，函數(shù)模型選取一元線性模型、二次多項(xiàng)式模型、對(duì)數(shù)模型、冪函數(shù)模型和指數(shù)模型5種形式，從中選取R2最大的模型作為該植被指數(shù)模型的最終形式，模型建立的結(jié)果如表2所示。

由表2可知，從模型類別來看，指數(shù)模型和二次多項(xiàng)式模型的擬合效果最好，其次為對(duì)數(shù)模型和冪函數(shù)模型，一元線性模型的擬合效果最差；從植被指數(shù)方面來看，DVI、EVI、SAVI模型的反演精度優(yōu)于其他植被指數(shù)，NDVI的擬合效果最差，原因在于冬小麥灌漿期植被密度達(dá)到了一定程度，NDVI不再隨小麥生長(zhǎng)繼續(xù)增長(zhǎng)，到了飽和期。DVI、EVI以及SAVI的二次多項(xiàng)式模型反演精度最高，模型決定系數(shù)R2分別為0.77、0.74、0.73。

2.2.2? 模型精度評(píng)定

為了驗(yàn)證模型的可靠性和普適性，利用12個(gè)驗(yàn)證集樣本對(duì)反演精度最好的3個(gè)模型進(jìn)行精度評(píng)定。將利用3個(gè)模型反演得到的LAI與驗(yàn)證集樣本LAI在SPSS軟件中進(jìn)行Pearson雙變量相關(guān)性分析，得出3種模型反演LAI與野外測(cè)量LAI的相關(guān)系數(shù)，其中DVI和SAVI的相關(guān)性最高，相關(guān)系數(shù)均為0.71，EVI的相關(guān)系數(shù)為0.69。

2.3? 冬小麥種植區(qū)域提取

遙感影像上含有不同種類的地物，論文的研究對(duì)象是冬小麥，因此需要在影像上提取出冬小麥種植區(qū)域。采用監(jiān)督分類法提取冬小麥種植區(qū)域，通過目視判讀，將遙感影像上的地物分為建設(shè)用地、林地、水、冬小麥、裸地共5類，采集5類樣本的感興趣區(qū)，感興趣區(qū)要均勻的分布在影像上而且有一定的代表性，并計(jì)算樣本的可分離性。然后進(jìn)行分類器選擇，論文中選擇支持向量機(jī)方法進(jìn)行監(jiān)督分類，并對(duì)監(jiān)督分類結(jié)果進(jìn)行后處理，最終得到研究區(qū)土地分類圖，如圖2所示。基于監(jiān)督分類結(jié)果將冬小麥種植區(qū)域提取出來，為了更直觀的顯示冬小麥的分布情況，將其分為冬小麥區(qū)域和非小麥區(qū)域并著色，結(jié)果如圖3所示。

2.4? 冬小麥長(zhǎng)勢(shì)監(jiān)測(cè)

基于DVI、EVI以及SAVI的植被指數(shù)影像，將建立的二次多項(xiàng)式模型應(yīng)用到Landsat8遙感影像上，得到研究區(qū)葉面積指數(shù)分布圖，將葉面積指數(shù)分布圖與冬小麥種植區(qū)域圖疊加，得到研究區(qū)冬小麥種植區(qū)域葉面積指數(shù)分布圖。結(jié)合研究區(qū)樣點(diǎn)地面實(shí)驗(yàn)數(shù)據(jù)和以往研究結(jié)果，將冬小麥的長(zhǎng)勢(shì)情況分為3個(gè)等級(jí)：長(zhǎng)勢(shì)一般、長(zhǎng)勢(shì)良好、長(zhǎng)勢(shì)過旺[16]，基于此分級(jí)標(biāo)準(zhǔn)，將研究區(qū)冬小麥葉面積指數(shù)分布圖進(jìn)行分級(jí)并著色，得到研究區(qū)灌漿期冬小麥長(zhǎng)勢(shì)監(jiān)測(cè)分布圖，如圖4所示。

為定量分析研究區(qū)冬小麥長(zhǎng)勢(shì)情況，計(jì)算三種模型下研究區(qū)冬小麥LAI的平均值，并統(tǒng)計(jì)三種模型下長(zhǎng)勢(shì)一般、良好和過旺等級(jí)的面積，如表3所示。

由圖4和表3可知，研究區(qū)冬小麥LAI的平均值為4.86，屬于長(zhǎng)勢(shì)良好等級(jí)，表明研究區(qū)冬小麥總體長(zhǎng)勢(shì)良好。從整個(gè)研究區(qū)看，長(zhǎng)勢(shì)一般區(qū)域占總種植面積的44.39%，主要分布在扶風(fēng)縣北部，原因是研究區(qū)西北部地區(qū)地勢(shì)高，種植條件差；長(zhǎng)勢(shì)良好的區(qū)域占總種植面積的40.13%，主要分布在扶風(fēng)縣中南部，武功縣中部，該地區(qū)水源充足，土地肥沃，靠近楊凌農(nóng)業(yè)高新技術(shù)產(chǎn)業(yè)示范區(qū)，受到示范區(qū)的影響，中南部地區(qū)的冬小麥種植有較好的農(nóng)田管理，大部分冬小麥長(zhǎng)勢(shì)良好；長(zhǎng)勢(shì)過旺的區(qū)域占總種植面積的15.5%，冬小麥長(zhǎng)勢(shì)過旺地區(qū)較為分散，主要分布在扶風(fēng)縣西部地區(qū)，武功縣西北部地區(qū)，出現(xiàn)這種情況可能是因?yàn)榉N植密度太大，施肥過多等，生長(zhǎng)后期會(huì)出現(xiàn)透風(fēng)性和透光性差。

3? 結(jié)? 論

為了對(duì)冬小麥長(zhǎng)勢(shì)進(jìn)行及時(shí)有效的監(jiān)測(cè)，文中選取陜西省武功縣與扶風(fēng)縣為研究區(qū)，綜合遙感影像和同步獲取的地面實(shí)驗(yàn)數(shù)據(jù)，構(gòu)建了LAI反演模型，對(duì)研究區(qū)冬小麥長(zhǎng)勢(shì)進(jìn)行遙感監(jiān)測(cè)，取得的主要結(jié)論如下：

（1）通過將8個(gè)植被指數(shù)與葉面積指數(shù)進(jìn)行相關(guān)性分析，8個(gè)植被指數(shù)與LAI相關(guān)性達(dá)到了顯著性檢驗(yàn)水平，其中DVI、SAVI和EVI與LAI相關(guān)性最高，分別為0.79、0.76和0.76，其次為OSAVI、GNDVI、NDVI和RVI，NRI與LAI相關(guān)性最低，相關(guān)系數(shù)為0.53。

（2）從模型建立形式看，二次多項(xiàng)式模型和指數(shù)模型的擬合效果最好，其次為對(duì)數(shù)模型和冪函數(shù)模型，一元線性模型的擬合效果最差。從植被指數(shù)方面來看，DVI、EVI、SAVI三種植被指數(shù)的擬合效果最好，其中基于DVI的二次多項(xiàng)式模型精度最高，模型R2為0.77，模型反演LAI與實(shí)測(cè)LAI有較好的相關(guān)性。

（3）研究區(qū)冬小麥總體長(zhǎng)勢(shì)良好，其中長(zhǎng)勢(shì)良好的區(qū)域占總種植面積的40.13%，主要分布在扶風(fēng)縣中南部，武功縣中部；長(zhǎng)勢(shì)一般的區(qū)域占44.37%，主要分布在扶風(fēng)縣北部地區(qū)；長(zhǎng)勢(shì)過旺的區(qū)域占15.5%，長(zhǎng)勢(shì)過旺地區(qū)較為分散，主要分布在扶風(fēng)縣西部地區(qū)，武功縣西北部地區(qū)。

參考文獻(xiàn)：

[1] 王紀(jì)華，趙春江，黃文江，等.農(nóng)業(yè)定量遙感基礎(chǔ)與應(yīng)用 [M].北京：科學(xué)出版社，2008.

[2] 劉占宇，黃敬峰，王福民，等.估算水稻葉面積指數(shù)的調(diào)節(jié)型歸一化植被指數(shù) [J].中國農(nóng)業(yè)科學(xué)，2008（10）：3350-3356.

[3] 張佳華，杜育璋，劉學(xué)鋒，等.基于高光譜數(shù)據(jù)和模型反演植被葉面積指數(shù)的進(jìn)展 [J].光譜學(xué)與光譜分析，2012，32（12）：3319-3323.

[4] 王毅.基于經(jīng)驗(yàn)?zāi)Ｐ偷娜~面積指數(shù)反演精度對(duì)比 [D].成都：成都理工大學(xué)，2017.

[5] 苗乃哲.基于LandsatTM數(shù)據(jù)的冬小麥不同生育期葉面積指數(shù)反演方法精度比較 [D].西安：西安科技大學(xué)，2012.

[6] 赫曉慧，馮坤，郭恒亮，等.基于PROSAIL模型和遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)模型的不同大豆種群葉面積指數(shù)反演比較 [J].河南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2021，55（4）：698-706.

[7] 王梟軒，盧小平，楊澤楠，等.基于PROSAIL結(jié)合VMG模型的冬小麥葉面積指數(shù)反演方法 [J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2022，53（6）：209-216.

[8] 包剛，覃志豪，周義，等.基于高光譜數(shù)據(jù)和RBF神經(jīng)網(wǎng)絡(luò)方法的草地葉面積指數(shù)反演 [J].國土資源遙感，2012（2）：7-11.

[9] 束美艷，顧曉鶴，孫林，等.基于新型植被指數(shù)的冬小麥LAI高光譜反演 [J].中國農(nóng)業(yè)科學(xué)，2018，51（18）：3486-3496.

[10] 孫華林，耿石英，王小燕，等.晚播條件下基于高光譜的小麥葉面積指數(shù)估算方法 [J].光譜學(xué)與光譜分析，2019，39（4）：1199-1206.

[11] MORIONDO M，MASELLI F，BINDI M. A simple model of regional wheat yield based on ndvi data [J].European Journal of Agronomy，2006，26（3）：266-274.

[12] 吾木提·艾山江，尼加提·卡斯木，買買提·沙吾提.基于波段組合優(yōu)化光譜指數(shù)的冬小麥LAI估算 [J].江蘇農(nóng)業(yè)科學(xué)，2022，50（13）：207-218.

[13] 李軍玲，彭記永.不同生育時(shí)期冬小麥葉面積指數(shù)地面高光譜遙感模型研究 [J].麥類作物學(xué)報(bào)，2018，38（8）：979-987.

[14] 劉婷，包廣道，張大偉，等.基于植被指數(shù)的GF-2影像防護(hù)林快速提取研究 [J].森林工程，2018，34（6）：13-19.

[15] 譚昌偉，王紀(jì)華，朱新開，等.基于Landsat TM影像的冬小麥拔節(jié)期主要長(zhǎng)勢(shì)參數(shù)遙感監(jiān)測(cè) [J].中國農(nóng)業(yè)科學(xué)，2011，44（7）：1358-1366.

[16] 李衛(wèi)國，李正金，楊澄.基于CBERS遙感的冬小麥長(zhǎng)勢(shì)分級(jí)監(jiān)測(cè) [J].中國農(nóng)業(yè)科技導(dǎo)報(bào)，2010，12（3）：79-83.

作者簡(jiǎn)介：蔡慶空（1986—），男，漢族，河南南召人，講師，博士，研究方向：資源環(huán)境遙感、高光譜遙感；李二俊（1986—），女，漢族，河南溫縣人，講師，碩士，研究方向：資源環(huán)境遙感；王果（1986—），男，漢族，河南內(nèi)鄉(xiāng)人，副教授，博士，研究方向：攝影測(cè)量與遙感；張迪（1987—），男，漢族，河南鄧州人，副教授，博士，研究方向：攝影測(cè)量與遙感；陳超（1989—），男，漢族，河南開封人，講師，博士，研究方向：攝影測(cè)量與遙感。

收稿日期：2022-10-18

基金項(xiàng)目：河南省高等學(xué)校重點(diǎn)科研項(xiàng)目（21A420003，23A420001，22A420003）；河南省重點(diǎn)研發(fā)與推廣項(xiàng)目（222102320155）