pkg load image

clear
clc

I = imread('Lena256color.jpg');
I2 = I;

for i = 1 : 1300
  myRandRow = randi(256);
  myRandCol = randi(256);
  I2(myRandRow,myRandCol,1) = 255;
  I2(myRandRow,myRandCol,2) = 255;
  I2(myRandRow,myRandCol,3) = 255;
endfor
for i = 1 : 1300
  myRandRow = randi(256);
  myRandCol = randi(256);
  I2(myRandRow,myRandCol,1) = 0;
  I2(myRandRow,myRandCol,2) = 0;
  I2(myRandRow,myRandCol,3) = 0;
endfor
for i = 1 : 1300
  myRandRow = randi(256);
  myRandCol = randi(256);
  I2(myRandRow,myRandCol,1) = 255;
  I2(myRandRow,myRandCol,2) = 0;
  I2(myRandRow,myRandCol,3) = 0;
endfor
for i = 1 : 1300
  myRandRow = randi(256);
  myRandCol = randi(256);
  I2(myRandRow,myRandCol,1) = 0;
  I2(myRandRow,myRandCol,2) = 255;
  I2(myRandRow,myRandCol,3) = 0;
endfor
for i = 1 : 1300
  myRandRow = randi(256);
  myRandCol = randi(256);
  I2(myRandRow,myRandCol,1) = 0;
  I2(myRandRow,myRandCol,2) = 0;
  I2(myRandRow,myRandCol,3) = 255;
endfor

subplot(1,2,1), imshow(I, [0 255]), title('Original')
subplot(1,2,2), imshow(I2, []), title('Noisy')