應(yīng)用方向：

高光譜成像技術(shù)在農(nóng)業(yè)領(lǐng)域具有廣泛應(yīng)用前景，主要體現(xiàn)在作物生長(zhǎng)監(jiān)測(cè)與精準(zhǔn)農(nóng)業(yè)管理方面。通過獲取作物在可見光—近紅外波段的連續(xù)光譜信息，可以實(shí)現(xiàn)對(duì)葉綠素含量等關(guān)鍵生理參數(shù)的無損反演，從而反映作物的營(yíng)養(yǎng)狀態(tài)、光合作用能力和生長(zhǎng)活力，為田間長(zhǎng)勢(shì)診斷、精準(zhǔn)施肥和產(chǎn)量預(yù)測(cè)提供科學(xué)依據(jù)。同時(shí)，結(jié)合無人機(jī)平臺(tái)的高空間分辨率和靈活觀測(cè)能力，高光譜成像能夠?qū)崿F(xiàn)農(nóng)田區(qū)域尺度的快速監(jiān)測(cè)與動(dòng)態(tài)分析。此外，通過融合地面高精度光譜數(shù)據(jù)與無人機(jī)高光譜數(shù)據(jù)，并結(jié)合機(jī)器學(xué)習(xí)或智能優(yōu)化算法進(jìn)行建模，可以進(jìn)一步提高生理參數(shù)反演精度，推動(dòng)高光譜遙感在智慧農(nóng)業(yè)和作物表型監(jiān)測(cè)中的應(yīng)用。

背景：

水稻是全球*重要的糧食作物之一，其生長(zhǎng)狀況直接關(guān)系到糧食**與農(nóng)業(yè)生產(chǎn)穩(wěn)定。其中，葉綠素含量是反映水稻光合能力、營(yíng)養(yǎng)狀況及生長(zhǎng)活力的重要生理指標(biāo)，因此對(duì)水稻葉綠素含量進(jìn)行準(zhǔn)確監(jiān)測(cè)，對(duì)于作物長(zhǎng)勢(shì)評(píng)估、精準(zhǔn)施肥管理以及產(chǎn)量預(yù)測(cè)具有重要意義。傳統(tǒng)的葉綠素檢測(cè)方法通常依賴化學(xué)分析或分光光度法，通過實(shí)驗(yàn)室提取和測(cè)量葉片樣品獲得葉綠素含量。雖然該類方法具有較高的測(cè)量精度，但需要對(duì)樣品進(jìn)行破壞性處理，且操作過程耗時(shí)、效率較低，難以滿足大尺度農(nóng)田快速監(jiān)測(cè)的需求。

隨著遙感技術(shù)的發(fā)展，高光譜遙感因能夠獲取連續(xù)且豐富的光譜信息，并揭示作物生理參數(shù)與光譜反射特性之間的關(guān)系，已成為作物生長(zhǎng)狀態(tài)監(jiān)測(cè)的重要技術(shù)手段。近年來，無人機(jī)搭載高光譜傳感器因具有機(jī)動(dòng)靈活、空間分辨率高和獲取成本相對(duì)較低等優(yōu)勢(shì)，在農(nóng)業(yè)監(jiān)測(cè)領(lǐng)域得到廣泛應(yīng)用。通過無人機(jī)高光譜數(shù)據(jù)，可以實(shí)現(xiàn)對(duì)水稻冠層光譜信息的大范圍獲取，并進(jìn)一步構(gòu)建葉綠素含量反演模型，從而實(shí)現(xiàn)作物生長(zhǎng)狀態(tài)的快速評(píng)估。

然而，無人機(jī)高光譜數(shù)據(jù)在實(shí)際應(yīng)用中仍存在一定局限性。受傳感器性能、大氣散射、氣溶膠吸收以及土壤和水體背景等因素影響，無人機(jī)獲取的光譜數(shù)據(jù)往往存在一定程度的誤差和噪聲，從而降低了葉綠素反演模型的精度和穩(wěn)定性。相比之下，地面光譜儀能夠在近距離條件下獲取高精度的光譜數(shù)據(jù)，在作物生理參數(shù)反演中具有較高可靠性，但其觀測(cè)范圍有限，難以滿足區(qū)域尺度監(jiān)測(cè)需求。因此，如何有效結(jié)合地面光譜數(shù)據(jù)的高精度優(yōu)勢(shì)與無人機(jī)高光譜數(shù)據(jù)的大范圍觀測(cè)能力，對(duì)無人機(jī)高光譜數(shù)據(jù)進(jìn)行優(yōu)化和校正，從而提高作物葉綠素含量反演精度，成為當(dāng)前農(nóng)業(yè)遙感研究中的一個(gè)重要問題。

