田 敏， 周 杰， 張 澤， 溫鵬飛， 呂 新

(1.石河子大學(xué)信息科學(xué)與技術(shù)學(xué)院，新疆石河子 832000；2.石河子大學(xué)農(nóng)學(xué)院/新疆生產(chǎn)建設(shè)兵團(tuán)綠洲生態(tài)農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室，新疆石河子 832000)

基于高光譜植被指數(shù)對(duì)棉花葉綠素含量的估算

田 敏1，2， 周 杰1， 張 澤2， 溫鵬飛2， 呂 新2

(1.石河子大學(xué)信息科學(xué)與技術(shù)學(xué)院，新疆石河子 832000；2.石河子大學(xué)農(nóng)學(xué)院/新疆生產(chǎn)建設(shè)兵團(tuán)綠洲生態(tài)農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室，新疆石河子 832000)

為了研究不同水氮組合條件下葉片葉綠素含量與光譜反射率之間的相關(guān)性，從而進(jìn)一步估算葉綠素含量，利用ASD Field SpecPro Fr型光譜儀室外測(cè)量棉花葉片的光譜反射率，同時(shí)采收棉花葉片獲得葉綠素含量值。計(jì)算光譜一階微分、原始光譜反射率組成的植被指數(shù)，以及對(duì)“三邊參數(shù)”與葉片葉綠素含量進(jìn)行相關(guān)性分析，結(jié)果表明：“三邊參數(shù)”中的紅邊內(nèi)最大一階微分值(Dr)與棉花葉綠素含量有很好的相關(guān)性，其決定系數(shù)r2為0.530 5**，估測(cè)模型決定系數(shù)r2為0.856 3**，均方根誤差RMSE為0.366；植被指數(shù)中BmSR705：(DR750-DR445)/(DR705-DR445)與棉花葉綠素含量有很好地相關(guān)性(r2=0.696 3**)，估測(cè)模型決定系數(shù)r2為0.815 7**，均方根誤差RMSE為0.278。因此利用特定的植被指數(shù)和“紅邊參數(shù)”能夠很好地預(yù)測(cè)葉綠素含量，從而為高光譜數(shù)據(jù)預(yù)測(cè)棉花葉片葉綠素含量提供理論基礎(chǔ)。

高光譜；植被指數(shù)；反射率；棉花；葉綠素含量；三邊參數(shù)

新疆是全國(guó)商品棉重要的生產(chǎn)基地，新疆獨(dú)特的機(jī)械化、規(guī)?；葪l件為高光譜分辨率遙感技術(shù)在棉花上的應(yīng)用提供了一個(gè)良好的應(yīng)用平臺(tái)[1]，而遙感技術(shù)的應(yīng)用為棉花生長(zhǎng)的快速監(jiān)測(cè)提供了便利。葉綠素含量與植被的光合能力、發(fā)育階段具有較好相關(guān)性，它通常是光合作用能力和植被發(fā)育階段(特別是衰老階段)的指示器[2]。當(dāng)前室內(nèi)化學(xué)分析葉綠素含量的方法既費(fèi)時(shí)、又費(fèi)力，使診斷得到的結(jié)果很難及時(shí)應(yīng)用到實(shí)際當(dāng)中。然而，高光譜遙感技術(shù)以其光譜分辨率高、波段連續(xù)性強(qiáng)、數(shù)據(jù)豐富、獲取快、無(wú)損監(jiān)測(cè)的特點(diǎn)，在獲取植被養(yǎng)分、葉綠素含量等研究中得到廣泛應(yīng)用[3]。本研究以棉花為研究對(duì)象，針對(duì)不同水氮處理?xiàng)l件下獲取棉花葉片葉綠素含量，針對(duì)不同棉花葉片的光譜敏感波段組合，構(gòu)建棉花葉片葉綠素光譜診斷模型，為快速獲取葉綠素含量數(shù)據(jù)提供依據(jù)。

1 材料與方法

1.1 研究區(qū)介紹

