吳 威，吳明清，陳桂云，于鎮(zhèn)偉，陳坤杰*

雞胴體表面稀釋污染物的連續(xù)投影-多元線性回歸-受試者特性分析檢測

吳 威，吳明清，陳桂云，于鎮(zhèn)偉，陳坤杰*

（南京農(nóng)業(yè)大學(xué)工學(xué)院，江蘇 南京 210031）

建立基于連續(xù)投影算法（successive projections algorithm，SPA）-多元線性回歸（multiple linear regression，MLR）-受試者特性（receiver operating characteristic，ROC）分析雞胴體表面稀釋污染物識(shí)別方法。首先采用高光譜成像系統(tǒng)獲取了20 個(gè)涂有稀釋污染物的雞胴體圖像，再利用SPA從1 232 個(gè)波長中提取出10 個(gè)特征波長，然后通過MLR方法構(gòu)建識(shí)別函數(shù)和特征波長光譜值之間的回歸模型，最后通過ROC分析方法，確定出能夠獲得高真陽性率（true positive rate，TPR）和低假陽性率（false positive rate，F(xiàn)PR）的最佳分類閾值，并據(jù)此對(duì)雞胴體表面稀釋后的污染物進(jìn)行檢測。結(jié)果顯示，利用SPA-MLR-ROC分類器檢測20 個(gè)污染雞胴體樣本，檢出污染區(qū)域的TPR達(dá)到98.08%，F(xiàn)PR僅為0.39%，高于波段比算法以及雙波段算法檢測相同樣本的準(zhǔn)確率，可見，SPA-MLRROC分類器方法檢測雞胴體表面稀釋污染物獲得了較好的效果，但由于樣本數(shù)量有限，還需要具有較大樣本量的研究來進(jìn)一步驗(yàn)證此方法檢測結(jié)果的穩(wěn)定性。

連續(xù)投影-多元線性回歸-受試者特性分析；分類器；污染物；高光譜圖像；雞胴體

在我國，雞肉是僅次于豬肉的第二大肉類消費(fèi)品。近20年來，雞肉產(chǎn)量以每年平均5%～6%的速度持續(xù)增長[1]，在這些規(guī)?；a(chǎn)加工過程中，雞胴體極易受到血液、膽汁、以及十二指腸、盲腸中的糞便等污染[2-6]。研究發(fā)現(xiàn)從商業(yè)屠宰廠提取的肉雞樣品中，95%～99%的雞胃中彎曲桿菌的檢測結(jié)果為陽性[7]，肉雞胃中的大腸桿菌含量為4.4（lg（CFU/g））[8]，這些含有致病菌的消化道內(nèi)容物在肉雞屠宰過程中可能會(huì)殘留在雞胴體表面。研究表明即使少量的盲腸糞便（5 mg）也會(huì)引起雞胴體表面彎曲桿菌數(shù)量上的顯著增加[9]。因此將受到污染的雞胴體檢測出來，并進(jìn)行適當(dāng)?shù)奶幚硎鞘直匾摹?/p>

人工及生化檢測無法實(shí)現(xiàn)污染物的在線自動(dòng)檢測。高光譜圖像同時(shí)具有光譜和圖像信息，并且光譜分辨率很高，可以提供識(shí)別物質(zhì)特征的豐富信息，因此被廣泛地應(yīng)用于雞肉檢測領(lǐng)域，如利用高光譜圖像技術(shù)可以檢測雞胸肉菌落總數(shù)[10]、預(yù)測新鮮雞肉彈性[11]、定量檢測雞肉中揮發(fā)性鹽基氮含量[12]、快速檢測雞肉表面脯氨酸含量[13]以及檢測雞胚胎性別[14]等。

國內(nèi)外的研究表明高光譜圖像技術(shù)可以成功地檢測出雞胴體表面的污染物，這些研究通常選取2～3 個(gè)特征波長，利用特征波長之間的波段比圖像進(jìn)行圖像處理[2-3,15-20]，從而實(shí)現(xiàn)污染物的識(shí)別，如文獻(xiàn)[2,15-19]均采用565/517 nm的波長比圖像進(jìn)行污染物的識(shí)別，文獻(xiàn)[3,20]在565/517 nm的波長比圖像處理的基礎(chǔ)上增加了802 nm波長處的圖像去除假陽性點(diǎn)。

既有研究的對(duì)象大多是從雞腸道中直接提取的糞便污染物[2-3,15-20]，但在實(shí)際的肉雞屠宰加工過程中，雞胴體表面的水分會(huì)將殘留在胴體表面的污染物稀釋，可見度很差，很難被人眼及普通的機(jī)器視覺識(shí)別。為了識(shí)別稀釋后的低可見污染物，陳坤杰等[21]提出一種雙波段的雞胴體表面稀釋污染物檢測方法，對(duì)雞胸腹局部區(qū)域殘留稀釋污染物的檢測真陽性率（true positive rate，TPR）可以達(dá)到81.6%。但這種方法檢測范圍過小，真陽性率也還有待提高。

本實(shí)驗(yàn)試圖構(gòu)建一種光譜學(xué)方法，通過光譜建模而不是圖像處理的方法精確提取低可見檢測目標(biāo)的特征，以實(shí)現(xiàn)雞胴體表面稀釋污染物的精確識(shí)別。首先通過連續(xù)投影算法（successive projections algorithm，SPA）選擇相關(guān)性最小的特征波長，然后利用多元線性回歸（multiple linear regression，MLR）構(gòu)建特征波長圖像中的反射率與識(shí)別函數(shù)之間的對(duì)應(yīng)關(guān)系，再通過受試者特性（receiver operating characteristic，ROC）分析方法選擇分類閾值，最后將識(shí)別效果與既有研究中普遍采用的波段比算法[2-3,15-20]，以及文獻(xiàn)[21]提出的雙波段檢測方法進(jìn)行對(duì)比。

1 材料與方法

1.1 材料

2015年12月從雨潤肉雞加工廠購買20 只45 d的白羽雞雞胴體，雞胴體為隨機(jī)挑選。在實(shí)驗(yàn)室人工將內(nèi)臟取出，并將腸道分離出來。每個(gè)雞胴體表面的污染物都是取自該胴體本身。

使用移液管從每個(gè)雞胴體分別抽取十二指腸和盲腸糞便各1 mL，從腹腔中抽取血液1 mL，從膽囊中抽取膽汁1 mL，放入不同的玻璃培養(yǎng)皿中。再分別抽取與污染物等體積的清水與污染物混合，分別制成稀釋的污染樣本，隨機(jī)涂布在雞胴體表面。

每個(gè)雞胴體表面涂布的污染區(qū)域數(shù)量最少6 個(gè)，最多16 個(gè)，20 個(gè)雞胴體表面共計(jì)有208 個(gè)污染區(qū)域。由于稀釋后的污染物流動(dòng)性強(qiáng)，無法使得所有污染區(qū)域面積保持一致，污染區(qū)域的面積在232 個(gè)像素和4 200 個(gè)像素之間。4 種污染區(qū)域（十二指腸糞便、盲腸糞便、膽汁和血液）的平均面積分別為1 501、1 354、1 630、1 421 個(gè)像素。

1.2 方法

1.2.1 高光譜圖像系統(tǒng)

實(shí)驗(yàn)采用的高光譜成像系統(tǒng)由圖像采集單元、照明單元和電動(dòng)平移臺(tái)構(gòu)成，如圖1所示。圖像采集單元包括一臺(tái)Raptor EM285CL高感度EMCCD相機(jī)、光譜范圍為358～1 021 nm的Imspector V10E線掃描成像光譜儀。照明單元為配有Ocean Optics HL-2000線性光導(dǎo)管組以及光纖和聚光鏡的鹵素光源，位于電動(dòng)平移臺(tái)上方220 mm處，垂直于電動(dòng)平移臺(tái)的軌道。電動(dòng)平移臺(tái)由一臺(tái)計(jì)算機(jī)控制，計(jì)算機(jī)中安裝高光譜系統(tǒng)自帶軟件Spectral Image（臺(tái)灣五鈴光學(xué)股份有限公司），用于控制曝光時(shí)間、電動(dòng)平移臺(tái)的速度、光源亮度等。

將雞胴體胸腹部向上，放在一個(gè)蒙著黑布的托盤上，再將托盤放在高光譜圖像系統(tǒng)電動(dòng)平移臺(tái)上，進(jìn)行高光譜圖像采集。設(shè)定電動(dòng)平移臺(tái)速率為1.3 mm/s，相機(jī)曝光時(shí)間為380 ms，進(jìn)行逐行掃描，采集得到20 個(gè)雞胴體圖像樣本，以raw格式存儲(chǔ)于計(jì)算機(jī)。每個(gè)高光譜圖像都包含空間信息（x，y）和光譜信息（z）。x方向像素為1 632；由于雞胴體大小不同，掃描時(shí)間不同，因此y方向的像素在600～1 000之間；z方向?yàn)?58～1 021 nm，共有1 232 個(gè)波長。圖像和數(shù)據(jù)的分析處理運(yùn)用了ENVI4.8可視化圖像處理平臺(tái)（美國Exelis VIS公司）和MATLAB R2009a矩陣實(shí)驗(yàn)室（美國MathWorks公司）軟件。

圖1 高光譜圖像系統(tǒng)的實(shí)景圖Fig. 1 Picture of the hyperspectral imaging system

1.2.2 訓(xùn)練集和驗(yàn)證集的選取

將獲取的高光譜圖像進(jìn)行黑白校正，校正方程[22]如公式（1）所示：

式中：I為校正后圖像；I0為從高光譜圖像系統(tǒng)中獲取的原始圖像；W為在與獲取樣本圖像相同拍攝條件下，拍攝標(biāo)準(zhǔn)白板獲得了白色參考圖像；為了去除相機(jī)傳感器暗電流的影響，在相機(jī)鏡頭關(guān)閉狀態(tài)下拍攝了暗參考圖像B。黑白校正后的圖像用于模型訓(xùn)練以及算法驗(yàn)證。

