嵌入式小波零树编码压缩(3)

（2）图像3级分解、10级编码、8级解码下的效果图。

ezw(Im64,3,10,8);

分解矩阵DecIm的部分数据

40.9531 87.2344 94.2188 161.2031 64.3125 40.4063 37.5781 65.6875

110.8281 129.3125 91.5625 164.3750 107.1719 50.5469 34.9688 75.2188

150.7813 134.5156 98.0625 164.9688 121.7344 67.5469 28.0938 77.1250

132.9219 117.7656 91.0938 155.3125 117.0781 53.0781 28.1563 83.3750

119.3281 118.5156 95.1875 137.1719 105.8750 30.3906 33.6094 79.6406

100.8906 102.9375 119.2031 132.8438 71.5000 28.6875 38.5313 75.3281

103.2813 85.3750 102.6250 107.4531 31.3594 28.2656 45.3438 79.7813

98.4844 120.0938 134.8438 73.3281 29.1719 33.6094 50.3906 86.

0781

解码矩阵DecodeMat的对应数据

40.7500 87.2500 94.2500 161.2500 64.2500 40.2500 37.7500 65.7500

110.7500 129.2500 91.7500 164.2500 107.2500 50.7500 34.7500 75.2500

150.7500 134.7500 98.2500 164.7500 121.7500 67.7500 28.2500 77.2500

132.7500 117.7500 91.2500 155.2500 117.2500 53.2500 28.2500 83.2500

119.2500 118.7500 95.2500 137.2500 105.7500 30.2500 33.7500 79.7500

100.7500 102.7500 119.2500 132.7500 71.7500 28.7500 38.7500 75.2500

103.2500 85.2500 102.7500 107.2500 31.2500 28.2500 45.2500 79.7500

98.2500 120.2500 134.7500 73.2500 29.2500 33.7500 50.2500 86.2500

运行效果：

EZW算法的过程详解和Matlab代码（1）构建扫描次序表（修正小波树结构）

收藏

前段时间，我们讨论了嵌入式小波零树算法的基本原理。（http://blog.csdn.net/chenyusiyuan/archive/2007/11/16/1888968.aspx）一个多星期过去了，我根据自己对算法的理解，编写出EZW算法的Matlab程序，可以实现图像的任意级别的小波分解和重构、以及任意精度的EZW编解码过程。下面，我们以一幅16*16的Lena局部图像经过3级小波分解后的小波数据为例详细的说明EZW算法的编解码过程，并给出相应的Matlab代码。原始图像如下

分解后的小波图像数据为：

一、在开始编码之前，首先要求出初始阈值T1：

MaxDecIm=max(max(abs(DecIm))); T=zeros(1,codeDim); T(1)=2^flor(log2(MaxDecIm));

二、然后是建立小波树结构，构建扫描次序表。这个扫描次序表非常重要，后面的编码、解码过程都要按照扫描次序表逐个处理数据矩阵的各个元素。构建过程如下：用（r，c）表示数据矩阵上各元素的位置。row、col作为全局变量，表示数据矩阵的行、列数。 1、小波树结构的特点：

（1）对于LL-N低频子带的点(r,c)，有3个孩子：(r,c+W)、(r+H,c)和(r+H,c+W)，其中W、H分别是LL-N子带的宽和高；

（2）第N~2高频子带(LH、HL、HH)的点都有4个孩子，即： tp=[2*r-1,2*c-1;2*r-1,2*c;2*r,2*c-1;2*r,2*c];

（3）第1高频子带的点没有孩子。根据小波树的这个特点，可编写如下“小波树”函数treeMat()，输入矩阵内任一点的位置（r，c），给出该点的子孙列表cp 。

function cp=treeMat(r,c) %这个函数是一个递归函数 global row col dim % dim是小波分解级数 HLL=row/2^dim; WLL=col/2^dim; if (r<=HLL)&&(c<=WLL)

tp1=[r,c+WLL;r+HLL,c;r+HLL,c+WLL];

cP=[tp1;treeMat(r,c+WLL);treeMat(r+HLL,c);treeMat(r+HLL,c+WLL)]; elseif (r>row/2)||(c>col/2) cP=[]; else

tp=[2*r-1,2*c-1;2*r-1,2*c;2*r,2*c-1;2*r,2*c]; tm1=[];tm2=[];tm3=[];tm4=[]; if (tp(4,1)<=row/2)&&(tp(4,2)<=col/2) t1=treeMat(tp(1,1),tp(1,2)); tm1=[tm1;t1];

t2=treeMat(tp(2,1),tp(2,2)); tm2=[tm2;t2];

t3=treeMat(tp(3,1),tp(3,2));

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库嵌入式小波零树编码压缩(3)在线全文阅读。