高分辨率遙感圖像SIFT和SURF算法匹配性能研究

齊冰潔，劉金國(guó)，張博研，左 洋，呂世良

(1.中國(guó)科學(xué)院 長(zhǎng)春光學(xué)精密機(jī)械與物理研究所，吉林 長(zhǎng)春 130033；2.中國(guó)科學(xué)院大學(xué)，北京 100049)

高分辨率遙感圖像SIFT和SURF算法匹配性能研究

齊冰潔，劉金國(guó)*，張博研，左 洋，呂世良

(1.中國(guó)科學(xué)院 長(zhǎng)春光學(xué)精密機(jī)械與物理研究所，吉林 長(zhǎng)春 130033；2.中國(guó)科學(xué)院大學(xué)，北京 100049)

遙感圖像匹配是圖像校正、拼接的基礎(chǔ)。由于遙感圖像特征相似度大，重疊區(qū)域小，遙感圖像對(duì)匹配算法的要求更高。本文首先從特征檢測(cè)、特征描述和特征匹配三個(gè)方面，比較了SIFT算法和SURF算法在計(jì)算速度和準(zhǔn)確度方面性能，然后研究了算法對(duì)遙感圖像重疊度、度量距離的要求，并針對(duì)SURF算法對(duì)特征方向誤差敏感的特點(diǎn)，提出一種oSURF算法；最后利用衛(wèi)星1A級(jí)條帶遙感圖像分析各個(gè)算法優(yōu)劣性。測(cè)試結(jié)果表明，相比于SIFT算法，SURF算法計(jì)算速度為SIFT的3倍，需要的圖像重疊寬度僅為1.25倍描述向量尺寸，而在保證同樣匹配率的情況下，SIFT算法則需要圖像重疊寬度為1.5倍描述向量尺寸。本文提出的 oSURF算法在保證計(jì)算速度的同時(shí)，準(zhǔn)確度相對(duì)于SURF算法提升5%～10%，因此，oSURF算法更適合1A級(jí)條帶遙感圖像的拼接。

遙感圖像；特征匹配；SIFT；SURF；重疊度

1 引 言

遙感圖像匹配是對(duì)不同傳感器、不同時(shí)相、不同視角的遙感圖像在其重疊部分計(jì)算同名點(diǎn)(Homologous Points)的過(guò)程。由于遙感圖像高分辨率和大視場(chǎng)的矛盾性，為了保證分辨率，衛(wèi)星相機(jī)只推掃條帶狀影像數(shù)據(jù)，再通過(guò)圖像拼接、融合等步驟實(shí)現(xiàn)大范圍圖像覆蓋率，而這些處理都是以圖像匹配為前提。圖像特征匹配包括檢測(cè)特征點(diǎn)、提取描述向量和向量匹配3個(gè)步驟。

本文選取計(jì)算機(jī)視覺(jué)領(lǐng)域中兩種經(jīng)典匹配算法：SIFT(Scale Invariant Feature Transform)和SURF(Speed Up Robust Features)進(jìn)行性能分析，SIFT算法是由David G.lowe[1-2]在1999年提出，該算法對(duì)圖像旋轉(zhuǎn)、縮放、光照變化都具有很好的魯棒性，但算法的計(jì)算速度較慢。SURF算法由Herbert Bay[3-4]在2006年提出，在保證算法魯棒性的基礎(chǔ)上對(duì)SIFT算法進(jìn)行了加速改進(jìn)。近幾年，Luo Juan、索春寶等人[5-9]利用Mikolajczyk和Schimd數(shù)據(jù)集，分析了兩種算法在普通圖像上的魯棒性。由于遙感圖像數(shù)據(jù)量大，一般為1G左右，導(dǎo)致計(jì)算量很大。而且遙感圖像特征相似度大，衛(wèi)星條帶圖像重疊區(qū)域較小，易受到非重疊區(qū)域影響，這些特點(diǎn)導(dǎo)致遙感圖像對(duì)匹配算法有很高要求。目前國(guó)內(nèi)張?jiān)粕?、羅楠等人[10-13]針對(duì)遙感圖像，對(duì)SIFT和SURF算法匹配性能做出了研究。南京理工大學(xué)的羅楠針對(duì)同一傳感器不同時(shí)相遙感圖像和不同傳感器遙感圖像分別分析兩種算法的效率，并提出一種基于擴(kuò)展SURF的描述符；葛盼盼在分析算法的匹配效率基礎(chǔ)上，比較了峰值信噪比(PNSR)、相似性度量(SSIM)和互信息量(MI)4個(gè)方面的匹配效果。但目前針對(duì)SIFT和SURF算法在遙感圖像邊緣處特性分析很少，比如特征鄰域尺寸和非重疊區(qū)域尺寸等影響，而且大部分實(shí)驗(yàn)數(shù)據(jù)均經(jīng)過(guò)深度處理，并不能充分反映遙感圖像特點(diǎn)。

