文章快速检索  
  高级检索
空-谱协同正则化稀疏超图嵌入的高光谱图像分类
黄鸿, 陈美利, 王丽华, 李政英     
重庆大学光电技术与系统教育部重点实验室, 重庆 400044
摘要:传统依据图嵌入的高光谱图像维数约简算法多数仅利用光谱信息表征像元间单一关系,忽视了数据间的多元几何结构。本文提出了一种面向高光谱图像分类的空-谱协同正则化稀疏超图嵌入算法(SSRSHE)。该算法首先利用稀疏表示揭示像元之间的相关性,自适应选择近邻,并构建稀疏本征超图和惩罚超图,以有效表征像元间的复杂多元关系,并进行正则化处理。然后利用遥感图像空间一致性原则,计算局部空间邻域散度来保持样本局部邻域结构,并引入样本总体散度来保持高光谱数据的整体结构。在低维嵌入空间中,尽可能使类内数据聚集、类间数据远离,提取鉴别特征用于分类。在Indian Pines和PaviaU高光谱遥感数据集上试验结果表明,本文算法总体分类精度分别达到86.7%和92.2%。相比传统光谱维数约简算法,该算法可有效改善高光谱图像地物分类性能。
关键词高光谱图像    维数约简    正则化稀疏超图模型    空-谱联合    分类    
Using spatial-spectral regularized hypergraph embedding for hyperspectral image classification
HUANG Hong, CHEN Meili, WANG Lihua, LI Zhengying     
Key Laboratory of Optoelectronic Technique System of the Ministry of Education, Chongqing University, Chongqing 400044, China
Abstract: In recent years, many graph embedding methods were developed for dimensionality reduction (DR) of hyperspectral image (HSI), while these methods only use spectral information to reveal a simple intrinsic relation and ignore complex spatial-spectral structure in HSI. A new DR method termed spatial-spectral regularized sparse hypergraph embedding (SSRSHE) is proposed for the HSI classification. SSRSHE explores sparse coefficients to adaptively select neighbors for constructing the regularized sparse intrinsic hypergraph and the regularized sparse penalty hypergraph. Based on the spatial consistency property of HSI, a local spatial neighborhood scatter is computed to preserve local structure, and a total scatter is computed for global structure of HSI. Then, the optimal discriminant projection is obtained by possessing better intrinsic data compactness and penalty pixels separability, which is beneficial for classification. The experimental results on Indian Pines and PaviaU hyperspectral data sets show that the overall classification accuracies respectively reach 86.7% and 92.2%. The proposed SSRSHE method can effectively improve classification performance compared with the traditional spectral DR algorithms.
Key words: hyperspectral image    dimensionality reduction    regularized sparse hypergraph    spatial-spectral features    image classification    

高光谱遥感图像通过从可见光到短波红外区域的密集光谱采样,可在数百个窄而连续的相邻光谱波段中提供空间场景,包含了丰富的空间、辐射和光谱信息,为地物精细分类提供了强有力的探测手段,目前已广泛应用于矿物勘探、环境监测、精准农业和目标识别等领域[1-2]。然而,高光谱数据具有数据量大、波段数多、波段间相关性强等特点,传统方法易导致“维数灾难”问题[3-4]。因此,如何减少波段数且尽量保留有用信息已成为高光谱遥感领域的研究热点问题。

维数约简是克服数据冗余的有效方法,可在降低数据维数的同时尽可能保留数据中的本征信息[5]。目前学者们提出了一系列的维数约简方法,如主成分分析(principal component analysis, PCA)[6]、线性判别分析(linear discriminant analysis, LDA)[7]、等距映射(isometric feature mapping, ISOMAP)[8]、局部线性嵌入(local linear embedding, LLE)[9]、邻域保持嵌入(neighborhood preserving embedding, NPE)[10]、拉普拉斯等距离映射(Laplacian eigenmaps, LE)[11]及局部保持投影(locality preserving projection, LPP)[12]。上述方法可统一在图嵌入框架(graph embedding, GE)[13-14]下,其异于如何定义本征图和惩罚图,但都为非监督方法,其分类性能受限。针对此问题,学者们通过将样本先验知识引入到图嵌入框架来改善分类性能,提出了边缘Fisher分析(marginal Fisher analysis, MFA)[15]和正则化局部判别嵌入(regularized local discriminant embedding, RLDE)[16]等监督学习方法,以提升分类精度。

然而,直接图嵌入方法只考虑数据间一元关系,在实际应用中高维数据通常具有复杂的多元几何结构[17-18]。为表征高维数据中的复杂结构,学者们试图引入超图学习来表示高光谱数据间的高阶关系。文献[19]提出了一种判别超-拉普拉斯投影(discriminant hyper-Laplacian projections, DHLP)方法,通过构造超图来获得超-拉普拉斯矩阵,实现维数约简。文献[20]提出了一种超图拉普拉斯联合稀疏化处理方法来分析像元的内在关系,以提取低维特征进行分类。

