引言

當(dāng)前，蔬菜中含有農(nóng)藥殘留的可能性越來(lái)越大，而相對(duì)成熟的農(nóng)藥殘留檢測(cè)技術(shù)如氣相色譜法等存在費(fèi)時(shí)費(fèi)力，破壞性大等明顯缺點(diǎn)，因此，急需一種快速無(wú)損的蔬菜農(nóng)藥殘留檢測(cè)技術(shù)。高光譜成像技術(shù)能夠區(qū)分地物光譜的微小差別，已有研究人員將其應(yīng)用于食品安全監(jiān)測(cè)領(lǐng)域，并且取得了一定的研究成果，在蔬菜農(nóng)藥殘留檢測(cè)方面具有較大的潛力。

異常檢測(cè)是高光譜遙感最常見的應(yīng)用之一，已經(jīng)在軍事目標(biāo)識(shí)別、偽裝目標(biāo)區(qū)分等方面取得了較好的成果。鑒于異常檢測(cè)在高光譜遙感領(lǐng)域目標(biāo)識(shí)別中的有效性，本文將其應(yīng)用于高光譜圖像，試圖探究其檢測(cè)蔬菜葉片上的農(nóng)藥點(diǎn)的可行性。

在本文中，我們首先搭建了一個(gè)高光譜成像平臺(tái)，獲取了包含農(nóng)藥點(diǎn)的菠菜葉片的高光譜圖像數(shù)據(jù)，對(duì)菠菜葉片的光譜特性和不同形態(tài)表現(xiàn)的農(nóng)藥點(diǎn)的光譜特性進(jìn)行了分析，確定了農(nóng)藥點(diǎn)光譜的特征波長(zhǎng)。

一、蔬菜農(nóng)藥點(diǎn)檢測(cè)的高光譜數(shù)據(jù)采集

1.1 蔬菜農(nóng)業(yè)點(diǎn)檢測(cè)的高光譜成像系統(tǒng)

面向蔬菜農(nóng)藥點(diǎn)檢測(cè)的高光譜成像系統(tǒng)由軟件和硬件兩部分組成。軟件部分主要是處理高光譜數(shù)據(jù)的軟件，常見的有ENVI、MATLAB等。硬件部分一般包括光源、分光器、相機(jī)、圖像采集系統(tǒng)和計(jì)算機(jī)等。根據(jù)成像方式的不同，目前在硬件設(shè)備搭建方面主要有兩種類型：一種是使用推掃成像方式實(shí)現(xiàn)，另一種是利用液晶可調(diào)濾光片加數(shù)碼相機(jī)拍照實(shí)現(xiàn)。

推掃成像方式中硬件設(shè)備主要包括光源、推掃的高光譜成像系統(tǒng)、鏡頭、傳送裝置和計(jì)算機(jī)等。光源一般采用鹵素?zé)艋蛘?a target="_blank">LED燈，能夠產(chǎn)生可見光到紅外光的連續(xù)光。推掃的高光譜成像系統(tǒng)用來(lái)記錄物體一個(gè)掃描線上的多個(gè)波段的入射光輻射能量，通常由前端的光譜儀和后端的光電傳感器兩部分組成。鏡頭用來(lái)調(diào)整光的范圍。傳送裝置用來(lái)輸送樣品，并實(shí)現(xiàn)在空間維度的掃描。計(jì)算機(jī)則用來(lái)存儲(chǔ)三維高光譜“立方體”數(shù)據(jù)，并實(shí)現(xiàn)數(shù)據(jù)的實(shí)時(shí)處理。圖1是該類設(shè)備的典型案例的示意圖。

圖 1 推掃型高光譜成像系統(tǒng)結(jié)構(gòu)原理圖

基于液晶可調(diào)濾光片(Liquid Crystal Tunable Filter, LCTF)實(shí)現(xiàn)的高光譜成像系統(tǒng)，不使用傳送帶，而是相機(jī)與樣品相對(duì)固定，利用LCTF和數(shù)碼相機(jī)，直接對(duì)樣品進(jìn)行拍照，通過(guò)連續(xù)采集一系列波長(zhǎng)下的二維圖像，得到三維高光譜數(shù)據(jù)，其結(jié)構(gòu)原理如圖2。該類系統(tǒng)其他構(gòu)成部件還包括光源、相機(jī)、計(jì)算機(jī)等。在文獻(xiàn)中，研究人員還利用這樣思路構(gòu)建了一個(gè)手持式的高光譜成像設(shè)備，實(shí)現(xiàn)對(duì)食品工廠食品設(shè)備的探測(cè)。

