楊青山， 張華

(1.武漢大學(xué)遙感信息工程學(xué)院,武漢 430079； 2.中國礦業(yè)大學(xué)(徐州)環(huán)境與測繪學(xué)院,徐州 221116)

?

融合像元形狀和光譜信息的高分遙感圖像分類新方法

楊青山1， 張華2

(1.武漢大學(xué)遙感信息工程學(xué)院,武漢 430079； 2.中國礦業(yè)大學(xué)(徐州)環(huán)境與測繪學(xué)院,徐州 221116)

在高空間分辨率(簡稱“高分”)遙感圖像分類中，由于存在“同譜異物”等現(xiàn)象，僅依靠光譜信息進(jìn)行分類的誤差較大。為提高圖像分類精度，提出一種融合像元形狀和光譜特征信息的高分多光譜遙感圖像分類新方法。首先利用像元及其鄰域的關(guān)系來描述其空間結(jié)構(gòu)，計(jì)算并提取像元同質(zhì)區(qū)域(pixel homogeneous regions,PHR)； 然后以所提取的同質(zhì)區(qū)域?yàn)榛A(chǔ)，分別計(jì)算中心像元的長/寬比(length-width ratio,LW)和面積/周長比(area-perimeter ratio,PAI)這2個(gè)像元形狀特征； 最后將歸一化后的像元形狀特征和光譜特征融合，并利用支持向量機(jī)分類方法進(jìn)行分類。以2個(gè)區(qū)域的QuickBird高分遙感圖像對(duì)該算法進(jìn)行驗(yàn)證，將實(shí)驗(yàn)結(jié)果與僅利用光譜信息分類和僅使用像元形狀指數(shù)(pixel shape index,PSI)分類的結(jié)果進(jìn)行比較。結(jié)果表明，所提出的方法得到的分類精度最高，該方法能有效地提高高分遙感圖像的分類精度。

像元同質(zhì)區(qū)域(PHR)； 像元形狀指數(shù)(PSI)； 閾值； 高空間分辨率遙感圖像

0 引言

傳統(tǒng)的基于像元的遙感數(shù)據(jù)監(jiān)督分類算法主要根據(jù)地物光譜數(shù)據(jù)，統(tǒng)計(jì)以像元為基本單元的圖像在不同光譜數(shù)據(jù)組合中的差別，對(duì)圖像進(jìn)行分類。但由于“同譜異物”現(xiàn)象和大量混合像元的存在，基于像元光譜特征的分類方法往往不能達(dá)到滿意的分類效果。近年來國內(nèi)外遙感研究人員提出了很多用于提高遙感圖像分類精度的方法來彌補(bǔ)僅依靠光譜特征進(jìn)行分類的不足。鄭淑丹等[1]基于分形理論和灰度共生矩陣(gray-level co-occurrence matrix,GLCM)2種方法提取不同的紋理特征，結(jié)合光譜信息對(duì)遙感圖像進(jìn)行最大似然監(jiān)督分類，得到較高的分類精度。李厚強(qiáng)等[2]將紋理特征和歸一化的色度和飽和度輸入反向傳播(back prepagation,BP)神經(jīng)網(wǎng)絡(luò)分類器進(jìn)行遙感圖像分類，較好地解決了光譜特征相似地物的區(qū)分問題。陳秋曉等[3]從圖像分割獲得的影像中提取光譜特征和拓?fù)涮卣鞯刃畔?，使影像易于識(shí)別和分類。張錦水等[4]基于支持向量機(jī)(support vector machine,SVM)的分類方法，利用復(fù)合光譜、紋理和結(jié)構(gòu)信息等多特征數(shù)據(jù)對(duì)IKONOS高空間分辨率(簡稱“高分”)圖像進(jìn)行分類，取得較高的分類精度。黃昕等[5]提出像元形狀指數(shù)(pixel shape index,PSI)，利用像元鄰域的光譜相似性探測其空間形狀特征，在高分遙感圖像分類中具有一定的優(yōu)勢。Shackelford等[6]提出一種長度-寬度提取算法，用于提取光譜特性相似的一組像元的長度與寬度信息，并得出該算法適合對(duì)不同長/寬比對(duì)象進(jìn)行分類的結(jié)論。Han等[7]提出的形狀尺寸指數(shù)(shape-size index, SSI)則是基于頻譜相似的同質(zhì)區(qū)域中的中心像元及其相鄰像元的特征信息提取的。國內(nèi)外學(xué)者已進(jìn)行的分類研究表明，從遙感圖像中可以提取大量的形狀、紋理和空間關(guān)系等特征信息，用于提高分類精度。