試驗(yàn)地位于新疆維吾爾自治區(qū)瑪納斯縣六戶地鎮(zhèn)創(chuàng)田地(86°07′36″E，44°39′15″N)，地處新疆腹地，屬溫帶大陸性氣候，冬季長(zhǎng)而嚴(yán)寒，夏季短而酷熱，晝夜溫差大，年平均氣溫7.2 ℃，全年無(wú)霜期170 d左右，年平均降水量170 mm。土壤質(zhì)地為中壤土，0～20 cm土層有機(jī)質(zhì)含量為19.09 g/kg，堿解氮含量為51.50 mg/kg，速效磷含量為19.89 mg/kg，速效鉀含量為212.88 mg/kg，田間持水量為30.1 g/m3，土壤容重為1.2 g/cm3。試驗(yàn)時(shí)間為2014—2015年，試驗(yàn)設(shè)4個(gè)氮(N)水平，即施氮量0(N0)、200(N200)、400(N400)、600(N600)kg/hm2，全部作為追肥；設(shè)3個(gè)水分水平，分別為3 000 (W3 000)、4 500(W4 500)、6 000(W6 000)m3/hm2。采用雙因素完全隨機(jī)分組，共12個(gè)處理，3次重復(fù)。供試棉花品種為新陸早50(XLZ 50)，種植模式為“1膜2管4行”，膜寬2.05 m，行距配置為30 cm+60 cm+30 cm，株距15 cm，種植密度 21.3 萬(wàn)株/hm2。試驗(yàn)田種植方式為棉花連作，南北走向，兩頭設(shè)置保護(hù)行。4月25日播種，并灌出苗水，6月份開(kāi)始第1次灌水，氮肥分7次隨水施肥，磷鉀肥各150 kg/hm2播前一次性施入。小區(qū)長(zhǎng)、寬分別為13、2.05 m?；瘜W(xué)調(diào)控與其他管理措施按高產(chǎn)栽培要求進(jìn)行。

1.2 高光譜數(shù)據(jù)獲取與校正

光譜數(shù)據(jù)采用美國(guó)ASD公司生產(chǎn)的ASD Field SpecPro Fr型光譜儀，光譜測(cè)量范圍為350～2 500 nm，其中350～1 000 nm 光譜采樣間隔為1 nm，在1 000～2 500 nm范圍內(nèi)光譜采樣間隔2 nm，重采樣間隔為1 nm[4]。2014—2015年分別在6月27日、7月10日、7月17日、8月5日、8月11日、8月29日選擇均勻且代表性的植株進(jìn)行棉花冠層光譜數(shù)據(jù)測(cè)量。同時(shí)采收葉片進(jìn)行室內(nèi)葉綠素含量測(cè)定。為了最大限度挖掘高光譜圖像數(shù)據(jù)中的光譜有效信息，去除噪音、樣品背景和雜散光影響等因素對(duì)模型預(yù)測(cè)性能的影響，需要對(duì)光譜數(shù)據(jù)進(jìn)行預(yù)處理[5]。經(jīng)過(guò)預(yù)處理后，將每個(gè)樣品測(cè)得的反射光譜曲線經(jīng)算術(shù)平均后則為該樣品的實(shí)際反射光譜數(shù)據(jù)。對(duì)光譜曲線作一階微分變換，獲得一階微分光譜數(shù)據(jù)。

1.3 葉片葉綠素含量測(cè)定

選取棉花第4完全展開(kāi)葉并剪下中間部分，剪碎混勻后稱取0.200 0 g，放入研缽中，加少量石英砂和碳酸鈣粉及2～3 mL 95%乙醇，研成勻漿，再加乙醇10 mL，繼續(xù)研磨至組織變白。然后將提取液倒入漏斗中，過(guò)濾到25 mL棕色容量瓶中，并用乙醇定容至25 mL，搖勻。提取液用可見(jiàn)光分光光度計(jì)分別測(cè)定645 nm和663 nm波長(zhǎng)處的吸光度，根據(jù)下列公式分別計(jì)算總?cè)~綠素值[6]。

葉綠素含量=C(20.29×D645 nm+8.04×D663 nm)×V/S。

式中：C為葉綠素濃度(mg/L)；V代表提取液體積；S是樣品的面積。

