自適應(yīng)距離和模糊拓?fù)鋬?yōu)化的模糊聚類SAR影像變化檢測(cè)

王建明，史文中，邵 攀,

1. 武漢大學(xué)遙感信息工程學(xué)院，湖北 武漢 430079； 2. 香港理工大學(xué)土地測(cè)量及地理資訊學(xué)系，香港 九龍

遙感變化檢測(cè)是對(duì)地觀測(cè)應(yīng)用中的關(guān)鍵技術(shù)[1]。因合成孔徑雷達(dá)(synthetic aperture radar，SAR)具有全天時(shí)、全天候和能夠穿透某些地物等成像優(yōu)勢(shì)[2]，SAR影像變化檢測(cè)技術(shù)受到了越來越多的關(guān)注[1-7]。

最常用的變化檢測(cè)一般基于“差分影像”[8]，該方法首先通過對(duì)比兩期影像生成差分影像，然后通過分析差分影像獲取地表變化信息。本文的重點(diǎn)是分析差分影像，其本質(zhì)是一個(gè)二分類問題[9]，即將差分影像劃分為未變化和變化兩類。

模糊C均值[10](fuzzy C-means，F(xiàn)CM)具有無須建立模型、性能穩(wěn)定以及對(duì)數(shù)據(jù)不確定性建模等優(yōu)勢(shì)[11]，能夠?qū)崿F(xiàn)對(duì)差分影像的有效分割，在遙感變化檢測(cè)領(lǐng)域有著廣泛的應(yīng)用[5,9,12-17]。文獻(xiàn)[13]利用FCM通過TM(thematic mapper)影像進(jìn)行火災(zāi)和湖水變化檢測(cè),取得了良好的檢測(cè)效果。但因FCM未考慮空間信息，其對(duì)噪聲比較敏感，為此，許多學(xué)者從不同角度將空間信息引入FCM中，提出許多改進(jìn)算法[5,11]。

為提高FCM的運(yùn)行效率和對(duì)噪聲的穩(wěn)健性，文獻(xiàn)[18]提出上下文敏感的快速模糊C均值算法(fast generalized FCM，F(xiàn)GFCM)。文獻(xiàn)[16]將FGFCM引入變化檢測(cè)中，提高了FCM對(duì)較小變化區(qū)域的檢測(cè)性能。為確保對(duì)噪聲不敏感同時(shí)，保留更多細(xì)節(jié)信息，文獻(xiàn)[11]通過定義一個(gè)全新的模糊因子Gkn增強(qiáng)了FCM的目標(biāo)函數(shù)，提出模糊局部信息C均值聚類算法(fuzzy local information C-means，F(xiàn)LICM)。文獻(xiàn)[14]將FLICM應(yīng)用到SAR影像變化檢測(cè)中，取得了良好的檢測(cè)效果。針對(duì)Gkn的不足，文獻(xiàn)[5]對(duì)其進(jìn)行優(yōu)化，提出增強(qiáng)的FLICM算法(reformulated FLICM，RFLICM)，實(shí)現(xiàn)了對(duì)差分影像空間信息的更有效利用。文獻(xiàn)[9]通過馬爾可夫隨機(jī)場(chǎng)模型將空間鄰域信息引入FCM中，增強(qiáng)了FCM的聚類效果，提高了變化檢測(cè)精度。文獻(xiàn)[12]提出穩(wěn)健半監(jiān)督模糊C均值(robust semi-supervised FCM，RSFCM)變化檢測(cè)技術(shù)，RSFCM同時(shí)考慮了差分影像的灰度值、類別知識(shí)和空間鄰域信息，提高了FCM聚類的變化檢測(cè)精度。

FCM及其變體都能實(shí)現(xiàn)對(duì)差分影像的有效分割，但仍存在兩個(gè)共有缺陷，影響了變化檢測(cè)精度的進(jìn)一步提高：①通常采用歐氏距離計(jì)算樣本點(diǎn)到聚類中心的距離，忽略了未變化類與變化類的形狀信息，影響了隸屬度函數(shù)的準(zhǔn)確程度；②通常采用最大隸屬度原則對(duì)模糊隸屬度函數(shù)進(jìn)行去模糊化操作[11]，但該原則未考慮像元隸屬度的分布差異：例如若像元A和B屬于變化類的隸屬度分別為0.98和0.52，屬于未變化類的隸屬度分別為0.02和0.48，則根據(jù)該原則，A和B都會(huì)被標(biāo)記為變化類，但像元B屬于變化和未變化類的隸屬度非常接近，其分類結(jié)果具有很強(qiáng)的不確定性。

