%**********************************************************
%**            Výpocet FFT spektrálnich obrazu           **
%**                Diskrétní transformace                **
%**        Leos Marsalek Leos.Marsalek@tiscali.cz        **
%**       Jan Skapa Jeniczek.s@seznam.cz (c) 2003        **
%**********************************************************
Gen_obr; % generovani obrazku, Ke stazeni na www.sendme.cz/goro
pic=pic1;           % vyber obrazku
ypic=rgb2ycbcr(pic);% prepocet do jineho barevneho prostoru.
obr=ypic(:,:,1);
OBR=fft2(obr);
velikost=size(obr);
Obr=OBR(1:velikost(1)/2,1:velikost(2)/2);
imwrite(abs(Obr)./100,'fft_ampl_pic1.tif','tif');
imwrite(angle(Obr)./2,'fft_phase_pic1.tif','tif');
imwrite(abs(OBR)./100,'fft_matlab_pic1.tif','tif');

pic=pic2;           % vyber obrazku
ypic=rgb2ycbcr(pic);% prepocet do jineho barevneho prostoru.
obr=ypic(:,:,1);
OBR=fft2(obr);
velikost=size(obr);
Obr=OBR(1:velikost(1)/2,1:velikost(2)/2);
imwrite(abs(Obr)./100,'fft_ampl_pic2.tif','tif');
imwrite(angle(Obr)./2,'fft_phase_pic2.tif','tif');
imwrite(abs(OBR)./100,'fft_matlab_pic2.tif','tif');

pic=pic3;           % vyber obrazku
ypic=rgb2ycbcr(pic);% prepocet do jineho barevneho prostoru.
obr=ypic(:,:,1);
OBR=fft2(obr);
velikost=size(obr);
Obr=OBR(1:velikost(1)/2,1:velikost(2)/2);
imwrite(abs(Obr)./100,'fft_ampl_pic3.tif','tif');
imwrite(angle(Obr)./2,'fft_phase_pic3.tif','tif');
imwrite(abs(OBR)./100,'fft_matlab_pic3.tif','tif');

pic=pic4;           % vyber obrazku
ypic=rgb2ycbcr(pic);% prepocet do jineho barevneho prostoru.
obr=ypic(:,:,1);
OBR=fft2(obr);
velikost=size(obr);
Obr=OBR(1:velikost(1)/2,1:velikost(2)/2);
imwrite(abs(Obr)./100,'fft_ampl_pic4.tif','tif');
imwrite(angle(Obr)./2,'fft_phase_pic4.tif','tif');
imwrite(abs(OBR)./100,'fft_matlab_pic4.tif','tif');

pic=pic5;           % vyber obrazku
ypic=rgb2ycbcr(pic);% prepocet do jineho barevneho prostoru.
obr=ypic(:,:,1);
OBR=fft2(obr);
velikost=size(obr);
Obr=OBR(1:velikost(1)/2,1:velikost(2)/2);
imwrite(abs(Obr)./100,'fft_ampl_pic5.tif','tif');
imwrite(angle(Obr)./2,'fft_phase_pic5.tif','tif');
imwrite(abs(OBR)./100,'fft_matlab_pic5.tif','tif');