1.   引言

极化合成孔径雷达(Polarimetric Synthetic Aperture Radar, PolSAR)具有全天时和几乎全天候的工作能力，通过收发极化状态正交的电磁波以获取目标的全极化散射信息^[1]。地物分类是农作物生长监控、农村与城市用地普查、环境监测等应用领域的共性基础问题，也是极化SAR图像理解与解译的重要应用方向。高精度的地物分类结果能够为上述应用领域提供可靠的信息支撑。

通常，提高极化SAR地物分类精度主要有两种途径^[2]。第1种途径专注于极化特征的挖掘与优选，通过精细化的极化散射机理建模与解译，从全极化信息中提取出对不同地物类别具有更强区分度的特征。常用的极化散射机理解译方法有基于特征值分解的方法和基于模型分解的方法。基于这些极化目标分解方法所得到的极化特征参数经常被用于极化SAR地物分类，例如Cloude-Pottier分解所得的极化熵/极化平均角/极化反熵(H/ $\alpha$/A)参数^[3]，Freeman-Durden分解^[4]、Yamaguchi分解^[5]和近年来提出的精细化极化目标分解^[6]所得的各散射机理的散射能量参数(如奇次散射、偶次散射、体散射、螺旋散射等)^[7]。第2种途径则从分类器入手，使用性能更好的分类器，以对现有的极化特征进行充分利用。常用的分类器包括C均值分类器、Wishart分类器、支持向量机(Support Vector Machine, SVM)分类器、随机森林分类器、神经网络分类器以及近来年在诸多领域取得成功应用的以卷积神经网络为代表的深度学习分类方法等^[8–11]。当然，对特征和分类器同时进行优化和优选也是提高极化SAR地物分类精度的有效途径。

在传统基于特征的极化SAR地物分类中，具有旋转不变特性的极化特征参数得到了广泛应用。例如，基于H/ $\alpha$/A和总散射能量SPAN的极化SAR地物分类就是一种常用的分类方法。然而，目标的极化响应与目标和SAR的相对几何关系密切相关。同一目标在不同方位取向下，其后向散射可以是显著不同的。同时，不同目标在某些特定方位取向下，其后向散射又是十分相似的。例如，具有不同方位取向的建筑物与森林等植被就是极化SAR图像解译的难点。这是诸多传统极化目标分解方法存在散射机理解译模糊的重要原因之一，同时也限制了基于旋转不变极化特征参数的传统分类方法所得精度的进一步提升。为避免这种解译模糊，一种思路是构建更精细化的目标散射模型和精细化的极化目标分解方法。而另一种思路则是挖掘利用目标方位取向与其后向散射机理之间的隐含关系。文献[12]提出的统一的极化矩阵旋转理论就是一种代表性的方法。该方法提出了在绕雷达视线的旋转域中理解目标散射特性的新思路，并导出了一系列旋转域极化特征。部分旋转域极化特征参数已经在农作物辨识^[13]、目标对比增强^[12]、人造目标提取^[14]等领域获得了成功应用。

由于这些旋转域极化特征包含有目标在旋转域中隐含的极化散射信息，且与其方位取向具有一定关系。若将它们与传统的旋转不变极化特征参数于H/ $\alpha$/A/SPAN联合作为地物分类特征集，则从极化特征挖掘的角度来看，两类不同的极化特征对于不同地物类别的区分能力势必会形成一定程度的互补，进而使分类精度得到进一步提升。基于这一思路，本文提出了一种结合旋转域极化特征与旋转不变特征H/ $\alpha$/A/SPAN的极化SAR地物分类方法。具体即基于不同地物类别样本集类间距最大的特征优选准则，以部分优选的旋转域极化特征参数与H/ $\alpha$/A/SPAN联合作为地物分类所用特征，并选用性能较为稳定的SVM^[15]作为分类器进行分类处理。由于该分类方法额外使用了目标在方位取向方面的隐含信息，故相较于仅使用旋转不变特征H/ $\alpha$/A/SPAN作为输入的SVM分类器^[10]，其能够达到更优的分类性能表现。

本文第2节简要介绍了统一的极化矩阵旋转理论及其所导出的旋转域极化特征参数；第3节提出结合旋转域极化特征的极化SAR地物分类方法；第4节基于AIRSAR和多时相UAVSAR实测数据开展了地物分类对比实验及分析；第5节总结本文方法并对后续研究工作进行展望。

