Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Particle–Hole Transformation in Strongly-Doped Iron-Based Superconductors
Next Article in Special Issue
Anomaly Detection Based on Mining Six Local Data Features and BP Neural Network
Previous Article in Journal
Numerical Simulation of PDEs by Local Meshless Differential Quadrature Collocation Method
Previous Article in Special Issue
Parametric Fault Diagnosis of Analog Circuits Based on a Semi-Supervised Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Edge Preserving Weighted Mean Filter for Removing Random-Valued Impulse Noise

1
Department of Electrical Engineering, University of Engineering and Technology, P.O. Box. 814, Peshawar 25120, Pakistan
2
School of Intelligent Mechatronics Engineering, Sejong University, Seoul 05006, Korea
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(3), 395; https://doi.org/10.3390/sym11030395
Submission received: 23 January 2019 / Revised: 13 March 2019 / Accepted: 14 March 2019 / Published: 18 March 2019
(This article belongs to the Special Issue Symmetry in Engineering Sciences)

Abstract

:
This paper proposes an adaptive noise detector and a new weighted mean filter to remove random-valued impulse noise from the images. Unlike other noise detectors, the proposed detector computes a new and adaptive threshold for each pixel. The detection accuracy is further improved by employing edge identification stage to ensure that the edge pixels are not incorrectly detected as noisy pixels. Thus, preserving the edges avoids faulty detection of noise. In the filtering stage, a new weighted mean filter is designed to filter only those pixels which are identified as noisy in the first stage. Different from other filters, the proposed filter divides the pixels into clusters of noisy and clean pixels and thus takes into only clean pixels to find the replacement of the noisy pixel. Simulation results show that the proposed method outperforms state-of-the-art noise detection methods in suppressing random valued impulse noise.

1. Introduction

