yasuo11.m

数字图像处理

kk=imread('ziji6.jpg');
X=rgb2gray(kk);
X=double(X);
figure();
imshow(X,[]);
axis square
wname='haar';lev=3;
[c,s]=wavedec2(X,lev,wname);
alpha=1.5;m=3.5*prod(s(1,:));
[thr,nkeep]=wdcbm2(c,s,alpha,m);
xd=wdencmp('lvd',c,s,wname,lev,thr,'h');
figure(1);
subplot(1,2,1);
imshow(X,[]);
title('原始信号');
subplot(1,2,2);
imshow(xd,[]);
title('压缩后的信号');