基于MATLAB的彩色图像平滑处理1.doc

1、1目录第一章、概述 21.1 图像平滑概述 21.2 图像平滑应用 .2第二章 、图像平滑方法 .52.1 空域低通滤波 52.1.1 均值滤波器 .52.1.2 中值滤波器 .62.2 频域低通滤波 7第三章、图像平滑处理与调试 .103.1 模拟噪声图像 .103.2 均值滤波法 123.3 中值滤波法 .153.4 频域低通滤波法 .17第四章、总结与体会 .20参考文献 .21第一章、概述21.1 图 像 平 滑 概 述图 像 平 滑 （ Smoothing） 的 主 要 目 的 是 减 少 图 像 噪 声 。 图 像噪 声 来 自 于 多 方 面 ， 有 来 自 于 系 统 外 部

2、的 干 扰 （ 如 电 磁 波 或 经 电 源窜 进 系 统 内 部 的 外 部 噪 声 ） ， 也 有 来 自 于 系 统 内 部 的 干 扰 （ 如 摄像 机 的 热 噪 声 ， 电 器 机 械 运 动 而 产 生 的 抖 动 噪 声 内 部 噪 声 ） 。实 际 获 得 的 图 像 都 因 受 到 干 扰 而 有 噪 声 ， 噪 声 产 生 的 原 因 决 定 了 噪声 分 布 的 特 性 及 与 图 像 信 号 的 关 系 。 减 少 噪 声 的 方 法 可 以 在 空 间 域或 在 频 率 域 处 理 。 在 空 间 域 中 进 行 时 ， 基 本 方 法 就 是 求 像 素 的 平

3、 均值 或 中 值 ； 在 频 域 中 则 运 用 低 通 滤 波 技 术 。图 像 中 的 噪 声 往 往 是 和 信 号 交 织 在 一 起 的 ， 尤 其 是 乘 性 噪 声 ，如 果 平 滑 不 当 ， 就 会 使 图 像 本 身 的 细 节 如 边 缘 轮 廓 ， 线 条 等 模 糊 不清 ， 从 而 使 图 像 降 质 。 图 像 平 滑 总 是 要 以 一 定 的 细 节 模 糊 为 代 价的 ， 因 此 如 何 尽 量 平 滑 掉 图 像 的 噪 声 ， 又 尽 量 保 持 图 像 的 细 节 ， 是图 像 平 滑 研 究 的 主 要 问 题 之 一 。1.2图 像 平 滑 应

4、 用图 像 平 滑 主 要 是 为 了 消 除 被 污 染 图 像 中 的 噪 声 ,这 是 遥感 图 像 处 理 研 究 的 最 基 本 内 容 之 一 ,被 广 泛 应 用 于 图 像 显 示 、传 输 、 分 析 、 动 画 制 作 、 媒 体 合 成 等 多 个 方 面 。 该 技 术 是 出 于 人 类视 觉 系 统 的 生 理 接 受 特 点 而 设 计 的 一 种 改 善 图 像 质 量 的 方 法 。 处 理对 象 是 在 图 像 生 成 、 传 输 、 处 理 、 显 示 等 过 程 中 受 到 多 种 因 素 扰 动形 成 的 加 噪 图 像 。 在 图 像 处 理 体 系

5、 中 ,图 像 平 滑 是 图 像 复 原 技术 针 对 “一 幅 图 像 中 唯 一 存 在 的 退 化 是 噪 声 ”时 的 特 例 。1.3噪声模型1.3.1噪声来源一幅图像可能会受到各种噪声的干扰，而数字图像的实质就是光电信息，因此图像噪声主要可能来源于以下几个方面：光电传感器噪声、大气层电磁暴、闪电等引起的强脉冲干扰、相片颗粒噪声和信道传输误差引起的噪声等。噪声的存在恶化图像质量，使图像模糊，更严重的甚至是图像的特征完全被淹没，以至于给图像识别和分析带来了困难。目前比较经典的去噪声的方法都或多或少给图像带来模糊，因此，探求一种既能去除噪声又不至于使图像模糊的方法，一直是图像增强处理中

