baseonptwaveletanalyse.m

一些个 人收集的小波变换的程序

%%%%%%%【基于峰值变换小波系数分析主程序】%%%%%%%%%%
clear all
disp('Running PTWT....')
%%%%%%%%%%%%%%%【读入图像数据】%%%%%%%%%%%%%%%%%%%
load woman;  %  下载图像
X=im2double(X);
f=X;         %  原始图像
%%%%%%%%%%%%%%%【二维小波分解】%%%%%%%%%%%%%%%%%%%
[sa,sb]=size(f);
[Ain]=hvanalysis(f,1,1,sa,sb,-1);
[Ain]=hvanalysis(Ain,1,1,sa,sb,1);
%%%%%%%%%%%%%%%%【峰值变换】%%%%%%%%%%%%%%%%%%%%
Ain=func_PeakTransform(Ain);           %调用func_PeakTransform函数，此函数包括前向峰值变换和后向峰值变换
%%%%%%%%%%%%%%%%%【编译码】%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%【编码】%%%%%%%%%%%%%%%%%%%%%%%
rate=input('请输入rate=');  
disp('Running.......');
Enc=SPIHT_Enc(Ain,1,rate);
%%%%%%%%%%%%%%%%%【译码】%%%%%%%%%%%%%%%%%%%%%%%%
Dec= func_SPIHT_Dec(Enc);
%%%%%%%%%%%%%%%【重构源图像】%%%%%%%%%%%%%%%%%%%
[Ain]=hvsynthesis(Dec,1,1,sa,sb,1);
[A]=hvsynthesis(Ain,1,1,sa,sb,-1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
construct_pic=A;  %  重建全部图像
%%%%%%%%%%%%%%%%%%【计算PSNR】%%%%%%%%%%%%%%%%
mse=sum(sum((construct_pic-f).^2))/sa/sb;
psnr=10*log10(255*255/mse)
%%%%%%%%%%%%%%%%%%%%%%【结果显示】%%%%%%%%%%%
% figure;colormap(map);subplot(2,1,1);
% image(f);                           %  源图像显示
% title('original pic');
% subplot(2,1,2);
% image(abs(construct_pic));          %  重构源图像显示
% title('reconstructed pic');
disp('%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%');