人臉識別MATLAB代碼

1、精品文檔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);%列像素for i=1:rowfor j=1:columnrr(i,j)=R1(i,j)/RGB(i,j);gg(i,j)=G1(i,j)/RG

2、B(i,j);endendrrr=mean(rr);r=mean(rrr);ggg=mean(gg);g=mean(ggg);2. 均值和協(xié)方差皮膚庫 1.jpg');r1,g1=rgb_RGB(t1);皮膚庫 2.jpg');r2,g2=rgb_RGB(t2);皮膚庫 3.jpg');r3,g3=rgb_RGB(t3);皮膚庫 4.jpg');r4,g4=rgb_RGB(t4);皮膚庫 5.jpg');r5,g5=rgb_RGB(t5);皮膚庫 6.jpg');r6,g6=rgb_RGB(t6);皮膚庫 7.jpg');r7,g7=rg

3、b_RGB(t7);皮膚庫 8.jpg');r8,g8=rgb_RGB(t8);。1歡迎下載精品文檔皮膚庫 9.jpg');r9,g9=rgb_RGB(t9);皮膚庫 10.jpg');r10,g10=rgb_RGB(t10);皮膚庫 11.jpg');r11,g11=rgb_RGB(t11);皮膚庫 12.jpg');r12,g12=rgb_RGB(t12);皮膚庫 13.jpg');r13,g13=rgb_RGB(t13);皮膚庫 14.jpg');r14,g14=rgb_RGB(t14);皮膚庫 15.jpg');r15,g

4、15=rgb_RGB(t15);皮膚庫 16.jpg');r16,g16=rgb_RGB(t16);皮膚庫 17.jpg');r17,g17=rgb_RGB(t17);皮膚庫 18.jpg');r18,g18=rgb_RGB(t18);皮膚庫 19.jpg');r19,g19=rgb_RGB(t19);皮膚庫 20.jpg');r20,g20=rgb_RGB(t20);皮膚庫 21.jpg');r21,g21=rgb_RGB(t21);皮膚庫 22.jpg');r22,g22=rgb_RGB(t22);皮膚庫 23.jpg');r

5、23,g23=rgb_RGB(t23);皮膚庫 24.jpg');r24,g24=rgb_RGB(t24);皮膚庫 25.jpg');r25,g25=rgb_RGB(t25);皮膚庫 26.jpg');r26,g26=rgb_RGB(t26);皮膚庫 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,

6、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);。2歡迎下載精品文檔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)

7、;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. 求偏轉(zhuǎn)角度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);b = b +

8、(j - xmean) * (i - ymean) * bw(i,j);c = c + (i - ymean)2 * bw(i,j);。3歡迎下載精品文檔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)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;bre

9、ak;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;。4歡迎下載精品文檔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)b

10、reak;end;end;6. 求起始坐標(biāo)function newcoord = checklimit(coord,maxval)newcoord = coord;if (newcoord<1)newcoord=1;end;。5歡迎下載精品文檔if (newcoord>maxval)newcoord=maxval;end;7. 模板匹配function ccorr, mfit, RectCoord = mobanpipei(mult,frontalmodel,ly,wx,cx, cy, angle)frontalmodel=rgb2gray(frontalmodel);model_

11、rot = imresize(frontalmodel,ly wx,'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,l u (r-l) (d-u);%選擇模板人臉區(qū)域modx,mody =center(bwmodel_rot);%求質(zhì)心morig, norig = size(bwmodel_rot);% 產(chǎn)生一個(gè)覆蓋了人臉模板的灰度圖像m

12、fit = 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);for i=

