以下是利用高光譜相機進行小麥生物量估算的一般流程：

1. 數(shù)據(jù)采集：使用搭載高光譜相機的無人機或地面平臺，在小麥生長的關(guān)鍵時期獲取高光譜圖像數(shù)據(jù)。確保在采集數(shù)據(jù)時，光照條件和相機參數(shù)保持一致，以獲得具有可比性的數(shù)據(jù)。

2. 數(shù)據(jù)預(yù)處理：對獲取的高光譜圖像進行預(yù)處理，包括噪聲去除、大氣校正、幾何校正等步驟。這一步驟旨在提高數(shù)據(jù)質(zhì)量，減少非生物因素對光譜的影響。

3. 特征選擇與提取：基于先驗知識或統(tǒng)計分析，選擇與小麥生物量相關(guān)的特征波長或植被指數(shù)。例如，歸一化差異植被指數(shù)(NDVI)是常用的植被指數(shù)，與生物量有較好的相關(guān)性。通過計算這些特征，可以量化小麥的生長狀況。

4. 建模與驗證：利用機器學(xué)習(xí)算法（如偏最小二乘回歸、支持向量回歸等）建立小麥生物量與提取的光譜特征之間的定量模型。使用一部分數(shù)據(jù)作為訓(xùn)練集來訓(xùn)練模型，并用另一部分獨立的數(shù)據(jù)作為驗證集來評估模型的性能。通過調(diào)整模型參數(shù)和算法，優(yōu)化模型的預(yù)測精度。

5. 生物量估算：將建立好的模型應(yīng)用于新的高光譜數(shù)據(jù)，實現(xiàn)對小麥生物量的預(yù)測和估算。這些預(yù)測結(jié)果可以提供關(guān)于小麥生長狀態(tài)、產(chǎn)量潛力的信息，為農(nóng)業(yè)生產(chǎn)管理提供決策支持。

需要注意的是，高光譜相機估算小麥生物量的準(zhǔn)確性受到多種因素的影響，如品種差異、生長環(huán)境、數(shù)據(jù)采集和處理方法等。因此，在實際應(yīng)用中，建議結(jié)合地面觀測數(shù)據(jù)和其他農(nóng)學(xué)參數(shù)，對模型進行修正和校準(zhǔn)，以提高估算結(jié)果的準(zhǔn)確性和可靠性。

相關(guān)實驗研究：

本研究應(yīng)用了400-1000nm的高光譜相機，可采用廣東賽斯拜克科技有限公司產(chǎn)品SP130M進行相關(guān)研究。光譜范圍在400-1000nm，波長分辨率優(yōu)于2.5nm，可達1200個光譜通道。采集速度全譜段可達128FPS，波段選擇后最高3300Hz（支持多區(qū)域波段選擇）。

1 材料與方法

1.1 試驗設(shè)計

為反映通常情況下區(qū)域性小麥的生長及品質(zhì)的特點，試驗區(qū)域為北京順義區(qū)和昌平區(qū)有代表性的農(nóng)田，該區(qū)氣候為典型的暖溫帶半濕潤大陸性季風(fēng)氣候，夏季炎熱多雨，冬季寒冷干燥，春、秋短促。年平均氣溫10～12℃,降雨量600多mm,降水季節(jié)分配很不均勻，全年降水的80%集中在夏季6、7、8三個月。

試驗于2003年進行。選擇冬小麥地20余塊，面積大于50畝，地塊平坦，長勢均勻。地塊由北京農(nóng)業(yè)技術(shù)推廣站提供，這些地塊屬于北京優(yōu)質(zhì)小麥訂單生產(chǎn)地。由農(nóng)民耕作，采用通常的肥水管理，不進行特殊處理，小麥品種是當(dāng)前推廣的高蛋白品種，如京冬8號、京9507等。小麥地塊分布于北京郊區(qū)(圖1),圖中對應(yīng)的紅點為選擇的小麥樣地。

1.2測試方法

