武倩雯，熊黑鋼,2，靳彥華，王莉鋒，王凱龍

(1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院，教育部綠洲生態(tài)重點實驗室， 新疆 烏魯木齊 830046；

2.北京聯(lián)合大學(xué)應(yīng)用文理學(xué)院， 北京 100083)

?

基于多個高光譜參數(shù)的玉米葉片葉綠素含量估測模型

武倩雯1，熊黑鋼1,2，靳彥華1，王莉鋒1，王凱龍1

(1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院，教育部綠洲生態(tài)重點實驗室， 新疆 烏魯木齊 830046；

2.北京聯(lián)合大學(xué)應(yīng)用文理學(xué)院， 北京 100083)

摘要：采用FieldSpecPro3光譜儀和SPAD-502葉綠素儀分別測定玉米葉片的光譜與其相對應(yīng)的葉綠素含量，通過分析紅邊位置、藍邊位置以及綠峰位置等高光譜參數(shù)與葉綠素含量的關(guān)系，建立葉綠素含量的單、雙和多變量光譜預(yù)測模型。結(jié)果表明： 在可見光區(qū)域，玉米葉綠素含量高，光譜反射率低，而進入近紅外區(qū)則剛好相反，葉綠素含量高，光譜反射率高； 紅邊位置、綠峰位置及藍邊位置各高光譜參數(shù)與葉綠素含量的相關(guān)性均達極顯著。其中紅邊位置與葉綠素含量的相關(guān)性最高，相關(guān)系數(shù)達0.84； 利用所選的3個高光譜參數(shù)分別建立的單、雙以及三變量模型，雖然大多數(shù)模型的精度R2大于0.71，但分析對比得出利用紅邊、藍邊及綠峰位置3個變量建立的模型具有最大模型精度R2、最小標(biāo)準(zhǔn)誤差(S)和均方根誤差(RMSE)，因此其模型預(yù)測能力較優(yōu)。

關(guān)鍵詞：玉米葉片；葉綠素含量；估測模型

葉綠素含量是作物生長過程中一個重要的生理指標(biāo)，對植被光合能力、發(fā)育階段以及營養(yǎng)狀況有重要的指示作用[1]，實時掌握玉米葉綠素含量及冠層光譜的動態(tài)變化，可以對玉米群體光合能力及營養(yǎng)狀況進行有效監(jiān)測，在玉米生長與產(chǎn)量預(yù)測以及精確診斷與動態(tài)調(diào)控中具有重要意義。

國內(nèi)外學(xué)者對建立作物葉綠素含量的最佳預(yù)測模型方面已有很多研究，并取得了一系列成果：通過對不同植被的光譜與葉綠素含量的研究，從而尋找它們各自的敏感波段[2]；分析作物不同生長期葉綠素含量與光譜的關(guān)系，建立不同生長階段的最適模型[3]；采用不同指數(shù)對葉綠素含量進行建模并且檢驗其模型精度[4]；利用反射光譜數(shù)據(jù)，特別是紅邊參數(shù)估計作物葉綠素[5]、氮素營養(yǎng)狀況[6]、葉面積指數(shù)[7]等；玉米、水稻的“紅邊”位置能很好地反映葉片葉綠素含量的信息[8]；通過研究植被光譜與葉綠素濃度的關(guān)系提出了光譜紅邊位置在植被葉綠素濃度估計中的重要作用[9]等。

這些研究大都是基于“紅邊”這一單變量建立作物葉綠素含量預(yù)測模型，而作物光譜的“峰谷”吸收特征是由葉片中相對應(yīng)的生化物質(zhì)引起的，這些化學(xué)成分是混合在一起的，彼此間加強或削弱了各自的吸收特征。所以，估測作物的某一生化成分時只用單一變量是不夠的，需要進行變量選擇和重組。由于玉米葉綠素含量不僅與紅邊位置相關(guān)性好[8]，與綠峰以及藍邊位置也具有較好的相關(guān)性[10]，并且這3個參數(shù)都是基于可見光波段，而此波段的反射率主要受葉片葉綠素含量的影響，因此它們能更好地反應(yīng)葉綠素含量信息。本文在前人研究的基礎(chǔ)上，分析多個玉米高光譜參數(shù)與其葉綠素含量之間的關(guān)系，建立預(yù)測模型。通過對比分析，選擇最優(yōu)的玉米葉片葉綠素含量預(yù)測模型。

1材料與方法

1.1研究區(qū)概況

