当前位置:首页 >> >>

A Modified Sigma Filter for Noise Reduction in Images_图文

A Modi?ed Sigma Filter for Noise Reduction in Images
RADU CIPRIAN BILCU1 , MARKKU VEHVILAINEN2
1,2 Multimedia

Technologies Laboratory, Nokia Research Center Visiokatu 1, FIN-33720, Tampere FINLAND

12

Abstract: - Noise reduction is a very important processing step in digital imaging and several different techniques were proposed in the open literature. Among them, sigma ?lter has been shown to be a good solution both in terms of ?ltering accuracy and computational complexity. However, the sigma ?lter does not preserve well small edges especially for high level of additive noise. In this paper, we introduce a new sigma ?lter for noise reduction in images. We propose here a new method that ?rst decomposes the input image in four components that are independently processed using a standard sigma ?lter. The output image is reconstructed from the ?ltered components. Comparative results between our approach and the sigma ?lter, on synthetic images and also on real images obtained with a camera phone are shown. Key-Words: - Image de-noising, Image Restoration, Sigma ?lter, Noise estimation, Mean Squared Error.

1 Introduction
Due to the fast growing of the mobile devices market, there is an increased demand for high performance and robust image processing algorithms. More and more mobile imaging devices, such as camera phones and PDA’s, incorporate several image processing algorithms and methods. Although, the manufacturing technologies of the CCD or CMOS sensors are very advanced, the snapped picture must always be digitally processed prior to storage or displaying. This is done in order to eliminate different distortions such as: noise, blur, geometrical and color distortions, to mention a few. Digital images obtained by a camera phone are used for various purposes such as: sending multimedia messages (MMS), storage, printing, etc. As a consequence, there are several places where image processing algorithms can be incorporated. For instance when the images are used for MMS, image processing take place in the mobile device. Due to their limited processing power, the implemented algorithms must have a low computational complexity. On the other hand, when the snapped images are to be stored or printed, post processing can be done on a PC. In such case, usually the computational complexity and performances of the implemented algorithms can be much higher. Noise reduction is a very important processing step

in all digital imaging applications. Moreover, even for the latest manufacturing technologies, of the camera sensors, the noise level is still high. As a consequence, image de-noising is and will always be an important research topic. Among many algorithms, that exist in the open literature, the sigma ?lter [6] is probably one of the simplest de-noising method. Due to its simplicity, this ?lter represents a good choice for implementation in mobile devices. However, the edge preservation performance of the sigma ?lter is not good, especially for small image details with variance close to the variance of the additive noise. In order to improve the detail preservation of the sigma ?lter, other more sophisticated approaches were proposed. For instance the fuzzy ?lter proposed in [5] uses some fuzzy estimates of the local derivative to perform directional ?ltering of the image. Although its good ?ltering performances this approach have the disadvantage of relative high complexity and a large number of parameters that must be setup. Another alternative, called hybrid sigma ?lter, was proposed in [2] for speckle noise reduction and showed improved performances compared with the Lee’s sigma ?lter. The hybrid sigma ?lter, however, does not address the problem of additive noise reduction which is the focus of our work. In this paper we propose a new modi?ed sigma ?lter for additive noise reduction in images. In our proposed method the input image is ?rst decomposed

Figure 1: The block diagram of sigma ?ltering with noise variance estimation.

in four components and a sigma ?lter is applied separately on each of them. The output image is then reconstructed from the four ?ltered components. The paper is organized as follows: in section 2 the standard sigma ?lter is described and its advantages and disadvantages are outlined. Based on the observations from section 2, in section 3 the new ?ltering scheme is introduced. In section 4 comparative results obtained with the sigma ?lter and the new proposed method are presented and section 5 concludes the paper.

2 Sigma Filter
In this section, we brie?y review the standard sigma ?lter for additive noise [6] and outline its advantages and disadvantages. We assume the following model for the input image: y(i, j) = x(i, j) + n(i, j). (1)

where y(i, j) is the observed image, x(i, j) is the original clean image and n(i, j) is a zero mean Gaussian distributed additive noise. The main idea of the sigma ?lter is based on the fact that for a Gaussian distributed variable with mean ? and variance σ 2 a percentage of 95.5% of its samples lies inside the range [? ? 2σ, ? + 2σ]. Applying this observation to the model from (1), for every pixel y(i, j) from the observed image, a local average is computed on those neighboring pixels that are inside the interval [y(i, j) ? 2σ, y(i, j) + 2σ]. Corresponding pixel of the output image f (i, j) is replaced with this local average. This ?ltering scheme is based on the assumption

that the pixel value y(i, j) is a good estimate of the local mean and there are two issues that must be addressed. The ?rst one is the selection of the neighboring pixels. Usually a rectangular M × M window centered at the current pixel is used for this (typically with M from 3 to 9). The second more important problem is the estimation of the additive noise 2 variance σn . In practical applications the level of the additive noise is unknown therefore, some noise estimation method must be applied to the input image, prior ?ltering. For implementation in mobile devices, such methods must have low computational complexity and good estimation performances. In this paper we will use the approach in [3] that showed good estimation performance at a low computational complexity. As a consequence, a block diagram for de-noising based on sigma ?lter is depicted in Fig. 1. The input image is ?rst passed through the noise estimation module and the estimated noise variance is then used in the sigma ?lter for de-noising. It is well known, that small details of the input image are not well preserved by the sigma ?lter. This is due to the fact that on the regions from the input image that have variance close to the noise variance almost all pixels from the local M × M window are used in the average process. This effect is in?uenced mainly by the method of noise estimation but it can also be noticed when the exact noise variance is known. For instance when the estimated noise variance is larger than the real noise level the blurring effect is evident. An immediate solution is to decrease the length of the selection range to [y(i, j) ? Γσ, y(i, j) + Γσ] with Γ < 2. This modi?cation reduces also the ?ltering capabilities of the sigma ?lter in smooth areas of the input image. Moreover, the average estimator is not a good choice for monotonically increasing/decreasing regions of the input image. In such regions, a better solution would be to use a higher order polynomial. In order to avoid the use of such polynomials to model the local monotonicity another alternative is to perform some transformation to the input image prior de-noising. Taking into account these observations, in the next section we will introduce a simple modi?cation that improves the performances of the sigma ?lter for regions of the input image that contain small details.

3 The Proposed Approach
In this section we introduce our modi?ed de-noising method based on the observations outlined in the pre-

Figure 2: The block diagram of our proposed approach.

Figure 3: The block diagram of the fast implementation of our method.

vious section. The block diagram of our proposed method is depicted in Fig. 2 where the blocks denoted as HP V and HP H perform a high-pass ?ltering on the vertical and respectively horizontal directions. The blocks denoted as LP V and LP H perform a low-pass ?ltering on vertical and horizontal directions and the block denotes as ADD recombine the ?ltered components to obtain the restored image. Noise estimation is performed by the blocks denoted as N E and the sigma ?lters are represented by the blocks Sigma. Our proposed ?ltering scheme can be described by the following steps: 1. Compute the horizontal differences between adjacent pixels of the input image: yHP H (i, j) = 1 (y(i, j) ? y(i, j ? 1)) 2 (2)

?uence the ?ltering accuracy and it can be discarded at this point. In such case, it must be included in the last step of the algorithm where the four components are combined to obtain the output image. 2. Compute the horizontal weighted sum of the adjacent pixels as follows: yLP H (i, j) = 1 (y(i, j) + y(i, j ? 1)) 2 (3)

This operation is the complementary of the one in (2) and it is done primarily with the scope to reconstruct the original image. In addition to that, the edges of yLP H (i, j) are reduced and the ?ltering of this component can be made using a wider interval (larger Γ). 3. Compute the vertical differences and sums between adjacent pixels similar to (2) and (3) respectively. yHP V (i, j) = yLP V (i, j) = 1 (y(i, j) ? y(i ? 1, j)) 2 1 (y(i, j) + y(i ? 1, j)) 2 (4) (5)

