李亮， 張云， 李勝， 應(yīng)國偉

(四川省第三測繪工程院，成都 610500)

融合空間關(guān)系的遙感圖像分類

李亮， 張云， 李勝， 應(yīng)國偉

(四川省第三測繪工程院，成都 610500)

針對光譜紋理特征分類方法的不足，提出了一種融合空間關(guān)系的遙感圖像分類方法。利用直方圖提取像斑特征，采用G統(tǒng)計量構(gòu)建單像斑概率，通過迭代統(tǒng)計方法計算地物類別鄰接概率，利用地物類別鄰接概率表達像斑鄰域概率，加權(quán)組合單像斑概率與像斑鄰域概率構(gòu)建像斑聯(lián)合概率，依據(jù)最大后驗概率準則獲取圖像分類結(jié)果。在QucikBird圖像上的試驗結(jié)果表明： 與傳統(tǒng)的光譜紋理分類方法相比，該方法能夠提高圖像分類的精度； 整體分類精度與Kappa系數(shù)分別提高了1.5%和2.1%。

空間關(guān)系； 像斑直方圖； G統(tǒng)計量； 地物類別鄰接概率； 圖像分類

0 引言

遙感圖像分類是通過地物特征及判別法則確定圖像中地物類別屬性的過程[1]，廣泛應(yīng)用于土地覆蓋分類[2]、專題信息提取[3]、土地利用變化檢測[4]等領(lǐng)域。高分辨率遙感圖像中不僅包括直接的光譜特征，還包括豐富的紋理[5-7]、形狀[8-9]、空間[10-11]等特征。僅利用光譜特征進行圖像分類的精度有限，融合多特征的影像分類方法能有效提高影像分類的精度[12]?？臻g特征用來表達地物在空間上的相關(guān)性，在高分辨率遙感圖像分類中廣泛應(yīng)用。常用空間特征表達方法可劃分為4類： 簡單空間約束法、地物類別鄰接矩陣法、多尺度空間特征法和馬爾科夫隨機場法。簡單空間約束法以地物空間位置上越相近，則地物類別越相似為空間約束進行分類[13]； 地物類別鄰接矩陣法利用相鄰地物的鄰接概率來表達空間特征[14]； 多尺度空間特征法對圖像進行多尺度分割，提取不同尺度下像斑的特征向量表達空間特征[15]； 馬爾科夫隨機場法中利用吉布斯分布來表達地物的空間關(guān)聯(lián)，吉布斯參數(shù)用來反映地物的空間特征[16-18]。4種方法中，馬爾科夫隨機場法的理論最為嚴密，但吉布斯參數(shù)的求解較為復雜； 其他3種方法的理論基礎(chǔ)較弱，但計算簡單，實用性強。

本文提出的融合空間關(guān)系的圖像分類方法，是一種基于地物類別鄰接矩陣的空間特征表達方法，能夠充分挖掘地物間的空間關(guān)系，提高遙感圖像分類的精度。采用G統(tǒng)計量構(gòu)建單像斑概率，利用地物類別鄰接矩陣計算像斑鄰域概率，依據(jù)單像斑概率自適應(yīng)確定像斑鄰域概率權(quán)重，加權(quán)組合單像斑概率與像斑鄰域概率構(gòu)建像斑的聯(lián)合概率，在最大后驗概率基礎(chǔ)上獲取圖像分類結(jié)果。

1 研究方法

對象，在影像上又稱之為像斑，是一系列光譜相似、空間相鄰的像元集合。本文方法是一種面向?qū)ο蟮谋O(jiān)督分類方法，圖1為方法流程圖。

圖1 流程圖

本文通過多尺度影像分割獲取像斑，選取訓練樣本像斑，利用直方圖提取像斑特征，采用G統(tǒng)計量度量像斑直方圖間的距離，構(gòu)建單像斑概率，依據(jù)最大后驗概率準則獲取初始分類結(jié)果； 遍歷所有像斑計算地物類別鄰接概率，利用地物類別鄰接概率構(gòu)建像斑鄰域概率，依據(jù)單像斑概率確定像斑鄰域概率權(quán)重，加權(quán)組合單像斑概率與像斑鄰域概率計算像斑聯(lián)合概率，采用最大后驗概率準則更新分類結(jié)果，迭代此過程，直到滿足收斂條件。