小麥冠層光譜測量使用ASD Fieldspec FR2500光譜儀。所有光譜測量是在天氣晴朗，無風(fēng)或風(fēng)速很小時測量的，測量時間為9:00～15:00。視場角25°,距地有150 cm,在視場范圍內(nèi)重復(fù)20次取平均，各測量前后利用參考板進行校正。2003年獲取地面光譜8次，時間分別為：3月30日，4月7日，4月15日，4月23日，5月1日，5月9日，5月17日和5月25日。

生物量測定方法：采用干重法，即破壞性取樣，烘干測定。

1.3研究方法

對在上述實驗區(qū)得到的光譜數(shù)據(jù)選取特征波段，計

算得到高光譜植被指數(shù)。本文選取了9種較常用的高光譜植被指數(shù)，分別是：1.RVI——比值植被指數(shù)RatioVegetation Index;2.NDVI——歸一化植被指數(shù)NormalizedDifference Vegetation Index; 3.NDVIgreen——綠波段植被指數(shù)Normalized Difference Vegetation Index in green;

4.SAVI-—調(diào)整土壤亮度的植被Soil-Adjusted Vegetation

Index;5.OSAVI——土壤調(diào)節(jié)植被指數(shù)OptimizedSoil-Adjusted Vegetation Index; 6.MSAVI——土壤調(diào)整植被指數(shù)Modified Soil Adjusted Vegetation Index;

7.MCARI——葉綠素吸收比值指數(shù)Modified ChlorophyllAbsorbtion Reflectance Index;8.TCARI改進型的葉綠素吸收比率指數(shù)Transformed Chlorophyll Absorption inReflectance Index;9.RDVI——重規(guī)一化植被指數(shù)Re-normalized Difference Vegetation Index.詳見表1。

本次試驗共采集了168組冠層高光譜植被指數(shù)和其相應(yīng)的生物量進行傳統(tǒng)回歸分析，隨機分成兩類，其中一類用于模型訓(xùn)練，包含112組數(shù)據(jù)：另一類用于檢驗?zāi)Ｐ蛿M合結(jié)果，包含56組數(shù)據(jù)。相關(guān)植被指數(shù)據(jù)相關(guān)的波譜數(shù)據(jù)運算所得。

誤差反向傳播(back-propagation簡稱BP)算法是目前人工神經(jīng)網(wǎng)絡(luò)(ANN--artificial neural network)中應(yīng)用最廣泛的學(xué)習(xí)算法，BP神經(jīng)網(wǎng)絡(luò)是一種采用有導(dǎo)師訓(xùn)練學(xué)習(xí)方式的前饋型網(wǎng)絡(luò)，它從導(dǎo)師樣本中學(xué)習(xí)并獲取知識，然后在導(dǎo)師訓(xùn)練模型引導(dǎo)下對數(shù)據(jù)庫中的信息進行

目標(biāo)查找、模式識別和屬性分類等。最基本的BP神經(jīng)網(wǎng)絡(luò)由輸入層、中間隱含層、輸出層3層構(gòu)成，這種網(wǎng)絡(luò)具有高度非線性映射能力和良好穩(wěn)健性模式識別特點。BP算法的基本思想是，學(xué)習(xí)過程由信號的正向傳播與誤差的反向傳播兩個過程組成。正向傳播時，樣本從輸入層傳入，經(jīng)隱含層處理后傳至輸出層。如果輸出層的實際輸出與期望輸出不符，則轉(zhuǎn)入誤差反向傳播，即將輸出誤差以某種形式通過隱含層向輸入層逐層反傳，并將誤差分攤給各層所有單元，進行權(quán)值修正。如此循環(huán)往復(fù)進行權(quán)值修正就是網(wǎng)絡(luò)的學(xué)習(xí)訓(xùn)練過程，此過程一直進行到網(wǎng)絡(luò)輸出誤差減少到可接受的程度[]。網(wǎng)絡(luò)訓(xùn)練過程中，需要指定循環(huán)次數(shù)，如果循環(huán)次數(shù)達到時，模式最大允許誤差尚未滿足要求，運算將停止，然后進行參數(shù)調(diào)整，重新訓(xùn)練。最后用訓(xùn)練好的網(wǎng)絡(luò)進行分類，