where i and j are the corresponding vertical and horizontal coordinates of the pixels. This operation is performed by the block denoted as HP H in Fig. 2. Computation of these differences transforms the horizontal monotonically increasing/decreasing regions of the input image into constant regions (see [3] and the references therein). Moreover this operation also preserves the horizontal edges from the input image. Transformation of the monotonic regions into constant regions makes the simple averaging, performed by the sigma ?lter, a better model.
1 The coef?cient 2 is introduced to preserve the dynamic range of yHP H (i, j). It does not in-

4. Apply a sigma ?lter separately on the four computed components yHP H , yLP H , yHP V and yLP V to obtain fHP H , fLP H , fHP V and fLP V respectively. This is done by the blocks Sigma in Fig. 2. The four sigma ?lters, in Fig. 2, necessitate estimation of the noise variance from the corresponding image component. This is done, in the implementation from Fig. 2, separately for

Original image

Input image

Original image

Input image

Proposed

Sigma

Proposed

Sigma

Figure 4: Parts of the: original image (up left), input 2 noisy image (σn = 100, up right), ?ltered image using the proposed algorithm (bottom left) and the result of sigma ?lter (bottom right).
Input image Sigma

Figure 5: Parts of the: original image (up left), input 2 noisy image (σn = 300, up right), ?ltered image using the proposed algorithm (bottom left) and the result of sigma ?lter (bottom right).
Proposed

Figure 6: Parts of the: input image (left), ?ltered image using sigma ?lter (middle) and the result of our proposed method (right). the four components by the blocks denoted as N E. A simpler and faster implementation in which the noise variance is estimated just once and most of the processing steps are done in parallel is discussed in the sequel. 5. Reconstruct the output image from the ?ltered components as follows: f (i, j) = 1 (fHP H (i, j) + fLP H (i, j)+ 2 +fHP V (i, j) + fLP V (i, j)) (6) in the case the four components are not scaled in (2)-(5), this coef?cient must be set to 1 in (6). 4 3.1 Fast implementation The block diagram from Fig. 2 might be too complicated for implementation into mobile devices with low processing power. From our experience, most of the processing time is spend on estimation of the noise variances from the four image components separately. Moreover, computing fHP H , fLP H , fHP V and fLP V independently, necessitate four scans of the complete input image. A solution to highly decrease the running time is

1 Again the factor 2 is included to preserve the dynamic range. We must emphasize here, that

Input 10.06 24.85 50.08 74.94 100.27 125.17 150.10 175.37 198.52 224.12 250.87 274.06 300.41

Lena Sigma Proposed 7.66 6.77 14.99 12.63 25.28 20.84 34.61 28.23 42.94 35.17 50.47 41.42 57.89 47.47 65.55 53.73 71.98 58.92 78.85 64.82 86.01 70.59 91.13 75.12 98.13 81.59

Input 9.98 25.04 50.49 75.76 100.70 125.18 150.02 174.21 200.05 225.41 249.65 276.04 301.45

Cameraman Sigma Proposed 6.0126 5.53 12.57 11.23 22.35 20.06 31.06 28.08 39.33 35.51 47.89 42.92 55.65 49.94 61.93 55.68 70.68 63.89 78.18 69.76 85.71 76.69 92.34 82.82 99.50 88.98

Input 9.95 24.96 50.03 74.90 99.94 124.85 150.09 174.16 199.18 226.37 248.86 274.88 298.93

Boats Sigma Proposed 10.15 9.73 18.27 16.25 29.90 25.58 40.40 34.01 49.88 41.78 59.25 48.97 68.11 56.37 76.21 62.70 83.89 69.44 92.64 76.39 99.08 81.67 107.39 88.67 114.01 94.22

Table 1: MSE of the compared algorithms for several input images and noise variances.

?rst to implement image decomposition in a single block. Doing this parallel processing it is possible to compute the four image components in just one scan of the input image. A great reduction in the processing time is then obtained if the noise variance is not estimated separately for the four image components but is done just once. This can be implemented at the beginning of the 2 algorithm by estimating the noise variance σn from the input image y(i, j). The noise variances of the 2 four components can be then obtained from σn taking into account the linear operations performed by (2)(5) as follows (see [3] and the references therein):
2 σn (7) 2 More than that, sigma ?ltering and reconstruction of the output image from the ?ltered components can be done in just one scan of the image to further reduce the processing time. The block diagram of the fast implementation described above is depicted in Fig. 3. 2 2 2 2 σHP H = σLP H = σHP V = σLP V =

