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