基于肤色的RGB多人脸检测

可以实现多人脸图像的检测，代码有注释，便于新手理解和学习。首先对人脸肤色的色素进行统计，如何在RGB空间下进行建立人脸肤色模型，最后使用MATLAB对其实现，准确度高，代码亲测可以实现。

%RGB （R:红、G:绿、B:蓝）空间进行肤色分割
I=imread（‘E:\tt.jpg‘）;
%function mm=rgb（I）
[widthheightdeep]=size（I）;
I=double（I）;
% rgbtoycbcr=[0.2990.5870.114;
%              -0.147-0.2890.436;
%              0.675-0.515-0.100];
% u=[156.560117.436];%#Cr#Cb
% C=[160.13012.143;
%       12.143299.457];
% d=inv（C）; 
%RGB to YCbCr
image=zeros（widthheight）;
for i=1:width 
    for j=1:height 
            R=I（ij1）;
            G=I（ij2）;
            B=I（ij3）;
            max=R;
            if（max                max=G;
            end   
            if（max                max=B;
            end
            min=R;
            if（min>G） 
                min=G;
            end  
            if（min>B）
                min=B;
            end
            if（R>98&&G>40&&B>20&&abs（R-G）>15&&R>G&&R>B&&max-min>15）
                image（ij）=255;
            else
                image（ij）=0;
            end               
    end
end
p=image;
se=strel（‘disk‘3）;
 se1=strel（‘disk‘1）;
 imf=imerode（pse）;
 imf=imdilate（imfse1）;
 [wh]=size（imf）;
 imf=bwareaopen（imfround（w*h/900））;
%  imf=bwareaopen（imf1000）;%开运算，去除像素点少于1000的区域
figure;
imshow（imf）;
[Lnum]=bwlabel（imf4）;     %连通区域标记 
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块区