with both algorithms are shown in Table 1. From these numerical values we clearly see that our proposed algorithm provide the lower MSE especially for high levels of the additive noise. To visually compare the performances of the two approaches, in Fig. 4 and Fig. 5, we show parts of the processed images. We note that our proposed method better preserve small details (for instance the ?ne details on Lena’s hat in Fig. 4 and the ?ne details at bottom of the boats image). Comparative results obtained with both algorithms on one image obtained with a Nokia cameraphone are shown in Fig. 6. Also here we can see better preservation of the ?ne details when the proposed de-noising method is applied.

5 Conclusions
In this paper we have introduced a new sigma ?lter for image de-noising. The new algorithm has improved performances in terms of MSE and also preserves better the ?ne details of the processed image as opposed with the standard sigma ?lter. A fast implementation with low computational complexity was also presented. The processing time of the new method is not much higher that that of the sigma ?lter and this makes it suitable for practical implementations on mobile devices. Further developments, of the new algorithm, with improved performances for signal dependent noise, are under consideration. References: [1] J. Astola, P. Kuosmanen, Fundamentals of Non-

4 Simulations and Results
In this section we show the comparative performances of our proposed approach and the standard sigma ?lter. To this end, we selected three images (lena, boats and cameraman) and we added zero mean Gaussian noise with different variances to them. The original images were represented on 8 bits (values in the range [0, 255]). The Mean Squared Error (MSE), between the ?ltered image and the original clean image, obtained

linear Digital Filtering, New York: CRC Press, 1997. [2] L. Alparone, S. Baronti, A. Garzelli, A Hybrid Sigma Filter for Unbiased and Edge-Preserving Speckle Reduction, in Proceedings of IEEE International Geoscience and Remote Sensing Symposium, IGARS ’95, Vol. 2, 1995, pp. 1409-1411. [3] R. C. Bilcu, M. Vehvilainen, A New Method for Noise Estimation in Images, in Proceedings of IEEE-EURASIP International Workshop on Nonlinear Signal and Image Processing, NSIP 2005, May 18-20, 2005, Sapporo, Japan, pp: 290-293. [4] A. Bosco, K. Findlater, S. Battiato, A. Castorina, A Noise Reduction Filter for Full-Frame Data Imaging Devices, in IEEE Transactions on Consumer Electronics, Vol. 49, No. 3, August, 2003, pp. 676-682. [5] D. Van De Ville, M. Nachtagael, D. Van der Weken, E. E. Kerre, W. Philips, I. Lemahieu, Noise Reduction by Fuzzy Image Filtering, in IEEE Transactions on Fuzzy Systems, Vol. 11, No. 4, August, 2003, pp. 429-436. [6] J. S. Lee, Digital Image Smoothing and the Sigma Filter, in Computer Graphics Image Processing, Vol. 24, 1983, pp. 255-269. [7] K. Deergha Rao, M.N.S. Swamy, E.I. Plotkin, Adaptive Filtering Approaches for Color Image and Video Restoration, in IEE Proceedings of Visual Image Signal Processing, Vol. 150, No. 3, June, 2003, pp. 168-177. [8] D.-H. Shin, R.-H Park, S. Yang, J.-H. Jung, Block-Based Noise Estimation Using Adaptive Gaussian Filtering, in IEEE Transactions on Consumer Electronics, Vol. 51, No. 1, February, 2005, pp. 218-226. [9] Z. Ma, H. R. Wu, B. Qiu, A Window Adaptive Hybrid Vector Filter for Color Image Restoration, in Proceedings of IEEE International Conference on Acoustics, Speech, and Signal Processing, ICASSP 2004, May 17-21, 2004, Vol. 3, pp. 205-208.


