郭樹旭 趙 靜 李雪妍

(吉林大學(xué)電子科學(xué)與工程學(xué)院 長春 130012)

基于中心-輪廓距離特征統(tǒng)計(jì)的形狀表示方法

郭樹旭 趙 靜 李雪妍*

(吉林大學(xué)電子科學(xué)與工程學(xué)院 長春 130012)

該文提出一種新的基于特征統(tǒng)計(jì)的形狀描述方法。通過對(duì)中心-輪廓距離(CCD)和傳統(tǒng)鏈碼(Chaincode)的聯(lián)合統(tǒng)計(jì)分析，使用中心-輪廓距離對(duì)形狀進(jìn)行層次分解，對(duì)各層的形狀映射部分的鏈碼描述進(jìn)行統(tǒng)計(jì)分析，從而形成中心-輪廓距離和鏈碼的聯(lián)合統(tǒng)計(jì)(JSCCDC)描述子。形狀之間的相似性可以用JSCCDC的城區(qū)距離來描述。實(shí)驗(yàn)結(jié)果表明，該表示方法兼具了形狀的全局特征和局部特征，相比于傳統(tǒng)的特征加權(quán)方法具有更優(yōu)越的性能，在形狀匹配和形狀檢索中具有較高的精度和可靠性。

模式識(shí)別；形狀表示；特征統(tǒng)計(jì)；鏈碼；中心-輪廓距離；匹配矩陣

1 引言

形狀匹配與分類是模式識(shí)別與計(jì)算機(jī)視覺研究的重要問題，被廣泛應(yīng)用在很多領(lǐng)域，如目標(biāo)識(shí)別、圖像檢索、人臉識(shí)別、醫(yī)學(xué)圖像診斷等，主要包括形狀表示、形狀匹配和度量學(xué)習(xí)3個(gè)模塊。對(duì)于一個(gè)二值形狀，先要提取其輪廓特征描述子，其區(qū)分能力的強(qiáng)弱將直接影響形狀識(shí)別的結(jié)果，大多數(shù)形狀匹配的工作都集中在此。匹配過程是找到一對(duì)不同形狀之間的整體及局部的對(duì)應(yīng)關(guān)系，對(duì)應(yīng)關(guān)系的準(zhǔn)確性也將會(huì)對(duì)之后得到的非相似度的區(qū)分能力產(chǎn)生直接影響。度量學(xué)習(xí)則是通過已知的數(shù)據(jù)庫形狀的上下文信息來改進(jìn)原有距離度量的方法，這一步驟能將原有距離的區(qū)分性能大幅度提升[1]。

形狀表示，又稱特征提取或形狀描述，是通過某種方法生成一個(gè)數(shù)值化的描述子來刻畫形狀特征的過程。輪廓特征描述子區(qū)分能力的強(qiáng)弱將直接影響形狀識(shí)別與分類的結(jié)果，它是有效完成形狀匹配任務(wù)的關(guān)鍵所在。形狀特征描述子主要可以分為基于輪廓和基于區(qū)域兩大類?；趨^(qū)域的方法是利用物體內(nèi)部區(qū)域(所有像素點(diǎn))的信息來表示形狀，而基于輪廓的方法主要是利用物體的邊界輪廓信息來表示形狀，與前者相比，其優(yōu)勢(shì)主要體現(xiàn)在對(duì)圖像低層特征的高識(shí)別度以及相對(duì)較小的計(jì)算量上。因此，該方法成為近年來形狀表示研究的主流。

