灰度投影的人眼识别matlab代码附带图像

灰度投影的人眼识别matlab代码（附带图像）

clear all;
[filenamepathname]=...
    uigetfile（{‘*.jpg‘:‘*.bmp‘:‘*.gif‘}‘         ‘）;%读入图像
str=[pathname filename];
l1=imread（str）;
X=rgb2gray（l1）;
X1=double（X）;
[mn]=size（X1）;
a=0;%垂直积分投影
f=0;d=0;PVD=0;
for j=1:n
    %for i=i0-（7-1）/2:i0+（7-1）/2
        %if i>=1&i<=512
    for i=1:m
     f=X1（ij）;
    d=d+f;
    %end
      %end
    end
    a=a+1;
    PVD（a）=d;
    d=0;
end
% plot（PVD）;
M=fix（n/2）;%向最接近零取整
[h1h2]=min（PVD（18:M））;
[h3h4]=min（PVD（1M+8:n））;
h5=h4+M;

a=0;%水平积分投影
f=0;d=0;PVD=0;
for i=1:m
    for j=1:n
     f=X1（ij）;
    d=d+f;
    end
    a=a+1;
    PVD（a）=d;
    d=0;
end
% % plot（PVD）;
N=fix（m/2）;

[w1w2]=max（PVD（11:N-10））;%w1是最小数，w2是坐标
[w3w4]=max（PVD（1N+1:m））;
w5=w4+N;
B=X1（w2:w5h2:h5）;
D=uint8（B）;
imshow（D）;


J=imadjust（D[0.150.35][01]）;%[0.20.3]对图像有影响， imadjust是线性灰度变换
BW=edge（J‘canny‘0.4）;%0.5对图像有影响 canny必须是在uint之后
figureimshow（BW）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      11382  2007-12-12 11:29  人眼定位\1\8_5.bmp

     文件        995  2010-05-25 09:36  人眼定位\1\face.m

     文件       2391  2010-05-24 15:52  人眼定位\1\face_1.asv

     文件       2406  2010-05-25 09:01  人眼定位\1\face_1.m

     文件        933  2008-02-15 15:49  人眼定位\1\hough_3.m

     文件      13454  2008-06-26 08:19  人眼定位\1\j2.bmp

     文件      13454  2008-06-26 08:19  人眼定位\1\j3.bmp

     文件      11382  2008-06-26 08:19  人眼定位\1\j4.bmp

    ..A.SH.     14336  2008-07-10 11:04  人眼定位\1\Thumbs.db

     目录          0  2010-05-25 09:36  人眼定位\1\新建文件夹

     目录          0  2010-05-25 09:36  人眼定位\1

     目录          0  2010-05-24 15:13  人眼定位

----------- ---------  ---------- -----  ----

                70733                    12