四步移相法提取相位（相位解包裹），初级问题

[sooAnderson]

最近在自学四步移相法提取相位，有个初级疑问：
有四幅移相图片a、b、c、d如下：完整解包裹代码如下：
程序运行后，主要是figure(7）和figure(8)没看太明白，;从figure(5）和figure(6）中可以看出被测物体的轮廓，可是从figure(7）和figure(8)中，怎么能看出物体的轮廓呢？figure(7）和figure(8)的解包是怎么解的？
由于是初学者，所以请您不吝赐教，谢谢。

　
clc;
clear all
clf
X1=imread('a.bmp');
X2=imread('b.bmp');
X3=imread('c.bmp');
X4=imread('d.bmp');
figure(1);I1=imresize(X1,1,'bilinear'); imshow(I1);
figure(2); I2=imresize(X2,1,'bilinear');imshow(I2);
figure(3); I3=imresize(X3,1,'bilinear'); imshow(I3);
figure(4); I4=imresize(X4,1,'bilinear'); imshow(I4);

[M N]=size(I1);
I1=double(I1);
I2=double(I2);
I3=double(I3);
I4=double(I4);

for j=1:615
for i=2:M
phase(i,j)=atan2(I2(i,j)-I4(i,j),I1(i,j)-I3(i,j));
end
end
%四步相移法计算出相位
figure(5);
imshow(mat2gray(phase));
figure(6);
imshow(phase);

n=zeros(M,615); %解包
n(1,1)=0;
for i=2:615
if abs(phase(1,i)-phase(1,i-1))<pi
n(1,i)=n(1,i-1);
elseif phase(1,i)-phase(1,i-1)<=-pi
n(1,i)=n(1,i-1)+1;
elseif phase(1,i)-phase(1,i-1)>=pi
n(1,i)=n(1,i-1)-1;
end
end

for i=2:576
for j=1:615
if abs(phase(i,j)-phase(i-1,j))<pi
n(i,j)=n(i-1,j);
elseif phase(i,j)-phase(i-1,j)<=-pi
n(i,j)=n(i-1,j)+1;
elseif phase(i,j)-phase(i-1,j)>=pi
n(i,j)=n(i-1,j)-1;
end
end
end

pphase=phase+2*pi.*n;
figure(7);
imshow(mat2gray(pphase));
figure(8);
surf(pphase(2:end-1,2:end-1));

[sooAnderson]

这里图片上传有点问题，可以参考下面的链接，谢谢。
http://zhidao.baidu.com/question/713095352075331125.html

[qq_30004127]

你说的figure5和6是截断相位（包裹相位），此时的相位通过atan2函数解得在在[-pi,pi]之间，所以后面就按点再次展开，获取每点相应的条纹级次n，带入pphase=phase+2*pi.*n;这就是最终相位了，也是解包裹相位

[qq_26702523]

楼主你好，请问你这里的解包裹用的什么方法？我最近也在学习四步相移法，谢谢指点。

[ZY520XUE]

有人还关注这个帖子吗

[吴帅杰]

你这个程序对吗？我在做三步相移取相位和解包裹，就是参考的你这个，我做出的图像不知道对不对

[weixin_43724679]

你好，我最近也在看这方面的。请问你为什么要放四张图呢，最后出来的这四张都是代表的同一张的包裹相位图和解包相位图吗？？

[qq_42920988]

引用 5 楼 吴帅杰 的回复:

你这个程序对吗？我在做三步相移取相位和解包裹，就是参考的你这个，我做出的图像不知道对不对

您好，我现在在做三步相移法取相位和解包裹，请问下您参考这个弄好了吗？

[weixin_41845437]

你好 问题解决了嘛 我最近也在学习这个 想请教一下