采用ENVI軟件中的感興趣區(qū)（region of interesting，ROI）工具手工選取訓(xùn)練集像素。任選10 個(gè)雞胴體高光譜圖像，隨機(jī)選取60 個(gè)污染區(qū)域（其中包括12 個(gè)盲腸糞便污染區(qū)域、17 個(gè)十二指腸糞便污染區(qū)域、18個(gè)血液污染區(qū)域、以及13 個(gè)膽汁污染區(qū)域）和57 個(gè)潔凈皮膚區(qū)域。污染區(qū)域和潔凈皮膚區(qū)域的平均面積分別為1 354 個(gè)像素和2 200 個(gè)像素。統(tǒng)計(jì)每一個(gè)區(qū)域中所有像素的平均光譜，將此117 個(gè)平均光譜作為訓(xùn)練集。訓(xùn)練集像素的光譜保存在一個(gè)矩陣中，矩陣中的前1 232 列代表樣本在1 232 個(gè)波長下的相對(duì)反射率，最后一列定義為識(shí)別函數(shù)值y。污染區(qū)域像素的y值指定為1，而非污染區(qū)域像素的y值指定為0；然后再利用ENVI的ROI工具，手工選取所有20 個(gè)高光譜圖像中的208 個(gè)污染區(qū)域，和胴體表面所有非污染區(qū)域作為驗(yàn)證集像素。提取每個(gè)驗(yàn)證集像素特征波長下的光譜作為驗(yàn)證集數(shù)據(jù)，檢驗(yàn)分類器的分類性能。

1.2.3 SPA-MLR-ROC分類器算法

1.2.3.1 SPA-MLR-ROC分類器算法流程

圖2 構(gòu)建及評(píng)估SPA-MLR- ROC分類器的流程圖Fig. 2 Flow chart of successive projections algorithm (SPA)-multivariable linear regression (MLR)-receiver operating characteristic curve (ROC) classifier algorithm and model evaluation

構(gòu)建及評(píng)估SPA-MLR-ROC分類器的整體流程如圖2所示。SPA選出特征波長，特征波長光譜與識(shí)別函數(shù)之間的對(duì)應(yīng)關(guān)系由MLR構(gòu)建，通過ROC方法選擇分類的最優(yōu)閾值，最后將SPA-MLR-ROC分類器識(shí)別效果與既有研究中的方法進(jìn)行比較。

1.2.3.2 SPA選擇特征波長

SPA是一種十分有效的變量選擇方法，可以從眾多的高光譜波長圖像中選取具有最小線性關(guān)系的波長，以此去除全波長模型的共線性[23]。

SPA包括兩個(gè)階段[24]，第1個(gè)階段是從任選的1 個(gè)波長開始，通過向剩余波長光譜矩陣連續(xù)投影，每次投影選出來的波長都與前面選出的波長的反射率之間具有最小的共線性。訓(xùn)練集中一共有117 個(gè)樣本光譜，因此選出的特征波長數(shù)量的最大值不能大于116。投影前任意選取的波長有1 232 種可能，因此會(huì)產(chǎn)生116×1 232 個(gè)備選波長鏈。

SPA的第2個(gè)過程是從所有的備選波長鏈中選出最優(yōu)集合。Galvao等[24]研究發(fā)現(xiàn)F快速檢測可以用來移除那些對(duì)預(yù)測結(jié)果作用不大的變量，以此增加模型的簡約性。計(jì)算包含不同數(shù)量波長的每種組合的均方根誤差（rootmean-square error，RMSE），如式（2）所示：

式中：xi為第i次投影選出的特征波長的光譜信息；x表示該組合所有波長下的平均光譜；n為此波長鏈的波長數(shù)量。選出比RMSE的最小值稍大的RMSE作為決定值，其對(duì)應(yīng)組合的波長數(shù)量的最小值即為通過SPA選出的特征波長的數(shù)量。

將訓(xùn)練集中污染區(qū)域的60 個(gè)平均光譜和非污染區(qū)域的57 個(gè)平均光譜輸入到SPA算法中，按照F檢驗(yàn)法進(jìn)行顯著性檢驗(yàn)，選擇顯著性水平α為0.25[24]，獲得通過SPA選出的特征波長。

1.2.3.3 MLR回歸建模

利用多元線性回歸模型構(gòu)建特征光譜與識(shí)別函數(shù)之間的回歸關(guān)系，如（3）式所示，

y=β0＋β1x1＋β1x2＋…βpxp（3）

式中：y為識(shí)別函數(shù)；xp為SPA選出的第p個(gè)特征波長下的相對(duì)反射率；βp為回歸系數(shù)；p為SPA選出的特征波長的數(shù)量。

1.2.3.4 ROC閾值選擇

