MatLab理想低通滤波及高通滤波实现高频增强实验报告_第1页
MatLab理想低通滤波及高通滤波实现高频增强实验报告_第2页
MatLab理想低通滤波及高通滤波实现高频增强实验报告_第3页
MatLab理想低通滤波及高通滤波实现高频增强实验报告_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

实验作业5：1、 用理想低通滤波器在频率域实现低通滤波程序代码：I=imread(d:/3.jpg);I=rgb2gray(I);figure(1),imshow(I);title(原图像);s=fftshift(fft2(I);figure(2);imshow(abs(s),);title(图像傅里叶变换所得频谱);a,b=size(s);a0=round(a/2);b0=round(b/2);d=10;for i=1:a for j=1:b distance=sqrt(i-a0)2+(j-b0)2); if distance=d h=1; else h=0; end; s(i,j)=h*s(i,j); end;end;s=uint8(real(ifft2(ifftshift(s);figure(3);imshow(s);title(低通滤波所得图像); 滤波d=102、用理想高通滤波器在频率域实现高频增强I=imread(d:/3.jpg);I=rgb2gray(I);figure(1),imshow(I);title(原图像);s=fftshift(fft2(I);figure(2);imshow(abs(s),);title(图像傅里叶变换所得频谱);figure(3);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);figure(4);imshow(s);title(高通滤波所得图像);figure(5);imshow(s+I);title(高通滤波所得高频增强图像); 心得体会：这次按照低通滤波器和高通滤波器的定义，按照低通滤波的过程，一步一步写，先是进行傅里叶变换，再对其频谱进行平移，使其中心位于中心，再对此时的频谱进行圆形滤波，刚开始纠结于公式，要怎么想出一个H（s）的滤波器的表达式，然后再进行相乘，后来没有想出来，就直接在每次循环里面直接进行乘法运算，如：for i=1:a for j=1:b distance=sqrt(i-a0)2+(j-b0)2); if distance=d h=1; else h=0; end; s(i,j)=h*s(i,j); end;end;后来低通滤波这么一写，高通滤波也就挺简单的了，没怎么想，把上面那个h=1和h=0的位置调换一下就变成了高通滤波器，至于高频增强，就改用了s(i,j)=(p+q*h)*s(i,
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论