作者信息：于豐華，沈陽農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院，博導(dǎo)

期刊來源：Computers and Electronics in Agriculture

研究?jī)?nèi)容

本研究基于單點(diǎn)地面光譜儀數(shù)據(jù)對(duì)無人機(jī)（UAV）高光譜數(shù)據(jù)進(jìn)行優(yōu)化，在可見光和近紅外區(qū)域開展光譜優(yōu)化。通過計(jì)算非成像高光譜反射率與不同區(qū)域積分后的無人機(jī)高光譜反射率的比值，獲得優(yōu)化系數(shù)并應(yīng)用于相應(yīng)波段，得到優(yōu)化后的高光譜反射率。在此基礎(chǔ)上，分別利用ELM、CPO-ELM和FLA-ELM構(gòu)建水稻葉片葉綠素含量反演模型。該方法既保留了無人機(jī)大面積觀測(cè)的優(yōu)勢(shì)，又通過單點(diǎn)數(shù)據(jù)提升了光譜精度，為水稻生長(zhǎng)狀況監(jiān)測(cè)提供了一種新的技術(shù)方案。本研究總體技術(shù)流程如圖1所示。

圖1 技術(shù)流程總圖

實(shí)驗(yàn)設(shè)計(jì)

本研究在遼寧省海城市耿莊鎮(zhèn)沈陽農(nóng)業(yè)大學(xué)海城實(shí)訓(xùn)基地試驗(yàn)田開展。試驗(yàn)水稻品種分別為"沈農(nóng)9816"和"金稻香103"。研究區(qū)概況如圖2所示，其中大田編號(hào)為1–8號(hào)，每塊面積660 m2；小田編號(hào)為1–15號(hào)，每塊面積40 m2。田間試驗(yàn)設(shè)置5個(gè)氮肥梯度處理：零氮、低氮、中氮、標(biāo)準(zhǔn)氮和高氮。其中N0不施氮肥；N3為220.5 kg/ha，基肥、蘗肥、穗肥施用比例為5:3:2；N1、N2、N4分別為N3的50%、75%和125%。所有小區(qū)田間管理措施一致。2023年數(shù)據(jù)集包括大田試驗(yàn)區(qū)1–8號(hào)地塊和小田試驗(yàn)區(qū)1–15號(hào)地塊的樣本；2024年數(shù)據(jù)集包括大田試驗(yàn)區(qū)1–4號(hào)地塊的樣本。

圖2 研究區(qū)概況

本研究采用深圳大疆生產(chǎn)的 M600PRO 六旋翼無人機(jī)，搭載GaiaSky-mini 內(nèi)置推掃式高光譜成像系統(tǒng)（江蘇雙利合譜科技有限公司，中國(guó)）作為核心數(shù)據(jù)采集設(shè)備。該高光譜系統(tǒng)的波長(zhǎng)范圍為 400–1000 nm，采樣間隔為 3 nm，光譜分辨率為 3.5 nm，無人機(jī)飛行高度為 100 m，共獲取 170 個(gè)有效波段。為減少天氣條件、太陽高度角及太陽輻射對(duì)光譜數(shù)據(jù)的影響，無人機(jī)高光譜數(shù)據(jù)采集選擇在晴朗、無云、無風(fēng)天氣條件下進(jìn)行。本研究將數(shù)據(jù)采集時(shí)間段設(shè)定為 12:00–14:00。在起飛前，對(duì)高光譜成像儀進(jìn)行了暗電流校正和白板校正。此外，在每個(gè)場(chǎng)景的高光譜數(shù)據(jù)采集區(qū)域內(nèi)放置了一塊尺寸為 1.5 m × 1.5 m、反射率為 60% 的漫反射校準(zhǔn)板，用于后續(xù)的反射率數(shù)據(jù)校正。數(shù)據(jù)采集完成后，利用機(jī)載高光譜成像系統(tǒng)配套的 SpectraView 軟件（雙利合譜，中國(guó)）對(duì)獲取的無人機(jī)高光譜遙感圖像進(jìn)行一系列預(yù)處理，包括鏡頭畸變校正、反射率定標(biāo)和大氣校正等步驟，*終生成水稻冠層的高光譜影像數(shù)據(jù)。

