untitled.m

嵌入式零树小波变换编码 EZW 的MATLAB 程序

clear
clc
X=imread('Girl.bmp');
% 对图像用小波进行层分解
[c,s]=wavedec2(X,2,'bior3.7');
%提取小波分解中的一层的低频和高频系数
ca1=appcoef2(c,s,'bior3.7',1); %提取低频系数
ch1=detcoef2('h',c,s,1); %提取高频系数
cv1=detcoef2('v',c,s,1);
cd1=detcoef2('d',c,s,1);
%小波重构
a1=wrcoef2('a',c,s,'bior3.7',1);
h1=wrcoef2('h',c,s,'bior3.7',1);
v1=wrcoef2('v',c,s,'bior3.7',1);
d1=wrcoef2('d',c,s,'bior3.7',1);
c1=[a1,h1,v1,d1];
ca1=appcoef2(c,s,'bior3.7',1);
% 保留小波分解第一层低频信息进行压缩
cal=wcodemat(ca1,400,'mat',0);
% 首先对第一层信息进行量化编码
% 改变图像高度
ca1=0.5*ca1;
ca2=appcoef2=(c,s,'bior3.7',2);
%%保留小波分解第二层低频信息进行压缩
ca2=wcodemat(ca2,400,'mat',0);
% 首先对第二层信息进行量化编码
%改变图像高度
ca2=0.25*ca2;
subplot(221)
% 显示原始图像
subplot(221)
imshow(X)
title('原始图像‘）
disp('原始图像的大小’）
whos('X')
subplot(222)
% 显示分频信息
c1=uint8(c1);
imshow(c1)
title('显示分频信息‘）
subplot(223)
disp('第一次压缩图像的大小’）
% 显示第一次压缩的图像
ca1=uint8(ca1);
whos('ca1')
imshow(ca1);
title('第一次压缩的图像‘）
disp('第二次压缩图像的大小’）
subplot(224)
%显示第二次压缩的图像
ca2=uint8(ca2);
imshow(ca2);
title('第二次压缩的图像‘）
whos('ca2')