Advanced
Enhancing Medical Images by New Fuzzy Membership Function Median Based Noise Detection and Filtering Technique
Enhancing Medical Images by New Fuzzy Membership Function Median Based Noise Detection and Filtering Technique
Journal of Electrical Engineering and Technology. 2015. Sep, 10(5): 2197-2204
Copyright © 2015, The Korean Institute of Electrical Engineers
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/ licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : August 06, 2014
  • Accepted : June 24, 2015
  • Published : September 01, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
G. Elaiyaraja
Corresponding Author: Dept. of Electronics and Communication Engineering, Adhiparasakthi College of Engineering, G.B.Nagar, Tamilnadu, India. (elaiyarg@gmail.com)
N. Kumaratharan
Dept. of Information and Technology, Srivenkateswara College of Engineering, Sriperumbudur, Tamilnadu, India. (kumaratharan@ rediffmail.com)

Abstract
In recent years, medical image diagnosis has growing significant momentous in the medicinal field. Brain and lung image of patient are distorted with salt and pepper noise is caused by moving the head and chest during scanning process of patients. Reconstruction of these images is a most significant field of diagnostic evaluation and is produced clearly through techniques such as linear or non-linear filtering. However, restored images are produced with smaller amount of noise reduction in the presence of huge magnitude of salt and pepper noises. To eliminate the high density of salt and pepper noises from the reproduction of images, a new efficient fuzzy based median filtering algorithm with a moderate elapsed time is proposed in this paper. Reproduction image results show enhanced performance for the proposed algorithm over other available noise reduction filtering techniques in terms of peak signal -to -noise ratio (PSNR), mean square error (MSE), root mean square error (RMSE), mean absolute error (MAE), image enhancement factor (IMF) and structural similarity (SSIM) value when tested on different medical images like magnetic resonance imaging (MRI) and computer tomography (CT) scan brain image and CT scan lung image. The introduced algorithm is switching filter that recognize the noise pixels and then corrects them by using median filter with fuzzy two-sided π- membership function for extracting the local information.
Keywords
1. Introduction
Medical images are distorted by salt and pepper noise is largely caused by the acquisition or transmission of digital images in a noisy channel [13] . Salt and pepper noise is one type of impulse noise which can degrade the image visualization, where the noisy pixels can take only the minimum and maximum pixel elements of image. The linear filter like mean filter and interrelated filters are not efficient for removal of salt and pepper noise. Non-linear filtering techniques like standard median filter (SMF) and Adaptive median filter (AMF) are widely used to remove salt and pepper noise [21] . However, the SMF, which restores each pixel with median pixel value in the filtered image, exhibits blurring of reconstructed image [3] . SMF is effective only at low noise densities. Several techniques have been implemented for removal of salt and pepper noise at high noise densities [1 - 9] .
Recent methods like efficient detail-preserving approach (EDPA) and efficient edge-preserving algorithm (EEPA) based on alpha-trimmed mean statistical estimator have been proposed for the removal of salt-and-pepper noise without degrading image fine information’s [2] . Then, adaptive impulse detector with center-weighted median (ACWM) filter efficiently to remove salt-and-pepper noise. These techniques simply perform well when an image is degraded with 50% of noise densities and small window size. On the other hand, Progressive Switching Median Filter(PSMF), decision based algorithm (DBA), modified decision based algorithm (MDBA), are one of the fastest and efficient algorithms of capable of salt-andpepper noise removal at high noise densities as 80% [18] . A most important disadvantage of this algorithm is striping problem at high noise densities. To overcome this drawback, modified decision based unsymmetric trimmed median filter (MDBUTMF) is used to remove salt-and-pepper noise at very high densities as 80-90%. In response to this noise problem, the switching median filter [6] , in which an impulse detector of employed to classify the center pixel of the filtering window is proposed. In the filtering algorithm, at high noise density, the processing pixel is replaced by the mean value of pixel elements within the window [5] . This causes blurring of fine image details or high computational time [10] . Examining all those problems exist in literature, we have proposed a new switching based algorithm, which aids to accurate medical image visualization. This algorithm gives better metric values than the existing median based algorithms [2 - 6] .
2. Proposed Algorithm
The modified fuzzy switching median based filter is a non-causal two-stage filter, one stage for finding the pixels are noise-free or noise pixels and second stage for noise reduction as shown in figure Fig. 1 . The first step to detect the two salt-and-pepper noise intensities. When a “corrupted pixel” is discovered, it is switched to the next filtering stage [11] . Otherwise, when a pixel is marked as “uncorrupted pixel” it will be left unchanged and the filtering process is spared to avoid replacing any fine details and textures in the noisy image [4] .
PPT Slide
Lager Image
Proposed medical image enhancing structure
- 2.1 Switching noise detection process
The filter initially estimates the two salt-and-pepper noise intensities of histogram of noise image shown in Fig. 2 , and then the detection stage starts by searching for the two local maximums, Lsalt=255 and Lpepper=0, representing the two salt-and-pepper noise intensities starting from both ends of noisy image histogram [17] . These two salt-and-pepper noise pixels will be used to recognize possible “corrupted pixels” in the image as shown in Fig. 3 . A binary noise values B ( m , n ) will be formed based on noise pixels as shown Fig. 4 .
PPT Slide
Lager Image
Histogram of Lung CT image distorted with 50% of Salt-and-pepper noise
PPT Slide
Lager Image
(7x7) Noisy image pixel values
PPT Slide
Lager Image
(7x7) Binary mask value of noisy image
PPT Slide
Lager Image
where X ( m , n ) is the pixel at location ( m , n ) with intensity X , B ( m , n )=1 indicates “uncorrupted pixels” to be obtained from the distorted image and B ( m , n )=0 indicates “corrupted pixels”.
- 2.2 Fuzzy noise cancellation
In this stage after creation of binary noise mask value, “corrupted pixels” marked with B(m, n)=0 is replaced by an estimated adjustment pixel value and then increase the window size W2s+1(m, n) with (2s+1)x(2s+1) of chosen noise image. Then, the number of “uncorrupted pixels” in the filtering window W2s+1(m, n) is calculated using
PPT Slide
Lager Image
where S value is choose the window size of image
After counting of uncorrupted pixel, the minimum number of “uncorrupted pixel” is less than one then expands the filtering window by one pixel [15] . This process is continued until the condition G 2s+1 (m, n) greater than or equal to one [18] . The uncorrupted pixels select the median pixel M(m, n), given by
PPT Slide
Lager Image
The corrupted pixels select the median pixel M(m, n), given by
PPT Slide
Lager Image
After finding the value of median pixel M(m, n), the local information in a 3 × 3 window is extracted by computing the absolute pixel difference α(m, n), given by
PPT Slide
Lager Image
Then the local information is described by the maximum absolute pixel difference in the 3×3 filtering window
PPT Slide
Lager Image
In filtering stage, the Fuzzy reasoning is applied to the extracted local information d(m, n). The adopted fuzzy set, which is defined by the Two-sided π-membership function(MF) f(m, n) with four predefined thresholds T1(5), T2(10), T3(245) and T4(250) as shown in Fig. 5 [19] .
PPT Slide
Lager Image
Two-sided π-membership function
The two-sided π-membership function [19] , given by
PPT Slide
Lager Image
Where S(d(m,n),T1,T2) is S-shaped open-right MF as shown in Fig.6 , given by
PPT Slide
Lager Image
S-Shaped membership function
PPT Slide
Lager Image
and Z(d(m, n), T3, T4) is Z-shaped open-left MF as shown in Fig. 7 , given by
PPT Slide
Lager Image
Z-Shaped membership function
PPT Slide
Lager Image
Where the local information d(m, n) is used as the fuzzy variable and the four thresholds T1, T2, T3 and T4 are set to 5,10,245 and 250 respectively for best performance [11] . The restored image Y(m, n) is acquired by linear combination of weighted processing pixel X(m,n) and median pixel M(m, n) of noisy image [15] .
PPT Slide
Lager Image
3. Simulation Results and Discussion
In this section, the performance of the proposed fuzzy based median filter was tested with medical images with noise levels are varied from 10% to 90% and compared with other existing salt-and–pepper noise removal filters. The restoration performances are quantitatively measured by several different kinds of metric value such as peak signal to noise ratio (PSNR), mean square error (MSE), RMSE, image enhancement factor (IMF), mean absolute error (MAE) and elapsed time [18] .
The PSNR is most commonly used as a quality measurement for lossy compressed images. The PSNR is the ratio of the maximal power of original image and the noise power of distorted image. It is represented in the logarithmic domain because the powers of signals are usually in a wide dynamic range [6] . Its formula is given by
PPT Slide
Lager Image
MSE is the most widely used image quality assessment based on error sensitivity [10] . It is computed by taking the average of squared intensity differences in every pixel of a reference image and a distorted image:
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
where I(m, n) is Original image, Y(m, n) is Denoised image X (m, n)-Noise image and (M x N)-size of the image.
To obtain more correct assessments, the structural simailarity (SSIM) metric which represents perceptual image quality based on the structural information [13] . SSIM is an objective image quality metric and is superior to traditional quantitative measures such as MSE and PSNR [20] .
A general form of SSIM is
PPT Slide
Lager Image
Where α>0, β>0 and γ >0 are parameters used to adjust the relative importance of the three components and m, n are image patches.
Where l(m,n) is luminance comparison
c(m, n) is contrast comparison
and s(m,n) is structural comparison
PPT Slide
Lager Image
Where C 1 , C2, C3 are constants and μ m , μ n m , σ n , σ mn are local sample means and standard deviations of m and n.
The MSE, RMSE, MAE, PSNR, IEF, elapsed time in seconds and SSIM value are calculated for the proposed algorithm by varying the noise density from 10% to 90% and comparisons of performance values of various existing median based filters for lung CT scan image are tabulated in Tables 1 - 7 .
From the Table 1 , shows that the MSE value of proposed algorithm is very low compared to other existing noise removal methods even noisy densities varied from 10% to 90%. The RMSE value in Table 2 shows that at high noise levels, RMSE values of existing noise removal filters are very high as compared to the proposed method. From Table 3 , we illustrate that the MAE value of proposed method is very low compared to other existing noise removal method.
MSE of Lung CT scan image for various filters at different noise densities
PPT Slide
Lager Image
MSE of Lung CT scan image for various filters at different noise densities
RMSE of Lung CT scan image for various filters at different noise densities
PPT Slide
Lager Image
RMSE of Lung CT scan image for various filters at different noise densities
MAE of Lung CT scan image for various filters at different noise densities
PPT Slide
Lager Image
MAE of Lung CT scan image for various filters at different noise densities
The PSNR, IEF and SSIM values of proposed method in Table 4 , Table 5 and Table 7 shows that at higher noise levels the PSNR, IEF and SSIM values of existing methods are very low as compared to the proposed method. From all the performance values of proposed method, it is clearly observed that the performance of the proposed algorithm is moderately significant than the existing algorithms at both low and high noise densities. The performance analysis against noise densities for lung CT scan image is also done with help of plot shown in Figs. 8 - 13 .
PSNR of Lung CT scan image for various filters at different noise densities
PPT Slide
Lager Image
PSNR of Lung CT scan image for various filters at different noise densities
IEF of Lung CT scan image for various filters at different noise densities
PPT Slide
Lager Image
IEF of Lung CT scan image for various filters at different noise densities
Elapsed time of Lung CT scan image for various filters at different noise densities
PPT Slide
Lager Image
Elapsed time of Lung CT scan image for various filters at different noise densities
SSIM of Lung CT scan image for various filters at different noise densities
PPT Slide
Lager Image
SSIM of Lung CT scan image for various filters at different noise densities
PPT Slide
Lager Image
Comparison of MSE value of different algorithms for Lung CT scan image
PPT Slide
Lager Image
Comparison of RMSE value of different algorithms for Lung CT scan image
PPT Slide
Lager Image
Comparison of MAE value of different algorithms for Lung CT scan image
PPT Slide
Lager Image
Comparison of PSNR value of different algorithms for Lung CT scan image
PPT Slide
Lager Image
Comparison of IEF value of different algorithms for Lung CT scan image
PPT Slide
Lager Image
Comparison of SSIM value of different algorithms for Lung CT scan image
The difference in detail preservation performance of the proposed method may be better observed by looking at the output images of the different methods [19]. For this purpose, the Lung CT scan image shown in Fig. 14(a) of corrupted by 80% noise and noise image shown in Fig. 14(b) is restored by using different median based filters. Fig. 14(c) - 14(h) shows the output images of the different methods. It is observed from this figure that performance of the SMF and AMF are very close to each other. Some noise mark are easily visible in the output images of these two filters. The output images of the DBA, PSMF and UTDBMF are almost indistinguishable from each other and they are significantly better than those of the SMF and AMF. They show very good noise removal performance but considerably blur the small details of the image. It is observed that the proposed algorithm yields much better preservation performance [15] .
PPT Slide
Lager Image
(a) Original Lung CT scan image; (b) Corrupted image with 70% of salt-and-pepper noise; (c) SMF image; (d) AMF image; (e) DBA image; (f) PSMF image; (g) UTDBMF image and (h) Proposed filtered image.
The PSNR value of proposed filter is high compared to existing methods as shown in Fig. 10 .
Fig. 15 show the reconstructed results of MRI scan brain corrupted with 80% of noise of various filters [17] . From the restored results the proposed filter produce excellent quality image for high noise density. From the Table 6 and Fig. 16 , we observe that the proposed filter provides better trade off among elapsed time.
PPT Slide
Lager Image
(a) Original MRI scan brain image; (b) 80% noise image; (c) Proposed filtered image; (d) Original CT scan brain image; (e) 80% noise image (f) Proposed filtered image
PPT Slide
Lager Image
Comparison of elapsed time value of different algorithms for Lung CT scan image
4. Conclusion
In this paper, a proposed fuzzy based median filtered algorithm has been implemented for enhancing the different medical images in terms of removal of salt-and-pepper noise from highly corrupted noisy image. The extracted local information for optimal performance is achieved by two-sided π-MF compare with other fuzzy membership function. Results show that the proposed algorithm improves the metric values in comparison with other existing algorithms in terms of higher PSNR, IEF and SSIM for noise densities from 10% to 90%. In addition, the moderately elapsed time and easy implementation for image denoising methods.
BIO
G. Elaiyaraja is currently working as an Assistant Professor in the Department of Electronics and Communication Engineering, Adhiparasakthi College of Engineering, G.B. Nagar, Vellore Tamilnadu, India since 2007.He received his bachelor degree from university of madras, Adhiyamaan of College of Engineering, Hosur in 1999 and Master degree from college of Engineering, Guindy, Anna University Chennai in 2004. He is currently doing his Ph.D under Anna University, Chennai, Tamilnadu, India. He has 14 years teaching experience. His research interest is digital image processing, nonlinear filters, and digital signal processing.
N.Kumaratharan is currently working as a Professor in the Department of Information and Technology, Sriperumbudur, Tamilnadu, India. He received his bachelor degree from University of Madras in 2001 and Master degree from college of Engineering, Guindy, in 2004. He received his Ph.D Degree in Electronics and Communication Engineering from Pondicherry University in 2010. His research interest includes wireless sensor networks, wireless and broadband communication, digital signal processing and digital image processing. He has published several papers in international journals and international Conferences. He is also member of review committees for many national and international journals.
References
wang Z. 1997 ” Impulse noise detection and removal using fuzzy Techniques,” ELEECTRONICS LETTER 33 (5) 378 - 379    DOI : 10.1049/el:19970257
Esakkirajan S. , Veerakumar T. , Subramanyan A. N. , PremChand C. H. 2011 ” Removal of High Density Salt and Pepper Noise Through Modified Decision Based Unsymmetric Trimmed Median Filter,” IEEE Signal Processing Letters 18 (5) 287 - 290    DOI : 10.1109/LSP.2011.2122333
Gonzalez R.C. , Woods R.E. 2002 Digital Image Processing Printice-Hall
Farbiz F. , Menhaj M. B. 2000 “Fuzzy Techniques in Image Processing,” Springer-Verlag New York
Vin Toh Kenny Kal , Mat Isa Nor Ashidi 2010 “Noise adaptive fuzzy switching median filter for salt-and pepper noise reduction,” IEEE signal processing letters 17 (3) 281 - 284    DOI : 10.1109/LSP.2009.2038769
Toh K. K. V. , Ibrahim H. , Mahyuddin M. N. 2008 “Salt-and-pepper noise detection and reduction using fuzzy switching median filter,” IEEE Trans. Consumer Electron. 54 (4) 1956 - 1961    DOI : 10.1109/TCE.2008.4711258
Toh Kenny Kal Vin , Mat Isa Nor Ashidi 2010 “Cluster-Based Adaptive Fuzzy Switching Median Filter for Universal Impulse Noise Reduction,” IEEE Transaction on consumer electronics 56 (4) 2560 - 2568    DOI : 10.1109/TCE.2010.5681141
Ng P. E. , Ma K.-K. 2006 “A switching median filter with boundary discriminative noise detection for extremely corrupted images,” IEEE Trans. Image Processing 15 (6) 1506 - 1516    DOI : 10.1109/TIP.2005.871129
Srinivasan K.S , Ebenezer D. 2007 “A new fast and efficient decision-based algorithm for removal of high-density impulse noises,” IEEE Signal Process. Lett. 14 (3) 189 - 192    DOI : 10.1109/LSP.2006.884018
Schulte S , De Witte V. , Nachtegael M , Van der Weken D. , Kerre EE 2007 “Fuzzy random impulse noise reduction method,” Fuzzy Sets and Systems 158 (3) 270 - 283    DOI : 10.1016/j.fss.2006.10.010
Akkoul Smail , Ledee Roger , Leconge Remy , Harba Rachid 2010 “A New Adaptive Switching Median Filter,” IEEE Signal Processing Letters 17 (6) 587 - 590    DOI : 10.1109/LSP.2010.2048646
Pushpavalli R. , Sivaradje G. 2012 “Nonlinear Filtering Technique for Preserving Edges and Fine Details on Digital Image,” International Journal of Electronics and Communication Engineering and Technology 3 (1) 29 - 40
Fabijanska A. , Sankowski D. 2010 “Noise adaptive switching median-based filter for impulse noise removal from extremely corrupted images,” IET image processing 5 (5) 472 - 480
Pratt W.K. 1978 Digital Image Processing Wiley
Zhang S. , Karim M. A. 2002 “A new impulse detector for switching median filters,” IEEE Signal Process. Lett. 9 (11) 360 - 363    DOI : 10.1109/LSP.2002.805310
Eng H.-L. , Ma K.-K. 2001 “Noise adaptive soft - switching median filter,” IEEE Trans. Image Processing 10 (2) 242 - 251    DOI : 10.1109/83.902289
Ng Pei-Eng , Ma Kai - Kuang 2006 “A Switching median filter with boundary Discriminative noise detection for extremely corrupted images,” IEEE Transactions on image Processing 15 (6) 1500 - 1516
Wang Zhou , Zhang David 1999 ” Progressive Switching Median Filter for the Removal of Impulse Noise from Highly Corrupted Images,” IEEE Transactions on Circuits and Systems — II: Analog and Digital Signal Processing 46 (1) 78 - 80
Jang J.S.R. , Sun C. T. , Mizutani E. 1997 Neuro-Fyzzy and Soft Computing Printice-Hall
Wang Z. , Bovik A. C. 2004 “Image quality assessment:From error visibility to structural similarity,” IEEE Trans. Image Processing 13 (4) 600 - 612    DOI : 10.1109/TIP.2003.819861
Ahmad M.Omair , Sundararajan Duraisamy 1987 ”A Fast Algorithm for Two-Dimensional Median Filtering.” IEEE Transactions on Circuits and System CAS-34 (11) 1364 - 1374