1.4 數(shù)據(jù)處理分析

采用Microsoft Excel 2003、SPSS 13.0和Origin 8.6分析數(shù)據(jù)和繪圖。

2 結(jié)果與分析

2.1 葉片葉綠素含量分析

2.1.1 葉片葉綠素含量統(tǒng)計(jì)分析 2014—2015年6個(gè)時(shí)期獲取的葉綠素值取2年的平均值，共計(jì)216個(gè)作為試驗(yàn)樣本。選取170個(gè)樣本組成高光譜數(shù)據(jù)模型，46個(gè)樣本作為驗(yàn)證數(shù)據(jù)集。棉花葉片整個(gè)生育時(shí)期葉綠素含量統(tǒng)計(jì)分析如表1所示。

表1 棉花葉綠素含量統(tǒng)計(jì)分析

2.1.2 不同水氮處理間葉片葉綠素含量分析 葉綠素含量的變化可能是由多種原因造成的，其他脅迫如缺氮、缺鐵以及病害等也能造成葉綠素含量的變化[7]。由圖1可知，不同水氮處理下棉花葉片葉綠素的含量總體上呈現(xiàn)波浪式分布，棉花出苗后50 d葉綠素含量較高，原因可能因?yàn)楝F(xiàn)蕾期時(shí)日照光線較強(qiáng)，長(zhǎng)波段有利于葉綠素a含量的增長(zhǎng)，從而葉綠素總量達(dá)到峰值；第2個(gè)峰值出現(xiàn)在棉花出苗后80 d左右(花鈴期)，這時(shí)正處于生殖生長(zhǎng)與營(yíng)養(yǎng)生長(zhǎng)的重要階段，光合作用強(qiáng)，承擔(dān)著光合作用吸收轉(zhuǎn)化的重要作用，因此葉綠素含量達(dá)到峰值。隨著生育期的發(fā)展葉片逐漸枯萎，養(yǎng)分減少，葉綠素含量降低。

2.2 不同水氮處理間“三邊參數(shù)”與葉片葉綠素含量估算模型

“三邊參數(shù)”能夠很好地反映綠色植被的光譜特征，對(duì)葉綠素含量變化較為敏感。利用已有的基于光譜反射率的“三邊參數(shù)”以及構(gòu)造的差值植被數(shù)據(jù)與歸一化植被指數(shù)，建立棉花葉片葉綠素含量與相關(guān)指數(shù)的相關(guān)性，比較估算效果。由表2結(jié)果分析，在紅邊(680～760 nm)內(nèi)最大一階微分值(Dr)、紅邊位置(WP_r)、紅邊內(nèi)一階微分的總和(SDr)相關(guān)性較高，r2分別為0.530 5、0.523 6、0.523 5，藍(lán)邊、黃邊所對(duì)應(yīng)的光譜位置與葉綠素之間決定系數(shù)較低。棉花葉片葉綠素含量與紅邊(680～760 nm)內(nèi)最大一階微分值(Dr)、紅邊位置(WP_r)、紅邊內(nèi)一階微分的總和(SDr)具有很好的相關(guān)性，因此利用這些光譜參數(shù)進(jìn)行棉花葉綠素含量估算具有很好的可行性。

挑選出相關(guān)性較好的模型進(jìn)行模型驗(yàn)證，通過(guò)均方根誤差(RMSE)，決定系數(shù)(r2)來(lái)驗(yàn)證模型的優(yōu)劣，由圖2可知，在紅邊(680～760 nm)內(nèi)最大一階微分值(Dr)、紅邊位置(WP_r)、紅邊內(nèi)一階微分的總和(SDr)對(duì)應(yīng)的決定系數(shù)(r2)分別為0.856 3、0.705 8、0.560 1，均方根誤差(RMSE)分別為0.366、0.531、0.524。可以看出紅邊內(nèi)一階微分最大值(Dr)擬合程度較好，模型精度較高，因此利用紅邊內(nèi)一階微分最大值(Dr)能夠預(yù)測(cè)棉花葉綠素含量。

表2 棉花葉片葉綠素含量與“三邊參數(shù)”相關(guān)性

