资源描述
创新性实验研究报告实验项目名称 matlab人脸识别、实验摘要本实验主要利用高斯肤色的归一化概率模型，将静态图像中的每一个像素限定在0到 1的范围内。利用人工阈值法将肤色与非肤色区域分开，形成二值图像。对上述图像进行 开闭运算，腐蚀膨胀等操作，弱化细节。通过脸部几何特征选取脸部区域，排除非脸部 区域。后利用脸部区域最大最小行列坐标得到最小外接矩形。进一步分析可得两眼，嘴 部坐标，以3点为基础画椭圆。最终圈定脸部区域。三、实验场地及仪器、设备和材料: 实验场地：J11-428实验室；仪器：PC机一台；材料：Matlab软件、待检测图片四、实验内容1、实验原理高斯肤色概率模型：由于统计表明不同人种的肤色区别主要受亮度信息影响，而受色度信息的影响较小， 所以直接考虑YCbCr空间的CbCr分量，映射为CbCr空间，在CbCr空间下，受亮度变化 的影响少，且是两维独立分布。通过实践，选取大量肤色样本进行统计，发现肤色在CbCr 空间的分布呈现良好的聚类特性。统计分布满足：77 Cb 127并且满足：133 Cr th)z(i,j)=1;elsez(i,j)=0;endendendfigure;imshow(z);title(Set threshold)%se=strel(square,3);f=imopen(z,se);f=imclose(f,se);%figure,imshow(f);%open and close processing; f=imfill(f,holes);%figure,imshow(f);%fill holes in the Img se1=strel(square,8);f=imerode(f,se1);f=imdilate(f,se1);%figure,imshow(f);%L,num=bwlabel(f,4);complextion probabilitythresholderosion and expansion;for i=1:num;%region loop;r,c=find(L=i);len=max(r)-min(r)+1;wid=max(c)-min(c)+1;area_sq=len*wid;area=size(r,1);for j=1:size(r,1)%pixel loop; if(len/wid2.4)|size(r,1)200|area/area_sq0.55L(r(j),c(j)=0;%not zero pixel =0;elsecontinue;endendendfigure;imshow(L);title(Face Region)%eliminate NOT facerengion;w=L&z;figure;imshow(w);title(Face and Features)%recover the features offacer c=find(L=0);r_min=min(r);r_max=max(r); c_min=min(c);c_max=max(c);figure;imshow(x);hold on %plot(round(.5*(c_max+c_min),round(.5*(r_max+r_min),+) flg=0;for i=round(.5*(r_min+r_max):-1:round(1/3*(r_max-r_min)+r_min)for j=round(.5*(c_max+c_min):-1:round(1/3*(c_max-c_min)+c_min) if(w(i,round(j)=0)flg=1;break;endendif(flg=1),break,endendw(i,j)if(w(i,j)=0)e1_x=j;e1_y=i;%plot(j,i,*)%Note left eyeflg=0;for i=round(.5*(r_min+r_max):-1:round(1/3*(r_max-r_min)+r_min)for j=round(.5*(c_max+c_min):round(2.2/3*(c_max-c_min)+c_min)if(w(i,j)=0)flg=1;break;endendif(flg=1)break,end;end w(i,j)e2_x=j;e2_y=i;%plot(j,i,*)%Note right eye for i=round(1/2*(r_max+r_min):round(r_min+2/3*(r_max-r_min)if(w(round(i),round(.5*(c_max+c_min)=0),break,endendw(i,round(.5*(c_max+c_min)m_x=.5*(c_max+c_min);m_y=i;%plot(round(.5*(c_max+c_min),i,*)%Note mouthox=(e1_x+e2_x+m_x)/3;oy=(e1_y+e2_y-e1_x+e2_x)/2;a=2*(e2_x-el_x);c=(el_y+e2_y)/2-m_y;b=(a*a-c*c)A0.5;t=0:.01:2*pi;x=ox+.65*b*cos(t);y=oy+.8*a*sin(t);plot(x,y,.)elserectangle(Position,c_min r_min c_max-c_min r_max-r_min,EdgeColor,r);title(Test Face Rec-ed)%test the face regionend
展开阅读全文