-
摘要: 极化合成孔径雷达(Polarimetric synthetic aperture radar,PolSAR)是一种多参数、多通道的微波成像系统,在农林业、地质、海洋和军事等领域有着广泛的应用前景.PolSAR图像的相干斑抑制和分类是数据解译的重要环节,已经成为遥感领域的研究热点.本文综述了现有PolSAR图像的相干斑噪声抑制和分类方法并进行展望.首先,简要介绍了PolSAR系统的主要进展和应用;然后,对PolSAR图像相干斑抑制的评价指标和方法进行综述并对几种代表性方法进行了实验对比;接下来,对PolSAR图像的特征进行分析归纳,分别对有监督、无监督和半监督的PolSAR分类方法进行总结并给出了几种有监督分类方法的实验比较;最后,对PolSAR图像相干斑抑制和分类方法未来可能的研究方向进行了思考和讨论.Abstract: Polarimetric synthetic aperture radar (PolSAR) is a multi-parameter multi-channel microwave imaging system and it has wide applications in many fields such as agroforestry, geology, oceanography and military. Speckle reduction and classification are important steps in the PolSAR image interpretation process and the related techniques have become a hot research topic in remote sensing. This paper reviews the existing PolSAR image speckle reduction and classification methods. Firstly, the main progresses and applications of PolSAR systems are briefly described. Then, the speckle reduction methods for PolSAR image are summarized and comparative experiments are conducted on several representative despeckling methods. Next, the features for PloSAR classification are discussed and the algorithms for supervised, unsupervised and semi-supervised classification are reviewed. Besides, the experiments on several supervised classification methods are given. Finally, the potential further research directions of the PolSAR image speckle reduction and classification methods are discussed.
-
合成孔径雷达(Synthetic aperture radar, SAR)是一种先进的主动式微波成像传感器, 具有较强的穿透性, 能全天时、全天候工作. SAR通过距离向的脉冲压缩技术和方位向的合成孔径原理, 可获取大面积的高分辨率地面图像[1].极化合成孔径雷达(Polarimetric synthetic aperture radar, PolSAR)是一种多参数、多通道的新型SAR, 通过收发电磁信号的不同极化方式实现对场景目标的极化散射特性的探测.电磁波的极化特性对目标的散射机理高度敏感, 而散射机理反映了散射体的材质、湿度、几何形状和粗糙度等性质.因此, 与传统的单极化SAR相比, PolSAR能够获取更丰富的地物信息和分类特征[2]. PolSAR已成为国内外雷达技术研究的热点, 是微波遥感成像发展的主要方向之一.
PolSAR的发展始于20世纪80年代, 目前很多国家已经发展了多种以航天飞机或卫星为载体的PolSAR系统.代表性的机载PolSAR系统有美国宇航局喷气推进实验室(NASA/JPL)研制的AIRSAR, 德国空间中心(DLR)的E-SAR和F-SAR, 加拿大遥感中心(CCRS)的CONVAIR-580 C/X-SAR以及法国的RAMSES和SETHI, 丹麦的EMISAR和日本的PISAR等; 代表性的星载PolSAR系统有美国的SIR-C/X SAR和Mini-RF, 欧洲的ENVISAT-ASAR, 德国的TerraSAR-X, 加拿大的RADARSAT-2, 意大利的Cosmo-SkyMed, 日本的ALOS2-PALSAR, 以色列的TecSAR和印度的Chandraayan1-MiniSAR等[2-3]. 2016年8月, 我国首颗1米分辨率C频段多极化合成孔径雷达卫星高分三号成功发射, 并于2017年1月正式投入使用, 实现了我国低地球轨道长寿命、高分辨率SAR卫星研制技术的里程碑式跨越.
PolSAR在农林业、地质、海洋和军事等领域具有广泛的应用和研究价值.在农业方面, PolSAR可用作农作物的种类识别、农作物生长状态监控和土地使用评估等.在林业方面, 利用PolSAR数据可对森林状态进行监控, 可及时发现火灾及过度伐木等问题, 另外可估计森林中的生物量.在军事方面, PolSAR可对桥梁、道路、飞机、车辆、油库等战略目标进行检测和识别.在地质方面, PolSAR可获取地质结构、矿物分布、地表粗糙度、地面覆盖程度、土壤湿度等信息, 还可以对地质灾害进行检测和评估.在海洋方面, PolSAR可用于冰川、海洋和海岸线的监控, 对海面溢油和舰船进行检测.在地形测绘方面, PolSAR可用来测定山脉、河流、城市等地面目标的形状和位置, 研究城市规划和变迁、湖泊分布及道路运输状况等[2].
本文综述了现有的PolSAR图像去噪和分类方法的研究进展, 并指出存在的问题和未来研究方向.第1节介绍PolSAR数据的形式.第2节对PolSAR数据的相干斑噪声抑制方法进行综述.第3节分析总结PolSAR图像的特征.第4节对有监督、无监督和半监督的PolSAR图像分类方法进行综述.第5节指出现有PolSAR数据去噪和分类方法存在的问题并探讨将来的研究方向.
1. PolSAR数据
PolSAR对目标的测量是通过收发电磁信号的不同极化方式来实现的, 获取的数据具有多通道、信息量大的性质.通常的极化散射体的数据表达方式有极化散射矩阵、极化协方差矩阵、极化相干矩阵、Stokes矩阵、Mueller矩阵及Kennaugh矩阵[2], 在此主要介绍前三种数据形式.
1.1 极化散射矩阵
PolSAR发射和接收的是极化电磁波, 可分解为水平极化波和垂直极化波.因此, 每个分辨率单元都可以用一个$ 2\times2 $的复散射矩阵$ S $表示, $ S $在线性极化基下可表达为
$ \begin{eqnarray*} {S = }\left[\begin{array}{cc} s_{hh} & s_{hv} \\ s_{vh} & s_{vv} \end{array}\right] \end{eqnarray*} $
其中, $ h $和$ v $分别表示水平和垂直极化, $ s_{hv} $是垂直发射水平接收的极化波的散射系数, 其他元素有相似定义. $ s_{hh} $和$ s_{vv} $被称为同极化分量, 对应入射场和散射场有相同的极化; $ s_{hv} $和$ s_{vh} $被称为交叉极化分量, 相应的入射场和散射场的极化是正交的.
1.2 极化协方差矩阵与极化相干矩阵
在介质互易、收发同置的情况下, 散射矩阵$ S $是对称的, 即$ s_{hv} = s_{vh} $.在Lexicographic矩阵基和Pauli spin矩阵基[2]的集合下, $ S $可分别分解为如下的散射向量: $ \boldsymbol{k}_L = [s_{hh}\; \sqrt{2}s_{hv}\; s_{vv}]^{\rm T} $和$ \boldsymbol{k}_P = $ $ [s_{hh}+s_{vv}\ \ s_{hh}-s_{vv}\ \ 2s_{hv}]^{\rm T}/\sqrt{2} $, 其中, $ \sqrt{2} $可使散射向量和散射矩阵的总功率保持一致.极化协方差矩阵和极化相干矩阵分别定义为
$ \begin{eqnarray*} C = {\rm E}(\boldsymbol{k}_{L}\boldsymbol{k}_{L}^{\rm H}), \quad T = {\rm E}(\boldsymbol{k}_{P}\boldsymbol{k}_{P}^{\rm H}) \end{eqnarray*} $
其中, $ {\rm E}(\cdot) $表示数学期望, 是假设随机散射介质各向同性下的空间统计平均, 上标$ {\rm H} $表示向量复共轭转置.矩阵$ C $和$ T $都是复Hermitan矩阵, 并且两者可以相互转化, 即$ T = VCV^{\rm T} $, 其中
$ \begin{eqnarray*} V = \frac{1}{\sqrt{2}} \left[\begin{array}{rrr} 1 &0 &1\\ 1 &0 &-1\\ 0 &\sqrt{2} &0 \end{array}\right] \end{eqnarray*} $
2. PolSAR数据相干斑抑制方法
2.1 相干斑噪声的产生机理
合成孔径雷达是一种相干成像系统[1].对于分布式目标, 可看做每个分辨率单元中有许多理想点目标, SAR接收的信号是这些理想点目标的回波的矢量和.因此交叠区会出现某些地方震动加强、某些地方震动减弱或完全抵消的现象, 这种现象称为干涉, 产生干涉现象的电磁波称为干涉波.由于PolSAR发射的是相干电磁波, 当照射到一个粗糙表面时, 返回的信号包含了一个分辨单元内部许多基本散射体的回波, 由于表面粗糙的原因, 各基本散射体与传感器之间的距离是不一样的.因此, 接收到的回波尽管在频率上是相干的, 但在相位上不再是相干的.如果回波相位一致, 那么接收到的是强信号, 如果回波相位不一致, 则接收到的是弱信号.因此, 散射回波系数会有很大的随机起伏, 进而导致图像上出现很多类似颗粒状的斑点, 这就是相干斑噪声[4].
相干斑噪声是所有基于相干原理的成像系统如合成孔径雷达、声纳、激光等固有的缺陷, 是系统噪声, 无法避免[1].相干斑噪声的存在降低了PolSAR图像的质量, 使图像变得很难解译, 严重影响了图像分类、目标检测和识别的精确性和有效性[2].因此, 相干斑抑制是PolSAR数据解译前的重要预处理.
由于PolSAR数据的协方差矩阵$ C $或相干矩阵$ T $中包含了极化通道间的统计相关性和相位差等信息, 因此PolSAR图像相干斑噪声抑制应该在整个$ C $或$ T $数据上执行, 而非使用散射矩阵$ S $数据[5].由于协方差矩阵的对角线元素是实数并且是乘性噪声, 而非对角线元素是复数并且是加性和乘性噪声的混合[6], 因此, PolSAR数据的去噪问题比光学图像或单通道SAR图像的去噪更难处理.
2.2 评价指标
Lee等[7]提出PolSAR数据相干斑抑制需要遵循以下原则: 1)协方差矩阵的每个元素应以相同的方式进行滤波处理以保持极化性质; 2)对各元素的滤波应相互独立以避免引入通道间的串扰; 3)尽量保持图像的分辨率及避免边缘模糊. Lopez-Martinez等[8]指出只要不引起信息损失并且能改善滤波图像质量, 可以松弛前两条原则.总的来说, 一个好的PolSAR去相干斑方法应该具备如下性质:在均匀区域相干斑得到充分的抑制; 很好保持图像的空间分辨率和极化性质并且不会引起图像的模糊和降低图像质量[2].下面介绍评估相干斑抑制方法性能的量化指标.
2.2.1 极化信息保持评估
对于模拟图像, 极化信息的保持程度可以通过下述指标进行评估, 首先计算相对误差$ |\eta-\hat{\eta}|/|\eta| $, 其中$ \eta $是真实值, 然后使用中值操作得到一个综合的性能指标[9], 这些值越小, 极化信息保持得越好.
1) 辐射度参数:使用Span数据$ |s_{hh}|^2+2|s_{hv}|^2 $ $ + $ $ |s_{vv}|^2 $计算.
2) 复相关参数:首先通过协方差矩阵中元素的不要超界比值, $ C_{12}/\sqrt{C_{11}C_{22}} $, $ C_{13}/\sqrt{C_{11}C_{33}} $, $ C_{23}/ $ $ \sqrt{C_{22}C_{33}} $计算相应的幅度和相位, 然后分别计算它们与对应的真实值的相对误差并取中值, 其度量了三个极化通道之间的相关性.
3) Cloude-Pottier参数:包括熵、各向异性和平均$ \alpha $角, 这些指标的计算方式可参考文献[2], 它们与相干矩阵的特征值和特征向量相关.
4) 同极化和交叉极化性质:计算方式为[10]
$ \begin{align} P = &\ [1, \cos(2\phi_i)\cos(2\psi_i), \\ &\ \cos(2\phi_i)\sin(2\psi_i), \sin(2\phi_i)]\; \times \\ &\ K [1, \cos(2\phi_j)\cos(2\psi_j), \end{align} $
$ \begin{align} &\ \cos(2\phi_j)\sin(2\psi_j), \sin(2\phi_j)]^{\rm T} \end{align} $
其中, $ \phi_i $和$ \psi_i $是椭圆率角和椭圆方向角, 描述了接收天线的极化性, $ \phi_j $和$ \psi_j $表示发射情况下相应的量, $ K $是Kennaugh矩阵, 可以从协方差矩阵中导出.同极化信号对应于$ \phi_i = \phi_j $和$ \psi_i = \psi_j $, 而交叉极化信号对应于$ \phi_i = -\phi_j $和$ \psi_i = \psi_j+\pi/2 $.先计算每一类的平均协方差矩阵, 在$ \phi_i\in[-\pi/4, \pi/4] $和$ \psi_i\in [-\pi/2, \pi/2] $的范围内计算这些参数, 然后求与真实数据的相对偏差, 最后取中值[9].
2.2.2 相干斑抑制评估
等效视数(Equivalent number of looks, ENL)度量了去噪图像在某一均匀区域的相干斑抑制程度, 其定义为[11]
$ \begin{eqnarray*} {\rm ENL} = \left(\frac{\mu}{\sigma}\right)^2 \end{eqnarray*} $
其中, $ \mu $为数据的均值, $ \sigma $为标准差, 都是在选出的均匀区域上计算的. ENL越大, 表明相干斑抑制得越充分.
2.2.3 边缘保持评估
基于比值平均的边缘保持程度(Edge-preservation degree based on the ratio of average, EPD-ROA)指标对乘性噪声模型度量边缘的保持程度具有很好的鲁棒性, 其定义为[12]
$ \begin{eqnarray*} {\text {EPD-ROA}} = \dfrac{\sum\limits_{i\in I}\left|\dfrac{E_{D1}(i)}{E_{D2}(i)}\right|}{\sum\limits_{i\in I}\left|\dfrac{E_{O1}(i)}{ E_{O2}(i)}\right|} \end{eqnarray*} $
其中, $ I $是一个灰度图像的指标集, $ E_{D1}(i) $和$ E_{D2}(i) $分别表示去噪后的图像沿着水平方向和垂直方向的相邻像素值.对于模拟数据与真实数据, $ E_{O1}(i) $和$ E_{O2}(i) $分别表示原始无噪图像与带噪图像沿水平和垂直方向的相邻像素. EPD-ROA的值越接近1, 说明边缘信息保持得越好.
近20年来, 研究人员针对PolSAR数据的相干斑抑制问题提出了多种方法, 包括空域滤波、非局部均值和变分方法等.其中有很多方法是从单极化SAR的噪声抑制方法发展而来, 甚至有不少方法都是将PolSAR数据噪声模型简化为多通道的单极化SAR噪声模型来处理, 或将其简单推广到PolSAR数据, 这在后续的各类方法的总结中会具体介绍.关于PolSAR数据相干斑抑制, 10年前的相关方法可参考文献[2, 13], 下面着重介绍近10年来提出的方法.
2.3 空域滤波方法
空域滤波是使用空域模板进行的一种邻域处理方法, 其机理是在待处理的图像中逐点地移动模板, 滤波器在该点的响应是通过事先定义的滤波器系数与滤波模板扫过区域的相应像素值的关系来计算[14].具体过程为寻找以某一像元为中心的邻域, 使用空域模板对该邻域进行事先设定的运算, 将滤波后的像素值作为中心像元的输出值.空域滤波可以分为线性滤波和非线性滤波, 线性滤波器将空域模板对应邻域内的像素值加权之和作为中心像素的相应输出, 包括均值滤波器、加权均值滤波器等; 非线性滤波器的滤波结果直接取决于考虑的邻域像素点的值, 包括中值滤波、最大最小值滤波器等.空域滤波方法往往执行简单、速度快, 是PolSAR相干斑抑制中最常用的方法.
Boxcar滤波是比较早期的空域滤波方法, 其对每个像素的邻域窗口内的协方差数据取平均作为中心像元的数据, 能有效去除各向同性区域内的噪声, 但会显著降低图像的空间分辨率并模糊图像中的边缘和目标点[2]. Lee等[7]提出的修正Lee滤波使用了边缘对齐的非方形窗口和最小均方误差滤波, 通过一组边缘检测模板找出同质区域, 并在该区域内进行Lee滤波.这两种方法都是从单极化SAR的斑点抑制方法发展而来, 修正Lee滤波能较好地保持边缘的纹理特征, 然而其使用span值来确定同质区域, 而span值未包含充分的像元散射特性信息, 因此该方法没有保持各像元主导的散射机制.为此, Lee等[15]提出了一种基于散射模型(Scattering model based filter, SMBF)的滤波方法, 基于极化协方差矩阵的Freeman-Durden分解得到的面散射、体散射和偶次散射的分量, 对像元进行无监督聚类, 然后进行类内滤波, 能在有效抑制相干斑的同时较好保持像元主散射机制. Vasile等[16]提出强度驱动的自适应邻域(Intensity-driven adaptive-neighborhood, IDAN)滤波, 使用区域增长技术, 相邻像素通过它们强度的相似性不断聚合, 从而对每个待滤波的像素建立一个自适应的邻域, 然后邻域内做简单平均或使用局部线性最小均方误差技术, 该方法比使用预定义窗口的方法更灵活, 并且更好地保持了空间分辨率, 但会产生一定的偏度[17]. Lopez-Martinez等[8]提出了一个基于模型的PolSAR (Model-based PolSAR, MBPolSAR)滤波, 基于协方差矩阵的乘性和加性混合噪声模型[6], 利用不同极化通道的相关结构, 将相同的协方差矩阵的对角线和非对角线元素分别处理, 能较好保持极化信息. Ding等[18]针对SMBF未考虑结构特征的问题, 提出了一种基于结构判别和混合四元素极化分解的斑点抑制方法, 分解得到的四元素除了面、体和偶次散射外还包含了它们的混合分量, 能有效地同时保持数据的极化性质和结构特征.前述的滤波方法往往受限于局部的邻域窗口, 双边滤波将窗口扩张到更大的尺寸. D$ ' $Hondt等[19]将传统的图像去噪中的双边滤波器拓展到PolSAR相干斑抑制中, 给出一种迭代的双边滤波方法, 充分利用数据的辐射信息和结构信息并基于K-L散度和测地线距离计算像素的相似性. Lang等[20]提出一种基于线条-边缘探测器和极化同质度量的自适应窗口的斑点滤波方法, 对于异质区域选择小的和非方形窗口以保持细节信息, 对同质区域选择大的方形窗口以尽可能抑制斑点噪声. Lee等[21]提出一种拓展的sigma滤波, 将单极化SAR斑点抑制的改进sigma滤波拓展到全PolSAR数据去噪中, 根据相干矩阵的对角线元素对数据进行分类, 对强目标不做滤波处理, 仅对分布式目标区域选取同质窗口进行sigma滤波, 能很好保持数据的散射特性.
上述方法都仅使用了与场景相关的空间信息, 而对于PolSAR时间序列数据, 时间信息的使用也非常重要. Alonso-González等[22]提出了一种基于二叉树(Binary partition trees, BPT)的PolSAR时间序列数据处理方法.首先根据目标性质构建BPT, 通过树修剪得到一组时空均匀区域, 在区域内进行数据处理, 能同时保持空间和时间域的分辨率而不会引入偏差和失真, 而BPT具有一个多层结构, 会带来较高的计算复杂度. Salehi等[23]提出了一种基于协方差矩阵所有元素的时间统计准则的多时间多维(Multitemporal multidimensional, MTMD)滤波方法, 能自适应地使用空间和时间域中具有相同散射特性的像素, 并能很好保持数据的极化信息.
2.4 非局部均值方法
非局部均值(Non-local means, NLM)方法是Buades等[24]于2005年最先提出, 基本思路是在以待估像素为中心的一个很大的搜索窗口中, 找出与该像素最相似的一些像素点, 相似性通过以这两个像素为中心的像素块进行衡量, 最后通过加权平均得到待估像素的值, 相似性越高的像素对应的权重也就越大.
针对PolSAR协方差数据斑点抑制问题, 相关的非局部均值方法有Chen等[25]提出的Pretest滤波, Zhong等[26]提出的NL-Lee滤波, Torres等[27]提出的随机散度滤波和Deledalle等[17]提出的NL-SAR方法.这些方法的最主要区别在于相似性像元块的度量方式不同. Pretest滤波[25]根据协方差矩阵的Wishart分布的似然比测试得到两个协方差矩阵$ X $和$ Y $的距离度量为
$ \begin{align*} d = &\ n(2q\ln (2)+\ln |X|\; + \end{align*} $
$ \begin{align*} &\ \ln |Y|-2\ln |X+Y|) \end{align*} $
其中, $ n $是等效视数, $ q $是矩阵的行数或列数, 通常$ q $ $ = $ $ 3 $, $ |\cdot| $表示矩阵的行列式.进一步, 可得到以$ X $和$ Y $为中心的两个像元块的距离度量为
$ \begin{align} D = \, &\sum\limits_{i = 1}^k d_i = n[2qk\ln (2)+\sum\limits_{i = 1}^k (\ln |X_i|\, +\\ &\ln |Y_i|-2\ln |X_i+Y_i|)] \end{align} $
其中, $ k $为像元块中像元的数目.两个像元间的权重设为$ W = \exp (-D/D_t) $, 其中$ D_t $为设定的阈值. Zhong等[26]将这种结构相似度量方法与Lee滤波的同质相似性度量相结合, 提出NL-Lee滤波.这种组合度量方法能改善非局部平均参数如像素块和搜索窗的尺寸的鲁棒性, 同时在相干斑抑制和细节保持之间能达到更好的平衡.随机散度滤波[27]基于统计中的假设检验以及随机散度判断两个样本是否来自同一分布, 从而计算出两个像元块的相似度, 然后通过P-value做软阈值操作得到权重, 该方法能有效增强极化熵并保持目标散射信息. NL-SAR滤波[17]给出了一个针对SAR、PolSAR和干涉SAR数据去相干斑的统一的非局部均值框架.其基于Wishart分布的广义似然比检验得到对应的两个协方差矩阵$ X $和$ Y $的距离度量为
$ \begin{eqnarray*} d = \frac{|X|^L|Y|^L}{|\frac{1}{2}(X+Y)|^{2L}} \end{eqnarray*} $
以$ X $和$ Y $为中心的两个像元块的距离为$ D = \sum_{i = 1}^k -\ln(d_i) $, 权重设为$ W = \exp (-D/h) $, 其中$ h $为滤波参数. Wang等[28]用有限混合模型刻画观测场景的变化, 并提出一种渐进的混合回归方法用于更好地估计先验目标的统计特性, 然后将导出的混合先验估计融合到非局部平均方法中, 从而将小图像块之间的局部变化模式纳入到相似性评估中.这种将有限混合模型回归引入到非局部平均的方法在平衡斑点抑制和特征保持上较有潜力, 能很好地保持显著的散射机制和小尺度的目标特征.近期的工作还有基于空间-波段信息和协方差矩阵的多波段PolSAR图像滤波方法[29]和基于恒虚警率边缘检测的自适应滤波方法[30].
这类方法在相干斑噪声抑制的效果上显著优于空域滤波方法, 并且能很好保持图像的重复性结构信息.然而, 由于需要进行图像块的相似性度量, 计算复杂度非常高; 另外, 在平坦区域噪声严重的情况下, 会试图保持某种结构而导致鬼影出现; 此外, 在相似结构缺乏的区域会出现块状效应.
2.5 变分方法
全变分(Total variation, TV)的概念首先由Jordan[31]于1881年引入, 是用以证明变差有界的不连续周期函数的傅里叶级数的收敛理论. 1992年Rudin等[32]提出了全变分正则化方法并用于图像去噪问题, 相应的模型记为ROF:
$ \begin{eqnarray*} \hat{X} = \arg\mathop{\min}_{X} \left\{ \|Y-X\|^2+ \lambda \|X\|_{\rm TV}\right\} \end{eqnarray*} $
其中, $ \|X\|_{\rm TV} $表示全变分正则化, 包括各项同性和各项异性两种定义形式.各项同性的TV算子更容易求解且更加常用, 其具体定义为
$ \|\nabla X\|_{2, 1} = \sum\limits_{i = 1}^N \sqrt{(\nabla_h X)_i^2+(\nabla_v X)_i^2} $
其中, $ \nabla $为一阶有限差分算子, $ \nabla_h X $和$ \nabla_v X $分别表示$ X $沿着水平方向和垂直方向的一阶有限差分.由于受噪声污染图像的总变分比无噪图像的总变分明显的大, 因而约束总变分就可以有效抑制噪声, 同时能很好保持图像的边缘特征.
近年来, 变分方法由于在图像去噪中的优势得到广泛的关注和研究.在PolSAR数据相干斑抑制方面, 2011年Li等[33]提出将PolSAR相干矩阵的对角线元素做对数变换, 然后直接使用ROF模型[32]做相干斑抑制的方法. 2013年Liu等[34]针对多通道SAR图像去相干斑提出了两个自适应的矢量变分模型, 如下:
$ \begin{align*} &\hat{X} = \arg\min\limits_{X} {VTV}(X)+\frac{1}{2}\sum\limits_{i, j} \lambda_{i, j} (X_{i, j}-Y_{i, j})^2 \\ &\hat{X} = \arg\min\limits_{X} {VTV}(X)+ \lambda \|X-Y\|^{p(Y, \nabla Y)} \end{align*} $
其中, $ Y $是观测数据, $ X $是待求数据, $ {VTV}(X) $是自然矢量全变差.上述两种方法都是将PolSAR数据噪声模型简化为多通道的单极化SAR噪声模型来处理, 仅考虑协方差矩阵或相干矩阵对角线上的实数数据, 而忽略了非对角线上的复数数据, 从而导致一些重要的极化性质的丢失, 例如极化通道间的相关性、Cloude-Pottier参数、同极化和交叉极化特征等.
为了解决上述问题, 相干斑抑制应该在整个协方差或相干矩阵的数据上执行[5], 这涉及到加性和乘性噪声的混合建模及复数优化问题.用于光学图像去加性噪声和单通道SAR图像去乘性噪声的变分模型不能直接推广用于PolSAR数据去噪.据此, 2015年Nie等[35]根据PolSAR协方差数据的统计特性和贝叶斯最大后验准则, 建立了Wishart全变分(Wishart total variation, WisTV)模型, 具体为
$ \begin{align} \hat{C} = &\ \arg\mathop{\min}_{C_i\in {H}_{++}^3} \sum\limits_{i\in I} \big[ \mathrm{tr}(C_i^{-1} Z_i) \\ &\ +\ln(\det(C_i))\big] +\lambda \| C \|_{\rm TV} \end{align} $
(1) 其中, $ Z_i $是观测的协方差矩阵, $ C_i $是待求的协方差矩阵, $ \lambda $是正则化参数, $ {H}_{++}^3 $表示$ 3\times 3 $ Hermitian正定矩阵的集合, $ \|C\|_{\rm TV} $表示四阶张量数据$ C $的全变差正则化. Nie等[36]针对WisTV方法会导致纹理等细节信息模糊的问题, 提出了基于Wishart分布的非局部变分(WisNLTV)模型.该模型是非凸的, 而且变量是高阶张量, 文中利用共轭函数思想和向前向后算子分裂方法求解.该方法在充分抑制相干斑噪声的同时能很好保持数据的极化性质和图像中的重复性结构信息.之后, Nie等针对文献[35-36]中的算法在优化过程中未考虑变量的正定约束而是直接把解往正定矩阵集做投影以及模型中的参数对数据敏感的问题, 提出基于sigmoid函数的数据归一化方法和基于矩阵正交分解的优化算法[37].
2017年, Deledalle等[38]针对多通道SAR图像斑点噪声抑制问题, 基于协方差矩阵的Wishart分布和矩阵对数变换, 提出了多通道对数高斯去噪(Multi-channel logarithm with Gaussian denoising, MuLoG)方法, 对应模型为
$ \begin{align} \hat{x} = &\arg\min\limits_{x} \!\sum\limits_{k = 1}^n \! \mathrm{tr}(\Omega (x_k)\!+{\rm e}^{\Omega (y_k)-\Omega (x_k)})\! +\sum\limits_{i = 1}^{D^2}\! \mathcal{R}(x^i) \end{align} $
(2) 其中, $ \mathcal{R}(x^i) $是正则化项, $ \tilde{\Sigma} = \Omega (x_k) = \mathcal{K}(A\Phi x+b) $是协方差矩阵$ \Sigma $经过矩阵对数变换后的矩阵, $ \mathcal{K} $是从$ \mathbb{R}^{D^2} $维的实值向量到$ D\times D $的Hermitian矩阵的映射, 对于PolSAR数据$ D = 3 $, $ (A, b) $是一个白化放射变换以去除通道间的相关性, $ \Phi $是一个平衡通道间噪声方差的归一化操作.因此, 估计的协方差矩阵$ \hat{\Sigma}_k = {\rm e}^{\Omega (\hat{x}_k)} $.通过交替方向乘子法, 可以将现有的高斯去噪子嵌入到算法中, 因此一组数据可以利用不同的算法做去噪处理, 然后通过比较选出去噪效果比较好的.这种方法也为相干斑抑制算法提供了一个参考基准.
除了上述三大类方法外, PolSAR数据的相干斑抑制方法还有Farage等[39]提出的变换域滤波器, 其是在小波变换域通过熵阈值收缩技术做去噪处理; Foucher等[40]提出的基于矩阵迹的偏微分方程滤波, 其等价于与方向高斯滤波的局部卷积; Sun等[41]提出的各项异性扩散滤波等. Ma等[42]近期发表的综述论文对基于偏微分方程和基于机器学习的相干斑抑制方法单独做了总结, 并对几种代表性方法做了系统的实验对比.
表 1总结比较了上述三大类PolSAR数据相干斑抑制方法的优缺点.其中, 空域滤波方法运算速度较快, 但重构的图像质量一般; 非局部均值方法能很好去噪并保持图像重复性结构特征, 但运算速度很慢; 变分方法能有效抑制噪声并保持图像边缘, 但会导致弱目标等细节信息的丢失.本节梳理了PolSAR数据相干斑抑制领域的研究现状, 并对一些代表性的方法展开介绍.需要说明的是, 上述的归类方式并非是严格的, 某些方法可能属于多个类别, 例如我们虽然将NL-Lee滤波归到非局部均值方法中, 但它实际上是结合了空域滤波和非局部均值滤波的思想, 类似地, WisNLTV则是结合了非局部均值和变分方法的思想.
表 1 各类PolSAR数据相干斑抑制方法比较表Table 1 Comparison of different types of PolSAR data speckle reduction methods相干斑抑制方法 代表性算法 优点 不足 空域滤波方法 修正的Lee、IDAN、改进的sigma滤波 操作简单、计算速度快 空间分辨率低、边缘模糊、细节信息丢失 非局部均值方法 Pretest、随机散度滤波、NL-SAR 去噪效果显著、结构特征保持得好 计算速度很慢、出现像素块状效应 变分方法 WisTV、WisNLTV、MuLoG算法 去噪显著、速度快、边缘保持得好 弱目标丢失、平坦区域出现阶梯效应 2.6 实验对比
本节从各类方法中挑选了几个代表性算法, 包括修正的Lee滤波[7]、SMBF滤波[15]、NL-SAR方法[17]和WisNLTV方法[36], 进行对比实验.其中, 修正的Lee滤波和SMBF滤波分别使用$ 5\times5 $和$ 7 $ $ \times $ $ 7 $的滤波窗口, 这两个滤波器都是使用PolSARpro V5.1软件执行的; NL-SAR和WisNLTV均使用作者提供的代码, NL-SAR的搜索窗口和像元块的尺寸分别设置为$ 25\times25 $和$ 9\times 9 $; WisNLTV的像元块尺寸设置为$ 9\times 9 $, 挑选的最近邻像元个数设置为10.实验使用的数据是由德国空间中心的机载ESAR系统获取的Oberfaffenhofen地区的全极化、L波段、单视数据.
图 1展示了不同方法对应的相干斑抑制结果的Pauli分解图, 可以看出这4种方法都不同程度地抑制了斑点噪声.但修正的Lee和SMBF滤波降低了图像的空间分辨率并且模糊了图像中的一些点和线特征.另外, 修正的Lee滤波会导致像素化效应, 见图 1 (b)中的平坦区域; 而SMBF滤波对黑色斑点的抑制不够充分. NL-SAR和WisNLTV滤波在空间分辨率保持方面要优于前两种方法, 因为它们都使用了中心像元的邻域像元块进行相似性度量, 而像元块可以较好空间结构信息.但NL-SAR在噪声严重的数据上表现欠佳, 并且会导致图像块效应, 如图 1 (d)中的线和点特征出现不同程度的模糊, 这是因为在这些结构附近很难找到相似性的像元块. WisNLTV较好地保持了图像中的结构, 如图 1 (e)中的点散射体和线特征.表 2展示了这4种方法在相干斑抑制和边缘保持方面的量化对比结果, 可以看出修正的Lee滤波得到了更加光滑的结果, 但边缘保持得不够好; WisNLTV方法在相干斑抑制上比其他方法要稍微弱一些, 这需要加大正则化参数进行调节, 但得到了更好的边缘保持结果.
表 2 不同PolSAR数据相干斑抑制方法的量化比较表Table 2 Quantitative comparison of different PolSAR data speckle reduction methods指标/方法 修正的Lee滤波 SMBF滤波 NL-SAR方法 WisNLTV方法 ENL 13.8698 10.9144 11.5435 10.7267 EPD-ROA水平方向 0.6863 0.7356 0.7505 0.7533 EPD-ROA垂直方向 0.7472 0.7957 0.8046 0.8109 3. PolSAR数据特征提取
特征提取是PolSAR图像分类的关键步骤, 本节将对PolSAR数据的特征进行整理总结, 并进行分类.
3.1 极化特征
极化特征包括基于原始极化数据的参数和基于极化目标分解的参数两类.原始参数是直接从PolSAR数据的元素或其变换中提取的特征[43], 包括:协方差矩阵和相干矩阵的对角线元素$ T_{11} $, $ T_{22} $, $ T_{33} $, $ C_{11} $, $ C_{22} $, $ C_{33} $, 这些参数包含了同极化后向散射系数; 交叉极化后向散射系数的幅度和相位$ |C_{12}| $, $ |C_{13}| $, $ |C_{23}| $, $ \angle C_{12} $, $ \angle C_{13} $, $ \angle C_{23} $, 其中$ |\cdot| $和$ \angle $分别表示取幅度和相位; 极化通道间比值$ 10 \ln({C_{33}}/{C_{11}}) $, $ 10 \ln({C_{22}}/{2C_{11}}) $, $ 10 \ln({C_{22}}/{2C_{33}}) $; 去极化率$ {C_{22}}/{(2C_{11}+2C_{33})} $; Span值$ C_{11}+C_{22}+C_{33} $.
极化目标分解是将极化观测数据分解成不同分量, 以揭示散射体的物理机制, 大致可分为四类[2]: 1)基于Mueller或Kennaugh矩阵的二分量分解方法, 例如Huynen分解[44], 其假定任意目标的Mueller矩阵可表示为单目标和$ N $目标的矩阵之和, 根据Mueller矩阵得到9个目标参数, 每个参数都包含目标的真实物理信息; 还有Yang分解[45]等. 2)基于协方差或相干矩阵的模型分解, 例如Freeman-Durden分解[46], 可将目标分解为由一阶Bragg表面散射得到的面散射、由二面角反射器得到的偶次散射和由植被冠层定向偶极子得到的体散射三类; Yamaguchi分解[47]是在Freeman-Durden分解基础上引入螺旋体散射, 常见于城市区域; 还有Freeman两元素分解[48]和Kusano等[49]提出的基于广义散射模型的相干矩阵分解方法. 3)基于协方差或相干矩阵的特征值或特征向量的分析, 例如Cloude-Pottier分解[50], 其根据相干矩阵的特征矢量分析定义了熵、平均$ \alpha $角和各项异性系数, 分别反映了散射随机程度、散射性质和次级散射机制的相对重要性; 还有van Zyl分解[51]等. 4)基于散射矩阵的相干分解, 例如Pauli分解[50], 其系数是$ |S_{hh} $$ + $ $ S_{vv} |^2 $, $ |S_{hh}-S_{vv} |^2 $和$ |2S_{hv} |^2 $, 其中Pauli基是完全正交的, 简单且具有一定抗噪能力, 但只能区分面散射和偶次散射; Krogager分解[52], 其系数为$ |S_{RL} | $, $ \min(|S_{LL} |, |S_{RR} |) $, $ \mathrm{abs}(|S_{LL} |-|S_{RR} |) $, 其中$ S_{RR} $ $ = {\rm i}S_{hv}+(S_{hh}-S_{vv} )/2 $, $ S_{LL} = {\rm i}S_{hv}-(S_{hh}-S_{vv})/2 $, $ S_{RL} = {\rm i}(S_{hh}+S_{vv})/2 $, 可区分大导电球、二面角和螺旋体散射; Cameron分解[53]是基于雷达目标的互易性和对称性做相干分解, 提供了散射体结构的物理信息, 可有效检测自然目标和人工目标, 但会对噪声敏感; 此外, 还有Touzi分解[54]等.
3.2 统计特征
根据第1节的介绍, PolSAR的每个分辨率单元都对应一个复散射矩阵$ S $, 在介质互易、收发同置情况下是对称的, 并可以退化成散射向量$ \boldsymbol{k}_L $和$ \boldsymbol{k}_P $.一般情况下, $ \boldsymbol{k}_L $可以用均值为零的多元复高斯分布进行建模[55]
$ \begin{eqnarray*} p_{\boldsymbol{k}_L} (\boldsymbol{k}) = \frac{1}{\pi^3 |C|} \exp(-\boldsymbol{k}^{\rm H} C^{-1} \boldsymbol{k}) \end{eqnarray*} $
其中, 复协方差矩阵$ C = {\rm E}(\boldsymbol{k}_{L}\boldsymbol{k}_{L}^{\rm H}) $, $ |C| $表示极化协方差矩阵的行列式.多视的极化协方差矩阵$ Z $是几个独立的单视协方差矩阵的平均, $ L $-视观测协方差矩阵可表示为$ Z = \frac{1}{L}\sum_{i = 1}^L {\boldsymbol{k}_{L_i}\boldsymbol{k}_{L_i}^{\rm H}} $, 其中$ \boldsymbol{k}_{L_i} $为第$ i $个单视样本的散射向量.注意到$ Z $是Hermitian矩阵, 并对任意的非零向量$ \mathit{\boldsymbol{u}} $, 都有$ \mathit{\boldsymbol{u}}^{\rm H}Z\mathit{\boldsymbol{u}}\leq 0 $, 因此$ Z $是半正定的.均值为$ C $的复Wishart分布[56]可以很好地描述$ Z $的统计性质, 其密度函数为
$ \begin{eqnarray*} p_z (Z) = \frac{L^{3L} |C|^{L-3} \exp (-L \mathrm{tr}(C^{-1} Z))}{K(L, 3) |C|^L } \end{eqnarray*} $
其中, $ \mathrm{tr}(\cdot) $表示矩阵的迹, $ K(L, 3) = \pi^3 \Gamma(L)\Gamma(L-1)\Gamma(L-2) $, $ \Gamma (\cdot) $表示gamma函数.同样地, $ L $-视观测相干矩阵$ W = \frac{1}{L}\sum_{i = 1}^{L} \boldsymbol{k}_{P_i}\boldsymbol{k}_{P_i}^{\rm H} $, 且有$ W = VZV^{\rm T} $, 由于$ V $是一个$ 3\times 3 $的满秩矩阵, 则$ W $服从均值为$ VCV^{\rm T} = T $的复Wishart分布, 其中$ T $ $ = $ $ {\rm E}(\boldsymbol{k}_{P}\boldsymbol{k}_{P}^{\rm H}) $.
然而Wishart分布是基于高斯模型的, 只能用于描述均匀区域的统计特性, 为更好刻画异质区域如城市或森林场景的统计性质, Freitas等[57]提出多元多视极化$ \mathcal{G} $分布, Khan等[58]给出了相应的单视极化$ \mathcal{G} $分布, Doulgeris等[59]提出K-Wishart分布.另外, Bombrun等[60]基于Fisher分布的纹理模型推导出多视协方差矩阵的Kummer U分布, 能区分具有不同纹理分布的区域. Frery等[61]在复Wishart分布基础上分析了Shannon, Rényi和限制Tsallis熵的解析表达.
3.3 纹理特征
纹理特征可以描述图像或其中小块区域的空间颜色分布和光强分布, 即刻画一个特定类的空间变化性.由于它是基于区域的特征, 因此对噪声有一定的鲁棒性. PolSAR图像的纹理特征一般是基于Span图像计算, 常用的[62]包括: 1)基于统计的特征. a)局部二元模式(Local binary pattern, LBP)特征, 其具有旋转不变性和多尺度特性; b)边缘柱状图描述子(Edge histogram descriptor, EHD), 将图像边缘分为5个方向, 即水平、垂直、45$ ^\circ $对角和135$ ^\circ $对角方向及各项异性方向, 然后计算直方图; c)灰度共生矩阵(Gray-level co-occurrence matrix, GLCM), 通常提取4个最有判别性的特征即能量、相关性、对比度和同质性, 能很好描述小的纹理和重复性结构模式, 如水和森林区域; d)序数共生矩阵(Ordinal co-occurrence matrix, OCM)是GLCM和序数描述子的结合, OCM在刻画人造目标区域这样的粗糙纹理和不规则模式方面要优于GLCM; e)多层局部模式柱状图(Multilevel local pattern histogram, MLPH)[63]专门用于刻画SAR图像纹理特征, 描述了以各种强度捕捉局部和全局结构信息的移动窗口中出现的明亮、黑暗和均匀模式的尺寸分布. 2)基于多尺度变换的特征. a) Gabor滤波, 能很好描述对应于空间尺度、空间位置及方向选择性的局部结构信息; b)小波变换, 利用联合的时间-尺度函数分析非平稳信号, 并通过伸缩和平移等运算功能对信号进行多分辨率细化分析, 具有良好的时频局部化、尺度变换和方向特性[64]; c) Wedgelet框架, 能以简单和最优的方式表示图像块内的直线边缘, 获取图像在不同尺度上的空间结构, 同时考虑了局部空间的复杂度[65].
3.4 颜色特征
由于原始的PolSAR数据未提供一种自然的颜色图像表示, 因此颜色特征在PolSAR图像分类中经常被忽略.而出于可视化的需要, 根据散射矩阵可创建伪彩色图像, 例如Pauli颜色编码图像.据此, 可以提取: 1)色调-饱和度-明度(Hue-saturation-value, HSV)空间的颜色柱状图, HSV是一种比较直观的颜色模型, 与人的视觉系统类似; 2)主导颜色描述子(Dominant color descriptor, DCD), 描述图像中任意不规则区域中处于支配地位的颜色信息, 即能用少量主导颜色充分描述图像的颜色内容; 3)颜色结构描述子(Color structure descriptor, CSD), 不仅能反映图像的色彩信息, 也能描述图像的局部颜色结构特征. Uhlmann等[62]证实了从伪彩色图像中提取的这些颜色特征能有效提升PolSAR数据分类的精度. Xing等[66]同时利用Pauli分解和Yamaguichi分解的伪彩色图, 并转化到CIELAB颜色空间用于聚类. Liu等[67]基于PolSAR数据的四分量分解模型, 计算了常用的4种颜色空间: YUV、RGB、HSI和CIELAB, 通过引入颜色熵从这些空间中量化选择颜色特征, 然后联合纹理特征和散射功率熵进行分割和聚类.
3.5 空间结构特征
如果仅使用单像素点的特征, 分类算法会对噪声比较敏感, 结果不够鲁棒.通过对PolSAR数据的空间结构信息的利用, 可以有效提高分类的精度和鲁棒性.目前主要有四类空间信息提取方法[68]: 1)基于超像素[69-70]和基于区域[71-72]的方法; 2)基于纹理特征提取的方法, 例如小波变换[73], Gabor滤波[62]和楔波分析[65]; 3)基于正则化准则的方法, 例如MRF[74]和流形正则化[75]; 4)基于语义同质区域提取的方法, 例如层次语义模型[76-77].与基于像素的分类方法相比, 这些方法考虑了图像中的局部空间结构信息, 可以抑制相干斑噪声对分类结果的影响, 使分类图更加准确和清晰, 并减少了孤立区域[78].
3.6 基于深度学习提取的特征
上面介绍的特征都是采用传统的人工指导的提取方法, 由于依赖手工设计, 这些特征往往比较简单, 而且手工设计出有效的特征是一个相当漫长的过程.另外, 它们主要依靠设计者的先验知识, 很难利用大数据的优势.目前出现了基于深度学习的特征提取方法, 与传统的提取方法思想完全不同.深度学习可以从大数据中自动学习特征的多层表示, 其中可以包含成千上万的特征, 这些学习出的特征的表达能力比手工设计出的特征更强.另外, 深度学习可以针对新的应用从训练数据中很快学习得到新的有效的特征表示.基于深度学习的方法也被广泛应用于PolSAR数据的特征提取中, 例如基于卷积神经网络(Convolutional neural networks, CNN)[79-80]、基于受限玻尔兹曼机的深度置信网络[81]以及基于多层自编码[82]的特征提取等.下文将详细介绍相关的方法.
4. PolSAR图像分类方法
PolSAR图像分类是数据解译的关键环节, 也是PolSAR数据处理的重大研究方向之一.分类的目的是通过PolSAR传感器获得的测量数据, 确定每个像素所属的类别. PolSAR图像分类精度的常用评价指标包括各子类的分类正确率、总体分类正确率以及混淆矩阵.其中, 总体分类正确率是被正确分类的像素总和除以总的像素数目; 混淆矩阵是分类结果和地面真实信息的比较阵列, 列代表地面检验为该类别的数目, 行代表通过分类算法得到的类别数目, 通过混淆矩阵可以计算各种精度统计值, 例如各子类和总体分类正确率、错分率、漏分率和Kappa系数等.目前的PolSAR分类方法主要包括有监督分类、无监督分类和半监督分类. 10年前的相关方法可参考文献[2-3], 本文重点梳理和归纳近十几年的研究成果.
4.1 有监督分类方法
有监督分类方法可以使用一定量的有标签样本训练分类器, 进而对无标签样本进行分类.有监督PolSAR图像分类方法包括特征提取和分类器设计两个阶段.针对特征提取方面, 第3节主要从特征本身的角度归纳总结了6类常用的特征及相应提取方法, 接下来从分类的角度梳理基于不同特征类型或多特征融合的有监督分类方法, 主要包括: 1)基于PolSAR数据统计特性的方法, 如分别使用复Wishart分布[83]、非高斯分布[59]和多个统计分布[84]的分类方法, 以及基于相干矩阵特征值统计分布的方法[85]. 2)基于极化分解特征的方法, 如Alberga等[86]比较了将相干参数和非相关参数作为输入的最小距离分类器、多层感知机(Multi-layer perceptron, MLP)和基于logistic回归(LR)分类器的效果, 结果显示MLP和LR分别在相干参数和非相干参数上的分类精度更高, 另外, 还有使用极化方向角[87]和使用基于散射机制统计特征[88]的PolSAR图像城市区域提取方法, 基于最优的Touzi分解参数的方法[89]以及基于同极化、交叉极化等多种极化性质融合的方法[90]等. 3)基于纹理特征的方法, 如基于小波变换纹理特征提取的方法[91], 结合小波极化信息、纹理基元和稀疏编码的方法[92], 结合纹理和极化信息的方法[64, 93], 结合时频分析技术和球面不变随机向量乘积模型的方法[94]. 4)融合颜色特征的分类方法, 鉴于之前的方法都忽视了颜色特征的使用, Uhlmann等[62]从PolSAR伪彩色图像中提取了颜色特征, 并与纹理、极化分解特征相结合, 实验证实了颜色特征的加入可有效提高分类精度. 5)基于线性或非线性降维的方法, 可以从冗余信息中提取更有判别性的低维特征, 能有效提升分类效果, 如Tu等[43]将极化分解参数堆叠成一个高维向量, 使用非线性流形降维技术即拉普拉斯特征图, 获取该高维向量在低维空间的本质特征.然而该方法没有使用样本的判别信息, 并且拉普拉斯特征图没有显示表达, 因此不能拓展到新样本上, 为此, Shi等[95]使用线性降维技术即基于流形学习的有监督图嵌入算法, 获取能保持判别信息的低维特征.然而这两种方法都是基于单像素的, 未考虑像素邻域的空间信息, 为此, Tao等[96]结合空间关系并利用Tucker张量分解和独立主成分分析技术得到低维本质特征集; Huang等[97]利用空间和极化信息构建三阶张量并提出张量局部判别嵌入方法实现降维. 6)基于多视角学习的方法, 上述不同属性的特征可以看做不同的视角, 目前的分类方法多数是基于单视角数据的, 也有基于某几个视角数据的, 但往往是排成一个向量进行处理, 忽视了特征间的不同属性和视角间的冗余性和互补性等关系, 为此, Nie等[98-100]提出基于多视角的在线学习算法, 利用协同正则化思想, 实现了多视角特征的融合; Huang等[101]结合线性判别嵌入和典型性相关分析提出多视角特征提取算法, 有效提升了分类精度.
在分类器设计方面, 包括: 1)基于贝叶斯最大似然的分类器, 如基于散射向量复高斯分布的分类器[102]和基于协方差矩阵的复Wishart分布的分类器[103]. 2)基于神经网络的分类器, 神经网络技术最初由Pottier等[104]引入到PolSAR图像分类领域, 后来各种新的神经网络方法又相继被应用到该领域, 如动态学习神经网络[105], 模糊神经网络[106], RBF神经网络[107], 四元数神经网络[108]和概率神经网络[109]. 3)基于支持向量机(Support vector machine, SVM)的分类器, 最初由Fukuda等[110]引入到PolSAR应用中, 由于SVM可以利用各种类型的特征并具有良好的分类效果和泛化能力, 在该领域得到了非常广泛的应用[92, 111-112]. 4)基于字典学习和稀疏表示的方法, 如Zhang等[113]使用训练样本的特征向量构建过完备字典并得到相应稀疏表示系数, 根据残差最小的原子判断最终类别, 还有基于监督非相干字典学习的目标检测方法[114]和基于多尺度压缩感知金字塔的分类方法[115]. 5)其他分类器, 如基于Boosting[116]、随机森林[117]、决策树[118]和扩散-反应系统[119]的分类方法. 6)基于多分类器融合的方法, 如Maghsoudi等[120]提出根据不同特征子集的SVM分类器的集成方法; Ma等[121]提出的利用像素间空间关系的Wishart、SVM、K-均值分类器融合的方法; Masjedi等[64]将MRF的能量函数结合到SVM的优化框架中, 同时使用了纹理、极化和上下文特征, 并组合了SVM、Wishart和MRF分类器, 使分类精度得到显著提升.这种多分类器融合的方法也是一种发展趋势.
上述介绍的分类器各有优缺点, 如贝叶斯分类器对大数据量训练时具有较高的速度, 并且结果容易解释, 但需要样本的独立性假设, 样本相关性较强时会影响效果; 神经网络学习能力强, 分类精度高, 对噪声数据的鲁棒性强, 但需要调整的参数较多, 学习过程比较长, 不能观察中间结果; SVM可以解决小样本和非线性问题, 很好处理高维数据集且泛化能力比较强, 但训练过程比较低效且合适的核函数的选择也是一个挑战; 稀疏表示分类器对数据缺损不敏感, 对测试集中的噪声比较鲁棒, 能达到较高的分类精度, 但对小样本问题, 分类效果一般; boosting方法很好地利用弱分类器进行级联, 比较简单, 不容易发生过拟合且不用做特征筛选, 但弱分类器数目不好确定且数据不平衡时会导致分类精度下降; 决策树可以可视化分析且容易提取规则, 运行速度也比较快, 但对缺失数据处理困难且容易出现过拟合问题.上面总结的分类算法都是离线学习方法, 不能适用于动态大规模PolSAR图像的分类问题, 而在线增量学习算法可以处理这种动态数据流问题, 其使用新来样本快速更新分类器, 目前这方面的工作包括: Kiranyaz等[122]提出的二分类器集成网络, 能自适应学习新样本、新类、新特征并不需要重新训练所有已知数据, 但速度较慢; Nie等[98-100]提出的在线多视角学习算法, 只使用新来样本更新分类器, 为多线性分类器组合, 执行速度非常快.
4.2 无监督分类方法
无监督方法不需要利用有标签样本进行预训练, 主要通过数据的不同特征分离出属于同一类型的像素.目前的方法主要包括两大类:基于目标统计或电磁散射特性的分类方法和基于机器学习理论的方法.
第一大类可划分为三类:基于极化目标分解的算法、基于数据统计特性的算法和基于两者结合的算法.基于极化目标分解的算法主要是根据不同地物的物理散射特性对目标进行分解进而进行分类, 如van Zyl[123]提出的基于入射波和散射波间的相位与自旋关系的方法; Cloude等[124]提出的基于$ H $/$ \alpha $的特征值分解方法, 该方法利用极化分解得到的极化熵$ H $和平均散射角$ \alpha $描述散射机理, 根据$ H $/$ \alpha $平面将地物分为八类; 另外, 还有基于目标散射向量去取向理论的分类方法[125], 利用目标散射分解的旋转不变参数进行分类的方法[54], 基于极化散射相似性的分割方法[126]以及基于交叉散射的方法[127].基于数据统计特性的算法是根据不同的数据统计分布和贝叶斯最大似然准则设计的分割算法, 如基于K-Wishart分布[128]和Kummer U分布[129]的方法.另外, Frery等[130]提出基于极化$ \mathcal{G} $分布和B样条的边缘检测方法, Silva等[131]提出基于Wishart分布最小随机距离的分割算法, 并给出了对应于基于Wishart分布的Kullback-Leibler, Bhattacharyya, Hellinger, Rényi和Chi-Square随机距离的测试统计量的解析表达式.基于两者结合的算法是在极化分解得到的初始分割结果的基础上, 利用数据的统计特性进一步进行迭代分类, 如Lee等提出的$ H/\alpha $-Wishart算法[132], 将$ H/\alpha $平面的分类结果与Wishart最大似然分类器结合起来; 还有用于多频全极化数据的$ H/\alpha/A $-Wishart算法[133]和引入模糊概念的$ H/\alpha $分类方法[134].这类方法克服了基于极化分解的分类结果杂乱及类间界限模糊的问题, 但分类后的类别丧失了散射机制的可解释性.为有效保持各类目标的物理散射特性, Lee等[135]基于Freeman分解和Wishart分类器构建了新的分类算法, 首先利用Freeman分解将地物分成三大类, 然后再将每大类分成30小类, 对各子类进行聚类和Wishart迭代分类, 此方法能很好地保持目标的主要散射特性, 并且拥有较优的收敛性以及稳定性.
基于机器学习理论的方法包括聚类算法、基于马尔科夫随机场(Markov random fields, MRF)、字典学习[136]、均值漂移(Mean shift)[67]、流形正则化[75]和二叉树[22]的方法等.由于之前的多数方法[132, 134]需要预先知道类数, 而很多情况下, 预设的类数可能并不准确.为此, Cao等[137]提出了基于$ Span/H/\alpha/A $的凝聚层次聚类算法, 根据$ Span/H/\alpha/A $获取初始聚类中心, 然后根据Wishart距离度量进行层次聚类, 聚类数目可通过数据的对数似然算法自动估计.之后, Dabboor等[138]提出基于Wishart-Chernoff距离的层次聚类方法, 该算法在分类过程中可以控制类别的自动合并.聚类方法还包括Kersten等[139]提出的基于模糊C均值的聚类算法, Ersahin等[72]提出的谱聚类算法, 其借助空间邻近关系, 基于区域设置的轮廓信息进行分割, 然后通过图划分进行分类, 结果优于Wishart分类器.基于MRF的PolSAR图像分割方法最初由Rignot等[74]提出, 由于具有明确的物理意义及良好的效果而得到广泛关注, 很多方法也相继被提出, 如结合张量空间聚类分析和MRF框架的算法[140], 结合区域增长和MRF边缘增强模型的算法[141]以及基于自动$ {\mathcal U} $分布和MRF的算法[142].针对传统MRF方法不能自动确定精确的聚类数目的问题, Bi等[143]结合有监督的softmax回归模型和马尔科夫随机场光滑约束设计能量函数, 将类别和分类器同时作为变量进行优化, 达到了较高的聚类精度; Song等[144]提出将Dirichlet过程混合模型和基于邻域的相似性度量融入到MRF框架中, 与其他基于MRF的方法相比, 该算法能自动识别聚类数目并获得更准确的边缘位置和更光滑的同质区域.除了上述提到的能自动估计PolSAR类别数目的方法外, 还有Liu等[70]提出的基于超像素的分类框架用于具有自适应类数的PolSAR图像中, 其使用超像素之间的成对相异度信息估计数据的类的数量和每个类的中心, 相应的分类结果更加准确鲁棒; Zhong等[145]在Hermitian正定的协方差矩阵数据上引入基于稀疏度相似性度量的无监督分类方法, 其中最优的聚类数目可通过自调整的谱聚类或同时聚类和模型选择的方法来确定; Kim等[146]使用四元数自编码器提取极化特征, 进而通过无监督的四元数自组织神经网络进行细分组, 然后执行标签化处理使得分组对应于人类可理解的地物概念, 不需要给定人类预设的土地分类情况, 能发现新的和更细的土地类别; Liu等[147]给出一种具有双重约束的空间变化的Wishart混合模型用于PolSAR图像无监督分类, 其基于像素块水平以获取空间信息, 通过变分推断算法实现参数的闭式形式的更新以及聚类数目的自动确定.
4.3 半监督分类方法
在实际应用中, 有标记样本往往难以获取, 需要有经验的专家做标注或进行物理实验, 非常昂贵和耗时.相比较而言, 无标记样本往往非常充足且比较容易获取.半监督学习方法结合有监督和无监督学习, 同时使用标记样本和大量的未标记样本进行分类训练.文献[148]结合无监督的确定性模拟退火方法和有监督的多层感知机方法构建了半监督学习方法. Niu等[149]提出了一种半监督的随机期望最大化(Stochastic estimation maximization, SEM)算法, 通过应用自适应马尔科夫随机场(MRF)有效开发了时空上下文信息, 此外, 提出了学习控制方案确保鲁棒的半监督映射过程以避免不希望的类合并. Uhlmann等[150]针对小样本问题提出一种基于集成分类器的自训练方法, 在小规模数据集上训练的分类器能达到有监督学习算法在显著更大训练集上得到的分类器相似的性能. Wei等[151]提出一种基于超图学习的半监督分类方法, 能有效提升分类精度, 但创建超图的过程往往比较耗时.文献[78]给出一种基于图的半监督学习框架, 首先做超像素分割, 然后利用Wishart距离做聚类, 进而构建空间图结构进行半监督分类, 在小样本情况下提升了分类精度. Liu等[152]提出一种有邻域约束的半监督降维算法, 该方法提取的特征不仅具有判别性而且保持了PolSAR数据的结构, 同时对噪声比较鲁棒, 但这种邻域约束也会导致分类边界上的像素难以区分. Hua等[153]提出一种基于改进的协同训练的半监督算法, 首先分析两视角数据的充分性和独立性, 然后执行带有样本选择策略的协同训练过程, 最后通过基于相似性原则和超像素算法的后处理方法改善样本分类的一致性. Hou等[154]针对有限的手动标记的样本量会限制自动分类方法的性能和手动标记的训练样本受噪声等影响导致质量较低的问题, 提出一种鲁棒的半监督概率图分类框架, 在半监督学习时使用邻域像素的结构信息, 在训练分类器时使用鲁棒的分类损失函数, 另外, 通过混合的生成/判别分类框架解决有限标记数据不准确问题, 其中标记和未标记样本同时被用于学习低质量像素的高级特征, 从而保证了方法的有效性. Zhang等[155]针对经典的半监督SVM需要估计正负样本的分布问题, 提出一种基于不同干扰因素的学习基和聚类评估手段的集成半监督SVM模型, 可以处理未知分布或不平衡数据情况下的半监督分类问题. Liu等[156]给出一种基于空间多属性图模型的分类方法, 根据PolSAR数据特征组的物理属性构造多个图, 然后自动优化每个图的权重, 并结合像素的空间信息用于标签传播和分类, 可以避免特征堆叠和维数灾难. Geng等[157]提出一种具有多个决策的超像素约束深度神经网络用于PolSAR图像分类, 通过提取超像素空间特征可降低相干斑噪声的影响, 另外, 这种半监督方法可有效处理有限训练样本的问题.与有监督学习相比, 半监督学习能有效利用无标记样本的信息, 从而在少量标记样本情况下获得更好的分类效果.
上述方法都需要先提取人工设计的特征然后利用分类或聚类等算法进行有监督、无监督或半监督学习, 而基于深度学习的方法通过深度神经网络的多个隐含层和非线性变换函数, 可以从数据中自动学习平移不变的空间特征.目前随着深度学习的发展, 具体实现也逐渐产生了有监督[158]、无监督[159-160]和半监督[161]等不同学习策略.基于深度学习的PolSAR数据分类方法也由于较好的性能正得到越来越多的关注. Zhou等[79]将PolSAR协方差或相干矩阵转化为6维实值向量, 然后利用四层卷积神经网络对数据做特征提取和分类. Liu等[81]提出了一种基于Wishart-Bernoulli受限玻尔兹曼机的深度置信网络, 但该网络模型是非任务主导的并且训练过程会比较耗时.为此, Jiao等[162]设计了一种Wishart深层堆叠网络, 能快速执行分类任务. Hou等[82]提出了基于多层自编码和超像素的PolSAR图像分类方法. Zhang等[163]采用堆叠稀疏自编码器学习PolSAR数据的深度空间稀疏特征. Xie等[164]将PolSAR数据的Wishart距离度量用于自编码和卷积自编码神经网络的训练中以实现无监督的特征提取. Liu等[165]提出一种基于邻域保持的深度神经网络用于PolSAR特征提取和分类, 其使用了空间邻域信息和超像素信息对每个像素进行加权, 可保持局部结构并减少了对样本数量的要求.上述这些方法大都仅考虑了PolSAR数据的幅度信息, 而未使用数据的相位信息.为此, Zhang等[80]提出了一种复值卷积神经网络用于PolSAR数据特征提取, 将卷积神经网络的输入输出层、卷积层、激活函数和pooling层都拓展到了复数域, 能同时处理PolSAR数据的幅度和相位信息.这类基于深度学习的分类方法对数据有很强的拟合能力和表达能力, 可以达到非常好的分类效果, 但也存在诸多问题, 例如需要大量的样本进行训练, 隐藏层的信息难以解释以及容易陷入局部最小值等.
图 2归纳总结了PolSAR图像的分类中常用的特征类型, 以及根据是否有人工指导划分的有监督、无监督及半监督的分类方法.总的来说, 有监督学习方法需要有标注样本进行模型训练, 充分利用了数据的类别先验信息; 无监督学习主要通过数据的不同特征分离出属于同一类型的样本, 不需要数据的标签信息, 因此很难达到有监督分类的精度; 半监督学习结合了前两者, 同时使用有标注样本和大量无标注样本训练分类器, 在小样本集上的结果能达到有监督学习在更大训练集上相似的效果.本节梳理了PolSAR数据分类领域的研究现状和发展脉络, 并对一些代表性的方法展开介绍, 为读者提供参考.需要说明的是本节对相关研究方法的归类并非是严格的, 某些方法可能属于多个类别, 如基于深度学习的方法既可用于有监督分类, 也可用于无监督和半监督分类中.
4.4 实验对比
本节使用几种有监督的分类算法进行对比实验, 包括有监督的Wishart分类器(Supervised Wishart classifier, SWC)[83]、SVM分类器、基于局部保持投影(Locality preserving projection, LPP)的方法[43]、实值卷积神经网络(RV-CNN)[79]以及复值卷积神经网络(CV-CNN)[80].前三种算法使用的特征是维数为45的极化特征, LPP先通过降维进一步提取特征再使用SVM做分类.这三种方法都是随机选取1 %的样本作为训练集, 其余作为测试集. RV-CNN和CV-CNN使用0.9 %的样本做训练, 0.1 %的样本做验证, 实验结果来自文献[80].
实验使用ESAR系统获取的Oberfaffenhofen地区的全极化数据, 尺寸为$ 1 300\times 1 200 $, 其Pauli分解图见图 3 (a).图 3 (b)展示了对应的真实分类图, 其中黑色、灰色和白色区域分别对应城市、林地和其他开放区域.图 3 (c) $ \sim $ 3 (g)展示了不同方法的分类图的对比结果, 表 3给出了这些方法的分类准确率的对比结果, 从表 3可以看出, SWC的分类精度最低, 其在城市区域有超过一半的像素被错分了.从图 3 (c)中也可以看到很多的黑色区域被错分成灰色和白色区域. LPP的结果要优于SVM, 因为SVM直接使用了45维的极化特征进行分类, 而LPP先提取更有判别性的低维本质特征再用SVM分类. RV-CNN的结果仅优于SWC, 因为它忽略了PolSAR数据的相位信息.与其他方法相比, CV-CNN的整体分类精度最高, 特别地, 其在城市区域的分类精度比其他方法高近6 %, 这也体现了深度神经网络对复杂环境特征提取的优势.
表 3 不同PolSAR图像分类方法的分类精度比较(%)Table 3 Classification accuracy of different PolSAR image classification methods (%)类别/方法 SWC SVM LPP RV-CNN CV-CNN 城市区域 49.4 78.7 82.5 85.7 91.3 林地区域 80.2 93.8 93.8 85.2 92.2 其他区域 95.6 95.9 96.5 93.4 94.6 总精度 81.2 91.2 92.5 89.9 93.4 5. 总结与展望
5.1 PolSAR相干斑抑制
针对PolSAR相干斑抑制问题, 本文在第2节总结了目前主要的三类方法:空域滤波、非局部均值和变分方法, 并对其中的几个代表性方法进行了实验比较.其中, 空域滤波方法操作简单、速度快, 但会降低图像空间分辨率并导致图像的边缘和细节模糊.非局部均值方法在相干斑抑制和分辨率保持上都优于空域滤波方法, 能显著抑制斑点噪声并很好保持图像的重复性结构特征, 但计算速度很慢而且会出现像素块状效应.变分方法在较好去噪的同时能很好保持图像的边缘特征, 并且计算速度相对较快, 但会导致图像中一些细节的模糊.总的来说, 每类方法都有各自的优缺点, 因此需要根据实际应用的需求选择特定的斑点抑制算法.
相干斑抑制作为一个预处理过程, 对PolSAR数据后续的解译影响很大.因此, 去相干斑算法需要能很好保持原数据的极化特征和图像的结构特征, 以保证后续分类或识别的精度.另外, 随着PolSAR数据的规模越来越大, 分辨率越来越高, 也需要去相干斑算法具有较高的处理效率.关于未来可能的研究方向, 我们认为可从以下几个方面进行考虑:
1) 设计快速且效果好的相干斑抑制算法.在目前的相干斑抑制算法中, 计算复杂度低的算法重建的图像分辨率往往不够高或细节保持的不好, 而重构图像质量高的算法往往计算复杂度也比较高而且需要调整参数.因此设计去相干斑效果好并且处理速度快的算法是实用化的关键.可以借助对PolSAR数据的噪声模型、统计分布或极化分解的最新研究成果, 结合现有的优势方法, 构建更有效的滤波算法; 或充分利用PolSAR分类的最新结果, 设计类内滤波的相干斑抑制算法, 对算法进行并行化处理以及设计参数自动选取策略, 提升计算效率.
2) 设计具有自适应去相干斑能力的算法.不同的PolSAR系统获取的图像往往分辨率不同, 噪声程度也不同, 即便同一PolSAR系统, 搭载不同波段的传感器获取的数据也差别很大.另外, 由于地形地貌千变万化, PolSAR数据对应的场景也往往各种各样.因此, 针对不同的传感器、不同分辨率、不同程度的噪声以及不同类型的场景, 设计能自适应地去相干斑的算法是未来的研究重点.不过设计这样通用且高效的算法往往比较困难, 可以针对特定程度的噪声或特定类型的场景, 充分利用先验信息, 设计高效的算法, 然后进行集成, 根据场景及应用需求自动选择特定方法.
3) 最新的图像处理或机器学习技术将为PolSAR数据斑点抑制提供新的途径.例如深度自编码和深度卷积神经网络等技术在光学图像的去噪中达到了非常好的效果.然而, PolSAR数据不像光学图像有真实数据做参考, 并受数据量的限制, 现有的这些新技术很难直接应用到PolSAR数据处理中.另外, 由于PolSAR数据的噪声是乘性和加性的混合, 并且涉及到复数据, 这也为最新技术的推广增加了困难.能否通过模拟的PolSAR数据进行网络训练或将在光学图像上训练的深度神经网络通过迁移学习用于PolSAR数据的去噪中, 还有待科研工作者的深入研究.
5.2 PolSAR分类
本文在第3节归纳总结了目前PolSAR数据的常用特征, 包括极化特征、统计特征、纹理特征、颜色特征、空间结构特征和基于深度学习提取的特征.在第4节总结了目前PolSAR图像的三大类分类方法:有监督学习、无监督学习和半监督学习, 并对每个大类中的方法进行了详细地归纳和分类.
随着PolSAR系统的发展和应用, 获取到的PolSAR数据的规模越来越大, 对海量PolSAR数据进行快速准确地分类是未来研究的重点问题.目前的分类方法主要集中于特征提取和分类器设计两个方面.在特征提取方面, 如何结合PolSAR数据的特点, 利用深度神经网络更有效地自动提取特征, 提取能有效反映散射体类型的新特征, 以及对冗余特征进行选择以得到低维的更有判别性的特征仍是研究的重点; 在分类器设计方面, 如何借鉴最新的分类技术或使用多分类器的集成学习等方式, 构造精度高、鲁棒性和泛化性强的分类算法也一直是研究的重点.除此之外, 关于未来可能的研究方向, 我们认为可从以下几个方面进行考虑:
1) 实现少量标注样本情况下的高精度分类. PolSAR图像不像光学数据那么直观, 因此判读困难, 需要有经验的专家进行标注, 这个过程耗时耗力.因此, 如何实现在只有少量标注样本情况下的高精度分类是研究的热点问题, 半监督学习和主动学习为此提供了解决途径.半监督学习充分利用有标注和大量无标注样本的信息, 能在小样本情况下获得不错的分类效果; 主动学习可以从大量样本中选择信息最丰富的样本进行标注, 以实现对尽可能少的样本进行标注而达到尽可能好的分类效果.目前, 针对PolSAR图像的这方面研究工作并不丰富, 但具有重要的应用价值, 值得深入研究.
2) 设计基于多视角信息融合的分类算法.多视角数据可以是从不同传感器获取的信息或多种属性的特征.目前的一些PolSAR系统可以提供多个波段的极化数据, 甚至一些平台可以搭载除PolSAR外的其他传感器, 从而可以获取大量的多视角数据.另外, PolSAR图像可以通过极化、纹理和颜色等多视角特征集进行刻画, 目前的PolSAR分类方法多数是基于单视角或某几个视角数据的, 但往往都是排成一个高维向量处理, 忽视了特征的不同属性和视角间的关系.如何充分利用这些视角间的冗余性、互补性和一致性, 设计有效的信息融合技术实现多种PolSAR数据源、多传感器数据或多种属性的特征的融合处理, 以提高分类精度也是未来研究的重要方向.
3) 设计在线增量算法用于PolSAR图像实时分类.在实际应用中, PolSAR数据并非是静态的, 而是动态数据流, 即PolSAR系统会不断获取新的数据, 因此样本随着时间的累计不断增长.目前的PolSAR分类方法基本上都是离线学习算法, 它们假定数据是静态的且在整个训练集上学习模型并用于后续的分类, 当有新的样本到来时不得不重新训练模型, 这将非常耗时.而在线增量学习算法能高效地从动态变化的大规模数据流中增量地学习分类模型, 并且无需使用之前的数据对模型重训练.在线增量算法能从新来样本中不断学习新信息更新分类器, 并且不会遗忘先前学出的分类模型, 具有很强的动态环境的自适应性, 这将为PolSAR图像的实时分类提供可能, 是未来值得探索的研究方向.
-
表 1 各类PolSAR数据相干斑抑制方法比较表
Table 1 Comparison of different types of PolSAR data speckle reduction methods
相干斑抑制方法 代表性算法 优点 不足 空域滤波方法 修正的Lee、IDAN、改进的sigma滤波 操作简单、计算速度快 空间分辨率低、边缘模糊、细节信息丢失 非局部均值方法 Pretest、随机散度滤波、NL-SAR 去噪效果显著、结构特征保持得好 计算速度很慢、出现像素块状效应 变分方法 WisTV、WisNLTV、MuLoG算法 去噪显著、速度快、边缘保持得好 弱目标丢失、平坦区域出现阶梯效应 表 2 不同PolSAR数据相干斑抑制方法的量化比较表
Table 2 Quantitative comparison of different PolSAR data speckle reduction methods
指标/方法 修正的Lee滤波 SMBF滤波 NL-SAR方法 WisNLTV方法 ENL 13.8698 10.9144 11.5435 10.7267 EPD-ROA水平方向 0.6863 0.7356 0.7505 0.7533 EPD-ROA垂直方向 0.7472 0.7957 0.8046 0.8109 表 3 不同PolSAR图像分类方法的分类精度比较(%)
Table 3 Classification accuracy of different PolSAR image classification methods (%)
类别/方法 SWC SVM LPP RV-CNN CV-CNN 城市区域 49.4 78.7 82.5 85.7 91.3 林地区域 80.2 93.8 93.8 85.2 92.2 其他区域 95.6 95.9 96.5 93.4 94.6 总精度 81.2 91.2 92.5 89.9 93.4 -
[1] Oliver C, Quegan S. Understanding Synthetic Aperture Radar Images. Raleigh, NC, USA: SciTech Publishing, 2004. [2] Lee J, Pottier E. Polarimetric Radar Imaging: From Basics to Applications. Florida, USA: CRC Press, 2009. [3] 周晓光, 匡纲要, 万建伟.极化SAR图像分类综述.信号处理, 2008, 24(5): 806-812 doi: 10.3969/j.issn.1003-0530.2008.05.023Zhou Xiao-Guang, Kuang Gang-Yao, Wan Jian-Wei. A review of polarimetric SAR image classification. Signal Processing, 2008, 24(5): 806-812 doi: 10.3969/j.issn.1003-0530.2008.05.023 [4] Goodman J W. Some fundamental properties of speckle. Journal of the Optical Society of America, 1976, 66(11): 1145-1150 doi: 10.1364/JOSA.66.001145 [5] Touzi R, Lopes A. The principle of speckle filtering in polarimetric SAR imagery. IEEE Transactions on Geoscience and Remote Sensing, 1994, 32(5): 1110-1114 doi: 10.1109/36.312901 [6] Lopez-Martinez C, Fabregas X. Polarimetric SAR speckle noise model. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41(10): 2232-2242 doi: 10.1109/TGRS.2003.815240 [7] Lee J S, Grunes M R, De Grandi G. Polarimetric SAR speckle filtering and its implication for classification. IEEE Transactions on Geoscience and Remote Sensing, 1999, 37(5): 2363-2373 doi: 10.1109/36.789635 [8] Lopez-Martinez C, Fabregas X. Model-based polarimetric SAR speckle filter. IEEE Transactions on Geoscience and Remote Sensing, 2008, 46(11): 3894-3907 doi: 10.1109/TGRS.2008.2002029 [9] Foucher S, López-Martínez C. Analysis, evaluation, and comparison of polarimetric SAR speckle filtering techniques. IEEE Transactions on Image Processing, 2014, 23(4): 1751-1764 doi: 10.1109/TIP.2014.2307437 [10] van Zyl J J, Zebker H A, Elachi C. Imaging radar polarization signatures: theory and observation. Radio Science, 1987, 22(4): 529-543 doi: 10.1029/RS022i004p00529 [11] Xie H, Pierce L E, Ulaby F T. SAR speckle reduction using wavelet denoising and Markov random field modeling. IEEE Transactions on Geoscience and Remote Sensing, 2002, 40(10): 2196-2212 doi: 10.1109/TGRS.2002.802473 [12] Feng H X, Hou B, Gong M G. SAR image despeckling based on local homogeneous-region segmentation by using pixel-relativity measurement. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(7): 2724-2737 doi: 10.1109/TGRS.2011.2107915 [13] 周晓光, 匡纲要, 万建伟.多极化SAR图像斑点抑制综述.中国图象图形学报, 2008, 13(3): 377-385 http://d.old.wanfangdata.com.cn/Periodical/zgtxtxxb-a200803001Zhou Xiao-Guang, Kuang Gang-Yao, Wan Jian-Wei. A review of polarimetric SAR speckle reduction. Journal of Image and Graphics, 2008, 13(3): 377-385 http://d.old.wanfangdata.com.cn/Periodical/zgtxtxxb-a200803001 [14] Gonzalez R C, Woods R E. Digital Image Processing (3rd edition). Upper Saddle River, NJ: Prentice-Hall, 2012. [15] Lee J S, Grunes M R, Schuler D L, Pottier E, Ferro-Famil L. Scattering-model-based speckle filtering of polarimetric SAR data. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(1): 176-187 doi: 10.1109/TGRS.2005.859338 [16] Vasile G, Trouve E, Lee J S, Buzuloiu V. Intensity-driven adaptive-neighborhood technique for polarimetric and interferometric SAR parameters estimation. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(6): 1609-1621 doi: 10.1109/TGRS.2005.864142 [17] Deledalle C A, Denis L, Tupin F, Reigber A, Jäger M. NL-SAR: a unified nonlocal framework for resolution-preserving (Pol)(In)SAR denoising. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(4): 2021-2038 doi: 10.1109/TGRS.2014.2352555 [18] Ding Z G, Zeng T, Dong F, Liu L S, Yang W F, Long T. An improved PolSAR image speckle reduction algorithm based on structural judgment and hybrid four-component polarimetric decomposition. IEEE Transactions on Geo- science and Remote Sensing, 2013, 51(8): 4438-4449 doi: 10.1109/TGRS.2013.2270036 [19] D$'$Hondt O, Guillaso S, Hellwich O. Iterative bilateral filtering of polarimetric SAR data. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2013, 6(3): 1628-1639 doi: 10.1109/JSTARS.2013.2256881 [20] Lang F K, Yang J, Li D R. Adaptive-window polarimetric SAR image speckle filtering based on a homogeneity measurement. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(10): 5435-5446 doi: 10.1109/TGRS.2015.2422737 [21] Lee J S, Ainsworth T L, Wang Y T, Chen K -S. Polarimetric SAR speckle filtering and the extended sigma filter. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(3): 1150-1160 doi: 10.1109/TGRS.2014.2335114 [22] Alonso-González A, López-Martínez C, Salembier P. PolSAR time series processing with binary partition trees. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(6): 3553-3567 doi: 10.1109/TGRS.2013.2273664 [23] Salehi M, Mohammadzadeh A, Maghsoudi Y. Adaptive speckle filtering for time series of polarimetric SAR images. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 10(6): 2841-2848 doi: 10.1109/JSTARS.2017.2698020 [24] Buades A, Coll B, Morel J M. A non-local algorithm for image denoising. In: Proceedings of the 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition. San Diego, CA, USA: IEEE, 2005. 60-65 https://www.iro.umontreal.ca/~mignotte/IFT6150/Articles/Buades-NonLocal.pdf [25] Chen J, Chen Y L, An W T, Cui Y, Yang J. Nonlocal filtering for polarimetric SAR data: a pretest approach. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(5): 1744-1754 doi: 10.1109/TGRS.2010.2087763 [26] Zhong H, Zhang J J, Liu G C. Robust polarimetric SAR despeckling based on nonlocal means and distributed lee filter. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(7): 4198-4210 doi: 10.1109/TGRS.2013.2280278 [27] Torres L, Sant$'$Anna S J S, Da Costa Freitas C, Frery A C. Speckle reduction in polarimetric SAR imagery with stochastic distances and nonlocal means. Pattern Recognition, 2014, 47(1): 141-157 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=b00fddd7cd8d22e065568b3eb6bd3c28 [28] Wang Y T, Ainsworth T L, Lee J S. Application of mixture regression for improved polarimetric SAR speckle filtering. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(1): 453-467 doi: 10.1109/TGRS.2016.2608501 [29] Ma X S, Wu P H, Shen H F. Multifrequency polarimetric SAR image despeckling by iterative nonlocal means based on a space-frequency information joint covariance matrix. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2019, 12(1): 274-284 doi: 10.1109/JSTARS.2018.2884194 [30] Sharma R, Panigrahi R K. CFAR-based adaptive PolSAR speckle filter. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2018, 11(12): 4895-4905 doi: 10.1109/JSTARS.2018.2877670 [31] Jordan C. Sur la série de Fourier. Comptes Rendus de l'Académie des Sciences Paris, 1881, 92(5): 228-230 http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ021418533/ [32] Rudin L I, Osher S, Fatemi E. Nonlinear total variation based noise removal algorithms. Physica D: Nonlinear Phenomena, 1992, 60(1-4): 259-268 doi: 10.1016/0167-2789(92)90242-F [33] Li H, Qin Z Y. PolSAR image speckle reduction algorithm based on TV-PDE. In: Proceedings of the SPIE 8002, MIPPR 2011: Multispectral Image Acquisition, Processing, and Analysis. Guilin, China: SPIE, 2011, 8002: Article No. 80020S http://adsabs.harvard.edu/abs/2011SPIE.8002E..0SL [34] Liu H Y, Yan F X, Zhu J B, Fang F M. Adaptive vectorial total variation models for multi-channel synthetic aperture radar images despeckling with fast algorithms. IET Image Processing, 2013, 7(9): 795-804 doi: 10.1049/iet-ipr.2013.0177 [35] Nie X L, Qiao H, Zhang B. A variational model for PolSAR data speckle reduction based on the Wishart distribution. IEEE Transactions on Image Processing, 2015, 24(4): 1209-1222 doi: 10.1109/TIP.2015.2396292 [36] Nie X L, Qiao H, Zhang B, Huang X Y. A nonlocal TV- based variational method for PolSAR data speckle reduction. IEEE Transactions on Image Processing, 2016, 25(6): 2620-2634 doi: 10.1109/TIP.2016.2552402 [37] Nie X L, Zhang B, Chen Y J, Qiao H. A new algorithm for optimizing TV-based PolSAR despeckling model. IEEE Signal Processing Letters, 2016, 23(10): 1409-1413 doi: 10.1109/LSP.2016.2602299 [38] Deledalle C -A, Denis L, Tabti S, Tupin F. MuLoG, or how to apply Gaussian denoisers to multi-channel SAR speckle reduction? IEEE Transactions on Image Processing, 2017, 26(9): 4389-4403 doi: 10.1109/TIP.2017.2713946 [39] Farage G, Foucher S, Benie G B. PolSAR image filtering based on feature detection using the wavelet transform. In: Proceedings of the 2007 IEEE International Geoscience and Remote Sensing Symposium IGARSS. Barcelona, Spain: IEEE, 2007, 2648-2652 https://www.researchgate.net/publication/224301188_PolSAR_image_filtering_based_on_feature_detection_using_the_wavelet_transform [40] Foucher S, Farage G, Bénié G B. Polarimetric SAR image filtering with trace-based partial differential equations. Canadian Journal of Remote Sensing, 2007, 33(3): 226- 236 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.5589/m07-023 [41] Sun S, Tian J W, Liu Z S, Cai N. Anisotropic diffusion for speckle filtering of polarimetric synthetic aperture radar imagery. Journal of Electronic Imaging, 2013, 22(1): Article No. 013003003-013 https://www.researchgate.net/publication/236116496_Anisotropic_Diffusion_for_Speckle_Filtering_of_Polarimetric_Synthetic_Aperture_Radar_Imagery [42] Ma X S, Wu P H, Wu Y L, Shen H F. A review on recent developments in fully polarimetric SAR image despeckling. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2018, 11(3): 743-758 doi: 10.1109/JSTARS.2017.2768059 [43] Tu S T, Chen J Y, Yang W, Sun H. Laplacian eigenmaps-based polarimetric dimensionality reduction for SAR image classification. IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(1): 170-179 doi: 10.1109/TGRS.2011.2168532 [44] Huynen J R. Phenomenological Theory of Radar Targets[Ph. D. dissertation], Technical University of Delft, The Netherlands, 1970. https://repository.tudelft.nl/islandora/object/uuid:e4a140a0-c175-45a7-ad41-29b28361b426?collection=research [45] Yang J, Peng Y N, Yamaguchi Y, Yamada H. On Huynen$'$s decomposition of a Kennaugh matrix. IEEE Geoscience and Remote Sensing Letters, 2006, 3(3): 369- 372 doi: 10.1109/LGRS.2006.873229 [46] Freeman A, Durden S L. A three-component scattering model for polarimetric SAR data. IEEE Transactions on Geoscience and Remote Sensing, 1998, 36(3): 963-973 doi: 10.1109/36.673687 [47] Yamaguchi Y, Moriyama T, Ishido M, Yamada H. Four-component scattering model for polarimetric SAR image decomposition. IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(8): 1699-1706 doi: 10.1109/TGRS.2005.852084 [48] Freeman A. Fitting a two-component scattering model to polarimetric SAR data from forests. IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(8): 2583-2592 doi: 10.1109/TGRS.2007.897929 [49] Kusano S, Takahashi K, Sato M. A new decomposition of a PolSAR coherency matrix using a generalized scattering model. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(8): 3933-3940 doi: 10.1109/JSTARS.2014.2367540 [50] Cloude S R, Pottier E. A review of target decomposition theorems in radar polarimetry. IEEE Transactions on Geoscience and Remote Sensing, 1996, 34(2): 498-518 doi: 10.1109/36.485127 [51] van Zyl J J. Application of cloude$'$s target decomposition theorem to polarimetric imaging radar data. In: Proceedings of the SPIE — the International Society for Optical Engineering. San Diego, CA, USA: SPIE, 1993. 184-191 [52] Krogager E. New decomposition of the radar target scattering matrix. Electronics Letters, 1990, 26(18): 1525- 1527 doi: 10.1049/el:19900979 [53] Cameron W L, Rais H. Conservative polarimetric scatterers and their role in incorrect extensions of the Cameron decomposition. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(12): 3506-3516 doi: 10.1109/TGRS.2006.879115 [54] Touzi R. Target scattering decomposition in terms of roll-invariant target parameters. IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(1): 73-84 doi: 10.1109/TGRS.2006.886176 [55] Goodman N R. Statistical analysis based on a certain multivariate complex Gaussian distribution (an introduction). The Annals of Mathematical Statistics, 1963, 34(1): 152- 177 doi: 10.1214-aoms-1177704250/ [56] Papoulis A. Probability, Random Variables, and Stochastic Processes. NewYork: McGraw-Hill College, 1965. [57] Freitas C C, Frery A C, Correia A H. The polarimetric G distribution for SAR data analysis. Environmetrics, 2005, 16(1): 13-31 doi: 10.1002/env.658 [58] Khan S, Guida R. On single-look multivariate G distribution for PolSAR data. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2012, 5(4): 1149-1163 doi: 10.1109/JSTARS.2012.2202557 [59] Doulgeris A P, Anfinsen S N, Eltoft T. Classification with a non-Gaussian model for PolSAR data. IEEE Transactions on Geoscience and Remote Sensing, 2008, 46(10): 2999- 3009 doi: 10.1109/TGRS.2008.923025 [60] Bombrun L, Beaulieu J -M. Fisher distribution for texture modeling of polarimetric SAR data. IEEE Geoscience and Remote Sensing Letters, 2008, 5(3): 512-516 doi: 10.1109/LGRS.2008.923262 [61] Frery A C, Cintra R J, Nascimento A D C. Entropy-based statistical analysis of PolSAR data. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(6): 3733-3743 doi: 10.1109/TGRS.2012.2222029 [62] Uhlmann S, Kiranyaz S. Integrating color features in polarimetric SAR image classification. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(4): 2197-2216 doi: 10.1109/TGRS.2013.2258675 [63] Dai D X, Yang W, Sun H. Multilevel local pattern histogram for SAR image classification. IEEE Geoscience and Remote Sensing Letters, 2011, 8(2): 225-229 doi: 10.1109/LGRS.2010.2058997 [64] Masjedi A, Zoej M J V, Maghsoudi Y. Classification of polarimetric SAR images based on modeling contextual information and using texture features. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(2): 932-943 doi: 10.1109/TGRS.2015.2469691 [65] Liu B, Zhang Z H, Liu X Z, Yu W X. Representation and spatially adaptive segmentation for PolSAR images based on wedgelet analysis. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(9): 4797-4809 doi: 10.1109/TGRS.2015.2410177 [66] Xing Y X, Zhang Y, Li N, Wang R, Hu G X. Improved superpixel-based polarimetric synthetic aperture radar image classification integrating color features. Journal of Applied Remote Sensing, 2016, 10(2): Article No. 026026- 026 https://www.researchgate.net/publication/303594157_Improved_superpixel-based_polarimetric_synthetic_aperture_radar_image_classification_integrating_color_features [67] Liu H Y, Wang S, Wang R F, Shi J F, Zhang E L, Yang S Y, Jiao L. A framework for classification of urban areas using polarimetric SAR images integrating color features and statistical model. Journal of Infrared and Millimeter Waves, 2016, 35(4): 398-406 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hwyhmb201604004 [68] Liu F, Shi J F, Jiao L C, Liu H Y, Yang S Y, Wu J, Hao H X, Yuan J L. Hierarchical semantic model and scattering mechanism based PolSAR image classification. Pattern Recognition, 2016, 59: 325-342 doi: 10.1016/j.patcog.2016.02.020 [69] Xiang D L, Ban Y F, Wang W, Su Y. Adaptive superpixel generation for polarimetric SAR images with local iterative clustering and SIRV model. IEEE Transcactions on Geoscience and Remote Sensing, 2017, 55(6): 3115-3131 doi: 10.1109/TGRS.2017.2662010 [70] Liu B, Hu H, Wang H Y, Wang K Z, Liu X Z, Yu W X. Superpixel-based classification with an adaptive number of classes for polarimetric SAR images. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(2): 907-924 doi: 10.1109/TGRS.2012.2203358 [71] Wu Y H, Ji K F, Yu W X, Su Y. Region-based classification of polarimetric SAR images using Wishart MRF. IEEE Geoscience and Remote Sensing Letters, 2008, 5(4): 668- 672 doi: 10.1109/LGRS.2008.2002263 [72] Ersahin K, Cumming I G, Ward R K. Segmentation and classification of polarimetric SAR data using spectral graph partitioning. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(1): 164-174 doi: 10.1109/TGRS.2009.2024303 [73] Du L J, Lee J S, Hoppel K, Mango S A. Segmentation of SAR images using the wavelet transform. International Journal of Imaging Systems and Technology, 1992, 4(4): 319-326 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1002/ima.1850040411 [74] Rignot E, Chellappa R. Segmentation of polarimetric synthetic aperture radar data. IEEE Transactions on Image Processing, 1992, 1(3): 281-300 doi: 10.1109/83.148603 [75] Ren B, Hou B, Zhao J, Jiao L C. Unsupervised classification of polarimetirc SAR image via improved manifold regularized low-rank representation with multiple features. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 10(2): 580-595 doi: 10.1109/JSTARS.2016.2573380 [76] Kim H, Hirose A. Unsupervised hierarchical land classification using self-organizing feature codebook for decimeter-resolution PolSAR. IEEE Transactions on Geoscience and Remote Sensing, 2019, 57(4): 1894-1905 doi: 10.1109/TGRS.2018.2870134 [77] 石俊飞, 刘芳, 林耀海, 刘璐.基于深度学习和层次语义模型的极化SAR分类.自动化学报, 2017, 43(2): 215-2265 http://www.aas.net.cn/CN/abstract/abstract19010.shtmlShi Jun-Fei, Liu Fang, Lin Yao-Hai, Liu Lu. Polarimetric SAR image classification based on deep learning and hierarchical semantic model. Acta Automatica Sinica, 2017, 43(2): 215-2265 http://www.aas.net.cn/CN/abstract/abstract19010.shtml [78] Liu H Y, Wang Y K, Yang S Y, Wang S, Feng J, Jiao L C. Large polarimetric SAR data semi-supervised classification with spatial-anchor graph. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(4): 1439-1458 doi: 10.1109/JSTARS.2016.2518675 [79] Zhou Y, Wang H P, Xu F, Jin Y Q. Polarimetric SAR image classification using deep convolutional neural networks. IEEE Geoscience and Remote Sensing Letters, 2016, 13(12): 1935-1939 doi: 10.1109/LGRS.2016.2618840 [80] Zhang Z M, Wang H P, Xu F, Jin Y Q. Complex-valued convolutional neural network and its application in polarimetricv SAR image classification. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(12): 7177-7188 doi: 10.1109/TGRS.2017.2743222 [81] Liu F, Jiao L C, Hou B, Yang S Y. Pol-SAR image classification based on Wishart DBN and local spatial information. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(6): 3292-3308 doi: 10.1109/TGRS.2016.2514504 [82] Hou B, Kou H D, Jiao L C. Classification of polarimetric SAR images using multilayer autoencoders and superpixels. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(7): 3072-3081 doi: 10.1109/JSTARS.2016.2553104 [83] Lee J S, Hoppel K W, Mango S A, Miller A R. Intensity and phase statistics of multilook polarimetric and interferometric SAR imagery. IEEE Transactions on Geoscience and Remote Sensing, 1994, 32(5): 1017-1028 doi: 10.1109/36.312890 [84] Frery A C, Correia A H, Da Freitas C. Classifying multifrequency fully polarimetric imagery with multiple sources of statistical evidence and contextual information. IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(10): 3098-3109 doi: 10.1109/TGRS.2007.903828 [85] Gou S P, Qiao X, Zhang X R, Wang W F, Du F F. Eigenvalue analysis-based approach for POL-SAR image classification. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(2): 805-818 doi: 10.1109/TGRS.2013.2244096 [86] Alberga V, Borghys D, Satalino G, Staykova D K, Borghgraef A, Lapierre F, Perneel C, et al. Comparison of algorithms for the classification of polarimetric SAR data. In: Proceedings of the 2009 SPIE 7477 — Image and Signal Processing for Remote Sensing XV. Berlin, Germany: SPIE, 2009. Article No. 7477-74771V https://www.researchgate.net/publication/228633764_Comparison_of_algorithms_for_the_classification_of_polarimetric_SAR_data [87] Kajimoto M, Susaki J. Urban-area extraction from polarimetric SAR images using polarization orientation angle. IEEE Geoscience and Remote Sensing Letters, 2013, 10(2): 337-341 doi: 10.1109/LGRS.2012.2207085 [88] Yang W, Yin X S, Song H, Liu Y, Xu X. Extraction of built-up areas from fully polarimetric SAR imagery via PU learning. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014, 7(4): 1207-1216 doi: 10.1109/JSTARS.2013.2289986 [89] Banerjee B, Bhattacharya A, Buddhiraju K M. A generic land-cover classification framework for polarimetric SAR images using the optimum Touzi decomposition parameter subset-an insight on mutual information-based feature selection techniques. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014, 7(4): 1167-1176 doi: 10.1109/JSTARS.2014.2304456 [90] Jafari M, Maghsoudi Y, Zoej M J V. A new method for land cover characterization and classification of polarimetric SAR data using polarimetric signatures. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(7): 3595-3607 doi: 10.1109/JSTARS.2014.2387374 [91] Fukuda S, Hirosawa H. A wavelet-based texture feature set applied to classification of multifrequency polarimetric SAR images. IEEE Transactions on Geoscience and Remote Sensing, 1999, 37(5): 2282-2286 doi: 10.1109/36.789624 [92] He C, Li S, Liao Z X, Liao M S. Texture classification of PolSAR data based on sparse coding of wavelet polarization textons. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(8): 4576-4590 doi: 10.1109/TGRS.2012.2236338 [93] Zhang L M, Zou B, Zhang J P, Zhang Y. Classification of polarimetric SAR image based on support vector machine using multiple-component scattering model and texture features. EURASIP Journal on Advances in Signal Processing, 2010(1): Article No. 960831 doi: 10.1155%2F2010%2F960831 [94] Wu W J, Guo H D, Li X W. Urban area SAR image man-made target extraction based on the product model and the time-frequency analysis. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(3): 943-952 doi: 10.1109/JSTARS.2014.2371064 [95] Shi L, Zhang L F, Yang J, Zhang L P, Li P X. Supervised graph embedding for polarimetric SAR image classification. IEEE Geoscience and Remote Sensing Letters, 2013, 10(2): 216-220 doi: 10.1109/LGRS.2012.2198612 [96] Tao M L, Zhou F, Liu Y, Zhang Z J. Tensorial independent component analysis-based feature extraction for polarimetric SAR data classification. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(5): 2481-2495 doi: 10.1109/TGRS.2014.2360943 [97] Huang X Y, Qiao H, Zhang B, Nie X L. Supervised polarimetric SAR image classification using tensor local discriminant embedding. IEEE Transactions on Image Processing, 2018, 27(6): 2966-2979 doi: 10.1109/TIP.2018.2815759 [98] Nie X L, Ding S G, Zhang B, Qiao H, Huang X Y. PolSAR data online classification based on multi-view learning. In: Proceedings of the 24th IEEE International Conference on Image Processing (ICIP). Beijing, China: IEEE, 2017. 2354-2358 https://www.researchgate.net/publication/323349225_Polsar_data_online_classification_based_on_multi-view_learning [99] Nie X L, Luo Y K, Qiao H, Zhang B, Jiang Z P. An incremental multi-view active learning algorithm for PolSAR data classification. In: Proceedings of the 24th International Conference on Pattern Recognition (ICPR). Beijing, China: IEEE, 2018. 2251-2255 https://nyuscholars.nyu.edu/en/publications/an-incremental-multi-view-active-learning-algorithm-for-polsar-da [100] Nie X L, Ding S G, Huang X Y, Qiao H, Zhang B, Jiang Z P. An online multi-view learning algorithm for PolSAR data real-time classification. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2019, 12(1): 302-320 doi: 10.1109/JSTARS.2018.2886821 [101] Huang X Y, Zhang B, Qiao H, Nie X L. Local discriminant canonical correlation analysis for supervised PolSAR image classification. IEEE Geoscience and Remote Sensing Letters, 2017, 14(11): 2102-2106 doi: 10.1109/LGRS.2017.2752800 [102] Kong J, Swartz A A, Yueh H A, Novak L M, Shin R T. Identification of terrain cover using the optimum polarimetric classifier. Journal of Electromagnetic Waves and Applications, 1988, 2(2): 171-194 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1163/156939387X00324 [103] Lee J S, Grunes M R, Kwok R. Classification of multi-look polarimetric SAR imagery based on complex Wishart distribution. International Journal of Remote Sensing, 1994, 15(11): 2299-2311 doi: 10.1080/01431169408954244 [104] Pottier E, Saillard J. On radar polarization target decomposition theorems with application to target classification, by using neural network method. In: Proceedings of the 7th International Conference on Antennas and Propagation (ICAP). York, UK, USA: IET, 1991. 265-268 [105] Chen K S, Huang W P, Tsay D H, Amar F. Classification of multifrequency polarimetric SAR imagery using a dynamic learning neural network. IEEE Transactions on Geoscience and Remote Sensing, 1996, 34(3): 814-820 doi: 10.1109/36.499786 [106] Tzeng Y C, Chen K S. A fuzzy neural network to SAR image classification. IEEE Transactions on Geoscience and remote Sensing, 1998, 36(1): 301-307 doi: 10.1109/36.655339 [107] Ince T, Kiranyaz S, Gabbouj M. Evolutionary RBF classifier for polarimetric SAR images. Expert Systems with Applications, 2012, 39(5): 4710-4717 doi: 10.1016/j.eswa.2011.09.082 [108] Shang F, Hirose A. Quaternion neural-network-based PolSAR land classification in Poincare-sphere-parameter space. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(9): 5693-5703 doi: 10.1109/TGRS.2013.2291940 [109] Antropov O, Rauste Y, Astola H, Praks J, Häme T, Hallikainen M T. Land cover and soil type mapping from spaceborne PolSAR data at L-band with probabilistic neural network. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(9): 5256-5270 doi: 10.1109/TGRS.2013.2287712 [110] Fukuda S, Hirosawa H. Support vector machine classification of land cover: application to polarimetric SAR data. In: Proceedings of the 2001 International Geoscience and Remote Sensing Symposium. Scanning the Present and Resolving the Future. Sydney, NSW, Australia: IEEE, 2001. 187-189 https://ieeexplore.ieee.org/document/976097 [111] Lardeux C, Frison P L, Tison C, Souyris J C, Stoll B, Fruneau B, Rudant J P, et al. Support vector machine for multifrequency SAR polarimetric data classification. IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(12): 4143-4152 doi: 10.1109/TGRS.2009.2023908 [112] Maghsoudi Y, Collins M J, Leckie D G. Radarsat-2 polarimetric SAR data for boreal forest classification using SVM and a wrapper feature selector. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2013, 6(3): 1531-1538 doi: 10.1109/JSTARS.2013.2259219 [113] Zhang L M, Sun L J, Zou B, Moon W M. Fully polarimetric SAR image classification via sparse representation and polarimetric features. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(8): 3923-3932 doi: 10.1109/JSTARS.2014.2359459 [114] 文伟, 王英华, 冯博, 刘宏伟.基于监督非相干字典学习的极化SAR图像舰船目标检测.自动化学报, 2015, 41(11): 1926- 1940 http://www.aas.net.cn/CN/abstract/abstract18767.shtmlWen Wei, Wang Ying-Hua, Feng Bo, Liu Hong-Wei. Supervised incoherent dictionary learning for ship detection with PolSAR images. Acta Automatica Sinica, 2015, 41(11): 1926-1940 http://www.aas.net.cn/CN/abstract/abstract18767.shtml [115] 何楚, 刘明, 冯倩, 邓新萍.基于多尺度压缩感知金字塔的极化干涉SAR图像分类.自动化学报, 2011, 37(7): 820-827 http://www.aas.net.cn/CN/abstract/abstract17450.shtmlHe Chu, Liu Ming, Feng Qian, Deng Xin-Ping. PolInSAR Image classification based on compressed sensing and multi-scale pyramid. Acta Automatica Sinica, 2011, 37(7): 820-827 http://www.aas.net.cn/CN/abstract/abstract17450.shtml [116] She X L, Yang J, Zhang W J. The boosting algorithm with application to polarimetric SAR image classification. In: Proceedings of the 1st Asian and Pacific Conference on Synthetic Aperture Radar. Huangshan, China: IEEE, 2007. 779-783 https://www.researchgate.net/publication/224300179_The_boosting_algorithm_with_application_to_polarimetric_SAR_image_classification [117] Zou T Y, Yang W, Dai D X, Sun H. Polarimetric SAR image classification using multifeatures combination and extremely randomized clustering forests. EURASIP Journal on Advances in Signal Processing, 2010, 2010(4): Article No. 4 doi: 10.1155/2010/465612 [118] Mishra P, Singh D. A statistical-measure-based adaptive land cover classification algorithm by efficient utilization of polarimetric SAR observables. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(5): 2889-2900 doi: 10.1109/TGRS.2013.2267548 [119] Gomez L, Alvarez L, Mazorra L, Frery A C. Fully PolSAR image classification using machine learning techniques and reaction-diffusion systems. Neurocomputing, 2017, 255: 52-60 doi: 10.1016/j.neucom.2016.08.140 [120] Maghsoudi Y, Collins M, Leckie D G. Polarimetric classification of boreal forest using nonparametric feature selection and multiple classifiers. International Journal of Applied Earth Observation and Geoinformation, 2012, 19: 139-150 doi: 10.1016/j.jag.2012.04.015 [121] Ma X S, Shen H F, Yang J, Zhang L P, Li P X. Polarimetric-spatial classification of SAR images based on the fusion of multiple classifiers. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014, 7(3): 961-971 doi: 10.1109/JSTARS.2013.2265331 [122] Kiranyaz S, Ince T, Uhlmann S, Gabbouj M. Collective network of binary classifier framework for polarimetric SAR image classification: an evolutionary approach. IEEE Transactions on Systems, Man, and Cybernetics, Part B (Cybernetics), 2012, 42(4): 1169-1186 doi: 10.1109/TSMCB.2012.2187891 [123] van Zyl J J. Unsupervised classification of scattering behavior using radar polarimetry data. IEEE Transactions on Geoscience and Remote Sensing, 1989, 27(1): 36-45 doi: 10.1109/36.20273 [124] Cloude S R, Pottier E. An entropy based classification scheme for land applications of polarimetric SAR. IEEE Transactions on Geoscience and Remote Sensing, 1997, 35(1): 68-78 doi: 10.1109/36.551935 [125] Xu F, Jin Y -Q. Deorientation theory of polarimetric scattering targets and application to terrain surface classification. IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(10): 2351-2364 doi: 10.1109/TGRS.2005.855064 [126] Chen Q, Kuang G Y, Li J, Sui L C, Li D G. Unsupervised land cover/land use classification using PolSAR imagery based on scattering similarity. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(3): 1817-1825 doi: 10.1109/TGRS.2012.2205389 [127] Xiang D L, Tang T, Ban Y F, Su Y, Kuang G Y. Unsupervised polarimetric SAR urban area classification based on model-based decomposition with cross scattering. ISPRS Journal of Photogrammetry and Remote Sensing, 2016, 116: 86-100 doi: 10.1016/j.isprsjprs.2016.03.009 [128] Doulgeris A P, Anfinsen S N, Eltoft T. Automated non-Gaussian clustering of polarimetric synthetic aperture radar images. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(10): 3665-3676 doi: 10.1109/TGRS.2011.2140120 [129] Bombrun L, Vasile G, Gay M, Totir F. Hierarchical segmentation of polarimetric SAR images using heterogeneous clutter models. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(2): 726-737 doi: 10.1109/TGRS.2010.2060730 [130] Frery A C, Jacobo-Berlles J, Gambini J, Mejail M E. Polarimetric SAR image segmentation with B-splines and a new statistical model. Multidimensional Systems and Signal Processing, 2010, 21(4): 319-342 doi: 10.1007/s11045-010-0113-4 [131] Silva W B, Freitas C C, Sant$'$Anna S J, Frery A C. Classification of segments in PolSAR imagery by minimum stochastic distances between Wishart distributions. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2013, 6(3): 1263-1273 doi: 10.1109/JSTARS.2013.2248132 [132] Lee J S, Grunes M R, Ainsworth T L, Du L J, Schuler D L, Cloude S R. Unsupervised classification using polarimetric decomposition and the complex Wishart classifier. IEEE Transactions on Geoscience and Remote Sensing, 1999, 37(5): 2249-2258 doi: 10.1109/36.789621 [133] Ferro-Famil L, Pottier E, Lee J -S. Unsupervised classification of multifrequency and fully polarimetric SAR images based on the H/A/AlphaK-Wishart classifier. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(11): 2332-2342 doi: 10.1109/36.964969 [134] Park S E, Moon W M. Unsupervised classification of scattering mechanisms in polarimetric SAR data using fuzzy logic in entropy and alpha plane. IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(8): 2652-2664 doi: 10.1109/TGRS.2007.897691 [135] Lee J S, Grunes M R, Pottier E, Ferro-Famil L. Unsupervised terrain classification preserving polarimetric scattering characteristics. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(4): 722-731 doi: 10.1109/TGRS.2003.819883 [136] Hou B, Chen C, Liu X J, Jiao L C. Multilevel distribution coding model-based dictionary learning for PolSAR image classification. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(11): 5262-5280 doi: 10.1109/JSTARS.2015.2460998 [137] Cao F, Hong W, Wu Y R, Pottier E. An unsupervised segmentation with an adaptive number of clusters using the SPAN/H/α K/A space and the complex Wishart clustering for fully polarimetric SAR data analysis. IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(11): 3454-3467 doi: 10.1109/TGRS.2007.907601 [138] Dabboor M, Collins M J, Karathanassi V, Braun A. An unsupervised classification approach for polarimetric SAR data based on the chernoff distance for complex wishart distribution. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(7): 4200-4213 doi: 10.1109/TGRS.2012.2227755 [139] Kersten P R, Lee J S, Ainsworth T L. Unsupervised classification of polarimetric synthetic aperture radar images using fuzzy clustering and EM clustering. IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(3): 519-527 doi: 10.1109/TGRS.2004.842108 [140] Wang Y H, Han C Z, Tupin F. PolSAR data segmentation by combining tensor space cluster analysis and Markovian framework. IEEE Geoscience and Remote Sensing Letters, 2010, 7(1): 210-214 doi: 10.1109/LGRS.2009.2031660 [141] Yu P, Qin A K, Clausi D A. Unsupervised polarimetric SAR image segmentation and classification using region growing with edge penalty. IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(4): 1302-1317 doi: 10.1109/TGRS.2011.2164085 [142] Doulgeris A P. An automatic U-distribution and Markov random field segmentation algorithm for PolSAR images. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(4): 1819-1827 doi: 10.1109/TGRS.2014.2349575 [143] Bi H X, Sun J, Xu Z B. Unsupervised PolSAR image classification using discriminative clustering. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(6): 3531- 3544 doi: 10.1109/TGRS.2017.2675906 [144] Song W Y, Li M, Zhang P, Wu Y, Jia L, An L. Unsupervised PolSAR image classification and segmentation using Dirichlet process mixture model and Markov random fields with similarity measure. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 10(8): 3556-3568 doi: 10.1109/JSTARS.2017.2684301 [145] Zhong N, Yang W, Cherian A, Yang X L, Xia G S, Liao M S. Unsupervised classification of polarimetric SAR images via Riemannian sparse coding. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(9): 5381-5390 doi: 10.1109/TGRS.2017.2707243 [146] Kim H, Hirose A. Unsupervised fine land classification using quaternion autoencoder-based polarization feature extraction and self-organizing mapping. IEEE Transactions on Geoscience and Remote Sensing, 2018, 56(3): 1839- 1851 doi: 10.1109/TGRS.2017.2768619 [147] Liu C, Liao W Z, Li H C, Fu K, Philips W. Unsupervised classification of multilook polarimetric SAR data using spatially variant wishart mixture model with double constraints. IEEE Transactions on Geoscience and Remote Sensing, 2018, 56(10): 5600-5613 https://ieeexplore.ieee.org/document/8345175 [148] Hänsch R, Hellwich O. Semi-supervised learning for classification of polarimetric SAR-data. In: Proceedings of the 2009 IEEE International Geoscience and Remote Sensing Symposium. Cape Town, South Africa: IEEE, 2009. Ⅲ-987-Ⅲ-990 https://ieeexplore.ieee.org/document/5417941 [149] Niu X, Ban Y F. An adaptive contextual SEM algorithm for urban land cover mapping using multitemporal high-resolution polarimetric SAR data. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2012, 5(4): 1129-1139 doi: 10.1109/JSTARS.2012.2201448 [150] Uhlmann S, Kiranyaz S, Gabbouj M. Semi-supervised learning for Ill-posed polarimetric SAR classification. Remote Sensing, 2014, 6(6): 4801-4830 doi: 10.3390/rs6064801 [151] Wei B H, Yu J, Wang C, Wu H Y, Li J. PolSAR image classification using a semi-supervised classifier based on hypergraph learning. Remote Sensing Letters, 2014, 5(4): 386-395 doi: 10.1080/2150704X.2014.912765 [152] Liu H Y, Zhu D X, Yang S Y, Hou B, Gou S P, Xiong T, Jiao L, et al. Semisupervised feature extraction with neighborhood constraints for polarimetric SAR classification. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(7): 3001-3015 doi: 10.1109/JSTARS.2016.2532922 [153] Hua W Q, Wang S, Liu H Y, Liu K, Guo Y H, Jiao L C. Semisupervised PolSAR image classification based on improved cotraining. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 10(11): 4971-4986 doi: 10.1109/JSTARS.2017.2728067 [154] Hou B, Wu Q, Wen Z D, Jiao L C. Robust semisupervised classification for PolSAR image with noisy labels. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(11): 6440-6455 doi: 10.1109/TGRS.2017.2728186 [155] Zhang D, Jiao L C, Bai X, Wang S, Hou B. A robust semi-supervised SVM via ensemble learning. Applied Soft Computing, 2018, 65: 632-643 doi: 10.1016/j.asoc.2018.01.038 [156] Liu H Y, Yang S Y, Gou S P, Liu S, Jiao L C. Terrain classification based on spatial multi-attribute graph using polarimetric SAR data. Applied Soft Computing, 2018, 68: 24-38 doi: 10.1016/j.asoc.2018.03.029 [157] Geng J, Ma X R, Fan J C, Wang H Y. Semisupervised classification of polarimetric SAR image via superpixel restrained deep neural network. IEEE Geoscience and Remote Sensing Letters, 2018, 15(1): 122-126 doi: 10.1109/LGRS.2017.2777450 [158] Lecun Y, Bottou L, Bengio Y, Haffner P. Gradient-based learning applied to document recognition. Proceedings of the IEEE, 1998, 86(11): 2278-2324 doi: 10.1109/5.726791 [159] Hinton G E, Osindero S, Teh Y W. A fast learning algorithm for deep belief nets. Neural Computation, 2006, 18(7): 1527-1554 doi: 10.1162/neco.2006.18.7.1527 [160] Bengio Y. Learning deep architectures for AI. Foundations and Trends in Machine Learning, 2009, 2(1): 1-127 http://d.old.wanfangdata.com.cn/OAPaper/oai_arXiv.org_1206.5538 [161] Rasmus A, Valpola H, Honkala M, Berglund M, Honkala M, Valpola H, Raiko T. Semi-supervised learning with ladder networks. In: Proceedings of the 28th International Conference on Neural Information Processing Systems Advances in Neural Information Processing Systems. Montreal, Canada: ACM, 2015. 3546-3554 https://arxiv.org/abs/1507.02672 [162] Jiao L C, Liu F. Wishart deep stacking network for fast PolSAR image classification. IEEE Transactions on Image Processing, 2016, 25(7): 3273-3286 doi: 10.1109/TIP.2016.2567069 [163] Zhang L, Ma W P, Zhang D. Stacked sparse autoencoder in PolSAR data classification using local spatial information. IEEE Geoscience and Remote Sensing Letters, 2016, 13(9): 1359-1363 doi: 10.1109/LGRS.2016.2586109 [164] Xie W, Jiao L C, Hou B, Ma W P, Zhao J, Zhang S Y, Liu F, et al. PolSAR image classification via Wishart-AE model or Wishart-CAE model. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 10(8): 3604-3615 doi: 10.1109/JSTARS.2017.2698076 [165] Liu H Y, Yang S Y, Gou S P, Zhu D X, Wang R F, Jiao L C. Polarimetric SAR feature extraction with neighborhood preservation-based deep learning. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 10(4): 1456-1466 doi: 10.1109/JSTARS.2016.2618891 期刊类型引用(7)
1. 韦雨岑,叶子毅,庾露. 基于神经网络的SAR图像超分辨率重建技术研究. 广西水利水电. 2024(02): 1-7+14 . 百度学术
2. 谢雯,滑文强,焦李成,王若男. 采用深度学习的极化SAR地物分类方法综述. 西安电子科技大学学报. 2023(03): 151-170 . 百度学术
3. 韩宾宾,韩萍,程争. 面向对象的PolSAR图像SVM分类. 中国民航大学学报. 2022(01): 21-26 . 百度学术
4. 崔浩贵,晏庆,张晓,杨飞,范路芳. 一种PolSAR图像G0分布参数估计新方法. 无线电工程. 2021(03): 179-184 . 百度学术
5. 邓少平,孙盛. 极化SAR影像深度学习分类方法对比研究. 测绘科学. 2021(07): 120-127 . 百度学术
6. 张继超,邹勇,宋伟东,张永红,李建飞. 联合对称不确定性ReliefF算法的PolSAR影像分类. 遥感信息. 2021(04): 20-27 . 百度学术
7. 鲁力,王洁,韩要昌,吴亚晖. 改进的DCGAN在SAR图像分类中的应用. 控制工程. 2020(03): 561-566 . 百度学术
其他类型引用(6)
-