干涉图条纹追踪Matlab

包含有程序代码和干涉图片。运行程序在干涉图中选取某一暗条纹中心一点，追踪该条纹所有中心点。用matlab实现。

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