基于高光譜成像的水果損傷分析研究

摘 要：高光譜成像技術(shù)包含圖像信息和光譜信息。本文利用高光譜成像技術(shù)檢測(cè)蘋(píng)果摔傷，主要采用主成分分析、波段比算法和支持向量機(jī)分析所采集的高光譜圖像數(shù)據(jù)。實(shí)驗(yàn)結(jié)果表明，波段比算法和主成分分析法分類(lèi)識(shí)別正確率為93.3%，與支持向量機(jī)相比更適用于蘋(píng)果摔傷的實(shí)時(shí)快速檢測(cè)。

關(guān)鍵詞：水果損傷；高光譜；波段比算法；主成分分析；支持向量機(jī)

中圖分類(lèi)號(hào)：TS255.7 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1003-5168（2018）10-0028-05

Damage Analysis of Fruits Based on Hyper spectral Imaging

HAN Haoran1，2 LI Meng1，2 DU Dewei1，2 PAN Mingkang1，2 WANG Xinye1，2

（1.College of Physics and Electronic， Yunnan Normal University，Kunming Yunnan 650500；2.Provincial Key Laboratory for Opto-electronic Information Technology，Kunming Yunnan 650500）

Abstract： Hyper-spectral imaging technology includes image information and spectral information. This paper used hyperspectral imaging technology to detect apple fall. In the process of experiment， principal component analysis， band ratio algorithm and support vector machine were used to analyze hyperspectral image data collected. The experimental results showed that the accuracy of band ratio algorithm and principal component analysis was 93. 3%， which was more suitable for real time and fast detection of apple fall than support vector machine.

Keywords： fruit injury；hyper-spectral；band ratio algorithm；principal component analysis； support vector machine

我國(guó)是水果生產(chǎn)和消費(fèi)大國(guó)，但目前的水果檢測(cè)與分類(lèi)分級(jí)技術(shù)還比較落后，主要靠人工進(jìn)行分類(lèi)分級(jí)，這種方法存在勞動(dòng)時(shí)間長(zhǎng)、容易疲勞、主觀性強(qiáng)等缺點(diǎn)。因此，迫切需要一種客觀、準(zhǔn)確、無(wú)損的檢測(cè)分類(lèi)技術(shù)。而融合了光譜信息和圖像信息的高光譜成像技術(shù)，可以更全面地反映農(nóng)產(chǎn)品的信息，因此，被廣泛應(yīng)用。高光譜成像技術(shù)，配合合適的光譜重建算法，其光譜維和空間維信息既可以檢測(cè)水果的內(nèi)在物理性質(zhì)和化學(xué)成分等，又能全面反映出農(nóng)產(chǎn)品的外部特征、表面缺陷[1]及污染[2]等情況。從農(nóng)產(chǎn)品的高光譜圖像中提取出相應(yīng)的特征光譜波長(zhǎng)，就能對(duì)農(nóng)產(chǎn)品的品質(zhì)進(jìn)行準(zhǔn)確快速地檢測(cè)。因此，在檢測(cè)水果的品質(zhì)時(shí)，可以從光譜信息或是圖像信息的角度分析其有效特征信息。由于高光譜圖像具有海量、高維的信息，因此，在光譜重建和后期光譜分析時(shí)需要對(duì)其進(jìn)行降維。主成分分析（principal component analysis，PCA）、波段比算法（Band Ratio Algorithm）和支持向量機(jī)法（Support Vector Machine Method）等是常用的降維方法，在高光譜圖像的處理中得到廣泛應(yīng)用。主成分分析法利用降維技術(shù)用少數(shù)幾個(gè)綜合變量來(lái)代替原始多個(gè)變量，這些綜合變量集中了原始變量的大部分信息；其通過(guò)計(jì)算綜合主成分函數(shù)得分，對(duì)客觀經(jīng)濟(jì)現(xiàn)象進(jìn)行科學(xué)評(píng)價(jià)。波段比算法可以使光照不均勻影響達(dá)到最小化，增強(qiáng)波段之間的差異性，提供一些單波段下無(wú)法得到的有效信息。支持向量機(jī)法選取感興趣區(qū)，抽取少量波段，運(yùn)算速度較快。本研究應(yīng)用高光譜成像技術(shù)，研究了無(wú)損檢測(cè)肉眼看不到的蘋(píng)果機(jī)械損傷的可能性，并用主成分分析、波段比運(yùn)算、監(jiān)督分類(lèi)的支持向量機(jī)三種方法對(duì)高光譜圖像進(jìn)行分析研究，確定了合適的檢測(cè)方法。

1 降維算法

1.1 主成分分析法

1.1.1 主成分分析的原理。主成分分析法是利用降維（線性變換）的思想，把多種指標(biāo)化簡(jiǎn)為少數(shù)的幾個(gè)綜合指標(biāo)（主成分），化簡(jiǎn)后的前幾個(gè)主成分就能反映原始數(shù)據(jù)的絕大部分信息，而且所含信息之間沒(méi)有重復(fù)。這種方法將復(fù)雜因素歸結(jié)為少數(shù)幾個(gè)主成分，使問(wèn)題得以簡(jiǎn)化，且數(shù)據(jù)信息更加科學(xué)有效[3]。主成分分析的特點(diǎn)是利用線性擬合的思路把分布在多個(gè)維度的高維數(shù)據(jù)投射到幾個(gè)軸上。如果每個(gè)樣本只有兩個(gè)數(shù)據(jù)變量，這種擬合就是：

[a1x1+a2x2=PC] （1）

其中，[x1]和[x2]分別是樣本的兩個(gè)變量，而[a1]和[a2]則被稱(chēng)為載荷，計(jì)算出的[PC]值就被稱(chēng)為主成分。

實(shí)際上，當(dāng)一個(gè)樣本只有兩個(gè)變量時(shí)，主成分分析本質(zhì)上就是一個(gè)線性回歸。式（1）本質(zhì)上就是一條直線。主成分坐標(biāo)旋轉(zhuǎn)圖[3]見(jiàn)圖1。

如果一個(gè)樣本有m個(gè)變量，那主成分就變?yōu)椋?/p>

[a1x1+a2x2+…+amxm=PC1] （2）

其中，[PC1]稱(chēng)為第一主成分，而且還可以獲得一系列與這個(gè)[PC]直線正交的其他直線，如：

[b1x1+b2x2+…+bmxm=PC2] （3）

[PC2]被稱(chēng)為第二主成分，以此類(lèi)推。若令[A=a1，a2，…，am]，[X=x1，x2，…，xm]，則

[ATX=PC] （4）

此時(shí)，向量[A]稱(chēng)為主成分的載荷（Loading），計(jì)算出的主成分的值[PC]稱(chēng)為得分（Score）。

1.1.2 主成分分析處理高光譜數(shù)據(jù)的計(jì)算步驟。設(shè)高光譜數(shù)據(jù)含有[P]個(gè)波段，每幅圖像包含[K=M×N]個(gè)像素，其中[M]表示行數(shù)，[N]表示列數(shù)，用矩陣[X=x1，x2，…，xp]表示高光譜圖像數(shù)據(jù)，其中[xi]（1≤i≤p）為各波段圖像，為[K×1]維的列向量，則X為[K×P]維矩陣[4]，主成分分析的計(jì)算步驟具體如下。

①對(duì)原始數(shù)據(jù)分波段進(jìn)行標(biāo)準(zhǔn)化處理，標(biāo)準(zhǔn)化后得到矩陣Y，Y矩陣為：

[yij=xij-xjVarxji=1，2，…，K，j=1，2，…，P] （5）

其中，[xj]和[Varxj]分別為：

[xj=1Kk=1Kxij] （6）

[Varxi=1K-1i=1Kxij-xj2j=1，2，…，K] （7）

②得到相關(guān)系數(shù)矩陣R，并計(jì)算 的特征值和特征向量。經(jīng)標(biāo)準(zhǔn)化處理后的數(shù)據(jù)構(gòu)成相關(guān)系數(shù)矩陣R（P×P方陣）：

[rij=1K-1t=1Kyityjti=1，2，3，…，P，j=1，2，3，…，P] （8）

相關(guān)系數(shù)矩陣R對(duì)應(yīng)的特征值[λ1≥λ2≥…≥λP]和特征向量[ai=ai1，ai2，ai3，…，aik]，其中[i=1，2，3，…，P]。