圖 2 基于液晶可調(diào)濾光片型高光譜成像系統(tǒng)示意圖

1.2 高光譜圖像獲取

獲取含有農(nóng)藥點(diǎn)的蔬菜高光譜圖像數(shù)據(jù)，首先需要確定蔬菜的種類和農(nóng)藥點(diǎn)的類型。

菠菜是一種一年生或二年生的草本植物，富含維生素A、C 及礦物質(zhì)，尤其是維生素A、C 和鐵的含量是所有蔬菜中最高的，對(duì)于胃腸障礙、痛風(fēng)、神經(jīng)類疾病、便秘、皮膚病和貧血有特殊食療效果,營(yíng)養(yǎng)價(jià)值豐富。隨著生長(zhǎng)環(huán)境的惡化，菠菜的病蟲害也越來(lái)越嚴(yán)重，絕大多數(shù)菠菜都需要多次噴灑農(nóng)藥后才能成熟上市，因此，一般菠菜葉上都含有農(nóng)藥殘留。考慮到菠菜春、秋、冬三季均可生長(zhǎng)，獲取相對(duì)容易，我們選擇了菠菜作為研究樣品。

毒死蜱，又名氯吡硫磷、氯蜱硫磷，，化學(xué)名稱為O,O-二乙基-O-(3,5,6-三氯-2-吡啶基)硫代磷酸，為白色結(jié)晶，是一種有機(jī)磷農(nóng)藥，對(duì)蔬菜、果樹、小麥中的多種咀嚼式和刺吸式口器害蟲防治效果較好，易與土壤中的有機(jī)質(zhì)結(jié)合，對(duì)土壤中的害蟲有特效。有機(jī)磷農(nóng)藥在人體內(nèi)易與膽堿酯酶發(fā)生化學(xué)反應(yīng)產(chǎn)生磷酰化膽堿酯酶，從而使膽堿酯酶喪失分解乙酰膽堿的作用，引發(fā)毒蕈堿樣等神經(jīng)類疾病。鑒于毒死蜱對(duì)較多害蟲具有較好的防治效果及對(duì)人體有較大的危害，我們?cè)趯?shí)驗(yàn)中選擇毒死蜱作為研究對(duì)象。

為了獲取含有農(nóng)藥的蔬菜高光譜圖像，我們搭建了一個(gè)簡(jiǎn)易的獲取高光譜圖像數(shù)據(jù)的實(shí)驗(yàn)平臺(tái)，主要的實(shí)驗(yàn)設(shè)備及材料如下：

高光譜成像儀：iSpecField-NIR/WNIR 便攜式地物光譜儀。

樣品：就近市場(chǎng)購(gòu)買的新鮮菠菜，有效成分為45%的毒死蜱農(nóng)藥。

光源：兩盞功率為90W 的白熾燈，一盞功率為90W 的熒光燈。由于我們所搭建的實(shí)驗(yàn)平臺(tái)比較簡(jiǎn)易，而專業(yè)的光源獲取又相對(duì)困難，熒光燈在可見光波段的輻射能量較強(qiáng)，白熾燈在近紅外波段的輻射能量較強(qiáng)，這兩種類型的光源比較容易獲取，因此，我們?cè)趯?shí)驗(yàn)中選用的光源為白熾燈和熒光燈。

將便攜式地物高光譜儀安裝在儀器自配的三角架上，將滴有10 至16 滴未經(jīng)稀疏的毒死蜱農(nóng)藥的新鮮菠菜葉經(jīng)風(fēng)干處理后，放置在事先鋪好的黑布上，菠菜葉旁邊放置相對(duì)反射率接近99%標(biāo)準(zhǔn)反射白板(由聚四氟乙烯制成)，光譜儀正對(duì)菠菜葉和白板拍攝。兩盞白熾燈分別旋轉(zhuǎn)在黑布的左右兩邊斜對(duì)樣品進(jìn)行照射，一盞熒光燈放置在高光譜儀的對(duì)面斜對(duì)樣品進(jìn)行照射。

