基于LBP特征的超像素空谱多核高光谱图像分类方法与流程

本发明涉及高光谱遥感图像处理
技术领域：
，尤其涉及一种基于lbp(localbinarypattern，局部二值模式)特征的超像素空谱多核高光谱图像分类方法。
背景技术：
：高光谱遥感的出现是遥感界的一场革命，它使本来在宽波段遥感中不可探测的物质，在高光谱遥感中能被探测。高光谱图像可以获得数百条相同的窄带光谱通道，并能提供更丰富的光谱信息，以支持各种地表覆盖材料的精细识别。因此，高光谱图像受到越来越多的关注，它被应用到分类、目标检测、异常检测、光谱解混等方面。而其中的高光谱图像分类任务又对地质勘探农作物检测、国防军事等领域起着实质性的重要作用，值得更加深入的研究。然而，在高光谱图像分类广泛应用的同时，它也面临着巨大的挑战，例如著名的休斯(hughes)现象。hughes现象是指在高光谱分析过程中，随着参与运算波段数目的增加，分类精度“先增后降”的现象。为了解决这一问题，特征提取被看作是高光谱图像处理的一个关键步骤。然而，由于光谱特征的空间变异性，高光谱图像特征提取是高光谱图像处理中最具挑战性的任务之一。最常用的高光谱图像特征提取方法是核方法，因其采用简单的线性加权方式，可以实现空间信息和光谱信息的联合学习与利用。核方法是指为了充分挖掘和利用高光谱数据蕴含的光谱信息和空间信息，构造新型核函数，实现样本相似性度量意义下不同异构特征融合，进而在svm框架下完成分类，以有效提高分类性能。从近年研究进展来看，针对高光谱核函数的设计经历了光谱加权核、空谱混合核、多核学习三个阶段。与仅考虑光谱信息的分类方法相比，结合空间信息和光谱信息的分类方法能够取得更好的分类效果。然而，由于空间核所采用的空间区域的大小和形状是固定的，常常导致高光谱图像的空间纹理无法得到充分的利用。超像素分割能够根据高光谱图像的空间结构自适应的改变区域的形状，并且利用像素之间特征的相似性将像素分组，用少量的超像素代替大量的像素来表达图片特征，能够很大程度上降低图像后处理的复杂度。基于超像素的复合核方法将超像素看作一个局部邻域，可以利用超像素来获得空间信息，从而避免了最佳空间邻域的选择。基于超像素的多核方法利用多个核函数来有效利用超像素的光谱空间信息，相比较基于超像素的单核方法，基于超像素的多核方法不仅利用了超像素内的空间信息，也利用了超像素间的空间信息，分类精度更高。然而，基于超像素的方法也存在缺陷：超像素内空间信息用像素的均值来表示，会丢失超像素内像素的边缘信息。技术实现要素：本发明提供了基于lbp特征的超像素空谱多核高光谱图像分类方法，以解决现有的基于超像素的方法的超像素内空间信息用像素的均值来表示，会丢失超像素内的边缘信息，导致图像分类的效率低、精度差等问题。本发明提供的基于lbp特征的超像素空谱多核高光谱图像分类方法，包括以下步骤：步骤1，采用熵率超像素分割算法对高光谱图像进行索引，生成带有超像素索引的高光谱图像；步骤2：利用所述带有超像素索引的高光谱图像获得三个核函数，所述三个核函数包括光谱核、基于超像素内的空间核和基于超像素间的空间核；步骤3：将得到的所述三个核函数进行融合，将获得的融合结果输入支持向量机分类器，完成训练后进行分类，得到分类结果图。进一步地，在一种实现方式中，所述步骤1包括：步骤1-1，使用主成分分析法对所述高光谱图像进行波段选择，获得三个各不相关的主成分图像；步骤1-2，对所述主成分图像应用熵率超像素分割算法，生成超像素分割图像；步骤1-3，将所述超像素分割图像与原始的高光谱图像结合，生成所述带有超像素索引的高光谱图像。进一步地，在一种实现方式中，所述步骤1-1包括：使用主成分分析法将所述高光谱图像线性投影到一组新的坐标空间中，选取所述高光谱图像的前三个主成分形成主成分图像，用于所述步骤1-2中的超像素分割。进一步地，在一种实现方式中，所述步骤1-2包括：步骤1-2-1，根据以下公式，结合所述主成分图像纹理的复杂程度选择超像素数：l＝lbase×rtexturertexture＝n/n其中，l是超像素数，lbase是基超像素数，rtexture是纹理比率，n表示经过纹理分析过滤后的主成分图像中的非零元素数量，n表示纹理分析过滤前主成分图像中非零元素的数量；步骤1-2-2，超像素分割首先在主成分图像上构造一个图g＝(v,e)；其中，v为基本图像像素对应的顶点集，e为相邻像素之间成对相似点的边缘集；步骤1-2-3，通过选择边的子集，将图分成多个连通的子图，每个所述子图对应一个超像素；步骤1-2-4，根据以下公式，在超像素分割的目标函数中加入熵率项h(a)和平衡项b(·)：其中，λ是为了控制h(a)和b(·)而引入的权重，λ≥0。进一步地，在一种实现方式中，所述步骤2包括：步骤2-1，采集所述高光谱图像中所有的光谱像素；单个超像素是由具有相似结构特征的相邻光谱像素组成，对所述光谱像素进行采集，获得所述高光谱图像的光谱特征，根据所述光谱特征计算获得光谱核；步骤2-2，利用局部二值模式算法，通过阈值标记法得出中心点像素与所述中心点像素的邻域像素之间的差分值，对所述主成分图像的局部纹理结构进行分析，实现超像素内空间核的提取；步骤2-3，采用加权平均算法提取所述超像素间的空间特征，在每个所述超像素中执行权值替换操作，将替换得到的所有所述超像素进行组合，得到所述超像素间的空间核。进一步地，在一种实现方式中，所述步骤2-1包括：步骤2-1-1，采集所述高光谱图像中所有的光谱像素，所述高光谱图像中所有的光谱像素构成高光谱图像的光谱特征；步骤2-1-2，利用采集到的所述光谱像素，计算光谱核：其中，为光谱核，为输入的n个光谱像素，σ为控制核测度光滑性的尺度参数，σ∈正实数。进一步地，在一种实现方式中，所述步骤2-2包括：步骤2-2-1，通过阈值标记法，根据中心点像素与中心点像素的邻域像素之间的差分值，对所述中心点像素的邻域像素进行标记；步骤2-2-2，完成标记后，利用局部二值模式算法，获取所有像素的lbp编码，实现对图像局部纹理结构的分析，提取超像素内的空间核所述局部二值模式算法的数学表达式为：其中，lbp(c、p)表示lbp编码，c表示中心点像素，p表示中心点像素周围的像素个数，ti表示相邻像素灰度值，tc表示中心点像素的灰度值，s(·)为符号函数。进一步地，在一种实现方式中，所述步骤2-2-1包括：在一个圆形区域内，将多个相邻的像素与中心点像素进行比较，如果任一所述相邻的像素的灰度值大于中心点像素的灰度值，则将所述相邻的像素标记为1，如果所述相邻的像素的灰度值小于中心点像素的灰度值，则将所述相邻的像素标记为0。进一步地，在一种实现方式中，所述步骤2-3包括：步骤2-3-1，采用加权平均算法提取超像素间的空间特征，即对于超像素xi，根据以下公式，计算加权平均像素：其中，为加权平均像素，为超像素xij的所有像素的平均值；其中，的权值ωi,j根据以下公式求出：其中，为超像素xi的所有像素的平均值，h为核尺度参数；步骤2-3-1，在每个所述超像素中都执行相同的权值替换操作；将替换得到的所有超像素进行组合，获得基于超像素间的空间核进一步地，在一种实现方式中，所述步骤3包括：步骤3-1，根据以下公式，将所述步骤2得到的三个核函数进行融合，形成复合核：其中，为光谱核，为基于超像素内的空间核，为基于超像素间的空间核；步骤3-2，将所述复合核输入支持向量机分类器，完成训练后进行分类，得到分类结果图。本发明提出了一种基于lbp特征的超像素空谱多核高光谱图像分类方法。首先，将经过主成分分析法降维的图像，用超像素分割生成带有超像素索引的高光谱图像。然后，采用加权平均滤波和lbp算法提取超像素间和超像素内的空间特征得到超像素间的空间核和超像素内的空间核，联合已提取到的光谱核进行融合。最后，将这个组合核输入到支持向量机分类器中生成分类结果图。本发明将lbp算法与超像素结合，用lbp算法来提取超像素内的边缘特征信息，可以有效的解决传统的多核方法用超像素内所有像素的平均值来代替超像素内所有像素值导致的像素边缘信息丢失以及边缘像素类别划分不准确的问题，实现分类精度与效率的提升。本发明与现有技术相比，其显著优点在于：(1)本发明在对高光谱图像完成主成分波段选择后，利用了熵率超像素算法生成二维超像素图像，可以大幅提高核提取的效率与精度。(2)本发明在采用超像素分割进行光谱核提取的基础上，采用lbp算法提取超像素内的空间核以及采用加权平均算法来提取超像素间的空间核。三种核函数各司其职，既能有效提高分类的精确度，又能提高分类的效率。(3)融合三种核函数的分类方法快速有效，方法巧妙，构思新颖，具有良好的应用前景。附图说明为了更清楚地说明本发明的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1是本发明实施例部分提供的基于lbp特征的超像素空谱多核高光谱图像分类方法的工作流程示意图；图2是本发明实施例部分提供的基于lbp特征的超像素空谱多核高光谱图像分类方法中对得到的超像素分割图与原始的高光谱图像结合生成的带有超像素索引的高光谱图像的工作流程示意图。具体实施方式为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。本发明实施例公开了基于lbp特征的超像素空谱多核高光谱图像分类方法，本方法应用于高光谱遥感图像的精确分类，是高光谱遥感影像预处理的关键步骤，对地质勘探、农作物检测、国防军事等应用领域起着实质性的关键作用。如图1所示，本实施例提供的基于lbp特征的超像素空谱多核高光谱图像分类方法，包括以下步骤：步骤1，采用熵率超像素分割算法对高光谱图像进行索引，生成带有超像素索引的高光谱图像；具体的，本步骤中，首先使用主成分分析法对高光谱图像进行波段选择。接着在基础图像上应用熵率超像素算法生成二维超像素图像。最后与原始高光谱图像结合生成带有超像素索引的高光谱图像。步骤2：利用所述带有超像素索引的高光谱图像获得三个核函数，所述三个核函数包括光谱核、基于超像素内的空间核和基于超像素间的空间核；具体的，本步骤中，首先根据每个超像素都是由一组相邻的光谱像素组成的原理，对高光谱图像中所有的光谱像素进行采集，光谱信息可以直接由采集的光谱像素构成，以此获取高光谱图像的光谱特征。其次利用lbp算法，通过阈值来标记中心点像素与其邻域像素之间的差别，从而对局部纹理结构进行分析，实现超像素内的空间核的提取。最后采用加权平均算法来提取超像素间的空间特征，在每个超像素中都执行权值替换操作，将替换得到的所有超像素进行组合，进而得到超像素间的空间核。步骤3：将得到的所述三个核函数进行融合，将获得的融合结果输入支持向量机分类器，完成训练后进行分类，得到分类结果图。本发明将lbp算法与超像素结合，用lbp算法来提取超像素内的边缘特征信息，可以有效的解决传统的多核方法用超像素内所有像素的平均值来代替超像素内所有像素值导致的像素边缘信息丢失以及边缘像素类别划分不准确的问题，实现分类精度与效率的提升。本实施例所述的一种基于lbp特征的超像素空谱多核高光谱图像分类方法中，所述步骤1包括：步骤1-1，使用主成分分析法对所述高光谱图像进行波段选择，获得三个各不相关的主成分图像；步骤1-2，对所述主成分图像应用熵率超像素分割算法，生成超像素分割图像；步骤1-3，将所述超像素分割图像与原始的高光谱图像结合，生成所述带有超像素索引的高光谱图像。本实施例所述的一种基于lbp特征的超像素空谱多核高光谱图像分类方法中，所述步骤1-1包括：使用主成分分析法将所述高光谱图像线性投影到一组新的坐标空间中，选取所述高光谱图像的前三个主成分形成主成分图像，用于所述步骤1-2中的超像素分割。本实施例所述的一种基于lbp特征的超像素空谱多核高光谱图像分类方法中，所述步骤1-2包括：步骤1-2-1，根据以下公式，结合所述主成分图像纹理的复杂程度选择超像素数：l＝lbase×rtexturertexture＝n/n其中，l是超像素数，lbase是基超像素数，rtexture是纹理比率，n表示经过纹理分析过滤后的主成分图像中的非零元素数量，n表示纹理分析过滤前主成分图像中非零元素的数量；步骤1-2-2，超像素分割首先在主成分图像上构造一个图g＝(v,e)；其中，v为基本图像像素对应的顶点集，e为相邻像素之间成对相似点的边缘集；步骤1-2-3，通过选择边的子集，将图分成多个连通的子图，每个所述子图对应一个超像素；步骤1-2-4，根据以下公式，在超像素分割的目标函数中加入熵率项h(a)和平衡项b(·)：其中，λ是为了控制h(a)和b(·)而引入的权重，λ≥0。如图2所示，本实施例所述的步骤2和步骤3包含对带有超像素索引的高光谱图像的核提取与融合。其中步骤2包含带有超像素索引的高光谱图像获得三个核函数实现三种不同类型的特征提取。本实施例所述的一种基于lbp特征的超像素空谱多核高光谱图像分类方法中，所述步骤2包括：步骤2-1，采集所述高光谱图像中所有的光谱像素；单个超像素是由具有相似结构特征的相邻光谱像素组成，对所述光谱像素进行采集，获得所述高光谱图像的光谱特征，根据所述光谱特征计算获得光谱核；步骤2-2，利用局部二值模式算法，通过阈值标记法得出中心点像素与所述中心点像素的邻域像素之间的差分值，对所述主成分图像的局部纹理结构进行分析，实现超像素内空间核的提取；步骤2-3，采用加权平均算法提取所述超像素间的空间特征，在每个所述超像素中执行权值替换操作，将替换得到的所有所述超像素进行组合，得到所述超像素间的空间核。本实施例所述的一种基于lbp特征的超像素空谱多核高光谱图像分类方法中，所述步骤2-1包括：步骤2-1-1，采集所述高光谱图像中所有的光谱像素，所述高光谱图像中所有的光谱像素构成高光谱图像的光谱特征；步骤2-1-2，利用采集到的所述光谱像素，计算光谱核：其中，为光谱核，为输入的n个光谱像素，σ为控制核测度光滑性的尺度参数，σ∈正实数。本实施例所述的一种基于lbp特征的超像素空谱多核高光谱图像分类方法中，所述步骤2-2包括：步骤2-2-1，通过阈值标记法，根据中心点像素与中心点像素的邻域像素之间的差分值，对所述中心点像素的邻域像素进行标记；步骤2-2-2，完成标记后，利用局部二值模式算法，获取所有像素的lbp编码，实现对图像局部纹理结构的分析，提取超像素内的空间核所述局部二值模式算法的数学表达式为：其中，lbp(c、p)表示lbp编码，c表示中心点像素，p表示中心点像素周围的像素个数，ti表示相邻像素灰度值，tc表示中心点像素的灰度值，s(·)为符号函数。本实施例所述的一种基于lbp特征的超像素空谱多核高光谱图像分类方法中，所述步骤2-2-1包括：在一个圆形区域内，将多个相邻的像素与中心点像素进行比较，如果任一所述相邻的像素的灰度值大于中心点像素的灰度值，则将所述相邻的像素标记为1，如果所述相邻的像素的灰度值小于中心点像素的灰度值，则将所述相邻的像素标记为0。本实施例所述的步骤2-2利用lbp算法中，通过阈值来标记中心点像素与其邻域像素之间的差别，从而对高光谱图像的局部纹理结构进行分析，此方法渐变且不受光照条件变化的影响。原始的高光谱图像经过主成分分析后保留了前3个主成分，在每个主成分上应用lbp算法提取特征，将得到的局部二值直方图的所有波段串联起来，形成空间特征向量。本实施例所述的一种基于lbp特征的超像素空谱多核高光谱图像分类方法中，所述步骤2-3包括：步骤2-3-1，采用加权平均算法提取超像素间的空间特征，即对于超像素xi，根据以下公式，计算加权平均像素：其中，为加权平均像素，为超像素xij的所有像素的平均值；其中，的权值ωi,j根据以下公式求出：其中，为超像素xi的所有像素的平均值，h为核尺度参数；步骤2-3-1，在每个所述超像素中都执行相同的权值替换操作；将替换得到的所有超像素进行组合，获得基于超像素间的空间核本实施例所述的一种基于lbp特征的超像素空谱多核高光谱图像分类方法中，所述步骤3包括：步骤3-1，根据以下公式，将所述步骤2得到的三个核函数进行融合，形成复合核：其中，为光谱核，为基于超像素内的空间核，为基于超像素间的空间核；步骤3-2，将所述复合核输入支持向量机分类器，完成训练后进行分类，得到分类结果图。将本发明在公开的indianpines高光谱图像数据集上进行实施，与近年知名的高光谱图像分类方法进行定量对比，其分类结果如表1所示。表1本发明在indianpines数据上进行实施所得分类结果类名svmckspatvlbpelmsc_mkrmkasmgssk本发明方法alfalfa0.52860.85100.98540.96100.95510.98780.9756corn-notill0.93520.97090.97030.94160.96600.97550.9865corn-mintill0.94050.98390.96560.97340.97310.99080.9969corn0.84330.99480.96890.97040.97100.96240.9948grass/pasture0.95080.95810.98330.96670.96440.98440.9837grass/tree0.98120.98490.98070.99820.98980.99110.9998grass/pasture-mowed0.16960.50000.93640.96400.96520.96520.9640hay-windrowed0.99071.00000.99500.98330.99370.99771.0000oats0.06670.05560.92800.91110.98330.90561.0000soybeans-notill0.89390.94370.92270.94100.96620.97770.9859soybeans-mintill0.96120.99040.98910.96980.99580.99140.9999soybeans-cleantill0.92180.97240.98310.94530.96680.98140.9919wheat0.99010.99530.99010.99460.99630.99010.9946woods0.97440.99440.99400.99380.99680.99380.9999bldg-grass-tree-drives0.81850.98450.99520.94670.99120.99091.0000stone-steeltowers0.89300.78260.92530.97140.98600.94190.9774mean_oa0.93750.97510.98160.96700.98230.98610.9946std_oa0.00260.00690.00450.00520.00300.00180.0017mean_kappa0.92860.87260.97890.96230.97990.98420.9939std_kappa0.00300.01350.00320.00600.00340.00200.0019本实施例中，总体分类精度(oa)和kappa系数以及它们的标准差被作为评价指标。表中的最优结果以粗体进行表示。可见，本发明提出的分类方法与其他的几种分类方法相比，精度提高了1％-6％，并且本发明的分类精度标准差与其他几种分类算法相比也是最小的。这说明在训练样本有限的情况下，本发明设计的算法不仅具有更好的分类性能，而且分类性能也更加稳定。本发明与现有技术相比，其显著优点在于：(1)本发明在对高光谱图像完成主成分波段选择后，利用了熵率超像素算法生成二维超像素图像，可以大幅提高核提取的效率与精度。(2)本发明在采用超像素分割进行光谱核提取的基础上，采用lbp算法提取超像素内的空间核以及采用加权平均算法来提取超像素间的空间核。三种核函数各司其职，既能有效提高分类的精确度，又能提高分类的效率。(3)融合三种核函数的分类方法快速有效，方法巧妙，构思新颖，具有良好的应用前景。具体实现中，本发明还提供一种计算机存储介质，其中，该计算机存储介质可存储有程序，该程序执行时可包括本发明提供的基于lbp特征的超像素空谱多核高光谱图像分类方法的各实施例中的部分或全部步骤。所述的存储介质可为磁碟、光盘、只读存储记忆体(read-onlymemory，rom)或随机存储记忆体(randomaccessmemory，ram)等。本领域的技术人员可以清楚地了解到本发明实施例中的技术可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本发明实施例中的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如rom/ram、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。本说明书中各个实施例之间相同相似的部分互相参见即可。以上所述的本发明实施方式并不构成对本发明保护范围的限定。当前第1页12