指数高通滤波器代码_高通滤波器matlab程序代码

%高斯低通滤波器

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

滤波函

数