人脸分割 matlab程序

人脸脸部的皮肤具有一定的规律（都落在某个YCbCr表示的一个椭圆区域内），符合此规律的辨识为人脸

程序：（自己写， 但结果颜色很奇怪）

function [] = detection()
clear;
close all;
l = imread('2-2.jpg');
figure(1);
imshow(l);
l=double(l);
k=rgb2ycbcr(l);
imshow(k);
f=zeros(267,400);


for i= 1:267
    for j=1:400
%Y=0.257*l(i,j,1)+0.504*l(i,j,2)+0.098*l(i,j,3);
%Cb=0.148*l(i,j,1)-0.291*l(i,j,2)+0.439*l(i,j,3);
%Cr=0.439*l(i,j,1)-0.368*l(i,j,2)-0.071*l(i,j,3);
Y=k(i,j,1);
Cb=k(i,j,2);
Cr=k(i,j,3);
%Wci(Y)
if(Y<125)
Wcb =23+(Y-16)*(46.97-23)/(125-16);
Wcr =20+(Y-16)*(38.76-20)/(125-16); 
aCb =108+(125-Y)*(118-108)/(125-16);%aCi(Y)(Ci(Y)平均数）
aCr =154-(125-Y)*(154-144)/(125-16);
end
if(Y>188)
Wcb =14+(235-Y)*(46.97-14)/(235-188);
Wcr =10+(235-Y)*(38.76-10)/(235-188);
aCb =108+(Y-188)*(118-108)/(235-188);
aCr =154-(Y-188)*(154-132)/(235-188);
end
%Cia(Y)代表Ci'(Y)
if (125<=Y<=88)
    Cba=Cb; 
    Cra=Cr;
end
aCbh=108+(188-188)*(118-108)/(125-16);
aCrh=154-(188-188)*(154-132)/(235-188);
if Y<125 
    Cba=(Cb-aCb)*46.97/Wcb+aCbh;
    Cra=(Cr-aCr)*38.76/Wcr+aCrh;
end
if Y>188
    Cba=(Cb-aCb)*46.97/Wcb+aCbh;
    Cra=(Cr-aCr)*38.76/Wcr+aCrh;
end
%椭圆判断
x = cos(2.53)*(Cba-109.38)+sin(2.53)*(Cra-152.02);
y = -sin(2.53)*(Cba-109.38)+cos(2.53)*(Cra-152.02);
if ((x-1.6)^2/25.39^2+(y-2.41)*(y-2.41)^2/14.03^2<=1)
    f(i,j) = 255; 
else
    f(i,j)=0;
end
    end
end
%figure(2);
%imshow(f);   

运行结果：

原图：

结果：