www.pudn.com > m_program.zip > nstdhaardemo.m

function nstdhaardemo(imgname) 
%  nstdhaardemo haar小波三级非标准规格化分解与重构演示程序 
%  nstdhaardemo(imgname) 
%  本程序作用:载入图像文件,显示图像3级非标准haar小波规格化分解与重构过程 
%  输入: imgname   ----要装载的图像名称(真彩色,灰度图,索引图) 
%  缺省图像为color256.png 
 
if nargin==0 
    imgname='color256.png'; 
end 
%  读入的X中含有被装载的图像信号,map中含有被装载的color 
[X,map] = imread(imgname); 
if ndims(X)==3 
    imgcolor=1; 
else 
    imgcolor=0; 
end 
X=double(X); 
h=size(X,1); 
 
%画出原始图像 
figure(1); 
subplot(4,4,1);imshow(uint8(X),map); 
title(' haar小波3级非标准规格化分解过程'); 
xlabel('原始图像'); 
axis square 
 
position=2; 
A=X; 
%依次序画出3级非标准规格化分解过程中对应的行变换,列变换图像 
for i=1:3 
    A011=decompose(A,imgcolor,1,1,h); 
    subplot(4,4,position);imshow(uint8(A011),map);axis square 
    A012=decompose(A011,imgcolor,1,2,h); 
    subplot(4,4,position+4);imshow(uint8(A012),map);axis square 
    xlabel(strcat(int2str(i),'级分解图像')); 
    A=A012; 
    h=h/2; 
    position=position+5; 
end 
xlabel('haar小波3级非标准规格化分解图象'); 
 
%画出3级分解图像 
figure(2); 
subplot(4,4,1);imshow(uint8(A),map); 
title('haar小波3级非标准规格化分解重构过程'); 
xlabel('3级分解图像'); 
axis square 
position=2; 
 
%依次序画出3级非标准规格化分解后重构过程中对应的行变换,列变换图像 
if imgcolor 
    h=size(A(:,:,1),1)/4; 
else 
    h=size(A,1)/4; 
end 
for i=1:3 
    RX11=reconstruct(A,imgcolor,1,1,h); 
    subplot(4,4,position);imshow(uint8(RX11),map);axis square 
    RX12=reconstruct(RX11,imgcolor,1,2,h); 
    subplot(4,4,position+4);imshow(uint8(RX12),map);axis square 
    xlabel(strcat(int2str(3-i),'级重构图像')); 
    A=RX12; 
    h=h*2; 
    position=position+5; 
end 
xlabel('haar小波3级非标准分解后重构图像'); 
figure(1); 
 
%-------------------------------------------------- 
%        内部程序 
%-------------------------------------------------- 
function C=decompose(A,imgcolor,level,roworcol,h) 
 
if imgcolor 
    Ar=A(:,:,1); 
    Ag=A(:,:,2); 
    Ab=A(:,:,3); 
    Cr=nstdhaardec2(Ar,level,roworcol,h); 
    Cg=nstdhaardec2(Ag,level,roworcol,h); 
    Cb=nstdhaardec2(Ab,level,roworcol,h); 
    C(:,:,1)=Cr; 
    C(:,:,2)=Cg; 
    C(:,:,3)=Cb; 
else 
    C=nstdhaardec2(A,level,roworcol,h); 
end 
 
%-------------------------------------------------- 
%        内部程序 
%-------------------------------------------------- 
function C=reconstruct(A,imgcolor,level,roworcol,h) 
 
if imgcolor 
    Ar=A(:,:,1); 
    Ag=A(:,:,2); 
    Ab=A(:,:,3); 
    Cr=nstdhaarrec2(Ar,level,roworcol,h); 
    Cg=nstdhaarrec2(Ag,level,roworcol,h); 
    Cb=nstdhaarrec2(Ab,level,roworcol,h); 
    C(:,:,1)=Cr; 
    C(:,:,2)=Cg; 
    C(:,:,3)=Cb; 
else 
    C=nstdhaarrec2(A,level,roworcol,h); 
end