%LIFTING dla Haar i Spline
clc
close all
clear all


%% wczytanie danych - rgb to gray, im2double
image = imresize(double(rgb2gray(imread('bikes.bmp'))),[64,64]);
[n,m] = size(image);
im = image;

%Liczba poziomów rozkładu
lP = 5;

%% Haar dla danej liczby poziomów
for i = 1 : lP
    [d, im] = haarTransform(im,1);
    [wynik{i}.detale_1, wynik{i}.image_1] = haarTransform(im,2);
    [wynik{i}.detale_3, wynik{i}.detale_2] = haarTransform(d,2);
    im = wynik{i}.image_1;
end

%% Wyrysowanie piramidy
piramida = [wynik{i}.image_1 wynik{i}.detale_1;wynik{i}.detale_2 wynik{i}.detale_3]; 
for i=lP-1:-1:1
       piramida = [piramida wynik{i}.detale_1;wynik{i}.detale_2 wynik{i}.detale_3];
end
imshow(piramida,[])
%% Zerowanie wartości
    wynik = clearLevel(wynik,[2],[0 0 1]);
%% Transformata odwrotna
obraz = [];
for i = lP:-1:1
    obraz= inversHaarTransform(wynik{i},obraz);
end
figure
imshowpair(obraz,image,'montage')