將117 個(gè)訓(xùn)練集區(qū)域中所有像素的特征波長光譜值代入MLR模型后，實(shí)際獲得的識(shí)別函數(shù)預(yù)測值是處于某一區(qū)間范圍內(nèi)的離散量，但可能取值的數(shù)量遠(yuǎn)遠(yuǎn)大于待分類數(shù)量，需要設(shè)置一個(gè)閾值T進(jìn)行歸類處理。識(shí)別函數(shù)預(yù)測值如果大于閾值T，則相應(yīng)的分類結(jié)果被指定為y＝1，即被檢測的像素被判斷處于污染區(qū)域。反之，如果識(shí)別函數(shù)預(yù)測值小于閾值T，則分類器給出的分類結(jié)果為y = 0，則該像素被認(rèn)為出自非污染區(qū)域。

檢測結(jié)果的TPR[25]定義為所有檢測出的污染區(qū)域數(shù)量與實(shí)際污染區(qū)域數(shù)量的比值。TPR是按區(qū)域定義的，因?yàn)闄z測目標(biāo)是稀釋污染物，污染區(qū)域有很多稀釋污染物和皮膚的混合像元，不可能將一個(gè)區(qū)域中的所有像素都檢測出來，因此如果一個(gè)區(qū)域中有大于50 個(gè)像素被檢出，這個(gè)區(qū)域就被標(biāo)記為污染區(qū)域；如果一個(gè)像素實(shí)際為非污染區(qū)域像素，但卻被分類器判定為污染區(qū)域像素，這樣的像素稱為假陽性點(diǎn)。假陽性點(diǎn)的數(shù)量與實(shí)際非污染區(qū)域總的像素的比值定義為假陽性率（false positive rate，F(xiàn)PR）。FPR定義在像素級(jí)，是由于非污染區(qū)域是一個(gè)整體區(qū)域，無法將FPR定義在區(qū)域級(jí)。

性能優(yōu)越的分類器或分類模型應(yīng)具有高的TPR以及低的FPR。通常一個(gè)閾值的選擇需要在高TPR與低FPR之間找到一個(gè)平衡點(diǎn)。ROC圖像可以綜合反映敏感性和特異性的雙重屬性[26]。ROC圖像以TPR為縱坐標(biāo)，以FPR為橫坐標(biāo)，通過設(shè)定多個(gè)不同的閾值繪制成曲線或分立點(diǎn)，以此直觀、形象地揭示TPR和FPR的相互關(guān)系。在ROC圖像上，最靠近左上角的點(diǎn)（TPR=1，F(xiàn)PR=0）對(duì)應(yīng)的閾值即為最優(yōu)閾值。

1.2.3.5 檢測結(jié)果可視化

為了直觀地了解分類器的性能，將分類器檢測出的陽性點(diǎn)標(biāo)記出來，而非污染區(qū)不加以標(biāo)記，這樣檢測結(jié)果的可視化圖中就可以直觀地觀察到檢測出的污染區(qū)域的數(shù)量和位置?？梢暬木唧w做法如下[27]：根據(jù)驗(yàn)證集中所有像素的特征光譜信息，應(yīng)用SPA-MLR-ROC分類器進(jìn)行分類，將檢測出的陽性像素的位置坐標(biāo)信息輸入到MATLAB的plot函數(shù)，并將其標(biāo)記出來。檢測結(jié)果圖中會(huì)有一些由于皮膚褶皺或者組織連接處的陰影所造成的假陽性點(diǎn)，這些假陽性點(diǎn)通常面積較小。在MATLAB中，通過區(qū)域標(biāo)記的方法[28]，標(biāo)記出檢測結(jié)果圖中所有的連通區(qū)域，保留所有面積大于50 個(gè)像素的連通區(qū)域，這樣便可以去除那些面積較小的假陽性點(diǎn)[29]。

1.2.4 用于對(duì)比的檢測算法

1.2.4.1 波段比算法

既有研究中多采用主成分分析選擇特征波長，按照主成分載荷最大的原則選出兩個(gè)特征波長λ1為517 nm和λ2為565 nm[2-3,5-6,16,18-20,30]，在565/517 nm的波段比圖像中采用圖像處理識(shí)別污染物。然后利用單線性回歸確定第3個(gè)特征波長λ3為803 nm，其中B為辨別函數(shù)（對(duì)于污染區(qū)域B為1，非污染區(qū)域B為0），通過517 nm和802 nm波段比圖像消除假陽性點(diǎn)[3,16,19]。

1.2.4.2 雙波段算法

前期研究[21]提出的基于高光譜成像技術(shù)的雙波段檢測算法也是用來檢測雞胴體表面低可見污染物的方法。雙波段檢測算法利用一個(gè)特征波長下的圖像進(jìn)行區(qū)域生長，生長后的區(qū)域即為雞胴體的胸腹部，將其作為掩膜。掩膜與污染物可分辨度最大的特征波長圖像進(jìn)行“與”操作，提取出最大面積的雞胴體待檢ROI，最后利用標(biāo)記法識(shí)別出ROI內(nèi)有污染物存在的雞胴體，可以檢測出雞胴體表面81.6%的污染區(qū)域。

