%------------------------------------------------------------------
% exa120602_db97, 例12.6.2
% 本程序利用提升方案實(shí)現(xiàn)db9/7小波濾波器，并用該小波進(jìn)行圖像的分解
%------------------------------------------------------------------
clear all;

format long;
pho0=double(imread('lena.bmp')); %讀取圖像
[h,v]=size(pho0);
pho1=zeros(h,v);
pho2=zeros(h,v);
hh=h;

subplot(2,2,1);colormap(gray);imagesc(pho0);axis off; %顯示原圖像
photemp=pho0;
for lever=1:3 %共進(jìn)行三級(jí)分解
h=floor(hh/2^(lever-1)); v=h;
for i=1:h %對(duì)每一行進(jìn)行分解    
[ca,cd]=lifting_db97(photemp(i,1:h));
pho1(i,1:h)=[ca,cd];
end

for i=1:v %對(duì)每一列進(jìn)行分解
[ca,cd]=lifting_db97(pho1(1:v,i));
pho2(1:v,i)=[ca,cd]';
end
subplot(2,2,lever+1);colormap(gray);imagesc(pho2);axis off; %顯示每一級(jí)分解后的圖像
photemp=pho2;
end