CN113920345B - 一种基于聚类多重流形测度学习的高光谱影像降维方法 - Google Patents
一种基于聚类多重流形测度学习的高光谱影像降维方法 Download PDFInfo
- Publication number
- CN113920345B CN113920345B CN202111053158.1A CN202111053158A CN113920345B CN 113920345 B CN113920345 B CN 113920345B CN 202111053158 A CN202111053158 A CN 202111053158A CN 113920345 B CN113920345 B CN 113920345B
- Authority
- CN
- China
- Prior art keywords
- manifold
- distance
- training
- matrix
- data
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 230000009467 reduction Effects 0.000 title claims abstract description 24
- 238000012549 training Methods 0.000 claims abstract description 53
- 238000005259 measurement Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 39
- 238000012545 processing Methods 0.000 claims description 7
- 230000005855 radiation Effects 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 7
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 101100379080 Emericella variicolor andB gene Proteins 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003331 infrared imaging Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- 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/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
Abstract
本发明提供了一种基于聚类多重流形测度学习的高光谱影像降维方法,在多重局部测度学习的基础上,引入流形学习。通过多重测度学习,将原本复杂的非线性数据分解称为多个相互关联数据,使得在训练过程中能够更好的把握数据结构中的细节部分,在有针对性的进行训练的同时,还在一定程度上避免了过拟合,能够很好的处理非线性数据。在上述基础上,利用数据标签信息的同时,还能够有效利用数据本身的结构信息,最大化的将有限的先验信息利用起来。在降维的同时,还使得最后求得的特征空间能够有效的将不同类区分开来,对影像后续的分类起到帮助。本发明的有益效果是:线性分线性数据都能进行处理,同时避免了分类计算中的过拟合,提高了分类精度。
Description
技术领域
本发明涉及遥感图像处理技术领域,尤其涉及一种基于聚类多重流形测度学习的高光谱影像降维方法。
背景技术
高光谱遥感影像处理在信息探测方面有着重要的地位,是遥感领域的重要课题。高光谱图像包含数百个光谱波段,提供了详细的空间结构和光谱信息。因此,通过高光谱影像可以区分地表覆盖类别之间的细微差异,广泛应用于大气监测、植被监测、灾害预警等领域等。随着传感器的不断发展,数据的维数也在不断提高,然而数据维数的提高并没有使得影像对于地物的分类和识别能力也得到提高,反而会在影像的处理过程中带来一系列的问题,使得影像分类精度下降。为了解决上述问题,就有研究提出了高光谱影像降维技术。通过高光谱降维技术,可以有效的降低数据的复杂程度,消除数据的冗余信息,避免了休斯现象(F.L.Luo,2017)。
高光谱数据降维,指的是通过处在低维空间中的数据来有效的表示处在原始高维空间中的数据,在减少数据维数同时,还能最大程度地保存有效信息,去除冗余数据,有助于提高高光谱数据在后续处理中的表现力。本质上,降维学习通过假设高斯模型寻找高维特征空间的低维表示,并在新的特征空间中区分类内和类间样本。测度学习可以利用先验样本信息最小化同一类别的距离,最大化不同类别的距离,实现不同类别的分离,满足降维的要求。例如,大间隔最近邻(LMNN)(K.Q.Weinberger and L.K.Saul,2009)通过训练使样本点周围的最近邻与样本点属于同一类,而不同类的样本点相距较远。信息理论测度学习(ITML)(J.V.Davis,B.Kulis,P.Jain,S.Sra,and I.S.Dhillon,2007)通过最小化将马氏测度学习转换为解决Bregman优化问题两个多元高斯分布之间的微分相对熵。以上算法都是基于全局测度学习,当面对异构分布的数据时,不足以获得令人满意的结果。因此有研究提出了一些多重局部测度学习算法。例如,聚类多重测度学习(CMML)(B.Nguyen,F.J.Ferri,C.Morell,andB.DeBaets,2019)。但是全局测度学习方法和多局部测度学习方法都忽略了数据的结构特征,不能利用数据本身的结构信息。这意味着该算法没有充分利用训练样本提供的信息。而高光谱数据中先验样本的获取成本十分昂贵,所以说如何最大化的利用有限的样本信息十分重要。
为了解决上述问题,本发明方法在多重局部测度学习的基础上,引入流形学习。流形旨在找出数据中嵌入高维空间中的低维流形结构,能够很好的利用数据本身的结构信息。所以本发明方法在多重局部测度学习的基础上,引入流形学习,结合多重局部测度和流形学习的优势,更好地将类与类区分开来。
发明内容
为了解决上述问题,本发明提供了一种基于聚类多重流形测度学习的高光谱影像降维方法,具体步骤如下:
S1:对获取的高光谱遥感图像进行处理,得到若干个训练样本和训练样本的类别;
S2:根据高光谱遥感图像X中各波段的像元辐射值和所述若干个训练样本,使用k-means算法将高光谱遥感图像聚类为Q个簇,每个训练样本只属于一个簇,每一个簇作为一个单独的训练单元进行训练;
S3:对簇中训练样本xi,i表示第i个训练样本,根据欧式距离选择最近的k个同类点视为邻近点,然后以训练样本为中心点联合其邻近点构建为一个流形Si,中心点xi的类别即为流形Si的类别;
S4:根据邻近点在重构样本点中的贡献,计算重构系数wiq;
S5:以中心点之间的欧式距离为标准,选择距离最近的同类别流形以及距离最近的异类别流形;
S6:结合步骤S4中计算得到的重构系数,计算流形Si与Sj之间的距离
S7:将流形Si以及它的同类别流形Sj,异类别流形Sl构建为一个三重约束;
S8:根据步骤S6中计算的流形距离,对步骤S7中的三重约束进行训练,得到局部测度矩阵
S9:将所述局部测度矩阵代入步骤S3中,采用步骤S4-8的方法,直至使用所有的训练样本都完成对三重约束的训练后,得到最终的局部测度矩阵M(c);
S10:根据所述局部测度矩阵M(c),对待分析高光谱遥感图像的原始数据进行特征空间变换,得到降维后的数据。
进一步地,步骤S1中,所述处理为,将获取的高光谱遥感图像读入读入大小为d×n的矩阵X中,矩阵X中各列元素为各波段对应的像元辐射值,其中{x1,x2,…,xn}∈Rd×n表示训练样本有n个,d为遥感图像的波段数,n为图像的像元数,将对应的标签读入大小为1×n的矩阵Y中,矩阵Y中各行元素为各个训练样本对应的类别标签,{y1,y2,…,yn}∈R1×n表示训练样本对应的类别标签有n个。
进一步地,计算重构系数wiq的公式如下:
其中,n为训练样本个数,k为邻近点个数,xiq为中心点xi的第q个邻近点,通过最小化该公式,得到最终的重构系数。
进一步地,流形距离的计算公式如下:
其中,代表中心点xi和中心点xj之间的马氏距离,θ为超参数,wiq为中心点xi与临近点q的重构系数,wjq为中心点xj与临近点q的重构系数,表示中心点xi和xjq之间的马氏距离,xjq表示中心点xj的第q个邻近点,表示中心点xj和xiq之间的马氏距离,xiq表示中心点xi的第q个邻近点。
进一步地,所述三重约束为:
其中,r表示三重约束,表示同类流形之间的距离,/>表示异类流形之间的距离,Si表示第i个流形,Sj表示第j个流形,Sl表示第l个流形。
进一步地,得到最终的局部测度矩阵M(c)的公式和限定条件如下:
M(0),...,M(Q)>0
其中M(0)代表全局测度矩阵,α,β为超参数,N(c)为三重约束r的个数,ξr为三重约束中的松弛变量,τ(c)代表第c个簇中所有三重约束的集合,Q表示簇的总个数,代表基于流形距离计算得到的损失值,M(0),...,M(Q)>0代表M(0),...,M(Q)均为正定矩阵。
本发明提供的技术方案带来的有益效果是:
1、通过多重测度学习,将原本复杂的非线性数据分解称为多个相互关联数据,使得在训练过程中能够更好的把握数据结构中的细节部分,在有针对性的进行训练的同时,还在一定程度上避免了过拟合,能够很好的处理非线性数据。
2、在测度学习的基础上,引入流形学习的思想,在利用数据标签信息的同时,还能够有效利用数据本身的结构信息,最大化的将有限的先验信息利用起来。在降维的同时,还使得最后求得的特征空间能够有效的将不同类区分开来,对影像后续的分类起到帮助。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例中一种基于聚类多重流形测度学习的高光谱影像降维方法的流程图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图详细说明本发明的具体实施方式。
本发明的实施例提供了一种基于聚类多重流形测度学习的高光谱影像降维方法,在多重局部测度学习的基础上,引入流形学习。本具体实施方式采用MATLAB平台实现,MATLAB遥感图像读写函数为实施基础。调用遥感图像读取函数,输入待读取遥感图像文件名,遥感图像就被读入大小为d×n的矩阵中,矩阵中各元素为各波段对应的像元辐射值,其中,d为遥感图像的波段数,n为遥感图像的像元数。调用遥感图像读写函数,将待分析高光谱遥感图像读入矩阵X,对应的标签读入矩阵Y。
请参考图1,图1是本发明实施例中一种基于聚类多重流形测度学习的高光谱影像降维方法的流程图,基于聚类多重流形测度学习进行目标探测,具体操作如下:
(1)根据各波段的像元辐射值和训练样本,调用MATLAB自带k-means算法将X聚类为Q个簇,簇与簇之间相互独立,即每个训练样本只属于一个簇。得到的每一个簇都被视为一个单独的训练单元进行训练;
(2)对簇中每一个训练样本xi,根据欧式距离选择最近的k个同类点视为邻近点,然后将训练样本及其邻近点构建为一个流形Si。欧式距离计算公式如下:
其中xi是X中第i个训练样本的像素辐射值。
(3)根据邻近点在重构样本点中的贡献,计算重构系数wiq,具体公式如下:
其中,n为训练样本个数,k为邻近点个数,xiq为样本点xi的第q个邻近点,通过最小化该公式,可以得到最终的重构系数wiq。每次训练,都是将一个样本点作为一个中心点。从公式中可以看到,邻近点xiq与中心点xi的距离越近,在重构过程中做出的贡献就越大,最终求得的权重系数也就越大。
(4)以流形中心点之间的欧式距离为标准,选择距离最近的同类别流形以及距离最近的异类别流形。
(5)结合步骤(3)中计算得到重构系数wiq,计算流形Si与Sj之间的距离
其中代表样本点xi和样本点xj之间的马氏距离,θ为超参数,用于平衡样本点与邻近点对于最终计算得到的流形距离的贡献。通过计算流形距离的方式,在传统测度学习的基础上,添加了数据本身的结构信息,使得训练过程能够更为充分的利用已有的信息。wiq为中心点xi与临近点q的重构系数,wjq为中心点xj与临近点q的重构系数,/>表示中心点xi和点xjq之间的马氏距离,/>表示中心点xi和点xiq之间的马氏距离,xjq表示中心点xj的第q个邻近点,xiq表示中心点xi的第q个邻近点。
(6)将流形Si以及它的同类别流形Sj,异类别流形Sl构建为一个三重约束。具体公式如下:其中,r表示三重约束,要求同类流形之间的距离要小于异类流形之间的距离,/>表示同类流形之间的距离,表示异类流形之间的距离,Si表示第i个流形,Sj表示第j个流形,Sl表示第l个流形。
(7)根据步骤(5)中得到的流形距离,对所述三重约束进行训练,得到局部测度矩阵然后将得到的局部测度矩阵/>代入步骤(5)中,重复步骤(5)-(7)。直至使用所有的训练样本都完成对三重约束的训练后,得到最终的局部测度矩阵M(c):
M(0),...,M(Q)>0
其中,M(0)代表全局测度矩阵,通过该矩阵可以将每个簇求得的局部测度矩阵关联起来,实现簇与簇之间的信息共享。α,β为超参数,用于平衡低秩约束部分和全局测度矩阵约束对于最后结果的影响。N(c)为三重约束的个数。ξr为三重约束中松弛变量。τ(c)代表第c个簇中所有三重约束的集合。
(8)通过随机梯度下降法求得最佳解,针对局部测度矩阵的求解,取步骤(7)中公式按如下操作进行求解:
M(c)>0
其中,M(0)代表全局测度矩阵,α,β为超参数,N(c)为三重约束r的个数,ξr为三重约束中的松弛变量,τ(c)代表第c个簇中所有三重约束的集合,Q表示簇的总个数,代表基于流形距离计算得到的损失值,M(0),...,M(Q)>0代表M(0),...,M(Q)均为正定矩阵。
将上面的公式转换为每个三重约束损失函数之和,如下所示:
(9)每次迭代过程是针对一组三重约束进行的,这里的三元组指的就是一组三重约束,在第i次迭代时,用基于单个三元组约束的近似替换步骤(8)中的函数,如下所示:
其中max(A,B)表示将A于B进行比较,最终输出较大值,若A>B,则max(A,B)=A。
(10)对步骤(9)中公式进行求导计算得到次梯度:
其中表示当括号内的条件成立时,结果值为1,反之则为0。/>代表第t次迭过程中计算得到的局部测度矩阵。
(11)通过步骤(10)中求得的次梯度,将更新到/>更新过程如下:
其中λmin是的最小特征值,umin是λmin对应的特征向量。当迭代完成后所求即为最终的局部测度矩阵。
(12)根据步骤(11)中得到的局部测度矩阵M(c),对原始数据进行特征空间变换,即可得到降维后的数据。
通过对比试验来验证本发明的效果,本试验采用的数据为1)Airborne Visible/Infrared Imaging Spectrometer(AVIRIS)传感器于1996年3月23日在佛罗里达州肯尼迪航天中心(KSC)获得的KSC数据集,图像大小为512像元×614像元,空间分辨率为18米,光谱分辨率为10nm,波段范围为400nm到2500nm,共176个可用波段;2)Center forAirborneLaser Mapping(NCALM)获取的休斯顿大学校园和邻近市区数据,图像大小为349像元×1905像元,空间分辨率为1.3米,波段范围为380nm到1050nm,共144个可用波段。3)Reflective Optics Spectrographic Imaging System(ROSIS)传感器获取的意大利的帕维亚城数据,图像大小为610像元×340像元,空间分辨率为1.3米,光谱分辨率为10nm,波段范围为430nm到860nm,共103个可用波段;分别采用经典的邻域保留嵌入(方法1)(X.F.He,D.Cai,S.C.Yan,H.J.Zhang,and S.O.C.,2005)、局部保持投影(方法2)(X.F.He andP.Niyogi,2004)、局部Fisher判别分析(方法3)(W.Li,S.Prasad,J.E.Fowler,andL.M.Bruce,2012)、大间隔最近邻方法(方法4)、信息论测度学习(方法5)、最大崩塌测度学习(方法6)(A.Globerson and S.T.Roweis,2006)、聚类多重测度学习(方法7)和本发明方法进行降维分类。
分类评价指标:1)Kappa系数:
Kappa系数是用来评价分类问题的权威评价指标。Kappa系数越大,精度越高。本试验中,选取的是方法1~7所能获得的最高Kappa系数来评价方法1~7的分类能力。
根据样本获取如表1所示的混淆矩阵:
表1混淆矩阵
表1中,Nnn表示未变化样本被检测成为未变化的数量;Nnc表示未变化样本被检测成为变化的数量;Ncn表示变化样本被检测成为未变化的数量;Ncc表示变化样本被检测成为变化的数量;Nnp为Nnc和Nnc之和,Ncp为Ncn和Ncc之和,Npn为Nnn和Ncp之和,Npc为Nnc和Ncc之和,N为总样本数。
Kappa系数的计算公式为:
1)整体精度:
整体精度(OverallAccuracy,OA)是用来评价分类问题的权威评价指标。整体精度越高,检测精度越高。OA的计算方法同样基于表1所示的混淆矩阵,整体精度OA计算公式为:
采用Kappa系数和整体精度评价经过方法1~7和本发明方法的降维后数据的分类能力,评价指标见表2。
表2对比试验结果
从表2可见,本发明方法能获得更高的整体精度和Kappa系数,表明经过本发明方法降维后的数据具有更好的分类结果。
由此可得出结论,与传统降维方法相比,本发明方法具有更好的降维效果,在完成数据降维的同时还能够有效的提升数据在分类中的表现力。本发明充分利用了有限的先验信息,考虑到了数据本身的结构特征,在训练的过程中能够更好的把握数据的细节部分。
本发明的有益效果是:
1、通过多重测度学习,将原本复杂的非线性数据分解称为多个相互关联数据,使得在训练过程中能够更好的把握数据结构中的细节部分,在有针对性的进行训练的同时,还在一定程度上避免了分类计算中过拟合,能够很好的处理非线性数据。
2、在测度学习的基础上,引入流形学习的思想,在利用数据标签信息的同时,还能够有效利用数据本身的结构信息,最大化的将有限的先验信息利用起来。在降维的同时,还使得最后求得的特征空间能够有效的将不同类区分开来,对影像后续的分类起到帮助提高了分类精度。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于聚类多重流形测度学习的高光谱影像降维方法,其特征在于:包括以下步骤:
S1:对获取的高光谱遥感图像进行处理,得到若干个训练样本和训练样本的类别;
S2:根据高光谱遥感图像中各波段的像元辐射值和若干个训练样本,使用k-means算法将高光谱遥感图像聚类为Q个簇,每个训练样本只属于一个簇,每一个簇作为一个单独的训练单元进行训练;
S3:对簇中训练样本xi,i表示第i个训练样本,根据欧式距离选择最近的k个同类点视为邻近点,然后以训练样本为中心点联合其邻近点构建为一个流形Si,中心点xi的类别即为流形Si的类别;
S4:根据邻近点在重构样本点中的贡献,计算重构系数wiq;
S5:以中心点之间的欧式距离为标准,选择距离最近的同类别流形以及距离最近的异类别流形;
S6:结合步骤S4中计算得到的重构系数,计算流形Si与Sj之间的距离
S7:将流形Si以及它的同类别流形Sj,异类别流形Sl构建为一个三重约束;
S8:根据步骤S6中计算的流形距离,对步骤S7中的三重约束进行训练,得到局部测度矩阵
S9:将所述局部测度矩阵代入步骤S3中,采用步骤S4-8的方法,直至使用所有的训练样本都完成对三重约束的训练后,得到最终的局部测度矩阵M(c);
S10:根据所述局部测度矩阵M(c),对待分析高光谱遥感图像的原始数据进行特征空间变换,得到降维后的数据。
2.如权利要求1所述的一种基于聚类多重流形测度学习的高光谱影像降维方法,其特征在于:步骤S1中,所述处理为,将获取的高光谱遥感图像读入大小为d×n的矩阵X中,矩阵X中各列元素为各波段对应的像元辐射值,其中{x1,x2,…,xn}∈Rd×n表示训练样本有n个,d为遥感图像的波段数,n为图像的像元数,将对应的标签读入大小为1×n的矩阵Y中,矩阵Y中各行元素为各个训练样本对应的类别标签,{y1,y2,…,yn}∈R1×n表示训练样本对应的类别标签有n个。
3.如权利要求1所述的一种基于聚类多重流形测度学习的高光谱影像降维方法,其特征在于:步骤S4中,计算重构系数wiq的公式如下:
其中,n为训练样本个数,k为邻近点个数,xiq为中心点xi的第q个邻近点,通过最小化该公式,得到最终的重构系数。
4.如权利要求1所述的一种基于聚类多重流形测度学习的高光谱影像降维方法,其特征在于:步骤S6中,流形距离的计算公式如下:
其中,代表中心点xi和中心点xj之间的马氏距离,θ为超参数,wiq为中心点xi与临近点q的重构系数,wjq为中心点xj与临近点q的重构系数,/>表示中心点xi和xjq之间的马氏距离,xjq表示中心点xj的第q个邻近点,/>表示中心点xj和xiq之间的马氏距离,xiq表示中心点xi的第q个邻近点。
5.如权利要求1所述的一种基于聚类多重流形测度学习的高光谱影像降维方法,其特征在于:步骤S7中,所述三重约束为:
其中,r表示三重约束,表示同类流形之间的距离,/>表示异类流形之间的距离,Si表示第i个流形,Sj表示第j个流形,Sl表示第l个流形。
6.如权利要求1所述的一种基于聚类多重流形测度学习的高光谱影像降维方法,其特征在于:步骤S9中,得到最终的局部测度矩阵M(c)的公式和限定条件如下:
其中,M(0)代表全局测度矩阵,α,β为超参数,N(c)为三重约束r的个数,ξr为三重约束中的松弛变量,τ(c)代表第c个簇中所有三重约束的集合,Q表示簇的总个数,代表基于流形距离计算得到的损失值,/>代表M(0),...,M(Q)均为正定矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111053158.1A CN113920345B (zh) | 2021-09-09 | 2021-09-09 | 一种基于聚类多重流形测度学习的高光谱影像降维方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111053158.1A CN113920345B (zh) | 2021-09-09 | 2021-09-09 | 一种基于聚类多重流形测度学习的高光谱影像降维方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113920345A CN113920345A (zh) | 2022-01-11 |
CN113920345B true CN113920345B (zh) | 2024-04-02 |
Family
ID=79234531
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111053158.1A Active CN113920345B (zh) | 2021-09-09 | 2021-09-09 | 一种基于聚类多重流形测度学习的高光谱影像降维方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113920345B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115861683B (zh) * | 2022-11-16 | 2024-01-16 | 西安科技大学 | 一种针对高光谱图像的快速降维方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104751191A (zh) * | 2015-04-23 | 2015-07-01 | 重庆大学 | 一种稀疏自适应半监督多流形学习的高光谱影像分类方法 |
CN105069482A (zh) * | 2015-08-21 | 2015-11-18 | 中国地质大学(武汉) | 基于空间正则化流形学习算法的高光谱遥感图像分类方法 |
CN109117712A (zh) * | 2018-06-26 | 2019-01-01 | 中国地质大学(武汉) | 一种最大边缘测度学习的高光谱影像目标探测方法及系统 |
CN110135479A (zh) * | 2019-04-29 | 2019-08-16 | 中国地质大学(武汉) | 基于随机森林测度学习的高光谱图像目标探测方法及系统 |
CN110378272A (zh) * | 2019-07-12 | 2019-10-25 | 河海大学 | 基于矩阵分块Isomap算法的高光谱遥感影像特征提取方法 |
CN112329654A (zh) * | 2020-11-10 | 2021-02-05 | 中国地震局地震预测研究所 | 基于多流形学习算法的高光谱影像数据分类方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160132754A1 (en) * | 2012-05-25 | 2016-05-12 | The Johns Hopkins University | Integrated real-time tracking system for normal and anomaly tracking and the methods therefor |
-
2021
- 2021-09-09 CN CN202111053158.1A patent/CN113920345B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104751191A (zh) * | 2015-04-23 | 2015-07-01 | 重庆大学 | 一种稀疏自适应半监督多流形学习的高光谱影像分类方法 |
CN105069482A (zh) * | 2015-08-21 | 2015-11-18 | 中国地质大学(武汉) | 基于空间正则化流形学习算法的高光谱遥感图像分类方法 |
CN109117712A (zh) * | 2018-06-26 | 2019-01-01 | 中国地质大学(武汉) | 一种最大边缘测度学习的高光谱影像目标探测方法及系统 |
CN110135479A (zh) * | 2019-04-29 | 2019-08-16 | 中国地质大学(武汉) | 基于随机森林测度学习的高光谱图像目标探测方法及系统 |
CN110378272A (zh) * | 2019-07-12 | 2019-10-25 | 河海大学 | 基于矩阵分块Isomap算法的高光谱遥感影像特征提取方法 |
CN112329654A (zh) * | 2020-11-10 | 2021-02-05 | 中国地震局地震预测研究所 | 基于多流形学习算法的高光谱影像数据分类方法及系统 |
Non-Patent Citations (2)
Title |
---|
基于局域自适应信息理论测度学习的高光谱目标探测方法;董燕妮;杜博;张乐飞;张良培;;武汉大学学报(信息科学版);20180709(08);全文 * |
空-谱协同流形重构的高光谱影像分类;黄鸿;陈美利;段宇乐;石光耀;;光学精密工程;20180715(07);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113920345A (zh) | 2022-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111860612B (zh) | 无监督高光谱图像隐低秩投影学习特征提取方法 | |
CN108830296B (zh) | 一种改进的基于深度学习的高分遥感影像分类方法 | |
CN110717354B (zh) | 基于半监督k-svd与多尺度稀疏表示的超像元分类方法 | |
CN110287869B (zh) | 基于深度学习的高分辨率遥感影像农作物分类方法 | |
CN109389080B (zh) | 基于半监督wgan-gp的高光谱图像分类方法 | |
CN110619263B (zh) | 低秩联合协同表示的高光谱遥感影像异常探测方法 | |
Zhang et al. | A modified stochastic neighbor embedding for multi-feature dimension reduction of remote sensing images | |
CN109190511B (zh) | 基于局部与结构约束低秩表示的高光谱分类方法 | |
CN105989597B (zh) | 基于像素选择过程的高光谱图像异常目标检测方法 | |
Rangarajan et al. | Detection of fusarium head blight in wheat using hyperspectral data and deep learning | |
Fang et al. | Confident learning-based domain adaptation for hyperspectral image classification | |
CN113920345B (zh) | 一种基于聚类多重流形测度学习的高光谱影像降维方法 | |
CN110399814B (zh) | 一种基于局部线性表示领域适应度量的人脸识别方法 | |
Qi et al. | Global-local three-dimensional convolutional transformer network for hyperspectral image classification | |
CN114937173A (zh) | 一种基于动态图卷积网络的高光谱图像快速分类方法 | |
Ling et al. | Semi-supervised learning via convolutional neural network for hyperspectral image classification | |
CN113139512B (zh) | 基于残差和注意力的深度网络高光谱影像分类方法 | |
Zhao et al. | Hyperspectral target detection method based on nonlocal self-similarity and rank-1 tensor | |
Dhandapani et al. | Multi-Channel Convolutional Neural Network for Prediction of Leaf Disease and Soil Properties | |
CN110717485B (zh) | 一种基于局部保留投影的高光谱图像稀疏表示分类方法 | |
CN110619370A (zh) | 一种高光谱图像超像素局部线性嵌入降维方法 | |
Wei et al. | The hyperspectral image clustering based on spatial information and spectral clustering | |
CN110188821B (zh) | 多参数边缘保持滤波和多特征学习的高光谱图像分类方法 | |
CN109145945B (zh) | 一种非局部加权联合稀疏表示的高光谱图像分类方法 | |
CN113902973A (zh) | 自编码器和低维流形建模的高光谱异常检测方法 |
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 |