二、蔬菜和農(nóng)藥光譜特性分析

2.1 實(shí)驗(yàn)數(shù)據(jù)選擇

在實(shí)驗(yàn)獲取的每幅圖像上，利用ENVI 4.7 軟件人工選取一塊合適大小的感興趣區(qū)域(ROI)，其選擇過(guò)程如圖4所示。左邊是反射率接近99%的標(biāo)準(zhǔn)漫反射白板，右邊是滴有毒死蜱的菠菜葉，其他部分是全黑的黑布。

實(shí)驗(yàn)中，我們用Matlab 2014a 版本軟件讀取高光譜圖像數(shù)據(jù)，發(fā)現(xiàn)農(nóng)藥點(diǎn)在高光譜圖像數(shù)據(jù)中的表現(xiàn)有3種形式，分別是全亮點(diǎn)、全暗點(diǎn)和部分亮部分暗點(diǎn)，如圖5所示，紅框顯示的是全亮農(nóng)藥點(diǎn)，黃框顯示的是全暗農(nóng)藥點(diǎn)，綠框顯示的是部分亮部分暗農(nóng)藥點(diǎn)。農(nóng)藥點(diǎn)之所以會(huì)顯示出這3 種形態(tài)，是因?yàn)槎舅莉甾r(nóng)藥是一種乳油狀的農(nóng)藥，風(fēng)干后，滴有農(nóng)藥的地方會(huì)在菠菜葉片上留下一個(gè)很小的油膜，而實(shí)驗(yàn)過(guò)程中使用了3盞燈，分別放置在蔬菜樣品的不同方位上，這樣位于樣品正上方的高光譜儀鏡頭接收到的油膜反射的三盞燈的光照強(qiáng)度也不一樣，根據(jù)高光譜儀鏡頭接收到的光的多少，就分別形成全亮、全暗和部分亮部分暗三種農(nóng)藥點(diǎn)。雖然有3種狀態(tài)類型的農(nóng)藥，但是總結(jié)起來(lái)就是兩類，亮農(nóng)藥點(diǎn)和暗農(nóng)藥點(diǎn)。

三、相對(duì)反射率的計(jì)算

為了提高樣品光譜的穩(wěn)定性，更好的分析不同表現(xiàn)形態(tài)農(nóng)藥點(diǎn)的光譜特性，我們挑選了亮、暗高光譜農(nóng)藥點(diǎn)圖像各25幅進(jìn)行數(shù)據(jù)分析，對(duì)每類農(nóng)藥點(diǎn)取平均光譜作為該類農(nóng)藥點(diǎn)的光譜;取50幅高光譜菠菜圖像的菠菜葉背景的平均光譜作為菠菜葉光譜。圖6顯示的是亮、暗農(nóng)藥點(diǎn)和菠菜葉的平均DN曲線。紅色線表示亮農(nóng)藥點(diǎn)的光譜DN曲線，綠色線表示暗農(nóng)藥點(diǎn)的光譜DN曲線，藍(lán)色線表示菠菜葉的光譜DN曲線。由圖可知樣品的原始光譜DN 值曲線在406.1nm、518.5nm和583.1nm處具有非常明顯的反射峰，需要說(shuō)明的是這3個(gè)反射峰并不是樣品的光譜特性，而是因?yàn)閷?shí)驗(yàn)中所使用的光源中有一盞是熒光燈，熒光燈在上述3個(gè)波長(zhǎng)處的輻射強(qiáng)度非常強(qiáng)，造成樣品的原始光譜DN值曲線在這3個(gè)波長(zhǎng)處有特別明顯的反射峰。除了這3個(gè)反射峰，樣品的原始光譜DN值曲線仍然有明顯的大小差異和波動(dòng)，說(shuō)明樣品的光譜受到入射光強(qiáng)度的影響和噪聲的影響很大。在此我們采用公式(2.1)對(duì)樣品的原始光譜進(jìn)行反射率計(jì)算，以減少入射光強(qiáng)度變化對(duì)樣本的影響。

