基于原型理论的极化SAR图像特征表达

黄晓菁　杨祥立　黄平平　杨　文*

①(武汉大学电子信息学院　武汉　430072)

②(内蒙古工业大学雷达技术研究所　呼和浩特　010051)

基于原型理论的极化SAR图像特征表达

黄晓菁①杨祥立①黄平平②杨文*①

①(武汉大学电子信息学院武汉430072)

②(内蒙古工业大学雷达技术研究所呼和浩特010051)

该文提出一种基于原型理论的极化SAR图像表达方法。该方法首先利用原型理论构建原型集，然后以正则化逻辑回归函数计算测试样本与每个原型集的相似度，最后通过集成投影获得图像的特征表达。在极化SAR数据上的非监督分类实验结果表明，该方法能够准确表达图像中各类地物的极化特性，达到较好的分类效果。

极化合成孔径雷达；特征表达；原型理论；非监督分类

引用格式：黄晓菁, 杨祥立, 黄平平, 等.基于原型理论的极化SAR图像特征表达[J].雷达学报, 2016, 5(2): 208–216.DOI: 10.12000/JR15071.

Reference format: Huang Xiaojing, Yang Xiangli, Huang Pingping, et al..Prototype theory based feature representation for PolSAR images[J].Journal of Radars, 2016, 5(2): 208–216.DOI: 10.12000/JR15071.

1　引言

全极化合成孔径雷达(Polarimetric Synthetic Aperture Radar, PolSAR)数据可提取地物目标的完整散射矩阵、几何结构的细节特性和地物目标的介电常数信息，极大地增强了成像雷达对目标信息的获取能力[1]。极化SAR遥感为对地观测提供了更丰富的场景和目标信息(如目标的几何形状等)，但同时也对极化图像的解译提出了新的挑战。其中，如何从极化SAR数据中有效地提取表达地物的特征信息是提高极化SAR图像解译性能的关键问题之一。

文献[2]中周晓光等人对极化SAR图像分类中应用的各种极化特征进行了详细的论述。目前，极化特征主要包括基于测量数据及其简单变换组合的特征和基于极化目标分解的特征。Cloude等人[3]将极化散射矩阵的分解方法分为3类：基于Mueller矩阵和斯托克斯矢量的方法、利用协方差或相干性矩阵的特征向量分析方法和采用相干散射矩阵模型的分解方法，并从不同角度分析了地物的散射特性。Freeman等人[4]提出三分量散射分解模型来刻画极化SAR数据中地物的散射机制，并取得了较好的建模效果。继Freeman之后，Yamaguchi等人[5]将螺旋散射功率组合到Freeman的散射机制中，增强了对密集建筑区的描述准确率。自此，通过极化分解的方式获取极化SAR数据中地物的极化特征在极化SAR图像处理中被广泛应用[6–9]。但是，对于复杂场景的分析，仅靠几种散射模型难以准确提取极化SAR数据中的地物特征信息。

近年来，原型理论(Prototype Theory)被广泛应用于图像分类领域[10,11]。K-prototype算法在数据聚类分析中也取得了很好的实验效果[12,13]。Crammer等人[14]的研究理论表明基于原型理论的分类器比最近邻分类器更准确。其中，最主要的原因是这种代表数据间本质差异的样本数据集––原型集，以图像子集间相互映射的方式表现了各类数据之间的差异，减少了传统分类器中易出现的模型过拟合情况。Martin等人[15]利用一定数量具有判别力的原型集来描述数据集，以提高分类表达与处理的效率。Dai等人[16]利用最大-最小采样方法构造原型集，并结合集成投影方法构建特征表达应用于半监督的图像分类，获得了满意的分类结果。另外，由于利用原型树或稀疏原型集能够准确学习并描述目标特征，达到快速匹配识别的效果，原型理论也被成功应用于目标跟踪与动作识别[17,18]等领域。

