当前位置:首页 >> >>

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....
2 Study of the Wiener Filter for Noise Reduction.pdf
Image De-noising for Dif... 暂无评价 4页 2.00元 Old Movies Noise Reducti...A Modified Sigma Filter ... 暂无评价 6页 免费 Abstract NOISE REDUCTION.....
Non-Iterative, Robust Monte Carlo Noise Reduction_....pdf
A Modified Sigma Filter ... 暂无评价 6页 免费 Adaptive soft-switching .....A GPU implementation of our algorithm runs in 500ms for a 512×512 image...
...PIECEWISE--LINEAR FILTER FOR IMAGE ENHANCEMENT_图文.pdf
PIECEWISE--LINEAR FILTER FOR IMAGE ENHANCEMENT_专业...in many application examples the Sigma lter both ...performing in the meantime a signi cant noise ...
sigma-delta学术报告.ppt
Σ-Δ ADCs are widely used in telecommunication ...Periodic noise suppression digit filter Typical ...modified-sinc filters for sigma-delta A/D ...
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 ...
...SIGMA-POINT PARTICLE FILTERS FOR SEQUENTIAL PROB....pdf
a bank of Sigma-Point Kalman Filters for the ...lter is called the Unscented Particle Filter and ...reduction in computational cost, making the GMSPPF...
污水零排放技术FET-WWTO-AM_图文.ppt
(414) 643-4150 (Sigma) rpeschel@thesigmagroup....modified for reason Anti-backsliding provisions ...occurs Congressional mandates WHAT’S IN A PERMIT...
Week_8_Six_Sigma (2)_图文.ppt
Design For Six Sigma (DFSS) ? Process Methods ...of this analysis is confirmed, or modified. In Six Sigma a symptom is defined as the outward...
CV模板.txt
02-07-2008 Last Modified: 02-07-2008 (07:59...(AGU) Phi Sigma Phi Phi Eta Sigma PUBLISHED ...for in such a way that recruiters find you ...
关于promodel的资料.doc
Filter Options now show all elements in the ...a Six Sigma analysis on in your model, run ...modified between scenarios for comparing the effect...
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 ...
O-1-2_图文.pdf
Both cell lines were cultured in Dulbecco’s Modified Eagle Medium (Life ...(Sigma-Aldrich, M2128), which is a colorimetric assay for measuring the ...
...of_Adenovirus_Specific_DNA_Sequence_图文.pdf
for the preparation of DNA modified electrode was illustrated in Scheme-I. ...EDC, NHS and dimercapto-succinic acid (DSA) were purchased from Sigma. ...
红花水提液对人移行细胞和鼠成纤维细胞的毒性作用_图文.pdf
filtered through a 0.2-mm filter to be ...s Modified Eagle’s Medium (SigmaAldrich, St ...in a humidified 5% CO2 atmosphere for 24, 48,...
改进型UKF滤波算法的卫星姿态估计_图文.pdf
a modified unscented Kalman filter ( U KF) algorit...noise covariance is also discussed in t he space...Reduced sigma point filters for t he propagation ...
...变换平方根的卫星姿态确定过滤器(IJCNIS-V3-N4-5)_图文.pdf
and the square root technique with modified Rodrigues...sigma points for an n-dimensional space, it is...filter in a more convenient way as in terms of...
铂_石墨烯氧还原电催化剂的共还原法制备及表征.pdf
Reduction Method for Oxygen Reduction Reactions WANG...modified Hummers method.18,19 In brief, the ...a mixed cellulose ester membrane filter (0.45 ...
Offspring from Nonstimulated Calves at an Age Young....pdf
g/ml gentamycine (Sigma Chemical Co., St. ...for 30 min at 56 C in a water bath, ...m millipore filter, and stored at 18 C until...
91 Pathway for Nitrate Assimilation in Corn (Zea ma....pdf
reduction. In corn, malate is a primary product...thioglycolate, and 2% (w/v) PVP-40 (Sigma ...spectrophotometer modified for direct illumination of...
更多相关标签: