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
| close all;
clear;
clc;
[x,y] = meshgrid(single(-5:0.01:5));
z = peaks(x,y)*1000;
% figure;
% imagesc(z);colorbar;colormap(parula(256));
z2 = imfilter(z, fspecial('gaussian', [601,601], 101));
z3 = real(fftshift(ifft2(fft2(padarray(fspecial('gaussian', [601,601], 101),[200,200])) .* fft2(z))));
z4 = real(fftshift(ifftn(fftn(fspecial('gaussian', [601,601], 101), size(z)) .* fftn(z, size(z)), size(z))));
% z2 = imfilter(z, ones(201)/201^2);
% z3 = real(fftshift(ifft2(fft2(padarray(ones(201)/201^2,[400,400])) .* fft2(z))));
% z2 = imfilter(z, fspecial('motion', 283, 45));
% z3 = real(fftshift(ifft2(fft2(padarray(fspecial('motion', 283, 45),[400,400])) .* fft2(z))));
figure;
imagesc(z2);colorbar;colormap(parula(256));
figure;
imagesc(z3);colorbar;colormap(parula(256));
figure;
imagesc(z4);colorbar;colormap(parula(256));
% figure;
% imagesc(z3-z2);colorbar;colormap(parula(256)); |
Partager