陸翔 章皖秋 鄭雅蘭 岳彩榮

?

特征選擇的全極化SAR影像面向?qū)ο笸恋馗采w分類

陸翔 章皖秋 鄭雅蘭 岳彩榮

（西南林業(yè)大學(xué)，昆明 650224）

全極化SAR數(shù)據(jù)信息豐富，僅利用單一的極化特征和基于像元的分類很難得到較好的分類效果。因此，提出了全極化數(shù)據(jù)特征優(yōu)選結(jié)合面向?qū)ο蠓椒ㄟM(jìn)行土地覆蓋分類。以云南西雙版納州勐臘縣和普洱市思茅區(qū)的Terra SAR-X的X波段全極化雷達(dá)數(shù)據(jù)為信息源，首先對(duì)全極化SAR數(shù)據(jù)進(jìn)行預(yù)處理，提取研究區(qū)Pauli RGB圖像后，利用影像分割技術(shù)對(duì)Pauli RGB圖像進(jìn)行分割，作為分類的基本單元；然后對(duì)SAR影像提取極化分解特征和Span影像的紋理特征，選取最優(yōu)特征集合；最后利用面向?qū)ο竽：诸惙椒ㄟM(jìn)行土地覆蓋分類，并采用實(shí)地調(diào)查數(shù)據(jù)對(duì)分類結(jié)果進(jìn)行了精度評(píng)價(jià)。試驗(yàn)結(jié)果表明，面向?qū)ο蠓椒梢院芎玫厝コ肼暤挠绊?，最?yōu)組合的特征波段使得分類結(jié)果更加精確。西雙版納州勐臘縣總體分類精度達(dá)到88.5%，普洱市思茅區(qū)總體分類精度達(dá)到86.8%，較之H/A/α-Wishart分類方法精度提高了40%以上。

全極化 合成孔徑雷達(dá) 影像分割 特征選擇 模糊分類 遙感應(yīng)用

0 引言

近年來(lái)，極化合成孔徑雷達(dá)（Polarimetric Synthetic Aperture Radar，PolSAR）技術(shù)在土地覆蓋分類、城市規(guī)劃、植被高度估測(cè)、環(huán)境監(jiān)測(cè)和探測(cè)軍事目標(biāo)等領(lǐng)域有著廣泛的應(yīng)用[1]。其中，遙感影像分類是Pol SAR的重要研究?jī)?nèi)容之一，尤其是SAR數(shù)據(jù)的極化特征，進(jìn)一步提升了遙感技術(shù)對(duì)土地覆蓋進(jìn)行分類的能力。

基于像元和對(duì)象分類是當(dāng)前進(jìn)行SAR影像土地覆蓋分類的兩種主要形式，且分類方法主要采用支持向量機(jī)分類方法、距離分類方法、模糊分類法、神經(jīng)網(wǎng)絡(luò)分類方法，SAR數(shù)據(jù)特征選取的不同，這些方法的分類效果也各不相同?；谙裨腟AR影像土地覆蓋研究是利用SAR影像的極化信息，從極化散射矩陣中提取與散射機(jī)理和特性密切相關(guān)的參數(shù)，并結(jié)合其他分類方法直接進(jìn)行分類。文獻(xiàn)[2]利用散射熵和平均散射角組成的/平面劃分8個(gè)區(qū)域進(jìn)行非監(jiān)督分類；文獻(xiàn)[3]利用非監(jiān)督分類結(jié)合基于極化目標(biāo)分解分類；文獻(xiàn)[4]利用Freeman分解結(jié)合復(fù)合Wishart分布的最大似然分類器，提出了可以保持極化散射特性的分類方法。近年來(lái)，基于像元SAR圖像分類方法大多是在上述方法的基礎(chǔ)上進(jìn)行優(yōu)化，在SAR圖像分類上取得了重大突破，但是也存在不足[5-6]。文獻(xiàn)[2]/的平面采用人工線性劃分的方式，導(dǎo)致地物細(xì)節(jié)信息丟失。文獻(xiàn)[3]的自動(dòng)迭代分類方法無(wú)法保持地物的散射特性，導(dǎo)致誤分錯(cuò)分現(xiàn)象嚴(yán)重。這些基于像元的SAR圖像分類方法在類別劃分上依靠經(jīng)驗(yàn)，沒(méi)有有力的理論支撐，且存在較多斑點(diǎn)噪聲，嚴(yán)重影響分類精度[7]。基于對(duì)象的SAR圖像分類研究依據(jù)SAR圖像本身提取出有利于分類的各類特征參數(shù)，再結(jié)合圖像分割技術(shù)和適當(dāng)?shù)姆诸惙椒ㄟM(jìn)行圖像分類。例如，文獻(xiàn)[1]利用圖譜分割技術(shù)獲得分類對(duì)象，然后提取基于極化信息和紋理信息的44種分類特征并用于分類；文獻(xiàn)[8]利用分割和特征提取結(jié)合改進(jìn)的分類器算法實(shí)現(xiàn)SAR圖像的分類；文獻(xiàn)[9]提出了一種基于統(tǒng)計(jì)區(qū)域合并分割和支持向量機(jī)的新型面向?qū)ο蠓诸惙椒?。以上基于面向?qū)ο蟮腟AR圖像分類中，分類特征的選擇是面向?qū)ο蠓诸惖闹匾獥l件之一?；谙裨蛯?duì)象的兩類方法中均把提取的所有特征未經(jīng)過(guò)甄選全部用于分類，這樣會(huì)造成數(shù)據(jù)冗余，甚至分類精度的降低。