基于輪廓的形狀表示方法大致可以分為4類：全局描述子，局部描述子，多尺度描述子和多方面描述子。早期描述子，如邊界長度、直徑、圓度等都屬于全局描述子。全局特征一般是平移，旋轉(zhuǎn)不變的，計(jì)算簡單，但僅僅只能表示形狀的大致特性，缺少細(xì)節(jié)描述，區(qū)分力不足。為此部分學(xué)者提出了局部形狀描述子的概念，如鏈碼和形狀顯著性描述符(Shape Salience Descriptor, SSD)[2]，此類描述子對(duì)形狀進(jìn)行了細(xì)致刻畫，但對(duì)噪聲非常敏感。針對(duì)這一問題，出現(xiàn)了多尺度描述子，如文獻(xiàn)[3]提出的多尺度分形(Multi-Scale Fractal, MS Fractal)維數(shù)算法和文獻(xiàn)[4]提出的用輪廓點(diǎn)控制尺度的算法。此外還有以文獻(xiàn)[5]提出的以形狀上下文為代表的多方面描述子，通過輪廓點(diǎn)的空間位置關(guān)系和分布來反映形狀特征。近年，還有部分學(xué)者提出了特征統(tǒng)計(jì)的方法，如同心離散圓簇描述法[6]。

單一特征描述子無法勝任復(fù)雜可變形狀的識(shí)別工作。如鏈碼是一種用曲線起始點(diǎn)的坐標(biāo)和邊界點(diǎn)方向代碼來描述曲線或邊界的方法，但其編碼只依賴于輪廓序列中相鄰點(diǎn)的相對(duì)位置信息，無法體現(xiàn)輪廓的全局特征，且碼串長，在傳輸?shù)倪^程中易受干擾。中心距離函數(shù)[7]是將輪廓線上的點(diǎn)到形狀的幾何中心的距離描述成中心角度的函數(shù)，這種方法一般能夠重構(gòu)被描述的形狀，但當(dāng)形狀的幾何中心位于形狀區(qū)域之外或者被描述形狀過于復(fù)雜時(shí)，就可能出現(xiàn)一個(gè)中心角對(duì)應(yīng)多中心距的情況，為了能將其轉(zhuǎn)化成1維函數(shù)，通常對(duì)這些多值有所取舍或者求取平均值，但無論怎樣取值，都造成了形狀信息的丟失從而不能重構(gòu)，無法滿足形狀描述的唯一性[8]。中心輪廓距離曲線(Centroid Contour Distance Curve, CCDC)[9]雖然克服了中心距離函數(shù)一個(gè)中心角對(duì)應(yīng)多個(gè)中心距的問題，但由于數(shù)字圖像量化定義的弧長精度，CCDC對(duì)形狀輪廓各部分的描述精度及所占帶寬隨著其像素的增加而遞增，因此抗噪性弱，且無法重構(gòu)圖像。

針對(duì)上述問題，本文從特征統(tǒng)計(jì)方法入手，融合了鏈碼在局部描述方面的優(yōu)勢(shì)和中心-輪廓距離(Centroid-Contour Distance, CCD)對(duì)全局特征刻畫的優(yōu)勢(shì)，提出基于中心-輪廓距離和鏈碼的聯(lián)合統(tǒng)計(jì)(Joint Statistical of Centroid-Contour Distance and Chaincode, JSCCDC)描述形狀的新方法：通過使用CCD對(duì)形狀進(jìn)行層次分解，對(duì)映射到每一層的形狀的鏈碼描述進(jìn)行統(tǒng)計(jì)分析，從而得到形狀的JSCCDC特征。實(shí)驗(yàn)證明JSCCDC能夠很好地表達(dá)形狀的全局與局部特征，相較于傳統(tǒng)的形狀表示方法更有效。

2 鏈碼的統(tǒng)計(jì)特征

2.1 統(tǒng)計(jì)鏈碼

由于鏈碼在形狀匹配中無法克服旋轉(zhuǎn)、噪聲和目標(biāo)尺度變化帶來的一系列問題，文獻(xiàn)[10]提出一種基于Freeman鏈碼的形狀描述方法：最小和統(tǒng)計(jì)方向碼(Minimum Summation Statistical Direction Code, MSSDC)。這種方法雖然在一定程度上實(shí)現(xiàn)了對(duì)形狀的固定或者相似角度的觀察，但由于噪聲或者柔性形變引起的變化，在進(jìn)行形狀匹配時(shí)，并不能保證以兩者最相似的視角進(jìn)行匹配。為此，本文將在2.2節(jié)提出匹配矩陣的概念。

