基于分數階微分的春小麥葉綠素含量高光譜估算

茹克亞·薩吾提，阿不都艾尼·阿不里，尼加提·卡斯木，李 虎，吾木提·艾山江，亞森江·喀哈爾，李曉航

(1.新疆大學資源與環境科學學院，新疆烏魯木齊 830046；2.新疆大學綠洲生態教育部重點實驗室，新疆烏魯木齊830046；3.新疆大學干旱生態環境研究所，新疆烏魯木齊 830046；4.滁州學院，安徽省高分辨率對地觀測系統數據產品與應用軟件研發中心，安徽省滁州市 239000)

葉綠素是綠色植物光合作用的重要色素，也是反映植被生長狀態的一個重要指標[1-2]。反射光譜與植物葉片葉綠素含量之間存在直接關系，依據此關系可利用遙感技術尤其是高光譜技術對農作物葉綠素含量進行高效、無損害、安全、快速估算，從而為農業區域生產和管理提供有效手段[3-4]。目前，小麥[5]、棉花[6]、玉米[7]、水稻[8]等不同農作物上已開展了相關研究。Croft等[9]利用Landsat-5、TM和其他多光譜傳感器數據運用物理建模方法來預測植物葉片葉綠素含量，并得到了較好的預測效果。Peng等[10]基于高光譜數據計算歸一化植被指數(normalized difference vegetation index，NDVI)、增強型植被指數(enhanced vegetation index，EVI)和樣本的比率(sample ratio,SR)，并對玉米與大豆葉綠素含量季節性變化進行了分析，發現初級生產量預測模型(gross primary productivity,GPP)在農作物葉片葉綠素含量預測方面具有一定潛力。劉 芬等[4]以水稻綠葉和黃葉為試驗材料，分析了高光譜指數與水稻葉綠素含量的關系，并發現用水稻綠葉RVI(363，675)和黃葉DVI(639，680)與葉綠素含量建立模型的估算精度最佳，表明利用高光譜技術結合一階微分光譜的方法可以監測水稻葉綠素含量。袁 杰等[11]研究表明，基于一階導數的反射率與葉綠素含量之間的關系，建立的線性模型的預測能力優于非線性模型。王 爍等[12]認為，RSI指數二次曲線回歸模型能夠很好地估算棉花葉綠素含量，且預測精度相對較高。

目前，在葉綠素含量高光譜估算模型研究中，主要利用原始高光譜、整數階微分、數學變換、植被指數計算等傳統方法[13-14]。雖然以上方法都能估算植被葉綠素含量，但是都有一定的局限性。利用高光譜數據直接對農作物葉綠素含量進行預測存在一定的難度。對于高光譜這類高維數據源，整數階微分變換會忽略其中的分數階微分信息，可能造成信息丟失，建模精度也會受到制約[15]。分數階微分在階數上對整數階微分進行了拓展，并在數字濾波、信號與圖像處理等領域已成熟運用[16]。鑒于此，本研究運用分數階微分處理方法，分析各階(0～2)微分光譜數據與春小麥抽穗期葉片葉綠素含量相關性變化趨勢，并選取相關性最強的分數階數，利用變量重要性投影法(VIP)篩選有效植被指數和運用偏最小二乘回歸(PLS)構建葉綠素含量估算模型，探討分數階微分預處理在小麥葉綠素含量高光譜估算中的可行性，以期為利用高光譜遙感反演春小麥葉綠素含量以及提高其估算精度提供參考。

1 材料與方法

1.1 研究區概況

阜康市位于新疆東天山山脈北坡，準噶爾盆地南部,地理坐標為87°46′～88°44′E，43°45′～45°30′N。境內氣候垂直地帶性分布特征顯著，具有新疆地區典型地理地貌特征(山地-綠洲-荒漠)，屬典型的溫帶大陸性荒漠氣候，四季分明，冬季嚴寒且夏季炎熱，降水稀少且分布不均[17]，年日照時數為2 932.7 h，平均氣溫13.5 ℃，年均降水量248.4 mm，年平均蒸發量1 553.2 mm，凍土深度202 mm，無霜期174 d。土地肥沃而廣闊，歷來為阜康的重要產糧區[18]。

研究區播種春小麥總面積為22 500 m2(150 m×150 m)，采樣區以1 m×1 m為樣方，并進行網格狀采樣，共165個采樣區(圖1)。為了減少測量誤差取均值進行計算。采樣時間為2017年6月初，正處于春小麥抽穗初期。

圖1 采樣區示意圖

1.2 高光譜數據采集及預處理

春小麥葉片光譜測量使用ASDFieldspecFR2500光譜儀，光譜范圍350～2 500 nm，其光譜分辨率分別為3 nm(350～1 000 nm)和10 nm(1 000～2 500 nm)，實際采樣間隔1 nm。在采樣區進行測量時，每個活體樣品由多個葉片組成，距植株葉片10 cm左右，測量時探頭的天頂角為15°；野外光譜測定時，各進行5次白板標定，重復測定10次，光譜曲線掃描設定時間為0.2 s，最后用其平均值作為該點的光譜反射值。測量光譜時間均在北京時間12:00-14:00，天氣晴朗，能見度好，無風或者風速很小的情況下進行。研究區獲取的野外高光譜數據利用View-Spectral軟件進行數據導出，為減少噪聲影響，剔除信噪比較低的邊緣波段(350～400 nm和2 401～2 500 nm)以及易受大氣中水汽影響的噪聲波段(1 800～ 2 000 nm和1 350～1 500 nm)后，利用Savitzky-Golay濾波對165個采樣區的401～2 400 nm 的反射光譜數據進行平滑去噪。

