PCA.zip_pca_闪电

tif1=double(imread('C:\image\tm_Kennedy\tm1.tif')); tif2=double(imread('C:\image\tm_Kennedy\tm2.tif')); tif3=double(imread('C:\image\tm_Kennedy\tm3.tif')); [z1,s1]=size(tif1); [z2,s2]=size(tif2); [z3,s3]=size(tif3); if(z1~=z2)|(z2~=z3)|(s1~=s2)|(s2~=s3) error('input images are not of same size'); end; s1=[tif1(:) tif2(:) tif3(:)]; s=cov(s1); [v,d]=eig(s); w=v(:,1)'; y=w*s1'; y=reshape(y,512,512); figure,imshow(y) patterns=imread('c:\image\LENA256.BMP'); patterns=double(patterns); [r,c]=size(patterns); dimension=200; if(r<dimension) disp('reqired dimension is larger than the data dimension.') disp('[will use dimension num2str(r)]') dimension=r; end %calculate cov matrix and the PCA matrixes m=mean(patterns')'; s=((patterns-m*ones(1,c))*(patterns-m*ones(1,c))'); [v,d]=eig(s); w=v(:,r-dimension+1:r)'; u=s*w'*inv(w*s*w'); %calculate new patterns uw=u*w; patternsy=w*patterns; figure,imshow(patternsy) patternsx=w'*patternsy; figure,imshow(patternsx)