数字图像-频域处理

关注

原创最新推荐文章于 2026-04-25 09:51:38 发布 852 阅读 AI 写同款· GEO 优化 ›

数字图像处理第三章——频域处理（上） 数字图像处理第二章数字图像处理---频域处理（一）二维离散傅里叶变换4.1线性空间滤波4.1.1 相关4.1.2 卷积4.1.3 imfilter函数4.2 非线性空间滤波（五）图像处理工具箱中标准的空间滤波器5.1线性空间滤波器5.2非线性空间滤波器 数字图像处理—频域处理（一）二维离散傅里叶变换邻域处理包含以下过程: 选择中心点( x，y) ; 仅对预定义的围绕点(x，y)的邻域内的... 阅读详情

冈萨雷斯第三版，频域处理代码示例：

% 低通滤波
clc, clear,

I = imread('blown_ic.tif');
[r,c] = size(I);
subplot(331),imshow(I);title('Origin');

Id = padarray(I,[r,c],'post');
subplot(332),imshow(Id);title('pad array');

for i = 1:2*r
for j = 1:2*c
hp(i,j) = (-1)^(i-1+j-1);
end
end
Ip = hp.*double(Id);
subplot(333),imshow(Ip,[]);title('shift array');

F = fft2(Ip);% padding of H
subplot(334),imshow(log(1+abs(F)),[]);title('spectrum of F');

H = lpfilter('gaussian',2*r,2*c,50,1);%
subplot(335),imshow(H,[]);title('H');

G = H.*F;
subplot(336),imshow(G);title('G=H.*F ');

gp = hp.*real(ifft2(G));% filtering
subplot(337),imshow(gp,[]);title('gp ');

g = gp(1:r,1:c);% cropped images
subplot(338),imshow(g,[]);title('g ');