引言

2019 年我國(guó)食用菌總產(chǎn)量達(dá)3934萬(wàn)t，總產(chǎn)值高達(dá)3127億元。雙孢蘑菇是食用菌最主要的栽培品種之一，在全球60多個(gè)國(guó)家和地區(qū)都有種植生產(chǎn)。雙孢蘑菇的短生長(zhǎng)周期、快繁殖、大數(shù)量等特點(diǎn)使得其經(jīng)濟(jì)價(jià)值廣受菇農(nóng)認(rèn)可。但在菇房種植過(guò)程中，由于機(jī)械化程度低、管理欠妥等問(wèn)題，雙孢蘑菇不可避免會(huì)遭受病害的侵襲。雙孢蘑菇真菌病疣孢霉病，又稱(chēng)褐腐病、濕泡病，是危害雙孢蘑菇的世界性病害。有害疣孢霉菌是該病的致病真菌，在我國(guó)長(zhǎng)江以南的種植區(qū)域廣為流行，嚴(yán)重影響著我國(guó)雙孢蘑菇菇房的生產(chǎn)效益。由于雙孢蘑菇疣孢霉病與寄主、病原菌、傳播媒介及菇房培養(yǎng)環(huán)境條件均有關(guān)系，且發(fā)病潛伏期長(zhǎng)、傳染快，存在使健康菇無(wú)法正常成長(zhǎng)的問(wèn)題，有經(jīng)驗(yàn)的菇農(nóng)最快也要在雙孢蘑菇子實(shí)體小菇期后才能鑒別染病與否；或者采用內(nèi)源轉(zhuǎn)錄間隔區(qū)（ITS）基因片段的聚合酶鏈?zhǔn)椒磻?yīng)（PCR）鑒定以及傳統(tǒng)的柯赫氏法則待菇體出現(xiàn)損耗后進(jìn)行檢測(cè)。這些方法往往因?yàn)闀r(shí)效性差導(dǎo)致病菇傳染速度大于遏制速度。因此，在早期對(duì)雙孢蘑菇疣孢霉病進(jìn)行快速、無(wú)損的檢測(cè)識(shí)別，能夠幫助菇農(nóng)及時(shí)發(fā)現(xiàn)菇房染病雙孢蘑菇，提醒菇農(nóng)及時(shí)采取止損措施，阻止病害繼續(xù)擴(kuò)散，從而保證菇房的雙孢蘑菇產(chǎn)量與品質(zhì)，保障菇農(nóng)的正 常經(jīng)濟(jì)收益。

國(guó)內(nèi)外學(xué)者采用無(wú)損檢測(cè)技術(shù)對(duì)雙孢蘑菇內(nèi)部品質(zhì)開(kāi)展了不同研究，高光譜成像技術(shù)能同時(shí)獲取雙孢蘑菇的光譜和圖像信息，具有信息量豐富、檢測(cè)便捷、無(wú)損等優(yōu)點(diǎn)，被廣泛應(yīng)用于農(nóng)作物病害檢測(cè)。例如 Parrag 等利用高光譜成像技術(shù)檢測(cè)雙孢蘑菇子實(shí)體成熟期的蜘蛛網(wǎng)病（由樹(shù)枝狀真菌引起），利用SVM建模成功分離出蜘蛛網(wǎng)病感染樣本，每組樣本的正確率均在80%以上；Gaston等利用可見(jiàn)-近紅外高光譜成像技術(shù)對(duì)雙孢蘑菇子實(shí)體成熟期菌蓋上的褐斑病 （細(xì)菌性病，無(wú)害）進(jìn)行識(shí)別研究，對(duì)正常、機(jī)械 損傷及在運(yùn)輸或儲(chǔ)存時(shí)被微生物腐敗后的雙孢蘑菇建立識(shí)別模型，使用決策樹(shù)將雙孢蘑菇分類(lèi)，分類(lèi)的準(zhǔn)確度能在95%以上。結(jié)果表明高光譜可以用于檢測(cè)并對(duì)蘑菇不同類(lèi)型的損傷進(jìn)行快速、自動(dòng)的分類(lèi)。

