数字图像处理实验报告.docVIP

实验报告书

系部名称

：

学生姓名

：

专业名称

：

班级

：

实验时间

：

实验一直方图均衡

一、实验目的

在学习图像直方图的概念、计算方法、性质和相关应用基础上，生成、绘制图像的直方图，并应用MATLAB编程实现图像直方图均衡化程序。

二、实验内容

（1）计算并绘制图像直方图；

（2）编程实现图像的直方图均衡化处理，显示均衡前后的直方图和图像；

三、实验运行结果

实验二频域图像增强

一、实验目的

1、频域图像增强

2、掌握基于频域的图像增强方法。

二、实验内容

（1）编程实现图像的理想低通和高通滤波；

（2）编程实现图像的巴特沃斯低通和高通滤波。

三、实验运行结果

四、实验中遇到的问题及解决方法

显示图像无法打开，最终查出来时图像格式弄错了。

五、思考题

分析为什么图像通过低通滤波器后变得模糊？为什么通过高通滤波器后得到锐化结果？

答：图像的精细结构及突变部分主要由高频成分起作用，故经低通滤波后图像的精细结构消失，变得模糊；经高通滤波后图像得到锐化。

六、实验心得体会

本实验中遇到很多问题及错误，例如图像打不开、处理后图像模糊等，都是经常容易发生的错误，最后实验几次，就能够逐一自己解决了。使自己对数字图像处理课程中的许多问题有了更实际和确切的深入了解。

七、程序清单

clc;

clear;

data4=imread(lena.gif);

subplot(3,2,1);

imshow(data4);

title(原图);

i=fft2(data4);

subplot(3,2,2);

i=fftshift(i);

z=log(abs(i));

x=0:1:255;

y=0:1:255;

[x,y]=meshgrid(x,y);

mesh(z);%以三维坐标显示该图像频谱图

title(原图频谱);

[n,m]=size(i);%对该图进行低通滤波

fork=1:1:n

forl=1:1:m

if(k^2+l^2)=190^2%选取D=190

result(k,l)=0;

elseresult(k,l)=i(k,l);

end

end

end

subplot(3,2,4);

z=log(abs(result));%三维方式显示低通滤波后的频谱图

x=0:1:255;

y=0:1:255;

[x,y]=meshgrid(x,y);

mesh(z);

title(理想低通滤波后的频谱);

subplot(3,2,3);%新建图像显示窗口

result=fftshift(result);%滤波后的数据去中心化

b=ifft2(result);%逆傅里叶变换

imshow(uint8(abs(b)));

title(理想低通滤波后的图像);

subplot(3,2,6);%新建图像显示窗口

%[n,m]=size(c);%对原图进行高通滤波

fork=1:1:n

forl=1:1:m

if(k^2+l^2)=190^2%选取D=190

result(k,l)=0;

elseresult(k,l)=i(k,l);

end

end

end

z=log(abs(result));

x=0:1:255;%三维方式显示高通滤波前的频谱图

y=0:1:255;

[x,y]=meshgrid(x,y);

mesh(z);

title(理想高通滤波后的频谱);

subplot(3,2,5);

result=fftshift(result);%滤波后的数据去中心化

d=ifft2(result);%逆傅里叶变换

imshow(uint8(abs(d)));

title(理想高通滤波后的图像);

%频域增强（巴特沃斯原型）

%二阶巴特沃斯（Butterworth）低通滤波器

%clc;

%clear;

Figure;

J1=imread(