これは、Ubuntu Online、Fedora Online、Windows オンライン エミュレーター、MAC OS オンライン エミュレーターなどの複数の無料オンライン ワークステーションの XNUMX つを使用して、OnWorks 無料ホスティング プロバイダーで実行できるコマンド grdgradientgmt です。
プログラム:
NAME
grdgradient - グリッドから方向導関数または勾配を計算します。
SYNOPSIS
勾配勾配 in_grdfile out_grdfile [ アジム[/アジム2]] [[[a][c][o][n]] [
[s|p]アジム/エレブ[/アンビエント/拡散する/スペキュラ/輝く]] [ フラグ ] [[e][t][amp] [/シグマ[/オフセット]]
] [ 地域 ] [ スロープファイル ] [[レベル]] [ -fg ] [ -n]
注意: オプションフラグと関連する引数の間にスペースを入れることはできません。
DESCRIPTION
勾配勾配 指定された方向の方向導関数を計算するために使用できます (-A),
または方向 (-S) [そして大きさ(-D)] データのベクトル勾配。
出力の最初/最後の行/列の推定値は、境界条件によって異なります (「
-L).
REQUIRED 議論
in_grdfile
方向導関数の計算元となる 2 次元グリッド ファイル。 (「グリッド ファイル形式」を参照してください)
下)。
-Gout_grdfile
方向導関数の出力グリッド ファイルの名前。 (「グリッド ファイル形式」を参照してください)
下)。
オプション 議論
-Aアジム[/アジム2]
方向導関数の方位角方向。 アジム X、Y の角度です
北 (+y 方向) から時計回りに正の角度で測定された平面
東(+x方向)。 方向導関数の負の値、
-[dz/dx*sin(アジム) + dz/dy*cos(アジム)]、 見つかった; 否定は正の値を生成します
z(x,y) の傾きが下り坂の場合 アジム 方向、正しい感覚
画像の照明をシェーディングする (「 画像 グリッドビュー) 光源による
から輝く X、Y 平面の上 アジム 方向。 オプションで、XNUMX つ指定します
方位角、 -Aアジム/アジム2、この場合、これらの各方向の勾配は次のようになります。
計算され、大きさが大きい方が保持されます。 これは役に立ちます
XNUMX 方向の線状構造でデータを照明する。例: -A0/270
北(上)と西(左)から照らします。
-D[a][c][o][n]
データの正 (上り勾配) の勾配の方向を見つけます。 代わりに
アスペクト (下り坂の方向) を見つけます。 -ダ。 デフォルトでは、方向は次のとおりです。
北から時計回りに測定すると、 アジム in -A その上。 追加 c 従来のものを使用する
正の x (東) 方向から反時計回りに測定したデカルト角度。
追加 o 方向 (0 ~ 180) ではなく方向 (0 ~ 360) を報告します。 追加 n
すべての角度に 90 度を追加します (たとえば、サーフェスに局所的なストライクを与えるため)。
-E[s|p]アジム/エレブ[/アンビエント/拡散する/スペキュラ/輝く]
で使用するのに適切なランバーシアン放射輝度を計算します。 画像 グリッドビューを選択します。
ランバート反射は、すべての光を反射する理想的な表面を想定しています。
衝撃が当たると、表面はどの方向から見ても同じように明るく見えます。 アジム
標高 は光ベクトルの方位角と仰角です。 オプションで、供給 アンビエント
拡散する スペキュラ 輝く 反射率のプロパティを制御するパラメータです
表面の。 デフォルト値は次のとおりです。 0.55/0.6/0.4/10 一部の値を残すには
そのままの場合は、新しい値として = を指定します。 例えば -E60/30/=/0.5 を設定します アジム 標高
拡散する 60、30、0.5 に変更し、他の反射率パラメータはそのままにします
手付かずの。 追加 s より単純なランバーシアン アルゴリズムを使用します。 このフォームでは、
指定する必要があるのは、方位角と仰角のパラメータだけです。 追加 p 使用します
Peucker 区分線形近似 (より単純だが高速なアルゴリズム。この場合)
アジム 標高 315 度および 45 度に配線されています。 これは、たとえあなたが
他の値を指定しても無視されます。)
-Lフラグ 境界条件 フラグ 多分 x or y or xy データが範囲内で周期的であることを示します
x または y またはその両方、または フラグ 多分 g 地理的条件を示します (x と y は
経度と緯度)。 [デフォルトでは「自然な」条件が使用されます(二次偏微分法線)
エッジまではゼロです)。
-ネット][amp] [/シグマ[/オフセット]]
正規化。 [デフォルト: 正規化なし] 実際の勾配 g オフセットされており、
正規化された勾配を生成するためにスケーリングされる gn 最大出力の大きさは amp.
If amp 指定されていない、デフォルト amp = 1. の場合 オフセット は与えられず、に設定されます
の平均 g. -N 収量 gn = amp *(g - オフセット)/最大(腹筋(g - オフセット))。 -ネ
累積ラプラス分布を使用して正規化すると、次の結果が得られます。 gn = amp * (1.0 -
exp(平方根(2) * (g - オフセット)/ シグマ)) どこ シグマ の L1 ノルムを使用して推定されます。
(g - オフセット)与えられない場合。 -Nt 累積コーシーを使用して正規化します
分配利回り gn = (2 * amp /ピ)*アタン((g - オフセット)/ シグマ) どこ シグマ
の L2 ノルムを使用して推定されます (g - オフセット)与えられない場合。
-NS[単位]xmin/xmax/イミン/ワイマックス[r] (もっと ...)
関心のある領域を指定します。 の使用 -R オプションはサブセクションを選択します
in_grdfile グリッド。 このサブセクションがグリッドの境界を超える場合、
共通領域が抽出されます。
-Sスロープファイル
勾配ベクトルのスカラーの大きさを含む出力グリッド ファイルの名前。 必要 -D
しかし、作ります -G 任意。
-V [レベル] (もっと ...)
詳細レベル[c]を選択します。
-fg 地理グリッド(経度、緯度の寸法)はメートルに変換されます
現在の楕円体パラメータを使用した「フラットアース」近似を介して。
-n[b|c|l|n][+a][+bBC][+c][+tしきい値] (もっと ...)
グリッドの補間モードを選択します。
-^ or ただ -
コマンドの構文に関する短いメッセージを出力してから終了します(注:Windowsの場合)
ただ使う -).
-+ or ただ +
任意の説明を含む広範な使用法(ヘルプ)メッセージを印刷します
モジュール固有のオプション(GMT共通オプションは除く)が終了します。
-? or いいえ 引数
オプションの説明を含む完全な使用法(ヘルプ)メッセージを印刷してから、
終了します。
- バージョン
GMTバージョンを印刷して終了します。
--show-datadir
GMT共有ディレクトリへのフルパスを出力して終了します。
GRID 距離 単位
グリッドに水平単位としてメーターがない場合は、追加します +u単位 入力ファイルへ
指定された単位からメートルに変換する名前。 グリッドが地理的である場合は、変換します
供給によるメートルまでの距離 -fg を代わりにお使いください。
ヒント
何かわからない場合は -N 強度ファイルの作成に使用するオプション 画像 or
グリッドビュー、最初の試みとしては、 -ネ0.6.
通常、視覚化の目的には 255 階調で十分です。 ディスクを 75% 節約できます
出力ファイル名に =nb/a を追加してスペースを追加する out_grdfile.
大規模なデータ セットのサブ領域の照明付きマップをいくつか作成したい場合、
すべてのマップで照明効果を一貫させる必要がある場合は、 -N オプション
同じ値を指定します シグマ オフセット 〜へ 勾配勾配 マップごとに。 良い推測は次のとおりです
オフセット = 0と シグマ によって発見された グリッド情報 -L2 or -L1 正規化されていない勾配グリッドに適用されます。
単に必要な場合は、 x- または y-グリッドの派生、使用 算数.
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次元のグリッド。
例
geoid.nc のデータを、正規化された勾配を使用して照明するためのファイルを作成するには、
北と西の方向の光源を模倣した範囲 [-0.6,0.6]:
gmt grdgradient geoid.nc -A0/270 -Ggradients.nc=nb/a -Ne0.6 -V
ファイル topo.nc 内の海底構造の方位角を見つけるには、次のようにします。
gmt grdgradient topo.nc -Dno -Gazimuths.nc -V
参考文献
Horn、BKP、Hill-Shading and the Reflectance Map、IEEE Proceedings、Vol. 69、いいえ。
1、1981 年 14 月、47 ~ XNUMX ページ。 (http://people.csail.mit.edu/bkph/papers/Hill-Shading.pdf)
onworks.net サービスを使用してオンラインで grdgradientgmt を使用する