13、starty:endy,for j=startx:endx,mfit(i,j) = model_rot(i-starty+1,j-startx+1);。6歡迎下載精品文檔end;end;ccorr = corr2(mfit,mult)%計(jì)算相關(guān)度l,r,u,d = bianjie(bwmodel_rot);sx = startx+l;sy = starty+u;RectCoord = sx sy (r-1) (d-u); %產(chǎn)生矩形坐標(biāo)8. 主程序clear;fname,pname=uigetfile('*.jpg''*.bmp''*.tif'&

14、#39;*.gif','Pleasechoose a color picture.');%返回打開的圖片名與圖片路徑名u,v=size(fname);y=fname(v);%圖片格式代表值switch ycase 0errordlg('YouShould Load Image FileFirst.','Warning.');case'g''G''p''P''f''F'%圖片格式若是 JPG/jpg 、BMP/bmp、TIF/tif或者 GIF/

15、gif，才打開I=cat(2,pname,fname);Ori_Face=imread(I);subplot(2,3,1),imshow(Ori_Face);otherwiseerrordlg('You Should Load Image File First.','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);。7歡迎下載精品文

16、檔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ù)for i=1:rowfor j=1:columnif RGB(i,j)=0rr(i,j)=0;gg(i,j)=0;elserr(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)

17、9;);%皮膚概率服從高斯分布endendendsubplot(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 = ;for threshold = 0.55:-0.1:0.05two_value = zeros(i,j);two_value(find(image_low>threshold) = 1;cha

18、nge = sum(sum(two_value - previousSkin2);changelist = changelist change;previousSkin2 = two_value;。8歡迎下載精品文檔endC, 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);%顯示二值圖像我的照片人臉模板 .jpg');%

19、讀入人臉模板照片F(xiàn)aceCoord=;imsourcegray=rgb2gray(Ori_Face);%將原照片轉(zhuǎn)換為灰度圖像L,N=bwlabel(two_value,8);%標(biāo)注二值圖像中連接的部分,L 為數(shù)據(jù)矩陣， N為顆粒的個(gè)數(shù)for i=1:N,x,y=find(bwlabel(two_value)=i);%尋找矩陣中標(biāo)號為i 的行和列的下標(biāo)bwsegment = bwselect(two_value,y,x,8);%選擇出第 i 個(gè)顆粒numholes = 1-bweuler(bwsegment,4);%計(jì)算此區(qū)域的空洞數(shù)if (numholes >= 1)%若此區(qū)域至少包含

20、一個(gè)洞，則將其選出進(jìn)行下一步運(yùn)算RectCoord = -1;m n = 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, 'b

21、ilinear');bw = bwfill(bw,'holes');。9歡迎下載精品文檔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, RectC

22、oord =mobanpipei(justface,frontalmodel,ly,wx, cx,cy, angle);endif (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);%顯示覆蓋了

23、模板臉的灰度圖像imsourcegray = final_image;subplot(2,3,6);imshow(Ori_Face);%顯示檢測效果圖end;if (RectCoord = -1)FaceCoord = FaceCoord; RectCoord;endend。10歡迎下載精品文檔endend% 在認(rèn)為是人臉的區(qū)域畫矩形numfaces x = size(FaceCoord); for i=1:numfaces,hd = rectangle('Position',FaceCoord(i,:);set(hd, 'edgecolor', 'y&

24、#39;);end人臉檢測是人臉識別、人機(jī)交互、智能視覺監(jiān)控等工作的前提。近年來,在模式識別與計(jì)算機(jī)視覺領(lǐng)域 , 人臉檢測己經(jīng)成為一個(gè)受到普遍重視、研究十分活躍的方向。本文針對復(fù)雜背景下的彩色正面人臉圖像, 將膚色分割、模板匹配與候選人臉圖像塊篩選結(jié)合起來, 構(gòu)建了人臉檢測實(shí)驗(yàn)系統(tǒng) , 并用自制的人臉圖像數(shù)據(jù)庫在該系統(tǒng)下進(jìn)行了一系列的實(shí)驗(yàn)統(tǒng)計(jì)。本文首先介紹了人臉檢測技術(shù)研究的背景和現(xiàn)狀,闡明人臉檢測技術(shù)發(fā)展的重要意義, 對目前常用的一些檢測算法進(jìn)行了總結(jié), 然后著重闡述了基于膚色分割和模板匹配驗(yàn)證的人臉檢測算法。 膚色是人臉重要特征, 在通過膚色采樣統(tǒng)計(jì)和聚類分析后, 確立一種在 YCbCr 空間下的基于高斯模型的膚色分割方法。在YCbCr 色彩空間中建立膚色分布的高斯模型, 得到膚色概率似然圖像 , 在最佳動態(tài)閾值選取算法下完成膚色區(qū)域的分割。采用數(shù)學(xué)形態(tài)學(xué)和一些先驗(yàn)知識對所得到的膚色區(qū)域進(jìn)行人臉特征篩選, 進(jìn)一步剔除非人臉膚色區(qū)域, 減少候選人臉數(shù)量, 簡化后續(xù)檢測過程的處