For image analysis, de-noising is an important pre-processing step. Digital images are oftentimes corrupted by impulse noise during acquisition, transmission and impaired camera sensors [1], which degrades image features such as edges, sharpness, depth, etc. These degradations can severely hamper the performance of some post-processing steps such as image segmentation, edge detection, feature extraction, target detection, recognition and classification. Therefore, the removal of impulse noise from these images is one of the most fundamental problems in the field of digital image processing. It refers to the removal of noise from corrupted images and the preservation of useful information such as edges and discontinuities. Unlike Gaussian noise, Impulse noise does not corrupt every pixel in the image. It corrupts only certain number of pixels based on noise density. Impulse noise is of two types: Salt & Pepper Noise (SPN) and Random Valued Impulse Noise (RVIN) [2]. SPN corrupts image with two fixed extreme Values while RVIN takes any arbitrary value in range [ n m i n , n m a x ] . For 8-bit images, n m i n = 0 and n m a x = 255. Therefore, detection of RVIN is much more complex and difficult as compare to SPN [3]. For simplicity, let x i , j and c i , j be the intensities of pixels of noisy and original images, respectively. Then, the noise model for RVIN [4] can be described as:
x i , j = c i , j , for p = 1 p 0 , n i , j , for p = p 0 ,
where n i , j shows gray level of uniformly distributed noise in range [ n m i n , n m a x ] while p 0 is the probability of noise.
To suppress impulse noise, different filters have been used. However, most of these filters have a disadvantage of losing some important information such as edges [5]. These are simple low pass filters where, for every pixel, some operation is performed on all the neighboring pixels in the predefined window and the result replaces the central pixel. The simplest filter to remove impulse noise is median filter (MF) [6]. This filter replaces every pixel in the image by the median value of the given window, irrespective of whether the pixel is noisy or non-noisy. Though MF is computationally simple and less expensive, but as the clean pixels are also replaced by median value of the selected window, this degrades the quality of image, and also produces edge jitters and streaking. Similarly, some patches are also created at higher noise density [7].
To make a trade-off between edges preservation and complexity, multiple variants of MF [8] have been proposed. Some proposed variants are Tri-State Median filter (TSMF) [9], Recursive Weighted Median Filter (RWMF) [10], the Multi-state Median Filter [11] the Central Weighted Median Filter (CWMF) [12], the Rank-Order Mean Filter and the Stack Filter [13]. However, these filters still degrade the quality of images as they replace all pixels in the image without considering whether the test pixel is noisy or not [14]. To overcome this drawback, several filtering schemes have been developed which are integrated with noise detectors. In these techniques, image de-noising is performed in two stages. In the first stage, the noisy pixels are detected by approximation from neighborhood pixels. In the second stage, only those pixels which are detected as noisy are cleaned while remaining pixels are left unchanged. The performance of detection stage heavily depends on the estimation of proper threshold value. A prior threshold is used by most of the techniques, where the predefined threshold is used throughout the image. Some notable examples of the filtering techniques which use a fixed threshold for noise detection are: adaptive center-weighted median (ACWM) filter [15], the optimal direction median filter [2], the two-pass switching rank-ordered arithmetic mean (TSRAM) [16] filter, Progressive Switching Median Filter (PSMF) [17], Luo-Iterative Median Filter (Luo) [18], Directional Weighted Median (DWM) filter [19], the contrast Enhancement-based Filter (CEF) [20] and robust outlyingness ratio non-local mean (ROR-NLM) [21], etc. The ACWM supresses the RVIN by using the difference between the output of CWM and the current pixel. In CEF, exponential function is used to enlarge the differences, which are then summed to detect the noisy pixels. The DWM only considers pixels along four certain directions to detect noise. Once the noise is detected, the noisy pixel is replaced by weighted median value in the optimal direction. In ROR-NLM [21], the weights of non-local means filter (NLM) [22] were used as noise detectors on initial denoised image. These filters usually perform well at lower noise density. However, at higher noise density, these filters have limited performance in terms of image details and edges preservation. This is because the statistics of the image are not uniform for the whole image, so the threshold value suitable for specific region may not be adapted well to other region in same image. In addition, the threshold value adjusted for a particular image may not work well for other image.
In order to reduce the shortcomings of the fixed threshold scheme, an Adaptive Switching Median Filter (ASWM) [23] was proposed. This filter automatically defines a new threshold for every pixel based on the local standard deviation, and does not require a priori knowledge for the threshold selection. Though this filter provides better detection results, when the noise is not too much high. However, due to random noise, some intensity values of noisy pixels may be very different from other neighbor pixels inside a current sliding window—in which case, the standard deviation of current window is very large, which can result a very large threshold. Thus, the higher standard deviation can cause inaccurate calculation of the local threshold. Secondly, the detected noisy pixel is replaced with simple median value of the whole current window. As the current sliding window also contains noisy pixels, the median value of such noisy window is thus not always accurate. This leads to incorrect replacement of a noisy pixel. Thus, the image details are not preserved. In order to overcome the shortcomings of this technique, we propose an adaptive two staged technique, where, for each pixel processing, a more accurate threshold is calculated. Hence, the main contributions of our proposed technique are:
  • We propose an adaptive noise detector and a new weighted mean filter to remove random valued impulse noise from the images.
  • In contrast to other state-of-the-art noise detectors, the proposed detector computes a novel adaptive threshold for each pixel. The accuracy of the detection scheme is further enhanced by employing additional steps to avoid faulty detection of edge pixels, which are noisy pixels. This step indeed ensures the preservation of the edges avoiding faulty detection of noise.
  • In the filtering stage, a novel weighted mean filter is designed to filter only those pixels which are identified as noisy in the detection stage. In contrast to the existing filtering techniques, the proposed filter first divides the pixels into clusters of clean and noisy pixels. Then, to design the filter, only clean pixels are utilized to find the replacement of the noisy pixel.
Extensive simulation results are presented that show that the proposed RVIN suppression mechanism is superior compared to the state-of-the-art schemes.

2. Proposed Noise Detection Scheme

In this section, we present the proposed efficient noise detector for the detection of noisy pixels. Before the detection process, window size is adjusted for the whole image. The window size is selected adaptively as it is set automatically based on the noise density. To estimate the noise density, we adopt the idea of the universal threshold [24]. In this method, wavelet transform is used to estimate the noise density σ of the noisy signal using Robust Median Estimator as:
σ = median ( H F ) 0.6745 ,
where H F is the high frequency sub-band. In this technique, the noisy image or signal is divided into sub-bands of low and high frequencies. The edges, sharp details and sudden changes are found in the high frequencies [25]. Thus, mostly, the high frequency components contain noise. Therefore, only a high frequency sub-band ( H F ) is used to estimate the noise of the signal, as illustrated in Equation (2). For example, if σ comes to be 50, the image has 50% noise and so on. The size of the window depends on σ , such that, for low σ , a smaller window is selected, while a larger window is selected for high σ .
Once the window is adjusted, then the detection process starts. Figure 1 illustrates the block diagram of the proposed noise detector. Each block given in Figure 1 is then elaborated with example in the following subsections.

2.1. Calculation of Separating Threshold

