KM-SVM法的SAR圖像無監(jiān)督變化檢測(cè)

田　淞， 宋建社， 張雄美， 任偉龍

(1. 第二炮兵工程大學(xué)七系， 陜西 西安 710025；

2. 重慶通信學(xué)院二系， 重慶 400035)

?

KM-SVM法的SAR圖像無監(jiān)督變化檢測(cè)

田淞1,2， 宋建社1， 張雄美1， 任偉龍1

(1. 第二炮兵工程大學(xué)七系， 陜西 西安 710025；

2. 重慶通信學(xué)院二系， 重慶 400035)

摘要：針對(duì)合成孔徑雷達(dá)(synthetic aperture radar, SAR)圖像變化檢測(cè)獲得有標(biāo)記樣本的數(shù)量十分有限且困難，傳統(tǒng)方法檢測(cè)率低等問題，提出了一種基于原始特征空間的K均值和支持向量機(jī)(K-means and support vector machine, KM-SVM)法SAR圖像無監(jiān)督變化檢測(cè)。首先，不需要任何先驗(yàn)信息的條件下，利用K-means聚類方法獲取差異圖像的分類閾值；其次，利用閾值，引入偏移量，自動(dòng)選取偽訓(xùn)練集和無標(biāo)簽集，并用偽訓(xùn)練集定義SVM的初始決策超平面；最后，用基于統(tǒng)計(jì)特征的半監(jiān)督學(xué)習(xí)算法和支持向量機(jī)相結(jié)合對(duì)圖像進(jìn)行變化類與非變化類的分類。實(shí)驗(yàn)結(jié)果表明：該算法優(yōu)于基于混合高斯分布模型的KI法和基于廣義高斯分布模型的KI法，能保持較好的分類、泛化能力和較穩(wěn)定的檢測(cè)精度。這些結(jié)果表明了文中方法的有效性。

關(guān)鍵詞：合成孔徑雷達(dá)圖像; 變化檢測(cè); 半監(jiān)督學(xué)習(xí); K-均值聚類; 支持向量機(jī)

0引言

變化檢測(cè)是指給定同一地區(qū)或同一目標(biāo)在不同時(shí)間獲得的配準(zhǔn)后的遙感圖像，采用圖像處理的方法檢測(cè)出該地區(qū)(或目標(biāo))地物的變化情況，通過對(duì)不同時(shí)相遙感圖像的綜合分析獲取變化信息，是實(shí)現(xiàn)遙感監(jiān)測(cè)的重要手段。合成孔徑雷達(dá)(synthetic aperture radar，SAR)因具有全天候、全天時(shí)、寬幅成像、多極化、可變觀測(cè)角、穿透能力強(qiáng)和高分辨率等特點(diǎn)，SAR圖像變化檢測(cè)已成為當(dāng)前遙感監(jiān)測(cè)的前沿技術(shù)之一。目前，SAR圖像變化檢測(cè)技術(shù)已廣泛應(yīng)用于自然災(zāi)害監(jiān)測(cè)、植被和森林變化、地面目標(biāo)監(jiān)視[1-7]等領(lǐng)域。

SAR圖像變化檢測(cè)的典型方法有：比較后分類，分類后比較和聯(lián)合分類，其中比較后分類方法具有簡(jiǎn)單、直觀、易于理解等優(yōu)點(diǎn)被廣泛關(guān)注，如基于混合高斯分布模型的EM法[8]、基于混合高斯分布模型的KI法(KI threshold selection criterion based on Gaussian model, GM-KI)、基于廣義高斯分布模型的EM法[9]、基于廣義高斯分布模型的KI法(KI threshold selection criterion based on general Gaussian model, GGM-KI)[10]、基于廣義Gamma模型的KI法[11]，還有采用瑞利-高斯模型[12]、非高斯分布模型[13]方法等。這些方法通常需要假定差異圖像上變化類與無變化類灰度值的統(tǒng)計(jì)分布模型，統(tǒng)計(jì)分布不同，選取的閾值也就不同。當(dāng)假設(shè)的統(tǒng)計(jì)分布模型能恰當(dāng)?shù)姆从砈AR圖像不同類別分布時(shí)，這些方法就能有效提取出變化信息；當(dāng)圖像不同類別差異較小或圖像受噪聲影響嚴(yán)重，使這些先驗(yàn)信息不恰當(dāng)時(shí)，這些方法通常難以獲得滿意的效果，特別是當(dāng)變化類和無變化類的灰度值分布是相互交疊時(shí)，基本不能實(shí)現(xiàn)變化信息的有效獲取。

