基于高光譜成像技術(shù)的不同時(shí)期番茄植株的快速無(wú)損判別研究

馬 玲，馬 倩，李亞嬌，張祎洋，王 靜，馬思艷，馬 燕，杜明華，吳龍國(guó)，2*

（1.寧夏大學(xué) 農(nóng)學(xué)院，寧夏 銀川 750001；2.寧夏現(xiàn)代設(shè)施園藝工程技術(shù)研究中心，寧夏 銀川 750001）

番茄作為我國(guó)設(shè)施栽培面積最大的園藝作物，已成為許多地區(qū)的主要經(jīng)濟(jì)支柱作物，寧夏番茄因顏色粉紅、果實(shí)扁圓、含糖度高、貨架期長(zhǎng)、營(yíng)養(yǎng)豐富等特點(diǎn)而深受消費(fèi)者喜愛(ài)。但寧夏降水量小且蒸發(fā)量大，水資源極為短缺，嚴(yán)重制約了設(shè)施蔬菜產(chǎn)業(yè)的發(fā)展。在種植過(guò)程中，適宜濃度的微咸水灌溉不僅能提供充足的水分，保障番茄作物的正常生長(zhǎng)，還可提高番茄果實(shí)的營(yíng)養(yǎng)價(jià)值［1］。

雖然微咸水灌溉可以緩解淡水資源危機(jī)，為作物生長(zhǎng)提供水分，但是持續(xù)微咸水灌溉極易引起作物根層土壤鹽分累積，導(dǎo)致次生鹽漬化，最終影響作物的產(chǎn)量［2－3］。持續(xù)微咸水灌溉可抑制作物生長(zhǎng)，主要?dú)w因于微咸水給作物生長(zhǎng)帶來(lái)的負(fù)效應(yīng)［4］。這種負(fù)效應(yīng)首先會(huì)導(dǎo)致土壤次生鹽漬化，直接或間接對(duì)土壤物理、化學(xué)和生物學(xué)性狀造成負(fù)面影響，如通氣性變差及抗侵蝕能力變?nèi)酰?］、有機(jī)質(zhì)和養(yǎng)分的有效性下降［5］、微生物和酶活性下降，進(jìn)而對(duì)微生物調(diào)控的各類(lèi)過(guò)程（如固氮、溶磷等養(yǎng)分轉(zhuǎn)換及協(xié)助作物抗逆）造成不利影響［6－7］。因此，發(fā)展快速監(jiān)測(cè)作物鹽脅迫的方法是助力寧夏設(shè)施蔬菜產(chǎn)業(yè)發(fā)展的重要任務(wù)。

李莎等［8］發(fā)現(xiàn)灌溉水礦化度是影響作物根層含鹽量的主要因素，楊樹(shù)青等［9］的研究結(jié)果表明對(duì)于不同作物的灌溉水利用效率礦化度的計(jì)算不同，柴付軍等［10］指出，長(zhǎng)期多次向土壤中澆水會(huì)對(duì)其返鹽產(chǎn)生影響。有學(xué)者向土壤澆灌礦化度微咸水，經(jīng)過(guò)長(zhǎng)時(shí)間試驗(yàn)，發(fā)現(xiàn)未對(duì)土壤產(chǎn)生次生鹽漬化效應(yīng)［11－13］。上述研究主要關(guān)注微咸水灌溉方法及其水土環(huán)境效應(yīng)等方面，但對(duì)作物本身的關(guān)注較少（集中在對(duì)作物產(chǎn)量和品質(zhì)的影響等）。所以，有必要從作物受到逆境威脅時(shí)對(duì)生理的參數(shù)響應(yīng)和調(diào)控的視角，探討微咸水灌溉下糧食作物的產(chǎn)量形成機(jī)制及對(duì)微咸水灌溉的調(diào)控技術(shù)，為合理開(kāi)發(fā)利用微咸水給出相應(yīng)的指導(dǎo)。

