%------------------------------------------------------------------
% exa120602_db53, 例12.6.2
% 本程序利用提升方案實現db5/3小波濾波器，并用該小波進行圖像的分解
%------------------------------------------------------------------
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 %共進行三級分解
h=floor(hh/2^(lever-1)); v=h;
for i=1:h %對每一行進行分解    
[ca,cd]=lifting_db53(photemp(i,1:h));
pho1(i,1:h)=[ca,cd];
end

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