Matlab实现简单的人脸识别程序

clear all;

close all;

F = imread(‘face9.jpg’);

F = rgb2gray(F); % 将图像转化为灰度图

level = graythresh(F);

BW = im2bw(F,level); % 将图像二值化

[n1 n2] = size(BW);

H = floor(n1/10); % 将图像分割成10个快

W = floor(n2/10);

h1 = 1;

h2 = H;

s = H*W; % 求得块的面积

h_i = H;

w_i = W;

for i = 1:10

w1 = 1;w2 =W; % 对应列初始化

for j = 1:10

if (w1 <= w_i || w2 >= 9*W) || (h1 <= h_i || h2 >= n1 - h_i) % 判断是否在图片四周的区域

loc = find(BW(h1:h2,w1:w2) == 0)

[num,~] = size(loc);

pr = num*100/s; % 计算灰度为黑色的像素占的比例

if pr < 100

BW(h1:h2,w1:w2) = 0;

end

% figure,imshow(BW);

% hold on

end

w1 = w1+W; % 跳到下一个块对应的列

w2 = w2+W;

end

h1 = h1+H; % 跳到下一个块对应的行

h2 = h2+H;

end

L = bwlabel(BW,8); % 区域标记

BB = regionprops(L,’BoundingBox’); % 得到包围的矩形框

Bd = cat(1,BB.BoundingBox);

[Bd_i,~] = size(Bd);

mx = 0;

for k = 1:Bd_i

p = Bd(k, 3)*Bd(k, 4);

% 如果满足面积块大，而且宽/高<1.8，则跳到Bd的该行

if p > mx && (Bd(k,3)/Bd(k,4)) < 1.8

mx = p;

j = k;

end

end

imshow(F); hold on; % 显示图片

% 画出矩形框

rectangle(‘Position’, Bd(j, :), …

‘EdgeColor’, ‘r’, ‘LineWidth’, 3);

title(‘标记图像’, ‘FontWeight’, ‘Bold’);