CN106469316A - 基于超像素级信息融合的高光谱图像的分类方法及系统 - Google Patents
基于超像素级信息融合的高光谱图像的分类方法及系统 Download PDFInfo
- Publication number
- CN106469316A CN106469316A CN201610810465.2A CN201610810465A CN106469316A CN 106469316 A CN106469316 A CN 106469316A CN 201610810465 A CN201610810465 A CN 201610810465A CN 106469316 A CN106469316 A CN 106469316A
- Authority
- CN
- China
- Prior art keywords
- super
- pixel
- characteristic
- dimension
- spectrum image
- 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.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 142
- 238000000034 method Methods 0.000 title claims abstract description 96
- 230000004927 fusion Effects 0.000 title claims abstract description 18
- 230000009467 reduction Effects 0.000 claims abstract description 36
- 230000011218 segmentation Effects 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000000605 extraction Methods 0.000 claims abstract description 13
- 238000013507 mapping Methods 0.000 claims description 15
- 238000012549 training Methods 0.000 claims description 14
- 238000012360 testing method Methods 0.000 claims description 12
- 238000012706 support-vector machine Methods 0.000 claims description 11
- 238000013480 data collection Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 abstract description 2
- 230000003595 spectral effect Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 241001466077 Salina Species 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 235000008331 Pinus X rigitaeda Nutrition 0.000 description 1
- 235000011613 Pinus brutia Nutrition 0.000 description 1
- 241000018646 Pinus brutia Species 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000000686 essence Substances 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical compound C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/7715—Feature extraction, e.g. by transforming the feature space, e.g. multi-dimensional scaling [MDS]; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/254—Fusion techniques of classification results, e.g. of results related to same input data
Abstract
本发明适用于图像处理领域,提供了一种基于超像素级信息融合的高光谱图像的分类方法及系统,旨在解决现有的分类方法计算复杂度高、分类精度较低以及波段间存在较大冗余的问题。所述方法包括下述步骤:滤波器生成步骤、Gabor特征提取步骤、高光谱图像超像素分割步骤、超像素特征数据计算步骤、超像素空间坐标计算步骤、数据降维步骤、高光谱特征数据生成步骤以及多任务支持向量机分类步骤。
Description
技术领域
本发明属于图像处理领域,尤其涉及一种基于超像素级信息融合的高光谱图像的分类方法及系统。
背景技术
本发明涉及到的是一种利用高光谱图像对地面物质进行分类的技术。高光谱图像是由遥感传感器在电磁波谱的可见光,近红外,中红外和热红外波段范围内,从地面上感兴趣的物质上获取的多光谱影像数据。高光谱图像光谱分辨率和维数的大量增加使得更精确和精细的分类成为可能。然而,高光谱特性和分类研究中主要存在以下两个难点:一是高维使得提高小样本分类的精度及其困难,即所谓的“维数灾难”现象;二是波段维数高使得计算量极具增大,而且波段间的强相关性增加了冗余性,如果不进行有效处理,会对结果有影响。
传统的分类方法(K近邻,支持向量机,基于稀疏表示的分类)在高光谱上只利用光谱数据直接进行分类,无法满足实际的分类效果。为了解决以上难点,谱域-空域结合的高光谱图像分类技术应运而生。鉴于高光谱图像具备空间、光谱信息,像素之间隐藏着丰富的有助于分类的信号变化信息,一个有效的方法就是希望能够提取出更强鉴别力的空间-光谱组合特征,从而提高分类精度。
当前,除了结合谱域-空域来提高分类性能,多任务的分类方法也被广泛采纳。多任务的分类方法主要表现在两个方面,一个是联合多种分类器进行的分类,该种方法认为不同的分类器具有不同的决策性能,综合在一起使用可以减少分类精度的方差,从而提高分类系统的性能。另一个方面则是联合多种特征进行的分类,该种方法通过各种特征的不同描述,自适应的补足特征之间的差异,使得对类别的判别能力更强。
然而,单纯的使用上述方法无法排除高光谱图像波段之间的冗余性,当前,解决以上问题的主要方法就是使用数据降维。现有的降维方法可以分为两类:一类是基于变换的方法,如主成分分析(PCA)、正交子空间投影(OSP)、正则分析(CA)、离散小波变换(DWT),薛定谔特征映射(SE)等。基于变换的降维方式,其优点是可以经过若干变换直接将高维数据降低到低维甚至一维,并且在降维过程中还可以融合空间与光谱信息,得到分辨力更强的特征,从而为分类精度的提高提供了可能;缺点就是改变了图像原有的特性。另一类是基于非变换的,如波段选择,数据源划分等。基于非变换的降维方式是在考察图像整体特点之后对图像进行选择和划分的,它的优点是保持了图像原有的特性;缺点就是波段选择后的图像仍然无法满足实际的分类效果。
目前,多任务的分类技术已经取得了一定的科研成果。采用了三维Gabor特征提取的方法选择性的得到多种高光谱数据特征,然后分别使用稀疏表示进行编码并重构得到各特征的重构误差,再将各特征的重构误差以线性加权的方式进行融合。最后,通过融合后的重构误差进行分类。该方法在分类阶段对各Gabor特征进行了融合,提高了分类精度,但提取出的Gabor特征仍然具备极大的冗余性,而且使用稀疏表示的方法进行多任务分类使得计算复杂度极高,短时内难以完成大量像素的分类。
除了多任务的分类技术,采用空谱结合的数据降维的方法也取得了一定的科研成果。使用超像素作用于薛定谔特征映射方法来对高光谱数据进行降维,加快了降维速度,但却由于特征单一,使用超像素级的降维难以达到小样本的分类精度要求。
发明内容
本发明所要解决的技术问题在于提供一种基于超像素级信息融合的高光谱图像的分类方法及系统,旨在解决现有的分类方法计算复杂度高、分类精度较低以及波段间存在较大冗余的问题。
为解决上述技术问题,本发明是这样实现的,本发明提供了一种基于超像素级信息融合的高光谱图像的分类方法,包括下述步骤:
滤波器生成步骤:生成若干个二维的Gabor滤波器;
Gabor特征提取步骤:将每一个所述Gabor滤波器分别与高光谱图像中的每一个波段进行卷积运算,并对卷积运算结果进行取幅值运算,得到若干个Gabor特征块;
高光谱图像超像素分割步骤:对所述高光谱图像进行超像素分割,得到若干个超像素;
超像素特征数据计算步骤:将每一个所述超像素与每一个所述Gabor特征块分别进行均值计算,得到若干个第一维度的超像素特征数据;
超像素空间坐标计算步骤:将每一个所述超像素与所述高光谱图像分别进行坐标均值计算,得到1个空间坐标数据集合;
数据降维步骤:将每一个所述第一维度的超像素特征数据分别与所述空间坐标数据集合使用空谱结合的薛定谔特征映射方法进行特征降维,从第一维度降到第二维度,得到若干个第二维度的超像素特征数据;
高光谱特征数据生成步骤:将每一个所述第二维度的超像素特征数据分别与所述空间坐标数据集合使用自然邻点插值法进行重构,得到若干个三维高光谱特征数据;
多任务支持向量机分类步骤:对所述每一个三维高光谱特征数据分别进行多任务的支持向量机分类。
进一步地,所述Gabor特征提取步骤包括:
将每一个所述二维Gabor滤波器与所述高光谱图像的每一个波段分别进行卷积运算,并对卷积运算结果根据下述公式进行取幅值运算,得到若干个Gabor特征块:
其中,表示若干个二维Gabor滤波器集合,表示第t个Gabor滤波器,(x,y)表示在二维平面进行卷积运算时对应的二元坐标变量,R表示所述高光谱图像,其中λ表示所述高光谱图像的每一波段,l表示高光谱图像的宽度,m表示高光谱图像的长度,B表示高光谱图像的波段数,即高光谱图像的高度,l×m×B表示三维,{Mt,t=1,2,..,X}表示若干个Gabor特征块集合,Mt表示第t个Gabor特征块,其中t表示个数,X为正整数。
进一步地,所述超像素特征数据计算步骤包括:{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt分别进行均值计算,得到n×B维的超像素特征数据Nt,最终共得到若干个n×B维的超像素特征数据集合{Nt,t=1,2,..,X};
在{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt进行均值计算时,S1对应Mt进行均值计算得到第1个B维的向量,S2对应Mt进行均值计算得到第2个B维的向量,S3对应Mt进行均值计算得到第3个B维的向量,依次类推至Sn对应Mt进行均值计算得到第n个B维的向量,最终共得到n个B维的向量,即n×B维的超像素特征数据Nt;
其中,{Si,i=1,2,..,n}表示对高光谱图像进行超像素分割得到的若干个超像素集合,Si表示第i个超像素,n表示超像素的个数;B维表示所述第一维度,{Nt,t=1,2,..,X}表示若干个n×B维的超像素特征数据集合,Nt表示第t个超像素特征数据;
其中,所述每一个超像素Si均包含若干个像素;
所述超像素空间坐标计算步骤包括:将所述每一个超像素Si分别与所述高光谱图像R进行坐标均值计算,得到1个n×2维的空间坐标数据集合C;其中,C表示空间坐标数据集合,
所述数据降维步骤包括:每一个超像素特征数据Nt分别与空间坐标数据集合C使用空谱结合的薛定谔特征映射的方法进行特征降维,从B维降到K维,得到K维的超像素特征数据Dt,最终共得到若干个K维的超像素特征数据集合{Dt,t=1,2,...,X};
其中,{Dt,t=1,2,..,X}表示维度为K的超像素特征数据集合,K表示所述第二维度,Dt表示第t个K维的超像素特征数据;所述每一个超像素特征数据Dt所对应的空间坐标数据集合仍是C。
进一步地,所述高光谱特征数据生成步骤包括:每一个超像素特征数据Dt分别和空间坐标数据集合C使用自然邻点插值法进行重构,对应原高光谱图像补足所有空间坐标的像素值,得到三维高光谱特征数据Gt,最终共得到若干个三维高光谱特征数据集合{Gt,t=1,2,..,X};
其中,{Gt,t=1,2,..,X}表示若干个三维的高光谱特征数据集合,l×m×K表示三维,Gt表示第t个三维的高光谱特征数据;
所述多任务支持向量机分类包括:将每一个高光谱特征数据Gt划分为训练数据G1t和测试数据G2t,g表示一个原始高光谱图像测试数据样本,其中g∈R,gt表示Gt当中与g同一位置坐标的像素特征数据,gt维度为K,{gt∈G2t,t=1,2,..,X}表示X个维度为K的像素特征数据集合,那么对g的分类过程如下:
(1)对训练数据{G1t,t=1,2,..,X}使用概率输出的支持向量机方法进行模型训练,得到概率输出模型{Modelt,t=1,2,..,X};
(2)使用概率输出模型Modelt对数据gt进行类别概率输出,得到gt属于每一个类别的概率{Pt(i),i=1,2,..,C},其中C为类别总数;
(3)高光谱图像样本g的类别预测公式为:
本发明还提供了一种基于超像素级信息融合的高光谱图像的分类系统,所述系统包括:
滤波器生成模块:用于生成若干个二维的Gabor滤波器;
Gabor特征提取模块:用于将每一个所述Gabor滤波器分别与高光谱图像中的每一个波段进行卷积运算,并对卷积运算结果进行取幅值运算,得到若干个Gabor特征块;
高光谱图像超像素分割模块:用于对所述高光谱图像进行超像素分割,得到若干个超像素;
超像素特征数据计算模块:用于将每一个所述超像素与每一个所述Gabor特征块分别进行均值计算,得到若干个第一维度的超像素特征数据;
超像素空间坐标计算模块:用于将每一个所述超像素与所述高光谱图像分别进行坐标均值计算,得到1个空间坐标数据集合;
数据降维模块:用于将每一个所述第一维度的超像素特征数据分别与所述空间坐标数据集合使用空谱结合的薛定谔特征映射方法进行特征降维,从第一维度降到第二维度,得到若干个第二维度的超像素特征数据;
高光谱特征数据生成模块:用于将每一个所述第二维度的超像素特征数据分别与所述空间坐标数据集合使用自然邻点插值法进行重构,得到若干个三维高光谱特征数据;
多任务支持向量机分类模块:用于对所述每一个三维高光谱特征数据分别进行多任务的支持向量机分类。
进一步地,所述Gabor特征提取模块具体用于:
将每一个所述二维Gabor滤波器与所述高光谱图像的每一个波段分别进行卷积运算,并对卷积运算结果根据下述公式进行取幅值运算,得到若干个Gabor特征块:
其中,表示若干个二维Gabor滤波器集合,表示第t个Gabor滤波器,(x,y)表示在二维平面进行卷积运算时对应的二元坐标变量,R表示所述高光谱图像,其中λ表示所述高光谱图像的每一波段,l表示高光谱图像的宽度,m表示高光谱图像的长度,B表示高光谱图像的波段数,即高光谱图像的高度,l×m×B表示三维,{Mt,t=1,2,..,X}表示若干个Gabor特征块集合,Mt表示第t个Gabor特征块,其中t表示个数,X为正整数。
进一步地,所述超像素特征数据计算模块具体用于:{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt分别进行均值计算,得到n×B维的超像素特征数据Nt,最终共得到若干个n×B维的超像素特征数据集合{Nt,t=1,2,..,X};
在{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt进行均值计算时,S1对应Mt进行均值计算得到第1个B维的向量,S2对应Mt进行均值计算得到第2个B维的向量,S3对应Mt进行均值计算得到第3个B维的向量,依次类推至Sn对应Mt进行均值计算得到第n个B维的向量,最终共得到n个B维的向量,即n×B维的超像素特征数据Nt;
其中,{Si,i=1,2,..,n}表示对高光谱图像进行超像素分割得到的若干个超像素集合,Si表示第i个超像素,n表示超像素的个数;B维表示所述第一维度,{Nt,t=1,2,..,X}表示若干个n×B维的超像素特征数据集合,Nt表示第t个超像素特征数据;
其中,所述每一个超像素Si均包含若干个像素;
所述超像素空间坐标计算模块具体用于:将所述每一个超像素Si分别与所述高光谱图像R进行坐标均值计算,得到1个n×2维的空间坐标数据集合C;其中,C表示空间坐标数据集合,
所述数据降维模块具体用于:每一个超像素特征数据Nt分别与空间坐标数据集合C使用空谱结合的薛定谔特征映射的方法进行特征降维,从B维降到K维,得到K维的超像素特征数据Dt,最终共得到若干个K维的超像素特征数据集合{Dt,t=1,2,..,X};
其中,{Dt,t=1,2,..,X}表示维度为K的超像素特征数据集合,K表示所述第二维度,Dt表示第t个K维的超像素特征数据;所述每一个超像素特征数据Dt所对应的空间坐标数据集合仍是C。
进一步地,所述高光谱特征数据生成模块具体用于:每一个超像素特征数据Dt分别和空间坐标数据集合C使用自然邻点插值法进行重构,对应原高光谱图像补足所有空间坐标的像素值,得到三维高光谱特征数据Gt,最终共得到若干个三维高光谱特征数据集合{Gt,t=1,2,..,X};
其中,{Gt,t=1,2,..,X}表示若干个三维的高光谱特征数据集合,l×m×K表示三维,Gt表示第t个三维的高光谱特征数据;
所述多任务支持向量机分类模块具体用于:将每一个高光谱特征数据Gt划分为训练数据G1t和测试数据G2t,g表示一个原始高光谱图像测试数据样本,其中g∈R,gt表示Gt当中与g同一位置坐标的像素特征数据,gt维度为K,{gt∈G2t,t=1,2,..,X}表示X个维度为K的像素特征数据集合,那么对g的分类过程如下:
(1)对训练数据{G1t,t=1,2,..,X}使用概率输出的支持向量机方法进行模型训练,得到概率输出模型{Modelt,t=1,2,..,X};
(2)使用概率输出模型Modelt对数据gt进行类别概率输出,得到gt属于每一个类别的概率{Pt(i),i=1,2,..,C},其中C为类别总数;
(3)高光谱图像样本g的类别预测公式为:
本发明与现有技术相比,有益效果在于:
针对现有的多任务稀疏表示分类方法计算复杂且计算量大的缺点,本发明采用基于多任务支持向量机分类的方法,大大降低了计算的复杂度;
针对现有的光谱数据在超像素级使用空谱结合的薛定谔降维后小样本分类精度不高的问题,本发明采用基于二维Gabor的多种特征结合超像素的空谱结合的薛定谔特征降维方法,具有更高的分类精度。
针对现有的波段间存在极大的冗余性问题,本发明使用的Gabor特征块包含了更加丰富的局部变化信息,使用数据降维的方法减少了波段之间的冗余信息。
附图说明
图1是本发明实施例提供的基于超像素级信息融合的高光谱图像的分类方法流程示意图;
图2是本发明实施例提供的不同频率和方向的滤波器示意图;
图3是本发明实施例提供的Gabor特征获取示意图;
图4是本发明实施例提供的Gabor滤波器对高光谱某一波段的图像进行卷积的结果示意图;
图5是本发明实施例提供的使用SLIC方法将高光谱图像分割成超像素的平面示意图;
图6是多任务支持向量机分类示意图;
图7是本发明实施例提供的基于超像素级信息融合的高光谱图像的分类系统示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明实施例提供了一种基于超像素级信息融合的高光谱图像的分类方法,如图1所示,包括下述步骤:
步骤S101:生成若干个二维的Gabor滤波器。
在本发明实施例中,采用了以下方法生成Gabor滤波器,生成Gabor滤波器的公式为:
其中,x′=xcosθv+ysinθv,y′=-xsinθv+ycosθv,a=0.9589fu,b=1.1866fu,fu表示Gabor滤波器的频率,θv表示Gabor滤波器的方向,(x,y)表示Gabor滤波器对应的二元变量。
本发明实施例设计了4个fu=[0.03589,0.09473,0.25,0.6577]和6个θv=[0,40,80,120,160,180],按照上述生成Gabor滤波器的公式计算共生成了24个二维的Gabor滤波器其中,表示24个二维Gabor滤波器集合,表示第t个Gabor滤波器。如图2所示,为本发明实施例生成的24个不同频率和方向的滤波器示意图。
在本发明实施例中,步骤S101采用了上述方法生成Gabor滤波器,但并不能限定本发明,步骤S101中的Gabor滤波器,可以采用比如Log-Gabor等不同的形式,生成滤波器的频率、角度等参数、个数也可以调整,本发明所采用的Gabor响应也可以通过幅值、相位等其他形式代替。
步骤S102:将每一个所述Gabor滤波器分别与高光谱图像中的每一个波段进行卷积运算,并对卷积运算结果进行取幅值运算,得到若干个Gabor特征块;
本发明实施例根据步骤S101的结果,将每一个所述二维Gabor滤波器与所述高光谱图像的每一个波段分别进行卷积运算,并对卷积运算结果根据下述公式进行取幅值运算,得到24个Gabor特征块:
其中,表示24个二维Gabor滤波器集合,表示第t个Gabor滤波器,(x,y)表示在二维平面进行卷积运算时对应的二元坐标变量,R表示所述高光谱图像,其中λ表示所述高光谱图像的每一个波段,l表示高光谱图像的宽度,m表示高光谱图像的长度,B表示高光谱图像的波段数,即高光谱图像的高度,l×m×B表示三维,{Mt,t=1,2,..,24}表示24个Gabor特征块集合,Mt表示第t个Gabor特征块。
如图3和图4所示,为本发明实施例提供的在步骤S102得到的Gabor特征获取示意图和Gabor滤波器对高光谱某一波段的图像进行卷积的结果示意图。
步骤S103:对所述高光谱图像进行超像素分割,得到若干个超像素;
本发明实施例使用SLIC(simple linear iterative clustering,简单线性迭代聚类)的超像素分割方法对高光谱图像进行分割,得到分割图,该分割图包含n个超像素{Si,i=1,2,..,n}。
在本发明实施例中使用SLIC方法对高光谱图像进行分割,并不能限定本发明,步骤S103还可以使用基于熵率的超像素分割等方法实现。如图5所示,为本发明实施例提供的使用SLIC方法将高光谱图像分割成超像素的平面示意图。
步骤S104:将每一个所述超像素与每一个所述Gabor特征块分别进行均值计算,得到若干个第一维度的超像素特征数据;
本发明实施例中,{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt分别进行均值计算,得到n×B维的超像素特征数据Nt,最终共得到24个n×B维的超像素特征数据集合{Nt,t=1,2,..,24}。
在{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt进行均值计算时,S1对应Mt进行均值计算得到第1个B维的向量,S2对应Mt进行均值计算得到第2个B维的向量,S3对应Mt进行均值计算得到第3个B维的向量,依次类推至Sn对应Mt进行均值计算得到第n个B维的向量,最终共得到n个B维的向量,即n×B维的超像素特征数据Nt。
其中,{Si,i=1,2,..,n}表示对高光谱图像进行超像素分割得到的若干个超像素集合,Si表示第i个超像素,n表示超像素的个数;B维表示所述第一维度,{Nt,t=1,2,..,24}表示24个n×B维的超像素特征数据集合,Nt表示第t个超像素特征数据;其中,所述每一个超像素Si均包含若干个像素。
步骤S105:将每一个所述超像素与所述高光谱图像分别进行坐标均值计算,得到1个空间坐标数据集合。
本发明实施例将所述每一个超像素Si分别与所述高光谱图像R进行坐标均值计算,得到1个n×2维的空间坐标数据集合C;其中,C表示空间坐标数据集合,
步骤S106:将每一个所述第一维度的超像素特征数据分别与所述空间坐标数据集合使用空谱结合的薛定谔特征映射方法进行特征降维,从第一维度降到第二维度,得到若干个第二维度的超像素特征数据。
本发明实施例中,将每一个超像素特征数据Nt分别与空间坐标数据集合C使用空谱结合的薛定谔特征映射的方法进行特征降维,从B维降到K维,得到K维的超像素特征数据Dt,最终共得到24个K维的超像素特征数据集合{Dt,t=1,2,..,24};
其中,{Dt,t=1,2,..,24}表示维度为K的超像素特征数据集合,K表示所述第二维度,Dt表示第t个K维的超像素特征数据;所述每一个超像素特征数据Dt所对应的空间坐标数据集合仍是C。
步骤S107:将每一个所述第二维度的超像素特征数据分别与所述空间坐标数据集合使用自然邻点插值法进行重构,得到若干个三维高光谱特征数据。
本发明实施例中,每一个超像素特征数据Dt分别和空间坐标数据集合C使用自然邻点插值法进行重构,对应原高光谱图像补足所有空间坐标的像素值,得到三维高光谱特征数据Gt,最终共得到24个三维高光谱特征数据集合{Gt,t=1,2,...,24}。
其中,{Gt,t=1,2,..,24}表示24个三维的高光谱特征数据集合,l×m×K表示三维,Gt表示第t个三维的高光谱特征数据。
步骤S108:对所述每一个三维高光谱特征数据分别进行多任务的支持向量机分类。
本发明实施例中,将每一个高光谱特征数据Gt划分为训练数据G1t和测试数据G2t,g表示一个原始高光谱图像测试数据样本,其中g∈R,gt表示Gt当中与g同一位置坐标的像素特征数据,gt维度为K,{gt∈G2t,t=1,2,..,24}表示24个维度为K的像素特征数据集合,那么对g的分类过程如下:
(1)对训练数据{G1t,t=1,2,..,24}使用概率输出的支持向量机方法进行模型训练,得到概率输出模型{Modelt,t=1,2,..,24};
(2)使用概率输出模型Modelt对数据gt进行类别概率输出,得到gt属于每一个类别的概率{Pt(i),i=1,2,..,C},其中C为类别总数;
(3)高光谱图像样本g的类别预测公式为:
如图6所示,为本发明实施例提供的在步骤S108得到的多任务支持向量机分类示意图。
如图7所示,本发明实施例还提供了一种基于超像素级信息融合的高光谱图像的分类系统,所述系统包括:
滤波器生成模块701:用于生成若干个二维的Gabor滤波器;
Gabor特征提取模块702:用于将每一个所述Gabor滤波器分别与高光谱图像中的每一个波段进行卷积运算,并对卷积运算结果进行取幅值运算,得到若干个Gabor特征块;
高光谱图像超像素分割模块703:用于对所述高光谱图像进行超像素分割,得到若干个超像素;
超像素特征数据计算模块704:用于将每一个所述超像素与每一个所述Gabor特征块分别进行均值计算,得到若干个第一维度的超像素特征数据;
超像素空间坐标计算模块705:用于将每一个所述超像素与所述高光谱图像分别进行坐标均值计算,得到1个空间坐标数据集合;
数据降维模块706:用于将每一个所述第一维度的超像素特征数据分别与所述空间坐标数据集合使用空谱结合的薛定谔特征映射方法进行特征降维,从第一维度降到第二维度,得到若干个第二维度的超像素特征数据;
高光谱特征数据生成模块707:用于将每一个所述第二维度的超像素特征数据分别与所述空间坐标数据集合使用自然邻点插值法进行重构,得到若干个三维高光谱特征数据;
多任务支持向量机分类模块708:用于对所述每一个三维高光谱特征数据分别进行多任务的支持向量机分类。
进一步地,所述Gabor特征提取模块702具体用于:将每一个所述二维Gabor滤波器与所述高光谱图像的每一波段分别进行卷积运算,并对卷积运算结果根据下述公式进行取幅值运算,得到若干个Gabor特征块:
其中,表示若干个二维Gabor滤波器集合,表示第t个Gabor滤波器,(x,y)表示在二维平面进行卷积运算时对应的二元坐标变量,R表示所述高光谱图像,其中λ表示所述高光谱图像的每一波段,l表示高光谱图像的宽度,m表示高光谱图像的长度,B表示高光谱图像的波段数,即高光谱图像的高度,l×m×B表示三维,{Mt,t=1,2,..,X}表示若干个Gabor特征块集合,Mt表示第t个Gabor特征块,其中t表示个数,X为正整数。
进一步地,所述超像素特征数据计算模块704具体用于:{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt分别进行均值计算,得到n×B维的超像素特征数据Nt,最终共得到若干个n×B维的超像素特征数据集合{Nt,t=1,2,..,X}。
在{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt进行均值计算时,S1对应Mt进行均值计算得到第1个B维的向量,S2对应Mt进行均值计算得到第2个B维的向量,S3对应Mt进行均值计算得到第3个B维的向量,依次类推至Sn对应Mt进行均值计算得到第n个B维的向量,最终共得到n个B维的向量,即n×B维的超像素特征数据Nt。
其中,{Si,i=1,2,..,n}表示对高光谱图像进行超像素分割得到的若干个超像素集合,Si表示第i个超像素,n表示超像素的个数;B维表示所述第一维度,{Nt,t=1,2,..,X}表示若干个n×B维的超像素特征数据集合,Nt表示第t个超像素特征数据;其中,所述每一个超像素Si均包含若干个像素。
所述超像素空间坐标计算模块705具体用于:将所述每一个超像素Si分别与所述高光谱图像R进行坐标均值计算,得到1个n×2维的空间坐标数据集合C;其中,C表示空间坐标数据集合,
所述数据降维模块706具体用于:每一个超像素特征数据Nt分别与空间坐标数据集合C使用空谱结合的薛定谔特征映射的方法进行特征降维,从B维降到K维,得到K维的超像素特征数据Dt,最终共得到若干个K维的超像素特征数据集合{Dt,t=1,2,..,X}。
其中,{Dt,t=1,2,..,X}表示维度为K的超像素特征数据集合,K表示所述第二维度,Dt表示第t个K维的超像素特征数据;所述每一个超像素特征数据Dt所对应的空间坐标数据集合仍是C。
进一步地,所述高光谱特征数据生成模块707具体用于:每一个超像素特征数据Dt分别和空间坐标数据集合C使用自然邻点插值法进行重构,对应原高光谱图像补足所有空间坐标的像素值,得到三维高光谱特征数据Gt,最终共得到若干个三维高光谱特征数据集合{Gt,t=1,2,..,X};其中,{Gt,t=1,2,..,X}表示若干个三维的高光谱特征数据集合,l×m×K表示三维,Gt表示第t个三维的高光谱特征数据。
所述多任务支持向量机分类模块708具体用于:将每一个高光谱特征数据Gt划分为训练数据G1t和测试数据G2t,g表示一个原始高光谱图像测试数据样本,其中g∈R,gt表示Gt当中与g同一位置坐标的像素特征数据,gt维度为K,{gt∈G2t,t=1,2,..,X}表示X个维度为K的像素特征数据集合,那么对g的分类过程如下:
(1)对训练数据{G1t,t=1,2,..,X}使用概率输出的支持向量机方法进行模型训练,得到概率输出模型{Modelt,t=1,2,..,X};
(2)使用概率输出模型Modelt对数据gt进行类别概率输出,得到gt属于每一个类别的概率{Pt(i),i=1,2,..,C},其中C为类别总数;
(3)高光谱图像样本g的类别预测公式为:
在本发明实施例中,采用了三个真实的高光谱数据集。
第一个数据集是Indian Pines,是由AVIRIS高光谱传感器从美国印第安纳州的一块测试地上获得的,该图像大小为145*145共21025个像素,一共有224个波段,在实际应用时去掉了4个零波段和35个杂波段,剩余185个波段。图像的空间分辨率约为20m.该数据包括16个地物类别,共10249个有标记样本点。
第二个数据是Salinas,该数据是由AVIRIS传感器在加利福尼亚萨利纳斯山谷的上空采集,总共有512*217个样本,其中地物共54129个样本,包含16类地物,其余为背景,由于污染去掉了20个谱段,剩余204个。
第三个数据是PaviaU,该数据是由ROSIS传感器从意大利北部的帕维亚上空得到的,其空间分辨率为1.3m每像素,大小为610*340,共103个波段,包含9类地物,总共207400个样本,其中地物42776个,背景164624个。
以每类15个训练样本的PaviaU数据为例,本发明能够达到91.75%的精度,传统的支持向量机核方法精度为70.11%,超像素的薛定谔特征映射特征提取加支持向量机核的方法精度为85.25%,形态学特征提取加支持向量机核的分类方法为精度为81.18%,Gabor特征提取加多任务的稀疏表示方法精度为83.00%。对比表明,本发明的方法在分类精度上远远超过传统的分类方法。
综上所述,本发明采用基于多任务支持向量机分类的方法,大大降低了计算的复杂度;本发明采用基于二维Gabor的多种特征结合超像素的空谱结合的薛定谔特征降维方法,具有更高的分类精度;本发明使用的Gabor特征块包含了更加丰富的局部变化信息,使用数据降维的方法减少了波段之间的冗余信息。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于超像素级信息融合的高光谱图像的分类方法,其特征在于,所述分类方法包括下述步骤:
滤波器生成步骤:生成若干个二维的Gabor滤波器;
Gabor特征提取步骤:将每一个所述Gabor滤波器分别与高光谱图像中的每一个波段进行卷积运算,并对卷积运算结果进行取幅值运算,得到若干个Gabor特征块;
高光谱图像超像素分割步骤:对所述高光谱图像进行超像素分割,得到若干个超像素;
超像素特征数据计算步骤:将每一个所述超像素与每一个所述Gabor特征块分别进行均值计算,得到若干个第一维度的超像素特征数据;
超像素空间坐标计算步骤:将每一个所述超像素与所述高光谱图像分别进行坐标均值计算,得到1个空间坐标数据集合;
数据降维步骤:将每一个所述第一维度的超像素特征数据分别与所述空间坐标数据集合使用空谱结合的薛定谔特征映射方法进行特征降维,从第一维度降到第二维度,得到若干个第二维度的超像素特征数据;
高光谱特征数据生成步骤:将每一个所述第二维度的超像素特征数据分别与所述空间坐标数据集合使用自然邻点插值法进行重构,得到若干个三维高光谱特征数据;
多任务支持向量机分类步骤:对所述每一个三维高光谱特征数据分别进行多任务的支持向量机分类。
2.如权利要求1所述的基于超像素级信息融合的高光谱图像的分类方法,其特征在于,所述Gabor特征提取步骤包括:
将每一个所述二维Gabor滤波器与所述高光谱图像的每一个波段分别进行卷积运算,并对卷积运算结果根据下述公式进行取幅值运算,得到若干个Gabor特征块:
其中,表示若干个二维Gabor滤波器集合,表示第t个Gabor滤波器,(x,y)表示在二维平面进行卷积运算时对应的二元坐标变量,R表示所述高光谱图像,其中λ表示所述高光谱图像的每一波段,l表示高光谱图像的宽度,m表示高光谱图像的长度,B表示高光谱图像的波段数,即高光谱图像的高度,l×m×B表示三维,{Mt,t=1,2,..,X}表示若干个Gabor特征块集合,Mt表示第t个Gabor特征块,其中t表示个数,X为正整数。
3.如权利要求2所述的基于超像素级信息融合的高光谱图像的分类方法,其特征在于:
所述超像素特征数据计算步骤包括:{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt分别进行均值计算,得到n×B维的超像素特征数据Nt,最终共得到若干个n×B维的超像素特征数据集合{Nt,t=1,2,..,X};
在{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt进行均值计算时,S1对应Mt进行均值计算得到第1个B维的向量,S2对应Mt进行均值计算得到第2个B维的向量,S3对应Mt进行均值计算得到第3个B维的向量,依次类推至Sn对应Mt进行均值计算得到第n个B维的向量,最终共得到n个B维的向量,即n×B维的超像素特征数据Nt;
其中,{Si,i=1,2,..,n}表示对高光谱图像进行超像素分割得到的若干个超像素集合,Si表示第i个超像素,n表示超像素的个数;B维表示所述第一维度,{Nt,t=1,2,..,X}表示若干个n×B维的超像素特征数据集合,Nt表示第t个超像素特征数据;
其中,所述每一个超像素Si均包含若干个像素;
所述超像素空间坐标计算步骤包括:将所述每一个超像素Si分别与所述高光谱图像R进行坐标均值计算,得到1个n×2维的空间坐标数据集合C;其中,C表示空间坐标数据集合,
所述数据降维步骤包括:每一个超像素特征数据Nt分别与空间坐标数据集合C使用空谱结合的薛定谔特征映射的方法进行特征降维,从B维降到K维,得到K维的超像素特征数据Dt,最终共得到若干个K维的超像素特征数据集合{Dt,t=1,2,..,X};
其中,{Dt,t=1,2,..,X}表示维度为K的超像素特征数据集合,K表示所述第二维度,Dt表示第t个K维的超像素特征数据;所述每一个超像素特征数据Dt所对应的空间坐标数据集合仍是C。
4.如权利要求3所述的基于超像素级信息融合的高光谱图像的分类方法,其特征在于:
所述高光谱特征数据生成步骤包括:每一个超像素特征数据Dt分别和空间坐标数据集合C使用自然邻点插值法进行重构,对应原高光谱图像补足所有空间坐标的像素值,得到三维高光谱特征数据Gt,最终共得到若干个三维高光谱特征数据集合{Gt,t=1,2,..,X};
其中,{Gt,t=1,2,..,X}表示若干个三维的高光谱特征数据集合,l×m×K表示三维,Gt表示第t个三维的高光谱特征数据;
所述多任务支持向量机分类包括:将每一个高光谱特征数据Gt划分为训练数据G1t和测试数据G2t,g表示一个原始高光谱图像测试数据样本,其中g∈R,gt表示Gt当中与g同一位置坐标的像素特征数据,gt维度为K,{gt∈G2t,t=1,2,..,X}表示X个维度为K的像素特征数据集合,那么对g的分类过程如下:
(1)对训练数据{G1t,t=1,2,..,X}使用概率输出的支持向量机方法进行模型训练,得到概率输出模型{Modelt,t=1,2,..,X};
(2)使用概率输出模型Modelt对数据gt进行类别概率输出,得到gt属于每一个类别的概率{Pt(i),i=1,2,..,C},其中C为类别总数;
(3)高光谱图像样本g的类别预测公式为:
5.一种基于超像素级信息融合的高光谱图像的分类系统,其特征在于,所述系统包括:
滤波器生成模块:用于生成若干个二维的Gabor滤波器;
Gabor特征提取模块:用于将每一个所述Gabor滤波器分别与高光谱图像中的每一个波段进行卷积运算,并对卷积运算结果进行取幅值运算,得到若干个Gabor特征块;
高光谱图像超像素分割模块:用于对所述高光谱图像进行超像素分割,得到若干个超像素;
超像素特征数据计算模块:用于将每一个所述超像素与每一个所述Gabor特征块分别进行均值计算,得到若干个第一维度的超像素特征数据;
超像素空间坐标计算模块:用于将每一个所述超像素与所述高光谱图像分别进行坐标均值计算,得到1个空间坐标数据集合;
数据降维模块:用于将每一个所述第一维度的超像素特征数据分别与所述空间坐标数据集合使用空谱结合的薛定谔特征映射方法进行特征降维,从第一维度降到第二维度,得到若干个第二维度的超像素特征数据;
高光谱特征数据生成模块:用于将每一个所述第二维度的超像素特征数据分别与所述空间坐标数据集合使用自然邻点插值法进行重构,得到若干个三维高光谱特征数据;
多任务支持向量机分类模块:用于对所述每一个三维高光谱特征数据分别进行多任务的支持向量机分类。
6.如权利要求5所述的基于超像素级信息融合的高光谱图像的分类系统,其特征在于,所述Gabor特征提取模块具体用于:
将每一个所述二维Gabor滤波器与所述高光谱图像的每一波段分别进行卷积运算,并对卷积运算结果根据下述公式进行取幅值运算,得到若干个Gabor特征块:
其中,表示若干个二维Gabor滤波器集合,表示第t个Gabor滤波器,(x,y)表示在二维平面进行卷积运算时对应的二元坐标变量,R表示所述高光谱图像,其中λ表示所述高光谱图像的每一波段,l表示高光谱图像的宽度,m表示高光谱图像的长度,B表示高光谱图像的波段数,即高光谱图像的高度,l×m×B表示三维,{Mt,t=1,2,..,X}表示若干个Gabor特征块集合,Mt表示第t个Gabor特征块,其中t表示个数,X为正整数。
7.如权利要求6所述的基于超像素级信息融合的高光谱图像的分类系统,其特征在于:
所述超像素特征数据计算模块具体用于:{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt分别进行均值计算,得到n×B维的超像素特征数据Nt,最终共得到若干个n×B维的超像素特征数据集合{Nt,t=1,2,..,X};
在{Si,i=1,2,..,n}集合对应每一个Gabor特征块Mt进行均值计算时,S1对应Mt进行均值计算得到第1个B维的向量,S2对应Mt进行均值计算得到第2个B维的向量,S3对应Mt进行均值计算得到第3个B维的向量,依次类推至Sn对应Mt进行均值计算得到第n个B维的向量,最终共得到n个B维的向量,即n×B维的超像素特征数据Nt;
其中,{Si,i=1,2,..,n}表示对高光谱图像进行超像素分割得到的若干个超像素集合,Si表示第i个超像素,n表示超像素的个数;B维表示所述第一维度,{Nt,t=1,2,..,X}表示若干个n×B维的超像素特征数据集合,Nt表示第t个超像素特征数据;
其中,所述每一个超像素Si均包含若干个像素;
所述超像素空间坐标计算模块具体用于:将所述每一个超像素Si分别与所述高光谱图像R进行坐标均值计算,得到1个n×2维的空间坐标数据集合C;其中,C表示空间坐标数据集合,
所述数据降维模块具体用于:每一个超像素特征数据Nt分别与空间坐标数据集合C使用空谱结合的薛定谔特征映射的方法进行特征降维,从B维降到K维,得到K维的超像素特征数据Dt,最终共得到若干个K维的超像素特征数据集合{Dt,t=1,2,..,X};
其中,{Dt,t=1,2,..,X}表示维度为K的超像素特征数据集合,K表示所述第二维度,Dt表示第t个K维的超像素特征数据;所述每一个超像素特征数据Dt所对应的空间坐标数据集合仍是C。
8.如权利要求7所述的基于超像素级信息融合的高光谱图像的分类系统,其特征在于:
所述高光谱特征数据生成模块具体用于:每一个超像素特征数据Dt分别和空间坐标数据集合C使用自然邻点插值法进行重构,对应原高光谱图像补足所有空间坐标的像素值,得到三维高光谱特征数据Gt,最终共得到若干个三维高光谱特征数据集合{Gt,t=1,2,..,X};
其中,{Gt,t=1,2,..,X}表示若干个三维的高光谱特征数据集合,l×m×K表示三维,Gt表示第t个三维的高光谱特征数据;
所述多任务支持向量机分类模块具体用于:将每一个高光谱特征数据Gt划分为训练数据G1t和测试数据G2t,g表示一个原始高光谱图像测试数据样本,其中g∈R,gt表示Gt当中与g同一位置坐标的像素特征数据,gt维度为K,{gt∈G2t,t=1,2,..,X}表示X个维度为K的像素特征数据集合,那么对g的分类过程如下:
(1)对训练数据{G1t,t=1,2,..,X}使用概率输出的支持向量机方法进行模型训练,得到概率输出模型{Modelt,t=1,2,..,X};
(2)使用概率输出模型Modelt对数据gt进行类别概率输出,得到gt属于每一个类别的概率{Pt(i),i=1,2,..,C},其中C为类别总数;
(3)高光谱图像样本g的类别预测公式为:
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610810465.2A CN106469316B (zh) | 2016-09-07 | 2016-09-07 | 基于超像素级信息融合的高光谱图像的分类方法及系统 |
PCT/CN2016/104661 WO2018045626A1 (zh) | 2016-09-07 | 2016-11-04 | 基于超像素级信息融合的高光谱图像的分类方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610810465.2A CN106469316B (zh) | 2016-09-07 | 2016-09-07 | 基于超像素级信息融合的高光谱图像的分类方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106469316A true CN106469316A (zh) | 2017-03-01 |
CN106469316B CN106469316B (zh) | 2020-02-21 |
Family
ID=58230161
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610810465.2A Active CN106469316B (zh) | 2016-09-07 | 2016-09-07 | 基于超像素级信息融合的高光谱图像的分类方法及系统 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN106469316B (zh) |
WO (1) | WO2018045626A1 (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107194936A (zh) * | 2017-05-24 | 2017-09-22 | 哈尔滨工业大学 | 基于超像素联合稀疏表示的高光谱图像目标检测方法 |
CN107633216A (zh) * | 2017-09-07 | 2018-01-26 | 深圳大学 | 高光谱遥感图像的特征提取方法及装置 |
CN108345892A (zh) * | 2018-01-03 | 2018-07-31 | 深圳大学 | 一种立体图像显著性的检测方法、装置、设备及存储介质 |
CN108664986A (zh) * | 2018-01-16 | 2018-10-16 | 北京工商大学 | 基于lp范数正则化的多任务学习图像分类方法及系统 |
CN109101977A (zh) * | 2018-06-27 | 2018-12-28 | 深圳大学 | 一种基于无人机的数据处理的方法及装置 |
CN109472199A (zh) * | 2018-09-29 | 2019-03-15 | 深圳大学 | 一种图像融合分类的方法及装置 |
CN109492593A (zh) * | 2018-11-16 | 2019-03-19 | 西安电子科技大学 | 基于主成分分析网络和空间坐标的高光谱图像分类方法 |
CN109711269A (zh) * | 2018-12-03 | 2019-05-03 | 辽宁工程技术大学 | 一种基于3d卷积谱空特征融合的高光谱图像分类算法 |
CN110084809A (zh) * | 2019-05-06 | 2019-08-02 | 成都医云科技有限公司 | 糖尿病视网膜疾病数据处理方法、装置及电子设备 |
CN110189328A (zh) * | 2019-06-11 | 2019-08-30 | 北华航天工业学院 | 一种卫星遥感图像处理系统及其处理方法 |
CN111079544A (zh) * | 2019-11-21 | 2020-04-28 | 湖南理工学院 | 一种基于加权联合最近邻的多任务稀疏表示检测方法 |
CN111695636A (zh) * | 2020-06-15 | 2020-09-22 | 北京师范大学 | 一种基于图神经网络的高光谱图像分类方法 |
WO2021022712A1 (zh) * | 2019-08-02 | 2021-02-11 | 平安科技(深圳)有限公司 | 图像识别模型训练方法、装置以及计算机设备 |
CN115761518A (zh) * | 2023-01-10 | 2023-03-07 | 云南瀚哲科技有限公司 | 一种基于遥感图像数据的作物分类方法 |
CN117315381A (zh) * | 2023-11-30 | 2023-12-29 | 昆明理工大学 | 一种基于二阶有偏随机游走的高光谱图像分类方法 |
Families Citing this family (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108985357B (zh) * | 2018-06-29 | 2023-10-20 | 湖南理工学院 | 基于图像特征的集合经验模态分解的高光谱图像分类方法 |
CN110879953A (zh) * | 2018-09-06 | 2020-03-13 | 华南农业大学 | 植物类别识别方法及系统 |
CN109598284A (zh) * | 2018-10-23 | 2019-04-09 | 广东交通职业技术学院 | 一种基于大间隔分布和空间特征的高光谱图像分类方法 |
CN110008844B (zh) * | 2019-03-12 | 2023-07-21 | 华南理工大学 | 一种融合slic算法的kcf长期手势跟踪方法 |
CN110188821B (zh) * | 2019-05-30 | 2023-04-07 | 安徽机电职业技术学院 | 多参数边缘保持滤波和多特征学习的高光谱图像分类方法 |
CN110222773B (zh) * | 2019-06-10 | 2023-03-24 | 西北工业大学 | 基于不对称分解卷积网络的高光谱图像小样本分类方法 |
CN110490270B (zh) * | 2019-08-27 | 2022-11-04 | 大连海事大学 | 一种基于空间信息自适应处理的高光谱图像分类方法 |
CN110516614A (zh) * | 2019-08-29 | 2019-11-29 | 大连海事大学 | 一种基于子空间投影的高光谱图像分类方法 |
CN110516754B (zh) * | 2019-08-30 | 2022-11-01 | 大连海事大学 | 一种基于多尺度超像元分割的高光谱图像分类方法 |
CN110807387A (zh) * | 2019-10-25 | 2020-02-18 | 山东师范大学 | 一种基于高光谱图像特征的对象分类方法及系统 |
CN110852371B (zh) * | 2019-11-06 | 2023-04-18 | 辽宁工程技术大学 | 基于累积变异比的集成超限学习机高光谱图像分类方法 |
CN111079797B (zh) * | 2019-11-25 | 2022-02-25 | 贝壳技术有限公司 | 一种图像分类的方法、装置和存储介质 |
CN111144463B (zh) * | 2019-12-17 | 2024-02-02 | 中国地质大学(武汉) | 一种基于残差子空间聚类网络的高光谱图像聚类方法 |
CN111274869B (zh) * | 2020-01-07 | 2023-04-07 | 中国地质大学(武汉) | 基于并行注意力机制残差网进行高光谱图像分类的方法 |
CN111242056B (zh) * | 2020-01-16 | 2023-03-10 | 陕西师范大学 | 一种最大-最小距离嵌入的无监督高光谱图像分类方法 |
CN111401428A (zh) * | 2020-03-12 | 2020-07-10 | Oppo广东移动通信有限公司 | 一种图像分类方法、装置、电子设备和存储介质 |
CN111460966B (zh) * | 2020-03-27 | 2024-02-02 | 中国地质大学(武汉) | 基于度量学习和近邻增强的高光谱遥感图像分类方法 |
CN111476287A (zh) * | 2020-04-02 | 2020-07-31 | 中国人民解放军战略支援部队信息工程大学 | 一种高光谱影像小样本分类方法及装置 |
CN111639587B (zh) * | 2020-05-27 | 2023-03-24 | 西安电子科技大学 | 基于多尺度谱空卷积神经网络的高光谱图像分类方法 |
CN111723731B (zh) * | 2020-06-18 | 2023-09-29 | 西安电子科技大学 | 基于空谱卷积核的高光谱图像分类方法、存储介质及设备 |
CN111860654B (zh) * | 2020-07-22 | 2024-02-02 | 河南大学 | 一种基于循环神经网络的高光谱图像分类方法 |
CN111967182B (zh) * | 2020-07-24 | 2024-04-02 | 天津大学 | 一种用于光谱分析的基于混合标记的高光谱建模方法 |
CN111967516B (zh) * | 2020-08-14 | 2024-02-06 | 西安电子科技大学 | 一种逐像素分类方法、存储介质及分类设备 |
CN112052758B (zh) * | 2020-08-25 | 2023-05-23 | 西安电子科技大学 | 基于注意力机制和循环神经网络的高光谱图像分类方法 |
CN112101271A (zh) * | 2020-09-23 | 2020-12-18 | 台州学院 | 一种高光谱遥感影像分类方法及装置 |
CN112819959B (zh) * | 2021-01-22 | 2022-03-04 | 哈尔滨工业大学 | 高光谱图像和激光雷达数据本征高光谱点云生成方法 |
CN112837293B (zh) * | 2021-02-05 | 2023-02-14 | 中国科学院西安光学精密机械研究所 | 基于高斯函数典型关联分析的高光谱图像变化检测方法 |
CN112884040B (zh) * | 2021-02-19 | 2024-04-30 | 北京小米松果电子有限公司 | 训练样本数据的优化方法、系统、存储介质及电子设备 |
CN113065518A (zh) * | 2021-04-25 | 2021-07-02 | 中国电子科技集团公司第二十八研究所 | 基于lbp特征的超像素空谱多核高光谱图像分类方法 |
CN113205143A (zh) * | 2021-05-10 | 2021-08-03 | 郑州轻工业大学 | 耦合空谱特征的多尺度超像素高光谱遥感图像分类方法 |
CN113139513B (zh) * | 2021-05-14 | 2023-11-10 | 辽宁工程技术大学 | 基于超像素轮廓和改进pso-elm的空谱主动学习高光谱分类方法 |
CN113255698B (zh) * | 2021-06-03 | 2022-11-08 | 青岛星科瑞升信息科技有限公司 | 用于高光谱影像空间特征提取的超像素级自适应ssa方法 |
CN113343900A (zh) * | 2021-06-28 | 2021-09-03 | 中国电子科技集团公司第二十八研究所 | 基于cnn与超像素结合的组合核遥感影像目标检测方法 |
CN113723255B (zh) * | 2021-08-24 | 2023-09-01 | 中国地质大学(武汉) | 一种高光谱影像分类方法和存储介质 |
CN113963207B (zh) * | 2021-10-21 | 2024-03-29 | 江南大学 | 基于空谱信息特征引导融合网络的高光谱图像分类方法 |
CN113989525B (zh) * | 2021-12-24 | 2022-03-29 | 湖南大学 | 自适应随机块卷积核网络的高光谱中药材鉴别方法 |
CN114187479A (zh) * | 2021-12-28 | 2022-03-15 | 河南大学 | 一种基于空谱特征联合的高光谱图像分类方法 |
CN114332534B (zh) * | 2021-12-29 | 2024-03-29 | 山东省科学院海洋仪器仪表研究所 | 一种高光谱图像小样本分类方法 |
CN114663790B (zh) * | 2022-05-24 | 2022-08-26 | 济宁德信测绘有限公司 | 一种智能遥感测绘方法及系统 |
CN114882291B (zh) * | 2022-05-31 | 2023-06-06 | 南京林业大学 | 基于高光谱图像像素块机器学习的籽棉地膜识别分类方法 |
CN115795225B (zh) * | 2022-12-09 | 2024-01-23 | 四川威斯派克科技有限公司 | 一种近红外光谱校正集的筛选方法及装置 |
CN116030352B (zh) * | 2023-03-29 | 2023-07-25 | 山东锋士信息技术有限公司 | 融合多尺度分割和超像素分割的长时序土地利用分类方法 |
CN116486160B (zh) * | 2023-04-25 | 2023-12-19 | 北京卫星信息工程研究所 | 基于光谱重建的高光谱遥感图像分类方法、设备及介质 |
CN116522690B (zh) * | 2023-07-03 | 2023-09-19 | 航天宏图信息技术股份有限公司 | 新一代海洋卫星水色水温扫描仪科学数据仿真方法及装置 |
CN116758361B (zh) * | 2023-08-22 | 2023-10-27 | 中国铁路设计集团有限公司 | 基于空间与光谱联合特征的工程岩组遥感分类方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6952499B1 (en) * | 1998-10-19 | 2005-10-04 | The United States Of America As Represented By The Secretary Of The Air Force | Aspire (autonomous spatial pattern identification and recognition algorithm) |
US20060251324A1 (en) * | 2004-09-20 | 2006-11-09 | Bachmann Charles M | Method for image data processing |
CN102819568A (zh) * | 2012-07-18 | 2012-12-12 | 哈尔滨工程大学 | 基于地形采样点位置的海底地形数据建立方法 |
CN104036289A (zh) * | 2014-06-05 | 2014-09-10 | 哈尔滨工程大学 | 一种基于空间-光谱特征和稀疏表达的高光谱图像分类方法 |
CN104809471A (zh) * | 2015-04-27 | 2015-07-29 | 哈尔滨工程大学 | 一种基于空间光谱信息的高光谱图像残差融合分类方法 |
CN105069478A (zh) * | 2015-08-19 | 2015-11-18 | 西安电子科技大学 | 基于超像素张量稀疏编码的高光谱遥感地物分类方法 |
CN105512661A (zh) * | 2015-11-25 | 2016-04-20 | 中国人民解放军信息工程大学 | 一种基于多模态特征融合的遥感影像分类方法 |
CN105654098A (zh) * | 2016-03-23 | 2016-06-08 | 深圳大学 | 一种高光谱遥感图像的分类方法及其系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8374404B2 (en) * | 2009-02-13 | 2013-02-12 | Raytheon Company | Iris recognition using hyper-spectral signatures |
CN104392454B (zh) * | 2014-12-03 | 2017-07-07 | 复旦大学 | 高光谱遥感图像空谱结合分类框架下的基于地物类别隶属度评分的合并方法 |
CN105184314B (zh) * | 2015-08-25 | 2018-05-25 | 西安电子科技大学 | 基于像素聚类的wrapper式高光谱波段选择方法 |
CN105740799B (zh) * | 2016-01-27 | 2018-02-16 | 深圳大学 | 基于三维Gabor特征选择的高光谱遥感图像分类方法及系统 |
-
2016
- 2016-09-07 CN CN201610810465.2A patent/CN106469316B/zh active Active
- 2016-11-04 WO PCT/CN2016/104661 patent/WO2018045626A1/zh active Application Filing
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6952499B1 (en) * | 1998-10-19 | 2005-10-04 | The United States Of America As Represented By The Secretary Of The Air Force | Aspire (autonomous spatial pattern identification and recognition algorithm) |
US20060251324A1 (en) * | 2004-09-20 | 2006-11-09 | Bachmann Charles M | Method for image data processing |
CN102819568A (zh) * | 2012-07-18 | 2012-12-12 | 哈尔滨工程大学 | 基于地形采样点位置的海底地形数据建立方法 |
CN104036289A (zh) * | 2014-06-05 | 2014-09-10 | 哈尔滨工程大学 | 一种基于空间-光谱特征和稀疏表达的高光谱图像分类方法 |
CN104809471A (zh) * | 2015-04-27 | 2015-07-29 | 哈尔滨工程大学 | 一种基于空间光谱信息的高光谱图像残差融合分类方法 |
CN105069478A (zh) * | 2015-08-19 | 2015-11-18 | 西安电子科技大学 | 基于超像素张量稀疏编码的高光谱遥感地物分类方法 |
CN105512661A (zh) * | 2015-11-25 | 2016-04-20 | 中国人民解放军信息工程大学 | 一种基于多模态特征融合的遥感影像分类方法 |
CN105654098A (zh) * | 2016-03-23 | 2016-06-08 | 深圳大学 | 一种高光谱遥感图像的分类方法及其系统 |
Non-Patent Citations (3)
Title |
---|
LEYUAN FANG等: ""Classification of Hyperspectral Images by Exploiting Spectral-Spatial Information of Superpixel via Multiple kernels"", 《IEEE TRANSACTION ON GEOSCIENCE AND REMOTE SENSING》 * |
SHUZHEN ZHANG等: ""Spectral-spatial classification of hyperspectral images via multiscale superpixels based sparse representation"", 《2016 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM》 * |
冯逍等: ""三维Gabor滤波器与支持向量机的高光谱遥感图像分类"", 《光谱学与光谱分析》 * |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107194936A (zh) * | 2017-05-24 | 2017-09-22 | 哈尔滨工业大学 | 基于超像素联合稀疏表示的高光谱图像目标检测方法 |
CN107194936B (zh) * | 2017-05-24 | 2021-02-02 | 哈尔滨工业大学 | 基于超像素联合稀疏表示的高光谱图像目标检测方法 |
CN107633216A (zh) * | 2017-09-07 | 2018-01-26 | 深圳大学 | 高光谱遥感图像的特征提取方法及装置 |
CN107633216B (zh) * | 2017-09-07 | 2021-02-23 | 深圳大学 | 高光谱遥感图像的三维表面空谱联合特征编码方法及装置 |
CN108345892A (zh) * | 2018-01-03 | 2018-07-31 | 深圳大学 | 一种立体图像显著性的检测方法、装置、设备及存储介质 |
CN108664986B (zh) * | 2018-01-16 | 2020-09-04 | 北京工商大学 | 基于lp范数正则化的多任务学习图像分类方法及系统 |
CN108664986A (zh) * | 2018-01-16 | 2018-10-16 | 北京工商大学 | 基于lp范数正则化的多任务学习图像分类方法及系统 |
CN109101977A (zh) * | 2018-06-27 | 2018-12-28 | 深圳大学 | 一种基于无人机的数据处理的方法及装置 |
CN109101977B (zh) * | 2018-06-27 | 2022-09-20 | 深圳大学 | 一种基于无人机的数据处理的方法及装置 |
CN109472199A (zh) * | 2018-09-29 | 2019-03-15 | 深圳大学 | 一种图像融合分类的方法及装置 |
CN109472199B (zh) * | 2018-09-29 | 2022-02-22 | 深圳大学 | 一种图像融合分类的方法及装置 |
CN109492593A (zh) * | 2018-11-16 | 2019-03-19 | 西安电子科技大学 | 基于主成分分析网络和空间坐标的高光谱图像分类方法 |
CN109492593B (zh) * | 2018-11-16 | 2021-09-10 | 西安电子科技大学 | 基于主成分分析网络和空间坐标的高光谱图像分类方法 |
CN109711269A (zh) * | 2018-12-03 | 2019-05-03 | 辽宁工程技术大学 | 一种基于3d卷积谱空特征融合的高光谱图像分类算法 |
CN110084809A (zh) * | 2019-05-06 | 2019-08-02 | 成都医云科技有限公司 | 糖尿病视网膜疾病数据处理方法、装置及电子设备 |
CN110189328A (zh) * | 2019-06-11 | 2019-08-30 | 北华航天工业学院 | 一种卫星遥感图像处理系统及其处理方法 |
CN110189328B (zh) * | 2019-06-11 | 2021-02-23 | 北华航天工业学院 | 一种卫星遥感图像处理系统及其处理方法 |
WO2021022712A1 (zh) * | 2019-08-02 | 2021-02-11 | 平安科技(深圳)有限公司 | 图像识别模型训练方法、装置以及计算机设备 |
CN111079544A (zh) * | 2019-11-21 | 2020-04-28 | 湖南理工学院 | 一种基于加权联合最近邻的多任务稀疏表示检测方法 |
CN111695636A (zh) * | 2020-06-15 | 2020-09-22 | 北京师范大学 | 一种基于图神经网络的高光谱图像分类方法 |
CN111695636B (zh) * | 2020-06-15 | 2023-07-14 | 北京师范大学 | 一种基于图神经网络的高光谱图像分类方法 |
CN115761518A (zh) * | 2023-01-10 | 2023-03-07 | 云南瀚哲科技有限公司 | 一种基于遥感图像数据的作物分类方法 |
CN117315381A (zh) * | 2023-11-30 | 2023-12-29 | 昆明理工大学 | 一种基于二阶有偏随机游走的高光谱图像分类方法 |
CN117315381B (zh) * | 2023-11-30 | 2024-02-09 | 昆明理工大学 | 一种基于二阶有偏随机游走的高光谱图像分类方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106469316B (zh) | 2020-02-21 |
WO2018045626A1 (zh) | 2018-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106469316A (zh) | 基于超像素级信息融合的高光谱图像的分类方法及系统 | |
Zhao et al. | Superpixel-based multiple local CNN for panchromatic and multispectral image classification | |
CN106778821B (zh) | 基于slic和改进的cnn的极化sar图像分类方法 | |
CN107239751B (zh) | 基于非下采样轮廓波全卷积网络的高分辨sar图像分类方法 | |
Xu et al. | A 3D convolutional neural network method for land cover classification using LiDAR and multi-temporal Landsat imagery | |
CN110533077B (zh) | 用于高光谱图像分类的形状自适应卷积深度神经网络方法 | |
CN110084159A (zh) | 基于联合多级空谱信息cnn的高光谱图像分类方法 | |
CN106503739A (zh) | 联合光谱和纹理特征的高光谱遥感影像svm分类方法及系统 | |
CN103886336B (zh) | 基于稀疏自动编码器的极化sar图像分类方法 | |
Meher et al. | Wavelet-feature-based classifiers for multispectral remote-sensing images | |
Prasad et al. | Information fusion in the redundant-wavelet-transform domain for noise-robust hyperspectral classification | |
CN102289673A (zh) | 基于偏最小二乘法的高光谱遥感图像波段选择方法 | |
Huang et al. | Automatic extraction of urban impervious surfaces based on deep learning and multi-source remote sensing data | |
CN107944470A (zh) | 基于轮廓波fcn‑crf的sar图像分类方法 | |
CN107292258A (zh) | 基于双边加权调制与滤波的高光谱图像低秩表示聚类方法 | |
CN105976376A (zh) | 一种基于部件模型的高分辨sar图像目标检测方法 | |
CN111563577B (zh) | 基于Unet的跳跃层分频和多尺度鉴别的本征图像分解方法 | |
Zhang et al. | Hyperspectral image classification using an unsupervised neuro-fuzzy system | |
Shi et al. | F 3 Net: Fast Fourier filter network for hyperspectral image classification | |
CN113935433A (zh) | 基于深度光谱空间倒残差网络的高光谱图像分类方法 | |
CN111460966B (zh) | 基于度量学习和近邻增强的高光谱遥感图像分类方法 | |
Shivakumar et al. | Classification of landsat 8 imagery using Kohonen’s self organizing maps and learning vector quantization | |
Jenicka | Land Cover Classification of Remotely Sensed Images | |
Pan et al. | Enhanced FCN for farmland extraction from remote sensing image | |
Kumar et al. | Evaluating the close range hyperspectral data for feature identification and mapping |
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 |