2 結(jié)果與分析

2.1 構(gòu)建的SPA-MLR-ROC模型

連續(xù)投影過程中，各個(gè)波長鏈的RMSE中最小值為0.256 99，比其稍大的RMSE的波長鏈中波長的數(shù)量為10，波長鏈中的波長為369、371、376、393、413、451、579、596、937、1 021 nm，即為通過SPA選出的特征波長。

通過多元線性回歸模型構(gòu)建的識(shí)別函數(shù)y與特征波長光譜之間的線性回歸關(guān)系如公式（4）所示：

式中：x為各個(gè)特征波長下的相對(duì)反射率。

圖3 不同閾值產(chǎn)生的ROC圖Fig. 3 ROC plot created with different thresholds.

將選出的117 個(gè)區(qū)域中所有像素的特征光譜值代入回歸方程，獲得的識(shí)別函數(shù)值范圍在-1.5～1.5之間。采用ROC分析方法，在-1.5～1.5之間每隔0.1選擇一個(gè)閾值，計(jì)算每一個(gè)閾值對(duì)應(yīng)的TPR 和 FPR，繪制ROC圖（圖3），可見，最優(yōu)閾值出現(xiàn)在T為1.3時(shí)，此時(shí)分類器獲得最優(yōu)性能（TPR=100%和FPR=0.43%）。

2.2 SPA-MLR-ROC分類器檢測結(jié)果

任選一個(gè)涂有污染物的雞胴體（圖4a），利用SPAMLR-ROC分類器檢測此胴體高光譜圖像中所有像素點(diǎn)，其檢測結(jié)果見圖4b?？梢?，所有污染區(qū)域均被檢測出來，圖像右下角的皮膚淤血也被檢測出來。

圖4 涂有污染物的雞胴體的彩色圖像（a）和通過SPA-MLR-ROC分類器獲得的檢測結(jié)果可視化圖（b）Fig. 4 Original image of contaminated carcass and visual results obtained by SPA-MLR-ROC classifier

SPA-MLR-ROC分類器檢測所有20 個(gè)雞胴體表面的208 個(gè)污染區(qū)域，結(jié)果顯示204 個(gè)污染區(qū)域被檢出，污染區(qū)域的TPR為98.08%，在像素級(jí)別FPR為0.385%。

2.3 既有算法的檢測結(jié)果

2.3.1 波段比算法的檢測結(jié)果

對(duì)訓(xùn)練集數(shù)據(jù)進(jìn)行主成分分析，分析結(jié)果顯示前3 個(gè)主成分占總體光譜變量的97.41%，提取出載荷最大的3 個(gè)波長：400、452 nm和679 nm。3 個(gè)特征波長相互組合，構(gòu)成波段比圖像，其中污染區(qū)域識(shí)別度最高的是679/452 nm。單線性回歸選取出去除噪聲的特征波長為516 nm。利用波段比算法檢測圖4a中雞胴體的檢測結(jié)果如圖5a所示。

圖5 兩個(gè)既有算法的檢測結(jié)果圖Fig. 5 Results obtained by two existing algorithms

利用波段比算法檢測所有20 個(gè)雞胴體，污染區(qū)域檢測的TPR為80.15%，像素級(jí)別FPR為0.042%。2.3.2 雙波段算法的檢測結(jié)果

按文獻(xiàn)[21]提出的可識(shí)別度定義，計(jì)算所有波長下圖像的可識(shí)別度，得到可識(shí)別度最大的波長為400 nm；掩膜選擇參數(shù)最大值出現(xiàn)在675 nm。利用雙波段算法檢測圖4a中雞胴體的檢測結(jié)果如圖5b所示。

對(duì)于20 個(gè)驗(yàn)證集中的雞胴體，雙波段算法檢測污染區(qū)域的TPR為82.37%，像素級(jí)別FPR為0.895%。

3 討 論

性能優(yōu)越的檢測算法應(yīng)該具有高的TPR，以及低的FPR。因此用3 種算法（SPA-MLR-ROC分類器、波段比算法以及雙波段檢測算法）分別檢測驗(yàn)證集中所有的雞胴體樣本，比較三者的檢測性能。

