Superpixel-Based Mixed Noise Estimation for Hyperspectral Images Using Multiple Linear Regression

HSIs (hyperspectral images) obtained by new-generation hyperspectral sensors contain both electronic noise and photon noise with comparable power. Therefore, both the SI (signal-independent) component and the SD (signal-dependent) component have to be considered. In this paper, a superpixel-based noise estimation algorithm using MLR (multiple linear regression) is proposed for the above mixed noise to estimate the noise standard deviation of both SI component and SD component. First, superpixel segmentation is performed on the first principal component obtained by MNF (minimum noise fraction)-based dimensionality reduction to generate non-overlapping regions with similar pixels. Then, MLR is performed to remove the spectral correlation, and a system of linear equations with respect to noise variances is established according to the local sample statistics calculated within each superpixel. By solving the equations in terms of the least-squares method, the noise variances are determined. The experimental results show that the proposed algorithm provides more accurate local sample statistics, and yields a more accurate noise estimation than the other state-of-the-art algorithms for simulated HSIs. The results of the real-life data also verify the effectiveness of the proposed algorithm.


Introduction
Emerging from the development of hyperspectral remote-sensing technology, HSIs (hyperspectral images) have both high spatial resolution and spectral resolution. The high resolution of HSIs makes them widely used in agriculture, forestry, geological survey, environmental detection, military reconnaissance and other fields. However, in the process of the acquisition and transmission of HSIs, noise is inevitably introduced which has many negative effects on target detection, spectral matching, above-ground estimation [1], classification [2] and other subsequent applications of hyperspectral remote-sensing images. For example, to realize a full-spectrum denoising of HSIs with a high SNR (signal-to-noise ratio), it is necessary to estimate the standard deviation of the noise for the post processing PCA (principal component analysis) and dual-tree complex wavelet transform thresholding [3]. Therefore, quantitative analysis of the noise in HSIs is of great significance to the post-processing of HSI data.
The noise of HSIs can generally be divided into two types: systematic noise and random noise. The systematic noise is generated by the imperfect calibration of detectors, such as the stripes, which can be effectively removed by the state-of-the-art method [4][5][6][7]. The random noise is composed of electronic noise and photon noise. The electronic noise is an additive noise generated by electronic circuitry, which is independent of the signal, while the photon noise is generated by the Poisson-distributed number of photons, which is dependent on the signal. Because of the complex process of noise generation and its statistical randomness, it is difficult to estimate the noise level of HSIs accurately.
In former remote sensors, the electronic noise is dominant. The SI (signal-independent) electronic noise contributes most of the noise. Thus the hypothesis of the additive and stationary noise model is widely accepted. Under this hypothesis, the photon noise is supposed to be negligible. Therefore, in previous work [8][9][10][11][12], the noise model is based on the SI-additive Gaussian white noise. In this model, the noise variances vary with the wavelength, while the noise variance of each band is constant. Gao estimates the noise variance for each band by taking account into the homogeneous regions of the image [8]. Exploiting the spectral and spatial correlation for noise estimation, Roger et al., perform MLR (multiple linear regression) in small blocks of the same size, and the average of the variances of the residuals for all blocks is taken as the estimated variance for each band [9]. Since this estimation is sensitive to a heterogeneous subset in a small block, Aiazzi et al., focus on estimating the information conveyed to a user by HSIs, and propose an estimation algorithm under the assumption that the variance of the observed signal measured on homogeneous areas will be equal to the variance of the noise [10]. Gao et al., present a noise estimation method based on MLR in homogeneous regions detected by an object-seeking algorithm [11]. Unlike the above spatial decorrelation in homogeneous regions, Xu et al., estimate the variance of the noise in the wavelet domain. The above noise estimation methods achieve good performance for the SI noise, but are not able to estimate the variances of the SD (signal-dependent) noise [12].
With the development of modern electronic devices, an increasing number of studies have shown that the SI electronic noise can be suppressed well. When the power of the electronic noise and the power of the photon noise become competitive, the SD photon noise component cannot be ignored as in previous works. Instead of an SI noise model, a more appropriate and practical Poissonian-Gaussian noise model is built, where the noise is composed of Poissonian SD noise and Gaussian SI noise [13]. The Poissonian process is treated as a special heteroskedastic Gaussian approximation, and the noise parameters are estimated by wavelet transform and maximum likelihood estimation. Yang et al., employ a particle filter to detect homogenous blocks grouped together to estimate noise levels by solving sparse representation under training [14]. Rakhshanfar et al., classify image patches according to their intensity and variance to find homogenous regions that represent the noise [15], and then clusters of connected patches are weighted and ranked to approximate the peak noise variance and noise level function. Zhang et al., detect homogeneous regions in wavelet transformed blocks, and combine them together to create a larger sample set for the variance estimation of mixed Poissonian-Gaussian noise [16]. Li et al., select homogenous blocks via local gray statistic entropy [17]. Haar wavelet-based local median absolute deviation and maximum likelihood estimation are applied to the homogenous blocks to estimate the noise parameters. Although the above proposed methods can estimate the noise level accurately, they are designed for the two-dimensional natural images taken by ordinary digital cameras, and do not consider the spectral correlation in HSIs.
Alparone et al., propose an estimation method for the mixed noise model of HSIs collected by new-generation imaging spectrometers [18]. The local expectation and variance pairs in homogenous regions are viewed as scatter-points clustered along a straight line, whose slope and intercept measure the variance of the SD noise and the SI noise, respectively. This process is performed on the HSIs band by band. To consider the significant spectral correlation in HSIs, MLR is performed on the HSIs to separate the noise and signal realizations, and the variances of the SD noise and SI noise are estimated by a maximum-likelihood approach in the noise realizations [19]. To make use of both the spectral correlation and the spatial correlation in HSIs, several mixed noise estimation methods are proposed. Fu et al., detect the homogenous blocks by data masking followed by MLR on each block, and the noise parameters are calculated by the maximum-likelihood estimation approach [20]. Sun divides HSIs into non-overlapping blocks to explore local statistics, and employs MLR to remove the strong Remote Sens. 2020, 12, 1324 3 of 27 spectral correlation. By applying the statistic properties of the local statistics and MLR on the HSI, a block-based system of linear equations with respect to the SD and SI noise variances is established and solved [21].
Currently, the existing parameter estimation methods for the mixed noise of HSIs are not mature and developed enough, especially for HSIs with rich details such as edges and textures. These methods tend to confuse the details of HSIs with noise, which leads to an over-estimation of the noise level, especially when the noise level is high. In order to overcome the above shortcomings, segmentation-based methods are studied to find homogenous regions, which can efficiently reduce the influence of the edges and textures on noise estimation. Although segmentation and classification have been applied in some previous works [13][14][15], these methods are sensitive to noise, and incorrect segmentation will decrease estimation accuracy.
Based on the above analysis, we propose a superpixel-based mixed noise estimation method for HSIs using MLR. The superpixel segmentation technique is introduced to generate the homogeneous regions; in order to reduce the influence of noise on the segmentation effect, the MNF (minimum noise fraction) is employed to obtain the PC (principle component) with the highest quality. Then, MLR is performed to remove the spectral correlation. The system of linear equations with respect to the noise variances is established according to the local sample statistics calculated within each superpixel. By solving the equations based on the least-squares method, the noise variances can be estimated.
The rest of the paper is organized as follows. Section 2 gives the mixed noise model and describes a two-stage scheme for the estimation algorithm. Section 3 describes the MNF-based superpixel segmentation particularly for noisy HSIs. Section 4 presents the mixed noise estimation on superpixels. Section 5 gives the simulated and real-life data experimental results and discussion, and Section 6 offers conclusions.