基于上述分析，提出了K均值和支持向量機(jī)(K-mean and support vector machine, KM-SVM)法的SAR圖像無監(jiān)督變化檢測(cè)方法。該方法利用K-means聚類算法獲取差異圖像在原始特征空間上的偽訓(xùn)練集和無標(biāo)簽集；在此基礎(chǔ)上，利用偽訓(xùn)練集和無標(biāo)簽集數(shù)據(jù)信息共同訓(xùn)練SVM[14-15]后，用SVM完成變化檢測(cè)，無需任何先驗(yàn)信息實(shí)現(xiàn)最優(yōu)變化檢測(cè)的目的。

1基于KM-SVM的SAR圖像變化檢測(cè)方法

設(shè)I1和I2分別為t1，t2時(shí)刻，同一區(qū)域，經(jīng)配準(zhǔn)后的SAR圖像，大小均為P×Q。分別對(duì)圖像I1和I2進(jìn)行預(yù)處理，預(yù)處理包括對(duì)SAR圖像的輻射、幾何定標(biāo)等處理。ωu，ωc分別為圖像像素?zé)o變化類和變化類集合，即Ω={ωu,ωc}。算法分為3步：①對(duì)數(shù)比值差異圖構(gòu)造；②K-means聚類算法獲取偽訓(xùn)練集和無標(biāo)簽集；③使用混合樣本訓(xùn)練SVM，利用無標(biāo)簽集樣本信息，得到一個(gè)泛化性能更好的分類器。

1.1對(duì)數(shù)比值差異圖構(gòu)造

差異圖構(gòu)造應(yīng)使Ω中的2類差異增強(qiáng)，便于檢測(cè)分類。SAR圖像中存在相干斑噪(specklenoise，SN)，它是由SAR的工作機(jī)制造成的，為了減低其對(duì)后續(xù)處理的影響。本文采用對(duì)數(shù)比值算子構(gòu)造差異圖像XDI，即

(1)

采用對(duì)數(shù)比值算子構(gòu)造差異圖有3個(gè)主要優(yōu)點(diǎn)：一是比值差異圖像的分布僅與圖像I1和I2強(qiáng)度的相對(duì)變化有關(guān)，與像素強(qiáng)度的大小無關(guān)。無論區(qū)域內(nèi)像素的強(qiáng)度強(qiáng)或弱，變化檢測(cè)的尺度是相同的；二是比值算子具有更強(qiáng)的魯棒性。因?yàn)镾AR成像過程中產(chǎn)生的乘性輻射誤差，能夠在重復(fù)軌道成像中精確重現(xiàn)，比值算子可以消去這種誤差。三是對(duì)比值算子取對(duì)數(shù)，可以減小比值差異圖像的波動(dòng)范圍，便于變化檢測(cè)處理。

1.2偽訓(xùn)練集和無標(biāo)簽集的生成

在SAR圖像變化檢測(cè)過程中，先驗(yàn)信息的獲取是困難的；另外由于相干斑噪的影響，獲得的先驗(yàn)信息存在較大誤差，而聚類分析是一種無監(jiān)督方法，不需要任何先驗(yàn)信息。K-means是一種簡(jiǎn)單、常用的聚類方法。首先選擇2個(gè)樣本作為初始聚類中心，按照聚類域中所有樣本到聚類中心的距離平方和最小的原則，使其他樣本向各個(gè)中心聚集，然后判斷分類是否合理，不合理則繼續(xù)迭代，直到聚類中心不再發(fā)生改變?yōu)橹埂５玫?個(gè)聚類中心值{k1,k2}，從而獲得閾值T(即T=[max{k1,k2}+min{k1,k2}]/2)。T是無變化類ωu和變化類ωc的粗分類閾值，靠近閾值T附近的像素處在2種類別的臨界狀態(tài)，受到誤差影響，一般分類精度較低，直接分類效果較差。但閾值T是正確定義偽訓(xùn)練集Xp的一個(gè)參考點(diǎn)，如圖1所示，首先把不確定其標(biāo)簽的區(qū)域定義為無標(biāo)簽集Xu，有

