




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、精品文档 各类滤波器的MATLAB程序 一、理想低通滤波器 IA=imread(le na.bmp); f1,f2=freqspace(size(IA),meshgrid); Hd=on es(size(IA); r=sqrt(f1.A2+f2.A2); Hd(r0.2)=0; Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure s
2、urf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); o o 二、理想高通滤波器 IA=imread(le na.bmp); f1,f2=freqspace(size(IA),meshgrid); Hd=on es(size(IA); r=sqrt(f1.A2+f2.A2); Hd(r0.2)=0; Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); la=real(ifft2(Ya); figure subplot(2,2,1),imshow(
3、ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); o o 三、Butterworth低通滤波器 IA=imread(le na.bmp); f1,f2=freqspace(size(IA),meshgrid); D=0.3; r=f1.A2+f2.A2; n=4; for i=1:size(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=1/(tA n+1
4、); end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); 0 0 四、Butterworth高通滤波器 IA=imread(le na.bmp); f1,f2=freqspa
5、ce(size(IA),meshgrid); D=0.3; r=f1.A2+f2.A2; n=4; for i=1:size(IA,1) for j=1:size(IA,2) t=(D*D)/r(i,j); Hd(i,j)=1/(tA n+1); end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf
6、(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); 0 0 五、高斯低通滤波器 IA=imread(le na.bmp);IB=imread(babarra.bmp); f1,f2=freqspace(size(IA),meshgrid); D=100/size(IA,1); r=f1.A2+f2.A2; Hd=on es(size(IA); for i=1:size(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=exp(-t); end end Y=fft2(doub
7、le(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); 0 0 六、高斯高通滤波器 IA=imread(le na.bmp);IB=imread(babarra.bmp); f1,f2=freqspace(size(
8、IA),meshgrid); %D=100/size(IA,1); D=0.3; r=f1.A2+f2.A2; for i=1:size(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=1-exp(-t); end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure
9、 surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); 0 0 七、梯形低通滤波器 IA=imread(le na.bmp);IB=imread(babarra.bmp); f1,f2=freqspace(size(IA),meshgrid); %D=100/size(IA,1); D0=0.1; D1=0.4; r=sqrt(f1.A2+f2.A2); Hd=zeros(size(IA); Hd(r=D0 end end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd
10、; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight in g,pho ng); 精品文档 0 0 八、梯形高通滤波器 IA=imread(le na.bmp);IB=imread(babarra.bmp); f1,f2=freqspace(size(IA),meshgrid); %D=100/size
11、(IA,1); D0=0.1; D1=0.4; r=sqrt(f1.A2+f2.A2); Hd=on es(size(IA); Hd(r=D0 end end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(ui nt8(IA); subplot(2,2,2),imshow(u in t8(Ia); figure surf(Hd,Facecolor,i nterp,Edgecolor, non e,Facelight
12、 in g,pho ng); 0 0 九、用其他方法编写的理想低通、理想高通、Butterworth低通、同态滤波程序 1、理想低通 i1=imread(le na.bmp); i2=i mno ise(i1,salt f=double(i2); k=fft2(f); g=fftshift(k); N1,N2=size(g); d0=50; uO=floor(N1/2)+1; vO=floor(N2/2)+1; for i=1:N1 for j=1:N2 d=sqrt(i-u0)A2+(j-v0)A2); if d=d0 h=1; else h=0; end y(i,j)=g(i,j)*h;
13、end end y=ifftshift(y); E1=ifft2(y); E2=real(E1); figure subplot(2,2,1),imshow(ui nt8(i1); subplot(2,2,2),imshow(ui nt8(i2); subplot(2,2,3),imshow(ui nt8(E2); 2、理想高通 i1=imread(le na.bmp); i2=i mno ise(i1,salt f=double(i2); k=fft2(f); g=fftshift(k); N1,N2=size(g); n=2; d0=10; uO=floor(N1/2)+1; vO=flo
14、or(N2/2)+1; for i=1:N1 for j=1:N2 d=sqrt(i-u0)A2+(j-v0)A2); if d=d0 h=0; else h=1; end y(i,j)=g(i,j)*h; end end y=ifftshift(y); E1=ifft2(y); E2=real(E1); figure subplot(2,2,1),imshow(ui nt8(i1); subplot(2,2,2),imshow(ui nt8(i2); subplot(2,2,3),imshow(ui nt8(E2); -.h Ji :“ .J1 3、Butterworth 低通 i1=imr
15、ead(le na.bmp); i2=i mno ise(i1,salt f=double(i2); k=fft2(f); g=fftshift(k); N1,N2=size(g); n=2; d0=50; uO=floor(N1/2)+1; vO=floor(N2/2)+1; for i=1:N1 for j=1:N2 d=sqrt(i-u0)A2+(j-v0)A2); h=1/(1+(d/d0)A(2* n); y(i,j)=g(i,j)*h; end end y=ifftshift(y); E1=ifft2(y); E2=real(E1); figure subplot(2,2,1),imshow(ui nt8(i1); subplot(2,2,2),imshow(ui nt8(i2); subplot(2,2,3),imshow(ui nt8(E2); 4、同态滤波 l=rgb2gray(imread(fabric00.bmp); M,N=size(l); T=double(l); L=l
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 商场内的商铺租赁合同
- 停车场承包合同
- 技术培训委托合同书
- 草籽草坪采购合同
- 房屋独家代理销售合同
- 车库转让合同协议书
- 医疗美容手术项目合同协议书
- 高层管理团队建设活动方案
- 上海餐饮商铺租赁合同
- 奶茶店转让合同(新标准版)8篇
- 第二十一章会阴部美容手术讲解
- 2025年度iPhone手机租赁与虚拟现实体验合同3篇
- 2025年度消防工程安全防护措施设计固定总价合同范本3篇
- 苏北四市(徐州、宿迁、淮安、连云港)2025届高三第一次调研考试(一模)语文试卷(含答案)
- 食品企业危机管理应对方案
- 2024年济南广播电视台招聘工作人员笔试真题
- 市场消防安全课件
- 名师工作室建设课件
- 2025-2025学年度人教版小学五年级美术下册教学计划
- 《电力建设工程施工安全管理导则》(NB∕T 10096-2018)
- 《黄色新闻的泛滥》课件
评论
0/150
提交评论