文章詳情
基于連續(xù)小波變換與高光譜成像技術(shù)的霉變花生精準(zhǔn)識(shí)別研究（上）
日期：2025-07-19 13:57
瀏覽次數(shù)：115
摘要：
應(yīng)用方向：
本研究結(jié)合高光譜成像技術(shù)和連續(xù)小波變換（CWT），聚焦于霉變花生的精準(zhǔn)識(shí)別，通過(guò)提取光譜敏感特征區(qū)分健康與霉變樣本。這一方法展示了高光譜技術(shù)在食品質(zhì)量與**檢測(cè)中的重要應(yīng)用價(jià)值，尤其在霉變與**感染監(jiān)測(cè)、農(nóng)產(chǎn)品自動(dòng)分選和在線(xiàn)監(jiān)控等方面具有廣泛潛力。相比傳統(tǒng)檢測(cè)方法，高光譜技術(shù)實(shí)現(xiàn)了非接觸、快速、綠色環(huán)保的檢測(cè)方式，為食品**監(jiān)管及智能農(nóng)業(yè)提供了高效解決方案。
背景：
花生作為全球廣泛種植和消費(fèi)的重要油料作物，具有很高的營(yíng)養(yǎng)價(jià)值。然而，由于其特殊的組成成分，花生極易受到**（...
應(yīng)用方向：
本研究結(jié)合高光譜成像技術(shù)和連續(xù)小波變換（CWT），聚焦于霉變花生的精準(zhǔn)識(shí)別，通過(guò)提取光譜敏感特征區(qū)分健康與霉變樣本。這一方法展示了高光譜技術(shù)在食品質(zhì)量與**檢測(cè)中的重要應(yīng)用價(jià)值，尤其在霉變與**感染監(jiān)測(cè)、農(nóng)產(chǎn)品自動(dòng)分選和在線(xiàn)監(jiān)控等方面具有廣泛潛力。相比傳統(tǒng)檢測(cè)方法，高光譜技術(shù)實(shí)現(xiàn)了非接觸、快速、綠色環(huán)保的檢測(cè)方式，為食品**監(jiān)管及智能農(nóng)業(yè)提供了高效解決方案。
背景：
花生作為全球廣泛種植和消費(fèi)的重要油料作物，具有很高的營(yíng)養(yǎng)價(jià)值。然而，由于其特殊的組成成分，花生極易受到**（如黃曲霉）的感染，進(jìn)而產(chǎn)生強(qiáng)致癌性的黃曲霉**（AFB1）。這種**不僅對(duì)人類(lèi)和牲畜健康構(gòu)成嚴(yán)重威脅，還可能導(dǎo)致農(nóng)民、加工商及分銷(xiāo)商的經(jīng)濟(jì)損失。
現(xiàn)有的黃曲霉**檢測(cè)方法（如薄層色譜、高效液相色譜等）盡管準(zhǔn)確，但存在耗時(shí)、操作復(fù)雜、依賴(lài)專(zhuān)業(yè)人員等局限性，不適用于自動(dòng)化檢測(cè)場(chǎng)景。而近紅外光譜及高光譜成像技術(shù)因其非接觸性及高效性，逐漸成為食品質(zhì)量評(píng)估領(lǐng)域的重要工具。然而，高光譜成像數(shù)據(jù)的高維特性給實(shí)時(shí)檢測(cè)系統(tǒng)帶來(lái)了計(jì)算負(fù)擔(dān)，特征降維成為后續(xù)數(shù)據(jù)處理的必要步驟。
連續(xù)小波變換（CWT）作為一種有效的特征提取方法，可在不同波長(zhǎng)和尺度上分解光譜數(shù)據(jù)，從而識(shí)別出細(xì)微的光譜特征。盡管CWT在某些領(lǐng)域已顯示出其優(yōu)越性，但其在谷物和油料作物**感染檢測(cè)中的應(yīng)用仍較少。本研究旨在結(jié)合CWT和高光譜成像技術(shù)，開(kāi)發(fā)一種高效、可靠的霉變花生識(shí)別方法。具體的研究目標(biāo)為：（1）探索霉變花生的光譜特性；（2）確定適用于霉變花生檢測(cè)的*佳CWT特征；（3）比較CWT提取的小波特征與傳統(tǒng)波段選擇方法的分類(lèi)性能。
實(shí)驗(yàn)設(shè)計(jì)
１.1材料與方法
（1）花生樣品制備
從市場(chǎng)上購(gòu)買(mǎi)了三種花生品種：花育（HY）、四粒紅（SLH）和小白沙（XBS）。手工挑選出完好無(wú)損的健康花生，并將其分成兩部分：一部分妥善保存作為健康樣本，另一部分用于獲取霉變花生樣本。為了獲得霉變花生樣本，將花生放入容器中，并置于37°C的恒溫箱中，相對(duì)濕度設(shè)定在85~90%，持續(xù)10天以促進(jìn)曲霉菌屬的快速繁殖。從第11天起，溫度調(diào)整至30°C，相對(duì)濕度保持不變。分別在第20天和第30天取出花生作為不同霉變階段的霉變樣本。然后，對(duì)于HY品種的花生，隨機(jī)選取8個(gè)健康、霉變20天和霉變30天的樣本，并使用黃曲霉**B1（AFB1）快速檢測(cè)試紙進(jìn)行評(píng)估。對(duì)于另外兩種花生品種，每種隨機(jī)選取10個(gè)健康、霉變20天和霉變30天的樣本，并使用AFB1快速檢測(cè)試紙進(jìn)行評(píng)估。AFB1快速檢測(cè)試紙*低檢測(cè)限為5 ppb。該產(chǎn)品采用競(jìng)爭(zhēng)性抑制金**層析法原理。當(dāng)樣本溶液中的 AFB1 濃度超過(guò)檢測(cè)限時(shí)，檢測(cè)線(xiàn)不會(huì)顯現(xiàn)顏色反應(yīng)，結(jié)果判定為陽(yáng)性；當(dāng)樣本溶液中的 AFB1 濃度低于檢測(cè)限時(shí)，檢測(cè)線(xiàn)呈現(xiàn)紫色反應(yīng)，結(jié)果判定為陰性。在實(shí)驗(yàn)過(guò)程中，健康花生樣本經(jīng) AFB1 快速檢測(cè)試紙檢測(cè)，其 AFB1 含量均低于 5 ppb，據(jù)此假設(shè)所有健康花生樣本均為健康。相反，霉變花生樣本經(jīng)檢測(cè)發(fā)現(xiàn) AFB1 含量均超過(guò) 5 ppb，因此假設(shè)所有霉變花生均為霉變。為排除水分含量對(duì)實(shí)驗(yàn)結(jié)果的干擾，健康與霉變的花生仁分別在 60°C 的干燥箱中干燥 24 小時(shí)。
對(duì)于每種花生品種，我們獲取了七張高光譜圖像（三張用于訓(xùn)練，四張用于測(cè)試）。圖1展示了用于獲取這七張高光譜圖像的花生樣本（HY）的照片。前三張圖像顯示了用作訓(xùn)練數(shù)據(jù)的花生樣本，**張展示了健康的花生樣本，**張展示了第20天獲得的霉變花生樣本，*后一張展示了第30天獲得的霉變花生樣本。后四張圖像展示了用作測(cè)試數(shù)據(jù)的花生樣本；其中前兩張是健康花生和第20天獲得的霉變花生的混合物，另外兩張是健康花生和第30天獲得的霉變花生的混合物。*終，共成像了547個(gè)花生樣本（HY：154個(gè)，XBS：175個(gè)，SLH：218個(gè)），包括252個(gè)健康花生樣本和295個(gè)霉變花生樣本。在這些樣本中，232個(gè)花生樣本，包括79個(gè)健康樣本和153個(gè)霉變樣本，被用作訓(xùn)練數(shù)據(jù)；315個(gè)花生樣本，包括173個(gè)健康樣本和142個(gè)霉變樣本，被用作測(cè)試數(shù)據(jù)。
圖1.用于獲取七張高光譜圖像的花生樣本（HY）照片：頂部三張圖像用于訓(xùn)練數(shù)據(jù)集，底部四張圖像用于測(cè)試數(shù)據(jù)集
（2）高光譜成像系統(tǒng)與圖像采集
高光譜圖像采集使用的是GaiaSorter（江蘇雙利合譜科技有限公司）設(shè)備。它由四個(gè)部分組成：一個(gè)光譜成像系統(tǒng)、一個(gè)照明系統(tǒng)、一個(gè)傳送臺(tái)和一個(gè)計(jì)算機(jī)。光譜成像系統(tǒng)由一個(gè)高光譜相機(jī)（Image-λ-N25E-HS）組成，其光譜范圍為920至2530納米，該相機(jī)連接到一個(gè)標(biāo)準(zhǔn)的C口變焦鏡頭（F/2, f = 22.5 mm, HSIA-OLES22）。照明系統(tǒng)由四個(gè)200瓦的溴鎢燈（HSIALS-T-200W）組成。傳送臺(tái)由一個(gè)樣品臺(tái)和集成的電動(dòng)精密傳送臺(tái)（HSIAT500）構(gòu)成。計(jì)算機(jī)安裝了SpecVIEW系統(tǒng)控制軟件。
花生樣本被放置在一個(gè)10厘米×10厘米的黑色托盤(pán)中以獲取高光譜圖像，如圖1所示。所獲得的圖像共有288個(gè)波段，光譜分辨率為5.6納米。此外，為了校正所采集的高光譜圖像，我們收集了一張白參考圖像，并通過(guò)完全關(guān)閉相機(jī)光圈記錄了一張黑參考圖像。
（3）高光譜圖像預(yù)處理與光譜提取
圖像采集完成后，將所得到的高光譜圖像使用黑參考圖像和白參考圖像進(jìn)行校正。另外，在捕獲高光譜圖像時(shí)，隨機(jī)噪聲是不可避免的，因此在進(jìn)行進(jìn)一步數(shù)據(jù)處理之前需要對(duì)數(shù)據(jù)進(jìn)行了濾波。在本研究中，采用了五點(diǎn)平滑濾波器來(lái)處理高光譜數(shù)據(jù)。
非花生像素不僅在后續(xù)數(shù)據(jù)處理中無(wú)用，甚至可能對(duì)分析結(jié)果產(chǎn)生干擾，因此被標(biāo)記為背景并賦予零值。為了提取緊密相連且無(wú)法通過(guò)簡(jiǎn)單閾值分割的花生區(qū)域，采用了標(biāo)記控制的水洗算法生成掩模圖像。隨后，利用該掩模圖像有效去除背景信息。在此基礎(chǔ)上，進(jìn)一步進(jìn)行了預(yù)處理操作，包括圖像裁剪以及低信噪比波段的移除。*終，篩選出波長(zhǎng)范圍為1000-2486nm的 265個(gè)波段用于分析
在圖像空間分割上，采用了區(qū)域生長(zhǎng)算法，其本質(zhì)是將具有相似屬性的像素組合在一起。對(duì)于每個(gè)區(qū)域，指定一個(gè)種子點(diǎn)作為生長(zhǎng)的起點(diǎn)；然后，比較周?chē)南袼攸c(diǎn)以確定它們是否包含在區(qū)域內(nèi)。在本研究中，應(yīng)用了四鄰域區(qū)域生長(zhǎng)算法到掩模圖像上，以連接同一核區(qū)域內(nèi)的像素，并為它們分配一個(gè)序列號(hào)。然后，每個(gè)花生仁都可以通過(guò)序列號(hào)輕松追蹤，這有利于實(shí)現(xiàn)花生分選的工業(yè)自動(dòng)化。
對(duì)于健康的花生樣本，我們從每種花生的訓(xùn)練圖像中隨機(jī)選取了100個(gè)像素。對(duì)于發(fā)霉的花生，首先從訓(xùn)練圖像中手動(dòng)選擇發(fā)霉的像素作為感興趣區(qū)域（ROI）；然后，從每種花生的ROI中隨機(jī)選取50個(gè)像素。因此，共計(jì)獲得了600個(gè)花生像素樣本，其中包括300個(gè)健康花生像素和300個(gè)不同程度**感染的發(fā)霉花生像素。所選取的花生像素樣本如圖2所示。
圖2.選定花生像素樣本的分布情況
（4）方法
特征波長(zhǎng)提取：
連續(xù)小波變換（CWT）是一種強(qiáng)大的方法，用于在不同尺度和分辨率下檢測(cè)和分析高光譜數(shù)據(jù)中的微弱信號(hào)。本質(zhì)上，CWT是一種線(xiàn)性操作，通過(guò)使用母小波函數(shù)在不同的波長(zhǎng)和尺度上，可以將高光譜反射光譜轉(zhuǎn)換為一系列系數(shù)。在本研究中，由于吸收特征的形狀類(lèi)似于準(zhǔn)高斯函數(shù)，因此選擇了墨西哥帽作為母小波基。為了減輕計(jì)算負(fù)擔(dān)，僅保留了二進(jìn)制尺度（21, 22, 23, …, 和 27）上的小波功率；因?yàn)榭偛ǘ螖?shù)為267，所以大于27 = 128的尺度上的分解分量不再攜帶有意義的光譜信息。
杰弗里斯-馬圖西塔（J-M）距離是一種靈活直觀(guān)的特征可分性指數(shù)，用于衡量?jī)蓚€(gè)不同類(lèi)別之間的特征，并廣泛用于指導(dǎo)特征選擇。杰弗里斯-馬圖西塔距離（J-M距離）的范圍是0到2，它提供了兩個(gè)類(lèi)別之間可分性的一般度量。J-M距離越大，意味著兩個(gè)類(lèi)別之間分離的概率越高，反之亦然。在本研究中，J-M距離被用來(lái)評(píng)估不同尺度和波長(zhǎng)下的小波系數(shù)，以進(jìn)行小波特征（WF）的選擇。
為了識(shí)別發(fā)霉花生仁，需要一種特征選擇方法來(lái)識(shí)別*重要的特征，因?yàn)樵S多WFs由于連續(xù)分解是冗余的，因此要選擇有意義的小波特征。首先，使用Mexican Hat作為母小波，通過(guò)CWT對(duì)每個(gè)花生仁的原始高光譜反射數(shù)據(jù)進(jìn)行分解，分解尺度為21、22、23，……，27。這樣，原始光譜就被轉(zhuǎn)換成了在不同波長(zhǎng)和尺度上的一組小波系數(shù)。其次，通過(guò)計(jì)算健康和發(fā)霉花生樣本之間的J-M距離，構(gòu)建了J-M距離尺度圖。*后，根據(jù)J-M距離降序排列特征，并應(yīng)用一個(gè)閾值J-M距離來(lái)劃分前1%的特征。由閾值劃分的特征在J-M距離尺度圖上形成了一個(gè)分散的特征區(qū)域。理論上，所有在特征區(qū)域內(nèi)的特征都是被選擇的小波特征。然而，它們攜帶了冗余的光譜信息，因?yàn)橥粎^(qū)域內(nèi)的特征是在連續(xù)的波長(zhǎng)位置和尺度上產(chǎn)生的。因此，對(duì)于每個(gè)特征區(qū)域，使用具有*大J-M距離的特征來(lái)代表該特征區(qū)域捕獲的光譜信息。
SPA是一種前向變量選擇算法，通常用于減少模型中的變量數(shù)量，提高建模的速度和效率。它從單一波長(zhǎng)開(kāi)始，然后在每次迭代中加入一個(gè)新的波長(zhǎng)，直到達(dá)到指定的波長(zhǎng)數(shù)量。新選擇的變量是在之前選擇的變量正交子空間上具有*大投影值的變量。因此，SPA提取的有用變量子集具有*小的共線(xiàn)性。*優(yōu)的變量數(shù)量通過(guò)均方根誤差（RMSE）來(lái)確定。在我們的研究中，SPA被用來(lái)選擇*優(yōu)的波段以與WFs進(jìn)行比較。
分類(lèi)模型：
我們將健康花生標(biāo)記為1，發(fā)霉花生標(biāo)記為2。為了評(píng)估WFs以及SPA選擇的*優(yōu)波段的性能，我們使用了兩種分類(lèi)模型：偏*小二乘判別分析（PLS-DA）和支持向量機(jī)（SVM）。通過(guò)計(jì)算健康花生和發(fā)霉花生的整體分類(lèi)準(zhǔn)確率來(lái)評(píng)估WFs的性能。在PLSR中，由于預(yù)測(cè)響應(yīng)很少直接產(chǎn)生等于1和2的二元結(jié)果，而是接近1或2的結(jié)果，因此使用了1.5作為分類(lèi)的閾值。此外，PLS-DA模型中理想的潛在變量數(shù)量是通過(guò)五折交叉驗(yàn)證下RMSE的值來(lái)確定的。在SVM中，徑向基函數(shù)（RBF）核被用作SVM分類(lèi)框架的核。在評(píng)估方法中，通過(guò)計(jì)算靈敏度、特異性和總體準(zhǔn)確性（OA）來(lái)評(píng)估分類(lèi)模型的性能。
為了更清晰地說(shuō)明圖像處理步驟，圖3展示了一個(gè)流程圖。虛線(xiàn)框中的部分專(zhuān)門(mén)針對(duì)訓(xùn)練數(shù)據(jù)，用于確定*佳波段和特征權(quán)重。驗(yàn)證圖像則跳過(guò)此部分，直接使用與訓(xùn)練圖像相同的*佳波段和特征權(quán)重。
圖3.圖像處理流程圖
１.2.結(jié)果與討論
（1）花生的光譜反射率
圖4展示了健康和霉變花生的反射光譜，這些光譜是通過(guò)計(jì)算不同類(lèi)別樣本光譜的平均值得出的。在健康花生的反射光譜中可以看到六個(gè)顯著的吸收峰，分別位于大約1208納米、1472納米、1747納米、1938納米、2145納米和2329納米（用黑點(diǎn)標(biāo)記），主要可以歸因于與花生中的水分、蛋白質(zhì)和油脂成分相關(guān)的CH、NH和OH的吸收。在1000-1500納米的波長(zhǎng)范圍內(nèi)，健康和霉變花生的光譜曲線(xiàn)形狀明顯不同。當(dāng)花**生霉變時(shí)，1118納米處的反射率降低，導(dǎo)致吸收峰減弱。在1365-2486納米的波長(zhǎng)區(qū)域，霉變花生的反射值高于健康花生，這可能歸因于**污染引起的散射和吸收特性。**的侵入可以使花生仁胚乳變得多孔，這可能導(dǎo)致霉變花生比健康花生散射更多的光線(xiàn)，這種散射會(huì)導(dǎo)致在反射模式下吸收更少的近紅外輻射。
圖4.健康與霉變花生的光譜響應(yīng)：黑點(diǎn)表示健康花生的吸收峰，垂直虛線(xiàn)表示由SPA選擇的*佳波長(zhǎng)變量。