現(xiàn)已有很多基于高光譜成像技術(shù)對(duì)雙孢蘑菇子實(shí)體成熟期時(shí)含水率、新鮮度和 損傷情況的預(yù)測(cè)研究，但目前針對(duì)雙孢蘑菇病害研究都是在子實(shí)體生長(zhǎng)的中后期，原基期至小菇期時(shí)的早期雙孢蘑菇可能處在菌絲生長(zhǎng)階段，其菌蓋無(wú)法用肉眼辨別。故采用高光譜成像技術(shù)獲取早期 雙孢蘑菇各個(gè)生長(zhǎng)周期的病害光譜信息，以實(shí)現(xiàn)雙孢蘑菇病害的早期診斷。本研究擬基于高光譜成像技術(shù)，在雙孢蘑菇生長(zhǎng)早期 （小菇期之前）對(duì)雙孢蘑菇疣孢霉病進(jìn)行檢測(cè)識(shí)別，能夠發(fā)現(xiàn)染病雙孢蘑菇的同時(shí)，還能夠減少農(nóng)藥使用、農(nóng)藥超標(biāo)、土壤污染等環(huán)境問(wèn)題。為消除高光譜圖像數(shù)據(jù)在采集過(guò)程中，可能存在光譜曲線(xiàn) 錯(cuò)誤旋轉(zhuǎn)、基線(xiàn)漂移和圖像變形等不利因素造成的噪聲信號(hào)大而過(guò)擬合，本研究對(duì)采集的高光譜數(shù)據(jù)進(jìn)行預(yù)處理，并選擇最優(yōu)預(yù)處理算法。對(duì)比不同建模方法結(jié)合預(yù)處理方法，通過(guò)定性評(píng)價(jià)標(biāo)準(zhǔn)選擇最終的病害鑒別模型，為開(kāi)發(fā)雙孢蘑菇疣孢霉病早期 鑒別設(shè)備提供參考。

2、材料與方法

2.1 實(shí)驗(yàn)材料

雙孢蘑菇菌菌種為W192，疣孢霉菌種為My.p0012，均由福建省農(nóng)業(yè)科學(xué)院食用菌研究所提供。

2.2 方法

2.2.1 雙孢蘑菇的培養(yǎng)

將播種、裝袋后的雙孢蘑菇放入人工氣候培養(yǎng)箱中培養(yǎng)，初始設(shè)置溫度22 ℃，相對(duì)濕度90%，無(wú)光照。待菌絲能夠長(zhǎng)滿(mǎn)90%栽培 料時(shí)（約兩周）進(jìn)行覆土處理，以福建省農(nóng)業(yè)科學(xué)院食用菌研究所的土料覆蓋栽培料約2.5cm為覆土層高。參考郭倩等研究的成果，此時(shí)調(diào)整人工氣候培養(yǎng)箱內(nèi)的環(huán)境為溫度20℃，濕度 85%，無(wú)光照。待菌絲長(zhǎng)到2/3覆土層時(shí)（約覆土后7～10d）， 調(diào)整人工氣候培養(yǎng)箱內(nèi)的環(huán)境為溫度20℃,相對(duì)濕度90%，無(wú)光照，等待出菇。培養(yǎng)期間每天早晚通風(fēng)1h。

2.2.2 雙孢蘑菇染病接種

將有害疣孢霉菌在馬鈴薯葡萄糖瓊脂培養(yǎng)基（PDA）上,于25℃ 恒溫暗箱中培養(yǎng)7d。于無(wú)菌操作臺(tái)內(nèi)用超純水將PDA培養(yǎng)基上的有害疣孢霉菌孢子洗下，制作成孢子懸液，并用血球計(jì)數(shù)板檢查 有害疣孢霉菌孢子懸液為1.0×105個(gè)mL?1。待雙孢蘑菇菌絲長(zhǎng)到2/3覆土層時(shí)，將5mL 有害疣孢霉菌孢子懸液均勻噴灑在覆土層表面。

2.2.3 病菇鑒定

在無(wú)菌臺(tái)內(nèi)將病菇上的褐色液滴和病菇切片分別放置在PDA平板培養(yǎng)基內(nèi)并封口后，于25 ℃恒溫培養(yǎng)箱中培養(yǎng)。5d后，在顯微鏡 下挑出菌絲并純化，觀(guān)察是否有發(fā)黃、凋零現(xiàn)象。也可通過(guò)菌蓋大小來(lái)判斷，處在同一時(shí)期的雙孢蘑菇，與健康樣本的菌蓋直徑相比，染病樣本的菌蓋直徑較短。

2.2.4 高光譜圖像獲取

