利用高光譜成像技術對紅辣椒根腐病的檢測與分析

背景：中(zhong)國(guo)是(shi)農(nong)業(ye)生(sheng)產(chan)大(da)國(guo)，農(nong)業(ye)在(zai)國(guo)民(min)經(jing)濟(ji)的(de)組(zu)成(cheng)和(he)發(fa)展(zhan)中(zhong)的(de)地(di)位(wei)至(zhi)關(guan)重(zhong)要(yao)。蔬(shu)菜(cai)種(zhong)植(zhi)是(shi)農(nong)業(ye)生(sheng)產(chan)種(zhong)植(zhi)業(ye)的(de)重(zhong)要(yao)分(fen)支(zhi)，是(shi)農(nong)民(min)重(zhong)要(yao)的(de)收(shou)入(ru)來(lai)源(yuan)之(zhi)一(yi)。自(zi)從(cong)中(zhong)國(guo)加(jia)入(ru)WTO之後
，中國的蔬菜生產產業在國際上具有更大的競爭優勢。辣椒作為產業鏈最為完整的蔬菜品類，在我國有400多年的種植曆史，因其能夠分泌辣椒素而味辣，維生素C占蔬菜首位
，從而在我國調味品種占有重要位置。近年來
，我國辣椒種植規模不斷增長
，產量不斷上升，辣椒種植麵積約220萬hm2
，產值達2500億元
，占全國蔬菜總值的11.36%，辣椒產業對農民收入貢獻率達到1.14%。近年來
，受新冠肺炎疫情、水澇、幹旱等自然災害以及病蟲害影響導致辣椒單產下降，利潤率下降。其中由於栽培模式固定
、重茬麵積大
、管理粗放、抗kang禦yu旱han澇lao災zai害hai的de能neng力li不bu強qiang，使shi得de辣la椒jiao根gen腐fu病bing在zai各ge個ge地di區qu呈cheng現xian逐zhu年nian加jia重zhong趨qu勢shi，造zao成cheng辣la椒jiao減jian產chan或huo絕jue產chan
，嚴yan重zhong製zhi約yue著zhe辣la椒jiao產chan業ye的de發fa展zhan，給gei農nong戶hu造zao成cheng嚴yan重zhong的de經jing濟ji損sun失shi。辣la椒jiao根gen腐fu病bing是shi典dian型xing的de土tu傳chuan病bing害hai，可ke以yi引yin起qi辣la椒jiao根gen係xi和he近jin地di麵mian莖jing稈gan木mu質zhi部bu

、renpibuhuaisi
，daozhilajiaobunengshusongshuifenheyingyangwuzhi
，fabinghoudixiabufengenxihejingganquanbuhuaisi
，dishangbugankusiwang。gaibinghaizaiwoguogelajiaochanqujunyoufasheng，yibaodaodelajiaogenfubingzhibingjunyoulianbaoshu
、疫黴屬和腐黴屬等，其中鐮孢菌是引起辣椒根腐病的主要病原菌，病原菌以厚垣孢子、junhehuojunsitizaiturangzhongyuedong，chengweiyinianzhuyaochuqinranyuan
，bingjunconggenjingbuhuogenbushangkouqinru，tongguoyushuihuoguangaishuichuanbohemanyan。yinci，youlianbaojunyinqilajiaogenfubingfashengzhouqichang，zaiyoumiaoqidaochengzhuqijunnengfasheng，daozhilajiaojianchanzaochengyanzhongjingjisunshi。lajiaogenfubingjuyoujiqiangdeyinbixing，fabingqianqirouyannanyifaxian。duiyulajiaogenfubingdengturangchuanbodebinghaizhili，yibandeyufangfangfaweidingqijianzhazhibingbingyuantizaiturangzhongdefenbuhekongzhiturangzhongbingyuantidemidu
，bingtongguoshengwuhehuaxuekongzhijilaizuduanqichuanbo。qizhongdingqijianzhaxuyaohaofeizhilirenyuandaliangshijianqieduizhilirenyuandezhuanyeyaoqiujiaogao。yinci，kaifayizhongnenggoutongguolajiaozhizhudedishangbufen（葉或莖）在未顯現症狀的情況下檢測出辣椒根腐病的方法（新思路）是辣椒種植農戶的迫切需求，對辣椒種植和育種產業有重要意義。