2.3 葉綠素含量與光譜指數(shù)的相關(guān)估算模型

總結(jié)分析近年來(lái)葉綠素含量與光譜指數(shù)相關(guān)性較好的植被指數(shù)，為了準(zhǔn)確地得到植被指數(shù)以及與葉綠素含量相關(guān)性較高的植被指數(shù)形式，采用原始波段反射率值、一階微分反射率值進(jìn)行波段優(yōu)選情況下組合，組成一般比值植被指數(shù)、一階微分比值植被指數(shù)、歸一化比值植被指數(shù)、其他歸一化比值植被指數(shù)，從而分析植被指數(shù)與葉綠素含量的相關(guān)性。由表3可見(jiàn)，在選擇的10個(gè)植被指數(shù)中均與葉綠素含量呈現(xiàn)出極顯著正相關(guān)關(guān)系，決定系數(shù)(r2)均大于0.36。其中3個(gè)植被指數(shù)相關(guān)性較好，分別是：PSSRb：R800/R635，所對(duì)應(yīng)表達(dá)式為y=0.768 3x+5.026 9(r2=0.664 4*)；PSNDb：(R800-R635)/(R800+R635)，對(duì)應(yīng)表達(dá)式為y=33.852x-11.312(r2=0.617 9**)；BmSR705：(DR750-DR445)/(DR705-DR445)，相關(guān)關(guān)系表達(dá)式為y=29.869x-6.984 7(r2=0.696 3**)。說(shuō)明波段800 nm和635 nm組成的比值植被指數(shù)與歸一化植被指數(shù)能夠很好地預(yù)測(cè)葉片葉綠素含量，另外750、705、445 nm組成的基值校正簡(jiǎn)單比值指數(shù)能夠很好地預(yù)測(cè)葉片葉綠素含量，并且具有很好的相關(guān)性。

挑選出相關(guān)性較好的模型進(jìn)行模型驗(yàn)證，通過(guò)均方根誤差(RMSE)，決定系數(shù)(r2)來(lái)驗(yàn)證模型的優(yōu)劣，由圖3可知，PSSRb：R800/R635、PSNDb：(R800-R635)/(R800+R635)、BmSR705：(DR750-DR445)/(DR705-DR445) 預(yù)測(cè)值與實(shí)測(cè)值擬合程度較高，決定系數(shù)分別為0.760 5、0.746 0、0.815 7，均方根誤差(RMSE)分別為0.327、0.343、0.278。因此利用組成的植被指數(shù)能夠很好地預(yù)測(cè)棉花葉綠素含量。

表3 特定植被指數(shù)與棉花葉片葉綠素含量相關(guān)性

注：PSND為色素歸一化差值指數(shù)(Pigment Specific ND)；PSSR為色素簡(jiǎn)單比值指數(shù)(Pigment Specific SR)；SRPI為簡(jiǎn)單比值色素指數(shù)(SR Pigment Index)；BmND為基值校正歸一化差值指數(shù)(Base modified ND)；BmSR為基值校正簡(jiǎn)單比值指數(shù)(Base modified SR)；“**”表示相關(guān)極顯著。

3 結(jié)論與討論

“三邊參數(shù)”與棉花葉綠素含量具有很好的相關(guān)性，其中紅邊(680～760 nm)內(nèi)最大一階微分值(Dr)、紅邊位置(WP_r)、紅邊內(nèi)一階微分的總和(SDr)與棉花葉綠素含量具有很好的相關(guān)性，通過(guò)預(yù)測(cè)值與實(shí)測(cè)值可知，紅邊(680～760 nm)內(nèi)最大一階微分值(Dr)決定系數(shù)r2=0.856 3，RMSE=0.366，因此紅邊(680～760 nm)內(nèi)最大一階微分值(Dr)能夠很好地預(yù)測(cè)棉花葉綠素含量。

