舒 田，陳智虎，劉春艷，許元紅，趙澤英*

(1.貴州省農(nóng)業(yè)科技信息研究所，貴州 貴陽(yáng) 550006；2.貴州師范大學(xué) 喀斯特研究院/國(guó)家喀斯特石漠化防治工程技術(shù)研究中心，貴州 貴陽(yáng) 550001)

0 引言

【研究意義】火龍果因其營(yíng)養(yǎng)豐富、功能獨(dú)特，含有一般植物少有的植物性白蛋白和花青素，以及豐富的維生素和水溶性膳食纖維，具有減肥美容、降低膽固醇、預(yù)防便秘和大腸癌等功效而深受大眾青睞，已成為近年來(lái)我國(guó)種植面積增長(zhǎng)最快的特色水果之一。截至2018年，貴州省火龍果種植面積達(dá)到5 867 hm2，產(chǎn)量達(dá)4.6萬(wàn)t，位居全國(guó)第3位。作為我國(guó)特色優(yōu)勢(shì)農(nóng)產(chǎn)品，在農(nóng)業(yè)增效、農(nóng)民增收、產(chǎn)業(yè)扶貧等方面發(fā)揮了重要作用[1]，已成為貧困地區(qū)富民興村的搖錢(qián)樹(shù)和致富果。高光譜具有光譜分辨率高、光譜連續(xù)性強(qiáng)、光譜信息量大等優(yōu)勢(shì)，可實(shí)時(shí)記錄植物器官生長(zhǎng)發(fā)育中的細(xì)微變化[2]?；诟吖庾V分析植被植株冠層與葉片的光譜信息，可以進(jìn)行作物種類(lèi)識(shí)別和產(chǎn)量預(yù)測(cè)估算[3]，還可以評(píng)估作物的生長(zhǎng)狀況和受病害程度[4]。利用高光譜遙感技術(shù)對(duì)作物進(jìn)行識(shí)別提取、產(chǎn)量預(yù)測(cè)、病蟲(chóng)害監(jiān)測(cè)等已成為當(dāng)前學(xué)術(shù)界的研究熱點(diǎn)?！厩叭搜芯窟M(jìn)展】閆苗等[5]研究卷積神經(jīng)網(wǎng)絡(luò)模型在高光譜遙感地物多分類(lèi)識(shí)別中可行性以及識(shí)別效果，證實(shí)卷積神經(jīng)網(wǎng)絡(luò)在高光譜遙感地物識(shí)別中具有較高的識(shí)別率。雷彤等[6]研究蘋(píng)果樹(shù)果期冠層反射光譜特性和敏感波段，篩選最佳光譜參數(shù)并構(gòu)建蘋(píng)果產(chǎn)量的最佳估測(cè)模型。劉穎等[2]利用Normalize預(yù)處理的高光譜數(shù)據(jù)的偏最小二乘法(PLS)模型對(duì)高紡錘形蘋(píng)果樹(shù)冠層單位面積花量的預(yù)測(cè)效果最優(yōu)。齊浩等[7]研究核桃、棗、香梨、蘋(píng)果、杏5種果樹(shù)，采用不同步長(zhǎng)間隔的平滑濾值及5種數(shù)據(jù)變換方式處理，結(jié)果表明一階微分變換的樹(shù)種識(shí)別精度高達(dá)99.3%。劉建雄等[8]利用光譜特征參數(shù)與葉面積指數(shù)的相關(guān)性，建立茶尺蠖危害程度擬合模型，對(duì)于產(chǎn)業(yè)化和規(guī)?；蟮牟鑸@茶尺蠖危害進(jìn)行監(jiān)測(cè)具有現(xiàn)實(shí)指導(dǎo)意義。張衡等[9]基于高光譜數(shù)據(jù)對(duì)馬尾松早期階段是否感染松萎蔫病進(jìn)行良好監(jiān)測(cè)。周偉等[10]對(duì)三江源區(qū)毒雜草草種和優(yōu)良牧草種進(jìn)行高光譜識(shí)別分析，并通過(guò)馬氏距離法提取其特征波段，為該地區(qū)植被群落生長(zhǎng)監(jiān)測(cè)和畜牧業(yè)可持續(xù)發(fā)展提供了科學(xué)依據(jù)?！狙芯壳腥朦c(diǎn)】許元紅等[11-12]分別利用無(wú)人機(jī)高光譜和可見(jiàn)光波段對(duì)火龍果種植株數(shù)進(jìn)行識(shí)別提取，目前對(duì)火龍果植株各部位的高光譜遙感識(shí)別鮮有研究，而高光譜對(duì)火龍果植株各部位的精準(zhǔn)識(shí)別與分類(lèi)對(duì)于監(jiān)測(cè)火龍果植株患病與否、病害分類(lèi)、產(chǎn)量估測(cè)、理化參數(shù)含量估算以及肥水精準(zhǔn)管理具有重要意義，同時(shí)也為果樹(shù)營(yíng)養(yǎng)和果實(shí)品質(zhì)無(wú)損檢測(cè)等提供了可能[13-15]?！緮M解決的關(guān)鍵問(wèn)題】充分利用高光譜優(yōu)勢(shì)實(shí)測(cè)火龍果植株冠層、果、枝、花的光譜數(shù)據(jù)，通過(guò)原始光譜反射率、光譜不同變換形式、不同類(lèi)型植被指數(shù)進(jìn)行識(shí)別能力分析，然后利用主成分分析(PCA)并提取特征敏感波段(1 054～1 089 nm和812～825 nm)區(qū)域的光譜反射值作為神經(jīng)網(wǎng)絡(luò)輸入層，構(gòu)建三層BP神經(jīng)網(wǎng)絡(luò)模型開(kāi)展火龍果植株各部位的識(shí)別分析，旨在為火龍果長(zhǎng)勢(shì)監(jiān)測(cè)、估產(chǎn)、植株病害高光譜遙感診斷提供可靠依據(jù)。

