CN109102477B - 一种基于非凸低秩稀疏约束的高光谱遥感图像恢复方法 - Google Patents
一种基于非凸低秩稀疏约束的高光谱遥感图像恢复方法 Download PDFInfo
- Publication number
- CN109102477B CN109102477B CN201811018402.9A CN201811018402A CN109102477B CN 109102477 B CN109102477 B CN 109102477B CN 201811018402 A CN201811018402 A CN 201811018402A CN 109102477 B CN109102477 B CN 109102477B
- Authority
- CN
- China
- Prior art keywords
- remote sensing
- sensing image
- hyperspectral remote
- image
- noise
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000011084 recovery Methods 0.000 title claims abstract description 54
- 239000011159 matrix material Substances 0.000 claims abstract description 46
- 238000005457 optimization Methods 0.000 claims description 8
- 238000000354 decomposition reaction Methods 0.000 claims description 5
- 238000012887 quadratic function Methods 0.000 claims description 5
- 238000002939 conjugate gradient method Methods 0.000 claims description 4
- 230000015556 catabolic process Effects 0.000 claims description 3
- HHHKFGXWKKUNCY-FHWLQOOXSA-N cilazapril Chemical compound C([C@@H](C(=O)OCC)N[C@@H]1C(N2[C@@H](CCCN2CCC1)C(O)=O)=O)CC1=CC=CC=C1 HHHKFGXWKKUNCY-FHWLQOOXSA-N 0.000 claims description 3
- 238000006731 degradation reaction Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 22
- 238000005516 engineering process Methods 0.000 abstract description 6
- 238000012545 processing Methods 0.000 abstract description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000000701 chemical imaging Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000001595 contractor effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 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
- 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
-
- 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/20021—Dividing image into blocks, subimages or windows
-
- 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/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
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
基于非凸低秩稀疏约束的高光谱遥感图像恢复方法,属于遥感图像处理中的高光谱遥感图像处理领域,为解决现有的高光谱遥感图像恢复技术不能有效去除噪声的问题,提高图像恢复质量,包括以下步骤:输入高光谱遥感图像;初始化权重系数矩阵、迭代次数和收敛阈值,初始化子图像尺寸和扫描步长,划分子块;建立图像恢复模型;引入辅助变量和正则项系数,利用最大值最小值算法迭代求解;判断恢复结果是否满足收敛条件;达到迭代次数获得符合要求的高光谱恢复图像,否则返回相应步骤继续迭代操作;计算权重系数矩阵,给每个子块分配合适的权重;还原高光谱遥感图像,得到最终恢复的高光谱遥感图像。去噪效果明显,并且保留图像细节。
Description
技术领域
本发明属于遥感图像处理中的高光谱图像处理领域,具体涉及一种基于非凸低秩稀疏约束模型的高光谱遥感图像恢复方法。
发明背景
高光谱遥感成像技术结合光谱分析和光学成像技术,探测目标的二维几何空间及一维光谱信息,获取高分辨率的连续、窄波段的图像数据。目前,高光谱成像技术发展迅速,因其蕴含丰富的地物光谱信息而被广泛应用于环境研究、地质勘探和军事监视等领域。但是由于传感器的物理缺陷、光子效应、传输损失和校准误差等原因,实际中获得的高光谱遥感图像常常会出现各种类型的退化,比如:噪音污染、条纹影响和数据缺失等。这不仅严重降低了图像的质量,还在很大程度上影响了分类、目标检测和解混等高光谱遥感图像的进一步应用。因此,高光谱遥感图像的恢复是一个重要的亟待解决的问题。
目前,多种不同高光谱遥感图像的恢复方法被应用于高光谱遥感图像的预处理工作中。近几年有学者提出利用高光谱遥感图像数据的相关性构建稀疏和低秩约束优化方程恢复高光谱遥感图像的方法。比如利用利用高光谱遥感图像低秩特性的低秩矩阵恢复(Low-Rank Matrix Recovery,LRMR)方法、噪声自适应迭代低秩矩阵分解(Noise-AdjustedIterative Low-Rank Matrix Approximation,NAILRMA)以及利用高光谱遥感图像稀疏特性的全变分正则化低秩矩阵分解(Total-Variation-Regularized Low-Rank MatrixApproximation,LRTV)方法等。但是此类算法多以低秩约束或稀疏约束的凸松弛l1范数进行约束优化模型的构建,但是l1范数由于收缩效应有时存在估计偏差,影响图像恢复质量;此外,当图像数据矩阵维度增加时,收敛速度常不能达到要求。
发明内容
本发明的目的旨在针对上述已有技术的不足,提出一种基于非凸低秩和非凸稀疏约束的高光谱遥感图像恢复方法,在去除高光谱遥感图像中噪声的同时保留图像细节,以提高高光谱遥感图像的恢复质量,满足后续应用要求。
为实现上述目的,本发明所采取的技术方案是:
一种基于非凸低秩稀疏约束的高光谱遥感图像恢复方法,所述方法包含以下步骤:
(1)输入含有混合噪声的高光谱遥感图像Y,Y的大小为M×N×B,其中M和N分别代表各波段高光谱遥感图像的长和宽,B代表高光谱遥感图像波段数;
定义恢复后的图像X*,初始化变量X(0)=Y,X*是X(0)经过一定迭代次数后的恢复图像,初始化X*的权重系数矩阵Weight,初始化迭代次数和收敛阈值;
定义子块的尺寸p和扫描步长s,并将高光谱遥感图像Y按照预设的尺寸p和扫描步长s划分为若干个部分重叠的子块;
(2)对于选取的第i个子块Yi,将受噪声污染的高光谱遥感图像分解为其中为没有噪声污染的干净的高光谱遥感图像,Ji为包括高斯噪声、条带噪声、脉冲噪声的混合噪声;根据求解逆问题并引入非凸低秩约束正则项LR(Xi)和非凸稀疏约束正则项S(Xi),构建高光谱遥感图像的每个子块的恢复模型:
式(2)是利用高光谱遥感图像的低秩性和稀疏性的先验信息建立的图像恢复模型,其中,||·||表示矩阵二范数,α1和α2是用于平衡数据一致性项||Xi-Yi||2和正则项LR(Xi)和S(Xi)的两个正则项系数;Xi表示第i个子块的当前迭代次数的恢复图像;
(3)分别在LR(Xi)和S(Xi)中引入辅助变量R、Tg,利用交替迭代最小化算法分别迭代求解恢复图像X*和辅助变量R、Tg,过程如下:
首先利用含有辅助变量R、Tg的X*的二次函数最大化惩罚项,构造恢复模型的代价函数,然后利用迭代交替最小化方法最小化代价函数,求解得到X*、R、Tg;
(4)判断当前恢复图像是否满足收敛条件,若满足,进入步骤(5),否则进入步骤(3),继续利用交替迭代最小化算法迭代求解恢复图像X*和辅助变量R、Tg;
(6)计算X*的权重系数矩阵Weight;
本发明的有益效果是:
本发明的技术要点为:输入高光谱遥感图像;初始化权重系数矩阵、迭代次数和收敛阈值,初始化子图像尺寸和扫描步长,划分子块;建立图像恢复模型;引入辅助变量和正则项系数,利用最大值最小值算法迭代求解;判断恢复结果是否满足收敛条件;达到迭代次数获得符合要求的高光谱恢复图像,否则返回相应步骤继续迭代操作;计算权重系数矩阵,给每个子块分配合适的权重;还原高光谱遥感图像,得到最终恢复的高光谱遥感图像。本发明解决了现有的高光谱遥感图像恢复技术不能有效去除噪声的问题,提高了图像恢复质量。
利用本发明提出的基于非凸低秩稀疏约束的方法恢复高光谱遥感图像,更加充分地挖掘了高光谱遥感图像潜在的低秩性和稀疏性,图像的恢复质量比较高,去噪效果明显,并且有效地保留了图像的边缘细节,能够更好的满足高光谱遥感图像分类、目标检测等后续应用的需求。实验结果表明,与低秩矩阵恢复方法(LRMR),全变分正则化低秩矩阵分解方法(LRTV),噪声自适应迭代低秩矩阵逼近方法(NAILRMA)相比,本发明方法能在相同的噪声水平下获得更高的图像恢复质量,去噪效果明显,并且保留图像细节,有效实现了高光谱遥感图像恢复的目的。
附图说明
图1为本发明方法流程图(基于非凸低秩稀疏约束的高光谱遥感图像恢复方法流程图);
图2为EO-1Hyperion Australia数据集恢复效果图,第52波段恢复效果。
图3为HYDICE Washington DC Mall数据集恢复效果图,第1波段恢复效果。
具体实施方法
下面结合附图和实施实例对本发明进行详细说明。
如图1所示,本发明的具体实施步骤如下:
一种基于非凸低秩稀疏约束的高光谱遥感图像恢复方法,所述方法包含以下步骤:
(1)输入含有混合噪声的高光谱遥感图像Y,Y的大小为M×N×B,其中M和N分别代表各波段高光谱遥感图像的长和宽,B代表高光谱遥感图像波段数;
定义恢复后的图像X*,初始化变量X(0)=Y,X*是X(0)经过一定迭代次数后的恢复图像,初始化X*的权重系数矩阵Weight,初始化迭代次数和收敛阈值;
定义子块的尺寸p和扫描步长s,并将高光谱遥感图像Y按照预设的尺寸p和扫描步长s划分为若干个部分重叠的子块;
(2)对于选取的第i个子块Yi,将受噪声污染的高光谱遥感图像分解为其中为没有噪声污染的干净的高光谱遥感图像,Ji为包括高斯噪声、条带噪声、脉冲噪声的混合噪声;根据求解逆问题并引入非凸低秩约束正则项LR(Xi)和非凸稀疏约束正则项S(Xi),构建高光谱遥感图像的每个子块的恢复模型:
式(2)是利用高光谱遥感图像的低秩性和稀疏性的先验信息建立的图像恢复模型,其中,||·||表示矩阵二范数,α1和α2是用于平衡数据一致性项||Xi-Yi||2和正则项LR(Xi)和S(Xi)的两个正则项系数;Xi表示第i个子块的当前迭代次数的恢复图像;
(3)分别在LR(Xi)和S(Xi)中引入辅助变量R、Tg,利用交替迭代最小化算法分别迭代求解恢复图像X*和辅助变量R、Tg,过程如下:
首先利用含有辅助变量R、Tg的X*的二次函数最大化惩罚项,构造恢复模型的代价函数,然后利用迭代交替最小化方法最小化代价函数,求解得到X*、R、Tg;
(4)判断当前恢复图像是否满足收敛条件,若满足,进入步骤(5),否则进入步骤(3),继续利用交替迭代最小化算法迭代求解恢复图像X*和辅助变量R、Tg;
(6)计算X*的权重系数矩阵Weight;
在步骤一中,初始化X*的权重系数矩阵Weight为大小是M×N×B的零矩阵;定义子块的尺寸为p×p×B,扫描步长为s,其中1≤s≤p;按照预设的子块尺寸和扫描步长把高光谱遥感图像划分为若干个部分重叠的子块,共计个;p×p×B表示子块的长度×宽度×波段数。
步骤(2)中包含以下操作步骤:
步骤2-2,基于秩函数的非凸近似函数的Schatten-p范数给出LR(Xi)的表达式,具体形式如下:
其中,0<p1<1,σj是Xi的第j个奇异值;
S(Xi)为Xi的非凸稀疏约束正则项,表示为Xi的TV表达式的lp2范数,0<p2<1,具体形式如下:
其中,TV(Xi)为Xi的全变分(Total Variation,TV)表达式,具体形式如下:
其中,Dx,Dy,Dz分别表示Xi在x,y,z轴方向上的差分算子,x,y,z轴所在坐标系以图像Xi的第一个波段的第一个像素点为坐标原点,x和y表示高光谱遥感图像的空间轴,z轴表示光谱轴。
步骤(3)中包含以下操作步骤:
步骤3-1,采用最大值最小值的方法来求解式(2),引入两个辅助变量R和Tg,利用二次函数最大化惩罚项,得到下式:
其中,||·||F表示弗罗贝尼乌斯范数;0<p1<1,Tg中g=1,2,3时分别表示在x,y,z轴方向上T的值,β1、β2分别是LR(Xi)和S(Xi)的收敛系数;
将式(6)、(7)代入式(2),得到需要求解的优化问题如下:
式中,C为恢复优化问题的代价函数,当C取到最小值时,即得到变量(Xi,R,{Tg})的最优解,C具体表示为:
步骤3-2,求解上述优化问题,依次迭代更新变量(Xi,R,Tg),第n+1次迭代格式如下:
其中,α1,α2为固定大小的系数,β1,β2为随着迭代不断增大的系数,Tg,n+1、Rn+1、Xi,n+1分别表示变量Tg经过n+1次迭代后的值、变量Rn经过n+1次迭代后的值和变量Xi,n经过n+1次迭代后的值;
步骤3-3,求解辅助变量Tg,n+1:
首先求中间变量P:
Tg,n+1可以通过下式求解:
步骤3-4,求解辅助变量Rn+1:
首先对Xi,n奇异值分解:
Xi,n=UΣVT(15)
其中,U是奇异值分解的右奇异矩阵,V是左奇异矩阵,Σ是Xi,n的奇异值构成的对角矩阵,VT表示左奇异矩阵V的转置;
Rn+1可以通过下式求解:
其中,(·)+表示取矩阵大于0的部分;
步骤3-5,利用共轭梯度法来求解变量Xi经过n+1次迭代后的值Xi,n+1:
式(13)的欧拉-拉格朗日方程:
其中,I表示单位矩阵(恒等运算符);
式(17)的右边表达式的值在每次迭代的过程中都是常数,式(17)可视为一个线性方程组,可用共轭梯度法进行求解Xi,n+1。
步骤(4)中,所述收敛条件为:
其中,Threshold为收敛阈值,Cn+1和Cn分别代表代价函数C经过n+1次和n次迭代求解后的值。
在步骤(6)中,X*的权重系数矩阵Weight的求解过程:
高光谱遥感图像上像素点每参与迭代求解一次,即在权重系数矩阵Weight相应位置加1,迭代完成后即得到最终的权重系数矩阵。
本发明将通过以下实验进行有效性的说明
1.实验环境:
实验室用台式机参数:CPU为Inter(R)Core(TM)i5-4590,主频为3.30GHz,内存为8G,操作***为Win10 64位***,实验软件为MatlabR2014a。
2.实验结果及结果分析
本发明实验所用的图像数据为网络上公开的高光谱遥感图像数据集:Hyperspectral Digital Imagery Collection Experiment(HYDICE)Washington DCMall,the Earth Observing-1(EO-1)Hyperion Australia dataset。为了验证方法的有效性,本发明与近几年提出的高光谱遥感图像恢复方法进行对比:低秩矩阵恢复方法(Low-Rank Matrix Recovery,LRMR),全变分正则化低秩矩阵分解方法(Total-Variation-Regularized Low-Rank Matrix Approximation,LRTV),噪声自适应迭代低秩矩阵逼近方法(Noise-Adjusted Iterative Low-Rank Matrix Approximation,NAILRMA)。本发明的方法:非凸低秩和非凸稀疏高光谱恢复方法(Nonconvex Low-Rank and Nonconvex Total-Variation-Regularized Hyperspectral Image Restoration,NonLRTV)。
实例一:本实例选择遥感领域典型数据集——EO-1Hyperion Australia数据集,原始图像的尺寸是3858×256×242,这里因为篇幅限制,且去除可见光近红外和短波红外范围内的重叠波段后,选取尺寸为400×200×150,将其每个像素值减去最小像素值后,再除以最大与最小像素值之差,得到归一化到[0,1]后的数据。使用本发明提出的方法进行图像恢复处理,并用上述提到的LRMR、LRTV、NAILRMA进行对比实验。实验结果如图2所示,其中图2(a)为第52波段原始的高光谱遥感图像,图2(b)为LRMR方法效果,图2(c)为LRTV方法效果,图2(d)为NAILRMA方法效果,图2(e)为本发明方法效果。可以看出,图2(a)中含有大量噪声,而图2(b)、图2(c)和图2(d)都有明显的噪声,虽然去除部分噪声,但是去除效果有限。图2(e)中本发明的去噪效果明显,不仅很好地去除了噪声,而且很大程度上保留了丰富的细节。
实例二:本实例选择遥感领域典型数据集——Hyperspectral Digital ImageryCollection Experiment(HYDICE)Washington DC Mall数据集,原始图像的尺寸是1208×307×191,因为篇幅限制,选取图像尺寸为256×256×11。与实例一相同,对原始图像数据进行归一化处理,人为地对所有波段随机添加均值为25dB的高斯噪声。然后用本发明提出的方法进行恢复处理,并与其他方法进行对比。实验结果如图3所示,其中图3(a)为其中一张清晰的高光谱遥感图像数据,图3(b)为添加高斯噪声和条带噪声后的图像,图3(c)为LRMR方法效果(SNR=10.8786dB,MSSIM=0.66377),图3(d)为LRTV方法效果(SNR=10.6603dB,MSSIM=0.74259),图3(e)为NAILRMA方法效果(SNR=15.5521dB,MSSIM=0.83878),图3(f)为本发明方法效果(SNR=17.087dB,MSSIM=0.89488)。可以看出,图3(c)和图3(e)仍有明显噪声,图3(d)虽然没有明显噪声,但是过于模糊,丢失了许多细节。图3(f)中本发明的去噪效果明显,能够在有效去除大部分噪声的同时,保留绝大部分的细节,恢复效果最好,并且评价指标SNR、MSSIM的值较前三种方法均有较大程度的提高。
在相同的噪声水平下,利用本发明中提出的方法恢复出的高光谱遥感图像具有较高的信噪比,并且与清晰的未加噪声的高光谱遥感图像的平均结构相似度最高。可见,相比于其他三种现已提出的方法,本发明所提方法可以获得更好的恢复结果,与为加噪声的原数据图像最为接近,恢复效果最好。
Claims (3)
1.一种基于非凸低秩稀疏约束的高光谱遥感图像恢复方法,其特征在于,所述方法包含以下步骤:
(1)输入含有混合噪声的高光谱遥感图像Y,Y的大小为M×N×B,其中M和N分别代表各波段高光谱遥感图像的长和宽,B代表高光谱遥感图像波段数;
定义恢复后的图像X*,初始化变量X(0)=Y,X*是X(0)经过一定迭代次数后的恢复图像,初始化X*的权重系数矩阵Weight,初始化迭代次数和收敛阈值;
定义子块的尺寸p和扫描步长s,并将高光谱遥感图像Y按照预设的尺寸p和扫描步长s划分为若干个部分重叠的子块;
(2)对于选取的第i个子块Yi,将受噪声污染的高光谱遥感图像分解为其中为没有噪声污染的干净的高光谱遥感图像,Ji为包括高斯噪声、条带噪声、脉冲噪声的混合噪声;根据求解逆问题并引入非凸低秩约束正则项LR(Xi)和非凸稀疏约束正则项S(Xi),构建高光谱遥感图像的每个子块的恢复模型:
式(2)是利用高光谱遥感图像的低秩性和稀疏性的先验信息建立的图像恢复模型,其中,||·||表示矩阵二范数,α1和α2是用于平衡数据一致性项||Xi-Yi||2和正则项LR(Xi)和S(Xi)的两个正则项系数;Xi表示第i个子块的当前迭代次数的恢复图像;
(3)分别在LR(Xi)和S(Xi)中引入辅助变量R、Tg,利用交替迭代最小化算法分别迭代求解恢复图像X*和辅助变量R、Tg,过程如下:
首先利用含有辅助变量R、Tg的X*的二次函数最大化惩罚项,构造恢复模型的代价函数,然后利用迭代交替最小化方法最小化代价函数,求解得到X*、R、Tg;
(4)判断当前恢复图像是否满足收敛条件,若满足,进入步骤(5),否则进入步骤(3),继续利用交替迭代最小化算法迭代求解恢复图像X*和辅助变量R、Tg;
(6)计算X*的权重系数矩阵Weight;
步骤(1)中,初始化X*的权重系数矩阵Weight为大小是M×N×B的零矩阵;定义子块的尺寸为p×p×B,扫描步长为s,其中1≤s≤p;按照预设的子块尺寸和扫描步长把高光谱遥感图像划分为若干个部分重叠的子块,共计个;p×p×B表示子块的长度×宽度×波段数;
步骤(2)中包含以下操作步骤:
步骤2-2,基于秩函数的非凸近似函数的Schatten-p范数给出LR(Xi)的表达式,具体形式如下:
其中,0<p1<1,σj是Xi的第j个奇异值;
其中,TV(Xi)为Xi的全变分(Total Variation,TV)表达式,具体形式如下:
其中,Dx,Dy,Dz分别表示Xi在x,y,z轴方向上的差分算子,x,y,z轴所在坐标系以图像Xi的第一个波段的第一个像素点为坐标原点,x和y表示高光谱遥感图像的空间轴,z轴表示光谱轴;
步骤(3)中包含以下操作步骤:
步骤3-1,采用最大值最小值的方法来求解式(2),引入两个辅助变量R和Tg,利用二次函数最大化惩罚项,得到下式:
其中,||·||F表示弗罗贝尼乌斯范数;0<p1<1,Tg中g=1,2,3时分别表示在x,y,z轴方向上T的值,β1、β2分别是LR(Xi)和S(Xi)的收敛系数;
将式(6)、(7)代入式(2),得到需要求解的优化问题如下:
式中,C为恢复优化问题的代价函数,当C取到最小值时,即得到变量(Xi,R,{Tg})的最优解,C具体表示为:
步骤3-2,求解上述优化问题,依次迭代更新变量(Xi,R,Tg),第n+1次迭代格式如下:
其中,β1,β2为随着迭代不断增大的系数,Tg,n+1、Rn+1、Xi,n+1分别表示变量Tg经过n+1次迭代后的值、变量Rn经过n+1次迭代后的值和变量Xi,n经过n+1次迭代后的值;
步骤3-3,求解辅助变量Tg,n+1:
首先求中间变量P:
Tg,n+1通过下式求解:
步骤3-4,求解辅助变量Rn+1:
首先对Xi,n奇异值分解:
Xi,n=UΣVT (15)
其中,U是奇异值分解的右奇异矩阵,V是左奇异矩阵,Σ是Xi,n的奇异值构成的对角矩阵,VT表示左奇异矩阵V的转置;
Rn+1通过下式求解:
其中,(·)+表示取矩阵大于0的部分;
步骤3-5,利用共轭梯度法来求解变量Xi经过n+1次迭代后的值Xi,n+1:
式(13)的欧拉-拉格朗日方程:
其中,I表示单位矩阵;
式(17)的右边表达式的值在每次迭代的过程中都是常数,式(17)视为一个线性方程组,用共轭梯度法进行求解Xi,n+1。
3.根据权利要求2所述的一种基于非凸低秩稀疏约束的高光谱遥感图像恢复方法,其特征在于,
在步骤(6)中,X*的权重系数矩阵Weight的求解过程:
高光谱遥感图像上像素点每参与迭代求解一次,即在权重系数矩阵Weight相应位置加1,迭代完成后即得到最终的权重系数矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811018402.9A CN109102477B (zh) | 2018-08-31 | 2018-08-31 | 一种基于非凸低秩稀疏约束的高光谱遥感图像恢复方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811018402.9A CN109102477B (zh) | 2018-08-31 | 2018-08-31 | 一种基于非凸低秩稀疏约束的高光谱遥感图像恢复方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109102477A CN109102477A (zh) | 2018-12-28 |
CN109102477B true CN109102477B (zh) | 2021-08-24 |
Family
ID=64864823
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811018402.9A Expired - Fee Related CN109102477B (zh) | 2018-08-31 | 2018-08-31 | 一种基于非凸低秩稀疏约束的高光谱遥感图像恢复方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109102477B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109978783A (zh) * | 2019-03-19 | 2019-07-05 | 上海交通大学 | 一种彩色图像修复方法 |
CN110335201A (zh) * | 2019-03-27 | 2019-10-15 | 浙江工业大学 | 结合Moreau增强TV和局部低秩矩阵恢复的高光谱图像去噪方法 |
CN110310234A (zh) * | 2019-05-17 | 2019-10-08 | 浙江工业大学 | 基于低秩逼近和香农全差分的人脸图像去噪方法及装置 |
CN110120026B (zh) * | 2019-05-23 | 2022-04-05 | 东北大学秦皇岛分校 | 基于Schatten Capped p范数的数据恢复方法 |
CN110443169B (zh) * | 2019-07-24 | 2022-10-21 | 广东工业大学 | 一种边缘保留判别分析的人脸识别方法 |
CN110728641A (zh) * | 2019-10-12 | 2020-01-24 | 浙江工业大学 | 一种遥感图像脉冲噪声去除方法及装置 |
CN111028172A (zh) * | 2019-12-10 | 2020-04-17 | 浙江工业大学 | 基于无参数的非凸低秩矩阵逼近的高光谱图像去噪方法 |
CN111260571B (zh) * | 2020-01-11 | 2022-09-23 | 浙江工业大学 | 一种基于非凸低秩低梯度的深度图像修复方法 |
CN111292266B (zh) * | 2020-01-22 | 2022-07-05 | 武汉大学 | 基于双低秩矩阵分解的gf-5遥感影像混合噪声去除方法 |
CN111598797B (zh) * | 2020-04-27 | 2023-08-29 | 浙江工业大学 | 一种基于非凸张量范数和全变分的高光谱图像恢复方法 |
CN111738937B (zh) * | 2020-05-28 | 2022-08-12 | 南京邮电大学 | 一种基于张量分解的总变分遥感图像去条带方法 |
CN111951183B (zh) * | 2020-07-14 | 2024-03-29 | 浙江工业大学 | 基于近端交替惩罚算法的低秩全变分高光谱图像修复方法 |
CN111985086B (zh) * | 2020-07-24 | 2024-04-09 | 西安理工大学 | 一种融合先验信息和稀疏约束的社团检测方法 |
US11788888B2 (en) * | 2020-09-08 | 2023-10-17 | Shenzhen Hypernano Optics Technology Co., Ltd | Method and device for restoring and reconstructing a spectrum of a light source based on a hyperspectral image |
CN112634167B (zh) * | 2020-12-29 | 2023-09-01 | 南京理工大学 | 一种全变差协同范数约束迭代投影的高光谱图像滤波方法 |
CN112950500B (zh) * | 2021-02-25 | 2022-08-26 | 桂林电子科技大学 | 一种基于边缘检测低秩全变分模型的高光谱去噪方法 |
CN112862729B (zh) * | 2021-03-30 | 2022-09-06 | 湖北地信科技集团股份有限公司 | 一种基于特征曲线指导的遥感图像去噪方法 |
CN113598716B (zh) * | 2021-08-24 | 2022-04-15 | 南方医科大学 | 一种光声与磁共振联合成像方法 |
CN113706418B (zh) * | 2021-09-07 | 2024-07-23 | 哈尔滨工业大学 | 一种基于光谱分离的长波红外遥感图像恢复方法 |
CN116012264B (zh) * | 2023-03-27 | 2023-06-13 | 山东省工业技术研究院 | 一种基于稀疏约束的图像恢复方法 |
CN116433534B (zh) * | 2023-06-09 | 2023-08-22 | 四川工程职业技术学院 | 一种高光谱图像修复方法、装置、存储介质及电子设备 |
CN117478365B (zh) * | 2023-10-17 | 2024-05-10 | 南京邮电大学 | 一种计及攻击的能源互联网用电数据软恢复方法、设备及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105551000A (zh) * | 2015-12-11 | 2016-05-04 | 中国科学院遥感与数字地球研究所 | 基于参考影像结构约束和非凸低秩约束的遥感图像重建方法 |
CN106408530A (zh) * | 2016-09-07 | 2017-02-15 | 厦门大学 | 基于稀疏与低秩矩阵逼近的高光谱图像恢复方法 |
CN106709881A (zh) * | 2016-12-14 | 2017-05-24 | 上海增容数据科技有限公司 | 一种基于非凸低秩矩阵分解的高光谱图像去噪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10235600B2 (en) * | 2015-06-22 | 2019-03-19 | The Johns Hopkins University | System and method for structured low-rank matrix factorization: optimality, algorithm, and applications to image processing |
-
2018
- 2018-08-31 CN CN201811018402.9A patent/CN109102477B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105551000A (zh) * | 2015-12-11 | 2016-05-04 | 中国科学院遥感与数字地球研究所 | 基于参考影像结构约束和非凸低秩约束的遥感图像重建方法 |
CN106408530A (zh) * | 2016-09-07 | 2017-02-15 | 厦门大学 | 基于稀疏与低秩矩阵逼近的高光谱图像恢复方法 |
CN106709881A (zh) * | 2016-12-14 | 2017-05-24 | 上海增容数据科技有限公司 | 一种基于非凸低秩矩阵分解的高光谱图像去噪方法 |
Non-Patent Citations (3)
Title |
---|
Hyperspectral Image Restoration via Iteratively Regularized Weighted Schatten p-Norm Minimization;Yuan Xie等;《IEEE Transactions on Geoscience and Remote Sensing》;20160415;全文 * |
基于低秩和稀疏性先验知识的压缩感知图像重构;陈长伟;《计算机应用研究》;20170331;第34卷(第3期);全文 * |
非凸稀疏低秩约束的高光谱解混方法;孔繁锵等;《西安电子科技大学学报》;20161231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109102477A (zh) | 2018-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109102477B (zh) | 一种基于非凸低秩稀疏约束的高光谱遥感图像恢复方法 | |
CN109741256B (zh) | 基于稀疏表示和深度学习的图像超分辨率重建方法 | |
CN106709881B (zh) | 一种基于非凸低秩矩阵分解的高光谱图像去噪方法 | |
CN105513026B (zh) | 一种基于图像非局部相似的压缩感知重构方法 | |
Wang et al. | Low rank constraint and spatial spectral total variation for hyperspectral image mixed denoising | |
CN108133465B (zh) | 基于空谱加权tv的非凸低秩松弛的高光谱图像恢复方法 | |
CN112069919A (zh) | 基于非凸低秩矩阵逼近和全变分正则化的高光谱图像去噪方法 | |
CN111127354B (zh) | 一种基于多尺度字典学习的单图像去雨方法 | |
Zhang et al. | Kernel Wiener filtering model with low-rank approximation for image denoising | |
CN112734763B (zh) | 基于卷积与k-svd字典联合稀疏编码的图像分解方法 | |
Huang et al. | Spatial-spectral weighted nuclear norm minimization for hyperspectral image denoising | |
Yang et al. | An improved weighted nuclear norm minimization method for image denoising | |
Han et al. | Edge-preserving filtering-based dehazing for remote sensing images | |
Chen et al. | Hyperspectral image denoising via texture-preserved total variation regularizer | |
Wang et al. | Weighted Schatten p-norm minimization for impulse noise removal with TV regularization and its application to medical images | |
Fu et al. | Hyperspectral image denoising via robust subspace estimation and group sparsity constraint | |
CN112837220B (zh) | 一种提高红外图像分辨率的方法及其用途 | |
CN113421198A (zh) | 一种基于子空间的非局部低秩张量分解的高光谱图像去噪方法 | |
Niresi et al. | Robust hyperspectral inpainting via low-rank regularized untrained convolutional neural network | |
Cui et al. | Compressed sensing based on L1 and TGV regularization for low-light-level images denoising | |
Chen et al. | A guidable nonlocal low-rank approximation model for hyperspectral image denoising | |
Luo et al. | Low-rank decomposition on transformed feature maps domain for image denoising | |
CN112862729B (zh) | 一种基于特征曲线指导的遥感图像去噪方法 | |
CN111598797B (zh) | 一种基于非凸张量范数和全变分的高光谱图像恢复方法 | |
CN114359064A (zh) | 基于双重梯度约束的高光谱图像恢复方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210824 |