孫翹楚，肖創(chuàng)柏，郭 犇，蘇 峰

(1.北京工業(yè)大學(xué)計(jì)算機(jī)學(xué)院，北京100124，sunqiaochu@hotmail.com;2.中國(guó)農(nóng)業(yè)銀行總行軟件開(kāi)發(fā)中心，北京100073)

近年來(lái)，紅外光譜技術(shù)在食品鑒別方面的應(yīng)用研究發(fā)展迅速.紅外光譜技術(shù)的優(yōu)勢(shì)在于樣品幾乎不需要分離處理，從而使得鑒別過(guò)程更加迅速［1］.同時(shí)根據(jù)樣品成分的紅外吸收率，紅外光譜技術(shù)能夠提供整體信息來(lái)鑒別樣品.紅外光譜根據(jù)不同的波數(shù)范圍分為近紅外(12 900～4 000 cm－1)、中紅外區(qū)(4 000～650 cm－1)和遠(yuǎn)紅外區(qū)(650～10 cm－1)3個(gè)光區(qū).其中，中紅外區(qū)是應(yīng)用最早和最廣的一個(gè)區(qū)，大多數(shù)有機(jī)化合物和無(wú)機(jī)離子的基頻振動(dòng)吸收都落在中紅外區(qū)內(nèi).化合物分子中所含的化學(xué)鍵或官能團(tuán)不同，或官能團(tuán)所處的化學(xué)環(huán)境不同，其振動(dòng)能級(jí)從基態(tài)躍遷到激發(fā)態(tài)所需能量不同會(huì)吸收不同波長(zhǎng)的紅外光，形成不同的紅外光譜圖，而相關(guān)吸收峰的高度與其含量成正比.因此，紅外光譜法被廣泛用于有機(jī)物的定性和定量研究.由于傅里葉變換紅外光譜法(Fourier transform infrared spectroscopy，F(xiàn)TIR)具有整體、快速、無(wú)損鑒定復(fù)雜體系多種組份的特點(diǎn)，目前己在食品鑒別領(lǐng)域獲得了廣泛應(yīng)用［2－3］.

劉沭華等［4］通過(guò)對(duì)幾種中藥的傅里葉變換紅外光譜進(jìn)行求一階導(dǎo)數(shù)對(duì)其進(jìn)行區(qū)分.但是在本文初期實(shí)驗(yàn)中，這種方法并不適合食用植物油的紅外光譜.對(duì)其求一階導(dǎo)數(shù)并不足以對(duì)其區(qū)分.同時(shí)對(duì)各個(gè)光譜的指紋區(qū)進(jìn)行了積分操作，提取的特征差別細(xì)微，不足以作為區(qū)分的標(biāo)準(zhǔn).賈素貞等［5］通過(guò)對(duì)幾種食用植物油的二維傅里葉變換紅外光譜進(jìn)行特征提取，對(duì)其進(jìn)行區(qū)分.這種方法可以很好的放大一維光譜的差異，但是計(jì)算復(fù)雜，對(duì)設(shè)備要求較高.本文提出一種利用分形理論提取食用植物油紅外光譜特征的算法，將此作為一種主要特征引入食用油紅外光譜模式識(shí)別中，在實(shí)驗(yàn)中取得了較好效果.

1 基于分形理論的紅外光譜特征提取

1.1 盒維數(shù)定義

分形學(xué)的概念的實(shí)質(zhì)是除傳統(tǒng)幾何學(xué)以外的不規(guī)則幾何體在標(biāo)度變換下的自相似性.

若F是分形集，則其應(yīng)具有以下性質(zhì)［6］:

1)F具有精細(xì)結(jié)構(gòu)，即有任意小比例細(xì)節(jié);

2)F極不規(guī)則，以致其整體與局部都無(wú)法用傳統(tǒng)幾何語(yǔ)言描述;

3)F通常具有某種自相似性;

4)F的分形維數(shù)大于其拓?fù)渚S數(shù);

5)F可以用非常簡(jiǎn)單的方法定義，可以由迭代產(chǎn)生.

一般而言若研究對(duì)象符合上述性質(zhì)中的全部或大多數(shù)，即使有某個(gè)性質(zhì)例外，也并不影響稱其為分形.

計(jì)盒維數(shù)或稱盒維數(shù)(Box-counting or Box dimension)是應(yīng)用最廣泛的維數(shù)之一，它的普遍應(yīng)用主要由于這種盒維數(shù)的數(shù)學(xué)計(jì)算以及經(jīng)驗(yàn)估計(jì)相對(duì)容易.

盒維數(shù)的定義為:設(shè)F為Rn中任一非空子集，記N(F，δ)表示最大直徑為δ且能覆蓋F的集合的最小數(shù)，則F的上下盒維數(shù)分別定義為