為克服FCM及其變體的上述缺陷，提出一種基于自適應(yīng)距離(adaptive distance，AD)[19]和模糊拓?fù)?fuzzy topology，F(xiàn)T)[20-21]理論的SAR影像變化檢測(cè)技術(shù)框架(adaptive distance and fuzzy topology based change detection,FATCD)：FATCD首先利用自適應(yīng)距離改進(jìn)模糊聚類，再利用改進(jìn)后算法對(duì)差分影像進(jìn)行軟分類，獲得更加準(zhǔn)確的模糊隸屬度函數(shù)；而后根據(jù)像元隸屬度的分布差異，利用模糊拓?fù)鋵⒉罘钟跋駝澐譃椴淮_定性較低的內(nèi)部和具有較高不確定性的模糊邊界區(qū)域；最后分別利用最大隸屬度原則和模糊拓?fù)涞闹芜B通性對(duì)內(nèi)部和模糊邊界像元分配類別。

1 FATCD變化檢測(cè)

假設(shè)X1和X2是同一地區(qū)不同時(shí)刻的，已經(jīng)過校正和配準(zhǔn)的兩幅SAR影像，采用對(duì)乘性相干斑噪聲具有穩(wěn)健性的“對(duì)數(shù)比”算子[1]來生成差分影像XD

XD(i,j)=|log(X1(i,j)/X2(i,j))|=

|log(X1(i,j))-log(X2(i,j))|

(1)

式中，XD(i,j)、X1(i,j)和X2(i,j)分別表示相應(yīng)影像上第i行、第j列的像元。得到差分影像后，利用所提出的FATCD(參見圖1)將其劃分為未變化和變化兩類。分別用符號(hào)wu和wc表示未變化類和變化類集合。

FATCD技術(shù)主要包括兩個(gè)階段：①計(jì)算模糊隸屬度函數(shù)：首先利用自適應(yīng)距離增強(qiáng)模糊聚類算法，使其能夠自適應(yīng)差分影像的統(tǒng)計(jì)特性，從而得到未變化類wu和變化類wc更加準(zhǔn)確的模糊隸屬度函數(shù)Mu和Mc。②生成變化檢測(cè)圖：根據(jù)Mu和Mc的取值不同，利用模糊拓?fù)鋵⒉罘钟跋駝澐譃榫哂休^高隸屬度的內(nèi)部和隸屬度較低的模糊邊界兩部分：內(nèi)部像元隸屬度較高，直接利用最大隸屬度原則判斷其類別；而邊界像元具有很強(qiáng)的不確定性，利用模糊拓?fù)渲芜B通性對(duì)其分配類別。

FATCD是一個(gè)優(yōu)化模糊聚類變化檢測(cè)的技術(shù)框架，在步驟①中可采用不同的FCM類型算法，如引言中提到的模糊聚類算法，給定其中一個(gè)，F(xiàn)ATCD就能夠誘導(dǎo)出一個(gè)相應(yīng)的基于模糊拓?fù)浜妥赃m應(yīng)距離的模糊聚類SAR影像變化檢測(cè)算法。本文以FCM和FLICM為例來研究和探討所提出的FATCD技術(shù)框架，給出兩種誘導(dǎo)算法FatFCM和FatFLICM。

1.1 計(jì)算模糊隸屬度函數(shù)

FCM和FLICM通過優(yōu)化目標(biāo)函數(shù)來求解每個(gè)樣本點(diǎn)到所有聚類中心的最優(yōu)隸屬度，從而實(shí)現(xiàn)對(duì)數(shù)據(jù)集的模糊劃分。公式(2)和(3)分別給出了FCM和FLICM的目標(biāo)函數(shù)

(2)

(3)

式中，m為加權(quán)指數(shù)，用來控制隸屬度函數(shù)的模糊度，本文中m=2；N為待聚類數(shù)據(jù)集的樣本個(gè)數(shù)，本文中指差分影像的像元個(gè)數(shù)；C為類別數(shù)目，本文中包括變化和未變化兩類；vk為第k個(gè)聚類中心；yn為第n個(gè)樣本；ukn為樣本yn屬于第k個(gè)類別的隸屬度；d(yn,vk)為樣本yn到聚類中心vk的歐氏距離；Gkn是FLICM聚類為利用空間信息引入的模糊因子[11]

(4)

式中，N(yn)是樣本yn的空間鄰域；yj是yn的鄰域像元；dnj是樣本yj和yn的空間距離；d(yj,vk)是樣本yj到聚類中心vk的歐氏距離。

從公式(2)、(3)和(4)可知，樣本點(diǎn)到聚類中心的距離公式d(yn,vk)是FCM和FLICM目標(biāo)函數(shù)的重要組成部分；因此其性能很大程度上決定著最終模糊隸屬度函數(shù)的準(zhǔn)確性。包括FCM和FLICM在內(nèi)的現(xiàn)有模糊聚類變化檢測(cè)中，通常采用歐氏距離d(yn,vk)=‖yn-vk‖來計(jì)算樣本點(diǎn)到聚類中心的距離。

