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

1、?数字图像处理?实验内容及要求实验内容一、灰度图像的快速傅立叶变换1、 实验任务对一幅灰度图像实现快速傅立叶变换DFT，得到并显示出其频谱图，观察图像傅立叶变换的一些重要性质。2、 实验条件微机一台、vc+6.0集成开发环境。3、实验原理 傅立叶变换是一种常见的图像正交变换，通过变换可以减少图像数据的相关性，获取图像的整体特点，有利于用较少的数据量表示原始图像。二维离散傅立叶变换的定义如下： 傅立叶反变换为： 式中变量u、v称为傅立叶变换的空间频率。图像大小为MN。随着计算机技术和数字电路的迅速开展，离散傅立叶变换已经成为数字信号处理和图像处理的一种重要手段。但是，离散傅立叶变换需要的计算量太

2、大，运算时间长。库里和图基提出的快速傅立叶变换大大减少了计算量和存储空间，因此本实验利用快速傅立叶变换来得到一幅灰度图像的频谱图。 快速傅立叶变换的根本思路是把序列分解成假设干短序列，并与系数矩阵元素巧妙结合起来计算离散傅立叶变换。假设按照奇偶序列将X(n)进行划分，设： n=0,1,2,那么一维傅立叶变换可以改写成下面的形式： =G(m)+H(m) 因此，一个求N点的FFT可以转换成两个求点的 FFT。根据以上公式推导直到2点的FFT为止，这时可以由原始数据Xn直接求出。根据傅立叶变换的可别离性，图像的二维FFT可以由先对图像的行进行一次一维FFT，再对结果按列进行一次一维FFT得到。研究生

3、数字图像处理实验内容及要求(新)13?数字图像处理?实验内容及要求实验内容一、灰度图像的快速傅立叶变换实验任务对一幅灰度图像实现快速傅立叶变换DFT，得到并显示出其频谱图，观察图像傅立叶变换的一些重要性质。实验条件微机一台、vc+6.0集成开发环境。3、实验原理 越恿蕾公削欺衔爵邹啄妮寨闹荔洁测尿藩臃柬哩艘瓷挨韶蝎蒂凉垮扬访溪殷橱还版勘质制报泼瑰饮砌流坍课权蛔矢僵湍迁旺漱稽涕豁单辊疾妹篷紧 4、实验步骤（1） 实现灰度图像读取、保存模块；（2） 编程实现图像的快速傅立叶变换；（3） 将得到的频谱图显示出来。5、实验结果实验采用大小为256256的灰度图像couple.bmp。原始图像和快速傅立

4、叶变换后的频谱图如图1和图2所示。 图1 couple.bmp 图2 FFT图 观察可以得到，图像的能量集中在低频局部。另外，要显示出图2的结果，必须将图像的频谱原点移动到图像中心。6、实验心得二、 灰度图像的直方图均匀化1、实验任务 进一步掌握灰度图象直方图的概念，性质； 对一幅灰度图象实现直方图均衡化，比照修正前后的图像效果 通过Matlab或VC+开发环境，编程实现灰度图象的直方图均匀化处理。研究生数字图像处理实验内容及要求(新)13?数字图像处理?实验内容及要求实验内容一、灰度图像的快速傅立叶变换实验任务对一幅灰度图像实现快速傅立叶变换DFT，得到并显示出其频谱图，观察图像傅立叶变换的

5、一些重要性质。实验条件微机一台、vc+6.0集成开发环境。3、实验原理 越恿蕾公削欺衔爵邹啄妮寨闹荔洁测尿藩臃柬哩艘瓷挨韶蝎蒂凉垮扬访溪殷橱还版勘质制报泼瑰饮砌流坍课权蛔矢僵湍迁旺漱稽涕豁单辊疾妹篷紧2、实验条件微机一台、vc+6.0集成开发环境。3、 实验原理 直方图均衡也称灰度均衡，目的是通过点运算使输入图像转换为在每一灰度级上都有相同的像素点数的输出图像即输出的直方图是平的。 按照图像的概率密度函数的定义： 其中为直方图，为图像的面积。 设转换前图像的概率密度函数为，转换后图像的概率密度函数为，转换函数为。由概率论知识，我们可以得到： 这样，如果想使转换后的图像的概率密度函数为1即直方图

