dbwt.m

2d wavelet transform

function dbwt(level)

h = [0   0.48296 0.83652 0.22414 -0.12941]
g = [0   -0.12941    -0.22414    0.83652 -0.48296]

x = imread('AG.jpg');
y = rgb2gray(x);
z = double(y(268:523,123:378));
%imshow(z);
[m,n] = size(z)

for l = 1:level

for i = 1:m
    th = conv(z(i,:),h); 
    sampth = th(1:2:n);
    tg = conv(z(i,:),g);
    samptg = tg(1:2:n);
    transH(i,1:n/2)=sampth;
    transH(i,(n/2+1):n)=samptg;
end

for j = 1:n
    th = conv(transH(:,j),h); 
    sampth = th(1:2:m);
    tg = conv(transH(:,j),g);
    samptg = tg(1:2:m);
    transH(1:m/2,j)=sampth;
    transH((m/2+1):m,j)=samptg;
end
    z = transH;
    [m,n] = size(z);
end

transH = transH./max(max(transH));

imshow(transH);