grdfftgmt - クラウドでオンライン

プログラム：

NAME

grdfft - 波数 (または周波数) ドメインのグリッドで数学演算を行う

SYNOPSIS

grdft イングリッド [ ingrid2 ] アウトファイル [ 方位角 ] [ zレベル ] [[階段|g]] [[[r|x|y][w[k]]] [
[r|x|y]paramsは ] [[階段|g]] [
[f|q|s|nx / ny][+a|d|h|l][+e|n|m][+t幅][+w[サフィックス]] [+z[p]]] [ 階段 ] [[レベル]] [ -fg
]

注意： オプションフラグと関連する引数の間にスペースを入れることはできません。

DESCRIPTION

grdft 2-D 順方向高速フーリエ変換を使用して、XNUMX つ以上を実行します
変換して空間に戻す前の周波数領域での数学的演算
ドメイン。 新しい値を出力に書き込む前にデータをスケーリングするオプションが提供されています
ファイル。 グリッドの水平方向の寸法は、メートル単位であると想定されます。 地理的
グリッドは、 -fg 度をメートルにスケーリングするオプション。 あなたが持っている場合
km 単位の寸法のグリッドを使用して、これをメートルに変更できます。 グリッド編集 またはスケーリング
出力 算数.

REQUIRED 議論

イングリッド 操作対象の 2 次元バイナリ グリッド ファイル。 (以下のグリッドファイル形式を参照してください)。 為に
クロススペクトル操作、XNUMX 番目のグリッド ファイルも提供 イングルド2.

-Gアウトファイル
出力グリッド ファイルまたは 1 次元スペクトル テーブルの名前を指定します ( -E） （見る
以下のグリッドファイル形式)。

オプション 議論

-A方位角
で方向導関数をとります。 方位角 時計回りの角度で測定された方向
北から。

-Czレベル
上向き（ zレベル > 0) または下向き ( zレベル < 0) フィールドを続行 zレベル
メートル。

