Hyperspectral image classification method based on neighborhood spectra and probability cooperative representation

高光谱图像具有光谱分辨率高、波段数量多、光谱信息丰富等特点，被广泛应用于军事^[1]、农业^[2]、矿业^[3]、渔业^[4]、林业^[5]等领域。在高光谱图像的众多应用中，分类是前提，其目标是给图像中的每个像元赋以唯一的类别标识。高光谱图像在为地物分类和识别提供了细致丰富的光谱信息的同时，其数据的高维性与谱间存在的大量冗余信息也给分类带来了困难。如果直接用高光谱数据对地物进行分类，容易产生Hughes现象，且计算量大、耗时长。为了克服这些困难，经常采用数据降维技术。

对于数据降维，学者们提出了许多方法，这些方法通常依据某种准则，将高维的谱向量投影到低维的特征空间，然后在低维特征空间进行分类。经典的线性降维方法有主成分分析(principal component analysis, PCA)和线性判别式分析(linear discriminant analysis, LDA)。PCA和LDA都是在高光谱数据满足全局线性的假设下通过寻找不同的线性模型进行降维，但高光谱数据往往是非线性的，因此，采用线性的PCA和LDA无法揭示其非线性结构。为了在提取的特征中揭示高光谱数据的非线性结构，研究者通过核函数将低维的线性空间映射到高维的非线性空间，然后进行数据降维。这些方法包括核主成分分析(kernel principal component analysis，KPCA)^[6-7]、核Fisher判别式分析^[8]、广义判别式分析^[9]和支持向量机(support vector machines，SVM)。流形学习算法是另一种揭示样本非线性结构的方法。与基于核函数的方法相比较，流形学习算法通过邻接图的嵌入能够有效保留数据的局部结构信息^[10-13]。

近来，稀疏表示分类器(sparse representation cla-ssification, SRC)被广泛应用^[14-16]。与其它的分类器相比较，SRC不需要训练过程，而是用所有的训练样本组成一个字典，并通过稀疏最小化重构误差实现对测试样本的分类。由于对每一个测试样本都需要进行复杂的编码过程，因此，在数据量较大时，SRC的算法复杂度是巨大的。ZHANG等人^[17]的研究表明，在对测试样本进行分类时，起主要作用的是所有样本对测试样本的协同表示而不是稀疏性约束，并提出一种基于l₂范数编码的协同表示分类方法(collaborative representation based classification, CRC)。与基于l₁范数编码的SRC相比较，这种基于l₂范数的编码方法能够得到解析解，因此具有更小的计算复杂度，且在人脸识别与高光谱分类方面具有较高的识别性能^[18-21]。最近，CAI等人从概率的角度出发，分析了CRC的机理并提出概率协同表示分类方法(probabilistic collaborative representation based classifier, ProCRC)^[22]，在MNIST数据库与AR数据库上取得了较好的识别结果。

上述方法在分类时仅仅考虑了像元的光谱信息，而从空间的角度来看，空间相邻的像元最有可能属于同一类，即具有空间一致性。为了充分利用像元的空间信息，近年来，研究者提出了许多空谱联合的高光谱图像的分类方法。这些方法在相似性度量和特征提取时考虑了像元的空间相关性，取得了优于纯光谱特征的分类精度^[23-24]。因此，如何充分利用空间信息与光谱信息进行像元的分类，是高光谱图像分类的研究热点。

现有的大多数基于空谱信息的高光谱分类方法仅利用的是单个像元的光谱数据，而事实上高光谱遥感数据是通过光谱技术与成像技术相结合得到的，最终的光谱数据中可能存在高光谱传感器噪声、无光照时仪器自身带的暗电流、大气校正时的偏差等；在反射率计算时，混合像元、地形、光照等问题也会影响到最终的数据。由此可见，每个像元的光谱数据易受到噪声污染且不是孤立的，而是与其邻域其它像元紧密相关的。为了降低单个像元光谱数据对噪声的敏感性并充分利用近邻像元的光谱相似特点，本文中对每个像元的光谱数据用半径为r的圆形邻域内像元光谱数据生成的邻域谱代替，这样可有效克服噪声影响并充分利用空间信息。生成邻域谱后，再采用概率协同表示分类器进行分类。

以坐标(i, j)为中心点的像元，图 1中给出了其半径为r的P个像元构成的邻域，表示为G(r, P)。P个邻域像元中，第m(m=1, …, P)个邻域点的坐标位置可通过下式计算：

Figure 1.  Example of neighborhood pixel sampling

当邻域点未落在像元中心位置时，其谱向量可通过双线性插值计算得到，如图 2所示。假定像元(i, j)的某个邻域像元为v，距离v最近的4个像元分别是v₁, v₂, v₃和v₄，且它们的谱向量分别为s₁, s₂, s₃和s₄，各点之间的距离为1个像素。假定像元v到点v₁和点v₃的水平距离为a，到点v₂的垂直距离为b，则v的谱向量g可由下式计算得到：

Figure 2.  v-point spectral vector interpolation

在计算得到像元(i, j)所有P个领域点的谱向量g_i(i=1, …, P)后，若像元(i, j)谱向量为s，则其邻域谱$\mathit{\boldsymbol{\hat s}}$可由下式得到：

假定第i类训练样本表示为矩阵X_i=[x_{i, 1}, x_{i, 2}, …, x_{i, ni}]，向量x_{i, j}∈R^d是第i类的第j个样本，则训练基共有N=c×n个样本，所有的训练样本集可表示为矩阵X=[X₁, X₂, …, X_c]，R^d表示d子维空间，c是样本的类别数量，l_X={1, 2, …, c}表示样本标签集合，并将全体训练样本张成的子空间S称为协同子空间。S空间中的样本向量x可表示为$\mathit{\boldsymbol{x}} = \mathit{\boldsymbol{X\alpha }} = \sum\limits_{i = 1}^c {{\mathit{\boldsymbol{X}}_i}} {\mathit{\boldsymbol{\alpha }}_i}$，其中向量α=[α₁, α₂, …, α_c]，α_i是对应于X_i的编码。若l(x)表示样本x的标签，则样本x可由训练集中的样本协同地表示为$\mathit{\boldsymbol{x}} = \mathit{\boldsymbol{X\alpha }} = \sum\limits_{i = 1}^c {{\mathit{\boldsymbol{X}}_i}} {\mathit{\alpha }_i}$，其属于第i类的概率可表示为^[22]：

式中，ε是一个常数。

对于待测样本y，其属于第i类的概率可表示为^[22]：

式中, λ和γ是常数。α可由下式计算得到^[22]：

式中，W_x为对角矩阵，其对角元素为：

式中，X_{i, j}表示矩阵X的第i行第j列元素。

在得到解向量α后，可由(5)式得到样本y属于第i(i=1, …, c)类的概率：

由(8)式可知，对所有的类，${\left\| {\mathit{\boldsymbol{y}} - \mathit{\boldsymbol{X}}\mathit{\boldsymbol{\widehat \alpha }}} \right\|_2}^2 + \lambda {\left\| {\mathit{\boldsymbol{\widehat \alpha }}} \right\|_2}^2$是相同的，因此，样本y属于第i类的概率可表示为：

这样，样本y的类别归属由下式表示的规则决定：

实验中选取美国印第安纳西北地区的Indian Pines高光谱图像和美国加州的萨利纳斯山谷的高光谱图像Salinas scene，利用总体识别率和kappa值作为分类性能的评价指标，并和主成分分析(PCA)、支持向量机(SVM)、基于稀疏表示的分类器(SRC)、基于协同表示的分类器(CRC)方法进行性能比较。

Indian Pines图像大小为145pixel×145pixel，波长范围为400nm~2500nm，光谱分辨率为10nm。原始波段数为224个，去除坏波段和和水体吸收波段，用剩余的200个波段进行分类实验，图 3显示了真实的地物分布情况。

Figure 3.  Indian Pines data and its ground truth value

图 4中给出的图像位置为80行50列~90行70列像元的原始谱数据与生成的领域谱数据的比较。其中图 4a为原始数据，图 4b为邻域谱数据(图中波带数和幅值均无单位)。对比两个图可以看出，像元的邻域谱较原始谱具有更好的聚类性，尤其在第40~第100波段之间的数据，明显优于原始的谱特征。

Figure 4.  Comparison between clustering of pixels original spectrum and neighborhood spectrum

在分类实验中，选择样本数大于400的地物进行训练，因此，参与分类实验的样本共有9种。每一类地物随机选择50个样本作为训练样本，其余作为测试样本进行实验。表 1中给出了4种算法的性能指标。从表 1可知，本文中算法的分类精度与kappa系数明显优于其余几种算法，具有较好的分类性能。

Salinas scene图像大小为512pixel×217pixel，空间分辨率为3.7m，用去掉20个水吸收波段后剩余的204个波段进行实验。本次实验中采用了Salinas scene图像的一个大小为86pixel×83pixel的子图像，其地面真实数据如图 5所示。

Figure 5.  Salinas scene data and its ground truth value

数据库2上的分类结果如表 2所示。从表 2可知，在Salinas scene图像上，除了PCA算法，其余4种算法都表现出良好的识别性能，但本文中算法在识别性能与kappa系数上都有一定的优势。

在Indian Pines与Salinas scene数据库上的实验结果表明，随着样本数量与类别数的增加，5种算法的识别性能都有所下降。本文中算法由于充分利用了近邻像元的光谱相似特性，得到的像元光谱特征更稳定，因此获得了更好的识别性能。

高光谱图像的空间一致性表明在同一邻域的像元很大概率上属于同类地物的，其光谱数据具有相似性，然而在高光谱数据的生成过程很容易受到由各种因素带影响，容易造成同一邻域同类像元谱特征的差异，为了克服这个困难，本文中通过采用插值方法生成像元的邻域谱，充分利用了高光谱图像空间一致性特点。与原始的谱特征相比较，邻域谱的类内聚类性更好，更利于分类。其次，概率协同表示分类方法是最近提出的一种分类方法，研究表明其具有较高的分类精度。本文中结合邻域谱与概率协同表示分类器，将其应用于高光谱遥感图像的分类中，实验结果表明，所提出的基于邻域谱与概率协同表示方法的高光谱分类方法是一种高性能的分类方法。