专利名称：基于像素空间信息的高光谱图像分割方法
技术领域：
本发明属于图像处理技术领域，特别是一种图像分割方法，该方法可以对复 杂地貌进行分析，从相近地貌中发现感兴趣的目标并大体划定地貌间的边界。
背景技术：
图像分割是指对自然图像或SAR图像进行处理，将图像中不同目标分离出来。 高光谱图像数据由高光谱遥感产生的，维数很高，它的每一维数据都是一幅灰度 图像。使用传统的图像分割方法无法对其进行整体处理，对数据的每一维进行处 理后再组合起来工作量太大，效果也很难保证。国外的一些学者提出了使用支持 向量机对高光谱图像进行分割，效果也较好。但支持向量机是一种有监督学习方 法，在进行训练时必须知道训练样本的类别属性，这一弱点限制了它在这一领域 的应用，因为在实际工作中往往不知道图像中各部分的类别。流形学习算法可以 在不知道数据类别的情况下降低数据维数并保持数据内部几何关系。这一特性使 得流形学习在高光谱图像分割领域有广泛的应用。2002年澳大利亚的一些学者 提出将流形学习中的等距映射算法应用到高光谱图像分割领域。在他们的工作中 用等距映射算法对图像进行降维，然后使用贝叶斯网络进行类别划分完成图像分 割工作。尽管他们的工作中只是挑选了像素特征差别大即高光谱图像中的每个像 素点整体的灰度信息的地貌进行试验性分割，但效果还是不很理想。
发明内容
本发明的在于克服上述已有技术的不足，提供一种基于像素空间信息的高光 谱图像分割方法，以提高图像分割效果。
实现本发明目的技术方案是将像素空间信息引入等距映射算法，将像素空间 信息与像素特征信息组合起来，找到一个从像素空间分布到像素特征分布的渐变 过程，以此找到地貌间的边界将高光谱图像分割，具体步骤如下
为实现上述目的，本发明包括如下步骤-
(1)对高光谱图像数据的像素特征信息和像素空间信息分别进行欧氏距离求 解，获得像素特征欧氏距离矩阵和像素空间欧氏距离矩阵，并对这两个矩阵进行 归一化；(2) 将归一化的像素特征欧氏距离矩阵和像素空间欧氏距离矩阵进行加权， 并通过调节加权参数改变二者的大小，使加权后的两个矩阵相加后形成联合相异 度矩阵，反复调整加权参数，得到多组联合相异度矩阵；
(3) 使用等距映射算法分别对每一组的联合相异度矩阵进行降维，得到多组 映射结果；
(4) 对多组映射结果分别进行欧氏距离求解，并统计像素点空间近邻变化情 况，将与近邻距离持续变大的点判定为边界点，完成高光谱图像的初步分割；
(5) 对图像初步分割结果中的边界点进行如下类别校正，得到最后的图像分 割结果，
5a)根据像素空间欧氏距离矩阵和像素特征欧氏距离矩阵提取每个边界点与 其像素空间八近邻点的像素特征欧氏距离，形成该边界点的距离集合；
5b)按每个边界点的八近邻点类别统计其距离集合，形成x个欧氏距离向量；
5c)计算每一个欧氏距离向量的中值，将x个中值中最小的那个向量代表的 类别设定为该边界点的最终类别。
本发明与现有技术相比具有以下优点-
1、 本发明由于将像素空间信息引入高光谱图像分割领域，增加了处理图像 时的信息量，与其它流形技术不同之处在于本方法依据从像素空间分布到像素特 征分布的渐进变化过程分割图像，相比于其它的无监督学习方法理论上更加完 备，效果也更好；
2、 本发明由于采用图像初次分割后，对边界点进行类别校正，因而减少了 错分的边界点，使最终得到的分割结果更接近于图像的地貌的真实分布。
图1是本发明的高光谱图像分割流程图； 图2是本发明待处理的高光谱数据类别示意图； 图3是使用等距映射算法对像素特征信息进行处理的映射结果图； 图4是本发明的高光谱图像分割结果图。
具体实施例方式
参照图1,本发明的具体实现步骤如下-
步骤1:求解高光谱图像数据的像素特征信息和像素空间信息的欧氏距离， 并进行归一化。1.1) 选取高光谱图像中图3所示的区域，该区域的数据包含200个波段，每 个波段都是一幅灰度图像，由此可以将这些数据视为200个灰度特征矩阵，按波 段顺序将每个灰度特征矩阵按照列分解，然后按顺序连接起来形成一个大的列向 量，由此得到了 200个列向量；
1.2) 将这200个列向量组合起来形成像素特征矩阵data,该特征矩阵data的 每一行都表示一个像素点的像素特征矢量，计算每两行像素特征矢量的欧氏距 离，得到像素特征欧氏距离矩阵Y1;
1.3) 按照1. 2)步中像素点的顺序获取每个像素点在图像上的空间位置矢量, 将所有点的空间位置矢量组成像素位置矩阵axis,该空间矩阵axis的每一行都表 示一个像素点的空间位置矢量，计算每两行像素空间位置矢量的欧氏距离，得到 像素空间欧氏距离矩阵Y2;
1.4) 得到像素特征欧氏距离矩阵Yl和像素空间欧氏距离矩阵Y2后，找到这 两个矩阵中的最大值，用Y1和Y2除以其本身的最大值，完成对这两个矩阵的 归一化。
步骤2:将归一化的像素特征欧氏距离矩阵和像素空间欧氏距离矩阵进行加 权，得到多组联合相异度矩阵。
得到归一化的像素特征欧氏距离矩阵Yl和像素空间欧氏距离矩阵Y2后， 设定一个加权参数;i,用A乘以像素特征欧氏距离矩阵Yl，用l-义乘以像素空 间欧氏距离矩阵Y2，将加权后的两个矩阵相加，得到联合相异度矩阵D，矩阵 D的具体形式是D-;^Y1 + (1-X"Y2，通过改变义的大小可以得到多组合相异 度矩阵，在本实验仿真中;i从0增加到0.5， 一次增加0.05，即保持两个矩阵的 加权参数合为l,先设定像素特征欧氏距离矩阵的加权参数初始值为0，然后将 该参数每次增加0.05，直到该参数增加到0.5，像素空间欧氏距离矩阵的加权参 数则相应由l减少到0.5。
步骤3:使用等距映射算法分别对每一组的联合相异度矩阵进行降维，得到 多组映射结果。
步骤4:对多组映射结果分别进行欧氏距离求解，并统计像素点空间近邻变 化情况，将与其近邻距离持续变大的点判定为边界点，完成高光谱图像的初步分 割。
4.1)将每个映射结果视为一个矩阵，矩阵的每一行表示一个像素点的映射坐
6标，计算每两行映射坐标的欧氏距离，得到每组映射结果的像素欧氏距离矩阵L;
4.2) 根据像素空间欧氏距离矩阵Y2取得每个像素点八近邻变化情况，即得 到像素空间欧氏距离矩阵Y2每一行中最小的八个非0值的列号，这些列号就是 该像素点的八近邻在像素欧氏距离矩阵L中的列号，将这些列号的数值组成该 行所代表的点的列表向量In，将这些列表向量In组成八近邻列表矩阵index, 该矩阵的每一行都代表一个像素点的列表向量；
4.3) 按照八近邻列表矩阵index在L中找到每个点的八近邻距离，将其中最 大三个相加得到一个距离和值，对每组映射结果进行相同处理，得到多个距离和 值，将这些距离和值按顺序组成该点的偏移行向量；
4.4) 将所有点的偏移行向量组成偏移矩阵dist，根据参数调整的顺序用dist 的后一个列向量减去前一个列向量，得到近邻变化矩阵L'， L'的每一行表示一 个像素点的近邻变化情况，即近邻变化行向量，统计近邻变化矩阵L'每一行中大 于0的数值个数，将这个个数作为该点的近邻变化值，由此得到变化描述列向量 S;
4.5) 针对dist设定一个阈值，査找dist每一行中大于这个阈值的数值个数， 将这个个数再加到变化描述列向量S的相应行的数字上；
4.6) 针对变化描述列向量S设定一个阈值m,将S各行中大于这个阈值m的
行所表示的点标记为边界点；
4.7) 将由4.6)步中标记的边界点隔开的数据块标记为不同的类别，得到初 步的图像分割结果W。
步骤5:对图像初步分割结果中的边界点进行类别校正，得到最后的图像分 割结果。
5.1) 根据像素空间欧氏距离矩阵和像素特征欧氏距离矩阵提取每个边界点与 其像素空间八近邻点的像素特征欧氏距离，即计算每个被标记为边界点的像素点 与其八近邻列表In中记录的像素点间的像素特征欧氏距离，形成该边界点的特 征距离集合；
5.2) 按每个边界点的八近邻点类别统计其距离集合，即依据步骤4得到的标 记类别统计5.1)得到的八个欧氏距离，如果近邻点属于同一类就将它们与该边 界点的欧氏距离值放在同一向量中，设该边界点的八近邻共分为x类，由此形成 x个向量；
75.3)计算5.2)中得到的每个向量的中值，将x个中值中中值最小的那个向 量的类别作为该边界点的最终类别，更新步骤4中得到的分割结果W得到更新 后的分割结果W'，对所有的边界点进行以上的处理，将最终得到的W'作为最后 的分割结果。
本发明的效果通过以下仿真进一步说明。
1. 仿真条件与内容
我们使用的实验数据是典型的AVIRIS高光谱图像。原始图像取自1992年6 月拍摄的美国印第安纳州西北部印第安遥感试验区的一部分，地貌类别共计16 类。该图像来源于Purdue大学遥感图像处理实验室。原始数据的波段是220，每 个波段图像的大小为145 x 145，每个像素16bits。在具体的实验中我们选取了噪 声污染少的200个波段。本发明中由于图像尺寸相对于等距映射算法来说过于巨 大，我们只能选取图像中一部分演示本方法的效果。
实验使用的软件平台为MATLAB7.0。
2. 仿真结果
图3是使用现有等距映射算法对像素特征信息进行处理的映射结果图，图4 是使用本发明的高光谱图像分割结果图。从图3和图4两幅图可以看出，直接使 用等距映射算法各类样本点间有很严重的混杂现象，很难进行类别划分，误分 率会很高；本发明中的方法能够将不同的地貌分割出来，混杂现象基本不存在， 分类效果要比直接使用等距映射算法好的多。
权利要求
1、一种高光谱图像分割方法，其特征在于包含步骤如下(1)对高光谱图像数据的像素特征信息和像素空间信息分别进行欧氏距离求解，获得像素特征欧氏距离矩阵和像素空间欧氏距离矩阵，并对这两个矩阵进行归一化；(2)将归一化的像素特征欧氏距离矩阵和像素空间欧氏距离矩阵进行加权，并通过调节加权参数改变二者的大小，使加权后的两个矩阵相加后形成联合相异度矩阵，反复调整加权参数，得到多组联合相异度矩阵；(3)使用等距映射算法分别对每一组的联合相异度矩阵进行降维，得到多组映射结果；(4)对多组映射结果分别进行欧氏距离求解，并统计像素点空间近邻变化情况，将与近邻距离持续变大的点判定为边界点，完成高光谱图像的初步分割；(5)对图像初步分割结果中的边界点进行如下类别校正，得到最后的图像分割结果，5a)根据像素空间欧氏距离矩阵和像素特征欧氏距离矩阵提取每个边界点与其像素空间八近邻点的像素特征欧氏距离，形成该边界点的特征距离集合；5b)按每个边界点的八近邻点类别统计其距离集合，形成x个欧氏距离向量；5c)计算每一个欧氏距离向量的中值，将x个中值中最小的那个向量代表的类别设定为该边界点的最终类别。
2、 根据权利要求1所述的高光谱图像分割方法，其特征在于所述的反复调 整加权参数，得到多组联合相异度矩阵，是保持两个矩阵的加权参数合为1，设 定像素特征欧氏距离矩阵的加权参数初始值为0，然后将该参数每次增加0. 05， 直到该参数增加到0.5，像素空间欧氏距离矩阵的加权参数则相应由1减少到 0.5。
3、 根据权利要求1所述的高光谱图像分割方法，特征在于所述的统计像素点空间近邻变化情况，按如下步骤进行3a)根据像素空间欧氏距离矩阵取得每个像素点八近邻变化情况，将每个点的八近邻变化最大三个相加得到一个距离和值，对每组映射结果进行相同处理， 得到多个距离和值，形成该点的和值集合；3b)将每一个点的和值集合组成偏移行向量，并根据参数调整的顺序用后一 个和值减去前一个和值，得到近邻变化行向量，这两个行向量表示该点的近邻变 化情况，统计近邻变化行向量中大于O的数值个数，得到近邻变化值。全文摘要
本发明公开了一种基于像素空间信息的高光谱图像分割方法，主要是解决现有方法中不能将相近地貌很好的分割出来的不足。其分割步骤是求解高光谱数据的像素特征和像素空间欧氏距离矩阵，并进行归一化；对像素特征和像素空间欧氏距离矩阵进行加权，将加权后的两个矩阵相加形成联合相异度矩阵，调节加权参数，得到多组联合相异度矩阵；使用等距映射算法对每一组的联合相异度矩阵进行降维，得到多组映射结果；对每组映射结果进行统计分析，完成高光谱图像的初步分割；对初步划分的边界点进行类别校正，得到最后的图像分割结果。本发明能够有效地找到高光谱图像中的不同地貌的细微差别，可用于军事目标识别，矿藏勘探，环境状况分析。
文档编号G06T7/00GK101667291SQ20091002405
公开日2010年3月10日 申请日期2009年9月25日 优先权日2009年9月25日
发明者周伟达, 周宏杰, 莉 张, 哲 晏, 焦李成 申请人:西安电子科技大学