一种基于自适应秩校正的高光谱图像去噪方法、系统及介质与流程

本发明涉及高光谱图像去噪技术，具体涉及一种基于自适应秩校正的高光谱图像去噪方法、系统及介质。

背景技术：

高光谱图像(hsi)包含丰富的光谱信息，在环境监测、城市规划、地质勘探和农林业生产等领域应用广泛。在成像和传输过程中，因受到多种复杂因素影响，获取的高光谱图像包含各种类型的混合噪声，如高斯噪声、脉冲噪声、条带、死线等。这些复杂的混合噪声不仅影响图像质量，且会对后续应用造成很大困扰。因此，提出一种有效的高光谱图像去噪方法对于高光谱图像处理分析和后续应用意义重大。

传统2d图像去噪方法，如非局部均值、k-svd、块匹配三维滤波(bm3d)等，采用逐波段处理方式对高光谱图像的每个波段逐一进行去噪。然而，逐波段处理方式忽略了光谱相关性，去噪效果不佳。为解决该问题，研究人员综合利用空间信息和光谱信息进行去噪。卢等人提出一种光谱-空间自适应稀疏表示高光谱图像去噪方法。另外，通过将高光谱图像视为3d张量，采用基于张量分解的方法对高光谱图像进行去噪，如多维维纳滤波与tucker分解结合方法，parafac方法等。通过充分利用空间相似性和光谱相关性，彭等人提出了一种可分解非局部张量字典学习(tdl)方法。总而言之，上述方法主要用于去除高斯噪声及混合高斯泊松噪声。

但是由于实际中，高光谱图像通常受到高斯噪声、脉冲噪声、死线和条带等多种噪声污染。因而，开发一种有效的高光谱图像去噪方法用以消除多种噪声就显得十分重要。由于干净高光谱图像hsi具有低秩特性，故而研究人员通常将此问题转化为低秩矩阵恢复(lrmr)问题。直接的矩阵秩最小化问题是一个非确定性多项式难(np-hard)问题，因此现有某些方法采用核范数作为凸代理来逼近矩阵秩，或采用秩逼近矩阵秩。非凸go分解(godec)算法被应用来更好地去除混合噪声。此外，通过利用非独立同分布(non-i.i.d.)混合高斯噪声(nmogs)假设，结合低秩矩阵分解模型，现有的某些方法提出非独立同分布混合高斯模型，该方法具有去除多种噪声的能力。核范数会过度惩罚大的奇异值，从而使获得的解偏离原始解。而基于秩逼近的低秩矩阵恢复方法，秩信息固定且需要提前定义，这在现实问题中不具有灵活性。

综上所述，现有大多数基于低秩矩阵恢复的高光谱去噪方法利用核范数逼近矩阵秩导致过度惩罚大的奇异值和秩r信息固定且需提前定义使得方法不灵活等问题。

技术实现要素：

本发明要解决的技术问题：为解决现有大多数基于低秩矩阵恢复的高光谱去噪方法利用核范数逼近矩阵秩导致过度惩罚大的奇异值和秩r信息固定且需提前定义使得方法不灵活等问题，提供一种基于自适应秩校正(arc)的高光谱图像去噪方法、系统及介质，本发明能够去除高光谱图像中的高斯噪声、脉冲噪声、死线和条带等复杂噪声。

为了解决上述技术问题，本发明采用的技术方案为：

一种基于自适应秩校正的高光谱图像去噪方法，实施步骤包括：

1)对退化高光谱图像矩阵化后采用非精确增广拉格朗日乘子算法进行初始估计，得到初始估计的干净高光谱图像

2)对退化高光谱图像初始估计的干净高光谱图像分别进行分块和矩阵化处理，分别得到矩阵化分块退化高光谱图像y、矩阵化分块初始估计的干净高光谱图像x；

3)针对各个矩阵化分块退化高光谱图像y、矩阵化分块初始估计的干净高光谱图像x建立基于自适应秩校正的高光谱图像去噪模型，所述基于自适应秩校正的高光谱图像去噪模型通过自适应秩校正惩罚函数以实现自适应地抵消核范数对大的奇异值的惩罚，且基于l2,1范数消除稀疏噪声；对所述基于自适应秩校正的高光谱图像去噪模型采用交替方向乘子算法进行求解得到分块干净图像x；

4)将所有的分块干净图像x组合获得最终估计的干净高光谱图像。