Once the window size is defined, the first step of the proposed noise detection process is the selection of separating threshold γ s . The algorithm to find separating threshold is illustrated with the help of example in Figure 2 and elaborated in the following steps.
  • Let x i , j be the pixel under observation at location ( i , j ) to know whether it is noise or clean, then it is considered as a central pixel. Window of the size found in Equation (2) is set around the central pixel.
  • The corresponding pixel values of of the current window are sorted in ascending order in vector s . The outliers are removed from s by adopting the “three sigma rule” with median absolute deviation (MAD) [21]. The MAD for s is found as:
    MAD = median s median ( s ) .
    Vector s is searched, and, if the difference between any value in s and the median value is greater than 3 × MAD , it is considered as an outlier and is not taken into account for further calculation. After removal of outliers, let s n be the resultant vector, where n denotes the size of s n .
  • Compute the absolute difference d k between every two adjacent pixels in s n as:
    d k = s k s k + 1 ,
    where k = 1 , 2 , , n 1 .
  • The separating threshold γ s is found as:
    γ s = m + e 1 σ / 100 ,
    where
    m = 1 n 1 k = 1 n 1 d ( k ) .
We can observe that the term e 1 σ / 100 depends upon the noise density σ . It can be observed that, for lower noise density, this term is large and, for higher noise density, this term is small. This term adjusts the threshold to ensure that clean pixels are not considered noisy pixels at very low noise density. Similarly, it also adjusts the threshold to make sure noisy pixels are not grouped into clean pixels, as this threshold will be used in grouping discussed in the following subsection.

2.2. Separation of Pixels into Clusters

Using the separating threshold γ s , in this section, we divide the pixels in the sorted vector s n in groups as follows:
  • Initialize from the first pixel in sorted vector s n . If the difference between current pixel s i and next pixel s i + 1 is less than or equal to γ s , keep the pixel s i + 1 in the current group together with s i .
  • However, if the difference is greater than γ s , finish the current group and start a new group, such that the first element of the new group will be s i + 1 . Apply the same procedure until the last element of vector s n .
  • The clean and noisy groups are decided based on the size of the group. The first two largest clusters are considered as clean clusters while all other clusters are considered as noisy and are discarded. The clean clusters are denoted as g 1 and g 2 . Cluster g 1 is used in the detection phase while g 1 and g 2 are used in filtering.

2.3. Noise Detection

Using the results of the group stage, in this section, the noise detection stage is presented. The noise detector is:
x i , j = Clean , if x i , j μ γ , Noisy , Otherwise ,
where threshold γ is
γ = 3 × MAD ( g 1 ) ,
and g 1 is the largest group discussed in the grouping stage.
Similarly, x i , j is the central pixel and μ is the mean value of g 1 . If the x i , j is detected as noisy, filtering is required to replace the x i , j with the filtered output. On the other hand, if x i , j is a clean pixel, it is left unchanged.

2.4. Edge Pixel Identification

In this section, we present a mechanism to avoid the faulty detection of edge pixels as noisy pixels. This step is performed only when the pixel is identified as a noisy pixel. This will make the detection scheme more robust by verification of the noisy pixels to avoid situations where the edge pixel is wrongly considered as a noisy pixel. This process is performed as follows:
  • Keeping x i , j as a reference pixel, pixels present along four directions are taken into account. These four directions θ A , θ B , θ C and θ D are illustrated in Figure 3.
  • θ A contains the pixels from the direction θ A excluding x i , j . Similarly, θ B , θ C and θ D contain the pixels from the directions θ B , θ C and θ D , respectively, excluding x i , j .
  • Find the absolute difference between central pixel x i , j with all pixels along direction θ A . The difference vector is defined as d A = x i , j θ A .
  • d A is searched to find the values less than or equal to m. Indexes of pixels corresponding to these values are identified in θ A . These pixels are stacked in vector v . The vector v is considered to be an edge if total number of pixels in v is at least 50% of direction θ A . Step 4 is only performed if this condition holds true. Otherwise, discard the v and go to step 5.
  • Take the standard deviation of central pixel x i , j with v . The pixel is identified as an edge if standard deviation of e and x i , j is less than or equal to threshold γ :
    x i , j = Edge , if std ( v , x i , j ) γ , Go to next Direction , Otherwise .
  • If x i , j is not identified an edge in direction θ A , then the three other directions θ B , θ C and θ D are checked for an edge identification by applying the five steps above on each direction.
If x i , j is not identified an edge in all four directions, it is considered a noisy pixel and filtering is required to de-noise it.

3. Illustration

