Image Denoising Using Non Linear Filter: Mr. Vijay R. Tripathi
Image Denoising Using Non Linear Filter: Mr. Vijay R. Tripathi
Image Denoising Using Non Linear Filter: Mr. Vijay R. Tripathi
<
Where T is a predefined threshold value. Once a pixel (i, j)
is detected as an impulse, the value of
( )
( , )
n
i j
x is subsequently modified
( 1) ( ) ( 1)
( , ) ( , ) ( , )
( )
( , )
( 1) ( ) ( 1)
( , ) ( , ) ( , )
,
,
n n n
i j i j i j
n
i j
n n n
i j i j i j
m if f f
x
x if f f
=
=
=
(4.3)
Suppose the impulse detection procedure is stopped after
the N
D
th iteration, then two output images-
( )
( , )
D
N
i j
x and
( )
( , )
D
N
i j
f are obtained, but only
( )
( , )
D
N
i j
f is useful for our noise
filtering algorithm.
2. Noise Filtering
Like the impulse detection procedure, the noise
filtering procedure also generates a gray scale image
sequence,
(0) (1) (2) ( )
( , ) ( , ) ( , ) ( , )
{ , , ,.... .....}
n
i j i j i j i j
y y y y and a binary
flag image sequence
(0) (1) ( )
( , ) ( , ) ( , )
{ , ,.... .....}
n
i j i j i j
g g g . In the
gray scale image sequence, we still use
(0)
( , ) i j
y to denote the
pixel value at position (i, j) in the noisy image to be filtered
and use
( )
( , )
n
i j
y to represent the pixel value at position (i, j)
in the image after the nth iteration. In a binary flag
image
( )
( , )
n
i j
g , the value
( )
( , )
n
i j
g =0 means the pixel (i, j) is
good and
( )
( , )
n
i j
g = 1 means it is an impulse that should be
filtered. A difference between the impulse detection and
noise-filtering procedures is that the initial flag image
(0)
( , ) i j
g of the noise-filtering procedure is not a blank image,
but the impulse detection result
( )
( , )
D
N
i j
f , i.e.,
(0)
( , ) i j
g =
( )
( , )
D
N
i j
f .
In the nth iteration (n = 1; 2; .), for each
pixel
( 1)
( , )
n
i j
y
, we also first find its median value
( 1)
( , )
n
i j
m
of a
W
F
W
F
(W
F
is an odd integer and not smaller than 3)
window centered about it. However, unlike that in the
impulse detection procedure, the median value here is
selected from only good pixels with
( 1)
( , )
n
i j
g
= 0 in the
window.
Let M denote the number of all the pixels with
( 1)
( , )
n
i j
g
= 0 in the W
F
W
F
window. If M is odd, then
( 1)
( , )
n
i j
m
=
( 1) ( 1)
( , ) ( , )
{ 0, ( , ) }
n n
i j i j F F
median y g i j W W
= e
The value of
( )
( , )
n
i j
y is modified only when the pixel (i, j) is
an impulse and M is greater than 0:
( 1) ( 1)
( , ) ( , )
( )
( , )
( 1)
( , )
, 1; 0
,
n n
i j i j
n
i j
n
i j
m if g M
y
y else
= >
Once an impulse pixel is modified, it is considered as a
good pixel in the subsequent iterations
( 1) ( ) ( 1)
( , ) ( , ) ( , )
( )
( , )
( ) ( 1)
( , ) ( , )
,
0,
n n n
i j i j i j
n
i j
n n
i j i j
g if y y
g
if y m
=
=
The procedure stops after the N
F
th iteration when all of the
impulse pixels have been modified, i.e.,
( , )
( , )
F
N
i j
i j
g
=0
Then we obtain the image
( )
( , )
{ }
F
N
i j
y which is our restored
output image.
Table 1.WF(3 * 3) Median Filter,CWM Filter&PWM Filter
Noise PSNR
With
Noise
PSNR After Filtering
Median
Filter
CWM
Filter
PWM
Filter
10%
15 12 13.5 29
20%
12 12 13 27
30%
11 12 13 26
40%
9.5 12 13 24
International Journal of Modern Engineering Research (IJMER)
www.ijmer.com Vol.2, Issue.6, Nov-Dec. 2012 pp-4543-4546 ISSN: 2249-6645
www.ijmer.com 4545 | Page
The Experimental results are shown in
Table.2 for WF (5 * 5) PWM Filter
Table 2. WF (5 * 5) PWM Filter
Noise
PSNR with
Noisy
PSNR after
Filtering
10 15 21
20 12 19
30 11 19
40 9.5 19
The Experimental results are shown in
Table.3 for WF (7 * 7) PWM Filter
Table 3. WF (7 * 7) PWM Filter
Noise
PSNR with
Noisy
PSNR after
Filtering
10 15 14
20 12 7. 3
30 11 6
40 9.5 5.9
.
(a)
(b)
(c)
(d)
(e)
<f>
(g)
(h)
<i>
International Journal of Modern Engineering Research (IJMER)
www.ijmer.com Vol.2, Issue.6, Nov-Dec. 2012 pp-4543-4546 ISSN: 2249-6645
www.ijmer.com 4546 | Page
Figure 2. a, b, c, d are noisy images of Lena (512512)
corrupted by salt and pepper noise with noise density of
10%, 20%, 30%, 40% respectively and corresponding
restored image by PWM are in e, f, g, h for WF (3 * 3)
III. Conclusion
In this entire dissertation work, two different non-linear
filters are implemented and extensive experiments are
performed to obtain the results with various parameters to
assess the performance of each filter. The plot of PSNR for
these two filters is given below. The Table.1, 2 &3 below
shows the PSNR value obtain using Lena Image of size 512
x 512.
From the PSNR value mention in the simulation result it is
very clear that PWM Filter shows better performance in
suppressing impulsive noise compare to above mention
filter in suppressing impulse noise when noise exceeds from
10 % to 40 %.
& Secondly extensive experimental result show that if we
increase window size i.e. 5 * 5 & 7 * 7,we find that by
increasing the window size image get more & more
corrupted & filter is not able to suppress impulsive noise
effectively compare to when window size in filter was (3 *
3) in filter. Though simulation time required is less, which
is given in table below.
Table4: Average Run Time in Sec
Therefore from the above table it is very cleared
that as window size increases image get more & more
blurred & distorted though it requires less simulation time
compare to that when window size in filter was (3 * 3). So
mostly window size of WF (3 * 3) is preferred compare to
that of WF (5 * 5) & WF (7 * 7).
Reference
[1] R.C. Gonzalez and R.E. Woods Digital Image
Processing Second Edition.
[2] A. C. Bovik, T. S. Huang, and D. C. Munson, A
generalization of median filtering using Linear
combinations of order statistics, IEEE Tran Acoust.,
Speech, Signal Process., volASSP-31, no.6, pp.1342
1350, Dec. 1983.
[3] J. Astola and P. Kuosmanen, Fundamental of
Nonlinear Filtering, Boca Raton, CRC Press, 1997.
[4] T. Sun and Y. Neuvo, Detail-preserving median
based filters in image processing, Pattern Recognit.
Lett., vol. 15, pp. 341347, Apr. 1994
[5] I. Pitas and A. N. Venetsanopoulos, Nonlinear Digital
Filters: Principles and Applications, Boston, Kluwer,
1990.
[6] B. I. Justus son, Median filtering: Statistical
properties, Two-Dimensional Digital Signal
Processing II, T. S. Huang Ed., New York;
Springer Verlag, 1981.
[7] S-J. KO and Y. H. Lee, Center-weighted
median filters and their applications to image
enhancement, IEEE Trans. Circuits and Syst.,
vol. 38, pp. 984-993, Sept. 1991.
[8] J.-H. Wang, Rescanned minmax center-
weighted filters for image restoration, Proc.
Inst. Elect. Eng., vol. 146, no. 2, pp. 101107,
1999.
[9] Tri-State Median Filter for Image Image Denoising,
T Chen, K-K.Ma, and L-H Chen IEEE
TRANSACTIONS ON IMAGE PROCESSING,
VOL. 8, NO. 12, DECEMBER 1999.
[10] T. Kasparis, N. S. Tzannes, and Q., Detail-
preserving adaptive conditional median filters, J.
Electron. Imag, vol. 1, no. 14, pp. 358364, 1992.
Filter Used PSM
Window Size Time In
Sec
WF (3 * 3) 16sec
WF (5 * 5) 13sec
WF (7 * 7)
10sec