实验三图像空域平滑和频域平滑

1、重庆大学光电工程学院数字图像处理实验指导书实验三第1页共8页重庆大学光电工程学院数字图像处理实验指导书实验三实验三图像空域平滑和频域平滑一、实验目的1）掌握图像典型噪声的基本特点；2）掌握图像空域高斯平滑的基本方法；3）掌握图像空域中值滤波的基本方法；4）掌握图像频域高斯平滑的基本方法；5）掌握根据图像特点进行平滑滤波的基本原理和方法。6）通过编程，上机调试程序，进一步提高编程能力及使用计算机解决问题 的能力。二、实验原理1. 图像的典型噪声噪声可理解为影响传感器对所接收图像源信息进行理解或分析的各种因素。 噪声一般是不可预测的随机信号，只能用概率统计的方法去认识。噪声对图像的 输入、采集和处

2、理的各种环节，以及输出结果的全过程均有影响。因此，去噪己 经成为图像处理中极其重要的手段，也是图像处理领域研究的一个重点。本节对 典型噪声的来源和性质进行简要介绍。对图像信号而言，灰度图像可视为二维亮度分布. ：，噪声可看作是对亮 度的干扰，用.:表示。噪声具有随机性，因而需要用随机过程来描述，即要 求知道其分布函数和密度函数。在许多情况下，这些函数很难测定和描述，甚至 无法得到，所以常用均值、方差、相关函数等统计特征来描述噪声，如噪声的总 功率描述为一；噪声的交流功率可由方差 心心上述-逍加3，疔描述；噪声的直流功率可用均值的平方； .表示。图像噪声的描述与建模方式主要有以下几种：（1）白噪

3、声（White Noise）：它具有常量的功率谱。白噪声的一个特例是高 斯噪声（Gaussian Noise，其直方图曲线服从一维高斯型分布：(1)其中，为均值与标准差。MATLAB产生函数为J = IMNOISE（l,gaussian,M,V），M，V为均值与方差，高斯噪声示例如图1所示。输入图像高斯噪声,均值Q.1方差00图1高斯噪声示例(2)椒盐噪声(Pepper & Salt Noise): 种在图像中产生黑色、白色点的 脉冲噪声。该噪声在图像中显现较为明显，对图像分害IJ、边缘检测、特征提取等 后续处理具有严重破坏作用。MATLAB 产生函数为 J = IMNOISE(I, alt

4、&pepperD), D为噪声密度。椒盐噪声示例如图 2所示。输入图像檄盐噪声，噪声密度01图2椒盐噪声示例(3)乘性噪声(Multiplicative Noise):乘性噪声与加性噪声对应，加性噪 声一般指热噪声、散弹噪声等，它们与信号的关系是相加，不管有没有信号，噪 声都存在。而乘性噪声一般由信道不理想引起，它们与信号的关系是相乘，依附于信号而存在。一般图像处理中把加性随机性噪声看成是系统的背景噪声；而乘性随机性噪声看成系统的时变性(如衰落或多普勒)或非线性所造成的。具体而 言乘性噪声是指一幅图像被斑点噪声像素破坏，而且这些像素的亮度与其邻域的 亮度显著不同，公式可表达为- -o MATL

5、AB产生函数为J=IMNOISE(l,speckle,V)，其中V为噪声方差。乘性噪声示例如图 3所示。输入图像乘性噪声，噪声方差0茨图3乘性噪声示例（4）量化噪声（Quatization Noise）:在量化级别不同时出现的噪声。例如, 将图像的亮度级别减少到1/16的时候会出现伪轮廓。量化噪声可从直方图分布 予以分析。量化噪声示例如图4所示。输入图像灰度级别2站量化噪声，灰度级别怡图4量化噪声示例2. 图像空域高斯滤波原理图像平滑是为了消除或减弱噪声对图像的干扰，有时也称图像滤波。图像平滑总体上讲可分为空域平滑（滤波）和频域平滑（滤波）两大类，根据滤波器特 点也可分为线性滤波和非线性滤波。

