人脸检测

function facedetec()
　　sum=0;
　　% clc,close all;
　　iptsetpref('ImshowBorder', 'tight')
　　f=imread('faces.jpg');
　　Ori_Face=f;
　　copy=f;
　　img=f;
　　% figure,imshow(f);
　　f=rgb2ycbcr(f);%rgb转换为ycbcr
　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
　　% figure,imshow(f);
　　f_cb=f(:,:,2);
　　f_cr=f(:,:,3);
　　f= (f_cb>=100) & (f_cb<=127) & (f_cr>=138) &(f_cr<=170) ;
　　%皮肤颜色在ycbcr色度空间的分布范围为:100<=cb<=127,138<=cr<=170
　　figure(1); imshow(f);
　　%K=medfilt2(f,[5 5]);%5×5中值滤波，去除噪声
　　%figure,imshow(f);
　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%利用开闭运算的方法消除噪声
　　se=strel('square',3);
　　f=imopen(f,se);
　　f=imclose(f,se);
　　figure(2),imshow(f);
　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%填孔处理
　　f=imfill(f,'holes');
　　figure(3),imshow(f);
　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%重构法
　　fe=imerode(f,ones(8,7));
　　fo=imopen(f,ones(8,7));
　　f=imreconstruct(fe,f);
　　figure(4),imshow(f);
　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相关匹配,作用不大,可注释掉
　　%w=ones(3);%w为模板大小.
　　%f=imfilter(f,w);
　　%figure(5),imshow(f);
　　% [M,N]=size(f);
　　% f=fft2(f);
　　% w=conj(fft2(w,M,N));
　　% g=real(ifft2(w.*f));
　　% gs=gscale(g);
　　% figure,imshow(gs);
　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
　　%人脸候选区筛选时,由于头部部分重合,以及头部与其他
　　%部分,例如衣服等的连接,对筛选造成了困难,故应先利用
　　%闭运算操作,断开连接,再进行处理.
　　se1=strel('square',8);
　　f=imerode(f,se1);
　　f=imdilate(f,se1);
　　figure(6),imshow(f);
　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%区域连通
　　[L,num]=bwlabeln(f,4);
　　for i=1:num;
　　[r,c]=find(L==i);
　　r_temp=max(r)-min(r);
　　c_temp=max(c)-min(c);
　　temp=size(r);
　　sum=sum+temp(1);
　　area_sq=r_temp*c_temp;
　　area=size(find(L==i),1);
　　ratio=area/area_sq;
　　if (r_temp/c_temp<0.8)|(r_temp/c_temp>2)|temp(1)>2000|temp(1)<200|ratio<0.6
　　%脸部区域<250的去掉,一般为手或其他干扰.
　　%利用脸部宽长比的大概上下限来确定一个模板范围.
　　%矩形面积area_sq=目标区长度*宽度,目标区面积为area,若area/area_sq<0.6,认为不是人脸区,删除之.
　　for j=1:temp(1);
　　L(r(j),c(j))=0;
　　end
　　else
　　continue;
　　end
　　end
　　L=bwperim(L,8);%边缘检测,检测出人脸的边缘区域
　　L=uint8(L);
　　z=find(L(:)>0);
　　L(z)=255;
　　% ave=sum/num;%ave=696.
　　figure(7),imshow(L);
　　L_r=L;
　　L_g=L;
　　L_b=L;
　　L_rgb=cat(3,L_r,L_g,L_b);%在原图上加框
　　% figure,imshow(L_rgb);
　　img1_r=min(L_r+img(:,:,1),255);
　　img1_g=min(L_g+img(:,:,2),255);
　　img1_b=min(L_b+img(:,:,3),255);
　　img1=cat(3,img1_r,img1_g,img1_b);
　　figure(8),imshow(img1);