close all;
clear all;

kw = imread('kw.bmp');
load maskiPP.mat;

kwConv2R1 = conv2(kw, R1, 'same');
kwConv2R1_128 = kwConv2R1 + 128;
kwConv2R1Modul = abs(kwConv2R1);

figure(1)
subplot(1,3,1)
imshow(kw);
subplot(1,3,2)
imshow(kwConv2R1_128, []);
subplot(1,3,3);
imshow(kwConv2R1Modul, []);

kwConv2R2 = conv2(kw, R2, 'same');
kwConv2R2_128 = kwConv2R2 + 128;
kwConv2R2Modul = abs(kwConv2R2);

figure(2)
subplot(1,3,1)
imshow(kw);
subplot(1,3,2)
imshow(kwConv2R2_128, []);
subplot(1,3,3);
imshow(kwConv2R2Modul, []);


kwConv2P1 = conv2(kw, P1, 'same');
kwConv2P1_128 = kwConv2P1 + 128;
kwConv2P1Modul = abs(kwConv2P1);

figure(3)
subplot(1,3,1)
imshow(kw);
subplot(1,3,2)
imshow(kwConv2P1_128, []);
subplot(1,3,3);
imshow(kwConv2P1Modul, []);

kwConv2P2 = conv2(kw, P2, 'same');
kwConv2P2_128 = kwConv2P2 + 128;
kwConv2P2Modul = abs(kwConv2P2);

figure(4)
subplot(1,3,1)
imshow(kw);
subplot(1,3,2)
imshow(kwConv2P2_128, []);
subplot(1,3,3);
imshow(kwConv2P2Modul, []);

kwConv2S1 = conv2(kw, S1, 'same');
kwConv2S1_128 = kwConv2S1 + 128;
kwConv2S1Modul = abs(kwConv2S1);

figure(5)
subplot(1,3,1)
imshow(kw);
subplot(1,3,2)
imshow(kwConv2S1_128, []);
subplot(1,3,3);
imshow(kwConv2S1Modul, []);



kwConv2S2 = conv2(kw, S2, 'same');
kwConv2S2_128 = kwConv2S2 + 128;
kwConv2S2Modul = abs(kwConv2S2);

figure(6)
subplot(1,3,1)
imshow(kw);
subplot(1,3,2)
imshow(kwConv2S2_128, []);
subplot(1,3,3);
imshow(kwConv2S2Modul, []);


kwConv2FK = sqrt(kwConv2S1.^2 + kwConv2S2.^2);
kwConv2FK_128 = kwConv2FK + 128;
kwConv2FKModul = abs(kwConv2FK);

figure(7)
subplot(1,3,1)
imshow(kw);
subplot(1,3,2)
imshow(kwConv2FK_128, []);
subplot(1,3,3);
imshow(kwConv2FKModul, []);


kwConv2FK2 = abs(kwConv2S1) + abs(kwConv2S2);
kwConv2FK2_128 = kwConv2FK2 + 128;
kwConv2FK2Modul = abs(kwConv2FK2);

figure(8)
subplot(1,3,1)
imshow(kw);
subplot(1,3,2)
imshow(kwConv2FK2_128, []);
subplot(1,3,3);
imshow(kwConv2FK2Modul, []);

jet = imread('jet.bmp')
jetConv2_1 = conv2(jet,S1,'same');
jetConv2_2 = conv2(jet,S2,'same');
jetConv2FK = sqrt(jetConv2_1.^2 + jetConv2_2.^2);
jetConv2FK_128 = jetConv2FK + 128;
jetConv2FKabs = abs(jetConv2FK);


figure(9)
subplot(1,3,1)
imshow(jet);
subplot(1,3,2)
imshow(jetConv2FK_128, []);
subplot(1,3,3);
imshow(jetConv2FKabs, []);

jetConv2FK2 = abs(jetConv2_1) + abs(jetConv2_2);
jetConv2FK2_128 = jetConv2FK2 + 128;
jetConv2FK2abs = abs(jetConv2FK2);

figure(10)
subplot(1,3,1)
imshow(jet);
subplot(1,3,2)
imshow(jetConv2FK2_128, []);
subplot(1,3,3);
imshow(jetConv2FK2abs, []);