根據(jù)衛(wèi)星圖像后期處理程度逐漸升高，可以分為0級(jí)、1A級(jí)、1B級(jí)、2級(jí)、3級(jí)和4級(jí)。1A級(jí)圖像是傳感器采集的原始圖像進(jìn)行初步均衡化輻射較正后的數(shù)據(jù)產(chǎn)品，該級(jí)數(shù)據(jù)可以充分反映遙感圖像特點(diǎn)。本文利用衛(wèi)星1A級(jí)遙感數(shù)據(jù)，對(duì)SIFT和SURF算法在遙感圖像匹配速度、準(zhǔn)確度以及非重疊區(qū)域影響等性能進(jìn)行對(duì)比分析。

2 實(shí)驗(yàn)測(cè)試與性能分析

本節(jié)在相同計(jì)算機(jī)配置下，從特征點(diǎn)檢測(cè)、特征向量描述和特征點(diǎn)匹配三方面，比較了SIFT和SURF算法計(jì)算速度、準(zhǔn)確度及重疊度等方面性能。

實(shí)驗(yàn)圖像為美國(guó)某海岸1A級(jí)遙感數(shù)據(jù)，圖1中(a)、(b)為兩條邊緣重疊的衛(wèi)星條帶圖像，本文采用的圖像重疊寬度為采集條帶寬度的1.23%。由于匹配前圖像重疊寬度未知，因此根據(jù)衛(wèi)星推掃特點(diǎn)，將特征檢測(cè)寬度限定為100像素(圖中框內(nèi)為檢測(cè)區(qū)域，為1 000×100像素)。

2.1 特征點(diǎn)檢測(cè)

首先進(jìn)行特征點(diǎn)檢測(cè)，本節(jié)從檢測(cè)時(shí)間和特征點(diǎn)數(shù)量?jī)煞矫姹容^算法性能。為了便于比較和顯示，本文對(duì)算法閾值進(jìn)行設(shè)定以確保檢測(cè)到特征點(diǎn)數(shù)量相當(dāng)，SIFT算法檢測(cè)閾值設(shè)定為0.071 2，SURF算法檢測(cè)閾值為0.000 6，檢測(cè)結(jié)果如圖2。

圖2 SIFT和SURF算法特征點(diǎn)檢測(cè)效果 Fig.2 Detection result of SIFT and SURF algorithm

表1為不同算法檢測(cè)性能比較，可以看出SURF算法檢測(cè)特征點(diǎn)速度相對(duì)SIFT算法有明顯提升，約是SIFT的3倍。

表1 SIFT和SURF算法特征檢測(cè)結(jié)果

SIFT算法是通過(guò)不同尺寸的高斯模板進(jìn)行濾波和降采樣形成尺度空間，再將相鄰層相減得到高斯差分DOG尺度空間[14]，通過(guò)非極大值抑制法尋找特征點(diǎn)。SURF算法尺度空間的建立不改變圖像的大小，而是改變?yōu)V波模板尺寸，通過(guò)Hessian矩陣行列式值來(lái)判斷特征點(diǎn)。SIFT和SURF算法通過(guò)建立尺度空間來(lái)實(shí)現(xiàn)尺度不變性，而兩種算法的檢測(cè)時(shí)間差也主要在尺度空間形成過(guò)程中產(chǎn)生，原因有兩個(gè)方面：

(1)SURF算法通過(guò)框狀濾波器代替高斯二階偏導(dǎo)，結(jié)合Paul Viola在2001年提出的積分圖像(Intergral Image)[15]的概念形成一種近似卷積模板，將SIFT算法中高斯濾波卷積運(yùn)算簡(jiǎn)化成框狀濾波器4個(gè)頂點(diǎn)積分值加減法運(yùn)算。這種近似卷積模板的計(jì)算量獨(dú)立于模板尺寸，使得SURF算法運(yùn)算量大大降低。