光譜技術(shù)作為一種新型的無(wú)損探測(cè)技術(shù)，可以利用物體的光譜特性推測(cè)其化學(xué)成分和特征結(jié)構(gòu)，在較大范圍上克服了常規(guī)探測(cè)技術(shù)所存在的缺陷。高光譜成像技術(shù)則由于具備非常廣泛的高光譜反應(yīng)范圍和納米量級(jí)的高光譜分辨率，能夠精確鑒定地物光譜的細(xì)微差別，進(jìn)而揭示地物組成與構(gòu)造等方面的重大差異。Magnus等［14］將相關(guān)優(yōu)化偏移（COW）算法在一個(gè)工業(yè)批處理過(guò)程的數(shù)據(jù)集上進(jìn)行了測(cè)試。結(jié)果表明，COW可用于在線(xiàn)實(shí)時(shí)對(duì)齊，是一種同步時(shí)變過(guò)程數(shù)據(jù)的實(shí)用工具。李嵐?jié)龋?5］考察了7個(gè)鎘脅迫梯度下菊苣葉片的高光譜響應(yīng)特征，王松磊等［16］通過(guò)可見(jiàn)近紅外高光譜成像技術(shù)構(gòu)建了快速檢測(cè)番茄葉片含水率的模型，鄭小慎等［17］認(rèn)為利用高光譜技術(shù)可以對(duì)作物葉綠素含量作出合理評(píng)價(jià)，孫紅等［18］通過(guò)高光譜成像檢測(cè)了土豆葉片的含水量，石吉勇等［19］提出了一種以葉綠素葉面分布特征診斷黃瓜中氮（N）、鎂（Mg）元素虧缺的高光譜檢測(cè)方法。文獻(xiàn)［20］構(gòu)建了一種新的多變量選擇方法——變量組合集群分析（VCPA）法提取特征波長(zhǎng)，通過(guò)指數(shù)遞減函數(shù)消除貢獻(xiàn)較小的變量，縮小了變量空間，其過(guò)程類(lèi)似于達(dá)爾文進(jìn)化論中的“適者生存”原則［21］。目前VCPA法已與其他建模方法相結(jié)合廣泛應(yīng)用于食品質(zhì)量檢測(cè)［22］、作物蛋白質(zhì)含量檢測(cè)［23］、醫(yī)藥［24］以及植物病害診斷［25］等多個(gè)領(lǐng)域。上述研究體現(xiàn)了高光譜成像技術(shù)在植物葉片物質(zhì)含量檢測(cè)方面的優(yōu)勢(shì)。

本文基于高光譜成像技術(shù)對(duì)番茄植株葉片進(jìn)行定性研究，通過(guò)對(duì)番茄植株進(jìn)行不同微咸水處理，采摘果實(shí)成熟期的葉片，利用高光譜成像技術(shù)提取番茄葉片特征光譜參數(shù)，采用偏最小二乘回歸算法（PLSR）分析光譜特征參數(shù)與不同鹽處理植株之間的相關(guān)關(guān)系，篩選并建立了番茄植株的鹽脅迫診斷模型，以期為番茄植株的在線(xiàn)原位監(jiān)測(cè)提供理論依據(jù)。

1 實(shí)驗(yàn)部分

1.1 材 料

供試番茄品種選用“博美2號(hào)”，購(gòu)于寧夏賀蘭縣天源種業(yè)公司。

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

于2020年7月至10月在寧夏大學(xué)國(guó)家實(shí)訓(xùn)基地玻璃溫室中，采用基質(zhì)盆栽的方式進(jìn)行種植，在番茄苗長(zhǎng)到四葉一心時(shí)移植，每盆移栽一株幼苗，待緩苗一周后，對(duì)番茄苗進(jìn)行不同鹽脅迫處理，設(shè)計(jì)A：50 mmol/L NaCl溶液；B：100 mmol/L NaCl溶液；C：150 mmol/L NaCl溶液；D：200 mmol/L NaCl溶液；E：250 mmol/L NaCl溶液共5個(gè)鹽分濃度進(jìn)行處理，對(duì)照組（CK）：0 mmol/L NaCl溶液。每個(gè)處理重復(fù)3次。每隔3 d，每株番茄植株按照500 mL灌溉量進(jìn)行灌溉。采摘果實(shí)成熟期的上、中、下層葉片，每層采摘2片，共采集192片樣本用于模型建立。

