function [y,a,b] = quant(x,q,s)
%        [y,a,b] = quant(x,q,s)
% y = sortie quantifiee
% a,b: affinite y = round(a*x+b)
% pas de quantification = q/log(variance)
% s = parametre optionnel,'v' pour voir les
% sous-images, 'c' pour les compresser
%------------------------------------------------

[nl,nc] = size(x);
xl = reshape(x,1,nl*nc);
mn = mean(xl);
pas = q/max(log(cov(xl)),log(1.1));
a = 1/pas;
b = -mn/pas;
y = round(a*x + b);
%y = a*x + b; % pour test PR

if (s == 'v')
  y = (y-b)/a;
  y = 5*y+100;
end
