结构光法的四步相移法求解相位差（matlab代码）

%读入第一组正弦照片

x1=imread('C:\Users\Administrator\Desktop\picture\1.png');
x2=imread('C:\Users\Administrator\Desktop\picture\2.png');
x3=imread('C:\Users\Administrator\Desktop\picture\3.png');
x4=imread('C:\Users\Administrator\Desktop\picture\4.png');

%实数化
x1=double(x1);
x2=double(x2);
x3=double(x3);
x4=double(x4);

%求解包裹相位
x5=atan2(x3-x1,x4-x2);

%在相应范围显示包裹相位图片
figure(1),imshow(x5,[])

%求解绝对相位
x6=unwrap(x5,pi,2);

%在相应范围显示包裹相位图片
figure(2),imshow(x6,[])

%读入第二组光线发生畸变的光栅图片
I1=imread('C:\Users\Administrator\Desktop\picture\a.png');
I2=imread('C:\Users\Administrator\Desktop\picture\b.png');
I3=imread('C:\Users\Administrator\Desktop\picture\c.png');
I4=imread('C:\Users\Administrator\Desktop\picture\d.png');