高通滤波-低通滤波matlab程序.doc

《高通滤波-低通滤波matlab程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、高通滤波matlab程序：clearall;I1=imread('E:matlab图片硬币.jpg');subplot(1,2,1);imshow(I1);title('原始图像');m=double(I1);f=fft2(m);f=fftshift(f);[N1,N2]=size(f);%返回矩阵的行和列n1=round(N1/2);n2=round(N2/2);n=2;d0=50;%滤波器截止频率,滤波半径fori=1:N1forj=1:N2d=sqrt((i-n1)^2+(j-n2)^2);%计算高通滤波转换函数ifd<=d0h=0+0.5;elseh=1+0.5;en

2、dy(i,j)=h*f(i,j);endendy=ifftshift(y);A=ifft2(y);B=uint8(real(A));subplot(1,2,2);imshow(B);title('高通滤波');低通滤波matlab程序：clearall;I1=imread('E:matlab图片硬币.jpg');subplot(1,2,1);imshow(I1);title('原始图像');m=double(I1);f=fft2(m);f=fftshift(f);[N1,N2]=size(f);%返回矩阵的行和列n1=round(N1/2);n2=round(N2/2);n=

3、2;d0=50;%滤波器截止频率,滤波半径fori=1:N1forj=1:N2d=sqrt((i-n1)^2+(j-n2)^2);%计算低通滤波转换函数ifd<=d0h=1;elseh=0;endy(i,j)=h*f(i,j);endendy=ifftshift(y);A=ifft2(y);B=uint8(real(A));subplot(1,2,2);imshow(B);title('低通滤波');