相关文章:
SNR Improvement by Photon Noise Filtering in Ocean ....pdf
paper provides a new photon noise filtering which combines Anscombe transformation [6] and modified Sigma filter [13] for photon noise reduction in images....
imatest -Noise_in_photographic_images.pdf
A Modified Sigma Filte... 暂无评价 6页 免费 A...film and pixel level variations in digital images...Noise appears relatively constant except for the ...
2 Study of the Wiener Filter for Noise Reduction.pdf
2财富值 LOW DELAY FILTER FOR ... 暂无评价 4页 免费 A Modified Sigma Filter ... 暂无评价 6页 免费 Abstract NOISE REDUCTION... 暂无评价 2页 免费喜...
Non-Iterative, Robust Monte Carlo Noise Reduction_....pdf
A Modified Sigma Filter ... 暂无评价 6页 免费 Adaptive soft-switching .....[1995] found that most Monte Carlo noise appeared in images as inter-...
Noise Reduction Prefiltering for Video Compression.pdf
in 8x8 block units, so that filter can be weakened or turned off for ...a wide range sigma of white noise and tau and gamma of the filter to ...
Oscillators_SJTC_图文.pdf
? 1. Narrowed Loop Filter Bandwidth ? Reduction ...Spreading of Quantization Noise 0 Sigma Delta A/...Modified Colpitts Oscillators RF cut (Yc=Gc+iBc)...
直升机设计-Rotorcraft Design 4_图文.ppt
Weight (Modified Configuration) Modified GW: (...sigma Noise 4 13.21 ft .563 ft -9 deg VR-...in 10^ -5 3 @RISK Trial Version For ...
A phenol biosensor based on modified core-shell at ....pdf
A phenol biosensor based on modified core-shell ...(APTES) were purchased from SigmaAldrich. The...and ground together in an agate mortar for about...
AN3494.pdf
A/D, Analog to Digital, A to D FILTER CIRCUITS...(also called sigma delta), typically used for ...This reduction in the requirements is the same ...
基于修正罗德里格参数的小卫星SSUKF姿态确定_图文.pdf
A AttitudeEstimationUsingModifiedRodriguesParameters ...decreasedthroughSSUKFwhilethefilter performancewas...UHLMANNJK.Reducedsigmapointthe F.Anewapproachfor...
【图文】面向天基光学监视的空间目标 TLE拟合与跟踪方法.pdf
Elements( TLE) conversion were modified for the bearing 2 only passive ...in the procedure of the tracking filter, and sigma points kalman fil2 ...
...nanodimensional structures: challenges a_图文.pdf
We developed a method for producing GaN nanofibers...(Sigma Aldrich) in a homogeneous mixture of ...A modified JEOL 6400 was used to perform the ...
M. Multiple communication protocols for software de....pdf
modified with an interpolation filter and a Phase...for a single protocol standard in software and ...Sigma-delta converters have been limited to 100 ...
A graphene modified anode to improve the performanc....pdf
A graphene modified anode to improve the ...a desired property for application in the MFCs. ...Sources 196 (2011) 54025407 5403 Sigma...
...onto zeolite by means of a modified_图文.pdf
of a modified_环境科学/食品科学_工程科技_专业...ammonium hydroxide were purchased from Sigma Aldrich...By indexing images in Fig 7, and analyzing the...
Nanoscopic LiquidLiquid Interface Arrays 支持材料_图文.pdf
catalyst (10 wt% Pd) were purchased from Sigma...modified Hummers’ method.2 In a typical ...Samples for TEM images were prepared by depositing...
一种新型节省芯片面积的数字插值滤波器的设计与实现_英文_.pdf
noise , and a 1bit D/ A for o ut2of2band...1 Block diagram of sigma2delta DAC In mixed ...Such a filter ’ s main ad2 vantage lies in ...
...for a Laser Patterned Neuronal Culture_图文.pdf
once the cells attach to the modified substrate....allows for patterning cells in a process known as...(PLL; Mw = 150,000; Sigma; St. Louis, MO)...
...analyses of the exopolysaccharides produced by a.pdf
10.0 and agar 15.0 in filter sterilized ...production were grown in a modified SWN medium. ...(Sigma) using 50 mM Phosphate buffer pH 7.0 as...
用MSP430F11x1 的Sigma-Delta 型毫伏表.pdf
12-bit sigma-delta type A/D converter using the...The feedback loop can be modified in length to...for the 1-bit DAC low-pass filter is important...
更多相关标签: