CN107392140B - 基于低秩稀疏分解和增量学习的极化sar地物分类方法 - Google Patents

基于低秩稀疏分解和增量学习的极化sar地物分类方法 Download PDF

Info

Publication number
CN107392140B
CN107392140B CN201710589213.6A CN201710589213A CN107392140B CN 107392140 B CN107392140 B CN 107392140B CN 201710589213 A CN201710589213 A CN 201710589213A CN 107392140 B CN107392140 B CN 107392140B
Authority
CN
China
Prior art keywords
matrix
low
rank
sample
value
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
CN201710589213.6A
Other languages
English (en)
Other versions
CN107392140A (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.)
Suzhou Wenjie Sensing Technology Co ltd
Original Assignee
Suzhou Wenjie Sensing Technology Co ltd
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 Suzhou Wenjie Sensing Technology Co ltd filed Critical Suzhou Wenjie Sensing Technology Co ltd
Priority to CN201710589213.6A priority Critical patent/CN107392140B/zh
Publication of CN107392140A publication Critical patent/CN107392140A/zh
Application granted granted Critical
Publication of CN107392140B publication Critical patent/CN107392140B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/13Satellite images
    • 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/213Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2136Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on sparsity criteria, e.g. with an overcomplete basis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/22Matching criteria, e.g. proximity measures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2415Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate

Landscapes

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

Abstract

本发明公开了一种基于低秩稀疏分解和增量学习的极化SAR地物分类方法。其方案是:1.输入待分类的极化SAR图像极化相干矩阵进行统计特征提取;2.进行低秩稀疏分解得到用于分类的低秩特征矩阵;3.初始化训练样本矩阵;4.构造一种相似性度量公式,给新增未标记样本标上类标签,实现增量式的循坏迭代标记分类。本发明将低秩稀疏分解、增量学习与极化SAR数据的统计特征相结合,在减少噪声干扰、充分保持样本相似性和细节信息的同时,实现增量式的样本标记分类,极大的提高了分类精度。

Description

基于低秩稀疏分解和增量学习的极化SAR地物分类方法
技术领域:
本发明属于图像处理和遥感技术领域,涉及极化SAR图像的地物分类,具体是一种基于低秩稀疏分解和增量学习的极化SAR地物分类方法,可用于对具有区域一致性的极化SAR图像进行分类。
背景技术:
极化合成孔径雷达(Polarimetric Synthetic Aperture Radar,极化SAR)是一种多通道、多参数的雷达成像系统,是通过实时的测量目标在每个分辨单元内的散射回波来获取目标的多种极化信息的,如极化散射矩阵、极化相干矩阵或者Stokes矩阵等。它是微波技术领域的先进测控技术。与SAR图像相比,极化SAR图像能够提供更多的目标地物信息,在目标检测与识别、地物分类等多方面有广泛的应用。极化SAR图像分类是图像解译的重要内容,旨在根据象元的极化信息,将包含相同或相似性质的象元分为一类,不同性质的象元分为多类。由于极化SAR图像包含全方位的极化信息,在加大目标辨识度的同时也增加了极化SAR图像分类的难度,因而提高极化SAR图像分类水平具有重要的现实意义。
目前的极化SAR图像分类方法主要是利用极化信息和统计特性进行分类的。典型的基于目标分解的分类方法有:Freeman分解和H/α分解。Freeman分解是把极化信息的协方差矩阵C分解为体散射、偶次散射和表面散射三种散射机制模型。H/α分解则是在基于相干矩阵T的基础上,根据散射熵和散射角提出的分解方法。除此以外,根据是否知道数据的先验信息,分为监督分类和无监督分类。监督分类是利用已知的训练样本和真实的地物分类结果等进行分类,常用的有贝叶斯分类器、SVM等。无监督分类则是在没有先验信息的前提下,利用极化数据的各种统计特征完成的分类,常用的有复Wishart分类器。监督分类精度高,但训练过程复杂,无监督分类算法速率快,但分类精度不高。从处理的数据的层次上来说,可以分为基于像素的分类和基于区域的分类。基于像素的分类方法的边缘信息保持良好,但易受相干斑噪声的影响。基于区域的分类方法对同类区域上效果良好,但边缘信息容易丢失。
增量学习是指能够不断的从新样本中学习新的知识同时保存之前学习到的知识的学习系统。它是一个连续的、动态的学习过程。对于增量学习,RobiPolokar曾提出几个条件:
(1)增量学习算法应该能够从新样本中学习新知识;(2)增量学习算法不需要处理之前已经学习过的样本,每次只有一个训练样本被看到和学习;(3)增量学习算法应该充分利用历史的训练结果。增量学习降低了时间和空间的复杂度,能够有效的解决图像分类中训练复杂度高、算法效率低等问题,是近年来的研究热点。
发明内容:
本发明的目的在于解决传统的基于像素的分类方法中易受相干斑噪声影响,细节保持不理想的缺陷,提出了一种基于低秩稀疏分解和增量学习的极化SAR地物分类方法,提高分类精度。
本发明的技术方案是这样实现的:
一.技术思路
本发明旨在解决目前极化SAR分类方法中对边缘信息保持不理想、分类精度不高的学习问题,将增量思想运用到极化SAR地物分类中,利用构造一种相似度权值因子W对新增的未标记样本逐个标上类标签,实现增量式的极化SAR特征数据标记分类,从而大幅度的提高分类精度。
二.实现方案
根据上述思路,本发明的实现步骤包括如下:
(1)输入待分类的极化SAR图像的极化相干矩阵T,并对极化相干矩阵T进行统计特征的提取,得到相应的特征数据矩阵;
(2)对特征数据矩阵进行低秩稀疏分解,得到低秩特征矩阵X;
(3)选取训练样本占标记样本比例为5%的训练样本矩阵Y;
(4)构造一种相似度权值度量因子W,对于新增的未标记样本,给其标上类标签,实现增量式的极化SAR特征数据标记分类;
(4a)计算第i行j列未标记样本xij与训练样本yh的相似度权值度量因子Wh,得到相似度距离矩阵W:
W=[W1,W2,...,Wh,...,WN]
其中N表示训练矩阵Y中的样本总数;
(4b)在相似度权值矩阵W上求得相似度权值因子最大值Wmax,并得到xij取得相似度权值因子最大值Wmax时相对应的训练样本的样本标签值lz
(4c)引入条件阈值τ,τ∈[0.970,0.995]。τ的初始值取0.980,变化步长为0.005。每遍历一遍未标记样本,τ的值振荡式增加或减少0.005;同时设定终止阈值B为图像数据大小的千分之一;
(4d)将相似度权值因子最大值Wmax与条件阈值τ比较,若Wmax>τ,则另lij=lz,其中lij表示未标记样本xij的样本标签值,同时将样本xij作为第N+1个训练样本放入到训练样本矩阵Y中,N=N+1。
(4e)计算未标记样本数目m=M-N,M表示原始未标记样本数目与训练样本数目之和,将未标记样本数目m与迭代终止阈值B做比较:
若m<B,则对预测标签矩阵L1进行局部邻域的微调,得到最终的分类结果;
若m≥B,则返回步骤(4a)。
本发明利用低秩稀疏分解对极化SAR特征数据进行处理,不但很好的去除了图像中的噪声,而且在保证稀疏性约束的同时保证了极化SAR特征数据的空间约束性。本发明将增量思想运用到极化SAR图像分类中,充分利用了历史的标记结果,对数据实现了增量式的循环迭代标记分类,实现了数据的可持续性学习,提高了算法的分类精度。
附图说明:
图1是本发明的实现流程图。
图2是本发明仿真实验采用的荷兰Flevoland1地区极化SAR图像的地物真实图;
图3是本发明和对比方法H/α分解及层次聚类方法和K最近邻分类方法对图2中的Flevoland1图像的分类结果显示图。
具体实施方式:
下面结合附图和所采用的极化SAR图像对本发明进行详细的描述。
参照图1,本发明的实现步骤如下:
步骤1.输入待分类的极化SAR图像的极化相干矩阵T,并对极化相干矩阵T进行统计特征的提取,得到相应的特征数据矩阵。
本步骤的具体实现如下:
(1a)极化相干矩阵T是一个复矩阵,其每一个样本点表示为:
Figure BDA0001354317720000031
其中Ti,j表示极化相干矩阵T中第i行j列的样本,
Figure BDA0001354317720000032
Figure BDA0001354317720000033
分别表示t21、t31和t32的共轭复数;
(1b)利用极化SAR数据的相干矩阵T的特征值分解得到每个像素点的散射功率熵H、散射角α和各向异性度A:
H=-(P1log3P1+P2log3P2+P3log3P3)
Figure BDA0001354317720000041
其中,Pt表示第t种散射机制发生的概率,Pt∈[0,1],αt表示第t种散射机制的类型,αi=0表示目标散射体发生表面散射,αi=π/4表示目标散射体发生偶极子散射(体散射),αi=π/2表示目标散射体发生二面角散射(偶次散射),λt表示极化相干矩阵T的第t个特征值,t∈[1,3];
(1c)利用Freeman分解得到每个像素点的三种散射矩阵:体散射Pv、偶次散射Pd、表面散射Ps
步骤2.对特征数据矩阵进行低秩稀疏分解,得到低秩特征矩阵X。
本步骤的具体实现如下:
(2a)对极化SAR特征数据进行低秩稀疏分解,有X=L+S+G,X表示输入的特征矩阵A,α,H,Pd,Ps,Pv,L表示特征数据矩阵的低秩部分,即特征数据矩阵中的所有样本可以由有限数量的基线性表示,S表示特征数据矩阵的稀疏部分,即其不能由在低秩中获得的基近似的来解释,G表示特征数据矩阵中的噪声部分,即不能由低秩中的基和稀疏值来解释的随机噪声。优化模型如下:
Figure BDA0001354317720000042
通过简单交替最小化的方法:在每次迭代中,首先优化L与固定S,然后优化S与固定 L。rank(L)表示矩阵L的秩,r是矩阵L的最大的秩,card(S)表示矩阵S的基数,k是矩阵S的最大基数。
Figure BDA0001354317720000046
表示矩阵X-L-S即噪声矩阵G的F范数的平方,其中F范数是一个矩阵中所有数值模的平方和的开方;
(2b)利用上述极化SAR特征数据的低秩稀疏分解的结果,则低秩特征矩阵X中的第i 行j列像素点xij可以表示为:
Figure BDA0001354317720000043
其中
Figure BDA0001354317720000044
Figure BDA0001354317720000045
分别表示各向异性度矩阵A、散射角α、散射功率熵H、偶次散射Pd、表面散射Ps和体散射Pv经过低秩稀疏分解后的相应低秩部分在第i行j列处的值。
步骤3.选取训练样本占标记样本比例为5%的训练样本矩阵Y。
本步骤的具体实现如下:
在各个类别中按照5%的训练样本比例选取训练样本构成训练样本矩阵yh
Y=[y1,...,yh,...,yN]
其中N表示训练样本的总数。
步骤4.构造一种相似度权值度量因子W,对于新增的未标记样本,给其标上类标签,实现增量式的极化SAR特征数据标记分类。
本步骤的具体实现如下:
(4a)计算第i行j列未标记样本xij与训练样本yh的相似度权值度量因子Wh,得到相似度距离矩阵W:
W=[W1,W2,...,Wh,...,WN]
其中N表示训练矩阵Y中的样本总数;
(4a1)计算第i行j列未标记样本xij与训练样本yh的相似度权值度量因子Wh
Figure BDA0001354317720000051
Figure BDA0001354317720000052
表示第i行j列未标记样本xij的第k个特征,yh表示训练样本矩阵Y中的类标签值为第h个训练样本;
Figure BDA0001354317720000053
表示训练样本yh的第k个特征,k∈{1,2,3,4,5,6};λ表示值参数,
Figure BDA0001354317720000054
(4a2)对于第i行j列未标记样本xij与训练矩阵相似度矩阵Y中的每一个训练样本求相似度权值因子,得到相似度权值矩阵W=[W1,W2,...,Wh,...,WN],N表示训练矩阵Y中的样本总数。
(4b)在相似度权值矩阵W上求得相似度权值因子最大值Wmax,并得到xij取得相似度权值因子最大值Wmax时相对应的训练样本的样本标签值lz
(4b1)在相似度权值矩阵W上求得相似度权值因子最大值Wmax
Wmax=max(W)
(4c)引入条件阈值τ,τ∈[0.970,0.995]。τ的初始值取0.980,变化步长为0.005。每遍历一遍未标记样本,τ的值振荡式增加或减少0.005;同时设定终止阈值B为图像数据大小的千分之一;
(4d)将相似度权值因子最大值Wmax与条件阈值τ比较,若Wmax>τ,则另lij=lz,其中lij表示未标记样本xij的样本标签值,同时将样本xij作为第N+1个训练样本放入到训练样本矩阵Y中,N=N+1。
(4e)计算未标记样本数目m=M-N,M表示原始未标记样本数目与训练样本数目之和,将未标记样本数目m与迭代终止阈值B做比较:
若m<B,则对预测标签矩阵L1进行局部邻域的微调,得到最终的分类结果;
若m≥B,则返回步骤(4a)。
步骤5.对预测标签矩阵L1进行局部邻域的微调。
(5e1)针对预测标签矩阵L1中的标签值为0的元素lp,q构建调整矩阵N1
Figure BDA0001354317720000061
其中lp,q表示预测标签矩阵L1中第p行q列的预测标签值;lp,q+1表示预测标签矩阵L1中第p行q+1列的预测标签优化值;lp+1,q表示预测标签矩阵L1第p+1行q列的预测标签优化值;lp+1,q+1表示预测标签矩阵L1中第p+1行q+1列的预测标签优化值;
(5e2)将第p行q+1列的预测标签值lp,q+1、第p+1行q列的预测标签值lp+1,q和第p+1行q+1列的预测标签值lp+1,q+1出现次数最多的值赋给标签值为0的元素lp,q,得到微调后的预测标签优化矩阵L2即为最终的分类结果。
本发明的效果可以通过以下实验进一步说明。
1.实验数据
本发明所用实验数据为Flevoland地区真实的极化SAR图像的大小为300×270的子地区 Flevoland1。Flevoland1地区的极化SAR图像是由NASA/JPL ARISAR获取的荷兰Flevoland1地区的L波段数据,它的真实地物图如图2所示,其包括马铃薯、甜菜、裸地、大麦、小麦和豌豆6类不同的地物,用不同的颜色显示。
2.本发明使用的对比实验方法如下所述:
对比方法1是基于H/α分解及层次聚类的分类方法。
对比方法2是K最近邻分类方法。
3.实验内容
实验1,对图2中Flevoland1地区大小为300×270的极化SAR图像,用本发明和对比分类方法进行分类,分类精度见表1,分类结果如图3,其中图(3a)表示H/α分解及层次聚类方法的分类结果,图(3b)表示K近邻分类方法的分类结果,图(3c)表示本发明的分类结果;
表1对荷兰Flevoland1地区,不同方法分类正确率:
Figure BDA0001354317720000071
4.实验结果分析
从表1中可以看出,在Flevoland1地区上本发明的分类精度明显高于H/α分解及层次聚类方法和K近邻分类方法的分类精度,分类精度平均高出11%以上。
由图3可以看出,在区域性一致的极化SAR地物的分类结果上图3(c)明显优于图3(a) 和图3(b),即本发明方法明显优于对比方法。
综上所述,本发明提出的基于低秩稀疏分解和增量学习的极化SAR地物分类方法,充分利用了极化SAR图像的空间相关性,减少了噪声干扰,将增量学习运用到极化SAR分类中,实现了数据的可持续性学习,提高了算法的分类精度。所以,本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。

Claims (5)

1.一种基于低秩稀疏分解和增量学习的极化SAR地物分类方法,包括:
(1)输入待分类的极化SAR图像的极化相干矩阵T,并对极化相干矩阵T进行统计特征的提取,得到相应的特征数据矩阵;
(2)对特征数据矩阵进行低秩稀疏分解,得到低秩特征矩阵X;
(3)选取训练样本占标记样本比例为5%的训练样本矩阵Y;
(4)构造一种相似度权值度量因子W,对于新增的未标记样本,给其标上类标签,实现增量式的极化SAR特征数据标记分类;
(4a)计算第i行j列未标记样本xij与训练样本yh的相似度权值度量因子Wh,得到相似度距离矩阵W:
W=[W1,W2,...,Wh,...,WN]
其中N表示训练矩阵Y中的样本总数;
(4b)在相似度权值矩阵W上求得相似度权值因子最大值Wmax,并得到xij取得相似度权值因子最大值Wmax时相对应的训练样本的样本标签值lz
(4c)引入条件阈值τ,τ∈[0.970,0.995];τ的初始值取0.980,变化步长为0.005,每遍历一遍未标记样本,τ的值振荡式增加或减少0.005;同时设定终止阈值B为图像数据大小的千分之一;
(4d)将相似度权值因子最大值Wmax与条件阈值τ比较,若Wmax>τ,则另lij=lz,其中lij表示未标记样本xij的样本标签值,同时将样本xij作为第N+1个训练样本放入到训练样本矩阵Y中,N=N+1;
(4e)计算未标记样本数目m=M-N,M表示原始未标记样本数目与训练样本数目之和,将未标记样本数目m与迭代终止阈值B做比较;
若m<B,则对预测标签矩阵L1进行局部邻域的微调,得到最终的分类结果;
若m≥B,则返回步骤(4a)。
2.根据权利要求1中所述的基于低秩稀疏分解和增量学习的极化SAR地物分类方法,所述步骤(1)中的对极化相干矩阵T进行统计特征的提取,得到相应的特征数据矩阵,按如下步骤进行:
(1a)利用极化SAR数据的相干矩阵T的特征值分解得到每个像素点的散射功率熵H、散射角α和各向异性度A:
H=-(P1log3P1+P2log3P2+P3log3P3)
Figure FDA0002787499470000021
其中,Pt表示第t种散射机制发生的概率,Pt∈[0,1],αt表示第t种散射机制的类型,αi=0表示目标散射体发生表面散射,αi=π/4为体散射,表示目标散射体发生偶极子散射;αi=π/2为偶次散射,表示目标散射体发生二面角散射;λt表示极化相干矩阵T的第t个特征值,t∈[1,3];
(1b)利用Freeman分解得到每个像素点的三种散射矩阵:体散射Pv、偶次散射Pd、表面散射Ps
3.根据权利要求1所述的基于低秩稀疏分解和增量学习的极化SAR地物分类方法,其特征在于:所述步骤(2)中对特征数据矩阵进行低秩稀疏分解,得到低秩特征矩阵X,按如下步骤进行:
(2a)对极化SAR特征数据进行低秩稀疏分解,有X=L+S+G,X表示输入的特征矩阵A,α,H,Pd,Ps,Pv,L表示特征数据矩阵的低秩部分,即特征数据矩阵中的所有样本可以由有限数量的基线性表示,S表示特征数据矩阵的稀疏部分,即其不能由在低秩中获得的基近似的来解释,G表示特征数据矩阵中的噪声部分,即不能由低秩中的基和稀疏值来解释的随机噪声,优化模型如下:
Figure FDA0002787499470000022
通过简单交替最小化的方法:在每次迭代中,首先优化L与固定S,然后优化S与固定L,rank(L)表示矩阵L的秩,r是矩阵L的最大的秩,card(S)表示矩阵S的基数,k是矩阵S的最大基数,
Figure FDA0002787499470000023
表示矩阵X-L-S即噪声矩阵G的F范数的平方,其中F范数是一个矩阵中所有数值模的平方和的开方;
(2b)利用上述极化SAR特征数据的低秩稀疏分解的结果,则低秩特征矩阵X中的第i行j列像素点xij可以表示为:
Figure FDA0002787499470000031
其中
Figure FDA0002787499470000032
Figure FDA0002787499470000033
分别表示各向异性度矩阵A、散射角α、散射功率熵H、偶次散射Pd、表面散射Ps和体散射Pv经过低秩稀疏分解后的相应低秩部分在第i行j列处的值。
4.根据权利要求1所述的基于低秩稀疏分解和增量学习的极化SAR地物分类方法,其特征在于:所述步骤(4a)中计算第i行j列未标记样本xij与训练样本yh的相似度权值度量因子Wh,得到相似度距离矩阵W,按如下步骤进行:
(4a1)计算第i行j列未标记样本xij与训练样本yh的相似度权值度量因子Wh
Figure FDA0002787499470000034
Figure FDA0002787499470000035
表示第i行j列未标记样本xij的第k个特征,yh表示训练样本矩阵Y中的类标签值为第h个训练样本;
Figure FDA0002787499470000036
表示训练样本yh的第k个特征,k∈{1,2,3,4,5,6};λ表示值参数,
Figure FDA0002787499470000037
(4a2)对于第i行j列未标记样本xij与训练矩阵相似度矩阵Y中的每一个训练样本求相似度权值因子,得到相似度权值矩阵W=[W1,W2,...,Wh,...,WN],N表示训练矩阵Y中的样本总数。
5.根据权利要求1所述的基于低秩稀疏分解和增量学习的极化SAR地物分类方法,其特征在于:所述步骤(4e)中对预测标签矩阵L1进行局部邻域的微调,按如下步骤进行:
(4e1)针对预测标签矩阵L1中的标签值为0的元素lp,q构建调整矩阵N1
Figure FDA0002787499470000038
其中lp,q表示预测标签矩阵L1中第p行q列的预测标签值;lp,q+1表示预测标签矩阵L1中第p行q+1列的预测标签优化值;lp+1,q表示预测标签矩阵L1第p+1行q列的预测标签优化值;lp+1,q+1表示预测标签矩阵L1中第p+1行q+1列的预测标签优化值;
(4e2)将第p行q+1列的预测标签值lp,q+1、第p+1行q列的预测标签值lp+1,q和第p+1行q+1列的预测标签值lp+1,q+1出现次数最多的值赋给标签值为0的元素lp,q,得到微调后的预测标签优化矩阵L2即为最终的分类结果。
CN201710589213.6A 2017-07-19 2017-07-19 基于低秩稀疏分解和增量学习的极化sar地物分类方法 Active CN107392140B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710589213.6A CN107392140B (zh) 2017-07-19 2017-07-19 基于低秩稀疏分解和增量学习的极化sar地物分类方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710589213.6A CN107392140B (zh) 2017-07-19 2017-07-19 基于低秩稀疏分解和增量学习的极化sar地物分类方法

