一、引 言

氮素是作物生長發(fā)育所必需的營養(yǎng)元素，與光合作用有著密切的關(guān)系，直接影響作物長勢(shì)和產(chǎn)量品質(zhì)的形成，而作物氮素營養(yǎng)主要來源于土壤氮素含量。土壤供氮狀況對(duì)作物增產(chǎn)起著至關(guān)重要的作用，在一定范圍內(nèi)增加氮肥施用量有利于小麥產(chǎn)量的提高，但施氮量過高會(huì)導(dǎo)致氮肥利用率的降低，施肥經(jīng)濟(jì)效益下降，因此，提高土壤氮素利用率已成為當(dāng)前精準(zhǔn)農(nóng)業(yè)的主要目標(biāo)之一。實(shí)時(shí)、準(zhǔn)確地獲取小麥葉片氮含量和土壤氮含量，兼顧二者的監(jiān)測(cè)，可以有效地做到合理施用氮肥，提高氮肥利用率及作物產(chǎn)量、品質(zhì)的形成。

高光譜遙感技術(shù)以實(shí)時(shí)、快速和非破壞性等優(yōu)勢(shì)成為當(dāng)前精準(zhǔn)農(nóng)業(yè)技術(shù)之一，在作物葉片氮含量和土壤供氮水平的監(jiān)測(cè)方面表現(xiàn)出良好的應(yīng)用前景。因此，本研究探討了利用高光譜遙感監(jiān)測(cè)小麥葉片氮含量及土壤對(duì)小麥供氮水平的研究。

本研究以2年3點(diǎn)不同施氮水平下2個(gè)小麥品種的田間試驗(yàn)數(shù)據(jù)為對(duì)象，運(yùn)用植被指數(shù)法和 PLSR，分析小麥冠層反射光譜與葉片氮含量、土壤氮含量的定量關(guān)系，比較并選擇小麥葉片氮含量、土壤氮含量的最佳預(yù)測(cè)模型。預(yù)期結(jié)果將為遙感技術(shù)對(duì)小麥葉片氮素營養(yǎng)及土壤氮素水平的診斷，以及合理施氮管理提供理論依據(jù)。

二、材料與方法

騰格里沙漠位于阿拉善地區(qū)東南部，是中國第四大沙漠。騰格里沙漠屬于典型的大陸性干旱氣候。沙漠內(nèi)部沙丘、湖盆、山地、殘丘及平地等交錯(cuò)分布，其中沙丘占71%，湖盆草灘占7%，山地、殘丘及平地占22%。2011年8月下旬，蘭州大學(xué)資源環(huán)境學(xué)院騰格里沙漠科學(xué)考察隊(duì)發(fā)現(xiàn)，騰格里沙漠的主要植物有芨芨草、沙冬青、沙蒿、白刺、老鴰頭、油蒿、駱駝蓬、鵝絨藤、鹽爪爪、沙蔥、梭梭、地錦、沙生針茅、蟲實(shí)、紅砂、霸王、蘆葦?shù)取?/p>

2.1 試驗(yàn)設(shè)計(jì)

本研究共涉及 3 個(gè)不同的田間試驗(yàn)，分別于2011—2013年小麥生長季節(jié)進(jìn)行。

試驗(yàn) 1：2011-2012 年度在山西省曲沃縣試驗(yàn)田進(jìn)行。供試土壤類型為褐土性土，土壤呈中性、微堿性反應(yīng)，礦物質(zhì)、有機(jī)質(zhì)積累較多，腐殖質(zhì)層較厚，肥力較高。土壤耕層有機(jī)質(zhì)含量10.75 g/kg，全氮含量0.49 g/kg，堿解氮27.77 mg/kg，速效磷 15.46 mg/kg。供試品種為‘運(yùn)麥 20410’。該試驗(yàn)數(shù)據(jù)用于建立模型。試驗(yàn)采用單因素隨機(jī)區(qū)組排列，3次重復(fù) 。設(shè)5個(gè)施氮水平，分別為0、75、150、225、300 kg/hm2 純氮，以尿素(N 46%)作氮肥，分40%基肥、60%返青后期施入。磷鉀肥以基肥一次性施入，過磷酸鈣(P2O516.5%)作磷肥，施磷量(P2O5) 100 kg/hm2 ，氯化鉀(K2O61.5%)作鉀肥，施鉀量(K2O) 100 kg/hm2。小麥返青后一次性追肥，田間光譜采樣時(shí)期分別為選擇返青期、拔節(jié)期、孕穗期和灌漿期。

