matlab理想低通滤波

以下是Matlab实现理想低通滤波的代码和步骤：

1.读取图像并加入噪声

I = rgb2gray(imread('fig.png'));
IA = imnoise(I,'gaussian'); % 加入高斯白噪声

2.进行理想低通滤波

[f1,f2] = freqspace(size(IA),'meshgrid');
D = 0.05; % 我们选择截止频率为0.05
r = f1.^2 + f2.^2;
Hd = double(r <= D^2); % 构造理想低通滤波器
Y = fft2(double(IA));
Y = fftshift(Y);
Ya = Y .* Hd;
Ya = ifftshift(Ya);
Ia = real(ifft2(Ya));

3.显示图像

subplot(221),imshow(uint8(I)); title('原图像');
subplot(222),imshow(IA); title('加入噪声后的图像');
subplot(223),imshow(uint8(Ia)); title('理想低通滤波处理');