 Research
 Open Access
 Published:
An adaptive decision based interpolation scheme for the removal of high density salt and pepper noise in images
EURASIP Journal on Image and Video Processing volume 2017, Article number: 67 (2017)
Abstract
An adaptive decision based inverse distance weighted interpolation (DBIDWI) algorithm for the elimination of highdensity salt and pepper noise in images is proposed. The pixel is initially checked for salt and pepper noise. If classified as noisy pixel, replace it with an inverse distance weighted interpolation value. This interpolation estimates the values of corrupted pixels using the distance and values nearby nonnoisy pixels in vicinity. Inverse distance weighted interpolation uses the contribution of nonnoisy pixel to the interpolated value. The window size is varied adaptively depending upon the nonnoisy content of the current processing window. The algorithm is tested on various images and found to exhibit good results both in terms of quantitative (PSNR, MSE, SSIM, Pratt’s FOM) and qualitative (visually) at high noise densities. The algorithm performs very well in restoring an image corrupted by highdensity salt and pepper noise by preserving fine details of an image.
Introduction
The term interpolation comes from the topic of resampling described by Hanumantharaju et al. By definition, resampling [1] is a process of transforming a discrete image that is indicated at one particular set of coordinate locations to a new set of coordinate points. In specific, the process of interpolation refers to estimating unknown pixel values by using known pixels, such that the estimated pixel is close to the original pixel. The details usually denotes coordinates, color, gray level, or density with the image having any dimensions. Image interpolation finds many leads in the field of image analysis such as zooming, given by Hanumantharaju et al. [1] and Min et al., [2] resolution enhancement, image impainting, image warping, given by Wing et al. [3] and Bender et al., [4] and geometric transformations. This paper gives a framework for the removal of salt and pepper noise using interpolation methods. Salt and pepper noise is often induced in an image due to faulty camera sensors or due to transmission errors. Filters are the obvious choice for the noise removal. Over the decades, many filters were formulated for salt and pepper noise (SPN) removal, of which nonlinear filters were good in removing the noise by preserving the edges. Hence, various nonlinear filters were proposed. The standard median filter [5] (SMF) eliminates salt and pepper noise by preserving fine details of an image as described by Huang et al. The main drawback of the SMF is that all the pixels of the images are replaced by median, irrespective whether the pixel is noisy or not. Hence, an adaptive window [6] for increasing noise densities was proposed by Hwang et al. The size of the window affects the performance of the image. A smaller window size does not have sufficient information for noise removal, and on the contrary, a larger window blurs the image. Several special filters such as centerweighted median filter [7] by Ko et al., weighted median filter [8] given by Brownrigg et al., and tristate median filter [9] proposed by Chen et al. were introduced to eliminate fixedvalued salt and pepper noise. These filters either do not restore the pixels or preserve edges for high noise densities. Progressive switching median filter [10] given by Wang et al. uses switching scheme for detection of impulses, thus only proportion of all the images was filtered. This algorithm flutters at high noise densities. Eng et al. proposed a switching median filter [11] which identifies the corrupted pixel by using a decision tree structure and replace it with any of the three filters (all pass, standard median, and a weighted fuzzy filter). The performance of the algorithm was poor at high noise densities. Raymond et al. introduced an algorithm [12] that works in two phases. Initially, the algorithm uses adaptive filter for estimating the pixels if it is noisy or not. In the second phase, a regulation method was applied to the noisy pixels which resulted in noise suppression with preservation of image details. Decision based algorithm [13] was proposed by Srinivasan et al. to eliminate high density impulse noise. This algorithm replaces the corrupted pixel with median value or the preprocessed neighbor. This repeated replacement of preprocessed neighborhood results in streaking. The effect of streaking was minimized using an improved decision based algorithm [14] given by Madhu et al. that used mean of preprocessed pixel which resulted in reduced streaks. A cascaded filters [15] were proposed by Balasubramanian et al. for high density salt and pepper noise. This algorithm initially applies decision based median filter as first stage, and asymmetrical trimmed median and midpoint were applied on the later stage. This algorithm exhibits smoothing effect at high noise densities. A novel decision based asymmetric trimmed median filter [16] was proposed by Aiswarya et al. which used reduced computation for calculation of median by asymmetrical trimming. This algorithm exhibits fading at high noise densities. The performance of asymmetrical trimmed median [17] by Esakkiarajan et al. was improved by finding mean of the window, which was entirely corrupted by salt and pepper noise. The above algorithm was improved further, given by VeeraKumar et al. [18], by replacing the corrupted pixel with trimmed global mean if the entire window is either 0 or 255. A decision based asymmetrical trimmed variants [19] proposed by Vasanth et al. (2012) replaced the corrupted pixel with asymmetrical trimmed midpoint based on the content of the current processing window. An adaptive cardinal Bspline algorithm [20] (ACBSA) given by Syamala Jeyashree et al. exploited the interpolation property of Bspline to estimate a pixel value to replace the corrupted pixel. The adaptive algorithm identifies the corrupted pixel and operates based on the number of nonnoisy pixels in the current processing window. At high noise densities, the edges do smudge. The Recursive Spline Interpolation Filter (RSIF) given by Veerakumar et al. [21] uses vicinity based noisefree pixels and preprocessed nonnoisy output. The algorithm too exhibits fading effect at high noise densities. A decision based neighborhoodreferred unsymmetrical trimmed variants [22] (DBNRUTVF), Vasanth et al. (2014) uses mean of the four neighbors or asymmetrical trimmed median or asymmetrical trimmed midpoint or the global trimmed mean for the replacement of the corrupted pixels. Hence, few nonlinear operations such as median, spline interpolation, nonlinear mean, asymmetrical trimmed median or asymmetrical trimmed midpoint, polynomial, bilinear, bicubic, etc. were used for the elimination of salt and pepper noise in images and videos. It was observed from the literatures that at high noise densities, the algorithm flatters or induce few ambiguities during noise removal. Peixuan et al. increased the adaptive size of the window [23] by increasing the window size until the two subsequent windows have the same outlier values. The corrupted pixel is replaced by weighted mean value of the current processing window. The adaptive window grows up to 39. At high noise densities, the restored image is attenuated. Hence, a proper noise detection and correction algorithm has to be used for a better result. The organization of the paper is as follows: Section 2 deals with existing interpolation techniques and decision based inverse distance weighted interpolation filter. Section 3 briefs the proposed algorithm. Section 4 gives the methodology of the interpolation technique. Section 5 gives the simulation results and discussions. Section 6 deals with the conclusion of the proposed work.
Interpolation Techniques
An interpolation algorithm is subdivided into adaptive and nonadaptive. The adaptive algorithm interpolation will work based on the user interest in images (information content (edges) or smooth texture). QImage, Photozoom pro, and genuine fractals are few renowned adaptive algorithms. But in the case of nonadaptive algorithm, all the pixels are treated equally. Nearest neighbor, linear, sinc, bilinear, Lanczos, and bicubic are few nonadaptive interpolating algorithms given by Hanumantharaju et al. The nearest neighbor is a simple interpolation method. This method chooses a pixel that is very close to the evaluated pixel. This method is simple but damages the straight edges and produce aliasing and blurring effects in images. Bilinear interpolation chooses closest four neighborhood of known pixel values in the vicinity of the unknown pixels. For the final interpolated value, weighted average of the closest four pixel is chosen. This interpolation results in smoothened images but computationally slower. Bicubic interpolation employs 16 closest neighbors of the known pixels, since these 16 pixels are at different distances from the estimated pixel. While in calculation, it was found that higher weighting is given to the pixels that are very close. When compared to the nearest neighborhood and bilinear method, bicubic interpolation retains the fine details well. Sinc interpolation can be usable for spatial convolution if it uses a smaller window and proper truncation. In order to approximate sinc functions, the following approximations were made such as nearest neighbors, bilinear, quadratic approximations, and Bspline approximations. These methods also induce strong blurring effects. After performing interpolation, sharp edge details and high local contrast are more affected. Bspline interpolation produces good results in terms of structural similarity in comparison to the initial image. Interpolation in the third order polynomial was found useful in many applications. The basic criteria for a good interpolation method is geometric invariance, contrast invariance, no noise, edge preservation, no aliasing (jagged or staircase edge), should preserve texture, should not over smoothen, and sensitive to specified (power weight) parameters ([24, 25]). In recent years, interpolation technique is used for the removal of salt and pepper noise.
Inverse distance weighted interpolation
This interpolation estimates the pixel values at corrupted pixel locations in an image using the distance and values nearby known pixels in confined vicinity. Inverse distance weighted interpolation (IDWI) can be used for data extraction or smoothing. IDWI reduces the contribution of a known pixel to the interpolated value. In this interpolation technique, weight of each sample pixel is a reciprocal of the distance. Hence, the impact of influencing pixel value with respect to other pixel diminishes with distance from the estimated pixel. IDW interpolator is good and efficient when we have a fairy good number of pixels in the closer vicinity. The weights and estimation formulae are given in Eqs. 1 and 2. The key factors influencing the effectiveness of the interpolating technique [26] are the power parameter and number of uncorrupted pixels (n) given in Isaaks et al. The choice of the spatial kernel size [27] is chosen arbitrary by Webster et al. There are two factors influencing the weights. They are number of uncorrupted pixel in the current processing window “n” and the power parameter “p”. If the number of nonnoisy pixels inside the current processing window “n” is more, weights will have less impact and uncorrupted pixel will have more impact on the estimated value. On the contradictory, if the n value decreases, then the weights will have more impact and the uncorrupted pixels will have less impact on the estimated values. Hence, for the inverse distance weighted interpolator to yield good results, “n” should contain at least three nonnoisy pixels for the filtering process. The impact of power and choice of its value are discussed in Section 2.2. The important reason for selecting inverse distance weighted interpolation is because it is an exact interpolator. The weight obtained diminishes with increase in distance. Hence, a closer distance between pixel values yield better interpolated results. Weighted functions can be controlled by the user. Inverse distance were used for weight generation and these weights combine with nonnoisy pixels to form the interpolated value which is used for salt and pepper noise removal in images
d _{ i } is the distance between nonnoisy pixel (g) and the pixel to be processed, p is a power parameter, and n represents the number of nonnoisy pixels in the current processing window used for the estimation and v gives the calculated weights.
Impact of “p” parameter
The impact of “p” parameter plays a vital role in restoration of images corrupted by salt and pepper noise. In general, adaptive window blurs the image. The blurring of the adaptive interpolated image increases as the power parameter increases. When the power p is equal to zero, then IDWI acts as moving average interpolator. When the power p is increased from 0 to 1, then IDWI acts as weighted mean interpolator, given by Brus et al. and Hosseini et al., and when p is greater than 1, IDWI acts as weighted average interpolator but the impact of uncorrupted pixel will be less for interpolation, given by Abdou et al. But the value of “p” will vary from image to image. Hence, for choosing the value of “p,” the edge information of different image is taken into account. Exhaustive experiments were conducted on different images by varying the value of “p.” It was found that value of p that lie between 0.8 and 2 yield very low mean square errors for all the images. Hence, for any images, the value of p will be between 0.8 and 2. If the image information content is less, then higher p value will be chosen. If the image information content is greater, then a very less p value will be chosen. If p is greater, the effect of pixels will be less. If p is small, the weights are evenly distributed among the neighboring data points.
Selection of window size
For the IDWI interpolation to work, we require at least three nonnoisy pixels in any given window. Hence, selection of window size plays an important role in this adaptive technique. For processing any given pixel in an image, a 3 × 3 window is used initially. As we require at least three nonnoisy pixels for the IDWI calculation, the size of the window increases depending on the presence of nonnoisy pixel in the current processing window (i.e., if the initial window size is 3 × 3 and we do not have three nonnoisy pixels, then the window size increases). For example, if the number of nonnoisy pixel is 3, then 3 × 3 window is used for calculation of weights. If number of nonnoisy pixel is 2 then the window size is increased by 2 (5 × 5 for this case) from the initial size. In this case, if the current processing window (3 × 3) has only one nonnoisy pixel, then the window size is increased by 4 (7 × 7 for this case). If there is no nonnoisy pixel in the current processing window (3 × 3), then the window size is increased by 6 (9 × 9 for this case). After increasing the window size, the algorithm counts the number of nonnoisy pixel in the new window. For example, if a 3 × 3 window has one nonnoisy pixel, then the window size is increased to 7 × 7. If there is only one nonnoisy pixel in the new 7 × 7 window, then window size will be increased by 4 (11 × 11 for this case). This process is repeated for the entire image.
Minimum number of pixels required for computation
For high density noise removal, the corrupted image will have a very less number of nonnoisy pixels sparsely located. The anatomy of interpolation is to find a suitable value that is closer to original pixel, using the meager nonnoisy available in a current processing window. If more number of nonnoisy pixels is employed, then the impact of the pixel in restoring the original pixel diminishes, thereby leading to blurring. Hence, lesser number of pixels will have greater impact in the information preservation after removing the noise. If we have only two nonnoisy pixels on either side of an edge, then interpolated value lies between the two values (which will lead to blurring). Hence, opting for two nonnoisy pixels is not reliable for estimation. So a minimum of three nonnoisy pixels for the calculation of weights in IDWI is chosen.
Proposed algorithm
The inverse distance weighted interpolation filter is presented to eliminate high density salt and pepper noise. Initially obtain the nonnoisy pixels from the image by checking each pixel of the image if it is noisy or not. To start with, the size of the window is assumed to be 3 × 3. The size of the window increases depending upon the number of noisy pixel in the current processing window as shown in Fig. 1. This algorithm works on a basic assumption that at least three uncorrupted pixels are required to interpolate the new pixel. In a current processing window, if the uncorrupted pixels are less than three, then the size of the current processing window increases as shown in Fig. 1. Let K be the number of uncorrupted pixels in the current processing window as shown in Fig. 1.
Methodology of the proposed algorithm
Initially, the proposed algorithm detects the salt and pepper noise in an image. The pixels in the image are termed as noisy, if the values take either 0 or 255. If the pixel holds other than 0 or 255, then it is termed as nonnoisy pixels. The initial size of the proposed algorithm is 3 × 3. The algorithm scans for at most three nonnoisy pixels inside the current processing window failing which the size of the window increases. The methodology of the proposed algorithm is briefed in the section as follows in three different cases. The pixel to be processed is denoted with a black background. The variable “k” refers to the number of nonnoisy pixels in the current processing window. The variable v refers to the calculated weights from the inverse distance, and the variable g (i) gives the value of the nonnoisy pixels inside the current processing window. The new value is estimated based on the equations furnished in Eqs. 1 and 2.
Case A: Pixel was found to be noisy, and hence, case A has three possible operations for estimating the new pixel.
In this example, the pixel to be processed p (x, y) is 255. It is considered as noisy and hence either case 1, 2, or 3 will be performed based on the nonnoisy pixel in the current processing window.
Case 1: Initially, the window size is assumed to be 3 × 3. Find the number of nonnoisy pixels (k) in the current processing window. If k is greater than or equal to three, then find the weights from their positions (−1, 0, 1).
In this example, the number of nonnoisy pixels in the current processing window is k = 4. The algorithm requires at least three pixels for interpolation. Hence, interpolate the values using inverse distance weighted interpolation filter. Find the inverse distance of the nonnoisy pixels using its pixel position (−1, 0, 1) with respect to the position of the pixel to be processed. The pixel position of 122 in the first row is (1, −1). The inverse distance is calculated using these above said positions as calculated in d (1).