1 數(shù)據(jù)與方法

1.1 數(shù)據(jù)來(lái)源

1.1.1 數(shù)據(jù)采集地概況 高光譜數(shù)據(jù)采集于貴州省花江示范區(qū)火龍果種植基地，中心位置為E 105°39′50″,N25°40′32″，海拔430～500 m?；瘕埞N植行距為3 m，株距2 m。基地火龍果植株生長(zhǎng)形態(tài)呈星狀，棱邊緣呈波浪狀或鋸齒狀，深綠色或綠色。枝條粗而長(zhǎng)，長(zhǎng)為30～150 cm，寬3～8 cm?；ǔ事┒窢?，花長(zhǎng)25～30 cm，直徑15～25 cm。果呈長(zhǎng)球形，深紅色或粉紅色。

1.1.2 采集方法 測(cè)量?jī)x器采用美國(guó)ASD公司的地物光譜儀(Field Spec 3)，波譜范圍為350～1 025 nm、350～2 500 nm，光譜采樣間隔分別為1.377 nm、1.377 nm@350～1 050 nm、2 nm@1 000～2 500 nm；其光譜分辨率分別為3nm@700nm、10 nm@1 400 nm、2 100 nm、30 nm@1 400 nm、2 100 nm。測(cè)定時(shí)間為2021年8月下旬，天空晴朗無(wú)云、無(wú)風(fēng)(或微風(fēng))、空氣濕度小，測(cè)定時(shí)段在12:00—13:00。測(cè)定時(shí)儀器探頭與作物冠層頂部相差0.5～1 m，垂直向下，每隔 0.5 h對(duì)儀器進(jìn)行優(yōu)化和白板校正。一共采集測(cè)量35處，每處測(cè)量點(diǎn)分別采集火龍果植株冠層、果、枝和花的光譜3次。由于反射率受到大氣、水汽以及植物自身因素的影響，光譜曲線在首尾兩端存在較大噪聲，數(shù)據(jù)顯得異常。因此，剔除異樣光譜曲線，只保留400～1 350 nm波段的反射光譜數(shù)據(jù)進(jìn)行分析處理，光譜反射率以ASCLL形式導(dǎo)出，便于后續(xù)處理分析。

1.2 研究方法

1.2.1 原始光譜不同形式變換 原始反射光譜經(jīng)預(yù)處理后求導(dǎo)，可以減弱土壤背景對(duì)目標(biāo)光譜的影響，也可以消除部分大氣效應(yīng)噪音，提高信噪比，從而可靠地反映植被內(nèi)部生化組成及含量等信息[16]。因此，以原始光譜的一階微分(First order Differential，F(xiàn)D，波長(zhǎng)間隔為3 nm、5 nm、7 nm、9nm、11 nm、13 nm、15 nm)、倒數(shù)的對(duì)數(shù)[lg(1/λ)]、倒數(shù)對(duì)數(shù)的一階微分[ FD(lg(1/λ))](波長(zhǎng)間隔為3 nm、5 nm、7 nm、9 nm、11 nm、13 nm、15 nm)等不同變換形式進(jìn)行分析。