采集單點(diǎn)高光譜數(shù)據(jù)。波長(zhǎng)范圍為 325–1075 nm，光譜測(cè)量分辨率為 3 nm，光譜采樣間隔為 1 nm。所有高光譜測(cè)量均在預(yù)先設(shè)定的樣點(diǎn)處進(jìn)行，這些樣點(diǎn)位于水稻移栽行的正北方向。在測(cè)量過程中，光譜儀的視場(chǎng)角設(shè)置為 25°，并將光譜儀置于距水稻植株 30 cm 處進(jìn)行觀測(cè)。為避免操作者的陰影投射到觀測(cè)目標(biāo)上，傳感器探頭保持垂直向下的方向進(jìn)行測(cè)量。為了減少高光譜測(cè)量誤差，光譜儀每 10 min 進(jìn)行一次標(biāo)準(zhǔn)白板校準(zhǔn)和背景噪聲數(shù)據(jù)校準(zhǔn)。

用分光光度計(jì)測(cè)定葉片的葉綠素含量。具體過程為：首先采集與光譜測(cè)量對(duì)應(yīng)的水稻葉片樣品，稱取約 0.1 g 新鮮葉片，并加入石英砂進(jìn)行充分研磨，使組織完*破碎。隨后使用 95% 乙醇溶液提取葉片中的葉綠素，并將提取液轉(zhuǎn)移至 50 mL 容量瓶中定容、混勻。將提取液置于 1 cm 光程的比色皿中，經(jīng)過避光處理約24 h 后，利用 UV–Vis 分光光度計(jì)在 649 nm 和 665 nm 波長(zhǎng)處測(cè)定溶液的吸光度。*后利用測(cè)得的吸光度值計(jì)算葉片中的葉綠素含量。

研究方法

本研究利用分析軟件工具提取采樣點(diǎn)。具體流程如下：首先確定各采樣網(wǎng)格中心點(diǎn)在真北方向上的坐標(biāo)，以此點(diǎn)為中心建立5 × 5像元的感興趣區(qū)域（ROI），提取該區(qū)域內(nèi)所有像元的平均光譜值。采用三次插值法，將無人機(jī)高光譜反射率數(shù)據(jù)重采樣為1 nm的波段間隔，以匹配非成像高光譜數(shù)據(jù)的波段數(shù)量和位置。鑒于稻田地物環(huán)境復(fù)雜，無人機(jī)采集的水稻冠層光譜反射率可能受到地物、水體等因素干擾，產(chǎn)生不同程度的噪聲。本研究采用高斯濾波算法對(duì)無人機(jī)高光譜反射率進(jìn)行平滑處理，高斯濾波器標(biāo)準(zhǔn)差設(shè)為11（圖3）。

圖3 無人機(jī)高光譜預(yù)處理結(jié)果。（a）感興趣區(qū)域光譜提??；（b）平滑處理前后對(duì)比圖；（c）平滑處理前后局部放大圖；（d）全樣本預(yù)處理結(jié)果

為確保無人機(jī)高光譜與非成像高光譜傳感器的光譜范圍一致，本研究選取400–1000 nm波長(zhǎng)范圍內(nèi)的非成像高光譜反射率數(shù)據(jù)。由于地面光譜儀各部件在運(yùn)行過程中易產(chǎn)生噪聲，采用高斯濾波算法對(duì)非成像高光譜反射率數(shù)據(jù)進(jìn)行平滑處理，高斯濾波器標(biāo)準(zhǔn)差設(shè)為11。預(yù)處理結(jié)果如圖4（a）所示。

以 700 nm 為界將光譜劃分為可見光區(qū)（400–700 nm）和近紅外區(qū)（701–1000 nm）。隨后分別計(jì)算單點(diǎn)光譜與 UAV 光譜在兩個(gè)波段內(nèi)的積分比值，得到可見光區(qū)優(yōu)化系數(shù)和近紅外區(qū)優(yōu)化系數(shù)，并將該系數(shù)分別作用于對(duì)應(yīng)波段的 UAV 反射率，從而獲得優(yōu)化后的光譜數(shù)據(jù)?？紤]到分區(qū)校正可能在 700 nm 附近產(chǎn)生光譜不連續(xù)現(xiàn)象，本文進(jìn)一步在 680–720 nm 范圍內(nèi)構(gòu)建過渡帶，通過三次樣條插值對(duì)該區(qū)域光譜進(jìn)行平滑重建，以保證光譜曲線在可見光與近紅外波段之間的連續(xù)性和穩(wěn)定性。

