王貴錦, 劉 博, 何 貝, 張樹君, 孟 龍
1.清華大學(xué)電子工程系,北京100084
2.山東明佳包裝檢測科技有限公司,山東泰安271022
聚對苯二甲酸(polyethylene terephthalate,PET)瓶的自動缺陷檢測是現(xiàn)代工業(yè)中高速生產(chǎn)的需要,也是圖像處理和計算機(jī)視覺的重要組成部分[1].傳統(tǒng)的人工檢測存在準(zhǔn)確率低、速度慢以及成本高等問題,無法滿足現(xiàn)代生產(chǎn)的需求[2];新興的自動檢測則結(jié)合圖像處理和機(jī)器學(xué)習(xí)算法,對采集得到的PET圖像進(jìn)行分類識別,從而顯著提高了準(zhǔn)確率,減少了檢測時間,大大降低了生產(chǎn)成本[3],但由于PET瓶的瓶型、相機(jī)拍攝條件以及成像位置的差異,現(xiàn)有的檢測算法魯棒性較差.PET瓶的缺陷檢測主要包括瓶蓋、液位以及噴碼等部分,其中瓶蓋檢測最重要.本文主要圍繞PET瓶瓶蓋的快速、高精度的檢測展開.
現(xiàn)有的PET瓶瓶蓋檢測算法主要是圍繞持胚環(huán)定位展開的:通過確定持胚環(huán)位置,提取待檢測瓶蓋和模板瓶蓋的特征差異來進(jìn)行分類.文獻(xiàn)[4]通過對瓶蓋圖像進(jìn)行邊緣檢測、Hough變換等預(yù)處理,利用支持向量機(jī)(support vector machine,SVM)分類算法對定位出的瓶蓋進(jìn)行檢測.文獻(xiàn)[5]設(shè)計了一套基于現(xiàn)場可編程門陣列(field programmable gate array,FPGA)的嵌入式瓶蓋檢測系統(tǒng),根據(jù)持胚環(huán)和瓶蓋頂部的斜率之差來檢測瓶蓋是否正常.文獻(xiàn)[6]對瓶蓋區(qū)域進(jìn)行預(yù)處理和分割;在此基礎(chǔ)上,通過提取瓶蓋特征與模板圖像匹配檢測瓶蓋的質(zhì)量.德國Heuft公司[6]開發(fā)的PET瓶瓶蓋檢測算法根據(jù)持胚環(huán)在瓶頸處突出的特性來定位持胚環(huán),再通過模板匹配完成瓶蓋的分類.文獻(xiàn)[7]利用邊緣檢測算法精確估計出持胚環(huán)和瓶蓋頂部的直線方程,通過距離的衡量完成瓶蓋的快速檢測.文獻(xiàn)[8-9]在工控機(jī)和硬件平臺上分別給出了在線PET檢測設(shè)備的設(shè)計方法,主要依靠了二值化并進(jìn)行水平方向的投影,利用持胚環(huán)的定位來完成缺陷識別.但在實際生產(chǎn)中,高速的生產(chǎn)線運(yùn)動會引起瓶體傾斜,而清洗PET瓶會使得瓶體附著大量水珠,加上壓蓋時可能產(chǎn)生的大量毛刺,顯著降低了持胚環(huán)的定位精度,減少了檢測率;同時,現(xiàn)有的檢測算法無法滿足現(xiàn)代高速生產(chǎn)的需求.
針對瓶體傾斜、存在的毛刺和水珠等問題,本文提出了一種結(jié)合對稱匹配的PET瓶瓶蓋的實時檢測算法.根據(jù)PET瓶瓶體的對稱性,通過水平翻轉(zhuǎn)前后瓶體圖像的配準(zhǔn)完成瓶體對稱軸的估計,避免了瓶身傾斜對于持胚環(huán)定位精度的影響.同時,利用對稱軸兩側(cè)點(diǎn)對的相關(guān)性,消除了毛刺和附著水珠的影響,大大提高了檢測精度.仿真和高速生產(chǎn)線上的實驗證明了本算法的準(zhǔn)確性和實時性.
PET瓶瓶蓋存在的缺陷主要包括4種:無蓋、高蓋、歪蓋以及安全環(huán)斷裂,如圖1中的(a)~(d)所示.圖1(e)中矩形框標(biāo)出的部分就是持胚環(huán),它被認(rèn)為是瓶蓋中最穩(wěn)定的區(qū)域[7],因為存在缺陷的持胚環(huán)已經(jīng)在PET瓶瓶胚成型階段被剔除.現(xiàn)有的算法主要通過持胚環(huán)的定位來檢測瓶蓋,但會受到瓶身傾斜、存在的毛刺和水珠的影響.本檢測算法的關(guān)鍵在于利用PET瓶瓶體的對稱性穩(wěn)定地定位出持胚環(huán),從而完成瓶蓋分類.
圖1 瓶蓋缺陷示意圖Figure 1 Diagram of the defects of PET cap
本算法的流程如圖2所示.首先,對圖像進(jìn)行局部動態(tài)二值化,完成PET瓶瓶體和背景區(qū)域的圖像分割;其次,利用對稱匹配算法估計出瓶體對稱軸的位置;再次,對瓶身偏移角進(jìn)行補(bǔ)償,并消除毛刺和水珠的影響;最后,完成持胚環(huán)的定位以及PET瓶瓶蓋的檢測.
圖2 算法流程圖Figure 2 Flowchart of the algorithm
軸對稱的物體存在一條以上的對稱軸,如圖3中的左側(cè)虛線所示.在瓶體圖像中,存在一條直線y=kx+b,使得翻轉(zhuǎn)前后的瓶體關(guān)于該直線重合.同時,由于瓶體對稱,而毛刺、水珠可以認(rèn)為是隨機(jī)分布在瓶體周圍的,故可以通過布爾運(yùn)算消除其影響.因此,可以采用對稱匹配(symmetry matching,SM)算法來估計瓶體的對稱軸,補(bǔ)償瓶體的傾斜角度,消除毛刺、水珠的影響.對稱匹配算法包括全局[11-12]與局部[13-16]兩類,由于瓶蓋圖像背景單一,瓶蓋區(qū)域涵蓋了整個圖像的信息,且采用尺度不變特征變換(scaleinvariant feature transform,SIFT)、最大穩(wěn)定極值區(qū)域(maximally stable extremal regions,MSER)等旋轉(zhuǎn)尺度不變子的局部對稱匹配速度較慢,難以滿足實時檢測PET瓶瓶蓋的需求,于是本文設(shè)計了一種全局對稱匹配算法來完成對稱軸的精確估計.
圖3 對稱圖像水平翻轉(zhuǎn)前后的對應(yīng)點(diǎn)變換關(guān)系Figure 3 Corresponding points between the original PET cap image and the horizontally f lipped one
如果將整個圖像按照中垂線x=W/2(其中W為圖像寬度)進(jìn)行水平翻轉(zhuǎn)(如圖3所示,左側(cè)為原圖而右側(cè)為水平翻轉(zhuǎn)后的圖像),那么水平翻轉(zhuǎn)前后的圖像間的對應(yīng)點(diǎn)滿足如下條件:橫坐標(biāo)之和為W且縱坐標(biāo)相同.假設(shè)在原對稱圖像中(x1,y1)和(W-x2,y2)關(guān)于y=kx+b對稱,那么翻轉(zhuǎn)變換前后圖像的配準(zhǔn)等價于求解(x1,y1)和(x2,y2)的變換關(guān)系.因此,根據(jù)對稱軸的性質(zhì)有
求解x1和y1,那么由于斜率k=tanβ=sinβ/cosβ,旋轉(zhuǎn)角α=2β-π,代入式(2)并利用三角函數(shù)的性質(zhì),可將對應(yīng)點(diǎn)的變換關(guān)系轉(zhuǎn)化為
在PET瓶的實際生產(chǎn)中,由于瓶身傾斜較小,旋轉(zhuǎn)角控制在-10°~10°之間,因此近似有sinα≈α,cosα≈1.即
式中,Tx=αb-W,表示x方向的偏移量.不難看出,翻轉(zhuǎn)變換前后的圖像之間僅存在2個自由度的變換[20].與Lucas-Kanade[20]配準(zhǔn)算法類似,定義翻轉(zhuǎn)變換前后圖像I1和I2之間的誤差為
對式(5)進(jìn)行Taylor展開可以得到
利用誤差函數(shù)E對兩個未知參數(shù)進(jìn)行求導(dǎo),并賦為0,那么有
式中,c=I2(x,y)-I1(x,y),其物理意義為翻轉(zhuǎn)變換前后的圖像誤差;矢量A如式(8)所示,對應(yīng)于Lucas-Kanade[20]配準(zhǔn)算法中的Hessian矩陣
根據(jù)最速下降法可以逐步修正翻轉(zhuǎn)后的圖像I2,使其與模板圖像I1匹配,進(jìn)而完成變換系數(shù)的迭代估計.
本算法的圖像采集設(shè)備放置于生產(chǎn)線上相對封閉的處理箱中,且打光方式為背光源模式,因此捕獲的圖像中瓶蓋區(qū)域和背景區(qū)域的分離度高且穩(wěn)定[10].在透光率低的地方,比如瓶蓋部分,像素為黑色;在透光率高的地方,比如瓶體部分和背景部分,像素為白色.如果要分割整個PET瓶和背景區(qū)域,需要考慮瓶體的完整性,避免瓶體內(nèi)部被劃入背景部分.
本文的瓶體分割算法主要利用二值化算法進(jìn)行,但PET瓶瓶壁較薄,直接進(jìn)行固定閾值二值化會導(dǎo)致瓶壁部分被劃為背景區(qū)域[11].本文涉及了局部動態(tài)二值化算法對瓶蓋進(jìn)行處理,如式(9)所示:
式中,Ik為像素點(diǎn)k處的灰度值,Rd表示以點(diǎn)k為中心的正方形鄰域且邊長為d,mean(·)為均值函數(shù),γ為比例系數(shù),Bk為二值化后的灰度值.下面從背景區(qū)域出發(fā),利用4鄰域的連通域分析即可完成瓶體區(qū)域的分割.如圖4所示,本文的局部動態(tài)二值化算法能夠自適應(yīng)地修正局部閾值,避免固定閾值對于瓶體等透光率較高區(qū)域的影響.
圖4 瓶蓋圖像的二值化結(jié)果對比Figure 4 Comparison of the results of different thresholding algorithms to cap image
根據(jù)第2節(jié)介紹的對稱匹配算法,可以利用圖像配準(zhǔn)估計出兩個變換參數(shù):旋轉(zhuǎn)角α和x軸方向的平移量Tx,計算出瓶體對稱軸所在直線的參數(shù),即斜率k≈-2/α以及截距b≈-(W-Tx)/α,進(jìn)而可以得到瓶體對稱軸的傾斜角θ=arctan k.通過對該角的補(bǔ)償即可保證瓶體圖像的水平,避免了生產(chǎn)線上由于鏈道運(yùn)動造成的瓶體傾斜的影響.
由于PET瓶瓶體關(guān)于對稱軸重合,瓶體左右兩側(cè)應(yīng)該是完全相同的;但由于毛刺、水珠的影響,該條件并不能完全滿足.因此,本算法將角度補(bǔ)償后的PET圖像關(guān)于對稱軸進(jìn)行水平翻轉(zhuǎn),再對兩側(cè)圖像進(jìn)行布爾運(yùn)算,如式(10)所示:
式中,OR(i,j)表示或運(yùn)算,當(dāng)且僅當(dāng)i和j都為0時結(jié)果為0.也就是說,對于每個像素點(diǎn),只有其本身和關(guān)于對稱軸對稱的像素點(diǎn)都屬于瓶體區(qū)域時,在結(jié)果圖像中才歸為瓶體區(qū)域.
如圖5所示,對稱軸的估計可以消除瓶身傾斜的影響.這一方面避免影響后續(xù)瓶蓋分類的計算精度;另一方面大大減少了毛刺和水珠,避免影響后續(xù)持胚環(huán)的定位精度.圖5(a)中的毛刺較大,這是無法通過簡單的形態(tài)學(xué)操作來完成的.
圖5 對稱軸補(bǔ)償?shù)慕Y(jié)果示意圖Figure 5 Effect of rotational compensation
持胚環(huán)是PET瓶瓶頸到瓶蓋第1個寬度最大的位置.考慮到瓶蓋上下抖動帶來的持胚環(huán)位置的變化,本文利用豎直方向的投影進(jìn)行估計.首先,對瓶蓋圖像進(jìn)行黑色像素點(diǎn)的豎直投影;其次,對投影向量進(jìn)行中值濾波,以防止殘留噪聲點(diǎn)和毛刺對于定位的影響;最后,定位從瓶頸到瓶蓋的第1個局部最大的位置作為持胚環(huán)的所在位置.
以持胚環(huán)的中點(diǎn)作為坐標(biāo)原點(diǎn),持胚環(huán)所在的直線為x軸,將待檢測瓶蓋和模板瓶蓋圖像投影到同一個坐標(biāo)系中.通過衡量待檢測瓶蓋和模板瓶蓋之間的相關(guān)度,即可完成瓶蓋的分類.本文采用級聯(lián)弱分類器得到強(qiáng)分類器的策略[12],在瓶蓋圖像上構(gòu)建若干矩形框,并將每個矩形框中白色像素占所有像素的百分比作為一個弱分類特征.強(qiáng)分類器C的表達(dá)式為
式中,Ci表示第i個弱分類器.本文選取直方圖特征作為弱分類器,其定義如式(12)所示:
式中,si表示待檢測瓶蓋圖像在第i個分類器中白色像素的百分比,表示模板瓶蓋圖像在該分類器中白色像素的百分比,δ表示百分比的正常波動范圍,min(·)和max(·)分別表示最小值函數(shù)與最大值函數(shù).
本文實驗的硬件配置為Intel Core(TM)2 Quad CPU處理器、2.5 GHz主頻以及3.2 5G內(nèi)存,編譯環(huán)境為VC2008,通過640×480、8 bpp的網(wǎng)絡(luò)相機(jī)抓拍瓶蓋圖像.本算法的主要參數(shù)設(shè)置如下:d=3,γ=0.85,δ=0.05.
首先需要評估對稱匹配算法的精度和相應(yīng)的計算時間.實驗選取了4種具有代表意義的PET瓶,包括輕量蓋的冰露礦泉水、大蓋的達(dá)能脈動、普通蓋的可樂以及加重蓋的營養(yǎng)快線,其中每種瓶型保留100個PET瓶.將PET瓶水平放置,利用靜態(tài)圖像采集模式獲得瓶蓋圖像,此時瓶體對稱軸偏離中垂線0°.對于每個PET瓶采集到的瓶蓋圖像,以[-10°,10°]、間隔0.1°進(jìn)行旋轉(zhuǎn)變換,共產(chǎn)生201張樣本圖像,于是可以利用對稱匹配算法對旋轉(zhuǎn)角度進(jìn)行估計.設(shè)第k個瓶子在角度αi時的估計值為,則估計誤差根據(jù)式(13),利用均值和標(biāo)準(zhǔn)差來衡量估計誤差和時間
圖6 對稱匹配算法的性能Figure 6 Performanceof thesymmetric matching algorithm
在上述構(gòu)建的不同旋轉(zhuǎn)角度的情況下,將本文算法與無對稱匹配的PET瓶瓶蓋檢測算法進(jìn)行比較,其中每個角度都選擇200張正常蓋、100張歪蓋、50張高蓋、50張安全環(huán)斷裂和50張無蓋(共250張不正常蓋).定義FN(false negative)表示正常蓋檢測為不正常的比率,即誤檢率;FP(false positive)表示不正常蓋檢測為正常的比率,即漏檢率,得到的檢測結(jié)果如圖7所示.當(dāng)旋轉(zhuǎn)角度較小時,有/無對稱匹配的算法都可以得到較好的檢測結(jié)果,特別是有對稱匹配的算法;當(dāng)旋轉(zhuǎn)角度增大時,有對稱匹配的算法能夠保持檢測結(jié)果的穩(wěn)定性,而無對稱匹配的算法會使得很多正常瓶蓋被誤判為不正常瓶蓋,顯著提高了誤檢率FN.除此之外還可以看出,誤檢率FN通常高于漏檢率FP,這是由于正常蓋情況相對單一而不正常蓋情況相對復(fù)雜.
在生產(chǎn)線的實時測試中,將本文算法與無對稱匹配的PET瓶瓶蓋檢測算法和文獻(xiàn)[7]的結(jié)果進(jìn)行比較.實驗包括20000瓶的靜態(tài)測試(PET瓶靜止放置在相機(jī)前)以及20000瓶的動態(tài)測試(生產(chǎn)線上實時采集).PET的測試用瓶包括輕量蓋的冰露礦泉水、大蓋的達(dá)能脈動、普通蓋的可樂以及加重蓋的營養(yǎng)快線等4種具有代表意義的瓶型.與靜態(tài)圖像相比,動態(tài)圖像存在瓶體傾斜、毛刺和水珠的影響.對于每一類圖像,分別選取5000張正常蓋、2000張歪蓋、1000張高蓋、1000張安全環(huán)斷裂以及1000張無蓋(共5000張不正常蓋).實驗結(jié)果如表1和2所示,其中最后一行給出了不同算法的平均誤檢率和漏檢率.
從表1和2中不難得出以下2點(diǎn)結(jié)論:
1)對于靜態(tài)采集的瓶蓋圖像,有/無對稱匹配的檢測算法和文獻(xiàn)[7]都能夠獲得較好的檢測結(jié)果;但當(dāng)生產(chǎn)線處于實時采集過程中時,由于瓶體傾斜、毛刺和水珠的影響,無對稱匹配的檢測算法和文獻(xiàn)[7]的性能顯著下降,而本文中的有對稱匹配算法依然能夠保持較高的檢測率.
圖7 有/無對稱匹配算法的檢測結(jié)果比較Figure 7 Result comparison of the defect inspection system with/without symmetric matching algorithm
表1 有/無對稱匹配算法和文獻(xiàn)[7]的誤檢(FN)結(jié)果比較Table 1 False negative result comparison of the defect inspection system with/without symmetric matching algorithm and reference[7] %
表2 有/無對稱匹配算法和文獻(xiàn)[7]的漏檢(FP)結(jié)果比較Table 2 False positive result comparison of the defect inspection system with/without symmetric matching algorithm and reference[7] %
2)在檢測過程中,3種算法的漏檢率均高于誤檢率,這是由于正常蓋的特征分布單一而不正常蓋的特征分布多樣化引起的.瓶體傾斜、毛刺和水珠的作用都會使得持胚環(huán)定位錯誤,從而造成瓶蓋圖像直方圖投影不準(zhǔn)確,使得瓶蓋特征在正常蓋的特征分布閾值之外,導(dǎo)致誤檢率提高.
3)本文算法的平均檢測速度為27.4 ms/幀,能滿足高速自動化生產(chǎn)需求.
表3比較了有/無對稱匹配的檢測算法和文獻(xiàn)[7]等3種算法的計算時間.首先,本算法需要迭代優(yōu)化翻轉(zhuǎn)前后的圖像誤差,因此會產(chǎn)生較高的計算復(fù)雜度;而無對稱匹配的PET檢測算法的速度則是優(yōu)于文獻(xiàn)[7]的.其次,3種算法都能達(dá)到30瓶/s,能夠滿足現(xiàn)有各種高速生產(chǎn)線的要求.最后,由于動態(tài)測試中存在瓶體傾斜、毛刺和水珠等問題,因此本算法需要更高的時間復(fù)雜度用于迭代配準(zhǔn),但平均速度仍然高達(dá)26ms/幀左右.
表3 有/無對稱匹配算法和文獻(xiàn)[7]的平均檢測時間比較Table 3 Processing time comparison of the defect inspection system with/without symmetric matching algorithm and reference[7] ms
本文提出了一種結(jié)合對稱匹配的PET瓶瓶蓋的檢測算法.現(xiàn)有的瓶蓋檢測算法依賴于持胚環(huán)的定位,主要面臨瓶體傾斜、瓶頸毛刺以及附著水珠的問題.本算法利用PET瓶瓶體的對稱性,通過配準(zhǔn)水平翻轉(zhuǎn)前后的圖像定位得到瓶體對稱軸,補(bǔ)償了瓶體傾斜的影響,通過瓶蓋圖像對稱相關(guān)的特性去除了毛刺和附著水珠對持胚環(huán)檢測的影響.實驗證明,本算法具有很高的準(zhǔn)確性和魯棒性,能滿足高速生產(chǎn)的檢測需求.
[1]MALAMAS E N,PETRAKIS E G M,ZERVAKIS M,PETIT L,LEGAT J D.A survey on industrial vision systems,applications and tools[J].Image and Vision Computing,2003,21(2):171-188.
[2]DUANF,WANGY,LIUH,LIY.A machine vision inspector for beer bottle[J].Engineering Applications of Artif icial Intelligence,2007,20(7):1013-1021.
[3]孫懷遠(yuǎn),周夫之,陸尚卿.機(jī)器視覺技術(shù)在藥液灌裝量檢測中的應(yīng)用[J].計算機(jī)應(yīng)用與軟件,2009,6:247-248.
SUNHuaiyuan,ZHOUFuzhi,LUShangqing.Applying machine vision technology in f illing quantity detection of liquid medicine[J].Computer Applications and Software,2009,6:247-248.(in Chinese)
[4]LIU H,WANG Y,DUAN F.An empty bottle intelligent inspector based on support vector machines and fuzzy theory[C]//Conference Proceedings of Sixth World Congress on Intelligent Control and Automation,2006:9739-9743.
[5]鄒振興,謝云,林奇鴻,陳炳成.基于SOPC技術(shù)的PET瓶缺陷檢測系統(tǒng)設(shè)計[J].現(xiàn)代電子技術(shù),2009,10:58-61.
ZOU Zhenxing,XIE Yun,LIN Qihong,CHEN Bingcheng.Design of PET bottle defect detection system based on SOPC[J].Modern Electronics Technique,2009,10:58-61.(in Chinese)
[6]YAZDIL,PRABUw ONOA S,GOLKAR E.Feature extraction algorithm for f ill level and cap inspection in bottling machine[C]//International Conference on Pattern Analysis and Intelligent Robotics,2011,1:47-52.
[7]鄭云峰,王貴錦,何貝,李艷麗,林行剛.基于直線擬合算法的PET瓶瓶蓋檢測方法[J].計算機(jī)應(yīng)用研究,2011,28(11):4398-4400.
ZHENGYunfeng,WANGGuijin,HEBei,LIYanli,LIN Xinggang.PET bottle cap inspection method based on linear f itting algorithm[J].Application Research of Computers,2011,28(11):4398-4400.(in Chinese)
[8]張樹君,陳大千,辛瑩瑩,國磊,姜輝.基于圖像的PET瓶瓶蓋和液位的檢測設(shè)備[J].輕工機(jī)械,2013(2):71-73.
ZHANG Shujun,CHEN Daqian,XIN Yingying,GUO Lei,JIANG Hui.Detection device of PET bottle and liquid based on image[J].Light Industry Machinery,2013,31(2):71-74.(in Chinese)
[9]鄭海成,謝云.PET瓶封蓋質(zhì)量在線檢測系統(tǒng)的設(shè)計與實現(xiàn)[J].電子技術(shù)應(yīng)用,2012,38(8):20-22.
ZHENGHaicheng,XIEYun.Design and implementation of PET bottle package defect detection system based on dual core processor[J].Application of Electronic Technique,2012,38(8):20-22.(in Chinese)
[10]CARRASCO M,PIZARRO L,MERY D.Image acquisition and automated inspection of wine bottlenecks by tracking in multiple views[C]//Iscgav'08:Proceedings of the 8th Wseas International Conference on Signal Processing,Computational Geometry and Artif icial Vision,2008:84-89,235.
[11]NGH F.Automatic thresholding for defect detection[J].Pattern Recognition Letters,2006,27(14):1644-1649.
[12]章品正,王健宏.一種應(yīng)用機(jī)器學(xué)習(xí)的車牌定位方法[J].應(yīng)用科學(xué)學(xué)報,2011,29(2):147-152.
ZHANG Pinzheng,WANG Jianhong.Vehicle license plate location based on machine learning[J].Journal of Applied Sciences,2011,29(2):147-152.(in Chinese)
[13]TZIMIROPOULOSG,ARGYRIOUV,STATHAKIT.Symmetry detection using frequency domain motion estimation techniques[C]//International Conference on Acoustics,Speech,and Signal Processing,2008:861-864.
[14]KELLER Y,SHKOLNISKY Y.A signal processing approach to symmetry detection[J].IEEE Transaction on Image Processing,2006,15(8):2198-2207.
[15]LOYG,EKLUNDHJ O.Detecting symmetry and symmetric constellations of features[C]//European Conference on Computer Vision,2006:508-521.
[16]CHOM,LEEK M.Bilateral symmetry detection via symmetry-growing[C]//Proceedings of the British Machine Vision Conference,2009:1-8.
[17]HE Bei,WANG Guijin,YIN Xuanwu,LIU Bo,LIN Xinggang.Self-clustering symmetry detection[J].IEICE Transactions on Information and Systems,2012,95(9):2359-2362.
[18]HAUAGGED C,SNAVELYN.Imagematching using local symmetry features[C]//IEEE Computer Society Conference on Computer Vision and Pattern Recognition,2012:206-213.
[19]ZITOVA B,FLUSSER J.Image registration methods:a survey[J].Image and Vision Computing,2003,21(11):977-1000.
[20]BAKER S,MATTHEw SI.Lucas-Kanade 20 years on:a unifying framework[J].International Journal of Computer Vision,2004,56(3):221-255.