Mixed Noise Model for Hyperspectral Images (HSI)
The mixed noise model for a given pixel in an image can be described as follows: where g is the observed noisy signal, and f is the noise-free signal. u is a stationary random process independent of f , with zero-mean and variance σ 2 u . Then the noise term v = f γ u is dependent on signal f , which is the SD noise. γ ≤ 1 is held for a great variety of images collected by several different imaging systems. For hyperspectral sensors, γ is equal to 0.5 [18][19][20][21]. w is a Gaussian white process independent of f , with zero-mean and variance σ 2 w , which is the SI noise. According to [13], the noise level function is defined as: where σ 2 T is the total variance of the noise f u + w, and I is the average intensity of the noise-free signal f .
Meanwhile, a total variance of the noise can be calculated by: where M, N are the height and width of the image, respectively, P N is the power of the noise term f u + w. Correspondingly, a total standard deviation σ T is defined as: and it holds that: The HSIs are regarded as a stack of images of L bands. The dataset in each band can be viewed as a static image, contaminated by mixed noise containing the SD noise and the SI noise. The SD noise variance σ 2 u and the SI noise variance σ 2 w are the two parameters to be estimated. Note that in different bands, the SD noise variance σ 2 u and the SI noise variance σ 2 w are different and vary with the wavelength.

Algorithm of the Noise Estimation
The estimation of mixed noise for HSI is incorporated into a two-stage scheme.
Step 1 Pre-processing to find the homogenous region. To reduce the effect of noise on segmentation, MNF-based dimensionality reduction is employed, and then superpixel segmentation is performed on the PC with the highest SNR to find the homogenous regions, wherein a superpixel is viewed as a homogenous region.
Step 2 Superpixel-based noise estimation. MLR is performed on the noisy data. For each band, the predicted values are obtained as linear combination of the pixels in the most adjacent bands. The coefficients of the linear combination for each band is obtained by the least-squares method. Using these coefficients and the local statistics in each superpixel, a system of linear equations are established with respect to the noise variances. By solving the above linear equations, the estimated parameters can be determined.

Superpixel Segmentation
The superpixel segmentation technique has been widely utilized to generate homogeneous regions for HSIs without any prior information, and has been a popular way to exploit the spatial and spectral correlations of HSIs. Pixels within each region have similar spatial and spectral signatures; this kind of stationary property is very helpful for the following noise estimation. Superpixel segmentation is considered in this paper to accurately estimate the noise levels from homogeneous regions in noisy HSIs. This method can segment an image into small regions that can fully adapt to the local structures. Note that the existing superpixel segmentation algorithms are widely used in a gray image or color image, however, unlike a gray or color image, HSI cannot be directly used for superpixel segmentation since the HSIs tested in this paper generally contain hundreds of bands with heavy noise. If the superpixel segmentation algorithm is performed on one of the bands directly, only a small amount of the local information for the HSI will be used. The presence of noise may also reduce the segmentation performance. In order to overcome this problem, a dimensionality reduction technique needs to be introduced for noisy HSIs. Both the PCA and MNF are widely used for dimensionality reduction [22]. The PCA orders all PCs according to their variances, and each succeeding component has the highest variance possible under the constraint that it is orthogonal to the preceding components. However, the noise variance may contribute to the variance of the succeeding components and, therefore, the variance is not perfectly able to measure image quality. Compared to variance, the SNR is a more effective evaluation criterion introduced by MNF. MNF can reorder all principal components according to SNR, the first principal component has higher image quality and is more conducive to subsequent processing. The MNF algorithm is designed based on the assumption that the noise is SI [23], and it is also effective on SD noise. Guangchun Luo et. al. report that their experiments demonstrate that MNF-based denoiseing methods achieve higher signal-to-noise ratios than PCA-based methods for SD noise [24]. Therefore, the MNF algorithm is suitable in the case of mixed SI and SD noise. In this paper, MNF-based dimensionality reduction is introduced to obtain the desired image with a high SNR for the following superpixel segmentation.

MNF-Based Dimensionality Reduction
The basic idea of MNF is to perform a noise-whitening process on the original HSIs, followed by applying PCA to the noise-whitened data. Because the noise variance is normalized after noise-whitening, the PCs are sorted according to their SNR values. Let P be the noise-whitening matrix, and let Σ n and Σ K be the noise covariance matrix and the covariance matrix of the noisy HSIs, respectively. Then we have: where I is an identity matrix. The covariance matrix after noise whitening can by obtained by the following operation: where Σ MNF is the covariance matrix after noise whitening. Then eigenvector U can be obtained by using eigenvalue decomposition for Σ MNF . Finally, the PCs achieved by MNF can be given as: It can be seen from the above calculation process that the key problem to perform MNF is estimating an accurate noise covariance matrix Σ n . However, this process is difficult. In [23,25] the estimated noise covariance matrix is given as follows: where η 2 k L k=1 is the diagonal elements of Σ −1 k , and Σ −1 k is the inverse matrix of Σ k . Once the MNF is performed on noisy HSIs, the first PC with the highest SNR will be retained for the following superpixel segmentation. Compared with arbitrary band of noisy HSIs, the first PC has three main advantages. First, it contains most information of HSIs. Second, the noise in the first PC may be suppressed. Finally, computational complexity can be significantly reduced for the superpixel segmentation process. Therefore, by using the first PC obtained by the MNF, the following superpixel segmentation is expected to achieve good performance. Algorithm 1 is the pseudo code of the MNF on HSI:

Superpixel Segmentation
Superpixels can be regarded as a set of non-overlapping regions with similar texture, brightness and other characteristics. The similarity of features between pixels is used to regroup pixels, which is a useful preprocessing step for the following noise estimation. In this paper, the superpixel segmentation technique is employed to obtain a large number of regions with adaptive shapes and sizes. Liu et al., propose a widely used ERS (entropy rate superpixel) algorithm based on a new objective function for the superpixel segmentation of natural images [26]. The experimental results show that the ERS significantly outperforms the state of the art with respect to the standard metrics. Therefore, we use the ERS algorithm for the superpixel segmentation of noisy HSIs.
The ERS algorithm can be regarded as a clustering problem, which considers the segmentation problem as a maximization optimization on a graph. First, a graph G = (V, E) is established by using the retained PC obtained by MNF, where V is the vertex set representing the pixels, and E is the edge set representing the similarity between adjacent pixels. The goal of the ERS algorithm is to select a subset of edges A such that the resulting graph, G = (V, A), contains exactly K connected subgraphs. Note that the objective function combines the entropy rate and the balancing function and, therefore, favors compact, homogeneous, and balanced superpixels. The segmentation is achieved via optimizing the following objective function with respect to the edge set: where Ψ(A) is the entropy rate of the random walk on G, Φ(A) is a balance term that reduces unbalanced superpixels, λ ≥ 0 is the weight of the balancing term, and N A is the number of connected components in the graph G,. Note that a novel greedy optimization scheme algorithm has been developed in [26] to solve the above optimization problem.

The Mixed Noise Estimation on Superpixels
The proposed algorithm to estimate the mixed noise is developed from [21]. There is an important assumption that the HSI can be divided into several non-overlapping patches with homogeneous regions, in which the local statistics can be calculated correctly. However, in [21], the blocks sized 7 × 7 are taken as the non-overlapping homogeneous regions, which lead to an incorrect estimation when different ground objects are in one block. To overcome the shortcomings of using a block as a homogenous region, the superpixel is taken for a homogenous region in this paper. This approah is more consistent with the relevant assumption. Therefore, HSI is divided into K superpixels Λ 1 , Λ 2 , . . . , Λ K via superpixel segmentation, and each superpixel is viewed as a homogeneous region for one ground object. There are L bands in HSI, and M × N (height × width) pixels in one band. Let g i,l denote intensity of the observed noisy data at position i (1 ≤ i ≤ M × N) in the l-th (1 ≤ l ≤ L) band. By superpixel segmentation, g i,l is assigned into the k-th superpixel, and then g i,l is denoted as g i,l , whose variances are σ 2 u l and σ 2 w l , respectively. σ 2 u l and σ 2 w l are the parameters to be estimated. Therefore, the noise model is as follows: The MRL model [27] is employed to remove the significant spectral correlation in the most adjacent bands. For 1 ≤ i ≤ M × N, the predicted value g i,l of g i,l is computed as follows: l = 1, a l g i,l−1 + b l g i,l+1 + c l , 1 < l < L, a L g i,L−2 + b L g i,L−1 + c L , l = L.
if l = 1, then and if l = L, then To achieve the minimum value of the following objective function, by using the least-squares method, we have: In the k-th superpixel Λ k , the predicted value g (k) i,l of g (k) i,l is computed as follows: where a l , b l and c l are computed by the MRL model.
i,l is computed as follows: Under the assumption that the ground objects are homogeneous in one superpixel, in the k-th i,l is zero, that is to say, where l . That is to say, Remote Sens. 2020, 12, 1324 8 of 27 Since i,l is a constant independent of i. Hence, we have: Since u i,l+1 and constant 1 are independent respectively, by Equations (11), (12), (18)- (22) and the propositions of variance, for 1 ≤ l ≤ L, it holds that, where σ 2 u k and σ 2 w k (1 ≤ l ≤ L) are the unknowns to be estimated. The matrix A is denoted as a partitioned matrix, where where (25) and (26), the linear system (23) is in the following form,

By Equations
where A is a KL × 2L matrix, x is a 2L × 1 vector, and d is a KL × 1 vector.
In practical computation, the theoretical expectation f Hence, the mixed noise parameter estimation turns to the solution of the linear system given by Equation (27). Note that A is a rectangular matrix. The least-squares solution can be determined by the existing algebraic method as follows:x where the elements in vectorx are the estimation of noise variances. The process for estimating the mixed noise is described in Algorithm 2.
Algorithm 2 Superpixel-Based Mixed Noise Estimation for HSI Using Multiple Linear Regression (MLR).
1. Input H: Noisy HSI of the size M × N × L (height × width × band) Input K: the number of superpixels Output:σ u l ,σ w l , 1 ≤ l ≤ L % The standard deviations of the SD and SI noise 2. Perform MNF on H by Algorithm 1, return PC 1 % Retain the first principal 3. Superpixel segmentation on PC 1 , and obtain K superpixels 4. Perform MLR on the two most adjacent bands, and obtain a l , b l and c l , 1 ≤ l ≤ L 5. Calculate A k by Equation (24) and Equation (25) for the k-th superpixel, 1 ≤ k ≤ K 6. Calculate d k by Equation (26) Solving linear equations by the least-squares method 9. returnσ u l ,σ w l , 1 ≤ l ≤ L % Output the standard deviations of the SD and SI noise

The Simulated Data Experiments
In the simulated data experiment, hyperspectral dataset Pavia Centre is taken, which was collected by the ROSIS (Reflective Optics Imaging Spectrometer) in July 2002, in Pavia, a northern city of Italy. The dataset contains 102 bands at the range of 0.43-0.86 µm, with spatial resolution of 1.3m and image size of 1096 × 1096. After removing the atmospheric absorption bands, the remaining 80 bands are considered as noise-free images. In the experiment, a dataset sized 256 × 256 × 80 is used as reference image to add simulated noise.
SNR is defined as SNR = 10 log 10 (P X /P N ), (29) where P X denotes the power of noise-free signal, and P N denotes the power of noise. Let x k i,j denote the noise-free signal in the i-th row, the j-th column, and the k-th band, and let n k i,j denote noise signal in the corresponding position. The SNR of the image in the k-th band is defined as: The relative error in the k-th band is defined as: whereσ k is the estimated standard deviation, i.e., the square root of variance, and σ k is the synthesized standard deviation in the simulated noise. The relative error for a whole HSI dataset is the average value of the relative errors of all the L bands. The absolute error of the k-th band is defined as: The absolute error for a whole HSI dataset is the average value of the absolute errors of all the L bands.
Let α denote the ratio of the SD noise power to the SI noise power, three cases are discussed in the experiment. Case 1: α = 3 : 1, the SD noise is dominant. Case 2: α = 1 : 1, the contribution of SD noise and SI noise is the same. Case 3: α = 1 : 3, the SI noise is dominant.
For abbreviation in the experiments, the proposed algorithm is denoted by the acronym superpixel-based MLRLS (MLR and least-squares method), the algorithm proposed in [21] is denoted by the block-based MLRLS, and the algorithm proposed in [19] is denoted by MLRML (MLR and maximum likelihood method), hereinafter.