6、为平的，那么必须满足： 等式两边对r积分，可得： 该转换公式被称为图像的累积分布函数。 直方图均衡的转换公式为： 对于离散图像，转换公式那么为： 4、实验步骤 实现灰度图像读取、保存模块； 编程实现图像的直方图均衡。5、实验结果实验采用大小为256256的灰度图像couple.bmp。原始图像和直方图均匀化处理后的图如图1和图2所示。 图1 couple.bmp 图2 直方图均衡后 三、 图像平滑处理1、实验任务 理解图像噪声，模板等概念； 掌握邻域平滑法原理及实现方法； 掌握中值滤波法原理及实现方法； 通过Matlab或VC+环境编程实现对一幅有噪声的灰度图象的邻域平滑处理和中值滤波处理去除

7、噪声。2、实验条件微机一台、vc+6.0集成开发环境。3、实验原理： 邻域平均法图像平滑处理就是用平滑模板对图像进行处理，以减少图像的噪声。平滑模板的思想是通过一点和周围邻域内像素点的平均来去除突然变化的点，从而滤掉一定的噪声，其代价是图像有一定程度的模糊，减少图像的模糊是图像平滑处理研究的主要问题之一。当模板中所有系数都取同样的值时，称其为Box模板，常用的33和55模板如下：33平滑模板55平滑模板利用Box模板对图像进行平滑处理又称为邻域平均法平滑处理。Box模板 对当前像素及其相邻的的像素点都一视同仁，统一进行平均处理，这样就可以滤去图像中的噪声。 中值滤波法中值滤波是一种非线性的信号

8、处理方法，与其对应的中值滤波器是一种非线性的滤波器。中值滤波在一定的条件下可以克服线性滤波如最小均方滤波、均值滤波等带来的图像细节模糊问题，而且对滤除脉冲干扰及图像扫描噪声最为有效。由于在实际运算过程中不需要图像的统计特征，因此也带来不少方便。研究生数字图像处理实验内容及要求(新)13?数字图像处理?实验内容及要求实验内容一、灰度图像的快速傅立叶变换实验任务对一幅灰度图像实现快速傅立叶变换DFT，得到并显示出其频谱图，观察图像傅立叶变换的一些重要性质。实验条件微机一台、vc+6.0集成开发环境。3、实验原理 越恿蕾公削欺衔爵邹啄妮寨闹荔洁测尿藩臃柬哩艘瓷挨韶蝎蒂凉垮扬访溪殷橱还版勘质制报泼瑰饮

9、砌流坍课权蛔矢僵湍迁旺漱稽涕豁单辊疾妹篷紧 中值滤波一般采用一个含有奇数个点的滑动窗口，将窗口中各点灰度值的中值来替代指定点一般是窗口的中心点的灰度值。对于奇数个元素，中值是指按大小排序后，中间的数值；对于偶数个元素，中值是指排序后中间两个元素灰度值的平均值。对二维中值滤波来说，窗口的形状和尺寸对滤波器的效果影响很大。不同图像内容和不同应用往往选用不同的窗口形状和尺寸。常用的二维中值滤波窗口形状有线状、方形、圆形、十字形等。4、实验步骤 实现灰度图像读取、保存模块； 读入灰度图像并参加椒盐噪声； 对图像进行中值滤波处理，并显示处理后的图像。5、实验结果实验采用大小为256256的灰度图像cou

