一种高光谱遥感数据的降维方法及系统与流程
本发明涉及数据降维技术领域，具体涉及一种高光谱遥感数据的降维方法及系统。
背景技术
高光谱遥感数据通常包含几十到几百个波段，能以较高的可信度识别地物目标。高光谱遥感采集得到的图像其中包含了丰富的空间、辐射和光谱信息，为地物的精确分类带来了机遇。但是同时高光谱遥感图像光谱分辨率的提高是基于增加数据量和数据维数作来实现的，这一点也为高光谱遥感图像进行地物分类和目标识别等应用带来了很大的难题。高光谱数据的“维数灾难”能够使分类识别不能正常进行，阻碍了高光谱遥感技术的应用推广。同时，繁多的数据量和高维度也给高光谱图像处理带来了巨大的计算压力。一般情况下，在遥感分类中，样本数目越多，分类精度就会越高；于是增加样本数目成为提高分类精度的一个可能。但是在高光谱遥感中，由于样本数据点获取的难度比较大而且代价昂贵，高光谱遥感影像分类精度一直受限于样本数目。另一方面，研究也发现，在不同维数的情况下，并不是样本数目越多，分类精度就越高。
地物分类是当前高光谱数据遥感图像处理领域研究的热点。高光谱数据的光谱分辨率可达到纳米数量级，光谱图像光谱范围窄、谱间相关性强，信息冗余较大，由此给传输、存储和处理带来了很多新问题。如高维数据与计算复杂度之间的矛盾、高维特征，与小样本之间的矛盾等。当利用高光谱遥感图像对地物分类时，若利用全部波段进行分类，计算量太大，分类算法的复杂度会随着光谱维数的增加而增加。
技术实现要素：
为解决上述现有技术中在高光谱遥感中，由于样本数据点获取的难度比较大而且代价昂贵，高光谱遥感影像分类精度一直受限于样本数目的问题，本发明的目的是提供一种高光谱遥感数据的降维方法及系统，利用一种基于主成分分析的高光谱遥感数据降维方法可以很大程度上降低高光谱数据的冗余度，获得精度较高的分类结果。
本发明的目的是采用下述技术方案实现的：
本发明提供一种高光谱遥感数据的降维方法，其改进之处在于：
根据相邻波段的高光谱遥感数据间的相关系数获取高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线；
对所述数据子集的相关系数曲线主成分变换后对所述高光谱遥感数据进行降维处理。
进一步地：所述高光谱遥感数据通过星载高光谱传感器拍摄以及航空飞机/无人机所搭载的高光谱仪获取。
进一步地：所述根据相邻波段的高光谱遥感数据间的相关系数获取高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线，包括：
计算相邻波段的高光谱遥感数据间的相关系数；
根据所述相邻波段的高光谱遥感数据间的相关系数构造相邻波段的高光谱遥感数据间的相关系数的相关矢量；
根据相关矢量获得高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线。
进一步地：所述计算相邻波段的高光谱遥感数据间的相关系数的公式为：
式中，rij为波段i和波段j间的相关系数，xip和xjp分别为波段i和波段j内的第p个像元的辐射亮度值；
进一步地：所述根据所述相邻波段的高光谱遥感数据的相关系数构造相邻波段的高光谱遥感数据的相关系数相关矢量的表达式为：
rntr＝(r12,r23,…,r(i,i+1),…,r(n-1,n))t
式中：r12,r23,…,r(i,i+1),…,r(n-1,n)分别表示波段1和波段2的相关系数，波段2和波段3的相关系数，....，波段i和波段i+1的相关系数，...，波段n-1和波段n的相关系数。
进一步地：所述根据相关矢量获得高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线，包括：
从所述相关系数矢量中提取k个局部相关极小值；
根据所述局部相关极小值划分相关矢量为k+1个数据子集；
将数据子集中的相关系数作为纵坐标，将所述数据子集中对应的高光谱遥感数据所处的波段作为横坐标，构建数据子集的相关系数曲线。
进一步地：所述对所述数据子集的相关系数曲线主成分变换后对所述高光谱遥感数据进行降维处理，包括：
设定的相关系数阈值，比较所述数据子集的相关系数曲线上的相关系数和所述相关系数阈值；
将所述数据子集的相关系数曲线上的相关系数大于所述相关系数阈值的相关系数判定为同一主成分；
将属于同一主成分的相关系数所对应的高光谱遥感数据进行合并映射，并组成新的高光谱遥感数据波段。
本发明还提供一种高光谱遥感数据的降维系统，其改进之处在于：
计算模块，用于根据相邻波段的高光谱遥感数据间的相关系数获取高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线；
降维模块，用于对所述数据子集的相关系数曲线主成分变换后对所述高光谱遥感数据进行降维处理。
进一步地：所述高光谱遥感数据通过星载高光谱传感器拍摄以及航空飞机/无人机所搭载的高光谱仪获取。
进一步地：所述计算模块，包括：
算术子模块，用于计算高光谱遥感数据中相邻光谱波段的相关系数；
第一构造子模块，用于根据所述高光谱遥感数据中相邻光谱波段的相关系数构造高光谱遥感数据中相邻光谱波段的相关系数矢量；
获得子模块，用于根据相关矢量获得高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线。
进一步地：所述获得子模块，包括：
提取单元，用于从所述相关系数矢量中提取k个局部相关极小值；
划分单元，用于根据所述局部相关极小值划分相关矢量为k+1个数据子集；
获得单元，用于将数据子集中的相关系数作为纵坐标，将所述数据子集中对应的高光谱遥感数据所处的波段作为横坐标，构建数据子集的相关系数曲线。
进一步地：所述降维模块，包括：
设定子模块，用于设定的相关系数阈值，比较所述数据子集的相关系数曲线上的相关系数和所述相关系数阈值；
判定子模块，用于将所述数据子集的相关系数曲线上的相关系数大于所述相关系数阈值的相关系数判定为同一主成分；
合并子模块，用于将属于同一主成分的相关系数所对应的高光谱遥感数据进行合并映射，并组成新的高光谱遥感数据波段。
与最接近的现有技术相比，本发明提供的技术方案具有的有益效果是：
本发明提供一种高光谱遥感数据的降维方法，获取高光谱遥感数据；对高光谱遥感数据进行计算，得到高光谱遥感数据相邻波段的相关系数曲线；对所述相关系数曲线进行主成分变换后对所述高光谱遥感数据进行降维处理。本发明利用一种基于主成分分析的高光谱遥感数据降维方法不仅可以很大程度上降低高光谱数据的冗余度，并且可以提高遥感解译效率，获得精度较高的分类结果，对于高光谱遥感解译技术的发展有很大的意义。
附图说明
图1是本发明提供的一种高光谱遥感数据的降维方法的流程图；
图2是本发明提供的主成分变换特征提取原理图；
图3是本发明提供的aviris高光谱遥感图像示意图；(a)为aviris高光谱影像示意图；(b)为地面真实数据图；
图4是本发明提供的具体实施例的高光谱数据不同波段间的近邻可传递曲线图；
图5是本发明提供的具体实施例的aviris图像相关系数矩阵示意图；
图6是本发明提供的主成分分析变换实验-原始图像最大似然分类结果图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
以下描述和附图充分地示出本发明的具体实施方案，以使本领域的技术人员能够实践它们。其他实施方案可以包括结构的、逻辑的、电气的、过程的以及其他的改变。实施例仅代表可能的变化。除非明确要求，否则单独的组件和功能是可选的，并且操作的顺序可以变化。一些实施方案的部分和特征可以被包括在或替换其他实施方案的部分和特征。本发明的实施方案的范围包括权利要求书的整个范围，以及权利要求书的所有可获得的等同物。在本文中，本发明的这些实施方案可以被单独地或总地用术语“发明”来表示，这仅仅是为了方便，并且如果事实上公开了超过一个的发明，不是要自动地限制该应用的范围为任何单个发明或发明构思。
本发明提供一种高光谱遥感数据的降维方法，其流程图如图1所示，包括下述步骤：
根据相邻波段的高光谱遥感数据间的相关系数获取高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线；
对所述数据子集的相关系数曲线主成分变换后对所述高光谱遥感数据进行降维处理。
进一步地：所述高光谱遥感数据通过星载高光谱传感器拍摄以及航空飞机/无人机所搭载的高光谱仪获取。
进一步地：所述根据相邻波段的高光谱遥感数据间的相关系数获取高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线，包括：
计算相邻波段的高光谱遥感数据间的相关系数；
根据所述相邻波段的高光谱遥感数据间的相关系数构造相邻波段的高光谱遥感数据间的相关系数的相关矢量；
根据相关矢量获得高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线。
进一步地：所述计算相邻波段的高光谱遥感数据间的相关系数的公式为：
式中，rij为波段i和波段j间的相关系数，xip和xjp分别为波段i和波段j内的第p个像元的辐射亮度值；
进一步地：所述根据所述相邻波段的高光谱遥感数据的相关系数构造相邻波段的高光谱遥感数据的相关系数相关矢量的表达式为：
rntr＝(r12,r23,…,r(i,i+1),…,r(n-1,n))t
式中：r12,r23,…,r(i,i+1),…,r(n-1,n)分别表示波段1和波段2的相关系数，波段2和波段3的相关系数，....，波段i和波段i+1的相关系数，...，波段n-1和波段n的相关系数。
进一步地：所述根据相关矢量获得高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线，包括：
从所述相关系数矢量中提取k个局部相关极小值；
根据所述局部相关极小值划分相关矢量为k+1个数据子集；
将数据子集中的相关系数作为纵坐标，将所述数据子集中对应的高光谱遥感数据所处的波段作为横坐标，构建数据子集的相关系数曲线。
进一步地：所述对所述数据子集的相关系数曲线主成分变换后对所述高光谱遥感数据进行降维处理，包括：
设定的相关系数阈值，比较所述数据子集的相关系数曲线上的相关系数和所述相关系数阈值；
将所述数据子集的相关系数曲线上的相关系数大于所述相关系数阈值的相关系数判定为同一主成分；
将属于同一主成分的相关系数所对应的高光谱遥感数据进行合并映射，并组成新的高光谱遥感数据波段。
主成分分析(pca)同时也叫k-l变换，它将多个变量通过线性变换选出较少的变量。理论依据是将原变量重新组合成一组新的无关的几个综合变量，根据实际需要从中取出较少的几个能尽可能多地反映原变量信息。这种方法在数学变换中保持变量的总方差不变，同时使第一变量具有最大的方差，称为第一主成分，第二变量的方差仅小于第一主成分，并且和第一变量并不相关，称为第二主成分；依次类推，第三个变量就有第三个主分。
在高光谱遥感数据上，利用主成分分析方法提取特征，是按照一定的准则将原始特征由高维空间映射到较低的维数空间提取的。特征应尽可能地保留针对具体应用的有价值的信息。经过特征提取后的光谱特征空间中，新的光谱向量反映特定地物某一性状的一个光谱参量。特征选择的过程如图2所示，其中x是测量空间，y是特征空间，变换y＝f(x1，x2，…..，xn)是一个线性或者非线性的特征提取器或者转换方程，它将高维的空间x投影到一个低维的空间y。
实施例二、
本发明实施例研究选取美国aviris高光谱数据92av3c。92av3c是公开的基准高光谱图像数据。如图3(a)和(b)所示。在目前已经发表的学术文章中，已经有大量的国内外学者在高光谱图像处理的研究上采用了该数据，因此有利于重复实验和不同的算法之间的比较。该数据采集自美国印第安纳州的一块印度松树测试地，图像大小为145×145共21025个像素，一共有220光谱波段，图像的空间分辨率约为20m，光谱分辨率约等于10nm，数据采集于1992年6月。aviris高光谱传感器搭载er-2(u2)飞机在约6万米高空测量。该数据及真实标记图由美国普度大学的davidlandgrebe教授提供。第50，27，17波段作为rgb通道的假彩色图像。
(1)计算高光谱数据的各个相邻光谱波段的互相关系数，构造互相关系数矢量。对得到的互相关系数矢量r进行处理，从中提取若干个局部相关的极小值，然后根据这些提取的局部极小值(k个)。将原始高光谱数据空间s划分为k+1个适合的数据子集，把计算统计的相关系数导入到excel表中获得高光谱数据相邻波段的相关系数曲线。对于某一类特定的光谱成像仪获取的高光谱图像，则会具有相似的相关系数曲线。具体实施例的高光谱数据不同波段间的近邻可传递曲线图如图4所示。
(2)主成分变换可以降低图像相关性，图像波段相关性很强，主成分变换就可以把信息压缩到很少波段。图像相关系数矩阵如图5所示，在aviris图像的相关系数矩阵图中，灰度值越高的地方表示相关系数越高，对角线上都是各个波段与它本身的相关系数，值全为1，显示为白色，由图5可知，aviris图像波段之间相关性特别强，分块明显，因此图像的冗余信息比较多。经过主成分变换后，高光谱图像能量容易集中到很少的几个波段。
(3)进行主成分变换后选取的个数为vm>t(贡献率的值t＝0.99)，从主成分的统计情况得出，只需要选择取前4个主成分包含了原始图像的大部分信息。
(4)验证图像降维后信息损失，分别对原始图像和选取的主成分进行贝叶斯最大似然有监督分类。图3(b)是真实地物图，如图6所示是前4个主成分最大似然分类结果，可以看出，在分类图中，在边缘地区，由于主成分并未包含原始数据所有的信息，所以有少量像元被误分。主成分变换在高光谱图像降维中具有很好的效果，因此，它是普遍被使用的降维工具。
实施例三、
基于同样的发明构思，本发明还提供一种高光谱遥感数据的降维系统，包括：
计算模块，用于根据相邻波段的高光谱遥感数据间的相关系数获取高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线；
降维模块，用于对所述数据子集的相关系数曲线主成分变换后对所述高光谱遥感数据进行降维处理。
进一步地：所述高光谱遥感数据通过星载高光谱传感器拍摄以及航空飞机/无人机所搭载的高光谱仪获取。
进一步地：所述计算模块，包括：
算术子模块，用于计算高光谱遥感数据中相邻光谱波段的相关系数；
第一构造子模块，用于根据所述高光谱遥感数据中相邻光谱波段的相关系数构造高光谱遥感数据中相邻光谱波段的相关系数矢量；
获得子模块，用于根据相关矢量获得高光谱遥感数据的数据子集以及所述数据子集对应的相关系数曲线。
进一步地：所述获得子模块，包括：
提取单元，用于从所述相关系数矢量中提取k个局部相关极小值；
划分单元，用于根据所述局部相关极小值划分相关矢量为k+1个数据子集；
获得单元，用于将数据子集中的相关系数作为纵坐标，将所述数据子集中对应的高光谱遥感数据所处的波段作为横坐标，构建数据子集的相关系数曲线。
进一步地：所述降维模块，包括：
设定子模块，用于设定的相关系数阈值，比较所述数据子集的相关系数曲线上的相关系数和所述相关系数阈值；
判定子模块，用于将所述数据子集的相关系数曲线上的相关系数大于所述相关系数阈值的相关系数判定为同一主成分；
合并子模块，用于将属于同一主成分的相关系数所对应的高光谱遥感数据进行合并映射，并组成新的高光谱遥感数据波段。
本发明针对海量、高维信息，进行了数据的降维操作，降低数据信息的冗余，从而提高数据处理速度和提高分类器的性能。提出了有效的高光谱降维方法，能够提高遥感图像解译效率。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员依然可以对本发明的具体实施方式进行修改或者等同替换，这些未脱离本发明精神和范围的任何修改或者等同替换，均在申请待批的本发明的权利要求保护范围之内。
- 还没有人留言评论。精彩留言会获得点赞！