由于數(shù)據(jù)類型的限制，以上研究多集中于單極化、雙極化信息的分類，而對(duì)全極化分類的研究相對(duì)不多。在高分辨率的支持下，SAR圖像的紋理特征又可以作為分類的重要條件。因此，本文針對(duì)基于像元和基于對(duì)象分類的缺點(diǎn)，結(jié)合本文的全極化和高分辨率的Terra SAR數(shù)據(jù)，提出了特征波段的選擇與面向?qū)ο蠼Y(jié)合用于SAR影像地物分類的方法。首先對(duì)提取的Pauli RGB圖像進(jìn)行分割得到分類對(duì)象，然后提取出原始影像的極化分解特征和Span影像的紋理特征，并采用特征空間優(yōu)化得到最優(yōu)分類特征組合，利用模糊分類器得到分類結(jié)果。結(jié)果表明本文較高的分類精度可以滿足勐臘研究區(qū)和普洱研究區(qū)的土地變化監(jiān)測(cè)、森林動(dòng)態(tài)監(jiān)測(cè)、城鎮(zhèn)變化監(jiān)測(cè)等的需求。

1 研究方法

1.1 研究區(qū)和數(shù)據(jù)

選取了面積較小但是地類豐富、具有典型代表的兩塊試驗(yàn)區(qū)，見(jiàn)圖1。勐臘縣位于北緯21°09′~22°23′，東經(jīng)101°05′~101°50′之間。地處北熱帶濕潤(rùn)季風(fēng)氣候區(qū)，植被以常綠闊葉林和熱帶雨林為主，研究區(qū)內(nèi)地形以山地為主。普洱市位于北緯22°02′~24°50′，東經(jīng)99°09′~102°19′，屬亞熱帶季風(fēng)氣候，全年氣候溫和，垂直氣候明顯，全市以山地為主。天然植被以亞熱帶常綠闊葉林、思茅松林為主，也是普洱茶的主要種植區(qū)。

圖1 研究區(qū)示例

Terra SAR-X空間分辨率達(dá)到2.1m，能夠獲取高品質(zhì)的X波段數(shù)據(jù)。目前，Terra SAR-X已被廣泛應(yīng)用于變化檢測(cè)、林業(yè)、農(nóng)業(yè)、國(guó)土、水文以及城市規(guī)劃等領(lǐng)域。本文研究區(qū)的數(shù)據(jù)是全極化數(shù)據(jù)，（見(jiàn)表1），在西雙版納及普洱這種多云雨天氣的地區(qū)具有無(wú)可代替的優(yōu)勢(shì)。

表1 研究區(qū)Terra SAR數(shù)據(jù)參數(shù)

Tab.1 Terra SAR data in the study area

1.2 影像分割

以云南西雙版納州勐臘縣和普洱市思茅區(qū)的Terra SAR-X的X波段全極化雷達(dá)數(shù)據(jù)為信息源，進(jìn)行了基于特征選擇的極化SAR影像面向?qū)ο笸恋馗采w分類試驗(yàn)。

