基于三重低秩模型的高光谱图像去噪方法与流程

本发明涉及一种解决高光谱图像中条纹、死线、椒盐噪声等类型噪声的去噪算法，属于遥感图像处理领域。

背景技术：

高光谱图像是一种具有数十、甚至上百个连续光谱波段的遥感图像，能提供目标区域的丰富的光谱信息，在地物分类和环境监测等领域具有广泛应用。然而在成像过程中，由于受到大气干扰、物理设备自身技术上的局限性以及传输带宽等因素影响，获取的高光谱图像往往受到各种噪声的影响，降低了高光谱图像的质量，对高光谱图像后续应用产生不利的影响。高光谱图像去噪是一种能有效降低噪声的图像处理技术，受到了人们的广泛关注。

常见高光谱图像去噪算法包括小波变换算法、优化重构算法和稀疏低秩算法等。小波变换算法主要是通过对高光谱图像进行小波变换，将其分解成不同尺度上的高频和低频，噪声主要集中在高频子带，从而实现对图像信息与噪声之间的分离。全变分约束是优化重构类算法常用的先验约束，其主要约束图像中相邻像素之间的差，但全变分约束产生阶跃效应。稀疏低秩算法是目前较流行的一类高光谱图像去噪技术，比如鲁棒主成分分析。通过低秩约束刻画无噪声高光谱图像中光谱间的关联性，并通过稀疏特性约束噪声，从而实现从含噪声的高光谱图像中重构出低秩无噪声图像。后续很多基于稀疏低秩模型的高光谱去噪算法都由鲁棒主成分分析算法衍生而来，比如截断核范数、图约束低秩模型以及双重鲁棒主成分分析等。现有的低秩模型主要从核范数的定义、低秩重构误差以及图像分布等方面提升低秩去噪模型的性能。然而，现有的鲁棒主成分分析模型及其改进形式对低秩部分的自相关性考虑还不够充分。

技术实现要素：

技术问题：本发明目的是提供一种新的高光谱图像低秩去噪算法，分别通过对低秩部分、低秩部分与线性变换域噪声图像之间的约束以及低秩部分的自相关性来进行三重低秩约束，从而去除高光谱图像中例如椒盐噪声、死线噪声等稀疏噪声。

技术方案：本发明的方法首先从含噪声的高光谱图像中依次提取图像块，再对提取的所有图像块分别进行三重低秩优化重构模型(即包括低秩部分约束、低秩部分与线性变换域噪声图像之间的约束以及低秩部分自相关性约束)，去除稀疏噪声，最后通过对重构图像块的重叠像素进行加权平均得到去噪的高光谱图像。

由于高光谱传感器受传输设备技术上的局限性以及大气折射等因素的影响高光谱图像往往含有很多噪声，比如死线噪声、条纹噪声等为此，本发明考虑一种下面的噪声模型

x＝l+e(1)

其中表示含噪声的高光谱图像，为干净图像，为死线、条纹等稀疏噪声，m×n是高光谱图像的空间尺寸，b是高光谱图像的波段数。

对于上述含噪高光谱图像的处理本发明提供了一种三重低秩约束模型，包括低秩部分约束、低秩部分与线性变换域噪声图像之间的约束以及低秩部分自相关性约束。

本发明具体步骤包括：

1.从噪声高光谱图像x中依次提取m×n×b的3-d图像块并变换为mn×b的矩阵xi，j作为该算法模型的输入。

2.通过以下优化算法针对每一个xi，j进行优化重构：

其中ai，j是原始数据xi，j的线性变换矩阵，bi，j是低秩部分li，j的自相关变换矩阵。||li，j||*

是表示关于li，j的核范数，||ei，j||1则是关于ei，j的l1范数，||·||f表示矩阵的f-范数，λi(i＝1，2，3，4，5)是平衡系数。

3.利用增广拉格朗日函数法对模型(2)进行求解。

首先，引入辅助变量将问题(2)等价变换为：

s.t.xi，j＝li，j+ei，j，ai，j＝ni，j，mi，j＝li，j-axi，j，pi，j＝bi，j，ki，j＝li，j(3)

问题(3)的增广拉格朗日函数为：

其中μ为惩罚参数，y1，y2，y3，y4，y5为拉格朗日乘子。

然后，关于增广拉格朗日函数(4)的最小优化问题通过下面的子问题进行交替迭代优化。

1)li，j子优化问题

其中，

ω(γ)＝us(∑γ)v^*是关于矩阵γ的奇异值阈值算子，∑γ＝diag({ηi}1≤i≤γ)是奇异值矩阵，u和v分别为左右正交矩阵，sε(x)＝sgn(x)*max(|x|-ε，0)是软阈值算子。

2)ei，j子优化问题

其中

3)ni，j子优化问题

4)mi，j子优化问题

5)bi，j子优化问题

6)pi，j子优化问题

7)ki，j子优化问题

8)ai，j子优化问题

其中

9)乘子更新

μ＝min(ρμ，μmax)(14)

其中ρ＞1μmax是常数。

10)交替迭代循环1)-9)，判断终止条件{||xi，j-li，j-ei，j||∞＜ε，||mi，j+ai，jxi，j-li，j||∞＜ε，||ai，j-ni，j||∞＜ε，||bi，j-pi，j||∞＜ε，||li，j-ki，j||∞＜ε}。如果满足终止条件则输出低秩部分li，j。

4.最后通过对重构图像块li，j的重叠像素进行加权平均得到去噪的高光谱图像l。

有益效果：本发明设计了一种三重低秩约束模型，包括低秩部分约束、低秩部分与线性变换域噪声图像之间的约束以及低秩部分自相关性约束。其中，低秩部分与线性变换域噪声图像之间的约束可以提高模型的鲁棒性以及图像去噪的精度，低秩部分自相关性约束可以使得该去噪模型相比较于其他去噪模型对死线噪声的去除效果更佳优越。

附图说明

图1为本发明一种基于三重低秩模型的高光谱图像去噪算法流程图。

具体实施方式

下面结合附图对本发明作进一步的说明。

本发明的优化算法包含的步骤为：

步骤1：输入含有稀疏噪声的高光谱图像终止误差ε，正则化系数λ1、λ2、λ3、λ4、λ5。进入步骤2；

步骤2：利从高光谱图像x中依次提取尺寸为mn×b的图像块

xi，j(i∈[1，m-m+1]，，j∈[1，n-n+1])。进入步骤3；

步骤3：对参数进行初始化li，j＝0，ei，j＝0，mi，j＝0，ki，j＝0，

y1＝0，y2＝0，y3＝0，y4＝0，y5＝0，μ＝0.1，

μmax＝10⁶，ai，j＝0，ni，j＝0，bi，j＝0，pi，j＝0，ρ＞1，

进入步骤4；

步骤4：更新li，j＝ω(γ)，

进入步骤5；

步骤5：更新进入步骤6；

步骤6：更新进入步骤7；

步骤7：更新进入步骤8；

步骤8：更新bi，j＝(μpi，j-y4+2λ3λ5k^tl)(2λ3λ5ki，j^tki，j+μi)^-1。进入步骤9；

步骤9：更新进入步骤10；

步骤10：更新ki，j＝(μl+y5+2λ3λ4li，jb^t)(2λ3λ4bi，jbi，j^t+μi)^-1。进入步骤11；

步骤11：更新进入步骤12；

步骤12：更新拉格朗日乘子

进入步骤13；

步骤13：更新惩罚系数μ＝min(ρμ，μmax)，进入步骤14；

步骤14：判断终止条件||xi，j-li，j-ei，j||∞＜ε，||ai，j-ni，j||∞＜ε，

||mi，j+ai，jxi，j-li，j||∞＜ε，||bi，j-pi，j||∞＜ε，||li，j-ki，j||∞＜ε，

则停止迭代，否则转向步骤4。进入步骤15；

步骤15：通过对重构图像块li，j的重叠像素进行加权平均得到去噪的高光谱图像l。

进入步骤16；

步骤16：输出去噪的高光谱图像l。

本发明测试了indianpine(145×145×224)和washingtondc(256×256×191)高光谱图像。设置图像块大小m×n为20×20，提取图像块步长为8×8(即：从横向和纵向移动8个像素)，参数λi(i＝1，2，3，4，5)分别设置为λ1＝0.125、λ2＝0.1、λ3＝0.1、λ4＝0.0003、λ5＝0.0003，终止误差为ε＝10^-6。实验测试了三种不同类型的模拟噪声，其分别定义为：

噪声1：每个波段添加强度为30％的椒盐噪声；

噪声2：每个波段添加强度为20％的椒盐噪声，并在波段126-145中添加死线噪声(宽度为1至3个像素)；

噪声3：每个波段添加强度为30％的椒盐噪声，并在波段126-145中添加条纹噪声(宽度为1至3个像素)；

本专利采用两种图像质量指标：平均峰值信噪比(mpsnr)和平均结构相似度(mssim)定义分别如下所示：

其中，psnri和ssimi表示ith波段的psnr和ssim值。

其中，表示x最大值的平方，x表示原始图像，l表示去噪图像，μx是x的平均值，μl是l的平均值，是x的方差，是l的方差，σxl是xl的协方差，c1＝(k1γ)²，c2＝(k2γ)²，γ是像素值的动态范围，k1＝0.01，k2＝0.03。

表1给出了indianpine图像和washingtondc图像在各种类型噪声下的去噪结果。结果表明本发明提出三重低秩去噪算法能有效地去除高光谱图像中的椒盐、条纹以及死线噪声。

表1：indianpine图像和washingtondc图像在各种噪声类型下的去噪结果

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