李 玲，宗 力，劉梅英，張 松

摘要：通過拍攝常用預混料有機載體玉米粉、全脂米糠粉、全麥粉、豆粕粉、膨化大豆粉顆粒不同放大倍數(shù)的掃描電子顯微鏡（SEM）圖像，利用基于周界的周長-面積法的Matlab程序計算其周長-面積分形維數(shù)。結(jié)果表明，周長-面積雙對數(shù)擬合線的相關(guān)系數(shù)都在0.99以上，且不同有機載體分形維數(shù)的均值存在顯著差異。采用周長-面積法計算玉米粉、全脂米糠粉、全麥粉、豆粕粉等具有統(tǒng)計自相似結(jié)構(gòu)的預混料載體是可行的，用以描述幾何形體不規(guī)則程度的分形維數(shù)可以作為表征預混料有機載體的一個參量。

關(guān)鍵詞：周長-面積法；分形維數(shù)；有機載體；預混料

中圖分類號：TP391.41文獻標識碼：A文章編號：0439-8114（2012）20-4634-04

自從Mandelbrot引入分形和分形維數(shù)的概念以后，分形理論以其研究不規(guī)則復雜形體的優(yōu)勢而被廣泛應用于自然科學和社會科學中，并衍生出眾多的分形維數(shù)測定方法。根據(jù)分形維數(shù)測定方法的不同，分形維數(shù)又分為Hausdorff維數(shù)、相似維數(shù)、盒維數(shù)、容量維數(shù)、關(guān)聯(lián)維數(shù)、信息維數(shù)、譜維數(shù)、周長-面積分形維數(shù)等[1，2]。不同的分形對象所適應的分形維數(shù)測定方法亦有所不同，例如，具有分形特性的不規(guī)則島形圖案的面積、周長存在固定的對應關(guān)系，通過這種對應關(guān)系即可計算其分形維數(shù)，利用這種對應關(guān)系計算其分形維數(shù)的方法即為周長-面積法。周長-面積法在分形研究中的應用較廣泛，例如，溫芝元[3]以周長-面積法研究了3類水果的花萼面和側(cè)面分形維數(shù)；曹樂平[4]將周長、面積和周長-面積分形維數(shù)應用于柑橘品種機器識別，品種正確識別率分別達到95%以上；李金萍等[5]應用改進的周長-面積法對鋁氫氧化顆粒的分形維數(shù)進行了計算；譚佐軍等[6]指出貯藏物昆蟲周長-面積分形維數(shù)反映了貯藏物昆蟲形狀的規(guī)則性和復雜程度；龔國淑等[7]利用周長-面積法和曲線長度法計算了4種植物病斑的分形維數(shù)，結(jié)果表明周長-面積法測得的分形維數(shù)相關(guān)系數(shù)更高。

預混料載體用于承載微量元素、維生素等飼料中不可或缺的微量組分，在預混料的分形分析方面已經(jīng)有了一些初步的探索[8，9]，研究以預混料常用載體為研究對象，將分形維數(shù)計算的周長-面積法引入預混料的分形分析中，以表征不同預混料載體微觀結(jié)構(gòu)的復雜程度。

1.2.2載體顆粒掃描電子顯微鏡（SEM）圖像的拍攝采用30目和80目的標準分樣篩對物料進行手動篩分，取30～80目的顆粒作為載體樣品，稱取5.00 g篩分后的載體樣品置于MC01000225電熱恒溫鼓風干燥箱中，在80 ℃條件下干燥8 h后置于干燥皿中冷卻。

待樣品冷卻后用牙簽取部分樣品顆粒粘到取樣膠帶上，使得膠帶上的樣品薄而均勻地分散開，經(jīng)JEOLJFC-1600離子濺射儀噴金處理之后再通過JSM-6390/LV掃描電子顯微鏡觀察并提取樣品顆粒的電子顯微圖像，對于視野內(nèi)清晰、獨立的每個樣品顆粒依次提取6個不同放大倍數(shù)的圖像。

1.2.3Matlab程序的編寫與檢驗采用Matlab軟件，根據(jù)周長-面積法計算分形維數(shù)的原理編寫分形維數(shù)計算程序，程序框架如圖1所示，該程序包括圖像處理和周長、面積、分形維數(shù)計算等幾個部分，周長和面積皆以像素個數(shù)來表示。

采用徐永春等[11]用Matlab編寫的Koch雪花繪制程序繪制經(jīng)過5次迭代的Koch雪花曲線，利用nEO iMAGING圖像處理軟件對該Koch雪花曲線進行放大，按照圖像的放大倍數(shù)的大小順序依次保存為6個Jpg文件。采用編寫的分形維數(shù)計算程序計算Koch雪花曲線的分形維數(shù)，并與Koch曲線的理論分形維數(shù)比較。

1.2.4預混料常用有機載體周長-面積分形維數(shù)的測定基于周長-面積法計算有機載體的粉體顆粒的分形維數(shù)，使用nEO iMAGING圖像處理軟件提取出單獨的樣品顆粒，再用Matlab軟件讀取對每個粉體顆粒提取的6個放大倍數(shù)的圖像，進行圖像處理，分別計算其周長和面積，得到一系列的周長{P1，P2，P3，P4，P5，P6}和面積{A1，A2，A3，A4，A5，A6}，利用這些數(shù)據(jù)，采用最小二乘法分析、擬合并繪制周長、面積的雙對數(shù)擬合線，求得直線的斜率，則粉體顆粒的面積-周長分形維數(shù)即為斜率值的2倍。

