Math 3360 Mathematical Imaging Lecture 11 Types of

Math 3360: Mathematical Imaging Lecture 11: Types of noises Prof. Ronald Lok Ming Lui Department of Mathematics, The Chinese University of Hong Kong

Class schedules Lecture 1: Introduction to Image Processing Lecture 2: Basic idea of image transformation Lecture 3: Image decomposition & Stacking operator Lecture 4: Singular Value Decomposition for Image decomposition & Error analysis Lecture 5: Haar & Walsh Transform Lecture 6: Examples of Haar & Walsh Transform; R-Walsh transform Lecture 7: Discrete Fourier transform Lecture 8: Even Discrete Cosine Transform (JPEG) Lecture 9: EDCT + ODCT+ EDST + ODST; Introduction to Image enhancement Lecture 10: Introduction to Linear filtering & Statistical images Lecture 11: Image denoising: Linear filtering model in the spatial domain; Image denoising: Nonlinear filtering model in the spatial domain; Relationship with the convolution Lecture 12: Image denoising: Linear filtering in the frequency domain Image denoising: Anisotropic diffusion Lecture 13: Image denoising: Total variation (TV) or ROF model Lecture 14: Image denoising: ROF model part 2 Lecture 15 to Lecture 17: Image deblurring Lecture 18 to Lecture 21: Image segmentation Lecture 22 to Lecture 24: Image registration

Type of noises n Recap: Preliminary statistical knowledge: n n n Random variables; Random field; Probability density function; Expected value/Standard deviation; Joint Probability density function; Linear independence; Uncorrelated; Covariance; Autocorrelation; Cross-correlation; Cross covariance; Noise as random field etc… Please refer to Supplemental note 6 for details.

Type of noises n Impulse noise: n n Change value of an image pixel at random; The randomness follows the Poisson distribution = Probability of having pixels affected by the noise in a window of certain size Poisson distribution: Gaussian noise: n Noise at each pixel follows the Gaussian probability density function:

Type of noises n Additive noise: n n Multiplicative noise: n n Noise parameter for the probability density function at each pixel are the same (same mean and same standard derivation) Zero-mean noise: n n Noisy image = original (clean) image * noise Homogenous noise: n n Noisy image = original (clean) image + noise Mean at each pixel = 0 Biased noise: n Mean at some pixels are not zero

Type of noises n Independent noise: n n Uncorrected noise: n n n Let Xi = noise at pixel i (as random variable); E(Xi Xj) = E(Xi) E(Xj) for all i and j. White noise: n n The noise at each pixel (as random variables) are linearly independent Zero mean + Uncorrelated + additive idd noise: n n n Independent + identically distributed; Noise component at every pixel follows the SAME probability density function (identically distributed) For Gaussian distribution,

Gaussian noise n Example of Gaussian noises:

White noise n Example of white noises:

Image components

Noises as high frequency component Why noises are often considered as high frequency component? (a) Clean image spectrum and Noise spectrum (Noise dominates the high-frequency component); (b) Filtering of high-frequency component

Linear filter = Convolution n n Linear filtering of a (2 M+1)x(2 N+1) image I (defined on [-M, M]x[-N, N]) = CONVOLUTION OF I and H H is called the filter. Different filter can be used: n Mean filter n Gaussian filter n Laplcian filter Variation of these filters (Non-linear) n Median filter n Edge preserving mean filter

Linear filter

Type of filter

In Photoshop

Mean filter

Mean filter Impulse noise After mean filter

Mean filter Gaussian noise After mean filter

Mean filter Real image After mean filter

Gaussian filter Define a function using Gaussian function Definition of H

Gaussian filter Real image After mean filter

Gaussian filter Real image After mean filter

Gaussian filter Real image After Gaussian filter

Gaussian filter Real image After mean filter

Gaussian filter Real image After Gaussian filter

Laplace filter ) Laplace filter (High pass filter)

Laplace filter Original Laplace filter

Laplace filter Original Laplace filter

Laplace filter Original Laplace filter

Median filter n Median n n Nonlinear filter Take median within a local window

Median filter Real image After mean filter

Median filter Salt & Pepper Mean filter Median filter

Median filter Noisy image Median filter

Median filter

Median filter Noisy image Median filter

Median filter

Median filter Noisy image Can you guess what it is? Median filter