③以AT對(duì)原圖像進(jìn)行計(jì)算，得到主成分PC，計(jì)算公式見(jiàn)式（4）。

④給出一個(gè)數(shù)值T（如0.85），根據(jù)公式（7）計(jì)算出累計(jì)貢獻(xiàn)[Tm≥T]的主成分個(gè)數(shù)m。

[Tm=i=1mλii=1Pλi] （9）

⑤當(dāng)前，m個(gè)主成分大于等于[T]時(shí)，選取前m個(gè)主成分，代替原來(lái)的高光譜圖像X。

1.2 波段比算法

1.2.1 波段比算法原理。波段比算法又稱(chēng)為比值增強(qiáng)（Ratio Enhancement），是計(jì)算同一幅圖像不同波段間相應(yīng)像元的像元值之比的圖像增強(qiáng)方法。波段比就是用一個(gè)波段除以另一個(gè)波段生成一幅能提供相對(duì)波段強(qiáng)度的圖像。比值函數(shù)的數(shù)學(xué)表達(dá)式如式（10）所示[5]：

[BVi，j，r=BVi，j，m/BVi，j，n] （10）

式（10）中，[BVi，j，m]和[BVi，j，n]分別是第m和n波段相同位置像素[i，j]的亮度值。[BVi，j，r]是該位置下像素[i，j]的m和n波段的比率值。分母[BVi，j，n]若為0，輸出比值[BVi，j，r]將賦值為0。使用歸一化函數(shù)進(jìn)行進(jìn)一步處理，如式（11）所示[5]：

[BVi，j，o= 0 BVi，j，r=0IntBVi，j，r×127+1 BVi，j，r∈1255，1IntBVi，j，r2+128 BVi，j，r∈1，255] （11）

式（5）中，[BVi，j，o]表示像素[i，j]的輸出亮度值，[Int]表示取整運(yùn)算。

1.2.2 特征波長(zhǎng)選取。由于高光譜圖像數(shù)據(jù)量非常大，相鄰波段之間相關(guān)性很強(qiáng)，因此特征波段的選擇非常關(guān)鍵。波段之間的相關(guān)性越小，波段比值圖像的信息量就越大，所以必須尋找相關(guān)性較小的特征波長(zhǎng)。Sheffield指數(shù)（Sheffield Index，SI）[6-8]可以很好地用來(lái)確定和評(píng)價(jià)最佳波段，如式（12）所示：

[SI=Covp×p] （12）

其中，p為所選擇的波段數(shù)目，當(dāng)選擇p>2時(shí)，根據(jù)[SI]確定的波段數(shù)也大于2，仍需確定相關(guān)性最小的兩個(gè)波段，因此只需考慮P=2的情況。分別計(jì)算各種波段組合的[SI]值，[SI]值越大，相關(guān)性越小。對(duì)[SI]值進(jìn)行排序，就可以選到特征波長(zhǎng)[6]。

1.3 監(jiān)督分類(lèi)的支持向量機(jī)

1.3.1 支持向量機(jī)方法原理。支持向量機(jī)（Support Vector Machine，SVM）是Corinna Cortes和Vapnik等于1995年首先提出的，其在解決小樣本、非線性及高維模式識(shí)別中表現(xiàn)出特有的優(yōu)勢(shì)，并能推廣應(yīng)用到函數(shù)擬合等機(jī)器學(xué)習(xí)問(wèn)題中[9]。在機(jī)器學(xué)習(xí)中，支持向量機(jī)是與相關(guān)學(xué)習(xí)算法有關(guān)的監(jiān)督學(xué)習(xí)模型，可以分析數(shù)據(jù)、識(shí)別模式以及用于分類(lèi)和回歸分析。

SVM的關(guān)鍵在于核函數(shù)。低維空間向量集很難劃分，解決的方法是將其映射到高維空間。但是，帶來(lái)的問(wèn)題是計(jì)算變得更復(fù)雜。解決這個(gè)問(wèn)題的方法是引入核函數(shù)。只要選用適當(dāng)?shù)暮撕瘮?shù)，就可以得到高維空間分類(lèi)函數(shù)。在支持向量機(jī)理論中，采用不同的核函數(shù)將產(chǎn)生不同的支持向量機(jī)算法。常見(jiàn)的核函數(shù)有以下四種[10]。

