n = 512;
lw = 2;
fList(:,:,1) = rescale( load_image('regular3',n) );
fList(:,:,2) = rescale( load_image('regular2',n) );
fList(:,:,3) = rescale( load_image('regular1',n) );
% fList(:,:,3) = rescale( load_image('tv-image',n) );
% fList(:,:,4) = rescale( load_image('geometrical',n) );

% obrazy oryginalne
clf;
for i=1:3
    imageplot(fList(:,:,i),'', 2,2,i);
end

Jmin = 1;
options.h = compute_wavelet_filter('Daubechies',20);

err = [];
for i=1:size(fList,3)
    % 
    fW = perform_wavortho_transf(fList(:,:,i),Jmin,+1, options);
    e = compute_l2err(fW,1);
    err = [err e];
end
sel = 10:round(n*n/10);
Err = err(sel,:);
Err = Err ./ repmat(Err(1,:), [size(Err,1), 1]);

pause;
clf;
h = plot(log10(sel), log10(Err));
if using_matlab()  
    set(h, 'LineWidth', lw);
end
title('log_{10}( \epsilon^2[M]/||f||^2 )');
%legend('regular', 'phantom', 'lena', 'mandrill');
legend('regular3', 'regular2', 'tv-image', 'geometrical');
axis([log10(sel(1)) log10(sel(length(sel))) -7 0]);
if using_matlab()  
    box('on');
end