人脸识别.zip_人脸识别_图像识别_模式识别_特征识别_识别

clc; clear all; close all; x=imread('C:\Users\75654\Desktop\图像处理\1.jpg'); y=rgb2ycbcr(x); [a, b c]=size(y); cb=double(y(:,:,2)); cr=double(y(:,:,3)); for i=1:a; for j=1:b w=[cb(i,j) cr(i,j)]; m=[117.4316 148.5599]; n=[260.1301 12.1430;12.1430 150.4574]; p(i,j)=exp((-0.5)*(w-m)*inv(n)*(w-m)'); end end z=p./max(max(p));%---------------------------------------------------------complextion?probability? imshow(z) %figure;imshow(x);? th=0.5; for i=1:a for j=1:b if(z(i,j)>th) z(i,j)=1; else z(i,j)=0; end end end figure;imshow(z);title('Set?threshold')%-----------------------------------threshold? se=strel('square',3); f=imopen(z,se); f=imclose(f,se); %figure,imshow(f);open?and?close?processing;? f=imfill(f,'holes'); %figure,imshow(f);%fill?holes?in?the?Img? se1=strel('square',8); f=imerode(f,se1); f=imdilate(f,se1); %figure,imshow(f);%----------------------------------------------------erosion?and?expansion;? [L,num]=bwlabel(f,4); for i=1:num;%region?loop;? [r,c]=find(L==i); len=max(r)-min(r)+1; wid=max(c)-min(c)+1; area_sq=len*wid; area=size(r,1); for j=1:size(r,1)%pixel?loop;? if(len/wid<.8)||(len/wid>2.4)||size(r,1)<200||area/area_sq<0.55 L(r(j),c(j))=0;%not?zero?pixel?=0;? else continue; end end end figure;imshow(L);title('Face and Region')%----------------------------eliminate?NOT face rengion; w=L&z; figure;imshow(w);title('Face and Features')%-----------------------recover?the features of face [r, c]=find(L~=0); r_min=min(r);r_max=max(r); c_min=min(c);c_max=max(c); figure;imshow(x); hold on %plot(round(.5*(c_max+c_min)),round(.5*(r_max+r_min)),'+')? flg=0; for i=round(.5*(r_min+r_max)):-1:round(1/3*(r_max-r_min)+r_min) for j=round(.5*(c_max+c_min)):-1:round(1/3*(c_max-c_min)+c_min) if(w(i,round(j))==0) flg=1; break; end end if(flg==1),break,end end w(i,j); if(w(i,j)==0) e1_x=j;e1_y=i; %plot(j,i,'*')%--------------------------------------------------------Note?left?eye? flg=0; for i=round(.5*(r_min+r_max)):-1:round(1/3*(r_max-r_min)+r_min) for j=round(.5*(c_max+c_min)):round(2.2/3*(c_max-c_min)+c_min) if(w(i,j)==0),flg=1;break;end end if(flg==1),break,end; end w(i,j); e2_x=j;e2_y=i; %plot(j,i,'*')%--------------------------------------------------------Note?right?eye? for i=round(1/2*(r_max+r_min)):round(r_min+2/3*(r_max-r_min)) if(w(round(i),round(.5*(c_max+c_min)))==0),break,end end w(i,round(.5*(c_max+c_min))); m_x=.5*(c_max+c_min);m_y=i; %plot(round(.5*(c_max+c_min)),i,'*')%-----------------------------------------Note?mouth? ox=(e1_x+e2_x+m_x)/3;oy=(e1_y+e2_y-e1_x+e2_x)/2; a=2*(e2_x-e1_x);c=(e1_y+e2_y)/2-m_y;b=(a*a-c*c)^0.5; t=0:.01:2*pi; x=ox+.65*b*cos(t);y=oy+.8*a*sin(t); plot(x,y,'.') else rectangle('Position',[c_min r_min c_max-c_min r_max-r_min],'EdgeColor','r'); title('Test?Face?Rec-ed')%---------------------------------------------test?the?face?region? end