1.3 葉綠素含量測定

每個采樣區使用五點法取樣，其中每個采樣點選取3個植株，對倒一、倒二葉的葉綠素含量各測定3次，取五點算術平均值作為采樣區冠層葉綠素含量值，從試驗田的11個橫斷面的165個采樣區獲得最終測量結果。在實驗基地內的實驗室中剪取新鮮葉片0.2 g，分別加入少量無水乙醇及石英砂和碳酸鈣粉將其研磨成勻漿至組織發白，過濾后定容至25 mL，然后用UV2450-紫外可見分光光度計測定665 nm、649 nm波長下的消光度值，再依據公式計算葉綠素a和葉綠素b的濃度值以及春小麥葉片單位質量葉綠素含量[19]。將冠層葉綠素含量從高到低進行排序，等間隔抽取99個建模集樣本與66個驗證集樣本，由冠層葉綠素含量數據的主要描述性統計(表3)可知，數據離散程度不強，屬于弱變異強度。

1.4 分數階微分算法

分數階微分對整數階微分的概念進行了擴展，是專門研究任意階微分數學性質及應用的領域[20]。分數階微分已經成功應用于系統建模、信號濾波、模式識別、分形理論等方面，主要是因為在這些應用中，分數階模型比原有的整數階模型的描述結果更為精確[21-22]。數學階微分運算主要基于Grünwald-Letnikov(G-L)分數階微分一元函數差分表達式來實現。分數階微分的計算公式如下(具體算法程序采用Matlab-2014實現)：

(1)

式中，Γ為Gamma函數；N為數據長度；q為實數。

當q值為分數時，式(1)就給出分數導數結果。令δ=λ/N，整理為下式：

(2)

1.5 植被指數的計算

表1 植被指數計算公式表Table 1 Formulation of vegetation indices

1.6 變量重要性準則(VIP)

變量重要性準則(VIP)由Wold[25]首先提出。VIP值代表自變量對模型擬合的程度，VIP值小于1的自變量對因變量的貢獻較小，可以考慮剔除。因此，本研究自變量VIP的閾值為大于等于1。VIP值在XLSTAT 2015軟件中計算，其計算公式如下：

(6)

其中，Wjf是j變量和f分量的權重值，SSYf是f分量和J個變量的方差平方和，SSYtotal是解釋因變量的方差總平方和，F是分量總數。

1.7 估算模型建立與驗證

表2 模型評價范圍Table 2 Classification of performances of models in terms of r2 and RPD

2 結果與分析

2.1 春小麥葉綠素含量統計與光譜特征

建模集和驗證集中春小麥葉綠素含量的最大值分別為6.04和5.98 mg·g-1，最小值分別為3.51和2.87 mg·g-1，平均值分別為5.07和 5.0 mg·g-1，變異系數分別為9.3%和 10.6%；全部采樣區的平均值為5.04 mg·g-1，變異系數為9.7%，屬于中小等變異程度(表3)。

春小麥葉片不同葉綠素含量樣本的光譜曲線形態上基本相似(圖2a)。從光譜曲線看，在可見光波段(400～760 nm)、近紅波段(760～1 300 nm)范圍內具有明顯的反射峰和吸收谷。除 760～800 nm波段范圍外，不同葉綠素含量小麥的其余波段光譜反射率曲線較容易區分。

表3 采樣區小麥葉綠素含量統計特征Table 3 Statistical characteristics of the wheat chlorophyll of sampling area mg·g-1

a圖中彩色線代表不同葉綠素含量樣本。

Chromatic lines in chart a represent the samples with different chlorophyll contents.

圖2 春小麥葉綠素含量與光譜反射率變化曲線

Fig.2 Change curve of chlorophyll content and spectral reflectance of spring wheat

2.2 分數階微分對葉片光譜反射率顯著性的影響

2.3 高光譜植被指數選擇

圖3 相關系數通過0.01顯著性檢驗水平的波段數量

表4 基于4種變換在1.2階微分的植被指數與葉綠素含量相關性Table 4 Correlation coefficients between vegetation indicexes calculated by 1.2-order differential of transformed spectral reflectance and chlorophyll contents

**：P<0.01；*：P<0.05.

2.4 基于PLS模型估算春小麥葉綠素含量

表5 不同光植被指數反演葉綠素含量的PLSR模型Table 5 Chlorophyll content in PLSR models derived from spectral vegetation indexesindices

5個模型中，基于lgR的1.2階微分建立的模型精度最優。從5種模型的預測值和實測值散點圖(圖4)可以看出，基于lgR的1.2階微分建立的模型對春小麥葉綠素含量具有較好的定量反演能力，其RPD=2.17，r2=0.87，RMSE= 0.24mg·g-1。

圖4 基于有效植被指數的春小麥葉綠素含量估算模型精度檢驗

3 討 論