基于超像素和全卷積網絡的極化SAR圖像分類

陳彥橋，張小龍，陳金勇，高 峰，柴興華

(1.中國電子科技集團公司第五十四研究所，河北 石家莊 050081；2.中國電子科技集團公司航天信息應用技術重點實驗室，河北 石家莊 050081)

0 引言

極化合成孔徑雷達(Synthetic Aperture Rader，SAR)圖像分類已經成為了遙感處理領域最重要的應用之一[1]。這些年，一系列機載和星載平臺提供了大量的極化SAR數據，極化SAR圖像研究也變得越來越重要[2]。許多相關算法也被提出用于極化SAR圖像分類任務，在這些方法中，許多方法是基于物理散射機制的極化目標分解方法，包括Cameron分解[3]、Cloude分解[4]、Krogager分解[5]、Freeman分解[6]、Huynen分解[7]、Yamaguchi分解[8]和Van Zyl分解[9]等方法。除此之外，極化SAR相關學者已經證明極化協方差矩陣(C)和極化相干矩陣(T)都服從復Wishart分布[10]。因此，極化SAR的數據分布在極化SAR圖像分類中扮演著很重要的角色。需要補充的是，極化SAR數據的散射機制和統計特性也常常在極化SAR圖像分類方法中同時使用[11]。近些年，機器學習在極化SAR圖像中也變得越來越重要。直至今日，許多機器學習方法用于極化SAR圖像分類，包括馬爾科夫隨機場(Markov Random Field，MRF)[12]、K近鄰(K-Nearest Neighbor，KNN)[13]、神經網絡(Neural Network，NN)[14]、隨機森林分類器(Random Forest Classifier，RFC)[15]和支持向量機(Support Vector Machine，SVM)[16]。

近些年，深度學習在極化SAR圖像分類中得到了應用，包括深度置信網(Deep Belief Network，DBN)[17]、稀疏自編碼(Sparse Auto-Encoder，SAE)[18]、卷積神經網絡(Convolution Neural Network，CNN)[19]和全卷積神經網絡(Fully Convolution Network，FCN)[2]。在上述方法中，FCN是最適合極化SAR圖像分類的架構，因為FCN可以利用圖像空間信息，而且是端到端、點到點的架構。因此，本文使用FCN作為基礎框架。

在之前的基于FCN的極化SAR圖像分類方法中，并沒有使用超像素，因為FCN已經可以使用圖像的空間信息。然而，在FCN沒有取得較為理想的分類結果時，超像素仍然可以在一定程度上修正分類結果。本文首先使用FCN得到初步分類結果，然后利用超像素對分類結果進行修正，最后使用FCN的分類結果中置信度較高的部分修正分類結果。

1 特征提取

1.1 極化相干矩陣

Sinclair散射矩陣包含全極化信息，極化圖像的單像素點可以通過散射矩陣描述?；谏⑸渚仃?，極化相干矩陣[2]也包含全極化信息。滿足互易定理SHV=SVH的情況下，T矩陣可以用式(2)表示，T矩陣的每個元素的物理意義已經被Huynen證明：

(1)

(2)

需要補充的是，T矩陣也可以表示為：

(3)

式中，a=SHH+SVV；b=SHH-SVV；c=SHV；“*”代表共軛操作。

1.2 H/A/α分解方法

在H/A/α分解模型中，利用T矩陣的特征值分析，給出了目標平均散射矩陣參數的估計。T矩陣的特征值分解[2]可以表示為：

(4)

(5)

2 本文方法

2.1 超像素SLIC算法

基于顏色空間和距離空間，Achanta等人[20]提出了簡單線性迭代聚類(Simple Linear Iterative Cluster，SLIC)算法。SLIC算法通過改進K-means聚類算法生成超像素。假設圖像的像素點數目為N，超像素的數目為K，以步長為S的網格初始化聚類中心；在2S×2S的范圍內計算每個像素點與最近的超像素中心的相似度，并將像素點賦值給最相似的超像素，不斷迭代上述過程直到最終收斂。相似度計算方式可以表示為：

(6)

(7)

(8)

式中，dlab代表顏色空間；dxy代表空間距離；DS代表像素空間距離；lk，ak，bk分別代表第k個像素點在lab顏色空間的值；xk，yk代表第k個像素點的坐標值；S代表步長；m是一個超參數，用于平衡空間距離和顏色距離。

2.2 滑窗全卷積網絡

FCN是CNN模型的改進，用于解決密集預測問題，例如圖像分割。為了在輸出層恢復輸入層的分辨信息，在經典FCN框架中加入了上采樣層。FCN還在上采樣通道和下采樣通道使用跳躍連接解決由于下采樣引起的分辨信息損失。跳躍連接對于上采樣通道從下采樣層恢復細節信息很有幫助。然而，FCN并不能直接應用于極化SAR圖像分類，因為不同尺寸的圖像需要不同的FCN框架，而且較大尺寸的輸入圖像也往往會增大設計框架的難度。為了為所有極化SAR圖像設計一種通用架構，李陽陽等人[2]提出了滑窗全卷積網絡(Sliding Window Fully Convolution Networks，SFCN)。SFCN框架圖如圖1所示。

