发明名称 基于局部和非局部多特征语义高光谱图像分类方法
摘要 本发明公开了一种基于局部和非局部多特征语义高光谱图像分类方法。主要解决高光谱图像分类中正确率低,鲁棒性差,空间一致性弱的问题。其步骤包括:输入图像,提取图像多种特征;数据集切分为训练集和测试集;概率支持矢量机将所有样本的多种特征映射成相应的语义表示;构造局部以及非局部近邻集合;构建降噪马尔可夫场模型,进行语义融合和降噪处理;对语义表示迭代优化;利用语义表示求得所有样本的类别,完成高光谱图像准确分类。本发明采用了多特征融合,并对存在于图像中的空间信息充分挖掘和利用,在小样本情况下,得到了非常高的分类精度,并拥有良好的鲁棒性和空间一致性,用于军事探测、地图绘制、植被调查、矿物检测等方面。
申请公布号 CN106529508A 申请公布日期 2017.03.22
申请号 CN201611119573.1 申请日期 2016.12.07
申请人 西安电子科技大学 发明人 张向荣;焦李成;高泽宇;冯婕;白静;侯彪;马文萍;李阳阳
分类号 G06K9/00(2006.01)I;G06K9/46(2006.01)I;G06K9/62(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 程晓霞;王品华
主权项 一种基于局部和非局部多特征语义高光谱图像分类方法,其特征在于,包括有如下步骤:(1)输入图像,提取图像的多种特征:输入高光谱图像<img file="FDA0001173486970000011.GIF" wi="259" he="95" />M是高光谱图像中所有像素点的总个数,h<sub>q</sub>为一个列向量,代表像素点q在每个波段的反射值所构成的向量;对高光谱图像分别提取多种特征,包括有:原始光谱特征、Gabor纹理特征、差分形态学特征;该高光谱图像包含c类像素点,其中有N个有标记像素点,m个无标记像素点,图像的每个像素点为一个样本,每个样本由V个特征向量构成,分别代表该样本在不同特征描述子下的表述,V是特征类别的个数;(2)高光谱图像数据集分为训练集和测试集:用N个有标记像素点作为训练样本构成训练集<img file="FDA0001173486970000012.GIF" wi="691" he="95" />其对应的类别标记集为<img file="FDA0001173486970000013.GIF" wi="211" he="78" />用m个无标记像素点作为测试样本构成测试集<img file="FDA0001173486970000014.GIF" wi="723" he="103" />其中,x<sub>i</sub>表示训练集的第i个样本,y<sub>j</sub>表示测试集的第j个样本,l<sub>i</sub>是第i个训练样本所属的类别标号,D<sub>v</sub>表示第v类特征的维数,R表示实数域;(3)利用概率支持矢量机(SVM)将所有样本的多种特征映射成相应的语义表示:分别利用训练集<img file="FDA0001173486970000015.GIF" wi="658" he="94" />中所有样本的V个特征向量以及其对应的类别标记集<img file="FDA0001173486970000016.GIF" wi="207" he="76" />构建V个概率支撑矢量机分类器,该分类器的核函数为径向高斯核,核参数以及惩罚参数由多倍交叉验证得到;将测试集<img file="FDA0001173486970000017.GIF" wi="693" he="95" />中所有样本的V个特征向量分别输入到构建的V个对应分类器中,得到在不同特征描述子表述下,每个测试样本y<sub>j</sub>,j=1,2,…,m属于每个类别的概率,作为每个测试样本的语义表示<img file="FDA0001173486970000018.GIF" wi="531" he="94" />对于训练集中的每个样本x<sub>i</sub>,i=1,2,…,N,它属于本身类别l<sub>i</sub>的概率为1,而属于其他类别的概率为0,得到多种特征对应的多种语义表示<img file="FDA0001173486970000019.GIF" wi="521" he="95" />其中<img file="FDA00011734869700000110.GIF" wi="276" he="55" />的第l<sub>i</sub>行为1,其它行为0;从而得到该高光谱图像中所有样本的多种语义表示<img file="FDA0001173486970000021.GIF" wi="226" he="87" />(4)构造测试集中所有样本的局部以及非局部近邻集合:对于测试集中的每个样本y<sub>j</sub>,j=1,2,…,m,构造其局部自适应近邻集合B<sub>j</sub>和非局部相似结构近邻集合C<sub>j</sub>;(5)构建降噪马尔可夫场模型,进行测试样本的多种语义表示融合以及语义表示的降噪处理:对每个测试样本y<sub>j</sub>,j=1,2,…,m分别进行如下操作,将y<sub>j</sub>所对应的语义表示<img file="FDA0001173486970000022.GIF" wi="235" he="93" />局部自适应近邻集合B<sub>j</sub>中所有样本的语义表示<img file="FDA0001173486970000023.GIF" wi="365" he="93" />和非局部相似结构近邻集合C<sub>j</sub>中所有样本的语义表示<img file="FDA0001173486970000024.GIF" wi="403" he="95" />均输入到局部能量函数中,最小化该能量函数,得到测试样本y<sub>j</sub>的一阶降噪语义表示<img file="FDA0001173486970000025.GIF" wi="98" he="86" />与此同时保持训练集中样本的语义表示不变,得到该高光谱图像所有样本的一阶降噪语义表示<img file="FDA0001173486970000026.GIF" wi="99" he="71" />(6)对高光谱图像所有样本的一阶降噪语义表示<img file="FDA0001173486970000027.GIF" wi="74" he="62" />进一步迭代优化:设定最大迭代次数T<sub>max</sub>,t为当前迭代代数,对每个测试样本进行如下操作:将测试样本y<sub>j</sub>,以及集合B<sub>j</sub>和集合C<sub>j</sub>中所有样本的第t阶降噪语义表示<img file="FDA0001173486970000028.GIF" wi="594" he="86" />作为降噪马尔可夫场局部能量函数的输入,最小化该能量函数,得到测试样本y<sub>j</sub>的第(t+1)阶语义表示<img file="FDA0001173486970000029.GIF" wi="115" he="80" />与此同时继续保持训练集中样本的语义表示不变,进而得到该高光谱图像所有样本的第(t+1)阶降噪语义表示<img file="FDA00011734869700000210.GIF" wi="126" he="64" />重复迭代过程,直到t=T<sub>max</sub>‑1停止,得到该高光谱图像所有样本的第T<sub>max</sub>阶降噪语义表示,也就是最终的语义表示<img file="FDA00011734869700000211.GIF" wi="146" he="71" />(7)利用最终的语义表示<img file="FDA00011734869700000212.GIF" wi="147" he="63" />求得测试集中所有样本的类别:对于测试集中的每个样本y<sub>j</sub>,j=1,2,…,m,其最终的语义表示为<img file="FDA00011734869700000213.GIF" wi="134" he="76" />即测试样本y<sub>j</sub>属于每个类别的概率所组成的列向量,选择该向量中最大值元素所在位置的标号作为y<sub>j</sub>的类别<img file="FDA00011734869700000214.GIF" wi="53" he="69" />从而得到测试集类别预测集合<img file="FDA00011734869700000215.GIF" wi="211" he="94" />完成该高光谱图像的分类任务。
地址 710071 陕西省西安市雁塔区太白南路2号