《人脸识别MATLAB代码》由会员分享，可在线阅读，更多相关《人脸识别MATLAB代码（11页珍藏版）》请在金锄头文库上搜索。
1、1.色彩空间转换function 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);% 列像素for i=1:row for j=1:columnrr(i,j)=R1(i,j)/RGB(i,j); gg(i,j)=G1(i,j)/
2、RGB(i,j);endendrrr=mean(rr);r=mean(rrr);ggg=mean(gg);g=mean(ggg);2.均值和协方差tl=imread(D:matlab皮肤库l.jpg);rl,gl=rgb_RGB(tl); 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=rg
3、b_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);tlO=imread(D:matlab 皮肤库10.jpg);rl0,gl0=rgb_RGB(tl0);t11=imread(D:matlab 皮肤库11.jpg);r11,g11=rgb_RGB(t1
4、1);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,g1
5、7=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 皮
6、肤库 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,r
7、15,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.求质心function xmean, ymean = center(bw) bw=bwfill(bw,holes);area = bwarea(bw);m n =size(bw); bw=double(bw); xmean =0;
8、 ymean = 0; for i=1:m, for j=1:n, xmean = xmean + j*bw(i,j); ymean = ymean + i*bw(i,j); end; end;if(area=0)xmean=0; ymean=0;elsexmean = xmean/area;ymean = ymean/area; xmean = round(xmean);ymean = round(ymean);end4. 求偏转角度function theta = orient(bw,xmean,ymean) m n =size(bw);bw=double(bw);a = 0; b = 0
9、;c = 0;for i=1:m,for j=1:n,a = a + (j - xmean)人2 * bw(i,j);b = b + (j - xmean) * (i - ymean) * bw(i,j); c = c + (i - ymean)A2 * bw(i,j);end;end;b = 2 * b;theta = atan(b/(a-c)/2;theta = theta*(180/pi); % 从幅度转换到角度5. 找区域边界function left, right, up, down = bianjie(A)m n = size(A);left = -1;right = -1;up
10、= -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;end;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
11、,if (A(i,j) = 0)down = i;break;end;end;if (down = -1)break;end;end;6. 求起始坐标function newcoord = checklimit(coord,maxval) newcoord = coord;if (newcoordmaxval) newcoord=maxval;end;7.模板匹配function ccorr, mfit, RectCoord = mobanpipei(mult, frontalmodel,ly,wx,cx, cy, angle)frontalmodel=rgb2gray(frontalmode
12、l);% 旋转模板% 求边界坐标% 选择模板人脸区域% 求质心model_rot = imresize(frontalmodel,ly wx,bilinear); % 调整模板大小model_rot = imrotate(model_rot,angle,bilinear); l,r,u,d = bianjie(model_rot);bwmodel_rot=imcrop(model_rot,l u (r-l) (d-u); modx,mody =center(bwmodel_rot);morig, norig = size(bwmodel_rot);% 产生一个覆盖了人脸模板的灰度图像 mfit
13、 = zeros(size(mult);mfitbw = zeros(size(mult); 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); endy = checklimit(endy,limy);for
14、 i=starty:endy, for j=startx:endx, mfit(i,j) = model_rot(i-starty+1,j-startx+1);end;end;ccorr = corr2(mfit,mult) % 计算相关度l,r,u,d = bianjie(bwmodel_rot);sx = startx+l;sy = starty+u;RectCoord = sx sy (r-1) (d-u); % 产生矩形坐标8.主程序clear;fname,pname=uigetfile(*.jpg;*.bmp;*.tif;*.gif,Please choose a colorpicture.);% 返回打开的图片名与图片路径名u,v=size(fname);y=fname(v); % 图片格式代表值switch ycase 0errordlg(You Should Load Image File First.,Warning.);caseg;G;p;P;f;F; % 图片格式若是 JPG/jpg、BMP/bmp、
