%FACE RECOGNITION SCRIPT
%First takes the last selected image
%Calls the images from database
%Applies PCA to images
%Compares the last selected image with database images and says if it is
%found at the database or not.
clc
close all
if exist('img')
img2=double(img(:));
if(exist('fdata.dat')==2)
try
load('fdata.dat','-mat');
mtr=zeros(size(data{1,1},1),fnumber);
for ii=1:fnumber
mtr(:,ii)=double(data{ii,1}); %Calling the images from database
mtr2=double(mtr)/255;
avr = mean(mtr2')'; %taking the average of the training set
for i=1:fnumber
mtr2(:,i) = mtr2(:,i) - avr;
end
end
Lmat = mtr2'*mtr2 ; %Obtaining the L matrix
[V,D] = eig(Lmat); %eigen values(D) and eigenvectors(V) of the L matrix
%obtaining eigenfaces
V = mtr2*V*(abs(D))^-0.5 ; %This line is taken from the software karhunenloeve.m
%Available at
%http://www.mathworks.com/matlabcentral/fileexchange/loadFile.do?objectId=6995&objectType=file
%Author:Alex Chirokov
f1=V'*(img2-avr) ; % feature vector of the new image
fdata=zeros(max_class,max_class); %empty matrice
for ii=1:fnumber
imdata=double(data{ii,1});
classdata=data{ii,2};
cor=V'*(imdata-avr); % weight vector of the stored faces
fdata(:,classdata)=fdata(:,classdata)+cor;
end
dist=zeros(max_class,1);
for ii=1:(max_class)
dist(ii)=norm(f1-fdata(:,ii));
end
[minf,pminf]=min(dist);
k=minf/1000
hll = waitbar(0,'CHECKING THE DATABASE...');
for i=1:2450,
waitbar(i/100)
end
close(hll)
if k
disp(strcat('The corresponding ID number of the face is --> ',num2str(pminf)));
pause
ldinfo
return
end
if k>=6.1
disp('The corresponding face does not exist at your database')
disp('Note: Save more pictures to the group for the best recognition')
disp(strcat('The similar of this face has an ID, number--> ',num2str(pminf)));
pause
ldinfo
return
end
catch
clc
disp('IMAGE PROCESSING IS NOT AVAILABLE')
disp('*********************************')
disp('Possible Reasons: ')
disp(' ')
disp('1- Database is empty')
disp('2- The size of the selected image is not suitable for processing')
disp('3- The color or format is not matching with database')
disp(' ')
disp('Press any key to continue')
pause
clear all
bdrfacerec
end
else
clc
disp('IMAGE PROCESSING IS NOT AVAILABLE')
disp('*********************************')
disp('Possible Reasons: ')
disp(' ')
disp('1- Database is empty')
disp('2- The size of the selected image is not suitable for processing')
disp('3- The color or format is not matching with database')
disp(' ')
disp('Press any key to continue')
pause
clear all
bdrfacerec
end
end