①線性核函數(shù)（Linear Kernel）。計(jì)算公式為：

[kx，y=xTy+c] （13）

②多項(xiàng)式核函數(shù)（Polynomial Kernel）。計(jì)算公式為：

[kx，y=acTy+cd] （14）

③徑向基函數(shù)（Radial Basis Function）。計(jì)算公式為：

[kx，y=exp-γx-y2] （15）

④二層神經(jīng)網(wǎng)絡(luò)核函數(shù)（Sigmoid Kernel）。計(jì)算公式為：

[kx，y=tanhaxT+c] （16）

本研究采用多項(xiàng)式核函數(shù)和徑向基函數(shù)。

2 實(shí)驗(yàn)及數(shù)據(jù)處理

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

市場(chǎng)購(gòu)買(mǎi)不同品種的蘋(píng)果15個(gè)，目視檢測(cè)均完好無(wú)損。將蘋(píng)果置于室溫下1d后開(kāi)始實(shí)驗(yàn)。將蘋(píng)果在高1m左右摔落于硬質(zhì)地板，表面肉眼觀察不到損傷。摔后損傷區(qū)手感稍微發(fā)軟。

2.2 圖像采集

圖像采集采用北京卓立漢光儀器有限公司的Gaia Sorter高光譜分選儀，獲得一個(gè)包含圖像信息和光譜信息的三維數(shù)據(jù)立方體。在采集高光譜圖像之前，需要對(duì)系統(tǒng)各參數(shù)進(jìn)行調(diào)整、設(shè)置，以保證采集到清晰的圖像。CCD曝光時(shí)間設(shè)置為15ms，掃描寬度設(shè)為150mm，輸送速度設(shè)為0.6cm/s。為校正CCD暗電流等對(duì)圖像的影響，對(duì)圖像進(jìn)行黑白校準(zhǔn)處理。設(shè)W為掃描白板采集到的全白的圖像，B為蓋上鏡頭采集到的全黑圖像，I是原始的高光譜圖像，R是黑白校準(zhǔn)后的高光譜圖像。采用高光譜系統(tǒng)配套的SPECview軟件進(jìn)行處理。光譜黑白校正的計(jì)算公式為[6-12]：

[R=I-BW-B] （17）

2.3 數(shù)據(jù)分析

2.3.1 主成分分析方法。主成分分析法處理高光譜圖像就是為了把實(shí)驗(yàn)樣品中包含的海量信息用少數(shù)幾個(gè)變量來(lái)表示，對(duì)高光譜圖像信息進(jìn)行降維。采用ENVI5.0軟件選取760～1 000nm波長(zhǎng)進(jìn)行主成分分析，取前6個(gè)主成分。各主成分經(jīng)ENVI5.0軟件處理后的圖像如圖2所示。

從圖2可以清楚地看到，第1至第3主成分圖像，完全觀察不到損傷的痕跡，第4至第6主成分圖像，可以觀察到損傷部分，以第5主成分圖像顯示損傷效果比較好，圖像最清晰。

2.3.2 波段比算法。波段比運(yùn)算就是用兩個(gè)波段相除生成新的相對(duì)波段強(qiáng)度圖像的過(guò)程。其可以降低蘋(píng)果表面不平整帶來(lái)光線反射不均勻的影響，還可以增強(qiáng)波段之間的波譜差異，提供一些單波段無(wú)法得到的信息。從幾組正常和摔傷的蘋(píng)果中各選3組光譜數(shù)據(jù)，正常完好的蘋(píng)果某一部位的光譜反射曲線如圖3所示，同一部位損傷后的光譜反射圖像如圖4所示。

對(duì)比圖3和圖4可得出，損傷部分和正常部分的光譜反射在紅外部分差別明顯。正常蘋(píng)果的在紅外部分反射率斜率較大；兩個(gè)損傷部分的光譜反射在560nm附近有個(gè)峰值，與正常部分相近波段的光譜反射率有明顯區(qū)別。損傷部分在紅外區(qū)反射率起伏不大，較正常部分變化平緩?？傊?，蘋(píng)果正常部分的光譜反射率與摔傷部分的光譜反射率存在明顯差異，為通過(guò)特征波段進(jìn)行分類(lèi)識(shí)別提供了可能性。