1.3 高光譜成像

采用HyperSIS－NIR高光譜成像系統(tǒng)（900 ~ 1 700 nm，光譜分辨率5 nm，256個(gè)波段）對(duì)番茄葉片進(jìn)行光譜圖像收集，該系統(tǒng)由北京卓立漢光儀器有限公司生產(chǎn)的高光譜成像光譜儀、CCD相機(jī)、4個(gè)35 W鹵鎢燈、電控位移平臺(tái)、計(jì)算機(jī)和數(shù)據(jù)采集軟件組成。通過(guò)反復(fù)調(diào)試獲得最優(yōu)采集參數(shù)：物距為400 mm，位移平臺(tái)速度為15 mm/s，曝光時(shí)間為10 ms。

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

1.4.1 光譜預(yù)處理由于高光譜圖像是三維數(shù)據(jù)，信息量冗余，同時(shí)原始光譜存在基線(xiàn)漂移，因此需進(jìn)行預(yù)處理。光譜預(yù)處理的目的在于減弱因儀器噪音、暗電流、多重共線(xiàn)性以及背景等產(chǎn)生的干擾［26－27］，主要方法有多元散射校正（MSC）、標(biāo)準(zhǔn)正態(tài)化（SNV）、正交信號(hào)校正（OSC）、COW等。

1.4.2 特征波長(zhǎng)提取光譜儀采集的光譜數(shù)據(jù)往往數(shù)百甚至上千個(gè)，大量的數(shù)據(jù)會(huì)加大計(jì)算的難度，增加模型的復(fù)雜度，且光譜信息中包含大量的冗余與共線(xiàn)性數(shù)據(jù)，影響模型的效果。一般采用數(shù)據(jù)壓縮提取特征信息或提取特征波長(zhǎng)的方法提取有效信息，減少計(jì)算量，簡(jiǎn)化模型并提高模型的精度［28］。本研究采用間隔變量迭代空間收縮法（iVISSA）、間隔隨機(jī)蛙跳法（IRF）、遺傳偏最小二乘算法（GAPLS）、競(jìng)爭(zhēng)性自適應(yīng)加權(quán)算法（CARS）、變量組合集群分析（VCPA）等方法提取特征波長(zhǎng)。

1.4.3 模型構(gòu)建與評(píng)價(jià)為正確評(píng)估所建模型的準(zhǔn)確度和穩(wěn)定性，并確定樣品間的含水量間隔，預(yù)先將192個(gè)樣本分成校正集和預(yù)測(cè)集，選取3/4的樣本作為建模集，另外1/4樣本則作為預(yù)測(cè)集。以土壤的高反射率和含水量為主要輸入量，通過(guò)偏最小二乘回歸的方法構(gòu)建了土壤含水量高光譜預(yù)測(cè)模型，模型有效性由決定系數(shù)（R2）、校正均方根誤差（RMSEC）、預(yù)測(cè)均方根誤差（RMSEP）和相對(duì)分析誤差（RPD）值進(jìn)行綜合評(píng)估。

2 結(jié)果與討論

2.1 番茄植株不同高度葉片的光譜特征分析

為了分析不同濃度鹽脅迫下上、中、下層葉片的光譜特征，分別提取不同濃度鹽脅迫處理下每株番茄植株不同高度下的葉片平均光譜，結(jié)果如圖1所示。

圖1 不同鹽脅迫下葉片的平均光譜Fig.1 Average spectra of leaves under different salt stress