基于文獻(xiàn)[10]，對(duì)于給定的形狀鏈碼表示，本文給出N-方向的歸一化的統(tǒng)計(jì)鏈碼S形式：

其中，fk表示碼元k出現(xiàn)的頻率，根據(jù)定義可知，S滿足歸一化條件。

2.2 匹配矩陣

為了實(shí)現(xiàn)統(tǒng)計(jì)鏈碼的最佳匹配，本文定義了N階匹配矩陣M來匹配N-方向統(tǒng)計(jì)鏈碼：

匹配矩陣M的每行每列分別對(duì)應(yīng)著形狀旋轉(zhuǎn)和翻轉(zhuǎn)一定角度后的形態(tài)，角度精度由鏈碼的方向數(shù)N決定。應(yīng)用匹配矩陣的統(tǒng)計(jì)鏈碼匹配流程為：

(1)首先，給出兩個(gè)待匹配形狀A(yù)和B的統(tǒng)計(jì)鏈碼形式SA和SB；

(2)其次，計(jì)算基于SB的匹配矩陣：

(3)再次，根據(jù)式(3)計(jì)算SA在匹配矩陣上的投影系數(shù)矩陣Q。

找到投影系數(shù)矩陣Q的最大值并標(biāo)記其所在的行(第a行)和列(第b列)，如式(4)所示。

在式(3)中，Q是本文定義的基于匹配矩陣M的投影系數(shù)矩陣，以SA和MB構(gòu)造的投影矩陣為例，Q中的每個(gè)元素對(duì)應(yīng)的是形狀A(yù)與不同角度觀察到的形狀B的相似程度。

(4)根據(jù)最小二乘原理，本文認(rèn)為使投影系數(shù)矩陣Q取得最大值的MB的行或者列是與待匹配形狀A(yù)最相近的視角觀察到的形態(tài)的統(tǒng)計(jì)鏈碼表示，記為

式(5)中，a和b是式(4)中標(biāo)記的Q的最大值的行列位置，WN(b)定義如式(6)所示。

(5)最后，使用L2范數(shù)計(jì)算SA和SB/A的距離，如式(7)所示，Ds的值越大，則表示形狀A(yù)和B的相似度越低，顯然，Ds(A,A)=0。

3 中心-輪廓距離(CCD)的特征統(tǒng)計(jì)

3.1 基于CCD的傳統(tǒng)描述

傳統(tǒng)的基于CCD的形狀描述方法主要有中心距離函數(shù)和中心輪廓距離曲線(CCDC)。中心距離函數(shù)將輪廓上的點(diǎn)到形狀幾何中心的距離描述成中心角度的函數(shù)，量化時(shí)會(huì)舍去部分點(diǎn)的信息；CCDC用序列號(hào)取代中心角的方案解決了輪廓點(diǎn)取舍的問題，卻也相當(dāng)于根據(jù)輪廓的周長對(duì)其進(jìn)行加權(quán)，這與人們平時(shí)的視覺處理方式不一致，降低了抗噪性，給匹配帶來了困難。

3.2 極半徑——統(tǒng)計(jì)CCD

由于CCD的尺度不變特征，對(duì)其進(jìn)行角度化處理可能損失部分點(diǎn)的信息，如中心距離函數(shù)，而CCDC卻由于序列的等間隔性天然地根據(jù)周長對(duì)CCD進(jìn)行了加權(quán)而引入了匹配過程難以復(fù)原的尺度信息，需要對(duì)CCDC做分段切割處理，也相當(dāng)于是一種角度化處理的過程。為了克服上述針對(duì)CCD角度化引起的問題，本文提出幅度化的概念，即對(duì)CCD的幅值進(jìn)行統(tǒng)計(jì)分類的過程。