統(tǒng)計(jì)3 種檢測方法對(duì)每只雞檢測的TPR和FPR，20 個(gè)雞胴體檢測結(jié)果的統(tǒng)計(jì)分析如表1所示。由表可見：1）SPA-MLR-ROC分類器的區(qū)域TPR遠(yuǎn)高于波段比算法和雙波段檢測算法，原因是SPA選取的10 個(gè)特征波長保留了更多的污染區(qū)域特征光譜信息，使得污染區(qū)域TPR大大提高；并且特征光譜信息能最大程度地區(qū)分污染區(qū)域和皮膚區(qū)域，使得SPA-MLR-ROC分類器可以將整個(gè)雞胴體作為檢測范圍，將胸腹部以外（如雞腿部）的污染區(qū)域也檢測出來，而不必?fù)?dān)心雞胸腹部與翅膀以及雞腿連接處的陰影造成的假陽性點(diǎn)。對(duì)比圖4和圖5b，可見，SPAMLR-ROC分類器檢出了雙波段檢測算法沒有檢出的雞腿部的十二指腸糞便和盲腸糞便污染，以及雞身上的淤血（圖5b黑色背景區(qū)域中用虛線橢圓標(biāo)出的3 個(gè)區(qū)域）；SPA-MLR-ROC分類器也檢出了波段比算法沒有檢出的雞腿部的十二指腸糞便，以及胴體胸腹部的膽汁污染區(qū)以及淤血區(qū)域（圖5a中虛線橢圓標(biāo)出的5 個(gè)區(qū)域）。2）SPA-MLR-ROC分類器的像素級(jí)FPR僅有0.385%，不到雙波段算法FPR的一半。雖然SPA-MLR-ROC分類器的FPR稍大于波段比算法，但是SPA-MLR-ROC分類器的TPR也遠(yuǎn)大于波段比算法，而作為一種識(shí)別算法，高的TPR尤為重要。3）SPA-MLR-ROC分類器的檢測結(jié)果中，TPR和FPR的標(biāo)準(zhǔn)偏差都小于波段比算法以及雙波段算法，證明了SPA-MLR-ROC分類器檢測性能更加穩(wěn)定。

綜上所述，SPA-MLR-ROC分類器算法與雙波段檢測算法以及波段比算法相比，不僅具有高TPR同時(shí)具有低FPR，因此是一個(gè)性能穩(wěn)定、檢測效果好的分類器。

表1 SPA-MLR-ROC分類器與波段比算法以及雙波段檢測算法在TPR和FPR方面的檢測結(jié)果對(duì)比Table 1 Comparison in detection results of the SPA-MLR-ROC classifier and the band ratio algorithm as well as the dual-band algorithm in terms of true positive rate (TPR) and false positive rate (FPR)

4 結(jié) 論

本實(shí)驗(yàn)提出了一種用于精確提取雞胴體表面污染物特征的SPA-MLR-ROC分類器方法，研究了SPAMLR-ROC分類器檢測雞胴體表面稀釋污染物的可行性。研究發(fā)現(xiàn)SPA從1 232 個(gè)高光譜波長中選取的10 個(gè)（369、371、376、393、413、451、579、596、937、1 021 nm）最優(yōu)的特征波長，能夠精確提取污染物的特征。SPA-MLR-ROC分類器能夠根據(jù)高光譜圖像中的每一個(gè)像素的10 個(gè)特征光譜信息，給出每一個(gè)像素的識(shí)別結(jié)果。SPA-MLR-ROC分類器對(duì)20 個(gè)雞胴體表面的污染區(qū)域的TPR可以達(dá)到98.08%，而檢出的FPR僅有0.385%。相比波段比算法大大提高了檢測TPR；相比雙波段檢測算法，擴(kuò)大了檢測范圍、提高了檢測精度。由此可見，采用SPA-MLR-ROC分類器可以更好地實(shí)現(xiàn)對(duì)雞胴體表面低可見污染物的精確檢測。本實(shí)驗(yàn)注重方法研究，其檢測結(jié)果的穩(wěn)定性還需要更多的雞胴體樣本加以驗(yàn)證，同時(shí)還需要繼續(xù)研究在實(shí)際生產(chǎn)線上利用多光譜相機(jī)代替高光譜相機(jī)后的檢測速度與檢測效果。

[1] 夏芳. 白羽肉雞行業(yè)迎來復(fù)蘇[J]. 農(nóng)產(chǎn)品市場周刊, 2016(3)∶ 15-17.DOI∶10.3969/j.issn.1009-8070.2016.03.009.

[2] PARK B, LAWRENCE K C, WINDHAM W R, et al. Hyperspectral imaging for detecting fecal and ingesta contaminants on poultry carcasses[J]. Transactions of the Asae, 2002, 45(6)∶ 2017-2026.DOI∶10.13031/2013.11413.

[3] HEITSCHMIDT G, PARK B, LAWRENCE K, et al. Improved hyperspectral imaging system for fecal detection on poultry carcasses[J]. Transactions of the Asabe, 2007, 50(4)∶ 1427-1432.DOI∶10.13031/2013.23615.

[4] WINDHAM W R, HEITSCHMIDT G W, LAWRENCE K C, et al.Effect of spectrally mixed hyperspectral image pixels on detection of cecal contaminated broiler carcasses[J]. International Journal of Poultry Science, 2007, 6(12)∶ 955-959. DOI∶10.3923/ijps.2007.955.959.

[5] 趙進(jìn)輝, 涂冬成, 歐陽靜怡, 等. 利用高光譜圖像技術(shù)檢測雞胴體內(nèi)部糞便污染物[J]. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào), 2011, 33(3)∶ 573-577.DOI∶10.3969/J.ISSN.1000-2286.2011.03.030.