表1 施肥方案設(shè)計(jì)

試驗(yàn)2：2011-2012年度，在山西省臨汾市堯都區(qū)驗(yàn)田各處理曲沃縣處理的方法相同。選擇在冬小麥返青期、拔節(jié)期、孕穗期和灌漿期測(cè)定小麥冠層反射光譜，同步采集葉片、土壤樣本并測(cè)定氮素指標(biāo)，該試驗(yàn)數(shù)據(jù)用于模型驗(yàn)證。

試驗(yàn) 3：2012-2013 年度，設(shè)在太原市東陽試驗(yàn)田，供試土壤為黃土母質(zhì)發(fā)育而成的石灰性褐土，土壤肥力水平中等，土壤耕層有機(jī)質(zhì)含量20.07 g/kg，堿解氮49.3 mg/kg，有效磷18.79 mg/kg，速效鉀247.3 mg/kg。試驗(yàn)采用單因素隨機(jī)區(qū)組排列，3次重復(fù)。供試小麥品種為國審‘晉太170’。試驗(yàn)設(shè)計(jì)、施肥量、田間光譜采集時(shí)間及測(cè)定方法以及田間管理同于試驗(yàn)1。該試驗(yàn)數(shù)據(jù)用于建立模型。

2.2 測(cè)定方法

1.2.1 冠層光譜的測(cè)定

冬小麥的冠層光譜數(shù)據(jù)測(cè)量可采用萊森光學(xué)生產(chǎn)的iSpecField-WNIR系列地物光譜儀，測(cè)量波段范圍350~2500 nm。冠層光譜測(cè)定選擇在天氣晴朗、無風(fēng)或風(fēng)速很小時(shí)進(jìn)行，測(cè)定時(shí)間為 10:00-14:00。測(cè)量時(shí)傳感器探頭垂直向下，距冠層頂垂直高度約1.5 m，地面視場(chǎng)范圍直徑0.44 m。測(cè)量時(shí)間為10:00-14:00。每小區(qū)重復(fù)測(cè)量3次，取平均值作為該觀測(cè)點(diǎn)的光譜反射值。測(cè)量過程中及時(shí)進(jìn)行標(biāo)準(zhǔn)白板校正（所測(cè)得的目標(biāo)物體光譜為無量綱的相對(duì)反射率）。

1.2.2 氮素的測(cè)定

在小麥返青期、拔節(jié)期、孕穗期和灌漿期測(cè)定小麥冠層反射光譜后，同步測(cè)定小麥葉片氮含量及土壤氮含量。具體方法為：在各小區(qū)隨機(jī)選取5 株代表性植株，將植株莖、葉分離，在 105℃下殺青30 min之后，75℃下烘干至恒重，稱干重，使用凱氏定氮法測(cè)定小麥葉片氮含量(%)；此外，按五點(diǎn)法對(duì)各小區(qū)耕層0~40 cm取土，每個(gè)小區(qū)取混合樣，使用凱氏定氮法測(cè)定并計(jì)算土壤氮含量(g/kg)。

2.3 數(shù)據(jù)分析方法

