1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55
| function GFD_test(I)
img=imread(I);
imshow(log(abs(fftshift(fft2(rgb2gray(img))))), []);
figure
subplot(1,3,1);
image(img),title('image origine'),grid
axis image off;
[img2,e,n]=logpolar(rgb2gray(img));
subplot(1,3,2)
image(img2),title('domaine polaire'),grid
axis image xy off
%pha=angle(fft_a2);
F = fft2(img2);
%subplot(1,4,3)
%imshow(real(F))
%axis image off
F=fftshift(F);
subplot(1,3,3)
tal1=size(F);
d1=max(F);
for j=1:tal1(1);
for k=1:tal1(2);
F(j,k)=(F(j,k)/d1(1))*255;
end;
end;
%imshow(uint8(F)),title('Spectrum shape'),grid;
imshow(log(abs(F)), []);
axis image xy off
%F=abs(fft2(img2));%imshow(real(fftshift(F)));
R=[];
for i=1:size(F)
R=[R F(i,:)];
end
R = R(:);
R=R/sum(R);
var = inputdlg('Entrez un nom pour la sauvegarde','Sauvegarde de descripteur de l''image.',1,{''});
assignin('base', [var{1}],R);
save ( var{1},'R')
end |
Partager