Elsevier

Displays

Volume 73, July 2022, 102200
Displays

Robust Low-Rank Analysis with Adaptive Weighted Tensor for Image Denoising

https://doi.org/10.1016/j.displa.2022.102200Get rights and content

Highlights

  • The adaptive weight tensor can effectively retain useful singular values and better preserve the low-rank properties of the unfolding matrix.

  • The proposed algorithm considers the spatial information and spectral information at the same time.

Abstract

In order to obtain better denoising results, this paper proposes the Robust Low-Rank Analysis with Adaptive Weighted Tensor (AWTD) method for image denoising tasks. On one hand, it uses the latest adaptive weight tensor, which obtains the low-rank approximation of the tensor by adding adaptive weights to the unfolding matrix of the tensor. The adaptive weight tensor can effectively retain useful singular values and better preserve the low-rank properties of the unfolding matrix. On the other hand, the proposed algorithm considers the spatial information and spectral information at the same time: for the RGB images, it retains the structural information inside the image patch and the connection between different channels (the spatial information of the image); for the hyperspectral images, it also retains the spectral information of the hyperspectral images. The experimental results show that the proposed method is superior to other test methods.

Introduction

The purpose of image denoising is to improve the quality of a given image due to noise. Various denoising technologies can effectively improve the image quality and increase the signal-to-noise ratio, to better reflect the original information of the image. With the development of computer technology and digital image processing technology, image denoising has been widely used in many fields such as medicine, public security, and military. Existing image denoising algorithms can be divided into traditional two-dimensional denoising methods and tensor-based multi-linear models.

There are many existing two-dimensional denoising algorithms, and they include local and non-local methods. The local methods [1], [2] use a designed filter to perform a convolution operation on the entire image, and they will lose the global structural information of the image. Correspondingly, there are many methods based on non-local algorithms [3], [4], [5], [6]. These algorithms use non-local means for image denoising and use the self-similarity of the image to remove noise. In addition, the method of rank minimization is also a popular research direction [7], [8], [9], [10], [11], [12]. The rank minimization problem is NP-hard, but in some cases, nuclear norm minimization can be used to approximate the original problem [13], [14]. Although the above methods have achieved good results in image denoising, they are only suitable for gray-scale image denoising tasks. Because they only consider the structural information inside the image patch and ignore the relationship between the different color channels in the RGB image, which limits their usefulness for extracting information from a multi-dimensional perspective.

The tensor-based multi-linear model can better deal with high-dimensional denoising tasks and can use the multi-linear structure to provide better understanding and higher accuracy. There are currently two popular directions of the tensor model: tensor decomposition and tensor low-rank. The tensor decomposition is usually divided into CANDECOMP/PARAFAC (CP) decomposition [15], [16] and Tucker decomposition (High-order SVD) [17], [18], as well as the newly proposed Ring decomposition model [19] and the special form of Tucker decomposition [20], [21] in recent years. The tensor decomposition is very similar to the Principal Component Decomposition (PCA) [22], [23] of the matrix, and PCA is sensitive to outliers and serious pollution (non-Gaussian noise), so the CP and Tucker decomposition also cannot successfully deal with these problems. Someone has proposed some algorithms to strengthen tensor decomposition to solve the problem of outliers and missing data [24], [25], [26]. Unfortunately, they lack global optimality guarantees.

In practice, whether it is caused by outliers or missing data, tensor data is often low-rank. This means that the main part of the data in practice is usually affected by a small number of potential factors. Therefore, the nuclear norm corresponding to the matrix is minimized, and the image denoising task can be accomplished by minimizing the tensor nuclear norm [27], [28], [29]. Although these methods have been successfully applied to color image denoising, they have not considered the adaptive weight of the tensor nuclear norm and therefore lack flexibility.

In summary, it can be concluded that on the one hand, the traditional two-dimensional methods are matrix-based and only consider the local information of the image patches or the structure information of the entire image patch. This ignores the relationship within the image block or between the entire image block and different channels (spatial information), which limits the scope of application of the algorithm in practical applications. On the other hand, in tensor-based multi-linear models, tensor decomposition cannot well overcome the problems of outliers and missing data. In the tensor low-rank models, the tensor nuclear norm minimization strategy of low-rank approximation does not take into account the adaptive weight of the nuclear norm, so it lacks flexibility.

In recent years, Kwok-Po Ng [30] proposed an adaptive weight tensor complement method for hyperspectral remote sensing image denoising (AWTC). By considering the contribution of spatial, spectral, and temporal information in each dimension to adaptively determine the weights, a new weighted tensor low-rank regularization model is constructed. In this paper, we apply the adaptive weight tensor proposed in the AWTC model to the High-order RPCA (HORPCA) [29] model in the tensor nuclear norm minimization strategy, and propose Robust Low-Rank Analysis with Adaptive Weighted Tensor for Image Denoising (AWTD). AWTD not only overcomes the problem of losing spatial information in the traditional two-dimensional method but also introduces an adaptive weight tensor, which is more flexible than the High-order RPCA (HORPCA) model.

In this paper, we make three following main contributions:

  • By introducing the adaptive weight tensor, a more flexible Robust Low-Rank Analysis with Adaptive Weighted Tensor for Image Denoising (AWTD) model is proposed to better complete the low-rank regularization of the tensor.

  • The results of different types of denoising experiments are shown, including RGB images with special and complex structures and hyperspectral remote sensing images.

  • The distribution of the weights of different unfolding matrices on different noise levels and different types of images is discussed. The sensitivity analysis of the parameters and the convergence analysis of the algorithm are given.