1.2.2 植被指數(shù) 植被對(duì)太陽(yáng)輻射的吸收與反射受到自身含水量和養(yǎng)分特征等因素的影響，而植被指數(shù)(Vegetation Index，VI)可結(jié)合不同波長(zhǎng)范圍的地物反射率以增強(qiáng)植被的某一特征[17]。在紅光與近紅外波段植被具有較強(qiáng)的吸收和反射特性，不同波段組合在一定條件下定量表明植被的生長(zhǎng)狀況。根據(jù)一般衛(wèi)星傳感器通道(近紅外、紅光、綠光)間波長(zhǎng)即830 nm、660 nm、560 nm，相應(yīng)獲取其反射率數(shù)據(jù)，選取DVI、RVI、NDVI、RDVI、SAVI、GNDVI、GRVI、VARI、NPCI、SIPI、WBI等11種常用植被指數(shù)對(duì)火龍果植株冠層、果、枝、花的光譜反射率進(jìn)行識(shí)別能力分析，各指數(shù)的計(jì)算方法見(jiàn)表1。

表1 不同植被指數(shù)及計(jì)算公式 Table 1 Different vegetation indexes and calculation formulas

1.2.3 主成分分析 由于高光譜數(shù)據(jù)繁多冗雜，若將所有光譜數(shù)據(jù)輸入模型，則不利于BP網(wǎng)絡(luò)模型的構(gòu)建。不僅計(jì)算處理能力下降，還增加運(yùn)算復(fù)雜性。因此，有效選取特征波段不僅可提高模型運(yùn)算速度，還具有解析準(zhǔn)確性。主成分分析(PCA)能在不丟失主要光譜信息前提下，將多個(gè)原始變量壓縮為幾個(gè)新變量的分析方法[18]。在主成分分析的基礎(chǔ)上獲取全波段(400～1 350 nm)中的敏感特征波段即最佳識(shí)別波段。載荷作為各個(gè)原始空間向量在各個(gè)新的向量空間(主成分)投影的單位向量，可用來(lái)表示這些原始向量分別在每個(gè)新的向量空間(主成分)上的相關(guān)性[19]。因此，將高光譜數(shù)據(jù)利用PCA降維后提取特征波段，將特征波段的反射值作為神經(jīng)網(wǎng)絡(luò)的輸入，然后構(gòu)建三層BP神經(jīng)網(wǎng)絡(luò)模型。為更好地構(gòu)建BP網(wǎng)絡(luò)模型，對(duì)采集的火龍果植株冠層、果、枝、花的光譜樣本數(shù)據(jù)按2∶1比例隨機(jī)抽取，分成建模集和預(yù)測(cè)集，即隨機(jī)抽樣選取70條實(shí)測(cè)的光譜樣本作為建模集，35條實(shí)測(cè)的光譜樣本作為預(yù)測(cè)集進(jìn)行主成分分析。調(diào)節(jié)神經(jīng)元激勵(lì)的傳遞函數(shù)sigmoid為0.9，輸入層節(jié)點(diǎn)數(shù)為47，輸出層節(jié)點(diǎn)數(shù)為1，通過(guò)多次試驗(yàn)計(jì)算確定隱含層節(jié)點(diǎn)數(shù)為35。允許誤差值為0.00001，最小訓(xùn)練速率為0.1，最大迭代次數(shù)為1000次。

1.3 數(shù)據(jù)處理

采用The Unscrambler V9.1和DPS7.05(Data Processing System)進(jìn)行數(shù)據(jù)分析處理。

2 結(jié)果與分析

2.1 火龍果植株不同部位的原始反射光譜

通過(guò)ViewSpec Pro對(duì)火龍果植株冠層、果、枝、花原始反射光譜進(jìn)行均值處理，獲得波長(zhǎng)在400～1 350 nm反射光譜曲線。從圖1看出，整個(gè)光譜區(qū)間內(nèi)，火龍果植株冠層、果、枝、花光譜反射率差異較大?？梢?jiàn)光400～630 nm花的反射率明顯高于植株冠層、果和枝，果的反射率曲線與其他部位明顯不一致且無(wú)明顯綠峰值；在640～680 nm范圍內(nèi)果實(shí)反射率最高，然后為花、冠層和枝，且均有“紅谷”出現(xiàn)。在680～1 150 nm波長(zhǎng)范圍內(nèi)，反射率差異顯著，反射率大小依次是果>花>冠層>枝；大于1 150 nm的波長(zhǎng)范圍，花與果實(shí)的反射率基本重疊并且最高，枝的反射率最低。由此可見(jiàn)，640～1 150 nm可作為識(shí)別火龍果植株不同部位的波長(zhǎng)范圍。

