clear; close all;clc;
img = imread('flower.tif');
img = imresize(img,[511 511]);
img = double(img);
img_fft = fftshift(fft2(img));
figure;mesh(real(img_fft));
filter = [1 1 1; 1 -8 1; 1 1 1];
filter_fft_pad = padarray(filter,[(size(img,1)-3)/2 (size(img,1)-3)/2]);
filter_fft = fftshift(fft2(filter_fft_pad));
figure;mesh(real(filter_fft_pad));
rt = fftshift(ifft2(fftshift(filter_fft.*img_fft)));
figure;
imshow(uint8(real(rt)));
%filter = rand(3,3);
% rt is almost same with img_conv
img_conv = conv2(img,filter,'same');
figure;
imshow(uint8(img_conv));
空域的卷积等于频率域的乘积例子
最新推荐文章于 2022-07-02 17:14:58 发布