天文图像上fft2的物理空间频率

Question

假设我有一张天文图像，它是512x512像素，像素比例为6弧分。我想用物理单位计算空间频率kx，ky（arcmin^{-1}），并将2D功率谱绘制成它们的函数。这就是我所做的

import numpy as np

PIXEL_SCALE = 6.0 # arcmin

N_PIXELS = 512
y, x = np.mgrid[0:N_PIXELS * PIXEL_SCALE:PIXEL_SCALE, \
   0:N_PIXELS * PIXEL_SCALE:PIXEL_SCALE]

k = 0.5
z = np.sin(2.0 * np.pi * k * y)

z_fft = np.fft.fftshift(np.fft.ifft2(np.fft.ifftshift(z)))

# FREQUENCIES
ky_fft = np.fft.fftshift(np.fft.fftfreq(y.shape[0], PIXEL_SCALE)) 
kx_fft = np.fft.fftshift(np.fft.fftfreq(x.shape[1], PIXEL_SCALE))

这是对的吗？如果我希望它们是弧度/弧分的单位，我想我还需要乘以2.0*np.pi在