CN107392140A - 基于低秩稀疏分解和增量学习的极化sar地物分类方法 - Google Patents
基于低秩稀疏分解和增量学习的极化sar地物分类方法 Download PDFInfo
- Publication number
- CN107392140A CN107392140A CN201710589213.6A CN201710589213A CN107392140A CN 107392140 A CN107392140 A CN 107392140A CN 201710589213 A CN201710589213 A CN 201710589213A CN 107392140 A CN107392140 A CN 107392140A
- Authority
- CN
- China
- Prior art keywords
- mrow
- matrix
- msub
- msubsup
- row
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2136—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on sparsity criteria, e.g. with an overcomplete basis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2415—Classification 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)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Probability & Statistics with Applications (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
本发明公开了一种基于低秩稀疏分解和增量学习的极化SAR地物分类方法。其方案是:1.输入待分类的极化SAR图像极化相干矩阵进行统计特征提取;2.进行低秩稀疏分解得到用于分类的低秩特征矩阵;3.初始化训练样本矩阵;4.构造一种相似性度量公式,给新增未标记样本标上类标签,实现增量式的循坏迭代标记分类。本发明将低秩稀疏分解、增量学习与极化SAR数据的统计特征相结合,在减少噪声干扰、充分保持样本相似性和细节信息的同时,实现增量式的样本标记分类,极大的提高了分类精度。
Description
技术领域:
本发明属于图像处理和遥感技术领域,涉及极化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是一个复矩阵,其每一个样本点表示为:
其中Ti,j表示极化相干矩阵T中第i行j列的样本,和分别表示t21、t31和t32的共轭复数;
(1b)利用极化SAR数据的相干矩阵T的特征值分解得到每个像素点的散射功率熵H、散射角α和各向异性度A:
H=-(P1log3P1+P2log3P2+P3log3P3)
其中,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表示特征数据矩阵中的噪声部分,即不能由低秩中的基和稀疏值来解释的随机噪声。优化模型如下:
通过简单交替最小化的方法:在每次迭代中,首先优化L与固定S,然后优化S与固定 L。rank(L)表示矩阵L的秩,r是矩阵L的最大的秩,card(S)表示矩阵S的基数,k是矩阵S的最大基数。表示矩阵X-L-S即噪声矩阵G的F范数的平方,其中F范数是一个矩阵中所有数值模的平方和的开方;
(2b)利用上述极化SAR特征数据的低秩稀疏分解的结果,则低秩特征矩阵X中的第i 行j列像素点xij可以表示为:
其中和分别表示各向异性度矩阵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:
表示第i行j列未标记样本xij的第k个特征,yh表示训练样本矩阵Y中的类标签值为第h个训练样本;表示训练样本yh的第k个特征,k∈{1,2,3,4,5,6};λ表示值参数,
(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:
其中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地区,不同方法分类正确率:
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)
<mrow>
<msub>
<mi>P</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mi>t</mi>
</msub>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>s</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</msubsup>
<msub>
<mi>&lambda;</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mi>a</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<msub>
<mi>P</mi>
<mi>t</mi>
</msub>
<msub>
<mi>&alpha;</mi>
<mi>t</mi>
</msub>
<mi>A</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>3</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>P</mi>
<mn>3</mn>
</msub>
</mrow>
</mfrac>
</mrow>
其中,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表示特征数据矩阵中的噪声部分,即不能由低秩中的基和稀疏值来解释的随机噪声。优化模型如下:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>L</mi>
<mo>.</mo>
<mi>S</mi>
</mrow>
</munder>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>X</mi>
<mo>-</mo>
<mi>L</mi>
<mo>-</mo>
<mi>S</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mi>F</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>n</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mi>r</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>r</mi>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mi>k</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
通过简单交替最小化的方法:在每次迭代中,首先优化L与固定S,然后优化S与固定L。rank(L)表示矩阵L的秩,r是矩阵L的最大的秩,card(S)表示矩阵S的基数,k是矩阵S的最大基数。表示矩阵X-L-S即噪声矩阵G的F范数的平方,其中F范数是一个矩阵中所有数值模的平方和的开方;
(2b)利用上述极化SAR特征数据的低秩稀疏分解的结果,则低秩特征矩阵X中的第i行j列像素点xij可以表示为:
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>L</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&alpha;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>L</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>H</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>L</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<msub>
<mi>d</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>L</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<msub>
<mi>s</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>L</mi>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<msub>
<mi>v</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>L</mi>
</msubsup>
<mo>&rsqb;</mo>
</mrow>
其中和分别表示各向异性度矩阵A、散射角α、散射功率熵H、偶次散射Pd、表面散射Ps和体散射Pv经过低秩稀疏分解后的相应低秩部分在第i行j列处的值。
4.根据权利要求1所述的基于低秩稀疏分解和增量学习的极化SAR地物分类方法,其特征在于:所述步骤(4a)中计算第i行j列未标记样本xij与训练样本yh的相似度权值度量因子Wh,得到相似度距离矩阵W,按如下步骤进行:
(4a1)计算第i行j列未标记样本xij与训练样本yh的相似度权值度量因子Wh:
<mrow>
<msub>
<mi>W</mi>
<mi>h</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mfrac>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>-</mo>
<mi>&lambda;</mi>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mi>&lambda;</mi>
<mo>|</mo>
<mo>|</mo>
<msubsup>
<mi>y</mi>
<mi>h</mi>
<mi>k</mi>
</msubsup>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<mn>2</mn>
<mo>|</mo>
<mo>|</mo>
<msubsup>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>k</mi>
</msubsup>
<msubsup>
<mi>y</mi>
<mi>h</mi>
<mi>k</mi>
</msubsup>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
</mrow>
表示第i行j列未标记样本xij的第k个特征,yh表示训练样本矩阵Y中的类标签值为第h个训练样本;表示训练样本yh的第k个特征,k∈{1,2,3,4,5,6};λ表示值参数,
(4a2)对于第i行j列未标记样本xij与训练矩阵相似度矩阵Y中的每一个训练样本求相似度权值因子,得到相似度权值矩阵W=[W1,W2,...,Wh,...,WN],N表示训练矩阵Y中的样本总数。
5.根据权利要求1所述的基于低秩稀疏分解和增量学习的极化SAR地物分类方法,其特征在于:所述步骤(4e)中对预测标签矩阵L1进行局部邻域的微调,按如下步骤进行:
(4e1)针对预测标签矩阵L1中的标签值为0的元素lp,q构建调整矩阵N1:
<mrow>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>q</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>q</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>p</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>q</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中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即为最终的分类结果。
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 true CN107392140A (zh) | 2017-11-24 |
CN107392140B 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) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108280466A (zh) * | 2018-01-12 | 2018-07-13 | 西安电子科技大学 | 基于加权核范数最小化的极化sar特征分类方法 |
CN108564006A (zh) * | 2018-03-26 | 2018-09-21 | 西安电子科技大学 | 基于自步学习卷积神经网络的极化sar地物分类方法 |
CN110210574A (zh) * | 2019-06-13 | 2019-09-06 | 中国科学院自动化研究所 | 合成孔径雷达图像解译方法、目标识别装置及设备 |
CN113344013A (zh) * | 2021-07-15 | 2021-09-03 | 西安电子科技大学 | 基于图低秩稀疏分解的极化sar图像特征挖掘方法 |
CN116797845A (zh) * | 2023-07-05 | 2023-09-22 | 中国科学院空天信息创新研究院 | 一种基于散射机制的无监督简缩极化分类方法 |
Citations (6)
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 | 重庆大学 | 一种快速局部约束低秩编码的图像分类方法 |
-
2017
- 2017-07-19 CN CN201710589213.6A patent/CN107392140B/zh active Active
Patent Citations (6)
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)
Title |
---|
季秀云: "基于内容的图像哈希检索算法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
王俊淑等: "融合光谱-空间信息的高光谱遥感影像增量分类算法", 《测绘学报》 * |
许凯等: "基于分布式低秩表示的子空间聚类算法", 《计算机研究与发展》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108280466A (zh) * | 2018-01-12 | 2018-07-13 | 西安电子科技大学 | 基于加权核范数最小化的极化sar特征分类方法 |
CN108280466B (zh) * | 2018-01-12 | 2021-10-29 | 西安电子科技大学 | 基于加权核范数最小化的极化sar特征分类方法 |
CN108564006A (zh) * | 2018-03-26 | 2018-09-21 | 西安电子科技大学 | 基于自步学习卷积神经网络的极化sar地物分类方法 |
CN108564006B (zh) * | 2018-03-26 | 2021-10-29 | 西安电子科技大学 | 基于自步学习卷积神经网络的极化sar地物分类方法 |
CN110210574A (zh) * | 2019-06-13 | 2019-09-06 | 中国科学院自动化研究所 | 合成孔径雷达图像解译方法、目标识别装置及设备 |
CN110210574B (zh) * | 2019-06-13 | 2022-02-18 | 中国科学院自动化研究所 | 合成孔径雷达图像解译方法、目标识别装置及设备 |
CN113344013A (zh) * | 2021-07-15 | 2021-09-03 | 西安电子科技大学 | 基于图低秩稀疏分解的极化sar图像特征挖掘方法 |
CN113344013B (zh) * | 2021-07-15 | 2023-06-30 | 西安电子科技大学 | 基于图低秩稀疏分解的极化sar图像特征挖掘方法 |
CN116797845A (zh) * | 2023-07-05 | 2023-09-22 | 中国科学院空天信息创新研究院 | 一种基于散射机制的无监督简缩极化分类方法 |
CN116797845B (zh) * | 2023-07-05 | 2024-01-26 | 中国科学院空天信息创新研究院 | 一种基于散射机制的无监督简缩极化分类方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107392140B (zh) | 2021-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107392140A (zh) | 基于低秩稀疏分解和增量学习的极化sar地物分类方法 | |
CN110135267B (zh) | 一种大场景sar图像细微目标检测方法 | |
CN103955702B (zh) | 基于深度rbf网络的sar图像地物分类方法 | |
CN106067042B (zh) | 基于半监督深度稀疏滤波网络的极化sar分类方法 | |
CN110598029A (zh) | 基于注意力转移机制的细粒度图像分类方法 | |
CN107944370B (zh) | 基于dccgan模型的极化sar图像分类方法 | |
CN108446582A (zh) | 基于纹理特征和仿射传播聚类算法的高光谱图像分类方法 | |
CN105184309A (zh) | 基于cnn和svm的极化sar图像分类 | |
CN107491734B (zh) | 基于多核融合与空间Wishart LapSVM的半监督极化SAR图像分类方法 | |
Wang et al. | Fruit classification model based on improved Darknet53 convolutional neural network | |
CN104850865A (zh) | 一种多特征迁移学习的实时压缩跟踪方法 | |
CN104732244A (zh) | 基于小波变换、多策略pso和svm集成的遥感图像分类方法 | |
CN105184298A (zh) | 一种快速局部约束低秩编码的图像分类方法 | |
CN104794730B (zh) | 基于超像素的sar图像分割方法 | |
CN104252625A (zh) | 样本自适应多特征加权的遥感图像分类方法 | |
CN107274387A (zh) | 基于进化多目标优化的高光谱遥感影像的端元提取方法 | |
CN103955709B (zh) | 基于加权合成核与tmf的极化sar图像分类方法 | |
CN106443632A (zh) | 基于标签保持多任务因子分析模型的雷达目标识别方法 | |
CN106446935A (zh) | 基于核稀疏表示和空间约束的极化sar图像分类方法 | |
CN103593853A (zh) | 基于联合稀疏表达的遥感影像多尺度面向对象分类方法 | |
CN109359525A (zh) | 基于稀疏低秩的判别谱聚类的极化sar图像分类方法 | |
CN107133648A (zh) | 基于自适应多尺度融合稀疏保持投影的一维距离像识别方法 | |
CN113723492A (zh) | 一种改进主动深度学习的高光谱图像半监督分类方法及装置 | |
CN109034213A (zh) | 基于相关熵原则的高光谱图像分类方法和*** | |
CN103714353A (zh) | 基于视觉先验模型的极化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 |