韓書霞 戚大偉

(東北林業(yè)大學(xué)，哈爾濱，150040)

采用計(jì)算機(jī)斷層掃描技術(shù)(CT)對(duì)原木進(jìn)行無(wú)損檢測(cè)，對(duì)取得的缺陷圖像采用分形理論進(jìn)行圖像處理，目前國(guó)內(nèi)外已有很多學(xué)者在研究［1］，這些研究從分形理論的分形維數(shù)方面已取得了非?？上驳某晒５珜?duì)多重分形頻譜的選擇卻各有不同。本文通過采用多重分形三種頻譜fh(α)、fg(α)、fl(α)的選擇，對(duì)原木缺陷CT圖像進(jìn)行了處理，找到了一種新的準(zhǔn)確確定原木內(nèi)部缺陷邊緣的方法。

1 多重分形分析理論與圖像分析

1.1 多重分形理論簡(jiǎn)介

多重分形也稱為標(biāo)度分形或復(fù)分形。常用來表示一個(gè)整體特征標(biāo)度指數(shù)(分形維數(shù))所不能完全描述的奇異性幾率分布，它是從信號(hào)的局部特征出發(fā)來研究其特征。多重分形分析的目的在于量化測(cè)度的奇異結(jié)構(gòu)，以及在尺度發(fā)生變化時(shí)為伴隨有不同范圍的冪定律的現(xiàn)象提供模型［2］，尤其適用于一些難以建模的不規(guī)則圖像處理和分析。由于多重分形分析具有良好的局部性和全局性，因此為圖像的邊緣檢測(cè)提供了一種新的途徑。

1.2 H?lder指數(shù)的應(yīng)用

設(shè)f(x)為f:R→R的函數(shù)，如果有∣f(x)－P(x－x0)∣＜C∣x－x0∣α，其中:α＞0，P為不超過 α 的多項(xiàng)式，x∈B(x0，δ)得閉球，C為一大于0的常數(shù)，令

α(x0)=sup(α:f∈Cαx0);

α(x)=sup(α:f∈Cαx)。

則稱α(x0)為函數(shù)f(x)在x0點(diǎn)的 H?lder指數(shù)，α(x)為f(x)在其定義域上的H?lder指數(shù)［3］，α是表征分形體某小區(qū)域的分維，稱為局部分維。因此H?lder指數(shù)α能夠很好地描述函數(shù)f(x)的局部特性。

1.3 多重分形頻譜與圖像分析

1.3.1 Hausdorff頻譜

多重分形用H?lder指數(shù)α表示分形體小區(qū)域的分維，由于小區(qū)域的數(shù)目很大，于是可得一個(gè)由不同α所組成的無(wú)窮序列構(gòu)成的譜f(α)，f(α)又稱為奇異譜。

設(shè)Eα為 H?lder指數(shù) α的點(diǎn)集，如果令fh(α)=dimh(Eα)，其中 dimh為的 Hausdorff維，那么稱fh(α)為Eα的Hausdorff頻譜，因此Hausdorff頻譜能夠很好地描述點(diǎn)集Eα的全局。

多重分形圖像分析過程，首先選擇一容量序列vn及參考測(cè)度μ［4］，計(jì)算出圖像的H?lder指數(shù)α(x)和多重分形頻譜fh(α(x));然后根據(jù)(α(x)，fh(α(x)))的值用幾何和概率的方法對(duì)圖像的每個(gè)像素點(diǎn)進(jìn)行分類和處理。因此，判斷某一像素點(diǎn)是否為邊緣，既取決于局部特征，又要考慮到圖像的整體性。

1.3.2 粗?；l譜和離差頻譜

當(dāng)n∈N，ε＞0，令

N(α)=，其中#()表示集合中元素的數(shù)目設(shè)為α的粗?；l譜和離差頻譜。

1.3.3 Legendre 頻譜

則存在一個(gè)凹函數(shù)τ，使得

設(shè)τ*為τ的Legendre多重分形頻譜，則τ的Legendre變換為fl(α)= τ*(α)=infq［qα－τ(q)］，稱fl(α)為 α 的 Legendre頻譜［5］。

1.4 多重分形頻譜的選擇

多重分形理論的中心問題是三種頻譜fh(α)、fg(α)、fl(α)的選擇，由以上分析可以看出fl(α)的計(jì)算最為簡(jiǎn)單，但會(huì)導(dǎo)致信息的丟失;fg(α)要逐點(diǎn)計(jì)算，使奇異點(diǎn)分布的概率描述，能夠很好的反映出頻譜的細(xì)節(jié)信息;fh(α)對(duì)應(yīng)于Hausdorff維，是奇異點(diǎn)的幾何描述，計(jì)算比較復(fù)雜。在某些特定的情況下fh(α)=fg(α)=fl(α)，定為f(α)，但三者之間一般滿足fh(α)≥fg(α)≥fl(α)的關(guān)系［3－4］。