所采集冠層光譜數(shù)據(jù)處理可采用SpecAnalysis地物光譜數(shù)據(jù)后處理分析軟件進(jìn)行處理。在MATLAB平臺(tái)下編程實(shí)現(xiàn)所有光譜參數(shù)的算法，選擇敏感光譜參數(shù)（表2）并建立線性和非線性（冪函數(shù)、指數(shù)函數(shù)和對(duì)數(shù)函數(shù)）回歸模型，從中選擇最佳預(yù)測(cè)模型。此外，采用偏最小二乘回歸(PLSR)重點(diǎn)分析350~2500 nm 光譜反射率與單因變量（葉片氮含量或土壤氮含量）之間的線性模型，具體分析過程基于 SAS 平臺(tái)下使用PROC PLS過程步，首先對(duì)光譜數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化變換，數(shù)據(jù)分析采用“舍一交叉驗(yàn)法”，并使用迭代的NIPALS算法。

表2 本研究所選高光譜參數(shù)及計(jì)算方法

注：Rλ為波長λnm處的冠層光譜反射率；RCi、RSi、REi分別為反射特征中心點(diǎn)、起點(diǎn)和結(jié)束點(diǎn)的光譜反射率；λCi、λSi、λEi分別為反射特征中心點(diǎn)、起點(diǎn)和結(jié)束點(diǎn)的波長。本文的Areai中的i取值為672 nm，NDi中i的取值為705 nm和672nm。

以試驗(yàn)1和試驗(yàn)3的數(shù)據(jù)為建模數(shù)據(jù)建立基于小麥冠層反射光譜的葉片氮含量及土壤氮含量的預(yù)測(cè)模型，進(jìn)而利用試驗(yàn)2為獨(dú)立資料數(shù)據(jù)對(duì)預(yù)測(cè)模型進(jìn)行驗(yàn)證，模型驗(yàn)證采用預(yù)測(cè)值和實(shí)測(cè)值的擬合決定系數(shù)(R2 )和相對(duì)均方根誤差(RRMSE)等指標(biāo)來評(píng)定模型預(yù)測(cè)精度的好壞，其中，相對(duì)均方根差(RRMSE) 計(jì)算公式如式(1)。

RRMSE值越小則模型的預(yù)測(cè)精度水平越高。

三、結(jié)果與分析

3.1冠層光譜與葉片氮含量及土壤氮含量的相關(guān)分析

利用小麥全生育期光譜反射率分別與葉片氮含量及土壤氮含量進(jìn)行相關(guān)性分析，圖1為小麥葉片氮含量與光譜反射率的相關(guān)性曲線，在可見光波段表現(xiàn)為顯著負(fù)相關(guān)(P<0.05)，其中367~719 nm光譜反射率與葉片氮含量達(dá)到極顯著負(fù)相關(guān)(P<0.01)，而大部分近紅外波段表現(xiàn)為顯著正相關(guān)(P<0.05)，其中738~1143 nm波段光譜反射率與葉片氮含量在全生育期都呈現(xiàn)極顯著正相關(guān)(P<0.01)，且在693 nm和792 nm波段分別達(dá)到最大負(fù)相關(guān)和負(fù)相關(guān)，相關(guān)系數(shù)分別為-0.863 和0.709（表3）。

圖1 小麥葉片氮含量與冠層光譜反射率及一階導(dǎo)數(shù)光譜的相關(guān)系數(shù)

表3 小麥冠層光譜反射率及一階導(dǎo)數(shù)光譜與葉片氮含量、土壤氮含量的相關(guān)系數(shù)

圖2為土壤氮含量與光譜反射率的統(tǒng)計(jì)相關(guān)性曲線，可見光波段表現(xiàn)為顯著負(fù)相關(guān)，432~717 nm光譜反射率與土壤氮含量達(dá)極顯著負(fù)相關(guān)，而大多數(shù)近紅外波段表現(xiàn)為顯著正相關(guān)，739~1141 nm 波段光譜反射率與土壤氮含量均呈顯著正相關(guān)，其中在796 nm和694 nm波段處分別達(dá)到最大負(fù)相關(guān)和正相關(guān)，相關(guān)系數(shù)分別為-0.647和0.721。

