nstdhaardec2.m

用matlab编写的压缩解压缩程序

function [a,lt] = nstdhaardec2(x,level,rorc,h)
%nstdhaardec2 二维haar小波非标准规格化分解程序(多级分解)
%作用:使用haar小波对每一行和每一列像素值进行小波变换
%[a,lt]=nstdhaardec2(x,level,rorc)
%输入:x      载入的二维图像像数值; 
%     level  小波分解次(级)数设定值(如果设定值超过最高可分解次数,按最高分解次数分解)
%     rorc   作行变换(1)或列变换(2),缺省为0,行列变换都做
%     h      分解的矩阵块大小,缺省为整个x矩阵的变换
%输出:a      分解后数值矩阵,大小与x相同,
%     lt     记录实际分解次数

a=double(x);
t=1;
sX=size(x);
if nargin==1
    level=1;
    h=sX(2);
    rorc=0;
else
    if nargin==2
        h=sX(2);
        rorc=0;
    else
        if nargin==3
            h=sX(2);
        end
    end            
end

lt=level;
while and(h>1,t<=level)
    if rorc==1
        for row=1:h
            a(row,:)=haardec(a(row,:),h);
        end
    else
        if rorc==2
            for col=1:h
                temp=haardec(a(:,col)',h);
                a(:,col)=temp';
            end
        else
            for row=1:h
                a(row,:)=haardec(a(row,:),h);
            end
            for col=1:h
                temp=haardec(a(:,col)',h);
                a(:,col)=temp';     
            end
        end
    end
    h=h/2;
    t=t+1;
end
if and(h<=1,lt~=t-1)
    lt=t-1;
end
%--------------------------------------------------
%        内部程序
%--------------------------------------------------
function y=haardec(c,h)
%haardec 1-D haar decompose program
%y=haardec(c,h)
%输入:c     载入的二维图像像数值; 
%     h     载入的要分解的矩阵大小
%输出:y     分解后数值矩阵,大小与x相同,

y=c;
sqrt2=sqrt(2);
h=h/2;
for i=1:h
    y(i)=(c(2*i-1)+c(2*i))/sqrt2;
    y(h+i)=(c(2*i-1)-c(2*i))/sqrt2;        
end