Image Denoising of Gaussian and Poisson Noise Based on Wavelet Thresholding
A dissertation submitted to the
Graduate School
of the University of Cincinnati
in partial fulfillment of the
requirements for the degree of
DOCTOR OF PHILOSOPHY (Ph.D.)
in the School of Electronic and Computing Systems
University of Cincinnati
Cincinnati OH 45221 USA
2013
by
Jin Quan
Bachelor of Science in Engineering, Tongji University
Shanghai, China, 2007
Committee chair: Dr. William G. Wee Abstract
Noise on images is generally undesirable and disturbing. It always plays a negative role on higher level processing tasks such as image registration and segmentation.
Thus,image denoising becomes a fundamental step necessarily required for better im- age understanding and interpretation. During the last couple of years, wavelet has been extensively employed in the application of suppressing noise and proven to be a successful tool which outperforms many conventional denoising filters due to its pre- ferred properties. Therefore, in this dissertation, the wavelet transform is applied to develop our denoising strategy.
Basically, two generic scenarios occur during the acquisition of images. First, when the detected intensities on the image are sufficiently high, the noise can be suitably modeled as following an additive independent Gaussian distribution. Second, when only a few photons are detected, this observed image is usually modeled as a Poisson process and the intensities to be estimated are assumed to be the underlying Poisson parameters. In this dissertation, these two scenarios are discussed respectively in Part
IandPartII.
In part I, we consider to reduce the typical additive white Gaussian noise (AWGN).
Our driving principle is to decrease the upper bound of the error restricted by the soft-thresholding strategy between the investigated image and noise-free image. Thus we develop a new context modeling method to group coefficients with similar statis-
ii tics and construct a smoothed version of the noisy image prior to the actual denoising operation. Then, we propose an optimized soft-thresholding denoising function with parameters derived from a modification of a closed form solution which has a more
flexible shape and is adaptively pointwise. Furthermore, we extend it to its overcom- plete representation by employing the “cycle spinning” method so that the property of shift invariance is achieved which leads to a boost of the denoising performance. By combining these strategies, the denoising results in our experiments confirm that the approach is very competitive to some state-of-the-art denoising methods in terms of quantitative measurements and computational simplicity.
In Part II, a new denoising method for Poisson noise corrupted images is proposed which is based on the variance stabilizing transformation (VST) with a new inverse.
The VST is used to approximately convert the Poisson noisy image into Gaussian distributed, so that the denoising methods aiming at Gaussian noise can be applied subsequently. The motivation for the improved inverse comes from a main drawback existing in the conventional VSTs such as the Anscombe transformation: its efficiency degrades significantly when the pixel intensities of the observed images are very low due to the biased errors generated by its inverse transformation. In order to correct the biased errors, we introduce a polynomial regression model based on weighted least squares as an alternate to its inverse. Moreover, we incorporate our developed wavelet thresholding strategy for Gaussian noise presented in Part I into the proposed method.
We also extend it to the overcomplete representation to suppress the Pseudo-Gibbs phenomena and therefore gains additional denoising effects. Experimental analysis indicates that this method is very competitive.
iii iv Acknowledgments
I would like to express my deepest gratitude and most sincere appreciation to my advisor Dr. William G. Wee for all his continuous guidance, constant encouragement and precious support during my stay at the University of Cincinnati. In every sense, he was always a fruitful source of inspiration from which I tremendously benefited. I am also very grateful to the members of the dissertation committee: Dr. Chia Y. Han, Dr.
Xuefu Zhou, Dr. Raj Bhatnagar and Dr. Ali Minai for kindly sharing their scientific knowledge, devoting priceless time and providing constructive advice and comments on my research. I learned a lot from enthusiastic discussions with them.
In addition, many thanks go to the former and current members of the Multimedia and Augmented Reality lab at the University of Cincinnati for their help, kindness and valuable suggestions.
Most of all, I owe my warmest thanks to my father Dr. Shuhai Quan and my mother Mrs. Yulan Ruan. Nothing could have been possible without their everlasting love, understanding and support.
v Table of Contents
Abstract ii
Acknowledgments v
List of Figures ix
List of Tables xiv
1 Introduction 1
1.1ImageDenoising...... 1
1.2ProblemStatement...... 2
1.3 Research Scope ...... 4
1.4Contributions...... 5
1.5DissertationOrganization...... 6
Part I Image Denoising for Additive Gaussian Noise 9
2 Background 10
2.1 Additive Gaussian Noise Model ...... 10
2.2ImageQualityEvaluation...... 11
2.2.1 Objective Image Quality Evaluation ...... 11
vi 2.2.2 Subjective Image Quality Evaluation ...... 13
2.3 Summary ...... 14
3 Literature Review 15
3.1SpatialDomainApproaches...... 15
3.1.1 LinearFilters...... 16
3.1.2 Nonlinear Filters ...... 17
3.2TransformedDomainApproaches...... 20
3.2.1 FourierTransformDenoising...... 20
3.2.2 WaveletTransformDenoising...... 21
3.2.3 Data-AdaptiveTransformDenoising...... 50
3.3OtherAlternativeApproaches...... 51
3.4 Summary ...... 53
4 The Proposed Gaussian Denoising Method: CWMT 55
4.1Motivation...... 55
4.2Overview...... 58
4.3 Denoising Operation 1–Improved Context Modeling ...... 59
4.4ExperimentalResultsofDenoisingOperation1 ...... 63
4.5 Denoising Operation 2–The Optimization of the Soft Thresholding Func-
tion...... 66
4.6ExperimentalResultsofDenoisingOperation2 ...... 69
4.7 Experimental Results of Combining Two Denoising Operations ..... 71
4.8TheStepsoftheDenoisingProposedMethod:CMWT...... 75
4.9 Summary ...... 79
5 Expansion to the Overcomplete Representation 81
5.1 Overview of Applying Overcomplete Expansion ...... 81
vii 5.2 Overcomplete Expansion Procedure ...... 83
5.3 Experimental Results for Overcomplete Expansion ...... 83
5.4 Summary ...... 91
Part II Image Denoising for Poisson Noise 92
6 Background on Poisson Denoising 93
6.1 Modeling of Low Intensity Images ...... 93
6.2PoissonNoiseModel...... 94
6.3RelatedWork...... 95
6.3.1 Variance Stabilization ...... 95
6.3.2 HypothesisTesting...... 98
6.3.3 WaveletFiltering...... 99
6.3.4 BayesianBasedApproach...... 100
6.4AnscombeTransformationandItsInversions...... 101
6.5 Summary ...... 104
7 The Proposed Poisson Denoising Method: CMWT-IAT 105
7.1InvestigationoftheBiasedErrors...... 106
7.2NewInverseTransformationforAnscombeTransformation...... 108
7.3CombinationoftheDenoisingMethodCMWTforAWGN...... 112
7.4 Summary ...... 113
8 Performance Evaluation 114
8.1ComparisonswithTwoConventionalInversions...... 115
8.2ComparisonswithSURE-LETUsingtheProposedInversion...... 115
8.3ComparisonswithState-of-the-ArtDenoisingMethods...... 117
8.4 Summary ...... 121
viii 9 Conclusion and Perspectives 125
9.1Conclusion...... 125
9.2Perspectives...... 126
Bibliography 130
ix List of Figures
1.1 (a) A noise-free Image Pepper,(b)Anoisyversionofit...... 2
3.1 Original Image Lena anditsFourierdecomposition...... 21
3.2Onescaleofwaveletdecomposition...... 22
3.3Somefamouswavelets...... 23
3.4 Original Image Lena and its first level decomposition by using db4 wavelet 24
3.5 Hard-thresholding function ...... 29
3.6 Soft-thresholding function ...... 30
3.7 Semisoft-thresholding function ...... 31
4.1Awaveletdenoisingflowchart...... 57
4.2 The parent-child relationship of a three level wavelet decomposition. . . 60
4.3 Subband designations ...... 61
4.4 Six standard testing images used on our experiments (a) Lena,(b)Boat,
(c) Goldhill,(d)Barbara,(e)Couple,(f)Man...... 64
4.5 Subbands of the 2D orthogonal wavelet transform ...... 65
4.6 Sensitivity of the denoising function with respect to variations of T on
Lena ...... 69
4.7 Sensitivity of the denoising function with respect to variations of T on
Boat ...... 70
x 4.8 Visual comparison between original soft-thresholding and the optimized
soft-thresholding functions on Image Goldhill. (a) Original Goldhill,
(b) Noisy version of noise level 30, (c) Denoised image by original soft-
thresholding function, (d) Denoised image by the optimized soft-thresholding
function...... 72
4.9 Visual comparison between original soft-thresholding and the optimized
soft-thresholding functions on Image Barbara. (a) Original Barbara,
(b) Noisy version of noise level 30, (c) Denoised image by original soft-
thresholding function, (d) Denoised image by the optimized soft-thresholding
function...... 73
4.10 Comparison of PSNR of CMWT and SURE-LET [1] on Lena ...... 74
4.11 Comparison of PSNR of CMWT and SURE-LET [1] on Goldhill .... 74
4.12 Comparison of PSNR of CMWT and SURE-LET [1] on Couple .... 75
4.13 Visual comparison on Image Lena. (a) Original image, (b)Noisy image
of noise level 20, (c) Denoised image by SURE-LET [1], (d) Denoised
imagebyCMWT...... 76
4.14 Visual comparison on Image Goldhill. (a) Original image, (b)Noisy im-
age of noise level 20, (c) Denoised image by SURE-LET [1], (d) Denoised
imagebyCMWT...... 77
4.15 Visual comparison on Image Couple. (a) Original image, (b)Noisy image
of noise level 20, (c) Denoised image by SURE-LET [1], (d) Denoised
imagebyCMWT...... 78
5.1 Relation between shifted times and PSNR of the denoised Image Lena .84
5.2 Relation between shifted times and PSNR of the denoised Image Goldhill 84
5.3 Comparison of PSNR (dB) with 3 other most efficient methods on Lena 85
5.4 Comparison of PSNR (dB) with 3 other most efficient methods on Goldhill 85
xi 5.5 Comparison of PSNR (dB) with 3 other most efficient methods on Couple 86
5.6 The Image Lena (a) Noise-free image, (b) Noisy image with noise level
of 50, (c) Denoised image using SURE-LET [62], (d) Denoised image
using BLS GSM [3], (e) Denoised image using BM3D [92], (f) Denoised
imageusingCMWT-OE...... 87
5.7 The Image Boat (a) Noise-free image, (b) Noisy image with noise level
of 50, (c) Denoised image using SURE-LET [62], (d) Denoised image
using BLS GSM [3], (e) Denoised image using BM3D [92], (f) Denoised
imageusingCMWT-OE...... 88
5.8 Wavelet domain images of Barbara at: (a) the first scale (b) the second
scale(scaledupbytwo)(c)thethirdscale(scaledupbyfour)...... 89
6.1GeneralPoissonnoisyimagedenoisingprocedure...... 96
7.1 Variance of Poisson distributed data sets ...... 107
7.2Varianceoftransformeddatasets...... 107
7.3BiasederrorsbetweenPoissonparametersandestimatedmeans.... 108
7.4 Curve fitting for Poisson parameter under 10 ...... 111
7.5CurvefittingforPoissonparameterfrom10to30...... 111
8.1 (a) The original Image Boat at peak intensity 30, (b) Poisson noise
corrupted image, (c) Image denoised with non-overcomplete SURE-LET
[1] and the proposed inversion, (d) Image denoised with CMWT-IAT. . 118
8.2 (a) Part of the original Image Man at peak intensity 30, (b) Poisson noise
corrupted image, (c) Image denoised with non-overcomplete SURE-LET
[1] and the proposed inversion, (d) Image denoised with CMWT-IAT. . 119
xii 8.3 (a) The original Image Goldhill at peak intensity 20, (b) Poisson noise
corrupted image, (c) Image denoised with overcomplete SURE-LET [62]
and the proposed inversion, (d) Image denoised with PURE-LET [6], (e)
Image denoised with BM3D and their unbiased inversion [7], (f) Image
denoised with our denoising method [109] and the proposed inversion. . 123
8.4 (a) The original Image Couple at peak intensity 10, (b) Poisson noise
corrupted image, (c) Image denoised with overcomplete SURE-LET [62]
and the proposed inversion, (d) Image denoised with PURE-LET [6], (e)
Image denoised with BM3D and their unbiased inversion [7], (f) Image
denoised with our denoising method [109] and the proposed inversion. . 124
xiii List of Tables
3.1 Comparison of several most famous thresholding denoising methods (PSNR)
...... 36
3.2 Comparison of the aforementioned image denoising methods(PSNR) . 45
4.1 Table for selecting A at different subbands and estimated noise levels . 63
4.2 Comparison of PSNR of 3 images from a reconstruction of Znew(i,j) and
the input PSNR ...... 63
4.3 Table for selecting Δa1 under different subbands and noise levels. . . . 68
4.4 Table for selecting Δa2 under different subbands and noise levels. . . . 68
4.5 Comparison of PSNR (dB) with Soft-thresholding (non-redundant) . . 71
4.6 Comparison of PSNR (dB) with the method SURE-LET(non-redundant)
...... 79
5.1 Comparison of several most famous thresholding denoising methods (PSNR)
...... 90
8.1 PSNR (dB) comparison of the arithmetical, asymptotical inverse Anscombe
transformationsandtheproposedinversetransformation...... 116
8.2 PSNR comparison of SURE-LET [1] and our denoising method [109]
combined with the proposed inversion for different images and peak
intensities...... 117
xiv 8.3 PSNR comparison of some of the best denoising methods for different
imagesandpeakintensities...... 121
xv Chapter 1
Introduction
1.1 Image Denoising
Images obtained from the real world are always mixed with noise. The noise brought in is derived from multiple sources. The imperfect instrument itself would produce a certain amount of noise when the image is taken. When transforming the optical signal into a digital signal, the pixel’s value at specific location is dependent to the number of photons the corresponding captor has received. So the instability of the number of receiving photons can cause the production of noise. Moreover, during image’s am- plification and transmission, additional perturbations can be introduced by electronic devices and transmission lines.
There are several different types of noise in digital images. For instance, shot noise is generated by the random way photons are emitted from a light source especially when the light intensity is limited and it is usually characterized by Poisson distribu- tion. Thermal noise, also known as dark current noise, is produced by thermal agitation of electrons at sensing sites and highly dependent on the sensor’s temperature and the exposure time. Images with impulsive noise, which is generally caused by the mal- functioning of elements in the camera sensors or timing errors in the data transmission
1 process, have bright pixels in dark areas and dark pixels in bright areas. And quanti- zation noise often happens due to the errors when an analog signal is converted to a number of discrete digital values.
Since noise seriously compromises the details of the image and hampers image un- derstanding and image analysis in scientific and commercial applications (see Figure
1.1), image denoising is extensively required. Thus it is highly necessary to use an appropriate and efficient denoising approach to eliminate or reduce noise while keeping the important image features when pre-processing images.
Figure 1.1: (a) A noise-free Image Pepper, (b) A noisy version of it
1.2 Problem Statement
Image denoising attempts to recover a noise-free image by eliminating or reducing the noise on the observed image. This processing can be modeled as obtaining an optimal estimate of the unknown noise-free image from the available noise-corrupted image. A large number of scientific literatures have emphasized on image denoising in the last decade and there is still existing a wide range of interest in the subject nowadays.
Although various algorithms and tools have been proposed, derived and improved, the
2 problem is that many denoising techniques always suffer over-softening the crucial im- age features as well as introducing artifacts. Thus the searching for an efficient image denoising method is still a challenging task.
Besides, the amount of noise usually depends on the signal intensity. Practitioners often consider it to be following a statistical distribution. Generally, when the magni- tude of the measured signal is sufficiently high, the noise is supposed to be independent of the original image that it corrupts, and modeled as an additive Gaussian random variable. On the other hand, when the magnitude of the observed signal is relatively low, it is often assumed to follow a Poisson distribution.
Thus, the general goal of this research is to design and implement an efficient im- age denoising method for Gaussian and Poisson noise, which can satisfy the following requirements:
• Competitive performance
The proposed algorithm should be competitive with other start-of-the-art denoising methods according to certain objective measurements, such as Peak Signal-to-Noise
Ratio (PSNR). It should also satisfy human visual assessment.
• Minimal human interaction
The human interaction should be minimized during the denoising process when ap- plying the proposed algorithm, in other words, the entire denoising process should be totally automatic.
• Low computational burden
The proposed algorithm should not require a very high computing capacity, a regular personal computer should satisfy the hardware requirement and be qualified to com- plete the whole process in a short period of time.
• Adequate reliability
The proposed algorithm should demonstrate consistent and repeatable experimental results regardless of the sources of images and how many times the denoising process
3 is performed.
• Wide targeted application
The major application of the algorithm is to denoise natural images corrupted by additive Gaussian noise. Under certain conditions, it should also be applicable to the denoising of images corrupted by other types of noise, such as Poisson noised corrupted images, which are commonly obtained in astronomy and biomedical research.
1.3 Research Scope
Although ultrasound images are more probably to be corrupted by speckle noise and the statistical distribution of the noise found in MRI images is more like Rician, most of the noise obtained during acquisition and transmission of the natural images are assumed to be additive white Gaussian noise (AWGN), so in particular, the scope of
Part I of this research is to focus on the suppression of AWGN on natural images.
Meanwhile, The image type focused on in the research is grayscale, though we make occasional references to papers on color image denoising. The images used for experimental purposes are all standard grayscale and natural testing images. These grayscale images contain 8 bit data which means the brightness levels are from 0∼255.
In addition, in order to broaden the range of applications of our denoising method, in Part II of the dissertation, we extend the work of this research to the denoising of non Gaussian-corrupted images, in particular, the denoising of low intensity Poisson- corrupted images. Poisson-corrupted images are commonly acquired in biomedical imaging and astronomy research when the intensity of the light source is limited and the exposure time is short. And noise on these low intensity images is often modeled as Poisson distribution.
4 1.4 Contributions
In this dissertation, we address the image denoising problem with focuses on the re- moval of additive Gaussian noise and Poisson noise. By researching on the wavelet transform and soft-thresholding strategy, we develop a very competitive image denois- ing method. Here we present the original contributions of this dissertation consisting of the following:
1. Many effective denoising methods involve sophisticated redundant transforms, which
carry a heavy computational burden. Our proposed method is based on orthogo-
nal wavelet transform which is simple and fast but retains many useful properties.
It is interesting to see that even with this non-redundant transform, our method
is comparable to the state-of-the-art denoising methods under redundant wavelet
transform framework.
2. A new statistical model in wavelet domain is used to smooth the image and reduce
the noise prior to thresholding. We fully exploit the parent-child relationship be-
tween wavelet coefficients and inspect the neighboring dependency. We also give
a detailed discussion about how to construct the estimate for each coefficient by
modifying the number of relevant coefficients involved with respect to the wavelet
decomposition levels in a proper way.
3. In order to achieve optimal performances, most of the existing denoising algorithms
require the optimization of several parameters by solving a nonlinear system, whereas
the parameters in our method are obtained in a practical manner from the noise
model. We do not to need to solve a complicated optimization equations to get
these parameters, but derive a modification of a closed-form solution based on an
unbiased estimate of MSE, which saves tremendous computing time.
5 4. The core principle driving us to develop the proposed denoising method is to min-
imize the mean-squared error (MSE). From a practical standpoint, although this
criterion cannot be optimized in real applications due to the unavailability of the
noise-free images, both our method for additive Gaussian denoising and its version
for Poisson denoising do not require any prior information about the noise-free im-
ages. We also do not assume a statistical modeling for the observed image unlike
many Bayesian based methods.
5. We extend the denoising method to perform wavelet overcomplete expansion, which
yields even fewer visual artifacts and better image quality. Although extra compu-
tational time it costs, if visual effect is the main concern in a specific application,
this expansion is worth implementing.
1.5 Dissertation Organization
This dissertation mainly consists of two parts.
The first part of this dissertation generally deals with the problem of the noise removal of additive white Gaussian noise (AWGN) on 2D images.
In Chapter 2, we investigate the Gaussian noise model with its important features and properties. Common measurements applied to evaluate the performance and the efficiency of denoising approaches are included. Both objective and subjective image quality assessments are discussed.
In Chapter 3, we provide an in-depth and comprehensive literature review of the related work on Gaussian image denoising. Several influential approaches are described and compared to each other in details.
In Chapter 4, first, we review the motivation of reducing the upper bound of soft- thresholding scheme so as to achieve largest denoising effect, which leads to our two- operation denoising method. Then, a low-complexity, but remarkably efficient denois-
6 ing method is provided. This method incorporates an improved context modeling into the optimization of parameterized thresholding functions in the form of derivatives of Gaussian (DOG). We also present the efficiency of each denoising operation with numerical and visual results and combine them together as a two-operation process by comparing with a start-of-the-art image denoising method–SURE-LET [1] under non-overcomplete wavelet transform.
In Chapter 5, we extend the approach to the overcomplete representation of the wavelet transform by using cycle spinning strategy [2] in order to suppress the Pseudo-
Gibbs phenomena introduced by the standard orthogonal wavelet transform. After applying cycle spinning strategy, we compare our method to the redundant versions of the state-of-the-art wavelet techniques such as [3, 4] in terms of Peak Signal-to-Noise
Ratio (PSNR), respectively, as well as a non-wavelet method (BM3D) [5].
The second part of this dissertation is basically dedicated to the problem of non-
Gaussian noise, in particular Poisson noise, reduction on 2D images.
In Chapter 6, we discuss the necessity and the rationale of modeling low intensity images as Poisson distributed and give the mathematical model for Poisson noise. Af- ter this, an extensive literature survey of Poisson denoising research work is included.
At the end of this chapter, the most widely used variance stabilizing transformation
(VST)–Anscombe transformation and its inversions are described in details.
In Chapter 7, we apply the Anscombe transformation to reduce the Poisson sig- nal dependence under the VST framework so that our previously developed denoising method for AWGN can be directly used. Moreover, to avoid the biased errors gen- erated by its conventional inversions, we devise a new inversion aiming to correcting these errors. This new inversion is based on a piecewise polynomial regression model in the sense of weighted least squares.
In Chapter 8, we compare our proposed new inverse transformation with two con- ventional inverse transformations. We also evaluate its performance combined with our
7 denoising method for AWGN proposed in Part I by comparing to some of the leading algorithms such as [1,6,7]. Numerical analysis and denoised images are presented and indicate it as competitive as some state-of-the-art denoisers.
Finally in Chapter 9, the conclusion and some perspectives are presented.
8 Part I
Image Denoising for Additive
Gaussian Noise
9 Chapter 2
Background
It is a challenging task to suppress the disturbing noise while preserving important features on images at the same time. In order to guarantee a successful denoising operation, an appropriate noise model should be set up prior to any further develop- ment. Hence, in this chapter, we investigate the commonly used Gaussian noise model with its basic characteristics. In addition, to assess the effects of a specific denoising method and compare it with other approaches, we introduce several popular criteria for performance evaluation.
2.1 Additive Gaussian Noise Model
As it is well known, the pixel intensity value at a specific location on an image is highly related to the number of photons obtained by the corresponding captor during a fixed period time. According to the central limit theorem, when the light source is stable, the number of photons received by a single captor fluctuates around its average. Thus, in many real applications, especially when the magnitude of observed signal is relevantly high, the the noise n(i, j) can be reasonably assumed as independent and Gaussian distributed with mean and a standard deviation σ. By assuming independent, it
10 simply means pixel intensity values at different locations should be independent random variables. The 1D probability density function (PDF) of the Gaussian distribution is defined by 1 (x−μ)2 G(x)=√ e 2σ2 (2.1) 2πσ
It can be also modeled as additive which means that each pixel intensity value in the noisy image is the sum of the underlying true intensity value and a random noise value following a Gaussian distribution. One can write
y(i)=x(i)+n(i) (2.2) where y(i) is the observed value at pixel location i, x(i) would be the true value which indicates that this value would be obtained if one averages each observed value throughout a very long period of time. Besides, n(i) is the noise.
2.2 Image Quality Evaluation
Choose a proper way to evaluate the performance of denoising approaches is crucial since it is directly connected to the validation of a denoising algorithm applied on real applications. Generally, there are two popular approaches to image quality assessment.
We discuss them in the following.
2.2.1 Objective Image Quality Evaluation
The most intuitive criterion is to calculate the standard deviation of the noise on de- noised image, evaluate it in accordance with image’s smoothness. Other than this, there are mainly five objective quality assessment measures.
• Mean squared error (MSE)
The mean squared error is widely used for quality measurement. It is a mathemat-
11 ical performance index that measures the similarity between the denoised image xˆ and the reference: original noiseless image x. Its expression is defined by