Adding noise into an image manually instead of using imnoise

48 次查看（过去 30 天）
Sufyan
Sufyan 2012-8-30
评论： Mintesnot 2022-9-29
I am to trying to understand the algorithms behind matlab way of adding noise into an image,
The algorithm which Matlab use to add Gaussian noise is this,
b = a + sqrt(p4)*randn(sizeA) + p3;
When I tried to implement this algorithm manually it worked successfully however it doesn't work unless i changed the image class to double. Why is that so ? Why should i suppose to change the class to double when adding gaussian noise ?
Here is my code,
I = imread('2.jpg');
J = rgb2gray(I);
p3= 0;
p4 = 0.05;
J = im2double(J);
b = J + sqrt(p4)*randn(size(J)) + p3;
imshow(b)
Just like Gaussian Noise i tried adding the Salt n Pepper noise manually, here is the algorithm Matlab use to add Salt n Pepper noise,
b = a; <-- Assign b to the input image
x = rand(sizeA); <--- Generate random pixels from the image pixels
d = find(x < p3/2); <--- Find the pixels whose values are less than half of the mean value
b(d) = 0 <-- Implement minimum saturation to them
d = find(x >= p3/2 & x < p3) <--- Find the pixels whose values are greater than half of the mean value & less than mean value
b(d) = 1; <-- Implement maximum saturation to them
and I implemented it as below,
I = imread('2.jpg');
J = rgb2gray(I);
p3= 0.5;
x = rand(size(J));
d = find(x < p3/2);
J(d) = 0; % Minimum value
d = find(x >= p3/2 & x < p3);
J(d) = 1; % Maximum (saturated) value
imshow(J)
but the output Image doesn't show any Salt n pepper noise in the Image , I wonder where the final image is actually stored in my code?
  3 个评论
显示 1更早的评论
sravani honey
sravani honey 2018-4-20
sir, i am doing my project on random impulse noise removal and i am adding random noise using noisyimage = (I + p*rand(size(I)))/(1+p) this function but it doesn't clear results on my project. can u please give any other command related to adding random impulse noise to the image?

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Walter Roberson
Walter Roberson 2012-8-30
Images read by imread() are often data class uint8(). The uint8() property will be retained in the rgb2gray() step. Your saturation assignments of 0 and 1 are, however, based upon the notion that the minimum and maximum are 0 and 1 instead of 0 and 255.
Also, the minimum value that a uint8 variable can be is 0, so there is no need to explicitly set 0 as the minimum.
  3 个评论
显示 1更早的评论
Walter Roberson
Walter Roberson 2012-8-31
MATLAB uses different value ranges for the different data classes of images. Images that are represented as uint8(), use the value range 0 to 255. Images that are represented as floating point, use the value range 0 to 1. Your difficulty is that you are trying to apply an algorithm that depends upon the 0 to 1 value range, to data that has been represented with the 0 to 255 value range.
It is not because of grayscale: you can have grayscale in uint8 or in floating point.
Please look at the routine im2double()

请先登录，再进行评论。

更多回答（2 个）

Image Analyst
Image Analyst 2012-8-31
Gaussian noise can have any magnitude. So how are you going to add 4.28435435 to an image that can have only integer values? If you could and then it was converted to an integer, you'd lose the fractional part of the number, which means it would no longer be Gaussian. So that's why you have to convert the image to floating point.
  3 个评论
显示 1更早的评论
Image Analyst
Image Analyst 2012-8-31
No. Gaussian noise can be any floating point value. However if you force it to be in a uint8 image you are not letting it be any value, and so it's not Gaussian noise anymore. So to prevent that, you have to convert your image to floating point (single or double).

请先登录，再进行评论。

Ameerah Omar
Ameerah Omar 2015-11-21
how i can add noise in image it come from Cauchy distribution???
  1 个评论
Image Analyst
Image Analyst 2015-11-21
Use the formula for the CDF and then use rand() to draw numbers from the distribution. Attached is an example where I do that for the Rayleigh distribution.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Images 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by