Selection of the Number of Superpixels
In this experiment, 30 dB noise is added to each band of HSI. According to the definition of SNR, the noise power in each band is proportional to the power of the noise-free signals in the corresponding band. Since the intensity varies with the wavelength, the noise level varies with the wavelength. The ratio of the SI noise power to the SD noise power is 1:1, that is to say, the power of SI noise and SD noise is competitive, none of which can be ignored when estimated. In Figure 1a,b, it is difficult for human eyes to distinguish the slight difference between the noise-free data and the noisy data of 30 dB noise. However, in many HSI post-processing methods that are sensitive to noise, 30 dB is high enough to reduce the performance significantly [28,29]. Therefore, we first choose 30 dB noise with α = 1 : 1 for the experiment. The first step of the algorithm is MNF. Figure 2 gives the top eight PCs obtained by MNF with the dataset Pavia Centre. The first PC has the strongest signal power, while the following PCs contain less signal power. Since the noise in HSI has been processed by the noise-whitening matrix, the first PC has the highest SNR compared to the other PCs. This advantage ensures the performance of the subsequent superpixel segmentation. The first step of the algorithm is MNF. Figure 2 gives the top eight PCs obtained by MNF with the dataset Pavia Centre. The first PC has the strongest signal power, while the following PCs contain less signal power. Since the noise in HSI has been processed by the noise-whitening matrix, the first PC has the highest SNR compared to the other PCs. This advantage ensures the performance of the subsequent superpixel segmentation.
The first step of the algorithm is MNF. Figure 2 gives the top eight PCs obtained by MNF with the dataset Pavia Centre. The first PC has the strongest signal power, while the following PCs contain less signal power. Since the noise in HSI has been processed by the noise-whitening matrix, the first PC has the highest SNR compared to the other PCs. This advantage ensures the performance of the subsequent superpixel segmentation. The superpixel segmentation results achieved by the ERS on PC1 are given in Figure 3, where the superpixel number K is 80, 600 and 1200. Note that the superpixel number has a direct impact on the final noise estimation result. Therefore, it is necessary to determine a reasonable number K. The superpixel segmentation results achieved by the ERS on PC1 are given in Figure 3, where the superpixel number K is 80, 600 and 1200. Note that the superpixel number has a direct impact on the final noise estimation result. Therefore, it is necessary to determine a reasonable number K.
Remote Sens. 2020, 12, x FOR PEER REVIEW 12 of 28 In the block-based MLRLS [21], the size of a block is 7 × 7. In this paper, more sizes N×N (height  width) each of them is viewed, where N ranges from 2 to 32. Meanwhile, in order to find a proper number K of the superpixels in the superpixel-based MLRLS, different numbers K are tested from 100 to 5000 with a step size of 100. The relative errors of the SD noise estimated by the block-based MLRLS and the superpixel-based MLRLS are given in Figure 4. The minimal relative error of the block-based MLRLS is 1.90%, with the length side of a block being 2 N  . When N is in a range from 2 to 7, the relative error varies stably in a narrow interval (1.90, 2.38). The blue curve in Figure 5 (b) shows the relative error of the SD noise estimated by the superpixel-based MLRLS, and the red line is the minimal relative error of the block-based MLRLS. When K is in a wide interval [900, 2600], the relative error fluctuates slightly along the red line, and is limited in a narrow interval (1.81, 2.07). The minimal relative error for the superpixel-based MLRLS is 1.81%, with 1800 K  . Figure 5 shows the relative errors of the SI noise estimated by the block-based MLRLS and the superpixel-based MLRLS. In Figure 5 (a), the minimal relative error of the block-based MLRLS is 2.73%, with 3 N  . When N is in the range from 2 to 7, the relative error varies stably and is limited in a narrow interval (2.73, 2.98). Figure 5 (b) shows the relative errors of the SI noise estimated by the superpixel-based MLRLS. The minimal relative error of the superpixel-based MLRLS is 2.22%, with 2300 K  . When K is in a wide interval [1000, 3600], all the values of the relative errors are below the red line. This means that the number of superpixels K taken in a large range can ensure In the block-based MLRLS [21], the size of a block is 7 × 7. In this paper, more sizes N × N (height × width) each of them is viewed, where N ranges from 2 to 32. Meanwhile, in order to find a proper number K of the superpixels in the superpixel-based MLRLS, different numbers K are tested from 100 to 5000 with a step size of 100. The relative errors of the SD noise estimated by the block-based MLRLS and the superpixel-based MLRLS are given in Figure 4. The minimal relative error of the block-based MLRLS is 1.90%, with the length side of a block being N = 2. When N is in a range from 2 to 7, the relative error varies stably in a narrow interval (1.90, 2.38). The blue curve in Figure 5b shows the relative error of the SD noise estimated by the superpixel-based MLRLS, and the red line is the minimal relative error of the block-based MLRLS. When K is in a wide interval [900, 2600], the relative error fluctuates slightly along the red line, and is limited in a narrow interval (1.81, 2.07). The minimal relative error for the superpixel-based MLRLS is 1.81%, with K = 1800.
Remote Sens. 2020, 12, x FOR PEER REVIEW 13 of 28 the superpixel-based MLRLS can still perform better than the block-based MLRLS with a wide selection of K when the ground truth noise is unknown.

Experiments for Different Cases
In this section, the algorithm is carried out on HSIs with different noise levels and different noise components. 25 dB, 30 dB and 35 dB noises with different  are added on the hyperspectral dataset Pavia Centre. Table 1 and Figure 7 show the relative errors of the superpixel-based MLRLS, the block-based MLRLS [21], and MLRML [19].
When SNR = 25 dB, the proposed algorithm does not improve much in terms of the relative errors. When

 
, the accuracy is slightly decreased. At low SNR, the high level noise affects the accuracy of superpixel segmentation and subsequently degrades the performance of the proposed algorithm. Therefore, when the noise level is high, the