本文针对极化SAR图像解译中的特征描述问题，尝试将原型理论引入极化SAR图像处理中，提出一种新的极化特征提取算法。不同于光学图像中原型理论的应用，极化SAR数据表现的极化特性增加了原型集构建的难度。特别是，在衡量图像区域相似度方面，适用光学图像的欧氏距离或马氏距离在极化SAR数据上难以取得令人满意的结果。本文算法在考虑极化SAR数据散射特性的同时，利用图像区域间基于协方差矩阵数据的相似度计算构建多组原型集，并以原型理论构建图像中匀质区域之间的映射关系。最后，用主成分分析方法进行特征降维，并通过非监督分类展示最终特征描述子的应用效果。

2　极化散射特征

极化SAR图像的解译通常依赖于极化数据中所含的丰富信息。极化SAR测量数据利用极化散射矩阵S表示，

它包含了目标的散射相位特性与极化特性，全面描述了目标的电磁散射特性。其协方差矩阵C完整描述了两个正交极化通道之间的关系：

其元素包含了雷达测量得到的全部极化信息，是进行极化SAR数据处理的基础。

表1　基本极化特征表达Tab.1　Fundamental polarimetric feature representation

以散射矩阵为主，基于SAR测量数据及其数学变换的方式可作为基本的极化特征提取方式。这些特征包括后向散射强度与强度比、相位差以及去极化比等[19](参见表1)。每种对极化SAR测量数据的数学处理都能够表达出地物固有的极化信息。但是在某些情况下，同类地物可能表现出不同的极化信息，不同类别地物也可能表现出相似的极化信息，即同物异谱或异物同谱。这使得纯粹依赖于对极化SAR测量数据进行数学处理的极化特征描述方法出现描述不准确或特征不具辨别力的情况。

基于极化目标分解的散射特征提取方法更清楚地表达了极化SAR数据中目标的散射机制。以目标散射机制的稳定性为依据，可将极化目标分解方法分为两种类型：相干分解与非相干分解。前者适用于具有确定的或稳态的散射特性的目标，散射回波是相干的，其针对Sinclair矩阵的分解；后者适用于散射特性不确定、需多视处理的目标，散射回波是部分或非相干的，其针对Mueller矩阵、相干矩阵和协方差矩阵的分解。Pauli分解[20]是一种经典的相干极化目标分解方法，它将目标的极化散射矩阵分解成3个基本矩阵的加权组合。而非相干极化目标分解中基本的分解方法是Cloude分解[3]，它以特征值分解方式将相干矩阵分解为3个独立矩阵。另一种常用的非相干分解方法是Freeman-Durden分解[4]，以3种散射类型描述了目标的散射机制(具体参数参见表2)。基于目标分解的极化特征各有侧重，但仍不足以全面描述整个雷达检测场景中的地物。因此，发展新的极化特征表达方法是极化SAR图像的解译的研究重点。

表2　基于目标分解的极化特征表达Tab.2　Target decomposition based polarimetric feature representation

3　基于原型理论的特征编码

在极化SAR图像的特征分析中，以匀质的图像区域作为分析对象，不仅可避免以像素点为单元分析导致的特征表达模糊，而且能提高地物极化特征表达的准确度与效率。文中利用简单线性迭代聚类算法[21](SLIC)将极化SAR图像分割为多个匀质区域，并借鉴文献[16]中的基本思想，以基于原型理论的特征编码模型对图像的原始特征进行编码，以期获取更全面、更具分辨力的图像极化特征。算法的完整流程图如图1所示。

图1　基于原型理论的极化SAR图像特征提取算法流程图Fig.1　Framework of feature extraction algorithm based on prototype theory for PolSAR image

3.1编码模型

基于原型理论的特征编码通过从图像中选取一部分具有代表性的样本集––原型集进行特征分析，从而将原始图像特征映射到另一个特征空间，得到更为准确的图像表达。这种仅针对图像中被选原型集的分析方式避免了对图像整体场景的细致分析，有效地降低了图像分析的复杂度。为提高编码后特征的辨别力与描述的全面性，选择合适的原型集是该编码模型的重点。