特定植被指數(shù)與葉綠素含量較高的相關(guān)性，選擇的10個(gè)植被指數(shù)中均與葉綠素含量呈現(xiàn)出極顯著正相關(guān)關(guān)系，決定系數(shù)(r2)均大于0.36。其中，BmSR705：(DR750-DR445)/(DR705-DR445)相關(guān)關(guān)系表達(dá)式為y=29.869x-6.984 7(r2=0.696 3**)。預(yù)測(cè)模型決定系數(shù)r2=0.815 7，RMSE=0.278，因此可以利用3個(gè)波段組成的基值校正簡(jiǎn)單比值指數(shù)，能夠很好地預(yù)測(cè)植被葉綠素含量。

最后，本研究應(yīng)用實(shí)測(cè)的光譜指數(shù)和葉綠素含量，雖然利用統(tǒng)計(jì)學(xué)的方法得到的模型具有一定的可靠性，但是試驗(yàn)得到的模型仍然具有一定的誤差。在下一步研究中應(yīng)當(dāng)結(jié)合不同品種、地域進(jìn)行試驗(yàn)，完善估測(cè)模型，充分發(fā)揮高光譜技術(shù)的優(yōu)越性。

[1]曹衛(wèi)彬，楊邦杰，裴志遠(yuǎn)，等. 新疆棉花遙感監(jiān)測(cè)系統(tǒng)建立的對(duì)策[J]. 石河子大學(xué)學(xué)報(bào)(自然科學(xué)版)，2003，7(2)：165-168.

[2]Minoha Co，Ltd. Chlorophyll SPAD-502 instruction manual[M]. Radiometric Instruments Operations，1989：17-21.

[3]劉 磊，沈潤(rùn)平，丁國(guó)香. 基于高光譜的土壤有機(jī)質(zhì)含量估算研究[J]. 光譜學(xué)與光譜分析，2011(3)：762-766.

[4]祁亞琴，呂 新，陳冠文，等. 基于高光譜數(shù)據(jù)提取棉花冠層特征信息的研究[J]. 棉花學(xué)報(bào)，2011，23(2)：167-171.

[5]褚小立. 化學(xué)計(jì)量學(xué)方法與分子光譜分析技術(shù)[M]. 北京：化學(xué)工業(yè)出版社，2011.

[6]李合生，孫 群，趙世杰，等. 植物生理生化實(shí)驗(yàn)原理和技術(shù)[M]. 北京：高等教育出版社，2000：131.

[7]林芬芳，丁曉東，付志鵬，等. 基于互信息理論的水稻磷素營(yíng)養(yǎng)高光譜診斷[J]. 光譜學(xué)與光譜分析，2009，29(9)：2467-2470.

[8]Vogelmann J E，Rock B N，Moss D M. Red edge spectral measurements from sugar maple leaves[J]. International Journal of Remote Sensing，1993，14(8)：1563-1575.

[9]Blackburn G A. Quantifying chlorophylls and carotenoids at leaf and canopy scales：An evaluation of some hyperspectral approaches[J]. Remote Sensing of Environment，1998，66(3)：273-285.

[10]le Maire G，F(xiàn)rancois C，Dufrene E. Towards Universal broad leaf chlorophyll indices using PROSPECT simulated database and hyperspectral reflectance measurements[J]. Remote Sensing of Environment，2004，89(1)：1-28.

10.15889/j.issn.1002-1302.2017.02.063

2015-12-10

國(guó)家自然科學(xué)基金(編號(hào)：61465011)；新疆生產(chǎn)建設(shè)兵團(tuán)成果轉(zhuǎn)化項(xiàng)目(編號(hào)：2014BD027)。

田 敏(1968—)，男，河南商丘人，博士，副教授，從事農(nóng)業(yè)信息技術(shù)與精準(zhǔn)農(nóng)業(yè)研究。E-mail：t_tianmin@sina.com。

呂 新，教授，主要從事農(nóng)業(yè)信息化和農(nóng)業(yè)生態(tài)研究。E-mail：lxshz@126.com。

S127；S562.01

A

1002-1302(2017)02-0216-04

田 敏，周 杰，張 澤，等. 基于高光譜植被指數(shù)對(duì)棉花葉綠素含量的估算[J]. 江蘇農(nóng)業(yè)科學(xué)，2017，45(2)：216-219.