clear;
[fname,pname]=uigetfile({'*.jpg';'*.bmp';'*.tif';'*.gif'},'Please choose a color picture...');
% 返回打开的图片名与图片路径名
[u,v]=size(fname);
y=fname(v); % 图片格式代表值
switch y
case 0
errordlg('You Should Load Image File First...','Warning...');
case{'g';'G';'p';'P';'f';'F'}; % 图片格式若是JPG/jpg、BMP/bmp、TIF/tif或者GIF/gif,才打开
I=cat(2,pname,fname);
Ori_Face=imread(I);
subplot(2,3,1),imshow(Ori_Face);
otherwise
errordlg('You Should Load Image File First...','Warning...');
end
R=Ori_Face(:,:,1);
G=Ori_Face(:,:,2);
B=Ori_Face(:,:,3);
R1=im2double(R); % 将uint8型转换成double型处理
G1=im2double(G);
B1=im2double(B);
RGB=R1+G1+B1;
m=[ 0.4144,0.3174]; % 均值
n=[0.0031,-0.0004;-0.0004,0.0003]; % 方差
row=size(Ori_Face,1); % 行像素数
column=size(Ori_Face,2); % 列像素数
for i=1:row
for j=1:column
if RGB(i,j)==0
rr(i,j)=0;gg(i,j)=0;
else
rr(i,j)=R1(i,j)/RGB(i,j); % rgb归一化
gg(i,j)=G1(i,j)/RGB(i,j);
x=[rr(i,j),gg(i,j)];
p(i,j)=exp((-0.5)*(x-m)*inv(n)*(x-m)'); % 皮肤概率服从高斯分布
end
end
end
subplot(2,3,2);imshow(p); % 显示皮肤灰度图像
low_pass=1/9*ones(3);
image_low=filter2(low_pass, p); % 低通滤波去噪声
subplot(2,3,3);imshow(image_low);
% 自适应阀值程序
previousSkin2 = zeros(i,j);
changelist = [];
for threshold = 0.55:-0.1:0.05
two_value = zeros(i,j);
two_value(find(image_low>threshold)) = 1;
change = sum(sum(two_value - previousSkin2));
changelist = [changelist change];
previousSkin2 = two_value;
end
[C, I] = min(changelist);
optimalThreshold = (7-I)*0.1
two_value = zeros(i,j);
two_value(find(image_low>optimalThreshold)) = 1; % 二值化
subplot(2,3,4);imshow(two_value); % 显示二值图像
frontalmodel=imread('D:\PFK\人脸模板.jpg'); % 读入人脸模板照片
FaceCoord=[];
imsourcegray=rgb2gray(Ori_Face); % 将原照片转换为灰度图像
[L,N]=bwlabel(two_value,8); % 标注二值图像中连接的部分,L为数据矩阵,N为颗粒的个数
for i=1:N,
[x,y]=find(bwlabel(two_value)==i); % 寻找矩阵中标号为i的行和列的下标
bwsegment = bwselect(two_value,y,x,8); % 选择出第i个颗粒
numholes = 1-bweuler(bwsegment,4); % 计算此区域的空洞数
if (numholes >= 1) % 若此区域至少包含一个洞,则将其选出进行下一步运算
RectCoord = -1;
[m n] = size(bwsegment);
[cx,cy]=center(bwsegment); % 求此区域的质心
bwnohole=bwfill(bwsegment,'holes'); % 将洞封住(将灰度值赋为1)
justface = uint8(double(bwnohole) .* double(imsourcegray));
% 只在原照片的灰度图像中保留该候选区域
angle = orient(bwsegment,cx,cy); % 求此区域的偏转角度
bw = imrotate(bwsegment, angle, 'bilinear');
bw = bwfill(bw,'holes');
[l,r,u,d] =bianjie(bw);
wx = (r - l +1); % 宽度
ly = (d - u + 1); % 高度
wratio = ly/wx % 高宽比
if ((0.8<=wratio)&(wratio<=2))
% 如果目标区域的高度/宽度比例大于0.8且小于2.0,则将其选出进行下一步运算
S=ly*wx; % 计算包含此区域矩形的面积
A=bwarea(bwsegment); % 计算此区域面积
if (A/S>0.35)
[ccorr,mfit, RectCoord] = mobanpipei(justface,frontalmodel,ly,wx, cx,cy, angle);
end
if (ccorr>=0.6)
mfitbw=(mfit>=1);
invbw = xor(mfitbw,ones(size(mfitbw)));
source_with_hole = uint8(double(invbw) .* double(imsourcegray));
final_image = uint8(double(source_with_hole) + double(mfit));
subplot(2,3,5);imshow(final_image); % 显示覆盖了模板脸的灰度图像
imsourcegray = final_image;
subplot(2,3,6);imshow(Ori_Face); % 显示检测效果图
end;
if (RectCoord ~= -1)
FaceCoord = [FaceCoord; RectCoord];
end
end
end
end
% 在认为是人脸的区域画矩形
[numfaces x] = size(FaceCoord);
for i=1:numfaces,
hd = rectangle('Position',FaceCoord(i,:));
set(hd, 'edgecolor', 'y');
end