此外，本文提出一種基于敏感域引導(dǎo)的隨機(jī)遍歷策略進(jìn)行特征波段選擇。在 400–700 nm 范圍內(nèi)，通過 Monte Carlo 隨機(jī)組合搜索生成候選波長(zhǎng)組合，并采用六波段組合進(jìn)行遍歷計(jì)算，通過約10,000 次隨機(jī)迭代建立候選波段庫，*終從中篩選出對(duì)葉綠素反演*具代表性的特征波段組合，從而減少冗余信息，提高模型建模效率和預(yù)測(cè)精度。

在模型的建立上，首先采用極限學(xué)習(xí)機(jī)（ELM）作為基礎(chǔ)模型。進(jìn)一步引入冠豪豬優(yōu)化算法（CPO）對(duì) ELM 進(jìn)行改進(jìn)，CPO通過全局搜索策略優(yōu)化參數(shù)組合。CPO算法的種群規(guī)模設(shè)置為30，*大迭代次數(shù)設(shè)置為250。此外，洪水優(yōu)化算法（FLA）是一種新提出的元啟發(fā)式優(yōu)化算法，本文將 FLA 算法與ELM相結(jié)合，構(gòu)建 FLA-ELM 模型，利用洪水優(yōu)化算法*有的多尺度搜索機(jī)制（包括梯度流探索與滲透效應(yīng)開發(fā)），以避免 ELM 中隨機(jī)參數(shù)引起的局部*優(yōu)問題。

本研究的樣本總量為 86 個(gè)，其中 2023 年樣本 57 個(gè)，2024 年樣本 29 個(gè)。所有樣本按照 8:2 的比例 隨機(jī)劃分為建模集和驗(yàn)證集，用于建立和驗(yàn)證反演模型。以 Monte Carlo 模擬生成的六個(gè)候選波長(zhǎng)組合作為輸入變量，以水稻葉綠素含量作為輸出變量，分別構(gòu)建 ELM、CPO-ELM 和 FLA-ELM 三種水稻葉綠素含量反演模型。在相同反演模型條件下，對(duì)原始 UAV 高光譜反射率數(shù)據(jù)與優(yōu)化后的 UAV 高光譜反射率數(shù)據(jù)的反演精度進(jìn)行對(duì)比分析。利用相關(guān)評(píng)價(jià)指標(biāo)（如 R2 和 RMSE）對(duì)模型預(yù)測(cè)性能進(jìn)行評(píng)估，以實(shí)現(xiàn)對(duì)水稻葉綠素含量的準(zhǔn)確估計(jì)。

結(jié)果

如圖4（b）所示，從總體形態(tài)和趨勢(shì)來看，單點(diǎn)高光譜數(shù)據(jù)與無人機(jī)高光譜的反射率曲線高度一致，主要波段光譜輪廓相似，均呈現(xiàn)出植被典型光譜特征，如紅邊處明顯躍升、近紅外波段高反射等。但兩者在反射率數(shù)值范圍和幅度上存在顯著差異：非成像高光譜反射率范圍更寬、整體數(shù)值更高，體現(xiàn)了其在標(biāo)定條件和儀器靈敏度方面的優(yōu)勢(shì)；而無人機(jī)高光譜反射率數(shù)值偏低且幅度變化較小，可能受到大氣干擾、光照條件不均或傳感器性能限制等因素影響。此外，在可見光區(qū)（400–700 nm），無人機(jī)高光譜數(shù)據(jù)整體反射率約為0.06–0.1，單點(diǎn)反射率略高，約為0.07–0.16；在550 nm葉綠素吸收谷處兩者差異達(dá)0.06。隨著波長(zhǎng)延伸至近紅外區(qū)，兩者差異急劇擴(kuò)大，非成像高光譜數(shù)據(jù)表現(xiàn)出更明顯的區(qū)分度，差異曲線在900 nm處達(dá)到0.25的峰值。綜上所述，兩種方法總體趨勢(shì)一致，但在反射率精度和細(xì)節(jié)表征方面仍存在差異。

圖4 （a）單點(diǎn)數(shù)據(jù)預(yù)處理結(jié)果；（b）單點(diǎn)與無人機(jī)高光譜成像反射率差異對(duì)比