如果上下維數(shù)相等，則稱這個(gè)共同的值F的計(jì)盒維數(shù)或者盒維數(shù)，定義為

盒維數(shù)存在一些等價(jià)定義，有時(shí)這些定義更適合應(yīng)用.通常為以下五種［7］:

1)覆蓋F的半徑為δ的最小閉球數(shù);

2)覆蓋F的邊長(zhǎng)為δ的最小立方塊數(shù);

3)相交于F的δ－網(wǎng)格塊數(shù);

4)覆蓋F的直徑至多為δ的最小集合數(shù);

5)球心在F中半徑為δ的互不相交球的最大數(shù).

本文利用定義(3)對(duì)紅外光譜的分形維數(shù)進(jìn)行計(jì)算.

1.2 盒維數(shù)的計(jì)算

設(shè)A∈H(X)，其中，(X，d)是度量空間，εn=cbn，0＜b＜1，c＞0和n=1，2，3，….

如果

則A具有分形維數(shù)D.在式(4)中用γ代替ε，計(jì)算分形維數(shù)的公式為

故可讓ln N(A，γ)與ln γ成為一對(duì)變量，則式(5)變?yōu)?/p>

于是D(A)就成了(ln γ，ln N(A，γ))坐標(biāo)系統(tǒng)中的斜率［8］.本文借以不同紅外光譜盒維數(shù)的不同，對(duì)其所代表的食用油進(jìn)行區(qū)分.

1.3 紅外光譜分形維數(shù)計(jì)算方法

紅外光譜在計(jì)算機(jī)中以圖像的格式存儲(chǔ)，文中使用Intel公司開(kāi)發(fā)的OpenCV開(kāi)源計(jì)算機(jī)視覺(jué)庫(kù)，實(shí)現(xiàn)了紅外光譜的分形維數(shù)計(jì)算.具體其分形維數(shù)計(jì)算過(guò)程如圖1所示.

圖1 紅外線光譜盒維數(shù)計(jì)算流程圖

1)首先讀入光譜圖像，對(duì)光譜圖像進(jìn)行二值化處理［9］.原始紅外光譜圖像為彩色圖像，進(jìn)行盒維數(shù)計(jì)算時(shí)，應(yīng)先將彩色圖像轉(zhuǎn)換成二值圖像(如圖2所示).利用OpenCV的on－trackbar()函數(shù)生成滑動(dòng)條，動(dòng)態(tài)設(shè)定閾值，閾值的取值以能顯示全部光譜圖像為準(zhǔn).然后利用cvThreshold()函數(shù)將讀入圖片二值化.包含光譜的像素點(diǎn)值為0.

圖2 二值化后的光譜圖像(局部)及其對(duì)應(yīng)的二值數(shù)據(jù)

2)將得到的二值圖像轉(zhuǎn)存為矩陣.運(yùn)用cv-Convert()函數(shù)將二值圖像轉(zhuǎn)存為與其同尺寸的矩陣I.

3)把得到的矩陣依次劃分為若干塊，每一塊的行數(shù)與列數(shù)均為 k，k取 1，2，4，…，2n，運(yùn)用cvGetReal2D()函數(shù)掃描每個(gè)劃分塊，把包含像素0的塊的個(gè)數(shù)記為Nk.這樣就得到k與Nk的對(duì)應(yīng)數(shù)組B［10］.

4)對(duì)數(shù)組B中的元素求對(duì)數(shù).并在雙對(duì)數(shù)坐標(biāo)平面中，以一元線性回歸用直線擬合數(shù)據(jù)點(diǎn)(ln k，ln Nk)，k=1，2，4，…，2n，所得到的斜率的負(fù)值D就是該紅外光譜的盒維數(shù).

2 實(shí)驗(yàn)結(jié)果及分析

2.1 算法正確性檢驗(yàn)

通過(guò)實(shí)驗(yàn)發(fā)現(xiàn)，分形維數(shù)的正確性與k的取值以及圖像的分辨率有關(guān).在實(shí)驗(yàn)中取1，2，4，8，16，32，64，128，256，這9個(gè)值對(duì)圖像計(jì)算分形.

根據(jù)盒維數(shù)的定義，所取k的尺寸越小，所計(jì)算出的分形維數(shù)越精確.但由于對(duì)存儲(chǔ)在計(jì)算機(jī)中的圖像進(jìn)行分形計(jì)算，只能分割到像素點(diǎn)尺寸.因此計(jì)算分形算法所得出的結(jié)果與實(shí)際分形維數(shù)之間存在誤差.而誤差的大小直接影響到特征提取的準(zhǔn)確性，若計(jì)算出的誤差大于光譜分形維數(shù)的差異，則無(wú)法利用此差異對(duì)光譜進(jìn)行區(qū)分.為驗(yàn)證分形維數(shù)計(jì)算算法的準(zhǔn)確性和可靠性，分別用直線、矩形、Koch曲線、Cantor三分集、Sierpinski三角墊和Sierpinski地毯等已知維度的分形圖形對(duì)算法進(jìn)行驗(yàn)證.結(jié)果如表1所示.