歐氏距離為樣本點(diǎn)間的直線距離，易于理解且計(jì)算方便，是使用最廣泛的距離度量。然而歐氏距離沒有充分考慮不同類別之間的分布差異，不能根據(jù)差分影像的統(tǒng)計(jì)特性自適應(yīng)優(yōu)化距離公式，從而使得FCM和FLICM在分析差分影像過程中，沒有充分利用未變化和變化類數(shù)據(jù)的形狀信息，一定程度上影響了隸屬度函數(shù)的準(zhǔn)確程度。為此，F(xiàn)ATCD借助自適應(yīng)距離來增強(qiáng)模糊聚類變化檢測(cè)。根據(jù)文獻(xiàn)[19]可知，自適應(yīng)距離公式有如下通用形式

(5)

式中，Ak為類別wk的距離誘導(dǎo)矩陣，k∈{u,c}。選取不同形式的矩陣Ak可得到不同的自適應(yīng)距離公式：若Ak為對(duì)角矩陣，則dAk(yn,vk)為加權(quán)歐氏距離；若Ak為類別wk協(xié)方差矩陣的逆矩陣，則dAk(yn,vk)為馬氏距離。

在本研究中，差分影像數(shù)據(jù)只包含一個(gè)屬性分量，使得未變化類和變化類的距離誘導(dǎo)矩陣Ak相對(duì)簡(jiǎn)單：Ak簡(jiǎn)化成一維矩陣，即一個(gè)常數(shù)。故無論Ak的取值如何變化，公式(5)總是加權(quán)歐氏距離。為避免人工設(shè)定參數(shù)，通過如下方式定義Ak：

為使dAk(yn,vk)能夠根據(jù)類別wk，k∈{u,c}的分布特征進(jìn)行自適應(yīng)調(diào)整，充分利用其形狀信息，本文使用類別wk的標(biāo)準(zhǔn)差來定義其距離誘導(dǎo)矩陣：假設(shè)σk為類別wk的標(biāo)準(zhǔn)差，則令A(yù)k=1/σk。通常情況下，類別wu與wc的標(biāo)準(zhǔn)差是未知的。本文通過如下方式進(jìn)行估算：首先利用FLICM對(duì)差分影像進(jìn)行聚類，而后根據(jù)聚類結(jié)果估算wu與wc的標(biāo)準(zhǔn)差σu和σc。

FATCD利用自適應(yīng)距離公式(5)增強(qiáng)FCM類算法，然后采用增強(qiáng)后算法對(duì)SAR差分影像進(jìn)行聚類，從而得到未變化類wu和變化類wc的隸屬度函數(shù)Mu和Mc。改進(jìn)后的FCM和FLICM的聚類過程與改進(jìn)前的聚類過程相似[10-11]。

1.2 生成變化檢測(cè)圖

得到未變化類和變化類的模糊隸屬度函數(shù)Mu和Mc后，本研究使用模糊拓?fù)鋬?yōu)化的去模糊化方案來生成變化檢測(cè)圖：首先利用模糊拓?fù)鋵⒉罘钟跋駝澐譃閮?nèi)部和模糊邊界區(qū)域；而后利用最大隸屬度原則和模糊拓?fù)涞闹芜B通性分別對(duì)內(nèi)部和模糊邊界像元分配類別。

1.2.1 分割差分影像

模糊拓?fù)涫菍?duì)傳統(tǒng)拓?fù)涞挠幸鏀U(kuò)展，是分析模糊集合的有力工具。為獲得未變化和變化類的不確定像元，本文將差分影像視作一個(gè)模糊拓?fù)淇臻g，該空間主要包含未變化集wu和變化集wc兩個(gè)模糊集合。模糊拓?fù)浼澳：負(fù)淇臻g的相關(guān)知識(shí)可參見文獻(xiàn)[20]，圖2(a)給出一個(gè)模糊拓?fù)淇臻g的簡(jiǎn)單示例。根據(jù)文獻(xiàn)[22]，可利用內(nèi)部算子來誘導(dǎo)模糊拓?fù)?，本文通過內(nèi)部算子Aα(·)

(6)

來誘導(dǎo)模糊拓?fù)?，并將Aα(·)的誘導(dǎo)拓?fù)溆涀鰽α-FT。接下來給出分割差分影像的方案：

圖2 模糊拓?fù)涫纠鼺ig.2 Fuzzy topology

(7)

(8)