實驗設計：樣品製備於2022年12月yue在zai山shan東dong省sheng濟ji寧ning市shi農nong業ye科ke學xue院yuan的de實shi驗yan室shi進jin行xing。為wei了le擴kuo大da樣yang本ben容rong量liang和he提ti高gao模mo型xing的de容rong錯cuo率lv
，方fang便bian將jiang模mo型xing移yi植zhi到dao其qi他ta辣la椒jiao品pin種zhong中zhong。試shi驗yan選xuan取qu兩liang種zhong辣la椒jiao品pin種zhong（山椒4號和滿山紅）種子進行幼苗培育。辣椒幼苗的培養步驟具體如下：將攝氏溫度55°的溫水倒入放好種子的燒杯中，使用玻璃棒連續攪拌10 min;倒掉燒杯中的溫水，將種子衝洗兩次後浸種8~10h;將種子均勻放在濕毛巾上（保持濕潤），將濕毛巾放在托盤中
，將托盤放置在光照培養箱下催芽（一般四天出芽）
；將發芽種子轉移到含育苗基質（珍珠岩
、多菌靈和有機營養土等）的育苗杯中，並放回光照培養箱下培養17天，每個育苗杯中成活1~3株。光照培養箱的條件設置為：亮光環境培養時間為10h，溫度=28°


，濕度=60%，光照強度為4500~10000Lux；暗環境培養時間為14h，溫度=23°

，濕度=60%。

通過人工接種的方式使辣椒幼苗感染根腐病病菌，接種步驟符合中華人民共和國農業行業標準 （NT/T2060 1-2011）。試驗所用根腐病病菌原液由濟寧市農業科學院培養和提供，使用灌根接種法接種根腐病病原菌。具體做法是：從兩個品種中分別挑選20株長勢一致、健康的辣椒幼苗進行病菌接種；接種前澆灌幼苗以保持土壤濕潤，使用玻璃棒在辣椒幼苗莖基部附近紮一個小穴
；使用移液槍將5ml孢子懸浮液（濃度為103個遊動孢子/ml）注入穴內。接種後

，將辣椒幼苗放回光照培養箱

，更改亮光環境培養時間為12h，暗環境培養時間為12h
，其他條件不變。

在(zai)實(shi)驗(yan)過(guo)程(cheng)中(zhong)，接(jie)種(zhong)辣(la)椒(jiao)幼(you)苗(miao)發(fa)病(bing)明(ming)顯(xian)。在(zai)接(jie)種(zhong)病(bing)菌(jun)第(di)一(yi)天(tian)肉(rou)眼(yan)無(wu)法(fa)觀(guan)察(cha)出(chu)症(zheng)狀(zhuang)。接(jie)種(zhong)後(hou)第(di)二(er)天(tian)，兩(liang)個(ge)品(pin)種(zhong)的(de)辣(la)椒(jiao)幼(you)苗(miao)都(dou)陸(lu)續(xu)出(chu)現(xian)根(gen)腐(fu)病(bing)的(de)典(dian)型(xing)症(zheng)狀(zhuang)，具(ju)體(ti)表(biao)現(xian)為(wei)葉(ye)片(pian)萎(wei)蔫(nian)且(qie)莖(jing)基(ji)部(bu)收(shou)縮(suo)
，莖(jing)部(bu)輕(qing)微(wei)倒(dao)伏(fu)。在(zai)接(jie)種(zhong)第(di)三(san)天(tian)
，辣(la)椒(jiao)幼(you)苗(miao)的(de)葉(ye)片(pian)萎(wei)縮(suo)嚴(yan)重(zhong)且(qie)莖(jing)基(ji)部(bu)呈(cheng)暗(an)褐(he)色(se)，造(zao)成(cheng)幼(you)苗(miao)折(zhe)倒(dao)枯(ku)死(si)。將(jiang)辣(la)椒(jiao)葉(ye)片(pian)樣(yang)本(ben)分(fen)為(wei)三(san)個(ge)階(jie)段(duan)
，接(jie)種(zhong)病(bing)菌(jun)前(qian)作(zuo)為(wei)健(jian)康(kang)期(qi)（0d），接種病菌第一天作為潛伏期（1d）
，接種病菌第二天作為發病期（2d）。

