研究生数字图像处理实验内容及要求(新)

《数字图像处理》实验内容及要求

实验内容

一、灰度图像的快速傅立叶变换

1、 实验任务

对一幅灰度图像实现快速傅立叶变换（DFT ），得到并显示出其频谱图，观察图像傅立叶变换的一些重要性质。 2、 实验条件

微机一台、vc++6.0集成开发环境。 3、实验原理

傅立叶变换是一种常见的图像正交变换，通过变换可以减少图像数据的相关性，获取图像的整体特点，有利于用较少的数据量表示原始图像。 二维离散傅立叶变换的定义如下：

11

2(

)00

(,)(,)ux vy

M N j M N

x y F u v f x y e

π---+===

∑∑

傅立叶反变换为：

112(

)00

1

(,)(,)ux vy M N j M N

u v f x y F u v e

MN

π--+===

∑∑

式中变量u 、v 称为傅立叶变换的空间频率。图像大小为M ×N 。随着计算机技术和数字电路的迅速发展，离散傅立叶变换已经成为数字信号处理和图像处理的一种重要手段。但是，离散傅立叶变换需要的计算量太大，运算

时间长。库里和图基提出的快速傅立叶变换大大减少了计算量和存储空间，因此本实验利用快速傅立叶变换来得到一幅灰度图像的频谱图。

快速傅立叶变换的基本思路是把序列分解成若干短序列，并与系数矩阵元素巧妙结合起来计算离散傅立叶变换。若按照奇偶序列将X(n)进行划分，设：

()(2)

()(21)

g n x n h n x n =??=+? （n=0,1,2,…,12N -）

则一维傅立叶变换可以改写成下面的形式：

1

0()()N mn

N

n X m x n W -==∑

11220

()()N N mn mn

N N n n g n W h n W --===+∑∑ 1122(2)(21)

(2)(21)N

N m n m n N N n n x n W x n W --+===++∑∑

11220

2

2

(2)(21)N

N mn mn m

N N N

n n x n W x n W W --===++∑∑

W H(m)

=G(m)+m N

N点的FFT。根据以上公因此，一个求N点的FFT可以转换成两个求

2

式推导直到2点的FFT为止，这时可以由原始数据X（n）直接求出。根据傅立叶变换的可分离性，图像的二维FFT可以由先对图像的行进行一次一维FFT，再对结果按列进行一次一维FFT得到。

4、实验步骤

（1）实现灰度图像读取、保存模块；

（2）编程实现图像的快速傅立叶变换；

（3）将得到的频谱图显示出来。

5、实验结果

实验采用大小为256×256的灰度图像couple.bmp。原始图像和快速傅立叶变换后的频谱图如图1和图2所示。

图1 couple.bmp 图2 FFT图

观察可以得到，图像的能量集中在低频部分。另外，要显示出图2的结果，必须将图像的频谱原点移动到图像中心。

6、实验心得

二、 灰度图像的直方图均匀化 1、实验任务

⑴ 进一步掌握灰度图象直方图的概念，性质；

⑵ 对一幅灰度图象实现直方图均衡化，对比修正前后的图像效果 ⑶ 通过Matlab 或VC++开发环境，编程实现灰度图象的直方图均匀化

处理。

2、实验条件

微机一台、vc++6.0集成开发环境。

3、 实验原理

直方图均衡也称灰度均衡，目的是通过点运算使输入图像转换为在每一灰度级上都有相同的像素点数的输出图像（即输出的直方图是平的）。 按照图像的概率密度函数的定义：

()A

D MAX

B A i i D

D f D H A ===∑0

1

()()p x H x A =

其中()H x 为直方图，

0A 为图像的面积。

设转换前图像的概率密度函数为

()r p r ，转换后图像的概率密度函数

为

()s p s ，转换函数为()s f r =。由概率论知识，我们可以得到：

()()

s r dr p s p r ds =

这样，如果想使转换后的图像的概率密度函数为1（即直方图为平的），

则必须满足：

() r ds

p r

dr

=

等式两边对r积分，可得：

00

1

()()()

r r

r

s f r P u du H u du

A

===

??

该转换公式被称为图像的累积分布函数。

直方图均衡的转换公式为：

()()

A

D

MAX

B A

D

D f D H u du

A

=

=?

对于离散图像，转换公式则为：

()A

D

MAX

B A i

i

D

D f D H

A

=

=

=∑

4、实验步骤

⑴实现灰度图像读取、保存模块；

⑵编程实现图像的直方图均衡。

5、实验结果

实验采用大小为256×256的灰度图像couple.bmp。原始图像和直方图均匀化处理后的图如图1和图2所示。

图1 couple.bmp 图2 直方图均衡后

三、图像平滑处理

1、实验任务

⑴理解图像噪声，模板等概念；

⑵掌握邻域平滑法原理及实现方法；

⑶掌握中值滤波法原理及实现方法；

⑷通过Matlab或VC++环境编程实现对一幅有噪声的灰度图象的邻域平

滑处理和中值滤波处理去除噪声。

2、实验条件

微机一台、vc++6.0集成开发环境。

3、实验原理：

⑴邻域平均法

图像平滑处理就是用平滑模板对图像进行处理，以减少图像的噪声。平滑模板的思想是通过一点和周围邻域内像素点的平均来去除突然变化的点，从而滤掉一定的噪声，其代价是图像有一定程度的模糊，减少图像的模糊是图像平滑处理研究的主要问题之一。当模板中所有系数都取同样的值时，称其为Box模板，常用的3×3和5×5模板如下：

