CN108388907B - 基于多视角学习的极化sar数据分类器实时更新方法 - Google Patents

基于多视角学习的极化sar数据分类器实时更新方法 Download PDF

Info

Publication number
CN108388907B
CN108388907B CN201711482208.1A CN201711482208A CN108388907B CN 108388907 B CN108388907 B CN 108388907B CN 201711482208 A CN201711482208 A CN 201711482208A CN 108388907 B CN108388907 B CN 108388907B
Authority
CN
China
Prior art keywords
view
classifier
time
classification model
classification
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711482208.1A
Other languages
English (en)
Other versions
CN108388907A (zh
Inventor
聂祥丽
黄夏渊
丁曙光
乔红
张波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Automation of Chinese Academy of Science
Original Assignee
Institute of Automation of Chinese Academy of Science
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Automation of Chinese Academy of Science filed Critical Institute of Automation of Chinese Academy of Science
Priority to CN201711482208.1A priority Critical patent/CN108388907B/zh
Publication of CN108388907A publication Critical patent/CN108388907A/zh
Application granted granted Critical
Publication of CN108388907B publication Critical patent/CN108388907B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/217Validation; Performance evaluation; Active pattern learning techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
  • Image Analysis (AREA)

Abstract

本发明属于机器学习算法和图像处理技术领域,具体涉及一种基于多视角学习的极化SAR数据分类器实时更新方法,旨在为了解决分类器不能实时更新或只能单独更新而忽略视角间一致性、互补性以致影响分类精度提升的问题,本方法包括:S1,基于t时刻的极化SAR图像,提取样本极化特征、颜色特征、纹理特征;S2,通过在线多视角分类模型,估计样本的地物类别标签;S3,根据真实地物类别标签计算损失,当损失大于零时通过拉格朗日乘子法求解所述在线多视角分类模型闭式解的方法,对分类器进行更新;S4,获取t+1时刻的极化SAR图像后,重复S1至S3,直至全部极化SAR图像处理完毕。本发明可以实现极化SAR数据在线实时分类,在线分类错误率更低。

Description

基于多视角学习的极化SAR数据分类器实时更新方法
技术领域
本发明属于机器学习算法和图像处理技术领域,具体涉及一种基于多视角学习的极化SAR数据分类器实时更新方法。
背景技术
极化SAR是一种先进的微波遥感工具,在不同的收发极化组合下,对地物目标的散射特性进行测量,能够获取目标的介电常数、物理特性、几何形状和取向等信息。相较于单通道SAR,极化SAR能够获取更丰富的地物信息和分类特征,是微波成像发展的主要方向之一。因此,极化SAR在地球资源普查、环境灾害监测、城市规划、军事侦察等领域都有着广阔的应用前景。随着极化SAR***的发展和应用,对海量极化SAR数据的在线分类技术的研究具有重要的理论价值和应用意义。
在过去的二十年里,极化SAR分类问题得到了广泛的关注,并涌现了包括有监督、半监督和无监督的一系列分类方法,如西安电子科技大学申请的专利“基于稀疏编码和小波自编码器的极化SAR图像分类方法”(专利申请号:201610407916.8,公开号:CN106096652 A)中国人民解放军国防科学技术大学申请的专利“一种结合旋转域极化零角特征的极化SAR地物分类方法”(专利申请号:201710088598.8,公开号:CN 106909939 A)。然而,现有的这些极化SAR分类方法都是离线学习算法,它们要求数据在训练开始时都是可用的,从训练数据中学习模型,只有当训练完成时,才能做分类预测。另外,当新来样本被错分时,训练的分类器不再进行更新,或者需要在整个新的数据集上重新进行训练。因此,这类方法环境自适应性不强,重训练也会比较耗时。除此之外,现有的大多数方法都是单视角分类算法,即仅使用了一种类型的特征或将多种特征简单的串联起来组合成一个向量特征,而忽视了各视角特征数据的不同属性和之间的关系,影响了分类精度。
针对上述问题,本发明提出一种在线多视角学习方法,并用于极化SAR数据的在线分类任务中。不同于离线学习,在线学习能高效地更新分类器并且不会重新使用之前所有数据。对于机载或星载的极化SAR***,数据往往是大规模的,并按照一个连续的序列不断被获取。通过引入在线学习,***能从数据流中增量式地学习一个模型,其对新增样本能高效更新分类器,对动态环境具有较强自适应性,对大规模数据具有很好的扩展性。因此,在线分类技术的研究对极化SAR的实际应用非常重要。
近年来,有很多种在线学习方法被提出,如感知机(perceptron)算法,在线梯度下降(online gradient descent,OGD)算法和被动攻击(passive-aggressive,PA)算法等。其中,PA算法最小化新分类器和先前分类器间的距离,并同时最小化新分类器在当前样本上的损失,(参见参考文献1:K.Crammer,O.Dekel,J.Keshet,S.Shalev-Shwartz,andY.Singer,“Online passive-aggressive algorithms,”J.Mach.Learn.Res.,vol.7,no.Mar,pp.551–585,2006.),由于其较好的效果和较低的计算复杂度而被广泛应用。然而,该方法只适用于单视角分类问题。Nguyen等人提出两视角PA算法(参见参考文献2:T.Nguyen,K.Chang,and S.Hui,“Two-view online learning,”in Proc.Pacific-AsiaConf.Knowl.Discov.Data Min.Springer,2012,pp.74–85.)和自适应两视角PA算法,记为AdaPA(参见参考文献3:T.Nguyen,K.Chang,and S.Hui,“Adaptive two-view onlinelearning for math topic classification,”in Proc.Joint European Conf.Mach.Learn.Knowl.Discov.Databases.Springer,2012,pp.794–809.),它们分别存在权重组合固定以致无法权衡视角间重要性,一致性约束项会加剧视角间的不一致性以及没有参数权衡不同视角的拟合项等问题,从而其效果不够理想。针对这些问题,我们提出了在线两视角PA算法(参见参考文献4:X.Nie,S.Ding,H.Qiao,B.Zhang,and X.Y.Huang,“PolSAR dataonline classification based on multi-view learning,”in Proc.Int.Conf.ImageProcess.(ICIP).IEEE,2017.),算法效果有了显著提升。然而,上述这三种方法只能用于两视角和二分类问题,对于任意数量的多视角和多分类问题并不适用。Wu等人提出一种在线多模态距离度量学习算法用于图像提取(参见参考文献5:P.Wu,S.C.Hoi,P.Zhao,C.Miao,and Z.-Y.Liu,“Online multi-modal distance metric learning with application toimage retrieval,”IEEE Trans.Knowl.Data Eng.,vol.28,no.2,pp.454–467,2016.),其可用于任意数量的多视角问题中,然而,其每个视角的分类器是单独更新的,没有利用视角间的一致性和互补性关系,从而影响了分类精度。针对上述问题,本发明提出的在线分类方法在建模时考虑了多视角之前的关系,适用于任意数量视角的二分类和多分类问题。
发明内容
为了解决现有技术中的上述问题,即为了解决分类器不能实时更新或只能单独更新而忽略视角间一致性、互补性以致影响分类精度提升的问题,本发明提出了一种基于多视角学习的极化SAR数据分类器实时更新方法,包括以下步骤:
步骤S1,基于t时刻的极化SAR图像,提取样本极化特征
Figure BDA0001534028680000031
颜色特征
Figure BDA0001534028680000032
纹理特征
Figure BDA0001534028680000033
三视角数据;
步骤S2,基于
Figure BDA0001534028680000034
通过在线多视角分类模型,估计样本的地物类别标签
Figure BDA0001534028680000035
步骤S3,根据真实地物类别标签yt计算损失lt,通过与设定损失阈值对比判断样本是否被正确表示;如果样本被错误表示,则通过拉格朗日乘子法求解所述在线多视角分类模型闭式解的方法,对所述在线多视角分类模型中的分类器进行更新;
步骤S4,获取t+1时刻的极化SAR图像后,重复步骤S1至步骤S3,直至全部极化SAR图像处理完毕。
进一步地,所述在线多视角分类模型为二分类任务的在线多视角分类模型;
所述二分类任务的在线多视角分类模型为:
估计的样本类别标签为
Figure BDA0001534028680000041
预测函数为
Figure BDA0001534028680000042
为权重参数且满足
Figure BDA0001534028680000043
损失函数为lt=max{0,1-ytft},yt为真实的样本标签;
Figure BDA0001534028680000044
其中,
Figure BDA0001534028680000045
为t时刻分类器的权重向量,m为提取样本的视角数量,λi为不同视角距离变差的均衡参数,di为视角间的耦合参数,c为一个正的惩罚参数,ξ为松弛变量,
Figure BDA0001534028680000046
为待求的t+1时刻分类器的权重向量,ft+1为t+1时刻的预测函数,
Figure BDA0001534028680000047
为t时刻视角为i的样本。
进一步地,所述“通过拉格朗日乘子法求解所述在线多视角分类模型闭式解的方法,对所述在线多视角分类模型中的分类器进行更新”,其方法为:
Figure BDA0001534028680000048
其中,
Figure BDA0001534028680000049
Figure BDA0001534028680000051
Figure BDA0001534028680000052
Figure BDA0001534028680000053
Figure BDA0001534028680000054
进一步地,所述在线多视角分类模型的初始化分类器权重
Figure BDA0001534028680000055
为一个随机的ni维的列向量,即
Figure BDA0001534028680000056
i为提取的第i个视角样本。
进一步地,所述在线多视角分类模型为多分类任务的在线多视角分类模型;
所述多分类任务的在线多视角分类模型为:
估计的样本类别标签为
Figure BDA0001534028680000057
预测函数为
Figure BDA0001534028680000058
为权重参数且满足
Figure BDA0001534028680000059
损失函数为
Figure BDA00015340286800000510
yt为真实的样本标签,
Figure BDA00015340286800000511
Figure BDA00015340286800000512
其中,
Figure BDA00015340286800000513
为t时刻分类器的权重矩阵,m为提取样本的视角数量,λi为不同视角距离变差的均衡参数,di为视角间的耦合参数,c为一个正的惩罚参数,ξ为松弛变量,
Figure BDA00015340286800000514
是t+1时刻分类器的权重矩阵,ft+1为t+1时刻的预测向量,
Figure BDA00015340286800000515
为第i视角在t时刻的样本;F为矩阵的Forbenius范数。
进一步地,所述“通过拉格朗日乘子法求解所述在线多视角分类模型闭式解的方法,对所述在线多视角分类模型中的分类器进行更新”,其方法为:
Figure BDA0001534028680000061
Figure BDA0001534028680000062
其中,
Figure BDA0001534028680000063
Figure BDA0001534028680000064
Figure BDA0001534028680000065
Figure BDA0001534028680000066
Figure BDA0001534028680000067
Figure BDA0001534028680000068
Figure BDA0001534028680000069
进一步地,所述在线多视角分类模型的初始化分类器权重
Figure BDA00015340286800000610
初始化为一个随机的K×ni的矩阵,
Figure BDA00015340286800000611
i为提取的第i个视角样本。
进一步地,通过交叉验证选择参数,选择所述在线多视角分类模型估计得类别标签错误率最小的一组;
通过交叉验证选择的参数包括:
不同视角距离变差的均衡参数λi,视角间的耦合参数di,惩罚参数c,权重参数r1,r2
进一步地,极化特征包括直接从获取的极化SAR数据及其变换中提取的原始特征和基于极化分解的特征;
颜色特征包括伪彩色图像元素、主导颜色权重及HSV图像和其直方图;
纹理特征包括局部二值模式直方图,灰度共生矩阵,Gabor和小波变换系数。
进一步地,输入参数的取值范围为:
不同视角距离变差的均衡参数λi包括λ1、λ2、λ3;λ1=1,λ2,λ3∈{1,1.5}:
视角间的耦合参数di包括d1、d2、d3;d1=d2=d3{0.001,0.01,0.1};
惩罚参数c∈{0.05,0.1,0.15};
权重参数r1,r2∈{0.3,0.4},且满足
Figure BDA0001534028680000071
本发明与现有技术相比较,具有以下优点:
(1)实现极化SAR数据在线分类
本发明提出的基于PA的在线学习方法更新步都有解析表达,因此能高效地更新分类器,克服了现有的极化SAR离线分类方法不对分类器做更新或需要使用全部数据重新训练分类器的问题,使得本发明能实现实时的分类,并对动态环境具有很强的自适应性以及对大规模数据有很强的扩展性。
(2)在线分类错误率更低
本发明从极化SAR数据中提取了极化、颜色和纹理特征,并将它们作为不同的视角,本发明充分利用了它们之间的一致性和互补性关系进行建模,克服了现有技术只使用某单一特征导致信息未充分开发或使用某几个特征串联成一个高维向量导致计算复杂度太高的不足,使得本发明利用了更完备的信息,从而在线分类的错误率会更低。
附图说明
图1是本发明一种实施例的
图2(a)是由ESAR获取的德国奥博珀法芬霍芬地区的单视极化SAR数据的Pauli分解伪彩色图像;
图2(b)是由ESAR获取的德国奥博珀法芬霍芬地区的与图2(a)对应的真实地物类别标记图;
图3是在线二分类完成后整体分类图的视觉对比结果;
图4是在线多分类完成后整体分类图的视觉对比结果。
具体实施方式
下面参照附图来描述本发明的优选实施方式。本领域技术人员应当理解的是,这些实施方式仅仅用于解释本发明的技术原理,并非旨在限制本发明的保护范围。
本发明的一种基于多视角学习的极化SAR数据分类器实时更新方法,包括以下步骤:
步骤S1,基于t时刻的极化SAR图像,提取样本极化特征
Figure BDA0001534028680000081
颜色特征
Figure BDA0001534028680000082
纹理特征
Figure BDA0001534028680000083
三视角数据;
步骤S2,基于
Figure BDA0001534028680000084
通过在线多视角分类模型,估计样本的地物类别标签
Figure BDA0001534028680000085
步骤S3,根据真实地物类别标签yt计算损失lt,通过与设定损失阈值对比判断样本是否被正确表示;如果样本被错误表示,则通过拉格朗日乘子法求解所述在线多视角分类模型闭式解的方法,对所述在线多视角分类模型中的分类器进行更新;
步骤S4,获取t+1时刻的极化SAR图像后,重复步骤S1至步骤S3,直至全部极化SAR图像处理完毕。
下面按照具体实施过程中从输入参数选择、参数初化到分类器更的流程对本发明技术方案进行详细说明,如图1所示,包括以下步骤:
步骤1:通过交叉验证选择参数,我们考虑的极化、颜色和纹理的特征,因此是三视角,即m=3,输入参数的具体范围如下:惩罚参数c∈{0.05,0.1,0.15},平衡参数λ1=1,λ2,λ3∈{1,1.5},耦合参数d1=d2=d3∈{0.001,0.01,0.1},权重参数r1,r2∈{0.3,0.4},且满足
Figure BDA0001534028680000091
这里给出了参数的选择范围,通过交叉验证的方式可以选出最优的参数取值,即错误率最小的一组参数。
步骤2:初始化分类器的权重,对二分类问题,
Figure BDA0001534028680000092
Figure BDA0001534028680000093
Figure BDA0001534028680000094
为一个随机的ni维的列向量;对多分类问题,
Figure BDA0001534028680000095
Figure BDA0001534028680000096
初始化为一个随机的K×ni的矩阵。
步骤3:获取极化SAR协方差数据,提取极化、颜色和纹理特征:极化特征包括直接从获取的极化SAR数据及其变换中提取的原始特征和基于极化分解的特征;颜色特征包括伪彩色图像元素,主导颜色权重及HSV图像和其直方图;纹理特征包括局部二值模式直方图,灰度共生矩阵,Gabor和小波变换系数等。本发明中具体使用的极化、颜色和纹理特征见表1,可以看出这三类特征的维数分别是n1=45,n2=34和n3=86,本发明将它们作为三个视角数据用于后续分类。
表1:用于在线分类的极化、颜色和纹理特征
Figure BDA0001534028680000097
Figure BDA0001534028680000101
步骤4:建立在线多视角分类模型,并根据分类函数预测样本的标签。
具体来讲,针对二分类任务的在线多视角分类模型如下:
Figure BDA0001534028680000102
s.t.l(w;(xt,yt))≤ξ;ξ≥0.其表示对二分类损失函数的松弛约束,其中松弛变量ξ必须是非负的。
其中,
Figure BDA0001534028680000103
为t时刻分类器的权重向量,m为提取样本的视角数量,λi是不同视角距离变差的均衡参数,di是视角间的耦合参数,c是一个正的惩罚参数,ξ为松弛变量,
Figure BDA0001534028680000104
为待求的t+1时刻分类器的权重向量,ft+1为t+1时刻的预测函数,
Figure BDA0001534028680000105
为t时刻视角为i的样本。
损失函数定义为hinge-loss损失lt=max{0,1-ytft}。预测函数定义为
Figure BDA0001534028680000106
其中ri∈(0,1)为权重参数且满足
Figure BDA0001534028680000107
则估计的样本的类别标签
Figure BDA0001534028680000108
针对多分类任务的在线多视角分类模型如下:
Figure BDA0001534028680000109
s.t.lMC(W;(xt,yt))≤ξ;ξ≥0。表示对多类分类损失函数的松弛约束,其中松弛变量ξ必须大于或等于0。
其中,λi,di,c是正参数,
Figure BDA0001534028680000111
为t时刻分类器的权重矩阵,
Figure BDA0001534028680000112
是t+1时刻分类器的权重矩阵,ft+1为t+1时刻的预测向量,F为矩阵的Forbenius范数。
预测函数
Figure BDA0001534028680000113
可知
Figure BDA0001534028680000114
lMC表示多类分类的损失函数,其定义为
Figure BDA0001534028680000115
yt为真实的样本标签,
Figure BDA0001534028680000116
多数情况下,估计的样本的类别标签
Figure BDA0001534028680000117
步骤5:根据样本对应的真实标签,计算损失lt,如果lt=0,表示样本被正确分类,不对分类器更新;如果lt>0,表示样本被错误分类,需要对当前分类器做更新。
针对二分类问题,损失lt=max{0,1-ytft},如果lt>0,对上一步骤中二分类优化问题,通过拉格朗日乘子法求解,可得到如下的闭式解,以此对分类器进行更新:
Figure BDA0001534028680000118
其中,
Figure BDA0001534028680000119
Figure BDA00015340286800001110
Figure BDA00015340286800001111
Figure BDA00015340286800001112
Figure BDA00015340286800001113
针对多分类问题,损失
Figure BDA0001534028680000121
如果lt>0,对上一步骤中多分类优化问题,通过拉格朗日乘子法求解,可得到如下的闭式解,以此对分类器进行更新:
Figure BDA0001534028680000122
Figure BDA0001534028680000123
其中,
Figure BDA0001534028680000124
Figure BDA0001534028680000125
Figure BDA0001534028680000126
Figure BDA0001534028680000127
Figure BDA0001534028680000128
Figure BDA0001534028680000129
Figure BDA00015340286800001210
步骤6:若还有新样本输入,返回到步骤1;如果所有样本都处理完毕,则计算整个在线学习过程的分类错误率,针对每个子类的错误率是子类中被错分的样本数比子类总样本数,总错误率是被错分的样本数比总样本数;画出最终的分类图,极化SAR图像中具有相同地物类别的像素,用同一种颜色表示,进而得到分类图。
下面分别结合二分类任务和多分类任务情况下的细分步骤进行进一步说明。
本发明一种实施例中,针对二分类任务的在线多视角学习算法包括以下步骤:
(1)通过交叉验证选择参数:惩罚参数c>0,平衡参数λi>0,耦合参数di>0和权重参数ri∈(0,1)且
Figure BDA00015340286800001211
(2)初始化分类器的权重向量
Figure BDA0001534028680000131
Figure BDA0001534028680000132
其中,m是视角的个数,ni是第i个视角的维数,可知
Figure BDA0001534028680000133
(3)在t时刻,接收极化SAR数据,提取极化
Figure BDA0001534028680000134
颜色
Figure BDA0001534028680000135
和纹理特征
Figure BDA0001534028680000136
分别作为不同视角,因此,接收的样本
Figure BDA0001534028680000137
Figure BDA0001534028680000138
(4)计算预测函数
Figure BDA0001534028680000139
(5)估计样本的类别标签
Figure BDA00015340286800001310
(6)接收正确的标签yt∈{+1,-1};
(7)计算损失lt=max{0,1-ytft};
(8)如果lt=0,表示样本被正确分类,不对分类器更新,直接进入下一轮迭代;如果lt>0,就对当前分类器做如下更新:
(9)根据下面公式计算辅助变量:
Figure BDA00015340286800001311
Figure BDA00015340286800001312
Figure BDA00015340286800001313
Figure BDA00015340286800001314
Figure BDA00015340286800001315
(10)更新分类器:
Figure BDA00015340286800001316
(11)若有新样本到来,t=t+1,返回执行第(3)步,否则算法终止。
本发明一种实施例中针对多分类任务的在线多视角学习算法包括以下步骤:
(1)通过交叉验证选择参数:c>0,λi>0,di>0和ri∈(0,1)且
Figure BDA00015340286800001317
其中m是视角的个数;
(2)初始化分类器的权重向量
Figure BDA0001534028680000141
Figure BDA0001534028680000142
其中K是类别数,ni是第i个视角的维数,可知
Figure BDA0001534028680000143
(3)在t时刻,接收极化SAR数据,提取极化
Figure BDA0001534028680000144
颜色
Figure BDA0001534028680000145
和纹理特征
Figure BDA0001534028680000146
分别作为不同视角,因此,接收的样本
Figure BDA0001534028680000147
Figure BDA0001534028680000148
(4)计算预测函数
Figure BDA0001534028680000149
可知
Figure BDA00015340286800001410
(5)估计样本的类别标签
Figure BDA00015340286800001411
(6)接收正确的标签yt∈Y=[1,2,...,K};
(7)计算
Figure BDA00015340286800001412
(8)计算损失
Figure BDA00015340286800001413
(9)如果lt=0,表示样本被正确分类,不对分类器更新,直接进入下一轮迭代,否则如果lt>0,就对当前分类器做如下更新:
(10)根据下面公式计算辅助变量:
Figure BDA00015340286800001414
Figure BDA00015340286800001415
Figure BDA00015340286800001416
Figure BDA00015340286800001417
Figure BDA00015340286800001418
Figure BDA00015340286800001419
Figure BDA00015340286800001420
(11)更新分类器:
Figure BDA00015340286800001421
Figure BDA00015340286800001422
(12)若有新样本到来,t=t+1,返回执行第(3)步,否则算法终止。
下面结合图2-4对本发明的效果做进一步的说明:
实验数据及条件:
本发明使用真实的极化SAR数据做测试实验,是由E-SAR传感器获取的德国奥博珀法芬霍芬地区的单视L波段数据,可从欧洲空间局网站下载到。该数据的Pauli分解伪彩色图像如图2(a)所示,其尺寸为1300*1200,图2(b)是其对应的真实地物分类图,其中的色块分别表示城区、林地、公路、农田、其它场景。
仿真实验中,使用的软件:MATLAB R2015b,处理器:Intel(R)Core(TM)i7-6700HQ,内存:20.0GB,操作***:64位Windows10。
实验内容及结果分析:
实验中分别考虑了二分类和多分类任务,它们分别对应城市场景提取和地物分类问题。为更好评估本发明提出的方法的效果,与PA分别在极化、颜色、纹理和它们的组合特征上的结果对比,分别记为PA_Pol,PA_Col,PA_Tex,PA_Cat。另外,也与参考文献3中提出的AdaPA方法和参考文献5中提出的OMDML方法的结果进行了比较,注意AdaPA方法只适用于二分类问题,所以在多分类任务中未与其比较。为了更好的对比这些方法,它们包含的参数通过交叉验证来选取,参数的选取范围设置如下:PA方法的攻击参数c∈[0.05,0.15];AdaPA方法的耦合参数d∈{0.001,0.01,0.1},权重参数r∈(0,1),惩罚参数c∈[0.01,0.15];OMDML方法的惩罚参数c∈[0.01,0.15],折扣参数β∈[0.8,1];本发明的方法,λ1=1,λ2,λ3∈{1,1.5},d1=d2=d3∈{0.001,0.01,0.1},r1,r2∈{0.3,0.4},c∈{0.05,0.1,0.15}。
图3展示了在线二分类完成后整体分类图的视觉对比结果,图3(a)是城区场景提取的真实类别标记图,非城区区域被标记为白色;图3(b)-(h)分别是PA_Pol,PA_Col,PA_Tex,PA_Cat,AdaPA,OMDML和本发明的分类结果图。表2给出了二分类情况下这些方法的分类错误率对比结果。从图3中可以看出,PA_Cat的分类结果要明显好于PA_Pol,PA_Col和PA_Tex的结果,因为PA_Cat中使用了更多的特征信息,这点也可以从表2中得到证实。另外,根据表2可知AdaPA和OMDML的结果中城市区域被错分为非城区的比率分别是43.95%和51.08%,要显著高于PA_Cat对城区的错分率29.72%,该结论也可以从图3(e)-(g)中看出,在(f)和(g)中有很多地方被错分成白色区域。从表2和图3中可得出结论,本发明提出的方法得到了最低的正样本(即城区)分类错误率22.89%和最低的总体分类错误率7.05%。
表2:二分类情况下的分类错误率对比结果
方法 PA_Pol PA_Col PA_Tex PA_Cat AdaPA OMDML 本方法
城区 0.4169 0.5449 0.4442 0.2972 0.4395 0.5108 0.2289
非城区 0.0777 0.1039 0.0847 0.0571 0.0346 0.0199 0.0396
总错误率 0.1333 0.1758 0.1435 0.0965 0.1007 0.1012 0.0705
图4是在线多分类完成后整体分类图的视觉对比结果,其对应的真实类别标签图见图2(b)。图4(a)-(f)分别是PA_Pol,PA_Col,PA_Tex,PA_Cat,OMDML和本发明的分类结果图。表3给出了该多分类情况下这些方法的分类错误率对比结果。可以看出相比于颜色和纹理特征,极化特征提供了更好的判别信息,因为PA_Pol的总体错误率要比PA_Col和PA_Tex的低10%,另外,PA_Col和PA_Tex在公路和农田区域有超过一半的像素被错分了。而PA_Cat的结果要显著优于PA_Pol,PA_Col和PA_Tex,并且PA_Cat改善了边界区域的判别性,如图4(d)所示。OMDML的总体分类精度比PA_Cat提升了6%,而且根据图4(e),不同区域的边界变得更加清晰。根据表3中的数值结果和图4中的视觉结果,可以看出本发明提出的方法使得大部分的样本被正确分类,与其它方法相比,达到了最低的总体错误率。
表3:多分类情况下的分类错误率对比
Figure BDA0001534028680000171
本领域技术人员应该能够意识到,结合本文中所公开的实施例描述的各示例的方法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明电子硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以电子硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。本领域技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
术语“包括”或者任何其它类似用语旨在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备/装置不仅包括那些要素,而且还包括没有明确列出的其它要素,或者还包括这些过程、方法、物品或者设备/装置所固有的要素。
至此,已经结合附图所示的优选实施方式描述了本发明的技术方案,但是,本领域技术人员容易理解的是,本发明的保护范围显然不局限于这些具体实施方式。在不偏离本发明的原理的前提下,本领域技术人员可以对相关技术特征做出等同的更改或替换,这些更改或替换之后的技术方案都将落入本发明的保护范围之内。

Claims (9)

1.一种基于多视角学习的极化SAR数据分类器实时更新方法,其特征在于,包括以下步骤:
步骤S1,基于t时刻的极化SAR图像,提取样本极化特征
Figure FDA0002733540270000011
颜色特征
Figure FDA0002733540270000012
纹理特征
Figure FDA0002733540270000013
三视角数据;
步骤S2,基于
Figure FDA0002733540270000014
通过在线多视角分类模型,估计样本的地物类别标签
Figure FDA0002733540270000015
步骤S3,根据真实地物类别标签yt计算损失lt,通过与设定损失阈值对比判断样本是否被正确表示;如果样本被错误表示,则通过拉格朗日乘子法求解所述在线多视角分类模型闭式解的方法,对所述在线多视角分类模型中的分类器进行更新;
步骤S4,获取t+1时刻的极化SAR图像后,重复步骤S1至步骤S3,直至全部极化SAR图像处理完毕;
其中,所述在线多视角分类模型包括二分类任务的在线多视角分类模型和多分类任务的在线多视角分类模型;
所述多分类任务的在线多视角分类模型为:
估计的样本类别标签为
Figure FDA0002733540270000016
预测函数为
Figure FDA0002733540270000017
ri∈(0,1)为权重参数且满足
Figure FDA0002733540270000018
损失函数为
Figure FDA0002733540270000019
yt为真实的样本类别标签,yt∈Y={1,2,...,K},K表示类别标签的总数,
Figure FDA00027335402700000110
Figure FDA00027335402700000111
Figure FDA0002733540270000021
Figure FDA0002733540270000022
为t时刻分类器的权重矩阵,m为提取样本的视角数量,λi为不同视角距离变差的均衡参数,di为视角间的耦合参数,c为一个正的惩罚参数,ξ为松弛变量,
Figure FDA0002733540270000023
是t+1时刻分类器的权重矩阵,ni是第i个视角的维数且i=1,...,m,ft+1为t+1时刻的预测函数,
Figure FDA0002733540270000024
为第i视角在t时刻的样本;F为矩阵的Forbenius范数。
2.根据权利要求1所述的基于多视角学习的极化SAR数据分类器实时更新方法,其特征在于,所述在线多视角分类模型为二分类任务的在线多视角分类模型;
所述二分类任务的在线多视角分类模型为:
估计的样本类别标签为
Figure FDA0002733540270000025
预测函数为
Figure FDA0002733540270000026
ri∈(0,1)为权重参数且满足
Figure FDA0002733540270000027
损失函数为lt=max{0,1-ytft},yt为真实的样本标签;
Figure FDA0002733540270000028
s.t.lt=max{0,1-ytft}≤ξ;ξ≥0
其中,
Figure FDA0002733540270000029
为t时刻分类器的权重向量,m为提取样本的视角数量,λi为不同视角距离变差的均衡参数,di为视角间的耦合参数,c为一个正的惩罚参数,ξ为松弛变量,
Figure FDA00027335402700000210
为待求的t+1时刻分类器的权重向量,ft+1为t+1时刻的预测函数,
Figure FDA00027335402700000211
为t时刻视角为i的样本。
3.根据权利要求2所述的基于多视角学习的极化SAR数据分类器实时更新方法,其特征在于,通过拉格朗日乘子法求解所述在线多视角分类模型闭式解的方法,对所述在二分类任务的线多视角分类模型中的分类器进行更新,其方法为:
Figure FDA0002733540270000031
其中,
Figure FDA0002733540270000032
Figure FDA0002733540270000033
Figure FDA0002733540270000034
Figure FDA0002733540270000035
Figure FDA0002733540270000036
4.根据权利要求3所述的基于多视角学习的极化SAR数据分类器实时更新方法,其特征在于,所述在线多视角分类模型的初始化分类器权重
Figure FDA0002733540270000037
为一个随机的ni维的列向量,
Figure FDA0002733540270000038
i为提取的第i个视角样本。
5.根据权利要求1所述的基于多视角学习的极化SAR数据分类器实时更新方法,其特征在于,通过拉格朗日乘子法求解所述在线多视角分类模型闭式解的方法,对所述在多分类任务的线多视角分类模型中的分类器进行更新,其方法为:
Figure FDA0002733540270000039
Figure FDA00027335402700000310
其中,
Figure FDA00027335402700000311
Figure FDA00027335402700000312
Figure FDA0002733540270000041
Figure FDA0002733540270000042
Figure FDA0002733540270000043
Figure FDA0002733540270000044
Figure FDA0002733540270000045
6.根据权利要求5所述的基于多视角学习的极化SAR数据分类器实时更新方法,其特征在于,所述在线多视角分类模型的初始化分类器权重
Figure FDA0002733540270000046
初始化为一个随机的K×ni的矩阵,
Figure FDA0002733540270000047
i为提取的第i个视角样本。
7.根据权利要求1-6任一项所述的基于多视角学习的极化SAR数据分类器实时更新方法,其特征在于,通过交叉验证选择参数,选择所述在线多视角分类模型估计得类别标签错误率最小的一组;
通过交叉验证选择的参数包括:
不同视角距离变差的均衡参数λi,视角间的耦合参数di,惩罚参数c,权重参数r1,r2
8.根据权利要求7所述的基于多视角学习的极化SAR数据分类器实时更新方法,其特征在于,
极化特征包括直接从获取的极化SAR数据及其变换中提取的原始特征和基于极化分解的特征;
颜色特征包括伪彩色图像元素、主导颜色权重及HSV图像和其直方图;
纹理特征包括局部二值模式直方图,灰度共生矩阵,Gabor和小波变换系数。
9.根据权利要求7所述的基于多视角学习的极化SAR数据分类器实时更新方法,其特征在于,输入参数的取值范围为:
不同视角距离变差的均衡参数λi包括λ1、λ2、λ3;λ1=1,λ2,λ3∈{1,1.5}:
视角间的耦合参数di包括d1、d2、d3;d1=d2=d3∈{0.001,0.01,0.1};
惩罚参数c∈{0.05,0.1,0.15};
权重参数r1,r2∈{0.3,0.4},且满足
Figure FDA0002733540270000051
CN201711482208.1A 2017-12-29 2017-12-29 基于多视角学习的极化sar数据分类器实时更新方法 Active CN108388907B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711482208.1A CN108388907B (zh) 2017-12-29 2017-12-29 基于多视角学习的极化sar数据分类器实时更新方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711482208.1A CN108388907B (zh) 2017-12-29 2017-12-29 基于多视角学习的极化sar数据分类器实时更新方法

Publications (2)

Publication Number Publication Date
CN108388907A CN108388907A (zh) 2018-08-10
CN108388907B true CN108388907B (zh) 2021-03-02

Family

ID=63076666

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711482208.1A Active CN108388907B (zh) 2017-12-29 2017-12-29 基于多视角学习的极化sar数据分类器实时更新方法

Country Status (1)

Country Link
CN (1) CN108388907B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110222778B (zh) * 2019-06-11 2021-04-13 中国科学院自动化研究所 基于深度森林的在线多视角分类方法、***、装置
CN110222793B (zh) * 2019-06-20 2021-06-22 中国科学院自动化研究所 基于多视角主动学习的在线半监督分类方法及***
CN110414623B (zh) * 2019-08-06 2022-06-24 广东工业大学 基于多视角学习的分类模型构建方法及装置
CN110766084B (zh) * 2019-10-28 2021-06-15 北京理工大学 基于cae与hl-cnn的小样本sar目标识别方法
CN111949794B (zh) * 2020-08-14 2024-06-28 扬州大学 一种用于文本多分类任务的在线主动机器学习方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103514456A (zh) * 2013-06-30 2014-01-15 安科智慧城市技术(中国)有限公司 基于压缩感知多核学习的图像分类方法及其装置
CN104463193A (zh) * 2014-11-04 2015-03-25 西安电子科技大学 基于深度稀疏ica的极化sar图像分类方法
CN105303198A (zh) * 2015-11-17 2016-02-03 福州大学 一种自定步长学习的遥感影像半监督分类方法
CN105512680A (zh) * 2015-12-02 2016-04-20 北京航空航天大学 一种基于深度神经网络的多视sar图像目标识别方法
CN107133649A (zh) * 2017-05-09 2017-09-05 西安电子科技大学 基于点‑域距离的增量式极化sar地物分类方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150316646A1 (en) * 2014-05-01 2015-11-05 Utah State University Synthetic aperture radar target modeling
US9846822B2 (en) * 2015-12-31 2017-12-19 Dropbox, Inc. Generating and utilizing normalized scores for classifying digital objects

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103514456A (zh) * 2013-06-30 2014-01-15 安科智慧城市技术(中国)有限公司 基于压缩感知多核学习的图像分类方法及其装置
CN104463193A (zh) * 2014-11-04 2015-03-25 西安电子科技大学 基于深度稀疏ica的极化sar图像分类方法
CN105303198A (zh) * 2015-11-17 2016-02-03 福州大学 一种自定步长学习的遥感影像半监督分类方法
CN105512680A (zh) * 2015-12-02 2016-04-20 北京航空航天大学 一种基于深度神经网络的多视sar图像目标识别方法
CN107133649A (zh) * 2017-05-09 2017-09-05 西安电子科技大学 基于点‑域距离的增量式极化sar地物分类方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
POLSAR DATA ONLINE CLASSIFICATION BASED ON MULTI-VIEW LEARNING;Xiangli Nie等;《IEEE Signal Processing Society》;20170915;第1页 *
基于极化SAR技术的舰船检测与识别算法研究;杨慧;《中国优秀硕士学位论文全文数据库》;20140915;第I136-469页 *

Also Published As

Publication number Publication date
CN108388907A (zh) 2018-08-10

Similar Documents

Publication Publication Date Title
CN108388907B (zh) 基于多视角学习的极化sar数据分类器实时更新方法
Makantasis et al. Tensor-based classification models for hyperspectral data analysis
Montazer et al. An improved radial basis function neural network for object image retrieval
Le Goff et al. Deep learning for cloud detection
CN108537102B (zh) 基于稀疏特征与条件随机场的高分辨sar图像分类方法
Ke et al. Adaptive change detection with significance test
Krylov et al. Supervised high-resolution dual-polarization SAR image classification by finite mixtures and copulas
Awad An Unsupervised Artificial Neural Network Method for Satellite Image Segmentation.
CN107491734B (zh) 基于多核融合与空间Wishart LapSVM的半监督极化SAR图像分类方法
CN109359525B (zh) 基于稀疏低秩的判别谱聚类的极化sar图像分类方法
CN103955709B (zh) 基于加权合成核与tmf的极化sar图像分类方法
CN117475236B (zh) 用于矿产资源勘探的数据处理***及其方法
Fu et al. Contextual online dictionary learning for hyperspectral image classification
Qin et al. Learning discriminative embedding for hyperspectral image clustering based on set-to-set and sample-to-sample distances
CN113723492A (zh) 一种改进主动深度学习的高光谱图像半监督分类方法及装置
CN108154511B (zh) 基于子模字典学习的sar图像分割方法
Li et al. Alter-cnn: An approach to learning from label proportions with application to ice-water classification
Alhichri et al. Clustering of hyperspectral images with an ensemble method based on fuzzy c-means and markov random fields
Du et al. Training SegNet for cropland classification of high resolution remote sensing images
CN117710745A (zh) 一种基于证据多视角核化的目标分类方法及装置
CN104050489B (zh) 一种基于多核最优化的合成孔径雷达自动目标识别方法
CN116129280B (zh) 一种遥感影像雪检测的方法
Bolton et al. Random set framework for context-based classification with hyperspectral imagery
Teke et al. Texture classification and retrieval using the random neural network model
Cao et al. Hyperspectral image classification based on three-dimensional adaptive sampling and improved iterative shrinkage-threshold algorithm

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant