CN109146984B - 一种基于粒子群优化的高光谱图像稀疏分解方法 - Google Patents
一种基于粒子群优化的高光谱图像稀疏分解方法 Download PDFInfo
- Publication number
- CN109146984B CN109146984B CN201810923137.2A CN201810923137A CN109146984B CN 109146984 B CN109146984 B CN 109146984B CN 201810923137 A CN201810923137 A CN 201810923137A CN 109146984 B CN109146984 B CN 109146984B
- Authority
- CN
- China
- Prior art keywords
- particle
- turning
- image
- particles
- sparse decomposition
- 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
Links
- 239000002245 particle Substances 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 25
- 238000005457 optimization Methods 0.000 title claims abstract description 10
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 16
- 238000012545 processing Methods 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 9
- 238000004088 simulation Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 241000488988 Monotropa Species 0.000 description 1
- 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
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000002474 experimental method Methods 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
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T9/00—Image coding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Biomedical Technology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Image Processing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于粒子群优化的高光谱图像稀疏分解方法,主要解决现有的正交匹配追踪稀疏分解算法计算复杂度高的问题。其技术关键是借助粒子群优化的思想,对正交匹配追踪算法的匹配过程进行改进,利用粒子表示冗余字典中的原子,依靠粒子群的快速搜索能力,找到能够对图像进行稀疏表示的最优原子,实现图像的稀疏分解。本发明方法能够在保证重构精度的条件下,提高稀疏分解的效率,且算法不需要事先产生冗余字典,减少对存储空间的占用,满足实时性要求。
Description
技术领域
本发明涉及图像处理技术领域,尤其涉及一种基于粒子群优化的高光谱图像稀疏分解方法。
背景技术
高光谱图像不但包含了被观测目标的空间分布信息,而且图像中的每个像元都有几十个甚至上百个窄波段的丰富光谱信息,具备“图谱合一”的性质。由于高光谱图像可以把反映物质性质的光谱特征和呈现物质几何空间信息的图像信息维系在一起,因此极大地提高了人类认知客观世界的能力,在遥感、军事、农业、医学等领域都被证明有着巨大的应用价值。
因高光谱图像的空间、谱间分辨率高的特点,使其包含了非常丰度的细节信息,导致其数据量巨大,给星载传感系统的传输和存储都带来较大的困难。为了解决这个问题,研究者引入压缩感知理论,对高光谱图像进行压缩处理。压缩感知理论将采样与压缩过程结合,直接采集数据的信息特性,可对稀疏信号进行处理,降低传感器的采样和计算成本。因高光谱图像波段多,图像分辨率高,如何设计计算复杂度低的稀疏分解算法成为研究的关键。
最常用的稀疏分解算法是基于全局搜索的贪婪算法,如正交匹配追踪(Orthogonal Matching Pursuit,OMP),OMP算法由于其收敛快且精度高被大量采用。但是,匹配追踪算法在处理基于冗余字典的稀疏分解问题时,由于冗余字典中原子过多,而匹配过程需遍历字典中所有原子,计算复杂度高,计算时间在现有计算条件下令人无法忍受。综上所述,现有技术存在的主要问题在于:稀疏分解过程的计算复杂度高,无法应用于高光谱图像的稀疏分解。
发明内容
有鉴于现有技术的上述缺陷,本发明所要解决的技术问题是提供一种基于粒子群优化的高光谱图像稀疏分解方法,利用粒子群对正交匹配追踪算法的匹配过程进行优化,快速搜索到最优原子,提高高光谱图像的稀疏分解速度。
为实现上述目的,本发明提供了一种基于粒子群优化的高光谱图像稀疏分解方法,包括以下步骤:
步骤1.高光谱图像总的波段数为J,波段序号为j,对每个波段图像进行分块处理,分块大小为B,分块个数为L,图像块号表示为l,则Xj,l表示第j个波段第l个图像块;设定最优原子个数为K,设定粒子群算法的种群个数为M,最大更新代数为T;
步骤2.令j=1,l=1;
步骤4.如果j>J,则稀疏分解过程结束;否则,转入步骤5;
步骤5.设定原子个数为k=1;初始化残差r0=Xj,l,最优原子索引集合为Λ0=[];
步骤6.利用粒子群算法搜索得到最优原子的索引;
Λk=Λk-1YGbest (1)
步骤8.根据公式(2)更新残差:
步骤11.令l=l+1,转入步骤3。
上述的一种基于粒子群优化的高光谱图像稀疏分解方法,所述步骤6具体为:
步骤6.1.令粒子标号为m=1;
其中,win表示高斯函数,n=1,2,...,B2;
步骤6.6.令粒子的更新代数为t=1;
步骤6.7.令粒子标号为m=1;
步骤6.8.根据公式(6)和公式(7)更新粒子的速度和位置:
其中,ls=1,2,3,4,w为惯性权重,c1和c2是非负的常数,称为加速度因子,r1和r2是分布于[0,1]区间的随机数;
其中,win表示高斯函数,n=1,2,...,B2;
步骤6.13.令t=t+1,如果t>T,输出Gbest,并转入步骤7;否则转入步骤6.6。
本发明的有益效果是:
第一,本发明由于采用粒子的位置代表原子,种群中粒子参数向量仅为一个四维向量,其占用的内存基本可以忽略不计;与正交匹配追踪算法需要事先产生冗余字典相比,本发明方法所需的计算内存较小;第二,本发明方法利用粒子群在整个参数空间中进行搜索,粒子位置的取值是连续的,与事先产生离散的冗余字典相比,本发明方法能生成更多类型的原子,更能稀疏的表示图像信号,从而提高了高光谱图像稀疏分解后的重构精度;第三,本发明方法利用粒子群搜索最优原子,仅需搜索少量的参数空间点,再由这些参数空间点产生原子,与正交匹配追踪算法的内积运算相比,本发明方法的运算量非常小,计算效率得到提升。
以下将结合附图对本发明的构思、具体结构及产生的技术效果作进一步说明,以充分地了解本发明的目的、特征和效果。
附图说明
图1为本发明的流程图;
图2为本发明仿真实验所用的四组高光谱原始图像的第40个波段图像的示意图;
图3为本发明方法得到的重构图像与已有的OMP方法得到的重构图像对比。
具体实施方式
一种基于粒子群优化的高光谱图像稀疏分解方法,其特征在于按如下步骤进行:
步骤1.高光谱图像总的波段数为J,波段序号为j,对每个波段图像进行分块处理,分块大小为B,分块个数为L,图像块号表示为l,则Xj,l表示第j个波段第l个图像块;设定最优原子个数为K,设定粒子群算法的种群个数为M,最大更新代数为T;
步骤2.令j=1,l=1;
步骤4.如果j>J,则稀疏分解过程结束;否则,转入步骤5;
步骤5.设定原子个数为k=1;初始化残差r0=Xj,l,最优原子索引集合为Λ0=[];
步骤6.利用粒子群算法搜索得到最优原子的索引:
步骤6.1.令粒子标号为m=1;
其中,win表示高斯函数,n=1,2,...,B2;
步骤6.6.令粒子的更新代数为t=1;
步骤6.7.令粒子标号为m=1;
步骤6.8.根据公式(3)和公式(4)更新粒子的速度和位置:
其中,ls=1,2,3,4,w为惯性权重,c1和c2是非负的常数,称为加速度因子,r1和r2是分布于[0,1]区间的随机数;
其中,win表示高斯函数,n=1,2,...,B2;
步骤6.13.令t=t+1,如果t>T,输出Gbest,并转入步骤7;否则转入步骤6.6;
Λk=Λk-1YGbest (7)
步骤8.根据公式(8)更新残差:
步骤11.令l=l+1,转入步骤3。
本发明的效果可以通过如下仿真实验具体说明:
1.仿真条件:
1)仿真实验中的四组高光谱图像分别Cupprite1场景、Cuprite2场景、IndianPines场景和Pavia University场景;Cupprite1场景、Cuprite2场景由AVIRIS采集得到,实验图像大小为256×256,波段数为188;Indian Pines场景由AVIRIS采集得到,实验图像大小为128×128,波段数为200;Pavia University场景由ROSIS采集得到,实验图像大小为256×256,波段数为103;四组高光谱数据集第40个波段的原始图像如2所示;
2)仿真实验所用的编程平台为Matlab R2012b;
3)仿真实验中,采用峰值信噪比(Peak Signal-to-Noise Ratio,PSNR)指标来评价实验结果,峰值信噪比PSNR的定义为:
2.仿真内容:
采用本发明方法和现有的OMP方法对四组高光谱图像进行稀疏分解,第40个波段图像的重构结果如图3所示;
从图3所显示的重构结果可以看出,本发明得到的重构图像比OMP方法得到的重构图像更接近原始图像。
3.峰值信噪比PSNR和计算时间对比
计算现有的OMP方法和本发明方法对四组高光谱图像进行稀疏分解得到的峰值信噪比PSNR以及计算时间,结果如表1所示。
表1四组高光谱图像稀疏表示的重构峰值信噪比(单位:dB)和计算时间(单位:s)对比
从表1可以看出,本发明方法重构的高光谱图像的峰值信噪比PSNR高于OMP方法的峰值信噪比,与OMP算法相比,本发明方法的运行时间能降低一个数量级,提高了计算效率。
以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术人员无需创造性劳动就可以根据本发明的构思做出诸多修改和变化。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。
Claims (1)
1.一种基于粒子群优化的高光谱图像稀疏分解方法,其特征在于,包括以下步骤:
步骤1.高光谱图像总的波段数为J,波段序号为j,对每个波段图像进行分块处理,分块大小为B,分块个数为L,图像块号表示为l,则Xj,l表示第j个波段第l个图像块;设定最优原子个数为K,设定粒子群算法的种群个数为M,最大更新代数为T;
步骤2.令j=1,l=1;
步骤4.如果j>J,则稀疏分解过程结束;否则,转入步骤5;
步骤5.设定原子个数为k=1;初始化残差r0=Xj,l,最优原子索引集合为Λ0=[];
步骤6.利用粒子群算法搜索得到最优原子的索引;
Λk=Λk-1∪Gbest(1)
步骤8.根据公式(2)更新残差:
步骤11.令l=l+1,转入步骤3;
其中,所述步骤6具体为:
步骤6.1.令粒子标号为m=1;
其中,win表示高斯函数,n=1,2,...,B2;
步骤6.6.令粒子的更新代数为t=1;
步骤6.7.令粒子标号为m=1;
步骤6.8.根据公式(6)和公式(7)更新粒子的速度和位置:
其中,ls=1,2,3,4,w为惯性权重,c1和c2是非负的常数,称为加速度因子,r1和r2是分布于[0,1]区间的随机数;
其中,win表示高斯函数,n=1,2,...,B2;
步骤6.13.令t=t+1,如果t>T,输出Gbest,并转入步骤7;否则转入步骤6.6。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810923137.2A CN109146984B (zh) | 2018-08-14 | 2018-08-14 | 一种基于粒子群优化的高光谱图像稀疏分解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810923137.2A CN109146984B (zh) | 2018-08-14 | 2018-08-14 | 一种基于粒子群优化的高光谱图像稀疏分解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109146984A CN109146984A (zh) | 2019-01-04 |
CN109146984B true CN109146984B (zh) | 2022-11-22 |
Family
ID=64792982
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810923137.2A Active CN109146984B (zh) | 2018-08-14 | 2018-08-14 | 一种基于粒子群优化的高光谱图像稀疏分解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109146984B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101739666A (zh) * | 2009-12-15 | 2010-06-16 | 西南交通大学 | 基于一维快速哈特莱变换和匹配追踪的图像稀疏分解快速算法 |
CN102013106A (zh) * | 2010-10-27 | 2011-04-13 | 西安电子科技大学 | 基于Curvelet冗余字典的图像稀疏表示方法 |
CN102609910A (zh) * | 2012-01-13 | 2012-07-25 | 西安电子科技大学 | 基于Ridgelet冗余字典的遗传进化图像重构方法 |
WO2018072351A1 (zh) * | 2016-10-20 | 2018-04-26 | 北京工业大学 | 一种基于粒子群优化算法对支持向量机的优化方法 |
-
2018
- 2018-08-14 CN CN201810923137.2A patent/CN109146984B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101739666A (zh) * | 2009-12-15 | 2010-06-16 | 西南交通大学 | 基于一维快速哈特莱变换和匹配追踪的图像稀疏分解快速算法 |
CN102013106A (zh) * | 2010-10-27 | 2011-04-13 | 西安电子科技大学 | 基于Curvelet冗余字典的图像稀疏表示方法 |
CN102609910A (zh) * | 2012-01-13 | 2012-07-25 | 西安电子科技大学 | 基于Ridgelet冗余字典的遗传进化图像重构方法 |
WO2018072351A1 (zh) * | 2016-10-20 | 2018-04-26 | 北京工业大学 | 一种基于粒子群优化算法对支持向量机的优化方法 |
Non-Patent Citations (3)
Title |
---|
一种优化稀疏分解的雷达目标识别方法;赵东波等;《现代电子技术》;20171201;第40卷(第23期);第1-5页 * |
基于粒子群优化的图像稀疏分解算法研究;王丽 等;《计算机仿真》;20151130;第32卷(第11期);第363-366页 * |
基于粒子群优化的稀疏分解最优匹配原子搜索算法;王春光等;《国防科技大学学报》;20080430;第30卷(第02期);第83-87页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109146984A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109087367B (zh) | 一种基于粒子群优化的高光谱图像快速压缩感知重构方法 | |
CN109683161B (zh) | 一种基于深度admm网络的逆合成孔径雷达成像的方法 | |
Wang et al. | Review of image fusion based on pulse-coupled neural network | |
CN112285712B (zh) | 一种提高sar图像中靠岸船只检测精度的方法 | |
Wang et al. | NAS-guided lightweight multiscale attention fusion network for hyperspectral image classification | |
CN108460391B (zh) | 基于生成对抗网络的高光谱图像无监督特征提取方法 | |
Wang et al. | TRC‐YOLO: A real‐time detection method for lightweight targets based on mobile devices | |
CN110119805B (zh) | 基于回声状态网络分类的卷积神经网络算法 | |
CN113902622B (zh) | 基于深度先验联合注意力的光谱超分辨率方法 | |
Jifara et al. | Hyperspectral image compression based on online learning spectral features dictionary | |
CN111709882A (zh) | 基于亚像素卷积与特征分割的超分辨率融合的计算方法 | |
Qi et al. | Global–local 3-D convolutional transformer network for hyperspectral image classification | |
CN115982736B (zh) | 用于计算机网络信息的数据加密方法及系统 | |
CN114612716A (zh) | 一种基于自适应解码器的目标检测方法及装置 | |
CN115222998A (zh) | 一种图像分类方法 | |
Shi et al. | F 3 Net: Fast Fourier filter network for hyperspectral image classification | |
Farooque et al. | Swin transformer with multiscale 3D atrous convolution for hyperspectral image classification | |
CN114170154A (zh) | 基于Transformer的遥感VHR图像变化检测方法 | |
CN109146984B (zh) | 一种基于粒子群优化的高光谱图像稀疏分解方法 | |
Singh et al. | A Pre-processing framework for spectral classification of hyperspectral images | |
CN115017366B (zh) | 基于多粒度语境化和多结构保存的无监督视频哈希检索方法 | |
Qian et al. | Structured sparse model based feature selection and classification for hyperspectral imagery | |
CN113435510B (zh) | 基于序列调整网络的逆合成孔径雷达图像分类方法 | |
CN114627370A (zh) | 一种基于transformer特征融合的高光谱影像分类方法 | |
CN114528918A (zh) | 一种基于二维卷积和lstm的高光谱图像分类方法及系统 |
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 |