lab 9 - poc

%1
close all;
clear all;

IMG =imread('dwieFale.bmp')

IMG_f= fft2(IMG)
IMG_optyczna = fftshift(IMG_f)

IMG_abs = abs(IMG_optyczna)
IMG_abs = log10(IMG_abs+1)
IMG_faza = angle(IMG_f.*(IMG_abs>0.0001))

IMG_F1 = fft(IMG,[],1)
IMG_F2 = fft(IMG_F1,[],1)
IMG_F3 = fft(IMG,[],2)
IMG_F4 = fft(IMG_F3,[],2)

IMG_w = IMG_F1 + IMG_F2;
IMG_w1 = IMG_F1 + IMG_F2;


F1 = fft(IMG,[],1)
F2 = fft(F1,[],2)

F3 = fft(IMG,[],2)
F4 = fft(F3,[],1)

IMG_roznica = F4-IMG_f
IMG_abs2 = abs(IMG_roznica)

IMG_roznica2 = F2-IMG_f
IMG_abs3 = abs(IMG_roznica2)

figure(1)
subplot(2,2,1)
imshow(IMG)
title('Obraz originalny')
subplot(2,2,2)
imshow(IMG_abs,[])
title('Logartym')
subplot(2,2,3)
imshow(IMG_faza,[])
title('Faza')
subplot(2,2,4)
imshow(IMG_w1,[])
title('Suma')

figure(2)
subplot(2,3,1)
imshow(F2)
title('fft2 - dwu')
subplot(2,3,2)
imshow(IMG_f)
title('fft - jedno')
subplot(2,3,3)
imshow(IMG_roznica,[])
title('Roznica')
subplot(2,3,4)
imshow(IMG_roznica2,[])
title('zmienona kolejnosc')
%2
close all;
clear all;

IMG =imread('dwieFale.bmp')
kwadrat = imread('kwadrat.bmp')

translacjaT = imread('kwadratT.bmp')
rotacja45 = imread('kwadrat45.bmp')
zmiana_rozmiaruS = imread('kwadratS.bmp')
kombinacja_liniowa45 = imread('kwadrat45.bmp')
kombinacja_liniowaKL = imread('kwadratKL.bmp')

translacja = imtranslate(kwadrat,[15,25])
translacjaT = imtranslate(translacjaT,[15,25])

rotacja = imrotate(kwadrat,15)
rotacja45 = imrotate(rotacja45,15)

zmiana_rozmiaru = imresize(kwadrat,0.5)
zmiana_rozmiaruS = imresize(zmiana_rozmiaruS,0.5)

kombinacja_liniowa = imlincomb(1.5,kwadrat);
kombinacja_liniowa45 = imlincomb(1.5,kombinacja_liniowa45);
kombinacja_liniowaKL = imlincomb(1.5,kombinacja_liniowaKL);


IMG_f= fft2(kwadrat)
IMG_op = fftshift(IMG_f)
IMG_abs = abs(IMG_op)
IMG_abs = log10(IMG_abs+1)
IMG_faza = angle(IMG_f.*(IMG_abs>0.0001))

IMG_f2= fft2(translacja)
IMG_op2 = fftshift(IMG_f2)
IMG_abs2 = abs(IMG_op2)
IMG_abs2 = log10(IMG_abs2+1)
IMG_faza2 = angle(IMG_f2.*(IMG_abs2>0.0001))

IMG_f3= fft2(rotacja)
IMG_op3 = fftshift(IMG_f3)
IMG_abs3 = abs(IMG_op3)
IMG_abs3 = log10(IMG_abs3+1)
IMG_faza3 = angle(IMG_f3.*(IMG_abs3>0.0001))

IMG_f4= fft2(zmiana_rozmiaru)
IMG_op4 = fftshift(IMG_f4)
IMG_abs4 = abs(IMG_op4)
IMG_abs4 = log10(IMG_abs4+1)
IMG_faza4 = angle(IMG_f4.*(IMG_abs4>0.0001))

IMG_f5= fft2(kombinacja_liniowa)
IMG_op5 = fftshift(IMG_f5)
IMG_abs5 = abs(IMG_op5)
IMG_abs5 = log10(IMG_abs5+1)
IMG_faza5 = angle(IMG_f5.*(IMG_abs5>0.0001))

figure(1)
subplot(2,3,1)
imshow(IMG_faza)

subplot(2,3,2)
imshow(IMG_faza2)

subplot(2,3,3)
imshow(IMG_faza3)

subplot(2,3,4)
imshow(IMG_faza4)

subplot(2,3,5)
imshow(IMG_faza5)

figure(2)
subplot(2,3,1)
imshow(IMG_abs)

subplot(2,3,2)
imshow(IMG_abs2)

subplot(2,3,3)
imshow(IMG_abs3)

subplot(2,3,4)
imshow(IMG_abs4)

subplot(2,3,5)
imshow(IMG_abs5)

%3
close all;
clear all;

IMG_kol =imread('kolo.bmp')

IMG_f= fft2(IMG_kol)
IMG_prz = ifftshift(IMG_f)
IMG_op = ifft2(IMG_prz)

IMG1 = abs(IMG_kol)
IMG2 = uint8(abs(IMG_op))
IMG3= IMG1+IMG2

figure(1)
subplot(1,3,1)
imshow(IMG_kol)

subplot(1,3,2)
imshow(IMG_op)

subplot(1,3,3)
imshow(IMG3)