Experiments for Different Cases
In this section, the algorithm is carried out on HSIs with different noise levels and different noise components. 25 dB, 30 dB and 35 dB noises with different  are added on the hyperspectral dataset Pavia Centre. Table 1 and Figure 7 show the relative errors of the superpixel-based MLRLS, the block-based MLRLS [21], and MLRML [19].
When SNR = 25 dB, the proposed algorithm does not improve much in terms of the relative errors. When

 
, the accuracy is slightly decreased. At low SNR, the high level noise affects the accuracy of superpixel segmentation and subsequently degrades the performance of the proposed algorithm. Therefore, when the noise level is high, the  Figure 5 shows the relative errors of the SI noise estimated by the block-based MLRLS and the superpixel-based MLRLS. In Figure 5a, the minimal relative error of the block-based MLRLS is 2.73%, with N = 3. When N is in the range from 2 to 7, the relative error varies stably and is limited in a narrow interval (2.73, 2.98). Figure 5b shows the relative errors of the SI noise estimated by the superpixel-based MLRLS. The minimal relative error of the superpixel-based MLRLS is 2.22%, with K = 2300. When K is in a wide interval [1000, 3600], all the values of the relative errors are below the red line. This means that the number of superpixels K taken in a large range can ensure that the relative error is smaller than that of the block-based MLRLS. Meanwhile, when K is in the interval (1300, 3300), the relative error fluctuates slightly in a narrow interval (2.22, 2.53).
Since the relative errors of both the SD noise and the SI noise do not always achieve the minimal value at the same point N or K, the average value of the relative errors of the SD noise and the SI noise is also tested for an overall optimization. Figure 6 shows the relative errors for the overall optimization taking both the SD and SI noise into account. The minimal relative error of the block-based MLRLS is 2.36%, with N = 2. When N is in a range of 2 to 7, the relative error varies stably in a narrow interval (2.36, 2.64). As N becomse large, the general trend of the relative error becomes large. The reason for this phenomenon is that when the size of a block becomes large, it will contain more types of ground objects in an identical block, which will lead to inhomogeneous regions in one block and violate to assumption of homogeneity. In Figure 6b, the minimal relative error of the superpixel-based MLRLS is 2.06%, with K = 1800. When K is in a wide interval [1000, 3300], all the relative errors of the superpixel-based MLRLS are smaller than the minimal relative error of the block-based MLRLS. Meanwhile, when K is in the range of (1300, 3300), the relative error fluctuates slightly in a narrow interval (2.06, 2.27).
By the comparison, the minimal relative error of the noise with K = 1800 estimated by the superpixel-based MLRLS is smaller than that of the block-based MLRLS with N = 2. However, the selection of the optimal length side of a block N in the block-based MLRLS and the optimal number of superpixels K in the superpixel-based MLRLS is dependent of the HSI dataset, and involves many factors such as the noise level, the edges, etc. The above experimental analysis indicates that the superpixel-based MLRLS can still perform better than the block-based MLRLS with a wide selection of K when the ground truth noise is unknown.

Experiments for Different Cases
In this section, the algorithm is carried out on HSIs with different noise levels and different noise components. 25 dB, 30 dB and 35 dB noises with different  are added on the hyperspectral dataset Pavia Centre. Table 1 and Figure 7 show the relative errors of the superpixel-based MLRLS, the block-based MLRLS [21], and MLRML [19].
When SNR = 25 dB, the proposed algorithm does not improve much in terms of the relative errors. When

 
, the accuracy is slightly decreased. At low SNR, the high level noise affects the accuracy of superpixel segmentation and subsequently degrades the performance of the proposed algorithm. Therefore, when the noise level is high, the

