Comparative Study of Noise Removal Algorithms for

13 downloads 0 Views 365KB Size Report
filter, Adaptive median filter, which are design using LabView. These filters are used ... But here we used labview for designing various ... obvious window is just the first few preceding and following entries .... Robotics Institution Superior Tech.
2015 International Conference on Computational Intelligence and Communication Networks

Comparative Study of Noise Removal Algorithms For Denoising Medical Image Using LabVIEW Sambit Satpathy, Mohan Chandra Pradhan, Subrat Sharma Department of embedded system design, Sambalpur University Institute of Information Technology, burla, india E-mail: [email protected], [email protected], [email protected] Medical imaging is the technique and process of creating visual representations of the interior of a body for clinical analysis and medical intervention. Medical imaging seeks to reveal internal structures hidden by the skin and bones, as well as to diagnose and treat disease. Medical imaging also establishes database of normal anatomy and physiology to make it possible to identify abnormalities

Abstract—Biomedical images are generally tainted by Gaussian noise, Salt & Pepper noise. Gaussian noises are additive where as salt & pepper noises are itself as separately occurring white & black pixels. This paper presents the study of five types of filters like Gaussian filter, Gabor filter, Box filter, Median filter, Adaptive median filter, which are design using LabView. These filters are used for removing of two type of noises like Gaussian noise, Salt & Pepper noise. Then respective PSNR value has been found out which is used for the comparative analysis of the filter.

G(x,y)

Keywords- Gaussian noise, Salt & Pepper noise, Gaussian filter, Gabor filter, Box filter, Median filter, Adaptive noise filter, PSNR ratio.

I.

A.

Denoising technique

Z(x,y) Fig-1 model of denoising III.

BASIC FILTERS

A. GAUSSIAN FILTER It is a filter whose impulse response is a Gaussian function. Image denoising is an important image processing task, both as a process itself, and as a component in other processes. Very many ways to denoise an image or a set of data exists. The main properties of a good image denoising model are that it will remove noise while preserving edges. Traditionally, linear models have been used. One common approach is to use a Gaussian filter. For filtering other than Gaussian filter we have used Gabor filter, box filter, median filter, adaptive median filter. Gaussian filters have the properties of having no overshoot to a step function input while minimizing the rise and fall time. This behavior is closely connected to the fact that the Gaussian filter has the minimum possible group delay. It is considered the ideal time domain filter, just as theism is the ideal frequency domain filter. Mathematically, a Gaussian filter modifies the input signal by convolution with a Gaussian function; mathematically it is represented by

ARCHITECTURE

MEDICAL IMAGE

Image denoising is an important image processing task, both as a process itself, and as a component in other processes. Very many ways to denoise an image or a set of data exists. The main properties of a good image denoising model are that it will remove noise while preserving edges. Traditionally, linear models have been used. One common approach is to use a Gaussian filter.

B. IMAGE DENOISING For filtering other than Gaussian filter we have used Gabor filter, box filter, median filter, adaptive median filter. possible , such as "box" or "cross" patterns. 978-1-5090-0076-0/15 $31.00 © 2015 IEEE DOI 10.1109/CICN.2015.67

W(x,y)

Fig-1 model of denoising

INTRODUCTION

Noise may be arised in creation, transmission, capturing process of images. Especially we considered medical images because these many importance in day-today life.In practical life various noises are added in medical images, like Gaussian noise, salt & pepper noise etc. In general we face several problems in MRI scan, x-ray, ultra sound, brain hemorrhage scanning. Here with real images various noises are added. So that images are not clearly visible. To overcome this process, we have to design various filters, through which noises are removed and we can see clear images. There are various soft wares available for filter designing, matlab gives convenient filter design tools. But here we used labview for designing various types of filters.Labview provides comprehensive tools that can be used to build any measurement or control application in significantly less time, it also gives idea about development environment for Finnovation, detection & accelerated results. II.

Linear operation

300