(1)本文給出形狀幾何中心的定義：

(io, jo)是目標(biāo)形狀輪廓的幾何中心，圖像大小是m×n, f(i,j)是邊界坐標(biāo)(i,j)灰度值，對(duì)于二值圖象，f(i,j)取值只可能是“0”或者“1”。

(2)首先，順時(shí)針找出目標(biāo)形狀的邊界坐標(biāo)，記做C={(i,j)|f(i,j)≠0}={(i1,j1),(i2,j2),…}。

(3)根據(jù)(io, jo)和C給出形狀的CCD串，為了區(qū)別于是其他基于CCD的描述子，本文定義使用的CCD為極半徑R：

(4)其次，把統(tǒng)計(jì)的思想應(yīng)用到R中，確定精度L后，給出分段步長：

(5)根據(jù)t給出統(tǒng)計(jì)極半徑的定義：

式(11)中，fr1+k?t,r2+k?t是極半徑R落在[r1+k·t,r2+k·t]區(qū)間的的頻率，即把形狀質(zhì)心映射到一個(gè)內(nèi)徑為r1+k·t，外徑為r2+k·t的圓環(huán)中心，統(tǒng)計(jì)圓環(huán)內(nèi)目標(biāo)形狀的弧長總和，再利用形狀的周長進(jìn)行歸一化得到極半徑在當(dāng)前區(qū)間的出現(xiàn)頻率。通過精度確定的步長t來定義圓環(huán)寬度，遍歷整個(gè)極半徑區(qū)間則得到定義的統(tǒng)計(jì)極半徑。

(6)最后，同樣采用L2范數(shù)來定義不同形狀的統(tǒng)計(jì)極半徑的距離：

4 CCD和鏈碼的聯(lián)合統(tǒng)計(jì)

針對(duì)統(tǒng)計(jì)鏈碼和統(tǒng)計(jì)極坐半徑，為了結(jié)合它們?cè)谛螤蠲枋龇矫娓髯缘膬?yōu)勢(shì)，本文提出了聯(lián)合統(tǒng)計(jì)。

4.1 鏈碼和極半徑的聯(lián)合統(tǒng)計(jì)描述

首先，本文在Freeman鏈碼的基礎(chǔ)之上進(jìn)行改進(jìn)，在原始鏈碼中加入了極半徑信息，獲得鏈碼-極半徑的聯(lián)合描述子，定義如下：

其中，F(xiàn)是對(duì)形狀A(yù)的輪廓進(jìn)行Freeman編碼的結(jié)果，R是對(duì)應(yīng)的編碼點(diǎn)的極半徑，可根據(jù)式(9)求得。

定義鏈碼和極半徑的聯(lián)合統(tǒng)計(jì)描述子(JSCCDC), JSCCDC在統(tǒng)計(jì)鏈碼的同時(shí)也統(tǒng)計(jì)極半徑空間位置，其形式為

元素RiFj表示的是在第i個(gè)極半徑空間中坐落的j方向的矢量線段的概率。N是N方向Freeman鏈碼，L是對(duì)極半徑區(qū)間的L等分。如圖1(a)所示的蝙蝠輪廓，首先設(shè)定參數(shù)L=8, N=8，然后根據(jù)式(8)標(biāo)記其幾何中心，并根據(jù)式(9)和值確定分段步長t，畫出等分圓環(huán)。最后統(tǒng)計(jì)落在每一個(gè)圓環(huán)內(nèi)的輪廓的鏈碼分布，并對(duì)整個(gè)輪廓進(jìn)行歸一化，即可得到式(14)中的PC矩陣，也就是形狀的JSCCDC描述，如圖1(b)所示。

4.2 JSCCDC匹配

圖1 JSCCDC描述子

