傅里叶变换 space

Question

我正在尝试制作一个简单的正弦波图像的动画。从逻辑上讲，这可以通过 Fourier-space 图像来实现，该图像具有恰好 1 个实部（像素）的功率，并逐渐改变像素对应的虚部（相位值）。下面的代码产生了一个非常逐渐向下移动的水平光栅。然而这种转变是缓慢的，一旦相位超过某个值，图像就不再具有动画效果。相位分量的单位和范围是什么？我在这里错过了什么？

import numpy as np, matplotlib as mpimg

f = (np.zeros([33, 33]) + 0j) # blank Fourier image
f[14, 16] = (1 + 0j) # a single component with power

steps = range(16) # number of images to produce
for i in range(len(steps)):
    c = np.fft.ifft2(np.fft.ifftshift(f, axes=[0, 1])) # grating image (cartesian)
    mpimg.imsave('grating %02i.png' % i, c, cmap='gray') # save image
    f[14, 16] += complex(0, 1/16.0) # increment phase

Answer 1

相位分量不仅仅是虚数分量。相位分量是 atan2(Im, Re)。要线性移动图像，您的增量可能应该将此相位修改 2*pi/nsteps。这意味着它必须同时修改实部和虚部。

傅里叶变换 space

translation in fourier space

python

numpy

fft