test.m

EZW implementation which will use wavelet transforms and do the compression decompression using EZW

clear all
close all
X=[63 -34 49 10 7 13 -12 7;
   -31 23 14 -13 3 4 6 -1;
   15 14 3 -12 5 -7 3 9;
   -9 -7 -14 8 4 -2 3 2;
   -5 9 -1 47 4 6 -2 2;
   3 0 -3 2 3 -2 0 4;
   2 -3 6 -4 3 6 3 6;
   5 11 5 6 0 3 -4 4];

[DD,L1,SS,L2]=ezw_encode(X,6);
[encode_x,h,sortindex]=huffman(DD);
encode_bin_x=change(encode_x);
dlmwrite('test_encode1.txt',encode_bin_x,' ');
dlmwrite('test_encode2.txt',SS,' ');

s=size(encode_bin_x,2)+size(SS,2)
[m,n]=size(X);
t=m*n*8
ratio=s/t

decode_x=ihuffman(encode_x,h,sortindex);
ezw_decode_x=ezw_decode(X,6,decode_x,L1,SS,L2)
E=X-ezw_decode_x