本文提出利用PSI，根據(jù)像元及其相鄰像元共同表征的形狀特征來提取影像的形狀特征信息； 再將歸一化后的形狀特征和光譜特征輸入SVM進(jìn)行高光譜遙感圖像分類； 詳細(xì)討論了像元形狀特征的提取方法以及長/寬比(length-width ratio,LW)和面積/周長比(area-perimeter ratio, PAI)以及融合二者的PL(PAI- LW)形狀特征對(duì)分類結(jié)果的影響； 并將實(shí)驗(yàn)結(jié)果與僅利用光譜信息分類和僅使用PSI分類的結(jié)果進(jìn)行比較。實(shí)驗(yàn)證明，在分類算法中融入對(duì)像元形狀特征的分析，能夠提高分類精度。

1 像元形狀特征

定義相鄰像元之間光譜相似的區(qū)域?yàn)橄裨|(zhì)區(qū)域(pixel homogeneous regions,PHR)，在PHR的基礎(chǔ)上提取2種像元形狀特征指標(biāo)，進(jìn)行圖像分類。

1.1 PHR確定

為了提高對(duì)空間特征的分類精度，利用中心像元與其8個(gè)相鄰像元之間的頻譜相似性計(jì)算PHR。

1.1.1 確定邊緣像元

邊緣像元的計(jì)算公式為

(1)

式中：edge(i)為估計(jì)概率值，表示第i個(gè)像元屬于邊緣像元的可能性；N為波段的個(gè)數(shù)；edged(i)為d波段的edge(i)值。edge(i)值用于判斷d波段中的第i個(gè)像元是否是邊緣像元： 當(dāng)edge(i)=1時(shí)，該像元是邊緣像元； 當(dāng)edge(i)=0時(shí)則不是。edge(i)的值用Canny濾波器獲得。

1.1.2 計(jì)算光譜相似性

光譜相似性的計(jì)算公式為

(2)

式中：PH(i)為第i個(gè)像元與其周圍像元的光譜相似性； arg min為使目標(biāo)函數(shù)取最小值時(shí)的變量值；Pd(i)和Pd(x)分別為d波段中的中心像元和與其相鄰像元的光譜值；x為與中心像元相鄰的8個(gè)像元。edge(i)值越大，該像元在此同質(zhì)區(qū)域中的可能性就越小。

1.1.3 形成中心像元的PHR

給定一個(gè)閾值T，若PH(i)

(3)

重復(fù)進(jìn)行光譜相似性計(jì)算和像元合并運(yùn)算，直到?jīng)]有新的像元需要合并，最終得到中心像元i的PHR(圖1)。

圖1 像元形狀特征計(jì)算示意圖

圖1中紅色像元是中心像元，灰色像元是PHR中除中心像元和邊緣像元之外的所有像元，綠色像元是PHR的邊緣像元。在所求出的PHR中進(jìn)行指標(biāo)提取的計(jì)算，遍歷圖像的所有像元，計(jì)算每1個(gè)像元的形狀特征，并進(jìn)行分類。

1.2 閾值設(shè)定

在構(gòu)造PHR的過程中，用來與PH(i)進(jìn)行比較的閾值通常都是人工設(shè)定的，可能會(huì)造成PHR提取結(jié)果有較大誤差。因此，本文采用自適應(yīng)閾值的方法設(shè)定閾值。在本質(zhì)上，圖像的邊緣像元可以表示整個(gè)圖像的局部不連續(xù)性，并可近似地表示整個(gè)圖像的特征，所以可利用PHR的邊緣像元獲得最佳閾值，即

(4)

式中：T(i)為第i個(gè)像元的PHR的閾值；Md為波段d中的像元數(shù)目；Xd(k)為邊緣像元在d波段中的光譜值。

1.3 像元形狀特征提取

