I=imread('t5.bmp');
n1=size(I,1);
n2=size(I,2);
XX=reshape(I,[n1*n2,3]);
imshow(I);
XX=im2double(XX);
covx=cov(XX);
[pc,variances,explained]=pcacov(covx);
myKL=XX*pc;
pic1=myKL(:,1);
pic11=reshape(pic1,[n1,n2]);
figure;imshow(pic11,[]);
imwrite(pic11,'1.jpg');
评论1