Publications (2)

Publication Number Publication Date
CN107392140A CN107392140A (zh) 2017-11-24
CN107392140B true CN107392140B (zh) 2021-01-08

Family

ID=60340040

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710589213.6A Active CN107392140B (zh) 2017-07-19 2017-07-19 基于低秩稀疏分解和增量学习的极化sar地物分类方法

Country Status (1)

Country Link
CN (1) CN107392140B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108280466B (zh) * 2018-01-12 2021-10-29 西安电子科技大学 基于加权核范数最小化的极化sar特征分类方法
CN108564006B (zh) * 2018-03-26 2021-10-29 西安电子科技大学 基于自步学习卷积神经网络的极化sar地物分类方法
CN110210574B (zh) * 2019-06-13 2022-02-18 中国科学院自动化研究所 合成孔径雷达图像解译方法、目标识别装置及设备
CN113344013B (zh) * 2021-07-15 2023-06-30 西安电子科技大学 基于图低秩稀疏分解的极化sar图像特征挖掘方法
CN116797845B (zh) * 2023-07-05 2024-01-26 中国科学院空天信息创新研究院 一种基于散射机制的无监督简缩极化分类方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103413146A (zh) * 2013-08-23 2013-11-27 西安电子科技大学 基于Freeman熵和自学习的极化SAR图像精细分类方法
CN103793713A (zh) * 2012-10-31 2014-05-14 中国科学院沈阳自动化研究所 低秩约束的在线自监督学习的场景分类方法
EP2767849A2 (en) * 2014-01-13 2014-08-20 Institute of Electronics, Chinese Academy of Sciences Method and apparatus for processing polarimetric synthetic aperture radar image
CN104361346A (zh) * 2014-10-21 2015-02-18 西安电子科技大学 基于k-svd和稀疏表示的极化sar图像分类方法
CN104504393A (zh) * 2014-12-04 2015-04-08 西安电子科技大学 基于集成学习的极化sar图像半监督分类方法
CN105184298A (zh) * 2015-08-27 2015-12-23 重庆大学 一种快速局部约束低秩编码的图像分类方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103793713A (zh) * 2012-10-31 2014-05-14 中国科学院沈阳自动化研究所 低秩约束的在线自监督学习的场景分类方法
CN103413146A (zh) * 2013-08-23 2013-11-27 西安电子科技大学 基于Freeman熵和自学习的极化SAR图像精细分类方法
EP2767849A2 (en) * 2014-01-13 2014-08-20 Institute of Electronics, Chinese Academy of Sciences Method and apparatus for processing polarimetric synthetic aperture radar image
CN104361346A (zh) * 2014-10-21 2015-02-18 西安电子科技大学 基于k-svd和稀疏表示的极化sar图像分类方法
CN104504393A (zh) * 2014-12-04 2015-04-08 西安电子科技大学 基于集成学习的极化sar图像半监督分类方法
CN105184298A (zh) * 2015-08-27 2015-12-23 重庆大学 一种快速局部约束低秩编码的图像分类方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于内容的图像哈希检索算法研究;季秀云;《中国优秀硕士学位论文全文数据库信息科技辑》;20150515(第5期);第22页 *
基于分布式低秩表示的子空间聚类算法;许凯等;《计算机研究与发展》;20160715;第1605-1611页 *
融合光谱-空间信息的高光谱遥感影像增量分类算法;王俊淑等;《测绘学报》;20150531;第44卷(第9期);第1003-1013页 *