由圖1可知：不同濃度鹽脅迫處理下的光譜曲線(xiàn)在1 473 nm處有吸收峰，主要是水分子OH—基團(tuán)的二倍頻吸收，可通過(guò)圖中光譜反射率的大小判斷含水率的高低［29］。對(duì)于CK，番茄上、中部葉片的反射率低于下部葉片的反射率，說(shuō)明下部葉片水分含量小于上、中部葉片，符合植物正常生長(zhǎng)的規(guī)律。如圖1A所示，植株在50 mmol/L NaCl的脅迫下，下部葉片的反射率均高于上、中部葉片，說(shuō)明下部葉片的含水率低，這與CK規(guī)律相一致，表明低鹽脅迫并未對(duì)植株上、中、下部葉片在1 473 nm處的光譜產(chǎn)生顯著影響；植株在100、200、250 mmol/L NaCl的脅迫下（圖1B、D、E），上部葉片反射率均高于中、下部葉片，說(shuō)明鹽脅迫處理使得上部葉片在1 473 nm處的光譜反射率增大，充分反映出鹽脅迫的影響，即番茄植株上部葉片最先反應(yīng)缺水；如圖1C，番茄植株中部葉片在1 473 nm處的反射率高于上、下部葉片，說(shuō)明中部葉片含水率低，表明150 mmol/L的NaCl溶液對(duì)番茄植株的影響相對(duì)較高。同一鹽濃度（200、250 mmol/L NaCl溶液）下，植株葉片水分分布呈現(xiàn)下部 > 中部 > 上部，可為今后番茄植株精準(zhǔn)灌溉提供指導(dǎo)。

2.2 不同鹽濃度下番茄植株葉片的光譜特征分析

為了進(jìn)一步分析鹽分對(duì)番茄葉片光譜的影響機(jī)制，對(duì)不同鹽濃度處理后番茄植株不同高度下葉片的光譜數(shù)據(jù)進(jìn)行分析，結(jié)果如圖2所示。

圖2 不同高度下番茄葉片的平均光譜Fig.2 Average spectra of tomato leaves at different heights

由圖2A可知，在1 473 nm處，上部葉片的光譜反射率呈現(xiàn)處理50 = 100 > 200 > CK > 150 > 250的趨勢(shì)，說(shuō)明微量的鹽脅迫導(dǎo)致植株增加了光合作用，避免了鹽脅迫對(duì)機(jī)體造成的損傷；當(dāng)濃度過(guò)大時(shí)，機(jī)體產(chǎn)生自我保護(hù)，氣孔關(guān)閉，從而積累水分，因此光譜反射率降低。中部葉片的光譜反射率呈現(xiàn)處理50 = CK > 100 > 150 > 200 > 250的趨勢(shì)（圖2B），說(shuō)明高濃度的鹽脅迫導(dǎo)致光譜反射率降低。下部葉片在1 473 nm處的光譜反射率呈現(xiàn)處理CK > 50 > 100 > 150 > 200 > 250的趨勢(shì)，且差異顯著（圖2C）。綜上，隨著鹽濃度的增大，中部和下部葉片的光譜反射率呈現(xiàn)降低的趨勢(shì)，而上部葉片則出現(xiàn)反射率先增加再降低的趨勢(shì)。

2.3 模型建立

為了建立不同鹽脅迫下番茄葉片的預(yù)測(cè)模型，采用光譜－理化值共生距離（Sample set partitioning based on joint x－y distance，SPXY）法對(duì)192個(gè)番茄葉片光譜數(shù)據(jù)進(jìn)行樣本劃分。選擇3/4的樣本作為校正集，1/4的樣本作為預(yù)測(cè)集，并保證預(yù)測(cè)集在校正集范圍內(nèi)。按照實(shí)驗(yàn)處理A、B、C、D、E、CK，依次將其賦值為1、2、3、4、5、0，用于模型建立。

2.3.1 光譜預(yù)處理建模為保障精確化學(xué)計(jì)量學(xué)模型的有效性，必須對(duì)原始光譜進(jìn)行預(yù)處理。本文分別選用MSC、SNV、OSC、COW進(jìn)行光譜預(yù)處理，采用PLSR進(jìn)行建模，結(jié)果如表1所示。其中RMSECV為交叉驗(yàn)證均方根誤差，其值越小表明模型效果越好。