(2)

同時(shí)，具有高概率屬于ωu和ωc的樣本定義為偽訓(xùn)練集Xp,有

(3)

式中，ε是用于確定偽訓(xùn)練集和無標(biāo)簽集的邊界調(diào)整系數(shù)。ε值越大，Xp中樣本數(shù)量越少，可靠性越高；反之，則Xp中樣本數(shù)量越多，可靠性越低。ε的設(shè)置應(yīng)保證Xp中樣本獲得正確標(biāo)簽的概率較高。偽訓(xùn)練集中無變化類和變化類標(biāo)簽分配如下：

(4)

一般情況下，ε的取值根據(jù)差異圖像的灰度值范圍確定。也可以定義更為復(fù)雜的ε取值方法，但考慮到計(jì)算復(fù)雜度和對(duì)最終變化檢測(cè)結(jié)果影響，根據(jù)經(jīng)驗(yàn)一般取ε=0.5。

圖1　偽訓(xùn)練集和無標(biāo)簽集生成示意圖

1.3SVM的半監(jiān)督學(xué)習(xí)分類檢測(cè)

這一部分的主要思想是在多時(shí)相SAR圖像的原始特征空間中定義決策函數(shù)，把無變化類和變化類樣本精確分類。偽訓(xùn)練集樣本的選取是基于平衡折中方法，不能完整表述無變化類和變化類的統(tǒng)計(jì)特征，因此利用無標(biāo)簽樣本信息可以進(jìn)一步提高SVM的分類和泛化能力[16]。

核函數(shù)選擇——在支持向量機(jī)中核函數(shù)起著非常重要的作用，它能夠解決非線性問題和維數(shù)災(zāi)難問題。核函數(shù)的作用是代替高維特征空間中的內(nèi)積運(yùn)算，從而避免復(fù)雜的高維運(yùn)算。統(tǒng)計(jì)學(xué)習(xí)理論指出，只要對(duì)稱函數(shù)K(x,y)滿足Mercer條件就能作為核函數(shù)[17]。當(dāng)缺乏先驗(yàn)信息或先驗(yàn)信息不可靠時(shí)，最好的方法是選用球面核[18]，即高斯徑向基函數(shù)核：

初始化—利用標(biāo)準(zhǔn)SVM作為代價(jià)函數(shù)和偽訓(xùn)練集Xp(包含偽變化類數(shù)據(jù)和偽無變化類數(shù)據(jù))獲得初始分類超平面。

(5)

式中，C1是正則化參數(shù)；φ(·)是非線性映射函數(shù)；n表示偽訓(xùn)練樣本數(shù)量；為處理不可分樣本，引入偽訓(xùn)練樣本的松弛變量ξp；根據(jù)標(biāo)準(zhǔn)SVM的符號(hào)要求，ωu和ωc的類標(biāo)簽分別用“-1”和“+1”表示，即yp∈{ωu,ωc}≡{-1,+1}。

半監(jiān)督學(xué)習(xí)—無標(biāo)簽集樣本通過半監(jiān)督學(xué)習(xí)后，SVM使用偽訓(xùn)練集和無標(biāo)簽集組成的混合樣本進(jìn)行訓(xùn)練，得到一個(gè)分類和泛化性能更好的分類器。文中采用的半監(jiān)督學(xué)習(xí)SVM是利用無標(biāo)簽樣本的均值特征，替代估計(jì)每一個(gè)無標(biāo)簽樣本標(biāo)簽，得到最優(yōu)分類超平面。

