i=rgb2gray(imread('test.jpg'));%%put your own test image in this folder
i=im2double(i);
gpre=i;
figure, imshow(mat2gray(gpre));
title('Gaussian pyramid');
figure, imshow(abs(fftshift(fft2(im2uint8(gpre)))),[20, 50000]),colormap jet;
title('Gaussian FFT Amplitude');
for index=1:4
gnew=gpre;
gnew=imgaussfilt(gnew,2);
gnew=downsample(downsample(gnew,2).',2).';
figure, imshow(mat2gray(gnew));
title('Gaussian pyramid');
figure, imshow(abs(fftshift(fft2(im2uint8(gnew)))),[20, 50000]),colormap jet;
title('Gaussian FFT Amplitude');
lappre=imresize(gnew,2);
lap=gpre-lappre;
figure, imshow(mat2gray(lap));
title('Laplace pyramid');
figure, imshow(abs(fftshift(fft2(im2uint8(lap)))),[20, 50000]),colormap jet;
title('Laplacian FFT Amplitude');
gpre=gnew;
end
评论0