%高斯低通滤波器
RGB
=
imread('132.jpg');
I0
=
rgb2gray(RGB);
subplot(2,3,1),imshow(I0);title('原图');
I1
=
imnoise(I0,'gaussian');
%对原图像加噪声
subplot(2,3,2),imshow(I1);title('加入噪声后')
%将灰度图像的二维不连续
Fourier
变换的零频率成分移到频谱的中
心
s=fftshift(fft2(I1));
subplot(2,3,3),imshow(log(1+abs(s)),[]);title('fftshift');
[M,N]=size(s);
%分别返回
s
的行数到
M
中,
列数到
N
中
%GLPF
滤波
d0=50;
%初始化
d0
n1=floor(M/2);
%对
M/2
进行取整
n2=floor(N/2);
%对
N/2
进行取整
for
i=1:M
for
j=1:N
d=sqrt((i-n1)^2+(j-n2)^2);
%点
(i,j)
到傅立叶
变换中心的距离
h(i,j)=1*exp(-1/2*(d^2/d0^2));
%GLPF
滤波函
数