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

プログラム：

NAME

dimfilter - 空間 (または時間) ドメインでの 2-D グリッド ファイルの方向フィルタリング

SYNOPSIS

ディムフィルタ 入力ファイル.nc 距離フラグ [モード] 出力ファイル.nc
[ コルズ ] [ 増加 ] [ 地域 ] [] [[[レベル]] [ -f]

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

DESCRIPTION

ディムフィルタ をフィルタリングします .nc 与えられた
円にフィルターをかける n_sectors、選択したプライマリ畳み込みの XNUMX つを適用するか、
非畳み込みフィルターを各セクターに適用し、
選択された二次フィルター。 デカルトまたは球面ジオメトリを使用して距離を計算します。
出力 .nc ファイルは、必要に応じて、入力のサブ領域として、および/または
新製品 -I増分します。 このようにして、入力データに「余分なスペース」を持たせることができます。
出力グリッドのエッジ効果はありません。 フィルタがローパスの場合、出力
入力よりもサンプリング頻度が低い場合があります。 -Q エラー分析モード専用です。
入力ファイルの列の総数が必要です。これには、フィルター処理された
深さ。 最後に、知っておくべきこと ディムフィルタ 他のように滑らかな出力を生成しません
空間フィルターは、最小の中央値を返すため、これを行います N の中央値 N セクター。 の
入力データにノイズがない場合を除き、出力は粗くなる可能性があります。 したがって、追加のフィルタリング
(例えば、ガウシアン グリッドフィルター) DiM フィルター処理されたデータの一般的に推奨されます。

REQUIRED 議論

入力ファイル.nc
フィルタリングするデータ グリッド。

-D距離フラグ
距離 フラグ グリッド (x,y) がフィルターにどのように関連するかを示します 幅、 次のように：

フラグ = 0: グリッド (x,y) と同じ単位 幅、デカルト距離。 フラグ = 1: グリッド
(x,y) 度、 幅 キロメートル単位のデカルト距離。 フラグ = 2: グリッド (x,y)
度で、 幅 km 単位、cos(middle y) でスケーリングされた dx、デカルト距離。

上記のオプションは、重み行列のみを計算できるため、最速です。
一度。 次の XNUMX つのオプションは、それぞれの重みを再計算するため遅くなります。
緯度。

フラグ = 3: 度単位のグリッド (x,y)、 幅 km、cosine(y) でスケーリングされた dx、デカルト
距離計算。

フラグ = 4: 度単位のグリッド (x,y)、 幅 km、球面距離の計算。

-F[モード]
プライマリ フィルタ タイプを設定します。 畳み込みフィルターと非畳み込みフィルターから選択します。
フィルター コードを追加し、その後に全直径を追加します。 幅. 利用可能な畳み込み
フィルターは次のとおりです。

(b）ボックスカー：すべての重量は同じです。

(c）コサインアーチ：重みはコサインアーチ曲線に従います。

(g）ガウス：重みはガウス関数によって与えられます。

非畳み込みフィルターは次のとおりです。

(m）中央値：中央値を返します。

(p）最尤確率（モード推定量）：最頻値を返します。 もっとあれば
複数のモードが見つかった場合は、それらの平均値を返します。 フィルタに-または+を追加します
モーダル値の最小値または最大値を返したい場合は幅。

-N
二次フィルタ タイプと蝶ネクタイ セクタの数を設定します。 n_sectors しなければなりません
0 より大きい整数である必要があります。 n_sectors が 1 に設定されている場合、XNUMX 次フィルターは
効果がありません。 使用可能な二次フィルターは次のとおりです。

(l) Lower: フィルタリングされたすべての値の最小値を返します。

(u) Upper: フィルタリングされたすべての値の最大値を返します。

(a) 平均: フィルタリングされたすべての値の平均を返します。

(m) 中央値: フィルタリングされたすべての値の中央値を返します。

(p) モード: フィルタリングされたすべての値のモードを返します。

-G出力ファイル.nc
出力ファイル.nc フィルタの出力です。

オプション 議論

-I x_inc [およびオプションで y_inc] は出力インクリメントです。 追加 m 分を示すために、
or c 秒を表示します。 新しい場合 x_inc, y_inc の整数倍ではない
古いもの (入力データ内) では、フィルタリングがかなり遅くなります。 [デフォルト: 同じ
入力として。

-R 西, 東, 南, 北 出力ポイントの領域を定義します。 [デフォルト：
入力と同じ]

-T の反対になるように、出力グリッドのノード登録を切り替えます。
入力グリッド [デフォルトでは、入力グリッドと同じレジストレーションが得られます]。

-Qコルズ コルズ 入力テキスト テーブル ファイルの列の総数です。 このモードでは、
複数の列で構成される深さを読み取ることを想定しています。 各列は
によって取得できるフィルター幅を持つフィルター処理されたグリッド grd2xyz -Z. 結果
中央値、MAD、および平均になります。 したがって、中央値の列は生成に使用されます
地域コンポーネントと MAD を含む列は、エラーを実行するために使用されます
解析。

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

-f [i | o]コリン情報 （もっと ...）
入力列または出力列、あるいはその両方のデータ型を指定します。

-^ 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次元のグリッド。

地理 そして タイム コーディネート

出力グリッドタイプがnetCDFの場合、座標には「経度」というラベルが付けられます。
入力データまたはグリッド（存在する場合）の属性に基づく「緯度」または「時間」、または
-f or -R オプション。 たとえば、両方 -f0x -f1t & -R90w / 90e / 0t / 3tは結果として
経度/時間グリッド。 x、y、またはz座標が時間の場合、グリッドに格納されます
のTIME_UNITおよびTIME_EPOCHで指定されたエポックからの相対時間として gmt.conf file
またはコマンドラインで。 加えて 単位 時間変数の属性は
このユニットとエポックの両方。

例

north_pacific_dbdb5.nc が東経 5 度から東経 140 度までの 260 分間の海底地形のファイルであるとします。
0N ～ 50N で、半径 300km (全周 600km) 内の値の中央値を見つけたい
150E から 250E および 10N から 40N の範囲で選択する出力ポイントの幅)、および
0.5度ごとに出力値が必要です。 中央値が偏らないようにするには
傾斜面では、フィルター円を 6 つのセクターに分割し、
6 つの中央値の中で最も低い値。 球面距離の計算を使用するには、次のものが必要です。

gmt dimfilter north_pacific_dbdb5.nc -Gfiltered_pa​​cific.nc -Fm600 -D4 \
-Nl6 -R150/250/10/40 -I0.5 -V

cape_verde.nc が 0.5W から 32W まで、15N から 8N までの XNUMX 分海底地形のファイルであるとします。
25N のうねりを定義するために、長さの小さいフィーチャを削除する必要があります。
27.5W から 20.5W および 12.5N から 19.5N の範囲で、出力値を
2分。 デカルト距離計算を使用するには、次のものが必要です。

gmt dimfilter cape_verde.nc -Gt.nc -Fm220 -Nl8 -D2 -R-27.5/-20.5/12.5/19.5 -I2m -V
gmt grdfilter t.nc -Gcape_swell.nc -Fg50 -D2 -V

特定の領域のフィルター幅の範囲を見つけ、その範囲をフィルター処理したとします。
フィルター幅の範囲を使用して与えられた深浅データ (例えば、 f100.nc f110.nc f120.nc
f130.nc)、フィルター幅の範囲を使用して地域の傾向を定義したい場合、および
各データ ポイントで中央絶対偏差 (MAD) 推定値を取得する場合。 次にあなた
する必要があります:

gmt grd2xyz f100.nc -Z > f100.d
gmt grd2xyz f110.nc -Z > f110.d
gmt grd2xyz f120.nc -Z > f120.d
gmt grd2xyz f130.nc -Z > f130.d
貼り付け f100.d f110.d f120.d f130.d > depths.d
gmt dimfilter depths.d -Q4 > output.z

制限

地理 (緯度、経度) グリッドを操作する場合、XNUMX つの畳み込みフィルター (boxcar、
コサイン アーチ、およびガウス) は、次の変動に対してフィルターの重みを適切に正規化します。
緯度でグリッドボックスのサイズを変更し、どのノードが必要かを正しく判断します
フィルター「円」が周期的な (0-360) 境界を越えるか、または含まれている場合の畳み込み
地理的な極。 ただし、メディアン フィルターやモード フィルターなどの空間フィルターでは、
したがって、デカルト グリッド (または非常に低い緯度) でのみ使用する必要があります。
このような空間フィルターを適用する場合は、データを等面積に投影する必要があります
射影し、結果のデカルト グリッドで dimfilter を実行します。

SCRIPT TEMPLATE

dim.template.sh は、完全な DiM をセットアップするために使用できるスケルトン シェル スクリプトです。
MAD分析を含む分析。

REFERENCE

Kim, S.-S.、および Wessel, P. (2008)、Regional-Residual の方向性メディアン フィルタリング
海底地形の分離、 ジオケム。 地球物理学。 ジオシスト。, 9、Q03005、
土井：10.1029/2007GC001850。

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