四步移相位解包程序

matlab代码，四步相移图进行四步相移测量的程序

clear all
clf
clc
I1=zeros（256256）;
I2=zeros（256256）;
I3=zeros（256256）;
I4=zeros（256256）;

figure（1）;        %生成四幅干涉图并显示 
a=-80:0.1:80;
b=meshgrid（a）;
I1=cos（b-pi）;
%Ia（::1）=I1;Ia（::2）=0;Ia（::3）=0;（red）

imshow（I1）;

figure（2）;
I2=cos（b-pi/2）;
imshow（I2）;

figure（3）;
I3=cos（b）;
imshow（I3）;

figure（4）;
I4=cos（b+pi/2）;
imshow（I4）;
for j=1:256
    for i=1:256
phase（ij）=atan2（I2（ij）-I4（ij）I1（ij）-I3（ij））;   
 end
    end

%四步相移法计算出相位
          
figure（5）;
imshow（mat2gray（phase））;

figure（6）;
imshow（phase）;

n=zeros（256256）;   %解包
n（11）=0;
for i=2:256
    if abs（phase（1i）-phase（1i-1））        n（1i）=n（1i-1）;
    elseif phase（1i）-phase（1i-1）<=-pi
        n（1i）=n（1i-1）+1;
    elseif phase（1i）-phase（1i-1）>=pi
        n（1i）=n（1i-1）-1;
    end
end

for i=1:256
    for j=2:256
        if abs（phase（ji）-phase（j-1i））            n（ji）=n（j-1i）;
        elseif phase（ji）-phase（j-1i）<=-pi
            n（ji）=n（j-1i）+1;
        elseif phase（ji）-phase（j-1i）>=pi
            n（ji）=n（j-1i）-1;
        end
    end
end    

pphase=phase+2*pi.*n;
figure（7）;
imshow（mat2gray（pphase））;
figure（8）;
surf（pphase（2:end-12:end-1））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    1327158  2006-07-04 21:55  四步相移法程序\a.bmp

     文件    1327158  2006-07-04 21:56  四步相移法程序\b.bmp

     文件    1327158  2006-07-04 21:58  四步相移法程序\c.bmp

     文件    1327158  2006-07-04 21:57  四步相移法程序\d.bmp

     文件       1249  2013-06-17 21:00  四步相移法程序\ma.m

     文件       1393  2013-06-18 11:27  四步相移法程序\ma2.m

    ..A.SH.     12800  2013-06-16 10:39  四步相移法程序\Thumbs.db

     目录          0  2013-06-17 16:58  四步相移法程序

    ..A..HR    919600  2007-06-02 08:31  wry.dat

----------- ---------  ---------- -----  ----

              6243674                    9