Utiliser SVD pour compresser une image dans MATLAB

je suis tout nouveau à MATLAB mais j'essaie de faire un peu de code de compression d'image pour les images en niveaux de gris.

Questions

Comment puis-je utiliser SVD pour couper des valeurs propres de faible valeur pour reconstruire une image compressée?

Travail/Tentatives de la mesure

Mon code pour l'instant est:

B=imread('images1.jpeg');   
B=rgb2gray(B);  
doubleB=double(B);  
%read the image and store it as matrix B, convert the image to a grayscale
photo and convert the matrix to a class 'double' for values 0-255  
[U,S,V]=svd(doubleB);

cela me permet de décomposer avec succès la matrice d'image avec des valeurs propres stockées dans la variable S.

Comment puis-je troncature S (qui est 167x301, classe double)? Disons que sur les 167 valeurs propres que je veux prendre seulement le top 100 (ou n'importe quel n vraiment), comment puis-je faire cela et reconstruire l'image compressée?

mise à Jour de code/pensées

au Lieu de mettre un tas de code dans la section des commentaires, c'est l'actuel projet que j'ai. J'ai réussi à créer l'image compressée en changeant manuellement N, mais j'aimerais faire 2 choses supplémentaires:

1- Afficher un panneau d'images pour différentes compressions (i/e, lancer une boucle pour N = 5,10,25, etc.)

2 - calculez D'une façon ou d'une autre la différence (erreur) entre chaque image et l'original et faites le graphe.

je suis horrible avec la compréhension des boucles et de la sortie, mais c'est ce que j'ai essayé:

B=imread('images1.jpeg');  
B=rgb2gray(B);  
doubleB=im2double(B);%  
%read the image and store it as matrix B, convert the image to a grayscale  
%photo and convert the image to a class 'double'  
[U,S,V]=svd(doubleB);   
C=S;  
for N=[5,10,25,50,100]  
C(N+1:end,:)=0;  
C(:,N+1:end)=0;  
D=U*C*V';  
%Use singular value decomposition on the image doubleB, create a new matrix  
%C (for Compression diagonal) and zero out all entries above N, (which in  
%this case is 100). Then construct a new image, D, by using the new  
%diagonal matrix C.  
imshow(D);  
error=C-D;  
end

évidemment il y a des erreurs parce que je ne reçois pas plusieurs images ou ne sais pas comment "graphe" la matrice d'erreur