The rest of this paper is organized as follows. In Section 2, we will give an introduction to mathematical symbols and review some definitions of matrices and tensors. We will first introduce the adaptive weight in Section 3, then we will propose our model: Robust Low-Rank Analysis with Adaptive Weighted Tensor for Image Denoising (AWTD). The experimental results and the sensitivity analysis of the parameters will be shown in Section 4. Finally, this paper will end with a conclusion in Section 5.

Section snippets

Related work

Referring to the practice in [31], we denote tensors by boldface Euler script letters, e.g., X, matrices by boldface capital letters, e.g., X, vectors by boldface lowercase letters, e.g., x. The mode-i unfolding(matricization) of the tensor X is the matrix denoted by X(i). The vectorization of X is denoted by vec(X).

Adaptive weighted tensor

The tensor trace norm was first proposed in [32] to solve the multiband structure of tensor data. A multiband image can be regarded as a third-order tensor, then its unfolding matrices in three dimensions are X(1), X(2), and X(3) respectively. Then the tensor trace norm can be described as: X=wii=13X(i).where wi is the weight corresponding to the ith dimension unfolding, and the weight wi has used the equal value in [32], i.e. 1/3. Although it has been successfully applied to tensor

Experimental results and analysis

In this section, we will give the comparison results of the visual and evaluation indicators, and give the sensitivity analysis of the parameters. The details are as follows.

Conclusion

In this article, we propose an AWTD strategy for the image denoising problem, which applies the adaptive weight tensor to the low-rank tensor model. It considers the internal structural information of the image block, the relationship between the different channels of the image block (spatial information) and spectral information. According to the low-rank strength of the unfolding matrix in different dimensions, an adaptive method is adopted to determine the weights. Finally, the ADAL method

CRediT authorship contribution statement

Lei Zhang: Conceptualization, Methodology, Software, Data curation, Writing – original draft, Data curation, Writing – original draft, Software, Validation, Writing – review & editing. Cong Liu: Supervision.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgment

The authors would like to appreciate all anonymous reviewers for their insightful comments and constructive suggestions to polish this paper in high quality. This research was supported by the National Natural Science Foundation of China (No. 61703278, 61673180).

References (50)

  • NejatiM. et al.

    Denoising by low-rank and sparse representations

    J. Vis. Commun. Image Represent.

    (2016)
  • GabayD. et al.

    A dual algorithm for the solution of nonlinear variational problems via finite element approximation

    Comput. Math. Appl.

    (1976)
  • PortillaJ. et al.

    Image denoising using scale mixtures of Gaussians in the wavelet domain

    IEEE Trans. Image Process.

    (2003)
  • Da CunhaA.L. et al.

    The nonsubsampled contourlet transform: theory, design, and applications

    IEEE transactions on image processing

    (2006)
  • DongW. et al.

    Nonlocal image restoration with bilateral variance estimation: a low-rank approach

    IEEE Trans. Image Process.

    (2013)
  • JinW. et al.

    Fast non-local algorithm for image denoising

  • PapyanV. et al.

    Multi-scale patch-based image restoration

    IEEE Trans. Image Process.

    (2016)
  • DabovK. et al.

    Image denoising by sparse 3-D transform-domain collaborative filtering

    IEEE Trans. Image Process.

    (2007)
  • GuS. et al.

    Weighted nuclear norm minimization with application to image denoising

  • UniversityS. et al.

    Robust principal component analysis?

    (2011)
  • PartridgeM. et al.

    Robust principal component analysis

  • NieF. et al.

    Low-Rank Matrix Recovery Via Efficient Schatten P-Norm Minimization

    (2012)
  • MohanK. et al.

    Iterative reweighted algorithms for matrix rank minimization

    J. Mach. Learn. Res.

    (2012)
  • CandesE. et al.

    Exact low-rank matrix completion via convex optimization

  • CaiJ. et al.

    A singular value thresholding algorithm for matrix completion

    SIAM J. Optim.

    (2018)
  • CarrollJ. et al.

    Analysis of individual differences in multidimensional scaling via an n-way generalization of “Eckart-Young” decomposition

    Psychometrika

    (1970)
  • HarshmanR.

    Foundations of the parafac procedure : Models and conditions for an “explanatory” multimodal factor analysis

    Ucla Work. Pap. Phonetics

    (1970)
  • TuckerL.

    Some mathematical notes on three-mode factor analysis

    Psychometrika

    (1966)
  • LathauwerL. et al.

    Multilinear singular value tensor decompositions

    SIAM J. Matrix Anal. Appl.

    (2000)
  • HeW. et al.

    Remote sensing image reconstruction using tensor ring completion and total variation

    IEEE Trans. Geoence Remote Sens.

    (2019)
  • BahriM. et al.

    Robust kronecker component analysis

    IEEE Trans. Pattern Anal. Mach. Intell.

    (2018)
  • BahriM. et al.

    Robust Kronecker-decomposable component analysis for low-rank modeling

    (2017)
  • JolliffeI.

    Principal component analysis

    J. Mar. Res.

    (2002)
  • AbdiH. et al.

    Principal component analysis

    Wiley Interdiscip. Rev. Comput. Stat.

    (2010)
  • EngelenS. et al.

    A fully robust PARAFAC method for analyzing fluorescence data

    Journal of Chemometrics: A Journal of the Chemometrics Society

    (2009)
  • This paper was recommended for publication by Jing Liu.

    View full text