数学建模论文作业(人脸识别程序代码)

《数学建模论文作业(人脸识别程序代码)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、附录附录1Matlab人脸定位程序%%%%%%%%%%%%%%%%R型聚类分析%%%%%%%%%%%%%%functionRer=Re()clearallcloseallclcnum_pic=1;path_pic=strcat('E:1_LYC0_WHUT4_lab',num2str(num_pic),'.jpg');I=imread(path_pic);figure;imshow(I);holdon;title('原始图片','FontWeight','Bold');O=rgb2ntsc(I);G=O(:,:,2);[mn]=size(G);

2、U=zeros(m,n);fori=1:mforj=1:nifG(i,j)>0.03&&G(i,j)<0.20U(i,j)=1;endendendfigure;imshow(U);holdon;title('二值化图片','FontWeight','Bold');sr=strel('disk',4);C=imclose(U,sr);C_morph=bwmorph(C,'erod',0.8);[n1,n2]=size(C_morph);r=floor(n1/100);c=floor(n2/100);x1=1;x2=r;s=r*c;C1=C_morph;fori

3、=1:100y1=1;y2=c;forj=1:100loc=find(C1(x1:x1+r-1,y1:y1+c-1)==0);[p,q]=size(loc);pr=p/s*100;ifpr<=157C1(x1:x1+r,y1:y1+c)=1;elseC1(x1:x1+r,y1:y1+c)=0;endy1=y1+c;y2=y2+c;endx1=x1+r;x2=x2+r;endfigure;imshow(C1);holdon;title('增强图片','FontWeight','Bold');[L,num]=bwlabel(C1,4);B=regionprops

4、(L,'area');Se=[B.Area];[Se1,Se2]=size(Se);stats=regionprops(L,'BoundingBox');Bd=cat(1,stats.BoundingBox);[n1,n2]=size(C_morph);mx=n1*n2;P=[];num_area=1;fori=1:Se2ifSe(1,i)/mx<=0.005

5、Se(1,i)/(Bd(i,3)*Bd(i,4))<=0.50;Se(1,i)=0;elseif(Bd(i,3)/Bd(i,4))<1.8;P(num_area)=i;num_area=num_are

6、a+1;endendnum_area=num_area-1;figure;imshow(I);holdon;title('标记图像','FontWeight','Bold');I_gray=rgb2gray(I);fori=1:num_areai_position=P(i);rectangle('Position',Bd(i_position,:),'EdgeColor','r','LineWidth',3);endfori=1:num_areai_position=P(i);x1=floor(Bd(i_position,1));ifx1==0;x1=1;e

7、ndx2=floor(Bd(i_position,1)+Bd(i_position,3));y1=floor(Bd(i_position,2));7ify1==0;y1=1;endy2=floor(Bd(i_position,2)+Bd(i_position,4));I_select=I_gray(y1:y2,x1:x2);I_select=imresize(I_select,[112,92]);imwrite(I_select,strcat('E:ORLA',num2str(num_pic),'',num2str(i),'.pgm'));endend

8、7附录2Matlab灰度投影程序functionprojection()clearallclccloseall%%读入图像数据I=imread('E:ORLA23.pgm');[mn]=size(I);figure;imshow(I);fory=1:nS(y)=sum(I(1:m,y));S(y)=S(y)/n;endy=1:n;figuresubplot(211),plot(y,S(y));title('垂直投影');forx=1:mS(x)=sum(I(x,:));S(x)=S(x)/m;endx=1:m;subplot(212),plot(x,S

9、(x));title('水平投影');end7附录3