半監(jiān)督SVM的優(yōu)化問題[19]為

(6)

式中，ξp和ξu為偽訓(xùn)練樣本和無標(biāo)簽樣本的松弛變量，以及相應(yīng)的懲罰系數(shù)C1和C2；n表示偽訓(xùn)練樣本數(shù)量;m表示參與訓(xùn)練的無標(biāo)簽樣本數(shù)量。為了構(gòu)造無標(biāo)簽集Xu中的2類樣本均值統(tǒng)計(jì)特征，則

(7)

(8)

即

(9)

因此,式(7)和式(6)是等價(jià)的。由于

(10)

得

(11)

(12)

(13)

式中，α=[α1,…，αN,αN+1,αN+2]T為L(zhǎng)agrange乘子，y=[y1,…，yN,1,-1]T為類標(biāo)簽集，Kij=(φi)T(φj)為核函數(shù)，其中非線性映射函數(shù)為

使用交替優(yōu)化算法求解式(13)，得到最終分類超平面。

2實(shí)驗(yàn)結(jié)果及分析

為了驗(yàn)證上述算法的有效性，進(jìn)行了相關(guān)實(shí)驗(yàn)，這里給出2組實(shí)驗(yàn)，驗(yàn)證本文算法的有效性。

2.1數(shù)據(jù)描述

實(shí)驗(yàn)1的SAR圖像如圖2所示。變化檢測(cè)原始圖像分別為1997年5月和8月加拿大渥太華地區(qū)的Radarsat-1SAR影像，如圖2(a)和2(b)所示。圖像尺寸均為350×290，灰度級(jí)為256，配準(zhǔn)誤差為1.5個(gè)像素左右。前后兩時(shí)相的變化信息主要是因?yàn)?月雨季來臨，洪水淹沒了部分陸地所致。圖2(c)為變化檢測(cè)參考圖，圖中的白色像素表示兩時(shí)相間發(fā)生變化的區(qū)域，像素個(gè)數(shù)為16 049。

圖2　渥太華地區(qū)水災(zāi)SAR圖像

實(shí)驗(yàn)2的SAR 圖像如圖3所示。變化檢測(cè)原始圖像是分別在1999年4月和5月瑞士Bern城市ERS-2的SAR圖像。前一時(shí)相的SAR圖像是在水災(zāi)發(fā)生前獲取的，如圖3(a)所示；后一時(shí)相的SAR圖像是在水災(zāi)發(fā)生后獲得的，圖像中陰暗部分為洪水淹沒區(qū)域，如圖3(b)所示。圖像大小均為 301×301像素，灰度級(jí)為256。圖3(c)為變化檢測(cè)參考圖，圖中白色像素表示兩時(shí)相間發(fā)生變化區(qū)域，像素個(gè)數(shù)為1 155。

其中支持向量機(jī)優(yōu)化后的參數(shù)選擇為C1=100，C2=0.1，高斯核寬度2σ2∈(0.1,1)。

圖3　Bern城市水災(zāi)SAR圖像

2.2實(shí)驗(yàn)結(jié)果

實(shí)驗(yàn) 1對(duì)數(shù)比算子生成的差異圖XDI，max(XDI)為4.060 4，min(XDI)為0。根據(jù)K-means聚類得到的閾值T=1.104 7，取ε為0.5，按照規(guī)則偽訓(xùn)練集Xp中xp≤0.552 35屬于無變化類ωu共71 457個(gè)樣本，xp≥1.657 05屬于變化類ωc共8 784個(gè)樣本，無標(biāo)簽集Xu樣本為21 259個(gè)。SVM是小樣本情況下的機(jī)器學(xué)習(xí)，其分類超平面僅與支持向量有關(guān)，因此對(duì)偽訓(xùn)練集和無標(biāo)簽集重采樣，測(cè)試KM-SVM方法的分類和泛化能力，以及利用無標(biāo)簽樣本的開放式半監(jiān)督學(xué)習(xí)能力。在學(xué)習(xí)過程中充分利用樣本均值統(tǒng)計(jì)特征，提高學(xué)習(xí)效率；并與GM-KI和GGM-KI法的檢測(cè)結(jié)果進(jìn)行比較，如表1和圖4所示。