1.1 像斑聯(lián)合概率

令地物類別數(shù)為m，對應(yīng)的類別集合為ω={ω1,ω2,…,ωm}，其中ωi(1≤i≤m)表示第i種地物類別。影像分割后被劃分為n個像斑的集合R={R1,R2,…,Rn}，其中Rk(1≤k≤n)表示第k個像斑； 像斑特征向量的集合為X={X1,X2,…,Xn}，其中Xk(1≤k≤n)表示像斑Rk的特征向量； 像斑地物類別的集合為C={C1,C2,…,Cn}，其中Ck(1≤k≤n)表示像斑Rk的地物類別，且Ck∈ω。令像斑鄰域系統(tǒng)集合為σ={σ1,σ2，…，σn}，其中σk(1≤k≤n)表示像斑Rk的鄰域像斑集合。鄰域系統(tǒng)具有如下特性： ①σk?R； ②Rkσk； ③若Rh∈σl，則Rl∈σh，1≤h，l≤n。

像斑的地物類別不僅依賴于像斑的特征，還受到鄰域像斑地物類別的影響。本文立足于像斑的特征向量與鄰域像斑的類別，確定像斑的地物類別。兼顧像斑自身特征及鄰域像斑類別，像斑Rk的聯(lián)合概率為

p(ωi|Xk,σk)=(1-λk)p(ωi|Xk)+λkp(ωi|σk),

(1)

式中：Xk和σk分別表示像斑Rk的特征向量與鄰域像斑；p(ωi|Xk,σk)為像斑聯(lián)合概率，表示像斑在特征向量為Xk、鄰域像斑為σk的條件下，地物類別為ωi的條件概率；p(ωi|Xk)為單像斑概率，表示像斑特征向量為Xk時，像斑地物類別為ωi的條件概率；p(ωi|σk)為表達空間特征的像斑鄰域概率，表示在像斑鄰域為σk時，當前像斑地物類別為ωi的概率；λk表示像斑鄰域概率的權(quán)重，0≤λk≤1，不同的像斑對應(yīng)不同的λk值。當λk=0時，像斑聯(lián)合概率即為單像斑概率，像斑的類別僅由像斑的特征決定，而與像斑的鄰域像斑無關(guān)； 當λk=1時，像斑聯(lián)合概率即為像斑鄰域概率，像斑的類別僅由鄰域像斑類別決定，而與像斑的特征無關(guān)。

1.2 單像斑概率

單像斑概率基于像斑自身的特征，衡量其歸屬某地物類別的概率。文中利用像斑直方圖提取像斑的特征。像斑直方圖是用來表達像斑內(nèi)灰度值及其出現(xiàn)頻數(shù)的二維圖，其橫坐標為灰度值g(0≤g≤L-1)，縱坐標為灰度值出現(xiàn)的頻數(shù)h(g)，其中L表示灰度級。像斑直方圖統(tǒng)計簡單，能夠較好地描述像斑內(nèi)部像元的灰度分布信息，在表達像斑光譜特征的同時，又能表達像斑紋理特征。為了簡便，通常將像斑直方圖歸一化，利用灰度值出現(xiàn)的頻率代替頻數(shù)。像斑直方圖同灰度級密切相關(guān)，灰度級過大，像斑直方圖能夠較好表達像斑的細節(jié)信息，但直方圖會過于稀疏，降低直方圖距離度量的準確性； 灰度級過小，像斑直方圖分布較為緊湊，但會丟失部分細節(jié)信息，降低像斑間的可區(qū)分性。

利用直方圖提取像斑特征后，像斑間的異質(zhì)性轉(zhuǎn)換為直方圖距離。度量直方圖距離的方法較多，如KL距離法、直方圖相交法及G統(tǒng)計量等。本文選擇G統(tǒng)計量度量直方圖的距離。G統(tǒng)計量計算簡單，且無需對直方圖的分布預先作任何假設(shè)[19-21]。令灰度級為L，2個像斑的歸一化直方圖分別為f和t，fj和tj分別表示灰度值j對應(yīng)的頻率，0≤j≤L-1。像斑直方圖間的G統(tǒng)計量為

(2)

像斑直方圖的累計頻率之和為1，則

(3)

(4)

(5)

將式(3)—(5)代入式(2)，即

(6)