JSCCDC是一種矩陣描述。本文借鑒圖像處理中像素間的城區(qū)距離(即模為1的距離)來描述兩個(gè)矩陣的相似性，即利用式(15)分別計(jì)算兩個(gè)矩陣對(duì)應(yīng)位置的差，取其絕對(duì)值和為兩個(gè)矩陣的距離。同時(shí)，也借鑒了文獻(xiàn)[8]中距離測(cè)量公式，即分別對(duì)矩陣的每行以L2范數(shù)來衡量其差異度，最后把各行差異度的平方和作為矩陣差異度，定義如式(16)。在5.1節(jié)的相似性評(píng)測(cè)實(shí)驗(yàn)中，為了和文獻(xiàn)[8]的方法比對(duì)，本文采用與其相同的距離測(cè)量，即式(16)，但時(shí)間復(fù)雜度較高，在檢索實(shí)驗(yàn)中使用城區(qū)距離來計(jì)算。

為了實(shí)現(xiàn)匹配的魯棒性，對(duì)矩陣PC的行向量引入式(2)的匹配矩陣，以實(shí)現(xiàn)JSCCDC描述子對(duì)平移、旋轉(zhuǎn)、翻轉(zhuǎn)和尺度變換的魯棒描述。

5 實(shí)驗(yàn)

5.1 相似性實(shí)驗(yàn)

本文選取了MPEG7 CE-Shape-1形狀數(shù)據(jù)庫中的Hammer圖像，對(duì)其進(jìn)行了如圖2所示的各種變換，并隨機(jī)從MPEG7 CE-Shape-1選取了圖3所示的其他形狀作為類間形狀對(duì)比，給出了不同方法下形狀的類間類內(nèi)距離測(cè)度。實(shí)驗(yàn)參數(shù)：N=8, L=8。

表1給出了使用JSCCDC結(jié)合匹配矩陣方法測(cè)得的圖2和圖3中圖例的類內(nèi)類間距離。通過類內(nèi)距離來看JSCCDC的表征能力，不僅能夠?qū)崿F(xiàn)對(duì)平移、旋轉(zhuǎn)和尺度變換的魯棒性，由于匹配矩陣的引入，還能還原形狀的翻轉(zhuǎn)變換。但由于數(shù)字圖像的量化問題，對(duì)形狀的旋轉(zhuǎn)和尺度變換還原還有一定的誤差。

結(jié)合類間距離再來看JSCCDC的分類能力：JSCCDC測(cè)得的Hammer形狀的類內(nèi)距離落在[0,0.0024]，狹義類間距離(類內(nèi)測(cè)試形狀與其他形狀的距離)落在[0.0255, 0.0376]，廣義類間距離(類間形狀的相互距離)落在[0.0092, 0.0408]，與類內(nèi)距離無交集，可以實(shí)現(xiàn)形狀分類。本文將分類能力定義為廣義類間距離與類內(nèi)距離的比值，由此得出JSCCDC對(duì)Hammer形狀的分類能力為47.6。

圖2 Hammer形狀的10種變形

圖3 節(jié)選自MPEG7的類間形狀示例

表1 JSCCDC測(cè)得的類內(nèi)類間距離

表2 不同算法相似度評(píng)測(cè)對(duì)比

表2是其他算法在相似度實(shí)驗(yàn)中和本文方法的各類指標(biāo)比較?？梢钥闯?，JSCCDC在相似性評(píng)測(cè)方面的性能是優(yōu)于很多傳統(tǒng)算法的。為了測(cè)試實(shí)驗(yàn)的普遍性，本文對(duì)圖3中的其他形狀均進(jìn)行了圖2所示的各種變換，代替Hammer形狀重復(fù)上述相似性實(shí)驗(yàn)過程，得到了一個(gè)均值統(tǒng)計(jì)，表2最后一列所示。

5.2 形狀檢索實(shí)驗(yàn)

