CN110570359B - 基于光谱和空间总变分最小限制的多层非负矩阵分解高光谱图像解混方法 - Google Patents
基于光谱和空间总变分最小限制的多层非负矩阵分解高光谱图像解混方法 Download PDFInfo
- Publication number
- CN110570359B CN110570359B CN201910320619.3A CN201910320619A CN110570359B CN 110570359 B CN110570359 B CN 110570359B CN 201910320619 A CN201910320619 A CN 201910320619A CN 110570359 B CN110570359 B CN 110570359B
- Authority
- CN
- China
- Prior art keywords
- total variation
- hyperspectral image
- matrix
- domain
- spectrum
- 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
- 239000011159 matrix material Substances 0.000 title claims abstract description 55
- 238000001228 spectrum Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 12
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 22
- 230000003595 spectral effect Effects 0.000 claims abstract description 22
- 239000013598 vector Substances 0.000 claims description 5
- 238000005457 optimization Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000009795 derivation Methods 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 abstract description 9
- 238000004088 simulation Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 description 7
- FDKCTEWMJWRPDS-UHFFFAOYSA-N dialuminum;trimagnesium;trisilicate Chemical compound [Mg+2].[Mg+2].[Mg+2].[Al+3].[Al+3].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-] FDKCTEWMJWRPDS-UHFFFAOYSA-N 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- NLYAJNPCOHFWQQ-UHFFFAOYSA-N kaolin Chemical compound O.O.O=[Al]O[Si](=O)O[Si](=O)O[Al]=O NLYAJNPCOHFWQQ-UHFFFAOYSA-N 0.000 description 2
- 229910052622 kaolinite Inorganic materials 0.000 description 2
- 229910001595 Buddingtonite Inorganic materials 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 229910052836 andradite Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011045 chalcedony Substances 0.000 description 1
- 238000000701 chemical imaging Methods 0.000 description 1
- BERDEBHAJNAUOM-UHFFFAOYSA-N copper(I) oxide Inorganic materials [Cu]O[Cu] BERDEBHAJNAUOM-UHFFFAOYSA-N 0.000 description 1
- LBJNMUFDOHXDFG-UHFFFAOYSA-N copper;hydrate Chemical compound O.[Cu].[Cu] LBJNMUFDOHXDFG-UHFFFAOYSA-N 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- YGANSGVIUGARFR-UHFFFAOYSA-N dipotassium dioxosilane oxo(oxoalumanyloxy)alumane oxygen(2-) Chemical compound [O--].[K+].[K+].O=[Si]=O.O=[Al]O[Al]=O YGANSGVIUGARFR-UHFFFAOYSA-N 0.000 description 1
- 229910052852 dumortierite Inorganic materials 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000011423 initialization method Methods 0.000 description 1
- 229910052901 montmorillonite Inorganic materials 0.000 description 1
- 229910052627 muscovite Inorganic materials 0.000 description 1
- 229910000273 nontronite Inorganic materials 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 229910052832 pyrope Inorganic materials 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10036—Multispectral image; Hyperspectral image
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了基于光谱和空间总变分最小限制的多层非负矩阵分解高光谱图像解混方法,设高光谱图像矩阵Y,端元光谱矩阵M,丰度矩阵R,随机噪声E,建立线性光谱混合模型将多层非负矩阵分解应用于线性混合模型:设计光谱域、空间域总变分函数;将光谱域和空间域总变分最小约束引入MLNMF,建立SSTV‑MLNMF目标函数;对得到的目标函数进行优化;选定实验参数,得到最终的解混结果。本发明将总变分最小的特点应用到了光谱域以及空间域,通过模拟高光谱图像和真实高光谱图像数据实验,验证了本算法的有效性,本方法比其它方法解混精度高。
Description
技术领域
本发明涉及图像信息处理技术领域,涉及一种光谱和空间总变分最小限制的 多层非负矩阵分解高光谱图像解混方法。
背景技术
由于含有丰富的波段信息,高光谱图像被广泛应用于卫星遥感、农作物观测 以及矿物勘测。然而,高光谱成像设备空间分辨率有限,并且离地观测距离较远, 使得在高光谱图像中存在混合像元,即一个像素中往往包含有多种地物类型的混 合光谱。高光谱图像解混便成为了将混合像元里所含地物类型的光谱(端元)分 辨出来并且计算地物类型所占比例(丰度)的重要技术手段。在众多的解混方法 中,非负矩阵分解因其简单、高效的解混性能受到了重度研究学者的关注。但是 由于非负矩阵分解算法的非凸性,其结果往往收敛于局部最小值,影响其最优解。 因此,根据高光谱图像的不同特性,多种不同的限制条件加入到非负矩阵分解算 法当中,如稀疏、流形、先验光谱知识等限制,以提高高光谱图像解混的精度。 传统的非负矩阵分解算法只有一层结构往往无法适应高度混合的高光谱图像,近年来多层非负矩阵分解很好的解决了这一问题。本算法就是基于多层非负矩阵分 解研究其光谱域和空间域总变分最小限制的方法。
发明内容
高光谱图像解混经过多年的发展已经取得了一定的成果,但其仍具有较大的 挑战性。
本发明的目的是为了更好解决基于多层非负矩阵分解目标函数非凸性的问 题,并结合高光谱图像光谱域以及空间域的特点,提出光谱域和空间域总变分最 小的约束条件,通过多层非负矩阵分解来提高高光谱图像解混精度,该方法称为 基于光谱和空间总变分最小限制的多层非负矩阵分解的高光谱图像解混方法 (spectral and spatial totalvariation regularized multilayer nonnegative matrix factorization),简称SSTV-MLNMF。
本发明采用的技术方案为基于光谱和空间总变分最小限制的多层非负矩阵 分解的高光谱图像解混方法,由以下步骤实现:
步骤一、设高光谱图像矩阵Y,端元光谱矩阵M,丰度矩阵R,随机噪声E, 建立线性光谱混合模型:
Y=MR+E (1)
步骤二、将多层非负矩阵分解(MLNMF)应用于线性混合模型:
Y=M1,…,MlRl+E (2)
其中M1,…,Ml为第l层;
步骤三、设计光谱域、空间域总变分函数;
步骤四、将光谱域和空间域总变分最小约束引入MLNMF,建立 SSTV-MLNMF目标函数;
步骤五、对得到的目标函数进行优化;
步骤六、选定实验参数,得到最终的解混结果。
本发明的有益效果为:
本发明提出的基于光谱和空间总变分最小限制的多层非负矩阵分解的高光 谱图像解混方法,将总变分最小的特点应用到了光谱域以及空间域,这种约束可 以在光谱域和空间域对得到的解混结果进行平滑;本发明在此约束基础之上,设 计了基于这两种约束的MLNMF目标函数,重新设计其迭代算法,最终得到其 解混结果。通过模拟高光谱图像和真实高光谱图像数据实验,验证了本算法的有 效性,通过对比SAD值和RMSE值,本发明所提出的高光谱图像解混方法比其 它方法解混精度提高5-10%。
附图说明
图1为本发明的流程图;
图2为本发明模拟数据所用UGGS不同光谱示例图;
图3(a)-图3(b)为本发明不同解混方法解混效果对比图;
图4(a)-图4(b)为本发明不同噪声强度时算法效果对比图;
图5(a)-图5(c)为本发明真实Cuprite数据X波段的示意图;
图6(a)-图6(l)为本发明Cuprite数据SSTV-MLNMF算法解混的丰度结 果;其中(a)Alunite;(b)Andradite;(c)Buddingtonite;(d)Chalcedony;(e) Dumortierite.(f)Kaolinite 1;(g)Kaolinite 2;(h)Montmorillonite;(i)Muscovite;(j) Nontronite;(k)Pyrope;(l)Sphene
具体实施方式
步骤一、设高光谱图像矩阵为Y,端元光谱矩阵为M,丰度矩阵为R,随机 噪声矩阵为E,则线性混合模型:
Y=MR+E (1)
其中,Y是H×N维的矩阵,M是H×P维的矩阵,R是P×N维矩阵,E是H×N 维矩阵。
步骤二、在建立了线性混合模型之后,应用多层非负矩阵分解(MLNMF) 来解决解混问题,其目标函数如下所示:
其中,Yl、Ml及Rl为矩阵Y、M、R的每一层,它们的关系如下:
此外,由于高光谱图像具有稀疏性,因此在MLNMF中我们也将L1/2稀疏性加 入其目标函数中,其目标函数为:
其中
步骤三、根据高光谱图像特性,分别设计光谱域和空间域总变分最小函数。 光谱域的总变分定义如下:
其中,Mj表示矩阵M中的第j列,也就是第j个端元;在这里,一维向量的总变 分定义如下:
空间域的总变分定义如下:
其中,Rj表示矩阵R中的第j行,也就是第j个端元的丰度值,F表示形状转换 操作,把N个像素的一维向量转化为m×n的图像,在这里m×n指的是高光谱图 像中图像的大小。二维图像的总变分定义如下:
步骤四、将光谱域和空间域总变分最小函数引入MLNMF,其目标函数为:
步骤五、建立优化模型
为了方便计算,上面的式子可以转化为如下形式:
类似的,也可以推出Ml的迭代方式如下:
算法SSTV-MLNMF实现高光谱解混
实现发明所提出的算法时,需要注意几个问题:首先是初始化矩阵M和R的 方法,而此初始化方法往往会影响最终的解混效果。在高光谱图像解混中,通常 有两种常见的方法初始化M和R矩阵,随机初始化和VCA-FCLS初始化。前者 随机选择0和1之间的值作为M和R的初始值,后者使用VCA来鉴定端元M 矩阵,然后使用FCLS来获取R的初始矩阵。本发明采用VCA-FCLS初始化方 法对矩阵M和R进行初始化。
第二个问题是在迭代的过程中保证丰度的非负性(Abundance NonnegativeConstraint,ANC)和丰度加和唯一性(Abundance Sum to one Constraint,ASC)。 在这里使用一个简单且有效的方法,当R迭代之后,矩阵和将取代原Y和M 作为输入,其定义为:
其中,δ调整了ASC的效果,在实验中选择δ=15。
最后,迭代停止标准也是一个重要问题。在这里选用两种方式作为迭代停止 的条件:
1.将目标函数的减小值如下所示:
|objnew-obj|≤ε (18)
作为一个停止条件,如果两次误差小于ε=10-3,我们就设置迭代结束。
2.最大迭代次数也被作为一个停止的标准,最大迭代次数被设置为3000。
仿真实验:
(一)性能指标
在本发明中,使用光谱角距离(Spectral Angle Distance,SAD)和均方根误差(Root Mean Square Error,RMSE)两个普遍采用的度量标准,来衡量端元及其丰度的准 确性。对于第p个端元,假设对应的向量是Mp,估计的光谱角距离(SAD), 定义为:
(二)模拟数据实验
为了产生模拟数据集,首先,从USGS光谱库中随机选择6个光谱,如图二 所示。其次,通过以下步骤来产生模拟数据集:(1)将一个z2×z2的图像分割成 z×z的块;(2)每一块以相同特征的光谱(在6个光谱中随机选择)填充;(3) 一个(z+1)×(z+1)的低通滤波器被用来形成混合像元;(4)对于丰度值大于 0.8的像元,每个端元的丰度值设为1/P。
在完成上述步骤之后,合成的图像中就不存在纯像元,低通滤波器不仅能产 生高度混合的像元,还能按照现实情况使丰度变化平稳。最后,高斯白噪声加入 到合成数据中,在这里信噪比(SNR)被定义为:
其中y和δ分别表示观测图像和图像噪声,E[]表示期望值。
不同方法对比实验:对VCA、L1/2-NMF、MLNMF及SSTV-MLNMF算法进行 解混效果比较。
抗噪性能实验:对VCA、L1/2-NMF、MLNMF及SSTV-MLNMF算法不同噪声 情况下(15、25、35dB、不含噪声)进行解混效果比较。
通过观察实验结果,可以得出本发明所提出的SSTV-MLNMF算法解混性能优于VCA、L1/2-NMF、MLNMF常用算法。
(三)真实数据实验
在真实数据实验中,采用已被广泛应用的Cuprite高光谱图像来对我们的算法 进行评价。在此高光谱数据中,其一共含有224个波段,其覆盖0.4-2.5μm。在 解混前,少量的低信噪比和水蒸气吸收谱带(1-2,104-113,148-167,221-224等 波段)被去除,实验中留下188条波段。图五展示了大小为250×191的原始图像 的第50、100、150波段。根据现有技术分析,图像中可估计出14种矿物。由于 不同化学成分的同一种矿物的某些光谱之间有细微的差别,所以解混的端元数量 减少到12个。图六阐明了SSTV-MLNMF解混的最终丰度值,其结果合理,表1 为四种算法SAD结果,其结果表明我们所提出的算法解混效果优于其它算法。
表1 Cuprite实验数据SAD值
(四)总结
本发明针对高光谱图像中光谱域以及空间域总变分最小的特点,提出一种基 于光谱和空间总变分最小限制的多层非负矩阵分解高光谱图像解混方法。通过模 拟数据与真实数据实验分析,表明本发明所提出的方法能有效利用光谱域以及空 间域总变分最小的特点,从而提高高光谱图像解混的精度。
Claims (4)
1.基于光谱和空间总变分最小限制的多层非负矩阵分解的高光谱图像解混方法,其特征在于:该方法由以下步骤实现,
步骤一、设高光谱图像矩阵Y,端元光谱矩阵M,丰度矩阵R,随机噪声E,建立线性光谱混合模型:
Y=MR+E (1)
步骤二、将多层非负矩阵分解应用于线性混合模型:
Y=M1,…,MlRl+E (2)
其中M1,…,Ml为第l层;
步骤三、设计光谱域、空间域总变分函数;
步骤四、将光谱域和空间域总变分最小约束引入MLNMF,建立SSTV-MLNMF目标函数;
步骤五、对得到的目标函数进行优化;
步骤六、选定实验参数,得到最终的解混结果;
在建立了线性混合模型之后,应用多层非负矩阵分解来解决解混问题,其目标函数如下所示:
其中,Yl、Ml及Rl为矩阵Y、M、R的每一层,它们的关系如下:
此外,由于高光谱图像具有稀疏性,因此在MLNMF中我们也将L1/2稀疏性加入其目标函数中,其目标函数为:
其中
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910320619.3A CN110570359B (zh) | 2019-04-20 | 2019-04-20 | 基于光谱和空间总变分最小限制的多层非负矩阵分解高光谱图像解混方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910320619.3A CN110570359B (zh) | 2019-04-20 | 2019-04-20 | 基于光谱和空间总变分最小限制的多层非负矩阵分解高光谱图像解混方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110570359A CN110570359A (zh) | 2019-12-13 |
CN110570359B true CN110570359B (zh) | 2023-01-20 |
Family
ID=68773435
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910320619.3A Active CN110570359B (zh) | 2019-04-20 | 2019-04-20 | 基于光谱和空间总变分最小限制的多层非负矩阵分解高光谱图像解混方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110570359B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113435366A (zh) * | 2021-06-30 | 2021-09-24 | 南京理工大学 | 一种小波域的多时高光谱图像贝叶斯解混方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105809105A (zh) * | 2016-02-06 | 2016-07-27 | 黑龙江科技大学 | 基于端元约束非负矩阵分解的高光谱图像解混方法 |
CN109085131A (zh) * | 2018-07-12 | 2018-12-25 | 重庆邮电大学 | 基于丰度稀疏和端元正交性约束nmf的高光谱解混方案 |
-
2019
- 2019-04-20 CN CN201910320619.3A patent/CN110570359B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105809105A (zh) * | 2016-02-06 | 2016-07-27 | 黑龙江科技大学 | 基于端元约束非负矩阵分解的高光谱图像解混方法 |
CN109085131A (zh) * | 2018-07-12 | 2018-12-25 | 重庆邮电大学 | 基于丰度稀疏和端元正交性约束nmf的高光谱解混方案 |
Non-Patent Citations (2)
Title |
---|
"An Endmember Dissimilarity Constrained Non-Negative Matrix Factorization Method for Hyperspectral Unmixing";Nan Wang 等;《 IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing》;20130221;第6卷(第2期);全文 * |
"Nonnegative-matrix-factorization-based hyperspectral unmixing with partially known endmembers";Lei Tong 等;《IEEE Transactions on Geoscience and Remote Sensing》;20160727;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110570359A (zh) | 2019-12-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111260576B (zh) | 一种基于去噪三维卷积自编码网络的高光谱解混算法 | |
CN109241843B (zh) | 空谱联合多约束优化非负矩阵解混方法 | |
CN110660038A (zh) | 一种基于生成对抗网络的多光谱图像与全色图像融合方法 | |
CN107609573B (zh) | 基于低秩分解和空谱约束的高光谱图像时变特征提取方法 | |
CN103745487B (zh) | 基于结构化稀疏先验的贝叶斯高光谱解混压缩感知方法 | |
CN111080555B (zh) | 一种基于三维拟递归神经网络的高光谱图像降噪方法 | |
CN104463223B (zh) | 基于空谱信息丰度约束的高光谱图像组稀疏解混方法 | |
CN103528968A (zh) | 基于迭代阈值法的反射率光谱重建方法 | |
CN111598786A (zh) | 一种基于深度去噪自编码网络的高光谱图像解混方法 | |
CN115311187B (zh) | 基于内部和外部先验的高光谱融合成像方法、系统及介质 | |
CN109785242B (zh) | 一种基于逐波段广义双线性模型的高光谱图像的解混方法 | |
CN112504975B (zh) | 一种基于约束非负矩阵分解的高光谱解混方法 | |
CN112699838B (zh) | 基于光谱诊断特征加权的高光谱混合像元非线性盲分解方法 | |
CN109671019B (zh) | 一种基于多目标优化算法和稀疏表达的遥感影像亚像元制图方法 | |
CN111161141A (zh) | 基于频带间注意力机制对抗学习高光谱单图超分辨率方法 | |
CN112712034A (zh) | 一种高光谱图像的解混方法及系统 | |
CN106067165B (zh) | 基于聚类化稀疏随机场的高光谱图像去噪方法 | |
CN110570359B (zh) | 基于光谱和空间总变分最小限制的多层非负矩阵分解高光谱图像解混方法 | |
CN115272093A (zh) | 一种基于空间结构信息约束的高光谱图像解混方法 | |
CN109727280B (zh) | 一种基于正交基的高光谱图像丰度估计方法 | |
CN108734672B (zh) | 基于光谱库裁剪与协同稀疏回归的高光谱数据解混方法 | |
Wang et al. | Hyperspectral unmixing via plug-and-play priors | |
CN117350373A (zh) | 一种基于局部自注意力机制的个性化联邦聚合算法 | |
CN104899850B (zh) | 一种基于加权联合稀疏回归的高光谱图像解混方法 | |
CN116626765A (zh) | 一种基于二维k-svd和卷积稀疏编码的多道地震反褶积方法 |
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 |