Next Article in Journal
SaltISNet3D: Interactive Salt Segmentation from 3D Seismic Images Using Deep Learning
Previous Article in Journal
D-Net: A Density-Based Convolutional Neural Network for Mobile LiDAR Point Clouds Classification in Urban Areas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Edge-Preserved Low-Rank Representation via Multi-Level Knowledge Incorporation for Remote Sensing Image Denoising

1
School of Electronic and Optical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China
2
School of Computer Science and Communication Engineering, Jiangsu University, Zhenjiang 212013, China
3
Department of Computer Science, University of Nigeria, Nsukka 410001, Nigeria
4
Jiangsu Yangjing Petrochemical Group Co., Ltd., Lianyungang 222000, China
5
Nanjing Les Information Technology Co., Ltd., Nanjing 210000, China
6
Key Laboratory of Digital Earth Science, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
*
Author to whom correspondence should be addressed.
Current address: Nanjing University of Science and Technology, Xiaolingwei 200, Nanjing 210094, China.
Remote Sens. 2023, 15(9), 2318; https://doi.org/10.3390/rs15092318
Submission received: 7 April 2023 / Revised: 25 April 2023 / Accepted: 25 April 2023 / Published: 27 April 2023
(This article belongs to the Section Remote Sensing Image Processing)

Abstract

:
The low-rank models have gained remarkable performance in the field of remote sensing image denoising. Nonetheless, the existing low-rank-based methods view residues as noise and simply discard them. This causes denoised results to lose many important details, especially the edges. In this paper, we propose a new denoising method named EPLRR-RSID, which focuses on edge preservation to improve the image quality of the details. Specifically, we considered the low-rank residues as a combination of useful edges and noisy components. In order to better learn the edge information from the low-rank representation (LRR), we designed multi-level knowledge to further distinguish the edge part and the noise part from the residues. Furthermore, a manifold learning framework was introduced in our proposed model to better obtain the edge information, as it can find the structural similarity of the edge part while suppressing the influence of the non-structural noise part. In this way, not only the low-rank part is better learned, but also the edge part is precisely preserved. Extensive experiments on synthetic and several real remote sensing datasets showed that EPLRR-RSID has superior advantages over the compared state-of-the-art (SOTA) approaches, with the mean edge protect index (MEPI) values reaching at least 0.9 and the best values in the no-reference index BRISQUE, which represents that our method improved the image quality by edge preserving.

1. Introduction

Remote sensing images have been widely used in a variety of areas, such as agricultural resource research, city planning, and others. However, additive noise is a critical problem for remote sensing images [1], which seriously affects their practical use and subsequent processing. Thus, remote sensing image denoising has been a hot issue. As a result, various remote sensing image denoising research works have been performed in the last few decades, especially for multispectral, hyperspectral, and SAR data [2].
Traditional approaches are mainly filter-based denoising methods with strong theoretical support. These filters can operate in either the spatial or frequency domains directly [3]. The Gaussian filter [4], bilateral filter [5], and guided image filter [6] are all classical spatial domain filters, which use image local information to help remove noise. Unlike these previous image filters, a new method called NL-means [7] has been proposed. It utilizes the redundant information of the original image to obtain a better result. Inspired by this, block-matching and 3D filtering (BM3D) [8] makes improvements on the basis of NL-means, which not only uses intra-fragment correlation, but also takes inter-correlation into account. The total variation (TV) algorithm [9] is an anisotropic model whose core is to use the gradient descent method to remove noise. The superiority of this method is that it smooths images while keeping the edges as much as possible. Due to the advantage of preserving the edge information, TV regulation has been added to many state-of-the-art (SOTA) denoising models.
Recently, low-rank representation (LRR) has gained much attention and has been successfully extended to the denoising field [10,11,12]. The low rank is an important feature, which indicates that high-dimensional data can be projected into a low-dimensional subspace. Based on this, several methods have been proposed for denoising. A method called NAILRMA [13] proposes an iterative regularization framework and tries to realize the purpose of denoising from the signal subspace. In NAILRMA, a noise-variance-based adaptive iteration factor selection method is used for each hyperspectral image (HSI) band. Meanwhile, some tensor-based low-rank denoising methods [14,15,16] have achieved satisfactory results.
However, such low-rank-based methods have a common problem, that is the over-smoothing caused by the loss of image edges. This drawback is fairly common due to the characteristic of LRR. The traditional LRR model filters out high-frequency edges as residues, which will be considered as only noise and be disregarded. In order to compensate for this defect, many researchers began to explore methods to enhance image edges. Some research has introduced TV regulation into low-rank models [17,18,19,20], trying to take advantage of the edge-preserving ability of the TV algorithm. A new method that makes use of low-rank regularization and a texture preserving prior [21] has then been proposed. Meanwhile, LRTL0 [22] uses the l 0 gradient regularization to realize edge preservation. In addition, more cutting-edge methods combining the latest developments and denoising models have been proposed. Liu et al. proposed a multi-source-information-based method [23], which employs priors to dictionary learning and edge feature predicting. It uses uncontaminated images as references to construct its own denoising object function, so that more details can be preserved through the similarities between the noisy image and clean reference image. A deep learning denoising network named RSIDNet [24] improves the usage of feature extraction through the combination of multiple modules, trying to reduce the loss of important details in the process of image denoising. The nonlocal 3D convolutional neural network (NL-3DCNN) [25] makes the best of the images’ nonlocal similarity and simultaneously utilizes the advantages of machine learning and deep learning. A large amount of novel approaches focus on tensor decomposition [26,27,28], which can avoid converting the original dimensions to keep as much information as possible. Furthermore, some filter-based methods [29,30,31] have been developed in recent years. They improve existing designs from a number of perspectives, including thresholding, image pre-processing, and building 3D dynamic kernels to address the shortcomings of traditional filters. Even though great efforts have been made in the edge-preserving field, the performances of existing methods still need to be improved to extract clearer edges by taking advantage of the edge information in residues.
In this paper, a novel edge-preserved denoising method for remote sensing images named EPLRR-RSID is proposed by concurrently describing the edge characteristics and LRR model. Most existing methods only care about whether the noise in remote sensing images has been effectively removed, which leads to the loss of edge features. Although a few advanced methods have considered the importance of edge preservation, their performances are not good enough due to the effect of inevitable noise. To overcome the aforementioned problems, we explored a more robust method to find edge features. We creatively considered the residues of LRR models as a mixture of high-frequency edge components and noisy components. Thus, we tried to extract clearer edges from the residues to remedy the shortcoming of traditional LRR methods. Unlike conventional methods, we considered the residues of LRR models as a mixture of high-frequency edge components and noisy components. Thus, we tried to extract clearer edges from the residues to remedy the shortcoming of traditional LRR methods. As shown in Figure 1, we first propose a new prior knowledge to better find structural edge information in the enhanced low-rank part of the residues. Then, manifold learning was integrated with our model to gain structural edge information discovered by the prior knowledge much better with affinity graph regularization. Through this framework, valuable edge information was extracted from those residues with smaller distances between them while the noisy part was suppressed due to its disorder. After the multi-level modeling, we can obtain the low-rank component, which contains blurred boundaries and edges, together with the distinct edge component, respectively. Thus, our proposed EPLRR-RSID can better preserve edge information while reducing noise, which provides high-quality denoised images with more edge details.
We list our main contributions as follows:
  • Different from conventional LRR-based denoising methods, we took a very different way to deal with low-rank residues, i.e., we did not simply discard them, but viewed those residues as a combination, which includes both edge information and the noisy part. By taking this different view, we tried to extract structural edge information in the residues to enrich the details of remote sensing images. To our best knowledge, this is the first time that the edge part in the residues has been valued and fully utilized.
  • In order to better learn the edge information from the residues, we designed a new multi-level prior knowledge regulation. By incorporating it, our proposed method can further separate the edge part from the noise part. Thus, the structural edge information can be further discovered by the manifold learning framework, which can find the structural similarity of the edge part in the enhanced low-rank part of the residues. This enabled our algorithm to preserve more high-frequency edge components.
  • Several experiments were conducted to prove the effectiveness of EPLRR-RSID. According to the experimental results, EPLRR-RSID had an outstanding advantage in preserving image edges, which made our method more competitive.