6、的难题，3至今尚在不断地探索。2.1.1噪声的分类所谓噪声，就是妨碍人的视觉器官或系统传感器对所接收的图像信息进行理解或分析的各种因素。一般噪声是不可预测的随机信号，它只能用概率统计的方法去认识。依据噪声产生的原因，我们可将经常影响图像质量的噪声源分为三类：2.1.1阻性元器件内部产生的高斯噪声，这类噪声是由于元器件中的电子随机热运动而造成的，很早就被人们成功的建模并研究，一般常用零均值高斯白噪声作为其模型。2.1.2光电转换过程中的泊松噪声,这类噪声是由光的统计本质和图像传感器中光电转换过程引起的。在弱光情况下，影响更为严重，常用具有泊松密度分布的随机变量作为这类噪声的模型。2.1.3感光过

7、程中产生的颗粒噪声，在显微镜下检查可发现，照片上光滑细致的影调在微观上其实呈现一个随机的颗粒性质。此外颗粒本身大小的不同以及每一颗粒暴光所需光子数目的不同，都会引入随机性。这些因素的外观表现称为颗粒性。对于多数应用，颗粒噪声可用高斯过程(白噪声)作为有效模型。根据噪声和信号的关系可以将其分为两种形式：(1)加性噪声：有的噪声与图像信号 g(x,y)无关，在这种情况下，含噪图像f(x,y)可表示为：f(x,y)= g(x,y)+n(x,y)信道噪声及扫描图像时产生的噪声都属于加性噪声。(2)乘性噪声：有的噪声与图像信号有关，这可以分为两种情况：一种是某像素点的噪声只与该像素点的图像信号有关，另一

8、种是某像素点的噪声与该点及其邻域的图像信号有关。如果噪声和信号成正比，则含噪图像 f(x,y)可以表示为：f(x,y)= g(x,y)+ n(x,y)g(x,y)图 加乘性噪声图另外，还可以根据噪声服从的分布对其进行分类，这时可以分为高斯噪声、4泊松噪声和颗粒噪声等。泊松分布噪声一般出现在照度非常小及用高倍电子线路放大的情况下，泊松噪声可以认为是椒盐噪声。其他的情况通常为加性高斯噪声。颗粒噪声可以认为是一种白噪声过程，在密度域中是高斯分布加性噪声，而在强度域中为乘性噪声。1.3.2 高斯噪声数字图像的噪声主要来源于图像的获取和传输过程。按其产生的原因可分为:光电子噪声、热噪声、KTC 噪声、量

9、化噪声和信道传输噪声等。按其是否独立于空间坐标以及和图像是否关联可分为加性噪声和乘性噪声。为了最大限度地减少噪声对图像的影响,人们从改善硬件质量和对受污图像进行处理两个方面做了许多的工作,文中主要考虑对受污图像进行处理的算法研究。为了对受污图像进行处理,人们对噪声进行了研究并建立了相应的数学模型。对噪声表述的数学建模主要考虑噪声的成因和分析受污图像上噪声的统计特性两个因素，这种噪声主要来源于电子电路噪声和低照明度或高温带来的传感器噪声,也称为正态噪声,是在实践中经常用到的噪声模型。高斯随机变量 z 的概率密度函数 ( PDF)由下式给出:其中, z 表示图像像元的灰度值; 表示 z 的期望；

10、表示 z 的标准差。图加高斯噪声图1.3.3椒盐噪声主要来源于成像过程中的短暂停留和数据传输中产生的错误。其 PDF 为如果 b a,灰度值 b 在图像中显示为一亮点,a 值显示为一暗点。如果 Pa 和 Pb 均不为零,在图像上的表现类似于随机分布图像上的胡椒和盐粉微粒,因此称为椒盐噪声。当 Pa 为零时,表现为“盐”噪声;当 Pb 为零时,表现为“胡椒”噪声。MATLAB图像处理工具箱提供的噪声添加函数imnoise,它可以对图像添加一些典型的噪声。其语法：5J=imnoise(I,type)J=imnoise(I,type,parameters)其功能是：返回对原图像I添加典型噪声的图像J

