黃鵬艷 范永良

(1.洛陽(yáng)理工學(xué)院土木工程學(xué)院，河南洛陽(yáng) 471000; 2.中航光電科技股份有限公司，河南洛陽(yáng) 471000)

0 引言

合成孔徑雷達(dá)(SAR)技術(shù)與系統(tǒng)已從單波段、單極化逐步發(fā)展到多波段、全極化［1］。全極化SAR(PolSAR)圖像中蘊(yùn)含的信息豐富，利用從圖像中獲取的特征向量對(duì)PolSAR圖像進(jìn)行分類(lèi)逐漸成為研究熱點(diǎn)，而最常見(jiàn)的方式是利用極化特征向量進(jìn)行分類(lèi)［2］。通過(guò)目標(biāo)分解的方法對(duì)圖像的散射數(shù)據(jù)提取參數(shù)構(gòu)成極化特征向量［3］。此外，PolSAR圖像中蘊(yùn)含著豐富的紋理信息［4］，紋理特征是一種重要的視覺(jué)特征，可以表現(xiàn)地物目標(biāo)在結(jié)構(gòu)排列和地表粗糙度上的一些特點(diǎn)，極化特性相似的地物可通過(guò)不同的紋理特征得到進(jìn)一步劃分。對(duì)獲取的極化SAR灰度圖提取紋理特征參數(shù)構(gòu)成紋理特征向量，與提取的極化特征向量相結(jié)合用于PolSAR圖像分類(lèi)。因此，本文研究加入紋理特征的PolSAR圖像分類(lèi)，給出了兩種特征向量的提取和組合過(guò)程，并采用支持向量機(jī)(SVM)方法進(jìn)行分類(lèi)［5，6］。本文以湖北咸寧地區(qū)的PolSAR數(shù)據(jù)為例進(jìn)行分類(lèi)實(shí)驗(yàn)，對(duì)在不同特征向量下分類(lèi)精度進(jìn)行分析，以證明紋理特征與極化特征結(jié)合用于極化SAR圖像分類(lèi)可以提高PolSAR圖像分類(lèi)精度。

1 算法原理

地物的極化特征是基于目標(biāo)散射屬性的特征，能與目標(biāo)的結(jié)構(gòu)形狀特征很好地聯(lián)系，可通過(guò)目標(biāo)分解［3］的方法得到用于Pol-SAR圖像分類(lèi)。極化特征是PolSAR圖像分類(lèi)最常使用的特征，但它也有一定的限制。極化信息是基于像素級(jí)別的描述，不能利用圖像的統(tǒng)計(jì)知識(shí)信息，限制了相似的地表類(lèi)別的區(qū)分能力，所以在反映圖像的空間分布特征和細(xì)節(jié)表現(xiàn)上能力有限。紋理特征是一種視覺(jué)特征［4］，可以表現(xiàn)地物目標(biāo)在結(jié)構(gòu)排列和地表粗糙度等表面信息。從圖像中提取紋理信息并將其與極化特征相結(jié)合用于圖像分類(lèi)中，使得圖像中蘊(yùn)含著的極化、紋理不同類(lèi)型的參數(shù)綜合發(fā)揮分類(lèi)功能，提供區(qū)別多種地物類(lèi)別的多方面分類(lèi)信息，提高分類(lèi)精度。

1.1 極化特征提取