研究區(qū)位于奇臺縣(43°25′～45°29′N，89°13′～91°22′E)，地處新疆維吾爾自治區(qū)東北部，天山山脈東段博格達山北麓，準(zhǔn)噶爾盆地東南緣。東邊與木壘哈薩克自治縣為鄰，西邊相接吉木薩爾縣，南連吐魯番市、鄯善縣，西北與富蘊縣、青河縣相交，東北與蒙古國接壤。南北長250 km，東西寬約45～150 km，總面積為1.81×104km2，邊界線長達131.47 km，是新疆的邊境縣之一。年平均氣溫5.5℃，7月份最高氣溫43℃，1月份極端最低氣溫-42.6℃。年平均降水176 mm，蒸發(fā)潛力2 141 mm，無霜期平均156 d，年平均風(fēng)速2.9 m·s-1。夏季炎熱，冬季寒冷， 四季分明，屬于典型的中溫帶大陸性干旱半干旱氣候。

2012年7月在平原區(qū)選擇3塊依靠井水灌溉的水澆地玉米，并用GPS進行現(xiàn)場定位。玉米行距50 cm，苗高大約2 m，葉寬大約10 cm，整體長勢較均勻，覆蓋度達95%～100%。品種是新玉23號，于4月中旬播種，9月底至10月初收獲。4月25日左右澆第一次水，在整個生長期共澆水6～7次，施肥1～2次，主要是底肥、農(nóng)家肥、尿素以及二胺分別約30 kg·667m-2。

1.2冠層光譜的采集及預(yù)處理

采用美國ASD公司的FieldSpecPro3光譜儀于2012年7月采集玉米光譜反射率數(shù)據(jù),其測量的波長范圍為350～2 500 nm。去除受外界噪聲影響較大的波段1 000～2 500 nm，取350～1 000 nm。這一波段包括了遙感常用的可見光和近紅外波段，可滿足本文分析要求。試驗選擇晴朗無云或少云的天氣，光譜儀探頭垂直向下，距玉米冠層頂部垂直高度約15 cm，視場角25°，每次采集玉米光譜前后進行白板校正以去除暗電流影響。選擇65個點進行光譜測定，每點重復(fù)測量5次取平均值作為該點光譜測量的結(jié)果。

采用移動平均法對實測光譜進行去噪處理[11]，即選取測定樣本某一點前后光譜曲線上一定范圍測定它的平均值。作為該點的值。公式：

(1)

1.3葉綠素含量的測定

日本Minolta Camera公司生產(chǎn)的SPAD-502葉綠素儀通過測量葉片在兩種波長范圍內(nèi)的透光系數(shù)來確定葉片葉綠素含量或 “綠色度”。測量光譜的同時，采用對植物無破壞性的SPAD-502葉綠素儀測定SPAD值(Soil plant analysis development)，以此代表玉米葉綠素含量[12-14]。在光譜儀探頭所能輻射的范圍內(nèi)，選取不同方向的玉米葉片各5片，每片葉隨機選取10個點測定其葉綠素含量，最終取5片葉的平均值。

1.4光譜參數(shù)選取

選取紅邊位置、藍邊位置和綠峰位置3個光譜參數(shù)，它們是描述植被色素狀態(tài)和健康狀況的重要指示波段。紅邊位置是在680～740 nm波長范圍內(nèi)的最大一階微分值對應(yīng)波長。當(dāng)綠色植被葉綠素含量高、生長旺盛時，“紅邊位置”會向長波方向偏移，即紅移，反之，當(dāng)綠色植被衰老或感染病蟲害，葉綠素含量減少時，“紅邊位置”將向短波方向偏移，即“藍移”[15-16]。490～530 nm內(nèi)的最大一階微分值對應(yīng)的波長為藍邊位置。在500～600 nm區(qū)間原始光譜出現(xiàn)峰值，將在此范圍內(nèi)光譜反射率的峰值對應(yīng)的波長作為綠峰位置[10]。

1.5精度評價

將總樣本隨機劃分為兩組，其中一組用于建立玉米葉片葉綠素含量預(yù)測模型，剩余一組用于模型精度檢驗。用模型精度(R2)、標(biāo)準(zhǔn)誤差(S)和均方根誤差(RMSE)對預(yù)測模型進行精度評價[17]。計算公式為：

(2)

(3)

(4)

2結(jié)果與分析

2.1不同葉綠素含量的玉米冠層光譜特征

