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

tm3=[tm3;t3];

t4=treeMat(tp(4,1),tp(4,2)); tm4=[tm4;t4]; end

cP=[tp;tm1;tm2;tm3;tm4]; end

示例，当row=8，col=8，dim=2时，LL-N低频子带的点(1,1) (1,2) (2,1) (2,2) 的子孙分布如下：

由这个小波树列表tree_p，我们可以进一步构建函数childMat()，给出矩阵数据Mat和矩阵任一点的位置（r，c），返回该点的子孙数据列表chMat。

function chMat=childMat(Mat,chRows,chCols) global row col dim

chPoint=treeMat(chRows,chCols); chMat=[];

[mRows,mCols]=size(chPoint); for iRows=1:mRows

chMat=[chMat;chPoint(iRows,1),chPoint(iRows,2),Mat(chPoint(iRows,1),chPoint(iRows,2))]; end

2、构建扫描次序表

本文EZW算法的扫描次序为Morton式，其特征是从（1，1）开始，每4个点组成一个“Z”型扫描单元，从微观到宏观上都是严格的“Z”型结构，可以用递归方法来构建扫描次序表。扫描次序表scanlist由两部分组成，一个是将数据矩阵Mat按照morton扫描次序转换成数据列表matlist，一个是按照扫描次序组成的矩阵各点位置的（r，c）列表lsorder。

function scanlist=morton(Mat) global row col

matlist=mat2list(Mat);

scanorder=listorder(row,col,1,1); scanlist=[]; for i=1:row*col

scanlist=[scanlist;i scanorder(i,:) matlist(i)]; end

function mls=mat2list(Mat) % 该函数为递归函数 [r,c]=size(Mat); if (r==2)&&(c==2)

mls=[Mat(1,1);Mat(1,2);Mat(2,1);Mat(2,2)]; else

M1=Mat(1:r/2,1:c/2); M2=Mat(1:r/2,c/2+1:c); M3=Mat(r/2+1:r,1:c/2);

M4=Mat(r/2+1:r,c/2+1:c); lt1=mat2list(M1); lt2=mat2list(M2); lt3=mat2list(M3); lt4=mat2list(M4); mls=[lt1;lt2;lt3;lt4]; end

function lsorder=listorder(mr,mc,pr,pc) % 该函数为递归函数 lso=[pr,pc;pr,pc+mc/2;pr+mr/2,pc;pr+mr/2,pc+mc/2]; mr=mr/2;mc=mc/2;

lm1=[];lm2=[];lm3=[];lm4=[]; if (mr>1)&&(mc>1)

ls1=listorder(mr,mc,lso(1,1),lso(1,2)); lm1=[lm1;ls1];

ls2=listorder(mr,mc,lso(2,1),lso(2,2)); lm2=[lm2;ls2];

ls3=listorder(mr,mc,lso(3,1),lso(3,2)); lm3=[lm3;ls3];

ls4=listorder(mr,mc,lso(4,1),lso(4,2)); lm4=[lm4;ls4]; end

lsorder=[lso;lm1;lm2;lm3;lm4]; len=length(lsorder);

lsorder=lsorder(len-mr*mc*4+1:len,1:2);

嵌入式小波零树（EZW）算法的过程详解和Matlab代码（2）编码过程（081227-补充quantifier程序） 收藏

首先根据实际需要确定扫描次数D=codeDim，然后进行D次主扫描和辅扫描。

1、初始化

首先获取扫描次序表scanlist，表中每行的数字依次是该元素的（1）扫描序号、（2）对应于数据矩阵中位置的行号r、（3）对应于数据矩阵中位置的列号c、（4）值Mat（r，c），即

scanlist (i) = [ i, r, c, Mat(r,c) ]

编码扫描的对象就是次序表scanlist，而非数据矩阵Mat。然后构建初始符号矩阵flagMat，其中每个元素均为字符“Z”。

2、主扫描

在第d次（1<= d <=codeDim）扫描中，输入阈值为T（d），首先依次对次序表的每一行scanlist(i)按照如下步骤扫描：

（1） 如果上一次扫描得到的flagMat中，处于位置（r，c）的符号为“O”，则表示相应的Mat（r，c）是重要的，不需要通过本次扫描再次编码，故作跳过处理。

（2） 将scanlist(i,4)（即Mat(r,c)）的绝对值与阈值T(d) 比较，如果是重要的系数，则令flagMat(r,c)的值为P或N，并存入主扫描表scancode中，并将Mat(r,c)的值以及相应的行、列数（r，c）存入重

要数列表imptvalue中，将符号P或N存入重要数符号列表imptflag中。

（3） 如果是不重要的系数，则首先检查flagMat(r,c)，若是符号“X”，则表示这个点是不重要的，作跳过处理不再扫描。

（4） 如果点（r，c）处于第一分解级，即r>row/2或c>col/2，没有子孙系数，且其系数是不重要的，则该点flagMat的符号为Z，存入主扫描表scancode中。

（5） 如果点（r，c）满足r

至此，扫描结束。然后，将扫描后得到的符号矩阵flagMat转存到scanflag中，再将flagMat中的字符P、N都置换为O，字符X和T都置换为Z，使得flagMat的元素只有字符O、Z两种，作为输出到下一次扫描中使用。

function [imptvalue,imptflag,scancode,scanflag,flaglist]=mainscan(Mat,scanlist,flaglist,imptvalue,imptflag,threshold) global row col scancode=[]; for i=1:row*col

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