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

整理文档很辛苦,赏杯茶钱您下走!

免费阅读已结束,点击下载阅读编辑剩下 ...

阅读已结束,您可以下载文档离线阅读编辑

资源描述

高通滤波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;endy(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=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('低通滤波');

1 / 2
下载文档,编辑使用

©2015-2020 m.777doc.com 三七文档.

备案号:鲁ICP备2024069028号-1 客服联系 QQ:2149211541

×
保存成功