% 
% Visualisation d'une image
% Decomposition sur la base de Haar
% Compression par seuillage
% Synthese
%
%  V. Perrier - Janvier 2000
%  Appelle WaveLab
%
%----------------
clc;clf;clear;
%  Visualisation de l'image
disp('Choix d image : Barton, Canaletto, Coifman, Daubechies, Fingerprint, Lincoln, Lenna, MRIScan, Phone')
     nom=input('Entrer le nom de l image :','s')
im=ReadImage(nom);
[n,J] = dyadlength(im);
% [n n]=size(im);
% jech=log(n)/log(2);
disp('Taille de l image '),n*n
figure(1)
imagesc(im)
     title('Original')
colormap('gray')

%
%--------------------
% Generation des filtres pour la decomposition

qmf = MakeONFilter('Haar',1) ;

%
%--------------------
% Calcul de la transformee en ondelettes et visualisation
L=J-4; % echelle la plus grossiere de la decomposition
wc = FWT2_PO(im,L,qmf);
disp('Affichage des coefficients')
     figure(2)
 clf
image(wc)
colormap('gray')

     disp('Appuyer sur une touche pour continuer')
pause
clf;


% 
%  Compression des coefficients
%  Seuillage "dur"
%
np=n*n;
rel=max(max(abs(wc)));

disp('COMPRESSION DE L IMAGE PAR SEUILLAGE DES COEFFICIENTS')
eps=input('entrer un seuil (relatif, compris entre 0 et 1) :')
eps=eps*rel;


for i=1:n;
for ii=1:n;
  if abs(wc(i,ii))<eps
    wc(i,ii)=0;
    np=np-1;
  end
end
end
%

% 
%  Synthese - algo de recomposition
%
disp('synthese apres seuillage - nbre de coef retenus')
disp( np) 
     yy = IWT2_PO(wc,L,qmf);
clf
  imagesc(yy);
title('Image filtree')
disp('Affichage de l image filtree')