一种基于无监督宽度学习的高光谱图像聚类方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN109685105A(43)申请公布日2019.04.26(21)申请号201811366436.7(22)申请日2018.11.16(71)申请人中国矿业大学地址221116江苏省徐州市大学路1号中国矿业大学科研院(72)发明人王雪松孔毅程玉虎(74)专利代理机构南京瑞弘专利商标事务所(普通合伙)32249代理人吴旭(51)Int.Cl.G06K9/62(2006.01)G06K9/00(2006.01)权利要求书3页说明书5页附图1页(54)发明名称一种基于无监督宽度学习的高光谱图像聚类方法(57)摘要本发明的一种基于无监督宽度学习的高光谱图像聚类方法，首先，为保留原始高光谱图像的本征流形结，利用图正则稀疏自动编码器微调输入到映射特征和映射特征到增量节点的权重；接着，无监督宽度学习的模型目标函数由图正则项和连接权重的l2范数项构成，并通过求解广义特征值分解问题，得到无监督宽度学习的模型权重以及无监督宽度学习的输出；最后，利用谱聚类得到聚类结果。CN109685105ACN109685105A权利要求书1/3页1.一种基于无监督宽度学习的高光谱图像聚类方法，其特征在于，包括如下步骤：步骤1：利用层次导向滤波对原始的高光谱图像进行滤波处理，得到高光谱图像的光谱-空间表示；步骤2：高光谱图像的光谱-空间表示作为输入，利用图正则稀疏自动编码器优化后的MEME权重Wi和W，先将输入通过Wi映射到映射特征，然后映射特征通过W映射到增量节点中；步骤3：建立由图正则项和连接权重的二范数正则项构成的无监督宽度学习的模型目标函数，求解模型目标函数，得到模型的输出层权重Wm；步骤4：根据输出层权重Wm，计算得到无监督宽度学习的输出向量步骤5：通过对输出向量进行谱聚类，得到高光谱图像的聚类结果。2.根据权利要求1所述的基于无监督宽度学习的高光谱图像聚类方法，其特征在于：所述步骤1中，利用层次导向滤波对原始的高光谱图像进行滤波处理的过程表示为：式中，为关于和的能量函数，和为基于输入高光谱图像和导向图像G的线性系数，上标p指第p个波段；ωv是像素v和近邻像素构成的大小为(2r+1)×(2r+1)的窗口，r为窗口半径，v是ωv中的一个像素；ε为控制参数；Gv为导向图像G的第v个像素；为输入高光谱图像的第v个像素的第p个波段；式(1)的解为：式中，|ω|为ωv中像素的个数，μv和σv分别表示导向图像G的均值和标准差，为ωv中高光谱图像的均值，为第p个波段的G的均值。3.根据权利要求1所述的基于无监督宽度学习的高光谱图像聚类方法，其特征在于：所n×m述步骤2中，给定高光谱图像的光谱-空间表示X＝{x1,…,xn}∈R，其中n表示样本个数，m表示样本维数，Rn×m表示n×m维空间；M假设所需的权重和特征分别为Wi和Zi，优化如下目标函数：式中，λ表示稀疏正则参数；根据流形假设，原始数据空间中相近邻的两个数据点xi和xj，在映射后的特征zi和zj也应相互靠近，则式(4)进一步写为：2CN109685105A权利要求书2/3页式中，α表示图正则稀疏自动编码器的图正则参数，tr(·)表示迹操作，L为拉普拉斯矩阵，wij为特征zi和zj的连接权重；n表示训练样本的个数；式(5)通过交替方向乘子法求解，由于l1范数||·||1为非凸函数，引入辅助变量O，则将式(5)写为：式(6)的拉格朗日表达式为：式中，ρ>0为一个常数，u为引入的拉格朗日乘子矩阵；M根据交替方向乘子法优化方法，Wi、O和u交替优化，同时固定其它两个变量，直到收敛M或者达到预先定义的迭代次数，从而得到所需的权重Wi；其中，所述优化方法包括如下步骤：MM步骤2.1：更新Wi，Wi的更新过程通过求解式(8)得到：M(k+1)M式中，Wi为k+1时刻的Wi，J(·)为代价函数；M通过计算J关于Wi的导数，并使导数为0，得到：(k)(k)式中，Zi为映射特征中的第i组映射特征，O为k时刻的O，u为k时刻的u，I为单位矩阵；步骤2.2：更新O，利用软阈值操作：式中，Sκ(g)为关于变量g的软阈值函数，g为自变量，κ为预先定义的阈值；则O的更新过程为：(k+1)M,(k+1)(k)O＝Sλ/ρ(Wi+u)(11)(k+1)式中，O为k+1时刻的O，Sλ/ρ为在参数λ和ρ下的软阈值操作；步骤2.3：更新u，u的更新表达式为：(k+1)(k)M,(k+1)(k+1)u＝u+(Wi-O)(12)3CN109685105A权利要求书3/3页式中，u(k+1)为k+1时刻的u；通过式(13)计算得到Zi：MZi＝XWi(13)从映射特征到增量节点的权重也利用图正则稀疏自动编码器进行微调，则所需的微调后的权重为WE，进而增量节点通过式(14)计算得