次のように、単純な正弦ベル アポダイゼーション関数を作成しました。これをデータに適用してから、2D FFT を実行し、結果の等高線図を作成します。
P = 90
apd_fn = np.sin(((np.pi - (np.pi/180.0) * P) / 1024) * np.arange(0, 1024, 1) + (np.pi/180.0) * P)**2
fid_mat = ws_list * apd_fn
jres_spec = abs(fftshift(fft2(fid_mat)))
max_val = (jres_spec[:,400:]).max()
figure()
CS=contour(jres_spec,arange(0,max_val,max_val/30))
xlabel('frequency (ppm)')
show()
x 次元でウィンドウ化すると、次のような 2D 分光画像が生成されます。
同じフィルターを y 次元に適用したいと思っていますが、少し精神的なブロックがあり、これを実装する方法を見つけるのに苦労しています。あなたが提供できるどんな助けも大歓迎です。
ありがとう。