This section illustrates the noise detection of the proposed method with example. Consider the image patch given in Figure 4 that is taken from a 40% noisy “lena” image. The size of the window depends on the noisy density σ such that:
Size of window = 7 × 7 , if σ 60 , 5 × 5 , otherwise .
To check whether the specified central pixel is noisy, the following steps are performed.
  • Using Equation (2), the noise density of the whole image is calculated. As the noisy density is 40%, therefore, according to Equation (10), the window size is selected 5 × 5.
  • Sort all the pixels inside the selected patch in ascending order as x = [17 64 64 65 65 65 65 66 66 67 68 68 69 69 83 84 112 112 133 143 152 199 227 236 252].
  • The extreme outliers are removed by using MAD given in Equation (3). The MAD value of vector x is 5, while the median value of x is 69. Therefore, according to the “Three Sigma Rule”, any value below 54 (69 − 15) and above 84 (69 + 15) is discarded and the vector x becomes [64 64 65 65 65 65 66 66 67 68 68 69 69 83 84].
  • Next, the deference between every two consecutive pixels is found to be d = [0 0 0 0 1 0 1 1 0 1 0 14 2].
  • The average value m of the vector d is equal to 1.4. Putting this value of m in Equation (5), γ s becomes 13.59.
  • Using the separating threshold, all pixels of x are divided into groups as [64 64 65 65 65 65 66 66 67 68 68 69 69], [83 84].
  • The largest group is [64 64 65 65 65 65 66 66 67 68 68 69 69], whose size is much greater compared to other group. Thus, this group is considered the clean group, while the other group is considered the noisy group.
  • Using Equation (8), threshold γ becomes 3. Similarly, mean value becomes μ = 66.
  • By employing Equation(9), any value greater than 69 and less than 63 is considered to be noisy.
In this example, the value of the central pixel is 67; therefore, it is a clean pixel and does not need filtering. In case the pixel comes out to the noisy pixel, then we need filtering to replace the noisy pixel with estimated clean pixel. In the next section, we present the filtering scheme.

4. Proposed Filtering Scheme

After the detection of the noise, the appropriate filter is required to replace the noisy pixel. For this, we design a novel Weighted Mean Filter (WMF), which is more robust compared to the existing median and mean filters. In a smooth region, the absolute difference among all neighbor pixels is very small and the whole region is flat [26], so, in most cases, only one large clean group is formed. On the other hand, the region that contains the edge is not as flat, in which case more than one cluster is required to replace the noisy pixel. Therefore, in the proposed filtering scheme, the second largest group g 2 is also considered. We remark that the g 2 is not required in detection stage. Instead, there is a separate section for edge detection. The two largest groups are assigned weights based on standard deviation and size of the group. Supposing that x i , j is the noisy pixel at location ( i , j ) , then the filtered pixel x ^ i , j is given as:
x ^ i , j = mean { w 1 g 1 , w 2 g 2 } ,
where ⋄ is the replication operator [27] that shows the repetition of any term. For example, if w 1 is 2, and w 2 is 1, then g 1 is repeated two times while g 2 is repeated 1 time, such as x ^ i , j = mean { g 1 , g 1 , g 2 } , as the size of g 1 is greater than g 2 . Therefore, to decide the weights, we also consider the standard deviation of both groups. The weights w 1 and w 2 are determined by a ratio R as:
R = s 2 s 1 × sd 1 sd 2 ,
where s 1 , s 2 represents the size, while sd 1 and sd 2 represent the standard deviation of g 1 and g 2 , respectively. This relation defines the weights of g 2 with respect to g 1 . As the higher standard deviation means less correlation, it is inversely proportional to the R, while the larger size shows larger priority; therefore, it is in direct relation to R. Based on the value of R, the following weights are assigned in Equation (11).
  • If the size of g 1 is very large compared to the size of g 2 or R ≤ 0.5, then w 1 = 1, w 2 = 0.
  • Similarly, when 0.5 R 1, the assigned weights are: w 1 = 1, w 2 = 1.
  • Finally, for 1 R 2, the assigned weights are: w 1 = 1, w 2 = 2.
The filtering is further improved by only considering the closest neighbor pixels of the central pixel from the selected groups, such that, if there are at least three pixels from the selected groups which are in a 3 × 3 window around the central pixel, then take the mean value of only these closest pixels.

5. Summary of the Proposed Denoising Scheme