分割差分影像的關(guān)鍵是確定閾值αk，k∈{u,c}。為此，提出一種自適應(yīng)的自動(dòng)搜尋αk最優(yōu)取值的算法。考慮到若搜尋步長(zhǎng)太大，則容易錯(cuò)過最優(yōu)解；若搜尋步長(zhǎng)太小，則會(huì)增加時(shí)間成本，同時(shí)考慮到αk∈(0.5,1)，本文根據(jù)經(jīng)驗(yàn)將搜尋步長(zhǎng)設(shè)定為0.05，那么可得到所提出算法搜尋最優(yōu)αk值的備選集合為C={ct|ct=0.5+0.05×t，t=1,2,…,9}。

模糊集合wk滿足：像元屬于類別wk的隸屬度越高，其屬于wk內(nèi)部的可能性越大，這要求αk的取值盡可能大；另一方面，αk取值越大，公式(7)得到的邊界像元數(shù)目越多，但當(dāng)?shù)玫降牟淮_定像元數(shù)目過多時(shí)，檢測(cè)精度會(huì)受到影響，這要求αk取值不能過大；根據(jù)經(jīng)驗(yàn)wk中邊界像元數(shù)目通常不超過其像元總數(shù)的10%。綜上給出如圖3(a)所示的搜尋算法:

(1) 計(jì)算比值集合R={Rt|Rt=Nt/N}，N和Nt分別表示集合{XD(i,j)|Mk(XD(i,j))>0.5}和{XD(i,j)|ct≥Mk(XD(i,j))>0.5}所包含的像元個(gè)數(shù)，其中ct為備選集合C的元素，t=1,2,…,9；并設(shè)定t=1；

(2) 判斷Rt是否滿足Rt>10%：若滿足，則令αk=ct-1，搜尋結(jié)束；否則轉(zhuǎn)入步驟③；

(3) 判斷ct是否滿足ct=0.95：若滿足，則令αk=0.95，搜尋結(jié)束；否則令t=t+1并轉(zhuǎn)入步驟②。

上述搜尋方案本質(zhì)上是搜尋滿足條件Rt≤10%的備選集合C的最大元素；若C的前8個(gè)元素均滿足Rt≤10%，則令αk=0.95。

圖3 流程及示例Fig.3 Flow chart and example

1.2.2 生成變化檢測(cè)圖

Mu(XD(i,j))>αu>0.5>1-Mu(XD(i,j))=Mc(XD(i,j))

(9)

如圖2(a)所示，模糊拓?fù)淇臻g包含兩種空間結(jié)構(gòu)：隸屬度空間(垂直方向)和像元的位置空間(x-y平面)。鑒于遙感影像的鄰域像元間存在很強(qiáng)的相關(guān)性，且真實(shí)世界中的變化通常是空間連通的區(qū)域，而不是孤立的點(diǎn)[12]，本研究借助已被標(biāo)記的內(nèi)部像元和模糊拓?fù)涞闹芜B通性(即x-y平面像元的位置連通)對(duì)邊界像元分配類別。具體步驟如下：

為敘述方便，用符號(hào)N8(i,j)表示像元XD(i,j)的八鄰域系統(tǒng)(參見圖3(b))。任給邊界像元XD(i,j)，假設(shè)Nu和Nc分別表示其八鄰域N8(i,j)中，已被標(biāo)記為未變化類和變化類的內(nèi)部像元個(gè)數(shù)，那么：①若Nu>Nc，則將邊界像元XD(i,j)標(biāo)記為未變化類；②若Nu

(10)

2 試驗(yàn)與分析

2.1 試驗(yàn)數(shù)據(jù)描述

第1組試驗(yàn)數(shù)據(jù)，即Bangladesh數(shù)據(jù)集，是由Envisat衛(wèi)星分別于2007年4月和2007年7月獲取的兩期Advanced SAR影像。該數(shù)據(jù)覆蓋了Bangladesh(孟加拉國)和印度部分地區(qū)。在影像獲取期間，由于雨季連續(xù)降雨造成了該地區(qū)的洪水泛濫。本文選取一個(gè)300×300像元的區(qū)域作為試驗(yàn)數(shù)據(jù)。前一期影像(圖4(a))反映了洪水發(fā)生前的情形，后一期影像(圖4(b))顯示了洪水淹沒的地區(qū)。圖4(c)為影像(a)和(b)的差分影像，圖4(d)為差分影像(c)的灰度直方圖。

