来自fft2的物理空间频率对天文图像的影响

Question

假设我的天文图像是512x512像素，像素比例为6弧分。我想以物理单位（arcmin ^ { - 1}）计算空间频率kx，ky，并绘制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))

这是对的吗？如果我想要它们是radians / arcmin的单位我想我也需要乘以2.0 * np.pi