全文预览已结束
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
实践二: 理想高通滤波器、Butterworth高通滤波器、高斯高通滤波器2.1.1 理想高通滤波器实践代码: I=imread(girl.bmp);subplot(221),imshow(I);title(原图像);s=fftshift(fft2(I);subplot(223),imshow(abs(s),);title(图像傅里叶变换所得频谱);subplot(224),imshow(log(abs(s),);title(图像傅里叶变换取对数所得频谱);a,b=size(s);a0=round(a/2);b0=round(b/2);d=10;p=0.2;q=0.5;for i=1:a for j=1:b distance=sqrt(i-a0)2+(j-b0)2); if distance=d h=0; else h=1; end; s(i,j)=(p+q*h)*s(i,j); end;end;s=uint8(real(ifft2(ifftshift(s);subplot(222),imshow(s);title(高通滤波所得图像);I=imread(girl.bmp);f1,f2=freqspace(size(I),meshgrid);Hd=ones(size(I);r=sqrt(f1.2+f2.2);Hd(r0.2)=0;figuresurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong); % 画三维曲面(色)图2.1.2 理想高通滤波器实践结果截图:2.2.1 Butterworth高通滤波器实践代码:I1=imread(flower.bmp);subplot(121),imshow(I1);title(原始图像);f=double(I1);g=fft2(f);g=fftshift(g);N1,N2=size(g);n=2;d0=5;n1=fix(N1/2);n2=fix(N2/2);for i=1:N1 for j=1:N2 d=sqrt(i-n1)2+(j-n2)2); if d=0 h=0; else h=1/(1+(d0/d)(2*n); end result(i,j)=h*g(i,j); endendresult=ifftshift(result);X2=ifft2(result);X3=uint8(real(X2);subplot(122),imshow(X3) ;title(Butterworth高通滤波);I1=imread(flower.bmp);f1,f2=freqspace(size(I1),meshgrid);D=0.3;r=f1.2+f2.2;n=4;for i=1:size(I1,1) for j=1:size(I1,2) t=(D*D)/r(i,j); Hd(i,j)=1/(tn+1); endendfiguresurf(Hd,Facecolor,interp,Edgecolor,none,Facelighting,phong); % 画三维曲面(色)图2.2.2 Butterworth高通滤波器实践结果截图:2.3.1 高斯高通滤波器实践代码:clear allIA=imread(girl.bmp); f1,f2=freqspace(size(IA),meshgrid);%D=100/size(IA,1);D=0.3;r=f1.2+f2.2;for i=1:size(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=1-exp(-t); endendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA);title(原始图像);subplot(2,2,2),imshow(uint8(Ia);title(高斯高通滤
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 抹灰、涂料技术交底
- 华北电力大学工程电磁场课件6月10日
- 2024年汽车导航仪项目发展计划
- 氧化还原滴定法
- 天然气灶松手熄火原因的方法
- 2024年生活垃圾焚烧发电项目建议书
- 2024年聚氨酯漆合作协议书
- 2024年网络隔离机(卡)项目合作计划书
- 2024年水稻新品种、新组合合作协议书
- 药品和医疗器械无菌包装全球前13强生产商排名及市场份额
- USB协议中文详解
- 父母资助声明-中英文模板(共2页)
- 新品电池组设计开发流程(共5页)
- 小学生身体红绿灯教学设计(共5页)
- 塑料检查井施工方法(共3页)
- 材料拉伸压缩实验报告
- 基于主题意义探究的初中英语写作教学实践
- 最新冷镦用碳素钢材料线材的化学成分及力学性能
- 《半导体集成电路》考试题目及参考答案
- 四电迁改及过渡方案及工艺
- 超声科主任个人述职报告5篇_超声科医生述职报告
评论
0/150
提交评论