%**************************************************************************
% główna funkcja do testowania metod nieliniowej aproksymacji w wielu
% bazach
%**************************************************************************

clear all; close all;

measurements = [ 1000 5000:5000:65000];
GraphError =   [0.01   0.02  0.05  0.1   0.2   0.5   1    3   5    10  20  50  100];

X = double(imread('target.pgm'));

XX = image_approx(X, 1, 8, 0.1);

figure, imshow(X,[]);

figure, imshow(XX,[]);

fprintf('PSNR: %f dB\n', psnr(X,XX));
fprintf('SSIM: %f\n', ssim(uint8(X),uint8(XX)));