数字图像处理实验报告doc文档格式.docx

数字图像处理实验报告doc文档格式.docx

《数字图像处理实验报告doc文档格式.docx》由会员分享，可在线阅读，更多相关《数字图像处理实验报告doc文档格式.docx（32页珍藏版）》请在冰豆网上搜索。

2.实验内容与要求

平滑空间滤波：

1）

读出一幅图像，给这幅图像分别加入椒盐噪声和高斯噪声后并与前一张图显示在同

一

图像窗口中。

2）对加入噪声图像选用不同的平滑（低通）模板做运算，对比不同模板所形成的效果，要求在同一窗口中显示。

3）使用函数imfilter时，分别采用不同的填充方法（或边界选项，如零填

充、’replicate’、’symmetric’、’circular

’）进行低通滤波，显示处理后

的图

像。

4）

运用for循环，将加有椒盐噪声的图像进行

10

次，20次均值滤波，查看其特点

显

示均值处理后的图像（提示:

利用fspecial

函数的’average’类型生成均值滤

波器）。

5）

对加入椒盐噪声的图像分别采用均值滤波法，

和中值滤波法对有噪声的图像做处理，

要

求在同一窗口中显示结果。

6）自己设计平滑空间滤波器，并将其对噪声图像进行处理，显示处理后的图像。

2.锐化空间滤波

读出一幅图像，采用

3×

3的拉普拉斯算子w=[1,1,1;

1

–81;

1,1,

1]

对其进行滤波。

2）

编写函数w=

genlaplacian（n）

，自动产生任一奇数尺寸

n的拉普拉斯算子，如

5

×

5的拉普拉斯算子