通過(guò)對(duì)比兩個(gè)光譜反射率，應(yīng)用波段比率法原理，對(duì)SI值進(jìn)行排序。損傷區(qū)在近紅外波段選取801nm和1 000nm做波段比運(yùn)算，處理后的圖像如圖5（c）所示。

圖5（a）是損傷水果在555nm附近的灰度圖像，圖5（b）損傷水果在700nm附近的灰度圖像，圖5（c）是經(jīng)過(guò)波段比法處理過(guò)的圖像。圖5（a）和圖5（b）比較亮的部分是鏡面反射，看不到損傷。根據(jù)人眼感知機(jī)理，人眼對(duì)555nm波長(zhǎng)最敏感，但在圖5（a）中卻看不到損傷情況；蘋(píng)果顏色多為紅色或黃色，但在圖5（b）700nm附近也察覺(jué)不到損傷；圖5（c）為經(jīng)過(guò)波段比法處理過(guò)的圖像，清楚地顯示出了中間近似圓形的損傷區(qū)，周?chē)鸀橥旰脜^(qū)域。通過(guò)圖5綜合分析，波段比運(yùn)算能在近紅外波段快速準(zhǔn)確檢測(cè)出蘋(píng)果原始圖像不易顯示的機(jī)械損傷。

2.3.3 監(jiān)督分類(lèi)的支持向量機(jī)法。支持向量機(jī)方法是采用ENVI5.0軟件，在損傷區(qū)與完好區(qū)選取感興趣區(qū)，采用多項(xiàng)式核函數(shù)和徑向基函數(shù)，分別進(jìn)行數(shù)據(jù)處理。因?yàn)槿庋劭床坏綋p傷，盡管損傷區(qū)在560nm有個(gè)峰值，不選為特征波段。從圖4與圖5可知，在近紅外區(qū)，正常蘋(píng)果與損傷區(qū)的光譜反射率差別較大，所以在近紅外區(qū)選取760、761、810、811、860、861、910、911、960nm和961nm作為特征波段。利用這些波段對(duì)圖像進(jìn)行處理，圖5（a）的蘋(píng)果處理后的圖像如圖6所示。圖6（a）是多項(xiàng)式核函數(shù)的圖像處理結(jié)果，圖6（b）是徑向基函數(shù)的圖像處理結(jié)果。圖6（a）分類(lèi)效果較主成分分析和波段比算法差，圖像不規(guī)則，因?yàn)樘O(píng)果是摔在硬的地板磚上的，損傷區(qū)應(yīng)該是接近圓形的規(guī)則圖像。圖6（b）沒(méi)有準(zhǔn)確反映出損傷部位，反映的只是光照較強(qiáng)的部分。

3 數(shù)據(jù)分析結(jié)果及對(duì)比評(píng)價(jià)

對(duì)15個(gè)樣本數(shù)據(jù)分別采用以上三種算法進(jìn)行分析，得出的分類(lèi)、篩選結(jié)果如表1所示。

主成分分析和波段比運(yùn)算有14個(gè)樣本能很好地顯示出損傷區(qū)，正確率為93.3%。對(duì)于同樣一組高光譜數(shù)據(jù)，監(jiān)督分類(lèi)的支持向量機(jī)（多項(xiàng)式核函數(shù)）方法只有11個(gè)在摔傷區(qū)有不規(guī)則投影，圖像質(zhì)量很差。

主成分分析選用的是波長(zhǎng)760nm到1 000nm進(jìn)行主成分分析，數(shù)據(jù)量大，檢測(cè)速度較慢。波段比算法選取801nm和1 000nm做波段比運(yùn)算，數(shù)據(jù)量小，計(jì)算速度快，損傷區(qū)檢測(cè)也比較準(zhǔn)確。監(jiān)督分類(lèi)的支持向量機(jī)效果較差。綜合比較，波段比運(yùn)算是最有前途的檢測(cè)水果機(jī)械損傷的無(wú)損檢測(cè)方法。