根據(jù)多重分形理論及上述分析可得到:①一般說來，粗糙區(qū)域的點(diǎn)的規(guī)律性較低，需要找到可以與輪廓線區(qū)分的方法，如果與某點(diǎn)的H?lder指數(shù)相應(yīng)的多重分形頻譜的數(shù)值為f(α)=1，則該點(diǎn)一定在輪廓上。②邊緣上的點(diǎn)除了有幾何學(xué)的特征外，還存在統(tǒng)計(jì)學(xué)特征:當(dāng)一個(gè)像素點(diǎn)是一個(gè)分辨率的圖像中任意選定的，通過它被找到的概率可以判斷其是否為邊緣上的點(diǎn)。例:若頻譜f(α)=1時(shí)，可以篩選出邊緣點(diǎn);若頻譜f(α)=1.5時(shí)，可以得到不規(guī)則的輪廓線;若頻譜f(α)接近2時(shí)，可以獲得光滑區(qū)域或粗糙區(qū)域的信息［6］。

1.5 圖像的多重分形邊緣檢測(cè)方法

根據(jù)所得出的多重分形頻譜，可以將圖像的像素點(diǎn)進(jìn)行分類，若像素點(diǎn)的多重分形頻譜f(α)=1.0，則此像素點(diǎn)為平滑邊緣點(diǎn);若f(α)=t(1.0＜t＜2.0)，則稱該像素點(diǎn)為t類邊緣平滑點(diǎn)。平滑邊緣點(diǎn)和奇異邊緣點(diǎn)就構(gòu)成所檢測(cè)圖像的邊緣。

2 方法驗(yàn)證

用計(jì)算機(jī)斷層掃描技術(shù)(CT)對(duì)木材進(jìn)行缺陷檢測(cè)，實(shí)驗(yàn)所用的圖像為原木CT腐朽圖像(見圖1)。圖2～圖4分別為多重分形頻譜f(α)=1.0、1.0≤f(α)≤1.1、1.2≤f(α)≤1.3邊緣檢測(cè)結(jié)果。從圖2～圖4可以看出，只要選擇適當(dāng)?shù)亩嘀胤中晤l譜閾值f(α)的值，使用多重分形頻譜分析的邊緣檢測(cè)方法能夠測(cè)出比較精細(xì)的邊緣，且閾值f(α)的范圍增大，檢測(cè)出的邊緣就更精確。

圖1 原木CT腐朽原圖像

圖2 f(α)=1.0邊緣檢測(cè)結(jié)果

3 結(jié)束語(yǔ)

通過實(shí)驗(yàn)證實(shí)，基于多重分形頻譜分析進(jìn)行的木材CT圖像的邊緣檢測(cè)是有效的，優(yōu)點(diǎn)是:該方法首先計(jì)算圖像各像素點(diǎn)的H?lder指數(shù)α(x，y)，然后估計(jì)出其多重分形頻譜f(α)的值，最后對(duì)圖像的像素點(diǎn)進(jìn)行分類;同時(shí)從實(shí)驗(yàn)的結(jié)果可以看出，多重分形頻譜的選擇直接影響邊緣檢測(cè)的結(jié)果，選擇適當(dāng)?shù)念l譜閾值，可以得到更多的邊緣細(xì)節(jié);該方法對(duì)于圖像邊緣檢測(cè)也具有良好的局部性。因此，多重分形頻譜分析為木材CT圖像邊緣檢測(cè)提供了一種新方法。

圖3 1.0≤f(α)≤1.1 邊緣檢測(cè)結(jié)果

圖4 1.2≤f(α)≤1.3 邊緣檢測(cè)結(jié)果

［1］Qi Dawei，Yu Lei，Han Shuxia，et al.Based on computed tomography multifractal analysis of wood defect［C］//IEEE International Conference on Control and Automation.Guangzhou，China:IEEE International Conference on Control and Automation，May 30 to June 1，2007:336－341.

［2］肯尼思·法爾科內(nèi).分形幾何中的技巧［M］.曾文曲，王向陽(yáng)，陸夷，譯.沈陽(yáng):東北大學(xué)出版社，1999.

［3］王耀南，王紹源.基于多重分形頻譜分析的圖像邊緣檢測(cè)［J］.湖南大學(xué)學(xué)報(bào):自然科學(xué)版，2004，31(4):28－33.

［4］Lévy Véhel J，VOJAKR.Multifractal analysis of choquet capacities:preliminary result［J］.Adv Appl Math，1998，20(1):1－43.

［5］戚大偉.基于分形理論的原木缺陷圖像處理［M］.哈爾濱:哈爾濱地圖出版社，2004.

［6］Canus C.Robust large deviation multifractal spectrum estimation［M］.Tangier:Proceedings of International Wavelets Conference，1998:739－740.