表1　渥太華地區(qū)水災(zāi)圖像變化檢測(cè)結(jié)果分析

圖4　渥太華地區(qū)變化檢測(cè)結(jié)果

實(shí)驗(yàn) 2對(duì)數(shù)比算子生成的差異圖XDI，max(XDI)為5.347 1，min(XDI)為0。根據(jù)K-means聚類得到的閾值T=1.649 1，取ε為0.5，按照規(guī)則偽訓(xùn)練集Xp中xp≤0.824 55屬于無變化類ωu共87 239個(gè)樣本，xp≥2.473 65屬于變化類ωc共624個(gè)樣本，無標(biāo)簽集Xu樣本為2 738個(gè)。實(shí)驗(yàn)結(jié)果如表2和圖5所示。

表2　Bern城市水災(zāi)圖像變化檢測(cè)結(jié)果分析

圖5　Bern城市變化檢測(cè)結(jié)果

2.3實(shí)驗(yàn)分析

利用典型變化檢測(cè)算法，包括GGM-KI算法和GM-KI算法，與KM-SVM算法相比較，表1和表2分別給出檢測(cè)結(jié)果。從表中數(shù)據(jù)可以得出以下幾點(diǎn)結(jié)論。

(1) GGM-KI算法要優(yōu)于GM-KI算法，這是因?yàn)閺V義高斯統(tǒng)計(jì)分布模型的適用范圍更廣。

(2) 對(duì)于GGM-KI算法和GM-KI算法，表1和圖4顯示的變化檢測(cè)結(jié)果較好，而表2和圖5的結(jié)果相對(duì)較差，是因?yàn)镚GM和GM模型與渥太華數(shù)據(jù)集數(shù)據(jù)模型更匹配。充分說明基于統(tǒng)計(jì)分布模型的檢測(cè)算法結(jié)果受統(tǒng)計(jì)分布模型選擇的影響較大。

(3) KM-SVM算法得到更優(yōu)的結(jié)果，并且檢測(cè)結(jié)果相對(duì)穩(wěn)定，這是因?yàn)镵M-SVM算法不需要模型選擇，并且SVM支持非線性分類，對(duì)于不同類別相互交疊的情況，也有較好的分類能力。

(4) 在2組實(shí)驗(yàn)中，KM-SVM算法的檢測(cè)精度分別提高了8%和13%，是因?yàn)槔昧藷o標(biāo)簽集的數(shù)據(jù)信息，能夠進(jìn)一步優(yōu)化分類決策超平面。

(5) 在運(yùn)算時(shí)間上，KM-SVM算法明顯比GGM-KI算法和GM-KI算法要長(zhǎng)，主要是半監(jiān)督學(xué)習(xí)過程運(yùn)算量較大，時(shí)間耗費(fèi)較多。

3結(jié)論

本文提出的KM-SVM法的SAR圖像變化檢測(cè)，通過半監(jiān)督學(xué)習(xí)，提取無標(biāo)簽數(shù)據(jù)信息，可以獲得較穩(wěn)定、較好的變化檢測(cè)結(jié)果，不會(huì)隨著SAR圖像數(shù)據(jù)源的變化而出現(xiàn)較大波動(dòng)。該方法是一種無監(jiān)督的方法，首先用K-means算法對(duì)對(duì)數(shù)比值差異圖進(jìn)行聚類，得到粗分類的閾值，對(duì)該閾值加上松弛量將像素點(diǎn)分成有標(biāo)簽的偽訓(xùn)練樣本和無標(biāo)簽樣本；然后，聯(lián)合利用偽訓(xùn)練樣本和無標(biāo)簽樣本訓(xùn)練SVM，此時(shí)，是用無標(biāo)簽樣本均值特征來修正由偽訓(xùn)練樣本訓(xùn)練得到的初始分類面，得到更優(yōu)的分類超平面，以提升檢測(cè)效果。實(shí)驗(yàn)結(jié)果證明，該方法可行、有效。下一步工作是利用圖像空間上下文信息進(jìn)一步提高變化檢測(cè)精度。

