yysuo.m

数字图像处理

kk=imread('ziji7.jpg');
X=rgb2gray(kk);
X=double(X);
figure();
axis square
wname='haar';lev=3;
[c,s]=wavedec2(X,lev,wname);
[thr,nkeep]=ddencmp('cmp','wv',X);
xd=wdencmp('gbl',c,s,wname,lev,thr,'s',1);
figure(1);
subplot(1,2,1);
imshow(X,[]);
title('原始信号');
subplot(1,2,2);
imshow(xd,[]);
title('压缩后的信号');