怎么求图像相位

[xiaofei26]

我用FFT 求相位 但是程序死机

相位＝ arctan (I 虚部/实部R) 也不知道对不对

double arctan(double x)
{
int i;
double r,e,f,sqr;
sqr=x*x;
r=0;
e=x;
i=1;
while(e/i>1e-15)
{
f=e/i;
r=(i%4==1)?r+f:r-f;
e=e*sqr;
i+=2;
}
return 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)；//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);
}
}

[xiaofei26]

上面的方法可能不对

有没有更好的思路

我现在激光图像采集 显示

主要想要 图像相位图 和干涉图

有没有 搞过这方面 给我思路 不胜感激

有代码更好 谢谢

[xiaofei26]

有没有 高手呀 给点思路呀

[菜牛]

不知道你的代码什么意思，不过求反余切有现成的函数：
atan, atanf, atan2, atan2f

[xiaofei26]

那是傅里叶的算法里的

不用管它 现在我就想求相位图 和干涉图代码