Also Published As

Publication number Publication date
CN107392140A (zh) 2017-11-24

Similar Documents

Publication Publication Date Title
CN107392140B (zh) 基于低秩稀疏分解和增量学习的极化sar地物分类方法
Ma et al. Infrared and visible image fusion methods and applications: A survey
WO2022001159A1 (zh) 无监督高光谱图像隐低秩投影学习特征提取方法
CN108446716B (zh) 基于FCN与稀疏-低秩子空间表示融合的PolSAR图像分类方法
Zhao et al. Detail-preserving smoothing classifier based on conditional random fields for high spatial resolution remote sensing imagery
CN110210545B (zh) 基于迁移学习的红外遥感水体分类器构建方法
Han et al. Combining 3D‐CNN and Squeeze‐and‐Excitation Networks for Remote Sensing Sea Ice Image Classification
CN102402685B (zh) 基于Gabor特征的三马尔可夫场SAR图像分割方法
Boggavarapu et al. A new framework for hyperspectral image classification using Gabor embedded patch based convolution neural network
CN107491734B (zh) 基于多核融合与空间Wishart LapSVM的半监督极化SAR图像分类方法
CN110598564B (zh) 基于OpenStreetMap的高空间分辨率遥感影像迁移学习分类方法
CN112990334A (zh) 基于改进原型网络的小样本sar图像目标识别方法
CN108446582A (zh) 基于纹理特征和仿射传播聚类算法的高光谱图像分类方法
CN109712149B (zh) 一种基于小波能量和模糊c均值的图像分割方法
CN106874862B (zh) 基于子模技术和半监督学习的人群计数方法
Huo et al. A batch-mode active learning algorithm using region-partitioning diversity for SVM classifier
CN113723492A (zh) 一种改进主动深度学习的高光谱图像半监督分类方法及装置
Yang et al. Unsupervised images segmentation via incremental dictionary learning based sparse representation
CN106611422A (zh) 基于素描结构的随机梯度贝叶斯sar图像分割方法
CN105160351A (zh) 基于锚点稀疏图的半监督高光谱分类方法
CN104123563B (zh) 一种基于Cloude分解的极化SAR图像分类方法
CN103456017B (zh) 基于种子集的半监督权重核模糊聚类的图像分割方法
CN110717531A (zh) 一种基于不确定性分析和贝叶斯融合的分类后变化类型检测方法
CN107464247B (zh) 一种基于g0分布的随机梯度变分贝叶斯sar图像分割方法
CN104680182B (zh) 基于nsct和判别字典学习的极化sar分类方法

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