%RGB?（R:紅、G:綠、B:藍(lán)）空間進(jìn)行膚色分割
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）;%開(kāi)運(yùn)算，去除像素點(diǎn)少于1000的區(qū)域
figure;
imshow（imf）;
[Lnum]=bwlabel（imf4）;?????%連通區(qū)域標(biāo)記?
B=zeros（size（imf））;
for?i=1:num
????Area（i）=bwarea（L==i）;%計(jì)算每個(gè)皮膚區(qū)域的面積
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?
????%計(jì)算各矩形區(qū)域面積
????Rect_Area（i）=（down（i）-up（i））*（right（i）-left（i））;
end
%計(jì)算各區(qū)域的填充率
Ratio=Area./Rect_Area;
for?i=1:num?
???if?Ratio（i）>=0.5%若相應(yīng)區(qū)域的填充率大于0.5則保留該區(qū)域
???????[xy]=find（L==i）;%第i塊區(qū)