基于膚色的人臉檢測matlab代碼

基于膚色的人臉檢測matlab代碼maincloseallclearallclc%輸入圖像名字img_name=input('請輸入圖像名字(圖像必須為RGB圖像，輸入0結(jié)束)：','s');%當輸入0時結(jié)束while~strcmp(img_name,'0')%進行人臉識別facedetection(img_name);img_name=input('請輸入圖像名字(圖像必須為RGB圖像，輸入0結(jié)束)：','s');endfacedetectionfunctionfacedetection(img_name)%讀取RGB圖像I=imread(img_name);%轉(zhuǎn)換為灰度圖像gray=rgb2gray(I);%將圖像轉(zhuǎn)化為YCbCr顏色空間YCbCr=rgb2ycbcr(I);%獲得圖像寬度和高度heigth=size(gray,1);width=size(gray,2);%根據(jù)膚色模型將圖像二值化fori=1:heigthforj=1:widthY=YCbCr(i,j,1);Cb=YCbCr(i,j,2);Cr=YCbCr(i,j,3);if(Y<80)gray(i,j)=0;elseif(skin(Y,Cb,Cr)==1)gray(i,j)=255;elsegray(i,j)=0;endendendend%二值圖像形態(tài)學處理SE=strel('arbitrary',eye(5));%gray=bwmorph(gray,'erode');%imopen先腐蝕再膨脹gray=imopen(gray,SE);%imclose先膨脹再腐蝕%gray=imclose(gray,SE);imshow(gray);%取出圖片中所有包含白色區(qū)域的最小矩形[L,num]=bwlabel(gray,8);STATS=regionprops(L,'BoundingBox');%存放經(jīng)過篩選以后得到的所有矩形塊n=1;result=zeros(n,4);figure,imshow(I);holdon;fori=1:numbox=STATS(i).BoundingBox;x=box(1);%矩形坐標xy=box(2);%矩形坐標yw=box(3);%矩形寬度wh=box(4);%矩形高度h%寬度和高度的比例ratio=h/w;ux=uint8(x);uy=uint8(y);ifux>1ux=ux-1;endifuy>1uy=uy-1;end%可能是人臉區(qū)域的矩形應(yīng)滿足以下條件：%1、高度和寬度必須都大于20，且矩形面積大于400%2、高度和寬度比率應(yīng)該在范圍（0.6,2）內(nèi)%3、函數(shù)findeye返回值為1ifw<20||h<20||w*h<400continueelseifratio<2&&ratio>0.6&&findeye(gray,ux,uy,w,h)==1%記錄可能為人臉的矩形區(qū)域result(n,:)=[uxuywh];n=n+1;endend%對可能是人臉的區(qū)域進行標記ifsize(result,1)==1&&result(1,1)>0rectangle('Position',[result(1,1),result(1,2),result(1,3),result(1,4)],'EdgeColor','r');else%如果滿足條件的矩形區(qū)域大于1則再根據(jù)其他信息進行篩選form=1:size(result,1)m1=result(m,1);m2=result(m,2);m3=result(m,3);m4=result(m,4);%標記最終的人臉區(qū)域ifm1+m3<width&&m2+m4<heigth<p="">rectangle('Position',[m1,m2,m3,m4],'EdgeColor','r');endendendfindeye%判斷二值圖像中是否含有可能是眼睛的塊%bImage----二值圖像%x---------矩形左上角頂點X坐標%y---------矩形左上角頂點Y坐標%w---------矩形寬度%h---------矩形長度%如果有則返回值eye等于1，否則為0functioneye=findeye(bImage,x,y,w,h)%根據(jù)矩形相關(guān)屬性得到二值圖像中矩形區(qū)域中的數(shù)據(jù)%存放矩形區(qū)域二值圖像信息part=zeros(h,w);%二值化fori=y:(y+h)forj=x:(x+w)ifbImage(i,j)==0part(i-y+1,j-x+1)=255;elsepart(i-y+1,j-x+1)=0;endendend[L,num]=bwlabel(part,8);%如果區(qū)域中有兩個以上的矩形則認為有眼睛ifnum<2eye=0;elseeye=1;endskin%AnilK.Jain提出的基于YCbCr顏色空間的膚色模型%根據(jù)當前點的CbCr值判斷是否為膚色functionresult=skin(Y,Cb,Cr)%參數(shù)%a=25.39;a=28;%b=14.03;b=18;ecx=1.60;ecy=2.41;sita=2.53;cx=109.38;cy=152.02;xishu=[cos(sita)sin(sita);-sin(sita)cos(sita)];%如果亮度大于230，則將長短軸同時擴大為原來的1.1倍if(Y>230)a=1.1*a;b=1.1*b;end%根據(jù)公式進行計算Cb=double(Cb);Cr=double(Cr);t=[(Cb-cx);(Cr-cy)];temp=xis