去除蜂窝状的噪声（matlab实现）

如下图所示，图片中有蜂窝状的噪声，影响视野，可以通过频域滤波进行消除。

%% 读取图像
i = rgb2gray(imread('图1.jpg'));
figure(1)
subplot(221),imshow(i,[]),title('原图')
i_fft = fftshift(fft2(i));
subplot(222),imshow(abs(log10(abs(i_fft)+1)),[]),title('原图对应频谱')
%% 生成频域高斯低通滤波
x = linspace(-23,24,480);
y = linspace(-31,32,640);
[X,Y] = meshgrid(y, x);
D=sqrt(X.^2+Y.^2);
D0=0.1*35;
H=exp(-(D.^2)/(2*(D0)^2));
%% 滤波结果
res_fft = i_fft.*H;
subplot(223),imshow(abs(log10(abs(res_fft)+1)),[]),title('滤波后频谱')
res = ifft2(fftshift(res_fft));
subplot(224),imshow(res,[]),title('滤波后结果')