11、,参数 type和parameters用于确定噪声的类型和相应的参数。三种典型的噪声：type=gaussian时,为 高 斯 噪 声；type=salt %读取图像I1=imnoise(I,gaussian); %加高斯噪声I2=imnoise(I,salt %加椒盐噪声I3=imnoise(I,speckle); %加乘性噪声subplot(221),imshow(I); %显示图像Isubplot(222),imshow(I1);subplot(223),imshow(I2);subplot(224),imshow(I3);运行结果如下：6(a)乘 乘 乘 (b)乘 乘 乘 乘 乘 乘

12、乘 乘 乘 乘(c)乘 乘 乘 乘 乘 乘 乘 乘 乘 乘 (d)乘 乘 乘 乘 乘 乘 乘 乘 乘 乘图 1 噪声污染的图像第 2 章、图像平滑方法2.1 空 域 低 通 滤 波将 空 间 域 模 板 用 于 图 像 处 理 ， 通 常 称 为 空 间 滤 波 ， 而 空 间 域 模板 称 为 空 间 滤 波 器 。 空 间 域 滤 波 按 线 性 和 非 线 性 特 点 有 ： 线 性 、 非线 性 平 滑 波 器 。线 性 平 滑 滤 波 器 包 括 领 域 平 均 法 （ 均 值 滤 波 器 ） ， 非 线 性平 滑 滤 波 器 有 中 值 滤 波 器 。72.1.1 均 值 滤 波

13、器对 一 些 图 像 进 行 线 性 滤 波 可 以 去 除 图 像 中 某 些 类 型 的 噪 声 ，如 采 用 邻 域 平 均 法 的 均 值 滤 波 器 就 非 常 适 用 于 去 除 通 过 扫 描 得 到 的图 像 中 的 颗 粒 噪 声 。 邻 域 平 均 法 是 空 间 域 平 滑 技 术 。 这 种 方法 的 基 本 思 想 是 ， 在 图 像 空 间 ， 假 定 有 一 副 NN 个 像 素 的 原始 图 像 f(x,y),用 领 域 内 几 个 像 素 的 平 均 值 去 代 替 图 像 中 的 每 一个 像 素 点 值 的 操 作 。 经 过 平 滑 处 理 后 得 到

14、一 副 图 像 g(x,y), 其 表 达 式 如 下 ： snmfMyxg),( ),(/1),(式 中 ： x,y=0,1,2,， N-1； s 为 （ x， y） 点 领 域 中 点 的坐 标 的 集 合 ， 但 不 包 括 （ x， y） 点 ； M 为 集 合 内 坐 标 点 的 总 数 。领 域 平 均 法 有 力 地 抑 制 了 噪 声 ， 但 随 着 领 域 的 增 大 ， 图像 的 模 糊 程 度 也 愈 加 严 重 。 为 了 尽 可 能 地 减 少 模 糊 失 真 ， 也 可 采用 阈 值 法 减 少 由 于 领 域 平 均 而 产 生 的 模 糊 效 应 。 其 公 式

15、 如 下 ： 其 他),( ),(/1),(),(/1),( ),(yxf TnmfMyxfnmfMgsn s式 中 ： T 为 规 定 的 非 负 阈 值 。上 述 方 法 也 可 称 为 算 术 均 值 滤 波 器 ， 除 此 之 外 还 可 以 采 用 几 何均 值 滤 波 器 、 谐 波 均 值 滤 波 器 和 逆 谐 波 均 值 滤 波 器 。 几 何 均 值 滤 波器 所 达 到 的 平 滑 度 可 以 与 算 术 均 值 滤 波 器 相 比 ， 但 在 滤 波 过 程 中 会丢 失 更 少 的 图 像 细 节 。 谐 波 均 值 滤 波 器 对 “盐 ”噪 声 效 果 更好 ，

16、但 是 不 适 用 于 “胡 椒 ”噪 声 。 它 善 于 处 理 像 高 斯 噪 声 那 样 的其 他 噪 声 。 逆 谐 波 均 值 滤 波 器 更 适 合 于 处 理 脉 冲 噪 声 ， 但 它 有 个 缺点 ， 就 是 必 须 要 知 道 噪 声 是 暗 噪 声 还 是 亮 噪 声 ， 以 便 于 选 择 合 适 的滤 波 器 阶 数 符 号 ， 如 果 阶 数 的 符 号 选 择 错 了 可 能 会 引 起 灾 难 性 的 后果 。82.1.2 中 值 滤 波 器中 值 滤 波 是 一 种 常 用 的 去 除 噪 声 的 非 线 性 平 滑 滤 波 处 理方 法 ， 其 基 本 思

17、想 用 图 像 像 素 点 的 领 域 灰 度 值 的 中 值 来 代 替 该 像素 点 的 灰 度 值 。 二 维 中 值 滤 波 可 以 用 下 式 表 示 ：ijij fMedy式 中 ： A 为 滤 波 窗 口 ； ijf 为 二 维 数 据 序 列 。 其 主 要 功 能 是 让周 围 象 素 灰 度 值 的 差 比 较 大 的 像 素 改 取 与 周 围 的 像 素 值 接 近 的 值 ，从 而 可 以 消 除 孤 立 的 噪 声 点 ， 所 以 中 值 滤 波 对 于 滤 除 图 像 的 椒 盐 噪声 非 常 有 效 。 中 值 滤 波 器 可 以 做 到 既 去 除 噪 声 又

18、 能 保 护 图 像 的 边 缘 ，从 而 获 得 较 满 意 的 复 原 效 果 ， 而 且 ， 在 实 际 运 算 过 程 中 不 需 要 图 像的 统 计 特 性 ， 这 也 带 来 不 少 方 便 ， 但 对 一 些 细 节 多 ， 特 别 是 点 、 线 、尖 顶 细 节 较 多 的 图 像 不 宜 采 用 中 值 滤 波 的 方 法 。 如 果 希 望 强 调 中间 点 或 距 中 间 点 最 近 的 几 个 点 的 作 用 ， 则 可 采 用 加 权 中 值 滤 波 。 其基 本 原 理 是 改 变 窗 口 中 变 量 的 个 数 ， 可 以 使 一 个 以 上 的 变 量 等

19、于 同一 点 的 值 ， 然 后 对 扩 张 后 的 数 字 集 求 中 值 。 这 种 方 法 比 简 单 中 值 滤波 性 能 更 好 地 从 受 噪 声 污 染 的 图 像 中 恢 复 出 阶 跃 边 缘 以 及 其 他 细 节 。另 有 一 种 可 以 处 理 具 有 更 大 概 率 的 冲 激 噪 声 的 是 自 适 应 中 值 滤 波 器 ，在 进 行 滤 波 处 理 时 ， 能 依 赖 一 定 条 件 而 改 变 领 域 的 大 小 。 其 优 点 是在 平 滑 非 冲 激 噪 声 时 可 以 保 存 细 节 ， 所 以 既 能 除 去 “椒 盐 ”噪 声 ， 平 滑 其 他 非

20、 冲 激 噪 声 ， 还 能 减 少 诸 如 物 体 边 界 细 化 或 粗 化 等失 真 。2.2 频 域 低 通 滤 波在 分 析 图 像 信 号 的 频 率 特 性 时 ， 对 于 一 副 图 像 ， 直 流 分 量 表 示了 图 像 的 平 均 灰 度 ， 大 面 积 的 背 景 区 域 和 缓 变 部 分 是 低 频 分 量 ， 其边 缘 、 细 节 、 跳 跃 部 分 以 及 颗 粒 噪 声 都 代 表 图 像 的 高 频 分 量 。 频 域低 通 滤 波 就 是 除 去 其 高 频 分 量 就 能 去 掉 噪 声 ， 从 而 使 图 像 得 到 平 滑 。利 用 卷 积 定 理

21、， 可 以 写 成 以 下 形 式 ：9G(u,v)=H(u,v)F(u,v)式 中 ， F(u,v)是 含 噪 图 像 的 傅 立 叶 变 换 ， G(u,v)是 平 滑 后 图像 的 傅 立 叶 变 换 ， H (u,v)是 传 递 函 数 。 利 用 H(u,v)使F(u,v)的 高 频 分 量 得 到 衰 减 ， 得 到 G(u,v)后 再 经 过 反 变 换 就 得到 所 希 望 的 图 像 g(u,v)了 。 低 通 滤 波 平 滑 图 像 的 系 统 框 图 如 下 所示 ：图 3-1 图 像 频 域 低 通 滤 波 流 程 框 图低 通 滤 波 法 又 分 为 以 下 几 种

22、：(1) 理 想 低 通 滤 波 器 (ILPF)一 个 理 想 的 低 通 滤 波 器 的 传 递 函 数 由 下 式 表 示 ： 0),(01),(DvuvuH式 中 D0是 一 个 规 定 的 非 负 的 量 ， 它 叫 做 理 想 低 通 滤 波器 的 截 止 频 率 。 D(u,v)代 表 从 频 率 平 面 的 原 点 到 （ u ,v） 点的 距 离 ， 即 ：2/1)(),vu理 想 低 通 滤 波 器 在 处 理 过 程 中 会 产 生 较 严 重 的 模 糊 和 振 铃 现 象 。(2) 巴 特 沃 思 低 通 滤 波 器 （ BLPF）n 阶 巴 待 沃 思 滤 波 器