可选地，步骤1)的详细步骤包括：

1.1)将退化高光谱图像矩阵化为为退化高光谱图像其中表示实数域，w、h和b分别表示退化高光谱图像的宽度、高度和光谱带数量；

1.2)采用非精确增广拉格朗日乘子算法ialm对退化高光谱图像进行初始估计得到初始估计的干净图像

1.3)通过将初始估计的干净图像的每一列重构为大小为w×h的矩阵，得到初始估计的干净高光谱图像

可选地，步骤2)的详细步骤包括：对退化高光谱图像，首先按照步长p以像素(u,v)为中心取大小为q×q×b的子立方体，然后将每个子立方体的每个波段转化为1d列向量来构造矩阵从而得到矩阵化分块退化高光谱图像y；对初始估计的干净高光谱图像首先按照步长p以像素(u,v)为中心取大小为q×q×b的子立方体，然后将每个子立方体的每个波段转化为1d列向量来构造矩阵从而得到矩阵化分块初始估计的干净高光谱图像x；其中，1≤u≤w,1≤v≤h，其中w、h和b分别表示退化高光谱图像的宽度、高度和光谱带数量，q表示子立方体块大小。

可选地，步骤3)中建立的基于自适应秩校正的高光谱图像去噪模型如下式所示；

上式中，||x||*表示矩阵x核范数，x为分块干净图像，<f(x),x>表示f(x)与x的欧式内积，f(x)表示谱算子，x表示矩阵化分块初始估计的干净高光谱图像，||s||2,1表示矩阵sl2,1范数，λ表示规则化参数，表示矩阵弗罗贝尼乌斯范数，y表示矩阵化分块退化高光谱图像，s表示矩阵化分块稀疏噪声，σ为与高斯噪声n标准差相关的常数。

可选地，步骤4)中将所有的分块干净图像x组合获得最终估计的干净高光谱图像时，还包括针对分块干净图像x的重叠区域取平均值的方式处理重叠区域像素的步骤。

此外，本发明还提供一种基于自适应秩校正的高光谱图像去噪系统，包括计算机设备，该计算机设备被编程或配置以执行所述基于自适应秩校正的高光谱图像去噪方法的步骤。

此外，本发明还提供一种基于自适应秩校正的高光谱图像去噪系统，包括计算机设备，其特征在于，该计算机设备的存储介质上存储有被编程或配置以执行所述基于自适应秩校正的高光谱图像去噪方法的计算机程序。

此外，本发明还提供一种计算机可读存储介质，该计算机可读存储介质上存储有被编程或配置以执行所述基于自适应秩校正的高光谱图像去噪方法的计算机程序。

此外，本发明还提供一种基于自适应秩校正的高光谱图像去噪系统，包括：

第一程序单元，用于对退化高光谱图像矩阵化后采用非精确增广拉格朗日乘子算法进行初始估计，得到初始估计的干净高光谱图像

第二程序单元，用于对退化高光谱图像初始估计的干净高光谱图像分别进行分块和矩阵化处理，分别得到矩阵化分块退化高光谱图像y、矩阵化分块初始估计的干净高光谱图像x；

第三程序单元，用于针对各个矩阵化分块退化高光谱图像y、矩阵化分块初始估计的干净高光谱图像x建立基于自适应秩校正的高光谱图像去噪模型，所述基于自适应秩校正的高光谱图像去噪模型通过自适应秩校正惩罚函数以实现自适应地抵消核范数对大的奇异值的惩罚，且基于l2,1范数消除稀疏噪声；对所述基于自适应秩校正的高光谱图像去噪模型采用交替方向乘子算法进行求解，得到分块干净图像x；

第四程序单元，用于将所有的分块干净图像x组合获得最终估计的干净高光谱图像。

和现有技术相比，本发明具有下述优点：本发明基于自适应秩校正的高光谱图像去噪方法提出的基于自适应秩校正的高光谱图像去噪模型通过自适应秩校正惩罚函数以实现自适应地抵消核范数对大的奇异值的惩罚，自适应秩校正惩罚函数不需要提前定义秩r信息，能够自适应地抵消核范数对大的奇异值的惩罚，从而能对秩进行自适应校正，准确捕捉干净hsi本质低秩信息；本发明基于自适应秩校正的高光谱图像去噪模型通过基于l2,1范数消除稀疏噪声，如死线和条带噪声；因此，本发明能够有效地去除混合噪声，还能有效保留高光谱图像的结构、细节信息，具有非常重要的应用价值。

