I1=double（l5）;
I2=double（l6）;
I3=double（l7）;
I4=double（l8）;
?[M?N]=size（I1）;


%參考平面的
Ii1=double（m5）;
Ii2=double（m6）;
Ii3=double（m7）;
Ii4=double（m8）;
?[Mi?Ni]=size（Ii1）;
for?i=1:M
?????I1（:i）=smooth（I1（:i）30）;
??????I2（:i）=smooth（I2（:i）30）;
???????I3（:i）=smooth（I3（:i）30）;
????????I4（:i）=smooth（I4（:i）30）;
????????Ii1（:i）=smooth（Ii1（:i）30）;
????????Ii2（:i）=smooth（Ii2（:i）30）;
????????Ii3（:i）=smooth（Ii3（:i）30）;
????????Ii4（:i）=smooth（Ii4（:i）30）;
?end
?
for?ji=1:512
????for?ii=2:Mi
phasei（ii