CN105740884A - 基于奇异值分解和邻域空间信息的高光谱影像分类方法 - Google Patents
基于奇异值分解和邻域空间信息的高光谱影像分类方法 Download PDFInfo
- Publication number
- CN105740884A CN105740884A CN201610044303.2A CN201610044303A CN105740884A CN 105740884 A CN105740884 A CN 105740884A CN 201610044303 A CN201610044303 A CN 201610044303A CN 105740884 A CN105740884 A CN 105740884A
- Authority
- CN
- China
- Prior art keywords
- classification
- matrix
- training sample
- category
- neighborhood
- 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
Classifications
-
- 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/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Image Analysis (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
本发明公开一种基于奇异值分解和邻域空间信息的高光谱影像分类方法,步骤是:输入每一类的训练样本矩阵,并对训练样本矩阵进行奇异值分解,得到第j类训练样本对应的右奇异矩阵;对于训练样本矩阵中的数据,利用最小二乘方法计算其对应每一类别的残差,并比较该数据对应每个类别的残差大小,将其划分到最小残差对应的类别中;如此重复,得到每一个训练样本的类别,与其原有类别做对比,通过迭代比较,得到使训练样本的分类正确率最高的参数,再对测试样本矩阵中每个数据进行分类,输出分类结果矩阵。此种分类方法可提高分类精度,缩短分类时间,提升高光谱影像分类的自动化水平。
Description
技术领域
本发明属于图像处理技术领域,特别涉及一种基于奇异值分解和邻域空间信息的高光谱影像分类方法。
背景技术
大多数高光谱影像监督分类方法需要足够多的训练样本来支持分类结果,如线性判别分析、朴素贝叶斯、神经网络等。样本获取往往需要耗费大量的时间和金钱。支持向量机方法在小样本的情况下也能获得较高的分类精度,但是支持向量机需要选择适合的核函数及其相应的参数,不同的参数对分类结果往往影响较大,而通过交叉验证等手段获得合适的参数往往需要消耗大量的训练时间。
发明内容
本发明的目的,在于提供一种基于奇异值分解和邻域空间信息的高光谱影像分类方法,其可提高分类精度,缩短分类时间,提升高光谱影像分类的自动化水平。
为了达成上述目的,本发明的解决方案是:
一种基于奇异值分解和邻域空间信息的高光谱影像分类方法,包括如下步骤:
步骤1,输入每一类的m×n训练样本矩阵其中,m为训练样本个数,n为高光谱向量维数,j表示其对应类别;同时输入l×n测试样本矩阵Dl×n,l为测试样本数量;
步骤2,利用公式Am×n=Um×rSr×rVr×n对训练样本矩阵进行奇异值分解,其中,Um×r和Vr×n分别为左、右奇异矩阵,Sr×r为对角线元素为奇异值的奇异值矩阵;设定初始值r,r为整数,其取值范围为0<r≤min(m,n),根据r的取值,从而计算得到第j类训练样本的
步骤3,对于第i个训练样本矩阵中的数据ai×n,利用步骤2得到的对应每个类别的采用最小二乘方法计算其对应每一类别的残差In×n为单位矩阵;比较该数据ai×n对应每个类别的残差大小,将其划分到最小残差对应的类别中;
步骤4,重复步骤3,得到每一个训练样本的类别,与其原有类别做对比,计算当前r对应得到训练样本的分类正确率,保留目前使正确率达到最大的rmax及各类然后将r加1,重复步骤2-4,在r的取值范围内,得到训练样本的分类正确率最高的rmax及各类从而完成训练过程并转至步骤5;
步骤5,对测试样本矩阵Di×n中第i个数据di×n进行分类,1≤i≤l,其中ej表示在每个类别中的残差;
步骤6,重复步骤5,对测试样本矩阵Dl×n中的每一条数据都进行分类,输出分类结果矩阵Rl×1。
上述步骤2中,r的取值为5,6,…,20,初始值为5。
上述步骤5的详细内容是:对于测试样本矩阵Dl×n中第i个数据di×n,1≤i≤l,利用计算其在每个类别中利用表示的残差,式中为步骤4得到的将其划分到残差最小的一类,k=3,表示方形邻域窗口的大小;表示每一行由di×n邻域像元组成的矩阵,表示邻域矩阵对应的系数向量,每一个系数为邻域像元的权重,(0,0)为中心像元,(i,j)表示邻域像元在方形邻域窗口中的位置,-k≤i,j≤k。
采用上述方案后,本发明在小样本的情况下,通过奇异值分解提取样本的主要信息,引入测试像元的邻域像元,依据最小二乘原则进行分类,提高分类精度的同时,缩短了分类的时间,通过自动化的参数设置,提升了高光谱影像分类的自动化水平,便于高光谱影像大规模处理应用。
具体实施方式
以下将结合具体实施例,对本发明的技术方案及有益效果进行详细说明。
本发明提供一种基于奇异值分解和邻域空间信息的高光谱影像分类方法,包括如下步骤:
步骤1,输入数据:输入每一类的m×n训练样本矩阵其中,m为训练样本个数,n为高光谱向量维数,j表示其对应类别(如水域、道路、树木等);同时输入l×n测试样本矩阵Dl×n,l为测试样本数量。
以下将进入迭代循环,找到最合适的r值,以得到最佳的分类方案。
步骤2,利用公式Am×n=Um×rSr×rVr×n对训练样本矩阵进行奇异值分解,其中,Um×r和Vr×n分别为左、右奇异矩阵,Sr×r为对角线元素为奇异值的奇异值矩阵;设定初始值r,r为整数,其取值范围为0<r≤min(m,n),根据r的取值,从而计算得到第j类训练样本的通常来说,由于前几个奇异值及其对应的左、右奇异向量就能表示绝大多数被分解矩阵的信息,因此在本实施例中,r取值5,6,…,20,初始值为5;
步骤3,对于第i个训练样本矩阵中的数据ai×n,利用步骤2得到的对应每个类别的采用最小二乘方法计算其对应每一类别的残差In×n为单位矩阵。比较该数据ai×n对应每个类别的残差大小,将其划分到最小残差对应的类别中;
步骤4,重复步骤3,得到每一个训练样本的类别,与其原有类别做对比,计算当前r对应得到训练样本的分类正确率,保留目前使正确率达到最大的rmax及各类然后将r加1,重复步骤2-4,在r的取值范围内,得到训练样本的分类正确率最高的rmax及各类从而完成训练过程并转至步骤5;
以下将利用得到的对测试样本矩阵每一条数据进行分类。
步骤5,对于测试样本矩阵Di×n中第i个数据di×n,1≤i≤l,测试样本及其邻域像元同时参与计算,利用计算其在每个类别中利用表示的残差,其中,k=3,表示方形邻域窗口的大小;式中为步骤4得到的将其划分到残差最小的一类,表示每一行由di×n邻域像元组成的矩阵,表示邻域矩阵对应的系数向量,每一个系数为邻域像元的权重,(0,0)为中心像元,(i,j)表示邻域像元在方形邻域窗口中的位置,-k≤i,j≤k;
步骤6,重复步骤5,对测试样本矩阵Dl×n中的每一条数据都进行分类,输出分类结果矩阵Rl×1。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (3)
1.一种基于奇异值分解和邻域空间信息的高光谱影像分类方法,其特征在于包括如下步骤:
步骤1,输入每一类的m×n训练样本矩阵其中,m为训练样本个数,n为高光谱向量维数,j表示其对应类别;同时输入l×n测试样本矩阵Dl×n,l为测试样本数量;
步骤2,利用公式Am×n=Um×rSr×rVr×n对训练样本矩阵进行奇异值分解,其中,Um×r和Vr×n分别为左、右奇异矩阵,Sr×r为对角线元素为奇异值的奇异值矩阵;设定初始值r,r为整数,其取值范围为0<r≤min(m,n),根据r的取值,从而计算得到第j类训练样本的
步骤3,对于第i个训练样本矩阵中的数据ai×n,利用步骤2得到的对应每个类别的采用最小二乘方法计算其对应每一类别的残差In×n为单位矩阵;比较该数据ai×n对应每个类别的残差大小,将其划分到最小残差对应的类别中;
步骤4,重复步骤3,得到每一个训练样本的类别,与其原有类别做对比,计算当前r对应得到训练样本的分类正确率,保留目前使正确率达到最大的rmax及各类然后将r加1,重复步骤2-4,在r的取值范围内,得到训练样本的分类正确率最高的rmax及各类从而完成训练过程并转至步骤5;
步骤5,对测试样本矩阵Dl×n中第i个数据di×n进行分类,1≤i≤l,其中ej表示在每个类别中的残差;
步骤6,重复步骤5,对测试样本矩阵Dl×n中的每一条数据都进行分类,输出分类结果矩阵Rl×1。
2.如权利要求1所述的基于奇异值分解和邻域空间信息的高光谱影像分类方法,其特征在于:所述步骤2中,r的取值为5,6,…,20,初始值为5。
3.如权利要求1所述的基于奇异值分解和邻域空间信息的高光谱影像分类方法,其特征在于:所述步骤5的详细内容是:对于测试样本矩阵Dl×r中第i个数据di×n,1≤i≤l,利用计算其在每个类别中利用表示的残差,式中为步骤4得到的将其划分到残差最小的一类,k=3,表示方形邻域窗口的大小;表示每一行由di×n邻域像元组成的矩阵,表示邻域矩阵对应的系数向量,每一个系数为邻域像元的权重,(0,0)为中心像元,(i,j)表示邻域像元在方形邻域窗口中的位置,-k≤i,j≤k。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610044303.2A CN105740884B (zh) | 2016-01-22 | 2016-01-22 | 基于奇异值分解和邻域空间信息的高光谱影像分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610044303.2A CN105740884B (zh) | 2016-01-22 | 2016-01-22 | 基于奇异值分解和邻域空间信息的高光谱影像分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105740884A true CN105740884A (zh) | 2016-07-06 |
CN105740884B CN105740884B (zh) | 2019-06-07 |
Family
ID=56246460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610044303.2A Expired - Fee Related CN105740884B (zh) | 2016-01-22 | 2016-01-22 | 基于奇异值分解和邻域空间信息的高光谱影像分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105740884B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108921088A (zh) * | 2018-06-29 | 2018-11-30 | 佛山市顺德区中山大学研究院 | 一种基于判别式目标方程的人脸识别方法 |
CN109271874A (zh) * | 2018-08-23 | 2019-01-25 | 广东工业大学 | 一种融合空间和光谱信息的高光谱图像特征提取方法 |
CN110738246A (zh) * | 2019-09-29 | 2020-01-31 | 深圳和而泰家居在线网络科技有限公司 | 产品分类方法、装置、计算设备及计算机存储介质 |
CN111191617A (zh) * | 2020-01-02 | 2020-05-22 | 武汉大学 | 一种基于分层结构的遥感场景分类方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6415233B1 (en) * | 1999-03-04 | 2002-07-02 | Sandia Corporation | Classical least squares multivariate spectral analysis |
CN104281855A (zh) * | 2014-09-30 | 2015-01-14 | 西安电子科技大学 | 基于多任务低秩的高光谱图像分类方法 |
CN104408478A (zh) * | 2014-11-14 | 2015-03-11 | 西安电子科技大学 | 一种基于分层稀疏判别特征学习的高光谱图像分类方法 |
CN104778482A (zh) * | 2015-05-05 | 2015-07-15 | 西安电子科技大学 | 基于张量半监督标度切维数约减的高光谱图像分类方法 |
-
2016
- 2016-01-22 CN CN201610044303.2A patent/CN105740884B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6415233B1 (en) * | 1999-03-04 | 2002-07-02 | Sandia Corporation | Classical least squares multivariate spectral analysis |
CN104281855A (zh) * | 2014-09-30 | 2015-01-14 | 西安电子科技大学 | 基于多任务低秩的高光谱图像分类方法 |
CN104408478A (zh) * | 2014-11-14 | 2015-03-11 | 西安电子科技大学 | 一种基于分层稀疏判别特征学习的高光谱图像分类方法 |
CN104778482A (zh) * | 2015-05-05 | 2015-07-15 | 西安电子科技大学 | 基于张量半监督标度切维数约减的高光谱图像分类方法 |
Non-Patent Citations (2)
Title |
---|
JIANI ZHANG ET AL.: "Randomized SVD Methods in Hyperspectral Imaging", 《JOURNAL OF ELECTRICAL AND COMPUTER ENGINEERING》 * |
何同弟: "高光谱图像的分类技术研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108921088A (zh) * | 2018-06-29 | 2018-11-30 | 佛山市顺德区中山大学研究院 | 一种基于判别式目标方程的人脸识别方法 |
CN108921088B (zh) * | 2018-06-29 | 2022-03-04 | 佛山市顺德区中山大学研究院 | 一种基于判别式目标方程的人脸识别方法 |
CN109271874A (zh) * | 2018-08-23 | 2019-01-25 | 广东工业大学 | 一种融合空间和光谱信息的高光谱图像特征提取方法 |
CN109271874B (zh) * | 2018-08-23 | 2022-02-11 | 广东工业大学 | 一种融合空间和光谱信息的高光谱图像特征提取方法 |
CN110738246A (zh) * | 2019-09-29 | 2020-01-31 | 深圳和而泰家居在线网络科技有限公司 | 产品分类方法、装置、计算设备及计算机存储介质 |
CN110738246B (zh) * | 2019-09-29 | 2024-01-19 | 广东中创智家科学研究有限公司 | 产品分类方法、装置、计算设备及计算机存储介质 |
CN111191617A (zh) * | 2020-01-02 | 2020-05-22 | 武汉大学 | 一种基于分层结构的遥感场景分类方法 |
CN111191617B (zh) * | 2020-01-02 | 2022-02-01 | 武汉大学 | 一种基于分层结构的遥感场景分类方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105740884B (zh) | 2019-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105243398B (zh) | 基于线性判别分析准则的改进卷积神经网络性能的方法 | |
CN105787516B (zh) | 一种基于空谱局部性低秩超图学习的高光谱图像分类方法 | |
CN105740884A (zh) | 基于奇异值分解和邻域空间信息的高光谱影像分类方法 | |
CN112560967B (zh) | 一种多源遥感图像分类方法、存储介质及计算设备 | |
CN111428818A (zh) | 基于神经通路激活状态的深度学习模型测试方法与装置 | |
CN110287983A (zh) | 基于最大相关熵深度神经网络单分类器异常检测方法 | |
CN108198147A (zh) | 一种基于判别字典学习的多源图像融合去噪的方法 | |
CN110675368A (zh) | 一种融合图像分割与分类的细胞图像语义分割方法 | |
CN104573699B (zh) | 基于中等场强磁共振解剖成像的实蝇识别方法 | |
CN104809471B (zh) | 一种基于空间光谱信息的高光谱图像残差融合分类方法 | |
CN112560966B (zh) | 基于散射图卷积网络的极化sar图像分类方法、介质及设备 | |
CN113793319B (zh) | 基于类别约束字典学习模型的织物图像瑕疵检测方法及系统 | |
CN110533063A (zh) | 一种基于卫星影像和gmdh神经网络的云量计算方法及装置 | |
Agarwal et al. | Assessment of optimizers impact on image recognition with convolutional neural network to adversarial datasets | |
Landa et al. | Accurate classification of fresh and charred grape seeds to the varietal level, using machine learning based classification method | |
CN115439654A (zh) | 动态约束下的弱监督农田地块精细化分割方法及系统 | |
Kar et al. | Benchmarking of deep learning algorithms for 3D instance segmentation of confocal image datasets | |
CN106908774B (zh) | 基于多尺度核稀疏保持投影的一维距离像识别方法 | |
Pang et al. | Convolutional neural network-based sub-pixel line-edged angle detection with applications in measurement | |
CN104866905A (zh) | 一种基于beta过程的非参数稀疏张量字典的学习方法 | |
CN115482227B (zh) | 机器视觉自适应成像环境调整方法 | |
CN115206455B (zh) | 基于深度神经网络的稀土元素组分含量预测方法及系统 | |
CN110866552A (zh) | 基于全卷积空间传播网络的高光谱图像分类方法 | |
CN108108758A (zh) | 面向工业大数据的多层增量特征提取方法 | |
US11062446B2 (en) | Analyzing method for perovskite structure using machine learning |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190607 Termination date: 20210122 |
|
CF01 | Termination of patent right due to non-payment of annual fee |