原型集的选取有两个约束条件：

(1)原型集内不同类别代表样本之间的差异性较大；

(2)同类别代表样本之间的差异性较小。

同时，各原型集内部及原型集之间均需保持一定的多样性。否则，仅处理几个或几组样本，即便这些样本都具有很高的代表性也不能确保其囊括图像的所有特征。虽然，这将消耗更多计算时间，但是相应得到的特征描述子质量也会有显著的提高。

在选取合适的原型集的基础上，原型理论中图像原始特征到另一个特征空间的映射方式也是影响最终特征表达效果的一个重要因素。在分类问题中，将一系列分类器单独决策结果以某种方式组合(一般为加权或非加权投票的方式)去决定新样本的分类结果比使用单个分类器进行决策得到的结果更好。这一思想可应用于图像的特征编码，即以集成多个子编码器的特征学习方式对图像原始特征进行映射以得到新的特征表达。本文正是利用这种集成映射的方式从多个不同的角度对基本极化特征进行分析，挖掘更深层的信息。

3.2基于原型理论的特征表达

为准确选取原型集，首先需要确定其中心点。中心点是原型集中各组样本集合的中心，确定了样本描述的对象。中心点的选择需满足条件(1)，使得各组样本尽量描述不同的对象或对象的不同方面。中心点的选取算法如算法1[16]所示。在超像素分割块中以随机游走的方式选取T个初始图像区域的集合S后，以匀质区域之间协方差信息的相似度为研究对象，多次迭代，选取相互之间差异性最大的一组中心点作为输出。协方差矩阵之间的相似度利用Bartlett距离公式[22]衡量。具体地，令为各匀质图像区域协方差信息的均值，，任意两个样本点的协方差信息之间的距离为：

其中q为协方差矩阵的维度，q=3。而初始类别数T的设定不依赖于图像场景中地物的类别数，可依据经验给出一个比图像中地物类别数略大的数值。本文实验中设定T=17。

算法 1：选取中心点算法 2：筛选原型集输入：N个匀质区域的协方差信息均值输入：中心点集合S与协方差信息均值b e初始化：=0，l(n)=0，n∈[1，N]初始化：，迭代计算：当迭代次数≤m时，D(n)=｛DB(n，j)，j∈｝距离计算：S=｛随机选择的T个图像区域｝，D(n，t)=minD(n)，l(n)=t，e= X i∈X j∈DB(i，j)t∈[1，T] t为标签，e>b e如果，选取原型样本：^e=eb=那么,　for t=1 to T当迭代次数>m，　st=M个最小D(n, t)对应的｛n｝=bt=[t···t···t]∈Rn输出：一组中心点集合S=｛t｝=｛t，t｝输出：一组原型集

其次，筛选原型集中的样本。原型集中同一类别样本描述了该类别的信息。为保证类别信息的准确清晰，样本筛选需要满足条件式(2)，具体如算法2所示。先通过Bartlett距离最小化划分匀质区域的类别，对于划分为同类别的样本，再选取离中心最近的M个样本作为原型集的类别样本。由于图像中每类地物所占的面积不同，用于描述它们的图像区域数量不同，因此不能保证每个类别中样本量相等。

最后，集成x在所有投影函数H上得到的映射向量｛y1··ynum··yNum｝，即可得到新的对象特征表达。特征矢量的集成方式有很多，为不丢失任一组投影结果，试验中选取的集成方式是向量串联。值得注意的是，由于最终投影向量的组成方式是串联，很可能导致极化特征描述维度过高。为准确有效地降低特征维度并尽可能保留特征信息的主要成分，在完成极化特征提取后可采用主成分分析算法进行降维。

图2　AIRSAR数据与本文特征的聚类结果Fig.2　Original AIRSAR image and clustering results of the proposed features

4　实验结果与分析

