CN104318243B - 基于稀疏表示和空谱拉普拉斯图的高光谱数据降维方法 - Google Patents
基于稀疏表示和空谱拉普拉斯图的高光谱数据降维方法 Download PDFInfo
- Publication number
- CN104318243B CN104318243B CN201410542949.4A CN201410542949A CN104318243B CN 104318243 B CN104318243 B CN 104318243B CN 201410542949 A CN201410542949 A CN 201410542949A CN 104318243 B CN104318243 B CN 104318243B
- Authority
- CN
- China
- Prior art keywords
- training sample
- data
- dimension
- point
- sample point
- 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
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/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2323—Non-hierarchical techniques based on graph theory, e.g. minimum spanning trees [MST] or graph cuts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/58—Extraction of image or video features relating to hyperspectral data
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Discrete Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种用于大规模高光谱数据的降维方法,主要用于解决传统流行学习信息单一和该类方法难以处理较大规模数据的问题。其实现步骤为:1.从大规模的高光谱数据中选择一定量的数据作为训练样本;2.对训练样本进行空谱拉普拉斯图的构造;3.对拉普拉斯矩阵进行特征分解得到训练样本的低维表示;4.利用训练样本及其低维表示构造高维字典和低维字典;5.计算剩余高光谱数据在高维字典上的稀疏表示系数;6.将该稀疏表示系数与低维字典进行相乘,得到剩余数据的低维表示;7.整合训练样本与剩余数据的低维表示得到完整的降维数据。本发明提升了流行降维的效果,可用于处理大规模的高光谱数据。
Description
技术领域
本发明属于数据处理技术领域,涉及高光谱数据的前期处理,主要目的是为了减少高光谱数据的维数,从而降低后期数据处理方法的计算复杂度,同时尽可能提升其性能。该方法可以被应用于大规模的高光谱数据聚类或分类中。
背景技术
数据降维处理在数据处理中起着很大的作用,许多维数过高的数据在处理之前都会进行降维处理,一方面可以降低计算量,另一方面还可以从原有特征中挖取更有用的特征,提升后期算法的处理效果。光谱数据随着成像设备的光谱分辨率的不断提高,数据的维数也是越来越高,数据降维必不可少,同时,随着设备的发展,空间分辨率也在不断提高,数据的规模也在不断的增大,如何处理大规模的高光谱数据也成为非常关键的一个问题。
现有的数据降维方法很多,常用的如:主成分分析PCA,线性判别分析LDA,局部保持投影LPP,拉普拉斯嵌入。主成分分析及线性判别分析方法简单实用,但是适合于线性数据,对于非线性的数据处理效果不是很好。以往研究表明,高光谱数据中存在着流形结构,线性方法并不能完全保持高光谱的数据背景。流形学习针对于非线性数据,利用图嵌入的方法捕捉数据的空间结构,将数据映射到具有同样空间结构的低维流行空间中,从而保持数据间的分布结构。
目前流形学习降维的方法有不少,如:
2000年Tenenbaum与Silva在《Science》上提出了ISOMAP,这个方法是利用非线性的局部变量信息学习数据集的全局集合结构,使用了测地距离来度量高维空间中的样本点距离,通过建立原数据的测地线距离与降维数据空间的空间距离的对等关系完成数据降维。该方法保证流形上的空间结构在低维流行空间中依然存在,但是当选择较大邻域时会出现短路现象。
2000年Roweis和Saul提出了局部线性嵌入法(Locally Linear Embedding,LLE),该方法的主要思想是具有低维子流形结构的数据集,原空间与低维空间中的点的邻域结构关系式不变的。该方法很好的保留了邻接点间的关系,使每个点的邻接权值保持不变,但对于等距流形,嵌入效果不是很好。
2003年M.Belkin与P.Niyogi提出了拉普拉斯特征映射LE,该方法的出发点为:高维空间中离得很近的点投影到低维空间中的像应该也会离得很近。该方法处理分类问题很好,但是权值计算所使用的heat kernel中的参数对嵌入结构有很大影响。
上述方法有两个统一的缺陷:(1)这些方法中很重要的步骤就是图的构造,当数据规模非常大的时候,图的存储及后期的计算都是非常困难的,一般的流形学习方法无法处理大规模数据;(2)普通的流形学习方法,并没有考虑到高光谱数据中所存在的空间结构,只是单纯的考虑其谱间的邻域关系,导致对高光谱数据降维效果不理想。
发明内容
本发明的目的在于克服上述已有技术的缺点,提出了一种基于稀疏表示和空谱拉普拉斯图的高光谱数据降维方法,以提高高光谱数据降维的效果,便于将流行学习能推广到大规模的高光谱数据中。
本发明的技术方案是:从大规模的高光谱数据中选择一定量的数据作为训练样本,对所选训练样本进行空谱拉普拉斯图的构造,对拉普拉斯矩阵进行特征分解得到训练样本的低维表示;利用高维训练样本及其低维表示构造高维字典和低维字典,将剩余的高光谱数据在高维字典上进行稀疏表示,得到对应的稀疏表示系数;将该稀疏表示系数与低维字典进行相乘,得到剩余高光谱数据的低维表示,整合训练样本与剩余高光谱数据的低维表示得到整体数据的低维表示。其具体步骤包括如下:
(1)从一幅高光谱图像数据I中选择n个数据点作为高维的训练样本,高光谱数据维数为p,n的数值由高光谱图像数据的规模确定,取整体数目的10%以上;
(2)对所选高维训练样本进行空谱拉普拉斯图G的构造:
(2a)构造谱间图G1:
使用谱信息散度SID作为训练样本点间的距离度量,计算第i个训练样本与其它训练样本间的距离,i=1,…,n,并对这些距离值进行由小到大排序,选择距离最小的N个样本作为第i个训练样本点的N近邻,N的值根据具体的实验数据进行设置;
根据第i个训练样本点的N近邻确定第i个训练样本点与其它训练样本点的连接关系:若第j个训练样本点在第i个训练样本点的N近邻中,则将第j个训练样本点与第i个训练样本点连接,并计算该连接边的权值反之,第j个训练样本点与第i个训练样本点不连接,Wij'=0,其中x,y分别为第i个训练样本点与第j个训练样本点所对应的光谱向量,参数t根据实际数据调试确定;
(2b)构造空间图G2:
比较第i个训练样本点与其它训练样本点的二维坐标,i=1,…,n,确定其它训练样本点是否在第i个训练样本点的K邻域中,若第j个训练样本点在第i个训练样本点的K邻域内,将第i个训练样本点与第j个训练样本点进行连接,反之第i个训练样本点与第j个训练样本点不连接,邻域参数K=11,该参数表示以第i个训练样本点为中心的11*11的邻域区域;
确定连接边的权值:将11*11的邻域划分为内邻域和外邻域,内邻域为以第i个训练样本点为中心的5*5的区域,外邻域为除去内邻域的剩余邻域区域;如果第j个训练样本点在第i个训练样本点的内邻域中,则连接边的权值为Wij”=1,如果第j个训练样本点在第i个训练样本点的外邻域中,则连接边的权值Wij”=0.8;若第i个训练样本点与第j个训练样本点间不存在连接,则Wij”=0;
(2c)将谱间图G1和空间图G2进行合并操作,保留这两个图中的所有连接边,得到空谱拉普拉斯图G,得到空谱拉普拉斯图G的权值矩阵为W,W=W'+W”,计算拉普拉斯矩阵L,L=D-W,其中D为由W的行或列求和得到的向量作为对角线元素的对角矩阵;
(3)对拉普拉斯矩阵L和对角矩阵D进行广义特征值分解,取最小r个特征值对应的特征向量作为训练样本所对应的低维表示TR;
(4)构造高维空间与低维空间的对偶字典:将n个p维的训练样本作为高维字典HD,将n个训练样本对应的r维表示TR作为低维字典LD,这两个字典的原子间存在一一对应的关系;
(5)对剩余高光谱数据进行稀疏表示求解,得到剩余高光谱数据在高维字典HD上的稀疏表示系数:Θ=[θ1,...,θs,...,θm],θs为第s个数据点的稀疏表示系数,s=1,...,m,m为剩余高光谱数据的个数;
(6)将剩余高光谱数据的稀疏表示系数Θ与低维字典LD相乘,得到剩余高光谱数据的r维表示RR=LD*Θ;
(7)结合训练样本的r维表示TR,得到整个高光谱数据的r维表示IR=[TR;RR]。
本发明具有如下优点:
1)本发明由于在构造谱间图时使用谱信息度量SID来度量光谱数据的相似度,能更加准确的描述光谱数据间的光谱域邻域结构;
2)本发明由于在构造空间谱时使用了分层邻域结构,使得空间域邻域结构更加精细;
3)本发明由于采用谱间图和空间图共同构成拉普拉斯图,所以能更好的表示高光谱数据的流行结构;
4)本发明由于使用稀疏表示的方法来模拟高维空间与低维空间的对应关系,从部分高光谱数据的低维表示中学习得到完整高光谱数据的低维表示,使得流行学习降维方法不再受到数据规模的影响,可将其应用于处理大规模的高光谱数据中。
实验证明,本发明通过构造空谱拉普拉斯图,提升了高光谱数据降维的效果,通过使用训练样本及其低维表示来表示高维空间和低维空间,利用稀疏表示学习得到剩余的高光谱数据的低维表示,打破了流行学习对数据规模的限制,可将其应用到更大规模的数据中。
附图说明
图1是本发明的整体实现流程图;
图2是本发明仿真所用到的数据的位置坐标图。
具体实施方式
参照图1,本发明的具体实现步骤如下:
步骤1,从一幅高光谱图像数据I中选择n个数据点作为高维的训练样本,高光谱数据维数为p,n的数值由高光谱图像数据的规模确定,取整体数目的10%以上。
步骤2,通过对训练样本进行分析,构造空谱拉普拉斯图G。
(2a)构造谱间图G1:
(2a.1)谱信息散度SID是一种光谱数据间的光谱相似度的度量,与一般的欧氏距离相比较,能够更好的捕获光谱数据间的相似性,因此使用谱信息散度SID作为谱间图的距离度量,使谱间图能更准确捕获训练样本点间的相似关系。谱信息散度SID的定义如下:
SID(x,y)=D(x||y)+D(y||x),
其中:x,y为光谱数据中的光谱向量,均为p维列向量,p等于光谱数据的光谱数目,y=(y1,...,yp)T,y所对应的概率向量为q=(q1,...,qi,...,qp)T,其中x=(x1,...,xp)T,x所对应的概率向量为e=(e1,...,ej,..,ep)T,其中上式中D(x||y)和D(y||x)分别由下式计算得到:
构造谱间图需要确定每个训练样本与其它训练样本间的关系,针对第i个训练样本,计算该训练样本与其它训练样本间的距离,并对这些距离值进行由小到大排序,选择距离最小的N个样本作为第i个训练样本点的N近邻,近邻参数N的值可以根据具体的实验数据进行设置,本实验中设置N=6;
(2a.2)根据第i个训练样本点的N近邻确定第i个训练样本点与其它训练样本点的连接关系:若第j个训练样本点在第i个训练样本点的N近邻中,则将第j个训练样本点与第i个训练样本点连接,并计算该连接边的权值反之,第j个训练样本点与第i个训练样本点不连接,Wij'=0,其中x,y分别为第i个训练样本点与第j个训练样本点所对应的光谱向量,参数t根据实际数据调试确定,本实例中设置t=0.01;
(2b)构造空间图G2:
(2b.1)构造空间图来表示训练样本点间的空间结构,由于每个高光谱数据都具有自己的空间坐标,可以通过比较光谱数据的空间坐标来分析它们之间的空间结构。比较第i个训练样本点与其它训练样本点的二维坐标,确定其它训练样本点是否在第i个训练样本点的K邻域中,若第j个训练样本点在第i个训练样本点的K邻域内,将第i个训练样本点与第j个训练样本点进行连接,反之第i个训练样本点与第j个训练样本点不连接,邻域参数K表示以第i个训练样本点为中心的K*K的邻域区域,该邻域参数K取值为奇数,如:3、7、9、11、21等,本实验中设置K=11;
(2b.2)使用邻域分层的方法来确定空间图的连接边权值,通过对空间邻域中的数据点进行更细的划分,将空间结构关系表现的更加准确:
将K*K的邻域划分为内邻域和外邻域,内邻域为以第i个训练样本点为中心的K1*K1的区域,K1<K,本实例中设置K1=5,外邻域为除去内邻域的剩余邻域区域;
如果第j个训练样本点在第i个训练样本点的内邻域中,则连接边的权值为Wij”=1,如果第j个训练样本点在第i个训练样本点的外邻域中,则连接边的权值Wij”=0.8;若第i个训练样本点与第j个训练样本点间不存在连接,则Wij”=0;
(2c)将谱间图G1和空间图G2进行合并操作,得到空谱拉普拉斯图G,该图G中不仅包含光谱域的信息还包含了空间域的信息,该空谱拉普拉斯图G的权值矩阵为:W=W'+W”,计算拉普拉斯矩阵:L=D-W,其中D为由W的行或列求和得到的向量作为对角线元素的对角矩阵。
步骤3,对拉普拉斯矩阵L和对角矩阵D进行广义特征值分解,对角矩阵D的逆矩阵存在,将L与D的广义特征值问题转化为D-1L的一般特征值问题,通过特征值分解得到n个特征值λ1,λ2,...,λn,n为方阵D-1L的行数,这n个特征值按照从小到大的顺序排列,即:λ1<λ2,...,<λn,以及对应的特征向量u1,u2,...,un,取最小的r个特征向量值对应的特征向量u1,u2,...,ur作为训练样本的r维表示TR,r表示降维后的数据维数,该参数可以根据实验数据设置,本实例中r=4。
步骤4,构造高维字典和低维字典,训练样本中的数据点作为高维字典HD的原子,训练样本的r维表示TR中的数据点作为低维字典LD的原子,高维字典与低维字典的原子间保持一一对应的关系,将高维字典原子看作高维空间的基原子,高维字典即代表了整个高维空间,同样,低维字典代表了整个低维空间。
步骤5,通过稀疏表示的方法确定剩余的高光谱数据在高维空间中的表示;剩余高光谱数据在高维字典HD上的稀疏表示系数:Θ=[θ1,...,θs,...,θm],θs为第s个数据点的稀疏表示系数,s=1,...,m,m为剩余高光谱数据的个数,通过最小化下式中的目标函数,得到解向量θ,使稀疏表示系数θs等于该解向量θ:
其中,xs为第s个数据点对应的光谱向量,||*||2为向量的2范数,||*||1为向量的1范数,β为模型调节参数,在本实例中设置β=0.1。
对上式中θ的求解,已有很多的成熟算法,最小绝对值收缩选择算子LASSO是其中使用非常广泛的一种求解方法,该方法是1996年由Robert Tibshirani提出,通过对表示系数中的某些系数原子进行收缩操作,并将其它的系数原子设置为0,从而保留下更重要的系数原子,在本实例中使用了SparseLab实验室工具包中的lasso函数来进行求解。
步骤6,将剩余高光谱数据的稀疏表示系数Θ与低维字典LD相乘,得到剩余高光谱数据的r维表示RR=LD*Θ,由于高维字典和低维字典的原子之间存在一对一的关系,因此,高维空间中的稀疏表示关系在低维空间中依然保持,可通过稀疏表示系数与低维字典计算得到剩余数据的低维表示。
步骤7,结合训练样本的r维表示TR,得到整个高光谱数据的r维表示IR=[TR;RR]。本发明的效果可以通过仿真实验具体说明:
1.实验条件
实验所用微机CPU为Intel i3 3.2GHz内存4GB,编程平台是MatlabR2010a。实验中采用的数据为高光谱图像数据,是1992年由AVIRIS传感器在印第安拉州拍摄的Indian_Pines高光谱图像数据,该图像尺寸为145×145,一共有220个波段,剔除噪声严重的20个波段,剩余200个波段。实验中所使用数据为原数据的部分数据,具体情况见表1,该实验数据的位置坐标图见图2,图中黑色的位置代表实验数据的空间位置。
表1
2.实验内容
使用本发明的方法在不同训练样本比率下对高光谱数据进行降维,然后再对降维后的数据进行K-mean聚类,计算聚类准确度ACC,训练样本的选择比率包含:10%,20%,30%,40%,K-mean聚类中的类别参数设置为4。
为了验证方法的有效性,对原始的高光谱数据和经PCA降维后的数据进行K-mean聚类作为对比实验;此外,为证明本发明中所用到的空谱拉普拉斯图对降维效果的影响,用以欧式距离作为距离度量的谱间N近邻图和使用不分层9*9邻域的空间图分别替换空谱拉普拉斯图进行实验。
聚类准确度ACC定义如下:
其中,cn为正确聚类的数据个数,n为训练样本的个数,m为剩余高光谱数据的个数。
3.实验结果
分别对原始数据、使用PCA方法和使用本发明方法对原始数据进行降维后的数据进行K-mean聚类,实验结果见表2。
表2
Method | Original | PCA | 10% | 20% | 30% | 40% |
ACC(%) | 68.1679 | 67.7714 | 75.3348 | 77.3998 | 78.4705 | 78.3117 |
表2中Original表示对原始数据进行K-mean聚类,PCA表示对原始数据进行PCA降维后进行K-mean聚类,10%、20%、30%、40%是流行学习所使用的训练样本比例,分别表示本发明的方法在相应训练样本比例下对原始数据进行降维,然后再进行K-mean聚类。
通过表2可以看出:本发明的方法虽然只是通过对部分高光谱数据的进行流行降维学习,可以得到比原始数据和使用PCA降维后的数据更好的聚类结果,由此可见,本发明的方法可通过对部分数据进行流行学习,实现对大规模高光谱数据的降维。
分别使用欧式距离作为度量的谱间图和使用不分层空间邻域的空间图替换本发明中的空谱拉普拉斯图对原数据进行降维,然后再进行K-mean聚类,实验结果见表3,
表3
Method | 10% | 20% | 30% | 40% |
SSLaplace | 75.3348 | 77.3998 | 78.4705 | 78.3117 |
G_s | 70.4935 | 71.3482 | 72.1639 | 73.6775 |
G_r | 71.8352 | 73.0829 | 73.7398 | 74.2538 |
表3中SSLaplace表示本方法中所使用的空谱拉普拉斯图,G_s表示使用欧式距离作为度量的谱间图,G_r表示邻域不分层的空间图。从表3中可以看出,本发明中使用的空谱拉普拉斯图与传统欧式距离谱间图、不分层空间图相比降维效果更好。
Claims (3)
1.一种基于稀疏表示和空谱拉普拉斯图的高光谱数据降维方法,包括以下步骤:
(1)从一幅高光谱图像数据I中选择n个数据点作为高维的训练样本,高光谱数据维数为p,n的数值由高光谱图像数据的规模确定,取整体数目的10%以上;
(2)对所选高维训练样本进行空谱拉普拉斯图G的构造:
(2a)构造谱间图G1:
使用谱信息散度SID作为训练样本点间的距离度量,计算第i个训练样本与其它训练样本间的距离,i=1,…,n,并对这些距离值进行由小到大排序,选择距离最小的N个样本作为第i个训练样本点的N近邻,N的值根据具体的实验数据进行设置;
根据第i个训练样本点的N近邻确定第i个训练样本点与其它训练样本点的连接关系:若第j个训练样本点在第i个训练样本点的N近邻中,则将第j个训练样本点与第i个训练样本点连接,并计算该连接边的权值反之,第j个训练样本点与第i个训练样本点不连接,W′ij=0,其中x,y分别为第i个训练样本点与第j个训练样本点所对应的光谱向量,参数t根据实际数据调试确定;
(2b)构造空间图G2:
比较第i个训练样本点与其它训练样本点的二维坐标,i=1,…,n,确定其它训练样本点是否在第i个训练样本点的K邻域中,若第j个训练样本点在第i个训练样本点的K邻域内,将第i个训练样本点与第j个训练样本点进行连接,反之第i个训练样本点与第j个训练样本点不连接,邻域参数K=11,该参数表示以第i个训练样本点为中心的11*11的邻域区域;
确定连接边的权值:将11*11的邻域划分为内邻域和外邻域,内邻域为以第i个训练样本点为中心的5*5的区域,外邻域为除去内邻域的剩余邻域区域;如果第j个训练样本点在第i个训练样本点的内邻域中,则连接边的权值为W″ij=1,如果第j个训练样本点在第i个训练样本点的外邻域中,则连接边的权值W″ij=0.8;若第i个训练样本点与第j个训练样本点间不存在连接,则W″ij=0;
(2c)将谱间图G1和空间图G2进行合并操作,保留这两个图中的所有连接边,得到空谱拉普拉斯图G,得到空谱拉普拉斯图G的权值矩阵为W,W=W'+W”,计算拉普拉斯矩阵L,L=D-W,其中D为由W的行或列求和得到的向量作为对角线元素的对角矩阵;
(3)对拉普拉斯矩阵L和对角矩阵D进行广义特征值分解,取最小r个特征值对应的特征向量作为训练样本所对应的低维表示TR;
(4)构造高维空间与低维空间的对偶字典:将n个p维的训练样本作为高维字典HD,将n个训练样本对应的r维表示TR作为低维字典LD,这两个字典的原子间存在一一对应的关系;
(5)对剩余高光谱数据进行稀疏表示求解,得到剩余高光谱数据在高维字典HD上的稀疏表示系数:Θ=[θ1,...,θs,...,θm],θs为第s个数据点的稀疏表示系数,s=1,...,m,m为剩余高光谱数据的个数;
(6)将剩余高光谱数据的稀疏表示系数Θ与低维字典LD相乘,得到剩余高光谱数据的r维表示RR=LD*Θ;
(7)结合训练样本的r维表示TR,得到整个高光谱数据的r维表示IR=[TR;RR]。
2.根据权利要求1所述的基于稀疏表示和空谱拉普拉斯图的高光谱数据降维方法,其中步骤(3)所述的对拉普拉斯矩阵L和对角矩阵D进行广义特征值分解,按如下步骤进行:
(3.1)将广义特征值问题转化为一般特征值问题:D-1Lu=λu,其中D-1为对角矩阵D的逆矩阵,λ为特征值,u为特征值λ对应的特征向量;
(3.2)对D-1L进行一般特征值分解得到n个特征值λ1,λ2,...,λn,n为方阵D-1L的行数,这n个特征值按照从小到大的顺序排列,即:λ1<λ2,...,<λn,以及对应的特征向量u1,u2,...,un,取最小的r个特征向量值对应的特征向量u1,u2,...,ur作为训练样本的r维表示TR,r表示降维后的数据维数,该参数可以根据实验数据设置。
3.根据权利要求1所述的基于稀疏表示和空谱拉普拉斯图的高光谱数据降维方法,其中步骤(5)所述的对剩余高光谱数据进行稀疏表示求解,是对每个数据点分别进行求解:
(5.1)设剩余高光谱数据在高维字典HD上的稀疏表示系数为:Θ=[θ1,...,θs,...,θm],θs为第s个数据点的稀疏表示系数,s=1,...,m,m为剩余高光谱数据的个数;
(5.2)最小化下式中的目标函数,得到对应的解向量θ,使稀疏表示系数θs等于该解向量θ:
其中,xs为第s个数据点对应的光谱向量,||.||2为向量的2范数,||.||1为向量的1范数,β为调节参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410542949.4A CN104318243B (zh) | 2014-10-14 | 2014-10-14 | 基于稀疏表示和空谱拉普拉斯图的高光谱数据降维方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410542949.4A CN104318243B (zh) | 2014-10-14 | 2014-10-14 | 基于稀疏表示和空谱拉普拉斯图的高光谱数据降维方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104318243A CN104318243A (zh) | 2015-01-28 |
CN104318243B true CN104318243B (zh) | 2017-09-26 |
Family
ID=52373472
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410542949.4A Active CN104318243B (zh) | 2014-10-14 | 2014-10-14 | 基于稀疏表示和空谱拉普拉斯图的高光谱数据降维方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104318243B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105574548B (zh) * | 2015-12-23 | 2019-04-26 | 北京化工大学 | 一种基于稀疏和低秩表示图的高光谱数据降维方法 |
CN105654517A (zh) * | 2016-02-22 | 2016-06-08 | 江苏信息职业技术学院 | 基于分层空间的rb化粒子滤波算法 |
CN106778032B (zh) * | 2016-12-14 | 2019-06-04 | 南京邮电大学 | 药物设计中配体分子海量特征筛选方法 |
CN107798345B (zh) * | 2017-10-20 | 2020-11-20 | 西北工业大学 | 基于块对角与低秩表示的高光谱伪装目标检测方法 |
CN109670418B (zh) * | 2018-12-04 | 2021-10-15 | 厦门理工学院 | 结合多源特征学习和组稀疏约束的无监督物体识别方法 |
CN109858531B (zh) * | 2019-01-14 | 2022-04-26 | 西北工业大学 | 一种基于图的高光谱遥感图像快速聚类算法 |
CN110580463B (zh) * | 2019-08-30 | 2021-07-16 | 武汉大学 | 单一光谱驱动的双类别稀疏表示高光谱图像目标检测方法 |
CN110648276B (zh) * | 2019-09-25 | 2023-03-31 | 重庆大学 | 基于流形映射与字典学习的高维图像数据降维方法 |
CN110929793A (zh) * | 2019-11-27 | 2020-03-27 | 谢国宇 | 用于生态环境监测的时空域模型建模方法及*** |
CN111079850B (zh) * | 2019-12-20 | 2023-09-05 | 烟台大学 | 波段显著度的深度空谱联合高光谱图像分类方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102938072A (zh) * | 2012-10-20 | 2013-02-20 | 复旦大学 | 一种基于分块低秩张量分析的高光谱图像降维和分类方法 |
EP2597596A2 (en) * | 2011-11-22 | 2013-05-29 | Raytheon Company | Spectral image dimensionality reduction system and method |
CN103413151A (zh) * | 2013-07-22 | 2013-11-27 | 西安电子科技大学 | 基于图正则低秩表示维数约简的高光谱图像分类方法 |
CN103996047A (zh) * | 2014-03-04 | 2014-08-20 | 西安电子科技大学 | 基于压缩谱聚类集成的高光谱图像分类方法 |
-
2014
- 2014-10-14 CN CN201410542949.4A patent/CN104318243B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2597596A2 (en) * | 2011-11-22 | 2013-05-29 | Raytheon Company | Spectral image dimensionality reduction system and method |
CN102938072A (zh) * | 2012-10-20 | 2013-02-20 | 复旦大学 | 一种基于分块低秩张量分析的高光谱图像降维和分类方法 |
CN103413151A (zh) * | 2013-07-22 | 2013-11-27 | 西安电子科技大学 | 基于图正则低秩表示维数约简的高光谱图像分类方法 |
CN103996047A (zh) * | 2014-03-04 | 2014-08-20 | 西安电子科技大学 | 基于压缩谱聚类集成的高光谱图像分类方法 |
Non-Patent Citations (2)
Title |
---|
Semisupervised Dual-Geometric Subspace;Shuyuan Yang等;《IEEE Transactions on Geoscience and Remote Sensing》;20140630;第52卷(第6期);第3866-3869页 * |
Semi-supervised Hyperspectral Image;Tatyana V等;《Geoscience and Remote Sensing Symposium》;20070630;第3587-3593页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104318243A (zh) | 2015-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104318243B (zh) | 基于稀疏表示和空谱拉普拉斯图的高光谱数据降维方法 | |
Yuan et al. | Factorization-based texture segmentation | |
Yuan et al. | Remote sensing image segmentation by combining spectral and texture features | |
WO2020103417A1 (zh) | 一种bmi评测方法、装置及计算机可读存储介质 | |
CN113362382A (zh) | 三维重建方法和三维重建装置 | |
CN108319957A (zh) | 一种基于超点图的大规模点云语义分割方法 | |
Li et al. | A multi-scale cucumber disease detection method in natural scenes based on YOLOv5 | |
KR102667737B1 (ko) | 특징점 포지셔닝 방법 및 장치 | |
CN104484658A (zh) | 一种基于多通道卷积神经网络的人脸性别识别方法及装置 | |
CN113191489B (zh) | 二值神经网络模型的训练方法、图像处理方法和装置 | |
CN107341505B (zh) | 一种基于图像显著性与Object Bank的场景分类方法 | |
CN112560967B (zh) | 一种多源遥感图像分类方法、存储介质及计算设备 | |
CN109726725B (zh) | 一种基于大间隔类间互异性多核学习的油画作者识别方法 | |
CN108460391A (zh) | 基于生成对抗网络的高光谱图像无监督特征提取方法 | |
CN107316005B (zh) | 基于稠密轨迹核协方差描述子的行为识别方法 | |
CN105654122B (zh) | 基于核函数匹配的空间金字塔物体识别方法 | |
CN110148103A (zh) | 基于联合优化的高光谱和多光谱图像融合方法、计算机可读存储介质、电子设备 | |
CN107862680B (zh) | 一种基于相关滤波器的目标跟踪优化方法 | |
CN107832786A (zh) | 一种基于字典学习的人脸识别分类方法 | |
CN109190511A (zh) | 基于局部与结构约束低秩表示的高光谱分类方法 | |
CN103646256A (zh) | 一种基于图像特征稀疏重构的图像分类方法 | |
Xu et al. | Robust PCANet on target recognition via the UUV optical vision system | |
CN116563096B (zh) | 用于图像配准的形变场的确定方法、装置以及电子设备 | |
CN111127490A (zh) | 一种基于循环残差U-Net网络的医学图像分割方法 | |
CN110443169B (zh) | 一种边缘保留判别分析的人脸识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |