2dpca.m

来自「经典的人脸识别算法：2Dpca的实现」· M 代码 · 共 57 行

num=104;
p=4;
q=4;
m=192/p;n=128/q;
G=zeros([m m]);
B=zeros([m n]);
temp=zeros([m n]);
T = dctmtx(8); 
mask = [1 1 1 1 0 0 0 0
        1 1 1 0 0 0 0 0
        1 1 0 0 0 0 0 0 
        1 0 0 0 0 0 0 0
        0 0 0 0 0 0 0 0
        0 0 0 0 0 0 0 0
        0 0 0 0 0 0 0 0
        0 0 0 0 0 0 0 0 ];
for i=1:num*2
    P=imread(strcat(num2str(i),'.jpg'));
    P=im2double(P);
    P = blkproc(P,[8 8], 'P1*x*P2',T,T');
    P = blkproc(P,[8 8],'P1.*x',mask) ;
    P = blkproc(P,[8 8] ,'P1*x*P2',T',T) ;
    I(:,:,i)=P;
end
 for i=1:num*2
     for  k=1:p
         for l=1:q
             temp=I(:,:,i);
             B=B+temp((k-1)*m+1:k*m,(l-1)*n+1:l*n);
         end
     end
 end
 B=B/(num*2*p*q);
 for i=1:num*2
     for k=1:p
         for l=1:q
             temp=I(:,:,i);
             G=G+(temp((k-1)*m+1:k*m,(l-1)*n+1:l*n)-B)*(temp((k-1)*m+1:k*m,(l-1)*n+1:l*n)-B)';
         end
     end
 end
G=G/(50*p*q);
[V,D]=eig(G);
X=V(:,48);
for i=1:2:2*num-1
    I(:,:,i)=(I(:,:,i)+I(:,:,i+1))/2;
end
for i=1:num
   W(:,:,i)=blkproc(I(:,:,2*i-1),[48 32] ,'(x'')*P1',X) ;
   fid=fopen(strcat('d:/',num2str(i),'.txt'),'w');
   fprintf(fid,'%12.8f ',W(:,:,i));
   status=fclose(fid);
end
fid=fopen('X.txt','w');
fprintf(fid,'%12.8f ',X);
status=fclose(fid);