針對(duì)全生育期一階導(dǎo)數(shù)光譜與葉片氮含量、土壤氮含量進(jìn)行相關(guān)性分析，分別提取對(duì)一階導(dǎo)數(shù)光譜與葉片氮含量、土壤氮含量的最大正相關(guān)和負(fù)相關(guān)的光譜信息（表3）。

F0.05為顯著相關(guān)；F0.01為極顯著相關(guān)

圖2 土壤氮含量與小麥冠層光譜反射率及一階導(dǎo)數(shù)光譜的相關(guān)系數(shù)

結(jié)果顯示，一階導(dǎo)數(shù)光譜與葉片氮含量、土壤氮含量的最大相關(guān)性均優(yōu)于原始冠層光譜?？赡茉蚴菍?duì)原始光譜的求導(dǎo)能有效地降低背景噪聲，從而增強(qiáng)冠層光譜反射率與目標(biāo)性狀的相關(guān)性。由圖1和圖2可以看出，小麥冠層反射光譜及一階導(dǎo)數(shù)光譜與葉片氮含量、土壤氮含量的相關(guān)性均表現(xiàn)較好，說明光譜信息中含有指示葉片氮含量及土壤氮含量的重要信息，可進(jìn)一步研究反射光譜與葉片氮含量、土壤氮含量之間的回歸關(guān)系。

3.2基于植被指數(shù)的葉片氮含量及土壤氮含量的估算模型

利用試驗(yàn)1和試驗(yàn)3的采集數(shù)據(jù)對(duì)小麥冠層光譜特征參數(shù)與葉片氮含量、土壤氮含量的回歸分析，建立了利用冠層光譜參量反演小麥葉片氮含量、土壤氮含量的預(yù)測(cè)模型，然后利用試驗(yàn)2資料數(shù)據(jù)進(jìn)行模型驗(yàn)證，通過模型間的比較分析，從中優(yōu)選出模型預(yù)測(cè)表現(xiàn)較好的模型（表4）。所選葉片氮含量預(yù)測(cè)模型多數(shù)為線 性 模 型 ，其 中 光 譜 參 數(shù) NDVI(810,560)、ND705、GNDVI、Area672、NDCI、PSNDb、mND705、ND672、SDr/SDb對(duì)小麥葉片氮含量的預(yù)測(cè)效果較好，決定系數(shù)(R2 )和標(biāo)準(zhǔn)誤(RMSE)變化范圍分別為 0.818~0.927，0.397~0.413，模型驗(yàn)證的預(yù)測(cè)精度及相對(duì)均方根誤差的變化范圍分別為0.785~0.808，0.124~0.135。

表4 小麥葉片氮含量和土壤氮含量與最佳光譜參數(shù)的回歸分析及其驗(yàn)證

根據(jù)模型決定系數(shù)最大和標(biāo)準(zhǔn)誤差最小的原則，基于ND705和GNDVI建立的葉片氮含量估算模型的預(yù)測(cè)表現(xiàn)最好（圖3），決定系數(shù)分別為0.827和0.826，均方根誤為0.397和0.398，相應(yīng)的模型驗(yàn)證也有較好的表現(xiàn)，預(yù)測(cè)精度分別為0.808和0.805，相對(duì)均方根誤分別為0.124和0.126。此外，所選光譜參數(shù)VOG2、VOG1、RI-2dB、mSR705、mND705、ND705 對(duì)土壤氮含量所建模型以非線性模型為主，決定系數(shù)和標(biāo)準(zhǔn)誤變化范圍分別為0.634~0.646，0.093~0.094，模型驗(yàn)證的預(yù)測(cè)精度及相對(duì)均方根誤差的變化范圍分別為0.607~0.625，0.127~0.131。根據(jù)模型決定系數(shù)最大和標(biāo)準(zhǔn)誤差最小的原則，基于VOG-2建立二項(xiàng)式回歸的土壤氮含量模型預(yù)測(cè)表現(xiàn)最好（圖4），從圖4可以看出，樣本觀測(cè)值基本散落在回歸曲線附近，決定系數(shù)和均方根誤分別為0.646和0.093，其模型驗(yàn)證表現(xiàn)較為穩(wěn)定，預(yù)測(cè)精度和相對(duì)均方根誤分別為0.625和0.127。