面向?qū)ο蟮腟AR影像土地覆蓋分類中，分割是決定分類精度的重要像因子，準(zhǔn)確的地物邊界分割是正確分類的必要條件。本文采用多尺度分割方法對(duì)研究區(qū)的Pauli RGB圖像進(jìn)行分割，得到分割對(duì)象，即具有相同統(tǒng)計(jì)特性像元組成的區(qū)域[10]，利用分割后的對(duì)象進(jìn)行分類可以有效避免噪聲對(duì)分類的影響。

分割時(shí)需要考慮兩個(gè)參數(shù)，一個(gè)是異質(zhì)度，一個(gè)是分割尺度。異質(zhì)度由顏色和形狀異質(zhì)度組成，二者權(quán)重之和為1。

式中color為顏色異質(zhì)度；為顏色異質(zhì)度的權(quán)重；shape為形狀異質(zhì)度；1–為形狀異質(zhì)度的權(quán)重，0≤≤1。compact為緊致度；s為緊致度的權(quán)重；smooth為光滑度；1–s為光滑度的權(quán)重，0≤s≤1。

根據(jù)影像的顏色和形狀特征設(shè)置分割參數(shù)，如果研究區(qū)內(nèi)地物形狀明顯則形狀的權(quán)重相應(yīng)的增大[11]。

式（1）中的color是兩個(gè)對(duì)象合并后得到的光譜異質(zhì)性值和合并前對(duì)象（obj1和obj2）的各自光譜異質(zhì)性值之和的差值

式中w為參與分割合并的波段權(quán)重，為參與分割的波段；merge，δ,merge分別為合并后的區(qū)域面積和光譜方差；δ,obj1，δ,obj2分別為兩個(gè)空間相鄰區(qū)域的光譜方差；obj1，obj2分別為兩個(gè)空間相鄰區(qū)域的面積。

形狀異質(zhì)度變化是用合并前后的差值來(lái)評(píng)估的，所以光滑度smooth和緊致度compact表示為

式中merge表示合并后的對(duì)象實(shí)際邊長(zhǎng)，obj1，obj2分別表示兩個(gè)空間相鄰對(duì)象的實(shí)際邊長(zhǎng)；merge為合并后對(duì)象外切矩形的周長(zhǎng)；obj1，obj2分別表示兩個(gè)空間相鄰對(duì)象的外切矩形的周長(zhǎng)。

1.3 影像的分類特征的選取

Terra SAR-X全極化數(shù)據(jù)具有豐富的極化特征信息，這些特征信息可間接反映地物的組成和結(jié)構(gòu)信息，是地物分類的重要依據(jù)。本文采用Freeman極化分解[12]、Yamaguchi極化分解[13]和Cloude-Pottier極化分解[14]三種極化分解方法提取了勐臘研究區(qū)和普洱研究區(qū)各10種常見(jiàn)的極化特征參數(shù)（見(jiàn)表2），作為地物覆蓋分類的備選分類特征，這些特征參數(shù)可以很好地表達(dá)地物的散射特性[15]。圖2、圖3分別為勐臘研究區(qū)和普洱研究區(qū)極化分解彩色合成圖?；叶裙采仃嚕℅ray Level Cooccurrence Matrix，GLCM）是提取紋理特征的常用方法[16]。紋理特征能夠反映影像中的同質(zhì)性區(qū)域，突出地物之間的差異性，在地物分類中起到重要作用。Terra SAR-X數(shù)據(jù)屬于短波SAR高分辨率數(shù)據(jù)，空間分辨率可達(dá)2.1m，具有豐富的紋理結(jié)構(gòu)特征。因?yàn)镾AR數(shù)據(jù)本身是復(fù)數(shù)形式，不能直接提取紋理，所以將Terra SAR數(shù)據(jù)轉(zhuǎn)換為Span（SAR數(shù)據(jù)的總功率圖）影像，再基于Span利用GLCM提取常用的6種紋理信息（見(jiàn)表2）。

表2 特征參數(shù)

Tab.2 Characteristic parameters

圖2 勐臘極化分解彩色合成圖

圖3 普洱極化分解彩色合成圖

1.4 最優(yōu)特征選擇

