基于無人機多光譜影像的完熟期玉米倒伏信息提取

李華森,夏晨真,張星宇,王 寅,3,張 月,3

(1.吉林農(nóng)業(yè)大學資源與環(huán)境學院,吉林 長春 130118;2.吉林省商品糧基地土壤資源可持續(xù)利用重點實驗室,吉林 長春 130118;3.秸稈綜合利用與黑土地保護教育部重點實驗室,吉林 長春 130118)

糧食是人類賴以生存的寶貴資源,糧食安全是保障國家長治久安的關(guān)鍵[1]。東北黑土區(qū)是我國重要的玉米生產(chǎn)基地,為我國糧食安全和糧食供給提供重要保障。玉米作為我國三大糧食作物之一,其莖稈較高且脆,易受大風、暴雨等極端天氣影響而發(fā)生倒伏,從而影響產(chǎn)量及籽粒品質(zhì)等。研究表明,倒伏是造成玉米減產(chǎn)的一個重要原因,一般能導致減產(chǎn)15%～30%,部分嚴重倒伏造成的損失可達到50%,甚至造成絕產(chǎn)[2-4]。因此,及時準確地獲取黑土區(qū)玉米倒伏發(fā)生的范圍和受災程度具有重要意義[5]。

通過傳統(tǒng)方法獲取作物的受災面積需要進行田間實地調(diào)查,耗費大量人力、物力和財力[6]。近年來遙感技術(shù)的快速發(fā)展,為大面積、高效提取作物倒伏信息提供了新的途徑[7-9]。如任紅玲等[10]基于HJ-1衛(wèi)星遙感影像,選取植被指數(shù)、紋理特征等,采用最大似然法對遙感影像進行監(jiān)督分類,最終確定玉米的倒伏面積。王立志等[11]利用多時相HJ-1B衛(wèi)星CCD多光譜影像分析了抽雄期玉米倒伏前后多種植被指數(shù)的變化,利用RVI比值植被指數(shù)進行數(shù)理統(tǒng)計及閾值劃分并最終確定玉米倒伏的受災范圍,總體分類精度為85.7%,Kappa系數(shù)為0.804。李宗南等[12]利用Worldview-2高分辨率影像,分析了灌漿期玉米倒伏地塊的紋理特征,從而識別出玉米的倒伏區(qū)域。但是,由于衛(wèi)星離地面距離較遠,且易受大氣等因素的干擾,所以利用衛(wèi)星平臺對作物倒伏的識別精度較低[13-15]。

近年來無人機遙感技術(shù)迅猛發(fā)展,與衛(wèi)星平臺相比,無人機技術(shù)具有諸多優(yōu)勢,例如受天氣影響小,能夠在多云情況下獲得清晰影像;分辨率高,可達厘米級,能實現(xiàn)對田塊級農(nóng)作物的精準識別;操作簡單,飛行高度低且時間靈活,能獲取中小尺度范圍內(nèi)高時空分辨率的影像數(shù)據(jù),有效彌補了衛(wèi)星遙感的不足,可滿足精準農(nóng)業(yè)的需求[16-17]。Liu等[18]利用無人機采集的熱紅外和RGB影像,構(gòu)建了水稻倒伏的識別模型,該模型的誤差率小于10%。毛智慧等[19]利用無人機獲取研究區(qū)的數(shù)字表面模型,并結(jié)合影像的色彩特征對圖像進行分類,從而得到玉米的倒伏信息,在以小區(qū)為單元的統(tǒng)計分析中誤差不到10%。Han等[20]使用無人機直接獲取了倒伏玉米的多光譜及可見光圖像,并通過提取紋理、植被指數(shù)、冠層結(jié)構(gòu)等參數(shù),構(gòu)建了2種用于倒伏信息提取的Logistic模型,最終利用2個模型判斷是否倒伏的概率分別為85%和90%。

