CN107220662B - 基于全局最优聚类的高光谱图像波段选择方法 - Google Patents
基于全局最优聚类的高光谱图像波段选择方法 Download PDFInfo
- Publication number
- CN107220662B CN107220662B CN201710342464.4A CN201710342464A CN107220662B CN 107220662 B CN107220662 B CN 107220662B CN 201710342464 A CN201710342464 A CN 201710342464A CN 107220662 B CN107220662 B CN 107220662B
- Authority
- CN
- China
- Prior art keywords
- wave band
- cluster
- wave
- band
- high spectrum
- 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
-
- 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/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
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)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明提供了一种基于全局最优聚类的高光谱图像波段选择方法。该方法使用聚类间差距和聚类内差距的比值作为待优化的目标函数,并用二分法和动态规划方法优化获得全局最优聚类结果,然后通过最小化波段线性重构误差的方法,最终完成具有代表性的高光谱图像波段选择。利用本发明方法进行高光谱图像波段选择,可以降低噪声波段被选择的概率,获得更高的分类精度。
Description
技术领域
本发明属图像处理技术领域,具体涉及一种基于全局最优聚类的高光谱图像波段选择方法。
背景技术
高光谱图像含有丰富的光谱信息,在地物目标只具有微小差异时仍然能够对其有较好的识别效果。然而对于高光谱图像,大量的光谱信息常常带来信息的冗余以及过大的计算量。波段选择,作为一种高光谱图像的降维方法,在近年来得到广泛的应用。文献“A.-UsóMartinez-Uso,F.Pla,J.M.Sotoca and P.-Sevilla,"Clustering-Based Hyperspectral Band Selection Using Information Measures,"in IEEETransactions on Geoscience and Remote Sensing,pp.4158-4171.2007.”中提出了一种基于聚类的波段选择方法,通过层次聚类的方法,基于信息论的一些度量指标,首先将高光谱图像的各个波段划分为若干个聚类,接着在每个聚类中选取一个最具有代表性的波段来表示整个聚类,以此达到降维的目的。该方法提出了一种用聚类方法来进行波段选择的思路,其不足在于:首先,它是一种不确定的方法,虽然能够降低聚类内部波段间的差异,但无法度量差异降低到何种程度;其次,在存在噪声波段的情况下,由于噪声波段与其他波段具有较大的差异性,更容易被分到一个单波段的聚类中,使得噪声波段更易于被选择,因而降低了降维的效果。
发明内容
为了克服现有方法的不足,本发明提出了一种基于全局最优聚类的高光谱图像波段选择方法。该方法的主要思想是:使用聚类间差距和聚类内差距的比值作为优化的目标函数,并用二分法和动态规划的方法取得全局最优解,然后通过最小化波段线性重构误差的方法,在每个聚类内选取具有代表性的波段,降低噪声波段被选择的概率。
一种基于全局最优聚类的高光谱图像波段选择方法,其特征在于包括以下步骤:
步骤1:按照pn=(po-M2)/(M1-M2)对高光谱图像进行归一化处理,其中,pn为归一化处理后图像的像素值,po为原高光谱图像的像素值,M1为原高光谱图像中像素值的最大值,M2为原高光谱图像中像素值的最小值;
步骤2:以聚类间差距和聚类内差距的比值为目标函数,使用二分法和动态规划方法优化该目标函数,得到一个全局最优的临界波段向量,以这些临界波段为边界将高光谱图像划分为C个连续波段的聚类;
其中,聚类内差距为聚类间差距为xj为高光谱图像第j个波段的波段向量,c为函数的参数,表示将高光谱图像的前sc个波段划分为c个聚类,0<c<C,C为待生成的聚类个数,为小于最终所要求选择波段数的正整数,si表示波段序号,i=0,…,c,满足0=s0<s1<...<sc≤L,L为高光谱图像的波段数,M(a,b)表示归一化后高光谱图像从a波段到b波段区间内所有波段的平均波段向量,如M(1,L)即表示归一化后高光谱图像所有波段的平均波段向量;
步骤3:对每个聚类,选择一个与类中心欧式距离最近的一个波段为该聚类的一个代表波段,所有代表波段构成已选波段集合;所述的类中心是指该聚类内所有波段的平均波段;
步骤4:在每个聚类内迭代地对剩余波段进行选择,直到所选择波段数量满足要求;
所述的迭代地对剩余波段进行选择的具体过程为:在第k次迭代中,k>0,对任一个尚未被选择的波段向量xi,设其所在的波段聚类为j,Bj,k为第k次迭代中波段聚类j内已经被选择的波段向量构成的矩阵,为在第k次迭代中聚类j内所有已选波段向量加上波段向量xi构成的矩阵,设为在第k次迭代中聚类j内所有尚未被选择的波段扣除波段xi得到的波段集合,并设该集合内所有波段向量的平均波段为分别用矩阵Bj,k和对进行线性重构,以两个重构误差的比值作为对波段向量xi的评价值,计算所有尚未被选择波段的评价值,并将评价值最高的波段加入已选波段集合。
本发明的有益效果是:首先,利用了高光谱图像所具有的物理意义,即相邻波段具有较大的相似性,设计了聚类间差距和聚类内差距比值的目标函数,缩小了聚类算法的解空间,可以得到为全局最优的聚类结果;其次,通过最小化波段线性重构误差的方法,在每个聚类内选取具有代表性的波段,可以降低噪声波段被选择的概率,获得更高的分类精度。
附图说明
图1是本发明的一种基于全局最优聚类的高光谱图像波段选择方法流程图
图2是采用不同方法进行高光谱图像波段选择后的分类精度结果图
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
如图1所示,本发明的一种基于全局最优聚类的高光谱图像波段选择方法,其实现步骤如下:
步骤1:归一化处理,即按照pn=(po-M2)/(M1-M2)对高光谱图像进行归一化处理,其中,pn为归一化处理后图像的像素值,po为原高光谱图像的像素值,M1为原高光谱图像中像素值的最大值,M2为原高光谱图像中像素值的最小值;
步骤2:根据高光谱图像的物理意义,即相邻波段具有较大的相似性,一个聚类内的波段应该连续,据此,本发明方法定义一个待优化目标函数,以缩小聚类算法的解空间,得到全局最优聚类。具体为:
以聚类间差距和聚类内差距的比值为目标函数,则待优化问题为:
其中,为聚类内差距,为聚类间差距,xj为高光谱图像第j个波段的波段向量,C为待生成的聚类个数,c为函数的参数,表示将数据集的前sc个波段划分为c个聚类,满足0<c<C si表示波段序号,i=0,…,c,满足0=s0<s1<...<sc≤L,L为高光谱图像的波段数,M(a,b)表示归一化后高光谱图像从a波段到b波段区间内所有波段的平均波段向量,如M(1,L)即表示归一化后高光谱图像所有波段的平均波段向量。
本发明采用二分法和动态规划方法求解待优化问题。先构造函数g(λ)如下:
可以证明,该函数的零点即为原问题的解,并且该函数为单调递减函数,于是可以用二分法来求得最优解λ*。
使用动态规划的方法求解g(λ):
定义P(i,j)表示将前i个波段划分为j个聚类能得到的Dinter(s)-λ·Dintra(s)的最大值,即:
那么有如下状态转移方程:
其中,Db(k+1,i)=||M(k+1,i)-M(1,L)||2,xz表示归一化后高光谱图像z波段向量,M(k+1,i)表示归一化后高光谱图像从k+1波段到i波段区间内所有波段的平均波段向量。
所有的P(i,j)求解完毕后,P(L,C)就是g(λ)的函数值,L为高光谱图像波段数,C为待生成的聚类个数。在状态转移的过程中,以F(i,j)记录最优解对应的波段划分方案:
根据F(i,j)的定义可知:F(L,C)=sC-1,F(sC-1,C-1)=sC-2…,F(s2,2)=s1,由此可以计算得到临界波段向量s=(s0,s1,…,sC)T。
以临界波段为边界将高光谱图像划分为C个连续波段的聚类。C为大于0的正整数,一般小于最终所要求选择的波段数。
步骤3:根据得到的聚类结果,在每个波段聚类中,首先计算该波段聚类所有波段向量的平均向量,即平均波段,再选择出该波段聚类中与平均波段的欧式距离最近的波段,加入已选波段集合。
步骤4:迭代地在每个波段聚类内对尚未被选择的波段进行选择,其过程如下:
步骤a:在第k次迭代过程中,对任一个尚未被选择的波段向量xi,设其所在的聚类为编号为j,设Bj,k为第k次迭代下,聚类j内已经被选择的波段向量构成的矩阵。设为在第k次迭代下,聚类j内所有已选波段向量加上波段向量xi构成的矩阵。再设为在第k次迭代下,聚类j内所有尚未被选择的波段扣除波段xi得到的波段集合,并设该集合内所有波段向量的平均波段为
步骤b:分别用矩阵Bj,k和对进行线性重构,得到两个重构误差:
步骤c:对所有尚未选择的波段i,计算两个重构误差的比值使得重构误差最小的波段,即为本次迭代选择的波段,将其加入已选波段。
当已选波段数量达到算法设定的数量时,停止上述迭代过程,完成高光谱图像波段选择。
本实施例是在中央处理器为i5-4590 3.30GHz CPU、内存16G、WINDOWS 7操作系统上,运用MATLAB软件进行的仿真。仿真中使用的数据为公开数据集Indian Pines。
图2为采用不同方法进行高光谱图像波段选择后,再利用SVM(支持向量机)分类器进行分类的分类精度结果图。其中,CBBS-KL(Clustering Based Band Selection–KL)为基于聚类的波段选择-KL散度方法,CBBS-MI(Clustering Based Band Selection–MI)为基于聚类的波段选择-互信息方法,GOC(Global Optimal Clustering)为全局最优聚类方法,即本发明方法,VGBS(Volume Gradient Band Selection)为体积梯度波段选择方法,UBS(Uniformed Band Selection)为等距波段选择方法,MTSP(Multi-task Sparsitypursuit)为多任务稀疏追踪方法,Total Bands为全波段。
可以看出,大部分情况下,利用本发明的方法进行波段选择后都比其他方法具有更高的分类精度。
Claims (1)
1.一种基于全局最优聚类的高光谱图像波段选择方法,其特征在于包括以下步骤:
步骤1:按照pn=(po-M2)/(M1-M2)对高光谱图像进行归一化处理,其中,pn为归一化处理后图像的像素值,po为原高光谱图像的像素值,M1为原高光谱图像中像素值的最大值,M2为原高光谱图像中像素值的最小值;
步骤2:以聚类间差距和聚类内差距的比值为目标函数,使用二分法和动态规划方法优化该目标函数,得到一个全局最优的临界波段向量,以这些临界波段为边界将高光谱图像划分为C个连续波段的聚类;
其中,聚类内差距为聚类间差距为xj为高光谱图像第j个波段的波段向量,c为函数的参数,表示将高光谱图像的前sc个波段划分为c个聚类,0<c<C,C为待生成的聚类个数,为小于最终所要求选择波段数的正整数,si表示波段序号,i=0,…,c,满足0=s0<s1<...<sc≤L,L为高光谱图像的波段数,M(a,b)表示归一化后高光谱图像从a波段到b波段区间内所有波段的平均波段向量,如M(1,L)即表示归一化后高光谱图像所有波段的平均波段向量;
步骤3:对每个聚类,选择一个与类中心欧式距离最近的一个波段为该聚类的一个代表波段,所有代表波段构成已选波段集合;所述的类中心是指该聚类内所有波段的平均波段;
步骤4:在每个聚类内迭代地对剩余波段进行选择,直到所选择波段数量满足要求;所述的迭代地对剩余波段进行选择的具体过程为:在第k次迭代中,k>0,对任一个尚未被选择的波段向量xi,设其所在的波段聚类为j,Bj,k为第k次迭代中波段聚类j内已经被选择的波段向量构成的矩阵,为在第k次迭代中聚类j内所有已选波段向量加上波段向量xi构成的矩阵,设为在第k次迭代中聚类j内所有尚未被选择的波段扣除波段xi得到的波段集合,并设该集合内所有波段向量的平均波段为分别用矩阵Bj,k和对进行线性重构,以两个重构误差的比值作为对波段向量xi的评价值,计算所有尚未被选择波段的评价值,并将评价值最高的波段加入已选波段集合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710342464.4A CN107220662B (zh) | 2017-05-16 | 2017-05-16 | 基于全局最优聚类的高光谱图像波段选择方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710342464.4A CN107220662B (zh) | 2017-05-16 | 2017-05-16 | 基于全局最优聚类的高光谱图像波段选择方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107220662A CN107220662A (zh) | 2017-09-29 |
CN107220662B true CN107220662B (zh) | 2019-07-19 |
Family
ID=59944002
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710342464.4A Active CN107220662B (zh) | 2017-05-16 | 2017-05-16 | 基于全局最优聚类的高光谱图像波段选择方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107220662B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107832793B (zh) * | 2017-11-08 | 2021-08-06 | 深圳大学 | 一种高光谱图像的分类方法及系统 |
CN108154094B (zh) * | 2017-12-14 | 2020-04-24 | 浙江工业大学 | 基于子区间划分的高光谱图像非监督波段选择方法 |
CN108491888B (zh) * | 2018-04-02 | 2020-07-14 | 安徽大学 | 基于形态学分析的环境监测高光谱数据谱段选择方法 |
CN108805061A (zh) * | 2018-05-30 | 2018-11-13 | 西北工业大学 | 基于局部自适应判别分析的高光谱图像分类方法 |
CN109460712B (zh) * | 2018-10-12 | 2024-01-12 | 西北工业大学 | 一种基于最优近邻重构的高光谱波段选择方法 |
CN109871768B (zh) * | 2019-01-18 | 2022-04-29 | 西北工业大学 | 基于共享最近邻的高光谱最优波段选择方法 |
CN112330622B (zh) * | 2020-10-30 | 2022-11-11 | 中国电子科技集团公司第五十四研究所 | 一种基于地物最大区分度的高光谱图像波段选择方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102903114A (zh) * | 2012-10-09 | 2013-01-30 | 河海大学 | 一种基于改进型层次聚类的高光谱遥感数据降维方法 |
CN102930533A (zh) * | 2012-10-09 | 2013-02-13 | 河海大学 | 一种基于改进k-均值聚类的半监督高光谱影像降维方法 |
CN103020912A (zh) * | 2012-12-18 | 2013-04-03 | 武汉大学 | 一种结合波段聚类和稀疏表达的遥感影像复原方法 |
CN105046276A (zh) * | 2015-07-14 | 2015-11-11 | 西安电子科技大学 | 基于低秩表示的高光谱图像波段选择方法 |
CN105512675A (zh) * | 2015-11-27 | 2016-04-20 | 中国石油大学(华东) | 一种基于记忆性多点交叉引力搜索的特征选择方法 |
-
2017
- 2017-05-16 CN CN201710342464.4A patent/CN107220662B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102903114A (zh) * | 2012-10-09 | 2013-01-30 | 河海大学 | 一种基于改进型层次聚类的高光谱遥感数据降维方法 |
CN102930533A (zh) * | 2012-10-09 | 2013-02-13 | 河海大学 | 一种基于改进k-均值聚类的半监督高光谱影像降维方法 |
CN103020912A (zh) * | 2012-12-18 | 2013-04-03 | 武汉大学 | 一种结合波段聚类和稀疏表达的遥感影像复原方法 |
CN105046276A (zh) * | 2015-07-14 | 2015-11-11 | 西安电子科技大学 | 基于低秩表示的高光谱图像波段选择方法 |
CN105512675A (zh) * | 2015-11-27 | 2016-04-20 | 中国石油大学(华东) | 一种基于记忆性多点交叉引力搜索的特征选择方法 |
Non-Patent Citations (4)
Title |
---|
A Novel Ranking-Based Clustering Approach for Hyperspectral Band Selection;Sen Jia 等;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;20160131;第54卷(第1期);88-102 |
Dual-Clustering-Based Hyperspectral Band Selection by Contextual Analysis;Yuan Yuan 等;《IEEE Transactions on Geoscience and Remote Sensing》;20160331;第54卷(第3期);1431-1445 |
Salient Band Selection for Hyperspectral Image Classification via Manifold Ranking;Qi Wang 等;《IEEE TRANSACTIONS ON NEURAL NETWORKS AND LEARNING SYSTEMS》;20160630;第27卷(第6期);1279-1289 |
基于谱聚类与类间可分性因子的高光谱波段选择;秦方普 等;《光谱学与光谱分析》;20150531;第35卷(第5期);1357-1364 |
Also Published As
Publication number | Publication date |
---|---|
CN107220662A (zh) | 2017-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107220662B (zh) | 基于全局最优聚类的高光谱图像波段选择方法 | |
US11341424B2 (en) | Method, apparatus and system for estimating causality among observed variables | |
Yuan et al. | Adaptive reverse graph learning for robust subspace learning | |
Wang et al. | CLUES: A non-parametric clustering method based on local shrinking | |
Yang et al. | A sparse singular value decomposition method for high-dimensional data | |
Shen et al. | Pcw-net: Pyramid combination and warping cost volume for stereo matching | |
Liu et al. | A weighted Lq adaptive least squares support vector machine classifiers–Robust and sparse approximation | |
CN108171012B (zh) | 一种基因分类方法与装置 | |
Zhu et al. | Multi-modal AD classification via self-paced latent correlation analysis | |
Kwon et al. | Tuning parameter selection for the adaptive Lasso in the autoregressive model | |
JPWO2015190593A1 (ja) | 情報処理方法、情報処理装置およびそのプログラム | |
Harris et al. | Feature clustering and ranking for selecting stable features from high dimensional remotely sensed data | |
Marjanovic et al. | On ${l} _ {q} $ optimization and sparse inverse covariance selection | |
CN109145975A (zh) | 一种矩阵范数导出正则化的多视图聚类机器学习方法 | |
Majumdar et al. | Improved group sparse classifier | |
Hoeltgen et al. | Clustering-based quantisation for PDE-based image compression | |
CN112800138B (zh) | 大数据分类方法及系统 | |
CN116343032A (zh) | 结合高斯回归混合模型和mrf高光谱函数数据的分类方法 | |
Husmeier et al. | Inter-time segment information sharing for non-homogeneous dynamic Bayesian networks | |
CN112148942A (zh) | 基于数据聚类的业务指标数据分类方法及装置 | |
Li et al. | A projection matrix design method for MSE reduction in adaptive compressive sensing | |
Yang et al. | High-dimensional trimmed estimators: A general framework for robust structured estimation | |
Wackers et al. | Error estimation for surrogate models with noisy small-sized training sets | |
CN114550906A (zh) | 基于多视图鲁棒表示的癌症亚型识别系统 | |
Shang et al. | A dynamic local cluster ratio-based band selection algorithm for hyperspectral images |
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 |