在同樣的分類器下特征的選取決定了分類精度，特征選擇的目的就是求得一組最有效的分類特征，所謂有效是指特征維數(shù)減少到同一水平時(shí)，分類性能最佳。因此需要有定量分析比較方法，判斷所得到的特征維數(shù)及所使用特征是否對(duì)分類最有利，這種用以定量檢驗(yàn)分類性能的準(zhǔn)則稱為類別可分離性判據(jù)，即類內(nèi)類間距離。

本文采用的特征選擇算法是在給定的初始特征集下（表2中的16類特征集合），利用選取的各個(gè)類別的樣本，尋找到不同類別的樣本之間的類內(nèi)類間距離最大的特征組合。在初始的特征集合中分析所有的特征組合，即1到16（特征個(gè)數(shù)）的所有組合，即在選擇1個(gè)特征或最多選擇16個(gè)特征時(shí)都可以用()計(jì)算出類內(nèi)類間距離，對(duì)每一個(gè)特征組合利用式（4）計(jì)算類內(nèi)類間距離()，()代表了樣本間的可分性，所以當(dāng)()最大時(shí)的特征組合即為最優(yōu)特征組合。

式中為總類別數(shù)；r,r分別為類及類中樣本數(shù)；p,p是相應(yīng)類別的先驗(yàn)概率。由式（4）計(jì)算得到的類內(nèi)類間距離（分離距離）。

1.5 模糊分類與H/A/α-Wishart分類

（1）模糊分類

本文在前面影像分割和特征提取的基礎(chǔ)上結(jié)合模糊分類方法下的最鄰近分類器進(jìn)行影像分類。通過(guò)提取SAR影像的極化特征和紋理特征信息，建立分類隸屬度規(guī)則集。面向?qū)ο蟮淖钹徑诸惙椒煞譃?個(gè)步驟：1）確定分類類別，在每一類中選取正確的訓(xùn)練樣本，通過(guò)特征優(yōu)選建立特征空間，設(shè)置函數(shù)斜率和最小隸屬度閾值，最后在特征空間中計(jì)算待分類影像對(duì)象與樣本之間的距離；2）將距離轉(zhuǎn)換為隸屬度（通常在0~1范圍的模糊值），如果待分類影像的隸屬度大于所設(shè)定的隸屬度閾值，則該對(duì)象被分為與它距離最近的樣本所屬的類別中[18]。例如：如果影像對(duì)象的最鄰近樣本對(duì)象屬于紅類樣本，則這個(gè)對(duì)象分為紅類樣本（圖4），所有對(duì)象的隸屬度值由式（5）計(jì)算。

由式（5）中的距離可以計(jì)算出一個(gè)多維的指數(shù)形式的隸屬度函數(shù)

式中()為對(duì)象的隸屬度；f為函數(shù)斜率。

()的減少量由參數(shù)決定，可以用不同的函數(shù)斜率來(lái)表達(dá)（見(jiàn)圖5）。

圖4 最鄰近分類原則

圖5 最鄰近隸屬度函數(shù)

本文使用函數(shù)的缺省值為0.2。最鄰近分類方法可以很好地評(píng)價(jià)對(duì)象特征間的相關(guān)性；可以有效解決二維空間中隸屬度函數(shù)的類描述的重疊問(wèn)題；可以快速處理分類的類層次結(jié)構(gòu)。

（2）H/A/α-Wishart分類

1996年Cloude和Pottier提出的Cloude-Pottier分解方法[14]，定義了散射熵、反熵和平均散射角三個(gè)物理量。Cloude-Pottier分解方法利用和組建一個(gè)–二維特征平面，并把該平面依據(jù)熵值高低和散射類型劃為8個(gè)有效區(qū)域，每個(gè)區(qū)域代表一類，可以依據(jù)和的值進(jìn)行簡(jiǎn)單的非監(jiān)督分類。

1999年Lee提出Wishart分布可以更好地描述自然分布的散射機(jī)制[3]。L Ferro-Famil提出將反熵引入分類，的值設(shè)定為0.5，反熵的值大于0.5為一類，小于0.5為另一類[19]。此時(shí)聚類中心由8類分為16類，然后進(jìn)行Wishart迭代分類，即H/A/α-Wishart分類。此方法是一種簡(jiǎn)單快速的非監(jiān)督分類方法，在迭代過(guò)程中無(wú)法保持地物的散射特性，造成誤分發(fā)生，導(dǎo)致分類精度低。