(i)
Find the inverse distance of nonnoisy pixel from the center pixel.

(ii)
Find the sum of the inverse distance.

(iii)
Divide all the inverse distance by sum which results in the weights of the nonnoisy pixels.

(iv)
Interpolate the values by multiplying the weights with the corresponding nonnoisy pixels in the current processing window. P (255) refers to the noisy pixel.
The noisy pixel is replaced by 122 as shown in Fig. 2.
Case 2: Initially, the window size is assumed to be 3 × 3 as shown in Fig. 3. Find the number of nonnoisy pixels (k) in the current processing window.
If k is equal to 2, then increase the window size by 2 to get the window size 5 × 5 and then find the weights from their positions of nonnoisy pixels. In this example, the number of nonnoisy pixels (k) in the current processing window is 2. So the window size was increased from 3 × 3 to 5 × 5 as shown in Fig. 4.
The expanded window consists of four numbers of nonnoisy pixels. Hence, estimate the processed pixel using inverse distance weighted interpolation filter as shown below. Find the inverse distance of the nonnoisy pixels using its pixel position (−2, −1, 0, 1, 2) with respect to the position of the pixel to be processed. The pixel position of 127 in the first row is (−2, −2). The inverse distance is calculated using these above said positions as calculated in d (1).

(i)
Find the inverse distance of nonnoisy pixel from the center pixel.