圖1 火龍果植株不同部位平均光譜反射率曲線 Fig.1 Reflectance curves of mean spectrum for different parts of pitaya plant

2.2 火龍果植株不同部位不同光譜數(shù)據(jù)變換的敏感波長(zhǎng)

從波長(zhǎng)間隔3 nm的一階微分光譜曲線圖(圖2)可看出，在400～680 nm范圍火龍果花的一階微分光譜曲線存在2個(gè)正峰，其他部位則僅1個(gè)正峰；680～1 350 nm范圍，果、枝、花及植株冠層波譜曲線走向基本一致，存在1個(gè)正峰和2個(gè)負(fù)谷?；诠庾V差異最大且易區(qū)分原則，特選取6個(gè)不同波長(zhǎng)處的一階微分(波長(zhǎng)間隔為3 nm)的最佳識(shí)別光譜，最佳識(shí)別光譜波段位置為516 nm、525 nm、686 nm、702 nm、724 nm、736 nm。同理，其余光譜變換形式的光譜波長(zhǎng)位置見(jiàn)表2，其中出現(xiàn)頻次最多的波長(zhǎng)即為識(shí)別火龍果各部位的敏感波長(zhǎng)(特征波長(zhǎng))，分別為520 nm、642 nm、664 nm、726 nm、950 nm、1 000 nm、1 130 nm、1 330 nm和1 345 nm。

圖2 火龍果植株不同部位的一階微分光譜曲線(波長(zhǎng)間隔為3 nm)Fig.2 First-order differential spectral curves of different parts of pitaya plant(with 3 nm wavelength interval)

表2 不同數(shù)據(jù)變換形式識(shí)別光譜的最佳波長(zhǎng)位置 Table 2 Optimal wavelength of spectrum recognition with different data transformation forms

2.3 火龍果植株不同部位的植被指數(shù)光譜識(shí)別

從表3看出，火龍果植株冠層、果、枝、花的RVI和GRVI的值均大于1，最大值分別為9.419和8.212，最小值分別為1.978和2.088，其余9種植被指數(shù)值均小于1。RVI和GRVI指數(shù)的方差與標(biāo)準(zhǔn)差均較大，遠(yuǎn)大于1；其他植被指數(shù)方差與標(biāo)準(zhǔn)差均小于1。結(jié)合舒田等[20]的研究可以得出，RVI對(duì)于火龍果植株冠層、枝和花識(shí)別能力最強(qiáng)，GRVI對(duì)于果的識(shí)別能力最強(qiáng)。

表3 火龍果植株不同部位植被指數(shù)、方差及標(biāo)準(zhǔn)差Table 3 Vegetation index,variance and standard deviation of different parts of pitaya plant

2.4 火龍果植株BP神經(jīng)網(wǎng)絡(luò)模型構(gòu)建與驗(yàn)證

2.4.1 主成分提取 表4反映火龍果植株各部位4個(gè)主成分分析的累積可信度，建模集的前2個(gè)主成分的累計(jì)可信度達(dá)94%，故選取前2個(gè)主成分作為特征敏感波段分析。

表4 火龍果各部位主成分分析的累計(jì)可信度Table 4 Accumulative credibility of different stems of pitaya plant by PCA

2.4.2 BP神經(jīng)網(wǎng)絡(luò)模型的輸入變量 從圖3可看出，主成分1與波長(zhǎng)1 054～1 087 nm范圍的相關(guān)性較大，即主成分1與1 054～1 087 nm波段范圍的反射值的相關(guān)性最強(qiáng)，主成分2與波長(zhǎng)812～825 nm波段范圍的反射值的相關(guān)性最強(qiáng)，故從400～1 350 nm范圍的950個(gè)光譜反射值中選出波長(zhǎng)在1054～1087 nm和812～825 nm范圍的共計(jì)47個(gè)反射值作為BP神經(jīng)網(wǎng)絡(luò)的輸入變量建立鑒別模型。

圖3 波長(zhǎng)與主成分相關(guān)性Fig.3 Correlation between wavelength and principal components

2.4.3 模型驗(yàn)證 在建模設(shè)置條件下不斷通過(guò)調(diào)整隱含層的節(jié)點(diǎn)數(shù)來(lái)優(yōu)化網(wǎng)絡(luò)模型，最終確定網(wǎng)絡(luò)模型結(jié)構(gòu)為47(輸入層)-35(隱含層)-1(輸出)共3層的BP神經(jīng)網(wǎng)絡(luò)模型，所建模型對(duì)70個(gè)建模樣本的擬合殘差值為0.000 325。通過(guò)對(duì)35個(gè)樣本的反射率實(shí)測(cè)值與BP神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)值(表5)進(jìn)行統(tǒng)計(jì)分析表明，所建模型識(shí)別準(zhǔn)確率為82.8%。雖然模型對(duì)樣本預(yù)測(cè)值準(zhǔn)確率不算太高，但將最佳識(shí)別波段縮小在近紅外的小區(qū)間內(nèi)。