(2)SIFT尺度空間通過(guò)降采樣改變圖像大小形成，需要層層遞推，而SURF尺度空間是改變盒子濾波尺寸，圖像的大小始終不變，因此不同尺度層可以同時(shí)建立。另一方面，由于圖像大小不變，SURF算法只需一個(gè)積分圖像，在濾波同時(shí)可以直接建立尺度空間，而不需要對(duì)圖像進(jìn)行降采樣處理，實(shí)現(xiàn)了計(jì)算量的優(yōu)化。

綜上所述，由于SURF算法采用了積分圖像，在檢測(cè)特征點(diǎn)相同的情況下，SURF算法的檢測(cè)速度比SIFT算法快3倍。

2.2 特征向量描述

檢測(cè)特征點(diǎn)之后，需要對(duì)特征點(diǎn)進(jìn)行局部信息描述，分為確定主方向和提取描述向量?jī)刹?。本?jié)從構(gòu)建向量速度、特征向量數(shù)目、描述鄰域尺寸三方面進(jìn)行比較和分析。本文中SURF算法特征向量采取64維和128維，由于SIFT算法受特征鄰域劃分規(guī)則所限，特征向量維度只有32和128，而32維度獨(dú)特性較差，故本文SIFT算法向量維度只選取128。表2是對(duì)SIFT和SURF算法檢測(cè)的特征點(diǎn)進(jìn)行描述的結(jié)果。

2.2.1 算法構(gòu)建向量速度分析

從表2可知，SIFT算法描述特征時(shí)間相對(duì)SURF算法更長(zhǎng)。由于SIFT算法需要通過(guò)計(jì)算特征點(diǎn)鄰域梯度信息構(gòu)建鄰域梯度直方圖來(lái)實(shí)現(xiàn)，而SURF算法直接利用積分圖結(jié)合haar小波濾波器生成的x、y方向上響應(yīng)實(shí)現(xiàn)，從而優(yōu)化了時(shí)間復(fù)雜度。另外，SURF算法的64維和128維向量雖然維度增加1倍，但是特征向量生成時(shí)間相差不大。

2.2.2 SIFT特征向量數(shù)目增加分析

由表2可知SIFT算法特征向量數(shù)目相對(duì)于檢測(cè)特征點(diǎn)數(shù)明顯增加，是由于SIFT算法不同于只有一個(gè)主方向的SURF算法，SIFT算法選取梯度方向直方圖中所有大于主方向0.8倍的方向作為特征方向，每一個(gè)特征方向作為一個(gè)特征點(diǎn)，因此向量數(shù)大大增加。這些輔助方向的存在有利于提高SIFT匹配穩(wěn)定性。

2.2.3 描述鄰域尺寸影響分析

由表2可知，特征描述后SURF算法特征數(shù)降低了，事實(shí)上SIFT算法特征點(diǎn)也有所降低，由于輔助方向的存在而沒(méi)有表現(xiàn)出來(lái)。因?yàn)镾IFT和SURF算法都屬于局部特征檢測(cè)算法[16-18]，即描述特征方向和構(gòu)建特征向量時(shí)需要特征點(diǎn)一定鄰域內(nèi)灰度信息。而遙感圖像重疊區(qū)域都集中在邊界，而且相對(duì)較小。因此，當(dāng)特征點(diǎn)與邊緣的距離小于鄰域半徑時(shí)，該特征點(diǎn)無(wú)法被描述。

表2 不同算法特征向量描述結(jié)果

圖3 檢測(cè)區(qū)域?qū)挾炔煌瑫r(shí)描述向量個(gè)數(shù)變化 Fig.3 Vector number change with the detection area

圖4 重疊區(qū)域不同時(shí)同名點(diǎn)概率變化 Fig.4 Probability change of Homologous Points with the overlap region