圖3 葉片氮含量與光譜參數(shù)ND705和GNDVI的關(guān)系散點(diǎn)圖

總之，上述所建小麥葉片氮含量和土壤氮含量遙感估算模型的精準(zhǔn)性和穩(wěn)定性均較好，而且模型驗(yàn)證結(jié)果表現(xiàn)也較好。說明利用冠層光譜估算小麥葉片氮含量及土壤氮含量具有可行性和可靠性，擬推薦上述所選模型為小麥葉片氮含量和土壤氮含量的最佳預(yù)測(cè)模型，可用于估算小麥氮素營養(yǎng)水平和土壤供氮狀況。

圖4 光譜參數(shù)VOG2與土壤氮含量的關(guān)系散點(diǎn)圖

四、討論

作物氮素營養(yǎng)與生長發(fā)育、光合作用及產(chǎn)量之間有密切的關(guān)系。葉片氮素是葉綠素分子的基本組分，而葉綠素是作物光合作用的主要色素，其含量的高低直接影響作物光合同化的物質(zhì)積累能力。利用遙感技術(shù)對(duì)作物葉片氮含量與營養(yǎng)豐缺的無損監(jiān)測(cè)，已成為合理施用氮肥，提高氮肥利用效率和作物產(chǎn)量、改善作物品質(zhì)的重要途徑。因此，根據(jù)不同氮肥水平下小麥全生育期氮素營養(yǎng)的高光譜響應(yīng)差異建立的預(yù)測(cè)模型，對(duì)于小麥氮素營養(yǎng)豐缺狀況的診斷及合理施肥管理具有重要意義。但是，小麥氮素營養(yǎng)主要來源于土壤供氮狀況，土壤氮含量及氮肥利用率也是指導(dǎo)氮肥合理施用的重要因素。盡管增施氮肥是提高作物產(chǎn)量的重要措施之一，但由于施肥方法、作物自身原因以及氮肥利用率高低等問題，導(dǎo)致不合理的增施氮肥，不僅浪費(fèi)資源，同時(shí)也對(duì)生態(tài)環(huán)境造成嚴(yán)重的破壞性。

因此，進(jìn)一步研究小麥冠層反射光譜對(duì)土壤供氮狀況的響應(yīng)差異，確立其估算模型，可以實(shí)時(shí)快速地診斷土壤對(duì)小麥的供氮狀況，為合理施氮及管理提供重要的理論依據(jù)。小麥葉片氮含量與冠層光譜及一階導(dǎo)數(shù)的相關(guān)性明顯地優(yōu)于土壤氮含量，說明冠層光譜所含信息主要來源于小麥冠層光譜結(jié)構(gòu)及其化學(xué)成分，而土壤的信息較少。此外，小麥葉片氮含量、土壤氮含量與一階導(dǎo)數(shù)光譜的相關(guān)性顯著地高于冠層原始光譜，說明一階導(dǎo)數(shù)光譜可有效地降低背景噪音，提高光譜指示目標(biāo)性狀的敏感性。冠層光譜與葉片氮含量在可見光部分表現(xiàn)出強(qiáng)烈的線性負(fù)相關(guān)，在近紅外部分表現(xiàn)出顯著正相關(guān)，而其與土壤氮含量的相關(guān)曲線的趨勢(shì)與前者呈現(xiàn)相反趨勢(shì)。究其原因，本研究所測(cè)光譜是以小麥冠層葉片的光譜反射率為主，由于小麥葉片對(duì)土壤的遮擋，使得所測(cè)光譜中僅含有少量反映土壤成份的信息。

