smildamain5.m

clc;
close all;
clear all;
TrainDatabasePath = uigetdir('C:\Users\Fathima\Desktop\thesis folder\orl', 'Select training database path');
M=40;
N=9;
%P=2;
M1=M*N;
gama=[];
num_lines= 1;
def = {'1','10'};
prompt={'Enter test image name (a number between 1 to 40):','enter a no:between 1 to 10'};
dlg_title = 'Input of PCA-Based Face Recognition System';
TestDatabasePath = uigetdir('C:\Users\Fathima\Desktop\thesis folder\orl', 'Select test database path');
answer = inputdlg(prompt,dlg_title,num_lines,def);
x = str2num(answer{1}); 
y = str2num(answer{2});
TestImage= strcat('\S',sprintf('%d',x),'\',int2str(y),'.pgm');
TestImage = strcat(TestDatabasePath,TestImage);
for i=1:M
    for j=1:N
%         if (i==x)&&(j==y)
%             continue;
%         end
        str = strcat('\S',sprintf('%d',i),'\',int2str(j),'.pgm');
        str = strcat(TrainDatabasePath,str);
         img = imread(str);
   
    temp=transpose(img);
    temp=reshape(temp,numel(temp),1);
   gama=[gama temp];
    end
end
mean_All= mean(gama,2);
A = [];  
for i = 1 : M1
    t = double(gama(:,i)) - mean_All;
    A = [A t];
end
L=A'*A;
[V D] = eig(L);
L_eig_vec = [];

for i=360:-1:320;
   
        L_eig_vec = [L_eig_vec V(:,i)];
   
end
Eigenfaces = A * L_eig_vec; 
 
 ProjectedImages_PCA= [];
 for i = 1 :M1
    temp = Eigenfaces'*A(:,i); % Projection of centered images into facespace
    ProjectedImages_PCA = [ProjectedImages_PCA temp]; 
 end
 P=M*N;
 Class_population=N;
 Class_number=M;
 m_PCA = mean(ProjectedImages_PCA,2); % Total mean in eigenspace
m = zeros(41,Class_number); 
Sw = zeros(41,41); % Initialization os Within Scatter Matrix
Sb = zeros(41,41); % Initialization of Between Scatter Matrix

 for i = 1 : Class_number
    m(:,i) = mean( ( ProjectedImages_PCA(:,((i-1)*Class_population+1):i*Class_population) ), 2 )';    
    
    S  = zeros(41,41); 
    for j = ( (i-1)*Class_population+1 ) : ( i*Class_population )
        S = S + (ProjectedImages_PCA(:,j)-m(:,i))*(ProjectedImages_PCA(:,j)-m(:,i))';
    end
    
    Sw = Sw + S; % Within Scatter Matrix
Sb = Sb + Class_number*((m(:,i)-m_PCA) * (m(:,i)-m_PCA)'); % Between Scatter Matrix
end
 [J_eig_vec, J_eig_val] = eig(Sb,Sw); % Cost function J = inv(Sw) * Sb
J_eig_vec = fliplr(J_eig_vec);

%%%%%%%%%%%%%%%%%%%%%%%% Eliminating zero eigens and sorting in descend order
for i = 1 : Class_number-1 
    V_Fisher(:,i) = J_eig_vec(:,i); % Largest (C-1) eigen vectors of matrix J
end
 
for i = 1 : Class_number*Class_population
    ProjectedImages_Fisher(:,i) = V_Fisher' * ProjectedImages_PCA(:,i);
end
 InputImage= imread(TestImage);
 figure;
 imshow(InputImage);title('Input Test Image');
 temp = InputImage(:,:,1);
[irow icol] = size(temp);
InImage = reshape(temp',irow*icol,1);
Difference = double(InImage)- mean_All; % Centered test image
ProjectedTestImage = V_Fisher' * Eigenfaces' * Difference; % Test image feature vectorED = [];
ED=[];
for i=1:M1
    
q = ProjectedImages_Fisher(:,i);
d= ProjectedTestImage - q;
d1=sum(d.*d);
ED=[ED d1];
end
[Euc_dist_min , Recognized_index] = min(ED);
n = Recognized_index;
for L=1:M
    if n<=9*L
        str =  strcat('Matched person is :  ',int2str(L));
        disp(str);
        break
    end
end

b= mod(n,9);
if b==0
    b=9;
end

outputname= strcat('\S',sprintf('%d',L),'\',int2str(b),'.pgm');
disp(outputname);
outputname = strcat(TestDatabasePath,outputname);
outputname=imread(outputname);
figure,imshow(outputname);
title('Matched Image');
INDEX=[];
%  for i=1:M1
% 
%      if ED(i)<=0.5
%          INDEX=[INDEX i];
%      end
%  end
 u=[];
 t=[];
 for i=1:10
[Euc_dist_min , Recognized_index] = min(ED);
n = Recognized_index;
u=[u n];
t=[t Euc_dist_min];
ED(n)=100
end
 INDEX=u;
 T=size(INDEX,2)
 if T==0
  display('Not authenticated'); 
  break
 end
 
 FT=[];
for i=1:10
fvi=u(i);
  FM=FINGERFIN5(fvi);
  FT=[FT FM];
end
 k=smipcamain5(INDEX,InputImage,gama);
 %[tem in]=min(k);
 f=k;
 u1=[];
 t1=[];
 for i=1:5
[Euc_dist_min , Recognized_index] = min(k);
n = Recognized_index;
u1=[u1 FT(n)];
t1=[t1 Euc_dist_min];
k(n)=100
end
%   fvi=INDEX(in);
%   FM=FINGERFIN5(fvi);
%   disp(FM);
%   d=mainedt5()
%   g=size(d);
%   for i=1:g
%       if d(g)==FM
%        str =  strcat('Matched person is :  ',int2str(FM));
%        disp(str);
%       display('Authenticated');
%       end
%   end