將用于建模樣本的光譜數(shù)據(jù)按葉綠素含量大小分為三組，得到不同葉綠素含量水平的玉米葉片反射光譜變化曲線(圖1)?？傮w來說，玉米葉片光譜曲線具有一般健康植被光譜的“峰和谷”特征。在450～500 nm與650～700 nm之間，藍、紅光波段形成兩個低反射谷，由于在此波段范圍內(nèi)光輻射被葉片中的葉綠素全部吸收進行光合作用，因而具有較低的反射率。從500 nm起葉片的吸收減少,形成一個小的反射峰。從680～740 nm，反射率有一個極大的增加過程，在740 nm以后曲線趨于平緩，成為高反射平臺區(qū)。在可見光波段，不同葉綠素含量的光譜反射率表現(xiàn)出的總體特點是：SPAD 35>SPAD 40>SPAD 46,而在近紅外區(qū)域，其呈現(xiàn)出的特點與可見光波段完全相反，即SPAD 46>SPAD 40>SPAD 35。這是因為在可見光波段，玉米葉片反射率主要取決于葉綠素含量的多少。葉綠素含量多，吸收率高，反射率就低。而近紅外區(qū)域玉米葉片光譜反射率一般受葉片內(nèi)部細胞結(jié)構(gòu)的影響，葉綠素含量高的葉片，其內(nèi)部細胞更為復(fù)雜，因而反射率高。

2.2玉米葉綠素含量與其紅邊、綠峰和藍邊位置的相關(guān)性

玉米紅邊、藍邊以及綠峰位置分別與葉綠素含量的相關(guān)分析表明，它們與葉綠素含量具有良好的相關(guān)性且均達極顯著水平(P<0.01)(表1)。其中，紅邊位置和藍邊位置分別與葉綠素含量呈極顯著正相關(guān)，相關(guān)系數(shù)分別為0.84和0.61；而綠峰位置與葉綠素含量呈極顯著負相關(guān)，r=-0.75。因此，可用這些高光譜參數(shù)建立玉米葉綠素含量估測模型。

注：** 極顯著水平(P<0.01)；* 顯著水平(P<0.05)

Note:“**”shows significant correlation atP<0.01; “*” significant correlation atP<0.05.

2.3玉米葉綠素含量估測模型的建立

利用回歸分析法，以葉綠素含量為因變量，分別以紅邊位置、藍邊位置、綠峰位置、紅邊和藍邊位置、紅邊和綠峰位置、藍邊和綠峰位置以及紅邊、藍邊和綠峰位置為自變量建立玉米葉綠素含量的預(yù)測模型(曲線模型R2與一次線性模型相當(dāng)，因此，只做一次線性模型)。

依據(jù)模型精確度R2最大和標(biāo)準(zhǔn)誤差最小的原則，可以看出基于3個變量建立的模型預(yù)測能力最好，其次是利用紅邊與綠峰位置、紅邊與藍邊位置、藍邊與綠峰位置建立的雙變量模型，最后是分別基于紅邊位置、綠峰以及藍邊位置建立的單變量模型(表2)。

2.4精度評價

利用25個樣本對玉米葉綠素含量估測模型進行精度評價。找出建模樣本的紅邊位置、藍邊位置和綠峰位置，分別將其代入不同的預(yù)測模型計算葉片葉綠素含量，并將結(jié)果與實際測定的玉米葉片葉綠素含量進行對比，以此來檢驗?zāi)Ｐ途?圖2)。

通過檢驗可知：在單變量預(yù)測模型中，由于紅邊位置與葉綠素含量關(guān)系密切，作物在生長過程中隨著葉綠素含量的升高，紅邊位置則向長波方向移動，反之亦然；因此，其模型精度不僅最高，而且不論是R2還是S均接近雙變量模型。雙變量模型預(yù)測能力優(yōu)于單變量模型。利用3個變量(紅邊、藍邊和綠峰位置)建立的預(yù)測模型最好，其模型精度R2=0.81最大，標(biāo)準(zhǔn)誤差和均方根誤差均最小(S=2.348、RMSE=1.973)，預(yù)測值與實測值的相關(guān)性達0.93。這是因為作物葉綠素含量不僅僅反映在紅邊位置，其受多種因素的共同影響，而紅邊位置、藍邊位置以及綠峰位置是描述植物色素狀態(tài)和健康狀況的重要指示波段[1,15-18]，三者結(jié)合能更好地反映作物的葉綠素含量情況，所建模型預(yù)測能力更高。

表2　玉米葉綠素含量估測模型

圖2不同模型實測值與預(yù)測值之間的關(guān)系

Fig.2Relationships between the measured and predicted values from different models

注：(a)紅邊位置；(b)藍邊位置；(c)綠峰位置；(d)紅邊位置和藍邊位置；(e)紅邊位置和綠峰位置；(f)藍邊位置和綠峰位置；(g)紅邊位置、綠峰位置和藍邊位置