Although the proposed scheme gives good results in single iterations, with second iterations, it gives the best performance and there is no need to apply further iterations as in other iterative methods. The larger threshold value is used in the first iteration while a smaller threshold value is used in the second iteration [2]. Therefore, if the threshold in the first iteration is 3 × M A D , in the second iteration, the threshold is set to 1.5 × M A D . We summarize our proposed technique in Algorithm 1 as:
Algorithm 1: Proposed algorithm for random valued impulse noise removal.
Input: The noisy image X, with size of L 1 × L 2
1 Set iteration i = 1.
2 Outer loop: for l 1 = 1 , L 1 , Inner loop: for l 2 = 1 , L 2
3 Set a sliding window centered at x i , j , the size of window is determined by σ .
4 Divide the window into groups of clean and noisy pixels using separating threshold computed in Equation (5).
5 Find mean μ and threshold γ from clean cluster, while discard the noisy clusters.
6 Using Equation (7) check x i , j for noise. If x i , j μ γ , go to step (2).
7 If the pixel is identified as noisy, check x i , j for an edge using Equation (9). Such that if std ( e , x i , j ) γ , go to step (2).
8 Filter x i , j as in Equation (11).
9 i = i + 1, if i > i m a x , stop iteration, else X ⟵ Y, go to next iteration.
Output: De-noised image Y

6. Experimental Results

To assess the capability the proposed noise removal procedure, in this section, comprehensive numerical results are compared with other state-of-the-art techniques.

6.1. Comparison of Noise Detection

