基于提升算法的二维53和97小波变换的MATLAB仿真DSP实现.pdfVIP

基于提升算法的二维5/3 和9/7 小波变换的MATLAB 仿真与DSP 实现 王靖琰，刘蒙 中国科学院上海应用物理研究所，上海 (201800) E-mail：wjycas@ 摘 要：本文讨论了基于提升算法的二维 5/3 和 9/7 小波的原理，对算法进行了MATLAB 仿真，并在浮点型DSP TMS320C6713B 上实现了图像的二维5/3、9/7 小波提升变换和逆变 换。实验结果证明了方法的有效性。 关键词：小波提升，二维9/7、5/3 小波，MATLAB ，TMS320C6713B 1．引言 随着人们对多媒体信息需求的日益增长，数码相机、移动电话、MP4 等多媒体信息处 理系统蓬勃发展。基于通用DSP 处理器的此类系统设计以灵活性强、扩展性好、可升级和 易维护的优点成为系统开发的首选方案 [1] 。 由于良好的时频局部特性和多分辨分析特性，小波已广泛应用于图像处理领域，并且被 吸收进新的一些国际标准中成为了标准算法。文中在MATLAB平台上对基于小波提升的二 维离散5/3和9/7小波变换算法进行了仿真，并在浮点型DSP TMS320C6713B上实现了算法， 该程序运算速度快，可充分利用硬件资源，特别适用于嵌入式系统的需求。 2 ．小波变换提升算法基本原理 1994 年 Sweldens 提出了小波的提升算法，有效地解决传统的基于 Mallat 的塔式分解小 波变换算法计算量大、对存储空间的要求高的问题，从算法方面提高了小波变换的实现效率 [2] 。 2.1 5/3小波提升格式 小波提升算法的基本思想是通过由基本小波(lazy wavelet)逐步构建出一个具有更加良 好性质的新小波，其实现步骤有 3 个：分解(split)、预测(predict)和更新(update) 。分解是将 数据分为偶数序列和奇数序列 2 个部分，预测是用分解的偶数序列预测奇数序列，得到的预 测误差为变换的高频分量，更新是由预测误差来更新偶数序列，得到变换的低频分量。在 J PEG2000 中，5/3 提升小波变换的算法为[3] ： ⎡x (2n ) +x (2n +2) ⎤ c(2n 1) x(2n 1) (1) + =+ −⎢⎣ 2 ⎥⎦ ⎡c(2n −1) +c(2n +1) +2⎤ d (2n) x(2n) =+⎢⎣ 4 ⎥⎦ (2) 由其正变换的反置即可得到逆变换的算法为 ⎡c(2n-1) + c(2n+1)+2⎤ x (2n) = d (2n) - ⎢⎣ 4 ⎥⎦ (3) ⎡x(2n)+x(2n+2) ⎤ x(2n+1)=c(2n)+ ⎢⎣ 2 ⎥⎦ (4) 从算式可以得出，提升算法是原位计算，即进行小波变换时在原位计算各个系数，计算 -1- 的系数可以直接替代原始数据而不需要附加数据存储空间。 2.2 9/7小波提升格式 Daubechies 9 /7 双正交小波基，具有线性相位，消失矩较大，能量集中性好等特性，在 图像处理领域有广泛的应用。图像经过9 /