2.   统一的极化矩阵旋转理论

极化SAR获得的目标全极化信息可以通过极化相干矩阵T表示。满足互易性原理时，极化相干矩阵T可以表示为：

其中， ${{{k}}_{\rm{P}}} = \displaystyle\frac{1}{{\sqrt 2 }}{\left[ {\begin{array}{*{20}{c}}{{S_{{\rm{HH}}}} + {S_{{\rm{VV}}}}}&{{S_{{\rm{HH}}}} - {S_{{\rm{VV}}}}}&{2{S_{{\rm{HV}}}}}\end{array}} \right]^{\rm{T}}}$为Pauli散射矢量。 ${S_{{\rm{HV}}}}$为以垂直极化天线发射并以水平极化天线接收条件下的散射系数， ${{{k}}_{\rm{P}}}$中其它元素可类似定义。 $\left\langle { \ } \right\rangle$表示集合平均。 ${T_{ij}}$则表示极化相干矩阵 ${{T}}$中第i行第j列所对应的元素。

将极化相干矩阵 ${{T}}$绕雷达视线进行旋转处理，则可得到旋转域中极化相干矩阵的表达式为：

其中，旋转矩阵为：

在旋转域中极化相干矩阵 ${{T}}\left( \theta \right)$的每个元素经过相应的数学变换即可被统一地由一个正弦函数进行表征^[12]：

其中，A为振荡幅度，B为振荡中心， $\omega$为角频率， ${\theta _{\rm{0}}}$为初始角度。文献[12]将这4类极化特征参数 $\left\{ {A,B,\omega ,{\theta _0}} \right\}$称为振荡参数集，其完整表征极化相干矩阵的各元素在旋转域中的特性。这样就可以导出一系列旋转域极化特征参数，如表1所示。其中， ${\rm{Angle}}\left\{ a \right\}$表示复数a的相位，相应取值范围为 $\left[ { - {{π}} ,{{π}} } \right]$。

表  1  旋转域极化特征参数^[12]

Table  1.  Polarimetric feature parameters derived from rotation domain^[12]

散射矩阵元素项	$A = \sqrt \bullet$	B	$\omega$	${\theta _0} = \frac{1}{\omega }{\rm{Angle}}\left\{ \bullet \right\}$
${\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$	${{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{13}}} \right]$	0	2	${\mathop{\rm Re}\nolimits} \left[ {{T_{13}}} \right] + j{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}} \right]$
${\mathop{\rm Re}\nolimits} \left[ {{T_{13}}\left( \theta \right)} \right]$	${{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{13}}} \right]$	0	2	$- {\mathop{\rm Re}\nolimits} \left[ {{T_{12}}} \right] + j{\mathop{\rm Re}\nolimits} \left[ {{T_{13}}} \right]$
${\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$	${{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{13}}} \right]$	0	2	${\mathop{\rm Im}\nolimits} \left[ {{T_{13}}} \right] + j{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}} \right]$
${\mathop{\rm Im}\nolimits} \left[ {{T_{13}}\left( \theta \right)} \right]$	${{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{13}}} \right]$	0	2	$- {\mathop{\rm Im}\nolimits} \left[ {{T_{12}}} \right] + j{\mathop{\rm Im}\nolimits} \left[ {{T_{13}}} \right]$
${\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$	$\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]$	0	4	$\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right) + j{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right]$
${T_{22}}\left( \theta \right)$	$\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]$	$\frac{1}{2}\left( {{T_{22}} + {T_{33}}} \right)$	4	${\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right] + j\frac{1}{2}\left( {{T_{22}} - {T_{33}}} \right)$
${T_{33}}\left( \theta \right)$	$\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]$	$\frac{1}{2}\left( {{T_{22}} + {T_{33}}} \right)$	4	$- {\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right] + j\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right)$
${\left| {{T_{12}}\left( \theta \right)} \right|^2}$	${\rm R}{{\rm e}^2}\left[ {{T_{12}}T_{13}^*} \right] + \frac{1}{4}{\left( {{{\left| {{T_{13}}} \right|}^2} - {{\left| {{T_{12}}} \right|}^2}} \right)^2}$	$\frac{1}{2}\left( {{{\left| {{T_{12}}} \right|}^2} + {{\left| {{T_{13}}} \right|}^2}} \right)$	4	${\mathop{\rm Re}\nolimits} \left[ {{T_{12}}T_{13}^*} \right] + j\frac{1}{2}\left( {{{\left| {{T_{12}}} \right|}^2} - {{\left| {{T_{13}}} \right|}^2}} \right)$
${\left| {{T_{13}}\left( \theta \right)} \right|^2}$	${\rm R}{{\rm e}^2}\left[ {{T_{12}}T_{13}^*} \right] + \frac{1}{4}{\left( {{{\left| {{T_{13}}} \right|}^2} - {{\left| {{T_{12}}} \right|}^2}} \right)^2}$	$\frac{1}{2}\left( {{{\left| {{T_{12}}} \right|}^2} + {{\left| {{T_{13}}} \right|}^2}} \right)$	4	$- {\mathop{\rm Re}\nolimits} \left[ {{T_{12}}T_{13}^*} \right] + j\frac{1}{2}\left( {{{\left| {{T_{13}}} \right|}^2} - {{\left| {{T_{12}}} \right|}^2}} \right)$
${\left| {{T_{23}}\left( \theta \right)} \right|^2}$	$\frac{1}{4}{\left\{ {\frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]} \right\}^2}$	$\begin{array}{l}\frac{1}{2}\left\{ {\frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2} + {\rm R}{{\rm e}^2}\left[ {{T_{23}}} \right]} \right\}\\ + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{23}}} \right]\end{array}$	8	$\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right){\rm Re}\left[ {{T_{23}}} \right] \\+ j\frac{1}{2}\left[ {{{{\mathop{\rm Re}\nolimits} }^2}\left[ {{T_{23}}} \right] - \frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2}} \right]$

下载: 导出CSV 

| 显示表格

基于上述振荡参数集，文献[12]还导出了一系列的极化角参数集，如极化零角参数、极化最大化角参数以及极化最小化角参数等。其中，极化零角参数的定义为在绕雷达视线的旋转域中使极化相干矩阵某元素取值为零的旋转角，即：

其中， ${\theta _{{\rm{null}}}}$即极化零角参数。由于表1中相互独立的5个初始角度 ${\theta _0}$分别为 ${\theta _0}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$, ${\theta _0}\_{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$, ${\theta _0}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$, ${\theta _0}\_{\left| {{T_{12}}\left( \theta \right)} \right|^2}$和 ${\theta _0}\_{\left| {{T_{23}}\left( \theta \right)} \right|^2}$，故相应的极化零角参数有 ${\theta _{{\rm{null}}}}\_{\rm Re}\left[\! {{T_{12}}\left( \theta \right)} \!\right]$, ${\theta _{{\rm{null}}}}\_{\mathop{\rm Im}\nolimits} \left[\! {{T_{12}}\left( \theta \right)} \!\right]$, ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ \!{{T_{23}}\left( \theta \right)} \!\right]$, ${\theta _{{\rm{null}}}}\_{\left| {{T_{12}}\left( \theta \right)} \right|^2}$和 ${\theta _{{\rm{null}}}}\_{\left| {{T_{23}}\left( \theta \right)} \right|^2}$。由文献[12]可知，各初始角度与其相应极化零角参数所包含的极化信息是相互等价的，且极化零角参数具有相对明确的物理意义，故在本文的后续部分均以极化零角参数代替相应的初始角度。

3.   结合旋转域极化特征的分类方法

3.1   旋转域极化特征的优选

