一种高光谱图像的降维方法及系统技术方案

本发明专利技术公开了一种高光谱图像的降维方法及系统，该方法包括：根据输入的高光谱图像样本集D＝{x

A dimension reduction method and system of hyperspectral image

【技术实现步骤摘要】
一种高光谱图像的降维方法及系统
本专利技术涉及高光谱图像处理
，特别涉及一种高光谱图像的降维方法及系统。
技术介绍
高光谱图像的降维处理是高光谱图像处理领域中一项重要的研究课题。众所周知，高光谱图像是一种三维图像数据立方体，它具有丰富的光谱信息和空间信息，可以更好地分析不同地物间的细微差异。但原始高光谱图像包含成百上千个波段，数据量大，其中含有大量的冗余信息，这为后续的目标异常检测等工作带来很大的困难和影响。目前，针对高光谱图像的降维处理问题，常规的处理方法是波段选择和特征提取。波段选择法选择出来的波段之间相关性较弱，能较好地保留图像中的信息量，但相比较特征提取法，处理后的图像携带原始信号的能量较少。现有算法还通过找到所有样本全局的最优解的思想来实现对高光谱图像进行降维，当数据量很大时或者样本维度很高时，计算量非常大，运算时间较长；对原高光谱图像中众多具有一定相关性的信息重新组合成一组新的相互无关的信息来代替原来的信息，但其得到的降维结果无法反映出样本点之间所隐藏的非线性性质，对保留主分量个数的估计依然比较困难。
技术实现思路
本专利技术针对上述现有技术中存在的问题，提出一种高光谱图像的降维方法及系统，能够通过较小的运算复杂度对高光谱图像进行降维，且处理效果理想。为解决上述技术问题，本专利技术是通过如下技术方案实现的：本专利技术提供一种高光谱图像的降维方法，其包括：S11：根据输入的高光谱图像样本集D＝{x1,x2,...xm}、最近邻数k以及降维到的维数d，计算和xi最近的k个最近邻(xi1,xi2,...xik)；S12：计算高光谱图像样本集D的局部协方差矩阵Zi；S13：根据所述S12中计算出的局部协方差矩阵Zi计算出其权重系数向量Wi；S14：根据所述S13计算出的权重系数向量Wi组成权重系数矩阵W，根据W计算得到矩阵：M＝(I-W)(I-W)T；S15：根据矩阵M计算矩阵M前d+1个特征值及这d+1个特征值对应的特征向量{y1,y2,...yd+1}，根据第二个特征向量到第d+1个特征向量，组成的矩阵即为输出低维高光谱图像样本集矩阵D′＝{y2,y3,...yd+1}。较佳地，所述S11中的计算和xi最近的k个最近邻(xi1,xi2,...xik)具体为：按欧式距离作为度量，计算和xi最近的k个最近邻(xi1,xi2,...xik)。较佳地，所述S12中的局部协方差矩阵Zi的公式为：Zi＝(xi-xj)(xi-xj)T。较佳地，所述S13中的权重系数向量的计算方法为：求最小化损失函数对应的m×k维度的权重系数W，将那些不在邻域位置的W的位置取值为0，将W扩充到m×m维度；所述最小化损失函数的公式为：较佳地，所述S13中的权重系数向量Wi的公式为：其中，1k表示k维全1向量。本专利技术还提供一种高光谱图像的降维系统，其包括：最近邻计算模块、局部协方差矩阵计算模块、权重系数向量计算模块、矩阵M计算模块以及低位高光谱图像样本集矩阵计算模块；其中，所述最近邻计算模块用于根据输入的高光谱图像样本集D＝{x1,x2,...xm}、最近邻数k以及降维到的维数d，计算和xi最近的k个最近邻(xi1,xi2,...xik)；所述局部协方差矩阵计算模块用于计算高光谱图像样本集D的局部协方差矩阵Zi；所述权重系数向量计算模块用于根据所述局部协方差矩阵计算模块计算出的局部协方差矩阵Zi计算出其权重系数向量Wi；所述矩阵M计算模块用于根据所述权重系数向量计算模块计算出的权重系数向量Wi组成权重系数矩阵W，根据W计算得到矩阵：M＝(I-W)(I-W)T；所述低位高光谱图像样本集矩阵计算模块用于根据矩阵M计算矩阵M前d+1个特征值及这d+1个特征值对应的特征向量{y1,y2,...yd+1}，根据第二个特征向量到第d+1个特征向量，组成的矩阵即为输出低维高光谱图像样本集矩阵D′＝{y2,y3,...yd+1}。较佳地，所述最近邻计算模块中的计算和xi最近的k个最近邻(xi1,xi2,...xik)具体为：按欧式距离作为度量，计算和xi最近的k个最近邻(xi1,xi2,...xik)。较佳地，所述局部协方差矩阵计算模块中的局部协方差矩阵Zi的公式为：Zi＝(xi-xj)(xi-xj)T。较佳地，所述权重系数向量计算模块中的权重系数向量的计算方法为：求最小化损失函数对应的m×k维度的权重系数W，将那些不在邻域位置的W的位置取值为0，将W扩充到m×m维度；所述最小化损失函数的公式为：较佳地，所述权重系数向量计算模块中的权重系数向量Wi的公式为：其中，1k表示k维全1向量。相较于现有技术，本专利技术具有以下优点：(1)本专利技术提供的高光谱图像的降维方法及系统，通过计算矩阵M的前d+1个特征值及这d+1个特征值对应的特征向量{y1,y2,...yd+1}，并将第二个特征向量到第d+1个特征向量组成矩阵，能够通过较小的运算复杂度对高光谱图像进行降维，且处理效果理想；(2)本专利技术提供的高光谱图像的降维方法及系统，通过计算局部协方差矩阵Zi以及其权重系数向量Wi，使得均方差损失函数最小，从而在降维后，仍能保持很好的线性关系。当然，实施本专利技术的任一产品并不一定需要同时达到以上所述的所有优点。附图说明下面结合附图对本专利技术的实施方式作进一步说明：图1为本专利技术的实施例的高光谱图像的降维方法的流程示意图；图2为本专利技术的实施例的AVIRIS高光谱图像中机场区域降维后的波段1图像；图3为本专利技术的实施例的AVIRIS高光谱图像中机场区域降维后的波段5图像；图4为本专利技术的实施例的AVIRIS高光谱图像中机场区域降维后的波段10图像；图5为本专利技术的实施例的HYDICE高光谱图像中Urban区域降维后的波段1图像；图6为本专利技术的实施例的HYDICE高光谱图像中Urban区域降维后的波段5图像；图7为本专利技术的实施例的HYDICE高光谱图像中Urban区域降维后的波段10图像。具体实施方式下面对本专利技术的实施例作详细说明，本实施例在以本专利技术技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本专利技术的保护范围不限于下述的实施例。如图1所示为本专利技术的实施例的高光谱图像的降维方法的流程示意图。请参考图1，本实施例的高光谱图像的降维方法包括以下步骤：S11：根据输入的高光谱图像样本集D＝{x1,x2,...xm}、最近邻数k以及降维到的维数d，计算和xi最近的k个最近邻(xi1,xi2,...xik)；S12：计算高光谱图像样本集D的局部协方差矩阵Zi；S13：根据S12中计算出的局部协方差矩阵Zi计算出其权重系数向量Wi；S14：根据S13计算出的权重系数向量Wi组成权重系数矩阵W，根据W计算得到矩阵：M＝(I本文档来自技高网...

【技术保护点】
1.一种高光谱图像的降维方法，其特征在于，包括：/nS11：根据输入的高光谱图像样本集D＝{x

【技术特征摘要】
1.一种高光谱图像的降维方法，其特征在于，包括：
S11：根据输入的高光谱图像样本集D＝{x1,x2,…xm}、最近邻数k以及降维到的维数d，计算和xi最近的k个最近邻(xi1,xi2,…xik)；
S12：计算高光谱图像样本集D的局部协方差矩阵Zi；
S13：根据所述S12中计算出的局部协方差矩阵Zi计算出其权重系数向量Wi；
S14：根据所述S13计算出的权重系数向量Wi组成权重系数矩阵W，根据W计算得到矩阵：M＝(I-W)(I-W)T；
S15：根据矩阵M计算矩阵M前d+1个特征值及这d+1个特征值对应的特征向量{y1,y2,…yd+1}，根据第二个特征向量到第d+1个特征向量，组成的矩阵即为输出低维高光谱图像样本集矩阵D′＝{y2,y3,…yd+1}。


2.根据权利要求1所述的高光谱图像的降维方法，其特征在于，所述S11中的计算和xi最近的k个最近邻(xi1,xi2,…xik)具体为：按欧式距离作为度量，计算和xi最近的k个最近邻(xi1,xi2,…xik)。


3.根据权利要求1所述的高光谱图像的降维方法，其特征在于，所述S12中的局部协方差矩阵Zi的公式为：Zi＝(xi-xj)(xi-xj)T。


4.根据权利要求1所述的高光谱图像的降维方法，其特征在于，所述S13中的权重系数向量的计算方法为：求最小化损失函数对应的m×k维度的权重系数W，将那些不在邻域位置的W的位置取值为0，将W扩充到m×m维度；所述最小化损失函数的公式为：


5.根据权利要求4所述的高光谱图像的降维方法，其特征在于，所述S13中的权重系数向量Wi的公式为：



其中，1k表示k维全1向量。


6.一种高光谱图像的降维系统，其特征在于，包括：最近邻计算模块、局部协方差矩阵计算模块、权重系数向量计算模块、矩阵M计算模块以及低位高光谱图像样本集矩...

【专利技术属性】
技术研发人员：张焕芹，张嘉伟，毛士杰，
申请(专利权)人：上海融军科技有限公司，
类型：发明
国别省市：上海;31