为验证本文方法的有效性，使用荷兰Flevoland地区的AIRSAR数据和德国Oberpfaffenhofen地区的ESAR数据进行极化特征分析，并与Lee等人提出的Freeman-Wishart非监督算法[24]进行对比，从定性与定量两方面评估实验结果。

4.1实验1

荷兰Flevoland地区获取的AIRSAR图像大小为1279×1024像素，包含地物种类较多。为对聚类结果进行细致分析，选取原数据中大小为340×200像素的区域进行极化特征分析(见图2(a))。

试验中，设定SLIC的最小分割块尺寸为5×5像素，得到2720个匀质区域。将这些匀质区域划分为T=17类，以d=23维原始极化特征(第2节)为信息筛选每类中距离中心点最近的M=60%的区域构成一组原型集。本组实验中，原型集组数设置为Num=70，正则化逻辑回归训练后，针对整幅图像映射得到17×70维极化特征。由于特征维度较高，利用PCA取精度为99%对上述极化特征进行降维，并利用K均值聚类算法对降维后的21维特征描述子进行聚类，聚类类别数为10类，结果见图2(b)。为说明实验中参数––原型集组数Num对本文中提出的极化特征的影响，图2(c)和图2(d)显示了其他参数(T与M)不变的情况下，原型集组数分别为Num=50与Num=100的聚类结果。经实验对比，映射得到的极化特征维数的改变不会对本文所提出的极化特征的表达特性有明显的影响。但是，当上述维度过低时，该极化特征的描述能力将会降低以致不能准确分辨地物；而当上述维度过高时，将会增加特征提取的计算量，而且特征描述过于精准也可能造成聚类结果中同质区域纯净度降低。算法中，初始类别数T和区域选取比例值M，与原型集组数Num值在对特征描述的影响上具有类似的特性，因此实验中，我们依据图像场景中地物特性以及图像处理经验选取相应参数值。

为比较基于原型理论的特征编码与PCA降维算法的特征映射效果，图3(b)给出了未经PCA降维的(17×70维)映射得到的极化特征的K均值聚类结果。对比图3(a)与图3(b)可见，除了降低特征的维度，PCA降维算法并未对映射得到的极化特征的表达能力造成明显影响，仅在极少部分区域使得表达结果更纯净(图3(a)与图2(b)相同)。

图3　AIRSAR数据与聚类结果Fig.3　Original AIRSAR image and clustering results

图4　纯净度评估子区域Fig.4　Selected regions for purity evaluation

为有效评估新的极化特征的聚类效果，本文将文中提出的极化特征聚类表达结果与Freeman-Wishart算法的非监督分类结果(聚为10类)进行对比。同时，为进一步说明基于原型理论的特征编码的有效性，定性和定量分析未编码的特征与编码后的特征在相同参数下的聚类结果。由于该地区缺乏完整的对应的地面参照资料，文中取图像中4块均质子区域(见图4)进行聚类纯净度[25]比较，

式中，Num表示图像聚类的类别数，nr表示第r类聚类的大小，为第i类数据中聚类为第r类的样本点数。纯净度常用于衡量聚类结果的准确度，纯净度值越高，聚类效果越好。两组对比实验结果显示在图3(c)、图3(d)以及表3中。对比3组聚类结果图可见，本文算法能够清晰准确并且稳定地表达图像的类别，而其他两种算法在多个匀质区域的表达中出现描述不均一的情况。从K均值聚类结果图以及纯净度分析中可知以原型理论对原始的极化特征进行编码，增强了原始极化特征的表达准确度以及类别分辨力。本文提出的极化特征表达在4个子区域上的聚类纯净度与原始特征相比均有所提高，最高提升了23%(区域4)。比较图3(a)与图3(d)中上部区域可见，编码后的特征能够准确区分两处地物而原始特征则将两个区域表达为同一类别。而Freeman-Wishart算法，由于其仅依赖于少量极化特征，影响了聚类结果的鲁棒性。这表现为不同子区域的聚类纯净度相差很大，最大差值接近50%(区域2与区域4)。对比之下，本文算法得到的聚类结果描述更加准确，纯净度整体约为90%，且在不同区域上的差异不超过30%。这些数据直观地表明了该组实验数据上本文算法的准确率与稳定性。