表1 經(jīng)不同預(yù)處理的PLSR模型對(duì)比Table 1 PLSR models of different pretreatment

由表1可知：經(jīng)MSC與SNV預(yù)處理所建立模型的決定系數(shù)R2C、R2P以及RPD值均小于原始光譜，RMSEP均大于原始光譜，說(shuō)明MSC與SNV預(yù)處理使得模型的有效性減弱。COW光譜預(yù)處理與原始光譜建模效果差異不明顯，且模型的穩(wěn)定性一致。OSC預(yù)處理所建立的模型的相關(guān)系數(shù)R2P和RPD值均高于其他預(yù)處理方法，RMSEP也比其他預(yù)處理所建模型低，可能是在使用OSC處理后，消除了光譜矩陣中的一些與鹽脅迫信息無(wú)關(guān)的噪聲信號(hào)，提高了模型的預(yù)測(cè)精度。綜上，OSC預(yù)處理后所建模型的有效性最佳，因此，優(yōu)選OSC處理的光譜進(jìn)行后續(xù)的光譜建模。

2.3.2 特征波長(zhǎng)提取為了降低數(shù)據(jù)處理量，通過(guò)Matlab 2020a軟件對(duì)原始數(shù)據(jù)進(jìn)行特征波長(zhǎng)提取，以減少無(wú)用信息，提高模型的預(yù)測(cè)能力和穩(wěn)定性。本文中采用iVISSA、IRF、GAPLS、CARS、VCPA提取特征波段，同時(shí)通過(guò)PLSR法構(gòu)建光譜模型，結(jié)果如表2所示。

表2 不同提取特征波長(zhǎng)的模型對(duì)比分析Table 2 Comparative analysis of models with different extracted characteristic wavelengths

由表2可知，由CARS和iVISSA提取特征波長(zhǎng)構(gòu)建的PLSR模型R2P值更好，且特征波長(zhǎng)相對(duì)較小。而GAPLS提取特征波長(zhǎng)構(gòu)建的PLSR模型的RPD值最低，說(shuō)明其所建模型的有效性最低。VCPA與IRF提取特征波長(zhǎng)所建模型的R2P、RPD差異性小。但VCPA傾向于選擇更少的變量，可在多次指數(shù)衰減函數(shù)（EDF）運(yùn)行后，從其余變量的所有組合中找到最優(yōu)變量子集。因此，VCPA可以在變量較少的數(shù)據(jù)集上獲得更好的預(yù)測(cè)性能［20］。綜合特征波長(zhǎng)數(shù)、R2、RMSEC、RMSEP和RPD值，優(yōu)選VCPA－PLSR模型用于無(wú)損快速診斷番茄植株的鹽脅迫程度，其預(yù)測(cè)模型的R2P與RMSEP分別為0.917、0.456，VCPA法提取的11個(gè)特征波長(zhǎng)為945、975、990、1 002、1 005、1 067、1 204、1 326、1 595、1 642、1 660 nm。該法為快速無(wú)損檢測(cè)番茄植株鹽脅迫程度提供了一種有前景的策略，預(yù)測(cè)效果如圖3所示。

圖3 VCPA－PLSR模型的預(yù)測(cè)效果Fig.3 Prediction effect of VCPA－PLSR model

3 結(jié) 論

本文通過(guò)應(yīng)用近紅外高光譜成像技術(shù)對(duì)不同鹽脅迫處理的番茄葉片進(jìn)行檢測(cè)，基于不同光譜預(yù)處理方法和特征波長(zhǎng)提取方法建立了鹽脅迫葉片的PLSR模型并進(jìn)行比較，優(yōu)選出最佳模型。結(jié)果顯示：預(yù)處理在一定程度上提高了數(shù)據(jù)的預(yù)期效果，對(duì)數(shù)據(jù)噪音波段有一定的改善作用，經(jīng)OSC預(yù)處理后所建模型的有效性最佳，VCPA法提取了11個(gè)特征波長(zhǎng)。所建立的番茄葉片定性模型具有較好的預(yù)測(cè)效果。