2結(jié)果與分析

經(jīng)檢驗發(fā)現(xiàn)，采用編寫的分形維數(shù)計算程序?qū)?次迭代的Koch雪花曲線（圖2）進行處理，處理流程如圖3所示，對周長和面積進行雙對數(shù)擬合，得到Koch雪花曲線的周長-面積雙對數(shù)擬合線（圖4），計算得到分形維數(shù)為1.267 3，相關(guān)系數(shù)達到0.999 9，與Koch雪花曲線理論的分形維數(shù)1.261 3十分接近，故采用該程序計算預混料常用有機載體的分形維數(shù)具有可行性。

采用上述分形維數(shù)計算程序?qū)?種預混料常用有機載體的周長-面積分形維數(shù)進行測定，預混料載體顆粒圖像處理流程以全脂米糠粉顆粒為例，該米糠粉顆粒不同放大倍數(shù)圖像如圖5所示，100倍放大倍數(shù)下圖像處理流程如圖6所示，提取顆粒在不同放大倍數(shù)下填充之后得到的面積和邊界，計算周長并進行雙對數(shù)擬合，得到全脂米糠粉顆粒的周長-面積雙對數(shù)擬合線如圖7所示。

預混料常用有機載體周長-面積分形維數(shù)的測定結(jié)果如表1所示，各樣品顆粒采用最小二乘法擬合的相關(guān)系數(shù)均達到0.99以上，表明這5種預混料常用載體具有統(tǒng)計自相似結(jié)構(gòu)，周長-面積法是測定這類物料分形維數(shù)的可行方法。

根據(jù)分形理論可知，分形維數(shù)的大小反映了結(jié)構(gòu)復雜程度的大小，周長-面積分形維數(shù)的測定結(jié)果顯示，不同載體的分形維數(shù)均值不同，同一載體不同顆粒的分形維數(shù)亦不同，表明不同載體顆粒微觀復雜程度不同，同種載體顆粒的微觀復雜程度也有一定的差異性，各載體分形維數(shù)均值呈現(xiàn)的關(guān)系為玉米粉>膨化大豆粉>全脂米糠粉>全麥粉>豆粕粉。

3討論

玉米粉的分形維數(shù)比全脂米糠粉的分形維數(shù)略高，這與廖娜[8]采用豪斯道夫維數(shù)和差分計盒維數(shù)測定的結(jié)果相一致，表明周長-面積法計算的分形維數(shù)和其他邊界分形維數(shù)在表征分形體邊界復雜程度上具有一定的聯(lián)系，周長-面積分形維數(shù)描述了預混料有機載體顆粒微觀形貌的復雜程度，可以作為表征預混料有機載體的一個參量，從而進一步研究預混料載體的分形維數(shù)與其流動性、承載能力等的關(guān)系。

自然界的自相似圖形都是統(tǒng)計自相似，同種物料的不同顆粒形狀亦千差萬別，故為了得到某種物料較準確的分形維數(shù)，在條件允許的情況下應盡可能多的采集樣本，取樣本顆粒分形維數(shù)的均值，或者用一個區(qū)間來表示某種物料的分形維數(shù)，樣本量偏少可能會降低分形維數(shù)的準確度。

參考文獻：

[1] 孫博玲.分形維數(shù)（Fractal dimension）及其測量方法[J].東北林業(yè)大學學報，2004，32（3）：116-119.

[2] 朱華，姬翠翠.分形理論及其應用[M]. 北京：科學出版社，2011.22-53.

[3] 溫芝元.基于數(shù)字圖像的水果分形維數(shù)特征[J].湖南農(nóng)業(yè)大學學報（自然科學版），2010，36（4）：478-482.

[4] 曹樂平.基于周長面積分形維數(shù)的柑橘品種機器識別[J].農(nóng)業(yè)工程學報，2010，26（2）：351-355.

[5] 李金萍，蓋國勝，鄭洲順，等.粉體顆粒的分形表征與分維計算[J]. 有色礦冶，2005，21（增刊）：92-94.

[6] 譚佐軍，盧軍，楊長舉，等.儲藏物昆蟲分形特性的研究[J].中國糧油學報，2008，23（2）：163-165.

[7] 龔國淑，張世熔.植物病害病斑形狀的分形研究[J].植物保護，2002，28（6）：9-13.

[8] 廖娜.有機粉狀載體表面分形維數(shù)與其承載能力關(guān)系研究[D].武漢：華中農(nóng)業(yè)大學，2008.

[9] 張龍.有機粉體分形表征與其流動性關(guān)系的研究[D].武漢：華中農(nóng)業(yè)大學，2010.

[10] MANDELBROT B B. Les objets fractals：forme，hasard et dimension[M]. Paris：Flammarion，1975.

[11] 徐永春，趙喜清.Koch雪花曲線的制作及其重要結(jié)論[J].長春師范學院學報，2003，22（1）：6-8.