5.2.1 MPEG-7數(shù)據(jù)庫 MPEG-7 CE-Shape-1數(shù)據(jù)庫是由文獻(xiàn)[11]在2000年發(fā)布的，共包含3個(gè)部分，其第2部分主要用來衡量基于相似性方法的檢索精度，是目前衡量形狀描述子的重要參考指標(biāo)之一。它包含70個(gè)類別、每類20幅不同形態(tài)共1400張形狀圖像，如圖4(a)和圖4(b)所示。使用Bull-eye方法度量檢索精度，對(duì)每一幅圖像檢索出最相似的40幅圖像，檢索率R可通過式(17)計(jì)算，其中，T(j)表示檢索出的與第j幅待檢形狀同一類形狀的個(gè)數(shù)。

表3是本文算法和部分經(jīng)典算法在MPEG-7數(shù)據(jù)庫上的檢索率。其中，F(xiàn)ourier是經(jīng)典的傅里葉描述子[12]，CS是一種基于輪廓顯著(Contour Salience, CS)的描述[13]，MSSDC是2.1節(jié)提到的最小統(tǒng)計(jì)和方向碼，MI是一種基于不變矩(Moment Invariant, MI )的描述[14]，MS Fractal是一種多尺度分形描述，SSD是一種基于局部距離函數(shù)的形狀顯著性描述子，SSD+GF是加入全局特征(Global Feature, GF)的SSD描述子[2]，JSCCDC+M是應(yīng)用匹配矩陣的JSCCDC描述。對(duì)比3和5以及9和10兩組數(shù)據(jù)，本文提出的匹配矩陣在統(tǒng)計(jì)鏈碼的匹配中實(shí)現(xiàn)了有效的對(duì)齊，檢索率都有所提高。雖然相較于加入全局特征的SSD描述子，JSCCDC在檢索率表現(xiàn)方面略遜一籌，但對(duì)比5和10以及8和11兩組數(shù)據(jù)，本文在統(tǒng)計(jì)鏈碼基礎(chǔ)之上加入全局特征而提出的JSCCDC描述子提高了統(tǒng)計(jì)鏈碼的檢索率絕對(duì)值達(dá)0.28，高達(dá)70%，而加入全局特征的SSD描述子在原有基礎(chǔ)之上提高了0.10，只有16%，可以看出本文提出的聯(lián)合統(tǒng)計(jì)的方法較傳統(tǒng)的特征引入法更能改善描述子性能，因此，JSCCDC這種多特征聯(lián)合統(tǒng)計(jì)的描述方法在形狀檢索工作中更具優(yōu)勢(shì)。

5.2.2 Kimia99數(shù)據(jù)庫 Kimia99數(shù)據(jù)庫由文獻(xiàn)[16]提出，共包含9類形狀，每一類11個(gè)子形狀，如圖4(c)所示。這些形狀中不僅同類形狀有局部的變形，還有局部的遮擋，且不同形狀甚至還有相似的全局特征。

Kimia99數(shù)據(jù)庫通用的評(píng)價(jià)標(biāo)準(zhǔn)通常以表格的形式呈現(xiàn)，數(shù)據(jù)庫中的每一個(gè)圖像都要作為模板對(duì)全庫進(jìn)行檢索。表4是本文方法和部分其他經(jīng)典算法的對(duì)比，其中SP加權(quán)法是本文提出的統(tǒng)計(jì)鏈碼和統(tǒng)計(jì)極半徑平均加權(quán)的方法?？梢钥闯鯦SCCDC的正確識(shí)別率不僅優(yōu)于其他經(jīng)典的基于單一特征的形狀描述子，而且效果優(yōu)于特征加權(quán)方法。

圖4 本文引用數(shù)據(jù)庫圖例

表3 不同算法在MPEG-7數(shù)據(jù)庫上的檢索率

5.3 時(shí)效性分析