Experiments for Different Cases
In this section, the algorithm is carried out on HSIs with different noise levels and different noise components. 25 dB, 30 dB and 35 dB noises with different α are added on the hyperspectral dataset Pavia Centre. Table 1 and Figure 7 show the relative errors of the superpixel-based MLRLS, the block-based MLRLS [21], and MLRML [19]. When SNR = 25 dB, the proposed algorithm does not improve much in terms of the relative errors. When α = 1 : 3, the overall relative error is only reduced by 0.13%, and when α = 3 : 1, the overall relative error is only reduced by 0.07%. When α = 1 : 1, the accuracy is slightly decreased. At low SNR, the high level noise affects the accuracy of superpixel segmentation and subsequently degrades the performance of the proposed algorithm. Therefore, when the noise level is high, the superiority of the superpixel-based MLRLS is not obvious over the block-based MLRLS. When SNR = 30 dB, the overall relative error of the superpixel-based MLRLS is 0.34% less than that of the block-based MLRLS, while when SNR = 35 dB, the overall relative error is 1.88% less than that of the block-based MLRLS, and the estimation accuracy is improved by 1.01%-2.89%. As the SNR becomes larger, the superpixel segmentation performs better. Therefore, Table 1 shows that the superpixel-based MLRLS outperforms the block-based MLRLS, especially for high-SNR images. In addition, the absolute errors are listed in Table 2 and Figure 8. These figures indicate that the absolute errors decrease with a higher SNR. Meanwhile, the MNF approach and using superpixel segmentation to find homogenous regions play an important role in the mixed noise estimation, especially for HSIs of high SNRs.   Figure 7. The optimal relative errors of the SD noise, the SI noise and the overall optimization.     A bar plot of Table 1 is given in Figure 7. As shown in Table 1 and Figure 7, the relative error becomes large as the SNR increases. However, this does not mean that a high SNR may lead to poor performance. In fact, a high SNR can ensure the performance of superpixel segmentation, which improves accuracy in terms of the absolute error. Since the standard deviations of the SD noise and the SI noise vary greatly in magnitude, to compare them at the same level, the relative error ε r = |σ − σ|/σ × 100% is introduced for normalization. In ε r = |σ − σ|/σ × 100%, if σ becomes smaller, ε r becomes larger. That is to say, σ is smaller with a higher SNR, accordingly, the relative error ε r becomes larger. To eliminate this misunderstanding that a high SNR leads to poor performance, the absolute errors |σ − σ| are also calculated in the experiments. Table 2 and Figure 8 show the absolute errors at different noise levels and in different cases. When the SNR is higher, the absolute errors are smaller. Figure 9 shows the standard deviation of the estimated SD and SI noise in each band compared with that of the synthetic noise, when α = 1 : 1, and SNR = 25 dB, 30 dB and 35 dB, respectively. When SNR = 35 dB, the curve of the standard deviation of the superpixel-based MLRLS is obviously closer to that of the synthetic noise than that of the block-based MLRLS. Next, we will discuss the selection of the side length N of a block in the block-based MLRLS and the number of the superpixels K in the superpixel-based MLRLS at different noise levels. Table 3 shows the side length N when the block-based MLRLS gets the optimal relative errors at different noise levels. It indicates that N fluctuates in a range from 2 to 8, and the average value of N is 4.04. Therefore, 4 N  is recommended for the block-based MLRLS performed on Furthermore, the proposed algorithm is compared with the MLRML method [19]. The MLRML only takes high spectral correlation into account, but ignores spatial correlation. When SNR = 25 dB, the noise level is high, the performance of the noise estimation achieved by the MLRML is acceptable. When the noise level is low, the MLRML fails to distinguish subtle difference between the noise and the details of the HSIs, which leads to an inaccurate noise estimation.
The experimental results with different noise levels and different noise components indicate that with the development of imaging quality using modern spectrometers, the proposed algorithm the superpixel-based MLRLS is more suitable for the noise estimation of modern HSI datasets.
Next, we will discuss the selection of the side length N of a block in the block-based MLRLS and the number of the superpixels K in the superpixel-based MLRLS at different noise levels. Table 3 shows the side length N when the block-based MLRLS gets the optimal relative errors at different noise levels. It indicates that N fluctuates in a range from 2 to 8, and the average value of N is 4.04. Therefore, N = 4 is recommended for the block-based MLRLS performed on the real-life data when the ground truth noise is unknown.  Table 4 and Figure 10 show the number of superpixels K when the superpixel-based MLRLS achieves the optimal relative errors at different noise levels. As can be seen, the optimal K is from 800 to 5700, and the average value of K is 2614.8, i.e., on average, 25 pixels in a superpixel. Meanwhile, when SNR = 25 dB, the average value of K is 1144.4, i.e., on average, 57 pixels in one superpixel; when SNR = 30 dB, the average value of K is 1988.9, i.e., on average, 33 pixels in one superpixel; when SNR = 35 dB, the average value of K is 4711.1, i.e., on average, 14 pixels in one superpixel. Thus, the selection of K is dependent on the noise level of the HSIs. Therefore, it is recommended that on average, 25 pixels in one superpixel is selected in the superpixel-based MLRLS. If the noise level is high, one can choose more pixels in a superpixel, i.e., a smaller K, and vice versa.  Figure 10. The number of superpixels K of the optimal relative error for the superpixel-based MLRLS. Table 5 and Figure 11 show the computational time of the two algorithms. The algorithms are compiled and run in Matlab R2015a. The computer processor is Intel(R) Core(TM) i7-10710U CPU Figure 10. The number of superpixels K of the optimal relative error for the superpixel-based MLRLS. Table 5 and Figure 11 show the computational time of the two algorithms. The algorithms are compiled and run in Matlab R2015a. The computer processor is Intel(R) Core(TM) i7-10710U CPU @1.10GHz and 1.61GHz, and the memory capacity is 16.0 GB. The average computational time of the superpixel-based MLRLS is 21.03 s, while the average computational time of the block-based MLRLS is 58.57 s. In the superpixel-based MLRLS, the computational time of MNF is 1.33 s, on average. The superpixel segmentation used in the proposed algorithm is based on a greedy optimization scheme, and is only performed on the first component. Therefore, the computational cost is not expensive. The superpixel segmentation takes only 0.32 s, on average. Most of the computational time is spent on calculating the matrix elements in Equations (24)- (26) and the least-squares method in Equation (28). The computational complexity is dependent on the number of superpixels or the number of blocks. Since the average number of superpixels K is 2614.8, while the average side length N of a block is 4.04, i.e., 16,222 blocks in the tested HSI data, the average computational time of the superpixel-based MLRLS is much lower than the block-based MLRLS.

Experiments on Low-Resolution Datasets
The proposed algorithm is based on the assumption that the HSIs have large spatial and spectral resolution. In this section, we discuss how lower spatial and spectral resolution affect the estimation performance of the proposed algorithm.

A. Lower Spatial Resolution
In this experiment, spatial subsampling is performed on the HSI dataset Pavia Centre to obtain a lower spatial resolution. The dataset is subsampled by 1:2 in spatial rows and 1:2 in columns, respectively. Thus, the subsampled dataset has a spatial resolution of 1:4 of the original image. The superpixel-based MLRLS and the block-based MLRLS are performed on the subsampled dataset. The estimated relative error is shown in Table 6 and Figure 12. This indicates that after spatial subsampling, the estimation errors of the two algorithms decrease compared to higher spatial resolution images. The estimated relative error of the superpixel-based MLRLS decreases by 1.26%-5.04%, with an average decrease of 2.33%. The estimated relative error of the block-based MLRLS decreases by 0.77%-4.23%, with an average decrease of 2.14%.  Figure 12. The optimal relative errors of the SD noise, the SI noise and the overall optimization for 1:4 spatial subsampling HSI.
relative error relative error relative error relative error relative error relative error relative error relative error relative error Figure 12. The optimal relative errors of the SD noise, the SI noise and the overall optimization for 1:4 spatial subsampling HSI.
Although the lower spatial resolution makes the estimation accuracy of the superpixel-based MLRLS decreased slightly more than that of the block-based MLRLS, the estimation accuracy of the superpixel-based MLRLS in most cases is still higher than that of the block-based MLRLS.

B. Lower Spectral Resolution
Spectral subsampling is also performed on the HSI Pavia Centre to obtain a lower spectral resolution. The dataset is subsampled by 1:2 in the spectrum. The superpixel-based MLRLS and the block-based MLRLS are performed on the subsampled dataset. The estimated relative error is shown in Table 7 and Figure 13. This indicates that after spectral subsampling, the estimation errors of the two algorithms decrease compared to higher spectral resolution images. The estimated relative error of the superpixel-based MLRLS decreases by 0.14%-23.84% (except for that of the SD noise when SNR = 25dB, α = 1 : 3), with an average decrease of 4.13%. The estimated relative error of the block-based MLRLS decreases by 0.43%-33.27%, with an average decrease of 8.14%. Table 7. The optimal relative errors of the SD noise, the SI noise and the overall optimization for 1:2 spectral subsampling HSI. Compared to the subsampling at a 1:4 spatial resolution, the subsampling at a 1:2 spectral resolution has a serious impact on the accuracy of noise estimation. This means that a high spectral correlation is necessary to ensure the accuracy of noise estimation, and a lower spectral resolution may lead to an unreliable result. The main reason for this is that the two algorithms are based on MLR, under the assumption of a high spectral correlation.  Figure 13. The optimal relative errors of the SD noise, the SI noise and the overall optimization for 1:2 spectral subsampling HSI.
relative error relative error relative error relative error relative error relative error relative error relative error relative error Figure 13. The optimal relative errors of the SD noise, the SI noise and the overall optimization for 1:2 spectral subsampling HSI.
Although both the algorithms are sensitive to spectral resolution, the performance of the superpixel-based MLRLS is still better than that of the block-based MLRLS by 0.33%-13.39%.

