orl人脸数据库在IBP特征提取下的识别matlab代码

clear; picname=dir('D:\Program Files\work\orl'); %put 400 images into picname(402*1 structure) % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % zhucenum = input('please input the Registered number in 10:'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% zhucenum=2; M=3 N=3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % k=10;t=1; % for i=1:40*zhucenum % for j=1:zhucenum % zhucename(t)=picname(j+k*(i-1)+2).name; %%put 40 images of 10 different people in zhucename(1*40 structure) % t=t+1; % end % end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% piccell=struct2cell(picname); %make it convenient to stract picture's name(cell) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% zhucename=cell(1,40*zhucenum); ceshiname=cell(1,400-40*zhucenum); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% k=10,t=1; for i=1:40 for j=1:zhucenum j+k*(i-1)+2 zhucename{1,t}=piccell{1,j+k*(i-1)+2}; t=t+1; end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% k=10,t=1; for i=1:40 for j=1+zhucenum+(i-1)*10:i*10 j ceshiname{1,t}=piccell{1,j+2}; t=t+1; end end % % k=10;t=1; % for i=1:41 % for j=1:2 % zhucename(t)=picname(j+k*(i-1)+2); % t=t+1; % end % end %%%%%%%%%%%%%%register one detective nine % k=10; % for i=1:40 % zhucename(i)=picname(1+k*(i-1)+2); %%put 40 images of 10 different people in zhucename(1*40 structure) % end %%%%%%%%%%%%% path='D:\Program Files\work\orl\'; pic=cell(1,400); for i=1:400 picpath=strcat(path,picname(i+2).name); pic{1,i}=imread(picpath); % read into an image pic{1,i}=rgb2gray(pic{1,i}); %reduce the dimension end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % prompt = {'Enter the number of lines:','Enter the number of columes:'}; % dlg_title = 'Input the size of the window'; % num_lines = 1; % answer = inputdlg(prompt,dlg_title,num_lines); % save test.mat answer;%%test是随便起的，只是用来说明clear allload test.matanswer(1) % M=str2num(cell2mat(answer(1))); % N=str2num(cell2mat(answer(2))); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% pich=cell(1,400); for i=1:400 pich{i}=lbp2h(pic{i},M,N);%% change pic{1,i}into k1*k2 segments end % save 400pich pich % save zhucenum %%%%%%%%%%%%%%%%%%%image discrimination%%%%%%%% % load('D:\Program Files\work\lbp\400pich.mat'); correct=0; %calculating right number % zhuce=cell(1,40*zhucenum); zhuceh=cell(1,40*zhucenum); ceshih=cell(1,400-40*zhucenum); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% k=10;t=1; for i=1:40 for j=1:zhucenum zhuceh{t}=pich{1,j+k*(i-1)};%%define zhuce's histogram t=t+1; end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% k=10,t=1; for i=1:40 for j=(1+zhucenum+(i-1)*10):i*10 j ceshih{t}=pich{j}; t=t+1; end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% for i=1:400-40*zhucenum i % if (mod(i,10)>zhucenum) for j=1:40*zhucenum he(j)=sum(min(ceshih{1,i},zhuceh{1,j})); %output the simvalue of two matrix end [matrix,num]=max(he); %%%%%find the index of the biggest value in sim[]%%% a=zhucename{num}(5:7); b=ceshiname{i}(5:7); if strcmp(a,b) correct=correct+1 end end correctrate=correct/(400-40*zhucenum)