上述方法仅利用了样本的光谱信息,却忽略了像元之间的空间位置关系,而研究表明空-谱联合维数约简方法可明显提高地物分类性能。文献[21]通过空间自适应方法提取影像的空间特征和光谱特征,在分类精度和计算效率上均取得了较好效果。文献[22]提出了一种空-谱协同嵌入方法(spatial-spectral coordination embedding, SSCE),利用样本空间块替代单个样本度量数据间相似性,降低异类地物被选为近邻的概率,从而改善地物分类效果。与此同时,空间信息也被引入超图模型中,文献[17]提出了一种融合空-谱信息的超图嵌入方法,利用像元空间邻域构造超边,能有效提取低维特征,但忽视了像元的类别信息。在文献[23]中,通过像元波段选取提取扩展形态学特征,并与光谱信息融合来构建超图模型,提取嵌入特征以提升地物分类性能。上述空-谱联合维数约简方法,或是忽略了像元间多元几何结构关系,或是在构造超图模型时没有充分利用样本标签信息,限制了分类性能的进一步提升。

针对上述问题,本文提出了一种空-谱协同正则化稀疏超图嵌入方法(spatial-spectral regularized sparse hypergraph embedding, SSRSHE)。该方法运用稀疏系数自适应揭示数据间近邻关系,并结合类别信息构建正则化稀疏超图,从而有效表征高光谱数据的多元几何结构。同时,融入图像的空间信息,构造局部空间邻域散度来表征样本局部邻域结构,同时定义样本总体散度矩阵来保证数据全局信息,提取有效鉴别特征,实现维数约简。在Indian Pines和PaviaU高光谱数据集上验证了本文算法的有效性。

1 本文算法

假设文中高光谱数据集Z=[z1, z2, …, zi, …, zn]∈Rd×n,其中d为波段数,n为样本数,类别标签集L=[l1, l2, …, li, …, ln],li∈{1, 2, …, u},其中u为样本类别数。低维嵌入特征可表示为Y=PTZ, YRτ×nτ(τ<<d)为嵌入维数,PRd×τ为映射矩阵。

1.1 图嵌入学习

为更好地理解维数约简算法,学者们提出了一种图嵌入框架(GE)来表示数据几何结构,并将PCA、LDA、ISOMAP、LLE、LE、NPE及LPP等算法统一到该框架中。在图嵌入框架下,需构建本征图和惩罚图两个无向图。本征图GI(V, WI)表征数据中需要保持的统计或几何性质,惩罚图GP(V, WP)描述数据中应避免的某种特性,其中V为顶点集,WIWP分别为图GIGP的权重矩阵,可通过简单法或热核函数来定义。

图嵌入框架意在低维空间中保留数据集的某些统计或几何属性,其低维嵌入特征可通过优化以下目标函数得到

(1)

式中,DI是对角矩阵,且DiiI=∑jwijILI=DI-WI为本征图GI的拉普拉斯矩阵; C为一常量矩阵;H为约束矩阵,可为单位阵以实现归一化处理,或为惩罚图GP的拉普拉斯矩阵,即H=LP=DP-WP, DiiP=∑jwijP

1.2 超图模型

直接图嵌入模型仅考虑了两点间一阶关系,而超图模型能有效表征数据间的多元特性[17]。超图模型可表示为GH=(VH, EH, WH),其中VH表示顶点集,EH为超边集,对应的相似权重矩阵是WH,以度量超边内各顶点间相关性。

为表示GH的内在关系,假设每一超边ei含有N(ei)个顶点,其权重表示为w(ei)∈EH,则关联矩阵H=[Hmn:h(em, vn)]∈R|EH|×|VH|、超边em的度d(em)和顶点vn的度d(vn)可分别定义为

(2)
(3)
(4)

式中,,即值等于超边内目标像素与其近邻点权重之和;h为热核系数。

综上,超图内每一超边由某一像元与其近邻点构成,揭示数据间内在多元关系。其对应的关联矩阵H,每行中的非零元素,描述每一超边内各点分布情况。超图通过多对顶点连通以表征邻域内顶点间多元结构,因而可更好地描述数据中多元关系。

1.3 SSRSHE算法

为表征高光谱数据中的多元几何结构关系,并联合像元的空间-光谱信息,本文提出了一种空-谱协同正则化稀疏超图嵌入(SSRSHE)方法。首先利用样本的稀疏系数来自适应性选择其近邻,构建稀疏本征超图和惩罚超图来揭示高光谱数据间的多元结构。同时,依据空间一致性原理构造局部空间邻域散度以保持像元局部空间近邻关系, 并采用样本总体散度来表征高光谱数据整体特性。在低维鉴别空间中,使类内数据尽可能聚集、类间数据尽可能发散,提取鉴别特征,提升地物分类性能。该算法的具体流程如图 1所示。

图 1 SSRSHE算法流程 Fig. 1 Flowchart of the proposed SSRSHE method

1.3.1 正则化稀疏超图模型构建

在构建超图时,首先需要选择合适的样本近邻点。目前的欧氏距离度量方法存在近邻点选取不准确及参数难以确定等问题,而稀疏表示具有自然鉴别力能自适应地揭示出数据的内在关系。某个样本可以由一个足够大的样本空间来近似线性表示,且表示系数大部分为零,只有极少数与该样本同类别数据对应的系数为非零,因此可反映数据的本征属性。

基于此,本文提出了一种正则化稀疏超图模型,首先通过稀疏表示[24]得到数据的稀疏系数矩阵,揭示数据内在关联特性,以自适应获取像元近邻。稀疏系数可通过以下l1范数求解

(5)