附图说明

图1为本发明实施例的自适应秩校正的高光谱图像去噪方法流程图。

图2为干净高光谱图像的低秩特性示意图。

图3为本发明实施例中当τ和ε不同取值时的标量函数φ(t)的函数图像。

图4为流行的5种现有去噪方法与本实施例方法的第48个波段去噪结果对比。

图5为流行的5种现有去噪方法与本实施例方法的第57个波段去噪结果对比。

图6为5种现有去噪方法与本实施例方法对频带数为108的去噪结果对比图。

图7为5种现有去噪方法与本实施例方法对频带数为150的去噪结果对比图。

图8为5种现有去噪方法与本实施例方法对频带数为103的去噪结果对比图。

图9为5种现有去噪方法与本实施例方法对频带数为138的去噪结果对比图。

具体实施方式

如图1所示，本实施例基于自适应秩校正的高光谱图像去噪方法的实施步骤包括：

1)对退化高光谱图像矩阵化后采用非精确增广拉格朗日乘子算法(ialm)进行初始估计，得到初始估计的干净高光谱图像

2)对退化高光谱图像初始估计的干净高光谱图像分别进行分块和矩阵化处理，分别得到矩阵化分块退化高光谱图像y、矩阵化分块初始估计的干净高光谱图像x；

3)针对各个矩阵化分块退化高光谱图像y、矩阵化分块初始估计的干净高光谱图像x建立基于自适应秩校正的高光谱图像去噪模型，基于自适应秩校正的高光谱图像去噪模型通过自适应秩校正惩罚函数以实现自适应地抵消核范数对大的奇异值的惩罚，且基于l2,1范数消除稀疏噪声；对基于自适应秩校正的高光谱图像去噪模型采用交替方向乘子算法(admm)进行求解得到分块干净图像x；

4)将所有的分块干净图像x组合获得最终估计的干净高光谱图像。

步骤1)采用非精确增广拉格朗日乘子算法(ialm)进行初始估计的目的在于得到初始估计的干净高光谱图像；本实施例中，步骤1)的详细步骤包括：

1.1)将退化高光谱图像矩阵化为为退化高光谱图像其中表示实数域，w、h和b分别表示退化高光谱图像的宽度、高度和光谱带数量；

1.2)采用非精确增广拉格朗日乘子算法ialm对退化高光谱图像进行初始估计得到初始估计的干净图像

1.3)通过将初始估计的干净图像的每一列重构为大小为w×h的矩阵，得到初始估计的干净高光谱图像

步骤2)的目的在于获取分块高光谱图像的矩阵表示形式，将去噪问题转化为数值优化问题。对于任一退化高光谱图像该退化图像一般可认为是由干净图像、稀疏噪声和高斯噪声三部分构成，即这里与有相同的大小，分别代表干净hsi、稀疏噪声、高斯噪声；w、h和b分别表示hsi的宽度、高度和光谱带数量。在退化高光谱图像中，对于以像素(u,v),(1≤u≤w,1≤v≤h)为中心，大小为q×q×b的子立方体，可通过将子立方体的每个波段转化为1d列向量来构造矩阵在初始估计的干净高光谱图像中，可通过同样的方式在相同位置构造矩阵用于后续秩校正项的构建。相应的高光谱图像退化模型可转换为矩阵形式如式(1)所示：

y＝x+s+n(1)

式(1)中，y表示矩阵化分块退化高光谱图像，x代表分块干净图像，s和n分别代表矩阵化分块稀疏噪声矩阵和矩阵化分块高斯噪声矩阵，矩阵x,s,n和y具有相同的大小。

本实施例中，步骤2)的详细步骤包括：对退化高光谱图像首先按照步长p以像素(u,v)为中心取大小为q×q×b的子立方体，然后将每个子立方体的每个波段转化为1d列向量来构造矩阵得到矩阵化分块退化高光谱图像y；对初始估计的干净高光谱图像首先按照步长p以像素(u,v)为中心取大小为q×q×b的子立方体，然后将每个子立方体的每个波段转化为1d列向量来构造矩阵从而得到矩阵化分块初始估计的干净高光谱图像x；其中，1≤u≤w,1≤v≤h，其中w、h和b分别表示退化高光谱图像的宽度、高度和光谱带数量，q表示子立方体块大小。

