大家好
小弟刚入门影像处理fft这个领域
最近碰到一个问题想请教版上的高手
就是我有一张512x512影像(灰阶) f(x,y)
我进行fft2运算后可以得到512x512个在频率域中的复数值 F(u,v)
透过实部和虚部平方和开根号可以得到 大小函数值(magnitude) |F(u,v)|
透过tan^-1(虚部/实部) 可以得到相位角(phase angle) phi(u,v)
我想要测试如果用a影像的大小函数与b影像(按照上述也可以得到b影像的
大小函数值与相位角)的相位角组合 令其为F'(u,v)
F'(u,v)=|F(u,v)|*cos[pih(u,v)]+j*|F(u,v)|*sin[phi(u,v)]
影像a 影像b的 影像a 影像b的
的大小 相位 的大小 相位
利用ifft2将新组合逆转换回影像空间
但是转回来为何是非常不合理的黑白线条呢
程式码如下
请板上程式神手指点迷津
ein = imread('einstein.jpg');
a = fft2(ein,512,512);
e_m = abs(a); %大小函数值
e_a = atan(imag(a)./real(a))*180/pi; %相位角(角度)
e_r=e_m*cos(m_a); %实部(相位角为另一张影像的相位角)
e_i=e_m*sin(m_a); %虚部(相位角为另一张影像的相位角)
c=e_r+e_i*j;
ein2=uint8(ifft2(c));
figure(1);
imshow(c);
谢谢