function?[eyeline]=eyefenge（grx）;
tempgrx=imread（‘1?（1）.bmp‘）;
[mn]=size（tempgrx）;
t1=graythresh（tempgrx）;?????????????????????%依據(jù)類內(nèi)方差最小法則求出該灰度圖像的最佳門限
bw=im2bw（tempgrxt1）;???????????????????????%將原灰度圖像二值化
????
%四向sobel算子的二值化
x=double（tempgrx）;
A=zeros（256256）;
for?i=2:255
????for?j=2:255
????????S1=abs（x（i-1j-1）+2*x（i-1j）+x（i-1j+1）-x（i+1j-1）-2*x（i+1j）-x（i+1j+1））;
????????S2=abs（x（i-1j-1）+2*x（ij-1）+x（i-1j-1）-x（i+1j+1）-2*x（ij+1）-x（i+1j+1））;
????????S3=abs（2*x（i-1j-1）+x（ij-1）+x（i-1j）-2*x（i+1j+1）-x（ij+1）-x（i+1j））;
????????S4=abs（2*x（i-1j+1）+x（ij+1）+x（i-1j）-2*x（i+1j-1）-x（ij-1）-x（i+1j））;
????????%A（ij）=S2;
????????gv=[S1S2S3S4];
????????A（ij）=max（gv）;
????????if?A（ij）>255
????????????A（ij）=255;
????????end
????end
end
A=255-A;
A=uint8（A）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3787??2010-12-06?21:13??eye.m

-----------?---------??----------?-----??----

?????????????????3787????????????????????1