表3　子区域聚类结果的纯净度(%)Tab.3　Purity(%)of the experimental results for sub image regions

图5　ESAR数据与聚类结果Fig.5　Original ESAR image and clustering results

4.2实验2

实验2提取德国柏林机场的ESAR全极化数据为研究区域，图像分辨率为3 m×3 m，大小为1200×1500像素。实验截图以及对应的地面参照分别显示在图5(a)与图5(b)中，是一幅大小为500×500像素的数据，地物类型包括林地、农田、建筑区和道路等。

本组试验中，设定SLIC的最小分割块尺寸为10×10像素，得到10000个匀质区域。其他参数为：原型集有Num=70组，每组类别数T=15，特征维数d=23，每类样本量为距离中心点最近的M=60%的区域。PCA降维后K均值的聚类结果显示在图5(e)和图5(h)中。文中选取道路与农田的聚类效果进行分析。

本文分别利用3种算法对实验区域进行5类与8类两组聚类实验，结果显示在图5(c)～图5(h)中。由图5(e)和图5(h)可见，无论是聚类为5类还是8类，本文算法能够准确表达道路和农田区域。而其他两种算法中在以5类表达图像时完全无法完成道路区域聚类，甚至在8类聚类时也不能准确描述道路区域，存在较为严重的错误分类现象。利用纯净度Purity和错误聚类率Pe分别对8类聚类中农田与道路区域的聚类结果进行评估。Pe则是聚类中属于某类而实际数据中不属于该类别的像素占实际数据中属于该类的像素的比例：

其中，Ci表示聚类为第i类地物的像素，TCi表示对比地面参照聚类中正确标注的像素，Gi表示地面参照中该类地物的像素。Pe值越低，聚类结果越好。评估结果显示在表4中。

表4　农田与道路区域聚类结果的定量分析(8类)Tab.4　Quantitative evaluation of the experimental results for farmland and road regions (8 classes)

定性和定量的分析均显示，在该组数据上，虽然Freeman-Wishart算法聚类结果的纯净度较高，均在90%左右，但是它易将不属于同类的数据划分为同一类别，导致错误聚类率高达77%。这样的聚类结果将给后续的极化SAR图像解译带来困难。本文提出的聚类算法在保证纯净度的同时，以更全面准确的极化特征进行表达，显著降低了错误聚类率。在两类地物上的聚类错误率不到Freeman-Wishart算法的1/3，这在很大程度上提高了此聚类结果的解译准确度。与原始极化特征的表达效果相比，本文算法编码后的极化特征在聚类纯净度与准确度上均有明显的提高。聚类结果的视觉观感与定量分析都说明了本文基于原型理论的特征提取算法的可行性与有效性。

5　结论

本文提出一种基于原型理论的特征编码方法，利用原型理论分析图像数据点之间极化信息的差异度映射关系，构建多组原型集，并以正则化逻辑回归函数进行训练得到多组投影函数，从而获得对地物目标特征信息的更准确、全面的表达。两组极化SAR数据的聚类实验的定性与定量评估分析均表明了该方法的有效性。

[1]Cloude S.Polarisation: Applications in Remote Sensing[M].London, U.K.: Oxford University Press, 2009.

[2]周晓光, 匡纲要, 万建伟.极化SAR图像分类综述[J].信号处理, 2008, 24(5): 806–812.Zhou Xiao-guang, Kuang Gang-yao, and Wan Jian-wei.A review of polarimetric SAR image classification[J].Signal Processing, 2008, 24(5): 806–812.

[3]Cloude S R and Pottier E.A review of target decomposition theorems in radar polarimetry[J].IEEE Transactions on Geoscience and Remote Sensing, 1996, 34(2): 498–518.