(ii)
Find the sum of the inverse distance.

(iii)
Divide all the inverse distance by sum which results in the weights of the nonnoisy pixels.

(iv)
Interpolate the values by multiplying the weights with the corresponding nonnoisy pixels in the current processing window. P (0) refers to the noisy pixel.
The corrupted pixel is replaced by 128 as shown in Fig. 4.
Case 3: Initially, the window size is assumed to be 3 × 3 as shown in Fig. 5. Find the number of nonnoisy pixels (k) in the current processing window.
If k is equal to 1, then increase the window size by 4 to get the window size 7 × 7 and then find the weights from their positions of the nonnoisy pixels. In this example, the number of nonnoisy pixels in the current processing window is k = 1. So the window size was increased from 3 × 3 to 7 × 7 as shown in Fig. 6.
The increased window will now consist of increased nonnoisy pixels. Now, the number of nonnoisy pixels increased to 7. Hence, interpolate the values using inverse distance weighted interpolation filter. Find the inverse distance of the nonnoisy pixels using its pixel position (−3, −2, −1, 0, 1, 2, 3) with respect to the position of the pixel to be processed. The pixel position of 129 in the second row is (1, −2). The inverse distance is calculated using these above said positions as calculated in d (1).

(i)
Find the inverse distance of nonnoisy pixel from the center pixel.