步骤3)基于自适应秩校正的高光谱图像去噪模型通过自适应秩校正惩罚函数以实现自适应地抵消核范数对大的奇异值的惩罚，且基于l2,1范数消除稀疏噪声，包括脉冲噪声、死线和条带等的本质稀疏结构。本实施例中，步骤3)中建立的基于自适应秩校正的高光谱图像去噪模型如式(2)所示；

式(2)中，||x||*表示矩阵核范数，x为分块干净图像，<f(x),x>表示f(x)与x的欧式内积，f(·)表示谱算子函数，x表示矩阵化分块初始估计的干净高光谱图像，||s||2，1表示矩阵l2,1范数，λ表示规则化参数，表示矩阵弗罗贝尼乌斯范数，y表示矩阵化分块退化高光谱图像，s表示矩阵化分块稀疏噪声，σ为与高斯噪声n标准差相关的常数。传统的基于秩-r逼近(rank(x)≤r)的低秩矩阵恢复方法中，秩r信息固定且需要提前定义r的上确界，使得此类在现实问题中不具有灵活性。本实施例基于自适应秩校正的高光谱图像去噪模型利用自适应秩校正惩罚函数||x||*-<f(x),x>来有效逼近矩阵秩，其中，||x||*代表矩阵x的核范数，定义为矩阵奇异值之和，x是对初始估计的干净高光谱图像进行分块操作并矩阵化得到，f(·)为谱算子函数。该惩罚函数能进行自适应秩校正，且能够自适应地抵消核范数对大的奇异值的惩罚。l0范数最小化为np难问题，将l0范数松弛为l2,1范数来描述稀疏噪声，包括脉冲噪声、死线和条带等的本质稀疏结构，其中l2,1范数为矩阵每一行的l2范数(向量各元素的平方和然后求平方根)之和。l2,1范数对于处理特定结构的稀疏噪声，如死线和条带噪声，非常有效。通过转化，得到基于自适应秩校正的高光谱图像去噪模型式(2)并获得分块后干净图像x，其中λ为一个正则化参数，参见式(2)可知，本实施例基于自适应秩校正的高光谱图像去噪模型的自适应秩校正惩罚函数为核范数减去基于初始估计的线性秩校正项，来有效逼近矩阵秩，该惩罚函数不需要提前定义秩r信息，且能够自适应地抵消核范数对较大奇异值的惩罚，具有自适应秩校正能力。同时，利用l2,1范数来逼近l0范数描述稀疏噪声(包括脉冲噪声、死线和条带等)的本质稀疏结构。利用交替方向乘子算法(admm)对得到的凸优化模型进行有效求解。与其它高光谱图像去噪方法相比可知，本实施例所提出的自适应秩校正方法不仅能同时去除高斯噪声、脉冲噪声、死线和条带等复杂噪声，还能有效保留高光谱图像的空间和光谱信息。

图2为干净高光谱图像的低秩先验特性的说明。图2显示干净图像的大部分奇异值都趋近于0，体现了干净图像具有低秩的先验特性。步骤3)中采用提出的自适应秩校正(arc)方法，对分块高光谱图像去噪，获得分块后干净图像。l0范数最小化问题是np难问题，且基于秩-r逼近(rank(x)≤r)的低秩矩阵恢复方法，秩r信息固定且需要提前定义r的上确界，使得此类在现实问题中不具有灵活性。利用自适应秩校正惩罚函数来有效逼近矩阵秩，利用l0范数松弛为l2,1范数来描述稀疏噪声，最终得到分块干净图像x。对于任何矩阵约束rank(x)≤r等价于式(3)所示函数表达式：

0＝σr+1(x)+…+σn(x)＝||x||*-||x||(r)(3)

式(3)中，σr+1(x)～σn(x)分别表示分块干净图像x的第r+1个奇异值～第n个奇异值，||x||*表示矩阵核范数，||x||(r)＝σ1(x)+…+σr(x)，n＝min(q²,b)。

基于当前的迭代值x^t，目标方程(3)获得下一步的迭代值x^t+1，通过解如下凸优化方程如式(4)所示；

式(4)中，x^t+1表示下一步的迭代值，||x||*表示矩阵核范数，x^t表示当前的迭代值，g^t是凸函数||x||(r)在x^t这一点的次梯度。