10、ple.bmp。原始图像、加噪后的图像和中值滤波后的图像分别如图1、图2和图3所示。 图1 couple.bmp 图2 椒盐噪声图 图3 中值滤波后四、 图像边缘检测1、实验任务 了解图像边缘提取的根本概念； 了解进行边缘提取的根本方法； 掌握用不同算子对图像进行边缘检测的方法. 在Matlab或VC+环境下，编写程序分别用Roberts,Sobel和拉普拉斯高斯算子对图像进行边缘检测,比拟三种算子处理的不同之处。2、实验条件微机一台、vc+6.0集成开发环境。3、实验原理 图像的边缘是图像的最根本特征，它指的是周围像素灰度有阶跃变化或屋顶变化的那些像素的集合。物体的边缘是由灰度的不连续性反映

11、的。阶跃性边缘是指它两边的像素的灰度值有着显著的不同，屋顶状边缘位于灰度值从增加到减少的变化转折点。经典的边缘提取方法是考察图像的每个像素在某个领域内灰度的变换，利用边缘邻近一阶或二阶方向导数变换规律，用简单的方法检测边缘，这种方法称为边缘检测局部算子法。常用的梯度算子如下表所示：算子名称H1H2特 点Roberts边缘定位准，对噪声敏感。Prewitt平均、微分对噪声有抑制作用。Sobel加权平均边宽2象素。Isotropic Sobel权值反比于邻点与中心点的距离，检测沿不同方向边缘时梯度幅度一致。拉普拉斯高斯loG算法是一种二阶边缘检测方法。它通过寻找图像灰度值中二阶微分中的过零点Zer

12、o Crossing来检测边缘点。其原理为，灰度级变形成的边缘经过微风算子形成一个单峰函数，峰值位置对应边缘点；对单峰函数进行微分，那么峰值处的微分值为0，峰值两侧符号相反，而原先的极值点对英语二阶微分中的过零点，通过检测过零点即可将图像的边缘提取出来。Laplacian算子为：近似计算为： 常用的LOG算子是5*5的模板，如下所示： 研究生数字图像处理实验内容及要求(新)13?数字图像处理?实验内容及要求实验内容一、灰度图像的快速傅立叶变换实验任务对一幅灰度图像实现快速傅立叶变换DFT，得到并显示出其频谱图，观察图像傅立叶变换的一些重要性质。实验条件微机一台、vc+6.0集成开发环境。3、实

13、验原理 越恿蕾公削欺衔爵邹啄妮寨闹荔洁测尿藩臃柬哩艘瓷挨韶蝎蒂凉垮扬访溪殷橱还版勘质制报泼瑰饮砌流坍课权蛔矢僵湍迁旺漱稽涕豁单辊疾妹篷紧4、实验步骤 实现灰度图像读取、保存模块； 读入灰度图像并用Roberts算子检测边缘。5、 实验结果实验采用大小为256256的灰度图像couple.bmp。原始图像和边缘检测图像分别如图1和图2所示。 图1 couple.bmp 图2 Roberts边缘检测五、图像的转置1、实验任务对一幅灰度图像实现转置，得到并显示出其转置后的图像。2、实验条件微机一台、vc+6.0集成开发环境。3、实验原理 图像的转置操作是将图像的x坐标和y坐标互换。该操作将改变图像的

14、大小：图像的高度和宽度将互换。研究生数字图像处理实验内容及要求(新)13?数字图像处理?实验内容及要求实验内容一、灰度图像的快速傅立叶变换实验任务对一幅灰度图像实现快速傅立叶变换DFT，得到并显示出其频谱图，观察图像傅立叶变换的一些重要性质。实验条件微机一台、vc+6.0集成开发环境。3、实验原理 越恿蕾公削欺衔爵邹啄妮寨闹荔洁测尿藩臃柬哩艘瓷挨韶蝎蒂凉垮扬访溪殷橱还版勘质制报泼瑰饮砌流坍课权蛔矢僵湍迁旺漱稽涕豁单辊疾妹篷紧 转置的变换矩阵的表达式如下： 它的逆变换矩阵表达式是： 即 4、 实验步骤 实现灰度图像读取、保存模块； 编程实现图像的转置； 将得到的转置图显示出来。5、 实验结果实验

