EZW压缩算法matlab

X=imread('cameraman.tif'); [cr,bpp] = wcompress('c',X,'cameraman.wtc','ezw','maxloop',10) Xc10 = wcompress('u','cameraman.wtc'); delete('cameraman.wtc') [cr,bpp] = wcompress('c',X,'cameraman.wtc','ezw','maxloop',11) Xc11 = wcompress('u','cameraman.wtc'); delete('cameraman.wtc') [cr,bpp] = wcompress('c',X,'cameraman.wtc','ezw','maxloop',12) Xc12 = wcompress('u','cameraman.wtc'); delete('cameraman.wtc') Xc10=uint8(Xc10); Xc11=uint8(Xc11); Xc12=uint8(Xc12); subplot(2,2,1), imagesc(X);colormap(gray); title('Original image') subplot(2,2,2), imagesc(Xc10);colormap(gray); title('compress image by ezw maxloop=10') subplot(2,2,3), imagesc(Xc11);colormap(gray); title('compress image by ezw maxloop=11') subplot(2,2,4), imagesc(Xc12);colormap(gray); title('compress image by ezw maxloop=12') PSNR10=calpsnr(X,Xc10); PSNR11=calpsnr(X,Xc11); PSNR12=calpsnr(X,Xc12); fprintf('SWT平稳小波滤波去噪PSNR=%d\n',PSNR10); fprintf('SWT平稳小波滤波去噪PSNR=%d\n',PSNR11); fprintf('SWT平稳小波滤波去噪PSNR=%d\n',PSNR12);