基于MATLAB的人脸识别源程序

1、色彩空间转换1.function r,g=rgb_RGB(Ori_Face)R=Ori_Face(:,:,1);G=Ori_Face(:,:,2);B=Ori_Face(:,:,3);将型转换成型double uint8 %R1=im2double(R);G1=im2double(G);B1=im2double(B);RGB=R1+G1+B1;行像素 row=size(Ori_Face,1);%列像素 %column=size(Ori_Face,2);for i=1:row for j=1:columnrr(i,j)=R1(i,j)/RGB(i,j);gg(i,j)=G1(i,j)/RGB(i

2、,j);endendrrr=mean(rr);r=mean(rrr);ggg=mean(gg);g=mean(ggg);均值和协方差2.学资学习网提供考研资料皮肤库学资学习网提供考研资料皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库皮肤库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,g

3、6,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.function xmean, ymean = center(bw) bw=bwfill(bw,'holes');area= bwarea(bw);m n =size(bw);bw=double(bw);xmean=0; ymean = 0;for i=1:m,for j=1:n,xmean= xmean+ j*bw(i,j);i*bw(i,j); + yme

4、an= ymean end;end;if(area=0) xmean=0;ymean=0;else xmean= xmean/area;ymean= ymean/area;xmean= round(xmean);ymean= round(ymean);end求偏转角度4.function theta = orient(bw,xmean,ymean)m n =size(bw);bw=double(bw);a= 0;b= 0;c= 0;for i=1:m,for j=1:n,a = a + (j - xmean)八2* bw(i,j);bw(i,j); * ymean) - (i * xmean)