實(shí)驗(yàn)條件：硬件為AMD Athlom(tm)ⅡX2 250 Processor, RAM= 4 GB；軟件為Win7, 64位操作系統(tǒng)，MATLAB；數(shù)據(jù)庫為MPEG-7(1400幅，gif格式)，Kimia99(99幅，pgm格式)。

5.1 節(jié)的相似性實(shí)驗(yàn)和5.2節(jié)的檢索實(shí)驗(yàn)證明了JSCCDC這種聯(lián)合統(tǒng)計(jì)描述子在形狀表示方面的優(yōu)越性，為了進(jìn)一步驗(yàn)證該算法是否具有實(shí)際應(yīng)用價(jià)值，特別設(shè)計(jì)了時(shí)效性實(shí)驗(yàn)，主要對(duì)比基于鏈碼的各種描述子在MPEG-7數(shù)據(jù)庫和Kimia99數(shù)據(jù)庫上的檢索時(shí)間和性能，為了體現(xiàn)基于輪廓描述法在處理時(shí)間上的優(yōu)越性，特別加入了基于區(qū)域的10階Zernike矩算法進(jìn)行對(duì)比。實(shí)驗(yàn)結(jié)果如表5所示。JSCCDC處理速度稍遜統(tǒng)計(jì)鏈碼，仍優(yōu)于MSSDC算法，但檢索性能卻較兩者有大幅度提高。

表4 部分算法在Kimia99數(shù)據(jù)庫上的檢索效果

6 結(jié)束語

本文提出了一種新的基于形狀中心-輪廓距離和鏈碼聯(lián)合統(tǒng)計(jì)的形狀描述和匹配方法。首先，它通過中心-輪廓距離對(duì)形狀進(jìn)行層次分解，然后統(tǒng)計(jì)每一層形狀映射部分的鏈碼描述，從而形成了中心-輪廓距離和鏈碼的聯(lián)合統(tǒng)計(jì)(JSCCDC)。這種描述方法通過對(duì)形狀的兩種簡單特征的聯(lián)合統(tǒng)計(jì)，既描述了形狀的局部特征，又包含了全局信息。經(jīng)典數(shù)據(jù)庫上測(cè)試結(jié)果表明，這種聯(lián)合統(tǒng)計(jì)的方法不僅優(yōu)于單一特征描述法，而且優(yōu)于加權(quán)方式結(jié)合的局部-全局特征描述。

表5 部分鏈碼算法在Kimia99和MPEG-7數(shù)據(jù)庫上性能對(duì)比

[1] 周瑜, 劉俊濤, 白翔. 形狀匹配方法研究與展望[J]. 自動(dòng)化學(xué)報(bào), 2012, 38(6): 889-910. Zhou Yu, Liu Jun-tao, and Bai Xiang. Research and perspective on shape matching[J]. Acta Automatica Sinica, 2012, 38(6): 889-910.

[2] Glauco V P, Marcos A B, and Celia A Z B. Image feature descriptor based on shape salience points[J]. Neurocomputing, 2013, 120(23): 156-163.

[3] Torres R S, Falc?o A X, and Costa L F. A graph-based approach for multiscale shape analysis[J]. Pattern Recognition, 2004, 37(6): 1163-1174.

[4] Alajlan N, Kamel M S, and Freeman G H. Geometry-based image retrieval in binary image databases[J]. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2008, 30(6): 1003-1013.

[5] Belongie S, Malik J, and Puzicha J. Shape matching and object recognition using shape contexts[J]. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2002, 24(4): 509-522.

[6] 孫景樂, 唐林波, 趙保軍, 等. 改進(jìn)的同心離散圓簇形狀描述方法[J]. 電子與信息學(xué)報(bào), 2013, 35(8): 1901-1906. Sun Jing-le, Tang Lin-bo, Zhao Bao-jun, et al.. An improved shape descriptor of cluster of concentric discrete circles[J]. Journal of Electronics & Information Technology, 2013, 35(8): 1901-1906.