[4]Freeman A and Durden S L.A three-component scattering model for polarimetric SAR data[J].IEEE Transactions on Geoscience and Remote Sensing, 1998, 36(3): 963–973.

[5]Yamaguchi Y, Moriyama T, Ishido M, et al..Fourcomponent scattering model for polarimetric SAR image decomposition[J].IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(8): 1699–1706.

[6]闫剑, 李洋, 尹嫱, 等.引入有取向二面角散射的Freeman-Durden分解[J].雷达学报, 2014, 3(5): 574–582.Yan Jian, Li Yang, Yin Qiang, et al..Freeman-Durden decomposition with oriented dihedral scattering[J].Journal of Radars, 2014, 3(5): 574–582.

[7]滑文强, 王爽, 侯彪.基于半监督学习的SVM-Wishart极化SAR 图像分类方法[J].雷达学报, 2015, 4(1): 93–98.Hua Wen-qiang, Wang Shuang, and Hou Biao.Semisupervised learning for classification of polarimetric SAR images based on SVM-Wishart[J].Journal of Radars, 2015, 4(1): 93–98.

[8]李洪忠, 程劲松, 王超, 等.基于相似性的POLSAR 占优散射归类及非监督聚类[J].电子与信息学报, 2012, 34(6): 1501–1505.Li Hong-zhong, Cheng Jin-song, Wang Chao, et al..POLSAR dominant scattering mechanism clustering and unsupervised classification based on similarity[J].Journal of Electronics & Information Technology, 2012, 34(6): 1501–1505.

[9]陈博, 王爽, 焦李成, 等.利用0–1矩阵分解集成的极化SAR图像分类[J].电子与信息学报, 2015, 37(6): 1495–1501.Chen Bo, Wang Shuang, Jiao Li-cheng, et al..Polarimetric SAR image classification via weighted ensemble based on 0–1 matrix decomposition[J].Journal of Electronics &Information Technology, 2015, 37(6): 1495–1501.

[10]Wohlhart P, Kostinger M, Donoser M, et al..Optimizing 1-Nearest prototype classifiers[C].2013 IEEE Conference on Computer Vision and Pattern Recognition, 2013: 460–467.

[11]Quattoni A, Collins M, and Darrell T.Transfer learning for image classification with sparse prototype representations[C].IEEE Conference on Computer Vision and Pattern Recognition, 2008: 1–8.

[12]王宇, 杨莉.模糊k-prototypes聚类算法的一种改进算法[J].大连理工大学学报, 2003, 43(6): 849–852.Wang Yu and Yang Li.An improved algorithm for fuzzy kprototypes algorithm[J].Journal of Dalian University of Technology, 2003, 43(6): 849–852.

[13]陈韡, 王雷, 蒋子云.基于K-prototypes的混合属性数据聚类算法[J].计算机应用, 2010, 30(8): 2003–2005.Chen Wei, Wang Lei, and Jiang Zi-yun.K-prototypes based clustering algorithm for data mixed with numeric and categorical values[J].Journal of Computer Application, 2010, 30(8): 2003–2005.

[14]Crammer K, Gilad-Bachrach R, Navot A, et al..Margin analysis of the LVQ algorithm[C].Advances in Neural Information Processing Systems, 2002: 462–469.

[15]Kostinger M, Wohlhart P, Roth P M, et al..Joint learning of discriminative prototypes and large margin nearest neighbor classifiers[C].IEEE International Conference on Computer Vision, 2013: 3112–3119.

[16]Dai D and Gool L V.Ensemble projection for semisupervised image classification[C].IEEE International Conference on Computer Vision, 2013: 2072–2079.

[17]Wang D.Online object tracking with sparse prototypes[J].IEEE Transactions on Image Processing, 2013, 22(1): 314–325.