圖1 SFCN框架圖Fig.1 Frame diagram of SFCN

SFCN的滑窗操作與CNN的滑窗操作類似，通過滑窗操作得到的圖像數目的計算方式為：

num=(ceil((Height-W)/S)+1)×

(ceil((Width-W)/S)+1)，

(9)

式中，ceil代表向上取整函數；Height，Width分別代表圖像的高度和寬度；W，S分別代表滑窗的尺寸和步長；num代表通過滑窗操作獲得圖像的數目。

2.3 超像素+SFCN

本文使用超像素+SFCN的方式對極化SAR圖像進行分類，具體步驟如下：

① 用超像素算法對整幅圖像進行分割；

② 用SFCN對整幅圖像進行分類，得到分類結果，標記為Result-SFCN；

③ 用超像素分割結果修正SFCN的分類結果，具體做法為：隸屬于同一超像素塊的所有像素點出現次數最多的分類類別賦值給該超像素塊所有像素點的分類類別，標記為Result-SLIC；

④ 根據信息熵，取Result-SFCN中置信度高的分類結果，修正Result-SLIC中對應區域的分類結果，并將最終分類結果標記為Result-Entropy，信息熵的計算方法為：

(10)

式中，N代表類別數目；i代表第i個類別；p(i)代表將x分為第i類的概率；h(x)越小，則x的分類結果的置信度越高。

3 實驗與分析

本文使用T矩陣(T11，T22，T33，Re(T12)，Re(T13)，Re(T23)，Im(T12)，Im(T13)，Im(T23))和H/A/α特征(H，A，α，λ1，λ2，λ3)作為原始特征，使用精致Lee濾波[1]抑制相干斑噪聲，使用平均正確率(Average Accuracy，AA)、總體正確率(Overall Accuracy，OA)和Kappa系數作為分類結果評價指標。

3.1 Xi’an圖像分類實驗

Xi’an圖像[2]是由RADARSAT-2衛星獲取的全極化C波段極化SAR圖像，位于中國西安市渭河，包含水域(Water)、草地(Grass)和建筑區(Building)3類地物。該圖像大小為512 pixel×512 pixel，分辨率為8 m×8 m。圖2給出了Xi’an圖像的Pauli偽彩色圖、類標圖和配色方案，該圖像的有標記樣本數目為237 416，本文隨機選取1%的有標記樣本作為訓練集，剩余的有標記樣本作為測試集。

圖2 Xi’an極化SAR圖像Fig.2 PolSAR image of Xi’an

該圖像的分類結果如圖3所示。SFCN得到了不錯的分類結果，而且SLIC也得到了不錯的分割結果，如圖3(c)所示，Result-SFCN經過超像素修正后，分類結果得到提高，如圖3(d)所示，Result-SLIC經過Result-SFCN中置信度較高的分類結果修正后，分類結果得到進一步提高。

圖3 Xi’an分類結果Fig.3 Classification results of Xi’an

表1給出了該圖像的分類正確率。本文方法的AA、OA以及Kappa系數都是所有方法中最高的。

表1 Xi’an圖像分類正確率Tab.1 Classification accuracy of Xi’an image

綜上所述，本文方法在Xi’an圖像的分類有效性得到證明。

3.2 Oberpfaffenhofen圖像分類實驗

Oberpfaffenhofen圖像[2]是由E-SAR獲取的全極化L波段的極化SAR圖像，位于德國Oberpfaffenhofen部分區域，包含開放區域(Open Areas)、森林(Wood Land)和建筑區(Built-up Areas)3類地物。該圖像大小為1 300 pixel×1 200 pixel，分辨率為3 m×2.2 m。圖4給出了該圖像的Pauli偽彩色圖、類標圖和配色方案，該圖像的有標記樣本數目為1 374 298，本文隨機選取0.1%有標記樣本作為訓練集，剩余的有標記樣本作為測試集。

圖4 Oberpfaffenhofen極化SAR圖像Fig.4 PolSAR image of Oberpfaffenhofen

圖5和表2分別給出了該圖像的分類結果和分類正確率。

圖5 Oberpfaffenhofen分類結果Fig.5 Classification results of Oberpfaffenhofen

表2 Oberpfaffenhofen圖像分類正確率

由圖5可以看出，SLIC得到了不錯的分割結果，本文方法得到了最好的分類結果。由表2可以看出，本文方法得到了最高的AA、OA和Kappa系數。

綜上所述，本文方法在Oberpfaffenhofen圖像的分類有效性得以證明。

4 結束語

本文提出了基于超像素和FCN的方法用于極化SAR圖像分類。得益于FCN端到端、點到點的分類架構，其在極化SAR圖像分類中很有前景。在此基礎上，使用超像素去修正分類結果，最后使用信息熵選取FCN分類結果中置信度較高的部分，進一步修正分類結果。本文使用T矩陣和H/A/α分解特征作為輸入特征，Xi’an和Oberpfaffenhofen兩幅極化SAR圖像用于測試方法的分類性能，與SFCN相比，本文方法得到了更好的分類結果。