これは、Ubuntu Online、Fedora Online、Windows オンライン エミュレーター、または MAC OS オンライン エミュレーターなどの複数の無料オンライン ワークステーションの XNUMX つを使用して、OnWorks 無料ホスティング プロバイダーで実行できるコマンド grdhisteqgmt です。
プログラム:
NAME
grdhisteq - グリッドのヒストグラム等化を実行します。
SYNOPSIS
グルディステク in_grdfile [ out_grdfile ] [ n_cells ] [[file]] [[[ノルム] ] [ ] 地域 [レベル]
注: オプションフラグと関連する引数の間にスペースを入れることはできません。
DESCRIPTION
グルディステク ユーザーは、特定のグリッド ファイルを分割するデータ値を見つけることができます。
等しい面積のパッチ。 一般的な使い方の XNUMX つは、 グルディステク 一種のヒストグラム均等化です
イメージの。 このアプリケーションでは、ユーザーは、
真ん中の山。 このファイルの通常のグレー シェーディング (grdimage/grdview を使用)
トポグラフィーからグレートーンへの線形マッピングでは、ほとんどの画像が非常に複雑になります。
濃い灰色で、山はほぼ白です。 使えるかもしれない グルディステク に書き込みます
データ範囲を分割するデータ値の ASCII リストを stdout またはファイルに保存します。
n_cells セグメント。それぞれのセグメントは画像内で等しい面積を持ちます。 使用する awk or メイクセプト XNUMXつ
この出力を取得して CPT ファイルを構築できます。 grdimage で CPT ファイルを使用すると、次のようになります。
すべてのグレーレベルが均等に発生する画像。 あるいは、次を参照してください。 grd2cpt.
XNUMX番目の一般的な使用法は、 グルディステク ある種の統計に基づいたグリッドを作成中です
累積分布関数の。 このアプリケーションでは、出力は比較的高くなります。
入力ファイルと同じ (x,y) 位置にありますが、値は次のように変更されます。
何らかの累積分布におけるそれらの位置を反映します。 一例として、
データの下位 10%: グリッドを取得し、実行します グルディステク を使用してグリッドを作成します n_cells = 10
次に、結果を輪郭化して 1 の輪郭をトレースします。 これには、下位 10% が含まれます。
元の値に関係なく、データ。 別の例は、次の出力のイコライズです。
勾配勾配。 シェーディングの目的では、データが滑らかに分布していることが望まれます。
ガウシアンなど。 走れば グルディステク からの出力で 勾配勾配 そしてグリッドファイルを作成します
ガウス オプションを指定して出力すると、値が分散されたグリッドが得られます。
ゼロ平均と単位分散のガウス分布に従います。 の場所
これらの値は入力の位置に対応します。 つまり最もネガティブな
出力値は、最も負の入力値の (x,y) 位置になります。
REQUIRED 議論
in_grdfile
イコライズ対象の 2D グリッド ファイル。 (以下の「グリッド ファイル形式」を参照してください)。
オプション 議論
-Cn_cells
データ範囲を何セル(または分割)にするかを設定します[16]。
-D レベル情報をダンプします file、ファイルが指定されていない場合は標準出力。
-Gout_grdfile
出力 2-D グリッド ファイルの名前。 と一緒に使用 -N のみ。 (以下の「グリッド ファイル形式」を参照してください)。
-NS[ノルム]
ガウス出力。 と一緒に使用する -G 標準の正規スコアを使用して出力グリッドを作成します。
追加 ノルム スコアを強制的に <-1,+1> の範囲に収めます [デフォルトは標準です
通常のスコア]。
-Q 二次強度スケーリングを使用します。 [デフォルトは線形]。
-NS[単位]xmin/xmax/イミン/ワイマックス[r] (もっと ...)
関心のある領域を指定します。 の使用 -R オプションはサブセクションを選択します
in_grdfile グリッド。 このサブセクションがグリッドの境界を超える場合、
共通領域が抽出されます。
-V [レベル] (もっと ...)
詳細レベル[c]を選択します。
-^ 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次元のグリッド。
例
ファイル heights.nc を 16 等分する高さの間隔を見つけるには
エリア:
gmt grdhisteq heights.nc -C16 -D > レベル.d
ファイル raw_intens.nc 内の不十分に分散された強度を、
画像 or グリッドビュー走る
gmt grdhisteq raw_intens.nc -Gsmooth_intens.nc -N -V
制限事項
あなたが使用している場合 グルディステク グラデーション シェーディングのガウス出力を作成するには 画像 or
グリッドビュー、次の点に注意してください: 出力は [-x, x] の範囲になります。
ここで、x は入力グリッド内のデータ数 (nx * ny) と累積値に基づきます。
ガウス分布関数 F(x)。 つまり、N = nx * ny とします。 次に、x が調整されるので、
F(x) = (N - 1 + 0.5)/N です。 標準的な正常値からの値の約 68% であるため、
分布が +/- 1 以内に収まる場合、これは出力グリッドにも当てはまります。 しかし、もし N が非常に
大きい場合、x が 4 より大きくなる可能性があります。 グリッドビュー プログラム
グラデーションを [-1, 1] の範囲にクリッピングすると、データの 68% の正しいシェーディングが得られます。
一方、それらの 16% は -1 にクリップされ、それらの 16% は +1 にクリップされます。 これもできれば
画像の大部分が明るすぎるか暗すぎる場合は、出力を取得する必要があります。 グルディステク および
を使用して再スケールします 算数 範囲を縮小するには、1.0 未満の値を乗算します。
したがって、画像の 68% 以上が [-1, 1] の範囲内に収まります。
あるいは、次のように正規化係数を指定します。 -N.
onworks.net サービスを使用してオンラインで grdhisteqgmt を使用する