[18]Lin Z, Jiang Z, and Davis L S.Recognizing actions by shape-motion prototype trees[C].2009 IEEE 12th International Conference on Computer Vision, 2009: 444–451.

[19]Molinier M, Laaksonen J, Rauste Y, et al..Detecting changes in polarimetric SAR data with content-based image retrieval[C].IEEE International on Geoscience and Remote Sensing Symposium, 2007: 2390–2393.

[20]Lee J S and Pottier E.Polarimetric Radar Imaging: From Basics to Applications[M].Boca Raton, FL: CRC Press, 2009.

[21]Achanta R, Shaji A, Smith K, et al..SLIC superpixels compared to state-of-the-art superpixel methods[J].IEEE Transactions on Pattern Analysis and Machine Intelligence, 2012, 34(11): 2274–2282.

[22]Kersten P R, Lee J S, and Ainsworth T L.Unsupervised classification of polarimetric synthetic aperture radar images using fuzzy clustering and EM clustering[J].IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(3): 519–527.

[23]Collins M, Schapire R E, and Singer Y.Logistic regression, adaboost and bregman distances[J].Machine Learning, 2002, 48(1): 253–285.

[24]Lee J S, Grunes M R, Pottier E, et al..Unsupervised terrain classification preserving polarimetric scattering characteristics[J].IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(4): 722–731.

[25]Zhao Y and Karypis G.Criterion functions for document clustering: experiments and analysis[R].Technical Report, 2001, 1: 40.

黄晓菁(1991–)，女，2013年获得武汉大学工学学士学位，现于武汉大学电子信息学院信号处理实验室攻读硕士学位，主要研究方向为合成孔径雷达图像解译。

E-mail: xiaojingH@whu.edu.cn

杨祥立(1991–)，男，2014年获得中南民族大学工学学士学位，现于武汉大学电子信息学院信号处理实验室攻读硕士学位，主要研究方向为合成孔径雷达图像解译。

E-mail: xiangliyang@whu.edu.cn

黄平平(1978–)，男，2007 年获内蒙古工业大学硕士学位，2010年获中国科学院电子学研究所博士学位，现任教于内蒙古工业大学信息工程学院，副教授，硕士生导师，研究方向为合成孔径雷达信号处理与图像应用。

E-mail: cimhwangpp@163.com

杨文(1976–)，男，1998年获武汉测绘科技大学工学学士学位，2001年获武汉大学工学硕士学位，2004年获武汉大学工学博士学位，现任教于武汉大学电子信息学院，教授，博士生导师，主要从事信号处理、图像分析和遥感信息处理方面的研究工作。

E-mail: yangwen@whu.edu.cn

Prototype Theory Based Feature Representation for PolSAR Images

Huang Xiaojing①Yang Xiangli①Huang Pingping②Yang Wen①

①(School of Electronic Information, Wuhan University, Wuhan 430072, China)
②(Radar Research Institute, Inner Mongolia University of Technology, Hohhot 010051, China)

This study presents a new feature representation approach for Polarimetric Synthetic Aperture Radar (PolSAR)image based on prototype theory.First, multiple prototype sets are generated using prototype theory.Then, regularized logistic regression is used to predict similarities between a test sample and each prototype set.Finally, the PolSAR image feature representation is obtained by ensemble projection.Experimental results of an unsupervised classification of PolSAR images show that our method can efficiently represent polarimetric signatures of different land covers and yield satisfactory classification results.

Polarimetric Synthetic Aperture Radar (PolSAR); Feature representation; Prototype theory; Unsupervised classification

TN957.52

A

2095-283X(2016)02-0208-09

10.12000/JR15071

2015-06-04；改回日期：2015-12-22；网络出版：2016-02-23

杨文yangwen@whu.edu.cn

国家自然科学基金(61271401, 61461040)，内蒙古自治区科技计划项目(20140155, 20131108)

Foundation Items: The National Natural Science Foundation of China (61271401, 61461040), The Projects of Inner Mongolia Science & Technology Plan (20140155, 20131108)