式中，G統(tǒng)計量的取值范圍為[0,4 ln 2]。2個像斑的直方圖完全相同時，G統(tǒng)計量取最小值0； 2個像斑的直方圖不相關(guān)時，G統(tǒng)計量取最大值4 ln 2。G統(tǒng)計量的值越大，像斑直方圖的距離越大，像斑間異質(zhì)性越大； 反之，G統(tǒng)計量的值越小，像斑直方圖的距離越小，像斑間異質(zhì)性越小。

利用G統(tǒng)計量可衡量待分類像斑與各訓練樣本像斑間的距離。受圖像中同物異譜現(xiàn)象的影響，同一地物類別下會存在較多子類； 選擇訓練樣本像斑時，同一地物類別需選取多個像斑。將待分類像斑與ωi類訓練樣本像斑中的最短距離作為其到ωi類的距離(圖2)，即

(7)

圖2 單像斑概率

利用待分類像斑到各地物類別的距離，采用倒數(shù)加權(quán)法建立像斑與地物類別距離到像斑對地物類別隸屬度的映射，即

(8)

式中：p(ωi|Xk)為單像斑概率，用來衡量基于像斑自身特征，像斑Rk對ωi類的隸屬度；G-1(k，ωi)為像斑Rk到ωi類距離的倒數(shù)。

1.3 像斑鄰域概率

式(1)中，p(ωi|σk)用來表達像斑類別受其鄰域像斑類別的影響。像斑的一階鄰域是與像斑直接相鄰的像斑，像斑的二階鄰域是與像斑一階鄰域像斑相鄰的像斑，依此類推像斑的更高階鄰域。圖3為鄰域系統(tǒng)的示意圖。

圖3 鄰域系統(tǒng)

當前像斑為A，一階鄰域像斑為A1，二階鄰域像斑為A1∪A2。本文在一階鄰域基礎(chǔ)上對空間特征進行表達。地物類別鄰接概率可以定量表達地物類別的空間關(guān)聯(lián)，2種地物類別的鄰接概率越大，則相鄰的狀態(tài)越穩(wěn)定。文中利用鄰域中所有像斑與當前像斑類別鄰接概率的平均值來表達像斑鄰域概率，即

(9)

式中：σk為像斑Rk的一階鄰域；n(σk)為鄰域像斑數(shù)；Ru和Cu分別為鄰域像斑及其地物類別；p(ωi|Cu)為鄰接概率，表示鄰域像斑類別為Cu時，當前像斑類別為ωi的條件概率。p(ωj|ωi)表示鄰域像斑類別為ωi時，當前像斑類別為ωj的鄰接概率，即

p(ωj|ωi)=n(ωi,ωj)/n(ωi),

(10)

式中：n(ωi)為ωi類像斑的鄰接像斑總數(shù)；n(ωi,ωj)為鄰接像斑類別為ωi、當前像斑類別為ωj的像斑對個數(shù)。

1.4 權(quán)重確定

光譜、紋理等特征是像斑的主要特征，是劃分像斑地物類別的重要依據(jù)?？臻g關(guān)系雖然能夠在一定程度上輔助像斑類別的劃分，但不能作為劃分像斑地物類別的主要依據(jù)。融合空間關(guān)系進行影像分類時，空間關(guān)系權(quán)重的確定應(yīng)滿足以下兩個條件： 其一，空間關(guān)系的權(quán)重應(yīng)小于非空間關(guān)系的權(quán)重； 其二，當非空間特征無法有效區(qū)分像斑類別時，空間關(guān)系的權(quán)重應(yīng)較大。為了滿足條件一，將非空間特征的權(quán)重值設(shè)定在[0.5,1]，空間關(guān)系的權(quán)重設(shè)定在[0,0.5]； 為了滿足條件二，利用單像斑概率的最大值與次最大值的差值來衡量像斑非空間特征的可分性。差值較小時，像斑利用非空間特征的可分性較差，容易產(chǎn)生誤分，此時應(yīng)設(shè)置較大空間關(guān)系權(quán)重； 差值較大時，像斑利用非空間特征的可分性較強，此時應(yīng)設(shè)置較小空間關(guān)系權(quán)重。圖4是空間關(guān)系權(quán)重示意圖。

圖4 空間權(quán)重示意圖