式(2.1)中， R 表示樣品的相對(duì)反射率， S I 表示樣品的原始DN 值， DI表示系統(tǒng)的暗電流強(qiáng)度， WI表示標(biāo)準(zhǔn)反射白板的DN 值。

在數(shù)據(jù)處理中，可以把黑布的DN 值看成是系統(tǒng)的暗電流，這樣我們就得到了樣品的相對(duì)反射率曲線，如圖2.9 所示。紅線表示亮農(nóng)藥點(diǎn)的相對(duì)反射率曲線，綠線表示暗農(nóng)藥點(diǎn)的光譜相對(duì)反射率曲線，藍(lán)線表示菠菜葉的相對(duì)反射率曲線。

對(duì)樣品的反射率進(jìn)行較正后，光譜曲線的波動(dòng)明顯減少。400nm 之前和940nm 之后的光譜數(shù)據(jù)受噪聲干擾嚴(yán)重，在實(shí)際數(shù)據(jù)處理中，將這兩段數(shù)據(jù)剔除，這樣，可用數(shù)據(jù)為400-940nm 合計(jì)226 個(gè)波段。

圖6 菠菜葉和毒死蜱農(nóng)藥的原始DN 值曲線圖

剔除噪聲較大的波段后，光譜相對(duì)反射率曲線仍然受到明顯的噪聲干擾，表現(xiàn)為相對(duì)反射率曲線起伏比較明顯，曲線不夠平滑，仍需對(duì)光譜數(shù)據(jù)進(jìn)行預(yù)處理以降低噪聲。

四、Savitzky-Golay(SG)濾波處理

目前，常見的光譜圖像數(shù)據(jù)預(yù)處理技術(shù)包括標(biāo)準(zhǔn)化(normalization)、平滑濾波(smoothing)、微分技術(shù)(derivative)、多元散射校正(MSC)、標(biāo)準(zhǔn)正態(tài)變量變換(SNV)、小波變換(WT)等，但是除平滑濾波外，大部分預(yù)處理方法都會(huì)改變光譜的形狀或者光譜的位置，不利于進(jìn)行光譜分析。

為了盡可能多的保持光譜曲線的原始特征，更好的進(jìn)行光譜特性分析，我們采用對(duì)光譜原始信息保留較好的平滑濾波方法對(duì)光譜曲線進(jìn)行降噪處理。平滑濾波法中比較常見的有滑動(dòng)平均濾波法、Savitzky-Golay(SG)濾波法等方法。

圖7 菠菜葉和毒死蜱農(nóng)藥的相對(duì)反射率曲線圖

滑動(dòng)平均濾波法是一種最簡(jiǎn)單的濾波方式，其表達(dá)式為

式(2.2)中，m 為濾波帶寬，2m+1 為滑動(dòng)窗口大小。在實(shí)際操作中，應(yīng)慎重選擇參數(shù)m，m 過(guò)小，則濾波效果不明顯，m 過(guò)大，否則會(huì)造成信號(hào)失真。

SG 濾波法是一種基于多項(xiàng)式擬合的平滑方法，又稱多項(xiàng)式平滑，波長(zhǎng)k 處光譜數(shù)值經(jīng)平滑后的平均值為：

式(2.3)中， hi為平滑系數(shù)， H 為歸一化系數(shù)，每一測(cè)量的光譜數(shù)值乘以平滑系數(shù) 是為了盡量減少平滑對(duì)有用信息的影響。SG 濾波法在時(shí)域內(nèi)利用最小二乘原理實(shí)現(xiàn)最佳擬合，在濾除噪聲的同時(shí)，還保持信號(hào)的形狀、寬度不變，在實(shí)際操作中應(yīng)用廣泛。