上述識別玉米倒伏的方法中大多為利用倒伏及未倒伏玉米間的光譜差異,利用監(jiān)督分類的方法并結(jié)合紋理、植被指數(shù)等參數(shù)進行倒伏信息的提取。但是,適用于遙感影像的分類方法有很多種,包括面向像元的監(jiān)督分類法(如支持向量機、最大似然法、隨機森林、平行六面體等)、非監(jiān)督分類法(如ISODATA、K-Means等)、面向?qū)ο蟮姆诸惙椒ǖ?究竟哪種方法更適合于對玉米倒伏信息的提取尚不清楚。Logistic回歸模型具有訓練快、效果好、可解釋性強、便于觀測樣本概率分數(shù)等優(yōu)點,已廣泛應(yīng)用于金融、貿(mào)易經(jīng)濟、臨床醫(yī)學等方面,但在玉米倒伏信息提取應(yīng)用中的研究較少。因此,本文采用了面向?qū)ο蠓诸惙椒?、面向像元的最大似然法、及多元Logistic回歸模型3種方法,并結(jié)合影像原始光譜波段、光譜指數(shù)、和紋理特征等進行玉米倒伏信息的識別。

本文基于研究區(qū)的無人機多光譜影像,利用從影像中提取的原始光譜特征、光譜指數(shù)及紋理特征,對比面向?qū)ο蠓诸惙椒ā⒆畲笏迫环ê投嘣狶ogistic回歸3種模型,進行完熟期玉米倒伏區(qū)域的定量識別,并利用在影像上目視選擇的樣本進行倒伏信息識別精度的驗證,以期尋求區(qū)分不同玉米倒伏狀態(tài)的敏感特征及特征組合,研究結(jié)果可為玉米倒伏后受災面積的精準提取及災后農(nóng)業(yè)生產(chǎn)、政府決策等提供數(shù)據(jù)基礎(chǔ)和技術(shù)支持。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于吉林省梨樹縣四棵樹鄉(xiāng)的王家橋村(43°14′48.9″N,124°29′09.9″E)(圖1,見203頁)。土壤類型為黑土,主要種植作物為玉米。該區(qū)域?qū)俦睖貛О霛駶櫞箨懶约撅L氣候,無霜期155 d,年均氣溫6.5℃,年日照時長2 541 h,作物生長期內(nèi)(5—10月)時長1 192 h。2020年玉米生育期內(nèi)降水量約為533.4 mm。日均氣溫20.80℃,>10℃積溫3 030℃。

1.2 試驗設(shè)計

通過設(shè)置不同的種植密度、耕作方式、養(yǎng)分投入量和來源、肥料施用運籌等田間管理模式對各試驗小區(qū)進行區(qū)分。試驗共設(shè)置4個處理,各處理包括4次重復,具體包括處理1:對照(CK),不施肥處理;處理2:農(nóng)民習慣(FP),通過進行農(nóng)戶問卷調(diào)查的方式,獲得當?shù)剞r(nóng)民主流的玉米田間管理方式;處理3:高產(chǎn)栽培(HY),以最大化提高玉米產(chǎn)量為目標,設(shè)置最高的種植密度并結(jié)合最大的養(yǎng)分投入和多次施肥;處理4:土壤-作物系統(tǒng)綜合管理(ISSM)[21],根據(jù)研究區(qū)的氣候條件和土壤狀況,通過合理的養(yǎng)分供應(yīng)和相對簡化的施肥方式,實現(xiàn)提升土壤肥力、作物產(chǎn)量及資源利用效率的管理方式。試驗小區(qū)的面積為120 m2(20 m×6 m),處理總面積為1 920 m2,小區(qū)設(shè)置見圖2。供試玉米品種為良玉99。在試驗開展期間,于4月下旬—5月上旬播種,9月下旬—10月初收獲。

CK:對照 Control;FP: 農(nóng)民習慣栽培 Farmers habits cultivation model;HY:高產(chǎn)栽培 High-yield cultivation;ISSM:土壤-作物系統(tǒng)綜合管理栽培 Integrated soil-crop system management model圖2 試驗小區(qū)設(shè)置Fig.2 Settings of the experimental plots

1.3 數(shù)據(jù)的獲取與處理