晞月原本聲音柔美，一哭起來(lái)愈加清婉悠亮，頗有一唱三嘆之效，十分哀戚。連遠(yuǎn)遠(yuǎn)站在外頭伺候的雜役小太監(jiān)們，亦不覺(jué)心酸起來(lái)。

2 試驗(yàn)結(jié)果

勐臘與普洱兩塊試驗(yàn)區(qū)試驗(yàn)過(guò)程完全一致。試驗(yàn)過(guò)程簡(jiǎn)述如下：1）在ENVI軟件下進(jìn)行研究區(qū)的裁剪，并進(jìn)行了圖像濾波（Lee濾波，3×3窗口大?。?；2）利用Polsarpro軟件提取極化信息和Span影像；3）利用eCognition進(jìn)行影像的多尺度分割，并在eCognition提取出紋理特征，選取訓(xùn)練樣本，然后對(duì)所有波段進(jìn)行特征空間優(yōu)化，最后在eCognition自帶的模糊分類器下進(jìn)行分類；4）H/A/α-Wishart分類在Polsarpro軟件下提取出熵、各項(xiàng)異性和平均散射角，在Wishart分類器下進(jìn)行迭代分類。

2.1 分割結(jié)果

考慮地物梯度、色彩、紋理、灰度級(jí)等特性，經(jīng)過(guò)多次分割試驗(yàn)，對(duì)兩個(gè)試驗(yàn)區(qū)的分割參數(shù)進(jìn)行設(shè)置，表3為最優(yōu)分割結(jié)果。圖6和圖7分別是兩個(gè)研究區(qū)分割結(jié)果，分割邊界為黑色線條。為驗(yàn)證分割邊界與實(shí)際地物的吻合度，在兩個(gè)試驗(yàn)區(qū)中用分割邊界疊加到光學(xué)影像對(duì)比顯示。

表3 分割參數(shù)設(shè)置

Tab.3 Segmentation parameter settings

圖6 勐臘研究區(qū)分割結(jié)果

圖7 普洱研究區(qū)分割結(jié)果

對(duì)比圖6（a）和圖6（b）可以看出，河流、裸地、農(nóng)田、林地的分割邊界較準(zhǔn)確，圖右下角建筑區(qū)的分割效果較差，這是由于建筑區(qū)組成復(fù)雜，例如，建筑物之間有樹(shù)木、道路以及空地等，影像的異質(zhì)性較大，使得局部分割結(jié)果不連續(xù)。由圖6（a）和圖6（b）的細(xì)節(jié)圖可以看出，SAR影像分割邊界與光學(xué)影像顯示的地物邊界吻合效果很好。

對(duì)比圖7（a）和圖7（b）可以看出，建筑用地、機(jī)場(chǎng)、水田、裸地、農(nóng)地、林地的邊界都有明顯的區(qū)分，由于分辨率的限制，道路未能分出，與建筑用地分為了一類。由圖7（a）和圖7（b）的細(xì)節(jié)圖可以看出，在機(jī)場(chǎng)的南端，建筑物和機(jī)場(chǎng)空地間的邊界清晰準(zhǔn)確，分割結(jié)果同樣符合光學(xué)影像的邊界。普洱研究區(qū)的地塊相較于勐臘研究區(qū)更小，采用的分割尺度比勐臘研究區(qū)的分割尺度要大，所以分割的更細(xì)，對(duì)象也更多。

2.2 特征優(yōu)選結(jié)果

勐臘和普洱研究區(qū)分離距離與特征量維數(shù)的相關(guān)關(guān)系見(jiàn)圖8、圖9。勐臘研究區(qū)通過(guò)對(duì)表2中16個(gè)備選分類特征的優(yōu)化，得到最佳維數(shù)為9，最大類內(nèi)類間分離距離0.302 2（圖8）。9類特征分別為：Freeman分解的偶次散射p，Yamaguchi分解的體散射f、螺旋體散射f，Cloude-Pottier分解的散射熵、反熵、平均散射角，灰度共生矩陣的差異性（Dissimilarity）、同質(zhì)性（Homogeneity）、均值（Mean）。普洱研究區(qū)通過(guò)對(duì)表2中16個(gè)備選分類特征的優(yōu)化，得到最佳維數(shù)為7，最大類內(nèi)類間距離為0.288 4（圖9）。7類特征分別為：Cloude-Pottier分解的散射熵、平均散射角，Yamaguchi分解的體散射f、螺旋體散射f，F(xiàn)reeman分解的體散射p，灰度共生矩陣的差異性（Dissimilarity）、同質(zhì)性（Homogeneity）。

圖8 勐臘研究區(qū)分離距離與特征量維數(shù)

Fig.8 Separation distance and feature dimension of Mengla study area

圖9 普洱研究區(qū)分離距離與特征量維數(shù)

2.3 基于面向?qū)ο蟮哪：诸惤Y(jié)果

圖10為勐臘研究區(qū)的分類樣本及分類結(jié)果，共6種地類：天然林地4個(gè)，橡膠林地3個(gè)，農(nóng)地4個(gè)，水體4個(gè)，建筑用地5個(gè)、裸地4個(gè)。圖11為普洱研究區(qū)的分類樣本及分類結(jié)果，共6種地類：機(jī)場(chǎng)3個(gè)、農(nóng)地3個(gè)、水田3個(gè)、林地3個(gè)、城區(qū)4個(gè)、裸地2個(gè)。試驗(yàn)驗(yàn)證數(shù)據(jù)采用目視解譯光學(xué)影像與實(shí)地考察相結(jié)合得到的ROI，并利用混淆矩陣進(jìn)行精度評(píng)價(jià)。勐臘及普洱研究區(qū)的精度評(píng)價(jià)結(jié)果分別如表4和表5所示。

圖10 勐臘研究區(qū)分類樣本及分類結(jié)果

表4 勐臘研究區(qū)精度評(píng)價(jià)表

Tab.4 The accuracy evaluation for Mengla study area

表5 普洱研究區(qū)精度評(píng)價(jià)表

Tab.5 The accuracy evaluation for Pu￠er study area

2.4 分類結(jié)果與分析

勐臘和普洱研究區(qū)總體分類精度分別達(dá)到了88.5%和86.9%，比H/A/α-Wishart分類方法提高了43.3%和44.6%。每一地類的分類效果也有很大的改善，在勐臘研究區(qū)中，相比H/A/α-Wishart分類方法在制圖精度方面：裸地提高了45.03%，水體提高了23.20%，建筑用地提高了36.64%，農(nóng)地提高了41.83%，橡膠林地提高了66.4%，天然林地提高了39.6%。普洱研究區(qū)的分類精度同樣有很大提高，詳見(jiàn)表5。此外，在勐臘研究區(qū)的分類結(jié)果中，H/A/α-Wishart分類方法將林區(qū)邊界的陰影部分錯(cuò)分為水體，這是由于單純的依靠散射機(jī)制不能有效地對(duì)陰影進(jìn)行分類，而本文方法則成功地解決了這個(gè)問(wèn)題。定性目視分析和定量精度評(píng)價(jià)表明，本文方法比H/A/α-Wishart分類方法有了較大提高，可以實(shí)現(xiàn)高精度的土地覆蓋分類。

3 結(jié)束語(yǔ)

以TerraSAR-X的X波段全極化數(shù)據(jù)為信息源，經(jīng)過(guò)微波遙感特征的提取和優(yōu)化結(jié)合面向?qū)ο蟮哪：诸惙椒?，開(kāi)展了云南西雙版納州勐臘縣和普洱市思茅區(qū)（均為局部）土地覆蓋分類研究，研究結(jié)果表明：

1）本文特征優(yōu)選方法在保證分類精度的同時(shí)，去除了冗余的分類特征，在準(zhǔn)確地物邊界分割和特征優(yōu)選結(jié)果下，結(jié)合模糊分類，有效避免了噪聲在分類中的影響，實(shí)現(xiàn)了SAR影像高精度的土地覆蓋分類；

2）兩個(gè)不同地區(qū)的試驗(yàn)結(jié)果中，散射熵、平均散射角、體散射、差異性、同質(zhì)性等微波遙感特征在土地覆蓋分類中發(fā)揮了重要作用，這一結(jié)果對(duì)同類型的微波遙感數(shù)據(jù)特征提取和優(yōu)化具有一定的參考意義；

3）研究中兩個(gè)不同氣候類型的研究區(qū)的分類精度都能達(dá)到85%以上，說(shuō)明本文的分類方法具有較高的可移植性，對(duì)其他地區(qū)開(kāi)展微波遙感土地覆蓋分類具有較高的借鑒價(jià)值，同樣具有很好的分類效果。

以下問(wèn)題有待進(jìn)一步開(kāi)展研究：1）如何確定最優(yōu)分割尺度，并在最優(yōu)分割結(jié)果的基礎(chǔ)上，用Terra SAR-X/TanDEM-X數(shù)據(jù)進(jìn)行極化干涉處理提取森林冠層高度和相位信息作為分類特征信息；2）如何提取反映對(duì)象形狀的特征參數(shù)并用于分類，以期得到更好的分類效果；3）研究發(fā)現(xiàn)特征波段的選取與分類樣本的依賴性較高，如何選取最優(yōu)的分類樣本是今后研究中有待深入研究的新問(wèn)題。

[1] 趙磊. 基于譜圖分割的極化SAR影像面向?qū)ο蠓诸惙椒ㄑ芯縖D]. 北京: 中國(guó)林業(yè)科學(xué)研究院, 2014. ZHAO Lei. Object-based Classification of Polarmetric SAR Images Based on Spectral Graph Partitioning[D]. Beijing: Chinese Academy of Forestry, 2014. (in Chinese)

[2] CLOUDE S R, POTTIER E. An Entropy Based Classification Scheme for Land Applications of Polarimetric SAR[J]. IEEE Transactions on Geoscience & Remote Sensing, 1997, 35(1): 68-78.

[3] LEE J S, GRUNES M R, AINSWORTH T L, et al. Unsupervised Classification Using Polarimetric Decomposition and the Complex Wishart Classifier[J]. IEEE Transactions on Geoscience & Remote Sensing, 2002, 37(5): 2249-2258.

[4] LEE J S, GRUNES M R, POTTIER E, et al. Unsupervised Terrain Classification Preserving Polarimetric Scattering Characteristics[J]. Geoscience & Remote Sensing IEEE Transactions on, 2004, 42(4): 722-731.

[5] PUTIGNANO C, SCHIAVON G, SOLIMINI D, et al. Unsupervised Classification of A Central Italy Landscape by Polarimetric L-band SAR Data[C]//IEEE International Geoscience and Remote Sensing Symposium. IEEE, 2005: 1291-1294.

[6] LUMSDON P, CLOUDE S R, WRIGHT G. Polarimetric Classification of Land Cover for Glen Affric Radar Project[J]. IEEE Proceedings–Radar Sonar and Navigation, 2005, 152(6): 404-412.

[7] 趙一博, 秦先祥, 鄒煥新. 基于目標(biāo)分解和SVM的極化SAR圖像分類方法[J]. 航天返回與遙感, 2013, 34(2): 50-56.ZHAO Yibo, QIN Xianxiang, ZOU Huanxin. Classification of Polarimetric SAR Image Based on Target Decomposition and SVM[J]. Spacecraft Recovery & Remote Sensing, 2013, 34(2): 50-56. (in Chinese)

[8] CHEN J, DU P, TAN K. An Improved Object Oriented Classification Scheme for Polarimetric SAR Image[J]. Journal of China University of Mining & Technology, 2015, 44(5): 944-951.

[9] LI H T, GU H Y. Object-oriented Classification of Polarimetric SAR Imagery Based on Statistical Region Merging and Support Vector Machine[C]//International Workshop on Earth Observation and Remote Sensing Applications. IEEE, 2008: 1-6.

[10] STEWART D, BLACKNELL D, BLAKE A, et al. Optimal Approach to SAR Image Segmentation and Classification[J]. Radar, Sonar and Navigation, IEEE Proceedings, 2000, 147(3): 134-142.

[11] 關(guān)元秀, 程曉陽(yáng). 高分辨率衛(wèi)星影像處理指南[M]. 北京: 科學(xué)出版社, 2008: 203-207. GUAN Yuanxiu, CHENG Xiaoyang. High Resolution Satellite Image Processing Guide[M]. Beijing: Science Press, 2008: 203-207. (in Chinese)

[12] FREEMAN A, DURDEN S L. A Three-component Scattering Model for Polarimetric SAR Data[J]. IEEE Transactions on Geoscience & Remote Sensing, 1998, 36(3): 963-973.

