CN109784392B - 一种基于综合置信的高光谱图像半监督分类方法 - Google Patents

一种基于综合置信的高光谱图像半监督分类方法 Download PDF

Info

Publication number
CN109784392B
CN109784392B CN201910012238.9A CN201910012238A CN109784392B CN 109784392 B CN109784392 B CN 109784392B CN 201910012238 A CN201910012238 A CN 201910012238A CN 109784392 B CN109784392 B CN 109784392B
Authority
CN
China
Prior art keywords
sample
training
semi
label
supervised
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910012238.9A
Other languages
English (en)
Other versions
CN109784392A (zh
Inventor
贺霖
余龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN201910012238.9A priority Critical patent/CN109784392B/zh
Publication of CN109784392A publication Critical patent/CN109784392A/zh
Application granted granted Critical
Publication of CN109784392B publication Critical patent/CN109784392B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Abstract

本发明公开了一种基于综合置信的高光谱图像半监督分类方法,包括:读入高光谱图像;计算图权重矩阵;8近邻连接稀疏化图权重矩阵;计算归一化图权重矩阵;获取初始训练集和候选集;设定协同训练迭代次数并开始训练过程;训练多项式逻辑回归分类器;使用多项式逻辑回归分类器得到候选集样本的预测标签;使用半监督图分类方法得到候选集样本的预测标签;选出两种预测标签一致的候选样本及对应的预测标签构成协定集,对应的置信系数构成综合置信集;筛选出综合置信系数高于99%的协定集样本及对应的预测标签,构成扩增集加入训练集;移除候选集中扩增集样本;判断训练是否达到设定次数,若未则继续迭代,若达到则停止,用半监督图分类高光谱图像。

Description

一种基于综合置信的高光谱图像半监督分类方法
技术领域
本发明涉及高维图像处理技术领域,具体涉及一种基于综合置信的高光谱图像半监督分类方法。
背景技术
高光谱图像是由高光谱传感器获取的、光谱分辨率达到纳米级的地物遥感图像。高光谱图像每个像素的光谱数据来自于不同波长的光在地物目标所处该像素位置上的反射率,光谱特征可以用于识别该像素所属的地物类别。高光谱图像有以下特点:大量的像素点、高维的光谱特征。同时,高光谱图像的“图谱合一”特性也反映了图像中空间信息较为丰富,光谱特征具有区分性。利用这些先验信息发展起来的高光谱图像分类技术相比全色图像、多光谱图像具有更加突出的地物识别能力,因此在城市规划、海洋探测、军事、天文等领域均受到广泛的关注。
另外,高光谱图像分类技术仍存在标记样本代价昂贵等问题,由于高光谱数据的光谱特征维度高、标记训练样本数量少,许多分类方法难以学习到有效的模型,也无法达到令人满意的分类精度。因此,针对数据量大、样本特征维度高、标记样本少的高光谱图像,设计合适的半监督学习方法从大量无标记样本中学习出更多的可分性信息,从而实现高光谱图像分类精度的提升,是高光谱图像分类领域较为重要的研究方向之一。
许多有监督分类方法如多项式逻辑回归虽然能够在不少分类问题中有不错的效果,但在半监督分类问题中,由于训练样本较少,对噪声较大的高光谱图像只能得到粗糙的分类结果;基于图的半监督分类方法虽然能够利用高光谱图像中的空间平滑性信息实现较为不错的分类结果,但该类方法的类边界区分能力十分弱。因此,已有的众多有监督、半监督学习方法在训练样本较少的情况下,对高光谱图像的分类能力有待提高。
发明内容
本发明的目的是针对现有技术的不足,提供了一种基于综合置信的高光谱图像半监督分类方法。首先利用高光谱数据建立了空间最近邻的图连接结构,然后结合空间信息和光谱信息以及多分类器综合置信筛选机制实现了训练样本的扩增,据此设计了一种基于综合置信和训练信息迭代扩增的协同训练自学习方法,最后利用半监督图分类方法对高光谱图像分类。所述基于综合置信的自学习方法十分适用于空间信息较多的高光谱图像数据,可以在无人为干预的情况下自学习出较为可靠的训练样本,极大地提升高光谱图像分类精度。
本发明的目的可以通过如下技术方案实现:
一种基于综合置信的高光谱图像半监督分类方法,所述方法包括以下步骤:
S1、读入三维高光谱图像立方体H(m,n,b),其中m、n表示空间像素点位置,b表示所处光谱波段位置;
S2、计算高光谱数据像素点之间的去样本均值的相关系数,用于构造图权重矩阵W,权重数值可以衡量像素点间的相似度;
S3、将步骤S2的图权重矩阵W中每个像素点与自身以及与空间8近邻外的像素点的相似度权重值置0,即只为每个像素点与其空间8近邻内的像素点间建立图连接,得到稀疏化的图权重矩阵
Figure GDA0002727099940000021
S4、计算归一化图权重矩阵
Figure GDA0002727099940000022
S5、获取初始训练集
Figure GDA0002727099940000023
和候选集
Figure GDA0002727099940000024
S6、设定协同训练的迭代次数M,并开始协同训练过程;
S7、使用训练集
Figure GDA0002727099940000025
训练多项式逻辑回归分类器,ξ∈[0,1,...,M-1]表示当前过程处于协同训练的第ξ+1次迭代;
S8、使用多项式逻辑回归分类器对候选集
Figure GDA0002727099940000026
中的样本分类,得到预测标签Y1和对应的置信度P;
S9、使用半监督图分类方法对候选集
Figure GDA0002727099940000027
中的样本分类,得到预测标签Y2
S10、取出候选集
Figure GDA0002727099940000028
中预测标签Y1和预测标签Y2相同的样本,以及对应的预测标签YA,构成协定集DA,从候选集置信度P中选出协定集样本对应的所有置信系数作为综合置信集PA
S11、筛选出协定集DA对应综合置信集PA中置信系数高于99%的样本,赋予预测标签YN,构成扩增集DN,加入训练集
Figure GDA0002727099940000029
中,得到更新的训练集
Figure GDA00027270999400000210
S12、将S11中扩增集DN的样本从候选集
Figure GDA00027270999400000211
中移除,得到更新的候选集
Figure GDA00027270999400000212
S13、迭代循环S7-S12步骤的高光谱图像协同训练过程,直至协同训练迭代达到设定值M次;
S14、使用半监督图分类方法和更新后的训练集
Figure GDA00027270999400000213
分类高光谱图像验证分类精度。
进一步地,步骤S2中的图权重矩阵W采用如下计算方式:
Figure GDA0002727099940000031
其中,vi代表第i个像素点的数据特征,va代表所有像素点的数据特征均值,wij是图权重矩阵W中对应第i行、第j列的元素,同时代表第i个像素点和第j个像素点的相似性权重值。
进一步地,步骤S4中的归一化图权重矩阵的计算采用
Figure GDA0002727099940000032
其中度矩阵D计算公式如下:
Figure GDA0002727099940000033
其中,wij是图权重矩阵W中第i行、第j列的元素,V是像素点总数。
进一步地,步骤S8中所述的初始候选集
Figure GDA0002727099940000034
由整个测试集内的样本构成,通过多项式逻辑回归分类器对整个候选集样本分类,得到第一种预测标签:
Figure GDA0002727099940000035
其中,γ代表候选集样本数,
Figure GDA0002727099940000036
代表多项式逻辑回归分类器对候选集
Figure GDA0002727099940000037
中第κ个样本的预测结果,κ∈[1,2,...,γ];计算多项式逻辑回归分类器对每个像素点的最大分类后验概率值,作为预测置信度P=[p1,p2,...,pγ],其中pκ代表候选集
Figure GDA0002727099940000038
中第κ个样本的预测置信度。
进一步地,步骤S9中的半监督图分类方法采用如下计算方式:
Figure GDA0002727099940000039
其中,α为训练参数,Y(0)由训练集
Figure GDA00027270999400000310
的样本标签
Figure GDA00027270999400000311
构成:
Figure GDA00027270999400000312
其中,i指代样本序号,c代表样本标签;
Figure GDA00027270999400000313
表示第i个样本的样本标签;[Γ12,...,Γδ]为训练集
Figure GDA00027270999400000314
中样本的样本序号;
使用半监督图分类方法对整个候选集
Figure GDA00027270999400000315
中的样本分类,得到第二种预测标签:
Figure GDA0002727099940000041
其中,γ代表候选集样本数,
Figure GDA0002727099940000042
代表半监督图分类方法对候选集
Figure GDA0002727099940000043
中第κ个样本的预测结果,κ表示候选集
Figure GDA0002727099940000044
中样本的序号,κ∈[1,2,...,γ]。
进一步地,步骤S10中的协定集DA由候选集
Figure GDA0002727099940000045
中预测标签Y1和预测标签Y2相同的样本构成:
Figure GDA0002727099940000046
其中,uκ表示
Figure GDA0002727099940000047
中的一个样本;
Figure GDA0002727099940000048
是预测标签Y1中的第κ个样本的预测标签,即多项式逻辑回归分类器对候选集
Figure GDA0002727099940000049
中第κ个样本的预测结果;
Figure GDA00027270999400000410
是预测标签Y2中的第κ个样本的预测标签,即半监督图分类方法对候选集
Figure GDA00027270999400000411
中第κ个样本的预测结果,κ∈[1,2,...,γ];
同时取出协定集样本对应的预测标签
Figure GDA00027270999400000412
以及协定集样本对应的所有置信系数作为综合置信集
Figure GDA00027270999400000413
pκ代表候选集
Figure GDA00027270999400000414
中第κ个样本的预测置信度;由于PA的确定是以两种不同的预测标签Y1和Y2为基础的,因此将其称为综合置信集。
进一步地,步骤S11中的扩增集DN由协定集DA中综合置信集PA>0.99的样本构成,扩增集DN的样本标签YN同样由YA中相应预测标签构成:
YN={YA(idx)|PA(idx)>0.99}
其中,YA(idx)、PA(idx)分别表示YA、PA中序号为idx的元素,而更新的训练集由之前的训练集和扩增集构成:
Figure GDA00027270999400000415
本发明与现有技术相比,具有如下优点和有益效果:
1、本发明提供的基于综合置信的高光谱图像半监督分类方法,使用了基于空间邻域相似性假设的图分类方法和基于光谱特征学习的多项式逻辑回归分类器,同时利用了空间信息和光谱信息进行自学习,从而增强分类器的分类性能。相比传统半监督图分类方法和多项式逻辑回归分类器,本发明方法结合了两种分类方法各自的优点,从而达到了增强的分类效果。
2、本发明提供的基于综合置信的高光谱图像半监督分类方法,在初始训练样本较少的情况下,通过选择多分类器决策一致的样本,自学习出更多的训练样本,增进了分类器的训练效果。再进一步筛选综合置信集中高置信系数的预测样本,使自学习的训练样本更加可靠。因此,在相同初始训练样本的条件下,相比传统方法,分类器拟合能力更强,分类精度更高,在半监督分类任务中更加有效。
附图说明
图1为本发明基于综合置信的高光谱图像半监督分类方法的工作流程图。
图2是对于AVIRIS Indian Pines图像,采用本发明的基于综合置信的高光谱图像半监督分类方法的分类精度随协同训练迭代次数增长的示意图;其中初始训练样本为对图像12类样本每类随机取6个的情况。
图3是对于AVIRIS Indian Pines图像,采用本发明的基于综合置信的高光谱图像半监督分类方法的分类精度在初始训练样本分别取每类2个、4个、6个、8个不同情况下的变化示意图;其中协同训练共迭代10次。
图4(a)是对于AVIRIS Indian Pines图像,采用半监督图分类方法的分类精度在初始训练样本分别取每类2个、4个、6个、8个不同情况下的变化示意图;图4(b)是对于AVIRIS Indian Pines图像,采用多项式逻辑回归分类器的分类精度在初始训练样本分别取每类2个、4个、6个、8个不同情况下的变化示意图;其中协同训练均迭代10次。
图5(a)-图5(d)是训练样本标签,图像中的像素点被分为12类,图中训练样本是对AVIRIS Indian Pines图像,通过在每类分别采样2个、4个、6个、8个像素点得到的。
图6(a)-图6(d)是测试样本标签,图像中的像素点被分为12类,四个图中抽取的训练样本是对AVIRIS Indian Pines图像,通过在每类分别采样2个、4个、6个、8个像素点得到的。
图7(a)-图7(d)是协同训练扩增后的训练样本标签,采用本发明提出的基于综合置信的自学习方法在图5(a)-图5(d)所示初始训练样本基础上对AVIRIS Indian Pines图像自学习得到的扩增训练样本,四个图中的初始训练样本是通过在每类分别采样2个、4个、6个、8个像素点得到的。
图8(a)-图8(d)是全体样本标签,采用本发明提出的基于综合置信的高光谱图像半监督分类方法得到对AVIRIS Indian Pines图像的分类结果,图像中的像素点被分为12类,四个图中抽取的训练样本是通过在每类分别采样2个、4个、6个、8个像素点得到的。
图9(a)-图9(d)是全体样本标签,采用半监督图分类方法得到对AVIRIS IndianPines图像的分类结果,图像中的像素点被分为12类,四个图中抽取的训练样本是通过在每类分别采样2个、4个、6个、8个像素点得到的。
图10(a)-图10(d)是全体样本标签,采用多项式逻辑回归分类器得到对AVIRISIndian Pines图像的分类结果,图像中的像素点被分为12类,四个图中抽取的训练样本是通过在每类分别采样2个、4个、6个、8个像素点得到的。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例:
本实施例提供了一种基于综合置信的高光谱图像半监督分类方法,所述方法的流程如图1所示,包括以下步骤:
S1、读入三维高光谱图像立方体H(m,n,b),其中m、n表示空间像素点位置,b表示所处光谱波段位置;
S2、计算高光谱数据像素点之间的去样本均值的相关系数,用于构造图权重矩阵W,权重数值可以衡量像素点间的相似度,计算方式如下:
Figure GDA0002727099940000061
其中,vi代表第i个像素点的数据特征,va代表所有像素点的数据特征均值,wij是图权重矩阵W中对应第i行、第j列的元素,同时代表第i个像素点和第j个像素点的相似性权重值;
S3、将步骤S2的图权重矩阵W中每个像素点与自身以及与空间8近邻外的像素点的相似度权重值置0,即只为每个像素点与其空间8近邻内的像素点间建立图连接,得到稀疏化的图权重矩阵
Figure GDA0002727099940000062
S4、计算归一化图权重矩阵
Figure GDA0002727099940000063
其中度矩阵D计算公式如下:
Figure GDA0002727099940000064
其中,wij是图权重矩阵W中第i行、第j列的元素,V是像素点总数;
S5、获取初始训练集
Figure GDA0002727099940000065
和候选集
Figure GDA0002727099940000066
初始候选集
Figure GDA0002727099940000067
由整个测试集内的样本构成;
S6、设定协同训练的迭代次数M,并开始协同训练过程;
S7、使用训练集
Figure GDA0002727099940000068
训练多项式逻辑回归分类器,其中ξ∈[0,1,...,M-1]表示当前过程处于协同训练的第ξ+1次迭代;
S8、使用多项式逻辑回归分类器对候选集
Figure GDA0002727099940000069
中的样本分类,得到第一种预测标签
Figure GDA00027270999400000610
计算多项式逻辑回归分类器对每个像素点的最大分类后验概率值,作为预测置信度P=[p1,p2,...,pγ];其中,γ代表候选集样本数,
Figure GDA0002727099940000071
代表多项式逻辑回归分类器对候选集
Figure GDA0002727099940000072
中第κ个样本的预测结果,pκ代表候选集
Figure GDA0002727099940000073
中第κ个样本的预测置信度,κ∈[1,2,...,γ];
S9、使用半监督图分类方法对整个候选集
Figure GDA0002727099940000074
中的样本分类:
Figure GDA0002727099940000075
得到第二种预测标签:
Figure GDA0002727099940000076
其中,γ代表候选集样本数,
Figure GDA0002727099940000077
代表半监督图分类方法对候选集
Figure GDA0002727099940000078
中第κ个样本的预测结果;α为训练参数,Y(0)由训练集
Figure GDA0002727099940000079
的样本标签
Figure GDA00027270999400000710
构成:
Figure GDA00027270999400000711
其中,i指代样本序号,c代表样本标签;
Figure GDA00027270999400000712
表示第i个样本的样本标签;[Γ12,...,Γδ]为训练集
Figure GDA00027270999400000713
中样本的样本序号;
S10、取出候选集
Figure GDA00027270999400000714
中预测标签Y1和预测标签Y2相同的样本,以及对应的预测标签
Figure GDA00027270999400000715
构成协定集
Figure GDA00027270999400000716
其中,uκ表示
Figure GDA00027270999400000717
中的一个样本;从候选集置信度P中选出协定集样本对应的所有置信系数作为综合置信集
Figure GDA00027270999400000718
由于PA的确定是以两种不同的预测标签Y1和Y2为基础的,因此将其称为综合置信集;S11、筛选出协定集DA对应综合置信集PA中置信系数高于99%的样本,赋予预测标签YN={YA(idx)|PA(idx)>0.99},其中,YA(idx)、PA(idx)分别表示YA、PA中序号为idx的元素,构成扩增集DN,加入训练集
Figure GDA00027270999400000719
中,得到更新的训练集
Figure GDA00027270999400000720
S12、将S11中扩增集DN的样本从候选集
Figure GDA00027270999400000721
中移除,得到更新的候选集
Figure GDA00027270999400000722
S13、迭代循环S7-S12步骤的高光谱图像协同训练过程,直至协同训练迭代达到设定值M次。
S14、使用半监督图分类方法和更新后的训练集
Figure GDA00027270999400000723
分类高光谱图像验证分类精度。
图2是对于AVIRIS Indian Pines图像,采用本实施例的基于综合置信的高光谱图像半监督分类方法的分类精度随协同训练迭代次数增长的示意图,其中初始训练样本为对图像12类样本每类随机取6个的情况;
图3是对于AVIRIS Indian Pines图像,采用本发明的基于综合置信的高光谱图像半监督分类方法的分类精度在初始训练样本分别取每类2个、4个、6个、8个不同情况下的变化示意图,其中协同训练共迭代10次;
图4(a)是对于AVIRIS Indian Pines图像,采用半监督图分类方法的分类精度在初始训练样本分别取每类2个、4个、6个、8个不同情况下的变化示意图;图4(b)是对于AVIRIS Indian Pines图像,采用多项式逻辑回归分类器的分类精度在初始训练样本分别取每类2个、4个、6个、8个不同情况下的变化示意图,其中协同训练均迭代10次;
图5(a)-图5(d)是训练样本标签,图像中的像素点被分为12类,图中训练样本是对AVIRIS Indian Pines图像,通过在每类分别采样2个、4个、6个、8个像素点得到的;
图6(a)-图6(d)是测试样本标签,图像中的像素点被分为12类,四个图中抽取的训练样本是对AVIRIS Indian Pines图像,通过在每类分别采样2个、4个、6个、8个像素点得到的;
图7(a)-图7(d)是协同训练扩增后的训练样本标签,采用本发明提出的基于综合置信的自学习方法在图5(a)-图5(d)所示初始训练样本基础上对AVIRIS Indian Pines图像自学习得到的扩增训练样本,四个图中的初始训练样本是通过在每类分别采样2个、4个、6个、8个像素点得到的。
图8(a)-图8(d)是全体样本标签,采用本发明提出的基于综合置信的高光谱图像半监督分类方法得到对AVIRIS Indian Pines图像的分类结果,图像中的像素点被分为12类,四个图中抽取的训练样本是通过在每类分别采样2个、4个、6个、8个像素点得到的;
图9(a)-图9(d)是全体样本标签,采用半监督图分类方法得到对AVIRIS IndianPines图像的分类结果,图像中的像素点被分为12类,四个图中抽取的训练样本是通过在每类分别采样2个、4个、6个、8个像素点得到的;
图10(a)-图10(d)是全体样本标签,采用多项式逻辑回归分类器得到对AVIRISIndian Pines图像的分类结果,图像中的像素点被分为12类,四个图中抽取的训练样本是通过在每类分别采样2个、4个、6个、8个像素点得到的;
由图8(a)-图8(d)得到表1利用基于综合置信的高光谱图像半监督分类方法以及图5(a)-图5(d)训练样本对AVIRIS Indian Pines高光谱图像测试样本分类,得到的准确率。其中,待分类的像素来自于12个类。
表1
每类训练样本个数 2 4 6 8
分类准确率(%) 67.27 76.74 83.50 86.75
由图9(a)-图9(d)得到表2利用半监督图分类方法以及图5(a)-图5(d)训练样本对AVIRIS Indian Pines高光谱图像测试样本分类,得到的准确率。其中,待分类的像素来自于12个类。
表2
每类训练样本个数 2 4 6 8
分类准确率(%) 59.75 64.24 68.99 74.70
由图10(a)-图10(d)得到表3利用多项式逻辑回归分类器以及图5(a)-图5(d)训练样本对AVIRIS Indian Pines高光谱图像测试样本分类,得到的准确率。其中,待分类的像素来自于12个类。
表3
每类训练样本个数 2 4 6 8
分类准确率(%) 37.06 52.90 56.49 57.94
本实施例所述方法的优势如下:
(1)更加符合高光谱图像数据特性:高光谱图像像素点是否属于同类不仅与光谱特征相似度有关,而且与空间位置的距离有关。光谱特征越相似的两个像素点越可能属于同一类,空间近邻的两个像素点也极有可能属于同一类。基于空间最近邻连接的半监督图分类方法虽然能够利用空间信息使空间邻接的像素产生图连接,从而得到空间上均匀成块分布的分类结果,但是这种方法缺乏了不同区域像素点之间的图连接,使得光谱特征极为接近但空间位置较远的像素点之间相关性不大。而基于光谱特征的判别式分类方法如多项式逻辑回归分类器却未利用空间信息,导致分类结果噪点明显。本发明提出的基于综合置信的高光谱图像半监督分类方法在距离初始训练样本较远的区域产生了更多较高综合置信度预测样本,并用于增强先验信息,相当于在原本只有空间近邻连接的图上新增了更多“远距离”连接。本发明提出的方法通过扩增更多训练样本的方式,将原本过分依赖空间信息的半监督图分类方法进行了增强,使改进后的图连接方式更加符合高光谱图像的数据特性。
(2)自学习挖掘出分类器的最大潜力:在半监督分类问题中,由于先验信息极少,连接测度的改进极难提升传统图分类方法的分类精度,也难以设计出优秀的图连接方式。因此,单纯靠半监督图分类方法,高光谱图像半监督分类任务仍十分受限于初始训练样本的数量和质量。本发明提出的基于综合置信的高光谱图像半监督分类方法利用协同训练的思想,通过半监督图分类方法和多项式逻辑回归分类器共同挑选新训练样本的扩增方式,实现了机器的自学习和分类效果的逐步提升。即使在初始训练样本极少的情况下,本发明方法仍旧能够找到更多益于分类的训练样本,挖掘出图分类方法和多项式逻辑回归分类器的最大潜力。
(3)协同训练效率较高:本发明的训练样本扩增过程采用两种互补的分类器,半监督图分类方法基于空间信息,预测的样本标签分布均匀成块,而多项式逻辑回归分类器基于光谱信息,预测的样本标签分布噪声较大但预测结果更符合样本的光谱特征分布。两种分类器预测意见一致的样本既能满足光谱特征分布,也能满足空间位置分布。再通过综合置信度做进一步筛选以保证扩增集的预测可靠性,这样的扩增集样本能够较为明显地改善半监督图分类方法的分类性能,进而提升协同训练过程的自学习效率。
以上所述,仅为本发明专利较佳的实施例,但本发明专利的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明专利所公开的范围内,根据本发明专利的技术方案及其发明专利构思加以等同替换或改变,都属于本发明专利的保护范围。

Claims (7)

1.一种基于综合置信的高光谱图像半监督分类方法,其特征在于,所述方法包括以下步骤:
S1、读入三维高光谱图像立方体H(m,n,b),其中m、n表示空间像素点位置,b表示所处光谱波段位置;
S2、计算高光谱数据像素点之间的去样本均值的相关系数,用于构造图权重矩阵W,权重数值能够衡量像素点间的相似度;
S3、将步骤S2的图权重矩阵W中每个像素点与自身以及与空间8近邻外的像素点的相似度权重值置0,即只为每个像素点与其空间8近邻内的像素点间建立图连接,得到稀疏化的图权重矩阵
Figure FDA0002727099930000011
S4、计算归一化图权重矩阵
Figure FDA0002727099930000012
S5、获取初始训练集
Figure FDA0002727099930000013
和候选集
Figure FDA0002727099930000014
S6、设定协同训练的迭代次数M,并开始协同训练过程;
S7、使用训练集
Figure FDA0002727099930000015
训练多项式逻辑回归分类器,ξ∈[0,1,...,M-1]表示当前过程处于协同训练的第ξ+1次迭代;
S8、使用多项式逻辑回归分类器对候选集
Figure FDA0002727099930000016
中的样本分类,得到预测标签Y1和对应的置信度P;
S9、使用半监督图分类方法对候选集
Figure FDA0002727099930000017
中的样本分类,得到预测标签Y2
S10、取出候选集
Figure FDA0002727099930000018
中预测标签Y1和预测标签Y2相同的样本,以及对应的预测标签YA,构成协定集DA,从候选集置信度P中选出协定集样本对应的所有置信系数作为综合置信集PA
S11、筛选出协定集DA对应综合置信集PA中置信系数高于99%的样本,赋予预测标签YN,构成扩增集DN,加入训练集
Figure FDA0002727099930000019
中,得到更新的训练集
Figure FDA00027270999300000110
S12、将S11中扩增集DN的样本从候选集
Figure FDA00027270999300000111
中移除,得到更新的候选集
Figure FDA00027270999300000112
S13、迭代循环S7-S12步骤的高光谱图像协同训练过程,直至协同训练迭代达到设定值M次;
S14、使用半监督图分类方法和更新后的训练集
Figure FDA00027270999300000113
分类高光谱图像验证分类精度。
2.根据权利要求1所述的基于综合置信的高光谱图像半监督分类方法,其特征在于,步骤S2中的图权重矩阵W采用如下计算方式:
Figure FDA0002727099930000021
其中,vi代表第i个像素点的数据特征,va代表所有像素点的数据特征均值,wij是图权重矩阵W中第i行、第j列的元素,同时代表第i个像素点和第j个像素点的相似性权重值。
3.根据权利要求1所述的基于综合置信的高光谱图像半监督分类方法,其特征在于,步骤S4中的归一化图权重矩阵的计算采用
Figure FDA0002727099930000022
其中度矩阵D计算公式如下:
Figure FDA0002727099930000023
其中,wij是图权重矩阵W中第i行、第j列的元素,V是像素点总数。
4.根据权利要求1所述的基于综合置信的高光谱图像半监督分类方法,其特征在于,步骤S8中的初始候选集
Figure FDA0002727099930000024
由整个测试集内的样本构成,通过多项式逻辑回归分类器对整个候选集样本分类,得到第一种预测标签:
Figure FDA0002727099930000025
其中,γ代表候选集样本数,
Figure FDA0002727099930000026
代表多项式逻辑回归分类器对候选集
Figure FDA0002727099930000027
中第κ个样本的预测结果,κ∈[1,2,...,γ];计算多项式逻辑回归分类器对每个像素点的最大分类后验概率值,作为预测置信度P=[p1,p2,...,pγ],其中pκ代表候选集
Figure FDA0002727099930000028
中第κ个样本的预测置信度。
5.根据权利要求1所述的基于综合置信的高光谱图像半监督分类方法,其特征在于,步骤S9中的半监督图分类方法采用如下计算方式:
Figure FDA0002727099930000029
其中,α为训练参数,Y(0)由训练集
Figure FDA00027270999300000210
的样本标签
Figure FDA00027270999300000211
构成:
Figure FDA00027270999300000212
其中,i指代样本序号,c代表样本标签;
Figure FDA0002727099930000031
表示第i个样本的样本标签;[Γ12,...,Γδ]为训练集
Figure FDA0002727099930000032
中样本的样本序号;
使用半监督图分类方法对整个候选集
Figure FDA0002727099930000033
中的样本分类,得到第二种预测标签:
Figure FDA0002727099930000034
其中,γ代表候选集样本数,
Figure FDA0002727099930000035
代表半监督图分类方法对候选集
Figure FDA0002727099930000036
中第κ个样本的预测结果,κ表示候选集
Figure FDA0002727099930000037
中样本的序号,κ∈[1,2,...,γ]。
6.根据权利要求1所述的基于综合置信的高光谱图像半监督分类方法,其特征在于,步骤S10中的协定集DA由候选集
Figure FDA0002727099930000038
中预测标签Y1和预测标签Y2相同的样本构成:
Figure FDA0002727099930000039
其中,uκ表示
Figure FDA00027270999300000310
中的一个样本;
Figure FDA00027270999300000311
是预测标签Y1中的第κ个样本的预测标签,即多项式逻辑回归分类器对候选集
Figure FDA00027270999300000312
中第κ个样本的预测结果;
Figure FDA00027270999300000313
是预测标签Y2中的第κ个样本的预测标签,即半监督图分类方法对候选集
Figure FDA00027270999300000314
中第κ个样本的预测结果,κ∈[1,2,...,γ];
同时取出协定集样本对应的预测标签
Figure FDA00027270999300000315
以及协定集样本对应的所有置信系数作为综合置信集
Figure FDA00027270999300000316
pκ代表候选集
Figure FDA00027270999300000317
中第κ个样本的预测置信度;由于PA的确定是以两种不同的预测标签Y1和Y2为基础的,因此将其称为综合置信集。
7.根据权利要求1所述的基于综合置信的高光谱图像半监督分类方法,其特征在于,步骤S11中的扩增集DN由协定集DA对应综合置信集PA>0.99的样本构成,扩增集DN的样本标签YN同样由YA中相应预测标签构成:
YN={YA(idx)|PA(idx)>0.99}
其中,YA(idx)、PA(idx)分别表示YA、PA中序号为idx的元素,而更新的训练集由之前的训练集和扩增集构成:
Figure FDA00027270999300000318
CN201910012238.9A 2019-01-07 2019-01-07 一种基于综合置信的高光谱图像半监督分类方法 Active CN109784392B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910012238.9A CN109784392B (zh) 2019-01-07 2019-01-07 一种基于综合置信的高光谱图像半监督分类方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910012238.9A CN109784392B (zh) 2019-01-07 2019-01-07 一种基于综合置信的高光谱图像半监督分类方法

Publications (2)

Publication Number Publication Date
CN109784392A CN109784392A (zh) 2019-05-21
CN109784392B true CN109784392B (zh) 2020-12-22

Family

ID=66499933

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910012238.9A Active CN109784392B (zh) 2019-01-07 2019-01-07 一种基于综合置信的高光谱图像半监督分类方法

Country Status (1)

Country Link
CN (1) CN109784392B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112183759B (zh) * 2019-07-04 2024-02-13 创新先进技术有限公司 模型训练方法、装置及系统
CN110321528B (zh) * 2019-07-11 2022-11-11 生态环境部南京环境科学研究所 一种基于半监督地理空间回归分析的高光谱影像土壤重金属浓度评估方法
CN110728187B (zh) * 2019-09-09 2022-03-04 武汉大学 一种基于容错性深度学习的遥感影像场景分类方法
CN110689064B (zh) * 2019-09-23 2023-06-09 云南电网有限责任公司电力科学研究院 一种图像半监督分类方法、装置及计算机可读存储介质
CN111144451B (zh) * 2019-12-10 2023-08-25 东软集团股份有限公司 一种图像分类模型的训练方法、装置及设备
CN111414942B (zh) * 2020-03-06 2022-05-03 重庆邮电大学 一种基于主动学习和卷积神经网络的遥感图像分类方法
CN111881989B (zh) * 2020-08-03 2022-09-09 齐齐哈尔大学 一种高光谱图像分类方法
CN112560960B (zh) * 2020-12-16 2024-08-13 北京影谱科技股份有限公司 一种高光谱图像分类方法、装置以及计算设备
CN112784774B (zh) * 2021-01-27 2022-08-23 山东农业大学 一种基于数据增强的小样本高光谱分类方法
CN113033641B (zh) * 2021-03-17 2022-12-16 华南理工大学 一种高维数据半监督分类方法
CN113869430A (zh) * 2021-09-29 2021-12-31 北京百度网讯科技有限公司 训练方法、图像识别方法、装置、电子设备以及存储介质
CN115410088B (zh) * 2022-10-10 2023-10-31 中国矿业大学 一种基于虚拟分类器的高光谱图像领域自适应方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107451616A (zh) * 2017-08-01 2017-12-08 西安电子科技大学 基于深度半监督迁移学习的多光谱遥感图像地物分类方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104376335A (zh) * 2014-11-13 2015-02-25 河南理工大学 一种基于信息熵的半监督高光谱遥感影像分类方法
CN107644235A (zh) * 2017-10-24 2018-01-30 广西师范大学 基于半监督学习的图像自动标注方法
CN108960282A (zh) * 2018-05-30 2018-12-07 昆明理工大学 一种基于半监督学习的在线服务信誉度量方法
CN109145950B (zh) * 2018-07-20 2021-09-21 华南理工大学 一种基于图信号采样的高光谱图像主动学习方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107451616A (zh) * 2017-08-01 2017-12-08 西安电子科技大学 基于深度半监督迁移学习的多光谱遥感图像地物分类方法

Also Published As

Publication number Publication date
CN109784392A (zh) 2019-05-21

Similar Documents

Publication Publication Date Title
CN109784392B (zh) 一种基于综合置信的高光谱图像半监督分类方法
CN106096561B (zh) 基于图像块深度学习特征的红外行人检测方法
CN110084159B (zh) 基于联合多级空谱信息cnn的高光谱图像分类方法
CN109684922B (zh) 一种基于卷积神经网络的多模型对成品菜的识别方法
CN111914907A (zh) 一种基于深度学习空谱联合网络的高光谱图像分类方法
CN109871875B (zh) 一种基于深度学习的建筑物变化检测方法
CN109766858A (zh) 结合双边滤波的三维卷积神经网络高光谱影像分类方法
CN110717526A (zh) 一种基于图卷积网络的无监督迁移学习方法
CN111639587B (zh) 基于多尺度谱空卷积神经网络的高光谱图像分类方法
CN109858557B (zh) 一种新的高光谱图像数据半监督分类方法
CN107590515A (zh) 基于熵率超像素分割的自编码器的高光谱图像分类方法
CN111160396B (zh) 一种基于多图结构的图卷积网络的高光谱图像分类方法
CN108229551B (zh) 一种基于紧凑字典稀疏表示的高光谱遥感图像分类方法
CN110414616B (zh) 一种利用空间关系的遥感图像字典学习分类方法
CN105160351B (zh) 基于锚点稀疏图的半监督高光谱分类方法
CN112464983A (zh) 一种用于苹果树叶病害图像分类的小样本学习方法
CN111639697B (zh) 基于非重复采样与原型网络的高光谱图像分类方法
CN108268890A (zh) 一种高光谱图像分类方法
CN112364747B (zh) 一种有限样本下的目标检测方法
CN106485238A (zh) 一种高光谱遥感图像特征提取和分类方法及其系统
CN108460400A (zh) 一种结合多种特征信息的高光谱图像分类方法
CN113052216A (zh) 基于双路图u-net卷积网络的溢油高光谱图像检测方法
Guo et al. Dual graph U-Nets for hyperspectral image classification
CN115240072A (zh) 一种基于多方向多尺度光谱-空间残差卷积神经网络的高光谱多类变化检测方法
CN109145950B (zh) 一种基于图信号采样的高光谱图像主动学习方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant