标准四步相移条纹

matlab实现标准四步相移条纹，可自己设定条纹的高度和宽度；最后得到包含相位信息的条纹包裹图。

%%初始化参数
wid=500;%图片宽度
height=1000;%图片高度
T=40;%条纹周期
a=100;%正弦幅值
b=120;%正弦y轴平移值
%四幅相位差的图片
I1=zeros（widheight）;
I2=zeros（widheight）;
I3=zeros（widheight）;
I4=zeros（widheight）;
%% 产生四步相移光栅条纹
for i=1:wid
    for k=1:height
        I1（ik）=a*cos（2*pi*k/T+0）+b;
        I2（ik）=a*cos（2*pi*k/T+pi/2）+b;
        I3（ik）=a*cos（2*pi*k/T+pi）+b;
        I4（ik）=a*cos（2*pi*k/T+1.5*pi）+b;
    end
end
figure（1）;%显示四步相移光栅
subplot（221）;imshow（I1[]）;title（‘初始相位0‘）;
subplot（222）;imshow（I2[]）;title（‘初始相位0.5pi‘）;
subplot（223）;imshow（I3[]）;title（‘初始相位pi‘）;
subplot（224）;imshow（I4[]）;title（‘初始相位1.5pi‘）;
%% 四步相移解相
I=atan（（I4-I2）./（I1-I3））+pi;%待修正的四步相移提取相位
figure（2）;
imshow（I[]）;title（‘提取相位图‘）;
 %        imwrite（I‘C:\Users\nec\Desktop\标栅提取相位图.bmp‘‘bmp‘）;%imwrite有问题，保存的图像是二值图
% 显示解出来的一列
figure（3）;
plot（1:heightI（3:））;grid on;title（‘相位图结果分析（第三列）‘）;
hold on;
figure（4）;mesh（I）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1045  2013-05-31 10:36  zhou3.m