これは、Ubuntu Online、Fedora Online、Windows オンライン エミュレーター、または MAC OS オンライン エミュレーターなどの複数の無料オンライン ワークステーションの XNUMX つを使用して、OnWorks 無料ホスティング プロバイダーで実行できるコマンド surfacegmt です。
プログラム：
NAME
サーフェス - 張力を調整可能な連続曲率スプラインを使用したグリッド テーブル データ
SYNOPSIS
表面 [ テーブル ] 出力ファイル.nc 増加 地域 [ アスペクト比 ] [ convergence_limit ] [
l下側 ] [ -ルーアッパー ] [ max_iterations ] [] [ 検索半径[m|s]] [[[i|b]張力係数
] [[レベル]] [ 過弛緩係数 ] [ -a] [ -bi] [ -du] [
-f] [ -h] [ -i] [ -:[i|o]]
注意： オプションフラグと関連する引数の間にスペースを入れることはできません。
DESCRIPTION
表面 標準入力からランダムな間隔の (x,y,z) トリプルを読み取ります [または テーブル』をプロデュース
以下を解くことにより、グリッド値 z(x,y) のバイナリ グリッド ファイルを作成します。
(1 - T) * L (L (z)) + T * L (z) = 0
ここで、T は 0 ～ 1 の張力係数、L はラプラシアン演算子を示します。 T = 0
SuperMISP および ISM と同等の「最小曲率」ソリューションを提供します。
パッケージ。 最小曲率は、望ましくない振動や誤った極大値または極大値を引き起こす可能性があります。
最小値 (Smith と Wessel、1990 年を参照)、これらを抑制するために T > 0 を使用することもできます。
効果。 経験上、潜在的なフィールド データと T には通常、T ~ 0.25 が適切であることが示唆されています。
急峻な地形データの場合は、より大きくする必要があります (T ~ 0.35)。 T = 1 は調和曲面を与えます (いいえ
最大値または最小値は、制御データ ポイントを除いて可能です)。 推奨されるのは、
ユーザーがデータを前処理する ブロック平均, ブロック中央値または ブロックモード 空間的なものを避けるために
エイリアスを作成し、冗長データを削除します。 下限および/または上限を課すことができます。
解決。 これらは、固定値、値を含むグリッド、または単純な形式で入力できます。
は入力データの最小値/最大値になります。 自然境界条件が適用されます。
エッジ（周期的な境界を適用する 360 度範囲の地理データを除く）
経度方向の様子。
REQUIRED 議論
-G出力ファイル.nc
出力ファイル名。 出力はバイナリ 2D です .nc ファイル。 最小のグリッドに注意してください。
次元は少なくとも 4 である必要があります。
-Iシンク[単位] [= | +] [/陰性[単位] [= | +]]
x_inc [およびオプションで y_inc]はグリッド間隔です。 オプションで、接尾辞を追加します
修飾子。 地理的 （度） 座標：追加 m アーク分を示すまたは s
アーク秒を示します。 ユニットのXNUMXつが e, f, k, M, n or u 追加されます
代わりに、増分はメートル、フィート、キロメートル、マイル、航海で与えられると想定されます
それぞれマイルまたは米国の測量フィートであり、同等のものに変換されます
地域の中緯度での経度（変換は
PROJ_ELLIPSOID）。 もしも /y_inc が与えられますが、0に設定すると、次のようにリセットされます。 x_inc;
それ以外の場合は、緯度に変換されます。 All 座標：もし = is
その後、対応する最大値を追加 x (東）または y (北）若干調整される場合があります
指定された増分に正確に適合するように[デフォルトでは、増分は調整される場合があります
与えられたドメインに少し合うように]。 最後に、増分を与える代わりに、
その 数 of ノード 追加して希望 + 提供された整数に
口論; 次に、増分はノードの数と
ドメイン。 結果の増分値は、選択したかどうかによって異なります。
グリッドライン登録またはピクセル登録グリッド。 詳細については、App-file-formatsを参照してください。
注： -Rグリッドファイル が使用される場合、グリッド間隔はすでに初期化されています。 使用する
-I 値を上書きします。
-NS[単位]xmin/xmax/イミン/ワイマックス[r] （もっと ...）
関心領域を指定します。
オプション 議論
テーブル XNUMXつ以上のASCII（またはバイナリ。を参照） -bi[ncols][type]）を保持するデータテーブルファイル
データ列の数。 テーブルが指定されていない場合は、標準入力から読み取ります。
-Aアスペクト比
アスペクト比。 必要に応じて、グリッドの異方性を方程式に追加できます。 入力
アスペクト比、ここで、dy = dx / アスペクト比 グリッドの寸法を関連付けます。 [デフォルト =
1 は等方性グリッドを想定しています。]
-Cconvergence_limit[％]
収束限界。 反復は、最大絶対値が得られた時点で収束したとみなされます。
グリッド値の変化が以下である convergence_limit。 (単位はデータzと同じ
単位）。 あるいは、% を追加して、rms 偏差のパーセンテージで制限を与えます。
[デフォルトは、データの二乗平均平方根偏差の 1e-4 にスケールされます。
最適な (最小二乗) 平面。]。 これは、次の時点での最終収束限界です。
望ましいグリッド間隔。 中間（粗い）グリッドの場合、効果的な収束
制限はグリッド間隔の乗数によって調整されます。
-Ll下側 -ルーアッパー
出力ソリューションに制限を課します。 l下側 下限を設定します。 下側 することができます
下限値、固定値を含むグリッド ファイルの名前。 d 最小に設定する
入力値、または u 制約なしの場合 [デフォルト]。 uアッパー 上限を設定し、
上限値、固定値を含むグリッド ファイルの名前です。 d に設定する
最大入力値、または u 制約なしの場合 [デフォルト]。 設定に使用されるグリッド ファイル
制限には NaN が含まれる場合があります。 NaN が存在する場合、ノードの制限は次のようにマスクされます。
NaN には制約がありません。
-Nmax_iterations
反復回数。 反復は次の場合に終了します convergence_limit に達するか、
反復回数が到達したとき max_iterations。 これが最後の反復です
必要なグリッド間隔で制限します。 中間（粗い）グリッドの場合は効果的です
反復制限は、グリッド間隔の乗数によって調整されます。 [デフォルトは500です。]
-Q 高度に複合的な最大公約数を持つグリッドの次元を提案します。 これ
表面がソリューション内でいくつかの中間ステップを使用できるようになり、結果がより速くなります
実行時間とより良い結果が得られます。 が提案するサイズ -Q を変更することで実現できます
-R および -I。 回復できます -R -I 後で使用して必要な場合 グリッドサンプル or
グリッドカット の出力について 表面.
-S検索半径[ミリ秒]
検索範囲。 入力 検索半径 X、Y データと同じ単位。 追加する m 示すために、
アーク分または s アーク秒間。 これは、実行前にグリッドを初期化するために使用されます。
最初の反復。 グリッド格子が素数でない限り、時間をかける価値はありません。
地域ステージを設けることはできません。 [デフォルト = 0.0 で、検索は行われません。]
-T[i|b]張力係数
張力係数[s]。 これらは 0 から 1 の間でなければなりません。張力は、
内部解 (スプリアス発振を抑制する上の方程式) と
境界条件 (境界条件に近づくと解が平坦化する傾向がある場合)
エッジ）。 両方の値にゼロを使用すると、自由な最小曲率サーフェスが得られます。
エッジ、つまり自然な双三次スプライン。 使用 -あなた張力係数 インテリアを設定する
緊張感と、 -Tb張力係数 境界張力を設定します。 先頭に付加しない場合 i or
b、両方が同じ値に設定されます。 [デフォルト = 0 で両方とも最小値が得られます
曲率の解決策。]
-V [レベル] （もっと ...）
詳細レベル[c]を選択します。 -V3 各反復後に収束を報告します。
-V 各地域グリッドが収束した後にのみレポートされます。
-Z過弛緩係数
過剰緩和要因。 このパラメータは、収束を加速するために使用されます。 それはです
1 と 2 の間の数値。値 1 は方程式を正確に反復し、
常に安定した収束を保証します。 値が大きいほど増分が過大評価されます
収束中に変化し、より迅速に解決に到達しますが、状況が悪化する可能性があります。
不安定。 この係数に大きな値を使用する場合は、次のことを監視することをお勧めします。
各反復では、 -Vl オプション。 [デフォルト = 1.4 はすぐに収束し、ほぼ
常に安定しています。]
-aコル=名[...] （もっと ...）
空間列の関連付けを設定する コル=名.
-bi [ncols][NS] （もっと ...）
ネイティブバイナリ入力を選択します。 [デフォルトは3入力列]。
-duデータなし （もっと ...）
等しい入力列を置き換えます データなし NaNで。
-f [i | o]コリン情報 （もっと ...）
入力列または出力列、あるいはその両方のデータ型を指定します。
-h [i | o] [n] [+ c] [+ d] [+ r発言] [+ rタイトル] （もっと ...）
ヘッダーレコードをスキップまたは生成します。 バイナリデータでは使用されません。
-iコルズ[l] [s階段] [oオフセット] [、...] （もっと ...）
入力列を選択します（0は最初の列です）。
-：[i | o] （もっと ...）
入力および/または出力の1番目と2番目の列を交換します。
-^ or ただ -
コマンドの構文に関する短いメッセージを出力してから終了します（注：Windowsの場合）
ただ使う -).
-+ or ただ +
任意の説明を含む広範な使用法（ヘルプ）メッセージを印刷します
モジュール固有のオプション（GMT共通オプションは除く）が終了します。
-? or いいえ 引数
オプションの説明を含む完全な使用法（ヘルプ）メッセージを印刷してから、
終了します。
- バージョン
GMTバージョンを印刷して終了します。
--show-datadir
GMT共有ディレクトリへのフルパスを出力して終了します。
GRID VALUES 精度
入力データの精度に関係なく、グリッド ファイルを作成する GMT プログラムは
内部ではグリッドを 4 バイト浮動小数点配列で保持します。 これはメモリを節約するために行われます
さらに、すべてではないにしてもほとんどの実データは 4 バイト浮動小数点を使用して保存できます。
価値観。 より高い精度のデータ (つまり、倍精度値) は、その精度を失います。
GMT がグリッド上で動作するか、新しいグリッドを書き出すと、精度が向上します。 損失を制限するには
データを処理するときは、精度を高める前にデータの正規化を常に考慮する必要があります。
処理。
例
hawaii_5x5.xyg の ASCII データから 5 分ごとに 5 分の重力ブロックをグリッド化するには、
張力係数 = 0.25、a convergence_limit = 0.1 ミリガル、結果をファイルに書き込みます
hawaii_grd.nc と呼ばれるファイルを作成し、各反復を監視して、次のことを試してください。
gmt サーフェス ハワイ_5x5.xyg -R198/208/18/25 -I5m -Ghawaii_grd.nc -T0.25 -C0.1 -Vl
onworks.net サービスを使用してオンラインで surfacegmt を使用する
