基于肤色模型的人脸检测MATLAB代码

基于ycbcr空间的MATLAB人脸检测代码，非常简单，注释和代码都写得很好，新手也可以很容易看懂。

clear all
close all
clc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 读入待检测图像
RGB= imread（‘renlian2.jpg‘）;
subplot（221）imshow（RGB）;title（‘原图‘）;

I=rgb2ycbcr（RGB）;        %颜色空间转换 
[abc]=size（I）; %得到图像的像素点个数
cb=double（I（::2））;
cr=double（I（::3））;

%肤色模型
for i=1:a
    for j=1:b
        w=[cb（ij）cr（ij）];
        m=[117.4316 148.5599];
        n=[260.1301 12.1430;12.1430 150.4574];
        p（ij）=exp（（-0.5）*（w-m）*inv（n）*（w-m）‘）;%算某象素点的概率
        if （p（ij）<0.5） 
            p（ij）=0;
        else 
            p（ij）=1;
        end
    end
end 

%去噪
SE = strel（‘square‘8）; 
imf=imopen（pSE）;         %开运算（即先腐蚀再膨胀），消除杂散点 
subplot（222）imshow（imf）;title（‘肤色模型分割效果图‘）;

%%%%%%%%%%%%%%%%%%%%%%%根据填充率去除手脚、胳膊等非人脸区域%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[Lnum]=bwlabel（imf8）;     %连通区域标记 
B=zeros（size（imf））;
for i=1:num
    Area（i）=bwarea（L==i）;%计算每个皮肤区域的面积
end
for i=1:num
    [rc] = find（L==i） ;
    left（i）=min（c）;
    right（i）=max（c）;
    up（i）=min（r）;
    down（i）=max（r）;
end

%计算各矩形区域面积
for i=1:num 
    Rect_Area（i）=（down（i）-up（i））*（right（i）-left（i））;
end

%计算各区域的填充率
Ratio=Area./Rect_Area;
for i=1:num 
   if Ratio（i）>=0.5%若相应区域的填充率大于0.5则保留该区域
       [xy]=find（L==i）;%第i块区域的坐标值
       B=B+bwselect（imfyx8）;%把填充率大于0.5皮肤区域叠加起来   
   end      
end
%%%%%%%%%%%%%%%%%%%%%%%%%%根据面积比来进一步除去一些较小的非人脸区域%%%%%%%%%%%%%%%%%%%%%%%%%%
[L1num1]=bwlabel（B8）;     %连通区域标记
B1=zeros（size（B））;
for i=1:num1
    Area（i）=bwarea（L1==i）;%计算每个皮肤区域的面积
end
maxarea=max（Area）;%取最大值
q=Area/maxarea;%每块区域的面积与最大区域面积的比值   
for i=1:num1 
   if q（i）>=0.3%若相应区域的面积比值大于0.3则保留该区域
       [xy]=find（L1==i）;%第i块区域的坐标值
       B1=B1+bwselect（Byx8）;%把面积比值大于0.3皮肤区域叠加起来   
   end      
end
%%%%%%%%%%%%%%%%%%%%%%%%%%根据肤色区域的长宽比来除去一些非人脸区域%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[L2num2]=bwlabel（B18）;     %连通区域标记 
B2=zeros（size（B1））;
for i=1:num2
    [rc] = find（L2==i）;  
    left（i）=min（c）;
    right（i）=max（c）;
    up（i）=min（r）;
    down（i）=max（r）;
end
for i=1:num2
    if （（down（i）-up（i））/（right（i）-left（i）））>0.8&（（down（i）-up（i））/（right（i）-left（i）））<2
        [xy]=find（L2==i）;
        B2=B2+bwselect（B1yx8）;%%%把满足长宽比在0.8到2的区域留下
    end
end
subplot（223）imshow（RGB）;title（‘人脸定位‘）;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%把人脸框出来%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[L3 num3]=bwlabel（B28）;     %连通区域标记
for i=1:num3
    [rc] = find（L3==i）;  
    left（i）=min（c）;
    right（i）=max（c）;
    up（i）=min（r）;
    down（i）=max（r）;
end

hold on;
for i=1:num3
    if（down（i）>（up（i）+（right（i）-left（i））*1.2））     %人脸长宽比限制
       down（i）=up（i）+（right（i）-left（i））*1.2;
    end
    x=[left（i）;left（i）;right（i）;right（i）;left（i）];
    y=[up（i）;down（i）;down（i）;up（i）;up（i）];
    plot（xy）;      %画框
end
hold off;

a = min（x）;b = min（y）;
w = max（x） - min（x）;
h = max（y） - min（y）;
rect  = [a b w h];
img_cut = imcrop（RGBrect）;
subplot（224）imshow（img_cut）;title（‘人脸‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        3234  2019-05-24 12:44  fusemoxing.m