基于相對(duì)形狀上下文和譜匹配方法的點(diǎn)模式匹配算法

2010-03-27 06:55:38孫即祥李智勇陳明生

電子與信息學(xué)報(bào) 2010年10期

趙鍵孫即祥李智勇陳明生

(國(guó)防科學(xué)技術(shù)大學(xué)電子科學(xué)與工程學(xué)院信息工程系長(zhǎng)沙 410073)

1 引言

點(diǎn)模式(或稱點(diǎn)集)匹配廣泛應(yīng)用于圖像配準(zhǔn)[1]、圖像分類[2]與檢索[3]、目標(biāo)識(shí)別[4]、形狀匹配[5]和立體視覺(jué)[6]等領(lǐng)域。目前，點(diǎn)模式匹配算法大致可以分為兩大類[7]，一是基于變換關(guān)系求解的算法，是通過(guò)估計(jì)點(diǎn)模式之間的空間變換參數(shù)，利用該參數(shù)恢復(fù)或模擬點(diǎn)模式間的變換，從而求解點(diǎn)模式匹配問(wèn)題，也稱之為基于變換參數(shù)估計(jì)的算法。這類算法主要有迭代最近點(diǎn)算法[8]，軟指派算法[9]等。二是基于匹配關(guān)系求解的算法，是通過(guò)提取點(diǎn)集中點(diǎn)的特征而后運(yùn)用匹配識(shí)別方法獲得點(diǎn)模式間的匹配關(guān)系，從而求解點(diǎn)模式匹配問(wèn)題，或更形象地稱為基于特征的匹配算法。這類算法主要有基于形狀上下文的方法[10]，基于不變量特征的方法[11]以及基于譜圖論的方法[12?16]等。

基于譜圖論的方法是一類利用鄰接矩陣或者與其密切相關(guān)的Laplacian矩陣的特征值和特征矢量來(lái)刻畫點(diǎn)集全局結(jié)構(gòu)的方法[12]。Scott和Longuethiggins[13]首次將譜方法應(yīng)用于點(diǎn)模式匹配中，即通過(guò)定義點(diǎn)集之間的親近矩陣，并對(duì)該矩陣進(jìn)行SVD分解，以獲得對(duì)應(yīng)關(guān)系。該方法可處理不同大小的點(diǎn)集，但對(duì)較大角度的旋轉(zhuǎn)效果不好。為了克服文獻(xiàn)[13]方法的不足，Shapiro和Brady[14]采用了點(diǎn)集內(nèi)部點(diǎn)親近矩陣來(lái)進(jìn)行匹配。其方法為，對(duì)待匹配的每一個(gè)點(diǎn)集，首先利用其內(nèi)部點(diǎn)之間距離的高斯函數(shù)值來(lái)構(gòu)造親近矩陣，然后分別計(jì)算親近矩陣的特征值和特征向量，其對(duì)應(yīng)關(guān)系可以通過(guò)比較點(diǎn)集親近矩陣的有序特征向量來(lái)獲得。若所需匹配的點(diǎn)集大小相同，文獻(xiàn)[14]的方法對(duì)隨機(jī)點(diǎn)抖動(dòng)和較小的仿射變換都可以給出相對(duì)較好的結(jié)果。