極化目標(biāo)分解是基于切合實(shí)際的物理約束解譯目標(biāo)的散射機(jī)制，不同的目標(biāo)分解方法從不同角度對(duì)基礎(chǔ)散射機(jī)制進(jìn)行物理解釋。Pauli分解是一種經(jīng)典的相干極化目標(biāo)分解方法，可得到|a|2，|b|2，|c|2特征參數(shù)，分別表示單次散射、二面角散射、不對(duì)稱分量的散射能量。Krogager分解得到ks，kd，kh特征參數(shù)，分別表示球、二面角、螺旋體的散射能量。Cloude-Pottier分解［7］使用最廣泛，可得到如下參數(shù):極化散射參數(shù)珔α，珔β，珋δ，珔γ;各項(xiàng)異性度A提供了相對(duì)較弱的兩個(gè)散射分量之間的大小關(guān)系;極化散射熵H描述各種不同散射類(lèi)型在統(tǒng)計(jì)意義上的無(wú)序性;λ1，λ2，λ3表示各散射機(jī)制的權(quán)重。Yamaguchi分解得到 Ps，Pd，Pv，Pc四個(gè)參數(shù)分別對(duì)應(yīng)表面散射、二次散射、體散射、螺旋體散射的功率，對(duì)人造目標(biāo)有獨(dú)特的識(shí)別能力。綜合利用Pauli，Krogager，Cloude-Pottier，Yamaguchi四種目標(biāo)分解方法得到的參數(shù)構(gòu)成極化特征向量，如式(1)所示。

1.2 紋理特征提取

提取紋理特征值的方法有多種，本文采用灰度共生矩陣法［8］。紋理特征參數(shù)較多，其中，對(duì)比度反映圖像局部灰度變化，水體和建筑物與周?chē)匚锘叶茸兓黠@，對(duì)比度圖像中表現(xiàn)為較亮區(qū)域;相異性表征紋理之間的差異性，能將地物的邊緣信息提取出來(lái);信息熵的大小表示圖像紋理的復(fù)雜度，建筑物和水體的紋理較為簡(jiǎn)單，熵值較小，熵圖像表現(xiàn)為較暗;協(xié)同性能識(shí)別局部的相似紋理信息;均值反映圖像灰度的均勻性，水域、道路、建筑物區(qū)域的灰度較為均勻，在均值圖像上表現(xiàn)為較高的亮度;一致性反映圖像局部紋理的一致性，水域的紋理一致性最為突出。分別用 CON，DIS，ENT，HOM，MEA，UNI來(lái)表示以上六個(gè)紋理特征參數(shù)。將六個(gè)紋理特征參數(shù)組成紋理特征向量如式(2)所示。

1.3 基于極化和紋理特征向量的SVM分類(lèi)

1)對(duì)PolSAR數(shù)據(jù)進(jìn)行濾波處理以降低相干斑噪聲的影響，同時(shí)保持圖像的邊緣信息和細(xì)節(jié)信息。

2)對(duì)濾波處理后的PolSAR數(shù)據(jù)進(jìn)行多種極化目標(biāo)分解，組成極化特征向量v。

3)提取紋理特征值，組成紋理特征向量vt。

4)極化特征向量與紋理特征向量進(jìn)行組合，構(gòu)成新的特征向量v+vt。

5)分別利用特征向量v和v+vt進(jìn)行分類(lèi)，并對(duì)比分類(lèi)精度。

2 分類(lèi)實(shí)驗(yàn)

2.1 實(shí)驗(yàn)數(shù)據(jù)

本文選用湖北咸寧地區(qū)在2013年4月12日采集的全極化數(shù)據(jù)進(jìn)行實(shí)驗(yàn)，數(shù)據(jù)采集傳感器為 Radarsat-2，為 C波段(37.5 mm～75 mm)成像的SLC數(shù)據(jù)，距離向分辨率為4.73 mm，方位向分辨率為4.96 mm。實(shí)驗(yàn)圖像為HH，HV，VH，VV四個(gè)極化通道的灰度圖，提取紋理參數(shù)需提供的總功率SPAN圖像。為對(duì)比分析分類(lèi)精度還提供同一研究區(qū)域的分辨率較高的光學(xué)影像圖。四個(gè)極化通道灰度圖像、總功率灰度圖、光學(xué)影像圖如圖1所示。

圖1 實(shí)驗(yàn)數(shù)據(jù)

在全極化SAR圖像中，純凈的水體介電常數(shù)非常高，幾乎不能被電磁波穿透，研究區(qū)域風(fēng)浪一般較小，粗糙度小，平靜的水面易發(fā)生鏡面反射，圖像上表現(xiàn)為較暗的區(qū)域。道路的特點(diǎn)為質(zhì)地堅(jiān)硬，表面平坦光滑，細(xì)長(zhǎng)型，散射機(jī)制為鏡面反射，標(biāo)準(zhǔn)高度差小，圖像上表現(xiàn)為較暗區(qū)域。建筑物為強(qiáng)散射體，表現(xiàn)為偶次散射，圖像上表現(xiàn)較亮。高矮不同的植被表現(xiàn)為體散射。裸地由于受到含水量和表面粗糙度的影響，散射過(guò)程較為復(fù)雜。本文中將試驗(yàn)區(qū)域的地物類(lèi)別分為五類(lèi)，即水域、道路、建筑、植被、裸地。

2.2 提取特征向量

在3×3窗口下采用極化LEE濾波對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行濾波去噪，之后進(jìn)行極化目標(biāo)分解，按照式(1)構(gòu)成極化特征向量v。對(duì)總功率SPAN灰度圖，按照式(2)構(gòu)成紋理特征向量vt。極化特征向量v與紋理特征向量vt結(jié)合構(gòu)成組合特征向量v+vt。

2.3 分類(lèi)結(jié)果及對(duì)比分析

使用SVM分類(lèi)方法，分別利用v和v+vt對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分類(lèi)，結(jié)果如圖2所示。

圖2 分類(lèi)結(jié)果

與分辨率較高的光學(xué)影像進(jìn)行對(duì)比可以看出，僅利用極化特征向量v的分類(lèi)結(jié)果(見(jiàn)圖2a))，五種地物類(lèi)別區(qū)域雖大致得到劃分，但區(qū)分度較差。加入紋理特征向量v+vt的分類(lèi)效果(如圖2b)所示)較好，地物區(qū)分明顯，水域邊緣、道路等細(xì)節(jié)信息清晰。對(duì)分類(lèi)結(jié)果進(jìn)行定量分析，計(jì)算混淆矩陣得到分類(lèi)精度，利用v和v+vt的分類(lèi)精度對(duì)比如圖3所示。

極化特征向量v的分類(lèi)精度為88.60%，加入紋理特征的特征向量v+vt的分類(lèi)精度為95.87%，精度增長(zhǎng)7.27%。對(duì)于道路，僅使用極化信息難以識(shí)別道路這種細(xì)長(zhǎng)的地物，分類(lèi)精度僅為72.92%，加入紋理特征參數(shù)后，道路的分類(lèi)精度提高為90.97%;植被散射過(guò)程復(fù)雜，裸地與道路、建筑物具有相似的散射性，僅利用極化特征分類(lèi)時(shí)植被和裸地的精度較低，分別為87.65%和89.36%，加入紋理特征向量后，植被和裸地的分類(lèi)精度提高為96.61%和97.08%;水域和建筑物是強(qiáng)散射體，僅利用極化信息已經(jīng)能將其很好地分類(lèi)識(shí)別，紋理特征的加入使其分類(lèi)精度小幅度提升。

圖3 分類(lèi)精度對(duì)比

3 結(jié)語(yǔ)

1)本文通過(guò)將紋理特征向量和極化特征向量相結(jié)合，進(jìn)行PolSAR圖像分類(lèi)方法的研究，闡述了算法原理并進(jìn)行了實(shí)驗(yàn)驗(yàn)證。實(shí)驗(yàn)結(jié)果表明，與僅利用極化特征向量相比，紋理特征向量的加入充分挖掘了隱藏在灰度圖像上的信息，改善了散射特性類(lèi)似的地物之間的區(qū)別能力，使得總分類(lèi)精度提升了7.27%。

2)紋理特征的加入對(duì)各類(lèi)地物的識(shí)別均有促進(jìn)作用。水域和建筑物這種強(qiáng)散射體僅利用極化特征的分類(lèi)精度已經(jīng)很高，紋理特征的加入使得這兩者的分類(lèi)精度小幅提升，而對(duì)散射特性不明顯的道路，植被，裸地，紋理特征的加入使得其分類(lèi)精度顯著增加。