文献[12]使用极化零角参数 ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ {{T_{12}}\left( \theta \right)} \right]$和 ${\theta _{{\rm{null}}}}\_{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$的组合能够成功辨识7类不同农作物，初步证实了极化零角参数集对于不同地物类别具有较好的区分能力。在此基础上，本文挖掘利用旋转域极化特征所蕴含目标在旋转域中的隐含信息，并将其应用于极化SAR地物分类。

在此之前，需要基于地物分类的应用背景对众多的旋转域极化特征进行优选处理。在文献[12]所导出的一系列旋转域极化特征之中，以不同地物类别样本集相互之间的“类间距最大化”为准则，进行相应的旋转域极化特征优选。具体步骤为：首先对各旋转域极化特征参数进行归一化处理；然后将不同的地物类别两两组合形成若干的地物类别对；接着针对各地物类别对，以其中两地物类别之间的类间距为标准，优选出使其取值达到最大的旋转域极化特征，则每个地物类别对均对应于一个优选的旋转域极化特征；最后，将各地物类别对的优选结果进行“取并集”处理，进而得到最终的优选结果。

文献[12]所导出相互独立的旋转域极化特征共有12个，分别为 ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ {{T_{12}}\left( \theta \right)} \right]$, ${\theta _{{\rm{null}}}}\_{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$, ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ {{T_{23}}\left( \theta \right)} \right]$, ${\theta _{{\rm{null}}}}\_{\left| {{T_{12}}\left( \theta \right)} \right|^2}$, ${\theta _{{\rm{null}}}}\_{\left| {{T_{23}}\left( \theta \right)} \right|^2}$, $A\_{\rm Re}\left[ {{T_{12}}\left( \theta \right)} \right]$, $A\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$, $A\_{T_{12}}\left( \theta \right)$, A_ ${T_{23}}\left( \theta \right)$, $B\_{T_{12}}\left( \theta \right)$, $B\_{T_{33}}\left( \theta \right)$, $B\_{T_{23}}\left( \theta \right)$。针对之后实验部分所使用的AIRSAR数据(15类地物，两两组合形成105个地物类别对；其它说明见4.1节)以及多时相UAVSAR数据(7类地物，两两组合形成21个地物类别对；4个数据获取日期；其它说明见4.2节)，上述特征优选流程所得结果如表2所示。

表  2  针对不同极化SAR实测数据的特征优选结果

Table  2.  Selected features for different PolSAR data