第2組試驗(yàn)數(shù)據(jù)稱作Ottawa數(shù)據(jù)集，是分別于1997年5月和1997年8月由Radarsat衛(wèi)星在加拿大Ottawa地區(qū)獲取的兩期SAR影像，分辨率為12 m，大小為290×350像元。兩期影像的變化主要是受雨季影響引起的。圖4(e)為正處在雨季中(5月)的影像，由于大量降水導(dǎo)致河水水位上漲，圖4(f)為雨季過后(8月)的影像，河水退去露出大片的陸地。圖4(g)為影像(e)和(f)的差分影像，圖4(h)為差分影像(g)的灰度直方圖。

通過公式(1)生成差分影像后，為降低噪聲干擾，采用3×3窗口濾波器對(duì)Ottawa數(shù)據(jù)的差分影像進(jìn)行中值濾波；對(duì)于Bangladesh數(shù)據(jù)，因其包含線狀目標(biāo)，未對(duì)其差分影像實(shí)施濾波操作。

圖4 試驗(yàn)數(shù)據(jù)集Fig.4 Experimental datasets

2.2 試驗(yàn)設(shè)計(jì)

為驗(yàn)證所提出技術(shù)框架FATCD的適用性，本文基于FCM和FLICM給出兩種FATCD增強(qiáng)的模糊聚類算法，即FatFCM和FatFLICM。為證明本文方法的有效性，根據(jù)其特點(diǎn)組織了如下對(duì)比試驗(yàn)：①與標(biāo)準(zhǔn)FCM和FLICM比較；②與基于信息熵的無分布假設(shè)的Kapur閾值技術(shù)[23]比較；③與上下文敏感的EMMRF(expectation maximization+Markov random field)[8]算法比較。

本文通過定性和定量?jī)煞N方式來評(píng)價(jià)變化檢測(cè)技術(shù)的性能。前者通過比較各算法生成的變化檢測(cè)圖實(shí)現(xiàn)；后者通過計(jì)算每一幅變化檢測(cè)圖的4個(gè)精度指標(biāo)[24]實(shí)現(xiàn)：漏檢(missed detections，MD，即未被檢測(cè)出的變化像元個(gè)數(shù))、虛警(false alarms，F(xiàn)A，即被誤檢為變化的未變化像元個(gè)數(shù))、總體錯(cuò)誤(overall error，OE，即漏檢+虛警)、Kappa系數(shù)[25]。

2.3 增強(qiáng)效果分析

FATCD框架首先通過自適應(yīng)距離優(yōu)化模糊聚類算法，而后利用模糊拓?fù)淅碚摳倪M(jìn)去模糊化操作，為清楚兩個(gè)步驟的具體效果，本部分以FatFCM為例來分析FATCD對(duì)FCM類算法的增強(qiáng)過程。表1列出了FCM、AFCM、FTFCM和FatFCM對(duì)兩組試驗(yàn)數(shù)據(jù)變化檢測(cè)結(jié)果的精度指標(biāo)。其中AFCM為只用自適應(yīng)距離(AD)增強(qiáng)FCM得到的聚類算法，F(xiàn)TFCM表示只利用模糊拓?fù)?FT)增強(qiáng)FCM得到的聚類算法。

表1FCM,AFCM,FTFCM和FatFCM檢測(cè)結(jié)果的定量分析指標(biāo)

Tab.1QuantitativeanalysisindicesofdetectionresultsfromFCM,AFCM,FTFCMandFatFCM

檢測(cè)算法BangladeshMDFAOEKappaFCM4764647700.7543AFCM30515531060.8496FTFCM38571038670.8072FatFCM222914923780.8884檢測(cè)算法OttawaMDFAOEKappaFCM233740227390.8934AFCM150295824600.9077FTFCM182838922170.9149FatFCM998101720150.9255

從表1可知：①對(duì)兩組數(shù)據(jù)，AFCM和FTFCM的檢測(cè)效果均優(yōu)于FCM，表明自適應(yīng)距離和模糊拓?fù)渚茉鰪?qiáng)FCM的性能；②FatFCM的檢測(cè)效果優(yōu)于AFCM和FTFCM，表明兩種理論綜合的增強(qiáng)效果優(yōu)于其中任何一種理論；③對(duì)Bangladesh數(shù)據(jù)，AFCM的檢測(cè)性能優(yōu)于FTFCM，而對(duì)Ottawa數(shù)據(jù)，F(xiàn)TFCM的檢測(cè)性能優(yōu)于AFCM，表明自適應(yīng)距離和模糊拓?fù)鋵?duì)FCM的增強(qiáng)效果與數(shù)據(jù)本身性質(zhì)相關(guān)。

2.4 試驗(yàn)結(jié)果與精度評(píng)價(jià)