6、空域滤波是直接对图像的数据做空间变换达 到滤波的目的。它是一种邻域运算，即输出图像中任何像素的值都是通过采用一 定的算法，根据输入图像中对用像素周围一定邻域内像素的值得来的。如果输出像素是输入像素邻域像素的线性组合则称为线性滤波（例如最常见的均值滤波和 高斯滤波），否则为非线性滤波（中值滤波等）。高斯滤波是利用高斯滤波器对图像进行滤波的过程， 这种线性平滑滤波器去除高斯噪声的效果很好，且在大多数情况下，对其它类型的噪声也有较好的效果。 高斯滤波器使用连续窗函数内像素加权和来实现滤波，它使用的高斯函数如图5所示。图5二维高斯函数实际使用高斯滤波器时要利用量化值，例如尺寸为3X 3的高斯滤波器(卷

7、第5页共8页重庆大学光电工程学院数字图像处理实验指导书实验三第#页共8页重庆大学光电工程学院数字图像处理实验指导书实验三积核)可简化为环=丄9 1G尺寸为5X 5的高斯滤波器可简化为2421484281GS448422 42 1滤波器元素取一.I.的原因是便于硬件实现。图像的空滤滤波公式可 表达为利用卷积公式有:! j :/；- - .(3)J1对于3X3的高斯滤波，核心代码可表示为O(x,y) = 0;for s = 1:3for t. = 1:3O(x.y) = O(x, y) + 爲貞艮2ty + f-2):endend由于滤波器模板尺寸因素，空域滤波需要考虑图像边界状况3. 图像空域中

8、值滤波原理中值滤波器是统计滤波器的一种，属于非线性的的空间滤波器。它将像素（中 值计算中包括原像素值）邻域内灰度的中值代替该像素的值。中值滤波器对处理 脉冲噪声非常有效，对于一定类型的随机噪声，中值滤波器也提供了优秀的去噪 能力，比小尺寸的线性平滑滤波器的模糊程度明显要低。中值滤波器的计算步骤：对于图像中每一像素首先设定其邻域像素 值集合（尺寸一般设为 1X 3，1X 5, 3X1，5X1, 3X 3, 5X 5，取奇数 是为了便于得到中值。），然后对其邻域集合元素进行排序，接着选取其中值像素 值，最后将该中值像素值赋给原图像素,.：0对于椒盐噪声，空域高斯滤波和 空域中值滤波的结果如图6所示

9、。输入图像橄盐噪声高斯滤波结果中值滤浪结果图6椒盐噪声的空域滤波结果示意4. 图像频域高斯滤波原理频率滤波需要先进行傅立叶变换至频域处理， 然后再反变换回空间域还原图 像。图像的二维傅立叶变换及逆变换公式为1A/-1 N-1FE)= MN工工心必如晋*即JI N-1=刀为F(u.硏严箭+卸u=0 t?=0这里为和数字信号处理统一，图像元素写作高斯滤波器的空域和频 域形式相同，从频域来看是一个低通滤波器，没有振铃现象。高斯传递函数在频 域可表示为其中 :为点I到中心距离，L为预设参数，即截止频率。高斯传递 函数及径向剖面曲线如图7所示。图7高斯传递函数径向剖面曲线空域滤波及对应的频域滤波可依次由

10、以下三个公式简单表示(用)=hx.G(u,v) = H(u,v)X(u,v) 曲皿=L (G(u.v)这里 ，分别代表待滤波图像在空域和频域的表达形式，v)分别代表已滤波图像在空域和频域的表达形式，/心期.仏础分 别代表滤波器冲击响应函数及传递函数。显然，空域卷积等同于频域乘积。使用 的MATLAB函数包括fft2, ifft2, fftshift,ifftshift 等。根据FFT运算特点，频域滤 波时对图像与滤波器尺寸的相容性必须充分考虑。图8示意了频域高斯滤波结果，注意为显示频域细节，对频幅图进行了一变换。输入图像高斯噪芦爵偉频幅图高斯传递函数滤液频帽图频域高斯滤液结果图8频域高斯滤波结

11、果三、主要实验仪器及材料微型计算机、Matlab6.5教学版四、实验内容1 知识准备认真复习以上基础理论，理解本实验所用到的实验原理。2 图像典型噪声分析1）利用imnoise函数为图像添加高斯噪声，改变参数，分析噪声图像特点。2）利用imnoise函数为图像添加椒盐噪声，改变参数，分析噪声图像特点。3）利用imnoise函数为图像添加乘性噪声，改变参数，分析噪声图像特点。3 图像空域高斯滤波对于图像中三类典型噪声（高斯噪声，椒盐噪声，乘性噪声，下同），禾I用 高斯滤波器进行滤波处理，改变卷积核尺寸，观察分析去噪实验结果。4 图像空域中值滤波对于图像中三类典型噪声，利用中值滤波器进行滤波处理，改变滤波器尺寸， 观察分析去噪实验结果。5 图像频域域高斯滤波将叠加了三类典型噪声之一的图像利用 FFT2方法变换到频域，利用频域高 斯滤波器进行滤波处理，然后变换回空域，改变滤波器参数，观察分析去噪实验 结果。6上机调试程序并打印或记录实验结果。7 完成实验报告。五、实验报告要求1）简述实验原理及目的。2）给出上述实验内容的实验