5、 - (j + b = b c = c + (i - ymean)八2* bw(i,j);end;end;b=2* b;theta= atan(b/(a-c)/2;从幅度转换到角度% = theta*(180/pi);theta找区域边界5.function left, right, up, down = bianjie(A)m n = size(A);left = -1;right = -1;up= -1;down = -1;for j=1:n, for i=1:m,if (A(i,j) =0)left = j;break;end;end;if (left = -1) break;end;e

6、nd;for j=n:-1:1,学资学习网提供考研资料for i=1:m,学资学习网提供考研资料if (A(i,j) =0)right =j;break;end;end;if (right = -1) break;end;end;for i=1:m,for j=1:n,if (A(i,j) =0) up=i;break;end;end;if (up = -1)break;end;end;for i=m:-1:1,for j=1:n,if (A(i,j) =0)down = i;break;end;end;if (down = -1) break;end;end;求起始坐标6.function

7、newcoord = checklimit(coord,maxval) coord; newcoord =(newcoord<1) if newcoord=1;end;(newcoord>maxval)ifnewcoord=maxval;end;模板匹配7.mobanpipei(mult,=RectCoordmfit,ccorr,functionfrontalmodel,ly,wx,cx, cy, angle)frontalmodel=rgb2gray(frontalmodel);调整模板大小%wx,'bilinear'); model_rot =imresize(

8、frontalmodel,ly旋转模板% imrotate(model_rot,angle,'bilinear');model_rot =求边界坐标% bianjie(model_rot);l,r,u,d =选择模板人脸bwmodel_rot=imcrop(model_rot,l u (r-l) (d-u);%区域求质心 =center(bwmodel_rot);%modx,modymorig, norig = size(bwmodel_rot);产生一个覆盖了人脸模板的灰度图像%mfit = zeros(size(mult);mfitbw = zeros(size(mult)

9、;limy, limx = size(mfit);计算原图像中人脸模板的坐标%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);checklimit(endy,limy); = endyfor i=starty:endy,for j=startx:endx,mfit(i,j) = model_r

10、ot(i-starty+1,j-startx+1);end;end;计算相关度corr2(mfit,mult)%ccorr =l,r,u,d = bianjie(bwmodel_rot);sx= startx+l;sy= starty+u;产生矩形坐标% (d-u); sx sy(r-1)RectCoord =主程序 8.clear;fname,pname=uigetfile('*.jpg''*.bmp''*.tif''*.gif','Pleasechoosea返回打开的图片名与图片路径名picture.');%c

11、oloru,v=size(fname);y=fname(v);%switch ycase0errordlg('You Should Load ImageFile First.','Warning.');图片格式若是%case'g''G''p''P''f''F'、 、或者，才打开GIF/gif JPG/jpgBMP/bmpTIF/tifI=cat(2,pname,fname);Ori_Face=imread(I);subplot(2,3,1),imshow(Ori_Fa

12、ce);otherwiseerrordlg('You Should Load ImageFile First.','Warning.');endR=Ori_Face(:,:,1);G=Ori_Face(:,:,2);B=Ori_Face(:,:,3);将型转换成型处理double uint8 %R1=im2double(R);G1=im2double(G);B1=im2double(B);RGB=R1+G1+B1;均值 0.4144,0.3174;%m=方差 n=0.0031,-0.0004;-0.0004,0.0003;%行像素数%row=size(Ori_F

13、ace,1);列像素数column=size(Ori_Face,2);% for i=1:rowfor j=1:columnif RGB(i,j)=0rr(i,j)=0;gg(i,j)=0;else归一化 rgbrr(i,j)=R1(i,j)/RGB(i,j); %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)');%布endend end显示皮肤灰度图像subplot(2,3,2);imshow(p);% low_pass=1/9*ones(3);低通滤波去

14、噪声%image_low=filter2(low_pass, p);subplot(2,3,3);imshow(image_low);previousSkin2 = zeros(i,j);changelist= ;for threshold= 0.55:-0.1:0.05 two_value = zeros(i,j);two_value(find(image_low>threshold) = 1;change= sum(sum(two_value- previousSkin2);changelist= changelist change;previousSkin2 = two_value

15、;endmin(changelist); = I C,optimalThreshold = (7-I)*0.1two_value = zeros(i,j);二值化 %1; two_value(find(image_low>optimalThreshold) =显示二值图像%subplot(2,3,4);imshow(two_value);读入人脸模板我的照片人脸模板.jpg');%照片FaceCoord=;将原照片转换为灰度imsourcegray=rgb2gray(Ori_Face);%图像标注二值图像中连接的部分为数,LL,N=bwlabel(two_value,8);%据矩

16、阵，为颗粒的个数Nfor i=1:N,寻找矩阵中标号为的i %x,y=find(bwlabel(two_value)=i);行和列的下标选择出第个颗粒% bwselect(two_value,y,x,8);ibwsegment=计算此区域的空洞数numholes= 1-bweuler(bwsegment,4);%若此区域至少包含一个洞，则将其选出% 1) if (numholes>=进行下一步运算RectCoord= -1;m n = size(bwsegment);求此区域的质心%cx,cy=center(bwsegment);将洞封住(将灰度值赋bwnohole=bwfill(bws

17、egment,'holes');%为) 1justface = uint8(double(bwnohole) .* double(imsourcegray);只在原照片的灰度图像中保留该候选区%域求此区域的偏转角度% orient(bwsegment,cx,cy);angle=bw = imrotate(bwsegment,angle,'bilinear');bw = bwfill(bw,'holes');l,r,u,d =bianjie(bw);宽度+1);% l wx = (r -高度ly = (d - u +% 1);高宽比 wratio

18、= ly/wx%if (0.8<=wratio)&(wratio<=2)如果目标区域的高度宽度比例大于且小于，则将2.0 0.8%/其选出进行下一步运算计算包含此区域矩形的面积%S=ly*wx;计算此区域面积A=bwarea(bwsegment);%if (A/S>0.35)ccorr,mfit, RectCoord = mobanpipei(justface,frontalmodel,ly,wx, cx,cy, angle);end(ccorr>=0.6) ifmfitbw=(mfit>=1);invbw = xor(mfitbw,ones(size(m

19、fitbw);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 end end end在认为是人脸的区域画矩形%numfacesx = size(FaceCoord);for i=1:numfaces, hd = rectangle('Position',FaceCoord(i,:);set(hd,'edgecolor','y');end人脸检测是人