參考文獻(xiàn)：

[1] Bouaraba A, Younsi A, Belhadj A A, et al. Robust techniques for coherent change detection using cosmo-skymed SAR images[J].ProgressinElectromagneticsResearchM, 2012, 22:219-232.

[2] Lincy P, Dr Ramamoorthy P. Synthetic aperture radar image change detection using fuzzy C-means clustering algorithm[J].InternationalJournalofAdvancedResearchinComputerandCommunicationEngineering, 2013, 2(3):1374-1379.

[3] Meghana P M K, Ravikumar A V. Robust technique for change detection in multitemporal synthetic aperture radar images[J].InternationalJournalofScientific&EngineeringResearch, 2013, 4(7):2228-2235.

[4] Radke R J, Andra S, Al-Kofahi O, et al. image change detection algorithms:a systematic survey[J].IEEETrans.onImageProcessing, 2005, 14(3):294-307.

[5] Fransson J E S, Walter F, Blennow K, et al. Detection of storm-damaged forested areas using airborne CARABAS-II VHF SAR image data[J].IEEETrans.onGeoscienceandRemoteSensing, 2002, 40(10):2170-2175.

[6] Bovolo F, Marin C, Bruzzone L. A hierarchical approach to change detection in very high resolution SAR images for surveillance applications[J].IEEETrans.onGeoscienceandRemoteSensing,2013,51(4):2042-2054.

[7] Huang S Q, Liu D Z, Hu M X, et al. Multi-temporal SAR image change detection technique based on wavelet transform[J].ActaGeodaeticaetCartographicaSinica, 2010, 39(2):180-186.(黃世奇, 劉代志, 胡明星，等. 基于小波變換的多時(shí)相SAR圖像變化檢測(cè)技術(shù)[J].測(cè)繪學(xué)報(bào), 2010, 39(2):180-186.)

[8] Bruzzone L, Fernandez P D. Automatic analysis of the difference image for unsupervised change detection[J].IEEETrans.onGeoscienceandRemoteSensing, 2000, 38(3):1171-1182.

[9] Bazi Y, Bruzzone L, Melgani F. Image thresholding based on the EM algorithm and the generalized distribution[J].PatternRecognition, 2007, 40(2):619-634.

[10] Hu Z L. An unsupervised change detection approach based on KI dual thresholds under the generalized Gauss model assumption in SAR images[J].ActaGeodaeticaetCartographicaSinica, 2013, 42(1):116-122.(胡召玲.廣義高斯模型及KI雙閾值法的SAR圖像非監(jiān)督變化檢測(cè)[J].測(cè)繪學(xué)報(bào),2013,42(1):116-122.)

[11] Gao C S, Zhang H, Wang C, et al. SAR change detection based on generalized Gamma distribution divergence and auto-threshold segmentation[J].JournalofRemoteSensing, 2010, 14(4):718-724.(高叢珊, 張紅, 王超，等．廣義Gamma模型及自適應(yīng)KI閾值分割的SAR圖像變化檢測(cè)[J]., 遙感學(xué)報(bào), 2010, 14(4):718-724.)

[12] Wang G T, Wang Y L, Jiao L C. Adaptive spatial neighborhood analysis and rayleigh-Gaussian distribution fitting for change detection in multi-temporal remote sensing images[J].JournalofRemoteSensing, 2009, 13(4):639-646.

[13] Moser G, Serpico S B. Generalized minimum-error thresholding for unsupervised change detection from SAR amplitude imagery[J].IEEETrans.onGeoscienceandRemoteSensing, 2006, 44(10):2972-2983.

[14] Zhang Y Y, Liu M H, Han T. Elevator-waiting people counting system based on mean shift segmentation and SVM classification[J].OpticsandPrecisionEngineering, 2013, 21(4):1079-1085.(張宇洋，劉滿華，韓韜. 基于Mean Shift圖像分割和支持向量機(jī)判決的候梯人數(shù)視覺檢測(cè)系統(tǒng)[J].光學(xué)精密工程, 2013, 21(4):1079-1085.)

