专利名称：融合空间信息及光谱信息的遥感图像的分类方法及装置的制作方法
技术领域：
本发明涉及遥感图像的分类方法及装置，更具体地说，涉及一种融合空间信息及光谱信息的遥感图像的分类方法及装置。
背景技术：
具有高光谱分辨率的高光谱图像具有广泛的应用领域，例如矿物探测、污染监测、精准农业等。其中分类是高光谱图像的一个非常重要的处理过程。传统的高光谱图像的分类主要是基于光谱信息，而很少考虑空间信息。普遍使用的融合空间信息及光谱信息的高光谱图像的分类方法包括马尔可夫随机场、形态学指针和形态学滤波。然而这些方法使用固定的窗口得到空间信息，导致了尺寸比例的选择问题。另一个分类方法是将空间信息集成到多核学习方法中，该方法也存在尺寸比例的选择问题，同时该方法的计算成本非常高。再有一种融合空间信息的方法是将逐 像素的分类结果和划分聚类(Partitional Clustering)得到的分割的结果融合起来,但这种方法的结果是不稳健的。
发明内容
本发明要解决的技术问题在于，针对现有技术的上述缺陷，提供一种稳健的、可行的、高精度的融合空间信息及光谱信息的遥感图像的分类方法及装置。本发明解决其技术问题所采用的技术方案是构造一种融合空间信息及光谱信息的遥感图像的分类方法，其包括以下步骤对所述遥感图像分别执行逐像素的支持向量机分类以得到分类图像，和基于区域的图像分割处理以得到分割图像；将所述分类图像和分割图像组合并执行大多数一致表决处理，以得到最终的分类图像。在本发明所述的融合空间信息及光谱信息的遥感图像的分类方法中，对所述遥感图像重复执行多次基于区域的图像分割处理后，得到所述分割图像。在本发明所述的融合空间信息及光谱信息的遥感图像的分类方法中，所述基于区域的图像分割处理为统计区域合并处理。在本发明所述的融合空间信息及光谱信息的遥感图像的分类方法中，对所述遥感图像重复执行多次统计区域合并处理后，得到所述分割图像。在本发明所述的融合空间信息及光谱信息的遥感图像的分类方法中，从所述遥感图像中选择三个光谱波段作为所述{R，G，B}颜色信道值，以执行所述统计区域合并处理。在本发明所述的融合空间信息及光谱信息的遥感图像的分类方法中，通过对所述遥感图像执行主成分变换，选择第一三主成分作为{R，G，B}颜色信道值。在本发明所述的融合空间信息及光谱信息的遥感图像的分类方法中，基于W分布，随机选择三个光谱波段作为所述{R，G，B}颜色信道值。
在本发明所述的融合空间信息及光谱信息的遥感图像的分类方法中，所述W分布是均匀分布、Wlda分布或WJ9分布。在本发明所述的融合空间信息及光谱信息的遥感图像的分类方法中，在执行所述统计区域合并处理之前，对每个选择的光谱波段执行线型对比度拉伸处理。根据本发明的另一个方面，提供一种融合空间信息及光谱信息的遥感图像的分类装置，其包括分类模块，用于对所述遥感图像执行逐像素的支持向量机分类以得到分类图像；分割模块，用于对所述遥感图像执行基于区域的图像分割处理以得到分割图像；融合空间信息及光谱信息分类模块，用于将所述分类图像和分割图像组合并执行大多数一致表决处理，以得到光谱空间分类图像。实施本发明的融合空间信息及光谱信息的遥感图像的分类方法及装置，具有以下有益效果由于将基于区域的图像分割结果和逐像素的分类结果相结合获得遥感图像的分类结果具有更好的稳健性，计算成本较低，易于执行及不产生尺寸比例的选择间题。解决了基于划分聚类(Partitional Clustering)的融合空间信息及光谱信息的高光谱图像的分类的不稳健的缺陷。
下面将结合附图及实施例对本发明作进一步说明，附图中图I是本发明融合空间信息及光谱信息的遥感图像的分类装置的结构示意图；图2是本发明融合空间信息及光谱信息的遥感图像的分类方法的流程图； 图3a是执行统计区域合并处理的分割图像的示意图；图3b是执行逐像素的支持向量机分类的示意图；图3c是将分割图像和分类图组合后的示意图；图3d是执行大多数一致表决后最终的分类图的示意图；图4a是本发明例子I中执行逐像素的支持向量机得到的分类图；图4b是本发明例子I中将图4a和图5a组合后的融合空间及光谱信息后得到的最终的分类图；图5a是本发明例子I中将PCA变换应用于SRM的分割图；图5b是本发明例子I中均匀分布的光谱波段选择应用于SRM处理的分割图；图5c是本发明例子I中将Wuia分布的光谱波段选择应用于SRM处理的分割图；图5d是本发明例子I中将Wjs分布的光谱波段选择应用于SRM处理的分割图；图6a是本发明例子I中执行200次将均匀分布的光谱波段选择SRM处理并执行大多数一致性表决后的融合空间及光谱信息后得到的最终的分类图；图6b是本发明例子I中执行200次将Wuia分布的光谱波段选择SRM处理并执行大多数一致性表决后的融合空间及光谱信息后得到的最终的分类图；图6c是本发明例子I中执行200次将Wjs分布的光谱波段选择SRM处理并执行大多数一致性表决后的融合空间及光谱信息后得到的最终的分类图；图7a是本发明例子2中执行逐像素的支持向量机分类得到的分类图；图7b是本发明例子2中将图4a和图5a组合后的融合空间及光谱信息后得到的最终的分类图；图8a是本发明例子2中将PCA变换应用于SRM的分割图；图8b是本发明例子2中均匀分布的光谱波段选择应用于SRM处理的分割图；图8c是本发明例子2中将Wuia分布的光谱波 段选择应用于SRM处理的分割图；图8d是本发明例子2中将Wjs分布的光谱波段选择 应用于SRM处理的分割图；图9a是本发明例子2中执行200次将均匀分布的光谱波段选择SRM处理并执行大多数一致性表决后的融合空间及光谱信息后得到的最终的分类图；图9b是本发明例子2中执行200次将Wuia分布的光谱波段选择SRM处理并执行大多数一致性表决后的融合空间及光谱信息后得到的最终的分类图；图9c是本发明例子2中执行200次将Wjs分布的光谱波段选择SRM处理并执行大多数一致性表决后的融合空间及光谱信息后得到的最终的分类图。
具体实施例方式为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。如图I所示，在图I是本发明融合空间信息及光谱信息的遥感图像的分类装置中，其主要包括分类模块11、分割模块12和光谱空间分类模块13 ;其中，分类模块11和分割模块12是对待处理的遥感图像进行第一级处理，该遥感图像可以是高光谱图像或多光谱图像，高光谱图像可以是例如ROSIS、AVIRIS。多光谱图像可以是例如SPOT、MODIS、LANDSET等等。该分类模块11对遥感图像进行第一级分类，也就是对遥感图像执行逐像素的支持向量机分类(pixel-wise support vector machine)以得到分类图像,该分类图像是初次分类图像；同时，分割模块12将分割该遥感图像，也就是对该遥感图像执行基于区域的图像分割处理以得到分割图像；与对遥感图像进行分割聚合处理相比，该分割模块12所得到的分割图像的结果更加稳健。为了使得，该分类模块11和分割模块12能得到更好的处理结果，可在遥感图像输入分类模块11和分割模块12之前进行预处理，该预处理主要为对于遥感图像，去除吸水性和低信噪比的光谱波段。当然也可包括辐射校正和几何校正。可以理解的，为了使得最终的分类结果更稳健，分割模块12可对遥感图像重复执行预设次数的基于区域的图像分割处理，例如可以是100次内的基于区域的图像分割处理，当然也可执行大于100次的基于区域的图像分割处理。一旦得到初次分类图像和分割图像后，光谱空间分类模块13就可执行光谱空间分类，即光谱空间分类模块13将分类图像和分割图像组合并执行大多数一致表决(majority voting)处理，以得到光谱空间分类图像。图2示出了本发明融合空间信息及光谱信息的遥感图像的分类方法流程，该方法流程基于图I所示的系统结构，具体过程如下S21 :对遥感图像分别执行逐像素的支持向量机分类以得到分类图像，和基于区域的图像分割处理以得到分割图像；可以理解的，在该步骤中的分割处理和分类处理为平行处理的过程。该遥感图像可以是高光谱图像或多光谱图像，高光谱图像可以是例如R0SIS、AVIRIS。多光谱图像可以是例如SP0T、M0DIS、LANDSET等等。为了使得，该分类和分割能得到更好的处理结果，可在遥感图像进行类和分割之前进行预处理，该预处理主要为对于遥感图像，去除吸水性和低信噪比的光谱波段。当然也可包括辐射校正和几何校正。可以理解的，为了使得最终的分类结果更稳健，可对遥感图像重复执行预设次数的基于区域的图像分割处理，例如可以是100次内的基于区域的图像分割处理，当然也可执行大于100次的基于区域的图像分割处理。优选的，该基于区域的图像分割处理为统计区域合并(SRM, StatisticalRegionMerging)处理。对经预处理的遥感图像重复执行多次统计区域合并处理后，由此可更快速和稳健地得到分割图像。对于遥感图像，可将其定义为I，由此该遥感图像I可包括I I I像素。该遥感图像I包括{R，G，B}颜色信道值，该{R，G，B}颜色信道值属于集合{0，2. . . g}，g = 255。将该遥感图像I的最优图像定义为I*。将该遥感图像的颜色信道采样自家族Q，取值范围为
，并分布在I*的每个像素。由此，Q的值将控制区域的数量，也就是Q的值越大，生成的区域数量越多。 在统计区域合并处理中有两点很重要，即合并条件和合并阶数。合并条件为吵，允)=[順
[false otherwise其中，i 1)= - Ir + -j— In —R和R'表示该遥感图像I的一组固定值，如果P(R，R' )= true，合并R和R'。合并阶数为选择函数f对分类该遥感图像I中的像素对。可将函数f选择为fa(p,p/ ) = |pa-pa'其中，Pe^P Pa'为像素信道值。为了改善统计区域合并处理的分割结果，本发明从遥感图像中选择三个光谱波段例如高光谱波段作为{R，G，B}颜色信道值，以执行统计区域合并处理。在一优选实施例中，通过对遥感图像执行主成分分析(PCA,principle component analysis)变换,选择第一三主成分(first three principle component)作为该{R,G,B}颜色信道值。该选择方法是静态的，其选择了三个预定的光谱波段。在另一优选实施例中，基于W分布，随机选择三个光谱波段作为{R，G，B}颜色信道值，该类选择方法为动态方法。具体处理时，该W分布是均匀分布、Wuia分布或WJ9分布。在采用均匀分布时，原始的遥感图像的光谱波段的权重是均等的。基于均匀分布，随机选择三个光谱波段作为{R，G，B}颜色信道值，以执行统计区域合并。在采用Wuia分布时，原始的遥感图像的光谱波段的权重是非均等的。可采用文献Martinez, A. M. andA. C. Kak, PCA versus LDA. IEEE Transactions on Pattern Analysis and MachineIntelligence, 2001. 23(2) p. 228-233.公开的Fisher LDA为每个光谱波段分配权重。通过下述公式计算每个光谱波段的权重
权利要求
1.一种融合空间信息及光谱信息的遥感图像的分类方法，其特征在于，包括以下步骤 对所述遥感图像分别执行逐像素的支持向量机分类以得到分类图像，和基于区域的图像分割处理以得到分割图像； 将所述分类图像和分割图像组合并执行大多数一致表决处理，以得到最终的分类图像。
2.根据权利要求I所述的融合空间信息及光谱信息的遥感图像的分类方法，其特征在于，对所述遥感图像重复执行多次基于区域的图像分割处理后，得到所述分割图像。
3.根据权利要求I所述的融合空间信息及光谱信息的遥感图像的分类方法，其特征在于，所述基于区域的图像分割处理为统计区域合并处理。
4.根据权利要求3所述的融合空间信息及光谱信息的遥感图像的分类方法，其特征在于，对所述遥感图像重复执行多次统计区域合并处理后，得到所述分割图像。
5.根据权利要求3所述的融合空间信息及光谱信息的遥感图像的分类方法，其特征在于，从所述遥感图像中选择三个光谱波段作为{R，G，B}颜色信道值，以执行所述统计区域合并处理。
6.根据权利要求5所述的融合空间信息及光谱信息的遥感图像的分类方法，其特征在于，通过对所述遥感图像执行主成分分析变换，选择第一三主成分作为所述{R，G，B}颜色信道值。
7.根据权利要求5所述的融合空间信息及光谱信息的遥感图像的分类方法，其特征在于，基于W分布,随机选择三个光谱波段作为所述{R，G, B}颜色信道值。
8.根据权利要求7所述的融合空间信息及光谱信息的遥感图像的分类方法，其特征在于，所述W分布是均匀分布、Wuia分布或WJ9分布。
9.根据权利要求5所述的融合空间信息及光谱信息的遥感图像的分类方法，其特征在于，在执行所述统计区域合并处理之前，对每个选择的光谱波段执行线型对比度拉伸处理。
10.一种融合空间信息及光谱信息的遥感图像的分类装置，其特征在于，包括 分类模块，用于对所述遥感图像执行逐像素的支持向量机分类以得到分类图像； 分割模块，用于对所述遥感图像执行基于区域的图像分割处理以得到分割图像； 融合空间信息及光谱信息的分类模块，用于将所述分类图像和分割图像组合并执行大多数一致表决处理，以得到最终的分类图像。
全文摘要
本发明涉及一种融合空间信息及光谱信息的遥感图像的分类方法及装置，该方法包括对所述遥感图像分别执行逐像素的支持向量机分类以得到分类图像，和基于区域的图像分割处理以得到分割图像；将所述分类图像和分割图像组合并执行大多数一致表决处理，以得到最终的分类图像。将基于区域的图像分割结果和逐像素的分类结果相结合获得遥感图像的分类结果具有更好的稳健性，计算成本较低，易于执行及不产生尺寸比例的选择问题。解决了传统的融合基于划分聚类(Partitional Clustering)分割的高光谱图像的融合空间及光谱信息分类的不稳健的缺陷。
文档编号G06K9/62GK102708373SQ201210003388
公开日2012年10月3日 申请日期2012年1月6日 优先权日2012年1月6日
发明者史文中, 苗则朗 申请人:香港理工大学