基于MATLAB的人脸识别源程序.doc

--------------------------校验:_____________-----------------------日期:[r,g]=rgb_RGB(Ori_Face)R=Ori_Face(:,:,1);G=Ori_Face(:,:,2);B=Ori_Face(:,:,3);R1=im2double(R);       %将uint8型转换成double型G1=im2double(G);B1=im2double(B);RGB=R1+G1+B1;row=size(Ori_Face,1);     %行像素column=size(Ori_Face,2);  %列像素fori=1:rowforj=1:column     rr(i,j)=R1(i,j)/RGB(i,j);     gg(i,j)=G1(i,j)/RGB(i,j);endendrrr=mean(rr);r=mean(rrr);ggg=mean(gg);g=mean(ggg);=imread('D:\matlab\皮肤库t2=imread('D:\matlab\皮肤库]=rgb_RGB(t2);t3=imread('D:\matlab\皮肤库t4=imread('D:\matlab\皮肤库t5=imread('D:\matlab\皮肤库t6=imread('D:\matlab\皮肤库t7=imread('D:\matlab\皮肤库\t8=imread('D:\matlab\皮肤库t9=imread('D:\matlab\皮肤库t10=imread('D:\matlab\皮肤库t11=imread('D:\matlab\皮肤库t12=imread('D:\matlab\皮肤库t13=imread('D:\matlab\皮肤库t14=imread('D:\matlab\皮肤库t15=imread('D:\matlab\皮肤库t16=imread('D:\matlab\皮肤库t17=imread('D:\matlab\皮肤库t18=imread('D:\matlab\皮肤库t19=imread('D:\matlab\皮肤库t20=imread('D:\matlab\皮肤库20]=rgb_RGB(t20);t21=imread('D:\matlab\皮肤库t22=imread('D:\matlab\皮肤库t23=imread('D:\matlab\皮肤库t24=imread('D:\matlab\皮肤库t25=imread('D:\matlab\皮肤库t26=imread('D:\matlab\皮肤库t27=imread('D:\matlab\皮肤库r=cat(1,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,r12,r13,r14,r15,r16,r17,r18,r19,r20,r21,r22,r23,r24,r25,r26,r27);g=cat(1,g1,g2,g3,g4,g5,g6,g7,g8,g9,g10,g11,g12,g13,g14,g15,g16,g17,g18,g19,g20,g21,g22,g23,g24,g25,g26,g27);m=mean([r,g])n=cov([r,g])[xmean,ymean]=center(bw)bw=bwfill(bw,'holes');area=bwarea(bw); [mn]=size(bw);bw=double(bw);xmean=0;ymean=0;fori=1:m,forj=1:n,xmean=xmean+j*bw(i,j);ymean=ymean+i*bw(i,j);end;end;if(area==0)   xmean=0;   ymean=0;else   xmean=xmean/area;   ymean=ymean/area;   xmean=round(xmean); ymean=round(ymean);[theta]=orient(bw,xmean,ymean)[mn]=size(bw);bw=double(bw);a=0;b=0;c=0;fori=1:m,forj=1:n,a=a+(j-xmean)^2*bw(i,j);b=b+(j-xmean)*(i-ymean)*bw(i,j);c=c+(i-ymean)^2*bw(i,j);end;end;b=2*b;theta=atan(b/(a-c))/2;theta=theta*(180/pi);  %