肤色检测matlab

matlab肤色检测, ycbcr效果一般

clear;close all;
rgb_image=imread（‘E:\9.jpg‘）;
z=~（~rgb_image）;
ycbcr_image=double（rgb2ycbcr（rgb_image））;
ycbcr_image=ycbcr_image.*z;
[rcz]=size（ycbcr_image）;
hsv_image=rgb2hsv（rgb_image）;
a=25.39;b=14.03;ecx=1.60;ecy=2.41;cx=109.38;cy=152.02;x=0;y=0;
sinr=sin（2.53）;cosr=cos（2.53）;
for i=1:1:r
    for j=1:1:c       
        imgy=ycbcr_image（ij1）;
        cb=ycbcr_image（ij2）;
        cr=ycbcr_image（ij3）;
        x=cosr*（cb-cx）+sinr*（cr-cy）;
        y=-sinr*（cb-cx）+cosr*（cr-cy）;
        if imgy<60 ycbcr_image（ij:）=0;
        elseif imgy>=60&&imgy<=230
            if （（x-ecx）^2）/（a^2）+（（y-ecy）^2）/（b^2）>1
                ycbcr_image（ij:）=0;
            end
        else 
            if （（x-ecx）^2）/（（1.2*a）^2）+（（y-ecy）^2）/（（1.2*b）^2）>1
                ycbcr_image（ij:）=0;
            end
        end          
    end
end
subplot（121）;imshow（rgb_image）;
subplot（122）;imshow（ycbcr_image）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    .CA....       937  2010-04-15 20:40  CbCrEllipse\CbCrEllipse.asv

    .CA....       938  2010-04-15 20:59  CbCrEllipse\CbCrEllipse.m

    .C.D...         0  2010-04-21 10:49  CbCrEllipse

----------- ---------  ---------- -----  ----

                 1875                    3