实测数据		优选所得旋转域极化特征(相应地物类别对的个数)
AIRSAR		${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]\left( {18} \right)$, ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]\left( {15} \right)$, ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]\left( {71} \right)$, $B\_{T_{33}}\left( \theta \right)\left( {1} \right)$
UAVSAR	6月17日	${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$(5), ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(12), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(4)
	6月22日	${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$(5), ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(14), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(2)
	7月03日	${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(3), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(18)
	7月17日	${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$(7), ${\theta _{{\rm{null}}}}\_{\rm Im}\left[ {{T_{12}}\left( \theta \right)} \right]$(5), ${\theta _{{\rm{null}}}}\_{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$(9)

下载: 导出CSV 

| 显示表格

综合考虑表2中的优选结果，并在追求较高地物分类精度的同时，将两组实测数据优选得到的旋转域极化特征进行统一，故本文优选部分的最终结果为3个极化零角参数，即 ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ {{T_{12}}\left( \theta \right)} \right]$, ${\theta _{{\rm{null}}}}\_{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$和 ${\theta _{{\rm{null}}}}\_{\rm Re}\left[ {{T_{23}}\left( \theta \right)} \right]$。

3.2   本文分类方法

为了将目标在旋转域中的隐含信息充分利用在极化SAR地物分类中，同时又发挥传统的旋转不变极化特征参数H/A/ $\alpha$/SPAN在极化散射机理解译方面的优点，本文提出了一种结合旋转域极化特征的极化SAR地物分类方法，其流程图如图1所示，相应的具体操作如下：

图  1  本文方法具体流程图

Figure  1.  Flowchart of proposed method

下载: 全尺寸图片 幻灯片

(1) 在进行Cloude-Pottier分解之前，需要对极化SAR数据进行相干斑滤波处理。本文采用新近提出的一种基于矩阵相似性检验的SimiTest自适应相干斑滤波方法^[16]对极化SAR数据进行滤波预处理。

(2) 基于滤波后的极化相干矩阵，计算总散射能量SPAN。

(3) 同样地，基于滤波后的极化相干矩阵，进行Cloude-Pottier分解，得到极化特征量H/ $\alpha$/A。

(4) 同时，将滤波后的极化相干矩阵绕雷达视线旋转，计算上述优选部分所得的3个极化零角参数。

(5) 对上述7个极化特征参数分别进行归一化处理，以作为地物分类特征集输入至SVM分类器。

(6) 通过SVM相应的训练与测试过程，实现对不同地物类别的分类处理。

4.   对比实验及分析

为了验证新极化特征(即3个旋转域极化零角参数)的引入对于传统地物分类方法性能的提升作用，在对极化相干矩阵中全部极化信息进行利用的前提之下，将本文方法与仅使用旋转不变特征H/A/ $\alpha$/SPAN作为SVM分类器输入的传统方法进行对比。首先使用AIRSAR数据15类地物的分类验证本文方法的分类性能，再使用多时相UAVSAR数据7类地物的分类进一步验证本文方法对多时相数据的稳健性。在对此两组数据分别进行SimiTest相干斑滤波^[16]时，所用滑窗大小均为15×15。对SVM分类器，各类地物样本的一半用于训练，另一半用于测试。

4.1   对AIRSAR数据的地物分类对比实验

本文首先使用NASA/JPL AIRSAR系统在荷兰Flevoland地区所获取的L波段全极化SAR数据进行地物分类实验。该数据方位向分辨率为12.1 m，距离向分辨率为6.6 m，所用区域大小为736×1010。SimiTest相干斑滤波后的Pauli RGB图如图2(a)所示。该区域的真值图如图2(b)所示，其中主要包含茎豆、豌豆、森林、苜蓿、小麦1、甜菜、土豆、裸地、草地、油菜籽、大麦、小麦2、小麦3、水域以及建筑物等15类地物。

图  2  AIRSAR数据

Figure  2.  AIRSAR data

下载: 全尺寸图片 幻灯片

使用传统方法和本文方法分别对滤波后的数据进行分类处理，所得结果如图3所示。

图  3  AIRSAR数据的分类结果

Figure  3.  Classification results of AIRSAR data

下载: 全尺寸图片 幻灯片

两种方法对AIRSAR数据15类地物分类处理所得精度如表3所示。通过比较可知，本文方法得到的总体分类精度为92.3%，优于传统方法91.1%的分类精度。且本文方法对草地77.3%的分类精度相较于传统方法的59.3%提升了18个百分点。另外，由于SVM分类器所用分类策略以总体分类精度的最大化为目标，无法保证单一地物类别的分类精度均达到最优。例如，本文方法在苜蓿、小麦1、裸地、大麦以及建筑物等5种地物类别区域所得分类精度均不及传统方法。针对其中分类精度差距最大(约8.3%)的裸地，由于其相应区域的主要散射机制为“面散射”，不同方位取向对其后向散射的影响较小，使用传统的旋转不变极化特征已经能较好地对其进行区分与辨识，本文方法额外引入的3个旋转域极化零角参数可能造成了分类信息的冗余，进而导致所得分类精度的较大幅度下降。

表  3  两种方法所得AIRSAR数据15类地物及总体的分类精度(%)

Table  3.  Classification accuracy of different terrains in AIRSAR data using two methods (%)

地物	传统方法	本文方法
茎豆	97.2	98.0
豌豆	93.7	96.9
森林	92.6	93.7
苜蓿	96.8	96.6
小麦1	88.7	85.9
甜菜	93.8	93.8
土豆	92.6	93.3
裸地	95.5	87.2
草地	59.3	77.3
油菜籽	83.9	88.0
大麦	92.6	91.5
小麦2	89.2	89.4
小麦3	94.3	95.9
水域	98.0	98.5
建筑物	84.9	83.2
总体精度	91.1	92.3

下载: 导出CSV 

| 显示表格

4.2   对多时相UAVSAR数据的地物分类对比实验

本文使用NASA/JPL UAVSAR系统在加拿大Manitoba地区所获取的多时相L波段全极化SAR数据进行地物分类实验。该数据方位向分辨率为7 m，距离向分辨率为5 m，所用区域大小为1325×1011。多时相极化SAR数据分别获取于6月17日、6月22日、7月3日以及7月17日。SimiTest相干斑滤波处理之后多时相极化SAR数据对应的Pauli RGB图如图4所示。该区域的主要地物类型是以谷物和油种产品为代表的混合型牧场农作物。相应的真值图如图5所示，其中主要包含阔叶林、草料、大豆、玉米、小麦、油菜籽以及燕麦等7类地物。

图  4  多时相UAVSAR数据滤波后Pauli RGB图

Figure  4.  Filtered Pauli RGB images of multi-temporal UAVSAR data

下载: 全尺寸图片 幻灯片

图  5  所用区域的真值图

Figure  5.  Gound truth of the multi-temporal data

下载: 全尺寸图片 幻灯片

使用传统方法和本文方法分别对滤波后的多时相极化SAR数据进行相互独立的分类处理，所得结果分别如图6和图7所示。

图  6  传统方法对多时相UAVSAR数据分类结果

Figure  6.  Classification results of multi-temporal UAVSAR data using conventional method

下载: 全尺寸图片 幻灯片

图  7  本文方法对多时相UAVSAR数据分类结果

Figure  7.  Classification results of multi-temporal UAVSAR data using proposed method

下载: 全尺寸图片 幻灯片

如图6(c)和图7(c)所示，基于7月3日获取的数据，传统方法将红色圆框内小麦与燕麦的绝大部分错分为了大豆，而本文方法在该区域的分类性能相较于前者有显著提升。又如图6(d)和图7(d)所示，基于7月17日获取的数据，传统方法将白色圆框内小麦的绝大部分错分为了大豆，而本文方法在该区域的分类精度相较于前者也有较大提升。

两种方法对多时相UAVSAR数据7类地物分类处理所得精度如表4所示。通过比较可知，对不同日期获取的数据，本文方法所得各类地物及总体的分类精度均优于或相当于传统方法。其中，对6月17日、6月22日、7月3日以及7月17日4个不同日期所获取的数据，本文方法得到的总体分类精度分别为94.98%, 95.12%, 95.99%以及96.78%，而传统方法所得总体分类精度则波动于80.87%至90.75%之间，出现约10%的起伏。具体就小麦和燕麦而言，本文方法得到的分类精度均分别保持在94%和92%以上，而传统方法所得相应分类精度则分别出现了约30%和23%的波动起伏。另外，本文方法95.72%的平均总体分类精度相较于传统方法的87.80%提升了约8个百分点。故本文方法较好的分类性能对于同一系统的多时相数据更具稳健性。

表  4  两种方法所得多时相UAVSAR数据7类地物及总体的分类精度 (%)

Table  4.  The classification accuracy of different terrains in multi-temporal UAVSAR data using two methods (%)

日期	方法	阔叶林	草料	大豆	玉米	小麦	油菜籽	燕麦	总体
6月17日	传统	98.47	62.24	92.64	96.12	93.63	91.70	86.37	90.19
	本文	98.49	81.65	96.76	98.19	96.08	92.25	96.32	94.98
6月22日	传统	98.05	61.38	94.14	97.30	97.89	93.82	77.29	90.75
	本文	97.96	72.60	96.86	98.18	97.07	96.84	95.13	95.12
7月3日	传统	97.41	54.38	90.45	98.89	68.75	98.81	63.46	80.87
	本文	97.77	76.68	98.12	99.08	96.95	98.93	94.22	95.99
7月17日	传统	96.86	64.51	97.38	99.78	84.76	92.19	82.98	89.39
	本文	97.27	93.15	99.31	99.58	94.73	99.71	92.16	96.78
平均	传统	97.70	60.63	93.65	98.02	86.26	94.13	77.53	87.80
	本文	97.87	81.02	97.76	98.76	96.21	96.93	94.46	95.72

下载: 导出CSV 

| 显示表格

另外，对于6月22日所获取数据中的阔叶林和小麦，以及7月17日所获取数据中的玉米，本文方法所得分类精度均略低于传统方法，且分类精度的差距均在1%以内。

在上述两组相互独立的对比实验所得结果中，本文方法所得分类精度均优于传统方法。故本文方法所表现出的较好分类性能对于不同系统的数据也具有较强稳健性。

5.   结论

目标方位取向对其后向散射响应的直接影响极易引起散射机理的解译模糊，进而限制仅使用旋转不变特征参数作为分类特征集的极化SAR地物分类所得精度。针对这一问题，本文将刻画目标旋转域隐含信息的旋转域极化特征用于极化SAR地物分类，并提出了一种结合旋转域极化特征和旋转不变特征H/A/ $\alpha$/SPAN的极化SAR地物分类方法，该方法将旋转域极化零角参数和H/A/ $\alpha$/SPAN联合作为分类特征集输入至SVM分类器。

将本文方法与仅使用旋转不变特征H/A/ $\alpha$/SPAN作为SVM分类器输入的传统方法进行比较：对AIRSAR数据15类地物分类而言，本文方法总体分类精度达到92.3%，优于传统方法的91.1%。对多时相UAVSAR数据7类地物分类而言，本文方法平均总体分类精度达到95.72%，显著优于传统方法的87.80%，表明本文方法对同一系统的多时相数据更具稳健性。这两组对比实验也表明本文方法较好的分类性能对于不同系统的数据具有较强稳健性。

通过对旋转域中目标极化散射信息的深入挖掘，能够为极化SAR图像的解译与应用提供一条新的可行途径。下一步将考虑旋转域极化特征与具有深度学习能力的卷积神经网络等分类器相结合，以实现更高的分类精度。另外，对极化特征参数更优的选择准则及相互融合也是我们未来将要深入研究讨论的内容。