%%?清理空間
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_