(i)
Find the sum of the inverse distance.

(ii)
Divide all the inverse distance by sum which results in the weights of the nonnoisy pixels.

(iii)
Interpolate the values by multiplying the weights with the corresponding nonnoisy pixels in the current processing window. P (0) refers to the noisy pixel.
The noisy pixel is replaced by 127 as shown in Fig. 6.
Case 4: Initially, the window size is assumed to be 3 × 3 as shown in Fig. 7.
Find the number of nonnoisy pixels (k) in the current processing window. If k is equal to 0, then increase the window size by 6 to get the window size 9 × 9 and then find the weights from their positions of nonnoisy pixels. In this example, the number of nonnoisy pixels in the current processing window is k = 0. So the window size was increased from 3 × 3 to 9 × 9.
The window size is increased to 7, resulting in increased number nonnoisy pixels to 4. Hence, estimate the new pixel values using inverse distance weighted interpolation filter. Find the inverse distance of the nonnoisy pixels using its pixel position (−4, −3, −2, −1, 0, 1, 2, 3, 4) with respect to the position of the pixel to be processed. The pixel position of 129 in the second row is (−3, −4). The inverse distance is calculated using these above said positions as calculated in d (1).

(i)
Find the inverse distance of nonnoisy pixel from the center pixel.

(ii)
Find the sum of the inverse distance.

(iii)
Divide all the inverse distance by sum which results in the weights of the nonnoisy pixels.

(iv)
Interpolate the values by multiplying the weights with the corresponding nonnoisy pixels in the current processing window. P (0) refers to the noisy pixel.
The noisy pixel is replaced by 124 as shown in Fig. 8.
Case B: Pixel was not found to be nonnoisy, and hence, the pixel is left unaltered.
In this case, the pixel to be processed is 173 which lie between 0 and 255. Hence, the pixel is termed as nonnoisy and left unaltered as shown in Fig. 9.
Simulation results and discussions
The inverse distance weighted interpolation filter is evaluated based on different metrics such as peak signaltonoise ratio (PSNR), image enhancement factor (IEF), and error rate which is given in Eqs. 3, 5, and 6, respectively. The metrics mean squared error (MSE) is given in Eq. 4. The structural similarity index metrics (SSIM) is calculated on various windows of an image. The measure between two windows x and y of common size MXN is given in Eq. 7.
Where r refers to original image, n gives the corrupted image, x denotes restored image, M × N is the size of processed image. C refers to number of restored pixels not equal with original pixels in restored and original image, respectively.
Where x refers to original image, y refers to restored image, μx is the average of x, μy is the average of y, σx standard deviation of x, σy is the standard deviation of y. C1 = (K _{1} L)^{2} and C2 = (K _{2} L)^{2} are two variables to stabilize the division with weak denominator; L is the dynamic range of the pixel values (for an 8bit image, it takes from 0 to 255), K1 = 0.01, and K2 = 0.03 by default [28]. The edge preservation performance of the proposed algorithm is evaluated using Pratt’s figure of merit (Pratt’s FOM) [29]. The figure of merit of Pratt, which calculates the alikeness between two edge images, is given in Eq. 8.
where KI and KB are the different points of edges in the restored image and original image, respectively, di is the distance between an edge pixel and the nearest edge pixel of the original and α is a constant and was used α = 1/9, optimal value established by Pratt in [29]. The algorithms used in this paper are derived from the references cited in the square brackets below. The existing algorithms used for the comparison are standard median filter of window size 3 (SMF (3 × 3) [30], adaptive median filter (AMF) (W _{max} = 39) [6], centerweighted median filter (CWF) [7], progressive switched median filter (PSMF) [10], modified decision based median filter (MDBMF) [16], alphatrimmed mean filter (ATMF) (trimming factor is 4) [30], decision based algorithm (DBA) [13], cascaded filters (CUTMF, CUDBMPF) [15], modified decision based unsymmetric trimmed median filter (MDBUTMF) [17], improved decision based median filter (IDBA) [14], noise adaptive fuzzy switching median (NAFSM) [31], modified decision based unsymmetrical trimmed median filter with global mean MDBUTMFGM [18], adaptive cardinal Bspline algorithm (ACBSA) [20], and adaptive weighted mean filter (AWMF) (W _{max} = 39) [30]. All the algorithms used in the paper were tested on Kodak natural image database hosted in University of Southern California and Signal and Image Processing Institute website [32]. The images used in the paper are based on the varying information content contained in it. Images such as Lena, cameraman, boat, peppers, and baboon were showcased in the paper as part of the illustration. Exhaustive experiments were conducted, and each algorithm were subjected to different images of the database. All the experiments of the quantitative analysis were done by varying the noise densities in the images from 10 to 90%. All the simulation was done in secondgeneration i32350 CPU with an operating frequency of 2.30 GHz with a 4GB RAM capability. Tables 1, 2, 3, 4, and 5 give the performance of various algorithms at different noise densities on Lena image corrupted by salt and pepper noise for peak signaltonoise ratio (PSNR), SSIM, IEF and Pratt’s FOM, respectively. The effectiveness of the algorithm in eliminating noise is decided by PSNR, IEF, and SSIM. The edge preservation capability is checked using Pratt’s FOM. The efficiency of the algorithm in correctly identifying noisy pixel from noise free is given by error rate. From Table 1, it is vivid that the algorithm exhibits excellent noise suppression capabilities by offering very good PSNR at high noise densities. Table 2 exhibits excellent structural preservation characteristics of the proposed algorithm when compared to other standard and existing algorithms. Table 3 shows a very low value for mean squared error for the proposed algorithm in comparison with other algorithms. Table 4 illustrates the edge preserving performance (Pratt’s FOM) of the proposed algorithm after the removal of salt and pepper noise. The Pratt’s FOM of the proposed algorithm is very high when compared with standard and existing algorithms. The proposed algorithm estimates the new pixel values based on the weights generated by inverse distance of the pixels in a small neighborhood with reference to noisy pixel. These weights or uncorrupted pixel of the current processing window has a strong influence on the new interpolated pixel. These uncorrupted neighbors preserve the originality of the corrupted pixel. This makes the proposed algorithm to have excellent quantitative results. From Tables 1, 2, and 3, we infer that algorithms such as SMF, PSMF, trimmed mean filter (TMF), and CWF fail to eliminate noise at high noise densities. Hence, for further analysis, the algorithms were not considered. Table 5 indicates the error rate of various weighted nonlinear filters while removing high density salt and pepper noise. Experiments were conducted to test the efficiency of various weighted algorithms in terms of error rate. It was found that the proposed algorithm induces less error when compared to the adaptive cardinal Bspline algorithm and adaptive weighted mean algorithm. The same trend continues even at high noise densities, making the proposed algorithm very good for salt and pepper noise removal. Table 6 gives the quantitative performance of the proposed algorithm on different images. Few additional measures were used to illustrate the performance of proposed algorithm such as cumulative probability of blur detection (CPBD) (no reference metrics) and normalized crosscorrelation (NCC) (similarity based metrics). The value of the CPBD indicates that even without the reference image (original image), the algorithm performs well. The NCC value indicates the similarity of the restored image with original image even after the removal of high density salt and pepper noise. Figure 10 shows the visual result of the proposed algorithm over the existing and standard algorithms on Lena image corrupted by 90% salt and pepper noise. Figure 11 briefs the qualitative performance of the proposed algorithm over the other algorithms on synthetic image corrupted by 90% salt and pepper noise. It was found from Fig. 10 that the proposed algorithm exhibits good visual quality when compared to other algorithms. Figure 11 is a synthetic image built using 21 different visual gray levels by a common eye. The information preservation (edges and fine details) of the images is quantified from various edge regions (line step and roof edge) of a synthetic image. The edge preservation of the proposed algorithm is explained by considering edges (information of the image) as line edge, ramp edge (semi constant region), step edge (constant region), and roof edge, respectively [22]. The proposed algorithm attenuates line, step, and roof edge, but preserves the semiconstant area (ramp edge). On careful observation from the synthetic images, it was found that the algorithms such as DBA, IDBA exhibits streaking, AMF attenuates edges completely, MDBUTMF also blurs the edges, and spline interpolation filter exhibits edge jittering. NAFSM does not remove salt and pepper noise completely. AWMF also eliminates step edge (first two regions of black and gray). The proposed algorithm showed better detail preserving capability than other algorithms because it attenuates the line edges (last two white regions) but preserves ramp edges (varying regions horizontally) even at high noise densities. At high noise densities, existing algorithms exhibits good noise suppression characteristics but induce blurring, smearing, streaking, and fading effect as seen in AMF, DBA, and MDBUTMF, respectively. The proposed algorithm gives excellent noise suppression capability without inducing any artifacts (blurring, smearing, streaking, and fading effect). Figure 12 gives the qualitative performance of various weighted algorithms corrupted by 90% salt and pepper noise on boat, pepper, and baboon image. It was found that proposed algorithm exhibits good noise removal capabilities. Existing weighted techniques such as AWMF attenuate edges at high noise densities. The maximum size of the window used in AMF and AWMF is 39. A larger window size blurs the image. In proposed algorithm, the window size is varied depending upon the nonnoisy content of the current processing window. Figure 13 gives the average PSNR and IEF of different algorithms on different images. Figure 14 illustrates the performance of the proposed algorithm on different images such as boat, Lena, and baboon for different images for the quantitative measure MSE and SSIM. Figure 15 gives the frequency spectrum of the different weighted methods on boat image corrupted by 90% salt and pepper noise. The proposed method is compared with the frequency response of ideal lowpass filter. An ideal lowpass filter should have a narrow central lobe. After filtering the restored image, it is subjected to discrete Fourier transform (DFT) and the frequency spectrum is plotted. It was found that for the proposed algorithm has a narrow central lobe in the spectrum (which is desired for any good filter). From the observation, other adaptive weighted methods such as AMF, NAFSM, and AWMF have a broader central lobe (in comparison to proposed algorithm). This makes the proposed weighted interpolation better than other weighted methods as their frequency response is closer to ideal lowpass filter.
Conclusion
A decision based inverse distance weighted interpolation filter for the elimination of high density salt and pepper noise in images is proposed. The proposed algorithm uses the inverse distance from nonnoisy pixels inside the current processing window to estimate the new pixel. This method requires at least three pixels for estimation of new pixels, failing which method uses an adaptive window based on the noisy content of the current processing window. The proposed algorithm shows a very good performance due to the usage of inverse distance of nonnoisy pixel in the current processing window for estimating new pixels. The interpolation technique is applied to an image only if the processed pixel is noisy. The proposed filter was applied on different images and found to show very high PSNR, SSIM, and IEF with good computational time supporting with excellent noise removal capability. The information preserving capability of the proposed algorithm is justified with an excellent Pratt’s FOM value even at high noise densities. The low error rate indicates the efficiency of the proposed algorithm in restoring the original image. The central lobe in the frequency spectrum of the restored image is narrow, thereby making the proposed algorithm suitable for noise removal. The proposed algorithm eliminates standard and existing algorithms in terms of excellent noise suppression and good detail preservation without leading to any ambiguity at very high noise densities.
References
 1.
M. Hanumantharaju, M. Ravishankar, D.R. Rameshbabu, S. Ramachandran, Color Image Enhancement Using Multiscale Retinex with Modified Color Restoration Technique, Second International Conference Proceedings on Emerging Applications of Information Technology (2011), p. 93–97
 2.
Z. Min, W. Jiechao, L. Zhiwei, L. Yonghua, An adaptive image zooming method with edge enhancement, 3rd International Conference Proceedings on Advanced Computer Theory and Engineering (2010), p. 608–611
 3.
W.S. Tam, C.W. Kok, W.C. Siu, A modified edge directed interpolation for images. Elec. Imag. J. 19(1), 1–20 (2010)
 4.
W. Bender, C. Rosenberg, Image enhancement using nonuniform sampling. Proc. SPIEInt. Soc. Opt. Eng. 1460, 59–70 (1991)
 5.
T.S. Huang, G.J. Yang, Fast two dimensional median filtering algorithm. IEEE Trans. Acoust. Speech Signal Process 27, 13–18 (1979)
 6.
H. Hwang, R.A. Haddad, Adaptive median filters: new algorithms and results. IEEE Trans. Image Process 4, 499–502 (1995)
 7.
S.J. KO, Y.H. Lee, Center weighted median filters and their application to image enhancement. IEEE Trans. Circuits Syst. 38(9), 984–993 (1991)
 8.
D.R.K. Brownrigg, The weighted median filter. Commun. ACM 27(8), 807–818 (1984)
 9.
T. Chen, K. Ma, Tristate median filter for image denoising. IEEE Trans. Image Process 8(12), 1834–1837 (1999)
 10.
Z. Wang, D. Zhang, Progressive switching median filter for removal of impulse noise from highly corrupted images. IEEE Trans. Circuits Syst. II 46, 78–80 (1999)
 11.
H. Eng, K. Ma, Noise adaptive softswitching median filter. IEEE Trans. Image Process 10(2), 242–251 (2001)
 12.
H. Raymond, C.W.H. Chan, M. Nikolova, Salt and pepper noise removal by median—type noise detectors and detail preserving regularization. IEEE Trans. Image Process 14(10), 1479–1485 (2005)
 13.
K.S. Srinivasan, D. Ebenezer, A new fast and efficient decisionbased algorithm for removal of highdensity impulse noises. IEEE Signal Process Lett. 14(3), 189–192 (2004)
 14.
M.S. Nair, K. Revathy, R Tatavarti, An Imroved Decision Based Algorithm for Impulse Noise Removal, Proceedings in Congress on Image and Signal Processing (2008), p. 426–431
 15.
S. Balasubramanian, S. Kalishwaran, R. Muthuraj, D. Ebenezer, V. Jayaraj, An efficient non linear cascade filtering algorithm for removal of high density salt and pepper noise in image and video sequence, Proceedings in International Conference on control, Automation, communication and Energy Conservation (2009), p. 1–6
 16.
Aiswarya K., Jayaraj V., Ebenezer D.A. New and Efficient Algorithm for the Removal of High Density Salt and Pepper Noise in Images and Videos, Proceedings in International conference on computer modelling and simulation (2010), p. 409–413
 17.
S. Esakkirajan, T. Veerakumar, A.N. Subramanyam, C.H. Prem Chand, Removal of high density salt and pepper noise through modified decision based unsymmetrial trimmed median filter. IEEE Signal Process Lett. 18(5), 287–290 (2011)
 18.
T. Veerakumar, S. Esakkirajan, I. Vennila, An approach to minimize very high density salt and pepper noise through trimmed global mean. Int. J. Compt. Appl. 39(12), 29–33 (2012)
 19.
K. Vasanth, V.J. Senthilkumar, A decision based asymmetricsal trimmed mid point algorithm for the removal of high density salt and pepper noise. App. Theo. Info Tech. J. 42(2), 553–563 (2012)
 20.
P. Syamala Jayasree, P. Raj, P. Kumar, R. Siddavatam, S.P. Ghrera, A fast novel algorithm for salt and pepper image noise cancellation using cardinal Bsplines. Int. J. Signal Image Video Process Springer 7(6), 1145–1157 (2012)
 21.
T. Veerakumar, S. Esakkirajan, I. Vennila, Recursive cubic spline interpolation filter approach for the removal of high density saltandpepper noise. Int. J. Signal Image Video Process Springer 8(1), 159–168 (2013)
 22.
K. Vasanth, V.J.S. Kuma, A decision based neighborhood referred unsymmetricsal trimmed variants filter for the removal of high density salt and pepper noise in images and videos. Int. J. Signal Image Video Process Springer 9(8), 1833–1841 (2014)
 23.
P. Zhang, F. Li, A new adaptive weighted mean filter for removing saltandpepper noise. IEEE Signal Process Lett. 21(10), 1280–1283 (2014)
 24.
D.J. Brus, J.J. De Gruijter, B.A. Marsman, R. Visschers, A.K. Bregt, A. Breeuwsma, J. Bouma, The performance of spatial interpolation methods and choropleth maps to estimate properties at points: a soil survey case study. Environmetrics 7(1), 1–16 (1996)
 25.
E. Hosseini, J. Gallichand, J. Caron, Comparison of several interpolators for smoothing hydraulic conductivity data in South West Iran. Trans. Am. Soc. Agric. Eng. 36(6), 1687–1693 (1993)
 26.
E.H. Isaaks, R.M. Srivastava, Applied Geostatistics (Oxford University Press, New York, 1989)
 27.
R. Webster, M. Oliver, Geostatistics for Environmental Scientists (Wiley, Chichester, 2001)
 28.
Z. Wang, A.C. Bovik, H.R. Sheikh, E.P. Simoncelli, Image quality assessment: from error measurement to structural similarity. IEEE. Trans. Image Process 13(4), 102–109 (2004)
 29.
I.A. Abdou, W. Pratt, Quantitative design and evaluation of enhancement/thresholding edge detectors. Proc. IEEE 67(5), 753–766 (1979)
 30.
J. Astola, P. Kusmanen, Fundamentals of non linear digital filtering (CRC press,Boca Raton, NewYork 1997)
 31.
K.K.V. Toh, N.A.M. Isa, Noise adaptive fuzzy switching median filter for saltandpepper noise reduction. IEEE Signal Process Lett. 17(3), 281–284 (2010)
 32.
SIPI Image Database, USE Signal and Image Processing Institute, http://sipi.usc.edu/database/. Accessed 15 May 2015
Acknowledgements
The authors would like to thank the Management of Vidya Jyothi Institute of Technology, Aziz nagar, Chilukur Road, Hyderabad75, Telengana, India for providing us with the infrastructure to complete the research work.
Funding
The research is done by the author himself, and there is no role of funding body involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
Data will not be shared, reason for not sharing the data and materials is that the work submitted for review is not completed. The research is still ongoing, and those data and materials are still required by the author and coauthors for further investigations.
Author information
Affiliations
Contributions
VK is involved in concept, design, acquisition of data, and initial and final analysis. GP is involved in detailed analysis of the code and manuscript drafting (after review). HK is involved in detailed interpretation and manuscript drafting (after revision). VS is involved in code development and initial analysis. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Authors’ information
Vasanth Kishore Babu received Electronics and Communication from Arulmigu Kalasalingam College of Engineering affiliated to Madurai Kamaraj University with firstclass distinction in 2003. He later got M.E degree in Applied Electronics from Sathyabama University in the year 2007 securing university Fifth rank. He completed his PhD from the same university in the year 2014. He acted as the project leader with the prestigious Sathyabama Nano Satellite for the 25yearold Sathyabama University. He currently works as a professor of the Department of Electronics and Communication Engineering at Vidya Jyothi Institute of Technology, Hyderabad, India. His current research interests include nonlinear filtering, image enhancement, VLSI signal processing for image processing algorithms, and spacequalified systems. He has published more than 50 research papers in reputed journals and conference proceedings.
Ganesan Packyanathan received his BE degree in Electronics and Communication Engineering from Madurai Kamaraj University, India, in 1998 and ME in Electronics and Control Engineering and PhD degree in Image Processing from Sathyabama University, India, in 2007 and 2015, respectively. Currently, he is a professor of the Department of Electronics and Communication Engineering at Vidya Jyothi Institute of Technology, Hyderabad, India. His research interests include image processing, soft computing, computer vision, and medical image analysis. He has published more than 50 papers in various journals and refereed conferences.
Harikrishna Kamatham obtained his B.Tech in Electronics and Communication Engineering from Jawaharlal Nehru Technological University Hyderabad and later moved to USA for his Master of Science degree in Electrical and Computer Engineering from Southern Illinois University Carbondale, IL, USA. He has completed his PhD degree in Telecommunication Engineering from SRM University, Chennai in 2013. His areas of interest include embedded systems, VLSI, and wireless communications. He is currently working in Vidya Jyothi Institute of Technology as a professor and head of ECE. He has published more than 20 papers in National and International Conferences and journals.
Vishnu Shankar is currently working as an analyst at Verizon Data Services, Chennai India. His research field includes computer vision for robots, nonlinear filtering, and mathematics for computer vision. He has won several laurels at Maths Olympiad at IIT.
Ethics approval and consent to participate
Not applicable.
Consent for publications
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Kishorebabu, V., Packyanathan, G., Kamatham, H. et al. An adaptive decision based interpolation scheme for the removal of high density salt and pepper noise in images. J Image Video Proc. 2017, 67 (2017). https://doi.org/10.1186/s1364001702150
Received:
Accepted:
Published:
Keywords
 Salt and pepper noise
 Interpolation technique
 Inverse distance weighted interpolation
 Edge preservation