ezw.m

ezw编码

%EZW编码主程序
Im=dwt;
%%选取阈值
MaxIm=double(max(max(abs(Im))));
T=zeros(1,5);%5级编码
T(1)=2^floor(log2(MaxIm));
for order=2:5
    T(order)=T(order-1)/2;
end
sublist=[];
sub_list=[];
[r,c]=size(Im);
A=mapping(r);
[m,n]=size(A);
global N;% 采用morton扫描，并建立全局矢量N
N=zeros(m*n,2);
     for i=1:m,
         for j=1:n,
             N(A(i,j),1)=i;
             N(A(i,j),2)=j;
         end
     end
global C;
C=zeros(m,n);% 建立0矩阵判断重要系数与否
disp('There are the EZW encoding outputs below')
order=1;
while order<=5,
[D,Im,sublist,sub_list] = mainpass(Im,T(order),sublist,sub_list);%主扫描
D
flag{order}=D;
significantlist{order}=sub_list;
T(order)=T(order)/2;
S = assistantpass(sublist,T(order));
S
SS{order}=S;
order=order+1;
end