matlab

clear;

zaklocenia=imread('zakИвcenia_1.jpg');
noise=imread('NOISE.TIF');

FL1 = [1 1 1; 1 1 1; 1 1 1];
FL2 = [1 1 1; 1 2 1; 1 1 1];
FL3 = [1 1 1; 1 4 1; 1 1 1];
Gauss = [1 2 1; 2 4 2; 1 2 1];

noiseD = (double(noise)/255);
zakloceniaD = (double(zaklocenia)/255);

FL1Z=colorFilter(zakloceniaD,FL1);
FL1N=filter2(FL1,noiseD)/(sum(sum(FL1)));
FL2Z=colorFilter(zakloceniaD,FL2);
FL2N=filter2(FL2,noiseD)/(sum(sum(FL2)));
FL3Z=colorFilter(zakloceniaD,FL3);
FL3N=filter2(FL3,noiseD)/(sum(sum(FL3)));
GaussZ=colorFilter(zakloceniaD,Gauss);
GaussN=filter2(Gauss,noiseD)/(sum(sum(Gauss)));

subplot(2,5,1);
imshow(zakloceniaD);
subplot(2,5,2);
imshow(noiseD);
subplot(2,5,3);
imshow(FL1Z);
subplot(2,5,4);
imshow(FL1N);
subplot(2,5,5);
imshow(FL2Z);
subplot(2,5,6);
imshow(FL2N);
subplot(2,5,7);
imshow(FL3Z);
subplot(2,5,8);
imshow(FL3N);
subplot(2,5,9);
imshow(GaussZ);
subplot(2,5,10);
imshow(GaussN);
function colorFilter = colorFilter(image,filter)
    filR=filter2(filter,image(:,:,1))/(sum(sum(filter)));
    filG=filter2(filter,image(:,:,2))/(sum(sum(filter)));
    filB=filter2(filter,image(:,:,3))/(sum(sum(filter)));
    colorFilter(:,:,1)=filR;
    colorFilter(:,:,2)=filG;
    colorFilter(:,:,3)=filB;
end