基于MATLAB的人臉識別源程序

精選優(yōu)質(zhì)文檔-----傾情為你奉上精選優(yōu)質(zhì)文檔-----傾情為你奉上專心---專注---專業(yè)專心---專注---專業(yè)精選優(yōu)質(zhì)文檔-----傾情為你奉上專心---專注---專業(yè)1.色彩空間轉(zhuǎn)換function[r,g]=rgb_RGB(Ori_Face)R=Ori_Face(:,:,1);G=Ori_Face(:,:,2);B=Ori_Face(:,:,3);R1=im2double(R);

%將uint8型轉(zhuǎn)換成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);2.均值和協(xié)方差t1=imread('D:\matlab\皮膚庫\1.jpg');[r1,g1]=rgb_RGB(t1);t2=imread('D:\matlab\皮膚庫\2.jpg');[r2,g2]=rgb_RGB(t2);t3=imread('D:\matlab\皮膚庫\3.jpg');[r3,g3]=rgb_RGB(t3);t4=imread('D:\matlab\皮膚庫\4.jpg');[r4,g4]=rgb_RGB(t4);t5=imread('D:\matlab\皮膚庫\5.jpg');[r5,g5]=rgb_RGB(t5);t6=imread('D:\matlab\皮膚庫\6.jpg');[r6,g6]=rgb_RGB(t6);t7=imread('D:\matlab\皮膚庫\7.jpg');[r7,g7]=rgb_RGB(t7);t8=imread('D:\matlab\皮膚庫\8.jpg');[r8,g8]=rgb_RGB(t8);t9=imread('D:\matlab\皮膚庫\9.jpg');[r9,g9]=rgb_RGB(t9);t10=imread('D:\matlab\皮膚庫\10.jpg');[r10,g10]=rgb_RGB(t10);t11=imread('D:\matlab\皮膚庫\11.jpg');[r11,g11]=rgb_RGB(t11);t12=imread('D:\matlab\皮膚庫\12.jpg');[r12,g12]=rgb_RGB(t12);t13=imread('D:\matlab\皮膚庫\13.jpg');[r13,g13]=rgb_RGB(t13);t14=imread('D:\matlab\皮膚庫\14.jpg');[r14,g14]=rgb_RGB(t14);t15=imread('D:\matlab\皮膚庫\15.jpg');[r15,g15]=rgb_RGB(t15);t16=imread('D:\matlab\皮膚庫\16.jpg');[r16,g16]=rgb_RGB(t16);t17=imread('D:\matlab\皮膚庫\17.jpg');[r17,g17]=rgb_RGB(t17);t18=imread('D:\matlab\皮膚庫\18.jpg');[r18,g18]=rgb_RGB(t18);t19=imread('D:\matlab\皮膚庫\19.jpg');[r19,g19]=rgb_RGB(t19);t20=imread('D:\matlab\皮膚庫\20.jpg');[r20,g20]=rgb_RGB(t20);t21=imread('D:\matlab\皮膚庫\21.jpg');[r21,g21]=rgb_RGB(t21);t22=imread('D:\matlab\皮膚庫\22.jpg');[r22,g22]=rgb_RGB(t22);t23=imread('D:\matlab\皮膚庫\23.jpg');[r23,g23]=rgb_RGB(t23);t24=imread('D:\matlab\皮膚庫\24.jpg');[r24,g24]=rgb_RGB(t24);t25=imread('D:\matlab\皮膚庫\25.jpg');[r25,g25]=rgb_RGB(t25);t26=imread('D:\matlab\皮膚庫\26.jpg');[r26,g26]=rgb_RGB(t26);t27=imread('D:\matlab\皮膚庫\27.jpg');[r27,g27]=rgb_RGB(t27);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])3.求質(zhì)心function[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);end4.求偏轉(zhuǎn)角度function[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);

%從幅度轉(zhuǎn)換到角度5.找區(qū)域邊界function[left,right,up,down]=bianjie(A)[mn]=size(A);left=-1;right=-1;up=-1;down=-1;forj=1:n,

fori=1:m,

if(A(i,j)~=0)

left=j;

break;

end;

end;

if(left~=-1)break;

end;end;forj=n:-1:1,fori=1:m,if(A(i,j)~=0)

right=j;

break;end;end;if(right~=-1)break;end;end;fori=1:m,forj=1:n,

if(A(i,j)~=0)

up=i;

break;

end;end;

if(up~=-1)

break;

end;end;fori=m:-1:1,

forj=1:n,

if(A(i,j)~=0)

down=i;

break;

end;

end;

if(down~=-1)

break;

end;end;6.求起始坐標(biāo)functionnewcoord=checklimit(coord,maxval)newcoord=coord;if(newcoord<1)

newcoord=1;end;if(newcoord>maxval)

newcoord=maxval;end;7.模板匹配function[ccorr,mfit,RectCoord]=mobanpipei(mult,frontalmodel,ly,wx,cx,cy,angle)frontalmodel=rgb2gray(frontalmodel);model_rot=imresize(frontalmodel,[lywx],'bilinear');%調(diào)整模板大小model_rot=imrotate(model_rot,angle,'bilinear');

%旋轉(zhuǎn)模板[l,r,u,d]=bianjie(model_rot);

%求邊界坐標(biāo)bwmodel_rot=imcrop(model_rot,[lu(r-l)(d-u)]);

%選擇模板人臉區(qū)域[modx,mody]=center(bwmodel_rot);

%求質(zhì)心[morig,norig]=size(bwmodel_rot);

%產(chǎn)生一個覆蓋了人臉模板的灰度圖像mfit=zeros(size(mult));mfitbw=zeros(size(mult));[limy,limx]=size(mfit);%計(jì)算原圖像中人臉模板的坐標(biāo)startx=cx-modx;starty=cy-mody;endx=startx+norig-1;endy=starty+morig-1;startx=checklimit(startx,limx);starty=checklimit(starty,limy);endx=checklimit(endx,limx);endy=checklimit(endy,limy);fori=starty:endy,forj=startx:endx,mfit(i,j)=model_rot(i-starty+1,j-startx+1);end;end;ccorr=corr2(mfit,mult)

%計(jì)算相關(guān)度[l,r,u,d]=bianjie(bwmodel_rot);sx=startx+l;sy=starty+u;RectCoord=[sxsy(r-1)(d-u)];%產(chǎn)生矩形坐標(biāo)8.主程序clear;[fname,pname]=uigetfile({'*.jpg';'*.bmp';'*.tif';'*.gif'},'Pleasechooseacolorpicture...');

%返回打開的圖片名與圖片路徑名[u,v]=size(fname);y=fname(v);

%圖片格式代表值switchy

case0

errordlg('YouShouldLoadImageFileFirst...','Warning...');

case{'g';'G';'p';'P';'f';'F'};

%圖片格式若是JPG/jpg、BMP/bmp、TIF/tif或者GIF/gif，才打開

I=cat(2,pname,fname);

Ori_Face=imread(I);

subplot(2,3,1),imshow(Ori_Face);otherwise

errordlg('YouShouldLoadImageFileFirst...','Warning...');endR=Ori_Face(:,:,1);G=Ori_Face(:,:,2);B=Ori_Face(:,:,3);R1=im2double(R);

%將uint8型轉(zhuǎn)換成double型處理G1=im2double(G);B1=im2double(B);RGB=R1+G1+B1;m=[0.4144,0.3174];

%均值n=[0.0031,-0.0004;-0.0004,0.0003];

%方差row=size(Ori_Face,1);

%行像素?cái)?shù)column=size(Ori_Face,2);

%列像素?cái)?shù)fori=1:row

forj=1:column

ifRGB(i,j)==0

rr(i,j)=0;gg(i,j)=0;

else

rr(i,j)=R1(i,j)/RGB(i,j);

%rgb歸一化

gg(i,j)=G1(i,j)/RGB(i,j);

x=[rr(i,j),gg(i,j)];

p(i,j)=exp((-0.5)*(x-m)*inv(n)*(x-m)');

%皮膚概率服從高斯分布

end

endendsubplot(2,3,2);imshow(p);

%顯示皮膚灰度圖像low_pass=1/9*ones(3);image_low=filter2(low_pass,p);

%低通濾波去噪聲subplot(2,3,3);imshow(image_low);%自適應(yīng)閥值程序previousSkin2=zeros(i,j);changelist=[];forthreshold=0.55:-0.1:0.05two_value=zeros(i,j);two_value(find(image_low>threshold))=1;change=sum(sum(two_value-previousSkin2));changelist=[changelistchange];previousSkin2=two_value;end[C,I]=min(changelist);optimalThreshold=(7-I)*0.1two_value=zeros(i,j);two_value(find(image_low>optimalThreshold))=1;

%二值化subplot(2,3,4);imshow(two_value);

%顯示二值圖像frontalmodel=imread('E:\我的照片\人臉模板.jpg');

%讀入人臉模板照片F(xiàn)aceCoord=[];imsourcegray=rgb2gray(Ori_Face);

%將原照片轉(zhuǎn)換為灰度圖像[L,N]=bwlabel(two_value,8);

%標(biāo)注二值圖像中連接的部分,L為數(shù)據(jù)矩陣，N為顆粒的個數(shù)fori=1:N,[x,y]=find(bwlabel(two_value)==i);

%尋找矩陣中標(biāo)號為i的行和列的下標(biāo)bwsegment=bwselect(two_value,y,x,8);

%選擇出第i個顆粒numholes=1-bweuler(bwsegment,4);

%計(jì)算此區(qū)域的空洞數(shù)if(numholes>=1)

%若此區(qū)域至少包含一個洞，則將其選出進(jìn)行下一步運(yùn)算RectCoord=-1;[mn]=size(bwsegment);[cx,cy]=center(bwsegment);

%求此區(qū)域的質(zhì)心bwnohole=bwfill(bwsegment,'holes');

%將洞封住（將灰度值賦為1）justface=uint8(double(bwnohole).*double(imsourcegray));

%只在原照片的灰度圖像中保留該候選區(qū)域angle=orient(bwsegment,cx,cy);

%求此區(qū)域的偏轉(zhuǎn)角度

bw=imrotate(bwsegment,angle,'bilinear');bw=bwfill(bw,'holes');[l,r,u,d]=bianjie(bw);wx=(r-l+1);

%寬度ly=(d-u+1);

%高度wratio=ly/wx

%高寬比

if((0.8<=wratio)&(wratio<=2))

%如果目標(biāo)區(qū)域的高度/寬度比例大于0.8且小于2.0，則將其選出進(jìn)行下一步運(yùn)算

S=ly*wx;

%計(jì)算包含此區(qū)域矩形的面積

A=bwarea(bwsegment);

%計(jì)算此區(qū)域面積

if(A/S>0.35)

[ccorr,mfit,RectCoord]=mobanpipei(justface,frontalmodel,ly,wx,cx,cy,angle);

end

if(ccorr>=0.6)

mfitbw=(mfit>=1);

invbw=xor(mfitbw,ones(size(mfitbw)));

source_with_hole=uint8(double(invbw).*double(imsourcegray));

final_image=uint8(double(source_with_hole)+double(mfit));

subplot(2,3,5);imshow(final_image);

%顯示覆蓋了模板臉的灰度圖像

imsourcegray=final_image;

subplot(2,3,6);imshow(Ori_Face);

%顯示檢測效果圖

end;

if(RectCoord~=-1)

FaceCoord=[FaceCoord;RectCoord];

end

endendend%在認(rèn)為