-NS[階段|g]
フィールドを微分します。つまり、d(フィールド)/dz を取ります。 これは乗算に相当します。
周波数領域の kr (kr は動径波数)。 にスケールを追加
(kr * を掛ける 階段） 代わりは。 または、追加 g あなたの
データはメートル単位のジオイド高であり、出力は mGal 単位の重力異常である必要があります。
[デフォルトはスケールなし]。

-E[r|x|y][w[k]]
動径方向のパワー スペクトルの推定 [r]。 場所 x or y 直後
-E 代わりに x または y 方向のスペクトルを計算します。 グリッドファイルはありません
作成した。 XNUMX つのグリッドが与えられた場合、f (つまり、周波数または波数), power[f],
とべき乗[f]の1つの標準偏差は、によって設定されたファイルに書き込まれます -G [標準出力]。 もしも
8 つのグリッドが与えられ、f と XNUMX つの量を書きます: Xpower[f]、Ypower[f]、コヒーレント
パワー[f]、ノイズ パワー[f]、位相[f]、アドミタンス[f]、ゲイン[f]、コヒーレンシ[f]。 各
数量の後に独自の 1-std dev エラー推定値が続くため、出力は 17 です。
列幅。 追加 w 周波数の代わりに波長を書きます。 あなたのグリッドが
さらに追加できる地理的 k メートル [デフォルト] から波長をスケーリングします。
km

-F[r|x|y]paramsは
データをフィルタリングします。 場所 x or y 直後 -F フィルタリングする x or y 方向のみ。
デフォルトは等方性 [r]。 コサイン テーパー バンドパス、ガウス バンドパスのいずれかを選択します。
バンドパス フィルター、またはバターワース バンドパス フィルター。

コサインテーパー:
XNUMXつの波長を指定 lc/lp/hp/hc 正しい単位で (参照 -fg) を設計する
バンドパス フィルター: より大きい波長 lc 以下 hc カットされます、
より大きい波長 lp 未満 hp 通過し、波長
間のコサイン テーパーになります。 例えば、 -F1000000 / 250000 / 50000 / 10000 -fg
> 1000 km および < 10 km の波長をカットし、通過します。
250 km から 50 km の間の波長。 ハイパスまたはローパス フィルターを作成するには、
ハイフン (-) を付ける hp/hc or lc/lp。 例えば、 -FX-/-/50/10 はローパスになります x,
通過波長 > 50 および拒絶波長 < 10。 -年度1000/250/-/-
ハイパスします y、通過波長 < 250、拒否波長 > 1000。

ガウシアン バンドパス:
追加 lo/hi、正しい単位の XNUMX つの波長 (参照 -fg) を設計する
バンドパスフィルター。 指定された波長では、ガウス フィルターの重みは次のようになります。
0.5 とする。 ハイパスまたはローパス フィルターを作成するには、ハイフン (-) を指定します。 hi
or lo 波長、それぞれ。 例えば、 -F-/30 は、
30 で半分の重みを持つガウス フィルター、 -F400/- はハイパス
データ。

バターワース バンドパス:
追加 lo/hi/注文、正しい単位の XNUMX つの波長 (参照 -fg）と
バンドパス フィルターを設計するためのフィルター次数 (整数)。 与えられた
バタワース フィルターの重みは 0.5 になります。 ハイパスを作るには
またはローパス フィルターには、ハイフン (-) を指定します。 hi or lo 波長、
それぞれ。 例えば、 -F-/30/2 は、2 次を使用してデータをローパスします。
バターワース フィルター、半分の重みを 30 に設定し、 -F400/-/2 はハイパス
データ。

-私[階段|g]
フィールドを積分します。つまり、integral_over_z (フィールド * dz) を計算します。 これは
周波数領域で kr で割ることに相当します (kr は動径波数)。
(kr * で除算するスケールを追加します。 階段） 代わりは。 または、追加 g 〜へ
データセットが mGal の重力異常であり、出力がジオイドであることを示します
メートル単位の高さ。 [デフォルトはスケールなし]。

-N [f | q | s |nx / ny] [+ a | [+ d | h | l] [+ e | n | m] [+ t幅] [+ w [サフィックス]] [+ z [p]]
FFTに適したグリッド寸法を選択または問い合わせ、オプションを設定
パラメーター。 FFT寸法を制御します。
-Nf FFTにデータの実際のディメンションを使用するように強制します。

-Nq より適切な寸法について問い合わせ、それらを報告してから続行します。

-NS オプションの寸法のリストが表示されたら、終了します。

-Nnx / ny 配列サイズでFFTを実行します nx / ny （> =グリッドファイルサイズである必要があります）。 デフォルト
FFTの速度と精度を最適化するデータ>=データを選択します。 FFTの場合
ディメンション>グリッドファイルのディメンション、データはゼロに拡張およびテーパーされます。

データのトレンド除去を制御する：線形トレンドを削除するための修飾子を追加します。
+d：データをトレンド除去します。つまり、最適な線形トレンドを削除します[デフォルト]。

+a：平均値のみを削除します。

+h：中間値、つまり0.5 *（最大+最小）のみを削除します。

+l：データはそのままにしておきます。

データの拡張とテーパーの制御：修飾子を使用して、拡張の方法を制御します
およびテーパーが実行されます。
+e エッジポイント対称を課すことによってグリッドを拡張します[デフォルト]、

+m エッジミラー対称性を課すことによってグリッドを拡張します

+n データ拡張をオフにします。

テーパリングは、データエッジからFFTグリッドエッジ[100％]まで実行されます。 変化する
このパーセンテージを介して +t幅。 時 +n 有効である場合、テーパーが適用されます
代わりに、拡張機能が利用できないため、データマージンに[0％]。

一時的な結果の書き込みを制御する：詳細な調査のために、
前方FFTに渡される中間グリッド。 これはおそらく
トレンド除去され、すべてのエッジに沿って点対称によって拡張され、テーパーが付けられます。 追加する
+w[サフィックス]そこから出力ファイル名が作成されます（つまり、 ingrid_prefix.ext)
[テーパー]、ここで EXT ファイル拡張子です。 最後に、複雑なグリッドを保存できます
を追加することにより、フォワードFFTによって生成されます +z。 デフォルトでは、実際の
架空のコンポーネント イングリッド_本物。EXT 　 イングリッド_imag。EXT。 追加する p 保存する
代わりに、ファイルに対する大きさと位相の極形式 イングリッド_mag。EXT 　
イングリッド_段階。EXT.

-S階段
各要素を掛ける 階段 空間領域 (周波数領域の後)
オペレーション）。 [デフォルトは 1.0]。

-V [レベル] （もっと ...）
詳細レベル[c]を選択します。

-fg 地理グリッド（経度、緯度の寸法）はメートルに変換されます
現在の楕円体パラメータを使用した「フラットアース」近似を介して。

-^ or ただ -
コマンドの構文に関する短いメッセージを出力してから終了します（注：Windowsの場合）
ただ使う -).

-+ or ただ +
任意の説明を含む広範な使用法（ヘルプ）メッセージを印刷します
モジュール固有のオプション（GMT共通オプションは除く）が終了します。

-? or いいえ 引数
オプションの説明を含む完全な使用法（ヘルプ）メッセージを印刷してから、
終了します。

- バージョン
GMTバージョンを印刷して終了します。

--show-datadir
GMT共有ディレクトリへのフルパスを出力して終了します。

GRID FILE 書式

デフォルトでは、GMTはCOARDSに苦情を申し立てるnetCDFで単精度浮動小数点数としてグリッドを書き出します。
ファイル形式。 ただし、GMTは他の多くの一般的に使用されるグリッドでグリッドファイルを生成できます
ファイル形式であり、浮動小数点を書き出す、いわゆるグリッドの「パッキング」も容易にします。
1バイトまたは2バイトの整数としてのデータ。 精度、スケール、オフセットを指定するには、ユーザーは次のことを行う必要があります
接尾辞を追加します =id[/階段/オフセット[/ナン]]、 どこ id グリッドのXNUMX文字の識別子です
タイプと精度、および 階段 　 オフセット オプションの倍率とオフセットは
すべてのグリッド値に適用され、 ナン 欠測データを示すために使用される値です。 万一に備えて
XNUMX人のキャラクター id のように提供されていません =/階段 よりも id=nf が想定されます。 いつ
グリッドを読み取ると、フォーマットは通常自動的に認識されます。 そうでない場合は、同じ接尾辞
入力グリッドファイル名に追加できます。 見る grdconvert およびセクションgrid-file-formatの
詳細については、GMTテクニカルリファレンスおよびクックブックを参照してください。

複数のグリッドを含むnetCDFファイルを読み取る場合、GMTはデフォルトで
そのファイルで見つけることができる最初の2次元グリッド。 GMTを別の読み物に誘導する
グリッドファイルの多次元変数、追加 ?変数名 ファイル名に、ここで
変数名 変数の名前です。 特別な意味から逃れる必要があるかもしれないことに注意してください
of ? シェルプログラムの前にバックスラッシュを置くか、シェルプログラムに
引用符または二重引用符の間のファイル名と接尾辞。 NS ?変数名 接尾辞も使用できます
出力グリッドの場合、デフォルトとは異なる変数名「z」を指定します。 見る
grdconvert GMTテクニカルのCFおよびグリッドファイル形式のセクション修飾子
特に3のスプライスの読み方に関する詳細については、リファレンスとクックブックを参照してください。
4次元または5次元のグリッド。

GRID 距離 単位

グリッドに水平単位としてメーターがない場合は、追加します +u単位 入力ファイルへ
指定された単位からメートルに変換する名前。 グリッドが地理的である場合は、変換します
供給によるメートルまでの距離 -fg を代わりにお使いください。

考慮事項

netCDFCOARDSグリッドは自動的に地理的として認識されます。 他のグリッドの場合
度をメートルに変換したい地理グリッド、選択 -fg。 データの場合
どちらかの極に近い場合は、グリッドファイルを長方形に投影することを検討する必要があります
を使用した座標系 グリッドプロジェクト

例

ファイル mag_0.nc の海面磁気異常をレベル 800 m まで上向きに継続するには
海抜：

gmt grdfft mag_0.nc -C800 -V -Gmag_800.nc

地理グリッド上の m 単位のジオイド高 (geoid.nc) を自由空気重力に変換するには
mGal の異常:

gmt grdfft geoid.nc -Dg -V -Ggrav.nc

mGal (faa.nc) の重力異常を垂直方向のたわみ (in
マイクロラジアン) 038 方向では、最初に重力を統合してジオイドを取得する必要があります。
方向導関数を取り、最後にラジアンをマイクロラジアンにスケーリングします。

gmt grdfft faa.nc -Ig -A38 -S1e6 -V -Gdefl_38.nc

重力異常の垂直方向の XNUMX 次導関数は、
分野。 2 回微分することにより、これらを mGal/m^XNUMX として計算できます。

gmt grdfftgravity.nc -D -D -V -Ggrav_2nd_derivative.nc

共同登録された海底地形と重力グリッドのクロススペクトル推定を計算するため、および
km 単位の波長の関数として結果を報告し、試してください。

gmt grdfft bathymetry.ncgravity.grd -Ewk -fg -V > cross_spectra.txt

トレンド除去、点対称反射、およびテーパリング後の FFT 前のグリッドを調べるには、
適用され、生のスペクトルの実数成分と虚数成分を保存するだけでなく、
topo.nc のデータを試してみてください

gmt grdfft topo.nc -N+w+z -fg -V

topo_taper.nc、topo_real.nc、および topo_imag.nc でデータのプロットを作成できるようになりました。

onworks.net サービスを使用してオンラインで grdfftgmt を使用する