在实际应用中，预先确定r值是一个挑战性的工作，尤其对那些基于分块的hsi去噪方法来说。另外，对于带噪声的低秩矩阵逼近问题，没有必要通过精确约束秩来得到目标解。因此，本发明提出一种惩罚函数来有效逼近秩函数，自适应进行秩校正。如果初始值与真实矩阵偏差不大，就会在一定程度上包含真实矩阵奇异值或者秩的一些信息。基于这个发现，给定一个与y具有相同大小的矩阵化分块初始估计的干净高光谱图像x，参见式(2)，本实施例提出的自适应秩校正惩罚函数定义为式(5)：

||x||*-<f(x),x〉(5)

式(5)中各符号定义参见式(2)。

这里f:是一个谱算子函数，定义见式(6)。

f(x)＝udiag(f(σ(x)))v^t(6)

式(6)中，f(·)表示谱算子函数,σ(x)表示矩阵x的奇异值，f(·)表示对称函数，diag用于构造一个对角矩阵，对角线上的元素等于f(σ(x))，u和v的列分别表示奇异值中的左、右奇异向量。

对称函数f:定义为式(7)：

式(7)中，fi(x)表示对称函数f在xi处的函数值，x表示任一向量x，xi表示x的第i个数值，||x||∞表示x的最大范数，φ表示标量函数。其中||x||∞＝maxi|xi|，由σ(x)降序排列可得||σ(x)||∞＝σ1(x)，σ(x)表示矩阵x的奇异值。

对于参数τ＞0和ε＞0，则有标量函数φ:形式如式(8)：

式(8)中，φ(t)为标量函数φ，ε^τ表示ε的τ次方，|t|^τ表示t绝对值的τ次方。

如果令f＝0，式(5)即为核范数。核范数对每个奇异值进行同等放缩。对图像而言，大的奇异值往往对应更多的结构信息，小的奇异值则对应噪声。因而，我们希望对大的奇异值放缩少或者对小的奇异值值放缩多。给定则有式(9)：

式(9)中，左侧各符号的定义参见式(2)，右侧各符号中，σi(x)表示x的第i个奇异值，φ(ti)表示函数φ在ti处的函数值。此处n＝min(q²,b)，q表示矩阵x的行数，b表示矩阵x的列数；表示矩阵x的第i个奇异值，表示矩阵x奇异值向量的最大范数。

在不同的τ＞0和ε＞0(τ表示输入参数，ε表示输入参数)的情况下：当自变量t在[0，1]上取值，图3显示了τ、ε在不同取值情况下标量函数φ的各种不同形状。当τ和ε取大于零的不同正数时，针对t∈[0,1]，φ(t)值的变化；(a)当ε＝0.1，τ取大于零的不同正数时，φ(t)值的变化；(b)当τ＝2，ε取大于零的不同正数时，φ(t)值的变化。可以看出，参数τ＞0主要控制标量函数φ的形状。对于任何满0＜ε＜1，如果0＜t＜ε，标量函数φ的函数值φ(t)趋于0，而如果ε＜t≤1，当τ→∞，则标量函数φ的函数值φ(t)趋于1。因此，如果ti接近0，有σi(x)-φ(ti)σi(x)→σi(x)。如果ti接近1，则σi(x)-φ(ti)σi(x)→0。因此，本发明提出的惩罚函数能够自适应地抵消核范数对大的奇异值的惩罚。

此外，本实施例基于自适应秩校正的高光谱图像去噪模型还将l0范数松弛为l2,1范数，来描述稀疏噪声，包括脉冲噪声、死线和条带等的本质稀疏结构。

对于任意的矩阵矩阵s的l2,1范数定义为：

式(10)中，||s||2,1表示矩阵s的l2,1范数，b表示退化高光谱图像的光谱带数量，q表示子立方体块大小，sij表示矩阵s第(i,j)个数值。l2,1范数对于处理特定结构的稀疏噪声，如死线和条带噪声，非常有效。基于上述推导过程，可得到步骤3)中建立的基于自适应秩校正的高光谱图像去噪模型如式(2)所示。

步骤3)中对基于自适应秩校正的高光谱图像去噪模型采用交替方向乘子算法(admm)进行求解，交替方向乘子算法(admm)为现有求解方法，目的是求解目标最小化方程，得到全局最优解。交替方向乘子算法(admm)利用目标函数的可分离性，将原问题分解为若干个子问题，然后交替地进行求解。交替方向乘子算法(admm)在处理大规模问题和多目标优化问题方面非常有效。基于增广拉格朗日函数的目标函数可表示为：

