matlab 肤色识别,MATLAB肤色检测.doc

MATLAB肤色检测

RGB=imread('1.jpg');

YCbCr=rgb2ycbcr(RGB); %??RGB??ê?×a?ˉ3éYCbCr??ê?

Y=YCbCr(:,:,1); %ìáè?áá?è

Cb=YCbCr(:,:,2);

Cr=YCbCr(:,:,3);

imshow(RGB);title('RGB'); %??ê??-í?

figure,imshow(YCbCr);title('YCbCr'); %??ê?YCbCrí?

I=RGB;

rows=size(YCbCr,1); %????DDêy

columns=size(YCbCr,2); %????áDêy

%ò????aó?í??2·?·¨?D?¨·?é???óòoí·?·?é???óò￡???DD?t?μ?ˉ

k=(2.53/180)*pi;

m=sin(k);n=cos(k);

cx=109.38;cy=152.02;ecx=1.60;ecy=2.41;a=25.39;b=14.03;

for i=1:rows

for j=1:columns

if Y(i,j)<80

I(i,j,:)=0;

elseif (Y(i,j)<=230&&Y(i,j)>=80)

x=(double(Cb(i,j))-cx)*n+(double(Cr(i,j))-cy)*m;

y=(double(Cr(i,j))-cy)*n-(double(Cb(i,j))-cx)*m;

if((x-ecx)^2/a^2+(y-ecy)^2/b^2)<=1

I(i,j,:)=255;

else I(i,j,:)=0;

end

elseif Y(i,j)>230

x=(double(Cb(i,j))-cx)*n+(double(Cr(i,j))-cy)*m;

y=(double(Cr(i,j))-cy)*n-(double(Cb(i,j))-cx)*m;

if ((x-ecx)^2/(1.1*a)^2+(y-ecy)^2/(1.1*b)^2)<=1

I(i,j,:)=255;

else I(i,j,:)=0;

end

end

end

end

figure,imshow(I); %??ê??t?μ?ˉoóμ?í???