引言

當前，蔬菜中含有農藥殘留的可能性越來越大，而相對成熟的農藥殘留檢測技術如氣相色譜法等存在費時費力，破壞性大等明顯缺點，因此，急需一種快速無損的蔬菜農藥殘留檢測技術。高光譜成像技術能夠區分地物光譜的微小差別，已有研究人員將其應用于食品安全監測領域，并且取得了一定的研究成果，在蔬菜農藥殘留檢測方面具有較大的潛力。

異常檢測是高光譜遙感最常見的應用之一，已經在軍事目標識別、偽裝目標區分等方面取得了較好的成果。鑒于異常檢測在高光譜遙感領域目標識別中的有效性，本文將其應用于高光譜圖像，試圖探究其檢測蔬菜葉片上的農藥點的可行性。

在本文中，我們首先搭建了一個高光譜成像平臺，獲取了包含農藥點的菠菜葉片的高光譜圖像數據，對菠菜葉片的光譜特性和不同形態表現的農藥點的光譜特性進行了分析，確定了農藥點光譜的特征波長。

一、蔬菜農藥點檢測的高光譜數據采集

1.1 蔬菜農業點檢測的高光譜成像系統

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

推掃成像方式中硬件設備主要包括光源、推掃的高光譜成像系統、鏡頭、傳送裝置和計算機等。光源一般采用鹵素燈或者LED燈，能夠產生可見光到紅外光的連續光。推掃的高光譜成像系統用來記錄物體一個掃描線上的多個波段的入射光輻射能量，通常由前端的光譜儀和后端的光電傳感器兩部分組成。鏡頭用來調整光的范圍。傳送裝置用來輸送樣品，并實現在空間維度的掃描。計算機則用來存儲三維高光譜“立方體”數據，并實現數據的實時處理。圖1是該類設備的典型案例的示意圖。

圖 1 推掃型高光譜成像系統結構原理圖

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

圖 2 基于液晶可調濾光片型高光譜成像系統示意圖

1.2 高光譜圖像獲取

獲取含有農藥點的蔬菜高光譜圖像數據，首先需要確定蔬菜的種類和農藥點的類型。

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

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

為了獲取含有農藥的蔬菜高光譜圖像，我們搭建了一個簡易的獲取高光譜圖像數據的實驗平臺，主要的實驗設備及材料如下：

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

樣品：就近市場購買的新鮮菠菜，有效成分為45%的毒死蜱農藥。

光源：兩盞功率為90W 的白熾燈，一盞功率為90W 的熒光燈。由于我們所搭建的實驗平臺比較簡易，而專業的光源獲取又相對困難，熒光燈在可見光波段的輻射能量較強，白熾燈在近紅外波段的輻射能量較強，這兩種類型的光源比較容易獲取，因此，我們在實驗中選用的光源為白熾燈和熒光燈。

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

二、蔬菜和農藥光譜特性分析

2.1 實驗數據選擇

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

實驗中，我們用Matlab 2014a 版本軟件讀取高光譜圖像數據，發現農藥點在高光譜圖像數據中的表現有3種形式，分別是全亮點、全暗點和部分亮部分暗點，如圖5所示，紅框顯示的是全亮農藥點，黃框顯示的是全暗農藥點，綠框顯示的是部分亮部分暗農藥點。農藥點之所以會顯示出這3 種形態，是因為毒死蜱農藥是一種乳油狀的農藥，風干后，滴有農藥的地方會在菠菜葉片上留下一個很小的油膜，而實驗過程中使用了3盞燈，分別放置在蔬菜樣品的不同方位上，這樣位于樣品正上方的高光譜儀鏡頭接收到的油膜反射的三盞燈的光照強度也不一樣，根據高光譜儀鏡頭接收到的光的多少，就分別形成全亮、全暗和部分亮部分暗三種農藥點。雖然有3種狀態類型的農藥，但是總結起來就是兩類，亮農藥點和暗農藥點。

三、相對反射率的計算

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

式(2.1)中， R 表示樣品的相對反射率， S I 表示樣品的原始DN 值， DI表示系統的暗電流強度， WI表示標準反射白板的DN 值。

在數據處理中，可以把黑布的DN 值看成是系統的暗電流，這樣我們就得到了樣品的相對反射率曲線，如圖2.9 所示。紅線表示亮農藥點的相對反射率曲線，綠線表示暗農藥點的光譜相對反射率曲線，藍線表示菠菜葉的相對反射率曲線。

對樣品的反射率進行較正后，光譜曲線的波動明顯減少。400nm 之前和940nm 之后的光譜數據受噪聲干擾嚴重，在實際數據處理中，將這兩段數據剔除，這樣，可用數據為400-940nm 合計226 個波段。

圖6 菠菜葉和毒死蜱農藥的原始DN 值曲線圖

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

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

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

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

圖7 菠菜葉和毒死蜱農藥的相對反射率曲線圖

滑動平均濾波法是一種最簡單的濾波方式，其表達式為

式(2.2)中，m 為濾波帶寬，2m+1 為滑動窗口大小。在實際操作中，應慎重選擇參數m，m 過小，則濾波效果不明顯，m 過大，否則會造成信號失真。

SG 濾波法是一種基于多項式擬合的平滑方法，又稱多項式平滑，波長k 處光譜數值經平滑后的平均值為：

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

SG 濾波法與滑動平均濾波法的基本思想相似，只是SG 濾波法是通過多項式來對滑動窗口內的數據進行多項式最小二乘擬合，而不是簡單的平均，SG 濾波實質上是一種加權平均法，更強調中心點的中心作用。在實際操作中，經過多次比較，我們選擇9 點SG 濾波方法對光譜曲線進行降噪處理，其處理結果如圖8所示。

圖8 經SG平滑處理的菠菜葉和毒死蜱農藥的光譜曲線圖

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

五、農藥點與葉片背景的光譜分析

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

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

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

審核編輯：湯梓紅