此外，本研究所光譜參數(shù) ND705 和GNDVI對(duì)葉片氮含量所建模型的決定系數(shù)分別達(dá)到0.827 和 0.826，光譜參數(shù) VOG2 與土壤氮含量所建二項(xiàng)式曲線回歸的決定系數(shù)達(dá)到0.646，說明敏感光譜參數(shù)均優(yōu)于冠層光譜與目標(biāo)性狀的敏感性。此外，所選光譜參數(shù)的構(gòu)成波段均落入相關(guān)性顯著(P<0.05)的波段范圍，進(jìn)一步證明冠層反射光譜與目標(biāo)性狀相關(guān)性顯著的波段范圍內(nèi)含有指示目標(biāo)性狀的重要信息。所選敏感植被指數(shù)僅包含了指示目標(biāo)性狀的2個(gè)特征波段，而其它波段中所蘊(yùn)含的指示目標(biāo)性狀的信息未能得以充分利用，因此深入挖掘冠層光譜中指示目標(biāo)性狀的重要信息應(yīng)盡量整合大多數(shù)光譜波段反射率的信息。光譜范圍350~1350 nm波段反射率與小麥葉片氮含量、土壤氮含量均有較強(qiáng)的相關(guān)性，因此，有必要進(jìn)一步深入挖掘光譜信息中指示小麥葉片氮含量和土壤氮含量的信息，明確該光譜范圍內(nèi)波段反射率與二者的定量關(guān)系。

五、結(jié)論

小麥冠層反射光譜與葉片氮含量的相關(guān)性分析表明，可見光波段367~719 nm光譜反射率與葉片氮含量達(dá)到極顯著負(fù)相關(guān)(P<0.01)，而近紅外波段738~1143 nm光譜反射率與葉片氮含量均呈現(xiàn)極顯著正相關(guān)(P<0.01)。冠層光譜與土壤氮含量的相關(guān)性分析表明可見光波段432~717 nm光譜反射率與土壤氮含量達(dá)極顯著負(fù)相關(guān)(P<0.01)，而近紅外波段739~1141 nm光譜反射率與土壤氮含量均呈顯著正相關(guān)(P<0.01)?；诠庾V參數(shù)ND705和GNDVI建立的葉片氮含量估算模型的決定系數(shù)分別為 0.827 和 0.826，標(biāo)準(zhǔn)誤分別為0.397和0.398?；诠庾V參數(shù)VOG2建立的土壤氮含量估算模型的決定系數(shù)和標(biāo)準(zhǔn)誤為分別為 0.646 和0.093。所建小麥葉片氮含量的PLSR估算模型的決定系數(shù)和均方根誤差分別為0.842和0.379，所建土壤氮含量 PLSR 模型的決定系數(shù)和均方根誤差分別為0.654和0.092。

推薦：

便攜式高光譜成像系統(tǒng)iSpecHyper-VS1000

專門用于公安刑偵、物證鑒定、醫(yī)學(xué)醫(yī)療、精準(zhǔn)農(nóng)業(yè)、礦物地質(zhì)勘探等領(lǐng)域的最新產(chǎn)品，主要優(yōu)勢(shì)具有體積小、幀率高、高光譜分辨率高、高像質(zhì)等性價(jià)比特點(diǎn)采用了透射光柵內(nèi)推掃原理高光譜成像，系統(tǒng)集成高性能數(shù)據(jù)采集與分析處理系統(tǒng)，高速USB3.0接口傳輸，全靶面高成像質(zhì)量光學(xué)設(shè)計(jì)，物鏡接口為標(biāo)準(zhǔn)C-Mount，可根據(jù)用戶需求更換物鏡。