對(duì)Bangladesh數(shù)據(jù)，圖5給出了不同算法的變化檢測(cè)圖。為方便比較，每幅變化檢測(cè)圖被劃分為四部分：黑色和白色分別表示被正確檢測(cè)的未變化和變化像元，紅色表示漏檢錯(cuò)誤，黃色表示虛警錯(cuò)誤。圖5(a)—(f)分別為Kapur、FCM、FLICM、EMMRF、FatFCM以及FatFLICM的變化檢測(cè)圖，圖5(g)為Bangladesh數(shù)據(jù)的變化參考圖，通過詳細(xì)對(duì)比兩期原始影像并人工數(shù)字化生成。

圖5 各算法對(duì)Bangladesh數(shù)據(jù)的變化檢測(cè)結(jié)果Fig.5 Change -detection results on Bangladesh dataset from different algorithms

從圖5可知：①6種檢測(cè)算法都能較為有效地檢測(cè)出實(shí)地大部分的洪水區(qū)域。相比于其他算法，本文方法FatFCM和FatFLICM獲得了最接近真實(shí)地表變化的變化檢測(cè)圖，在增加少量虛警的情況下，很大程度上解決了其相應(yīng)標(biāo)準(zhǔn)算法FCM和FLICM的漏檢問題；②FCM和FLICM幾乎沒有虛警，但漏檢了大片洪水區(qū)域，整體檢測(cè)效果不佳；③Kapur算法一定程度上緩解了漏檢問題，但仍有較多漏檢像元。④EMMRF漏檢錯(cuò)誤最少，但較多虛警誤差影響了其整體檢測(cè)效果。

為更加客觀的評(píng)價(jià)各變化檢測(cè)技術(shù)的性能，表2給出了圖5中各變化檢測(cè)圖的定量分析指標(biāo)。從表2可看出：

(1) 本文方法FatFCM和FatFLICM的整體檢測(cè)效果優(yōu)于其他檢測(cè)技術(shù)，其總體錯(cuò)誤分別為2378和2485像元，比其他技術(shù)降低約500～3500像元，Kappa系數(shù)分別為0.888 4和0.883 8，比其他檢測(cè)技術(shù)提高約2.7%～13%。

(2) 相比于標(biāo)準(zhǔn)FCM和FLICM，本文方法檢測(cè)效果顯著提高：總體錯(cuò)誤降低2000像元以上，Kappa系數(shù)提高11%以上。這是因?yàn)楸疚姆椒ㄔ谟?jì)算像元隸屬度時(shí)采用更靈活的自適應(yīng)距離，得到了更準(zhǔn)確的隸屬度函數(shù)，并利用模糊拓?fù)淅碚撛鰪?qiáng)了模糊聚類的去模糊化過程。

表2圖5中各變化檢測(cè)結(jié)果的定量分析指標(biāo)

Tab.2Quantitativeanalysisindicesofchange-detectionresultsshowninFig.5

檢測(cè)算法MDFAOEKappaKapur29035929620.8573EMMRF450548159310.7745FCM4764647700.7543FLICM4581945900.7652FatFCM222914923780.8884FatFLICM221726824850.8838

(3) 與基于信息熵的Kapur閾值技術(shù)相比，本文方法漏檢錯(cuò)誤降低約700像元，檢測(cè)效果有了較大改善。這主要是因?yàn)殡m然兩者均無分布假設(shè)，但本文方法在變化檢測(cè)過程中不僅通過自適應(yīng)距離結(jié)合了未變化類與變化類的形狀信息，且通過支撐連通性利用了差分影像的空間鄰域信息。

(4) 與EMMRF技術(shù)相比，本文方法的虛警顯著減少，Kappa系數(shù)提高約11%，這是因?yàn)楸疚姆椒ǎ簾o分布假設(shè)，避免了假設(shè)分布模型與真實(shí)分布不太符合的影響(圖4(d))；只對(duì)差分影像的模糊邊界像元運(yùn)用空間信息，有效地抑制了對(duì)變化邊界的過度平滑。

對(duì)Ottawa數(shù)據(jù)，圖6(a)—(f)展示了不同算法的變化檢測(cè)圖；圖6(g)為其變化參考圖，通過人工對(duì)比兩期原始影像生成。表3給出了圖6中每幅變化檢測(cè)圖的4個(gè)精度指標(biāo)。

圖6 各算法對(duì)Ottawa數(shù)據(jù)的變化檢測(cè)結(jié)果Fig.6 Change -detection results on Ottawa dataset from different algorithms

Tab.3Quantitativeanalysisindicesofchange-detectionresultsshowninFig.6

檢測(cè)算法MDFAOEKappaKapur927208530120.8917EMMRF509218926980.9042FCM233740227390.8934FLICM237822426020.8982FatFCM998101720150.9255FatFLICM563167122340.9196

