一个简单的人脸检测

<span style="font-size:18px;"><strong>I=imread('lena.jpg');
imshow(I)
% I=rgb2gray(I);
I=wiener2(I,[5 5]);
figure;
imshow(I);
bw=im2bw(I);
figure;
imshow(bw);
[m n]=size(bw);
r=floor(m/10);
c=floor(n/10);
x1=1;x2=r;
s=r*c;
%吧图像的边缘的十分之一区域进行去除
for i=1:10
    y1=1;y2=c;
    for j=1:10
        if(y2<=c | y2>=9*c)|(x1==1 | x2==r*10)
            bw(x1:x2,y1:y2)=0;
        end
        y1=y1+c;
        y2=y2+c;
    end
    x1=x1+r;
    x2=x2+r;
end                      %吧图像边缘十分之一的弄成黑色
figure;
imshow(bw);
L=bwlabel(bw,4);             %标注个联通区域
BB=regionprops(L,'BoundingBox');                 %计算包含这个连通区域的最小矩形的坐标
BB1=struct2cell(BB);                              %将个联通区域的坐标结构转换成元包
BB2=cell2mat(BB1);                                 %将元包转换成数组
[s1 s2]=size(BB2);
mx=0;
for k=3:4:s2-1
    p=BB2(1,k)*BB2(1,k+1);
    if p>mx & (BB2(1,k+1)/BB2(1,k))<2
        mx=p;
        j=k;
    end
end                                     %根据人脸区域长度与宽度的比小于2来确定

figure;
imshow(I);
hold on
rectangle('Position',[BB2(1,j-2) BB2(1,j-1) BB2(1,j) BB2(1,j+1)],'EdgeColor','r')
%用红矩形框画出人脸区域</strong></span>