式(11)中，λ表示拉格朗日乘子，μ表示惩罚参数，其余各符号定义参见式(2)。目标函数(11)可以通过优化其中一个变量，固定其余变量，迭代求解。特别地，在第t+1迭代，目标函数(11)中的变量更新如下：

λ^k+1＝λ^k+μ^k(y-x^k+1-s^k+1)(14)

μ^k+1＝min(t*μ^k,μmax)(15)

式(12)～(15)中，x^k+1表示分块干净图像x在第k+1次迭代中的值，表示自变量为x,s^k,λ^k,μ^k的增广拉格朗日函数，s^k+1表示矩阵化分块稀疏噪声s在第k+1次迭代中的值，表示自变量为x^k+1,s,λ^k,μ^k的增广拉格朗日函数，λ^k+1表示拉格朗日乘子λ在第k+1次迭代中的值，λ^k表示拉格朗日乘子λ在第k次迭代中的值，μ^k表示惩罚参数u在第k次迭代中的值，y表示矩阵化分块退化高光谱图像，μ^k+1表示惩罚参数u在第k+1次迭代中的值，μmax表示惩罚参数u的最大值，t表示步长，此处步长t取值为1.5。通过简单的推导，对于x^k+1与s^k+1的迭代可表示为式(16)～(17)：

式(16)～(17)中，||x||*表示分块干净图像x核范数，f(·)表示谱算子函数，||s||2,1表示矩阵化分块稀疏噪声s的l2,1范数，其余各符号与式(12)～(15)中相同。分别将子问题式(16)～(17)称为x子问题与s子问题。这两个子问题都有封闭解。关于x子问题和s子问题的解可分别由如下【引理1】和【引理2】得出。

【引理1】：给定一矩阵e^k计算如式(18)，它的奇异值分解(svd)可表示为式(19)，让α^k＝1/μ^k，其奇异值收缩算子服从式(20)，且有式(21)和式(22)；

式(18)中，e^k表示式(18)右边变量在第k次迭代中的值，y表示矩阵化分块退化高光谱图像，s^k表示矩阵化分块稀疏噪声s在第k次迭代中的值，f(·)表示谱算子函数，λ^k表示朗格朗日乘子λ在第k次迭代中的值，μ^k表示惩罚参数u在第k次迭代中的值；

e^k＝u^kσ^k(v^k)^t,σ^k＝diag(σ(e^k))(19)

式(19)中，矩阵e^k表示(18)式右边变量在第k次迭代中的值，矩阵σ^k对角线上的元素等于e^k的奇异值，u^k和v^k的列分别表示奇异值中的左、右奇异向量，σ(e^k)表示e^k的奇异值，diag(σ(e^k))用于构造一个对角矩阵，对角线上的元素等于e^k的奇异值。

式(20)中，x^k+1表示分块干净图像x在第k+1次迭代中的值，α^k＝1/μ^k，μ^k表示惩罚参数u在第k次迭代中的值，||x||*表示x的核范数，x表示分块干净图像，e^k表示式(18)右边变量在第k次迭代中的值，||x-e^k||f表示矩阵x-e^k的弗罗贝尼乌斯范数(f-范数)。

式(21)中，x^k+1表示分块干净图像x在第k+1次迭代中的值，矩阵σ^k为对角矩阵，对角线上的元素等于e^k的奇异值，u^k和v^k的列分别表示奇异值中的左、右奇异向量，参见式(22)。

式(22)中，e^k表示式(18)右边变量在第k次迭代中的值，α^k＝1/μ^k，μ^k表示惩罚参数u在第k次迭代中的值。

【引理2】：给定一矩阵w^k如式(23)，s子问题(17)有最优解，其第j列更新如式(24)。

式(23)中，y表示矩阵化分块退化高光谱图像，x^k+1表示分块干净图像x在第k+1次迭代中的值，λ^k表示拉格朗日乘子λ在第k+1次迭代中的值，μ^k表示惩罚参数μ在第k次迭代中的值。

式(24)中，s^k+1(:,j)表示矩阵化分块稀疏噪声s的第j列在第k+1次迭代中的值，||w^k(:,j)||2表示w^k的第j列的2-范数，λ表示规则化参数，μ^k表示惩罚参数μ在第k次迭代中的值。