Experiments of the Textured Datasets
In the above experiments, the superpixel-based MLRLS achieves better performance on the HSI data having well-distinguishable homogeneous regions. In this experiment, we will test the estimation on a textured HSI dataset and discuss what happens if the scene is highly textured when there are few recognizable homogeneous regions. The dataset was acquired over Yellow Stone by the AVIRIS (airborne visible/infrared imaging spectrometer), from JPL (Jet Propulsion Laboratory), NASA (National Aeronautics and Space Administration). One scene is 614 × 512 × 224 in size. A dataset sized 256 × 256 × 30 is taken from the 61st band to the 90th band as a reference dataset. Figure 14 shows its snapshot and superpixel segmentation.
Remote Sens. 2020, 12, x FOR PEER REVIEW 21 of 28 Compared to the subsampling at a 1:4 spatial resolution, the subsampling at a 1:2 spectral resolution has a serious impact on the accuracy of noise estimation. This means that a high spectral correlation is necessary to ensure the accuracy of noise estimation, and a lower spectral resolution may lead to an unreliable result. The main reason for this is that the two algorithms are based on MLR, under the assumption of a high spectral correlation.
Although both the algorithms are sensitive to spectral resolution, the performance of the superpixel-based MLRLS is still better than that of the block-based MLRLS by 0.33%-13.39%.

Experiments of the Textured Datasets
In the above experiments, the superpixel-based MLRLS achieves better performance on the HSI data having well-distinguishable homogeneous regions. In this experiment, we will test the estimation on a textured HSI dataset and discuss what happens if the scene is highly textured when there are few recognizable homogeneous regions. The dataset was acquired over Yellow Stone by the AVIRIS (airborne visible/infrared imaging spectrometer), from JPL (Jet Propulsion Laboratory), NASA (National Aeronautics and Space Administration). One scene is 614 × 512 × 224 in size. A dataset sized 256 × 256 × 30 is taken from the 61st band to the 90th band as a reference dataset. Figure  14 shows its snapshot and superpixel segmentation. The estimated relative errors are shown in Table 8 and Figure 15, and the estimated absolute errors are shown in Table 9 and Figure 16. The experimental results show that the proposed algorithm performs not as good as on the regular HSI datasets. This is because superpixel segmentation cannot accurately detect the boundary or contour of a highly textured HSI as a regular image. Therefore, the HSI might not have enough homogeneous regions to guarantee a reliable estimation of noise variances. Thus, when estimating the noise level of a highly textured HSI dataset, the block-based MLRLS may be a proper substitute for the superpixel-based MLRLS. Table 8. The optimal relative errors of the SD noise, the SI noise and the overall optimization for the textured HSI. The estimated relative errors are shown in Table 8 and Figure 15, and the estimated absolute errors are shown in Table 9 and Figure 16. The experimental results show that the proposed algorithm performs not as good as on the regular HSI datasets. This is because superpixel segmentation cannot accurately detect the boundary or contour of a highly textured HSI as a regular image. Therefore, the HSI might not have enough homogeneous regions to guarantee a reliable estimation of noise variances. Thus, when estimating the noise level of a highly textured HSI dataset, the block-based MLRLS may be a proper substitute for the superpixel-based MLRLS.  Figure 15. The optimal relative errors of the SD noise, the SI noise and the overall optimization for the textured HSI.
relative error relative error relative error relative error relative error relative error relative error relative error Relative error Figure 15. The optimal relative errors of the SD noise, the SI noise and the overall optimization for the textured HSI.

The Real-Life Data Experiment
The Dataset Low Altitude is used for a real-life data experiment. The dataset was acquired by AVIRIS from JPL, NASA. Experiment 1. Two sub-datasets of different ground objects are cut from the scene. The size of the selected datasets is 256 × 256 × 224. Figure 16 shows the 60th band of the two datasets. One dataset covers complex ground objects including fields, forests, river and roads. In the second dataset, forest covers most of the land, except for some field and river areas. Although the land-cover of the ground objects in two datasets are quite different, they are cut from the same scene and acquired in the same condition, so their noise levels should follow the same pattern, and the estimation results should not vary significantly. Figure 17 shows the 60th band of the sub-datasets and their superpixel segmentation.   Figure 16. The optimal absolute errors of the SD noise, the SI noise and the concerned optimization for the textured HSI.

The Real-Life Data Experiment
The Dataset Low Altitude is used for a real-life data experiment. The dataset was acquired by AVIRIS from JPL, NASA. Experiment 1: Two sub-datasets of different ground objects are cut from the scene. The size of the selected datasets is 256 × 256 × 224. Figure 16 shows the 60th band of the two datasets. One dataset covers complex ground objects including fields, forests, river and roads. In the second dataset, forest covers most of the land, except for some field and river areas. Although the land-cover of the ground objects in two datasets are quite different, they are cut from the same scene and acquired in the same condition, so their noise levels should follow the same pattern, and the estimation results should not vary significantly. Figure 17 shows the 60th band of the sub-datasets and their superpixel segmentation. The number of superpixels is 3000 in this experiment. Figure 18 gives the estimated standard deviations vs. the bands of the two datasets. It is observed that the noise levels in bands 1-6, 31-42, 107-112, 154-168 and 221-224 are extremely high, even several times higher than those of the other bands. Meanwhile, these bands are contaminated by water vapor absorption, and their spectral correlation is weak. This means that the algorithm is unreliable in the bands severely contaminated by water vapor absorption. This is consistent with the discussion in Section 5.1.3 that a lower spectral correlation will lead to an unreliable estimation. The number of superpixels is 3000 in this experiment. Figure 18 gives the estimated standard deviations vs. the bands of the two datasets. It is observed that the noise levels in bands 1-6, 31-42, 107-112, 154-168 and 221-224 are extremely high, even several times higher than those of the other bands. Meanwhile, these bands are contaminated by water vapor absorption, and their spectral correlation is weak. This means that the algorithm is unreliable in the bands severely contaminated by water vapor absorption. This is consistent with the discussion in Section 5.1.3 that a lower spectral correlation will lead to an unreliable estimation.

