gtstfilt の解説

97/02/01 堀之内

● 機能

時空間 Fourier 変換＆逆変換によるフィルター

● オプション 

 以下で「東進」などの表現については README-gtstspct を参照のこと。

 sdim [I] sdim番目の軸をFourier変換を適用する「空間」軸とする(default=1)
 ewd  [L] T／F : 「東進」／「西進」成分を残す (default=T)
 mmin [I] 波数mmin(正値)以上を残す。0以下を与えると1を仮定 (default=-999)
 mmax [I] 波数mmax(正値)以下を残す。0以下を与えると最大波数-1を仮定(def.=-999)
      (注：Fourier変換の性質上 最大波数=IMAX/2 は東西進を区別できないので除外)
 fmin [I] 振動数fmin(正値)以上を残す。0以下を与えると1を仮定 (default=-999)
          (振動数も小さい方から何番目といった整数値で指定)
 fmax [I] 振動数mmax(正値)以下を残す。0以下だと最大振動数-1を仮定(def.=-999)
 m4 [I,4個]
 f4 [I,4個]
          波数振動数の範囲を4辺形で設定する。
          波数・振動数の４隅 (m4(i),f4(i)),i=1,..,4 は反時計回りになっ
          ていることを仮定する. すべての値に負でない値を与えたときのみ
          有効(defaultはすべて-999)。m4=2,10,10,2 などとコンマで区切っ
          て与える。
 help[L]  オプションの表示と簡単な説明。
 out,apnd,item,unit,title,dset,edit,ettl,greset: 他のgtoolコマンドと同様


● 備考

mmin,mmax,fmin,fmax による範囲指定と (m4(i),f4(i)) による範囲指定は同
時に行える。よって最大で8角形の切断が可能。

計算に用いる FFT の都合上、空間軸時間軸ともグリッド数が 2,3,5 で因数分
解できる偶数でなければならない。

● プログラマーズガイド

現状では波数ゼロ成分等、必ず削られる部分がある。これを残せるようにした
ほうが気分がいいかも知れない。

波数の指定が EWD と MMIN の二本立てになっているが、波数に負の値を許す
ようにすれば EWD は廃止できる。そうしたほうが使いやすいだろう。
