close all
[A,map]=imread(’lena.gif’);
B=im2double(A,’indexed’);
imshow(B,map)
[u,s,v]=svd(B);
C=zeros(size(B));
for j=1:k
C=C+s(j,j)*u(:,j)*v(:,j).’;
end
C=floor(C);
imshow(C,map)
k=find(C<1);
C(k)=1;
set(gcf,’Unit’,’inches’,’Paperposition’,[0,0,2,1])
print -djpeg ’lenak.jpg'
http://online.redwoods.cc.ca.us/instruct/darnold/laproj/fall2001/adamdave/textwriteup.pdf
http://online.redwoods.cc.ca.us/instruct/darnold/laproj/fall2001/adamdave/textwriteup.pdf