




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 日文版劳动合同
- 游戏开发引擎授权使用及技术支持服务协议
- 三字经学习感悟
- Oleic-acid-Standard-生命科学试剂-MCE
- 小学生经典国学启蒙读后感
- 2-2-Dibutylpropane-1-3-diol-生命科学试剂-MCE
- 年终团队建设活动安排
- 年度销售部门工作总结报告新趋势与成效分析
- 公司融资计划及投资策略分析
- 试用期员工劳动合同
- 原发性血小板减少性紫癜患者的生活质量
- HG20202-2014 脱脂工程施工及验收规范
- 基层医疗机构公共卫生健康教育工作手册(各类表格模板)
- 2024版《保密法》培训课件
- 2024年内蒙古中考地理生物试卷(含答案)
- 广东省汕尾市汕尾市2024年中考一模英语试题(含答案)
- LY/T 3378-2024木蜡油地板
- 元宵节猜灯谜 11
- 施工现场视频监控系统施工方案
- 2024年江西电力职业技术学院单招职业适应性测试题库含答案
- (正式版)JTT 1495-2024 公路水运危险性较大工程安全专项施工方案审查规程
评论
0/150
提交评论