Note: (a) red edge; (b) blue edge; (c) green peak; (d) red edge and blue edge; (e) red edge and green peak; (f) blue edge and green peak; (d) red edge, blue edge and green peak.

3討論

葉綠素含量是植物營養(yǎng)脅迫、光合作用能力和生長狀況的良好指示指標(biāo)，直接控制著植被能量和物質(zhì)交換。前人在建立玉米葉綠素含量的估測模型方面已經(jīng)做了很多工作。采用不同指數(shù)對葉綠素含量進行建模并且檢驗其模型精度，對不同色素含量與其光譜特征的關(guān)系進行研究[4,19]；唐延林等[20]在玉米高光譜與紅邊特征的分析中，利用紅邊參數(shù)建立葉綠素含量預(yù)測模型等；吉海彥等[21]利用偏最小二乘方法建立了冬小麥葉片葉綠素和水分含量與反射光譜的定量分析模型等。但大多數(shù)研究都是經(jīng)過相關(guān)分析后選擇相關(guān)性較大的某一個指數(shù)或參數(shù)建立葉綠素含量的預(yù)測模型。由于作物的葉綠素含量受多種因素影響，所以只選用一種參數(shù)是遠遠不夠的。本文是在前人研究的基礎(chǔ)上選擇與葉綠素含量相關(guān)性較大的3個參數(shù)，分別基于它們建立葉綠素含量的多變量預(yù)測模型。本試驗的數(shù)據(jù)樣本是在干旱地區(qū)采集的，由于樣本有限及環(huán)境條件不同，所建模型能否適應(yīng)其他地區(qū)及不同地類種植的玉米還需進一步研究。

4結(jié)論

本文測定了玉米葉片光譜及其對應(yīng)的葉綠素含量，分析玉米葉片光譜紅邊、藍邊以及綠峰位置3個變量分別與葉綠素含量的相關(guān)性，并建立玉米葉片葉綠素含量監(jiān)測的單變量、雙變量和三變量模型。通過分析對比得出：利用紅邊、藍邊及綠峰位置3個變量建立的玉米葉片葉綠素含量預(yù)測模型最優(yōu)，其模型預(yù)測值與實測值相關(guān)系數(shù)達0.93,均方根誤差RMSE最小，為1.973。其次是雙變量模型，最后是單變量模型。

參 考 文 獻：

[1]Co Minolta. Chlorophyll SPAD-502 Instruction Manual[Z]. Radiometeric Instruments Operations,1989:17-21.

[2]劉克,趙文吉,郭逍宇,等.野鴨湖典型濕地植物光譜特征[J].生態(tài)學(xué)報,2010,30(21):5853-5861.

[3]王凱龍,熊黑鋼,張芳.干旱區(qū)冬小麥不同生長階段的光譜特征與葉綠素含量估測研究[J].干旱區(qū)資源與環(huán)境,2013,27(11):45-49.

[4]楊杰,田永超,等.水稻上部葉片葉綠素含量的高光譜估算模型[J].生態(tài)學(xué)報,2009,12(29):6561-6571.

[5]袁杰,王登偉,黃春燕,等.基于高光譜數(shù)據(jù)的棉花葉綠素密度定量提取研究[J].干旱地區(qū)農(nóng)業(yè)研究,2007,25(3):89-93.

[6]張雪紅,劉紹民,何蓓蓓.不同氮素水平下油菜高光譜特征分析[J].北京師范大學(xué)學(xué)報(自然科學(xué)版),2007,43(3):245-249.

[7]黃敬峰,王淵,王福民,等.油菜花紅邊特征及其葉面積質(zhì)素的高光譜估算模型[J].農(nóng)業(yè)工程學(xué)報,2006,22(8):22-26.