圖4中有林地D和居民地C共2個像斑。D像斑非空間特征的可區(qū)分性較好，單像斑概率的最大值與次最大值相差較大，應(yīng)設(shè)定較小的空間關(guān)系權(quán)重； C像斑利用非空間特征的可區(qū)分性較差，與道路、裸地容易混淆，單像斑概率的最大值與次最大值相差較小，應(yīng)設(shè)定較大的空間關(guān)系權(quán)重?？臻g關(guān)系權(quán)重λk的計算公式為：

(11)

(12)

式中，max(·)和secmax(·)分別表示最大值與次最大值。

1.5 像斑類別劃分

從式(1)、(9)可以看出，要獲得像斑的分類結(jié)果，需要已知地物類別鄰接概率，而要獲得地物類別鄰接概率，需要已知像斑的分類結(jié)果。為此，采用迭代的方法來對像斑類別進行劃分： ①提取像斑直方圖特征，采用G統(tǒng)計量構(gòu)建像斑歸屬各地物類別的單像斑概率，依據(jù)單像斑概率最大原則獲取初始影像分類結(jié)果； ②遍歷像斑分類結(jié)果，采用式(10)統(tǒng)計地物類別鄰接矩陣，利用式(9)計算像斑鄰域概率； ③依據(jù)式(11)、(12)，利用單像斑概率的最大值與次最大值確定像斑鄰域概率的權(quán)重，加權(quán)組合單像斑概率與像斑鄰域概率構(gòu)建像斑的聯(lián)合概率，依據(jù)聯(lián)合概率最大原則更新影像分類結(jié)果； ④重復執(zhí)行②、③2步驟，直到連續(xù)2次類別發(fā)生改變的像斑小于一定比例。

2 實驗及分析

實驗數(shù)據(jù)為2002年武漢市武昌城區(qū)的QuickBird遙感圖像(如圖5(a)所示)，包含藍、綠、紅及近紅外4個波段，空間分辨率為2.4 m，影像大小為1 650×1 582像元。采用多尺度分割算法對遙感圖像進行分割，分割后獲得732個像斑(圖5(b))。選擇道路、耕地、居民地、林地、裸地和湖泊的訓練樣本像斑對圖像進行分類。利用直方圖提取像斑特征，直方圖灰度級設(shè)定為16。從圖5(a)中可以看出，道路、居民地、裸地3者的反射率較強，具有較高的亮度值，在圖像上的光譜紋理特征相似，因此這3種地物在圖像分類時容易互相混淆。

(a) QuickBird圖像 (b) 圖像分割圖

圖5實驗數(shù)據(jù)

Fig.5Experimentaldata

2.1 像斑空間特征分析

經(jīng)過3次迭代后，圖像分類結(jié)果不再發(fā)生改變，最終的地物類別鄰接矩陣見表1。從表中可以看出，對角線元素的值均較大，這表明地物類別具有較強的空間聚集性，同類地物相互鄰接的概率較大。所有地物與居民地的鄰接概率均較大，道路、耕地、林地、裸地、湖泊與居民地的鄰接概率分別為0.41，0.27，0.36，0.44和0.24，這是因為圖像對應(yīng)區(qū)域為城區(qū)，其中居民地的覆蓋面積較大，其他地物均分布在居民地的內(nèi)部及四周。城區(qū)中裸地主要來自居民地的拆遷，裸地較多被居民地圍繞，所以居民地與裸地的鄰接概率較大(為0.17)。耕地與湖泊的鄰接概率為0.15，這是因為耕地普遍分布在湖泊的周圍。地物類別鄰接矩陣所表達的地物空間關(guān)聯(lián)關(guān)系與實際情況基本一致，這表明了其表達地物空間關(guān)系的有效性。

2.2 像斑分類結(jié)果分析

為了衡量像斑分類結(jié)果的精度，利用目視解譯的方法繪制了一幅標準分類圖(圖6(a))； 同時為了驗證引入空間關(guān)系的有效性，將單像斑概率分類結(jié)果(圖6(b))與本文分類結(jié)果(圖6(c))進行了比較。

(a) 標準結(jié)果 (b) 單像斑概率分類 (c) 本文方法分類

圖6圖像分類結(jié)果

Fig.6Imageclassificationmap

由圖6可知，圖6(b)中大量居民地被錯誤劃分為道路(黑色矩形區(qū)域內(nèi))與裸地(黑色橢圓形區(qū)域內(nèi))，部分湖泊被錯誤劃分為道路(黑色圓形區(qū)域內(nèi))。利用本文方法進行分類，圖中大量被誤分為道路和裸地的居民地均能得到正確劃分，同時誤分為道路的部分湖泊也得到正確劃分。這是因為道路、居民地和裸地3者的光譜紋理特征相似，對應(yīng)的空間關(guān)系權(quán)重較大，而在空間上，居民地與居民地相鄰的概率大于居民地與道路、裸地相鄰的概率。鄰域像斑類別為居民地時，當前像斑類別為居民地的概率較大，因此被誤分為裸地和道路的像斑能夠得到正確劃分。部分湖泊由于反射率較強，導致誤分為道路。在空間上，湖泊與湖泊鄰接的概率大于湖泊與道路鄰接的概率； 引入空間關(guān)系后，在鄰接像斑為湖泊的條件下，當前像斑能被正確劃分為湖泊。對比圖6(b)和6(c)可以看出，除道路、居民地和裸地3類地物外，其他3類地物的分類結(jié)果變化不大。這是因為其他3類地物間的光譜紋理可區(qū)分性較好，對應(yīng)像斑的空間關(guān)系權(quán)重較小，對像斑分類結(jié)果的影響較小。為了定量對比2種方法的分類精度，繪制混淆矩陣表2和表3，其中列代表真實類別，行代表檢測類別，單元格中的數(shù)字表示像元個數(shù)。

表2 單像斑概率法混淆矩陣

表3 本文方法混淆矩陣

對比表2和表3可以看出，本文方法分類精度較優(yōu)。用戶總體精度從90.3%提高至91.8%，Kappa系數(shù)從87.0%提高至89.1%。各地物類別中居民地分類精度明顯改善，制圖精度從89.6%提高至93.4%。道路與裸地的精度也有提高，道路的用戶精度從73.1%提高至83.4%，裸地的用戶精度從74.9%提高至79.1%。耕地、林地和湖泊3類地物的分類精度改善程度較小，這是因為3類地物的光譜紋理可區(qū)分性較好，對應(yīng)的空間關(guān)系權(quán)重較小，空間關(guān)系對圖像分類的貢獻較小，分類結(jié)果變化不大。

3 結(jié)論

1) 地物類別鄰接概率可以較好地表達地物類別的空間關(guān)聯(lián)性。地物類別的鄰接概率越大，則地物在空間上相鄰的狀態(tài)越穩(wěn)定。

2) 依據(jù)單像斑概率自適應(yīng)確定像斑鄰域概率權(quán)重的方法具有可行性。像斑非空間特征的可分性較差時，像斑鄰域概率設(shè)定較大值，否則設(shè)定較小值。

3) 引入空間關(guān)系后，能夠提高圖像分類精度?？臻g關(guān)系雖然是圖像分類時的輔助特征，但在非空間特征的可分性較差時，引入空間關(guān)系能夠改善圖像分類精度。

4) 本文在定量表達像斑鄰域概率時，僅以一階鄰域為基礎(chǔ)，在精確表達像斑空間關(guān)系方面不足，未來應(yīng)進一步研究在高階鄰域中的空間關(guān)系表達。

[1] 白穆,劉慧平,喬瑜,等.高分辨率遙感圖像分類方法在LUCC中的研究進展[J].國土資源遙感,2010(1):19-23.doi:10.6046/gtzyyg.2010.01.03. Bai M,Liu H P,Qiao Y,et al.New progress in the classification of high spatial resolution satellite images for LUCC[J].Remote Sensing for Land and Resources,2010(1):19-23.doi:10.6046/gtzyyg.2010.01.03.

[2] 楊耘,徐麗,顏佩麗.條件隨機場框架下基于隨機森林的城市土地利用/覆蓋遙感分類[J].國土資源遙感,2014,26(4):51-55.doi:10.6046/gtzyyg.2014.04.09. Yang Y,Xu L,Yan P L.Urban land use /cover classification of remote sensing using random forests under the framework of conditional random fields[J].Remote Sensing for Land and Resources,2014,26(4):51-55.doi:10.6046/gtzyyg.2014.04.09.

[3] 張樓香,阮仁宗,夏雙.洪澤湖濕地紋理特征參數(shù)分析[J].國土資源遙感,2015,27(1):75-80.doi:10.6046/gtzyyg.2015.01.12. Zhang L X,Ruan R Z,Xia S.Parameter analysis of image texture of wetland in the Hongze Lake[J].Remote Sensing for Land and Resources,2015,27(1):75-80.doi:10.6046/gtzyyg.2015.01.12.

[4] 楊勝,李敏,彭振國,等.一種新的多波段遙感影像變化檢測方法[J].中國圖象圖形學報,2009,14(4):572-578. Yang S,Li M,Peng Z G,et al.A novel multi-band remote sensing image change detection algorithm[J].Journal of Image and Graphics,2009,14(4):572-578.

[5] Ferro C J S,Warner T A.Scale and texture in digital image classification[J].Photogrammetric Engineering and Remote Sensing,2002,68(1):51-63.

[6] 王增茂,杜博,張良培,等.基于紋理特征和形態(tài)學特征融合的高光譜影像分類法[J].光子學報,2014,43(8):0810002-1-0810002-8. Wang Z M,Du B,Zhang L P,et al.Based on texture feature and extend morphological profile fusion for hyperspectral image classification[J].Acta Photonica Sinica,2014,43(8):0810002-1-0810002-8.

[7] 劉昌振,舒紅,張志,等.基于多尺度分割的高分遙感圖像變異函數(shù)紋理提取和分類[J].國土資源遙感,2015,27(4):47-53.doi:10.6046/gtzyyg.2015.04.08. Liu C Z,Shu H,Zhang Z,et al.Variogram texture extraction and classification of high resolution remote sensing images based on multi-resolution segmentation[J].Remote Sensing for Land and Resources,2015,27(4):47-53.doi:10.6046/gtzyyg.2015.04.08.

[8] 黃昕,張良培,李平湘.融合形狀和光譜的高空間分辨率遙感影像分類[J].遙感學報,2007,11(2):193-200. Huang X,Zhang L P,Li P X.Classification of high spatial resolution remotely sensed imagery based on the fusion of spectral and shape features[J].Journal of Remote Sensing,2007,11(2):193-200.

[9] Segl K,Roessner S,Heiden U,et al.Fusion of spectral and shape features for identification of urban surface cover types using reflective and thermal hyperspectral data[J].ISPRS Journal of Photogrammetry and Remote Sensing,2003,58(1/2):99-112.

[10]蔡曉斌,陳曉玲,王濤,等.基于圖斑空間關(guān)系的遙感專家分類方法研究[J].武漢大學學報(信息科學版),2006,31(4):321-324. Cai X B,Chen X L,Wang T,et al.Remote sensing expert classification method based on patch spatial relationship[J].Geomatics and Information Science of Wuhan University,2006,31(4):321-324.

[11]喬程,沈占鋒,吳寧,等.空間鄰接支持下的遙感影像分類[J].遙感學報,2011,15(1):88-99. Qiao C,Shen Z F,Wu N,et al.Remote sensing image classification method supported by spatial adjacency[J].Journal of Remote Sensing,2011,15(1):88-99.

[12]陳秋曉,駱劍承,周成虎,等.基于多特征的遙感影像分類方法[J].遙感學報,2004,8(3):239-245. Chen Q X,Luo J C,Zhou C H,et al.Classification of remotely sensed imagery using multi-features based approach[J].Journal of Remote Sensing,2004,8(3):239-245.

[13]趙紅蕊,閻廣建,鄧小煉,等.一種簡單加入空間關(guān)系的實用圖像分類方法[J].遙感學報,2003,7(5):358-363. Zhao H R,Yan G J,Deng X L,et al.A classification method based on spatial information[J].Journal of Remote Sensing,2003,7(5):358-363.

[14]王凱,舒寧,孔祥兵,等.一種多特征轉(zhuǎn)換的高光譜影像自適應(yīng)分類方法[J].武漢大學學報(信息科學版),2015,40(5):612-616. Wang K,Shu N,Kong X B,et al.A multi-feature conversion adaptive classification of hyperspectral image[J].Geomatics and Information Science of Wuhan University,2015,40(5):612-616.