圖3為SIFT和SURF兩種算法檢測(cè)范圍變化時(shí)檢測(cè)到特征點(diǎn)數(shù)量。由于構(gòu)建向量的鄰域比提取主方向鄰域大，因此以構(gòu)建向量鄰域作為主要影響參數(shù)。為了保證尺度不變性，特征描述是在檢測(cè)到特征點(diǎn)的尺度層上進(jìn)行的。無(wú)論是SIFT還是SURF算法特征鄰域大小都與金字塔層對(duì)應(yīng)sigma值成正比，sigma越大鄰域越大，式(1)和式(2)分別為本文SIFT和SURF算法特征向量鄰域大小，根據(jù)公式可知，SIFT和SURF算法最小特征向量鄰域分別為40和20個(gè)像素。由圖3可知，兩種算法均在超過(guò)最小鄰域范圍后才開(kāi)始出現(xiàn)特征向量。圖4為重疊尺寸不同時(shí)同名向量占向量總數(shù)的概率變化情況。圖像重疊100個(gè)像素時(shí)同名點(diǎn)概率都在90%以上，但隨著重疊區(qū)域尺寸下降，SIFT算法由于特征向量鄰域較大，同名點(diǎn)概率下降更快。為了保證后續(xù)的匹配質(zhì)量，同名點(diǎn)的概率至少要50%以上，因此SIFT算法重疊區(qū)域保證在最小特征鄰域的1.5倍以上，SURF算法重疊區(qū)域至少為最小鄰域的1.25倍。

(1)

(2)

其中，surf0=2，ssurf0=1.2。

針對(duì)描述向量鄰域影響，部分算法采取在邊界處降低檢測(cè)鄰域大小，但遙感圖像特征相似度高，降低特征鄰域意味著降低算子獨(dú)特性，會(huì)引起誤配率上升，此種方法并不可取。還有一些算法采取插值法增加邊界像素點(diǎn)，但這樣會(huì)增加1倍的計(jì)算量。

綜上所述，在計(jì)算速度方面， SURF算法無(wú)論是64維還是128維其特征向量構(gòu)建速度都遠(yuǎn)遠(yuǎn)大于SIFT算法。在特征向量數(shù)方面，由于特征檢測(cè)鄰域尺寸的影響，部分邊緣檢測(cè)點(diǎn)無(wú)法被描述，為了保證后續(xù)的匹配質(zhì)量，SIFT算法要求重疊區(qū)域在最小特征鄰域的1.5倍以上，而SURF算法重疊區(qū)域至少為最小鄰域的1.25倍。

2.3 特征點(diǎn)匹配

特征點(diǎn)檢測(cè)和描述后，根據(jù)描述向量間距離進(jìn)行特征點(diǎn)匹配，本節(jié)從匹配速度、準(zhǔn)確度、非重疊區(qū)域影響以及匹配距離四方面對(duì)兩種算法進(jìn)行了比較和分析。本文采用歐式距離作為衡量標(biāo)準(zhǔn)，結(jié)合最近鄰比次近鄰的方法匹配特征點(diǎn)，閾值選取0.7。本文選取Mikolajczyk和Schmid提出的匹配準(zhǔn)確度[19]作為判定標(biāo)準(zhǔn)。

(3)

其中，誤配點(diǎn)數(shù)通過(guò)RANSAC(Random Sample Consensus)[20]對(duì)極幾何約束模型算法計(jì)算，迭代次數(shù)設(shè)定為100，總特征點(diǎn)數(shù)為去除誤配前的匹配點(diǎn)數(shù)。

圖5 不同算法匹配效果圖 Fig.5 Matching result of different algorithms

2.3.1 算法匹配速度分析

從表3中可以看出，SIFT算法匹配時(shí)間為SURF算法的3倍左右，主要由于SIFT算法特征向量較多，單位向量的匹配時(shí)間相差并不大。而且SURF算法的不同維度特征向量匹配時(shí)間相差很小。從總時(shí)間來(lái)看，SIFT算法計(jì)算時(shí)間約為SURF算法的3倍，主要時(shí)間差集中在檢測(cè)和描述上，尤其是特征描述步驟。

表3 不同算法特征匹配結(jié)果比較

2.3.2 算法匹配精度分析

從匹配準(zhǔn)確性來(lái)看，SIFT算法的匹配準(zhǔn)確度最好，其次是SURF128維算法，SURF64維算法最差。由于遙感圖像拍攝距離遠(yuǎn)，圖像包含大量相似結(jié)構(gòu)，由于特征描述符僅利用特征點(diǎn)局部鄰域信息，對(duì)散落在這些相似結(jié)構(gòu)中的點(diǎn)極易發(fā)生誤匹配現(xiàn)象，進(jìn)而放大了匹配準(zhǔn)確度方面的差異。

2.3.3 受非重疊區(qū)域影響分析