這個試驗使用GaiaField-V10E便攜式可見光/近紅外成像係統（Jiangsu Dualix Spec­tral Imaging Technology Co., Ltd, Wuxi, Jiangsu, China）采集高光譜圖像信息。該係統包括光譜儀（GaiaField-V10E）
、成像鏡頭（HSIA-OL23）、專用光源（HSIA-LS-T-200 W）、標準白板（HSIA-CT-150 × 150）和專用電腦（(Surface Go 2 4425Y, Microsoft., USA）等組件。該設備可以在可見光和近紅外波段範圍內進行成像，波長覆蓋範圍為400~1000nm，光譜分辨率為2.8nm
，每個樣本可獲取696×697×256波段的高光譜圖像立方體，曝光時間設置為12毫秒，入射狹縫寬度為30微米，視場角為22°
，CCD像素為1392*1040。

接種根腐病病菌的辣椒植株為各品種20株

，共40株可用辣椒植株。從每一株植株中選取葉片最大,葉脈清晰的兩片葉子作為研究對象
，共標記了80篇葉片。連續在健康期、潛伏期和患病期對辣椒植株進行高光譜圖像拍攝。將所有辣椒葉片樣本按2:1的比例使用SPXY算法劃分訓練集與預測集，以建立適用範圍更加廣泛的辣椒根腐病早期診斷模型。結果為總樣本集辣椒葉片數據240個

，其中總訓練集辣椒葉片數據162個
，總預測集數據78個。其中健康期、潛伏期和發病期總樣本數均為80個，訓練集均為54個，預測集均為26個。

使用SPA和光譜指數方法進行數據降維
，使用PLS-DA,LSSVM,BP神經網絡等技術進行病害早期預測。

結論：圖1為辣椒在400 ~ 1000 nm範圍內的平均光譜標準差曲線及其對應的葉片狀態的圖片。結果表明兩類辣椒品種樣品的光譜曲線趨勢差異很小，差異主要集中在550-650和750-850的波長區間內。550-650的波長區間內的光譜信息通常用於表征植被特征如植被的生長情況和健康狀況等
，可以用於衡量葉綠素的含量。山椒4haodehuanbingqiyangbendeguangpufanshelvxiajiangqushigengjiamingxian，biaozhenghuanbingqiyangbenheqitayangbenchayixingda
，kenengzaishougenfubingxieposhishiqulegengduodeyelvsu。zai750-850的(de)波(bo)段(duan)區(qu)間(jian)內(nei)，滿(man)山(shan)紅(hong)的(de)潛(qian)伏(fu)期(qi)葉(ye)片(pian)樣(yang)本(ben)的(de)反(fan)射(she)率(lv)有(you)從(cong)低(di)到(dao)高(gao)增(zeng)加(jia)的(de)趨(qu)勢(shi)，潛(qian)伏(fu)期(qi)樣(yang)本(ben)和(he)健(jian)康(kang)期(qi)樣(yang)本(ben)的(de)光(guang)譜(pu)反(fan)射(she)率(lv)差(cha)異(yi)大(da)，和(he)患(huan)病(bing)期(qi)樣(yang)本(ben)的(de)差(cha)異(yi)小(xiao)，可(ke)能(neng)是(shi)潛(qian)伏(fu)期(qi)失(shi)水(shui)較(jiao)多(duo)


，使(shi)這(zhe)一(yi)品(pin)種(zhong)的(de)辣(la)椒(jiao)的(de)健(jian)康(kang)期(qi)葉(ye)片(pian)和(he)潛(qian)伏(fu)期(qi)葉(ye)片(pian)差(cha)別(bie)較(jiao)大(da)。山(shan)椒(jiao)4號(hao)的(de)潛(qian)伏(fu)期(qi)葉(ye)片(pian)樣(yang)本(ben)的(de)光(guang)譜(pu)反(fan)射(she)率(lv)則(ze)更(geng)加(jia)平(ping)緩(huan)。總(zong)體(ti)而(er)言(yan)
，兩(liang)個(ge)品(pin)種(zhong)的(de)辣(la)椒(jiao)光(guang)譜(pu)曲(qu)線(xian)趨(qu)勢(shi)基(ji)本(ben)相(xiang)似(si)，所(suo)以(yi)這(zhe)個(ge)論(lun)文(wen)把(ba)兩(liang)類(lei)辣(la)椒(jiao)品(pin)種(zhong)合(he)並(bing)在(zai)一(yi)起(qi)進(jin)行(xing)建(jian)模(mo)研(yan)究(jiu)。這(zhe)樣(yang)做(zuo)

，能(neng)增(zeng)加(jia)病(bing)害(hai)早(zao)期(qi)檢(jian)測(ce)模(mo)型(xing)的(de)可(ke)移(yi)植(zhi)性(xing)，使(shi)這(zhe)個(ge)研(yan)究(jiu)的(de)實(shi)驗(yan)結(jie)果(guo)可(ke)以(yi)應(ying)用(yong)到(dao)其(qi)他(ta)辣(la)椒(jiao)品(pin)種(zhong)中(zhong)。觀(guan)察(cha)圖(tu)2不難發現健康辣椒葉片和染病辣椒葉片在550nm、580~680 nm和760~1000nm波長範圍內的光譜反射率有明顯的區別。這種顯著性的差異證實了高光譜成像技術在識別葉片病變方麵的有效性（10）。進一步分析表明
，健康期辣椒葉片反射率最高
，發病期的辣椒葉片的光譜反射率明顯低於潛伏期辣椒葉片的光譜反射率。在580~680nm範圍內
，健康期葉片的反射率較低，而且隨著患病程度的加重，該範圍內的反射率逐漸減小。染病（包括潛伏期和發病期）deyepianzaijinhongwaiboduandeguangpufanshelvmingxiandiyujiankangyepiandefanshelv，qiesuizhelajiaoyepianhuanbingchengdudejiashen，jinhongwaiboduanfanweineidefanshelvzhujianjianxiao，zheshiyouyubingjunyinqigenjingxibaojiegoupohuaierdaozhideyepianhanshuiliangheyelvsuhanliangxiajiangzaochengde。

利用SPA從全光譜數據中挑選出*有效的特征波長，減少冗餘變量。如圖3（a）所示，設置特征波段數量範圍為10至32個，經過反複測試
，當變量個數為13時，使用的特征波長組合具有最小的校正均方根誤差（RMSE=0.2797）。最終篩選出的特征波段為：425、433
、437、453、552、639、670

、680
、690、715、741、819
、909nm，占全部波段的5.08%。

將所有辣椒樣本按健康期
、潛伏期和發病期三類進行賦值分類，將辣椒葉片樣本在400nm到1000nm任意兩波段下反射率計算出的NDSI值，利用皮爾森算法（Pearson correlation coefficient）與其對應的染病階段進行相關性分析，將相關係數繪製成分布矩陣，如圖4所示。

紅色區域顯示最高正相關值，藍色區域顯示最高負相關值。560~610nm及附近波長範圍構建出的光譜指數與根腐病顯示出極高的相關性。790~810nm及附近波長範圍構建出的光譜指數與根腐病具有較高的相關性，該區域相關係數絕對值明顯低於560~610nm波段區域，相關係數絕對值（|r|）最高值為0.725。其他區域構建出的光譜指數與根腐病的相關性普遍較低，|r|範圍區間在0到0.3。結果表明，與辣椒根腐病相關的敏感波段主要分布在560~610nm與560~610nm這個波段區域。

在高光譜分類中，BP神經網絡得到了最高的準確率。整個BP神經網絡模型由input，hidden layer、output layer和output組成，其中hidden layer的激活函數設置為tansig，output layer 的激活函數設置為purelin。將訓練參數中的最大迭代次數設置為10000
，目標訓練誤差設置為10-6，學習率設置為0.01。整個BP網絡使用levenberg-marquardt訓練算法
，validation checks設置為6，數據分類方式設置為random。最後經過多次調試，得出了模型的預測集訓練結果如圖5所示。

具體的訓練結果在混淆矩陣圖6所示，圖6（a）是SPA的分類結果
，其中29個健康樣本有3個被誤判為潛伏期，28個潛伏期樣本有2個被誤判為健康期，21個患病期樣本被誤判為潛伏期樣本。圖6（b）是光譜指數的分類結果，其中27個健康樣本有2個被誤判為潛伏期
，29個潛伏期樣本有5個被誤判為健康期，有1個被誤判為患病期。潛伏期和健康期葉片之間的誤判大部分是山椒4號品種的辣椒
，誤判數為5個（共7個），僅僅有兩個是滿山紅品種的辣椒
，或許滿山紅辣椒的潛伏期和健康期的差別比山椒4號更大一些

圖6