上述譜方法的顯著優(yōu)點(diǎn)是構(gòu)造簡(jiǎn)單、計(jì)算量小，但由于它們是精確點(diǎn)模式匹配算法，因此，若待匹配的兩個(gè)點(diǎn)集大小不同以及存在位置噪聲時(shí)性能較差。針對(duì)該問(wèn)題，Carcassoni和Hancock[15]給出了數(shù)種構(gòu)造親近矩陣的方法，并將基于改進(jìn)后親近矩陣的譜方法與EM算法框架相結(jié)合從而提高了算法對(duì)點(diǎn)集大小和位置噪聲的魯棒性。由于該方法是迭代的，它相對(duì)較慢且對(duì)初始值比較敏感。Carcassoni等人[16]提出了基于模式聚類的分層算法，首先利用點(diǎn)集親近矩陣尋找點(diǎn)集間模式聚類中心的對(duì)應(yīng)性，再以聚類中心為約束確定該類中各點(diǎn)之間的匹配概率。雖然該算法在一定程度上改善了點(diǎn)集大小不同時(shí)的匹配性能，但僅適用于具備聚類特征的點(diǎn)集之間的匹配問(wèn)題。上述的幾種譜方法均要對(duì)親近矩陣的有序特征向量所組成的模式矩陣進(jìn)行截?cái)嘀?，才能進(jìn)行不同大小點(diǎn)集的匹配。相比之下，Leordeanu和Hebert[17]提出的基于成對(duì)約束的譜匹配算法(Spectral Matching，SM)無(wú)需進(jìn)行模式矩陣的截?cái)嘁材軐?shí)現(xiàn)非精確點(diǎn)模式匹配(即不同大小點(diǎn)集的匹配)。該算法先構(gòu)造分配圖及其親近矩陣，再求解親近矩陣最大特征值所對(duì)應(yīng)的主特征向量作為正確匹配點(diǎn)對(duì)的指示向量，從而實(shí)現(xiàn)點(diǎn)集之間的匹配對(duì)應(yīng)。若點(diǎn)集之間滿足等距變換的關(guān)系，該算法具有較好的抗噪聲和抗出格點(diǎn)的能力，但是該算法卻不能應(yīng)用于相似變換或其他更復(fù)雜變換的場(chǎng)合。針對(duì)上述問(wèn)題，本文提出了基于相對(duì)形狀上下文和譜匹配方法的點(diǎn)模式匹配算法(Relative Shape Context & Spectral Matching，RSC-SM)。與文獻(xiàn)[17]所提出的SM算法中分配圖親近矩陣的構(gòu)造方法不同，本文提出了一種新的基于點(diǎn)集的不變特征相對(duì)形狀上下文，并以點(diǎn)集間相對(duì)形狀上下文的統(tǒng)計(jì)檢驗(yàn)匹配測(cè)度作為新的相容性度量，以該相容性度量作為分配圖中各對(duì)應(yīng)邊的權(quán)值，最后利用分配圖的主特征向量來(lái)實(shí)現(xiàn)點(diǎn)模式匹配。

急性心肌梗死是一種臨床上病死率極高的心臟疾病，威脅著患者的生命健康，抓住搶救時(shí)間做好急救和護(hù)理有助于贏得下一步治療時(shí)間，提高恢復(fù)效果[3] 。臨床護(hù)理路徑的護(hù)理模式與急性心肌梗死患者的疾病治療護(hù)理需要有著密切關(guān)聯(lián)，是整體護(hù)理當(dāng)中的一個(gè)重要組成部分，能夠使各項(xiàng)護(hù)理工作井然有序，減少資源浪費(fèi)，提升護(hù)理質(zhì)量。此次研究表明，臨床護(hù)理路徑的護(hù)理方案有著以下優(yōu)勢(shì)：搶救流程設(shè)置合理，能夠提升搶救成功率，具備預(yù)見(jiàn)性和有序性特征；醫(yī)護(hù)合作效率高，可提升醫(yī)療安全；護(hù)理流程規(guī)范合理，可減少護(hù)理差錯(cuò)，提升護(hù)理完整性與延續(xù)性；縮短住院時(shí)間，降低醫(yī)療費(fèi)用；提升患者及其家屬滿意度，優(yōu)化護(hù)患關(guān)系。

2 相對(duì)形狀上下文

得到廣泛應(yīng)用的形狀上下文[10]是一種用于刻畫目標(biāo)形狀的描述子，點(diǎn)集中每個(gè)點(diǎn)的形狀上下文就是點(diǎn)集中其它點(diǎn)相對(duì)于該點(diǎn)的角度以及對(duì)數(shù)距離的直方圖分布。具體方法是對(duì)于任意點(diǎn)pi提取特征時(shí)，以其為原點(diǎn)把點(diǎn)集分布的空間按對(duì)數(shù)極坐標(biāo)劃分為60個(gè)部分，其中對(duì)數(shù)距離lgρ分為5部分，角度θ分為12部分。而后記錄除pi點(diǎn)外的n?1個(gè)點(diǎn)在這60個(gè)部分的分布數(shù)目bik(k=1,2,…,60)，由該數(shù)目來(lái)作為pi的特征。形狀上下文描述子具有平移不變性；為達(dá)到旋轉(zhuǎn)不變性，可以將每個(gè)點(diǎn)處的切線方向作為對(duì)數(shù)極坐標(biāo)軸的正方向；為具備尺度不變性，可以先求點(diǎn)集中所有點(diǎn)對(duì)之間的平均距離，然后用該值來(lái)規(guī)則化所有的點(diǎn)對(duì)距離。由于形狀上下文描述子在距離上的不均勻分組方式及其較為豐富的形狀特征描述能力，因此該描述子具有對(duì)局部區(qū)域形變的不變性。但是，當(dāng)出格點(diǎn)較多和噪聲較嚴(yán)重時(shí)，形狀上下文描述子的上述不變性就會(huì)變得較差。為了解決這些局限性，本文提出了一種新的魯棒性較強(qiáng)的形狀特征描述子相對(duì)形狀上下文(Relative Shape Context，RSC)。

定義1 如圖1所示，在具有n個(gè)點(diǎn)的點(diǎn)集P中，任意點(diǎn)pi相對(duì)于其他任意點(diǎn)pj的相對(duì)形狀上下文Hpij(s)就是在以pi為原點(diǎn)，矢量為正軸的極坐標(biāo)系下，點(diǎn)集中其他n?2個(gè)點(diǎn)相對(duì)于該極坐標(biāo)系正方向的極角所構(gòu)成的直方圖分布，即

圖1 相對(duì)形狀上下文描述子

式中#代表集合的勢(shì)，∠pjpip代表矢量與之間的逆時(shí)針旋轉(zhuǎn)角度，設(shè)點(diǎn)集空間按角度劃分為s個(gè)區(qū)間，則每區(qū)間角度跨度為Δθs= 360°/s ，則第k個(gè)角度區(qū)間為bin(k)=[(k?1)Δθs,kΔθs]。

(4)從L中剔除所有與v?=(i, i ′)相矛盾的分配，在本文的一對(duì)一匹配約束條件下需要從L中剔除掉諸如(i, k)和(q, i′)的分配，其中k≠i′，q≠i；

利用藍(lán)墨云班課創(chuàng)建課程和班級(jí)從而生成二維碼和邀請(qǐng)碼，所在班級(jí)學(xué)生掃碼或輸入邀請(qǐng)碼進(jìn)入班級(jí)，完成教師布置的預(yù)習(xí)任務(wù)，通過(guò)學(xué)生完成任務(wù)的質(zhì)量，教師課前智能診斷，從而安排或者調(diào)整教學(xué)內(nèi)容。以醫(yī)護(hù)英語(yǔ)項(xiàng)目模塊At thesurgery為例，筆者提前在云班課上傳了英國(guó)和美國(guó)一些診所背景知識(shí)的介紹，推薦了TED網(wǎng)易公開課和美劇《急診室的故事》，讓學(xué)生通過(guò)課前預(yù)習(xí)記錄下至少十個(gè)和醫(yī)護(hù)英語(yǔ)相關(guān)的詞匯并上傳到云班級(jí)，這樣可以讓學(xué)生在課前對(duì)所學(xué)習(xí)的內(nèi)容有一個(gè)基本的了解，激發(fā)他們的學(xué)習(xí)興趣。

(2)求解M的主特征向量x?，令x為解向量，并將其初始化為n×1的零向量。令L為分配圖G中所有頂點(diǎn)(即所有可能的分配)所組成的集合；

式(3)中的w{ vs, vt}反映了P中點(diǎn)對(duì)(i, j)的歐式距離特征與Q中點(diǎn)對(duì)(i′, j′)的歐式距離特征之間的相似程度，由于歐式距離是正交不變量，只能在等距變換下(即P與Q之間僅存在旋轉(zhuǎn)和平移變換時(shí))保持不變，如果當(dāng)兩個(gè)點(diǎn)集間若存在尺度變換或者其他更為復(fù)雜的變換時(shí)，顯然文獻(xiàn)[17]中的權(quán)值將不再適用；當(dāng)點(diǎn)集中存在出格點(diǎn)時(shí)以及存在由于隨機(jī)抖動(dòng)所產(chǎn)生的位置噪聲時(shí)，歐式距離特征的魯棒性較差。因此，需要定義新的更加魯棒的權(quán)值。本文第2節(jié)中提出的相對(duì)形狀上下文特征具有相似變換不變性，即平移、旋轉(zhuǎn)和尺度變換的不變性，而且對(duì)于噪聲和出格點(diǎn)也有較強(qiáng)的魯棒性，因此可以考慮利用其來(lái)代替歐式距離特征重新定義分配圖中邊的權(quán)值。具體定義如下，當(dāng)vs=(i, i′)，vt=(j, j′)時(shí)，定義P中點(diǎn)對(duì)(i, j)所確定的相對(duì)形狀上下文特征與Q中點(diǎn)對(duì)(i′, j′)所確定的相對(duì)形狀上下文特征之間的相似程度(相容性度量)，即分配圖中邊的權(quán)值w{ vs, vt}為式(4)中，Cij;i′j′為pi相對(duì)于pj的相對(duì)形狀上下文與qi′相對(duì)于qj′的相對(duì)形狀上下文的匹配測(cè)度；Cji;j′i′為pj相對(duì)于pi的相對(duì)形狀上下文與qj′相對(duì)于qi′的相對(duì)形狀上下文的匹配測(cè)度。Cij;i′j′與Cji;j′i′可以直接由式(2)求得，這里不再贅述。式(4)中w{ vs, vt}的物理意義為：當(dāng)pi與qi′是正確匹配點(diǎn)對(duì)，且pj與qj′也是正確匹配點(diǎn)對(duì)時(shí)，Cij;i′j′=0且Cji;j′i′=0從而w{ vs, vt}=1，此時(shí)表明vs=(i, i′)與vt=(j, j')共存的相容性最大；若w{ vs, vt}＜1，則表示vs=(i, i′)與vt=(j, j′)共存的相容性減小。這些特性不僅滿足分配圖中邊的權(quán)值所需要達(dá)到的要求，而且使得權(quán)值能適用于相似變換的場(chǎng)合，同時(shí)也增強(qiáng)了對(duì)噪聲和出格點(diǎn)的魯棒性。

頭孢菌素類是臨床另一類廣泛使用的抗菌藥物。本研究中，無(wú)論是ICU還是普通病房，大腸埃希菌和肺炎克雷伯菌對(duì)頭孢菌素均存在不同程度的耐藥。但兩種病房中，肺炎克雷伯菌對(duì)頭孢曲松的耐藥率均低于文獻(xiàn)[10]報(bào)道，可能與本研究中碳青霉烯類耐藥的肺炎克雷伯菌所占比例較低有關(guān)。本研究中，大腸埃希菌ESBLs在兩種病房的檢出率約為60%，與王啟等[12]研究一致。ESBLs陽(yáng)性可導(dǎo)致大腸埃希菌對(duì)所有青霉素類和頭孢菌素均耐藥，而酶抑制劑合劑如頭孢哌酮/舒巴坦和哌拉西林/他唑巴坦對(duì)陰性桿菌均有較好的敏感性，可作為產(chǎn)ESBLs腸桿菌科細(xì)菌感染的經(jīng)驗(yàn)用藥[13]。

3 基于相對(duì)形狀上下文和譜匹配方法的點(diǎn)模式匹配算法

本節(jié)首先簡(jiǎn)單介紹點(diǎn)模式匹配問(wèn)題的數(shù)學(xué)表達(dá)，再將相對(duì)形狀上下文特征應(yīng)用于分配圖及其親近矩陣的構(gòu)造過(guò)程，結(jié)合分配圖譜匹配方法，提出一種新的魯棒的點(diǎn)模式匹配算法—基于相對(duì)形狀上下文和譜匹配方法的點(diǎn)模式匹配算法。

3.1 點(diǎn)模式匹配問(wèn)題的提出

設(shè)模板點(diǎn)集為P={pi|i=1,2,…,nP}，目標(biāo)點(diǎn)集為Q={qj|j=1,2,…,nQ}。解決點(diǎn)模式匹配問(wèn)題的本質(zhì)在于求解點(diǎn)集P與Q之間的正確匹配映射關(guān)系。設(shè)P與Q之間的正確匹配映射為m:pi→qm(i)，即表明點(diǎn)pi匹配于點(diǎn)qm(i)。為簡(jiǎn)便起見(jiàn)，可令m(i)=i′，此時(shí)映射m可表示為m: i→i′。在映射m下，如果Q中點(diǎn)qj存在與之相匹配的模板點(diǎn)pi，則稱qj為Q中的內(nèi)點(diǎn)(inliers)；若點(diǎn)qj在映射m下無(wú)與之相匹配的模板點(diǎn)，則稱qj為Q中的出格點(diǎn)(outliers)。分配(assignments)定義為由點(diǎn)集P與Q的點(diǎn)所組成的組合(i, i′)，其中i∈P，i′∈Q。此時(shí)m又可定義為滿足所有正確分配(i, i′)的集合：Cm={(i, i′)|m(i)=i′, i∈P, i′∈Q}。所有屬于集合Cm的P(或Q)中的點(diǎn)稱為P(或Q)的內(nèi)點(diǎn)，不屬于集合Cm中任一分配的點(diǎn)稱為出格點(diǎn)。不同的場(chǎng)合所要求的匹配映射約束條件也不同，本文主要研究的是一對(duì)一的匹配約束，即點(diǎn)集P中某個(gè)點(diǎn)至多只能與點(diǎn)集Q中某一個(gè)點(diǎn)相匹配。因此，本文所研究的點(diǎn)模式匹配問(wèn)題實(shí)質(zhì)上就是求解滿足一對(duì)一匹配約束條件的所有正確分配所組成的集合Cm。

3.2 分配圖及其親近矩陣的構(gòu)造

首先，將P與Q之間所有可能的分配所組成的集合作為圖G的頂點(diǎn)集合，即V(G)={v1, v2,…,vn}，其中vs=(i, i′)，i∈P，i′∈Q，n=nPnQ。圖G的邊集合為E(G)={e1, e2,…,em}，其中ek={vs, vt}，即ek為以vs和vt為端點(diǎn)的無(wú)向邊，設(shè)w{ vs, vt}為邊ek上的權(quán)，此時(shí)所構(gòu)造的加權(quán)無(wú)向圖G稱為分配圖。各邊權(quán)值由匹配約束條件和成對(duì)幾何約束條件來(lái)共同確定：

(1)分配圖G中的任意兩個(gè)頂點(diǎn)設(shè)為vs=(i, i′)，vt=(j, j′)，其中i, j∈P，i′, j′∈Q。如果i=j而i′≠j′時(shí)(或當(dāng)i′=j′而i≠j時(shí))，即vs和vt所代表的分配與一對(duì)一匹配約束條件相矛盾，這兩種分配不可能同時(shí)存在，則令w{ vs, vt}=0。

綜上所述，Tomita評(píng)分可有效判斷預(yù)后，作為治療決策的參考指標(biāo)，對(duì)于Tomita評(píng)分4～7分的患者應(yīng)用球囊擴(kuò)張椎體后凸成形術(shù)治療轉(zhuǎn)移瘤性椎體壓縮骨折患者，能夠明顯降低VAS評(píng)分、ODI指數(shù)，恢復(fù)椎體高度，減輕痛苦，明顯提高總體生活質(zhì)量，并且具有較高的安全性。

(2)如果滿足了上述匹配約束條件，則w{ vs, vt}＞0。權(quán)值w{ vs, vt}可以理解為衡量由P中點(diǎn)對(duì)(i, j)所確定的特征與Q中點(diǎn)對(duì)(i′, j′)所確定的相應(yīng)特征之間的相似程度。一般來(lái)說(shuō)，w{ vs, vt}越大則表明此時(shí)vs和vt所代表的分配相容程度越高，越有可能都是正確的分配；反之則表明vs和vt所代表的分配相容性低，均是正確分配的可能性低。因此，w{ vs, vt}可以理解為vs和vt所代表的分配的相容性度量。