2. Related Works

As we previously mentioned, LRR has laid the foundation for many low-rank-based denoising methods [11,32,33]. In the previous works, we noticed that many researchers tried to design rational image priors that could inspire some significant properties [14]. A low-rank tensor recovery method named LLRT [34] introduces an analysis-based hyper-Laplacian prior to address the problem of ringing artifacts. This prior can help the model better use the structure information of the images. A new kind of mixed prior was proven effective, which combines the external patch prior and the internal self-similarity prior [35]. A three-stage mixed denoising model [36] was creatively proposed, which employs a new global-guided nonlocal prior mechanism. These methods try to find the most-appropriate prior to reconstruct the noise-free image from the low-rank model.
Among these methods, a weighted group LRR model called WGLRR [37] has attracted our attention. According to WGLRR, it is quite hard to recover a noise-free image from a noisy image if there is no side information for a heavily corrupted image. Different from other methods that use prior knowledge, it introduces a local image rank-ordered absolute differences statistic called ROAD [38] to construct its prior knowledge P creatively. After obtaining the P matrix, the WGLRR model is written as
min Z , E Z * + λ E F 2 s . t . P X = P Z + E
where X is the image with the original noise, Z is the low-rank component, and E is the noisy component. P represents a weight matrix to constrain the fidelity of recovered clean patches. “∘” denotes the Hadamard product.
However, such LRR methods have an unavoidable problem: they will filter out high-frequency signals together with corrupt elements during the denoising procedure, which means useful edge information will be abandoned. This drawback is inherited from the characteristic of the traditional LRR model, which only focuses on how to preserve a low-rank structure while ignoring the importance of edge information.

3. Methodology

We define the objective function and provide an optimization strategy that can successfully tackle the problem. What needs to be declared in advance is that the whole model was completed under the NLM framework [39]. We divided the noisy image into several patches to create group matrices, and then, our method was applied to each group matrix.

3.1. Model Formulation

As mentioned earlier, most low-rank-based denoising methods cannot avoid the problem of over-smoothing, because they always filter out high-frequency edge information along with the noise. This problem leads to the reduction of remote sensing image quality and negatively impacts downstream tasks. Hence, we tried to build a denoising model that can preserve edge details to obtain clear and high-quality remote sensing images. We note that the traditional LRR model only cares about keeping the overall image structure, but neglects the edges by filtering them out as part of the residues. Motivated by this phenomenon, we tried to greatly utilize the edge information in the residues to enrich the edge details of the denoised images, which can achieve further improvement of the low-rank model.
Inspired by the previous works on prior knowledge, it can be clearly observed that prior knowledge is useful in the aspect of smoothing noisy parts. In (1), WGLRR uses the ROAD statistic to design its prior knowledge matrix, which can distinguish between clean pixels and noisy pixels. On this basis, we modified the prior knowledge matrix in WGLRR to be more suitable for remote sensing images and integrated it into our model. Therefore, our model is formulated as
min Z , R Z * + λ R 2 , 1 s . t . P X = P Z + R
where P represents the prior knowledge matrix we designed for remote sensing images. R represents the low-rank residues. X and Z still have the same meaning as (1).
However, we cannot directly extract edge information from R, because there is much noise in it, which will affect the subsequent learning effect of the edge subspace. Then, we further propose a more relaxed prior knowledge to help better extract edge regions in the residues. By using this prior knowledge, it can ensure that noise and clean pixels are separated and more edge information is preserved in the clean part. Therefore, we introduce a new constraint term on the basis of (2).
min Z , R , R 1 , R 2 Z * + λ R 2 , 1 s . t . P 1 X = P 1 Z + R P 2 R = P 2 R 1 + R 2
Through this new term, the original residue R is separated into two parts: edge information R 1 and the noise part after extraction R 2 . Note that we used different P matrices in (3), which is a special design for our multi-level model. Here, we set two different prior matrices for different purposes. P 1 was employed to obtain the low-rank part from the noisy image, while P 2 was employed to find more edges from the residues. As we anticipated, P 1 should be stricter to ensure the noise will be removed as much as possible. The constraint of P 2 was relatively loose, because we wanted to use P 2 to extract more useful information from the residues, even though some noise will also be introduced. Through setting the two P matrices, we constructed a multi-level model, where the low-rank part and the edge part are obtained, respectively, at the corresponding level.
In order to learn the edge information we need, a manifold learning framework was then incorporated in our model, which can further learn the structural information from R 1 , so as to extract the edge subspace. Specifically, we tried to find the robust projection V = ( v 1 , v 2 , . . . , v r ) from R 1 and make sure that the geometric structural information of the data will be seized concurrently along the way.
min V , Z , W , R , R 1 , R 2 1 2 i , j V T r 1 , i V T r 1 , j 2 2 W i j + λ 3 W F 2 s . t . P 1 X = P 1 Z + R P 2 R = P 2 R 1 + R 2 V T V = I W 1 = 1 , W i i = 0 , W 0
where 1 2 i , j V T r 1 , i V T r 1 , j 2 2 is a coupling term integrated into our model, because it is impossible to predict which part of the data will be influenced by noise in advance. In short, if two similar samples r 1 , i and r 1 , j have very small values, they will mainly contain edge information, which also indicates that the distance between them will be quite small. Therefore, it is possible for us to extract those edge details from residue R. In contrast, the distance will be large if the noise part exceeds the edge information part. As a result, V is adaptively learned with similar residues through decreasing the effect of distant pairs by 1 2 i , j V T r 1 , i V T r 1 , j 2 2 . Put differently, a dynamic affinity graph W is integrated by utilizing similar points in R i , which is useful for the reconstruction of edge information. To achieve this goal, we set W i j as non-zero when two residue samples were similar. Otherwise, we set W i j as zero, attempting to mitigate the effect of noise. To effectuate our algorithm, we added two nonnegative constraints W i i = 0 and W 0 . W 1 = 1 was added into the model as well to guarantee that the affinity graph will be normalized.
After all, considering that the purpose of our method is to preserve the edge information as much as possible, we further added a TV operator to our model, which is widely used due to its outstanding performance on edge preservation. We can obtain the complete model:
min V , Z , W , R , R 1 , R 2 1 2 i , j V T r 1 , i V T r 1 , j 2 2 W i j + λ 1 Z * + λ 2 R 2 , 1 + λ 3 W F 2 + λ 4 Z T V s . t . P 1 X = P 1 Z + R P 2 R = P 2 R 1 + R 2 V T V = I W 1 = 1 , W i i = 0 , W 0
In order to solve the problem of (5), we propose an efficient optimization algorithm. Besides, we also introduced another variable F to avoid Sylvester’s equation appearing in the solution of R ^ . First, the following auxiliary terms S = Z , R 1 = R ^ , Z = T , and G = R were considered to make (5) easier to solve. Next, the coupling term 1 2 i , j V T r 1 , i V T r 1 , j 2 2 W i j was reformulated as t r ( V T ( F L R ^ T ) V ) , where L = d i a g i W i j W means the graph Laplacian of W. Here, d i a g means the operation of making a diagonal matrix. In addition, considering the existence of Sylvester’s equation, we introduced another variable F during the solution of R ^ . Consequently, the whole model can be described as
min V , Z , S , W , R , R 1 , R 2 , R ^ , F , T , G t r ( V T ( F L R ^ T ) V ) + λ 1 S * + λ 2 R 2 , 1 + λ 3 W F 2 + λ 4 T T V s . t . P 1 X = P 1 Z + R P 2 R = P 2 R 1 + R 2 R 1 = R ^ , S = Z , T = Z , G = R V T V = I W 1 = 1 , W i i = 0 , W 0