本文在充分考慮相鄰像元的光譜相似性和相同形狀區(qū)域內(nèi)的像元具有相近特征值的情況下，提出了2個(gè)像元PHR的形狀特征值，即長/寬比(LW)和面積/周長比(PAI)。

1.3.1 LW

LW的計(jì)算公式為

LW=PLength/PWidth,

(5)

式中：PLength為與中心像元所在的PHR具有相同標(biāo)準(zhǔn)二階中心矩的橢圓的長軸長度；PWidth為與中心像元所在的PHR具有相同標(biāo)準(zhǔn)二階中心矩的橢圓的短軸長度(參見圖1)。

1.3.2 PAI

PAI的計(jì)算公式為

PAI=PPerimeter/PArea,

(6)

式中：PPerimeter為中心像元所在PHR的邊緣像元個(gè)數(shù)之和；PArea為中心像元所在PHR內(nèi)所有像元個(gè)數(shù)之和。當(dāng)PHR面積相同時(shí)，PAI越大，則區(qū)域越復(fù)雜。形狀特征提取的優(yōu)勢在于： ①考慮了邊緣像元； ②自適應(yīng)閾值用于獲得PHR； ③多邊形的特征在不確定的窗口內(nèi)。

2 實(shí)驗(yàn)方法與步驟

本文的實(shí)驗(yàn)方法是基于光譜和形狀特征，利用 SVM分類器進(jìn)行分類： ①獲得1景高分多光譜圖像，用Canny濾波器[8]提取邊緣信息，得到模糊邊緣圖； ②基于Canny濾波器和原始圖像，用所設(shè)定的閾值計(jì)算每1個(gè)像元的PHR； ③基于模糊邊緣圖、自適應(yīng)閾值、原始圖像以及基于PHR的像元形狀特征，獲得每1個(gè)像元的PHR； ④應(yīng)用SVM分類器[9]對(duì)圖像進(jìn)行基于光譜和形狀特征的分類。SVM分類器的核函數(shù)為徑向基函數(shù)(radial basis function, RBF)，采用交叉驗(yàn)證獲取最優(yōu)參數(shù)C和g的值[9]。PL為綜合了PAI和LW的分類結(jié)果。

實(shí)驗(yàn)步驟如下：

1)獲取模糊邊緣圖。為了在求取原始圖像的最小變化值和保存邊緣信息時(shí)避免人工的干擾，采用中值濾波器，首先用一個(gè)3×3鄰域窗口對(duì)圖像進(jìn)行運(yùn)算； 然后用 Canny算法獲得濾波后圖像中的邊緣信息，得到edge(i)值； 最后用式(1)計(jì)算得到模糊邊緣圖。

2)按照前文提到的方法設(shè)定一個(gè)合適的閾值。

3)獲取像元的PHR。利用獲取的模糊邊緣圖，根據(jù)設(shè)定的閾值選取像元的PHR。

4)獲得像元的形狀特征集。形狀特征包括長/寬比(LW)和面積/周長比(PAI)。

5)基于SVM分類。將提取的像元形狀特征和多光譜特征均標(biāo)準(zhǔn)化為[0,1]，使所提取的形狀特征和多光譜特征集成1個(gè)整體特征(PL)； 然后應(yīng)用SVM分類器對(duì)所集成的圖像進(jìn)行分類。

3 實(shí)驗(yàn)結(jié)果與分析

本文實(shí)驗(yàn)采用遙感圖像和地面參考數(shù)據(jù)作為分類數(shù)據(jù)。分類算法是廣泛使用的基于Matlab7.8、由LIBSVM提供的SVM庫。用2個(gè)實(shí)驗(yàn)論證本文提出的新方法。將像元形狀特征和光譜特征組合作為分類特征，輸入SVM分類器中進(jìn)行分類。

3.1 實(shí)驗(yàn)一

3.1.1 研究數(shù)據(jù)

實(shí)驗(yàn)采用分層隨機(jī)抽樣的方法，從地面參考數(shù)據(jù)中分別抽取各類別樣本數(shù)據(jù)總數(shù)的1%作為訓(xùn)練樣本數(shù)據(jù)集，其余數(shù)據(jù)作為檢驗(yàn)樣本數(shù)據(jù)集。原始遙感圖像、地面實(shí)況圖像、實(shí)驗(yàn)用樣本及提取的不同形狀特征圖像如圖2所示。

