matlab 识别调试,有关matlab的人脸识别程序，但调试是不成功

有关matlab的人脸识别程序，但调试是不成功，求高手帮忙指点修改.在此先谢了

1.色彩空间转换

function[r,g]=rgb_RGB(Ori_Face)

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;

row=size(Ori_Face,1);% 行像素

column=size(Ori_Face,2);% 列像素

fori=1:row

forj=1:column

rr(i,j)=R1(i,j)/RGB(i,j);

gg(i,j)=G1(i,j)/RGB(i,j);

end

end

rrr=mean(rr);

r=mean(rrr);

ggg=mean(gg);

g=mean(ggg);

2.均值和协方差

t1=imread('D:\matlab\皮肤库\1.jpg');[r1,g1]=rgb_RGB(t1);t2=imread('D:\matlab\皮肤库\2.jpg');[r2,g2]=rgb_RGB(t2);t3=imread('D:\matlab\皮肤库\3.jpg');[r3,g3]=rgb_RGB(t3);t4=imread('D:\matlab\皮肤库\4.jpg');[r4,g4]=rgb_RGB(t4);t5=imread('D:\matlab\皮肤库\5.jpg');[r5,g5]=rgb_RGB(t5);t6=imread('D:\matlab\皮肤库\6.jpg');[r6,g6]=rgb_RGB(t6);

t7=imread('D:\matlab\皮肤库\7.jpg');[r7,g7]=rgb_RGB(t7);

t8=imread('D:\matlab\皮肤库\8.jpg');[r8,g8]=rgb_RGB(t8);

t9=imread('D:\matlab\皮肤库\9.jpg');[r9,g9]=rgb_RGB(t9);

t10=imread('D:\matlab\皮肤库\10.jpg');[r10,g10]=rgb_RGB(t10);

t11=imread('D:\matlab\皮肤库\11.jpg');[r11,g11]=rgb_RGB(t11);

t12=imread('D:\matlab\皮肤库\12.jpg');[r12,g12]=rgb_RGB(t12);

t13=imread('D:\matlab\皮肤库\13.jpg');[r13,g13]=rgb_RGB(t13);

t14=imread('D:\matlab\皮肤库\14.jpg');[r14,g14]=rgb_RGB(t14);

t15=imread('D:\matlab\皮肤库\15.jpg');[r15,g15]=rgb_RGB(t15);

t16=imread('D:\matlab\皮肤库\16.jpg');[r16,g16]=rgb_RGB(t16);

t17=imread('D:\matlab\皮肤库\17.jpg');[r17,g17]=rgb_RGB(t17);

t18=imread('D:\matlab\皮肤库\18.jpg');[r18,g18]=rgb_RGB(t18);

t19=imread('D:\matlab\皮肤库\19.jpg');[r19,g19]=rgb_RGB(t19);

t20=imread('D:\matlab\皮肤库\20.jpg');[r20,g20]=rgb_RGB(t20);

t21=imread('D:\matlab\皮肤库\21.jpg');[r21,g21]=rgb_RGB(t21);

t22=imread('D:\matlab\皮肤库\22.jpg');[r22,g22]=rgb_RGB(t22);

t23=imread('D:\matlab\皮肤库\23.jpg');[r23,g23]=rgb_RGB(t23);

t24=imread('D:\matlab\皮肤库\24.jpg');[r24,g24]=rgb_RGB(t24);

t25=imread('D:\matlab\皮肤库\25.jpg');[r25,g25]=rgb_RGB(t25);

t26=imread('D:\matlab\皮肤库\26.jpg');[r26,g26]=rgb_RGB(t26);

t27=imread('D:\matlab\皮肤库\27.jpg');[r27,g27]=rgb_RGB(t27);

r=cat(1,r1,r2,r3,r4,r5,r6,r7,r8,r9,r10,r11,r12,r13,r14,r15,r16,r17,r18,r19,r20,r21,r22,r23,r24,r25,r26,r27);

g=cat(1,g1,g2,g3,g4,g5,g6,g7,g8,g9,g10,g11,g12,g13,g14,g15,g16,g17,g18,g19,g20,g21,g22,g23,g24,g25,g26,g27);

m=mean([r,g])

n=cov([r,g])

3.求质心

function[xmean, ymean] = center(bw)

bw=bwfill(bw,'holes');

area= bwarea(bw);

[mn] =size(bw);

bw=double(bw);

xmean=0; ymean = 0;

fori=1:m,

forj=1:n,

xmean= xmean + j*bw(i,j);

ymean= ymean + i*bw(i,j);

end;

end;

if(area==0)

xmean=0;

ymean=0;

else

xmean= xmean/area;

ymean= ymean/area;

xmean= round(xmean);

ymean= round(ymean);

end

4.求偏转角度

function[theta] = orient(bw,xmean,ymean)

[mn] =size(bw);

bw=double(bw);

a= 0;

b= 0;

c= 0;

fori=1:m,

forj=1:n,

a= a + (j - xmean)^2 * bw(i,j);

b= b + (j - xmean) * (i - ymean) * bw(i,j); c = c + (i - ymean)^2 * bw(i,j);

end;

end;

b= 2 * b;

theta= atan(b/(a-c))/2;

theta= theta*(180/pi); % 从幅度转换到角度

5.找区域边界

function[left, right, up, down] = bianjie(A)

[mn] = size(A);

left= -1;

right= -1;

up= -1;

down= -1;

forj=1:n,

fori=1:m,

if(A(i,j) ~= 0)

left= j;

break;

end;

end;

if(left ~= -1) break;

end;

end;

forj=n:-1:1,

fori=1:m,

if(A(i,j) ~= 0)

right= j;

break;

end;

end;

if(right ~= -1) break;

end;

end;

fori=1:m,

forj=1:n,

if(A(i,j) ~= 0)

up= i;

break;

end;

end;

if(up ~= -1)

break;

end;

end;

fori=m:-1:1,

forj=1:n,

if(A(i,j) ~= 0)

down= i;

break;

end;

end;

if(down ~= -1)

break;

end;

end;

6.求起始坐标

functionnewcoord = checklimit(coord,maxval) newcoord = coord;

if(newcoord<1)

newcoord=1;

end;

if(newcoord>maxval)

newcoord=maxval;

end;

7.模板匹配

function[ccorr, mfit, RectCoord] = mobanpipei(mult, frontalmodel,ly,wx,cx, cy, angle)frontalmodel=rgb2gray(frontalmodel);

model_rot= imresize(frontalmodel,[ly wx],'bilinear'); % 调整模板大小

model_rot= imrotate(model_rot,angle,'bilinear'); % 旋转模板

[l,r,u,d]= bianjie(model_rot); % 求边界坐标

bwmodel_rot=imcrop(model_rot,[lu (r-l) (d-u)]); % 选择模板人脸区域

[modx,mody]=center(bwmodel_rot); % 求质心

[morig,norig] = size(bwmodel_rot);

% 产生一个覆盖了人脸模板的灰度图像

mfit= zeros(size(mult));

mfitbw= zeros(size(mult));

[limy,limx] = size(mfit);

% 计算原图像中人脸模板的坐标

startx= cx-modx;

starty= cy-mody;

endx= startx + norig-1;

endy= starty + morig-1;

startx= checklimit(startx,limx);

starty= checklimit(starty,limy);

endx= checklimit(endx,limx);

endy= checklimit(endy,limy);

fori=starty:endy,

forj=startx:endx,

mfit(i,j)= model_rot(i-starty+1,j-startx+1);

end;

end;

ccorr= corr2(mfit,mult) % 计算相关度

[l,r,u,d]= bianjie(bwmodel_rot);

sx= startx+l;

sy= starty+u;

RectCoord= [sx sy (r-1) (d-u)]; % 产生矩形坐标

8.主程序

clear;

[fname,pname]=uigetfile({'*.jpg';'*.bmp';'*.tif';'*.gif'},'Pleasechoose a color picture...'); % 返回打开的图片名与图片路径名

[u,v]=size(fname);

y=fname(v);% 图片格式代表值

switchy

case0

errordlg('YouShould 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('YouShould 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);% 列像素数

fori=1:row

forj=1:column

ifRGB(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= [];

forthreshold = 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('E:\我的照片\人脸模板.jpg');% 读入人脸模板照片

FaceCoord=[];

imsourcegray=rgb2gray(Ori_Face);% 将原照片转换为灰度图像

[L,N]=bwlabel(two_value,8);% 标注二值图像中连接的部分,L为数据矩阵，N为颗粒的个数

fori=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;

[mn] = 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

% 在认为是人脸的区域画矩形

[numfacesx] = size(FaceCoord);

fori=1:numfaces,

hd= rectangle('Position',FaceCoord(i,:));

set(hd,'edgecolor', 'y');

end