四步相移法-解相

%% 清理空间
clc;
clear;
close all;

%% 读入并显示图像
X1=imread（‘o_1000000.bmp‘）;
X2=imread（‘o_1000001.bmp‘）;
X3=imread（‘o_1000002.bmp‘）;
X4=imread（‘o_1000003.bmp‘）;

figure（1）;I1=imresize（X11‘bilinear‘）;imshow（I1）;
figure（2）;I2=imresize（X21‘bilinear‘）;imshow（I2）;
figure（3）;I3=imresize（X31‘bilinear‘）;imshow（I3）;
figure（4）;I4=imresize（X41‘bilinear‘）;imshow（I4）;
 
%% 归一化【0-1】
I1=double（I1）;
I2=double（I2）;
I3=double（I3）;
I4=double（I4）;

%% 四步相移法，计算相位并显示
for j=1:615
    for i=2:576
        mask_0 = I2（ij）-I4（ij）;
        mask_1 = I1（ij）-I3（ij）;
        if mask_1 ~= 0
            phase（ij）=atan2（mask_