(a) 遙感圖像 (b) 地面實(shí)況(c) 訓(xùn)練樣本

(d) 檢驗(yàn)樣本(e) PAI (f) LW

圖2 實(shí)驗(yàn)一樣本及形狀特征影像

Fig.2 Samples and shape feature of image for experiment 1

3.1.2 分類結(jié)果分析

圖2(e)(f)是提取的2種像元形狀特征，將其與圖像的光譜信息作為分類的特征信息輸入SVM分類器，分別對(duì)不同的形狀特征數(shù)據(jù)進(jìn)行了分類，結(jié)果見圖3。

(a) 光譜分類 (b) PAI分類(c) LW分類

(d) PL分類(e) PSI分類

圖3 實(shí)驗(yàn)一SVM分類結(jié)果

Fig.3 Results of SVM classification for experiment 1

通過實(shí)驗(yàn)一可知，陰影、草地和樹木的光譜特性很相似，故僅使用光譜信息很難得到好的分類效果。要想提高分類精度，必須引入其他的分類特征信息。圖3(a)是僅使用光譜信息分類的分類結(jié)果； 圖3(b),(c),(d)和(e)分別是加入了PAI,LW,PL和PSI像元形狀特征后的分類結(jié)果。從表1可以看出，加入像元形狀特征后的分類精度得到了大幅度的提高，最高的(PL)分類總體精度比僅使用光譜信息分類提高了8.47%。不同的像元形狀特征對(duì)各類地物分類精度的提高程度也不同，從各個(gè)類的分類精度提高情況來看，雖然單一利用PAI和LW形狀特征對(duì)某些地物的分類精度比PSI的低，但是綜合了PAI和LW的PL形狀特征對(duì)各個(gè)類的分類精度的提高最為明顯，最終的總體分類精度達(dá)到了94.66%，Kappa系數(shù)達(dá)到了0.932 1，均比PSI的分類精度高。說明在分類圖像中存在各種復(fù)雜形狀的圖斑，很難用單一的形狀特征加以區(qū)分和識(shí)別，在實(shí)際分類中應(yīng)采用多種像元形狀特征組合來提高分類的精度。

表1 實(shí)驗(yàn)一的分類精度統(tǒng)計(jì)

3.2 實(shí)驗(yàn)二

3.2.1 研究數(shù)據(jù)

實(shí)驗(yàn)二同樣采用分層隨機(jī)抽樣的方法，從地面參考數(shù)據(jù)中分別抽取了各類別樣本數(shù)據(jù)總數(shù)的1%作為訓(xùn)練樣本數(shù)據(jù)集，其余數(shù)據(jù)作為驗(yàn)證樣本數(shù)據(jù)集。實(shí)驗(yàn)二樣本及不同形狀特征圖像如圖4所示。

(a) 遙感圖像 (b) 地面實(shí)況(c) 訓(xùn)練樣本

(d) 檢驗(yàn)樣本(e) PAI(f) LW

圖4 實(shí)驗(yàn)二樣本及形狀特征影像

Fig.4 Samples and shape feature images for experiment 2

3.2.2 分類結(jié)果分析

圖4(e)(f)是提取的2種像元形狀特征的圖像，將其和圖像的光譜特征信息作為分類的特征信息輸入SVM分類器，分別對(duì)不同的形狀特征數(shù)據(jù)進(jìn)行了分類，結(jié)果如圖5所示。

(a) 光譜分類(b) PAI分類(c) LW分類

(d) PL分類(e) PSI分類

圖5 實(shí)驗(yàn)二 SVM分類結(jié)果

Fig.5 Results of SVM classification for experiment 2