[6] 趙進(jìn)輝, 吁芳, 吳瑞梅, 等. 基于分段主成分分析與波段比的雞胴體表面糞便污染物檢測[J]. 激光與光電子學(xué)進(jìn)展, 2011(7)∶ 163-167.DOI∶10.3788/lop48.073001.

[7] MUSGROVE M T, BERRANG M E, BYRD J A, et al. Detection of Campylobacter spp. in ceca and crops with and without enrichment[J].Poultry Science, 2001, 80(6)∶ 825-828. DOI∶10.1093/PS/80.6.825.

[8] BERRANG M E, BUHR R J, CASON J A. Campylobacter recovery from external and internal organs of commercial broiler carcass prior to scalding[J]. Poultry Science, 2000, 79(2)∶ 286-290. DOI∶10.1093/PS/79.2.286.

[9] BERRANG M E, SMITH D P, WINDHAM W R, et al. Effect of intestinal content contamination on broiler carcass campylobacter counts[J]. Journal of Food Protection, 2004, 67(2)∶ 235-238.DOI∶10.4315/0362-028X-67.2.235.

[10] JIANG F C, LONG Y, TANG X Y, et al. Detection of the total viable counts in chicken based on visible/near-infrared spectroscopy[C]//Proceedings of spie-the international society for optical engineering.International society for optics and photonics, 2014, 9108∶ 91080S.DOI∶10.1117/12.2053201.

[11] XIONG Z J, SUN D W, DAI Q, et al. Application of visible hyperspectral imaging for prediction of springiness of fresh chicken meat[J]. Food Analytical Methods, 2015, 8(2)∶ 380-391. DOI∶10.1007/S12161-014-9853-3.

[12] 趙杰文, 惠喆, 黃林, 等. 高光譜成像技術(shù)檢測雞肉中揮發(fā)性鹽基氮含量[J]. 激光與光電子學(xué)進(jìn)展, 2013(7)∶ 158-164. DOI∶10.3788/LOP50.073003.

[13] XIONG Z J, SUN D W, XIE A G, et al. Potential of hyperspectral imaging for rapid prediction of hydroxyproline content in chicken meat[J]. Food Chemistry, 2015, 175∶ 417-422. DOI∶10.1016/j.foodchem.2014.11.161.

[14] G?HLER D, FISCHER B, MEISSNER S. In-ovo sexing of 14-day-old chicken embryos by pattern analysis in hyperspectral images (vis/Nir spectra)∶a non-destructive method for layer lines with gender-specific down feather color[J]. Poultry Science, 2017, 96(1)∶ 1-4. DOI∶10.3382/ps/pew282.

[15] PARK B, LAWRENCE K C, WINDHAM W R, et al. Multispectral imaging system for fecal and ingesta detection on poultry carcasses[J].Journal of Food Process Engineering, 2004, 27(5)∶ 311-327.DOI∶10.1111/j.1745-4530.2004.00464.x.

[16] WINDHAM W R, SMITH D P, BERRANG M E, et al. Effectiveness of hyperspectral imaging system for detecting cecal contaminated broiler carcasses[J]. International Journal of Poultry Science, 2005,4(9)∶ 657-662. DOI∶10.3923/ijps.2005.657.662.

[17] YOON S C, LAWRENCE K C, PARK B, et al. Optimization of fecal detection using hyperspectral imaging and kernel density estimation[J]. Transactions of the ASABE, 2007, 50(3)∶ 1063-1071.DOI∶10.13031/2013.23120.

[18] PARK B, YOON S-C, WINDHAM W R, et al. Line-scan hyperspectral imaging for real-time in-line poultry fecal detection[J].Sensing and Instrumentation for Food Quality and Safety, 2011, 5(1)∶25-32. DOI∶10.1007/s11694-011-9107-7.

[19] YOON S C, PARK B, LAWRENCE K C, et al. Line-scan hyperspectral imaging system for real-time inspection of poultry carcasses with fecal material and ingesta[J]. Computers and Electronics in Agriculture, 2011, 79(2)∶ 159-168. DOI∶10.1016/j.compag.2011.09.008.

[20] WINDHAM W R, HEITSCHMIDT G W, SMITH D P, et al. Detection of ingesta on pre-chilled broiler carcasses by hyperspectral imaging[J].International Journal of Poultry Science, 2005, 43(12)∶ 10157-10163.DOI∶10.3923/ijps.2005.959.964.

[21] 吳威, 陳桂云, 夏建春, 等. 雞胴體表面低可見污染物的雙波段檢測方法研究[J]. 光譜學(xué)與光譜分析, 2014, 34(12)∶ 3363-3367.DOI∶10.3964/j.issn.1000-0593(2014)12-3363-05.

[22] 鄒小波, 李志華, 石吉勇, 等. 高光譜成像技術(shù)檢測肴肉新鮮度[J]. 食品科學(xué), 2014, 35(8)∶ 89-93. DOI∶10.7506/spkx1002-6630-201408017.