w=[1

1

-24

3）

分别采用

5×

5，9×

9，15×

15

和25×

25大小的拉普拉斯算子对

blurry_moon.tif

进

行锐化滤波，并利用式g（x,y）

2f（x,y）完成图像的锐化增强，观察其有何

f（x,y）

不同，要求在同一窗口中显示。

4）采用不同的梯度算子对该幅图像进行锐化滤波，并比较其效果。

5）自己设计锐化空间滤波器，并将其对噪声图像进行处理，显示处理后的图像；

3.傅立叶变换

1）读出一幅图像，对其进行快速傅立叶变换，分别显示其幅度图像和相位图像。

仅对相位部分进行傅立叶反变换后查看结果图像。

2）仅对幅度部分进行傅立叶反变换后查看结果图像。

3）将图像的傅立叶变换F置为其共轭后进行反变换，比较新生成图像与原始图像的差异。

4.平滑频域滤波

1）设计理想低通滤波器、巴特沃斯低通滤波器和高斯低通滤波器，截至频率自选，分别给出各种滤波器的透视图。

2）读出一幅图像，分别采用理想低通滤波器、巴特沃斯低通滤波器和高斯低通滤波器

对其进行滤波（截至频率自选），再做反变换，观察不同的截止频率下采用不同低通

滤波器得到的图像与原图像的区别，特别注意振铃效应。

（提示:

1）在频率域滤波同

样要注意到填充问题；

2）注意到（-1）;

）

5.锐化频域滤波

1）设计理想高通滤波器、巴特沃斯高通滤波器和高斯高通滤波器，截至频率自选，分别给出各种滤波器的透视图。

2）读出一幅图像，分别采用理想高通滤波器、巴特沃斯高通滤波器和高斯高通滤波器

对其进行滤波（截至频率自选），再做反变换，观察不同的截止频率下采用不同高通滤波器得到的图像与原图像的区别。

3.实验具体实现

1.平滑空间滤波：

（1）.读出一幅图像，给这幅图像分别加入椒盐噪声和高斯噪声后并与前一张图显示在同一图像窗口中。

img=imread（'

lena.png'

figure,subplot（1,3,1）;

imshow（img）;

title（'

原始图像'

）;

img2=imnoise（img,'

salt&

pepper'

0.02）;

subplot（1,3,2）;

imshow（img2）;

title（'

椒盐噪声图像

'

img3=imnoise（img,'

gaussian'

subplot（1,3,3）,imshow（img3）;

高斯噪声图像'

实验结果如下：

（2）.对加入噪声图像选用不同的平滑（低通）模板做运算，对比不同模板所形成的效果，

要求在同一窗口中显示。

平滑滤波是低频增强的空间域滤波技术。

它的目的有两个，一是模

糊，二是消除噪声。

将空间域低通滤波按线性和非线性特点有：

线性、非线性平滑滤波器，

线性平滑滤波器包括均值滤波器，非线性的平滑滤波器有最大值滤波器，中值滤波器，最小值滤波器。

代码如下：

）img=rgb2gray（img）;

imshow（img）;

img2=imnoise（img,'

pepper'

subplot（1,3,2）;

椒盐噪声图像'

高斯噪声图像

%对椒盐噪声图像进行滤波处理

h=fspecial（'

average'

3）;

I1=filter2（h,img2）/255;

I2=medfilt2（img2,[33]）;

figure,subplot（2,2,1）,imshow（img）,title（'

原图像'

subplot（2,2,2）,imshow（img2）,title（'

椒盐噪声图'

subplot（2,2,3）,imshow（I1）,title（'

3*3

均值滤波图'

subplot（2,2,4）,imshow（I2）,title（'

中值滤波图'

%对高斯噪声图像进行滤波处理

G1=filter2（h,img3）/255;

G2=medfilt2（img3,[33]）;

subplot（2,2,2）,imshow（img3）,title（'

高斯噪声图'

subplot（2,2,3）,imshow（G1）,title（'

subplot（2,2,4）,imshow（G2）,title（'

（3）.使用函数

充、’replicate

imfilter’、’

时，分别采用不同

symmetric’、’circular

的填充方法（或边界选项，如零填

’）进行低通滤波，显示处理后的图像。

g=imfilter（f,w,filtering_mode,boundary_options,size_options），其中，f为输入图像，w为滤波掩模，g为滤波后图像。

motion'

50,45）;

%创建一个运动模糊滤波器

filteredimg=imfilter（img,h）;

boundaryReplicate=imfilter（img,h,'

replicate'

boundary0=imfilter（img,h,0）;

boundarysymmetric=imfilter（img,h,'

symmetric'

boundarycircular=imfilter（img,h,'

circular'

figure,subplot（3,2,1）,imshow（img）,title（'

OriginalImage'

subplot（3,2,2）,imshow（filteredimg）,title（'

MotionBlurredImage'

subplot（3,2,3）,imshow（boundaryReplicate）,title（'

Replicate'

subplot（3,2,4）,imshow（boundary0）,title（'

0-Padding'

subplot（3,2,5）,imshow（boundarysymmetric）,title（'

subplot（3,2,6）,imshow（boundarycircular）,title（'

（4）.运用for循环，将加有椒盐噪声的图像进行10次，20次均值滤波，查看其特点

利用fspecial函数的’average’类型生成均值滤波器）。

fori=1:

J1=imfilter（img2,h）;

end

forj=1:

20

J2=imfilter（img2,h）;

figure,subplot（1,3,1）,imshow（img2）,title（'

pepperNoise'

subplot（1,3,2）,imshow（J1）,title（'

10AverageFiltering'

subplot（1,3,3）,imshow（J2）,title（'

20AverageFiltering'

实验结果:

（5）.对加入椒盐噪声的图像分别采用均值滤波法，和中值滤波法对有噪声的图像做处理，要求在同一窗口中显示结果。

h1=fspecial（'

J=imfilter（img2,h1）;

J2=medfilt2（img2）;

subplot（1,3,2）,imshow（J）,title（'

AveragingFiltering'

MedianFiltering'

实验结果为：

（6）.自己设计平滑空间滤波器，并将其对噪声图像进行处理，显示处理后的图像。

[mn]=size（img2）;

figure,subplot（1,2,1）,imshow（img2）;

s=zeros（1,9）;

fori=2:

1:

m-1

forj=2:

n-1

h=1;