從實(shí)驗(yàn)二的分類效果來看，各種加入像元形狀特征后的分類效果要比僅使用光譜特征進(jìn)行分類的效果好。圖5(a)是僅使用光譜信息的分類結(jié)果，圖中水體與陰影、道路與房屋、裸地與房屋及草地與樹的錯(cuò)分現(xiàn)象比較嚴(yán)重，其原因是它們的光譜特性太相似，故僅利用光譜信息來區(qū)分是很難達(dá)得好的分類效果的，必須加入形狀特征等其他信息。圖5(b),(c),(d)和(e)分別是加入了PAI,LW,PL和PSI像元形狀特征后的分類結(jié)果。從表2可以看出，加入像元形狀特征后的分類精度得到了大幅度的提高，其中PL分類的總體分類精度達(dá)到了93.47%，在所有特征組合分類精度中最高，比僅利用光譜信息提高了18.61%，對(duì)樹木、裸地、道路、陰影、草地、房屋和水體分類精度分別提高了57.69%,50.70%,43.10%,39.65%,6.01%,5.33%和4.19%。PL對(duì)各個(gè)類的分類精度的提高是所有分類特征中最大的，再次說明了分類圖像中存在各種形狀復(fù)雜的圖斑，很難用單一的形狀特征來識(shí)別，在實(shí)際分類中應(yīng)該采用多種像元形狀特征組合來提高分類的精度。

表2 實(shí)驗(yàn)二的分類精度統(tǒng)計(jì)

4 結(jié)論

1)由于地物之間的光譜相似性，傳統(tǒng)的僅利用光譜信息的分類方法不能得到精度滿意的分類結(jié)果； 要想提高分類精度，需要更多考慮圖像中地物的結(jié)構(gòu)、形態(tài)、分布等空間特征信息的遙感圖像分類技術(shù)。

2)引入長/寬比(LW)和面積/周長比(PAI)這2個(gè)像元形狀特征，與影像的光譜信息組合作為SVM分類器的輸入數(shù)據(jù)，可提高遙感圖像的分類精度。2個(gè)實(shí)驗(yàn)的結(jié)果表明： 所引入的像元形狀特征充分考慮了相鄰像元間關(guān)系等形狀和空間位置特征，能大幅度提高遙感圖像的分類精度； 不同像元形狀特征對(duì)不同遙感圖像中不同地物的分類精度的提高程度也有所不同。

3)1種像元形狀特征難以適應(yīng)不同的地物類別，所以在實(shí)際分類過程中，應(yīng)將多種形狀特征組合用于遙感圖像分類，以得到相對(duì)滿意的分類結(jié)果； 對(duì)于多特征組合分類，一般需要進(jìn)行特征選擇和優(yōu)化。

4)本文中所涉及的形狀特征有限，還需研究更多的特征，以保證所得結(jié)論的普適性和通用性。在特征提取中，對(duì)尺度的影響及在像元形狀特征計(jì)算中的閾值選擇規(guī)則等問題尚須進(jìn)一步研究。

[1] 鄭淑丹,鄭江華,石明輝,等.基于分形和灰度共生矩陣紋理特征的種植型藥用植物遙感分類[J].遙感學(xué)報(bào),2014,18(4):868-886. Zheng S D,Zheng J H,Shi M H,et al.Classification of cultivated Chinese medicinal plants based on fractal theory and gray level co-occurrence matrix textures[J].Journal of Remote Sensing,2014,18(4):868-886.

[2] 李厚強(qiáng),劉政凱,林峰.基于分形理論的航空?qǐng)D像分類方法[J].遙感學(xué)報(bào),2001,5(5):353-357. Li H Q,Liu Z K,Lin F.Aerial image classification method based on fractal theory[J].Journal of Remote Sensing,2001,5(5):353-357.

[3] 陳秋曉,駱劍承,周成虎,等.基于多特征的遙感影像分類方法[J].遙感學(xué)報(bào),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.

[4] 張錦水,何春陽,潘耀忠,等.基于SVM的多源信息復(fù)合的高空間分辨率遙感數(shù)據(jù)分類研究[J].遙感學(xué)報(bào),2006,10(1):49-57. Zhang J S,He C Y,Pan Y Z,et al.The high spatial resolution RS image classification based on SVM method with the multi-source data[J].Journal of Remote Sensing,2006,10(1):49-57.

[5] 黃昕,張良培,李平湘.融合形狀和光譜的高空間分辨率遙感影像分類[J].遙感學(xué)報(bào),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.