式中,ε为稀疏误差;E是全为1的向量。在具体计算中,可通过将式(5)问题进一步转化为Lasso问题求解[25],即可得到稀疏系数矩阵S=[s1, s2, …, sn]T

图 2为基于稀疏系数自适应选取近邻构造超边示意图。因稀疏系数可反映数据间相似性,对应系数非零则表示像元间具有相关性,其值越大则属于同类近邻点可能性越大。因此相比欧氏度量,利用稀疏系数自适应选择近邻能更为有效反映数据内蕴信息。

图 2 基于稀疏系数的自适应选取近邻构造超边 Fig. 2 Construction of sparse hyperedge

根据样本稀疏特性和类别信息,构建稀疏本征超图Gw={Z, EHw, WHw},揭示类内数据的本征特性;同时,构建稀疏惩罚超图Gb={Z, EHb, WHb},表征异类数据间几何结构;其中Z为顶点集,EHwEHb分别是本征超边集和惩罚超边集,对应权重矩阵为WHwWHb。在稀疏本征超图Gw中,根据稀疏系数矩阵具有的自适应表征能力,本征超边eiw由样本zi与其对应稀疏系数非零点连接而成。则其超边权值wijw可定义为

(6)

式中,sij为稀疏系数矩阵S中第i行第j列元素,以表征像元zizj互相关程度;系数α(α>1)用于调节同类像元权重,强化同类样本贡献率,进而提升分类性能。由此,本征超边eiw的相似权重w(eiw)为

(7)

式中,N(eiw)表示本征超边eiw包含的顶点数目。

同时,稀疏本征超图的关联矩阵Hw=[Hijw:h(ejw, zi)]∈R|Ew|×|Z|计算如下

(8)

式中,。根据Hww(eiw),计算像元zi与超边eiw的度

(9)
(10)

同理,在稀疏惩罚超图Gb中,惩罚超边eib由样本zi与其对应稀疏系数非零且类别不同的点连接而成,其边权值可定义为

(11)

由式(11)可获取惩罚超边eib的权重w(eib),表述了zi与其不同类近邻样本间的相似特性

(12)

式中,N(eib)表示惩罚超边eib含有的顶点数。

其对应的关联矩阵Hb=[Hijb:h(ejb, zi)]∈R|Eb|×|Z|可表示为

(13)

由式(13)求解顶点zi和惩罚超边eib的度分别是

(14)
(15)

在低维映射空间中,为提取鉴别特征,应使同类数据尽可能聚集、不同类数据尽可能远离,因此目标函数可表示为

(16)
(17)

式中,Lw=DzwHwWw(Dew)-1(Hw)TLb=DzbHbWb(Deb)-1(Hb)T分别是GwGb对应的拉普拉斯矩阵;DzwDzb分别为θwθb对角阵;Dew表示本征超边度对角矩阵,Deb代表惩罚超边度对角矩阵。稀疏本征超图保持了同类样本的聚集性,而稀疏惩罚超图避免了非同类样本在低维嵌入空间内过于接近。

由式(16)、式(17)中的目标函数,可进一步转化为以下最优化问题

(18)

在训练样本较少的情况下,式(18)易受奇异点影响。故在此引入正则化项,则式(18)中的优化问题可拓展为

(19)

式中,η(0 < η < 1)表示正则化参数。正则项ZZT用于保持样本的多样性。将ZLwZT对角化,以改善式(19)问题求解的稳定性,即其对应的特征值在较大时可自适应减小,在极小或是零时增大。因此,式中分母项矩阵满足非奇异性。假如η=0,式(19)即为式(18);若η=1,以单位矩阵替代对角矩阵,式(19)则等效为PCA。

1.3.2 局部空间邻域散度和总体散度计算

鉴于高光谱图像空间一致性特点,即在空间局部邻域内近邻属于同类概率较大。以像元zi:(xi, yi)为中心作方形窗口δ(zi),(xi, yi)为zi在图像中的空间坐标位置,则窗口为γ×γ(γ是正奇数)的空间邻域像元集可记作

(20)

式中,c=(γ-1)/2,zim:(xm, ym)对应空间邻域里第m个像元点。δ(zi)共有γ×γ个像元。则空间邻域距离可定义为

(21)

式中,,度量像元zi与邻域内像元间相似权重,且。对于所有训练样本局部空间邻域散度矩阵

(22)

此外,为揭示影像数据多样性,保持数据的整体结构,定义总体散度矩阵

(23)

式中, z是训练样本均值。

1.3.3 空-谱协同低维嵌入

为在嵌入空间中提取低维空-谱鉴别特征,不仅要保持高光谱数据局部空间近邻结构,还需使超图中的类内数据聚集、类间数据远离。因此,式(19)、式(22)和式(23)可进一步转化为以下优化问题

(24)

式中,参数η, ξ∈[0, 1];Mw=ξ[(1-η)Nw+ηdiag(diag(Nw))]+(1-ξ)A,用于表征类内与局部数据紧致性,Nw=ZLwZTMb=ξ[(1-η)Nb+ηZZT]+(1-ξ)B,用于表示类间与全局数据发散度,Nb=ZLbZT

依据拉格朗日乘子法,式(24)可转换为以下广义特征值求解

(25)

将式(25)特征值降序排列,选取前τ个特征值对应的特征向量构成最优映射矩阵P=[p1, p2, …, pτ-1, pτ]。在低维空间里,测试样本ztest的空-谱协同特征是ytest=PTztest

2 试验结果与分析

为验证本文算法的有效性,在公开的Indian Pines和PaviaU高光谱数据集上进行分类试验,并与相关的维数约简算法进行了对比。

2.1 数据集

(1) Indian Pines数据集为美国宇航局在1992年利用AVIRIS传感器拍摄位于美国Indian州西北100 km2范围的高光谱遥感影像,其尺寸为145×145像素,共220个波段,空间分辨率为20 m,剔除受水气(噪声)影响的波段后,余下200个波段用于试验。该数据集主要包含16类地物,其假彩色图和真实地物图如图 3所示。

图 3 Indian Pines高光谱图像 Fig. 3 Indian Pines hyperspectral image

(2) PaviaU数据集为2002年采用ROSIS传感器拍摄的意大利北部的帕维亚大学周围的高光谱影像,其尺寸为610×340像素,空间分辨率为1.3 m,共有115个波段,去除受噪声影响严重的12个波段后,剩余103个波段用于对比试验。该数据集包括道路、砖块、屋顶和裸土等9类地物,图 4为其假彩色图和真实地物图。

图 4 PaviaU高光谱图像 Fig. 4 University of Pavia hyperspectral image

2.2 试验设置

在试验中,每次试验随机选取一定数目的样本用于训练,其余进行测试。鉴于在实际应用中,高光谱图像中存在部分地物类别样本数量非常少,例如在Indian Pines数据集中,Alfalfa(46)、Oats(20)、Stone-steel towers(93),括号中为对应的样本数。为避免出现某些类别选取训练样本所占比例过高或数量过少,在试验中设置如下:假设每类地物随机选取样本量为niNi表示某类地物的总样本数,若niNi/2,则niNi/2;若ni≤10,则定ni=10。通过采用各维数约简算法得到投影矩阵后,将所有样本投影到低维空间得到嵌入特征,并通过利用最近邻分类器(1-NN)进行分类。在每种试验条件下均进行10次重复试验,将总体分类精度(the overall accuracies,OAs)、平均分类精度(the average accuracies,AAs)及Kappa系数作为分类结果的评价指标。

试验中,将本文方法与PCA、LDA、MFA、LPP、RLDE、DHLP、SSCE、LPSNPE等维数约简算法进行比较,采用交叉验证方法获得各算法的最佳参数。SSCE在两个数据集中空间窗口均设置为5,SSCE和LPP的最近邻取5,DHLP中近邻数为9;RLDE和MFA的类内和类间近邻数分别为3、5,8、60。LDA的嵌入维数为u-1,u为类别数,其他算法的嵌入维数均设置为30。

为探索本文方法中参数ηξ,空间窗口γ对分类精度的影响,从数据集中每类地物中随机选取5个样本进行训练,其余样本作为测试样本。令α=10,ε=0.006,ηξ的取值范围均设置为{0, 0.01, 0.05, 0.1, 0.2, …, 0.9, 1},γ={3, 5, …, 39}。图 5为本文SSRSHE算法在不同ηξ值下的分类结果,图 6是本文SSRSHE算法在不同γ下的分类结果。

图 5 SSRSHE在不同ηξ参数值下的总体分类精度 Fig. 5 OAs of SSRSHE with different values of parameters η and ξ on Indian Pines and PaviaU data sets

图 6 SSRSHE在不同空间窗口γ下的总体分类精度 Fig. 6 OAs of SSRSHE with different size γ on different data sets

图 5可知,随着ξ的增加,其分类精度随之增加而后达到平稳,但是ξ值过大时,分类精度有所下降。这是因为在SSRSHE中,ξ用于平衡光谱信息和空间结构在特征提取中作用,ξ过小时未能有效利用超图所表征的像元间的复杂多元结构关系,过大时则忽略了空间结构,也不利于鉴别特征提取。与此同时,尽管试验中每类样本数量仅有5个,但是在同一η值下,分类结果比较稳定,有利于实际场景应用。为平衡光谱信息与空间信息对分类性能的影响,依据试验结果,本文在Indian Pines数据集设置ξ为0.3,η为0.7;对于PaviaU数据集,设置η=0.5及ξ=0.2。

图 6知,随着空间窗口γ变大,能利用的空间信息愈发丰富,分类精度随之增加;但γ过大时,空间窗口内包含来自于不同类数据的可能性增大,导致分类性能下降,且窗口过大,会导致计算复杂度增加。因此,综合考虑算法性能及计算效率,在Indian Pines数据集上设置γ=7,在PaviaU数据集上γ=15。

2.3 Indian Pines试验结果与分析

试验中,从Indian Pines数据集的每类地物里分别按照5、20、50、100、200样本数随机选取数据用于训练,剩余样本用于测试。采用各维数约简算法训练得到嵌入特征后,采用1-NN进行分类。表 1为在不同样本数量下不同算法的总体分类精度和Kappa系数值。

表 1 不同降维算法在Indian Pines数据集上的分类效果 Tab. 1 Classification with different numbers of training data via different DR methods on Indian Pines data set
5 20 50 100 200
总体分类精准度/(%) Kappa 总体分类精准度/(%) Kappa 总体分类精准度/(%) Kappa 总体分类精准度/(%) Kappa 总体分类精准度/(%) Kappa
RAW 43.6±2.8 0.372 54.9±1.7 0.495 60.1±1.4 0.552 63.6±0.9 0.588 66.9±0.6 0.622
PCA 43.4±2.7 0.370 54.9±1.6 0.495 60.2±1.2 0.553 63.9±0.8 0.591 67.0±0.6 0.622
LDA 32.5±4.8 0.253 51.6±1.9 0.459 64.4±1.2 0.599 71.0±0.5 0.672 74.4±0.7 0.706
LPP 43.6±3.7 0.371 54.5±1.8 0.491 59.7±1.2 0.546 62.7±1.0 0.578 65.8±0.5 0.609
MFA 44.1±4.0 0.377 57.1±1.6 0.520 66.8±1.9 0.625 70.8±1.1 0.669 72.0±1.0 0.680
RLDE 41.7±3.7 0.351 60.9±1.5 0.561 69.8±1.4 0.659 74.6±0.7 0.711 78.4±0.6 0.751
RSHE 48.6±3.4 0.422 63.2±1.9 0.587 71.0±1.7 0.672 77.1±0.9 0.739 80.0±0.7 0.770
DHLP 44.1±3.8 0.377 57.2±2.1 0.522 68.9±1.2 0.649 73.8±0.8 0.702 77.6±0.7 0.741
SSCE 30.2±4.5 0.230 69.7±1.0 0.658 76.3±0.9 0.730 79.1±0.5 0.760 82.9±0.6 0.801
LPSNPE 60.2±3.5 0.594 74.0±1.4 0.706 79.3±0.7 0.759 81.6±0.6 0.791 84.2±0.6 0.817
SSRSHE 65.6±2.3 0.615 74.8±1.2 0.706 80.0±1.0 0.765 82.9±1.0 0.803 86.7±1.0 0.829

表 1可得知,各种维数约简算法的分类性能都随着训练样本数目的增大而不断提高,这是由于随着训练数据量的增加,蕴含的信息就越丰富,有利于特征提取。DHLP、RSHE等超图方法的分类精度大多数情况下均优于传统图嵌入方法,表明利用数据间的多元几何结构特性可有效提高分类精度。与此同时,SSCE、LPSNPE等空-谱类方法,通过融合样本数据的空间信息,其分类性能要优于PCA、LDA、LPP、MFA、RLDE等仅利用了光谱信息的图嵌入方法。在各种训练条件下SSRSHE方法的分类性能均优于其他算法,因为它利用了超图框架来表示各样本邻域内顶点间的多元几何关系,因而可更好描述数据中复杂邻域结构。同时SSRSHE将样本类别信息融入超图框架,分别构建了稀疏本征超图和惩罚超图,能充分揭示数据间的复杂判别多元关系,提取出更有效的低维鉴别特征,进一步提升分类精度。

为进一步探索SSRSHE对每种地物的分类性能,从Indian Pines数据集每一类里随机选择3%的像元为训练样本,余下数据用于测试。表 2为不同维数约简算法对于每一种地物的总体分类精度、平均分类精度、Kappa系数及降维运行时间,其对应在整个数据集上的分类结果如图 7所示。

表 2 不同算法在Indian Pines数据集每类地物上的分类精度 Tab. 2 Classification accuracy of different types of features on Indian Pines data set by different algorithms
(%)
class train test RAW PCA LDA LPP MFA RLDE RSHE DHLP SSCE LPSNPE SSRSHE
1 10 36 41.67 41.67 77.78 36.11 50.00 61.11 47.22 63.89 55.56 77.78 94.44
2 143 1285 53.39 52.45 64.12 51.05 56.03 72.14 66.15 66.69 60.31 80.78 88.17
3 83 747 57.30 55.29 57.70 47.12 50.07 61.58 63.32 61.58 63.45 74.30 80.46
4 24 213 41.78 44.60 52.58 43.19 21.60 58.69 59.62 59.15 51.17 77.00 84.04
5 48 435 78.85 78.62 89.43 77.93 78.85 86.90 82.99 87.82 80.46 91.72 96.55
6 73 657 90.26 89.50 95.74 91.02 94.67 95.28 97.41 95.89 94.52 96.04 97.02
7 10 18 77.78 88.89 100 88.89 77.78 94.44 100 94.44 100 94.44 100
8 48 430 95.58 95.58 99.53 93.95 93.26 99.30 93.49 99.77 92.56 99.53 98.60
9 10 10 70.00 70.00 60.00 50.00 70.00 80.00 80.00 90.00 90.00 100 80.00
10 97 875 61.03 60.46 60.91 57.49 42.06 68.91 73.83 63.20 72.11 82.74 83.89
11 246 2209 69.76 69.85 71.89 69.62 58.85 79.36 82.66 79.22 74.02 85.92 89.50
12 59 534 39.33 37.45 65.36 32.02 47.38 67.42 60.30 62.73 50.56 87.83 83.71
13 21 184 88.04 88.04 97.83 88.04 94.57 97.28 95.65 98.37 94.57 98.91 100
14 127 1138 94.02 93.94 94.11 92.88 90.69 96.66 93.15 95.61 93.94 96.10 95.52
15 39 347 31.12 30.55 54.18 25.07 42.65 40.92 56.48 48.41 55.04 71.47 83.57
16 10 83 91.57 91.57 90.36 85.54 85.54 90.36 84.34 92.77 84.34 92.77 97.59
OA 68.33 67.88 74.44 65.91 64.03 78.27 78.17 77.00 74.06 87.65 89.78
AA 67.59 68.03 76.97 64.37 65.87 78.15 77.29 78.72 75.79 88.02 90.88
Kappa 0.638 0.633 0.707 0.609 0.589 0.751 0.750 0.736 0.704 0.858 0.884
DR time 0 0.01 0.01 0.15 0.20 3.28 11.87 5.86 432.8 15.15 27.87

