人脸特征的选取和定位

A=imread('e:\21.jpg');a=double(A);c=mean(mean(a));B=imread('e:\2.jpg');b=double(B);fori=11:222forj=11:214p=b((i-10):(i+10),(j-10):(j+10));d=mean(mean(b));s1=sum(sum(a.*p))-441*c*d;s2=sum(sum(a.^2))-441*c^2;s3=sum(sum(p.^2))-441*d^2;s=sqrt(s2*s3);r(i,j)=s1/s;endend说明：21.jpg是