//Lecture et décomposition de l'image

N=500*640;
a= read('im34.txt',-1,640);
[U,S,V]=svd(a);
ap = uint8(255*a);
figure();
ShowImage(ap,'Image originale');

//Suppression des valeurs singulières

valS=diag(S);
seuil=length(valS)*0.5;
//
for(i=seuil:length(valS))
    S(i,i)=0;
    end

//Reconstruction de la matrice et affichage de l'image

a = U*S*V';
ap=uint8(255*a);

figure();
ShowImage(ap,'Image Reconstruite'); 