4 結(jié)論

①本研究采用主成分分析法、波段比算法和監(jiān)督分類(lèi)支持向量機(jī)三種方法分析高光譜圖像，對(duì)肉眼看不到的蘋(píng)果機(jī)械損傷進(jìn)行識(shí)別，前兩種方法檢測(cè)正確率達(dá)到93.3%，支持向量機(jī)效果較差。

②主成分分析法與選取的波段、選用協(xié)方差矩陣或相關(guān)系數(shù)矩陣均有關(guān)。選用不同的波段及不同的矩陣，最后采用的主成分不一樣。

③實(shí)驗(yàn)結(jié)果表明，克服相鄰波段之間的相關(guān)性強(qiáng)的方法是SI。SI可以快速確定出特征波長(zhǎng)，從而快速實(shí)現(xiàn)高光譜數(shù)據(jù)降維。

④波段比算法可以有效克服光照不均勻?qū)悠樊a(chǎn)生的影響，提取光譜間的差異，迅速?gòu)母吖庾V數(shù)據(jù)中找出準(zhǔn)確信息，準(zhǔn)確地確定肉眼看不見(jiàn)的損傷區(qū)域的輪廓。通過(guò)選取合適的特征波段，就可以達(dá)到快速檢測(cè)出水果機(jī)械損傷的目的。

參考文獻(xiàn)：

[1]趙娟，彭彥坤，趙松瑋，等.基于高光譜技術(shù)檢測(cè)蘋(píng)果外觀缺陷[J].食品安全質(zhì)量檢測(cè)學(xué)報(bào)，2012（6）：681-684.

[2]索少增，劉翠玲，吳靜珠，等.高光譜圖像技術(shù)檢測(cè)梨表面農(nóng)藥殘留試驗(yàn)研究[J].北京工商大學(xué)學(xué)報(bào)（自然科學(xué)版），2011（6）：73-77.

[3]彭謹(jǐn).主成分分析入門(mén)[EB/OL].（2016-10-16）[2018-03-01].https：//www.cnblogs.com/SCUJIN/p/5965946.html.

[4]張明月，賀金蘭，田尉霞.主成分分析在高光譜圖像降維中的應(yīng)用[J].內(nèi)蒙古科技與經(jīng)濟(jì)，2017（15）：88-90.

[5]蔡健榮，王建黑，陳全勝，等.波段比算法結(jié)合高光譜圖像技術(shù)檢測(cè)柑橘果銹[J].農(nóng)業(yè)工程學(xué)報(bào)，2009（1）：127-131.

[6]Inoue Y， Pe?uelas J， Miyata A， et al. Normalized difference spectral indices for estimating photosynthetic efficiency and capacity at a canopy scale derived from hyperspectral and CO2 flux measurements in rice[J]. Remote Sensing of Environment， 2008（1）：156-172.

[7] Steve De Backer， Aleksandra Pi?urica， Bruno Huysmans， et al. Denoising of multicomponent images using wavelet least-squares estimators[J]. Image and Vision Computing， 2008（7）：1038-1051.

[8] Thomas Koenig，Lester Melie-Garcia， Maria Stein， et al. Establishing correlations of scalp field maps with other experimental variables using covariance analysis and resampling methods[J]. Clinical Neurophysiology， 2008（6）：1262-1270.

[9]百度百科.支持向量機(jī)[EB/OL].（2015-01-13）[2018-03-01].https：//baike.baidu.com/item/%E6%94%AF%E6%8C%81%E5%90%91%E9%87%8F%E6%9C%BA/9683835.

[10]常用的核函數(shù)[EB/OL].（2016-07-03）[2018-03-01].https：//blog.csdn.net/qq_27231343/article/details/51817866.

[11]趙杰文，劉劍華，陳全勝，等.利用高光譜圖像技術(shù)檢測(cè)水果輕微損傷[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2008（1）：106-109.

[12]陳全勝，趙杰文，蔡健榮，等.利用高光譜圖像技術(shù)評(píng)判茶葉的質(zhì)量等級(jí)[J].光學(xué)學(xué)報(bào)，2008（4）：669-674.