main.m

来自「EZW implementation which will use wavele」· M 代码 · 共 25 行

clear all
%close all

X=imread('lena.bmp');
%X=rgb2gray(m);
subplot(1,2,1);imshow(X,[]);title('原始图像');axis square;

de_x=haardec(X);
[DD,L1,SS,L2]=ezw_encode(de_x,6);
[encode_x,h,sortindex]=huffman(DD);
encode_bin_x=change(encode_x);
dlmwrite('encode1.txt',encode_bin_x,' ');%
dlmwrite('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(de_x,6,decode_x,L1,SS,L2);
re_x=haarrec(ezw_decode_x);
re_x=uint8(re_x);

subplot(1,2,2);imshow(re_x,[]);title('重构图像');axis square;