高光谱成像技术不仅能够精细描述真实场景,提供丰富的光谱空间特征,而且具有广
泛的应用,它可提高物体检测和分类的准确性。但在采集和传输高光谱图像过程中,不可
避免地受到各种噪声的破坏
[1]
,这不仅使得高光谱图像的质量大打折扣,而且对后续应用
(如图像分解、图像检测)产生不利影响。近年来,关于高光谱图像去噪的方法有很多。假
设观测的高维图像看作一个张量,记为 y,用数学表达式可以描述为
其中: x∈RI1×I2×⋯×IN 为干净图像; N∈RI1×I2×⋯×IN 是噪声项。高光谱图像(HSI)去噪
目的是从观察到的 y 中去估计 x。为了有效利用低秩模型的局部及非局部自相似性,Wu 等
[2]
通过有效地控制 CP 秩提出了一种加权秩-1 张量分解方法(称为 WTR1),用于高光谱图像
去噪。Xu 等
[3]
提出了一种基于塔克分解的高光谱图像恢复方法,通过在合适的窗口尺寸中
利用非局部相似性,从高维图像中提取的 3-D 全波段块分组来形成三阶张量。Yuan 等
[4]
提
出了一种基于空间光谱总变差的高光谱图像去除噪声算法,高光谱图像
[19-21]
的固有结构是
通过利用沿空间维的二维总变差和沿光谱维的一维总变差来进行优化的。上述方法都是将
高维图像降维,这样会破坏图像块的二维结构。基于此,Gao 等
[5]
提出一种高阶奇异值分
解的高光谱图像恢复方法,将几个相似的图像块堆叠成三维数组,并视为三阶张量,然后
进行高阶奇异值分解,针对分解得到的核心张量,设计了一种自适应硬阈值系数收缩的迭
代算法来求解模型。通过利用张量环(TR)潜在空间的低秩结构,Xiong 等
[6]
提出了一种新的
高光谱图像恢复方法,该模型在潜在 TR 因子上引入核范数正则化,利用奇异值分解(SVD)
的优化步骤执行的数据规模要小得多。为了使秩估计更加准确,Xue 等
[7]
提出了一种用于
非局部低秩正则化 CP 分解模型。沿波段全局相关性(GCS)和空间非局部自相似性(NSS)先
验信息在低秩正则化 CP 分解下有效地进行探索,该方法可减少张量秩估计偏差对降噪性
能的影响。
上述方法均是将高维张量转化为向量或矩阵,虽然易于处理,但是会导致信息丢失和
性能下降。在经典分解中张量秩对参数的反应很敏感,估计最优秩成为重要的一步,如果
把握不好度,就会出现过度去噪和无法将噪声去除的现象,并且其中涉及的参数比较多,
导致计算量大。针对上述缺陷,提出了一种基于张量环分解的非局部正则化高光谱图像去
除噪声算法(TRTD-NRM)。将张量环用于高光谱图像的混合噪声去除,且在非局部的基础
上进行,与以往张量分解相比,张量环分解无需估计最佳秩,且数据量较小,有利于后续
计算。数值实验结果表明,在主观和客观方面,本算法优于现有算法。
1. 张量环分解
张量环(TR)分解
[8]
的目的是通过一系列三阶张量的循环相乘来表示高阶(或高维)张
量。具体来说,令 T 为 n
1
×n
2
×…×n
d
的 d 阶张量,表示为 T∈ R
n
1
×n
2
×…×n
d
,TR 表示是将其
分解为一系列潜在张量 T∈ R
n
1
×n
2
×…×n
d
,k=1, 2, …, d,可用以下元素形式表示:
T(i1,i2,⋯,id)=Tr{F1(i1)F2(i2)⋯Fd(id)}=Tr{∏k=1dFk(ik)},