陶惠林 徐良驥 馮海寬,3 楊貴軍,4 代 陽 牛亞超
(1.北京農(nóng)業(yè)信息技術(shù)研究中心農(nóng)業(yè)農(nóng)村部農(nóng)業(yè)遙感機(jī)理與定量遙感重點實驗室, 北京 100097;2.安徽理工大學(xué)測繪學(xué)院, 淮南 232001; 3.國家農(nóng)業(yè)信息化工程技術(shù)研究中心, 北京 100097;4.北京市農(nóng)業(yè)物聯(lián)網(wǎng)工程技術(shù)研究中心, 北京 100097)
葉面積指數(shù)(Leaf area index,LAI)是反映作物長勢的重要參數(shù),與作物產(chǎn)量具有緊密聯(lián)系,對農(nóng)業(yè)生產(chǎn)管理具有重要作用[1-2]。準(zhǔn)確、高效地監(jiān)測LAI能夠提高長勢監(jiān)測效果和產(chǎn)量預(yù)測精度[3-5]。因此,對LAI進(jìn)行動態(tài)監(jiān)測顯得尤為重要。傳統(tǒng)測量LAI需要實地采集數(shù)據(jù),測量中會對作物造成損傷,同時需要消耗大量的人力物力,并且測量范圍有限[6]。
遙感技術(shù)憑借宏觀、適時和動態(tài)等特點在作物L(fēng)AI監(jiān)測中得到廣泛應(yīng)用。由于平臺的不同,采用遙感技術(shù)監(jiān)測LAI主要分為高空、低空和地面尺度。高空主要利用衛(wèi)星進(jìn)行遙感監(jiān)測,衛(wèi)星可以獲取大范圍的遙感數(shù)據(jù),對于較大區(qū)域的監(jiān)測具有較好效果,但衛(wèi)星運行周期長、空間分辨率低、容易受到云層的影響,使衛(wèi)星遙感監(jiān)測精度有限[7-10]。地面主要通過地物光譜儀獲取遙感數(shù)據(jù),在操作過程中由于平臺的高度限制,很難得到正射影像數(shù)據(jù)[11]。低空主要采用載人飛機(jī)和無人機(jī)獲取數(shù)據(jù),相比載人飛機(jī),無人機(jī)機(jī)動靈活、更加安全,其操作簡單、對起飛場地要求較低,因此,無人機(jī)遙感在農(nóng)業(yè)中得到了廣泛應(yīng)用[12-14]。目前,無人機(jī)攜帶的傳感器主要包括數(shù)碼相機(jī)、多光譜相機(jī)和高光譜相機(jī),數(shù)碼相機(jī)和多光譜相機(jī)獲取的影像波段較少、得到的光譜信息有限,而高光譜相機(jī)能夠獲取較多光譜信息,更適于監(jiān)測作物L(fēng)AI[15-19]。
采用無人機(jī)高光譜遙感技術(shù)監(jiān)測LAI需要使用植被指數(shù)。植被指數(shù)是由2個或多個光譜通過一定方式組合而成,能夠有效反映植被狀況,在遙感監(jiān)測中具有重要的意義。在無人機(jī)高光譜高效監(jiān)測LAI方面,研究者進(jìn)行了大量的研究工作[20-24]。但這些基于高光譜遙感的研究主要通過單個植被指數(shù)或多個植被指數(shù)估算LAI,而利用植被指數(shù)結(jié)合作物高度估算LAI的研究還很少。為了提高LAI估算精度,本文使用無人機(jī)高光譜遙感影像生成作物表面模型(CSM),提取出冬小麥植株高度(Hcsm),然后基于植被指數(shù)、植被指數(shù)結(jié)合Hcsm,使用偏最小二乘回歸(PLSR)方法構(gòu)建LAI估算模型,探究利用無人機(jī)高光譜和Hcsm估算LAI的方法,以期為提高作物L(fēng)AI估算精度提供一種新的研究手段。
田間試驗在北京市昌平區(qū)小湯山鎮(zhèn)國家精準(zhǔn)農(nóng)業(yè)研究示范基地(北緯40°10′48″~40°10′54″,東經(jīng)116°26′51″~116°26′53″)進(jìn)行。該區(qū)域?qū)儆谂瘻貛Ш桶霛駶櫞箨懶约撅L(fēng)氣候,降水多在夏季和秋季,年平均降水量約42 mm,年平均溫度約11.8℃,試驗田前茬作物為玉米,土壤為潮土類型,土質(zhì)比較肥沃。進(jìn)行小區(qū)試驗,采用了中麥175(ZM175,中國農(nóng)業(yè)科學(xué)院)和京麥9843(J9843,北京市農(nóng)業(yè)農(nóng)村局)2種冬小麥品種,試驗田共有48個小區(qū),每個小區(qū)面積為48 m2,其中16個小區(qū)為一個重復(fù)區(qū),每個重復(fù)區(qū)都進(jìn)行不同程度的氮肥和水分處理,其中氮肥設(shè)置了4種水平(N1: 0 kg/hm2;N2: 195 kg/hm2;N3: 390 kg/hm2;N4: 585 kg/hm2),水分設(shè)置了3種灌溉(W0:僅雨水;W1:灌溉量675 m3/hm2;W2:灌溉量1 012.5 m3/hm2),詳細(xì)的田間試驗設(shè)計如圖1所示。
于2015年4月21日(拔節(jié)期)、2015年4月26日(挑旗期)和2015年5月13日(開花期),選擇天氣晴朗、無風(fēng)無云的時間進(jìn)行無人機(jī)高光譜遙感作業(yè),飛行高度為80 m。試驗無人機(jī)搭載的傳感器為Cubert UHD185 Firefly型成像光譜儀(測量過程中保證光譜儀垂直向下),成像光譜儀由德國生產(chǎn),質(zhì)量470 g;由于是高光譜,有較多的波段,波長范圍450~950 nm;光譜分辨率為8 nm@532 nm,空間分辨率為1.92 cm;與地面高光譜不同,每個波段間隔4 nm;有125個光譜通道;成像速度為每秒拍攝5個高光譜圖像立方體。成像的光譜分辨率高,但空間分辨率并不高。為了去除土壤背景的影響,需對影像重采樣。因為成像光譜儀曝光時間具體取決于太陽光強(qiáng)度,因此需要進(jìn)行校正。在進(jìn)行無人機(jī)遙感作業(yè)前,在地面利用黑白板進(jìn)行輻射定標(biāo)。在進(jìn)行遙感作業(yè)時各時期的飛行航線保持一致。
無人機(jī)高光譜遙感數(shù)據(jù)獲取后需要進(jìn)行處理,數(shù)據(jù)處理主要包括兩部分:①高光譜影像的輻射校正和拼接。需要先將影像像元亮度(Digital number,DN)轉(zhuǎn)換為地表反射率[25],再通過Agisoft LLC公司生產(chǎn)的Agisoft PhotoScan 軟件對影像進(jìn)行拼接[26]。②提取冠層光譜反射率。為了避免試驗小區(qū)邊緣和田壟之間背景光譜的影響,使用ArcGIS軟件根據(jù)小區(qū)面積繪制出20個矢量,每個小區(qū)得到20個中心像元,對矢量進(jìn)行編號,結(jié)合IDL語言提取出所有矢量感興趣區(qū)像元值的光譜反射率,統(tǒng)計出各小區(qū)的平均光譜反射率,將獲取的矢量感興趣區(qū)的像元值平均光譜作為各試驗小區(qū)冬小麥冠層的平均光譜反射率。
地面數(shù)據(jù)獲取與無人機(jī)遙感作業(yè)同步進(jìn)行,獲取了實測的冬小麥LAI和株高數(shù)據(jù)。為了測量LAI,在每個小區(qū)長勢均勻區(qū)域采樣,選取20株植株,將樣本放入密封袋中帶回實驗室;經(jīng)過莖葉分離處理,通過美國CID生物科技公司生產(chǎn)的CI-203型激光葉面積儀,測定采樣樣本葉片面積,得到總的葉面積;獲得小區(qū)單位面積的單莖數(shù),乘以總面積得總莖數(shù),進(jìn)而得到LAI。冬小麥株高實測時間分別為2015年4月14日(拔節(jié)期)、2015年4月26日(挑旗期)和2015年5月13日(開花期),由于天氣原因,其中無人機(jī)獲取的拔節(jié)期影像時間為2015年4月21日。為了測量不同小區(qū)的冬小麥株高,在每個小區(qū)對角線的1/3和2/3處取4個測量點,用直尺測量,每個測量點取1株冬小麥,取平均值作為測量小區(qū)的冬小麥株高。
目前,植被指數(shù)的種類有很多,通過篩選不同的植被指數(shù),并根據(jù)已有研究成果,選取了一些在LAI監(jiān)測方面效果較好的植被指數(shù),如LCI[27]、NPCI[28]、MCARI[29]、TCARI[29]、PBI[30]、BGI[31]、TVI[32]、OSAVI[33]、NDVI[34]、SR[35]這10種植被指數(shù),用于構(gòu)建LAI估算模型,從而監(jiān)測LAI。
采用偏最小二乘(Partial least squares regression,PLSR)方法分析冬小麥生育期的遙感數(shù)據(jù),構(gòu)建出LAI估算模型。PLSR是將多元線性回歸、典型相關(guān)分析和主成分分析結(jié)合為一體,可以提供一種多對多的線性回歸建模方法,特別當(dāng)變量個數(shù)較多時,存在多重相關(guān)性,而觀測變量數(shù)據(jù)較少時,利用PLSR建立的模型具有傳統(tǒng)的典型回歸分析所沒有的優(yōu)點。利用最小化誤差的平方和構(gòu)建最佳的模型,以達(dá)到較好的預(yù)測效果。因此得到的分析結(jié)果,除了提供一個合理的模型外,還可以同時完成一些類似于主成分分析和典型相關(guān)分析的信息[36-38]。
通過測量獲取了每個生育期48組LAI數(shù)據(jù)集,采用32組數(shù)據(jù)集作為建模集,剩余16組數(shù)據(jù)驗證模型效果。構(gòu)建本文的LAI估算模型,為了評估建模和驗證的擬合效果和預(yù)測精度,選用決定系數(shù)(Coefficient of determination,R2)、均方根誤差(Root mean square error,RMSE)作為評價指標(biāo)[23]。
為了提取出基于無人機(jī)高光譜遙感數(shù)據(jù)下的冬小麥株高,使用作物表面模型(CSM)來提取冬小麥株高。步驟如下:
(1)通過無人機(jī)高光譜影像的拼接處理,生成數(shù)字表面模型(DSM),在影像的不同小區(qū)共選取500個土壤樣本點,利用ArcGIS軟件的ArcTooLbox工具,基于DSM獲取土壤樣本點的高程。
(2)利用ArcGIS軟件的kriging工具,基于獲取的土壤樣本點的高程生成數(shù)字高程模型(DEM)。
(3)將DSM減去土壤樣本點的DEM,并使用ArcGIS軟件的grid calculator工具,提取出CSM。
(4)最后,利用ArcGIS軟件的ROI工具計算出生育期試驗小區(qū)的冬小麥株高。
為了驗證提取株高的效果,將獲取的數(shù)據(jù)與實測株高進(jìn)行分析,如圖2所示。由圖2可知,實測冬小麥高度和基于無人機(jī)高光譜提取的冬小麥高度(Hcsm)之間有較高的擬合性,R2為0.95,說明提取得到的Hcsm有較高的預(yù)測效果。因此,用提取的Hcsm進(jìn)行LAI估算研究。
圖2 實測和預(yù)測的冬小麥高度Fig.2 Measured and predicted winter wheat height
將選取的植被指數(shù)以及基于高光譜遙感提取的Hcsm分別與不同生育期的LAI進(jìn)行相關(guān)性分析,得到相關(guān)系數(shù)絕對值如表1所示。由表1可知,大部分植被指數(shù)與LAI表現(xiàn)為0.01顯著水平,相關(guān)性較好。拔節(jié)期,MCARI、TCARI和TVI表現(xiàn)為無顯著相關(guān),這3個植被指數(shù)與LAI相關(guān)性較低,剩余植被指數(shù)均為0.01顯著水平,其中相關(guān)系數(shù)絕對值最高的為0.713,對應(yīng)的植被指數(shù)為PBI。挑旗期,僅MCARI為無顯著相關(guān),其余植被指數(shù)為極顯著相關(guān)(0.01顯著水平),LCI的相關(guān)系數(shù)絕對值最高,為0.736。開花期,相比前2個生育期而言,大部分植被指數(shù)與LAI的相關(guān)系數(shù)絕對值有所增加,表現(xiàn)出較高的相關(guān)性,其中除TCARI外,均達(dá)到0.01顯著水平,相關(guān)系數(shù)絕對值最高達(dá)到0.828,為植被指數(shù)PBI。對于提取的Hcsm,從拔節(jié)期到開花期,隨著冬小麥生育期推移,Hcsm與LAI的相關(guān)性表現(xiàn)一直增強(qiáng),相關(guān)系數(shù)絕對值變大,在開花期達(dá)到最高的相關(guān)性,為0.585。植被指數(shù)和Hcsm與LAI的相關(guān)性在不同生育期表現(xiàn)有差異,說明相關(guān)性受生育期影響,整體上,植被指數(shù)、Hcsm與LAI的相關(guān)性呈現(xiàn)增強(qiáng)趨勢。
表1 植被指數(shù)、Hcsm與LAI的相關(guān)系數(shù)絕對值Tab.1 Absolute value of correlation coefficient of vegetation indices, Hcsm and LAI
根據(jù)表1中10種植被指數(shù)與LAI的相關(guān)性,發(fā)現(xiàn)3個生育期相關(guān)性最強(qiáng)的分別為PBI、LCI、PBI,為了得到不同生育期的最優(yōu)植被指數(shù)模型,分別構(gòu)建基于PBI、LCI和PBI的線性回歸模型,同時也構(gòu)建了不同生育期的基于Hcsm的線性回歸模型,如表2和圖3所示。根據(jù)表2和圖3可知,對于最優(yōu)植被指數(shù),從拔節(jié)期到開花期,建模R2的范圍是0.55~0.65,RMSE的變化為0.62~0.72,驗證R2的范圍是0.62~0.77,RMSE的變化為0.78~0.58,這3個生育期表現(xiàn)出最優(yōu)植被指數(shù)估算LAI效果逐漸增強(qiáng)。對于Hcsm,從拔節(jié)期到挑旗期,建模R2的變化范圍為0.25~0.26,RMSE為0.80~1.43,驗證R2的變化范圍為0.23~0.61,RMSE為0.88~1.23,建模和驗證結(jié)果都顯示出最優(yōu)植被指數(shù)估算LAI能力越來越強(qiáng)。從挑旗期到開花期,建模R2和RMSE的范圍分別為0.26~0.30和1.43~1.02,驗證R2和RMSE的范圍分別為0.61~0.39和1.23~0.85,隨著生育期建模效果逐漸增強(qiáng),驗證效果降低,其中建模R2增加的幅度低于驗證降低的幅度,因此從挑旗期到開花期,估算LAI效果逐漸降低。
表2 不同生育期最優(yōu)植被指數(shù)、Hcsm與LAI的回歸關(guān)系Tab.2 Regression relationships between optimal vegetation index, Hcsm and LAI in different growth stages
圖3 基于最優(yōu)植被指數(shù)和Hcsm的LAI實測值與預(yù)測值Fig.3 Measured and predicted values of LAI based on optimal vegetation index and Hcsm
為了探究Hcsm對植被指數(shù)估算LAI的精度影響,將每個生育期的最優(yōu)植被指數(shù)結(jié)合Hcsm,得到3個生育期最優(yōu)植被指數(shù)結(jié)合Hcsm與LAI的關(guān)系,如表3和圖4所示。從拔節(jié)期到開花期,最優(yōu)植被指數(shù)結(jié)合Hcsm建模R2一直增加(0.57~0.69),RMSE為0.61~0.67,驗證結(jié)果和建模保持一致,驗證R2也一直變大(0.63~0.79),RMSE從0.77降為0.54,說明建模效果較好。對比僅基于最優(yōu)植被指數(shù)和Hcsm構(gòu)建的LAI估算模型,將最優(yōu)植被指數(shù)結(jié)合Hcsm構(gòu)建的估算LAI模型效果更佳,擬合性和精度更高。
由表1可知,在拔節(jié)期中MCARI、TCARI和TVI表現(xiàn)為無顯著相關(guān),同時在挑旗期和開花期中這3種植被指數(shù)與LAI之間表現(xiàn)出較低的相關(guān)性。因此為了構(gòu)建LAI估算模型,選取剩余的7種植被指數(shù)進(jìn)行模型構(gòu)建,即模型輸入因子為LCI、NPCI、PBI、BGI、OSAVI、NDVI、SR,并使用PLSR回歸方法構(gòu)建不同生育期的LAI估算模型。為了探究植被指數(shù)結(jié)合Hcsm估算LAI能力,將選取的植被指數(shù)與Hcsm共同作為模型因子,也使用PLSR方法建立LAI估算模型,如表4、5和圖5所示。
表3 不同生育期最優(yōu)植被指數(shù)結(jié)合Hcsm與LAI的回歸關(guān)系Tab.3 Regression relationships between optimal vegetation index combined with Hcsm and LAI in different growth stages
圖4 基于最優(yōu)植被指數(shù)結(jié)合Hcsm的LAI實測值與預(yù)測值Fig.4 Measured and predicted values of LAI based on optimal vegetation index combined with Hcsm
從表4、5和圖5可以看出,以植被指數(shù)為模型因子,拔節(jié)期到開花期,建模R2呈現(xiàn)上升趨勢,驗證R2和建模R2變化保持一致,也是逐漸增加,在開花期達(dá)到最佳估算效果(建模R2=0.70,RMSE為0.67;驗證R2=0.79,RMSE為0.53)。以植被指數(shù)結(jié)合Hcsm為模型因子,3個生育期中隨著冬小麥生長,建模和驗證R2均表現(xiàn)出一直增加,開花期的估算效果最佳(建模R2=0.73,RMSE為0.64;驗證R2=0.83,RMSE為0.49)。將植被指數(shù)結(jié)合Hcsm構(gòu)建LAI估算模型,3個生育期建模和驗證的R2均高于以植被指數(shù)為因子的建模和驗證R2,Hcsm比較明顯地提高了LAI估算精度,構(gòu)建的估算模型擬合性和穩(wěn)定性都較好。
表4 基于植被指數(shù)的冬小麥不同生育期LAI估算模型評估參數(shù)Tab.4 Estimation of LAI in different growth stages of winter wheat based on vegetation indices
表5 基于植被指數(shù)結(jié)合Hcsm的冬小麥不同生育期LAI估算模型評估參數(shù)Tab.5 Estimation of LAI in different growth stages of winter wheat based on vegetation indices combined with Hcsm
圖5 基于植被指數(shù)以及植被指數(shù)結(jié)合Hcsm的LAI實測值與預(yù)測值Fig.5 Measured and predicted values of LAI based on vegetation indices and vegetation indices combined with Hcsm
為了得到LAI的空間分布,使用每個生育期的基于植被指數(shù)結(jié)合Hcsm的最佳估算模型,生成各生育期的LAI空間分布圖,如圖6所示。從圖6可知,從拔節(jié)期到挑旗期,LAI逐漸增大,而從挑旗期到開花期,LAI開始下降,這個表現(xiàn)和冬小麥實際LAI分布是一致的。各小區(qū)的LAI有明顯差別,其中重復(fù)2區(qū)域的LAI較高,重復(fù)1和重復(fù)3區(qū)域LAI較低,并且在3個生育期中均表現(xiàn)為重復(fù)2區(qū)域LAI高于重復(fù)1和重復(fù)3區(qū)域。拔節(jié)期LAI范圍為0~5,挑旗期LAI范圍為0~9,開花期LAI范圍為0~5,不同生育期的LAI空間分布和實際測量結(jié)果是一致的,說明構(gòu)建的不同生育期的最佳LAI估算模型是可靠的。
圖6 冬小麥各生育期的LAI空間分布Fig.6 Spatial distributions of LAI in different growth stages of winter wheat
傳統(tǒng)上測量作物株高主要是通過實地調(diào)查,這種測量方式消耗較多的人力、物力和時間。為了解決這個問題,高效準(zhǔn)確地監(jiān)測作物株高就變得尤為重要,可以提高農(nóng)業(yè)管理者的工作效率[39-40]。本研究通過冬小麥3個生育期的無人機(jī)高光譜影像數(shù)據(jù),提取出了不同生育期的Hcsm,實測Hcsm和Hcsm擬合的R2達(dá)到0.95,證明了提取的Hcsm精度較高。這個結(jié)果和文獻(xiàn)[41-42]的研究一致。然而,基于無人機(jī)高光譜影像提取的Hcsm相比實測的株高整體上偏小。主要是由于傳感器獲取的是作物的冠層數(shù)據(jù),冠層結(jié)構(gòu)的三維點云是由攝影測量計算的,包含不止第1片葉子,還含有植株的其余葉子和裸土像元;同時構(gòu)建冠層三維點云,可能會除去植株葉子的空間信息。另外,由于拔節(jié)期株高獲取時間和無人機(jī)影像獲取時間不完全同步,導(dǎo)致實測株高與預(yù)測株高精度降低,同步獲取株高和無人機(jī)影像數(shù)據(jù)可以提高株高監(jiān)測精度。
將不同的植被指數(shù)、Hcsm分別與LAI相關(guān)性分析,結(jié)果表明隨著生育期推移,大部分植被指數(shù)與LAI的相關(guān)性逐漸增強(qiáng),主要由于作物的生育期與光譜信息的敏感性相關(guān),生育期后期植被指數(shù)與LAI敏感性更強(qiáng);Hcsm與LAI整個生育期相關(guān)性都很好,說明作物高度能夠反映生長狀況[43],Hcsm隨著生育期增加,從而影響LAI的準(zhǔn)確估算。拔節(jié)期、挑旗期和開花期的最優(yōu)植被指數(shù)分別為PBI、LCI、PBI,這與文獻(xiàn)[44]發(fā)現(xiàn)PBI能夠很好地估算LAI以及文獻(xiàn)[23]探究出LCI在預(yù)測LAI方面有較高的能力結(jié)論一致。將最優(yōu)植被指數(shù)結(jié)合Hcsm構(gòu)建的LAI估算模型精度更高,效果更好,證明了Hcsm能夠提高估算LAI能力。
為了探究PLSR對LAI估算模型的效果,選取了LCI、NPCI、PBI、BGI、OSAVI、NDVI、SR作為模型的輸入變量,利用PLSR估算LAI。相比僅基于植被指數(shù)估算LAI,PLSR可以明顯地提高LAI估算效果。而基于植被指數(shù)結(jié)合Hcsm,并使用PLSR構(gòu)建LAI估算模型,模型的R2進(jìn)一步變大,RMSE變小,此時估算LAI的效果更好。2種模型均說明了使用PLSR能夠提高模型精度。這與PLSR處理光譜數(shù)據(jù)能力有關(guān),PLSR利用了主成分分析,更好地使用了光譜信息,并且能夠解決多個變量之間的共線性問題,從而在作物參數(shù)監(jiān)測方面能夠取得較好的成效[45]。同時文獻(xiàn)[23]證明了使用PLSR結(jié)合植被指數(shù)估算作物參數(shù)精度優(yōu)于單個植被指數(shù),這個研究和本文得到的結(jié)果是一致的。
(1)由無人機(jī)高光譜提取得到的Hcsm與實測Hcsm之間具有較高的擬合性(R2=0.95),Hcsm有較高的精度。
(2)在不同生育期,大部分植被指數(shù)與LAI表現(xiàn)為0.01顯著相關(guān)水平,Hcsm與LAI也表現(xiàn)為0.01顯著相關(guān)水平。
(3)不同生育期的最優(yōu)植被指數(shù)(PBI、LCI、PBI)估算LAI效果逐漸增強(qiáng),同時Hcsm以及最優(yōu)植被指數(shù)結(jié)合Hcsm估算LAI的能力也越來越強(qiáng),其中基于最優(yōu)植被指數(shù)結(jié)合Hcsm估算LAI效果優(yōu)于僅基于最優(yōu)植被指數(shù)或Hcsm的LAI估算效果。
(4)以植被指數(shù)、植被指數(shù)結(jié)合Hcsm為模型輸入因子,采用PLSR方法構(gòu)建的不同生育期LAI估算模型均在開花期達(dá)到最佳估算效果(建模R2=0.70,RMSE為0.67;建模R2=0.73,RMSE為0.64),并且以植被指數(shù)結(jié)合Hcsm估算LAI的精度最佳。
(5)通過無人機(jī)高光譜遙感影像提取作物Hcsm、并結(jié)合植被指數(shù),明顯提高了LAI估算精度,對農(nóng)業(yè)生產(chǎn)管理具有重要意義。