高光譜成像系統(tǒng)（見(jiàn)圖1）主要包括內(nèi)置推掃高光譜成像儀、4 個(gè) 50W鹵鎢燈光源、樣品升降臺(tái)及平板電腦等。為保證采集的高光譜圖像清晰，避免失真現(xiàn)象，調(diào)試確定物鏡距離、CCD相機(jī)曝光時(shí)間和推掃移動(dòng)速度。經(jīng)過(guò)反復(fù)調(diào)試，最終確定物鏡高度為 30cm，CCD 相機(jī)曝光時(shí)間為70ms，相機(jī)推掃移動(dòng)速度為 1.30 mm·s?1。分別采集健康雙孢蘑菇早期子實(shí)體菌蓋的可見(jiàn)/近紅外高光譜圖像（401～1046 nm，分辨率2.8nm）后提取感 興趣區(qū)域（Region of interest，ROI）的平均光譜數(shù)據(jù)。由于感興趣區(qū)域?yàn)榫w周?chē)糠郑矢信d趣區(qū)域提取形狀選擇橢圓形，逐一選取圖像樣本的ROI后利用新區(qū)域（New Region）疊加待提取數(shù)據(jù)的樣本圖像區(qū)域，最后在401～1046 nm 波段下提取所有圖像感興趣區(qū)域的光譜數(shù)據(jù)。為清除由于傳感器暗電流等原因產(chǎn)生圖像中的噪音，將采集到的原始高光譜圖像利用標(biāo)準(zhǔn)反射白板和黑暗背景下高光譜圖像數(shù)據(jù)校正，如公式 1 所示

I = I0? B /W ? B (1)

式中：I——校正之后的高光譜圖像數(shù)據(jù)；I0——原始的高光譜圖像數(shù)據(jù)；B——黑暗背景下的高光譜圖像數(shù)據(jù)；W——標(biāo)準(zhǔn)反射白板的高光譜圖像數(shù)據(jù)。

2.2.5 診斷模型

隨機(jī)森林分類(lèi)器（RF）是一種集成分類(lèi)器，對(duì)于異常值和噪聲有良好的容忍度，隨機(jī)選擇的訓(xùn)練樣本子集和變量后，生成多個(gè)決策 樹(shù)，其具體執(zhí)行步驟見(jiàn)參考文獻(xiàn)。本研究隨機(jī)森林通過(guò)設(shè)定500 顆決策樹(shù)對(duì)分類(lèi)結(jié)果進(jìn)行投票，例如，當(dāng)對(duì)健康樣本的投票數(shù)高于對(duì)染病樣本的投票數(shù)時(shí)，隨機(jī)森林就將這個(gè)變量確定為健康樣本。支持向量機(jī)（SVM）的基本原理是將原始特征空間映射到更高維度的空間，使得數(shù)據(jù)線(xiàn)性可分，這種轉(zhuǎn)換是通過(guò)對(duì)原始數(shù)據(jù)應(yīng)用內(nèi)核函數(shù)來(lái)隱式執(zhí)行的。本研究采用徑向基函數(shù)（RBF）核函數(shù)結(jié)合 網(wǎng)格法優(yōu)化參數(shù)，使建模結(jié)果最佳。對(duì)于二分類(lèi)的分類(lèi)目標(biāo)是尋找最優(yōu)超平面，得到正確分類(lèi)結(jié)果并使分類(lèi)間隔最大。極限學(xué)習(xí)機(jī)（ELM）由Huang等提出，優(yōu)化了傳統(tǒng)的單隱含層前饋神經(jīng)網(wǎng)絡(luò)。不同于傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)，ELM 無(wú)需迭代調(diào)整參數(shù)，能夠自動(dòng)實(shí)現(xiàn)參數(shù)調(diào)優(yōu)。ELM隨機(jī)分配隱含層的權(quán)值和偏差值，通過(guò)解簡(jiǎn)單的方程式確定最終的輸出權(quán)值 ，計(jì)算速度快。

圖 1 可見(jiàn)/近紅外高光譜成像系統(tǒng)

可以由算法隨機(jī)設(shè)定輸入層和隱含層的連接權(quán)值以及隱含層的閾值，在設(shè)定完成后，不 需要反向調(diào)整。在保證學(xué)習(xí)精度的前提下，泛化性能滿(mǎn)足需求。通過(guò)定性評(píng)價(jià)指標(biāo)：樣本識(shí)別準(zhǔn)確度（見(jiàn)公式2），針對(duì)分類(lèi)問(wèn)題判斷所建模型對(duì)樣本分類(lèi)的效果。準(zhǔn)確度越高，該模型對(duì)輸入樣本的分類(lèi)完成的越好。

