CN114841875B - 一种基于图学习和降噪的高光谱图像解混方法 - Google Patents
一种基于图学习和降噪的高光谱图像解混方法 Download PDFInfo
- Publication number
- CN114841875B CN114841875B CN202210432150.4A CN202210432150A CN114841875B CN 114841875 B CN114841875 B CN 114841875B CN 202210432150 A CN202210432150 A CN 202210432150A CN 114841875 B CN114841875 B CN 114841875B
- Authority
- CN
- China
- Prior art keywords
- matrix
- unmixing
- hyperspectral image
- image
- hyperspectral
- 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
- 238000000034 method Methods 0.000 title claims abstract description 96
- 239000011159 matrix material Substances 0.000 claims abstract description 70
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 230000004927 fusion Effects 0.000 claims description 4
- 238000002156 mixing Methods 0.000 claims description 4
- 230000017105 transposition Effects 0.000 claims description 4
- 230000001737 promoting effect Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 229940050561 matrix product Drugs 0.000 claims 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 11
- 229910052802 copper Inorganic materials 0.000 description 11
- 239000010949 copper Substances 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 4
- BERDEBHAJNAUOM-UHFFFAOYSA-N copper(I) oxide Inorganic materials [Cu]O[Cu] BERDEBHAJNAUOM-UHFFFAOYSA-N 0.000 description 3
- LBJNMUFDOHXDFG-UHFFFAOYSA-N copper;hydrate Chemical compound O.[Cu].[Cu] LBJNMUFDOHXDFG-UHFFFAOYSA-N 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 239000010424 alunite Substances 0.000 description 2
- 229910052934 alunite Inorganic materials 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000011045 chalcedony Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- KPZTWMNLAFDTGF-UHFFFAOYSA-D trialuminum;potassium;hexahydroxide;disulfate Chemical compound [OH-].[OH-].[OH-].[OH-].[OH-].[OH-].[Al+3].[Al+3].[Al+3].[K+].[O-]S([O-])(=O)=O.[O-]S([O-])(=O)=O KPZTWMNLAFDTGF-UHFFFAOYSA-D 0.000 description 2
- 229910001595 Buddingtonite Inorganic materials 0.000 description 1
- 238000012899 de-mixing Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- 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
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于图学习和降噪的高光谱图像解混方法,所述方法包括如下步骤:步骤1:从观测高光谱图像中图学习;步骤2:高光谱图像的降噪处理;步骤3:基于降噪后的高光谱图像Y解混;步骤4:终止条件。本发明利用图学习的方法从含有噪声的观测的高光谱图像中学习得到一个空间Laplacian矩阵,通过空间Laplacian矩阵可自适应不同情况的高光谱图像。该Laplacian矩阵保留了实际的高光谱图像的像素间的真实空间结构关系的同时还能够用于去除高光谱图像中存在的噪声,然后再基于去噪后的高光谱图像进行解混处理,这可在一定程度上提升高光谱图像解混的性能,能够自适应不同情况的高光谱图像,更加鲁棒和准确。
Description
技术领域
本发明属于遥感技术的研究领域,涉及一种高光谱图像解混方法,具体涉及一种对降噪后的高光谱图像的像素进行解混处理并得到相应地表物质丰度表示的方法。
背景技术
在高光谱遥感中,高光谱图像记录了场景中的地表物质的光谱特征,它的每个像素的光谱包含数百个连续波段丰富的光谱信息。高光谱图像已广泛用于各个领域的许多应用,包括矿产勘探、环境监测、农业和军事应用等。在实践中,由于高光谱传感器的低空间分辨率和地表物质分布的高度复杂性,导致了高光谱图像包含高度的混合像素,其每个像素的光谱是来自多种不同材料的光谱特征的混合。混合像素的存在严重限制了高光谱图像的应用,例如,它会给高光谱图像分析在聚类和分类任务中带来巨大挑战。对高光谱的每个像素进行解混是解决这一问题的有效的方法,其已经成为了高光谱图像分析与进一步开发的必要过程。因而,需要对高光谱图像进行有效的解混处理。
在实际场景的高光谱图像中,由于光照条件、大气效应以及设备自身精度等因素的干扰,导致了实际的高光谱图像中存在大量的干扰信息,这些干扰信息可视作是高光谱图像中的噪声。目前已有的稀疏解混方法都是基于含有噪声的观测的高光谱图像进行解混,但由于观测图像受到了噪声的干扰,这会在很大程度上影响对高光谱图像的混合像素进行解混处理。
发明内容
本发明基于观测的高光谱遥感图像中存在噪声干扰这一动机,提供了一种基于图学习和降噪的高光谱图像解混方法。该方法基于降噪后的高光谱图像进行解混,去除了噪声对解混过程的干扰,可进一步提升高光谱图像的解混性能,可为相关的专业人员对高光谱图像的进一步分析与应用提供一定的理论分析依据。
本发明的目的是通过以下技术方案实现的:
一种基于图学习和降噪的高光谱图像解混方法,包括如下步骤:
步骤1:从观测高光谱图像中图学习
首先将3-D的观测高光谱图像S进行分割与融合,然后将预处理后的图像用于空间Laplacian矩阵的图学习;
步骤2:高光谱图像的降噪处理
基于步骤1图学习过程所学得的空间Laplacian矩阵L,对观测的含噪高光谱图像S按下式进行降噪处理:
Y=S(In+λlL)-1;
其中,In表示维数为n的单位阵,λl为正则参数,上角标-1表示矩阵的求逆运算,Y为降噪后的高光谱图像;
步骤3:基于降噪后的高光谱图像Y解混
步骤3.1、解混过程中变量X按下式进行计算:
X(k+1)=[Λ1Λ2]+;
其中,[z]+=max(z,0),Λ1=(ATA+μ(W(k))-2)-1,Λ2=ATY+μ(W(k))-1(V(k)-U(k)),上角标-2表示对矩阵所有元素进行平方后的求逆运算,上角标-1表示矩阵的求逆运算,上角标T表示矩阵的转置运算;
步骤3.2、解混过程内层循环的变量V按下式进行计算:
其中,H=W(k)X(k+1)+U(k),H(i,:)表示矩阵H的第i行;
步骤3.3、解混过程中Lagrangian乘子按下式进行更新:
U(k+1)=U(k)+(W(k))-1X(k+1)-V(k+1);
步骤3.4、解混过程中促进丰度稀疏性的权重矩阵W按下式计算:
其中,ε为解混过程设定的超参数,X(k+1)(i,:)表示在解混过程的第k+1次迭代时,估计矩阵X的第i行;
步骤4:终止条件
若解混过程中的迭代次数k=Qmax,则终止迭代并输出矩阵X(k+1);否则转到步骤3继续执行;当依据步骤4停止时,输出的X(k+1)就是基于降噪后的高光谱图像Y所估计出的丰度矩阵。
相比于现有技术,本发明具有如下优点:
1、本发明与目前广泛用于高光谱遥感图像解混的其他稀疏解混方法相比,所提出的解混方法具有更好的解混性能和很低的计算代价。
2、本发明基于降噪后的高光谱图像进行解混,通过使用观测的高光谱图像利用图学习的方法从含有噪声的观测的高光谱图像中学习得到一个空间Laplacian矩阵,通过学习得到的空间Laplacian矩阵可自适应不同情况的高光谱图像。该Laplacian矩阵保留了实际的高光谱图像的像素间的真实空间结构关系的同时还能够用于去除高光谱图像中存在的噪声,然后再基于去噪后的高光谱图像进行解混处理,这可在一定程度上提升高光谱图像解混的性能,能够自适应不同情况的高光谱图像,更加鲁棒和准确。因此,本发明的方法具有计算代价低、可靠性好以及精度相对较高的特点,适合推广使用。
3、本发明基于国际公开的光谱数据库,使用了被广泛应用于高光谱图像解混任务中的两个合成数据集DC1、DC2和一个真实的铜矿石(Cuprite)数据集进行实验,这些数据集可被用于定量分析高光谱解混方法的性能。基于这些数据集的实验表明,本发明中所提出的解混方法远好于当前其他的稀疏解混方法。
附图说明
图1为采用现有的SUnSAL、SUnSAL-TV、S2WSU、MUA-SLIC及SBGLSU五种稀疏解混方法与本发明中所述的解混方法,在输入信噪比为10dB时,对高光谱图像模拟数据集DC1进行解混时的性能表现图;
图2为采用现有的SUnSAL、SUnSAL-TV、S2WSU、MUA-SLIC及SBGLSU五种稀疏解混方法与本发明中所述的解混方法,在输入信噪比为10dB时,对高光谱图像模拟数据集DC2进行解混时的性能表现图;
图3为采用现有的SUnSAL、SUnSAL-TV、S2WSU、MUA-SLIC及SBGLSU五种稀疏解混方法与本发明中所述的解混方法,对真实的铜矿石(Cuprite)高光谱图像的数据集进行解混时的性能表现图;
图中,Ground Truth为用于衡量解混性能的真实参考丰度图,SUnSAL、SUnSAL-TV、S2WSU、MUA-SLIC以及SBGLSU为五种稀疏解混对比方法所生成的丰度图,GLDWSU为本发明所提出的方法所生成的丰度图。
具体实施方式
下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。
本发明提供了一种基于图学习和降噪的高光谱图像解混方法,所述方法包括如下步骤:
步骤1:从观测高光谱图像中图学习
为了同时减少高光谱图像谱域中的冗余信息和计算成本,首先将3-D的观测高光谱图像S进行分割与融合,然后将预处理后的图像更好的用于空间Laplacian矩阵的图学习,具体步骤如下:
步骤1.1:首先对观测图像S进行分割处理,将图像S的每个波段的图像视为一个子图像,则这些子图像可用集合表示为其中,si表示3-D观测高光谱图像的第i个波段的子图像,d为总的波段数;然后将d个子图像划分为由相邻波段所构成的多个子集/>其中,K是子集的数量,即:/>其中,Nk是由与第k个波段相邻的多个波段的子图像所构成的索引集合;最后,每个子集合Si由下式进行融合:
其中,|Nk|表示集合的大小,包含了每个子集Sk中的平均空间信息。
在预处理过程中,不同的数据集的参数设置如下:
模拟数据集DC1中子集数量K:K=5。
模拟数据集DC2中子集数量K:K=5。
真实的铜矿石数据集中子集数量K:K=20。
步骤1.2:基于融合的图像来学习空间Laplacian矩阵,在Laplacian矩阵的学习过程中,相关参数及输入量设置如下:
所有数据集在图学习过程中的初始化设置:其中,k=1,…,K。
模拟数据集DC1的相关参数设置:高斯滤波器的标准差:σ=1.2;图学习的正则参数:λl=0.003;超参数:ξ=0.001;图学习过程的最大迭代次数:Tmax=5。
模拟数据集DC2的相关参数设置:高斯滤波器的标准差:σ=0.5;图学习的正则参数:λl=0.008;超参数:ξ=0.01;图学习过程的最大迭代次数:Tmax=5。
铜矿石数据集的相关参数设置:高斯滤波器的标准差:σ=0.1;图学习的正则参数:λl=0.05;超参数:ξ=15;图学习过程的最大迭代次数:Tmax=10。
具体步骤如下:
步骤1.2.1、计算对角权重矩阵Wr和Wc按下式计算:
其中,Dr和Dc分别为水平和垂直两个方向上的离散差分算子,Gσ是标准差为σ的高斯滤波器,ξ为图学习过程设定的超参数,“*”为卷积算子,为基于元素运算的乘积算子,以及除法和绝对值都是基于元素的运算;
步骤1.2.2、计算空间Laplacian矩阵
其中,Dr和Dc分别为水平和垂直两个方向上的离散差分算子,上角标T为矩阵的转置运算;
步骤1.2.3、在图学习的第t次迭代过程中,变量按下式进行更新:
其中,In表示维数为n的单位阵,上角标-1表示矩阵求逆运算;
步骤1.2.4、图学习过程的终止条件:
若迭代次数t=Tmax,则终止迭代并按下式输出空间Laplacian矩阵L:
否则空间Laplacian矩阵的图学习转到步骤1.2.1继续执行;
步骤2:高光谱图像的降噪处理
基于步骤1图学习过程所学得的空间Laplacian矩阵L,对观测的含噪高光谱图像S按下式进行降噪处理:
Y=S(In+λlL)-1 (7);
其中,In表示维数为n的单位阵,λl为步骤1中所设定的正则参数,上角标-1表示矩阵的求逆运算,Y为降噪后的高光谱图像。
步骤3:基于降噪后的高光谱图像Y解混
所有降噪后的高光谱图像数据集的解混过程的相关参数及输入量设置如下:
模拟数据集DC1、DC2的相关参数设置:正则参数:λu=0.5;罚参数:μ=0.01;超参数:ε=0.0001;解混过程的最大迭代次数设置:Qmax=100。
真实铜矿石数据集的相关参数设置:正则参数:λu=0.1;罚参数:μ=0.6;超参数:ε=0.0001;解混过程的最大迭代次数设置:Qmax=200。
算法所要估计的丰度矩阵的初始化X(0)设置如下:
式(8)中的变量说明:“svd”表示奇异值分解运算,“*”为矩阵乘积运算,上角标T表示矩阵的转置运算,A为已知的谱库矩阵,Y为降噪的高光谱图像,“diag(·)”表示构建一个对角阵,或以向量形式返回一个矩阵对角线上的元素,μ为解混过程设定的罚参数。
解混过程的辅助变量V(0)的初适化设置:V(0)=Y;
解混过程的Lagrangian乘子初始化设置:U(0)为全零矩阵。
具体步骤如下:
步骤3.1、解混过程中变量X按下式进行计算:
X(k+1)=[Λ1Λ2]+ (9);
其中,[z]+=max(z,0),Λ1=(ATA+μ(W(k))-2)-1,
Λ2=ATY+μ(W(k))-1(V(k)-U(k)),上角标-2表示对矩阵所有元素进行平方后的求逆运算,上角标-1表示矩阵的求逆运算,上角标T表示矩阵的转置运算;
步骤3.2、解混过程内层循环的变量V按下式进行计算:
其中,H=W(k)X(k+1)+U(k),H(i,:)表示矩阵H的第i行;
步骤3.3、解混过程中Lagrangian乘子按下式进行更新:
U(k+1)=U(k)+(W(k))-1X(k+1)-V(k+1) (11);
步骤3.4、解混过程中促进丰度稀疏性的权重矩阵W按下式计算:
其中,ε为解混过程设定的超参数,X(k+1)(i,:)表示在解混过程的第k+1次迭代时,估计矩阵X的第i行,i=1,…,m。
步骤4:算法的终止条件
若解混过程中的迭代次数k=Qmax,则终止迭代并输出矩阵X(k+1);否则算法转到步骤3继续执行;当算法依据步骤4停止时,算法输出的X(k+1)就是基于降噪后的高光谱图像Y所估计出的丰度矩阵。
为验证本发明中所提出的解混方法的有效性,我们分别采用SUnSAL、SUnSAL-TV、S2WSU、MUA-SLIC以及SBGLSU五种稀疏解混方法和本发明中所述的解混方法(GLDWSU),使用了两个合成的模拟数据集DC1、DC2和一个真实的铜矿石(Cuprite)数据集进行解混处理。模拟数据集DC1的高光谱图像包含了224个波段,每个波段的图像由75×75的像素构成,包含25个以网格方式空间排列的正方形,其中每个正方形是一个同质区域,其像素以相同的丰度生成。模拟数据集DC2含有224个波段,每个波段的图像包含了100×100个像素,其通过使用以高斯随机场为中心的狄利克雷分布的相应丰度分数组成。真实的铜矿石(Cuprite)数据集包含了188个波段,每个波段的图像由250×191的像素构成。
我们在实验中,对于数据集DC1选择了第3端元、第4端元和第5端元的丰度图作为代表;对于数据集DC2选择了第3端元、第8端元和第9端元的丰度图作为代表;对于真实的铜矿石数据集选择了明矾石(Alunite)、布丁英石(Buddingtonite)以及玉髓(Chalcedony)这三种矿物作为代表,所提出的方法在每个数据集的解混性能分别如图1、图2和图3所示。从图中可以看出,本发明所提出的解混方法的性能表现都远好于同类型的其他稀疏解混方法,此实验充分证实了所提出方法的有效性。
此外,为了更进一步的说明本发明所提出方法的有效性,我们利用信号的重构误差(signal-reconstruction-error,SRE)来度量所提出解混方法的性能,该指标被广泛的应用于高光谱解混方法性能的度量中。信号重构误差SRE的具体定义方式如下:
其中,X和分别为真实的丰度图和估计的丰度图,通常情况,SRE值越大说明解混方法的性能越好。
在实验中,我们分别对模拟数据集DC1和DC2加入了信噪比为10dB的高斯噪声,以用于模拟含有噪声的高光谱图像的数据集。我们分别采用SUnSAL、SUnSAL-TV、S2WSU、MUA-SLIC以及SBGLSU五种稀疏解混方法和本发明中所述的解混方法,对加入10dB高斯噪声的模拟数据集DC1和DC2进行了解混处理,所提出的方法的SRE值和计算时间如表1和表2所示。从表1和表2可以看出,本发明所提出的解混方法的性能表现都远好于同类型的其他稀疏解混方法,同时还具有很低的计算代价。此实验更进一步的证实了本发明所提出方法的有效性。
表1各种方法在模拟数据集DC1、DC2上输入信噪比为10dB时的解混SRE值
Data | SREValue | SUnSAL | SUnSAL-TV | S2WSU | MUA-SLIC | SBGLSU | GLDWSU |
DC1 | SRE(dB) | 0.07 | 4.67 | 3.65 | 7.89 | 7.62 | 16.81 |
DC2 | SRE(dB) | 1.37 | 2.10 | 3.46 | 9.10 | 8.73 | 14.63 |
表2各种方法在所有高光谱数据集上的计算时间
Data | SUnSAL | SUnSAL-TV | S2WSU | MUA-SLIC | SBGLSU | GLDWSU |
DC1 | 2.46s | 57.53s | 13.56s | 1.94s | 26.32s | 2.78s |
DC2 | 3.93s | 88.92s | 21.98s | 2.51s | 58.93s | 4.60s |
RealCupritedata | 52.77s | 693.44s | 166.71s | 47.15s | 777.32s | 86.56s |
注:由于铜矿石数据是真实的高光谱图像数据,不存在真实的丰度值作为参考,因此没有列出关于铜矿石数据集的性能值的表格。
Claims (2)
1.一种基于图学习和降噪的高光谱图像解混方法,其特征在于所述方法包括如下步骤:
步骤1:从观测高光谱图像中图学习
首先将3-D的观测高光谱图像S进行分割与融合,然后将预处理后的图像用于空间Laplacian矩阵的图学习,具体步骤如下:
步骤1.1:首先对观测图像S进行分割处理,将图像S的每个波段的图像视为一个子图像,则这些子图像可用集合表示为其中,si表示3-D观测高光谱图像的第i个波段的子图像,d为总的波段数;然后将d个子图像划分为由相邻波段所构成的多个子集/>其中,K是子集的数量,即:/>其中,Nk是由与第k个波段相邻的多个波段的子图像所构成的索引集合;最后,每个子集合Si由下式进行融合:
其中,|Nk|表示集合的大小,包含了每个子集Sk中的平均空间信息;
步骤1.2:基于融合的图像来学习空间Laplacian矩阵;
步骤2:高光谱图像的降噪处理
基于步骤1图学习过程所学得的空间Laplacian矩阵L,对观测的含噪高光谱图像S按下式进行降噪处理:
Y=S(In+λlL)-1;
其中,In表示维数为n的单位阵,λl为正则参数,上角标-1表示矩阵的求逆运算,Y为降噪后的高光谱图像;
步骤3:基于降噪后的高光谱图像Y解混
步骤3.1、解混过程中变量X按下式进行计算:
X(k+1)=[Λ1Λ2]+;
其中,[z]+=max(z,0),Λ1=(ATA+μ(W(k))-2)-1,Λ2=ATY+μ(W(k))-1(V(k)-U(k)),上角标-2表示对矩阵所有元素进行平方后的求逆运算,上角标-1表示矩阵的求逆运算,上角标T表示矩阵的转置运算;
步骤3.2、解混过程内层循环的变量V按下式进行计算:
其中,H=W(k)X(k+1)+U(k),H(i,:)表示矩阵H的第i行;
步骤3.3、解混过程中Lagrangian乘子按下式进行更新:
U(k+1)=U(k)+(W(k))-1X(k+1)-V(k+1);
步骤3.4、解混过程中促进丰度稀疏性的权重矩阵W按下式计算:
其中,ε为解混过程设定的超参数,X(k+1)(i,:)表示在解混过程的第k+1次迭代时,估计矩阵X的第i行;
丰度矩阵的初始化X(0)设置如下:
[UF,SF]=svd(ATA)
X(0)=IF*AT*Y
其中,“svd”表示奇异值分解运算,“*”为矩阵乘积运算,上角标T表示矩阵的转置运算,A为已知的谱库矩阵,Y为降噪的高光谱图像,“diag(·)”表示构建一个对角阵,或以向量形式返回一个矩阵对角线上的元素,μ为解混过程设定的罚参数;
解混过程的辅助变量V(0)的初始化设置:V(0)=Y;
解混过程的Lagrangian乘子初始化设置:U(0)为全零矩阵;
步骤4:终止条件
若解混过程中的迭代次数k=Qmax,则终止迭代并输出矩阵X(k+1);否则转到步骤3继续执行;当依据步骤4停止时,输出的X(k+1)就是基于降噪后的高光谱图像Y所估计出的丰度矩阵。
2.根据权利要求1所述的基于图学习和降噪的高光谱图像解混方法,其特征在于所述步骤1.2的具体步骤如下:
步骤1.2.1、计算对角权重矩阵Wr和Wc按下式计算:
其中,Dr和Dc分别为水平和垂直两个方向上的离散差分算子,Gσ是标准差为σ的高斯滤波器,ξ为图学习过程设定的超参数,“*”为卷积算子,为基于元素运算的乘积算子;
步骤1.2.2、计算空间Laplacian矩阵
其中,Dr和Dc分别为水平和垂直两个方向上的离散差分算子,上角标T为矩阵的转置运算;
步骤1.2.3、在图学习的第t次迭代过程中,变量按下式进行更新:
其中,In表示维数为n的单位阵,上角标-1表示矩阵求逆运算;在Laplacian矩阵的学习过程中,所有数据集在图学习过程中的初始化设置:其中,k=1,…,K;
步骤1.2.4、图学习过程的终止条件:
若迭代次数t=Tmax,则终止迭代并按下式输出空间Laplacian矩阵L:
否则空间Laplacian矩阵的图学习转到步骤1.2.1继续执行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210432150.4A CN114841875B (zh) | 2022-04-22 | 2022-04-22 | 一种基于图学习和降噪的高光谱图像解混方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210432150.4A CN114841875B (zh) | 2022-04-22 | 2022-04-22 | 一种基于图学习和降噪的高光谱图像解混方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114841875A CN114841875A (zh) | 2022-08-02 |
CN114841875B true CN114841875B (zh) | 2023-08-11 |
Family
ID=82565600
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210432150.4A Active CN114841875B (zh) | 2022-04-22 | 2022-04-22 | 一种基于图学习和降噪的高光谱图像解混方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114841875B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115829886B (zh) * | 2022-12-21 | 2023-08-11 | 哈尔滨师范大学 | 基于端元自适应不相干和空间约束的盲高光谱解混方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104008574A (zh) * | 2014-06-16 | 2014-08-27 | 浙江大学 | 一种基于无限高斯混合模型的高光图图像解混方法 |
CN110070518A (zh) * | 2019-03-15 | 2019-07-30 | 南京航空航天大学 | 一种基于双路径支持下的高光谱图像超分辨率制图方法 |
CN110309851A (zh) * | 2019-05-16 | 2019-10-08 | 中国科学院西安光学精密机械研究所 | 基于子空间聚类约束的高光谱影像解混方法、计算机可读存储介质、电子设备 |
CN110363712A (zh) * | 2019-06-14 | 2019-10-22 | 江苏理工学院 | 一种稀疏对偶约束的高光谱图像解混方法 |
CN111260576A (zh) * | 2020-01-14 | 2020-06-09 | 哈尔滨工业大学 | 一种基于去噪三维卷积自编码网络的高光谱解混算法 |
CN113661381A (zh) * | 2019-01-31 | 2021-11-16 | 南加州大学 | 高光谱成像系统 |
-
2022
- 2022-04-22 CN CN202210432150.4A patent/CN114841875B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104008574A (zh) * | 2014-06-16 | 2014-08-27 | 浙江大学 | 一种基于无限高斯混合模型的高光图图像解混方法 |
CN113661381A (zh) * | 2019-01-31 | 2021-11-16 | 南加州大学 | 高光谱成像系统 |
CN110070518A (zh) * | 2019-03-15 | 2019-07-30 | 南京航空航天大学 | 一种基于双路径支持下的高光谱图像超分辨率制图方法 |
CN110309851A (zh) * | 2019-05-16 | 2019-10-08 | 中国科学院西安光学精密机械研究所 | 基于子空间聚类约束的高光谱影像解混方法、计算机可读存储介质、电子设备 |
CN110363712A (zh) * | 2019-06-14 | 2019-10-22 | 江苏理工学院 | 一种稀疏对偶约束的高光谱图像解混方法 |
CN111260576A (zh) * | 2020-01-14 | 2020-06-09 | 哈尔滨工业大学 | 一种基于去噪三维卷积自编码网络的高光谱解混算法 |
Also Published As
Publication number | Publication date |
---|---|
CN114841875A (zh) | 2022-08-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111260576B (zh) | 一种基于去噪三维卷积自编码网络的高光谱解混算法 | |
CN107705306B (zh) | 一种基于多特征矩阵低秩分解的织物疵点检测方法 | |
CN115272093B (zh) | 一种基于空间结构信息约束的高光谱图像解混方法 | |
CN103413151B (zh) | 基于图正则低秩表示维数约简的高光谱图像分类方法 | |
CN110660038A (zh) | 一种基于生成对抗网络的多光谱图像与全色图像融合方法 | |
CN108197650B (zh) | 局部相似性保持的高光谱图像极限学习机聚类方法 | |
CN104978573A (zh) | 一种应用于高光谱图像处理的非负矩阵分解方法 | |
CN108229551B (zh) | 一种基于紧凑字典稀疏表示的高光谱遥感图像分类方法 | |
CN107491793B (zh) | 一种基于稀疏散射全卷积的极化sar图像分类方法 | |
CN107609573A (zh) | 基于低秩分解和空谱约束的高光谱图像时变特征提取方法 | |
CN109671029A (zh) | 基于伽马范数最小化的图像去噪算法 | |
CN114841875B (zh) | 一种基于图学习和降噪的高光谱图像解混方法 | |
CN102682441A (zh) | 基于亚像元映射的高光谱图像超分辨重建方法 | |
CN104268561B (zh) | 基于结构先验低秩表示的高光谱图像解混方法 | |
CN107392863A (zh) | 基于亲和矩阵融合谱聚类方法的sar图像变化检测方法 | |
CN108256557B (zh) | 结合深度学习和邻域集成的高光谱图像分类方法 | |
CN103226825B (zh) | 基于低秩稀疏模型的遥感图像变化检测方法 | |
CN116403046A (zh) | 一种高光谱影像分类装置及方法 | |
CN110673138B (zh) | 一种基于奇异值分解和模糊c均值法的探地雷达图像处理方法 | |
CN115457311A (zh) | 一种基于自表达迁移学习的高光谱遥感影像波段选择方法 | |
CN107392211A (zh) | 基于视觉稀疏认知的显著目标检测方法 | |
CN107368846A (zh) | 基于小波变换和稀疏表示的高光谱图像分类方法 | |
CN107273919A (zh) | 一种基于置信度构造类属字典的高光谱无监督分类方法 | |
CN113421198A (zh) | 一种基于子空间的非局部低秩张量分解的高光谱图像去噪方法 | |
CN113392871A (zh) | 一种基于散射机制多通道扩张卷积神经网络的极化sar地物分类方法 |
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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Deng Shiwen Inventor after: Song Fuxin Inventor after: Kan Chao Inventor after: Li Jia Inventor after: Diao Jie Inventor before: Kan Chao Inventor before: Deng Shiwen Inventor before: Diao Jie Inventor before: Li Jia Inventor before: Song Fuxin |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |