フーリエ空間でグリッドを作成する

Question

私は、-2 arcsecから+2 arcsecまでの4x4 arcsecの正方形の画像を作成するコードを持っています。これは80x80グリッド上に作成されています。これに私は別のイメージを追加したい。 この第2の画像は、80×80グリッドのFFTによって作成され、したがってフーリエ空間ではじまります。 FFTの後、最初のイメージと同じ大きさのイメージを実空間に表示します。

フーリエ空間は尺度を表し、波数はk = 2pi/xと定義されています（この場合、numpy.fftはk = 1/xと考えるところの定義を使用します）。最小のk値と最小のk値を有するようにする。

その場合X_MAX = 2（第1の画像のx方向の寸法）とdim_x = 80（グリッド内の列の数）：

k_x、最大= 1 /（2 * X_MAX/dim_x）

k_x、分= 1 /（2 * X_MAX）

とk_xから実行フーリエ空間内のグリッドは、分

）、y方向についても同様maxは（k_xすることせ

私はこれを十分に説明して欲しいと思っていますが、私はリテラトゥでこれについて何らかの確認や説明を見つけることができませんでしたFFTのことについては本当に正しいかどうか知りたいと思っています。

ありがとうございます。

Answer 1

これは正しくありません。 k空間値は-N/2*omega_0から(N-1)/2*omega_0の範囲になります。ここで、omega_0はサンプル長の逆数で、2*pi/(max(x)-min(x))で与えられ、Nはサンプル数です。あなたの場合は、この行に沿って何かを得る：

N = len(x) 
dx = x[-1]-x[0] 
k = np.linspace(-N*pi/dx, (N+1)*pi/dx, N)