分析圖6和表3可得出與Bangladesh數(shù)據(jù)類似的結(jié)論，即對(duì)于Ottawa數(shù)據(jù)，本文方法的檢測(cè)效果仍是最好的：從定性的角度，圖6表明本文方法FatFCM和FatFLICM的變化檢測(cè)圖彩色(紅色和黃色)面積最小，與變化參考圖最接近。從定量的角度，表3也表明本文方法的檢測(cè)效果優(yōu)于其他檢測(cè)技術(shù)：其總體錯(cuò)誤最小，Kappa系數(shù)最高，能在總體錯(cuò)誤最小的情況下檢測(cè)出更多的變化區(qū)域。

兩組真實(shí)SAR影像數(shù)據(jù)的試驗(yàn)結(jié)果證明了所提出FATCD技術(shù)框架的有效性和對(duì)不同F(xiàn)CM類算法的適用性。另外，從圖4可知，兩組試驗(yàn)數(shù)據(jù)差分影像的灰度直方圖間存在著很大差異，這表明FATCD能夠適用于具有不同概率密度分布的差分影像。

3 結(jié) 論

本文提出一種改進(jìn)FCM及其變體的SAR影像變化檢測(cè)技術(shù)框架FATCD。FATCD通過自適應(yīng)距離優(yōu)化聚類過程中像元隸屬度的計(jì)算公式；通過模糊拓?fù)涓倪M(jìn)了去模糊化過程，為其取得較優(yōu)的檢測(cè)結(jié)果提供了理論保障。試驗(yàn)結(jié)果表明經(jīng)過FATCD技術(shù)增強(qiáng)后的FCM類型算法的變化檢測(cè)效果優(yōu)于Kapur、FCM、FLICM和EMMRF四種經(jīng)典變化檢測(cè)技術(shù)。

所提出的FATCD技術(shù)能夠較好地適用于中高分辨率的遙感影像，例如本文中使用的SAR影像。今后的工作重點(diǎn)是將FATCD應(yīng)用到更高分辨率的遙感影像。

參考文獻(xiàn)：

[1] 公茂果, 蘇臨之, 李豪, 等. 合成孔徑雷達(dá)影像變化檢測(cè)研究進(jìn)展[J]. 計(jì)算機(jī)研究與發(fā)展, 2016, 53(1): 123-137.

GONG Maoguo, SU Linzhi, LI Hao, et al. A Survey on Change Detection in Synthetic Aperture Radar Imagery[J]. Journal of Computer Research and Development, 2016, 53(1): 123-137.

[2] 胡召玲. 廣義高斯模型及KI雙閾值法的SAR圖像非監(jiān)督變化檢測(cè)[J]. 測(cè)繪學(xué)報(bào), 2013, 42(1) 116-122.

HU Zhaoling. An Unsupervised Change Detection Approach Based on KI Dual Thresholds under the Generalized Gauss Model Assumption in SAR Images[J]. Acta Geodaetica et Cartographica Sinica, 2013, 42(1): 116-122.

[3] 莊會(huì)富, 鄧喀中, 范洪冬. 紋理特征向量與最大化熵法相結(jié)合的SAR影像非監(jiān)督變化檢測(cè)[J]. 測(cè)繪學(xué)報(bào), 2016, 45(3): 339-346. DOI: 10.11947/j.AGCS.2016.20150022.

ZHUANG Huifu, DENG Kazhong, FAN Hongdong. SAR Images Unsupervised Change Detection Based on Combination of Texture Feature Vector with Maximum Entropy Principle[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(3): 339-346. DOI: 10.11947/j.AGCS.2016.20150022.

[4] 郝洪美, 張永紅, 石海燕, 等. 統(tǒng)計(jì)假設(shè)檢驗(yàn)方法在全極化SAR變化檢測(cè)中的應(yīng)用[J]. 遙感學(xué)報(bào), 2012, 16(3): 520-532.

HAO Hongmei, ZHANG Yonghong, SHI Haiyan, et al. Application of Test Statistic Method in Fully Polarimtric SAR Change Detection[J]. Journal of Remote Sensing, 2012, 16(3): 520-532.

[5] GONG Maoguo, ZHOU Zhiqiang, MA Jingjing. Change Detection in Synthetic Aperture Radar Images Based on Image Fusion and Fuzzy Clustering[J]. IEEE Transactions on Image Processing, 2012, 21(4): 2141-2151.

[6] 萬紅林, 焦李成, 王桂婷, 等. 在感興趣的區(qū)域?qū)用嫔线M(jìn)行SAR圖像變化檢測(cè)的方法研究[J]. 測(cè)繪學(xué)報(bào), 2012, 41(2): 239-245.

WAN Honglin, JIAO Licheng, WANG Guiting, et al. A Region-of-Interest Level Method for Change Detection in SAR Imagery[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(2): 239-245.

[7] 黃世奇, 劉代志, 胡明星, 等. 基于小波變換的多時(shí)相SAR圖像變化檢測(cè)技術(shù)[J]. 測(cè)繪學(xué)報(bào), 2010, 39(2): 180-186.

HUANG Shiqi, LIU Daizhi, HU Mingxing, et al. Multi-Temporal SAR Image Change Detection Technique Based on Wavelet Transform[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(2): 180-186.

[8] BRUZZONE L, PRIETO D F. Automatic Analysis of the Difference Image for Unsupervised Change Detection[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(3): 1171-1182.

[9] GONG Maoguo, SU Linzhi, JIA Meng, et al. Fuzzy Clustering with a Modified MRF Energy Function for Change Detection in Synthetic Aperture Radar Images[J]. IEEE Transactions on Fuzzy Systems, 2014, 22(1): 98-109.

[10] BEZDEK J C. Pattern Recognition with Fuzzy Objective Function Algorithms[M]. New York: Plenum Press, 1981: 65-80.

[11] KRINIDIS S, CHATZIS V. A Robust Fuzzy Local Information C-means Clustering Algorithm[J]. IEEE Transactions on Image Processing, 2010, 19(5): 1328-1337.

[12] SHAO Pan, SHI Wenzhong, HE Pengfei, et al. Novel Approach to Unsupervised Change Detection Based on a Robust Semi-Supervised FCM Clustering Algorithm[J]. Remote Sensing, 2016, 8(3): 264.

[13] GHOSH A, MISHRA N S, GHOSH S. Fuzzy Clustering Algorithms for Unsupervised Change Detection in Remote Sensing Images[J]. Information Sciences, 2011, 181(4): 699-715.

[14] MA Jingjing, GONG Maoguo, ZHOU Zhiqiang. Wavelet Fusion on Ratio Images for Change Detection in SAR Images[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(6): 1122-1126.

[15] 崔瑩, 熊博蒞, 蔣詠梅, 等. 結(jié)合結(jié)構(gòu)相似度的自適應(yīng)多尺度SAR圖像變化檢測(cè)[J]. 中國圖象圖形學(xué)報(bào), 2014, 19(10): 1507-1513.

CUI Ying, XIONG Boli, JIANG Yongmei, et al. Multi-Scale Approach Based on Structure Similarity for Change Detection in SAR Images[J]. Journal of Image and Graphics, 2014, 19(10): 1507-1513.

[16] MISHRA N S, GHOSH S, GHOSH A. Fuzzy Clustering Algorithms Incorporating Local Information for Change Detection in Remotely Sensed Images[J]. Applied Soft Computing, 2012, 12(8): 2683-2692.

[17] HAO Ming, ZHANG Hua, SHI Wenzhong, et al. Unsupervised Change Detection Using Fuzzy C-means and MRF from Remotely Sensed Images[J]. Remote Sensing Letters, 2013, 4(12): 1185-1194.

[18] CAI Weiling, CHEN Songcan, ZHANG Daoqiang. Fast and Robust Fuzzy C-means Clustering Algorithms Incorporating Local Information for Image Segmentation[J]. Pattern Recognition, 2007, 40(3): 825-838.

[19] BOUCHACHIA A, PEDRYCZ W. Enhancement of Fuzzy Clustering by Mechanisms of Partial Supervision[J]. Fuzzy Sets and Systems, 2006, 157(13): 1733-1759.

[20] LIU Yingming, LUO Maokang. Fuzzy Topology[M]. Singapore: World Scientific, 1997: 40-79.

[21] SHI Wenzhong, SHAO Pan, HAO Ming, et al. Fuzzy Topology-Based Method for Unsupervised Change Detection[J]. Remote Sensing Letters, 2016, 7(1): 81-90.

[22] LIU K, SHI Wenzhong. Computing the Fuzzy Topological Relations of Spatial Objects Based on Induced Fuzzy Topology[J]. International Journal of Geographical Information Science, 2006, 20(8): 857-883.

[23] PATRA S, GHOSH S, GHOSH A. Histogram Thresholding for Unsupervised Change Detection of Remote Sensing Images[J]. International Journal of Remote Sensing, 2011, 32(21): 6071-6089.

[24] BOVOLO F, BRUZZONE L, MARCONCINI M. A Novel Approach to Unsupervised Change Detection Based on a Semisupervised SVM and a Similarity Measure[J]. IEEE Transactions on Geoscience and Remote Sensing, 2008, 46(7): 2070-2082.

[25] CONGALTON R G. A Review of Assessing the Accuracy of Classifications of Remotely Sensed Data[J]. Remote Sensing of Environment, 1991, 37(1): 35-46.