小波变换算法应用

软软 件件 开开 发发 课程设计课程设计 题目 小波算法的设计题目 小波算法的设计 题目要求 将小波算法在 MATLAB 中实现 并将其应用于数字图像处理中 学学 院 数学学院院 数学学院 专业班级 应用数学专业班级 应用数学 09 209 2 班班 姓姓 名 李名 李 明明 学学 号号 2009631220096312 指导教师 邢燕 何蕾指导教师 邢燕 何蕾 2013 3 5 小波算法的设计小波算法的设计 1 小波变换背景小波变换背景 小波变换是当前应用数学中一个迅速发展的领域 是分析和处理非平稳信号的一种有 力 工具 它是以局部化函数所形成的小波基作为基底而展开的 具有许多特殊的性能和优点 小波分析是一种更合理的时频表示和子带多分辨分析 对它的研究开始于 20 世纪 80 年代 理论基础奠基于 20 世纪 80 年代末 经过十几年的发展 它已在信号处理与分析 地震信 号处理 信号奇异性监测和谱古迹 计算机视觉 语音信号处理 图像处理与分析 尤其 是图像编码等领域取得了突破性进展 成为一个研究开发的前沿热点 2 小波变换概念小波变换概念 小波变换是一窗口大小固定不变但其形状可改变的时频局部化分析方法 小波变换在 信号的高频部分 可以取得较好的时间分辨率 在信号的低频部分 可以取得较好的频率 分辨率 从而能有效地从信号 语音 图像等 中提取信息 设是平方可积分函数 即 则该连续函数的小波变换定义为 tf 2 RLtf dt a bt tf a baWTf 1 0 a 式中称为母小波 基本小波 生成的位移和尺度伸缩 其中 1 t a bt a ba t a 为尺度参数 b 为平移参数 连续小波变换有明确的物理意义 尺度参数 a 越大 则越宽 该函数的时间分 a t 辨率越低 前增加因子 a 1 是为了使不同的 a 下的能量相同 而 t ab t ab baWTf 在频域可以表示为 是幅频特性比 deF a baWT bj f 2 较集中的带通函数 小波变换具有表征分析信号 F 频域上局部性质的能力 采用不同 的 a 值做处理时 的中心频率和带宽都不同 但品质因数 中心频率 带宽 却不变 3 小波变换需求分析小波变换需求分析 小波分析的应用领域十分广泛 其中包括数学分析 信号分析 图像处理 量子 力学 理论物理 军事电子对抗与武器的智能化 计算机分类与识别 音乐与语言的人工 合成 医学成像与诊断 地震勘探数据处理和大型机械的故障诊断等方面 例如 在数学 方面 它已用于数值分析 构造快速数值方法 曲线曲面构造 微分方程求解 控制论等 领域 在信号分析方面 它主璧用于滤波 去噪声 压缩 传递等领域 在图像处理方面 它已应用于图像的压缩 分类 识别与诊断等领域 在医学成像方面 它已应用于减少 B 超 CT 核磁共振成像的时间 以及提高图像分辨率等领域 本课程设计要求将小波算法在 MATLAB 中实现 并将其应用于数字图像处理中 所以本 课程设计中主要将小波变换应用于图像处理的去噪 4 图像去噪常用函数图像去噪常用函数 通常处理的图像很多为索引图像 图像矩阵各元素表示的是调色板中的序号 而小波 分析是对数值进行分析的 因此要将索引图像进行编码 进行小波分析才有实际意义 MATLAB 提供了 wcodemat 函数来对图像进行编码 1 wcodemat 函数语法格式 Y wcodemat X NB 功能 对索引图像的数据矩阵 X 进行编码 Y 为编码返回值 NB 是最大编码值 决定了编码的范围是 0 NB 2 dwt2 函数实现二维离散小波变换 其语法格式 cA cH cV cD dwt2 X wname 其功能 使用指定的小波基函数 wname 对图像 X 进行二维离散小波变换 cA cH cV cD 分 别为图像分解的近似分量 水平分量 垂直分量和细节分量 3 idwt2 函数实现二维离散小波反变换 其语法格式 X idwt2 cA cH cV cD wname 其功能是 利用小波分解得到的 cA cH cV cD 分量进行二维离散小波反变换得到原始图像 wname 函数指定二维小波反变换采用的小波基函数 4 wavedec2 函数用于二维图像进行多层小波分解 其语法格式 C S wavedec X N wname 其功能 使用指定的小波基函数 wname 对图像 X 进行 N 层二维离散小波分解 5 利用小波变换进行图像去噪利用小波变换进行图像去噪 下面以二维小波分析用于图像去噪为例 对于二维图像信号的去噪方法同样适用于一 维信号 尤其对于几何图像更适合 二维模型可以表示为 jiejifjis 1 2 1 0 mji 其中 e 是标准偏差不变的高斯白噪声 幅度分布服从高斯分布 功率谱密度又是均匀分布 的噪声 二维信号用二维小波分析的去噪步骤如下 1 二维信号的小波分解 选择一个小波和小波分解的层次 N 然后计算信号 s 到第 N 层的分解 2 对高斯系当选进行阀值量化 对于从 1 到 N 的每一层 选择一个阀值 并对这 一层的高频系数进行软阀值量化处理 3 二维小波重构 根据小波分解的第 N 层的低频系数和经过修改的从第一层到第 N 层的各层高频系数计算二维信号的小波重构 值得注意的是 重点是如何选取阀值和阀值的量化 在 MATLAB 中的去噪函数有 ddendmp 和 wdencmp 等 MATLAB 运行结果 运行结果 6 6 结论结论 这次设计利用小波变换完成了对图像进去噪目的 基本上实现了设计的要求 图像去 噪是一个很有发展前途的研究领域 由于成像传感器噪声 相片颗粒噪声 图片在传输过 程中的通道传输误差等因素会使图片上出现一些随机的 离散的 孤立的像素点 这就是 图像噪声 图像噪声在视觉上通常与它们相邻的像素明显不同 例如黑区域中的白点 白 区域中的黑点等 与Fourier 变换相比 小波变换是空间 时间 和频率的局部变换 因而 能有效地从信号中提取信息 通过伸缩和平移等运算功能可对函数或信号进行多尺度的细 化分析 解决了 Fourier 变换不能解决的许多困难问题 因此利用小波变换进行图像的去 噪处理 有其独特的优势 此外 小波变换联系了应用数学 物理学 计算机科学 信号与信息处理 图像处理 地震勘探等多个学科 小波分析是一个新的数学分支 它是泛函分析 Fourier 分析 样 调分析 数值分析的完美结晶 小波分析是时间 尺度分析和多分辨分析的一种新技术 它在信号分析 语音合成 图像识别 计算机视觉 数据压缩 地震勘探 大气与海洋波 分析等方面的研究都取得了有科学意义和应用价值的成果 七 参考文献七 参考文献 1 谢平 王娜 林洪彬等编 数字信号处理 燕山大学 2007 年 3 月 2 余成波编 数字图像处理及 MATLAB 实现 重庆大学出版社 2007 年 8 月 3 薛年喜主编 MATLAB 在数字信号处理中的应用 清华大学出版社 2003 年 附录 源程序附录 源程序 load detfingr 装入图像 init 3718025452 下面进行噪声产生 randn seed init 设置一个种子 设置后下面的随机数是一定的 Xnoise X 18 randn size X 产生噪声 colormap map 显示原始图像及他的含噪声图像 subplot 2 2 1 image wcodemat X 192 title 原始图像X subplot 2 2 2 image wcodemat Xnoise 192 title 含噪声的图像Xnoise c s wavedec2 X 2 sym5 用sym5小波对图像信号进行二层的小波分解 下面进行图像的去噪处理 使用ddencmp函数来计算去噪的默认阀值和熵标准 使用wdencmp函