[23] ARAúJO M C U, SALDANHA T C B, GALV?O R K H, et al. The successive projections algorithm for variable selection in spectroscopic multicomponent analysis[J]. Chemometrics & Intelligent Laboratory Systems, 2001, 57(2)∶ 65-73. DOI∶10.1016/S0169-7439(01)00119-8.

[24] GALV?O R K H, ARAúJO M C U, FRAGOSO W D, et al. A variable elimination method to improve the parsimony of MLR models using the successive projections algorithm[J]. Chemometrics and Intelligent Laboratory Systems, 2008, 92(1)∶ 83-91. DOI∶10.1016/j.chemolab.2007.12.004.

[25] YUAN Y, FANG J W, WANG Q. Online anomaly detection in crowd scenes via structure analysis[J]. IEEE Transactions on Cybernetics,2015, 45(3)∶ 548-561. DOI∶10.1109/TCYB.2014.2330853.

[26] METZ C E. Basic principles of roc analysis[J]. Seminars in Nuclear Medicine, 1978, 8(4)∶ 283-298. DOI∶10.1016/S0001-2998(78)80014-2.

[27] WU D, SHI H, WANG S J, et al. Rapid prediction of moisture content of dehydrated prawns using online hyperspectral imaging system[J].Analytica Chimica Acta, 2012, 726(9)∶ 57-66. DOI∶10.1016/j.aca.2012.03.038.

[28] STEFANO L D, BULGARELLI A. A simple and efficient connected components labeling algorithm[C]//Image analysis and processing,1999∶ 322-327. DOI∶10.1109/ICIAP.1999.797615.

[29] WU W, CHEN G Y, WU M Q, et al. Detection of diluted contaminants on chicken carcasses using a two-dimensional scatter plot based on a two-dimensional hyperspectral correlation spectrum[J]. Applied Optics, 2017, 56(9)∶ 72-78. DOI∶10.1364/AO.56.000D72.

[30] PARK B, LAWRENCE K C, WINDHAM W R, et al. Assessment of hyperspectral imaging system for poultry safety inspection[J].Proceedings of Spie, 2003, 4879∶ 269-279. DOI∶10.1117/12.462394.

Successive Projections Algorithm-Multiple Linear Regression-Receiver Operating Characteristic Analysis for Diluted Contaminant Identification on Chicken Carcasses

WU Wei, WU Mingqing, CHEN Guiyun, YU Zhenwei, CHEN Kunjie*
(College of Engineering, Nanjing Agricultural University, Nanjing 210031, China)

This paper presents a method for the identification of diluted contaminants on the surface of chicken carcasses based on successive projections algorithm (SPA)-multiple linear regression (MLR)-receiver operating characteristic (ROC)classifier. Firstly, a total of 20 images of carcasses with diluted contaminants were acquired by hyperspectral imaging system, and 10 characteristic bands were extracted from the 1 232 bands by SPA. Then the MLR method was used to construct a regression model between the discriminant function and the characteristic spectral bands. Finally, the optimal classification threshold with high true positive rate (TPR) and low false positive rate (FPR) was determined by ROC analysis. Thus, the SPA-MLR-ROC classifier allowed the identification of the diluted contaminants. The results showed that the TPR of the SPA-MLR-ROC classifier was 98.08% and the FPR was only 0.39%. The detection accuracy was higher than that of the band ratio algorithm and the dual-band algorithm. Hence, the SPA-MLR-ROC classifier exhibited good performance for the detection of diluted contaminants on the surface of chicken carcass. However,because of the limited number of samples, further study using more samples is needed to verify the stability and feasibility of this method.

successive projections algorithm-multiple linear regression-receiver operating characteristic analysis; classifier;contaminants; hyperspectral imaging; chicken carcass

DOI∶10.7506/spkx1002-6630-201724040

S123；TP391.4

A

1002-6630（2017）24-0247-06

吳威, 吳明清, 陳桂云, 等. 雞胴體表面稀釋污染物的連續(xù)投影-多元線性回歸-受試者特性分析檢測[J]. 食品科學(xué), 2017,38(24)∶ 247-252. DOI∶10.7506/spkx1002-6630-201724040. http∶//www.spkx.net.cn

2017-04-26

江蘇省科技支撐計(jì)劃項(xiàng)目（BE2014708）；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)（KYZ201665）

吳威（1978—），女，講師，博士，研究方向?yàn)檗r(nóng)產(chǎn)品無損檢測。E-mail：wuwei@njau.edu.cn

*通信作者：陳坤杰（1963—），男，教授，博士，研究方向?yàn)檗r(nóng)產(chǎn)品檢測與加工技術(shù)。E-mail：kunjiechen@njau.edu.cn

WU Wei, WU Mingqing, CHEN Guiyun, et al. Successive projections algorithm-multiple linear regression-receiver operating characteristic analysis for diluted contaminant identification on chicken carcasses[J]. Food Science, 2017, 38(24)∶247-252. (in Chinese with English abstract) DOI∶10.7506/spkx1002-6630-201724040. http∶//www.spkx.net.cn