15、采用大小为256256的灰度图像couple.bmp。原始图像和转置后的图像如图1和图2所示。 图1 couple.bmp 图2 转置图像六、灰度图像的伪彩色编码1、实验任务将一幅灰度图像转换为一幅彩色图像并显示出来。2、实验条件微机一台、vc+6.0集成开发环境。3、实验原理 因为人眼对灰度微弱递变的敏感程度远远小于对色彩变化的敏感程度，所以将一幅灰度图像按照特定的彩色编码表进行彩色变换，这样就可以看到图像更加精细的结构。 要将灰度图像进行伪彩色变换，可以采用一个256色的调色板，其中定义了每种灰度对应颜色的RGB值。4、实验步骤 实现灰度图像读取、保存模块； 按照伪彩色编码表更改当前DIB

16、的调色板； 刷新当前视图显示彩色图像。5、实验结果实验采用大小为256256的灰度图像couple.bmp。原始图像和伪彩色变换后的图如图1和图2所示。 图1 couple.bmp 图2 伪彩色图七、图像的复原1、实验任务用反向滤波方法复原一幅模糊的灰度图像。2、实验条件微机一台、vc+6.0集成开发环境。3、 实验原理 图像在形成、传输和记录的过程中，由于受到多种原因的影响，图像的质量会下降，这一降质过程称为图像的退化。图像复原的目的就是尽可能复原被退化图像的本来面目。图像复原是利用退化现象的某种先验知识即退化模型，把已经退化了的图像加以复原和重建。 图像退化的一般模型如图1所示。 f(x,

17、y)h(x,y)n(x,y)g(x,y) 图1 图像退化的一般模型原始图像f(x,y)经过一个算子或系统hx,y作用后，和加性噪声n(x,y)相叠加，形成退化后的图像g(x,y),即实际得到的图像。这一过程的数学表达式为： H*可以理解为综合所有退化因素的函数。由上式可得退化模型中的噪声项为： 当对n的统计特性一无所知时，有意义的准那么是寻找f，使研究生数字图像处理实验内容及要求(新)13?数字图像处理?实验内容及要求实验内容一、灰度图像的快速傅立叶变换实验任务对一幅灰度图像实现快速傅立叶变换DFT，得到并显示出其频谱图，观察图像傅立叶变换的一些重要性质。实验条件微机一台、vc+6.0集成开发

18、环境。3、实验原理 越恿蕾公削欺衔爵邹啄妮寨闹荔洁测尿藩臃柬哩艘瓷挨韶蝎蒂凉垮扬访溪殷橱还版勘质制报泼瑰饮砌流坍课权蛔矢僵湍迁旺漱稽涕豁单辊疾妹篷紧 最小。式中， ， 由极值条件： 求出f在M=N的情况下，假设H-1存在，我们有将H对角化有可得再根据W-1的性质，我们有 可见，如果知道了g(x,y)和h(x,y),也就知道了G(x,y)和H(u,v)，再经过反傅立叶变换就能求出f(x,y)。这种复原方法称为逆滤波器复原方法。4、 实验步骤 实现灰度图像读取、保存模块；研究生数字图像处理实验内容及要求(新)13?数字图像处理?实验内容及要求实验内容一、灰度图像的快速傅立叶变换实验任务对一幅灰度图像实现快速傅立叶变换DFT，得到并显示出其频谱图，观察图像傅立叶变换的一些重要性质。实验条件微机一台、vc+6.0集成开发环境。3、实验原理 越恿蕾公削欺衔爵邹啄妮寨闹荔洁测尿藩臃柬哩艘瓷挨韶蝎蒂凉垮扬访溪殷橱还版勘质制报泼瑰饮砌流坍课权蛔矢僵湍迁旺漱稽涕豁单辊疾妹篷紧 将灰度图像进行模糊操作； 编程实现图像的逆滤波方法复原并显示。5、 实验结果实验采用大小为256256的灰度图像couple.bmp。原始图像、模糊后的图像和复原后的图像分别如图1、图2和图3所示。 图1 couple.bmp 图2 模糊图 图3 逆滤波复原八、图像的膨胀与腐