SG 濾波法與滑動(dòng)平均濾波法的基本思想相似，只是SG 濾波法是通過(guò)多項(xiàng)式來(lái)對(duì)滑動(dòng)窗口內(nèi)的數(shù)據(jù)進(jìn)行多項(xiàng)式最小二乘擬合，而不是簡(jiǎn)單的平均，SG 濾波實(shí)質(zhì)上是一種加權(quán)平均法，更強(qiáng)調(diào)中心點(diǎn)的中心作用。在實(shí)際操作中，經(jīng)過(guò)多次比較，我們選擇9 點(diǎn)SG 濾波方法對(duì)光譜曲線進(jìn)行降噪處理，其處理結(jié)果如圖8所示。

圖8 經(jīng)SG平滑處理的菠菜葉和毒死蜱農(nóng)藥的光譜曲線圖

對(duì)比圖7和圖8可知，經(jīng)過(guò)SG濾波處理后，菠菜葉和毒死蜱農(nóng)藥的光譜曲線得到明顯光滑，并且保持了原始光譜曲線的形狀，這說(shuō)明SG濾波處理去掉了很大一部分噪聲，有利于下步進(jìn)行光譜特性分析。

五、農(nóng)藥點(diǎn)與葉片背景的光譜分析

由圖8可以看出，菠菜葉片的光譜特性是典型的植被光譜特性。在可見光波段內(nèi)，葉綠素、葉黃素等色素是影響菠菜葉片光譜響應(yīng)的主要因素，尤其是葉綠素。在中心波長(zhǎng)分別為450nm的藍(lán)光波段和650nm的紅光波段，葉綠素吸收大部分入射能量，而位于這兩個(gè)葉綠素吸收帶之間，由于吸收作用較小，在520nm的綠光波段附近形成了一個(gè)反射峰，因此，新鮮的菠菜葉看起來(lái)是綠色的。在可見光與近紅外波段之間，大約在680nm附近，反射率急劇上升，形成“紅邊”效應(yīng)，這也是植被曲線最明顯的特征。

由于亮農(nóng)藥點(diǎn)的光譜受光照條件影響太大，不能代表農(nóng)藥的光譜曲線，而暗農(nóng)藥點(diǎn)的光譜則受光照條件的影響相對(duì)較小，可以代表農(nóng)藥的光譜曲線。毒死蜱農(nóng)藥的有效成分為45%，滴在菠菜葉片上的毒死蜱農(nóng)藥光譜曲線表現(xiàn)為菠菜葉和農(nóng)藥的混合光譜，其光譜曲線形狀與菠菜葉的光譜曲線形狀類似，但是與菠菜葉片的光譜曲線相比，仍然有一些明顯的區(qū)別。整體而言，毒死蜱農(nóng)藥的相對(duì)反射率比菠菜葉的低，尤其是在520nm和700-900nm之間，這種表現(xiàn)更為明顯。而在660-680nm之間，兩者的相對(duì)反射率最為接近。在520nm處，菠菜葉有一個(gè)比較明顯的反射峰，而毒死蜱農(nóng)藥在此處的峰卻相對(duì)平緩。在884.3nm和920.2nm處，毒死蜱農(nóng)藥各有一個(gè)小的吸收峰，菠菜葉在相應(yīng)位置沒(méi)有吸收峰。在571.1nm和870nm處，毒死蜱農(nóng)藥各有一個(gè)小的反射峰，菠菜葉則沒(méi)有。

根據(jù)亮、暗農(nóng)藥點(diǎn)的光譜特性及它們與菠菜葉的光譜特性的區(qū)別，可以分別提取亮、暗農(nóng)藥點(diǎn)的特征波長(zhǎng)。如圖6所示，亮農(nóng)藥點(diǎn)的特征波長(zhǎng)為408.5nm、418.1nm、449.1nm、461.1nm、475.5nm、518.5nm、540nm、554.4nm、573.5nm、583.1nm、597.4nm和831.7nm，共12個(gè)特征波長(zhǎng)，暗農(nóng)藥點(diǎn)的特征波長(zhǎng)為408.5nm、468.3nm、520.9nm、571.1nm、607nm、645.2nm、726.5nm、740.9nm、810.2nm、870nm、884.3nm和920.2nm，共12個(gè)特征波長(zhǎng)。特征波長(zhǎng)的提取為下步進(jìn)行異常檢測(cè)提供了數(shù)據(jù)支撐。

審核編輯：湯梓紅