function xq = quant_laplace(x,R); 
% Quantizzatore di Laplace

N = 2^R;
sigma = std(x(:));
x = x/sigma;  % dati normalizzati

% livelli di restituzione per R =1,2,3,4 bit/pixel
liv = cell(4,1);
liv{1} = [-0.7058 0.7058];      
liv{2} = [-1.8534 -0.4249 0.4249 1.8534];
liv{3} = [-2.9190 -1.5407 -0.7718 -0.2212 0.2212 0.7718 1.5407 2.9190];
liv{4} = [-3.9546 -2.6787 -1.8946 -1.3609 -0.9590 -0.6379 -0.3592 -0.1092 0.1092 0.3592 0.6379 0.9590 1.3609 1.8946 2.6787 3.9546];

y = liv{R};
t   = [min(x(:))-1 (y(1:N-1)+y(2:N))/2 max(x(:))+1];
xq  = 0;
for i=1:N,
   xq = xq+((t(i)<=x)&(x<t(i+1)))*y(i);
end;
xq = xq*sigma;
