数字图像处理实验作业4.docx

1、数字图像处理实验作业4数字图像处理 学院：电信学院 专业：信息工程 班级：2008级3班 姓名：朱志俊 学号：200830101124实验作业4：1、 用理想低通滤波器在频率域实现低通滤波（1） 程序：I=imread(E:大三下数字图像处理MATLAB图像处理水果2.jpg);m,n=size(I);figureimshow(uint8(I);%显示原图像title(原始图像);for i=1:m %对原图像矩阵上的每个元素乘与(-1)(x+y)后再进行傅里叶变换 for j=1:n I(i,j)=I(i,j)*(-1)(i-1+j-1); end end A=fft2(I);A1=log2

2、(1+abs(A);figureimshow(uint8(A1);title(原图像频谱);m,n=size(A);H=zeros(m,n); %构造理想低通滤波器，截止频率为50k1=round(m/2);k2=round(n/2)for i=1:m for j=1:n if sqrt(i-k1)2+(j-k2)2)=50 H(i,j)=1; else H(i,j)=0; end end endF=A.*H; %在频域对原始图像进行滤波处理F1=log2(1+abs(F);figureimshow(uint8(F1);title(低通滤波后图像频谱);f=ifft2(F) %傅里叶反变换并做

3、相应处理（取实部再乘以(-1)(x+y)）f=real(f); for i=1:m for j=1:n f(i,j)=f(i,j)*(-1)(i-1+j-1); end endfigureimshow(uint8(f);title(经理想低通滤波器处理后的图像)（2） 结果：（3） 体会：本实验采用理想低通滤波器对灰度图像进行处理。针对理想低通滤波器的设计是自行设定一个截止频率，采用循环对一个预先设定的零矩阵在满足一定的条件下进行赋值，得到频域理想低通滤波器：m,n=size(A);H=zeros(m,n); %构造理想低通滤波器，截止频率为50k1=round(m/2);k2=round(n

4、/2)for i=1:m for j=1:n if sqrt(i-k1)2+(j-k2)2)=20 H(i,j)=1; else H(i,j)=0; end end endF=A.*(H+ones(m,n)*1.7); %在频域对原始图像进行滤波处理， %并做相应的增强处理figure F1=log2(1+abs(F);imshow(uint8(F1);title(高通滤波后图像频谱); f=ifft2(F) %傅里叶反变换并做相应处理（取实部再乘以(-1)(x+y)）f=real(f); for i=1:m for j=1:n f(i,j)=f(i,j)*(-1)(i-1+j-1); end

5、 endfigureimshow(uint8(f);title(经理想高通滤波器处理后的图像)(2)结果：（3）体会： 将理想低通滤波器的代码稍加改动即可得到理想高通滤波器：理想低通滤波器 m,n=size(A);H=zeros(m,n); %构造理想低通滤波器，截止频率为50k1=round(m/2);k2=round(n/2)for i=1:m for j=1:n if sqrt(i-k1)2+(j-k2)2)=20（修改后） H(i,j)=1; else H(i,j)=0; end end end有了高通滤波器就可以对图像进行高通滤波处理，但为了保证处理后的图像显示效果，有必要做相应的增强处理,在过程中我才用以下代码对对图像做高通滤波增强处理： F=A.*(H+ones(m,n)*1.7);接着就可以顺利得到经高通滤波增强处理的图像。