測(cè)試圖像重疊區(qū)域?yàn)?0個(gè)像素，圖6中為匹配準(zhǔn)確度隨重疊區(qū)域占檢測(cè)區(qū)域比率變化的情況，非重疊區(qū)域的影響反映了算法的獨(dú)特性。由圖6可知不存在非重疊區(qū)域時(shí)3種算法準(zhǔn)確度都在95%以上，隨著非重疊區(qū)域比率的增加算法準(zhǔn)確度都有所下降。當(dāng)重疊區(qū)域占檢測(cè)區(qū)域比率為20%時(shí)，3種算法受到非重疊區(qū)域的影響加大，準(zhǔn)確度開(kāi)始急速下降。相比較而言，SIFT算法和SURF128算法受到非重疊區(qū)域的影響相對(duì)更小。

圖6 非重疊區(qū)域的影響 Fig.6 Influence of non-overlapping domain

2.3.4 特征向量匹配距離分析

從圖7中可以看出，SIFT算法和SURF64算法采用余弦距離匹配準(zhǔn)確度高于歐氏距離，而SURF128算法無(wú)法通過(guò)余弦距離進(jìn)行匹配。原因是余弦距離表示向量間夾角，當(dāng)SURF算法向量由64維增加到128維，將每個(gè)種子點(diǎn)的向量參數(shù)∑dx、 ∑|dx|分成了dy>0，dy<0兩種情況，∑dy、 ∑|dy|同理。由于增加的信息不直接反映方向信息，因此采用余弦距離計(jì)算無(wú)法匹配。

圖7 不同算法計(jì)算距離不同時(shí)正確率的變化情況 Fig.7 Matching ratio change with the calculated distance

綜上所述，從匹配時(shí)間和匹配準(zhǔn)確度兩方面綜合來(lái)看，SURF128算法的匹配性能最好。同時(shí)SURF128算法受非重疊區(qū)域影響相對(duì)較小，但進(jìn)行特征匹配時(shí)只能采用歐式距離。

3 oSURF算法

通過(guò)第二節(jié)分析可知，相對(duì)于SIFT算法和SURF64算法，SURF128算法綜合性能最好，但在特征描述和匹配過(guò)程中，SURF算法過(guò)分依賴局部區(qū)域提取的特征方向，特征方向上小的偏差會(huì)造成匹配準(zhǔn)確度上較大誤差。因此，本節(jié)在SURF128算法基礎(chǔ)上，對(duì)特征方向做出改進(jìn)，提出一種oSURF算法。

3.1 算法改進(jìn)原理

在以興趣點(diǎn)為中心，6s(s為興趣點(diǎn)對(duì)應(yīng)尺度)為半徑的鄰域內(nèi)，用尺寸為4s的haar小波濾波器對(duì)圓形鄰域內(nèi)所有點(diǎn)求x、y方向小波響應(yīng)，再以興趣點(diǎn)為中心對(duì)響應(yīng)進(jìn)行高斯加權(quán)。將圓形鄰域分成若干個(gè)60度扇形，每個(gè)扇形間隔15度，分別計(jì)算每個(gè)扇形中所有x、y方向haar響應(yīng)和，再構(gòu)成的矢量。

圖8 Haar小波濾波器 Fig.8 Haar wavelet filters

圖9 扇形窗口 Fig.9 Sector window

將所有扇形窗口對(duì)應(yīng)矢量按照大小排序，最大矢量的方向即為該特征點(diǎn)主方向。同時(shí)，當(dāng)存在另一個(gè)矢量大小超過(guò)最大矢量的80%時(shí)，則將這個(gè)矢量的方向作為該特征點(diǎn)輔助方向，存儲(chǔ)為新的特征點(diǎn)。雖然輔助方向會(huì)導(dǎo)致特征向量數(shù)目增加，但多個(gè)特征方向可以彌補(bǔ)由于特征方向失準(zhǔn)引起的匹配準(zhǔn)確率下降問(wèn)題。

3.2 算法性能分析

對(duì)oSURF算法進(jìn)行測(cè)試，采用與第2節(jié)相同的實(shí)驗(yàn)數(shù)據(jù)和條件，匹配結(jié)果如下：

圖10 oSURF算法匹配效果圖 Fig.10 Matching result of oSURF algorithm

表4 oSURF算法特征匹配結(jié)果比較