structure in the image. These filters smooth the data while keeping the small and sharp details. The median is just the middle value of all the values of the pixels in the neighborhood. Note that this is not the same as the mean instead; the median has half the values in the neighborhood larger and half smaller. The median is a stronger "central indicator" than the average. In particular, the median is hardly affected by a small number of discrepant values among the pixels in the neighborhood.

Where x= the distance from the origin in the horizontal axis Y= the distance from the origin in the vertical axis. σ = is the standard deviation of the Gaussian distribution.

IV.

B. GABOR FILTER

LABVIEW MODEL DESIGNING AND IMPLEMENTATION

A. GAUSSIAN FILTER

Gabor filter named after Dennis Gabor, is a filter used for edge detection. Frequency and orientation representations of Gabor filters are similar to those of the human visual system, and they have been found to be particularly appropriate for texture representation and discrimination. In the spatial domain, a 2D Gabor filter is a Gaussian kernel function modulated by a sinusoidal plane wave. Simple cells in the visual cortex of mammalian brains can be modeled by Gabor functions. Thus, image analysis with Gabor filters is thought to be similar to Perception in the human visual system. C. BOX FILTER Box filter is a spatial domain linear filter in which each pixel in The resulting image has a value equal to the average value of its neighboring pixels in the input image. It is a form of low-pass filter and also called box linear filter which can written in this way like 3*3 Matrix of 1/9 * determinant matrix Due to its property of using equal weights.

Fig-2 Gaussian function implementation

B. GABOR FILTER D. MEDIAN FILTER The median filter is a nonlinear filtering technique, often used to remove noise. Such noise reduction is a typical pre-processing step to improve the results of later processing for example, edge detection on an image. Median filtering is very widely used in digital image processing because, under certain conditions, it preserves edges while removing noise . The main idea of the median filter is to run through the signal entry by entry, replacing each entry with each entry with the median of neighboring entries. The pattern of neighbors is called the "window", which slides, entry by entry, over the entire signal. For 1D signals, the most obvious window is just the first few preceding and following entries, whereas for 2D signals such as images, more complex window patterns are possible , such as "box" or "cross" patterns.

E. ADAPTIVE MEDIAN FILTER The median filter is a nonlinear filtering technique, often used to remove noise. Such noise reduction is a typical pre-processing step to improve the results of later processing for example, edge detection on an image. Median filtering is very widely used in digital image processing because, under certain conditions, it preserves edges while removing noise . The main idea of the median filter is to run through the signal entry by entry, replacing each entry with We see various application of median filter, As i compare it with advanced standard median filtering, the Adaptive Median Filter performs vital processing to preserve detail and smooth non-impulsive noise. A basic benefit to this adaptive approach of median filtering is that repeated applications of this Adaptive Median Filter do not erode away edges or other small

Fig-3 Gabor function implementation

A. BOX FILTER(3*3)

301

C.

ADAPTIVE MEDIAN FILTER USING MATLAB SCRIPT

Fig-4 box filter (3*3) implementation

C. BOX FILTER (11*11) Fig-7 block diagram of adaptive median filtering using matlab script. V.

EXPERIMENT

In this work we have collected real time medical images which contain salt & pepper noise, Gaussian noise. Then we passing medical noisy images through Gaussian filter, Gabor filter, box filter, median filter, adaptive median filter. As a result we get denoised images. Below we taken two Gaussian noise medical images and denoised images by Gaussian filter, Gabor filter, box filter, median filter, adaptive median filter respectively.

A.

ORIGINAL NOISED IMAGES

Fig-5 box filter (11*11) implementation

B. MEDIAN FILTER USING IMAQ TOOL

Fig-8.1 image containing Gaussian noise

Fig-6 block diagram of median filter using imaq Tool

302

Fig-9.3 Box filter image

Fig-8.2 image containing salt & pepper noise B. DENOISED IMAGES(REMOVING GAUSSIAN NOISE)

Fig-9.4 Median filter image

Fig-9.1 Gaussian filter image

Fig-9.5 Adaptive median filter image C. DENOISED IMAGES (REMOVING SALT & PEPPER NOISE)

Fig-9.2 Gabor filter image

303

Fig-10.5 Adaptive median filter image Fig-10.1 Gaussian filter image VI.