選擇晴朗無云、無風或微風天氣,于2020年9月24日10∶00—15∶00時進行研究區(qū)無人機影像的采集。無人機平臺采用深圳市大疆創(chuàng)新科技有限公司生產(chǎn)的大疆Matrice 600 Pro六旋翼無人機,其最大載重為6 kg,最大水平飛行速度為65 km·h-1。傳感器選擇 parrot sequoia 農(nóng)業(yè)遙感專用4通道多光譜相機,該相機由多光譜傳感器和光照傳感器兩部分組成,包含4個120萬像素的單波段相機和1個1600萬像素的RGB相機。光譜儀鏡頭選擇聚焦25 mm,對應(yīng)的視場角約為13°。起飛前對無人機進行航線設(shè)定和白板校正,航高為100 m,空間分辨率為5 cm。影像采集時的航向重疊度80%,旁向重疊度60%。

影像采集后,利用 Pix4D Mapper 軟件進行影像拼接。幾何校正以大疆Matrice 600 Pro無人機拍攝的 RGB 正射影像為參考,均勻且隨機地在影像上選取部分點位進行幾何校正,確保校正后影像誤差小于0.5個像元。隨后進行輻射定標,輻射定標采用偽標準地物輻射糾正法,以地面1.5 m×1.5 m的白板為參考,利用ENVI 5.3軟件中的FLAASH模塊進行大氣校正,隨后得到研究區(qū)影像中各像元的反射率值。

通過實地測量以獲得研究區(qū)倒伏玉米、未倒伏玉米和半倒伏玉米的實際面積,分別為2 062.29、2 754.71、1 210.52 m2,用于后續(xù)驗證3種模型方法對玉米倒伏信息的提取精度。

1.4 玉米倒伏敏感特征的選取

本研究選取從無人機影像中提取的光譜指數(shù)和紋理特征等信息作為識別玉米倒伏信息的特征。其中光譜指數(shù)包括15種指數(shù),詳見表1。

紋理特征是指地物實體表面所共有的內(nèi)在特征,可反映區(qū)域內(nèi)的灰度變化,并且不依賴于地物的顏色與亮度[25]。這些量化值主要包括均值 (Mean)、方差 (Variance)、差異性(Dissimilarity)、熵 (Entropy)、二階矩 (Second Moment)和相關(guān)性 (Correlation)、逆差距 (Homogeneity)、對比度 (Contrast)。本研究采用灰度共生矩陣(Gray level co-occurrence matrix, GLCM)紋理分析法,基于二階概率統(tǒng)計方法(Co-occurrence measures)對這8種濾波特征在ENVI 5.3軟件中進行提取。

1.5 玉米倒伏信息的提取與精度驗證

1.5.1 面向?qū)ο蠓诸惙?面向?qū)ο蠓诸惙ㄊ侵覆捎梅指钏惴ㄔ诜诸惽皩⑼|(zhì)鄰近像元構(gòu)建為光譜、紋理、顏色特征等信息相似的同質(zhì)均勻?qū)ο?并將這個對象作為提取地物的最小單元,再通過每個對象上的特征信息對這些影像對象進行分類[26-27],進而實現(xiàn)對目標地物特征的識別。影像完成分割后,構(gòu)成影像的基本單元從像元變成了包含各種特征的對象,在一定程度上降低了混合像元對地物特征信息提取的干擾。本研究中的面向?qū)ο蠓诸惙椒ㄔ趀Cognition 9.0軟件中實現(xiàn)。

1.5.2 最大似然分類法 最大似然分類法是通過遙感影像的光譜統(tǒng)計特征,對類別樣本對應(yīng)的影像數(shù)據(jù)進行分類處理,并記錄影像數(shù)據(jù)隸屬各類別的數(shù)量。根據(jù)訓練樣本提取的特征向量,構(gòu)造出各類別的多維正態(tài)分布模型,即概率密度函數(shù)或者概率分布函數(shù)[28-30],最終通過最大似然分類方法進行判斷,得到較為準確的地物分類結(jié)果。

1.5.3 多元Logistic回歸模型 Logistic回歸模型是指研究某個無序或有序多分類的目標變量與有關(guān)因素的關(guān)系,它分為二分類和多分類兩種回歸模型[6]。其中多元Logistic回歸模型是一種因變量為多分類的Logistic回歸,其模型構(gòu)建簡單且推廣度較好。本文中將倒伏玉米設(shè)置為1,半倒伏玉米設(shè)置為2,未倒伏玉米設(shè)置為3,進行模型的構(gòu)建。