表5 火龍果植株35個(gè)樣本的反射率實(shí)測(cè)值與BP神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)值 Table 5 Predicted value of BP neural network and actual value of 35 samples of pitaya plant

3 討論

火龍果植株各部位(冠層、果、枝、花)的反射光譜存在較大差異，火龍果果實(shí)、花的光譜曲線在400～630 nm處與綠色植物葉片或冠層不太一致，沒(méi)有“綠峰”出現(xiàn)；可見(jiàn)光區(qū)域植被原始光譜反射率較低，經(jīng)過(guò)導(dǎo)數(shù)變換后限制了低頻背景光譜對(duì)目標(biāo)光譜的影響，經(jīng)過(guò)對(duì)數(shù)變換則可以增強(qiáng)可見(jiàn)光區(qū)域的光譜差異，還能減少因光照條件變換引起的乘性因素影響[21-22]。一階微分、倒數(shù)的對(duì)數(shù)和倒數(shù)對(duì)數(shù)的一階微分等數(shù)據(jù)變換形式對(duì)于增強(qiáng)作物在可見(jiàn)光區(qū)域的光譜信息、提升識(shí)別能力起到很大作用。植被指數(shù)是利用植被光譜數(shù)據(jù)線性和非線性組合構(gòu)建的光譜指數(shù)，近紅外波段是作物葉片健康狀況最靈敏的區(qū)域，指示著植物光合作用是否正常運(yùn)行，近紅外和紅光區(qū)域是構(gòu)建植被指數(shù)的最經(jīng)典波段，構(gòu)建不同植被指數(shù)可以間接反映作物生長(zhǎng)狀況和提高不同作物的識(shí)別能力。基于PCA降維操作剔除了貢獻(xiàn)率相對(duì)較小的主成分從而達(dá)到降維的目的，同時(shí)利用PCA構(gòu)建BP神經(jīng)網(wǎng)絡(luò)模型提取了火龍果各部位識(shí)別的最佳特征波段。未來(lái)將利用高光譜遙感影像-高光譜反射率開(kāi)展“圖-譜”的綜合判別分析，提升火龍果植株不同部位識(shí)別精度，同時(shí)針對(duì)病蟲(chóng)害植株開(kāi)展高光譜調(diào)查，從而為火龍果長(zhǎng)勢(shì)監(jiān)測(cè)、產(chǎn)量估測(cè)、植株病蟲(chóng)害高光譜診斷提供強(qiáng)有力的科學(xué)支撐。

4 結(jié)論

采集田間火龍植株冠層、果、枝、花等各部位反射光譜，通過(guò)原始光譜反射率、光譜不同變換形式、不同類(lèi)型植被指數(shù)進(jìn)行識(shí)別能力分析，利用主成分分析法構(gòu)建BP神經(jīng)網(wǎng)絡(luò)模型并提取出火龍果植株各部位識(shí)別的最佳特征波段，得出以下結(jié)論：一是可見(jiàn)光400～630 nm范圍內(nèi)，火龍果花的光譜反射率最高，冠層、枝、果差異不明顯，而640～1 140 nm光譜反射差異非常明顯；二是采用原始光譜不同變化形式得到區(qū)別火龍果冠層、果、枝、花最佳識(shí)別光譜波段位置為520 nm、642 nm、664 nm、726 nm、950 nm、1 000 nm、1 130 nm、1 330 nm和1 345 nm；三是通過(guò)不同植被指數(shù)計(jì)算，RVI對(duì)于火龍果植株冠層、枝和花識(shí)別能力最強(qiáng)，GRVI對(duì)于果識(shí)別能力最強(qiáng)，同時(shí)也驗(yàn)證了RVI、GRVI植被指數(shù)對(duì)不同作物識(shí)別能力強(qiáng)的結(jié)論；四是利用高光譜主成分降維后提取出1 054～1 087 nm和812～825 nm為火龍果植株各部位識(shí)別的最佳特征波段，通過(guò)構(gòu)建BP神經(jīng)網(wǎng)絡(luò)模型，識(shí)別準(zhǔn)確率達(dá)82.8%。