基于空?谱信息的高光谱图像半监督分类方法
【专利摘要】本发明公开了一种基于空?谱信息的高光谱图像半监督分类方法，该方法将高光谱图像中的光谱信息、空间信息联合作用于支持向量机分类器，采用自训练的半监督分类框架，并利用主动学习的方法作为半监督分类的样本选取策略，将半监督分类得到的初始分类结果进行按类分解得到的各类的二值图像作为边缘保留滤波器的输入图像，光谱信息的第一主成分分量作为滤波器的参照图像，利用边缘保留滤波器进行局部平滑，消除噪声，然后将像元按照所属概率最大的类别进行划分，完成分类过程。本发明联合光谱信息和空间信息提高了类别的可分性，利用自训练的半监督分类框架解决高光谱图像小样本分类问题，能有效地消除初始分类结果中的斑点状误差，提高分类精度。
【专利说明】
基于空-谱信息的高光谱图像半监督分类方法
技术领域
[0001 ]本发明属于信息技术领域，涉及模式识别、图像处理技术，具体是一种联合使用高 光谱图像光谱信息和空间信息的半监督分类方法。
【背景技术】
[0002 ]高光谱图像(光谱分辨率在101数量级范围内的光谱图像)地物分类是高光谱遥感 数据信息处理的重要内容，也是人们从遥感图像上获取有用信息的重要途径。高光谱图像 分类是基于图像像元的空间和光谱特性，对每个像元或比较匀质的像元组中所代表的不同 类别地物进行类别属性的确定和标注。
[0003] 近年来，随着航天航空技术和遥感传感器技术的不断发展，遥感对地观测系统采 集到图像的光谱分辨率也越来越高。这些丰富的光谱信息增加了对感兴趣物体实现更准确 识别的可能性。然而高光谱遥感图像光谱分辨率在提升的同时也使得数据维数以及数据量 不断增加，并且高光谱图像覆盖区域广，实地勘察困难且周期长，有标签样本的获取需要耗 费大量的人力物力，所以有标签样本的数量十分有限。
[0004] 采用监督分类方法对高光谱图像进行分类时，由于高光谱图像数据的高维度与有 标签样本数量之间的不平衡，会引起Hughes现象，严重影响了分类器的性能，不能得到很好 的分类效果。而半监督分类方法则是综合利用有标签样本和无标签样本中所包含的有用信 息，来改善分类器的性能，所以半监督分类方法在有标签样本数量少时也能取得较好的分 类效果。
[0005] 传统的分类方法通常只使用数据光谱维上的信息，却忽略了空间维所包含的信 息，而且单一的光谱信息对分类精度的提升效果是有限的。
【发明内容】
[0006] 针对上述现有技术中存在的问题，本发明的目的在于，提供一种基于空-谱信息的 高光谱图像半监督分类方法，有效地利用光谱图像的空间信息，以更好地解决高光谱图像 小样本分类问题。
[0007] 为了实现上述任务，本发明采用以下技术方案：
[0008] -种基于空-谱信息的高光谱图像半监督分类方法，包括以下步骤：
[0009] 步骤一，空-谱信息结合的高光谱图像半监督分类
[0010]利用主成分分析法和独立成分分析法提取高光谱图像的光谱信息，利用Gabor滤 波器组从光谱信息的第一主成分分量中提取高光谱图像的空间信息，将光谱信息和空间信 息联合作用于支持向量机分类器，采用自训练的半监督分类框架，并利用主动学习的方法 选取出包含信息最丰富的无标签样本扩充到初始训练样本中，以进行初始分类；
[0011] 步骤二，基于边缘保留滤波的分类结果空间优化
[0012] 将初始分类结果进行按类分解得到的各类的二值图像作为边缘保留滤波器的输 入图像以提供类别信息，光谱信息的第一主成分分量作为边缘保留滤波器的参照图像以提 供边缘信息，利用边缘保留滤波器在保留图像中原有边缘结构的同时进行局部平滑，从而 消除初始分类结果中的斑点状噪声；
[0013] 步骤三，对于待分类的某一像元，将其划分到所属概率最大的一类中，得到最终的 分类结果。
[0014] 进一步地，所述的步骤一中，Gabor滤波器组以及支持向量机分类器的参数设置过 程包括：
[0015] 步骤1-1，选取样本，并把所选样本随机分成测试样本集T和训练样本集，选取N个 类别进行分类，每类从训练样本集中选取s个样本，共n = s X N个有标签样本构成初始训练 样本集Dl = {(Xll，yil)，···，（Xln，yin)}，剩余的样本组成无标签样本集U;其中XII，· · .Xln表示 样本，yii，. . .，yin表示样本的标签；
[0016] 步骤1-2,设置Gabor滤波器的正弦平面波的最大频率为fmax，频率和尺度个数分别 为P、Q，高斯半径为〇,相位角为Φ ;
[0017]步骤1-3,利用M次交叉验证得到支持向量机分类器最优的高斯半径和惩罚参数。
[0018] 进一步地，所述的步骤一中进行光谱信息提取、空间信息提取并将光谱信息和空 间信息联合的具体过程包括：
[0019] 步骤1-4,光谱信息提取
[0020] 利用主成分分析法对高光谱图像进行降维，得到每个像元的光谱特征xwer，再对 得到的光谱特征利用独立成分分析法得到最终的光谱信息x keRk;其中R表示实数域，w、k表 示维数，Xk，X w分别是k，W维的实数向量；
[0021] 步骤1-5,空间信息提取
[0022] 利用Gabor滤波器组从光谱信息的第一主成分分量中提取空间信息，得到每个像 元的空间信息Xs e Rs; Xs、Rs分别指的是S维的实数向量和S维的实数空间;Gabor滤波器的核 函数为：
[0023]
[0024] 上式中，X' =xcos9+ysin9，y' =-xsin9+ycos9，χ和y是图像中像素坐标，f和Θ分别 是正弦平面波的频率和方向角，Φ是Gabor函数的相位角，σ和γ分别表示高斯半径和方向 角;利用具有不同方向和频率的Gabor滤波器来提取图像中不同的空间纹理信息，具体表示 为：
[0025]
[0026]上式中，Q和P表示频率和尺度个数，fmax表示正弦平面波的最大频率；
[0027]步骤1-6，光谱信息和空间信息联合
[0028] 将光谱信息Xk e Rk和空间信息Xk e Rk进行级联，得到(Xk，Xs)。
[0029] 进一步地，所述的步骤一中进行初始分类的具体过程包括：
[0030] 步骤1-7,利用训练样本集0!对支持向量机分类器进行训练，并用训练好的分类器 对无标签样本集U和测试样本集T中的样本进行分类，得到各样本的预测标签；用测试样本 集的分类精度作为支持向量机分类器的性能评价参数；
[0031] 步骤1-8,候选样本集的产生
[0032]将训练样本集D1中样本的二阶邻域样本点选取出来，去除冗余点，将它们与步骤 1-7中的预测标签和中心训练样本的标签比较，保留和中心训练样本的标签一致的点，得到 的候选样本集表示为Ud = {(xdi，ydi)，…，（xdh，ydh) }，其中，xdi, .. .Xdh表示候选样本，ydi, ...ydh表示样本的预测标签，h表示样本个数；
[0033]步骤1-9,简化候选样本集
[0034] 利用BT主动学习的方法简化候选样本集Ud，选出包含信息量最丰富的c个样本K， i表示第i次迭代;然后更新训练样本M ;
[0035] 步骤1-10,当迭代次数到达给定阈值S时停止，否则跳转至步骤1-7。
[0036]进一步地，所述的步骤1-9中BT主动学习方法的策略为：
[0037]
[0038] 式中，xB%选出的样本，U为无标签样本集，yi*为样本的标签，Xi表示样本点，ω表 示某一类别，N是类别总数，ω +表示样本^所属概率最大的那一类。
[0039] 进一步地，所述的步骤二的具体过程包括：
[0040] 步骤2-1，分解初始分类结果
[0041 ]把步骤一得到的初始分类结果c按类别分解成二值图像P= (P1，…，ΡΝ)，Ρ表示所有 类的二值图像的集合，Pm(m=l，-_，N)表示第m类的二值图像，用Pl,m表示P m中像素 i处的元 素，它的取值如下所示：
[0042]
[0043]上式中，m表示第m类地物的数字编号，i表示图像中的像素点，(^表示初始分类结 果c在像素点i的值；
[0044] 步骤2-2,边缘保留滤波
[0045] 把步骤2-1得到的二值图像作为边缘保留滤波器的输入，第一主成分分量作为边 缘保留滤波器的参照图像，对初始的分类结果进行局部的区域平滑操作，得到各类的滤波 结果，具体滤波过程表示如下：
[0046]
[0047]上式中，i表示图像中的像素点，j表示以像素点i为中心的局部窗口范围内的点，m 表示第m类地物的数字编号，表示第m类的二值图像中像素点j对应的值，I是边缘保留滤 波器的参照图像，久,"是第m类在像素点i处的滤波结果，W^(I)是权重，它的作用使滤波结 果中能够保留参照图像中显著的边缘信息，权重W^(I)由双边滤波器得到：
[0048]
[0049] 上式中，i，j表示像素点，Coi表示一个以i为中心大小为(2δ3+1)Χ(2δ3+1)的窗口， 瓦f是一个归一化因子，参数Sr表示当参照图像像素间差异增加时对应的权重减少的比例， 参数33表示窗口的大小
j 分别表示参考图像中像素点i和j对应的灰度值;经边缘保留滤波器后，得到各类的平滑结 果，其值是一个概率，表示像元属于该类的可能性。
[0050] 本发明与现有技术相比具有以下技术特点：
[0051] 1.本发明在对高光谱图像进行分类时，引入空间信息，可以获取更多对地物分类 有用的信息，从而对光谱信息的分类结果进行优化；
[0052] 2.本发明联合光谱信息和空间信息提高了类别的可分性，利用自训练的半监督分 类框架解决高光谱图像小样本分类问题，通过边缘保留滤波器的方式利用空间信息对初始 分类结果进行空间优化，能够有效地消除初始分类结果中的斑点状误差，提高分类精度。
【附图说明】
[0053]图1为本发明的总体框架图；
[0054]图2为自训练半监督分类方法流程图；
[0055] 图3(a)为地物分布Ground truth;
[0056] 图3 (b)为支持向量机分类器的分类结果；
[0057] 图3(c)为初始分类结果；
[0058]图3(d)为最终分类结果；
【具体实施方式】
[0059]遵从上述技术方案，如图1所示，本发明的具体过程如下：
[0060]步骤一，空-谱信息结合的高光谱图像半监督分类
[0061] 利用主成分分析法和独立成分分析法提取高光谱图像的光谱信息，利用Gabor滤 波器组从光谱信息的第一主成分分量中提取高光谱图像的空间信息，将光谱信息和空间信 息联合作用于支持向量机分类器，采用自训练的半监督分类框架，并利用主动学习的方法 选取出包含信息最丰富的无标签样本扩充到初始训练样本中，以进行初始分类;具体包括 以下步骤：
[0062] 本方案中首先需要选取训练样本集和测试样本集，以及确定Gabor滤波器和基于 概率模型的支持向量机(SVM)分类器的参数：
[0063] 步骤1-1，选取样本，并把所选样本随机分成测试样本集T和训练样本集，选取N个 类别进行分类，并用数字1到N来表示各个不同的类;每类从训练样本集中选取s个样本，共η =8父1'1个有标签样本构成初始训练样本集〇1={(111，711)，"_，（1111，7111)}，剩余的样本组成 无标签样本集U;其中?1，. . .Xln表示样本，yii，...，yin表示样本的标签；
[0064] 步骤1-2,设置Gabor滤波器的正弦平面波的最大频率为fmax，频率和尺度个数分别 为P、Q，高斯半径为〇,相位角为Φ ;
[0065]步骤1-3,利用M次交叉验证得到支持向量机分类器最优的高斯半径和惩罚参数。
[0066]样本集及初始参数设置好之后，进行光谱信息、空间信息的提取：
[0067]步骤1-4,光谱信息提取
[0068]利用主成分分析法对高光谱图像进行降维，得到每个像元的光谱特征xwer，再对 得到的光谱特征利用独立成分分析法得到最终的光谱信息xkeRk;其中R表示实数域，w、k表 示维数，x k，Xw分别是k，w维的实数向量；
[0069] 步骤1-5,空间信息提取
[0070] 利用Gabor滤波器组从光谱信息的第一主成分分量中提取空间信息，得到每个像 元的空间信息Xs e Rs; Xs、Rs分别指的是S维的实数向量和S维的实数空间;Gabor滤波器的核 函数为：
[0071]
[0072] 上式中，X' =xcos9+ysin9，y' =-xsin9+ycos9，χ和y是图像中像素坐标，f和Θ分别 是正弦平面波的频率和方向角，Φ是Gabor函数的相位角，σ和γ分别表示高斯半径和方向 角;利用具有不同方向和频率的Gabor滤波器来提取图像中不同的空间纹理信息，具体表示 为：
[0073]
[0074] 上式中，Q和P表示频率和尺度个数，fmax表示正弦平面波的最大频率；
[0075] 光谱信息以及空间信息提取后，将二者联合，以作用于支持向量机分类器：
[0076] 步骤1-6，光谱信息和空间信息联合
[0077] 将光谱信息Xk e Rk和空间信息xk e Rk进行级联，得到(xk，xs)。
[0078]进行初始分类的具体过程包括：
[0079] 步骤1-7,利用训练样本集0!对支持向量机分类器进行训练，并用训练好的分类器 对无标签样本集U和测试样本集T中的样本进行分类，得到各样本的预测标签；用测试样本 集的分类精度作为支持向量机分类器的性能评价参数；
[0080] 步骤1-8,候选样本集的产生
[0081] 将训练样本集D1中样本的二阶邻域样本点选取出来，去除冗余点（如背景点、测试 样本等），将它们与步骤1-7中的预测标签和中心训练样本的标签比较(预测标签指的用训 练好的分类器分类得到的标签，中心训练样本的标签是训练样本的标签，这个标签是已知 的，保留和中心训练样本的标签一致的点，得到的候选样本集表示为Ud= {(Xdl，ydl)，···， (xdh，ydh)}，其中，xdi，· · .Xdh表示候选样本，ydi，· · .ydh表示样本的预测标签，h表示样本个 数；
[0082]步骤1-9,简化候选样本集
[0083] 利用BT主动学习的方法简化候选样本集Ud，选出包含信息量最丰富的c个样本.， i表示第i次迭代;然后更新训练样本A = ·(舄，Μ}
[0084] 步骤1-9中BT主动学习方法的策略为：
[0085]
[0086] 式中，χΒΙ为选出的样本，U为无标签样本集，yi*为样本的标签，Xl表示样本点，ω表 示某一类别，N是类别总数，ω +表示样本^所属概率最大的那一类。
[0087] 步骤1-1至步骤1-9的重复迭代过程构成了采用自训练的半监督分类框架。
[0088] 步骤1-10,对于步骤1-9,当迭代次数到达给定阈值S时停止，否则跳转至步骤1-7。
[0089]步骤二，基于边缘保留滤波(EPF)的分类结果空间优化
[0090]将初始分类结果进行按类分解得到的各类的二值图像作为边缘保留滤波器的输 入图像以提供类别信息，光谱信息的第一主成分分量作为Gabor滤波器的参照图像以提供 边缘信息，利用边缘保留滤波器在保留图像中原有边缘结构的同时进行局部平滑，平滑的 到的值是一个概率，表示像素属于该类的可能性;从而消除初始分类结果中的斑点状噪声 (错误分类）；具体过程如下：
[0091 ]步骤2-1，分解初始分类结果
[0092]把步骤一得到的初始分类结果c按类别分解成二值图像P= (P1，…，Pn)，P表示所有 类的二值图像的集合，Pm(m=l，-_，N)表示第m类的二值图像，用Pl,m表示P m中像素 i处的元 素，它的取值如下所示：
[0093]
[0094] 上式中，m表示第m类地物的数字编号，i表示图像中的像素点，(^表示初始分类结 果c在像素点i的值；
[0095]步骤2-2,边缘保留滤波
[0096] 把步骤2-1得到的二值图像作为边缘保留滤波器的输入，第一主成分分量作为边 缘保留滤波器的参照图像，对初始的分类结果进行局部的区域平滑操作，得到各类的滤波 结果，具体滤波过程表示如下：
[0097]
[0098]上式中，i表示图像中的像素点，j表示以像素点i为中心的局部窗口范围内的点，m 表示第m类地物的数字编号，表示第m类的二值图像中像素点j对应的值（即步骤2-1中的 P1,m)，I是边缘保留滤波器的参照图像，氣m是第m类在像素点i处的滤波结果，W^(I)是权 重，它的作用使滤波结果中能够保留参照图像中显著的边缘信息，权重W^(I)由双边滤波 器得到：
[0099]
[0100] 上式中，i，j表示像素点，Wi表示一个以i为中心大小为(2δ3+1)Χ(2δ3+1)的窗口， K6是一个归一化因子，参数Sr表示当参照图像像素间差异增加时对应的权重减少的比例， 参数33表示窗口的大小
分别表示参考图像中像素点i和j对应的灰度值;经边缘保留滤波器滤波后，得到各类的平 滑结果，其值是一个概率，表示像元属于该类的可能性。
[0101]步骤三，取最大值
[0102] 对于某一像元，将其划分到所属概率最大的那一类中，得到最终的分类结果。
[0103] 以下是发明人给出的一个具体实施例：
[0104] 实验数据是1992年由AVIRIS在美国印第安纳州西北部的Indian Pines实验基地 拍摄得到。该数据包含224个波段的145X 145大小的图像，空间分辨率为20m，波长范围为 0.4-2.5μπι(水吸收波段104-108，150-163，220在实验前已去除，实际使用200个波段）。它包 含农田、森林以及其他多种原生植被。
[0105] 在本实施例中，方法设置的具体参数为：
[0106] 随机选取总样本的25 %作为测试样本集T，75 %作为训练样本集，选取N=9个类别 进行分类，用数字1到9来表示各个不同的类，每类从训练样本集中选取S = 15个样本，共η = 135个有标签样本构成初始训练样本集Di = {(Xii，yn)，…，（xii35，yii35)}，剩余的样本为无 标签样本集U;设置Gabor滤波器的正弦平面波的最大频率为f max = 0.25,频率和尺度个数分 别为P = 6、Q = 3,高斯半径为σ = 1，相位角为φ =〇;利用M = 5次交叉验证得到支持向量机分 类器最优的高斯半径和惩罚参数;双边滤波器参数心= 〇.2，δ3 = 3;利用PCA变换的第一主成 分作为边缘保留滤波器的参照图像，上一步的分类结果作为边缘保留滤波器的输入图像。
[0107] 用主成分分析法对高光谱图像进行降维，得到每个像元的光谱特征xw e Rw，W = 20， 再对得到的光谱特征利用独立成分分析法得到最终的光谱信息xkeRk，k=14。
[0108] 按照上述的参数设置，通过本发明方法得到的分类精度详情如表1所示，具体处理 过程中不同步骤的处理结果如图3所示。在图3中，图3(a)为地物分布Ground truth，图3(b) 为支持向量机分类器的分类结果，图3(c)为初始分类结果，图3(d)为最终分类结果;在表1 中，S 3C表示本发明步骤一的空-谱信息结合的半监督分类方法得到的初始分类结果，S3C_ EPF表示经步骤二的边缘保留滤波优化后的最终分类结果：
[0109] 表1不同方法得到的分类精度
[0112]实验结果中可以看出，通常监督分类方法对小样本分类问题的分类效果不是很 好，如Indian Pines数据的SVM分类精度都不够高。在Indian Pines数据的实验结果中，采 用半监督分类方法得到的结果与直接使用SVM分类的结果相比，分类精度得到了一定程度 上的提高。其中S 3C_BT方法得到的提升最大，Indian Pines数据提升了 13%，由此可以看出 本专利提出的分类方法的有效性。
【主权项】
1. 一种基于空-谱信息的高光谱图像半监督分类方法，其特征在于，包括以下步骤： 步骤一，空-谱信息结合的高光谱图像半监督分类 利用主成分分析法和独立成分分析法提取高光谱图像的光谱信息，利用Gabor滤波器 组从光谱信息的第一主成分分量中提取高光谱图像的空间信息，将光谱信息和空间信息联 合作用于支持向量机分类器，采用自训练的半监督分类框架，并利用主动学习的方法选取 出包含信息最丰富的无标签样本扩充到初始训练样本中，以进行初始分类； 步骤二，基于边缘保留滤波的分类结果空间优化 将初始分类结果进行按类分解得到的各类的二值图像作为边缘保留滤波器的输入图 像以提供类别信息，光谱信息的第一主成分分量作为边缘保留滤波器的参照图像以提供边 缘信息，利用边缘保留滤波器在保留图像中原有边缘结构的同时进行局部平滑，从而消除 初始分类结果中的斑点状噪声； 步骤三，对于待分类的某一像元，将其划分到所属概率最大的一类中，得到最终的分类 结果。2. 如权利要求1所述的基于空-谱信息的高光谱图像半监督分类方法，所述的步骤一 中，Gabor滤波器组以及支持向量机分类器的参数设置过程包括： 步骤1-1，选取样本，并把所选样本随机分成测试样本集T和训练样本集，选取N个类别 进行分类，每类从训练样本集中选取s个样本，共n = sXN个有标签样本构成初始训练样本 集Dl = {(Xll，yil)，…，（Xln，yin)}，剩余的样本组成无标签样本集U;其中XII，. . .Xln表示样 本，yn，...，yin表示样本的标签； 步骤1-2,设置Gabor滤波器的正弦平面波的最大频率为fmax，频率和尺度个数分别为P、 Q，高斯半径为〇,相位角为Φ ; 步骤1-3,利用Μ次交叉验证得到支持向量机分类器最优的高斯半径和惩罚参数。3. 如权利要求1所述的基于空-谱信息的高光谱图像半监督分类方法，其特征在于，所 述的步骤一中进行光谱信息提取、空间信息提取并将光谱信息和空间信息联合的具体过程 包括： 步骤1 _4，光谱fg息提取 利用主成分分析法对高光谱图像进行降维，得到每个像元的光谱特征xwer，再对得到 的光谱特征利用独立成分分析法得到最终的光谱信息xkeRk;其中R表示实数域，w、k表示维 数，x k，xw分别是k，w维的实数向量； 步骤1_5，空间彳目息提取 利用Gabor滤波器组从光谱信息的第一主成分分量中提取空间信息，得到每个像元的 空间信分别指的是s维的实数向量和s维的实数空间;Gabor滤波器的核函数 为：上式中，X' = xcos9+ysin9，y' = -xsin9+ycos9，χ和y是图像中像素坐标，f和Θ分别是正 弦平面波的频率和方向角，Φ是Gabor函数的相位角，σ和γ分别表示高斯半径和方向角;利 用具有不同方向和频率的Gabor滤波器来提取图像中不同的空间纹理信息，具体表示为： 上式中，Q和P表示频率和尺度个数，fmax表示正弦平面波的最大频率；步骤1-6，光谱信息和空间信息联合 将光谱信息xk e Rk和空间信息xk e Rk进行级联，得到(xk，xs)。4. 如权利要求1所述的基于空-谱信息的高光谱图像半监督分类方法，其特征在于，所 述的步骤一中进行初始分类的具体过程包括： 步骤1-7,利用训练样本集Ddi支持向量机分类器进行训练，并用训练好的分类器对无 标签样本集U和测试样本集T中的样本进行分类，得到各样本的预测标签；用测试样本集的 分类精度作为支持向量机分类器的性能评价参数； 步骤1-8,候选样本集的产生 将训练样本集D冲样本的二阶邻域样本点选取出来，去除冗余点，将它们与步骤1-7中 的预测标签和中心训练样本的标签比较，保留和中心训练样本的标签一致的点，得到的候 选样本集表示为Ud={(xdi，ydi)，…，（xdh，ydh)}，其中，xdi,.. .xdh表示候选样本，ydi,.. .ydh表 示样本的预测标签，h表示样本个数； 步骤1-9，简化候选样本集 利用BT主动学习的方法简化候选样本集Ud，选出包含信息量最丰富的c个样本t/!，i表 示第i次迭代;然后更新训练样本私=[ A， 步骤1-10,当迭代次数到达给定阈值S时停止，否则跳转至步骤1-7。5. 如权利要求4所述的基于空-谱信息的高光谱图像半监督分类方法，其特征在于，所 述的步骤1 -9中BT主动学习方法的策略为：式中，xB%选出的样本，U为无标签样本集，yi*为样本的标签，xi表示样本点，ω表示某 一类别，Ν是类别总数，ω+表示样本^所属概率最大的那一类。6. 如权利要求1所述的基于空-谱信息的高光谱图像半监督分类方法，其特征在于，所 述的步骤二的具体过程包括： 步骤2-1，分解初始分类结果 把步骤一得到的初始分类结果c按类别分解成二值图像P=(Pi，…，Ρν)，Ρ表示所有类的 二值图像的集合，Pm(m=l，-_，N)表示第m类的二值图像，用pi,m表示Pm中像素 i处的元素，它 的取值如下所示：上式中，m表示第m类地物的数字编号，i表示图像中的像素点，Cl表示初始分类结果c在 像素点i的值； 步骤2-2,边缘保留滤波 把步骤2-1得到的二值图像作为边缘保留滤波器的输入，第一主成分分量作为边缘保 留滤波器的参照图像，对初始的分类结果进行局部的区域平滑操作，得到各类的滤波结果， 具体滤波过程表示如下：上式中，i表示图像中的像素点，j表示以像素点i为中心的局部窗口范围内的点，m表示 第m类地物的数字编号，Pj,m表示第m类的二值图像中像素点j对应的值，I是边缘保留滤波器 的参照图像，氣 m是第m类在像素点i处的滤波结果，Wu(I)是权重，它的作用使滤波结果中 能够保留参照图像中显著的边缘信息，权重Wu(I)由双边滤波器得到：上式中，i，j表示像素点，ω,表示一个以i为中心大小为(2δ3+1)Χ(2δ3+1)的窗口，<是 一个归一化因子，参数心表示当参照图像像素间差异增加时对应的权重减少的比例，参数心 表不窗口的大小，，.Ii和Ij分别表 示参考图像中像素点i和j对应的灰度值;经边缘保留滤波器后，得到各类的平滑结果，其值 是一个概率，表示像元属于该类的可能性。
【文档编号】G06K9/62GK106056157SQ201610382887
【公开日】2016年10月26日
【申请日】2016年6月1日
【发明人】彭进业, 赵二龙, 刘胜杰, 王珺, 章勇勤, 李展, 祝轩, 周剑虹, 艾娜
【申请人】西北大学