[6] Shackelford A K,Davis C H.A hierarchical fuzzy classification approach for high-resolution multispectral data over urban areas[J].IEEE Transactions on Geoscience and Remote Sensing,2003,41(9):1920-1932.

[7] Han Y,Kim H,Choi J,et al.A shape-size index extraction for classification of high resolution multispectral satellite images[J].International Journal of Remote Sensing,2012,33(6):1682-1700.

[8] 薛麗霞,李濤,王佐成.一種自適應(yīng)的Canny邊緣檢測算法[J].計(jì)算機(jī)應(yīng)用研究,2010,27(9):3588-3590. Xue L X,Li T,Wang Z C.Adaptive Canny edge detection algorithm[J].Application Research of Computers,2010,27(9):3588-3590.

[9] 付紀(jì)剛.基于支持向量機(jī)的邊緣檢測算法研究[D].西安:西安科技大學(xué),2011. Fu J G.Research of Edge Detection Based on Support Vector Machine[D].Xi’an:Xi’an University of Science and Technology,2011.

[10]張華.遙感數(shù)據(jù)可靠性分類方法研究[D].徐州:中國礦業(yè)大學(xué),2012. Zhang H.Study on Reliable Classification Methods Based on Remotely Sensed Image[D].Xuzhou:China University of Mining and Technology,2012.

[11]Chang C C,Lin C J.LIBSVM:A library for support vector machines[EB/OL].2001.http://www.csie.ntu.edu.tw/～cjlin/libsvm/.

[12]Huang X,Zhang L P,Li P X.Classification and extraction of spatial features in urban areas using high-resolution multispectral imagery[J].IEEE Geoscience and Remote Sensing Letters,2007,4(2):260-264.

[13]Zhang L P,Huang X,Huang B,et al.A pixel shape index coupled with spectral information for classification of high spatial resolution remotely sensed imagery[J].IEEE Transactions on Geoscience and Remote Sensing,2006,44(10):2950-2961.

(責(zé)任編輯： 劉心季)

A new method for classification of high spatial resolution remotely sensed imagery based on fusion of shape and spectral information of pixels

YANG Qingshan1， ZHANG Hua2

(1.SchoolofRemoteSensingandInformationEngineering，WuhanUniversity，Wuhan430079，China； 2.SchoolofEnvironmentScienceandSpatialInformatics，ChinaUniversityofMiningandTechnology(Xuzhou)，Xuzhou221116，China)

In the classification of high spatial resolution remotely sensed imagery，due to the presence of the same object with different spectra, the dependence only on spectral information for classification is not enough. To improve the accuracy of classification, the authors proposed a novel spatial features extraction method for classification of the HSRMI. Firstly, neighborhood pixels’ spatial relationship was described and used to calculate and extract the pixel homogeneous regions (PHR). Then, based on the extracted PHR, the pixels’ shape index features, including length-width ratio(LW) and area-perimeter ratio(PAI), were extracted. Lastly, the pixel shape index features were normalized and combined with the spectral information to perform classification by using SVM classification method. Two different areas’ QuickBird images were used to test the performance of proposed method. The experimental results show that the proposed method has the highest performance compared with pixel shape index(PSI)and spectral information, and can improve the classification accuracy of high spatial resolution remotely sensed imagery.

pixel homogeneous regions(PHR)；pixel shape index(PSI)；threshold；high spatial resolution remotely sensed imagery

10.6046/gtzyyg.2016.04.10

楊青山,張華.融合像元形狀和光譜信息的高分遙感圖像分類新方法[J].國土資源遙感,2016,28(4)：64-70.(Yang Q S,Zhang H.A new method for classification of high spatial resolution remotely sensed imagery based on fusion of shape and spectral information of pixels[J].Remote Sensing for Land and Resources,2016,28(4)：64-70.)

2016-01-13；

2016-02-18

TP 751.1

A

1001-070X(2016)04-0064-07

楊青山(1993-),男,碩士研究生,主要研究方向?yàn)镚IS和遙感圖像處理與特征提取。Email： wilm_yang@foxmail.com。

張華(1979-)，男,博士,副教授,主要從事遙感數(shù)據(jù)不確定性、空間分析及GIS算法與應(yīng)用系統(tǒng)開發(fā)等方面的研究。Email: zhhua_79@163.com。