图 7 在Indian Pines数据集上,各降维算法对应的全分类结果 Fig. 7 Classification map of different DR methods on Indian Pines data set

表 2可发现,SSRSHE的分类性能表现最佳,在每类上的总体分类精度、平均分类精度、Kappa值均优于其他方法,且对比SSCE算法,其运行效率快,优势明显。这是因为SSRSHE算超图学习,充分揭示了数据间高阶关系,以及像元空间特征的有效利用,有效表征了影像内蕴特性,提取的嵌入特征更具鉴别力,更有助于地物分类。同时,从图 7可以看到,本文算法相比其他算法,在其分类结果图更趋于平滑,尤其在“Alfalfa”、“Soybeans-min”、“Stone-steel towers”等区域更明显。由此可见,本文算法基于空-谱信息与超图模型协同学习,实现有效鉴别特征提取,改善影像分类精度,确实具有一定实践意义。

2.4 PaviaU试验结果与分析

在试验中,从每种地物中随机选取5、20、50、100、200个样本用于训练,其余数据用来测试,采用最近邻分类器进行分类。表 3为在不同的训练样本数目下各维数约简算法对应的总体分类精度及Kappa值。

表 3 不同算法在PaviaU数据集上的分类效果 Tab. 3 Classification with different numbers of training data via different DR methods on PaviaU data set
算法 5 20 50 100 200
总体分类精准度/(%) Kappa 总体分类精准度/(%) Kappa 总体分类精准度/(%) Kappa 总体分类精准度/(%) Kappa 总体分类精准度/(%) Kappa
RAW 60.5±4.2 0.512 66.4±2.4 0.583 73.5±1.6 0.663 76.4±0.8 0.698 78.8±0.8 0.724
PCA 60.5±4.2 0.512 66.5±2.2 0.583 73.4±1.6 0.662 76.4±0.8 0.697 78.7±0.8 0.724
LDA 46.7±6.4 0.351 59.6±1.8 0.495 73.5±1.4 0.662 78.9±0.9 0.727 83.4±0.6 0.782
LPP 47.0±5.6 0.354 59.3±2.6 0.500 72.8±2.3 0.654 78.3±1.3 0.722 82.2±1.2 0.768
MFA 64.5±4.3 0.555 69.2±4.5 0.613 76.4±2.0 0.699 78.1±2.4 0.715 79.1±2.2 0.730
RLDE 64.4±3.2 0.555 74.6±2.7 0.677 77.9±2.2 0.718 82.1±1.0 0.770 84.8±1.0 0.802
RSHE 63.2±4.1 0.540 75.4±2.3 0.685 78.3±1.4 0.720 83.4±0.9 0.784 84.9±1.3 0.802
DHLP 56.8±8.0 0.471 62.2±3.6 0.530 70.8±2.1 0.629 77.5±2.7 0.711 80.2±1.5 0.742
SSCE 42.3±5.3 0.309 63.3±2.9 0.543 75.8±1.7 0.692 82.7±1.2 0.814 87.0±0.8 0.828
LPSNPE 68.0±4.2 0.606 80.0±2.2 0.747 86.3±1.3 0.822 87.9±0.9 0.842 89.9±0.6 0.877
SSRSHE 71.6±2.7 0.646 82.6±2.3 0.776 87.5±1.1 0.837 90.0±1.5 0.882 92.2±0.2 0.908

依据表 3,在大多数训练条件下,DHLP、RSHE等超图方法和SSCE、LPSNPE等空-谱联合方法的分类结果要优于直接图嵌入方法,这表明超图学习和空-谱融合信息均有利于高光谱数据鉴别特征提取,有效改善地物分类性能。本文提出的SSRSHE方法在各种试验条件下,均具有最佳分类性能,这是因为其不仅通过超图学习发现高光谱数据中复杂结构,且有效融入了空间信息,在低维空间中使同类信息聚集、非同类信息远离,提高了数据可分性,进而有效提高地物分类效果。

为进一步分析SSRSHE方法在每种地物上的分类性能,从每类地物里随机选择5%的像元组成训练样本集,其他部分为测试样本集。表 4反映了不同维数约简方法在每类地物的分类效果,图 8则为各方法对整个PaviaU遥感图像分类的结果图。由表 4可以看到,SSRSHE在大多数地物类别中的分类性能要优于其他方法,表明在影像地物分类过程中,SSRSHE算法可使同类数据的关联性,异物数据间奇异性增强,鉴别特征尤为突出,分类性能更佳。同时,在图 8中,本文方法在“Asphalt”,“Meadows”,“Gravel”等地物区域的分类结果较为光滑,误分点较少,且运行时间并没大幅度增加,表明联合空-谱特性与超图学习的SSRSHE算法的地物分类性能有明显提升,更适合实际应用场景。

表 4 不同算法在PaviaU数据集每种地物上的分类精度 Tab. 4 Classification accuracy of different types of features on PaviaU data set by different algorithms
(%)
class train test RAW PCA LDA LPP MFA RLDE RSHE DHLP SSCE LPSNPE SSRSHE
1 332 6299 85.62 85.62 87.68 87.82 82.76 90.19 87.46 63.89 89.73 90.20 91.19
2 933 17716 94.65 94.57 94.88 94.76 93.90 97.73 95.80 66.69 96.70 97.53 98.12
3 105 1994 65.15 64.64 63.34 67.00 61.84 74.77 69.71 61.58 72.37 77.28 78.69
4 154 2910 77.22 77.36 81.79 79.01 77.02 84.13 79.59 59.15 84.78 87.83 89.26
5 68 1277 98.83 98.83 98.84 99.30 99.77 99.53 98.98 87.82 99.37 99.77 99.77
6 252 4777 60.26 60.32 65.17 65.47 69.72 70.36 65.24 95.89 73.86 89.68 85.22
7 67 1263 75.30 75.30 66.67 75.69 71.26 80.36 82.58 94.44 88.60 86.06 90.18
8 185 3497 80.27 80.27 74.39 81.42 77.36 84.79 75.24 79.77 82.85 84.68 79.33
9 48 899 100 100 99.44 100 99.67 100 98.44 90.00 99.78 99.89 100
OA 84.92 84.88 85.40 86.27 84.73 89.70 87.31 78.00 89.60 91.30 92.59
AA 81.92 81.88 81.47 83.38 81.48 86.87 85.51 77.72 87.56 89.53 90.55
Kappa 0.797 0.796 0.804 0.815 0.796 0.861 0.825 0.736 0.861 0.883 0.902
DR time 0 0.08 0.03 0.81 1.67 6.80 8.21 5.86 896.2 7.29 15.19

图 8 在PaviaU数据集上,各降维算法对应的全分类结果 Fig. 8 Classification map of different DR methods on PaviaU data set

3 总结

针对传统图嵌入降维方法存在不能表征高光谱数据中的多元关系且未有效利用空间信息等问题,本文提出了一种空-谱协同正则化稀疏超图嵌入算法。本文算法利用稀疏系数实现自适应近邻选取,构建正则化稀疏超图模型来揭示高光谱数据间的多元几何结构。此外,考虑到保持样本的全局特性和局部邻域结构分别定义样本总体散度与局部空间邻域散度,实现空-谱鉴别特征提取。在Indian Pines和PaviaU高光谱数据集上试验结果表明,相比其他算法,在训练样本数较少时,SSRSHE地物分类性能仍有明显提升。但本文方法仅运用光谱信息构建超图,在下一步工作将考虑空-谱联合超图模型构建,以进一步提升地物分类效果。


参考文献
[1] 侯榜焕, 王锟, 姚敏立, 等. 面向高光谱图像分类的半监督空谱判别分析[J]. 测绘学报, 2017, 46(9): 1098–1106.
HOU Banghuan, WANG Kun, YAO Minli, et al. Semi-supervised spatial-spectral discriminant analysis for hyperspectral image classification[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(9): 1098–1106. DOI:10.11947/j.AGCS.2017.20170121
[2] 张良培, 武辰. 多时相遥感影像变化检测的现状与展望[J]. 测绘学报, 2017, 46(10): 1447–1459.
ZHANG Liangpei, WU Chen. Advance and future development of change detection for multi-temporal remote sensing imagery[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(10): 1447–1459. DOI:10.11947/j.AGCS.2017.20170340
[3] MENEZES J, POOJARY N. Dimensionality reduction and classification of hyperspetral images using DWT and DCCF[C]//Proceedings of the 3rd MEC International Conference on Big Data and Smart City. Muscat: IEEE, 2016.
[4] CHEN Mulin, WANG Qi, LI Xuelong. Discriminant analysis with graph learning for hyperspectral image classification[J]. Remote Sensing, 2018, 10(6): 836. DOI:10.3390/rs10060836
[5] FENG Fubiao, LI Wei, DU Qian, et al. Dimensionality reduction of hyperspectral image with graph-based discriminant analysis considering spectral similarity[J]. Remote Sensing, 2017, 9(4): 323.
[6] 罗甫林. 高光谱图像稀疏流形学习方法研究[J]. 测绘学报, 2017, 46(3): 400.
LUO Fulin. Sparse manifold learning for hyperspectral imagery[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(3): 400. DOI:10.11947/j.AGCS.2017.20160621
[7] WANG Hao, FAN Yuanyuan, FANG Baofu, et al. Generalized linear discriminant analysis based on Euclidean norm for gait recognition[J]. International Journal of Machine Learning and Cybernetics, 2018, 9(4): 569–576. DOI:10.1007/s13042-016-0540-0
[8] FAN Mingyu, ZHANG Xianqin, QIAO Hong, et al. Efficient isometric multi-manifold learning based on the self-organizing method[J]. Information Sciences, 2016(345): 325–339.
[9] 王东, 张强, 严亮. 一种融合聚类的监督局部线性嵌入算法研究[J]. 半导体光学, 2017, 38(3): 419–424.
WANG Dong, ZHANG Qiang, YAN Liang. Study on supervised local linear embedding algorithm based on fusion clustering[J]. Semiconductor Optoelectronics, 2017, 38(3): 419–424.
[10] WANG Qi, MENG Zhaotie, LI Xuelong. Locality adaptive discriminant analysis for spectral-spatial classification of hyperspectral images[J]. IEEE Geoscience and Remote Sensing Letters, 2017, 14(11): 2077–2081. DOI:10.1109/LGRS.2017.2751559
[11] JIANG Quansheng, ZHU Qixin, WANG Bangfu, et al. Nonlinear machine fault detection by semi-supervised Laplacian eigenmaps[J]. Journal of Mechanical Science and Technology, 2017, 31(8): 3697–3703. DOI:10.1007/s12206-017-0712-1
[12] DENG Yangjun, LI Hengchao, PAN Lei, et al. Modified tensor locality preserving projection for dimensionality reduction of hyperspectral images[J]. IEEE Geoscience and Remote Sensing Letters, 2018, 15(2): 277–281. DOI:10.1109/LGRS.2017.2786223
[13] KUMAR S, BHUYAN M K, LOVELL B C, et al. Hierarchical uncorrelated multiview discriminant locality preserving projection for multiview facial expression recognition[J]. Journal of Visual Communication and Image Representation, 2018(54): 171–181.
[14] TAŞKIN G, KAYA H, BRUZZONE L. Feature selection based on high dimensional model representation for hyperspectral images[J]. IEEE Transactions on Image Processing, 2017, 26(6): 2918–2928. DOI:10.1109/TIP.2017.2687128
[15] LUO Fulin, HUANG Hong, DUAN Yule, et al. Local geometric structure feature for dimensionality reduction of hyperspectral imagery[J]. Remote Sensing, 2017, 9(8): 790. DOI:10.3390/rs9080790
[16] ZHOU Yicong, PENG Jiangtao, CHEN C L P. Dimension reduction using spatial and spectral regularized local discriminant embedding for hyperspectral image classification[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(2): 1082–1095.
[17] YUAN Haoliang, TANG Yuanyan. Learning with hypergraph for hyperspectral image feature extraction[J]. IEEE Geoscience and Remote Sensing Letters, 2015, 12(8): 1695–1699. DOI:10.1109/LGRS.2015.2419713
[18] DU Weibao, QIANG Wenwen, LÜ Meng, et al. Semi-supervised dimension reduction based on hypergraph embedding for hyperspectral images[J]. International Journal of Remote Sensing, 2018, 39(6): 1696–1712. DOI:10.1080/01431161.2017.1415480
[19] HUANG Sheng, YANG Dan, GE Yongxin, et al. Discriminant hyper-Laplacian projections and its scalable extension for dimensionality reduction[J]. Neurocomputing, 2016, 173(2): 145–153.
[20] ZHANG Zhihong, BAI Lu, LIANG Yuanheng, et al. Joint hypergraph learning and sparse regression for feature selection[J]. Pattern Recognition, 2017(63): 291–309.
[21] WU Zebin, SHI Linlin, LI Jun, et al. GPU parallel implementation of spatially adaptive hyperspectral image classification[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2018, 11(4): 1131–1143. DOI:10.1109/JSTARS.2017.2755639
[22] 黄鸿, 郑新磊. 高光谱影像空-谱协同嵌入的地物分类算法[J]. 测绘学报, 2016, 45(8): 964–972.
HUANG Hong, ZHENG Xinlei. Hyperspectral image land cover classification algorithm based on spatial-spectral coordination embedding[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(8): 964–972. DOI:10.11947/j.AGCS.2016.20150654
[23] SUN Yubao, WANG Sujuan, LIU Qingshan, et al. Hypergraph embedding for spatial-spectral joint feature extraction in hyperspectral images[J]. Remote Sensing, 2017, 9(5): 506. DOI:10.3390/rs9050506
[24] SUN Tao, YIN Penghang, CHENG Lizhi, et al. Alternating direction method of multipliers with difference of convex functions[J]. Advances in Computational Mathematics, 2018, 44(3): 723–744. DOI:10.1007/s10444-017-9559-3
[25] SUN Yanfeng, ZHAO Jiangang, HU Yongli. Supervised sparsity preserving projections for face recognition[C]//Proceedings of SPIE 8009, Third International Conference on Digital Image Processing. Chengdu, China: SPIE, 2011: 357-366.
http://dx.doi.org/10.11947/j.AGCS.2019.20180469
中国科学技术协会主管、中国测绘地理信息学会主办。
0

文章信息

黄鸿,陈美利,王丽华,李政英
HUANG Hong, CHEN Meili, WANG Lihua, LI Zhengying
空-谱协同正则化稀疏超图嵌入的高光谱图像分类
Using spatial-spectral regularized hypergraph embedding for hyperspectral image classification
测绘学报,2019,48(6):676-687
Acta Geodaetica et Cartographica Sinica, 2019, 48(6): 676-687
http://dx.doi.org/10.11947/j.AGCS.2019.20180469

文章历史

收稿日期:2018-10-15
修回日期:2019-03-21

相关文章

工作空间