3.2.1 oSURF算法速度分析

由表2、表3和表4對(duì)比可知， 由于oSURF算法增加了輔助方向，特征描述向量相對(duì)SURF128增加了1倍，但總計(jì)算時(shí)間僅增加了0.345 s。綜上所述，oSURF算法保留了SURF算法計(jì)算速度快的優(yōu)點(diǎn)。

3.2.2 oSURF算法匹配準(zhǔn)確度分析

本節(jié)對(duì)oSURF算法與SIFT、SURF64、SURF128算法進(jìn)行匹配準(zhǔn)確度比較，實(shí)驗(yàn)圖像為10組衛(wèi)星條帶圖像，非重疊度在75%～85%之間，檢測(cè)寬度為200個(gè)像素點(diǎn)。圖11為不同算法匹配準(zhǔn)確度比較結(jié)果。

圖11 不同算法匹配結(jié)果 Fig.11 Matching result of different algorithms

通過(guò)2.3節(jié)實(shí)驗(yàn)分析可知，SURF算法雖然比SIFT算法計(jì)算速度提升了3倍，但犧牲了一定的匹配準(zhǔn)確度。由圖11可知，在10組條帶圖像匹配實(shí)驗(yàn)中，oSURF算法相對(duì)于SURF128算法匹配準(zhǔn)確度提升了5%～10%，和SIFT算法匹配準(zhǔn)確度相當(dāng)，恰恰彌補(bǔ)了SURF算法采用框狀濾波器近似帶來(lái)的準(zhǔn)確度損失。

通過(guò)速度和匹配準(zhǔn)確度分析可知，相對(duì)SURF算法，oSURF算法在保證匹配速度的前提下，使得匹配準(zhǔn)確度提高5%～10%。

4 結(jié) 論

本文針對(duì)遙感圖像數(shù)據(jù)，首先，從特征檢測(cè)、描述和匹配三個(gè)方面，對(duì)SIFT算法和SURF算法在重疊率和匹配準(zhǔn)確度方面進(jìn)行了比較；其次，分析了算法對(duì)遙感圖像重疊度、匹配度量距離的要求；最后，通過(guò)對(duì)SURF128算法特征方向的改進(jìn)，提出了oSURF算法，并與SIFT、SURF64、SURF128算法進(jìn)行了匹配準(zhǔn)確度比較。實(shí)驗(yàn)結(jié)果表明，相對(duì)于SIFT算法，SURF算法計(jì)算速度為SIFT的3倍，受非重疊區(qū)域影響更小，為保證50%以上同名點(diǎn)概率，SIFT算法要求圖像重疊區(qū)域不低于其描述向量鄰域大小的1.5倍,SURF算法僅要求為1.25倍；相對(duì)于SURF算法，在保證匹配速度前提下，oSURF算法準(zhǔn)確度提高了5%～10%。因此，oSURF算法更適合遙感圖像拼接。

[1] LOWE D G. Object recognition from local scale-invariant features[C]. International Conference on Computer Vision,Kerkyra,Greece:ICCV,1999:1150-1157 .

[2] LOWE D G. Distinctive image features from scale-invariant keypoints[J].InternationalJ.ComputeVision,2004,60(60):91-110.

[3] BAY H,TUYTELAARS T,GOOL L. SURF:speeded up robust features[J].ComputerVision&ImageUnderstanding,2006,110(3):404-417.

[4] BAY H,ESS A,TUYTELAARS T,etal.. Speeded-Up Robust Features(SURF)[J].ComputerVision&ImageUnderstanding,2008,110(3):346-359.

[5] LUO J,GWUN O. A Comparison of SIFT, PCA-SIFT and SURF[J].J.BusinessEducation,2009,3(4):143-152.

[6] 索春寶，楊東清，劉云鵬.多種角度比較SIFT、SURF、BRISK、ORB、FREAK算法[J].北京測(cè)繪，2014，(4):23-26. SUO CH B,YANG D Q,LIU Y P. Comparison of SIFT, SURF, BRISK, ORB and FREAK algorithms in various aspects[J].BeijingSurveyingandMapping,2014,(4):23-26.(in Chinese)

[7] 安維勝，余讓明，伍玉鈴.基于FAST和SURF的圖像配準(zhǔn)算法[J].計(jì)算機(jī)工程，2015，41(10):232-235. AN W SH,YU R M,WU Y L. Image registration algorithm based on FAST and SURF [J].ComputerEngineering,2015,41(10):232-235.(in Chinese)

[8] 王新華，黃瑋，歐陽(yáng)繼紅.多探測(cè)器拼接成像系統(tǒng)實(shí)時(shí)圖像配準(zhǔn)[J].中國(guó)光學(xué)，2015,8(2):211-219. WANG X H,HUANG W,OUYANG J H. Real-time image registration of the multi-detectors mosaic imaging system[J].ChineseOptics,2015,8(2):211-219.(in Chinese)

[9] 許佳佳.結(jié)合Harris與SIFT算子的圖像快速配準(zhǔn)算法[J].中國(guó)光學(xué)，2015,8(4):574-581. XU J J. Fast image registration method based on Harris and SIFT algorithm[J].ChineseOptics,2015,8(4):574-581.(in Chinese)

[10] 張?jiān)粕?，鄒崢嶸.基于改進(jìn)ORB算法的遙感圖像自動(dòng)配準(zhǔn)方法[J].國(guó)土資源遙感，2013，25(3):20-24. ZAHNG Y SH,ZOU ZH R. Automatic registration method for remote sensing images based on improved ORB algorithm[J].RemoteSensingforLandandResources,2013,25(3):20-24.(in Chinese)

[11] 羅楠，孫權(quán)森，耿蕾蕾,等.一種擴(kuò)展SURF描述符及其在遙感圖像配準(zhǔn)中的應(yīng)用[J].測(cè)繪學(xué)報(bào)，2013，42(3):383-388. LUO N,SUN Q S,GENG L L ,etal.. An extended SURF descriptor and its application in remote sensing images registration[J].ActaGeodaeticaetCartographicaSinica,2013,42(3):383-388.(in Chinese)

[12] 閆利，陳林.一種改進(jìn)的SURF及其在遙感影像匹配中的應(yīng)用[J].武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2013,38(7):770-773. YAN L,CHEN L. A modified SURF descriptor and its application in remote sensing images matching[J].GeomaticsandInformationScienceofWuhanUniversity,2013,38(7):770-773.(in Chinese)

[13] 葛盼盼，陳強(qiáng).基于SURF特征提取的遙感圖像自動(dòng)配準(zhǔn)[J].計(jì)算機(jī)系統(tǒng)應(yīng)用，2014，23(3):16-24. GE P P,CHEN Q. Remote sensing image automatic registration based on SURF feature extraction[J].ComputerSystems&Applications,2014,23(3):16-24.(in Chinese)

[14] LINDEBERG T. Scale-space theory:a basic tool for analyzing structures at different scales[J].J.AppliedStatistics,1994,21(1):225-270.

[15] VIOLA P,JONES M. Rapid object detection using a boosted cascade of simple features[C]. Proceedings of the 2001 IEEE Computer Society Conference on Computer Vision and Pattern Recognition,Kauai,HI,USA,2001:511-518

[16] 翟優(yōu)，曾巒，熊偉.不同局部鄰域劃分加速魯棒特征描述符的性能分析[J].光學(xué) 精密工程，2013，21(9):2395-2404. ZHAI Y,ZENG L,XIONG W. Performance analysis of SURF descriptor with different local region partitions[J].Opt.PrecisionEng.,2013,21(9):2395-2404.(in Chinese)

[17] 孫浩，王程，王潤(rùn)生.局部不變特征綜述[J].中國(guó)圖象圖形學(xué)報(bào)，2011，16(2):141-151. SUN H,WANG CH,WANG R SH. A review of local invariant features[J].J.ImageandGraphics,2011,16(2):141-151.(in Chinese)

[18] 聶海濤,龍科慧,馬軍,等.采用改進(jìn)尺度不變特征變換在多變背景下實(shí)現(xiàn)快速目標(biāo)識(shí)別[J].光學(xué) 精密工程，2015,23(8):2349-2356. NIE H T,LONG K H,MA J,etal.. Fast object recognition under multiple varying background using improved SIFT method[J].Opt.PrecisionEng.,2015,23(8):2349-2356.(in Chinese)

[19] MIKOLAJCZYK K,SCHMID C. A performance evaluation of local descriptors[J].IEEETransactionsonPatternAnalysisandMachineIntelligence,2005,27(10):1615-1630.

[20] FISCHLER M,BOLLES R.RandomSampleConsensus:AParadigmforModelFittingwithApplicationstoImageAnalysisandAutomatedCartography[M]. Burlington:Morgan Kaufmann Publishers Inc.,1987.

《光學(xué) 精密工程》(月刊)

● 中國(guó)光學(xué)開(kāi)拓者之一王大珩院士親自創(chuàng)辦的新中國(guó)歷史最悠久的光學(xué)期刊

● 現(xiàn)任主編為國(guó)家級(jí)有突出貢獻(xiàn)的青年科學(xué)家曹健林博士

● Benjamin J Eggleton, John Love等國(guó)際著名光學(xué)專家為本刊國(guó)際編委

《光學(xué) 精密工程》主要欄目有現(xiàn)代應(yīng)用光學(xué)(空間光學(xué)、纖維光學(xué)、信息光學(xué)、薄膜光學(xué)、光電技術(shù)及器件、光學(xué)工藝及設(shè)備、光電跟蹤與測(cè)量、激光技術(shù)及設(shè)備)；微納技術(shù)與精密機(jī)械(納米光學(xué)、精密機(jī)械)；信息科學(xué)(圖像處理、計(jì)算機(jī)應(yīng)用與軟件工程)等。

? 美國(guó)工程索引EI核心期刊

? 中國(guó)出版政府獎(jiǎng)期刊提名獎(jiǎng)

? 中國(guó)精品科技期刊

? 中文核心期刊

? 百種中國(guó)杰出學(xué)術(shù)期刊

? 中國(guó)最具國(guó)際影響力學(xué)術(shù)期刊

主管單位：中國(guó)科學(xué)院

主辦單位： 中國(guó)科學(xué)院長(zhǎng)春光學(xué)精密機(jī)械與物理研究所

中國(guó)儀器儀表學(xué)會(huì)

地址：長(zhǎng)春市東南湖大路3888號(hào) 郵編：130033

電話：0431-86176855 傳真：0431-84613409

電郵：gxjmgc@sina.com 網(wǎng)址：http://www.eope.net

定價(jià)：100.00元/冊(cè)

Research on matching performance of SIFT and SURF algorithms for high resolution remote sensing image

QI Bing-jie, LIU Jin-guo*, ZHANG Bo-yan, ZUO Yang, LYU Shi-liang

(1.ChangchunInstituteofOptics,FineMechanicsandPhysics，ChineseAcademyofSciences,Changchun130033,China；2.UniversityofChineseAcademyofSciences,Beijing100049,China)

Image matching is the basis of image rectification and mosaic. Because of higher features similarity and smaller overlapped area than ordinary images, the remote sensing images have higher requirements on matching algorithm in both performance and iteration speed. The performances in three aspects: feature detection, feature description and feature matching, are analyzed between the SIFT algorithm and the SURF algorithm in terms of speed and accuracy. The requirements of the degree of overlapping between remote sensing images and the matching distance of the genvector is discussed as well. In view of the characteristic that SURF algorithm is sensitive to the error in feature detection, oSURF algorithm is presented in this paper. Finally, the advantages and disadvantages of each algorithm are analyzed by using satellite remote sensing data of level 1A. The results show that iteration speed of SURF algorithm is three times faster than SIFT algorithm. Under the same matching rate, the width of overlapped area on image required in SURF algorithm is 1.25 times of the dimension of genvector but 1.5 times in SIFT, and the accuracy of oSURF algorithm is increased by 5%～10% compared with SURF algorithm in the same computation speed, which indicate that oSURF is more suitable for remote sensing image stitching.

remote sensing image;feature matching;SIFT;SURF;degree of overlapping

2017-01-18；

2017-03-27

國(guó)家自然科學(xué)基金資助項(xiàng)目(No.61405191) Supported by National Natural Science Foundation of China(No.61405191)

2095-1531(2017)03-0331-09

TP391

A

10.3788/CO.20171003.0331

齊冰潔(1991—)，女，遼寧開(kāi)原人，碩士研究生，主要從事遙感圖像處理方面的研究。E-mail: doudouqbj@126.com

劉金國(guó)(1968—)，男，吉林蛟河人，碩士，研究員，碩士生導(dǎo)師，主要從事空間遙感與圖像處理方面的研究。E-mail: liujg@ciomp.ac.cn

*Correspondingauthor，E-mail:liujg@ciomp.ac.cn