高通滤波器matlab程序代码

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

%高斯低通滤波器

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滤波函数

s(i,j)=h(i,j)*s(i,j);%GLPF滤波后的频域表示

end

end

s=ifftshift(s);%对s进行反FFT移动

%对s进行二维反离散的Fourier变换后,取复数的实部转化为无符号8位整数

s=uint8(real(ifft2(s)));

subplot(2,3,4),imshow(h);title('传递函数');%显示GHPF滤波器的传递函数

subplot(2,3,5),imshow(s);title('GLPF滤波(d0=50)');%显示GLPF滤波处理后的图像

相关文档
最新文档