步骤3)中对基于自适应秩校正的高光谱图像去噪模型采用交替方向乘子算法(admm)进行求解的详细步骤包括：

3.1)输入噪声污染的图像

3.2)使用采用非精确增广拉格朗日乘子算法(ialm)估计获得初始估计图像设置参数值s，λ，μ，δ＝10^-4，k＝0，其中s表示矩阵化分块稀疏噪声，λ表示拉格朗日乘子，μ表示惩罚参数，δ表示误差范围，k表示迭代次数；

3.3)判断是否满足收敛，当不收敛时，跳转执行下一步；否则，结束并退出；

3.4)根据式(18)计算矩阵e^k；

3.5)根据式(21)更新x^k+1；

3.6)根据式(23)计算矩阵w^k；

3.7)根据式(24)更新s^k+1；

3.8)根据式(14)更新λ^k+1；

3.9)根据式(15)更新μ^k+1；

3.10)跳转执行步骤3.3)。

步骤3.3)判断是否满足收敛的函数表达式为：

||y-x^k+1-s^k+1||f≤δ*||y||f(25)

式(25)中，y表示矩阵化分块退化高光谱图像，x^k+1表示分块干净图像x在第k+1次迭代中的值，s^k+1表示矩阵化分块稀疏噪声s在第k+1次迭代中的值，||y-x^k+1-s^k+1||f表示y-x^k+1-s^k+1的弗罗贝尼乌斯范数(f-范数)，δ表示误差范围，||y||f表示y的弗罗贝尼乌斯范数(f-范数)。

本实施例中，步骤4)中将所有的分块干净图像x组合获得最终估计的干净高光谱图像时，还包括针对分块干净图像x的重叠区域取平均值的方式处理重叠区域像素的步骤。

下文将通过实验对本实施例基于自适应秩校正的高光谱图像去噪方法进行验证。本实施例中利用washingtondcmall数据集中的一张波段数为191的256×256的图像进行模拟实验，添加6种不同的噪声情景，详情如下：1)独立同分布高斯噪声：测试数据集的每个频带添加零均值，标准差为0.05的高斯噪声。2)非独立同分布高斯噪声：测试数据集的每个频带添加零均值的高斯噪声。噪声的标准差是0到0.05之间的随机值。3)高斯+脉冲噪声：测试数据集的每个频带添加的高斯噪声如情形2)。另外，随机选取测试集中的30个波段添加脉冲噪声。添加的脉冲噪声比率在0.4到0.5之间随机取值。4)高斯+条带噪声：测试数据集的每个频带添加的高斯噪声如情形2)。随机选取测试集中的30个波段添加条带噪声。条带的数量在20到40之间随机取值。5)高斯+死线噪声：测试数据集的每个频带添加的高斯噪声如情形2)。随机选取测试集中的30个波段添加死线噪声。死线的数量在5到10之间随机取值。6)混合噪声：高斯噪声、脉冲噪声、条带、死线添加到测试数据集的每个频带。这四种噪声的添加方式如情形2)-5)所示。

作为本实施例基于自适应秩校正的高光谱图像去噪方法的对比，采用了5种流行的现有去噪算法,即四维块匹配滤波(block-matching4dfiltering，bm4d)，低秩矩阵恢复(low-rankmatrixrecovery，lrmr)，张量字典学习(tensordictionarylearning，tdl)，空间光谱全变差(spatio-spectraltotalvariation,sstv)，非独立同分布混合高斯(non-independentandidenticallydistributedmixtureofgaussians，nmog)，与本实施例基于自适应秩校正的高光谱图像去噪方法(arc)在以上的6种不同的噪声情形下进行比较，具体结果如下表1所示。评价指标有四种，分别是平均峰值信噪比(mpsnr)、平均结构相似度(mssim)、平均光谱角距离(msad)和计算时间(秒)。

表1：本实施例方法和5种现有去噪方法在四种去噪指标下的效果比较。

从表1可以看出，在独立同分布高斯噪声情景下，本实施例基于自适应秩校正的高光谱图像去噪方法(arc)略低于tdl；在其他复杂更加符合实际的噪声情形下，本发明所提出的arc在所有评估指标方面都优于其他方法。