COMPARISONS OF PSNR RATIO

Noise densit y

Noisy image s

0.05

14.264 9 12.602 5 10.550 8 9.7685 1 8.7635 7 8.1773 5

0.10 0.15 0.20 0.30 0.40

Fig-10.2 Gabor filter image

Denoise d image using Gaussia n filter 23.7703

Denoise d image using Gabor filter 23.7785

Denoise d image using Box filter 5.71139

Denoise d image using Median filter 21.2472

Denoise d image using AMF filter 18.3179

21.3328

21.3354

5.36576

18.5115

15.2925

19.7279

19.7281

5.11506

16.8979

14.5096

18.6728

18.6723

22.3983

15.7698

13.6797

17.1371

17.1356

15.6226

14.2295

12.6284

16.2774

16.2774

13.2595

13.2595

12.0707

Fig-11 PSNR Ratio for 2d b (Gaussian noise)

Fig-10.3 Box filter image

Noise density

Noisy images

17.7227

Denoised image using Gaussian filter 26.3213

Denoised image using Gabor filter 26.3414

Denoised image using Box filter 5.95379

Denoised image using Median filter 38.1054

Denoised image using AMF filter 42.6656

0.05 0.10

14.7031

24.6245

24.6352

5.73889

36.6994

40.4358

0.15

12.864

22.8758

22.8807

5.51272

33.6396

37.8844

0.20

11.637

21.4984

21.5004

5.32434

29.7598

36.8312

0.30

9.89726

19.1226

19.1225

25.3836

23.45

33.5118

0.40

8.64189

17.1717

17.1702

15.4086

18.5472

31.3005

Fig-12 PSNR Ratio for 2d b (Salt & pepper noise) VII.

RESULT & DISCUSSION

In this research work, we have used LabView for image processing, and found it suitable. Here, we implemented five filters i.e. Gaussian filter, Gabor filter, Box filter, median filter and adaptive median filter for removing various noises (Gaussian and salt & pepper noise) from medical images. PSNR value is used for the comparative analysis of these filters. Based on PSNR values with different noises and concentration of noise present in images we observed that, adaptive median filter gives better results as compared to others for impulsive noise (salt & pepper noise). It also found that in Gabor filter when center frequency tends to 0.02

Fig-10.4 Median filter image

304

Hz it gives the maximum results for Gaussian noise otherwise the PSNR value decreases. Again as compared to Gaussian filter in Gabor filter, we found that in maximum cases PSNR value is higher by taking σ = 2. The performance of box filter is good for particular noise level. But in maximum cases Gaussian performance is much better than the box filter (depending on level of noise).

of Imaging Sc. And tech, Delft Uni of Technology [4] Maria Isabel Ribeiro, “Gaussian probability Density Functions: Properties and Error Characterization” Institute for System and Robotics Institution Superior Tech. [5] Nezamoddin Kachouri, Paul Fieguth “A Gabor based technique for Image Denoising Saskatoon, Mau 2005 [6] Selective removal of impulse noise based on homogeneity level information, Gouchol Pok; Jyh Charn Liu; Nair, A.S; Image Processing, IEEE Transactions on , Volume: 12 Issue: 1, Jan. 2003: Page(s): 85 -92 [7] Impulse noise reduction in medical images with the use of switch mode fuzzy adaptive median filter Abdullah Toprak ; ˙Inan Güler; Digital Signal Processing 17 (2007) 711–723, ELSIVIER [8] LabVIEW Graphical Programming (2nded.), Gary W. Johnson, 1997, McGraw Hill

REFERENCES [1] R. C. Gonzalez and R. E. Woods, “Digital Image Processing”, Pearson Education, India, 2002. [2] Ying Shen . Weihua Zhu .“ Medical Image Processing using A Machine Vision-based Approach”. International journal of signal processing, Image processing and Pattern Recognition, Vol.6, No.3, June 2013. [3] Tuan Q. Pham Lucas J. Van Viet, “Separable Bilateral Filtering for Fast Video Preprocessing”, Quantitative Imaging Group Dept

.

305