Experiment 2:
In this experiment, the proposed algorithm the superpixel-based MLRLS is applied to the second sub-dataset of Low Altitude. Since there is no ground truth to evaluate the results of noise estimation for a real-life dataset, the denoising method in the [30,31] is chosen for comparison, indirectly verifying the reliability of the proposed algorithm. This denoising method is based on wavelet transform and MRL, accordingly, denoted by the acronym Wavelet+MRL. In fact, a denoising method cannot separately estimate the variances of the SD noise and SI noise. However, it can estimate the noise power by the difference between the noisy dataset and the denoised dataset. On the one hand, by Equation (5) On the other hand, Figure 18 demonstrates that although the ground object components of the two sub-datasets are quite different, their standard deviation curves are similar. In this paper, the Pearson correlation coefficient is used as the similarity metrics for the standard deviation vs. the band. The Pearson correlation coefficient measures how highly correlated two vectors are in a range of [−1, 1]. A correlation coefficient of 1 means that the two vectors are perfectly positively linearly related, while a correlation of 0 means that the two vectors do not have any linear relation. The coefficient of the estimated standard deviation vs. band of the SD noise is 0.9218, and that of the SI noise is 0.8449. Figure 19 gives a histogram of the differences between the estimated standard deviations of the two datasets. In most of the bands, the differences between the estimated standard deviations of the two datasets are very close to zero. Occasionally, a few values corresponding to unreliable estimation in the water vapor absorption bands deviate far from zero. Eliminating these bands, the coefficient of the estimated standard deviation vs. the band of the SD noise is 0.9828, and that of the SI noise is 0.9402. The estimation on the two datasets covered by different ground objects from one scene demonstrates two significantly similar standard deviation curves in Figure 18. This indicates that the estimated noise level is determined by the sensor, but does not rely on the ground objects. Therefore, the proposed algorithm embodies the intrinsic noise of the imaging instruments, and performs stably and reliably.

Experiment 2.
In this experiment, the proposed algorithm the superpixel-based MLRLS is applied to the second sub-dataset of Low Altitude. Since there is no ground truth to evaluate the results of noise estimation for a real-life dataset, the denoising method in the [30,31] is chosen for comparison, indirectly verifying the reliability of the proposed algorithm. This denoising method is based on wavelet transform and MLR, accordingly, denoted by the acronym Wavelet+MLR. In fact, a denoising method cannot separately estimate the variances of the SD noise and SI noise. However, it can estimate the noise power by the difference between the noisy dataset and the denoised dataset. On the one hand, by Equation (5), the total standard deviation σ T k can be calculated by the noise power. On the other hand, by Equation (5), the total standard deviation σ T k can be calculated by the standard deviations of the SD noise and the SI noise estimated by the superpixel-based MLRLS. If the total standard deviations σ T k estimated by the two methods match, the effectiveness and reliability of the superpixel-based MLRLS can be verified. Figure 20 shows the curves of the total standard deviation σ T k vs. the band estimated by the two methods. Figure 21 shows a histogram of the differences of the total standard deviations σ T k estimated by the two methods. It is observed that in most of the bands, the differences of the total standard deviation σ T k between the two methods are very close to zero, except for a few water vapor absorption bands. The Pearson correlation coefficient of the estimated total standard deviation σ T k is 0.9885, which is very close to 1. This indicates that the estimated total standard deviation σ T k from the two different methods match perfectly. Therefore, the proposed algorithm is effective and reliable.

Conclusions
This paper proposes a mixed noise variance estimation algorithm for HSIs developed from the the block-based MLRLS method. The proposed algorithm decouples the spectral correlation MLR in homogenous regions detected by MNF-based superpixel segmentation. Using this method, local statistics are calculated according to the residual image in each homogenous region to build a linear system, and the noise variances are obtained by solving the linear equations.
The main difference between the proposed method, the superpixel-based MLRLS and its predecessor, the block-based MLRLS is that the former calculates the local statistics for estimation on superpixels viewed as homogenous regions, and the latter calculates on blocks of a fixed size. In

Conclusions
This paper proposes a mixed noise variance estimation algorithm for HSIs developed from the the block-based MLRLS method. The proposed algorithm decouples the spectral correlation MLR in homogenous regions detected by MNF-based superpixel segmentation. Using this method, local statistics are calculated according to the residual image in each homogenous region to build a linear system, and the noise variances are obtained by solving the linear equations.
The main difference between the proposed method, the superpixel-based MLRLS and its predecessor, the block-based MLRLS is that the former calculates the local statistics for estimation on superpixels viewed as homogenous regions, and the latter calculates on blocks of a fixed size. In the experiments, the difference between segmenting an HSI image into fixed-size blocks or total standard deviation number of bands Figure 21. A histogram of the differences between the estimated standard deviations of the two datasets.

Conclusions
This paper proposes a mixed noise variance estimation algorithm for HSIs developed from the the block-based MLRLS method. The proposed algorithm decouples the spectral correlation MLR in homogenous regions detected by MNF-based superpixel segmentation. Using this method, local statistics are calculated according to the residual image in each homogenous region to build a linear system, and the noise variances are obtained by solving the linear equations.
The main difference between the proposed method, the superpixel-based MLRLS and its predecessor, the block-based MLRLS is that the former calculates the local statistics for estimation on superpixels viewed as homogenous regions, and the latter calculates on blocks of a fixed size. In the experiments, the difference between segmenting an HSI image into fixed-size blocks or superpixels lies in the performance of the relative error versus the side length of a block or the number of superpixels. In the block-based MLRLS, the relative error of estimation roughly increases with the number of blocks for SD noise as well as for SI noise. However, for the superpixel-based MLRLS there is a local minimum in the graph of relative error against the number of superpixels, both for SD and SI noise. The local minimum value for the superpixel-based MLRLS is less than the minimum error value for the block-based MLRLS. Therefore, the proposed method demonstrates a better performance than its predecessor.
The computational expense is also related to the side length of a block or the number of superpixels. Since the superpixel segmentation used in the proposed method is a greedy optimization and is performed only on the first component, its computational expense is acceptable. The majority of the computational time is spent on calculating the local statistics and solving the overdetermined equation with the least square methods. To achieve an optimal estimation, the number of blocks in the block-based MLRLS is far greater than the number of superpixels in the superpixel-based MLRLS, therefore, therefore, the superpixel-based MLRLS costs less time and is superior to the block-based MLRLS in terms of the computational expense.
For further discussion, HSIs at a lower spatial and spectral resolution are tested. After spatial or spectral subsampling, the estimation errors of the two algorithms decrease compared to higher-resolution images, and the lower spectral resolution has a greater impact on the estimation. However, at lower resolution, the superpixel-based MLRLS still performs better than the block-based MLRLS.
In the future, the proposed algorithm should be further improved to achieve better performance in noise estimation on HSIs at a low spectral resolution and with rich textures.