3×3平滑模板5×5平滑模板

利用Box模板对图像进行平滑处理又称为邻域平均法平滑处理。Box模板对当前像素及其相邻的的像素点都一视同仁，统一进行平均处理，这样就可以滤去图像中的噪声。

⑵中值滤波法

中值滤波是一种非线性的信号处理方法，与其对应的中值滤波器是一种非线性的滤波器。中值滤波在一定的条件下可以克服线性滤波如最小均方滤波、均值滤波等带来的图像细节模糊问题，而且对滤除脉冲干扰及图像扫描噪声最为有效。由于在实际运算过程中不需要图像的统计特征，因此也带来不少方便。

中值滤波一般采用一个含有奇数个点的滑动窗口，将窗口中各点灰度值的中值来替代指定点（一般是窗口的中心点）的灰度值。对于奇数个元素，中值是指按大小排序后，中间的数值；对于偶数个元素，中值是指排序后中间两个元素灰度值的平均值。对二维中值滤波来说，窗口的形状和尺寸对滤波器的效果影响很大。不同图像内容和不同应用往往选用不同的窗口形状和尺寸。常用的二维中值滤波窗口形状有线状、方形、圆形、十字形等。

4、实验步骤

⑴实现灰度图像读取、保存模块；

⑵读入灰度图像并加入椒盐噪声；

⑶对图像进行中值滤波处理，并显示处理后的图像。

5、实验结果

实验采用大小为256×256的灰度图像couple.bmp。原始图像、加噪后的图像和中值滤波后的图像分别如图1、图2和图3所示。

图1 couple.bmp 图2 椒盐噪声图图3 中值滤波后

四、图像边缘检测

1、实验任务

⑴了解图像边缘提取的基本概念；

⑵了解进行边缘提取的基本方法；

⑶掌握用不同算子对图像进行边缘检测的方法.

⑷在Matlab或VC++环境下，编写程序分别用Roberts,Sobel和拉普

拉斯高斯算子对图像进行边缘检测,比较三种算子处理的不同之处。

2、实验条件

微机一台、vc++6.0集成开发环境。

3、实验原理

图像的边缘是图像的最基本特征，它指的是周围像素灰度有阶跃变化或屋顶变化的那些像素的集合。物体的边缘是由灰度的不连续性反映的。

阶跃性边缘是指它两边的像素的灰度值有着显著的不同，屋顶状边缘位于灰度值从增加到减少的变化转折点。经典的边缘提取方法是考察图像的每个像

素在某个领域内灰度的变换，利用边缘邻近一阶或二阶方向导数变换规律，用简单的方法检测边缘，这种方法称为边缘检测局部算子法。

常用的梯度算子如下表所示：

算子名称H1 H2 特点

Roberts 边缘定位准，对噪声敏感。

Prewitt 平均、微分对噪声有抑制作用。

Sobel 加权平均边宽≥2象素。

Isotropic Sobel 权值反比于邻点与中心点的距离，检测沿不同方向边缘时梯度幅度一致。

拉普拉斯高斯（loG）算法是一种二阶边缘检测方法。它通过寻找图像灰度值中二阶微分中的过零点（Zero Crossing）来检测边缘点。其原理为，灰度级变形成的边缘经过微风算子形成一个单峰函数，峰值位置对应边缘点；对单峰

函数进行微分，则峰值处的微分值为0，峰值两侧符号相反，而原先的极值点对英语二阶微分中的过零点，通过检测过零点即可将图像的边缘提取出来。Laplacian算子为：

近似计算为：

常用的LOG算子是5*5的模板，如下所示：

4、实验步骤

⑴实现灰度图像读取、保存模块；

⑵读入灰度图像并用Roberts算子检测边缘。

5、实验结果

实验采用大小为256×256的灰度图像couple.bmp。原始图像和边缘检测图像分别如图1和图2所示。

图1 couple.bmp 图2 Roberts 边缘检测

五、图像的转置

1、实验任务

对一幅灰度图像实现转置，得到并显示出其转置后的图像。 2、实验条件

微机一台、vc++6.0集成开发环境。 3、实验原理

图像的转置操作是将图像的x 坐标和y 坐标互换。该操作将改变图像的大小：图像的高度和宽度将互换。 转置的变换矩阵的表达式如下：

1010011

000100

11x x y y ??????

?

? ?= ? ? ? ? ? ????

???

它的逆变换矩阵表达式是：

001010100110011x x y y ??????

? ???=

? ??? ? ?????????

即 01

01x y y x =??=?

4、实验步骤

⑴实现灰度图像读取、保存模块；

⑵编程实现图像的转置；

⑶将得到的转置图显示出来。

5、实验结果

实验采用大小为256×256的灰度图像couple.bmp。原始图像和转置后的图像如图1和图2所示。

图1 couple.bmp 图2 转置图像

六、灰度图像的伪彩色编码

1、实验任务

将一幅灰度图像转换为一幅彩色图像并显示出来。

2、实验条件

微机一台、vc++6.0集成开发环境。

3、实验原理

因为人眼对灰度微弱递变的敏感程度远远小于对色彩变化的敏感程度，所以将一幅灰度图像按照特定的彩色编码表进行彩色变换，这样就可以