图像处理-FFT

最新推荐文章于 2024-07-31 00:30:00 发布

fzk118gd147

最新推荐文章于 2024-07-31 00:30:00 发布

阅读量1k

点赞数

分类专栏：图像处理（机器视觉）

本文链接：https://blog.csdn.net/fzk118gd147/article/details/71525474

版权

图像处理（机器视觉）专栏收录该内容

2 篇文章 0 订阅

订阅专栏

im = im2double( rgb2gray( imread('3.png')));% “im” should be a gray-scale floating point image
[imh, imw] = size(im);
fftsize = 1024; % should be order of 2 (for speed) and include padding
im_fft = fft2(im, fftsize, fftsize); % 1) fft im with padding
hs = 30; % filter half-size
fil = fspecial('gaussian', hs*2+1, 10);
fil_fft = fft2(fil, fftsize, fftsize); % 2) fft fil, pad to same size as image
im_fil_fft = im_fft .* fil_fft; % 3) multiply fft images
im_fil = ifft2(im_fil_fft); % 4) inverse fft2
im_fil = im_fil(1+hs:size(im,1)+hs, 1+hs:size(im, 2)+hs); % 5) remove padding

subplot(1,2,1)
imshow(im)

subplot(1,2,2)
imshow(im_fil)

figure(2)
subplot(1,2,1)
imagesc((abs(fftshift(fil_fft))))

subplot(1,2,2)

imagesc(log(abs(fftshift(im_fft)))), axis image, colormap jet