[15] Zhou Z H . Disagreement-based semi-supervised learning[J].ActaAutomaticaSinica, 2013, 39(11):1871-1878.(周志華. 基于分歧的半監(jiān)督學(xué)習(xí)[J].自動(dòng)化學(xué)報(bào),2013,39(11):1871-1878.)

[16] Bruzzone L, Chi M, Marconcini M. A novel transductive SVM for semisupervised classification of remote-sensing images[J].IEEETrans.onGeoscienceandRemoteSensing, 2006, 44(11):3363-3373.

[17] Burges C J. A tutorial on support vector machines for pattern recognition[J].KnowledgeDiscoveryandDataMining, 1998, 2(2):121-167.

[18] Chapelle O, Vapnik V, Bousquet O, et al. Choosing multiple parameters for support vector machines[J].MachineLearning, 2002, 46(1):131-159.

[19] Bovolo F, Bruzzone L, Marconcini M. A novel approach to unsupervised change detection based on a semisupervised SVM and a similarity measure[J].IEEETrans.onGeoscienceandRemoteSensing, 2008, 46(7):2070-2082.

[20] Li Y F, Tsang I W, Kwok J T, et al. Tighter and convex maximum margin clustering[C]∥Proc.ofthe12thInternationalConferenceonArtificialIntelligenceandStatistics, 2009:344-351.

田淞(1975-)，男，講師，博士研究生，主要研究方向?yàn)檫b感圖像處理、模式識(shí)別。

E-mail：cqtianyi423@sina.com

宋建社(1954-)，男，教授，博士研究生導(dǎo)師，主要研究方向?yàn)楹铣煽讖嚼走_(dá)信號(hào)處理、圖像處理、信息建模和系統(tǒng)工程。

E-mail：songjshe@126.com

張雄美(1983-)，女，講師，博士，主要研究方向?yàn)檫b感圖像處理。

E-mail：zxw.ok@163.com

任偉龍(1988-)，男，博士研究生，主要研究方向?yàn)楹铣煽讖嚼走_(dá)圖像處理。

E-mail：renweilong350@163.com

網(wǎng)絡(luò)優(yōu)先出版地址：http://www.cnki.net/kcms/detail/11.2422.TN.20141121.0934.007.html

KM-SVM approach to unsupervised change detection in SAR images

TIAN Song1,2， SONG Jian-she1， ZHANG Xiong-mei1， REN Wei-long1

(1. 7thDepartment,SecondArtilleryEngineeringUniversity,Xi’an710025,China;

2. 2ndDepartment,ChongqingCommunicationInstitute,Chongqing400035,China)

Abstract:It is difficult to obtain training data in practical synthetic aperture radar (SAR)image change detection tasks, and the detection rates of classical methods are low. A K-means and support vector machine (KM-SVM) approach, which aims at extracting the change information in the original feature space without any training data, is proposed. First, the threshold of difference images is selected by the K-means clustering method. Second, an offset should be obtained for deriving an unlabeled set and a pseudo-training set that is necessary for initializing a binary support vector machine (SVM) classifier. Finally, according to a semi-supervised learning algorithm based on statistical characteristics, the SVM performs change detection by considering unlabeled data in definition of the decision boundary between unchanged and changed pixels. Experiment results show that, the proposed approach achieves better performance, higher generalization ability and more stable change detection precision, than the classical GM-KI, GGM-KI methods. These results prove the efficiency of the proposed approach.

Keywords:synthetic aperture radar (SAR)images; change detection; semi-supervised learning; K-means; support vector machine (SVM)

作者簡(jiǎn)介：

中圖分類號(hào)：TP 751

文獻(xiàn)標(biāo)志碼：ADOI:10.3969/j.issn.1001-506X.2015.05.10

基金項(xiàng)目：國家自然科學(xué)基金(61072141,61132008)資助課題

收稿日期：2014-05-04；修回日期：2014-10-23；網(wǎng)絡(luò)優(yōu)先出版日期：2014-11-21。