图4为流行的5种现有去噪方法与本实施例方法(arc)针对washingtondcmall数据集中在混合噪声情景下的第48个波段去噪结果对比；其中(a)混合噪声情景下的噪声图像，(b)bm4d去噪后的图像，(c)lrma去噪后的hsi，(d)tdl去噪后的图像，(e)sstv去噪后的图像，(f)nmog去噪后的结果，(g)arc去噪后的图像，及(h)原清洁图像。图5为流行的5种现有去噪方法与本实施例方法(arc)针对washingtondcmall数据集中在混合噪声情景下的第57个波段去噪结果对比；其中(a)混合噪声情景下的噪声图像，(b)bm4d去噪后的图像，(c)lrma去噪后的hsi，(d)tdl去噪后的图像，(e)sstv去噪后的图像，(f)nmog去噪后的结果，(g)arc去噪后的图像，及(h)原清洁图像。图6为5种现有去噪方法与本实施例方法(arc)在真实的遥感图像数据集avirisindianpines中频带数为108的去噪结果对比图；其中(a)被污染的噪声图像，(b)bm4d去噪后的图像，(c)lrma去噪后的hsi，(d)tdl去噪后的图像，(e)sstv去噪后的图像，(f)nmog去噪后的结果，及(g)arc去噪后的图像。图7为5种现有去噪方法与本实施例方法(arc)在真实的遥感图像数据集avirisindianpines中频带数为150的去噪结果对比图；其中(a)被污染的噪声图像，(b)bm4d去噪后的图像，(c)lrma去噪后的hsi，(d)tdl去噪后的图像，(e)sstv去噪后的图像，(f)nmog去噪后的结果，及(g)arc去噪后的图像。图8为5种现有去噪方法与本实施例方法(arc)在真实的遥感图像数据集hydiceurban中频带数为103的去噪结果对比图；其中(a)被污染的噪声图像，(b)bm4d去噪后的图像，(c)lrma去噪后的hsi，(d)tdl去噪后的图像，(e)sstv去噪后的图像，(f)nmog去噪后的结果，及(g)arc去噪后的图像。图9为5种现有去噪方法与本实施例方法(arc)在真实的遥感图像数据集hydiceurban中频带数为138的去噪结果对比图；其中(a)被污染的噪声图像，(b)bm4d去噪后的图像，(c)lrma去噪后的hsi，(d)tdl去噪后的图像，(e)sstv去噪后的图像，(f)nmog去噪后的结果，及(g)arc去噪后的图像。结合图4～图9也可以看出，本实施例基于自适应秩校正的高光谱图像去噪方法(arc)可以有效地去除混合噪声并且很好地保留真实的图像结构，具有重要的应用价值。

此外，本实施例还提供一种基于自适应秩校正的高光谱图像去噪系统，包括计算机设备，该计算机设备被编程或配置以执行本实施例前述基于自适应秩校正的高光谱图像去噪方法的步骤。此外，本实施例还提供一种基于自适应秩校正的高光谱图像去噪系统，包括计算机设备，该计算机设备的存储介质上存储有被编程或配置以执行本实施例前述基于自适应秩校正的高光谱图像去噪方法的计算机程序。此外，本实施例还提供一种计算机可读存储介质，该计算机可读存储介质上存储有被编程或配置以执行本实施例前述基于自适应秩校正的高光谱图像去噪方法的计算机程序。此外，本实施例还提供一种基于自适应秩校正的高光谱图像去噪系统，包括：第一程序单元，用于对退化高光谱图像矩阵化后采用非精确增广拉格朗日乘子算法进行初始估计，得到初始估计的干净高光谱图像第二程序单元，用于对退化高光谱图像初始估计的干净高光谱图像分别进行分块和矩阵化处理，分别得到矩阵化分块退化高光谱图像y、矩阵化分块初始估计的干净高光谱图像x；第三程序单元，用于针对各个矩阵化分块退化高光谱图像y、矩阵化分块初始估计的干净高光谱图像x建立基于自适应秩校正的高光谱图像去噪模型，基于自适应秩校正的高光谱图像去噪模型通过自适应秩校正惩罚函数以实现自适应地抵消核范数对大的奇异值的惩罚，且基于l2,1范数消除稀疏噪声；对基于自适应秩校正的高光谱图像去噪模型采用交替方向乘子算法进行求解，得到分块干净图像x；第四程序单元，用于将所有的分块干净图像x组合获得最终估计的干净高光谱图像。

以上所述仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。