[15]張倩,黃昕,張良培.多尺度同質(zhì)區(qū)域提取的高分辨率遙感影像分類研究[J].武漢大學學報(信息科學版),2011,36(1):117-121. Zhang Q,Huang X,Zhang L P.Multiscale image segmentation and classification with supervised ECHO of high spatial resolution remotely sensed imagery[J].Geomatics and Information Science of Wuhan University,2011,36(1):117-121.

[16]Zhao Y D,Zhang L P,Li P X,et al.Classification of high spatial resolution imagery using improved Gaussian Markov random-field-based texture features[J].IEEE Transactions on Geoscience and Remote Sensing,2007,45(5):1458-1468.

[17]龔龑,舒寧,王琰,等.遙感影像像斑空間關(guān)系分析的非規(guī)則無參數(shù)馬爾可夫隨機場模型[J].測繪學報,2013,42(1):101-107. Gong Y,Shu N,Wang Y,et al.Segments spatial relationship feature analysis of remote sensing imagery based on non-parameter irregular-MRF model[J].Acta Geodaetica et Cartographica Sinica,2013,42(1):101-107.

[18]譚熊,余旭初,張鵬強,等.基于MKSVM和MRF的高光譜影像分類方法[J].國土資源遙感,2015,27(3):42-46.doi:10.6046/gtzyyg.2015.03.08. Tan X,Yu X C,Zhang P Q,et al.Hyperspectral images classification based on MKSVM and MRF[J].Remote Sensing for Land and Resources,2015,27(3):42-46.doi:10.6046/gtzyyg.2015.03.08.

[19]Ojala T,Pietik?inen M.Unsupervised texture segmentation using feature distributions[J].Pattern Recognition,1999,32(3):477-486.

[20]Ojala T,Pietikainen M,Maenpaa T.Multiresolution gray-scale and rotation invariant texture classification with local binary patterns[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,2002,24(7):971-987.

[21]巫兆聰,胡忠文,張謙,等.結(jié)合光譜、紋理與形狀結(jié)構(gòu)信息的遙感影像分割方法[J].測繪學報,2013,42(1):44-50. Wu Z C,Hu Z W,Zhang Q,et al.On combining spectral, textural and shape features for remote sensing image segmentation[J].Acta Geodaetica et Cartographica Sinica,2013,42(1):44-50.

(責任編輯：李瑜)

Classificationofremotesensingimagesbasedonthefusionofspatialrelationship

LI Liang, ZHANG Yun, LI Sheng, YING Guowei

(TheThirdAcademyofEngineeringofSurveyingandMapping,Chengdu610500,China)

In order to overcome the disadvantages of the classification method based on spectral and texture features, the authors put forward a classification method based on the fusion of spatial relationship in this paper. Single object probability was built by G statistic after image object feature was extracted by histogram. The neighborhood object probability was described by land cover adjacency probability which was calculated by iterative statistics method. The joint probability of the object was built by the weighted combination of single object probability and neighborhood object probability. The classification result of the image was obtained according to the maximum a posteriori. The experimental results based on QuickBird image show that the proposed method can improve the classification accuracy compared with the traditional classifier using spectral and texture features. The overall classification accuracy and kappa coefficient are increased by 1.5% and 2.1%, respectively.

spatial relationship; object histogram; G statistics; land cover adjacency probability; image classification

10.6046/gtzyyg.2017.03.02

李亮，張云，李勝,等.融合空間關(guān)系的遙感圖像分類[J].國土資源遙感,2017,29(3):10-16.(Li L,Zhang Y,Li S,et al.Classification of remote sensing images based on the fusion of spatial relationship[J].Remote Sensing for Land and Resources,2017,29(3):10-16.)

2016-02-29；

2016-05-23

測繪地理信息公益性行業(yè)科研專項“衛(wèi)星遙感與地面?zhèn)鞲芯W(wǎng)一體化的湖泊流域地理國情監(jiān)測關(guān)鍵技術(shù)研究”(編號： 201512026)和數(shù)字制圖與國土信息應(yīng)用工程國家測繪地理信息局重點實驗室開放基金資助項目“基于遙感影像的矢量圖更新關(guān)鍵技術(shù)研究”(編號： DM2016SC04)共同資助。

李亮(1987-)，男，博士，工程師，主要從事遙感影像的智能化解譯研究。Email： liliang1987wuda@163.com。

TP 79

： A

： 1001-070X(2017)03-0010-07