人脸检测中的皮肤检测程序

close all clear all clc %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%读入图像 I=imread('41.jpg'); %读入图像 [Ix,Iy,Iz]=size(I); if Ix>400&Iy>300 I=imresize(I,[400,Iy*400/Ix],'nearest');%图像大小变换 end figure imshow(I) title('原始图像') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 皮肤检测 I=double(I); %整数转化为小数 [hue,s,v]=rgb2hsv(I); %RGB颜色空间转化为HSV颜色空间 YCBCR = rgb2ycbcr(I);; %RGB颜色空间转化为YCbCr颜色空间 cb=0.148*I(:,:,1)-0.291*I(:,:,2)+0.439*I(:,:,3)+128; cr=0.439*I(:,:,1)-0.368*I(:,:,2)-0.071*I(:,:,3)+128; %RGB颜色空间转化为YCrCb颜色空间 [w h]=size(I(:,:,1)); %图像大小 for i=1:w for j=1:h if 145<=cr(i,j)&cr(i,j)<=170&145<=cb(i,j)&cb(i,j)<=195&0.01<=hue(i,j)&hue(i,j)<=0.15 %分离出黄色皮肤区域 segment(i,j)=1; %皮肤区域 else segment(i,j)=0; end end end figure imshow(segment); skin=segment; skin=bwareaopen(skin,round(w*h/900));%去除少像素连通区域 se=strel('disk',5); skin=imdilate(skin,se); %膨胀作用 im(:,:,1)=I(:,:,1).*skin; im(:,:,2)=I(:,:,2).*skin; im(:,:,3)=I(:,:,3).*skin; figure imshow(uint8(im)); title('皮肤区域')