face.m

基于matlab人脸检测的方法很多

function facedetec()
sum=0;
 clc,close all;
iptsetpref('ImshowBorder', 'tight')
f=imread('E:\vv.bmp');
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; 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,imshow(f);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%填孔处理
f=imfill(f,'holes');
 figure,imshow(f);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%重构法
fe=imerode(f,ones(8,7));
fo=imopen(f,ones(8,7));
f=imreconstruct(fe,f);
 figure,imshow(f);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相关匹配,作用不大,可注释掉
%w=ones(3);%w为模板大小.
%f=imfilter(f,w);
%figure,imshow(f);
 %[M,N]=size(f);
 %f=fft2(f);
 %w=conj(fft2(w,M,N));
 %g=real(ifft2(w.*f));
 %gs=scale(g);
 %figure,imshow(gs);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%人脸候选区筛选时,由于头部部分重合,以及头部与其他
%部分,例如衣服等的连接,对筛选造成了困难,故应先利用
%闭运算操作,断开连接,再进行处理.
se1=strel('square',8);
f=imerode(f,se1);
f=imdilate(f,se1);
figure,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);