3.2. Design of Multi-Level Knowledge

In this section, we propose the special design of a novel multi-level knowledge regulation. We hope that our prior knowledge can distinguish between clean pixels and noise, so that the manifold learning introduced to our model can better learn the edge details we need. However, this requirement is usually put forward for optical images. Therefore, we attempted to achieve this goal through prior knowledge for remote sensing images. Such prior knowledge should simultaneously combine the characteristics of both optical and remote sensing images. We laid the foundation of our prior knowledge on the basis of existing methods. We have mentioned before that the ROAD statistic was introduced into a trilateral filter to reduce mixed noise. Thanks to that foundation, a new statistic called ROLD [40] has been introduced, which is able to identify contaminated pixels more accurately. ROLD makes an improvement to ROAD, which helps ROLD better distinguish between clean and noisy pixels. Then, Yu et al. [41] created a calculation for the uncorrupted probability of pixel q u .
P = exp ( 0.01 R O A D ( q u ) )
We combined (7) with the ROLD statistic to better distinguish the two types of pixels. Furthermore, we attempted to set an adjustable coefficient to meet our needs under multi-level modeling. Therefore, (7) can be written as
P = exp ( α R O L D ( q u ) )
where α represents the coefficient, helping gain various image features we need under distinct conditions. In our multi-level model, we set different α for different levels. At the first level, the α was large to easily find the difference between clean pixels and noisy pixels. This setting ensured that noise will be filtered out to the greatest extent. At the second level, the α was set to be small relatively so as to preserve the edge information we need. With different values of α , we can achieve different P matrices in (3).

3.3. Optimization Procedure

