基于形态学处理算法的人脸识别matlab仿真,含仿真操作录像

clc; clear; close all; warning off; addpath(genpath(pwd)); rng('default') i=imread('face8.jpg'); I=rgb2gray(i);%转灰度图像 figure,imshow(i) hold on; [xt, yt] = meshgrid(round(linspace(1, size(I, 1), 10)), ... round(linspace(1, size(I, 2), 10))); mesh(yt, xt, zeros(size(xt)), 'FaceColor', ... 'None', 'LineWidth', 3, ... 'EdgeColor', 'r'); BW=im2bw(I);%转二值图像 figure,imshow(BW) figure,imshow(BW) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%% minimisation of background portion [n1 n2]=size(BW); r=floor(n1/10);%%即对n1/10向负方向舍入如floor(3.2)=3,floor(3.6)=3,floor(-4.2)=-5 c=floor(n2/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) loc=find(BW(x1:x2, y1:y2)==0); [o p]=size(loc); pr=o*100/s; if pr<=100 BW(x1:x2, y1:y2)=0; r1=x1;r2=x2;s1=y1;s2=y2; pr1=0; end imshow(BW); end y1=y1+c; y2=y2+c; end x1=x1+r; x2=x2+r; end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%% detection of face object L = bwlabel(BW,8);%%%函数bwlabel是把四连通或八连通的区域连接起来 %regionprops统计被标记的区域的面积分布,显示区域总数.返回值是结构数组,其相应域定义了每一个区域相应属性下的度量. BB = regionprops(L, 'BoundingBox');%‘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)/BB2(1,k+1))<1.8 mx=p; j=k; end end figure,imshow(I); hold on; %rectangle函数功能：创建二维矩形对象。 rectangle('Position',[BB2(1,j-2),BB2(1,j-1),BB2(1,j),BB2(1,j+1)],'EdgeColor','r' )%从点(x,y)开始绘制一个宽w高h的矩形，