基于波段迁移的高光谱图像聚类方法
【专利摘要】本发明公开了一种基于波段迁移的高光谱图像聚类方法。主要解决现有基于高斯混合模型的聚类方法对于高光谱图像聚类不准确，速度较慢，稳定性较差的问题。其实现步骤是：先从高光谱图像中随机选出一些波段分别作为目标域，根据聚类一致性值，从目标域选出一些聚类不稳定样本；再根据欧氏距离原则从高光谱图像剩余波段中选取一些波段作为源域，利用EM算法估计源域中高斯混合模型的参数；最后将源域中的参数迁移至目标域来改善目标域中不稳定样本的聚类结果。本发明与传统的聚类方法相比，具有速度快、效果好等优点，而且可使用较少的波段达到较高的分类精度，大大的节约了计算成本，可用于矿产探测、资源调查及环境监测。
【专利说明】基于波段迁移的高光谱图像聚类方法
【技术领域】
[0001]本发明属于图像处理【技术领域】，涉及对高光谱遥感图像的处理，可用于矿产探测、资源调查及环境监测。
【背景技术】
[0002]伴随着遥感技术的发展和成像设备的进步，光学遥感技术进入了高光谱遥感时代，高光谱遥感技术的出现和发展将使得遥感成像设备以更快的速度、更高的分辨率、更大的信息量为人类提供海量的观测数据，将人们通过遥感技术观测和认识地物的能力带入了一个崭新的阶段。针对由成像光谱仪采集得到的高光谱遥感数据，如何对其中蕴含的海量信息进行充分而且有效的利用，给广大的科研工作者和技术人员提出的挑战并带来了机遇。通过对高光谱遥感数据进行有效的分类，能够更加深入的挖掘高光谱遥感图像在矿产探测、资源调查以及环境监测等方面的应用。由于实际应用中，有标签的数据难以获得并且代价巨大，于是无监督的分类在高光谱遥感图像中的应用显得尤为重要。
[0003]聚类，是将物理或抽象对象的集合组成为由类似的对象组成的多个类的过程。聚类是数据挖掘技术中的重要组成部分，根据样本间的相似度的度量标准将数据分成若干个簇，所组成的簇是一组数据对象的集合，且使同一簇中的样本相似度尽量高，不属于同一簇中的样本之间的相似度尽量低。在聚类的过程中，一些簇中的数据对象可以被看作为一个整体。一个聚类分析系统的输入是一组样本和一个度量样本间相似度或距离的标准，而输出则是簇集，即数据集的几个类，这些类构成一个分区或者分区结构。聚类是研究数据间逻辑上或物理上相互关系的技术，其分析结果不仅可以揭示数据间的内在联系与区别，还可以为进一步的数据分析与知识发现提供重要依据。作为统计学的重要研究内容之一，聚类分析具有坚实的理论基础，并形成了系统的方法学体系。
[0004]众所周知，高光谱图像数据波段服从高斯分布，用高斯混合模型对高光谱图像进行建模会得到良好的效果。高斯混合模型GMM是一个常用的描述混合密度分布的模型，即多个高斯分布的混合分布。高斯混合模型是一种半参数的密度估计方法，它融合了参数估计法和非参数估计法的优点，既不局限于特定的概率密度函数的形式，而且模型的复杂度仅与所求解的问题有关而与样本集合的大小无关。高斯混合模型的一个重要特性是，如果模型中的成员足够多，它能够以任意精度逼近任意的连续分布。所以高斯混合模型在许多领域得到了广泛的应用。高斯混合模型是一种统计模型，它能够很好地刻画参数空间中数据的空间分布及其特性。高斯混合模型中常用EM算法估计它的参数。然而，单纯的用基于高斯混合模型的EM算法对高光谱进行分类会有对初始值敏感，容易陷入局部最优的问题，并且由于高光谱中波段数目多，使得聚类速度缓慢。
【发明内容】
[0005]本发明的目的在于针对上述高光谱遥感图像中聚类的缺陷，提出一种基于高光谱图像波段迁移的聚类方法，以加快聚类速度，提高聚类精度。[0006]为实现上述目的，本发明包括如下步骤:
[0007](I)输入高光谱图像1_，1?是不同谱对应的11^11大小的图像』为高光谱图像的波段数目；
[0008](2)随机选出高光谱的10个波段数据作为目标域样本Dt，根据这10个波段数据，利用欧氏距离作为度量准则，从剩下的波段中再另选取出10个波段数据作为源域样本Ds ；
[0009](3)对目标域样本Dt用EM、K_means和FCM聚类算法分别聚3次，得到的9个聚类标签，计算这9个聚类标签的聚类一致性值Cl ；
[0010](4)设定一个阈值Q = 0.6，将聚类一致性值Cl与阈值Q进行比较，若Cl大于等于Q，则将样本标记为确定性样本Da，否则标记为非确定性样本Db，同时记录确定性样本Da的聚类标签Ia ；
[0011](5)对源域样本Ds的高斯混合模型用EM算法进行参数估计，得到源域中混合高斯模型的类混合权值参数a”均值μP协方差σP并将得到的这些参数迁移至目标域样本Dt，得到目标域样本Dt新的类混合权值参数α/、均值μ/、协方差；
[0012](6)根据新的类混合权值参数α /、均值μ/、协方差σ /，计算出非确定样本Db的聚类标签Ib ;
[0013](7)根据确定性 样本Da的聚类标签Ia和非确定性样本Db的聚类标签lb，得到目标域样本Dt的聚类标签1:1 = Ia U Ib ；
[0014]本发明与现有的技术相比具有以下优点:
[0015]1、本发明通过引入迁移学习的方法，充分利用了高光谱图像波段间的信息，减少了需要使用的波段数目，使聚类速度加快；同时通过引入其它波段的分布信息，提高了聚类的准确性；
[0016]2、本发明针对EM算法中对于初始值敏感，容易陷入局部最优等问题，引入聚类一致性值和迁移学习的方法，对聚类一些不稳定的样本重新进行聚类，大大的提高了聚类的精度和稳定性。
【专利附图】
【附图说明】
[0017]图1是本发明的总流程图；
[0018]图2是本发明所使用的高光谱遥感图像的灰度图；
[0019]图3是本发明中所使用高光谱图像的真实分布标签；
[0020]图4是本发明中对高光谱数据用GMM-EM聚类方法得到的聚类结果；
[0021]图5是本发明中对高光谱数据用基于MIFS算法聚类之后的结果；
[0022]图6是本发明中对高光谱数据用基于PCA算法聚类之后的结果；
[0023]图7是本发明中用基于波段迁移之后目标域的聚类结果。
【具体实施方式】
[0024]参照图1，本发明基于波段迁移的高光谱图像聚类方法，包括如下步骤:
[0025]步骤1:输入高光谱图像IKXN，尺是不同谱对应的mXn大小的图像，N为高光谱图像的波段数目，这幅高光谱遥感图像是美国加利福尼亚萨利纳斯山谷的AVIRIS遥感图像，大小为86X83，灰度图像如图2所示，该图总共有6类地物，真实分布标签如图3所示。[0026]步骤2:随机选出高光谱的10个波段数据作为目标域样本Dt，根据这10个波段数据，利用欧氏距离作为度量准则，从剩下的波段中再另选取出10个波段数据作为源域样本Ds。
[0027]步骤3:对目标域样本Dt用EM、K_means和FCM聚类算法分别聚3次，得到的9个聚类标签，计算这9个聚类标签的聚类一致性值Cl。
[0028]3a)利用EM算法对于目标域样本Dt进行3次聚类得到3个类别标签，再利用K-means算法对目标域样本Dt进行3次聚类得到3个类别标签，最后再利用FCM算法对目标域样本Dt进行3次聚类得到3个类别标签，总共得到9个类别标签；
[0029]3b)通过如下公式计算9个聚类标签的聚类一致性值Cl:
[0030]Cl = {Cl (X)}，
[0031]其中，
【权利要求】
1.一种基于波段迁移的高光谱图像聚类方法，包括如下步骤: (1)输入高光谱图像IKXN，R是不同谱对应的HiXn大小的图像，N为高光谱图像的波段数目； (2)随机选出高光谱的10个波段数据作为目标域样本Dt，根据这10个波段数据，利用欧氏距离作为度量准则，从剩下的波段中再另选取出10个波段数据作为源域样本Ds ； (3)对目标域样本Dt用EM、K-means和FCM聚类算法分别聚3次，得到的9个聚类标签，计算这9个聚类标签的聚类一致性值Cl ； (4)设定一个阈值Q= 0.6，将聚类一致性值Cl与阈值Q进行比较，若Cl大于等于Q，则将样本标记为确定性样本Da，否则标记为非确定性样本Db，同时记录确定性样本Da的聚类标签Ia ; (5)对源域样本Ds的高斯混合模型用EM算法进行参数估计，得到源域中混合高斯模型的类混合权值参数a”均值μ P协方差σ P并将得到的这些参数迁移至目标域样本Dt，得到目标域样本Dt新的类混合权值参数α/、均值μ/、协方差σ/ ； (6)根据新的类混合权值参数α/、均值μ/、协方差σ/，计算出非确定样本Db的聚类标签Ib ; (7)根据确定性样本Da的聚类标签Ia和非确定性样本Db的聚类标签lb，得到目标域样^Dt的聚类标签I:1 = Ia U lbo
2.根据权利要求1所述的方法，其中所述步骤(3)中计算9个聚类标签的聚类一致性值Cl，通过如下公式计算:
3.根据权利要求1所述的方法，其中步骤(5)所述的将得到的参数迁移至目标域样本Dt，得到目标域样本Dt新的类混合权值参数a/、均值μ/、协方差，按如下步骤进行: (5a)根据类混合权值参数α P均值μ j和协方差σ ρ计算源域中的隶属度Zi^j:
4.根据权利要求1所述的方法，其中步骤(6)所述的根据新的类混合权值参数a/、均值μ/、协方差σ/，计算出非确定样本Db的聚类标签lb，按如下步骤进行: (6a)计算目标域样本的隶属度Zi,/:
【文档编号】G06K9/62GK103745232SQ201410032062
【公开日】2014年4月23日 申请日期:2014年1月23日 优先权日:2014年1月23日
【发明者】缑水平, 刘芳, 张观侣, 马晶晶, 马文萍, 侯彪 申请人:西安电子科技大学