最后,本研究采用較為常用的混淆矩陣方法進行玉米倒伏信息提取結(jié)果的精度評價。為保證本次樣本選擇的有效性,本研究通過手動方式在研究區(qū)均勻選取400個樣本點,并同步應(yīng)用于對3種分類方法的精度驗證。

2 結(jié)果與分析

2.1 不同田間管理條件下的玉米倒伏狀態(tài)

實地測量倒伏面積時發(fā)現(xiàn),倒伏玉米多分布在CK和FP處理中,半倒伏玉米多分布在HY處理中,未倒伏玉米多分布在ISSM處理中。在極端天氣下玉米倒伏不僅與天氣情況相關(guān),也與不同的田間管理條件相關(guān)。當作物氮肥施用不足(如CK和FP處理),植株根系發(fā)育較差,面臨惡劣氣候時易發(fā)生倒伏現(xiàn)象;但氮肥施用量也不是越高越好,當其處于過量狀態(tài)時(如HY處理),會造成植株徒長,植株整體韌性較弱,同樣容易造成倒伏,因此適量施用氮肥有利于玉米抗倒伏。

2.2 不同玉米倒伏狀態(tài)的光譜特征差異

本研究通過計算每種類型的玉米倒伏狀態(tài)所對應(yīng)的像元在不同波段的光譜反射率均值,繪制出3種待識別類型(倒伏、半倒伏、未倒伏)在不同波段的光譜反射率平均值曲線。由圖3可知,倒伏玉米在各波段的反射率值均大于半倒伏、未倒伏玉米,倒伏玉米與半倒伏玉米之間的光譜反射率差異較大,半倒伏與未倒伏玉米之間的光譜反射率差異相對較小,半倒伏玉米與其他狀態(tài)玉米在紅光波段的差值較大。

2.3 不同玉米倒伏狀態(tài)的紋理特征差異

通過計算采樣區(qū)內(nèi)每種玉米倒伏形態(tài)的紋理特征均值,本文分析了不同倒伏形態(tài)間紋理特征的相對差異度。如表2所示,對比倒伏與未倒伏玉米的紋理特征,差異最大的是相異性指標,兩種倒伏形態(tài)的相對差值為0.28,比其他紋理特征平均高出0.13;對比倒伏與半倒伏玉米的紋理特征,差異最大的是對比度指標,兩種倒伏形態(tài)的相對差值為0.42,比其他紋理特征平均高出0.26;對比半倒伏與未倒伏玉米的紋理特征,差異最大的是對比度指標,兩種倒伏形態(tài)的相對差值為0.19,比其他紋理特征平均高出0.12;綜上,區(qū)分3種玉米倒伏形態(tài)的難點在于對未倒伏與半倒伏玉米形態(tài)的區(qū)分。

2.4 不同玉米倒伏狀態(tài)的光譜指數(shù)差異

光譜指數(shù)是兩個或多個波段反射率進行線性或非線性的組合,不同光譜指數(shù)均值在面對倒伏的趨勢上有所不同,見表3,本文中區(qū)分倒伏較敏感的指數(shù)有NDVI、GNDVI、SRVI等。NDVI從未倒伏到半倒伏、倒伏,其值越來越小,GNDVI與NDVI趨勢相同,越趨向于倒伏,其值就越小。SRVI值從未倒伏到半倒伏先增加,再到倒伏狀態(tài)時,反而減小至小于未倒伏狀態(tài)。

表3 不同倒伏狀態(tài)植被指數(shù)均值的趨勢Table 3 Trend of mean vegetation index in different lodging states

2.5 玉米倒伏信息的提取

2.5.1 面向?qū)ο蠓椒?在利用面向?qū)ο蠓椒ㄟM行玉米倒伏信息的提取過程中,多尺度分割中閾值如果過小會導致同一地物被過度分割,形成過多分類單元;如果將閾值設(shè)置過大會導致分類單元欠分割,不同的地物會被分割成同一對象[31]。經(jīng)過多次嘗試,反復目視比對得到如下分類參數(shù):固定緊密度和形狀指數(shù)的權(quán)重,得到最優(yōu)的分割尺度為100;固定分割尺度和緊密度的權(quán)重,得到最優(yōu)的形狀指數(shù)為0.3;固定分割尺度和形狀指數(shù)的權(quán)重,得到最優(yōu)的緊密度指數(shù)為0.5。因此,最終得到影像分割的最佳尺度為100,緊密度權(quán)重為0.5,形狀指數(shù)權(quán)重為0.3(圖4)。然而,分割后的影像仍然有部分區(qū)域與實際地物不符,主要是由于同一地物可能包含不同的亮度,易造成漏分和錯分的情況[26]。因此本研究還另外加入了15種光譜指數(shù)和8種紋理特征參與后續(xù)分割。后續(xù)分割過程中,經(jīng)多次目視試錯后確定將無人機影像的權(quán)重設(shè)為0.5,光譜指數(shù)的權(quán)重設(shè)為0.5,紋理特征的權(quán)重設(shè)為1,能夠達到較好的分割效果(圖5,見205頁)。最后,將無人機影像、8種紋理特征及15種光譜指數(shù)全部輸入到eCognition 9.0軟件中,對玉米的倒伏信息進行提取及精度驗證。在提取倒伏面積過程中發(fā)現(xiàn)對區(qū)分倒伏與未倒伏狀態(tài)較為敏感的光譜指數(shù)為NDVI,區(qū)分倒伏與未倒伏狀態(tài)較敏感的紋理特征為相異性。

圖4 不同參數(shù)設(shè)置下的影像分割結(jié)果(研究區(qū)局部)Fig.4 Image segmentation results under different parameters (Part of the study area)

2.5.2 最大似然法 在ENVI 5.3軟件中,對玉米倒伏、未倒伏、半倒伏3種類型分別選擇訓練感興趣區(qū)。其中,倒伏、未倒伏訓練感興趣區(qū)的數(shù)量均為60個,半倒伏訓練感興趣區(qū)數(shù)量為80個,共計200個樣本點。判斷樣本點的選取是否合格以及其是否能夠進行較為準確的分類,可以通過驗證樣本之間的可分離性來進行。如果兩個感興趣區(qū)組合的可分離性值在0～2.0之間,且小于1.8則需要重新選擇樣本;小于1則考慮將兩類樣本合并成一類[32-33],各樣本的分離性見表4,其中,半倒伏與倒伏狀態(tài)之間分離性較小,說明區(qū)分半倒伏與倒伏狀態(tài)較難,其余各類之間的分離性均較好。而后通過最大似然分類法進行分類,在選擇待分類影像后,進行參數(shù)設(shè)置:選中所有分類類別,設(shè)置似然度閾值為無,數(shù)據(jù)比例系數(shù)為255,輸出分類結(jié)果后進行精度評價,并生成混淆矩陣以驗證提取玉米倒伏信息的精度。

2.5.3 多元Logistic模型 將多光譜影像導入ENVI 5.3軟件中,隨機選取200個樣本點作為訓練集用于建模,隨后分別獲取訓練集200個點在紅、綠、藍3個波段的光譜均值,用SPSS對樣本點進行主成分分析,提取到1個主成分,其貢獻率為98.07%,計算式為:

F=0.992Rgreen+0.992Rbule+0.986Rred

式中,F為可見光波段均值的主成分提取值,Rgreen、Rblue、Rred分別為綠、藍、紅波段的均值。完成主成分分析后進行多元Logistic分析,因變量為倒伏狀態(tài),分別為倒伏玉米設(shè)置為1,半倒伏玉米設(shè)置為2,未倒伏玉米設(shè)置為3。自變量為主成分,在模型中勾選偽R方、模型擬合度信息等,參數(shù)設(shè)置置信區(qū)間為95%,最大迭代次數(shù)為100。結(jié)果顯示本模型擬合度的顯著性為0.021,說明模型具有統(tǒng)計意義,但偽R方表格中的值偏低,最大為0.55,擬合程度不佳。

2.6 提取結(jié)果的精度評價

利用3種方法提取的玉米倒伏信息精度對比如表5所示?；诿嫦?qū)ο蠓椒ǖ淖R別結(jié)果可知,半倒伏的漏分誤差較高,大部分被錯分的半倒伏對象被分為未倒伏狀態(tài),少部分被錯分為倒伏狀態(tài),分類總體精度能達到88.13%,Kappa系數(shù)能達到0.83,分類精度較高。在利用最大似然分類法的玉米倒伏精度評價過程中發(fā)現(xiàn),倒伏玉米多被錯分為半倒伏玉米,分類總體精度為73.51%,Kappa系數(shù)為0.65。在利用多元Logistic回歸模型的玉米倒伏精度評價中,倒伏玉米被錯分為半倒伏玉米,未倒伏玉米同樣多被錯分為半倒伏玉米,分類總體精度為77.64%,Kappa系數(shù)為0.71。

表5 3種方法提取的玉米倒伏信息精度對比Table 5 Comparison of accuracy of lodging information extracted by three methods

利用3種方法對玉米倒伏信息的識別情況與實際倒伏情況對比如表6所示。相比而言,面向?qū)ο蠓椒ǚ謩e提取的玉米倒伏、半倒伏和未倒伏面積與實際面積差距最小,總體精度為88.13%,Kappa系數(shù)為0.83,最大似然法與多元Logistic回歸模型總體精度分別為73.51%和77.64%?？傮w分析發(fā)現(xiàn)基于面向?qū)ο蠓ㄌ崛〉挠衩椎狗畔⒕茸罡?因此,本研究利用此方法獲得研究區(qū)玉米倒伏信息的最終識別結(jié)果圖(圖6)。

表6 3種方法提取玉米倒伏面積與實際面積的對比/m2Table 6 Comparison of the maize lodging extraction areas based on the three methods and actual lodging areas

圖6 最終分類圖Fig.6 Final classification results of maize lodging

3 討 論

本研究在實地測量玉米的倒伏面積時,倒伏玉米多分布在CK和FP處理中,半倒伏玉米多分布在HY處理中,未倒伏玉米多分布在ISSM處理中,由于CK處理沒有進行土壤旋耕和施用化肥及有機肥,植株間水肥競爭過于激烈,作物主根發(fā)育不良,扎根較淺,造成倒伏較嚴重,因此CK處理地塊中的倒伏玉米面積較大;HY處理中作物種植密度相對較大,且養(yǎng)分投入過多,因此植株的生長空間較小,透風透光環(huán)境較差。同時,由于氮肥施用過量造成植株徒長、根莖韌性較弱,在遭遇極端天氣的情況下,也易發(fā)生倒伏;ISSM處理中,由于根據(jù)區(qū)域氣候和土壤狀況設(shè)計了作物生產(chǎn)體系以及對應(yīng)的管理辦法,并通過合理密植、適量的養(yǎng)分供應(yīng)提升了作物產(chǎn)量以及資源利用效率,還提高了植株莖稈質(zhì)量,減少了植株倒伏情況的發(fā)生。

倒伏與半倒伏玉米在光譜特征和紋理特征方面的差異均較大,較易區(qū)分;半倒伏玉米與未倒伏玉米在上述兩方面特征上差異較小,所以本研究中分類的難點就在于對半倒伏與未倒伏玉米狀態(tài)之間的區(qū)分。本文中利用對倒伏敏感的光譜指數(shù)(NDVI、GNDVI、SRVI等)與紋理特征(相異性、對比度)的組合進行分類,所得到的分類結(jié)果較好,精準度較高。

相較于衛(wèi)星影像,無人機多光譜影像可提供更加精細的地物形狀及紋理特征等信息,且無人機具有影像分辨率高、靈活性強、時效性高等優(yōu)點,能夠為獲取倒伏作物的面積提供精準的數(shù)據(jù)源。本文基于從無人機多光譜影像中獲取的光譜特征、紋理特征及光譜指數(shù),并應(yīng)用面向?qū)ο蠓诸惙?、最大似然分類法和多元Logistic回歸模型3種方法,對玉米倒伏、半倒伏、未倒伏信息進行定量提取。其中,利用面向?qū)ο蠓诸惙▽τ衩椎狗畔⑻崛〉木茸罡?其總體精度可達到88.13%,Kappa系數(shù)為0.83。由于面向?qū)ο蟮姆诸惙椒蓪⒕哂邢嗨乒庾V及紋理特征的對象劃分到同一類別中,避免了面向像元分類方法中易出現(xiàn)的“同物異譜、同譜異物”現(xiàn)象,因此其對玉米倒伏信息的提取精度高于最大似然分類法的精度(總體精度=73.51%,Kappa系數(shù)=0.65);但是在利用面向?qū)ο蠓椒ㄟM行玉米倒伏信息的提取時,識別精度在很大程度上依賴于對象的分割尺度。而分割尺度多是由人工判斷,因此多尺度分割標準的界定主觀性較強。由于分割尺度會直接影響到地物信息的識別精度,所以后續(xù)的研究將考慮對分割結(jié)果的評價,以進一步提高玉米倒伏信息的識別精度。利用多元Logistic回歸模型提取倒伏信息的總體準確率為77.64%,其結(jié)果優(yōu)于最大似然法73.51%的準確率,但低于面向?qū)ο蠓诸惙椒?8.13%的準確率。原因分析如下:在進行多元Logistic回歸分析時,由于模型本身形式較為簡單,難以擬合數(shù)據(jù)的真實分布,從而影響了分類的精確度,分析結(jié)果的偽R方的數(shù)據(jù)為0.55,這意味著模型僅能解釋因變量55%的變化,對原始變量的變異解釋稍顯不足,還有部分信息無法解釋,因此模型的分類精度較低。

張新樂等[15]的研究總體精度為90.18%,Kappa系數(shù)為0.84,略高于本研究精度(總體精度=88.13%,Kappa系數(shù)=0.83)。該研究利用無人機多光譜影像識別完熟期玉米的倒伏信息,采用植被指數(shù)、光譜和紋理特征等不同組合并應(yīng)用最大似然法提取玉米倒伏面積,其通過多類紋理特征向量集提取的玉米倒伏面積較為準確。與之相比,本研究應(yīng)用了更多的光譜指數(shù)參與分類,且對比了3種不同分類方法對玉米倒伏信息的提取精度,同時更好地區(qū)分了不同倒伏狀態(tài)的差異。孫乾[2]的研究總體精度為86.61%,Kappa系數(shù)為0.83;其研究通過灌漿期玉米的多光譜影像提取了4種光譜指數(shù)、主成分以及紋理特征,并通過最大似然法對紋理特征與光譜指數(shù)組合進行分類,應(yīng)用特征組合方式分類玉米倒伏,相比單獨使用某種分類特征進行分類精度要高。本研究同樣通過結(jié)合紋理特征、光譜特征以及光譜指數(shù)組合進行分類,應(yīng)用紋理特征盡可能避免了應(yīng)用光譜特征識別的“同物異譜、同譜異物”現(xiàn)象;與之相比,本文應(yīng)用的光譜指數(shù)易于區(qū)分不同倒伏狀態(tài)之間的差異,分類精度要稍高一些。綜上,利用面向?qū)ο蟮姆椒ㄟM行田塊尺度玉米倒伏的識別是可行的。

4 結(jié) 論

1)基于無人機多光譜影像,能夠?qū)崿F(xiàn)對田塊尺度玉米倒伏信息的精準識別,且識別精度較利用衛(wèi)星影像識別倒伏信息的精度高。

2)對比3種方法對玉米倒伏信息的提取精度,面向?qū)ο蠓椒ǖ木茸罡?主要是由于其利用了從高分辨率的無人機影像中提取的光譜指數(shù)及紋理特征,避免了面向像元方法中的混合像元問題,增強了對玉米倒伏和半倒伏狀態(tài)的辨識能力。

3)區(qū)分玉米倒伏和未倒伏狀態(tài)的最佳特征為植被指數(shù)中的NDVI和紋理特征中的相異性,結(jié)合使用植被指數(shù)和紋理特征后,所獲得的分類結(jié)果較好,精度較高。