[8]吳長山,童慶禧,鄭蘭芬,等.水稻、玉米的光譜數(shù)據(jù)與葉綠素的相關(guān)分析[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報,2000,8(1):31-37.

[9]Horler D N H, Barber J, Barringer A R. Effects of heavy metals on the absorbance and reflectance spectra of plants[J]. International Journal of Remote Sensing，1980,1(2):121-136.

[10]錢育蓉,楊峰,李建龍,等.利用高光譜數(shù)據(jù)快速估算高羊茅牧草光合色素的研究[J].草業(yè)學(xué)報,2009，18(4):94-102.

[11]何挺,王靜.土壤氧化鐵光譜特征研究[J].地理與地理信息科學(xué),2006,3(2):30-34.

[12]李敏夏,張林森,李丙智,等.蘋果葉片高光譜特性與葉綠素含量和SPAD值的關(guān)系[J].西北林學(xué)院學(xué)報,2010,25(2):35-39.

[13]朱新開,盛海君,顧晶,等.應(yīng)用SPAD值預(yù)測小麥葉片葉綠素和氮含量的初步研究[J].麥類作物學(xué)報,2005,25(2):46-50.

[14]盧曉萍,楊丙賢,徐嬋娟,等.3種小檗科植物葉片SPAD值與葉綠素的相關(guān)性及通徑分析[J].浙江大學(xué)學(xué)報(農(nóng)業(yè)與生命科學(xué)版),2012,39(3):261-266.

[15]張金恒,王珂,王人潮.高光譜評價植被葉綠素含量的研究進展[J].上海交通大學(xué)學(xué)報(農(nóng)業(yè)科學(xué)版),2003,21(1):742-780.

[16]張永賀,郭嘯用,褚武道.等.基于紅邊位置的木荷葉片葉綠素含量估測模型研究[J].紅外與激光工程,2013,42(3):798-804.

[17]潘蓓,趙庚星,朱西存,等.利用高光譜植被指數(shù)估測蘋果樹冠層葉綠素含量[J].光譜學(xué)光譜分析,2013,33(8):2203-2206.

[18]黃春燕,王登偉,張煜星.基于棉花紅邊參數(shù)的葉綠素密度及葉面積指數(shù)的估算[J].農(nóng)業(yè)工程學(xué)報,2009,25(增刊2):137-141.

[19]張俊華,張佳寶.夏玉米光譜特征對其不同色素含量的響應(yīng)差異[J].西北農(nóng)業(yè)學(xué)報,2010,19(4):70-76.

[20]唐延林,王秀珍,等.玉米高光譜與紅邊特征的分析[J].山地農(nóng)業(yè)生物學(xué)報,2003,22(3):189-194.

[21]吉海彥,王鵬新,嚴(yán)泰來.冬小麥活體葉片葉綠素和水分含量與反射光譜的模型建立[J].光譜學(xué)與光譜分析,2007,27(3):514-516.

Prediction model on chlorophyll content in maize leaf based on

several high spectral parameters

WU Qian-wen1, XIONG Hei-gang1,2, JIN Yan-hua1, WANG Li-feng1, WANG Kai-long1

(1.CollegeofResources&EnvironmentScience,XinjiangUniversity,KeyLaboratoryof

OasisEcologyMinistryofEducation,Urumqi,Xinjiang830046,China;

2.CollegeofAppliedArts&Science,BeijingUnionUniversity,Beijing100083,China)

Abstract：In this research, FieldSpecPro3 spectroscopy and SPAD-502 chlorophyll meter were used to measure the spectrum of corn leaf chlorophyll content and its amount respectively. Through the analysis on the relationships between chlorophyll content and parameters of red, blue, and green edge positions and spectral peak positions, single, double and multivariate spectral prediction models about the chlorophyll content were established. The results showed that in the visible spectrum, the higher corn chlorophyll content, the lower reflectance spectrum. However, exact opposite was observed in the infrared spectrum. Spectral parameters of red edge, green peak and blue edge positions were in significant correlations with the chlorophyll content, reaching up to 0.84 between red edge position and the chlorophyll content. As a result, single, double and three-variable models using three high spectral parameters were established. Although the accuracy ofR2 were mostly greater than 0.71, it was further suggested that three-variable model had the best accuracy inR2, the minimal standard deviation (S) and root mean square error (RMSE), which might provide better prediction results than the other two models.

Keywords:maize leaf; chlorophyll content; prediction models

中圖分類號：S513

文獻標(biāo)志碼：A

通信作者：熊黑鋼(1956—)，男，湖南湘鄉(xiāng)人，博士，教授，博士生導(dǎo)師，主要從事資源與環(huán)境研究。 E-mail: heigang@sbuu.edu.cn。

作者簡介：武倩雯(1989—)，女，陜西寶雞人，碩士，主要研究方向為干旱區(qū)資源與環(huán)境。 E-mail:wqw200952199@163.com。

基金項目：國家自然科學(xué)基金(41171165)；北京市屬高等學(xué)校高層次人才引進與培養(yǎng)計劃項目(IDHT20130322)；北京聯(lián)合大學(xué)人才強校計劃資助項目(BPHR2012E01)

收稿日期：2014-12-10

doi:10.7606/j.issn.1000-7601.2016.01.31

文章編號：1000-7601(2016)01-0201-05