樣本識(shí)別準(zhǔn)確度 = 正確分類(lèi)樣本數(shù)輸入樣本總數(shù) ×100% (2）

3、結(jié)果與分析

3.1 平均光譜曲線(xiàn)分析

采集原基期、菇蕾期、幼菇期、小菇期的健 康、染病雙孢蘑菇高光譜圖像，各個(gè)生長(zhǎng)周期約1～2d。小菇期后染病雙孢蘑菇病害特征較為明顯， 菌蓋表面出現(xiàn)褐變等患病癥狀，因此將接種后7～ 11d 內(nèi)的樣本認(rèn)定為早期。獲取菌蓋部分的ROI高光譜圖像，最終獲得共400個(gè)高光譜圖像的ROI光譜數(shù)據(jù)樣本，具體見(jiàn)表 1。

表1樣本信息

本研究采用401～1046nm 波長(zhǎng)范圍，共360個(gè)波段做進(jìn)一步研究。健康和染病雙孢蘑菇高光譜平均光譜曲線(xiàn)如圖2所示，健康和染病雙孢蘑菇的高光譜曲線(xiàn)輪廓在原基期和菇蕾期較為相似，隨著蘑菇的生長(zhǎng)，在幼菇期后健康雙孢蘑菇能夠長(zhǎng)出光滑白凈的子實(shí)體，其對(duì)應(yīng)反射率逐漸提高，而染病雙孢蘑菇從小菇期開(kāi)始發(fā)病癥狀顯現(xiàn)，菌蓋表面由于有害疣孢霉菌的侵染產(chǎn)生不規(guī)則褐斑。

圖2健康和染病雙孢蘑菇高光譜平均光譜曲線(xiàn)

染病雙孢蘑菇光譜反射率與健康雙孢蘑菇光譜反射率的差異主要體現(xiàn)在401～450nm 的紫光波段附近和800～1046nm的近紅外波段范圍。400～500nm是紫光波段范圍，此時(shí)物質(zhì)反射率越高，則對(duì)紫光的吸收就越多，呈現(xiàn)出的黃色就越深，發(fā)病癥狀越顯著。由于水的O-H在960nm 處的合頻吸收，C-H 在800～900nm區(qū)域的第三倍頻吸收，N-H在1000～1100nm區(qū)域的第二倍頻在這個(gè)區(qū)域的吸收比較弱，表明有害疣孢霉菌對(duì)雙孢蘑菇的侵染還使得其內(nèi)部含水率發(fā)生了變化。

3.2 預(yù)處理

如圖3可見(jiàn)SG卷積一階求導(dǎo)處理后，消除干擾的同時(shí)在1000～1100nm范圍內(nèi)引入了一定噪聲；SG 卷積平滑處理后，光譜曲線(xiàn)較為平滑，噪聲信號(hào)有所降低；多元散射校正（MSC）處理后，光散射和基線(xiàn)漂移等問(wèn)題得到明消除。

3.3 鑒別模型

將預(yù)處理后的高光譜數(shù)據(jù)用于建立基于高光譜成像技術(shù)的雙孢蘑菇疣孢霉病早期檢測(cè)模型，將樣本按 2∶1的比例劃分為測(cè)試集和預(yù)測(cè)集，得到測(cè)試集樣本共267個(gè)，預(yù)測(cè)集樣本共133個(gè)，最終定性評(píng)價(jià)的結(jié)果如表2所示。其中，結(jié)果最好的為MSC-SVM 模型，測(cè)試集和預(yù)測(cè)集總體樣本識(shí)別準(zhǔn)確度為 92.21%和91.04%。在無(wú)預(yù)處理的全波段建模結(jié)果 中，測(cè)試集樣本總體識(shí)別準(zhǔn)確度在85%以上，預(yù)測(cè) 集樣本總體識(shí)別準(zhǔn)確度在87%以上。

表2不同預(yù)處理與不同建模方法結(jié)果

在經(jīng)過(guò)MSC法 對(duì)原始光譜數(shù)據(jù)進(jìn)行預(yù)處理后RF建模預(yù)測(cè)集準(zhǔn)確度由88.06%提升至89.34%，SVM 建模預(yù)測(cè)集準(zhǔn)確度由87.38%提升至91.04%，SG 卷積一階求導(dǎo)后ELM 建模預(yù)測(cè)集準(zhǔn)確度由89.62% 提升至90.93%。表明MSC能夠有效去除冗余信息，提高檢測(cè)準(zhǔn)確度。經(jīng)過(guò)MSC預(yù)處理后，SVM 建模的測(cè)試集和預(yù)測(cè)集染病樣本識(shí)別準(zhǔn)確度較RF和 ELM建模結(jié)果好，分別提至91.71%和 91.52%，這可為在線(xiàn)監(jiān)測(cè)設(shè)備的研究提供依據(jù)。模型預(yù)測(cè)集分類(lèi)結(jié)果如圖4～6所示。

圖3雙孢蘑菇高光譜曲線(xiàn)

圖4隨機(jī)森林分類(lèi)結(jié)果

圖5SVM分類(lèi)結(jié)果

圖6ELM分類(lèi)結(jié)果

由圖4可以觀(guān)察到，500顆決策樹(shù)對(duì)每個(gè)樣本的投票結(jié)果都往兩端聚集，僅有少量樣本投票結(jié)果在誤差框（虛線(xiàn)框）內(nèi)，表明樣本特征較明顯，隨機(jī)森林能夠?qū)颖咎卣魈崛?，并以此作為分?lèi)依據(jù)，得到良好的分類(lèi)結(jié)果。對(duì)比不同模型的預(yù)測(cè)集總體樣本識(shí)別準(zhǔn)確度，MSC預(yù)處理過(guò)后的隨機(jī)森林分類(lèi)模型比無(wú)預(yù)處理、SG卷積平滑和SG卷積一階求導(dǎo)的準(zhǔn)確度提升的百分比多，說(shuō)明MSC預(yù)處理能夠消除雙孢蘑菇菌蓋表面的顆粒大小及其分布不均所產(chǎn)生的散射影響。如圖5所示，在400個(gè)分類(lèi)樣本中，僅有少量樣本被用RBF核函數(shù)建模的SVM分類(lèi)模型誤判類(lèi)別。當(dāng)全波長(zhǎng)分類(lèi)模型的準(zhǔn)確率較RF分類(lèi)器提升的同時(shí)，SVM分類(lèi)模型測(cè)試集準(zhǔn)確度也維持在87% 以上。其中，無(wú)預(yù)處理的SVM診斷模型預(yù)測(cè)集總體樣本準(zhǔn)確度可達(dá)到87.38%，經(jīng)MSC預(yù)處理后SVM 診 斷模型的預(yù)測(cè)集總體樣本準(zhǔn)確度也由87.38% 提高到91.04%，較SG 卷積平滑和SG卷積一階求導(dǎo)提升的多，分別為90.38%和90.70%。SVM采用徑向基函數(shù)（RBF）核函數(shù)，用網(wǎng)格搜索法自動(dòng)選擇優(yōu)化參數(shù)，使得診斷模型準(zhǔn)確度最佳。

