求相位图 干涉图

xiaofei_net

我现在激光图像采集显示
主要想要图像相位图和干涉图 ，有没有搞过这方面给我思路，不胜感激 ，有代码更好   


我用FFT求相位   
但是程序死机
相位 arctan (I   虚部/实部R)         也不知道对不对      
for(y=0;y <h;y++)
{
for(x=0;x <w;x++)
{
m=arctan(FD[x*h+y].im/FD[x*h+y].re)；//m=arctan   (I   虚部/实部R)

if   (m> 255)
m=255;

Point((x <w/2?x+w/2:x-w/2),nHeight-1-(y <h/2?y+h/2:y-h/2))=(BYTE)(m);
}
}

[ 本帖最后由 zhangnan3509 于 2007-11-1 17:29 编辑 ]