In this section, we try to tackle the problem we propose in (6). An efficient optimization algorithm was applied in our method. Each variable will be optimized one after another. Firstly, we obtain the augmented Lagrange function of (6):
min V , Z , T , S , W , R , R ^ , R 1 , R 2 , F , G t r ( V T ( F L R ^ T ) V ) + λ 1 S * + λ 2 R 2 , 1 + λ 3 W F 2 + λ 4 T T V + t r ( Y 1 T ( P 1 X P 1 Z R ) ) + t r ( Y 2 T ( P 2 G P 2 R 1 R 2 ) ) + t r ( Y 3 T ( S Z ) ) + t r ( Y 4 T ( R ^ R 1 ) ) + t r ( Y 5 T ( F R ^ ) + t r ( Y 6 T ( T Z ) + t r ( Y 7 T ( G R ) ) + μ 2 ( P 1 X P 1 Z R F 2 + P 2 G P 2 R 1 R 2 F 2 + S Z F 2 + R ^ R 1 F 2 + F R ^ F 2 + T Z F 2 + G R F 2 ) s . t . V T V = I W 1 = 1 , W i i = 0 , W 0
where Y 1 , Y 2 , Y 3 , Y 4 , Y 5 , Y 6 , and Y 7 are Lagrange multiples, which are used to deal with constrained optimization problems.
V subproblem:
This subproblem can be solved by removing the unrelated terms from (9):
min V t r ( V T ( F L R ^ T ) V ) s . t . V T V = I
Denoting F L R ^ T = Q , we can obtain
min V t r ( V T Q V ) s . t . V T V = I
Notice that (11) is an eigenvalue decomposition problem. To solve it, the set of k eigenvectors of the topmost k smallest eigenvalues of Q is needed.
Z subproblem:
By excluding other unrelated variables, we can simplify the optimal problem of Z as follows:
min Z P 1 X P 1 Z R + Y 1 μ F 2 + S Z + Y 3 μ F 2 + T Z + Y 6 μ F 2
Setting the derivative Z = 0 , the optimization procedure of Z can be described as:
Z = [ ( P 1 X P 1 ) ( R P 1 ) + ( Y 1 μ P 1 ) + S + Y 3 μ + T + Y 6 μ ] ÷ ( P 1 P 1 + 2 I )
where “÷” denotes elementwise division of two matrices.
T subproblem:
min T λ 4 T T V + μ 2 T Z + Y 6 μ F 2
To solve this subproblem, an iterative method was developed [42]. Firstly, U is defined as the set of matrix pairs ( u 1 , u 2 ) , u 1 R ( m 1 ) × n , u 2 R m × ( n 1 ) . Then, we can obtain
T = U ( ( Z Y 6 μ ) 2 λ 4 μ L ( u 1 , u 2 ) )
where L represents the linear operation, which can change the elements to the size of m × n . Finally, the strategy of fast gradient projection (FGP) [43] was used to solve this subproblem.
S subproblem:
min S λ 1 S * + μ 2 S Z + Y 3 μ F 2
By utilizing the singular-value thresholding (SVT), the optimization procedure of S can be described as
S = Φ λ 1 μ ( Z Y 3 μ )
where Φ represents a soft-thresholding operator.
W subproblem:
min W 1 2 i , j V T r 1 , i V T r 1 , j 2 2 W i j + λ 3 W F 2 s . t . W 1 = 1 , W i i = 0 , W 0
It can be found that (18) is independent for each i; therefore, the problem of W i could be rewritten as follows:
min W i 1 2 i , j V T r 1 , i V T r 1 , j 2 2 W i j + λ 3 W i F 2 s . t . W 1 = 1 , W i i = 0 , W 0
For simplicity, we used A i j to symbolize V T r 1 , i V T r 1 , j ; as a result, we have
min W i j W i + A i λ 3 2 2 s . t . W 1 = 1 , W i i = 0 , W 0
Furthermore, it can be unconstrained as
min W i j W i + A i λ 3 2 2 η ( 1 T W i 1 ) ζ T W i
Next, setting the derivative W i = 0 , we have
W i + 1 A i λ 3 η 1 ζ = 0
W i ’s j-th entry can be obtained through the following formula.
W i j + A i j λ 3 η ζ i = 0
Based on the KKT conditions, W i j can be updated as follows:
W i j = ( A i j λ 3 + η ) +
R subproblem:
min R λ 2 R 2 , 1 + μ 2 P 1 X P 1 Z R + Y 1 μ F 2 + μ 2 G R + Y 7 μ F 2
Denote T = 1 2 ( P 1 X P 1 Z + G + Y 1 μ + Y 7 μ ) , τ = λ 2 μ .
The i- t h column of R is given as follows:
R ( : , i ) = t i τ t i t i , i f τ < t i 0 , o t h e r w i s e
R ^ subproblem:
min R ^ t r ( V T ( F L R ^ T ) V ) + μ R ^ R 1 + Y 4 μ F 2 + μ F R ^ + Y 5 μ F 2
Setting the derivative R ^ = 0 , then the final result of the optimal R ^ can be described as follows:
R ^ = 2 μ R 1 V V T F L 2 Y 4 + 2 μ F + 2 Y 5 4 μ
R 1 subproblem:
min R 1 P 2 G P 2 R 1 R 2 + Y 2 μ F 2 + R ^ R 1 + Y 4 μ F 2
Setting the derivative R 1 = 0 , then the final result of the optimal R 1 can be described as follows:
R 1 = ( P 2 G R 2 + R ^ + Y 2 μ + Y 4 μ ) ÷ ( P 2 + 1 )
R 2 subproblem:
min R 2 P 2 G P 2 R 1 R 2 + Y 2 μ F 2
Setting the derivative R 2 = 0 , then we can obtain the optimal R 2 .
R 2 = ( P 2 G ) ( P 2 R 1 ) + Y 2 μ
F subproblem:
min F t r ( V T ( F L R ^ T ) V ) + μ F R ^ + Y 5 μ F 2
Setting the derivative F = 0 , then F is obtained as follows:
F = 2 μ R ^ P P T R ^ L T 2 Y 5 2 μ
G subproblem:
min G P 2 G P 2 R 1 R 2 + Y 2 μ F 2 + G R + Y 7 μ F 2
Setting the derivative G = 0 , then G is obtained as follows:
G = ( P 2 P 2 R 1 + P 2 R 2 P 2 Y 2 μ + R Y 7 μ ) ÷ ( 1 + P 2 P 2 )
To sum up, Algorithm 1 outlines the initialization of each variable and the multipliers.
Algorithm 1 The proposed EPLRR-RSID algorithm.
1:
Input: Data X i R d × n , parameter λ 1 , λ 2 , λ 3 , λ 4
2:
Initialize:
Z = S = W = R = R 1 = R 2 = R ^ = F = T = G = 0 ;
Y 1 = Y 2 = Y 3 = Y 4 = Y 5 = Y 6 = Y 7 = 0 ;
the projection V with orthogonal column vector;
ρ = 1.1 ; μ = 10 6 ; μ max = 10 6 ; ε = 10 5 ;
3:
while not converged do
4:
   Update V , Z , T , S , W , R , R ^ , R 1 , R 2 , F , G by (11), (13), (15), (17), (24), (26), (28), (30), (32), (34), and (36).
5:
   Update the multipliers Y 1 , Y 2 , Y 3 , Y 4 , Y 5 , Y 6 , Y 7 .
Y 1 = Y 1 + μ ( P 1 X P 1 Z R ) ,
Y 2 = Y 2 + μ ( P 2 G P 2 R 1 R 2 ) ,
Y 3 = Y 3 + μ ( S Z ) ,
Y 4 = Y 4 + μ ( R ^ R 1 ) ,
Y 5 = Y 5 + μ ( F R ^ ) ,
Y 6 = Y 6 + μ ( T Z ) ,
Y 7 = Y 7 + μ ( G R )
6:
   Update μ by μ = min ( ρ μ , μ m a x ) .
7:
   Check the following convergence conditions:
P 1 X P 1 Z R < ε , P 2 G P 2 R 1 R 2 < ε , S Z < ε , R ^ R 1 < ε , F R ^ < ε , T Z < ε , G R < ε
8:
end while
9:
Output: V R d × r and Z R d × n .

4. Results and Analysis

In this section, both simulated and real image data experiments were undertaken to demonstrate the effectiveness of our method. The results of EPLRR-RSID were compared with following denoising methods: LRTA [44], ANLM3D [45], BM4D [46], TDL [47], ITSReg [48], SNLRSF [49], and RCTV [50].

4.1. Experimental Settings

In the experiments of this paper, there were some parameters that needed to be determined. For the parameters in each algorithm, we set the values according to the optimal settings found in their respective literature works. The settings of parameters used in our experiments are shown in Table 1.
In the case of EPLRR-RSID, the parameters α 1 and α 2 in prior knowledge matrices P 1 and P 2 have a significant impact on the denoising performance and edge preserving performance of the denoised result, as we mentioned earlier. The specific setting of α 1 and α 2 depends on the different datasets. We used [0.001, 0.01, 0.1, 1] as the candidate set for selection. Furthermore, the parameters λ 1 , λ 2 , λ 3 , and λ 4 need tuning to find a robust solution under which the best denoising results will be obtained. However, it was hard for us to find a good combination of these parameters due to the complexity of the different combinations. To the best of our current knowledge, there are no accurate and efficient solutions for this problem in the existing literature. Thus, a reasonable strategy was used from [51] to seek the optimal combination of the parameters. The specific determination of these four parameters is discussed in the following section.

4.2. Experiments on Simulated Data

In this section, we chose one HSI dataset, the Pavia city center dataset (Available: http://www.ehu.es/ccwintco/index.php/Hyperspectral4_4Remote4_4Sensing4_4Scenes accessed on 2 November 2022.)), which is shown in Figure 2. The Pavia city center dataset was collected by the Reflective Optics System Imaging Spectrometer (ROSIS-03), whose size is 1096 × 715 × 102. However, we only used the subimage of size 200 × 200 × 80 (only clean bands). In the simulated experiments, we chose four indexes to assess the final denoised result: the peak signal-to-noise ratio (PSNR), the structural similarity (SSIM), the erreur relative globale adimensionnelle de synthese (ERGAS), and the mean spectral angle distance (MSAD). For an observed image whose size is m × n × b , the definitions are as follows:
P S N R = 255 2 × m × n i = 1 M j = 1 N I g ( i , j ) I ^ g ( i , j ) 2
where I g denotes the original image and I ^ g denotes the final result after denoising.
S S I M = ( 2 E ( I g ) E ( I g ) ) × ( 2 σ I I ^ ) ( ( E ( I g ) ) 2 + ( E ( I ^ g ) ) 2 ) ( σ I 2 + σ I ^ 2 )
where E ( I g ) and E ( I ^ g ) represents the mean value of I g and I ^ g . σ I 2 and σ I ^ 2 represent the variance of I g and I ^ g . σ I I ^ is the estimated covariance of image I g and I ^ g , respectively.
E R G A S = 1 b i = 1 b M S E ( I g i , I ^ g i ) E ( I g i )
where I g i and I ^ g i represent the i t h band of the original image and denoised image.
M S A D = 1 m n i = 1 m n 180 π × arccos ( χ i ) T · ( χ ^ i ) χ i · χ ^ i
where χ i and χ ^ i are the i t h spectral signatures of the original image and denoised image, respectively.
For the simulated process, additive Gaussian noise with different variance was added to the testing images to generate the noisy observations, whose values were 0.02, 0.04, 0.06, 0.08, and 0.1, respectively.
Figure 3 shows the denoised results of the Pavia city center image, which is contaminated by Gaussian noise with a noise variance value of 0.06. Figure 4 presents the magnified results from Figure 3. It is evident that our proposed method outperformed the others and efficiently reduced the Gaussian noise. ANLM3D obviously suffered from the problem of over-smoothing, which led to the loss of significant details in the image. In fact, one can barely recognize the basic outline of the denoised image. The result of the LRTA method still retained the most noise. The other five algorithms basically completed the denoising task and preserved enough edge details. However, we can further compare their denoising performance from the magnified results in Figure 4. From the enlarged local area, there were still unremoved noise points in the homogeneous area in the results of TDL and ITSReg. Although BM4D and SNLRSF showed good denoising performance in these regions, please note that these two methods both removed some obvious and important bright spots, while EPLRR-RSID well preserved these features and had the closest result to the original image.
Table 2 presents the four index values of the different compared methods for the Pavia city center dataset with different variances of Gaussian noise. It is easy to observe that our method EPLRR-RSID was more competitive in almost all cases. One can see that SNLRSF, RCTV, and EPLRR-RSID all achieved better results. Among these three methods, RCTV improved the denoising effect via representative coefficient total variation, while SNLRSF attached importance to spatial nonlocal self-similarity. Our method EPLRR-RSID took both the TV term and nonlocal self-similarity into account and further incorporated the extracted edge part, thereby leading to the superiorities of our method on all indexes.

4.3. Experiments on Real HSIs

In our experiments, we adopted two real HSI datasets to further demonstrate the effectiveness of the EPLRR-RSID method. The details of these two datasets are as follows:
1
AVIRIS Indian Pines Dataset (Available: http://www.ehu.es/ccwintco/index.php/Hyperspectral$_$Remote$_$Sensing$_$Scenes (accessed on 15 November 2022): This dataset was acquired by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) sensor over the Indian Pines test site in Northwestern Indiana in 1992. The size of the data is 145 × 145 × 220, and the false-color image is presented in Figure 5a. Several bands contained a mixture of water absorption, impulse noise, and Gaussian noise.
2
HYDICE Urban Dataset (Available: http://www.erdc.usace.army.mil/Media/Fact-Sheets/Fact-Sheet-Article-View/Article/610433/hypercube/ (accessed on 15 November 2022): This dataset was acquired by the HYDICE sensor, whose size is 307 × 307 × 210. This dataset includes roads, roofs, grass, and trees. A subimage of size 200 × 200 × 210 was chosen in our experiments, and the false-color image is presented in Figure 5b.
Figure 6 provides the denoising results of different methods on the Indian Pines dataset. It is clear that EPLRR-RSID enhanced the image quality to the utmost extent and preserved more fine details than the other compared algorithms. Among all these methods, we can see that the results of ANLM3D, BM4D, TDL, and ITSReg contained evident blurry areas, which led them to lose many important details. On the other hand, SNLRSF, RCTV, and our proposed EPLRR-RSID method presented denoising images with clear edge details. However, the result of RCTV was not satisfying in homogeneous regions, while SNLRSF and EPLRR-RSID both achieved a balance between noise suppression and edge preservation. This is because SNLRSF and EPLRR-RSID both make the best use of the low-rank characteristic and nonlocal self-similarity of HSIs, further improving the denoising performance of SNLRSF and EPLRR-RSID. In addition, EPLRR-RSID adds the edge part extracted from the low-rank residues, which will greatly enhance edge features in denoised results, giving a distinct advantage to our method.
Figure 7 shows the denoising results of different methods on the HYDICE Urban dataset. As displayed in Figure 7b, LRTA still could not effectively remove noise. Although BM4D can effectively remove some noise and stripes, the problem of over-smoothing still exists, as we feared, in Figure 7d. Remarkably, the denoised image of our method had the best visual effect, showing significant details to the fullest extent. We also present the vertical mean profiles of Band 139 of the HYDICE Urban dataset in Figure 8 to show the different performance of various methods.
As revealed in Figure 8, the curves of all methods have the effect of smoothing compared with Figure 8a. Furthermore, it is worth noting that EPLRR-RSID achieved the smoothest vertical mean profile. Such a phenomenon representing the noise in the original images was removed. Because the curves of the vertical mean profiles are typically believed to be smooth, whereas the complex noise in the original image will inevitably cause the rapid fluctuations in the curves of noisy image, the reduction of rapid fluctuations means that noise in the original image was effectively suppressed after denoising. Therefore, the smoothest curve obtained by EPLRR-RSID revealed that our method had the ability to achieve the best denoising performance in the HYDICE Urban experiment.
We have mentioned that our proposed method has the ability of edge preservation; therefore, the Edge Preservation Index (EPI) was selected to assess this ability of EPLRR-RSID. Using the EPI to assess the image quality, it is necessary to select a homogeneous region whose center is at ( i , j ) in the observed image firstly. In this region, the calculation of the EPI can be defined as
E P I = i M j N I ^ g ( i , j ) 8 × E ( I ^ g ) i M j N I g ( i , j ) 8 × E ( I g )
where I g , I ^ g , E ( I g ) , and E ( I ^ g ) have the same meaning as (37) and (38).
The EPI values of each band for two real HSI datasets are shown in Figure 9. It is obvious that our proposed EPLRR-RSID can achieve the highest EPI values in the vast majority of cases, manifesting the advantage of our method in edge preservation. Meanwhile, we can see that SNLRSF, which employs the nonlocal low-rank factorization, could not obtain satisfactory EPI results. This phenomenon showed that conventional low-rank-based methods have poor performances on edge preservation, while EPLRR-RSID further extracted the edge details to offset this defect.
We calculated the mean EPI values to further evaluate the performance of each algorithm on various HSI datasets. As shown in Figure 10, we can see that EPLRR-RSID had obvious advantages on both datasets. The reason is that our method further extracts edge information from low-rank residues, which made EPLRR-RSID contain more edge details than the other algorithms. Through adding the edge part R 1 , the ability of edge preservation will be greatly strengthened, although some noise will be introduced inevitably.
For real HSIs, there are no reference images for image quality assessment. In this paper, we adopted a no-reference image quality assessment named BRISQUE [52] to evaluate the two denoised real HSI datasets. BRISQUE employs scene statistics of locally normalized luminance coefficients to evaluate possible losses of naturalness in the image caused by the presence of distortions. It has been proven that this index is comparable with the PSNR and SSIM or even better. The BRISQUE results of different comparison algorithm are shown in Table 3; a lower value represents higher image quality. Thereby, it is obvious that our method had some advantages over all the compared methods on the whole.

4.4. Parameter Sensitivity

In this section, we discuss the selection of four hyperparameters in our methods, along with a test of their impacts on denoising performance. The sensitivity analysis was undertaken on the simulated Pavia city center dataset, of which the noise variance was 0.1. The results are shown in Figure 11 and Figure 12.
As shown in Figure 11a, the MPSNR value of EPLRR-RSID was higher when λ 1 and λ 2 changed in the range of 0.01–10 and 1–10. When fixing λ 1 and λ 2 , we can see that the results maintained stability as the λ 3 was within 0.001–1 and λ 4 was within 0.001–0.1. Furthermore, the MSSIM value was comparatively stable, no matter whether the value of λ 2 was high or low, but the results of the MSSIM were clearly better when λ 1 changed in the range of 0.01–1. When fixing λ 1 and λ 2 , the changes of λ 3 and λ 4 had little effect on the denoising results. In general, the change of the four parameters had relatively little impact on the MSSIM. This is because the value of the MSSIM itself fluctuates within a quite small range, and the change of parameters will not influence edge extraction greatly. Therefore, there are no significant trend changes in Figure 12.

4.5. Complexity Analysis

Due to the NLM framework we used in EPLRR-RSID, we only discuss the computational complexity on each group of data. We mainly took multiplication operations, nuclear norm minimization, the inverse of a matrix operation, and eigendecomposition into consideration. For one patch of the image, the multiplication operation, nuclear norm minimization, and the inverse of a matrix operation all cost O ( n 3 ) , while the eigendecomposition cost O ( d 2 r ) . Therefore, the computational complexity of EPLRR-RSID in each iteration can be written as O ( ( r + 1 ) n 3 + d 2 r ) .

5. Conclusions

In this paper, a novel denoising algorithm for remote sensing images named EPLRR-RSID was proposed. We attempted to take another perspective of low-rank residuals and preserve those edge details. Therefore, our proposed EPLRR-RSID method extracts the edge information from the residuals by finding robust image projections under a manifold learning framework. As our greatest advantage, our method had excellent performance in edge preservation. We also integrated a multi-level prior knowledge regulation, which helped distinguish pixels of different categories. Our multi-level modeling was proven to be efficient since more structural information can be learned. In this way, our method had more advantages than the existing methods that focus on edge preservation, because the edge information we extracted from the residues contained considerable details. The experimental results showed that EPLRR-RSID is highly effective on synthetic data and real HSI datasets. Furthermore, it was proven that EPLRR-RSID is better than existing denoising methods in edge preservation, indicating that it has the ability to better preserve edge information. More details can be observed and found in the denoised remote sensing images, which help better complete downstream tasks, such as target detection, classification, and clustering. This new method can be used for remote sensing image pre-processing to enhance the performance of subsequent processes. However, there is still a drawback in our method, which is that the running time of EPLRR-RSID depends on the size of the image. When processing large real images, our method was not fast enough. Future works include extending this method into the deep learning field to improve the denoising performance and computational efficiency.

Author Contributions

Conceptualization and validation, X.F. and S.T.; methodology, S.T. and X.S.; writing—original draft preparation and editing, X.F.; writing—review, S.E.A. and J.W.; project administration, S.T., H.Z., and C.W.; software, Z.X., X.Z., and W.G.; supervision, C.W.; funding acquisition, S.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the Key Program of National Natural Science Foundations of China, Grant Number 41930110.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Editor and the anonymous Referees for their instructive comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, P.; Huang, F.; Li, G.; Liu, Z. Remote-sensing image denoising using partial differential equations and auxiliary images as priors. IEEE Geosci. Remote Sens. Lett. 2011, 9, 358–362. [Google Scholar] [CrossRef]
  2. Liu, S.; Liu, M.; Li, P.; Zhao, J.; Zhu, Z.; Wang, X. SAR image denoising via sparse representation in shearlet domain based on continuous cycle spinning. IEEE Trans. Geosci. Remote Sens. 2017, 55, 2985–2992. [Google Scholar] [CrossRef]
  3. Li, L.; Hu, J.; Wu, F.; Zhao, J. A research and strategy of remote sensing image denoising algorithms. In Proceedings of the Advances in Natural Computation, Fuzzy Systems and Knowledge Discovery: Volume 2; Springer: Berlin/Heidelberg, Germany, 2020; pp. 704–712. [Google Scholar]
  4. Aurich, V.; Weule, J. Non-linear Gaussian filters performing edge preserving diffusion. In Proceedings of the Mustererkennung 1995: Verstehen Akustischer und Visueller Informationen; Springer: Berlin/Heidelberg, Germany, 1995; pp. 538–545. [Google Scholar]
  5. Paris, S.; Kornprobst, P.; Tumblin, J.; Durand, F. Bilateral filtering: Theory and applications. Found. Trends Comput. Graph. Vis. 2009, 4, 1–73. [Google Scholar] [CrossRef]
  6. He, K.; Sun, J.; Tang, X. Guided image filtering. IEEE Trans. Pattern Anal. Mach. Intell. 2012, 35, 1397–1409. [Google Scholar] [CrossRef] [PubMed]
  7. Buades, A.; Coll, B.; Morel, J.M. A nonlocal algorithm for image denoising. In Proceedings of the 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR’05), San Diego, CA, USA, 20–25 June 2005; IEEE: Piscataway, NJ, USA, 2005; Volume 2, pp. 60–65. [Google Scholar]
  8. Dabov, K.; Foi, A.; Katkovnik, V.; Egiazarian, K. Image denoising by sparse 3-D transform-domain collaborative filtering. IEEE Trans. Image Process. 2007, 16, 2080–2095. [Google Scholar] [CrossRef] [PubMed]
  9. Rudin, L.I.; Osher, S.; Fatemi, E. Nonlinear total variation based noise removal algorithms. Phys. D Nonlinear Phenom. 1992, 60, 259–268. [Google Scholar] [CrossRef]
  10. Wright, J.; Ganesh, A.; Rao, S.; Peng, Y.; Ma, Y. Robust principal component analysis: Exact recovery of corrupted low-rank matrices via convex optimization. Adv. Neural Inf. Process. Syst. 2009, 22. [Google Scholar] [CrossRef]
  11. Guan, D.; Xiang, D.; Tang, X.; Kuang, G. SAR image despeckling based on nonlocal low-rank regularization. IEEE Trans. Geosci. Remote Sens. 2018, 57, 3472–3489. [Google Scholar] [CrossRef]
  12. Chen, G.; Li, G.; Liu, Y.; Zhang, X.P.; Zhang, L. SAR image despeckling based on combination of fractional-order total variation and nonlocal low rank regularization. IEEE Trans. Geosci. Remote Sens. 2019, 58, 2056–2070. [Google Scholar] [CrossRef]
  13. He, W.; Zhang, H.; Zhang, L.; Shen, H. Hyperspectral image denoising via noise-adjusted iterative low-rank matrix approximation. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 3050–3061. [Google Scholar] [CrossRef]
  14. Ma, T.H.; Xu, Z.; Meng, D. Remote sensing image denoising via low-rank tensor approximation and robust noise modeling. Remote Sens. 2020, 12, 1278. [Google Scholar] [CrossRef]
  15. Zhang, H.; Liu, L.; He, W.; Zhang, L. Hyperspectral image denoising with total variation regularization and nonlocal low-rank tensor decomposition. IEEE Trans. Geosci. Remote Sens. 2019, 58, 3071–3084. [Google Scholar] [CrossRef]
  16. Li, C.; Ma, Y.; Huang, J.; Mei, X.; Ma, J. Hyperspectral image denoising using the robust low-rank tensor recovery. JOSA A 2015, 32, 1604–1612. [Google Scholar] [CrossRef]
  17. Yang, J.H.; Zhao, X.L.; Ma, T.H.; Chen, Y.; Huang, T.Z.; Ding, M. Remote sensing images destriping using unidirectional hybrid total variation and nonconvex low-rank regularization. J. Comput. Appl. Math. 2020, 363, 124–144. [Google Scholar] [CrossRef]
  18. Cao, W.; Chang, Y.; Han, G.; Li, J. Destriping remote sensing image via low-rank approximation and nonlocal total variation. IEEE Geosci. Remote Sens. Lett. 2018, 15, 848–852. [Google Scholar] [CrossRef]
  19. Sun, L.; Zhan, T.; Wu, Z.; Xiao, L.; Jeon, B. Hyperspectral mixed denoising via spectral difference-induced total variation and low-rank approximation. Remote Sens. 2018, 10, 1956. [Google Scholar] [CrossRef]
  20. Fan, H.; Li, C.; Guo, Y.; Kuang, G.; Ma, J. Spatial–spectral total variation regularized low-rank tensor decomposition for hyperspectral image denoising. IEEE Trans. Geosci. Remote Sens. 2018, 56, 6196–6213. [Google Scholar] [CrossRef]
  21. Zhang, M.; Desrosiers, C. Structure preserving image denoising based on low-rank reconstruction and gradient histograms. Comput. Vis. Image Underst. 2018, 171, 48–60. [Google Scholar] [CrossRef]
  22. Wang, M.; Wang, Q.; Chanussot, J. L 0 Gradient Regularized Low-Rank Tensor Model for Hyperspectral Image Denoising. In Proceedings of the 2019 10th Workshop on Hyperspectral Imaging and Signal Processing: Evolution in Remote Sensing (WHISPERS), Amsterdam, The Netherlands, 24–26 September 2019; IEEE: Piscataway, NJ, USA, 2019; pp. 1–6. [Google Scholar]
  23. Liu, P.; Wang, M.; Wang, L.; Han, W. Remote-sensing image denoising with multi-sourced information. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2019, 12, 660–674. [Google Scholar] [CrossRef]
  24. Han, L.; Zhao, Y.; Lv, H.; Zhang, Y.; Liu, H.; Bi, G. Remote sensing image denoising based on deep and shallow feature fusion and attention mechanism. Remote Sens. 2022, 14, 1243. [Google Scholar] [CrossRef]
  25. Wang, Z.; Ng, M.K.; Zhuang, L.; Gao, L.; Zhang, B. Nonlocal self-similarity-based hyperspectral remote sensing image denoising with 3-d convolutional neural network. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–17. [Google Scholar] [CrossRef]
  26. Zeng, H.; Xie, X.; Cui, H.; Yin, H.; Ning, J. Hyperspectral image restoration via global L 1-2 spatial–spectral total variation regularized local low-rank tensor recovery. IEEE Trans. Geosci. Remote Sens. 2020, 59, 3309–3325. [Google Scholar] [CrossRef]
  27. Wang, M.; Wang, Q.; Chanussot, J. Tensor low-rank constraint and l_0 total variation for hyperspectral image mixed noise removal. IEEE J. Sel. Top. Signal Process. 2021, 15, 718–733. [Google Scholar] [CrossRef]
  28. Sun, L.; He, C. Hyperspectral Image Mixed Denoising Using Difference Continuity-Regularized Nonlocal Tensor Subspace Low-Rank Learning. IEEE Geosci. Remote Sens. Lett. 2021, 19, 1–5. [Google Scholar] [CrossRef]
  29. Lei, S.; Lu, M.; Lin, J.; Zhou, X.; Yang, X. Remote sensing image denoising based on improved semi-soft threshold. Signal Image Video Process. 2021, 15, 73–81. [Google Scholar] [CrossRef]
  30. Das, S.; Chakravortty, S. Spectral-spatial 3D dynamic trimmed median filter for removal of impulse noise in remotely sensed images. Multimed. Tools Appl. 2022, 82, 15945–15982. [Google Scholar] [CrossRef]
  31. Pandey, S.; Miri, R.; Sinha, G.; Raja, R. AFD filter and E2N2 classifier for improving visualization of crop image and crop classification in remote sensing image. Int. J. Remote Sens. 2022, 43, 5848–5873. [Google Scholar] [CrossRef]
  32. Geng, J.; Fan, J.; Ma, X.; Wang, H.; Cao, K. An iterative low-rank representation for SAR image despeckling. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; IEEE: Piscataway, NJ, USA, 2016; pp. 72–75. [Google Scholar]
  33. Xue, J.; Zhao, Y.; Liao, W.; Kong, S.G. Joint spatial and spectral low-rank regularization for hyperspectral image denoising. IEEE Trans. Geosci. Remote Sens. 2017, 56, 1940–1958. [Google Scholar] [CrossRef]
  34. Chang, Y.; Yan, L.; Zhong, S. Hyper-laplacian regularized unidirectional low-rank tensor recovery for multispectral image denoising. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Hawaii, USA, 21–26 July 2017; pp. 4260–4268. [Google Scholar]
  35. Chen, F.; Zhang, L.; Yu, H. External patch prior guided internal clustering for image denoising. In Proceedings of the IEEE International Conference on Computer Vision, Santiago, Chile, 7–13 December 2015; pp. 603–611. [Google Scholar]
  36. Xie, M.; Liu, X.; Yang, X. Novel hybrid low-rank tensor approximation for hyperspectral image mixed denoising based on global-guided-nonlocal prior mechanism. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–17. [Google Scholar] [CrossRef]
  37. Tang, C.; Cao, L.; Chen, J.; Zheng, X. Speckle noise reduction for optical coherence tomography images via nonlocal weighted group low-rank representation. Laser Phys. Lett. 2017, 14, 056002. [Google Scholar] [CrossRef]
  38. Garnett, R.; Huegerich, T.; Chui, C.; He, W. A universal noise removal algorithm with an impulse detector. IEEE Trans. Image Process. 2005, 14, 1747–1754. [Google Scholar] [CrossRef]
  39. Guo, Q.; Zhang, C.; Zhang, Y.; Liu, H. An efficient SVD-based method for image denoising. IEEE Trans. Circuits Syst. Video Technol. 2015, 26, 868–880. [Google Scholar] [CrossRef]
  40. Dong, Y.; Chan, R.H.; Xu, S. A detection statistic for random-valued impulse noise. IEEE Trans. Image Process. 2007, 16, 1112–1120. [Google Scholar] [CrossRef]
  41. Yu, H.; Gao, J.; Li, A. Probability-based nonlocal means filter for speckle noise suppression in optical coherence tomography images. Opt. Lett. 2016, 41, 994–997. [Google Scholar] [CrossRef]
  42. Wang, H.; Cen, Y.; He, Z.; He, Z.; Zhao, R.; Zhang, F. Reweighted low-rank matrix analysis with structural smoothness for image denoising. IEEE Trans. Image Process. 2017, 27, 1777–1792. [Google Scholar] [CrossRef] [PubMed]
  43. Beck, A.; Teboulle, M. Fast gradient-based algorithms for constrained total variation image denoising and deblurring problems. IEEE Trans. Image Process. 2009, 18, 2419–2434. [Google Scholar] [CrossRef]
  44. Renard, N.; Bourennane, S.; Blanc-Talon, J. Denoising and dimensionality reduction using multilinear tools for hyperspectral images. IEEE Geosci. Remote Sens. Lett. 2008, 5, 138–142. [Google Scholar] [CrossRef]
  45. Manjón, J.V.; Coupé, P.; Martí-Bonmatí, L.; Collins, D.L.; Robles, M. Adaptive nonlocal means denoising of MR images with spatially varying noise levels. J. Magn. Reson. Imaging 2010, 31, 192–203. [Google Scholar] [CrossRef]
  46. Maggioni, M.; Katkovnik, V.; Egiazarian, K.; Foi, A. Nonlocal transform-domain filter for volumetric data denoising and reconstruction. IEEE Trans. Image Process. 2012, 22, 119–133. [Google Scholar] [CrossRef]
  47. Peng, Y.; Meng, D.; Xu, Z.; Gao, C.; Yang, Y.; Zhang, B. Decomposable nonlocal tensor dictionary learning for multispectral image denoising. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Columbus, OH, USA, 23–28 June 2014; pp. 2949–2956. [Google Scholar]
  48. Xie, Q.; Zhao, Q.; Meng, D.; Xu, Z.; Gu, S.; Zuo, W.; Zhang, L. Multispectral images denoising by intrinsic tensor sparsity regularization. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Las Vegas, NV, USA, 27–30 June 2016; pp. 1692–1700. [Google Scholar]
  49. Cao, C.; Yu, J.; Zhou, C.; Hu, K.; Xiao, F.; Gao, X. Hyperspectral image denoising via subspace-based nonlocal low-rank and sparse factorization. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2019, 12, 973–988. [Google Scholar] [CrossRef]
  50. Peng, J.; Wang, H.; Cao, X.; Liu, X.; Rui, X.; Meng, D. Fast Noise Removal in Hyperspectral Images via Representative Coefficient Total Variation. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–17. [Google Scholar] [CrossRef]
  51. Wen, J.; Xu, Y.; Liu, H. Incomplete multiview spectral clustering with adaptive graph learning. IEEE Trans. Cybern. 2018, 50, 1418–1429. [Google Scholar] [CrossRef] [PubMed]
  52. Mittal, A.; Moorthy, A.K.; Bovik, A.C. No-reference image quality assessment in the spatial domain. IEEE Trans. Image Process. 2012, 21, 4695–4708. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The framework of EPLRR-RSID. As illustrated using the above corrupted scenario, image edges are preserved by incorporating the multi-level prior knowledge matrices. Through the two-level model, we can extract clear edges from low-rank residues.
Figure 1. The framework of EPLRR-RSID. As illustrated using the above corrupted scenario, image edges are preserved by incorporating the multi-level prior knowledge matrices. Through the two-level model, we can extract clear edges from low-rank residues.
Remotesensing 15 02318 g001
Figure 2. Pavia city center dataset (R: 80, G: 42, B: 8).
Figure 2. Pavia city center dataset (R: 80, G: 42, B: 8).
Remotesensing 15 02318 g002
Figure 3. Denoised results in the simulated experiments under the noise variance 0.06: (a) The original image of Band 59. (b) noisy band. (c) LRTA. (d) ANLM3D. (e) BM4D. (f) TDL. (g) ITSReg. (h) SNLRSF. (i) RCTV. (j) EPLRR-RSID.
Figure 3. Denoised results in the simulated experiments under the noise variance 0.06: (a) The original image of Band 59. (b) noisy band. (c) LRTA. (d) ANLM3D. (e) BM4D. (f) TDL. (g) ITSReg. (h) SNLRSF. (i) RCTV. (j) EPLRR-RSID.
Remotesensing 15 02318 g003
Figure 4. Magnified results of Figure 3: (a) The original image of Band 59. (b) Noisy band. (c) LRTA. (d) ANLM3D. (e) BM4D. (f) TDL. (g) ITSReg. (h) SNLRSF. (i) RCTV. (j) EPLRR-RSID.
Figure 4. Magnified results of Figure 3: (a) The original image of Band 59. (b) Noisy band. (c) LRTA. (d) ANLM3D. (e) BM4D. (f) TDL. (g) ITSReg. (h) SNLRSF. (i) RCTV. (j) EPLRR-RSID.
Remotesensing 15 02318 g004
Figure 5. Two real HSI datasets. (a) AVIRIS Indian Pines (R: 29, G: 42, B: 89). (b) HYDICE Urban (R: 2, G: 109, B: 207).
Figure 5. Two real HSI datasets. (a) AVIRIS Indian Pines (R: 29, G: 42, B: 89). (b) HYDICE Urban (R: 2, G: 109, B: 207).
Remotesensing 15 02318 g005
Figure 6. Denoised results of Indian Pines dataset. (a) The original image of Band 3. (b) LRTA. (c) ANLM3D. (d) BM4D. (e) TDL. (f) ITSReg. (g) SNLRSF. (h) RCTV. (i) EPLRR-RSID.
Figure 6. Denoised results of Indian Pines dataset. (a) The original image of Band 3. (b) LRTA. (c) ANLM3D. (d) BM4D. (e) TDL. (f) ITSReg. (g) SNLRSF. (h) RCTV. (i) EPLRR-RSID.
Remotesensing 15 02318 g006
Figure 7. Denoised results of HYDICE Urban dataset. (a) The original image of Band 139. (b) LRTA. (c) ANLM3D. (d) BM4D. (e) TDL. (f) ITSReg. (g) SNLRSF. (h) RCTV. (i) EPLRR-RSID.
Figure 7. Denoised results of HYDICE Urban dataset. (a) The original image of Band 139. (b) LRTA. (c) ANLM3D. (d) BM4D. (e) TDL. (f) ITSReg. (g) SNLRSF. (h) RCTV. (i) EPLRR-RSID.
Remotesensing 15 02318 g007
Figure 8. Vertical mean profiles of HYDICE Urban dataset. (a) Original. (b) LRTA. (c) ANLM3D. (d) BM4D. (e) TDL. (f) ITSReg. (g) SNLRSF. (h) RCTV. (i) EPLRR-RSID. The selected band in this experiment is Band 139.
Figure 8. Vertical mean profiles of HYDICE Urban dataset. (a) Original. (b) LRTA. (c) ANLM3D. (d) BM4D. (e) TDL. (f) ITSReg. (g) SNLRSF. (h) RCTV. (i) EPLRR-RSID. The selected band in this experiment is Band 139.
Remotesensing 15 02318 g008
Figure 9. EPI values of each band for two real HSI datasets. (a) Indian Pines dataset. (b) HYDICE Urban dataset.
Figure 9. EPI values of each band for two real HSI datasets. (a) Indian Pines dataset. (b) HYDICE Urban dataset.
Remotesensing 15 02318 g009
Figure 10. Mean EPI values of different methods on real HSI datasets. The corresponding values are also marked on the figure.
Figure 10. Mean EPI values of different methods on real HSI datasets. The corresponding values are also marked on the figure.
Remotesensing 15 02318 g010
Figure 11. The parameter sensitivity study in terms of MPSNR. (a) Fix λ 3 , λ 4 ( λ 3 = 1, λ 4 = 1) to adjust λ 1 and λ 2 . (b) Fix λ 1 and λ 2 ( λ 1 = 0.01, λ 2 = 1) to adjust λ 3 and λ 4 .
Figure 11. The parameter sensitivity study in terms of MPSNR. (a) Fix λ 3 , λ 4 ( λ 3 = 1, λ 4 = 1) to adjust λ 1 and λ 2 . (b) Fix λ 1 and λ 2 ( λ 1 = 0.01, λ 2 = 1) to adjust λ 3 and λ 4 .
Remotesensing 15 02318 g011
Figure 12. The parameter sensitivity study in terms of MSSIM. (a) Fix λ 3 , λ 4 ( λ 3 = 1, λ 4 = 1) to adjust λ 1 and λ 2 . (b) Fix λ 1 and λ 2 ( λ 1 = 1, λ 2 = 1) to adjust λ 3 and λ 4 .
Figure 12. The parameter sensitivity study in terms of MSSIM. (a) Fix λ 3 , λ 4 ( λ 3 = 1, λ 4 = 1) to adjust λ 1 and λ 2 . (b) Fix λ 1 and λ 2 ( λ 1 = 1, λ 2 = 1) to adjust λ 3 and λ 4 .
Remotesensing 15 02318 g012
Table 1. The settings of parameters used in our experiments.
Table 1. The settings of parameters used in our experiments.
DefinitionNotationValue
Patch size P s 9 × 9
The number of similar patches in each group matrix L s 81
The parameters in prior knowledge matrices α 1 , α 2 [0.001, 0.01, 0.1, 1]
Hyperparameters λ 1 , λ 2 , λ 3 , λ 4 [0.0001, 0.01, …, 1000]
Table 2. Numerical evaluation for the simulated Pavia city center image.
Table 2. Numerical evaluation for the simulated Pavia city center image.
Noise VarianceEvaluation IndexLRTAANLM3DBM4DTDLITSRegSNLRSFRCTVEPLRR-RSID
0.02MPSNR42.051540.938643.249445.502143.494746.395837.347748.7693
MSSIM0.98750.98080.99170.97080.99500.99240.99300.9943
ERGAS29.835233.858225.363621.042624.845818.431352.407018.0475
MSAD0.06590.05320.05330.03880.03850.03580.07190.0336
0.04MPSNR38.654933.991538.983541.487339.965842.666436.335145.8265
MSSIM0.97530.91890.97860.98780.98330.99080.96430.9927
ERGAS43.477773.906041.300931.601137.190329.005357.077327.3846
MSAD0.08220.07290.07840.05280.04950.05120.07730.0428
0.06MPSNR36.220430.809136.530438.849337.837940.266834.687743.1284
MSSIM0.95880.85530.96300.97880.97380.98470.94900.9851
ERGAS57.5013105.628354.698342.553247.448438.555569.200735.1563
MSAD0.10030.08450.09560.06260.05750.06380.08750.0532
0.08MPSNR34.329029.070734.801037.034236.461038.473233.944040.2722
MSSIM0.93900.80280.94590.96920.96450.97770.93920.9806
ERGAS71.1996128.525566.693452.332055.465547.374675.493943.2758
MSAD0.11780.09110.10950.07190.06370.07460.09290.0639
0.1MPSNR33.099528.003633.469935.589735.558236.654933.109537.9364
MSSIM0.92250.76240.92780.95900.95720.96310.92900.9718
ERGAS81.6941145.046477.698261.801461.564658.286381.903756.7236
MSAD0.12430.09580.12180.07870.06870.07950.09820.0675
The best values are highlighted in bold.
Table 3. BRISQUE value of different denoising methods on the two real HSI datasets.
Table 3. BRISQUE value of different denoising methods on the two real HSI datasets.
LRTAANLM3DBM4DTDLITSRegSNLRSFRCTVEPLRR-RSID
Indian Pines114.97115.82115.94114.84116.1695.47116.2391.60
HYDICE Urban96.1199.26101.2696.16100.5697.56101.1486.44
The best values are highlighted in bold.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Feng, X.; Tian, S.; Abhadiomhen, S.E.; Xu, Z.; Shen, X.; Wang, J.; Zhang, X.; Gao, W.; Zhang, H.; Wang, C. Edge-Preserved Low-Rank Representation via Multi-Level Knowledge Incorporation for Remote Sensing Image Denoising. Remote Sens. 2023, 15, 2318. https://doi.org/10.3390/rs15092318

AMA Style

Feng X, Tian S, Abhadiomhen SE, Xu Z, Shen X, Wang J, Zhang X, Gao W, Zhang H, Wang C. Edge-Preserved Low-Rank Representation via Multi-Level Knowledge Incorporation for Remote Sensing Image Denoising. Remote Sensing. 2023; 15(9):2318. https://doi.org/10.3390/rs15092318

Chicago/Turabian Style

Feng, Xiaolin, Sirui Tian, Stanley Ebhohimhen Abhadiomhen, Zhiyong Xu, Xiangjun Shen, Jing Wang, Xinming Zhang, Wenyun Gao, Hong Zhang, and Chao Wang. 2023. "Edge-Preserved Low-Rank Representation via Multi-Level Knowledge Incorporation for Remote Sensing Image Denoising" Remote Sensing 15, no. 9: 2318. https://doi.org/10.3390/rs15092318

APA Style

Feng, X., Tian, S., Abhadiomhen, S. E., Xu, Z., Shen, X., Wang, J., Zhang, X., Gao, W., Zhang, H., & Wang, C. (2023). Edge-Preserved Low-Rank Representation via Multi-Level Knowledge Incorporation for Remote Sensing Image Denoising. Remote Sensing, 15(9), 2318. https://doi.org/10.3390/rs15092318

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop