clear?all;
clc;
I=imread（‘IFA2.bmp‘）;
imshow（I）;
hold?on;
[mn]=ginput（1）;
n=round（n）;
m=round（m）;
plot（mn‘g+‘‘color‘[1?1?0]）;
hold?on;
[GmagGdir]=imgradient（I）;
[rowcol]=size（I）;
%--------向左下搜索------------
n1=n;
m1=m;
X1（1）=m;
Y1（1）=n;
z1=2;
while（std2（I（n1-3:n1+3m1-3:m1+3））>2.8）?%n1>9&&m1>9&&n1????for?i=n1-9:n1+9
????????flag=0;
????????for?j=m1-9:m1+9
????????????if（Gmag（ij）<10&&i>j+（n1-m1））
????????????????flag=1;
????????????????X1（z1）=j;Y1（z1）=i;z1=z1+1;
????????????????n1=i;
????????????????m1=j;
????????????????break;
????????????end
????????end
????????if（flag==1）
????????????break;
????????end
????end
end
plot（X1（1:z1-1）Y1（1:z1-1）‘g‘）;
hold?on;
%-------------向右上搜索----------------
n2=n;
m2=m;
X2（1）=m;
Y2（1）=n;
z2=2;
while（std2（I（n2-3:n2+3m2-3:m2+3））>2.8）?%n2>9&&m2>9&&n2????for?i=n2-9:n2+9
????????flag=0;
????????for?j=m2-9:m2+9
????????????if（Gmag（ij）<10&&i????????????????flag=1;
????????????????X2（z2）=j;Y2（z2）=i;z2=z2+1;
????????????????n2=i;
????????????????m2=j;
????????????????break;
????????????end
????????end
????????if（flag==1）
????????????break;
????????end
????end
end
plot（X2（1:z2-1）Y2（1:z2-1）‘g‘）;
hold?on;

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

?????文件???????1279??2015-11-04?11:06??fringetrace\FringeTrace.m

?????文件?????395576??2015-10-29?21:28??fringetrace\IFA2.bmp

?????目錄??????????0??2015-11-16?17:23??fringetrace

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

???????????????396855????????????????????3