表1 分形維數(shù)計(jì)算結(jié)果

由實(shí)驗(yàn)結(jié)果可知，誤差在小數(shù)點(diǎn)后三位到四位之間，而光譜分形維數(shù)的差異在小數(shù)點(diǎn)后兩位.因此，誤差不會(huì)對(duì)特征提取造成影響.

2.2 實(shí)驗(yàn)分析

采用本文的算法，本文對(duì)市面上常見(jiàn)的福臨門(mén)、金龍魚(yú)、歐尚等品牌的幾種食用油的紅外光譜的分形維數(shù)進(jìn)行了計(jì)算.通過(guò)對(duì)每種食用油各30個(gè)樣本的盒維數(shù)計(jì)算(抽樣結(jié)果如表2所示，其中，k為窗口大小，Nk為盒子數(shù))，可知食用油的盒維數(shù)總體在1.293 1～1.315 3之間.可利用此特征對(duì)食用油同其他液體進(jìn)行區(qū)分.另外通過(guò)對(duì)玉米油和芝麻油的紅外光譜分形維數(shù)的計(jì)算得到，玉米油的紅外線光譜的分形維數(shù)在1.293 1～1.297 2之間，而芝麻油的紅外光譜分形維數(shù)在1.311 3～1.315 3之間，分形維數(shù)差別明顯，因此可利用此特征對(duì)兩種植物油的真?zhèn)芜M(jìn)行鑒別.

表2 食用植物油紅外光譜盒維數(shù)計(jì)算結(jié)果

3 結(jié)論

1)分形維數(shù)能有效地表現(xiàn)紅外光譜的復(fù)雜度以及自相似性特征.

2)分形維數(shù)的準(zhǔn)確性與源圖像的分辨率有關(guān)，在實(shí)驗(yàn)中采用的光譜圖像分辨率為1280× 880.隨著儀器設(shè)備的不斷進(jìn)步，更大分辨率的光譜圖像可以得到更為準(zhǔn)確的分形維數(shù)特征.

3)通過(guò)引入其它特征，如譜峰高度、光譜導(dǎo)數(shù)、面積等特征，可以對(duì)其它食用油進(jìn)行更好的區(qū)分.

4)通過(guò)對(duì)市面上幾種常見(jiàn)食用油的紅外線光譜的分形維數(shù)的計(jì)算表明，利用分形維數(shù)可以對(duì)玉米油以及芝麻油的真?zhèn)芜M(jìn)行較好的區(qū)分.

［1］Van de VOORT F R，GHETLER A，GARCIA-GONZA'LEZ D L.Perspectives on quantitative mid-ftir spectroscopy in relation to edible oil and lubricant analysis:Evolution and integration of analytical methodologies［J］.Food Analytical Methods，2008，3(1):153－163.

［2］MILLS B L，Van de VOORT F R，KAKUDA Y.The quantitative analysis of fat and protein in meat by transmission infrared analysis［J］.Meat Science，1984，11(4):253－262.

［3］Van de VOORT F R，SEDMAN J，ISMAIL A A.A rapid FTIR quality control method for determining fat and moisture in high-fat products［J］.Food Chemistry，1993，48(2):213－221.

［4］劉沭華，張學(xué)工，周群，等.近紅外漫反射光譜法和模式識(shí)別技術(shù)鑒別中藥材產(chǎn)地［J］.光譜學(xué)與光譜分析，2006，26(4):629－632.

［5］賈素貞，鮑猛，聶傳平.幾種食用油的二維紅外相關(guān)分析［J］.濟(jì)南大學(xué)學(xué)報(bào)，2007，21(2):124－129.

［6］BRIGGS J.Fractals:The Patterns of Chaos［M］.London:Thames and Hudson，1992:148.

［7］FALCONER K.Techniques in Fractal Geometry［M］.NJ，USA:John Willey and Sons，1997:59－61.

［8］張平，王新柯，李海濤，等.基于分形理論的太赫茲光譜識(shí)別［J］.量子電子學(xué)報(bào)，2007，24(6):672－676.

［9］文志英.分形幾何的數(shù)學(xué)基礎(chǔ)［M］.上海:上?？萍冀逃霭嫔?，2002.

［10］CHAUDHURI B B，SARKAR N.Texture segmentation using fractal dimension［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，1995，17(1): 72－77.

［11］王韶娟，曾國(guó)蓀.分形維數(shù)的一個(gè)并行算法［J］.計(jì)算機(jī)應(yīng)用與軟件，2005，22(10):19－20.