結(jié)論

本研究應(yīng)用可見(jiàn)/近紅外高光譜成像技術(shù)，以雙孢蘑菇疣孢霉病為研究對(duì)象，對(duì)該病害早期檢測(cè)方法進(jìn)行研究，探討了MSC、SG卷積一階求導(dǎo)和SG卷積平滑3種預(yù)處理方法對(duì)建模結(jié)果的影響，得知經(jīng)過(guò) MSC預(yù)處理后的建模結(jié)果中，無(wú)關(guān)信息被有效消除，光譜信息的信噪比提高，為提升雙孢蘑菇疣孢霉病早期診斷模型的準(zhǔn)確度奠定基礎(chǔ)。對(duì)比預(yù)處理后的全波段高光譜圖像數(shù)據(jù)經(jīng)過(guò)RF、SVM和ELM 等3種建模方法對(duì)染病雙孢蘑菇的鑒別效果影響，結(jié)果表明，3種建模方法均有效提升診斷模型的準(zhǔn)確度。其中，經(jīng)過(guò)MSC預(yù)處理和SVM建模方法組合的模型，測(cè)試集和預(yù)測(cè)集總體樣本識(shí)別準(zhǔn)確度為 92.21%和91.04%，優(yōu)于其他預(yù)處理和建模方法的 組合。該研究結(jié)果為進(jìn)一步開(kāi)發(fā)雙孢蘑菇疣孢霉病 的快速無(wú)損檢測(cè)設(shè)備提供了理論基礎(chǔ)。

歡迎關(guān)注公眾號(hào)：萊森光學(xué)，了解更多光譜知識(shí)。

萊森光學(xué)（深圳）有限公司是一家提供光機(jī)電一體化集成解決方案的高科技公司，我們專(zhuān)注于光譜傳感和光電應(yīng)用系統(tǒng)的研發(fā)、生產(chǎn)和銷(xiāo)售。

審核編輯黃宇