[7] Bernier T and Landry J A. A new method for representing and matching shapes of natural objects[J]. Pattern Recognition, 2003, 36(8): 1711-1723.

[8] 王斌, 舒華忠, 施朝健, 等. 一種基于輪廓線的形狀描述與匹配方法[J]. 電子與信息學(xué)報(bào), 2008, 30(4): 949-952.

Wang Bin, Shu Hua-zhong, Shi Chao-jian, et al.. A contour-based shape description and matching method[J]. Journal of Electronics & Information Technology, 2008, 30(4): 949-952.

[9] Wang Zhi-yong, Chi Zhe-ru, and Feng Da-gan. Shape based leaf image retrieval[J]. IEE Proceedings-Vision, Image and Signal Processing, 2003, 150(1): 34-43.

[10] 王小玲, 謝康林. 一種新的方向碼描述的圖像檢索方法[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2006, 38(9): 1545-1548.

Wang Xiao-ling and Xie Kang-lin. Novel shape-based image retrieval using direction code[J]. Journal of Harbin Institute of Technology, 2006, 38(9): 1545-1548.

[11] Latecki L J, Lakaemper R, and Eckhatdt T. Shape descriptors for non-rigid shapes with a single closed contour[C]. Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Hilton Head Island, USA, 2000: 424-429.

[12] Zhang Deng-sheng and Lu Guo-jun. Shape-based image retrieval using generic Fourier descriptor[J]. Signal Processing: Image Communication, 2002, 17(10): 825-848.

[13] Torres R S and Falc?o A X. Contour salience descriptors for effective image retrieval and analysis[J]. Image Vision Computing, 2007, 25(1): 3-13.

[14] Liao S X and Pawlak M. On image analysis by moments[J]. IEEE Transactions on Pattern Analysis & Machine Intelligence, 1996, 18(3): 254-266.

[15] Nacéra L, Slimane L, Farouk L, et al.. Curve normalization for shape retrieval[J]. Signal Processing: Image Communication, 2014, 29(4): 556-571.

[16] Sebastian T B, Klein P N, and Kimia B B. Recognition of shapes by editing their shock graphs[J]. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2004, 26(5): 550-571.

[17] 王斌. 一種用于形狀描述的拱高半徑復(fù)函數(shù)[J]. 電子學(xué)報(bào), 2011, 39(4): 831-836. Wang Bin. Shape description using arc-height radius complex function[J]. Acta Electronica Sinica, 2011, 39(4): 831-836.

郭樹旭： 男，1959年生，博士，教授，研究方向?yàn)閳D像處理與信號(hào)分析.

趙 靜： 女，1989年生，碩士，研究方向?yàn)閳D像處理與模式識(shí)別.

李雪妍： 女，1980年生，博士，講師，研究方向?yàn)閳D像理解與模式識(shí)別.

Research on Shape Representation Based on Statistical Features of Centroid-contour Distance

Guo Shu-xu Zhao Jing Li Xue-yan
(College of Electronic Science and Engineering, Jilin University, Changchun 130012, China)

This paper proposes a novel shape representation method based on statistical features. According to the joint analysis on Centroid-Contour Distance (CCD) and chaincode, the silhouette is decomposed into several levels based on CCD. And then, the chaincode describing laying in each level is analyzed to extract the Joint Statistical of Centroid-Contour Distance and Chaincode (JSCCDC) descriptor for the silhouette. The similarity between different shapes can be measured by the city-block distance. Experiment results show that the proposed method describes both global and local features. Compared with traditional feature weighting method, JSCCDC is more accurate and reliable for shape matching and retrieval.

Pattern recognition; Shape representation; Feature statistics; Chaincode; Centroid-Contour Distance (CCD); Matching matrix

TP391.4

： A

：1009-5896(2015)06-1365-07

10.11999/JEIT140960

2014-07-21收到，2015-01-15改回

*通信作者：李雪妍 leexy@jlu.edu.cn