獲得各類目標(biāo)的信息分布特征。

決定系數(shù)R可以判斷實測值和預(yù)測值的線性擬合度，但是由于光譜植被指數(shù)和作物氮素含量間存在的非線性關(guān)系，它無法提供有關(guān)光譜預(yù)測能力的信息。因此，本文使用均方根誤差(RMSE)和T值檢驗，來選擇最佳結(jié)果。其表達式如下：

本文數(shù)據(jù)處理和神經(jīng)網(wǎng)絡(luò)訓(xùn)練均采用Matlab 7.0(Mathworks,USA)。

2結(jié)果與分析

2.1回歸分析

圖2清楚顯示光譜與生物量的關(guān)系，在660～760nm波段，相關(guān)系數(shù)的絕對值達到峰值，為0.68,這與其他學(xué)者的研究結(jié)果類似415],在近紅外波段，相關(guān)系數(shù)較低；780nm以上波長的光譜與生物量表現(xiàn)出顯著相關(guān)。

比較各種回歸模型，各植被指數(shù)與生物量的對數(shù)回歸模型具有較高的決定系數(shù)，各決定系數(shù)結(jié)果如表2所示，其中NDVI、NDVIgreen、SAVI、OSAVI和MSAVI與氮素含量具有較高的相關(guān)系數(shù)(R2>0.5)。這主要是由于NDVI、NDVIgreen、SAVI、OSAVI和MSAV主要受氮素含量影響，對高氮素水平更敏感，而氮素含量與生物量有很大的相關(guān)性，這從它們的表達式中也可以看出。

2.2神經(jīng)網(wǎng)絡(luò)分析

本文采用BP神經(jīng)網(wǎng)絡(luò)含有1個輸入層、2個隱含層和1個輸出層。1個輸入層即為每類植被指數(shù)；第一個隱含層包含16個神經(jīng)元，第二個隱含層包含8個神經(jīng)元；1個輸出層即為植株氮素含量。為了更好地訓(xùn)練網(wǎng)絡(luò)，首

先對輸入向量進行歸一化處理，初始訓(xùn)練速率為0.05,初始權(quán)重和閥值為任意值。網(wǎng)絡(luò)訓(xùn)練結(jié)果見表3。

圖3中橫坐標(biāo)為植株生物量實測值，縱坐標(biāo)為各模型估測值。當(dāng)實測值大于0.01時，估測值達到一個平臺值

(約為0.04)。

2.3 BP神經(jīng)網(wǎng)絡(luò)與統(tǒng)計回歸模型比較

結(jié)果顯示，神經(jīng)網(wǎng)絡(luò)方法較傳統(tǒng)的回歸模型有更佳

的擬合效果(見表3)。使用神經(jīng)網(wǎng)絡(luò)方法可以得到更高的R2和T值，以及較低的RMSE。從表3中看出，對于TCARI,采用神經(jīng)網(wǎng)絡(luò)方法，R2比對數(shù)回歸方法提高了

74.45%,RMSE減小了24.53%,而在NDVlgreen中T

值提高了99.8%之多。對于RMSE和T值改進最為顯著的植被指數(shù)有SAVI、MSAVI和MCARI。因此，人工神經(jīng)網(wǎng)絡(luò)模型可以改進作物光譜植被指數(shù)和生物量的相關(guān)性，預(yù)測范圍也有所擴展。

3結(jié) 論

小麥生物量和其高光譜植被指數(shù)的量化為改進遙感作物診斷提供更好的平臺，并可以擴展應(yīng)用于作物其它生化變量的估算。由于植被指數(shù)和作物生物量的復(fù)雜非線性關(guān)系，人工神經(jīng)網(wǎng)絡(luò)方法應(yīng)用于作物的生物量診斷，取得了較好的效果，比較傳統(tǒng)的回歸模型，本文采用的BP神經(jīng)網(wǎng)絡(luò)模型表現(xiàn)出明顯的估算優(yōu)勢，顯著改進了傳統(tǒng)的分析方法。