[13] YAMAGUCHI Y, MORIYAMA T, ISHIDO M, et al. Four-component Scattering Model for Polarimetric SAR Image Decomposition[J]. IEEE Transactions on Geoscience & Remote Sensing, 2005, 43(8): 1699-1706.

[14] CLOUDE S R, POTTIER E. A Review of Target Decomposition Theorems in Radar Polarimetry[J]. IEEE Transactions on Geoscience & Remote Sensing, 1996, 34(2): 498-518.

[15] YANG H, SHAO Y, YANG H, et al. Full Polarimetric SAR Classification Based on Yamaguchi Decomposition Model and Scattering Parameters[C]//IEEE International Conference on Progress in Informatics and Computing, 2010: 1104-1108.

[16] GAO C C, HUI X W. GLCM-based Texture Feature Extraction[J]. Computer Systems & Applications, 2010, 19(6): 195-198.

[17] 李海峰, 張建州, 游志勝. 基于類距離的可分離性判據(jù)[J]. 計(jì)算機(jī)工程與應(yīng)用, 2003, 39(26): 97-99. LI Haifeng, ZHANG Jianzhou, YOU Zhisheng. A Separative Criterion Based on Class Distance[J]. Computer Engineering and Applications, 2003, 39(26): 97-99. (in Chinese)

[18] 陳金麗. 面向?qū)ο蟮淖钹徑惴ㄑ芯颗c實(shí)現(xiàn)[D]. 北京: 中國(guó)地質(zhì)大學(xué), 2009. CHEN Jinli. Research and Realization of Object Oriented Nearest Neighbour Classification Algorithm[D]. Beijing: China University of Geosciences(Beijing), 2009. (in Chinese)

[19] FERRO-FAMIL L, POTTIER E, LEE J S. Unsupervised Classification of Multifrequency and Fully Polarimetric SAR Images Based on The H/A/Alpha-Wishart Classifier[J]. IEEE Transactions on Geoscience & Remote Sensing, 2002, 39(11): 2332-234.

（編輯：王麗霞）

Object-oriented Land Cover Classification Based on Feature Selection in Quad-polarimetric SAR Images

LU Xiang ZHANG Wanqiu ZHENG Yalan YUE Cairong

（Southwest Forestry University, Kunming 650224, China）

As the quad-polarimetric SAR data has plentiful information, it is difficult to obtain good classification results just by single polarimetric feature and pixel-based classification. Based on the X-band quad-polarimetric radar data of Terra SAR-X in Mengla of Xishuangbanna and Simao Pu'er city, Yunnan province, the object-oriented land cover classification experiments based on feature selection are carried out. Firstly, the quad-polarimetric SAR data is preprocessed. The Pauli RGB image is extracted from the study area data, and then is segmented as the basic unit of classification using the image segmentation technique. Then, the polarimetric features and the texture features of Span images are extracted from the SAR images, and the optimal feature sets are selected. Finally, the object-oriented fuzzy classification method is used into classification experiments, and the classification results are evaluated by field survey data. The experimental results show that object-oriented method can effectively remove the impact of noise, and the optimal combination of feature bands makes the classification results more accurate. The overall classification accuracy is up to 88.5% for Mengla in Xishuangbanna, and is 86.8% in Simao of Pu'er city. Compared with the H/A/alpha -Wishart classification method, the accuracy is improved by more than 40%.

quad-polarimetric; synthetic aperture radar; image segmentation; feature selection; fuzzy classification; remote sensing application

TP751.1

A

1009-8518(2018)02-0093-11

10.3969/j.issn.1009-8518.2018.02.011

2017-08-28

陸翔，男，1992年生，2015年獲滁州學(xué)院地理信息系統(tǒng)專業(yè)學(xué)士學(xué)位，現(xiàn)在西南林業(yè)大學(xué)地圖學(xué)與地理信息系統(tǒng)專業(yè)攻讀碩士學(xué)位。研究方向是微波遙感應(yīng)用。E-mail：1821332789@qq.com。岳彩榮，男，1964年生，博士，教授，博士生導(dǎo)師，研究方向?yàn)橘Y源環(huán)境遙感，微波遙感。E-mail：cryue@163.com。