We compare our proposed noise detection scheme with different state-of-the-art techniques such as (ACWM) [15], the (DWM) [19], (TSM) [9] (PSMF [17], (Luo) [18], (ROR-NLM) [21], SBF [28] and (ASWM) [23].
In Table 1, the noise detection results of noisy “lena” image with noise densities 40%, 50% and 60% are given. The table lists three results of the noise detection: the number of undetected noisy pixels (False Negative), the number of clean pixels that are falsely detected as noisy (False Positive) and the Total number (False Negative + False Positive).
From Table 1, it is clear that the proposed detection scheme has generated less ‘Total’ number than all of the above techniques especially at higher noise density. At lower noise density, only ROR-NLM produces a lower total number, but they produce much higher false negative terms, indicating that they miss most of the noisy pixels, and most of noise patches are left unchanged. Thus, the proposed detection technique shows superiority over other detection schemes. It is because a noise detector is considered to be efficient and robust, if it reduces both false positive and false negative numbers. When the noise density increases, the superiority of the proposed method is clear as it generates far less ‘Total’ terms than all other methods. This means that our technique is more efficient and robust, even at high noise density.

6.2. Comparison of Image Restoration

The de-noising performance of the proposed method is tested on different standard images such as Lena, pepper, Bridge, and Boat with different noise densities. All these images are 8-bit gray scale and have the same size of 512 × 512 pixels. It should be noted that the size of the window is defined only once based on the noise density using Equation (2), which is used on all pixels of the current image. In our case, for noise density less than 60%, window size of 5 × 5 is used, while window size of 7 × 7 is used otherwise. Peak Signal to Noise Ratio (PSNR) is used as quantitative measurement to compare the proposed method with different well known techniques. The PSNR is defined as:
PSNR ( dB ) = 10 log 10 255 2 MSE ,
where MSE is mean squared error and is defined as:
MSE = 1 M N i = 1 m j = 1 n O i , j R i , j 2 ,
where O is the original clean image while R is the output restored image. Table 2 and Table 3 illustrate comparison of the proposed method with other state-of-the-art methods.
Table 2 and Table 3 list the PSNR values of different test images that are contaminated by RVIN of noise densities ranges from 40% to 60%. From both tables, it is clear that our method has performed better as compared to all other listed methods in terms of PSNR. Figure 5 and Figure 6 show de-noised results of Lena and Pepper images, respectively. From these figures, it is obvious that TSM, ACWM, PSM and Luo have created artifacts and damaged image details. While DWM, ASWM and ROR-NLM have performed better at lower noise densities, but their performance is degraded at higher noise densities. This is due to the fact that these filters replace the noisy pixel by considering all the window pixels. On the other hand, the proposed method only considers clean pixels to replace the noisy pixel. As the noise density increases, the number of corrupted pixels in the window are also increased, which increases the error when the central noisy pixel is replaced by mean or median value of noisy neighbor pixels. Finally, in Figure 7, we present performance results of different filters in restoring Lena image corrupted by different noise densities. These results confirm the results given in Table 2 and Table 3 and clearly show that the proposed method outperforms other state of the art techniques.

7. Conclusions

In this paper, a new algorithm is presented to remove random valued impulse noise. The proposed approach is based on calculation of new and adaptive thresholds for every pixel. The pixels are divided into noisy and clean clusters. Then, only a clean cluster is used to calculate the threshold. A new edge identification step is proposed to ensure that a clean pixel is not considered as a noisy pixel falsely. If the pixel is identified as noisy, a novel weighted mean filter is designed to filter the noisy pixel. Instead of taking the mean of all pixels in the selected window, the proposed filter considers only clean pixels to replace the noisy pixel. Simulation results show that the proposed method outperforms state-of-the-art de-noising techniques, both visually and quantitatively (PSNR), even when the noise density is high.

Author Contributions

Conceptualization, S.A., I.K. and B.M.L.; Data curation, N.I. and I.K.; Formal analysis, S.A., I.K. and N.I.; Funding acquisition, I.K. and B.M.L.; Investigation, S.A., I.K. and N.I.; Methodology, S.A., I.K., B.M.L. and N.I.; Project administration, N.I. and I.K.; Software, S.A., I.K. and N.I.; Writing—original draft, S.A., I.K. and B.M.L.; Writing—review & editing, S.A., I.K., B.M.L. and N.I.

Funding

This work was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (Grant No.: NRF-2017R1D1A1B03028350).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gao, G.; Liu, Y.; Labate, D. A two-stage shearlet-based approach for the removal of random-valued impulse noise in images. J. Vis. Commun. Image Represent. 2015, 32, 83–94. [Google Scholar] [CrossRef]
  2. Awad, A.S. Standard deviation for obtaining the optimal direction in the removal of impulse noise. IEEE Signal Process. Lett. 2011, 18, 407–410. [Google Scholar] [CrossRef]
  3. Ilke, T. A new method to remove random-valued impulse noise in images. AEU-Int. J. Electron. Commun. 2013, 67, 771–779. [Google Scholar]
  4. Chan, R.H.; Hu, C.; Nikolova, M. An iterative procedure for removing random-valued impulse noise. IEEE Signal Process. Lett. 2004, 11, 921–924. [Google Scholar] [CrossRef]
  5. Jayasree, S.; Bodduna, K.; Pattnaik, P.K.; Siddavatam, R. An expeditious cum efficient algorithm for salt-and-pepper noise removal and edge-detail preservation using cardinal spline interpolation. J. Vis. Commun. Image Represent. 2014, 25, 1349–1365. [Google Scholar] [CrossRef]
  6. Gonzalez, R.C.; Woods, R.E. Digital Image Processing; Prentice Hall: Upper Saddle River, NJ, USA, 2002. [Google Scholar]
  7. Astola, J.; Kuosmaneen, P. Fundamentals of Nonlinear Digital Filtering; CRC: BocaRaton, FL, USA, 1997. [Google Scholar]
  8. Singh, N.; Thilagavathy, T.; Lakshmipriya, R.T.; Umamaheswari, O. Some studies on detection and filtering algorithms for the removal of random valued impulse noise. IET Image Process. 2017, 11, 953–963. [Google Scholar] [CrossRef]
  9. Chen, T.; Ma, K.; Chen, L. Tri-state median filter for image denoising. IEEE Trans. Image Process. 1999, 8, 1834–1838. [Google Scholar] [CrossRef] [PubMed]
  10. Arce, G.R.; Paredes, J.L. Recursive weighted median filters admitting negative weights and their optimization. IEEE Trans. Signal Process. 2000, 48, 768–779. [Google Scholar] [CrossRef]
  11. Ari, N.; Heinonen, P.; Neuvo, Y. A new class of detail-preserving filters for image processing. IEEE Trans. Pattern Anal. Mach. Intell. 1987, 1, 74–90. [Google Scholar]
  12. Ko, S.-J.; Lee, Y.H. Center weighted median filters and their applications to image enhancement. IEEE Trans. Circuits Syst. 1991, 38, 984–993. [Google Scholar] [CrossRef]
  13. Coyle, E.J.; Lin, J.; Gabbouj, M. Optimal stack filtering and the estimation and structural approaches to image processing. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 2037–2066. [Google Scholar] [CrossRef]
  14. Vikas, G.; Chaurasia, V.; Shandilya, M. Random-valued impulse noise removal using adaptive dual threshold median filter. J. Vis. Commun. Image Represent. 2015, 26, 296–304. [Google Scholar]
  15. Chen, T.; Hong, R.W. Adaptive impulse detection using center-weighted median filters. IEEE Signal Process. Lett. 2001, 8, 1–3. [Google Scholar] [CrossRef]
  16. Lin, T.-C. Switching-based filter based on Dempster’s combination rule for image processing. Inf. Sci. 2010, 180, 4892–4908. [Google Scholar] [CrossRef]
  17. Zhou, W.; Zhang, D. Progressive switching median filter for the removal of impulse noise from highly corrupted images. IEEE Trans. Circuits Syst. II Analog Digit. Signal Process. 1999, 46, 78–80. [Google Scholar] [CrossRef]
  18. Luo, W. A new efficient impulse detection algorithm for the removal of impulse noise. IEICE Trans. Fund. Electron. Commun. Comput. Sci. 2005, 88, 2579–2586. [Google Scholar] [CrossRef]
  19. Dong, Y.; Xu, S. A new directional weighted median filter for removal of random-valued impulse noise. IEEE Signal Process. Lett. 2007, 14, 193–196. [Google Scholar] [CrossRef]
  20. Ghanekar, U.; Singh, A.K.; Pandey, R. A contrast enhancement-based filter forremoval of random valued impulse noise. IEEE Signal Process. Lett. 2010, 1, 47–50. [Google Scholar]
  21. Bo, X.; Yin, Z. A universal denoising framework with a new impulse detector and nonlocal means. IEEE Trans. Image Process. 2012, 21, 1663–1675. [Google Scholar]
  22. Antoni, B.; Coll, B.; Morel, J.-M. A non-local 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. [Google Scholar]
  23. Akkoul, S.; Ledee, R.; Leconge, R.; Harba, R. A new adaptive switching median filter. IEEE Signal Process. Lett. 2010, 17, 587–590. [Google Scholar] [CrossRef]
  24. Johnstone, I.M.; Silverman, B.W. Wavelet threshold estimators for data with correlated noise. J. R. Stat. Soc. Ser. B (Stat. Methodol.) 1997, 59, 319–351. [Google Scholar] [CrossRef]
  25. Cihan, T.; Akinlar, C. Edge drawing: a combined real-time edge and segment detector. J. Vis. Commun. Image Represent. 2012, 23, 862–872. [Google Scholar]
  26. Chen, J.; Zhan, Y.; Cao, H.; Wu, X. Adaptive probability filter for removing salt and pepper noises. IET Image Process. 2018, 12, 863–871. [Google Scholar] [CrossRef]
  27. Bovik, A.C. Handbook of Image and Video Processing; Academic Press: Cambridge, MA, USA, 2010. [Google Scholar]
  28. Lin, C.H.; Jia, S.T.; Ching, T.C. Switching bilateral filter with a texture/noise detector for universal noise removal. IEEE Signal Process. Soc. 2010, 19, 2307–2320. [Google Scholar]
Figure 1. Proposed noise detection scheme.
Figure 1. Proposed noise detection scheme.
Symmetry 11 00395 g001
Figure 2. Separating threshold calculation.
Figure 2. Separating threshold calculation.
Symmetry 11 00395 g002
Figure 3. Edges in different directions (a) θ A ; (b) θ B ; (c) θ C ; (d) θ D .
Figure 3. Edges in different directions (a) θ A ; (b) θ B ; (c) θ C ; (d) θ D .
Symmetry 11 00395 g003
Figure 4. A patch from Lena image: (a) original image; (b) 40% noisy image.
Figure 4. A patch from Lena image: (a) original image; (b) 40% noisy image.
Symmetry 11 00395 g004
Figure 5. Results of different filters in restoring 40 % corrupted Lena image: (a) adaptive switching median filter (ASWM); (b) tri-state median filter (TSM); (c) progressive switching median filter (PSM); (d) luo-iterative median filter (Luo); (e) directional weighted median filter (DWM); (f) Adaptive Switching Median Filter (ASWM); (g) Switching Bilateral Filter (SBF); (h) Robust Outlyingness Ratio Non-Local Mean (ROR-NLM); (i) proposed method.
Figure 5. Results of different filters in restoring 40 % corrupted Lena image: (a) adaptive switching median filter (ASWM); (b) tri-state median filter (TSM); (c) progressive switching median filter (PSM); (d) luo-iterative median filter (Luo); (e) directional weighted median filter (DWM); (f) Adaptive Switching Median Filter (ASWM); (g) Switching Bilateral Filter (SBF); (h) Robust Outlyingness Ratio Non-Local Mean (ROR-NLM); (i) proposed method.
Symmetry 11 00395 g005
Figure 6. Results of different filters in restoring 60 % corrupted pepper image: (a) adaptive switching median filter (ASWM); (b) tri-state median filter (TSM); (c) progressive switching median filter (PSM); (d) luo-iterative median filter (Luo); (e) directional weighted median filter (DWM); (f) adaptive switching median filter (ASWM); (g) switching bilateral filter (SBF); (h) robust outlyingness ratio non-local mean (ROR-NLM); (i) proposed method.
Figure 6. Results of different filters in restoring 60 % corrupted pepper image: (a) adaptive switching median filter (ASWM); (b) tri-state median filter (TSM); (c) progressive switching median filter (PSM); (d) luo-iterative median filter (Luo); (e) directional weighted median filter (DWM); (f) adaptive switching median filter (ASWM); (g) switching bilateral filter (SBF); (h) robust outlyingness ratio non-local mean (ROR-NLM); (i) proposed method.
Symmetry 11 00395 g006
Figure 7. Results of different filters in restoring Lena image corrupted by different noise levels.
Figure 7. Results of different filters in restoring Lena image corrupted by different noise levels.
Symmetry 11 00395 g007
Table 1. Detection comparison of different filters at different noise densities.
Table 1. Detection comparison of different filters at different noise densities.
Methods 40% 50% 60%
FNFPTotalFNFPTotalFNFPTotal
ACWM14,590236716,85721,897370625,60330,198652636,724
Luo’s14,679183116,51021,665306824,73333,987278036,767
TSM18,921520124,12223,921621830,13928,123890337,026
PSM18,672498223,65423,762612329,88527,987839336,380
ASWM748911,56419,05311,77912,78624,56519,98216,48236,464
DWM11,786893120,71715,321875624,07715,72814,81630,544
SD-OOD13,50010,67524,17511,98715,82727,81417,82118,26136,082
ROR-NLM12,890332816,21815,297348718,78421,827780829,635
Propoposed10,908797318,88111,668961321,24113,571976023,331
Table 2. Peak Signal-to-Noise Ratio (dB) values of different filters for Lena and Pepper image corrupted by random valued impulse noise of different noise densities.
Table 2. Peak Signal-to-Noise Ratio (dB) values of different filters for Lena and Pepper image corrupted by random valued impulse noise of different noise densities.
MethodsLenaPepper
(Noise Density)40%50%60%40%50%60%
TSM24.9122.2218.9824.0222.5617.78
ACWM28.1224.6420.3228.1226.0921.52
LOU’S29.6225.8222.6928.4126.3323.89
PSM29.1225.3124.2328.2325.4123.83
ASWM30.9128.4226.2529.0227.1225.33
DWM30.9228.9126.5129.1327.7925.48
SBF30.1227.1223.1328.8927.0224.67
ROR-NLM31.4229.2125.6129.6127.8925.45
PROPOSED31.7730.0128.0329.7528.1126.62
Table 3. Peak Signal -to- Noise Ratio (dB) values of different filters for Bridge and Boat image corrupted by random valued impulse noise of different noise densities.
Table 3. Peak Signal -to- Noise Ratio (dB) values of different filters for Bridge and Boat image corrupted by random valued impulse noise of different noise densities.
MethodsBridgeBoat
(Noise Density)40%50%60%40%50%60%
TSM21.5519.7217.2622.8920.7217.55
ACWM23.7222.1919.1225.3223.5621.45
LOU’S23.8422.7819.1725.4523.7821.61
PSM23.6521.9119.3924.9822.9820.87
ASWM24.0222.6921.1226.7625.2623.23
DWM24.0722.7221.1926.8325.3123.57
SBF22.1221.3120.1525.3324.8822.67
ROR-NLM24.2722.9121.2127.2325.4324.21
PROPOSED24.3523.0821.7527.8526.6124.87

Share and Cite

MDPI and ACS Style

Iqbal, N.; Ali, S.; Khan, I.; Lee, B.M. Adaptive Edge Preserving Weighted Mean Filter for Removing Random-Valued Impulse Noise. Symmetry 2019, 11, 395. https://doi.org/10.3390/sym11030395

AMA Style

Iqbal N, Ali S, Khan I, Lee BM. Adaptive Edge Preserving Weighted Mean Filter for Removing Random-Valued Impulse Noise. Symmetry. 2019; 11(3):395. https://doi.org/10.3390/sym11030395

Chicago/Turabian Style

Iqbal, Nasar, Sadiq Ali, Imran Khan, and Byung Moo Lee. 2019. "Adaptive Edge Preserving Weighted Mean Filter for Removing Random-Valued Impulse Noise" Symmetry 11, no. 3: 395. https://doi.org/10.3390/sym11030395

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