clear?all;
I1=imread（‘H:\新建文件夾\8_5.bmp‘）;%讀入圖像
X1=double（I1）;
[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對圖像有影響?cann

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-10-16?15:58??人臉圖像歸一化\
?????目錄???????????0??2012-10-16?15:58??人臉圖像歸一化\1\
?????文件???????11382??2007-12-12?11:29??人臉圖像歸一化\1\8_5.bmp
?????文件????????2304??2008-07-04?06:56??人臉圖像歸一化\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