肤色模型人脸识别matlab代码

clear all
close all
clc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 读入待检测图像
RGB= imread（‘1.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）=b