23、的 传 递 系 数 为 nvuDvuH20),(/1),(DLPF与 ILPF不 同 ， 它 的 通 带 与 阻 带 之 间 没 有 明 显 的 不 连 续 性 ，10因 此 它 没 有 “振 铃 ”现 象 发 生 ， 模 糊 程 度 减 少 ， 但 从 它 的 传 递 函数 特 性 曲 线 H(u,v)可 以 看 出 ， 在 它 的 尾 部 保 留 有 较 多 的 高 频 ， 所以 对 噪 声 的 平 滑 效 果 还 不 如 ILPF（ 理 想 低 通 滤 波 器 ） 。(3) 指 数 滤 波 器 (ELPF)其 传 递 函 数 表 示 为 ：nvuDevuH),(/0),(由 于 ELPF

24、 具 有 比 较 平 滑 的 过 渡 形 ， 为 此 平 滑 后 的 图 像 没 有“振 铃 ”现 象 ， 而 ELPF 与 BLPF 相 比 它 具 有 更 快 的 衰 减 特 性 ，所 以 经 ELPF 滤 波 的 图 像 比 BLPF 处 理 的 图 像 稍 微 模 糊 一 些 。(4) 梯 形 滤 波 器 （ TLPF)梯 形 滤 波 器 的 传 递 函 数 介 于 理 想 低 通 滤 波 器 和 具 有 平滑 过 渡 带 的 低 通 滤 波 器 之 间 ， 它 的 传 递 函 数 为 ：0101 1),(0/),( ,),( DvuDvuH式 中 ： D0 为 梯 形 低 通 滤 波

25、 器 截 止 频 率 ， D0、 D1 须 满 足D0D1， 它 的 性 能 介 于 ILPF 和 BLPF 之 间 ， 对 图 像 有 一 定 的 模糊 和 振 铃 效 应 。第 三 章 、 图 像 平 滑 处 理 与 调 试本 课 程 设 计 中 程 序 运 行 的 环 境 是 windows平 台 ， 并 选 用MATLAB作 为 编 程 开 发 工 具 ， MATLAB是 一 种 向 量 语 言 ， 它 非 常 适合 于 进 行 图 像 处 理 。3.1 模 拟 噪 声 图 像图 像 增 强 操 作 主 要 是 针 对 图 像 的 各 种 噪 声 而 言 的 ， 为 了 说 明 图像

26、处 理 中 的 滤 波 方 法 和 用 途 ， 需 要 模 拟 数 字 图 像 的 各 种 噪 声 来 分 析滤 波 效 果 。11MATLAB图 像 处 理 工 具 箱 提 供 的 噪 声 添 加 函 数 imnoise,它可 以 对 图 像 添 加 一 些 典 型 的 噪 声 。 其 语 法 ：J=imnoise(I,type)J=imnoise(I,type,parameters)其 功 能 是 ： 返 回 对 原 图 像 I添 加 典 型 噪 声 的 图 像 J,参 数type和 parameters用 于 确 定 噪 声 的 类 型 和 相 应 的 参 数 。三 种 典 型 的 噪

27、声 ：type=gaussian时,为 高 斯 噪 声；type=salt %读 取 图 像I1=imnoise(I,gaussian); %加 高 斯 噪 声I2=imnoise(I,salt %加 椒 盐 噪 声I3=imnoise(I,speckle); %加 乘 性 噪 声subplot(221),imshow(I); %显 示 图 像 Ititle(原 图 像 );subplot(222),imshow(I1);title(高 斯 噪 声 污 染 的 图 像 );subplot(223),imshow(I2);title(椒 盐 噪 声 污 染 的 图 像 );subplot(224

28、),imshow(I3);title(乘 性 噪 声 污 染 的 图 像 ); 12运 行 结 果 如 下 ：图 3-1 噪声污染的图像3.2 均 值 滤 波 法在 MATLAB 图 像 处 理 工 具 箱 中 ， 提 供 了 imfilter 函 数 用于 实 现 均 值 滤 波 ， imfilter 的 语 法 格 式 为 ：B=imfilter(A,H)其 功 能 是 ， 用 H 模 板 对 图 像 A 进 行 均 值 滤 波 ，取 平 均 值 滤 波 模 版 为H1=1/91 1 1;1 1 1;1 1 1；H2=1/251 1 1 1 1； 1 1 1 1 1； 1 1 1 1 1；

29、 1 1 1 1 1； 1 1 1 1 1；分 别 以 这 两 个 平 均 值 滤 波 算 子 对 图 3-1中 的 四 幅 图 像 进 行 滤波 操 作 。13取 H1,程 序 如 下 ：I=imread(lena.bmp);I1=imnoise(I,gaussian);I2=imnoise(I,salt I3=imnoise(I,speckle);H1=ones(3,3)/9; %33领 域 模 板J=imfilter(I,H1); %领 域 平 均J1=imfilter(I1,H1);J2=imfilter(I2,H1);J3=imfilter(I3,H1);subplot(221),i

30、mshow(J);subplot(222),imshow(J1);subplot(223),imshow(J2);subplot(224),imshow(J3);运 行 结 果 如 图 3-2取 H2， 程 序 如 下 ：RGB=imread(h:赵 云 -真 三 国 六 .jpg);I=rgb2gray(RGB);I1=imnoise(I,gaussian);I2=imnoise(I,salt I3=imnoise(I,speckle);H2=ones(5,5)/25; %55领 域 模 板J=imfilter(I,H2); %领 域 平 均J1=imfilter(I1,H2);J2=imf

31、ilter(I2,H2);J3=imfilter(I3,H2);subplot(221),imshow(J);14title(原 图 像 滤 波 后 );subplot(222),imshow(J1);title(高 斯 污 染 图 像 滤 波 后 );subplot(223),imshow(J2);title(椒 盐 污 染 图 像 滤 波 后 );subplot(224),imshow(J3);title(乘 法 污 染 图 像 滤 波 后 ); 运 行 结 果 如 图 3-3：图3-2 图3-1中图像经过平均值算子H1滤波后图像15图 3-3 图 3-1 中图像经过平均值算子 H2 滤波

32、后图像比 较 处 理 后 的 图 像 结 果 可 知 ， 领 域 平 均 处 理 后 ， 图 像 的 噪 声 得到 了 抑 制 ， 但 图 像 变 得 相 对 模 糊 ， 对 高 斯 噪 声 的 平 滑 效 果 比 较 好 。领 域 平 均 法 的 平 滑 效 果 与 所 选 用 的 模 板 大 小 有 关 ， 模 板 尺 寸 越 大 ，则 图 像 的 模 糊 程 度 越 大 。 此 时 ， 消 除 噪 声 的 效 果 将 增 强 ， 但 同 时 所得 到 的 图 像 将 变 得 更 模 糊 ， 图 像 细 节 的 锐 化 程 度 逐 步 减 弱 。3.3 中 值 滤 波 法MATLAB图 像

33、 处 理 工 具 箱 提 供 了 medfilt2函 数 用 于 中 值 滤 波 。其 语 法 格 式 为 ：B=medfilt2(A)其 功 能 为 ： 用 33的 滤 波 窗 口 对 图 像 A进 行 中 值 滤 波 ；B=medfilt2(A,m,n)其 功 能 是 ： 用 大 小 为 mn的 窗 口 对 图 像 A进 行 中 值 滤 波 ；B=medfilt2(A,indexed,.)其 功 能 为 ： 对 索 引 图 像 A进 行 中 值 滤 波 ；16可 运 行 以 下 程 序 实 现 ：RGB=imread(h:赵 云 -真 三 国 六 .jpg);I=rgb2gray(RGB)

34、;I1=imnoise(I,gaussian);I2=imnoise(I,salt I3=imnoise(I,speckle);J1=medfilt2(I1,3,3); %33中 值 滤 波 模 板J2=medfilt2(I2,3,3);J3=medfilt2(I3,3,3);J4=medfilt2(I1,5,5); %55中 值 滤 波 模 板J5=medfilt2(I2,5,5);J6=medfilt2(I3,5,5);figure,subplot(121),imshow(J1);title(高 斯 3*3中 值 滤 波 );subplot(122),imshow(J2);title(椒

35、盐 3*3中 值 滤 波 );figure,subplot(121),imshow(J3);title(乘 法 3*3中 值 滤 波 );subplot(122),imshow(J4);title(高 斯 5*5中 值 滤 波 );figure,subplot(121),imshow(J5);title(椒 盐 5*5中 值 滤 波 );subplot(122),imshow(J6);title(乘 法 5*5中 值 滤 波 );运 行 结 果 如 下 ：1718图 3-4 受到高斯、椒盐及乘法噪声污染的图像经不同模版的中值滤波后的图像由 图 3-4可 知 ， 此 方 法 能 够 非 常 好

36、地 将 椒 盐 噪 声 去 除 掉 ， 可 见中 值 滤 波 方 法 对 于 椒 盐 噪 声 或 脉 冲 式 干 扰 具 有 很 强 的 滤 除 作 用 ， 但对 于 高 斯 和 乘 性 噪 声 效 果 不 佳 。 与 图 3-2、 图 3-3相 比 ， 当 噪声 为 椒 盐 噪 声 时 ， 中 值 滤 波 器 的 效 果 比 均 值 滤 波 好 。 因 为 这些 干 扰 值 与 其 邻 近 像 素 的 灰 度 值 有 很 大 的 差 异 ， 经 过 排 序 后 取 中 值的 结 果 就 将 此 干 扰 强 制 变 成 与 其 邻 近 的 某 些 像 素 值 一 样 ， 从 而 达 到去 除

37、干 扰 的 效 果 。 但 是 由 于 中 值 滤 波 方 法 在 处 理 过 程 中 会 带 来 图 像模 糊 ， 所 以 对 于 细 节 丰 富 ， 特 别 是 点 、 线 和 尖 顶 细 节 较 多 的 图 像 不适 用 。3.4 频 域 低 通 滤 波 法频 域 低 通 滤 波 处 理 噪 声 图 像 的 方 法 如 下 :首 先 构 建 二 维 滤 波 器 d;19f1,f2=freqspace(25,meshgrid);Hd=zeros(25,25);d=sqrt(f1.2+f2.2)0.5; %0.5为 截 止 半 径 大 小Hd(d)=1;h=fsamp2(Hd);figure

38、,freqz2(h,64,64); 图 3-5 用频率采样法构建的二维滤波器然 后 用 所 构 建 的 二 维 滤 波 器 对 以 上 图 像 进 行 滤 波RGB=imread(h:赵 云 -真 三 国 六 .jpg);I=rgb2gray(RGB);I1=imnoise(I,gaussian);I2=imnoise(I,salt I3=imnoise(I,speckle);J=imfilter(I,h,replicate);J1=imfilter(I1,h,replicate);J2=imfilter(I2,h,replicate);J3=imfilter(I3,h,replicate);

39、subplot(221),imshow(J);title(原 图 像 滤 波 后 );20subplot(222),imshow(J1);title(高 斯 污 染 图 像 滤 波 后 );subplot(223),imshow(J2);title(椒 盐 污 染 图 像 污 染 后 );subplot(224),imshow(J3);title(乘 法 污 染 图 像 滤 波 后 );其 运 行 结 果 如 下 图 :图3-6 图3-1中图像经过二维滤波器 h滤波后图像频域低通滤波具有更好的选择性，对噪声在一定范围内也可以起到抑制作用，同时也对图像的边缘细节和高频信息分量有更好的保持作用。使

40、得图像在轮廓上显得更清晰。21第 四 章 、 总 结 与 体 会图像平滑是为了人类视觉系统的生理接受特点而设计的改善图像质量的方法,有很强的主观性。不同的平滑算法适应于特定类型的噪声模型,实际应用中应该根据实际图像中包含的噪声情况灵活地选取适当的平滑算法。文中对噪声模型及分类，图像平滑处理的两大方面即空间域和频率域,以及两种处理方向里的几种处理方法进行了介绍，并对一种常用的简单平滑算法进行了分析。本文结合不同的噪声类型，分析讨论了均值滤波、中值滤波和频域低通滤波3种平滑滤波技术在去几种典型噪声方面性能的优劣和视觉效果的对比，对在图像平滑处理方面提供了一定的参考。通过本次设计使我加深了对图像处理

41、尤其是图像平滑的认识及其应用，同时感谢在此设计中帮助我的老师和同学参 考 文 献1 杨 帆 ， 等 数 字 图 像 处 理 与 分 析 北 京 ： 北 京 航 空 航 天 大 学 出 版 社 ，2007.102 闫 娟 .数 字 图 像 的 平 滑 处 理 方 法 研 究 .软 件 导 刊 ， 2009.013 平 丽 .图 像 平 滑 处 理 方 法 的 比 较 研 究 .信 息 技 术 ， 2010.014 梁 一 江 .图 像 平 滑 处 理 方 法 初 探 及 简 单 的 算 法 介 绍 .才 智 ，2009.045李 骞 ,陈占伟.图像边缘检测新技术及其应用.许昌学院学报,2006.36求是科技编著.MATLAB7.0从入门到精通.人民邮电出版社,20067程正兴.小波分析算法与应用,西安交通常大学出版社,19978章毓晋.图像处理和分析 ，清华大学出版社，19999王宇生等.一种基于积分变换的边缘检测方法,中国图形图像学报,200510王晓明,顾晓东,刘健.基于张量的图像边缘检测及滤波,中国图形图像学报,20082211张英琦,何洪涛， 一种基于数学形态学的灰度图象边缘检测方法 ，西南交通大学学报，2006