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

プログラム：

NAME

img2grd - メルカトル形式または地理形式で img ファイルのサブセットを抽出します

SYNOPSIS

img2grd 画像ファイル グリッドファイル 地域 type [] [[ミンラット/マックスラット]] [] [ 分 ] [] [ ナビゲーション
] [[階段]] [[[レベル]] [ マックスロン ] [ -n]

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

DESCRIPTION

img2grd img 形式のファイルを読み取り、サブセットを抽出して、グリッド ファイルに書き込みます。 の -M
オプションは、img ファイルの球面メルカトル図法を使用するかどうかを指定します。
保存するか、メルカトル図法を元に戻して地理グリッドを書き込む必要があるかどうか。 もしも
地理的グリッドが選択されている場合は、正確な地理的グリッドへのリサンプリングをリクエストすることもできます。 -R 与えられた。

REQUIRED 議論

画像ファイル
海洋重力フィールドや海底地形フィールドなどのメルカトル img 形式ファイル
Sandwell と Smith による衛星高度計データから推定。 ユーザーが設定した場合
環境変数 $ GMT_DATADIRをタップし、その後、 img2grd 見つけようとします 画像ファイル in
$ GMT_DATADIR; それ以外の場合は開こうとします 画像ファイル 直接。

-Gグリッドファイル
グリッドファイル 出力グリッド ファイルの名前です。

-NS[単位]西/東/南/北[/ズミン/zmax][NS]
西, 東, 南, 北 関心領域を指定すると、指定できます
それらは小数度または[+-] dd：mm [：ss.xxx] [W | E | S | N]形式です。 追加する r 低い場合
w / e / s / nの代わりに、左右のマップ​​座標が表示されます。 二つ
速記 -Rg & -Rd グローバルドメインの略（経度で0/360および-180 / + 180
それぞれ、緯度で-90 / + 90）。 または、グリッドを作成する場合は、
Rコデロン/ラット/nx/nyここで、 コード L、C、Rの2文字の組み合わせです（左の場合、
中央、または右）およびT、M、B（上、中央、または下）。 たとえば、左下のBL。
これは、長方形領域のどのポイントが LON/ラット 座標参照
to、およびグリッドの寸法 nx & ny を介してグリッド間隔で -I 作成に使用されます
対応する地域。 または、既存のグリッドファイルの名前を指定します
と -R 設定（および該当する場合はグリッド間隔）がグリッドからコピーされます。
使い方 -R単位 選択したものと互換性のある投影（デカルト）座標を期待します -J
逆に投影して、実際の長方形の地理的領域を決定します。 にとって
透視図 （-p）、オプションで/を追加ズミン/zmax。 透視図の場合
(-p）、z範囲（ズミン, zmax）を追加して、XNUMX次元を示すことができます。 この
を使用する場合にのみ実行する必要があります -Jz オプション、のみを使用する場合ではありません -p オプションを選択します。
後者の場合、平面の斜視図がプロットされます。
寸法。

オプション 議論

-C 投影中心を基準とした x および y メルカトル座標を設定します [デフォルトは
グリッドの左下隅を基準に]。 必要 -M.

-NS[ミンラット/マックスラット]
拡張緯度範囲 -80.738/+80.738 を使用します。 または、追加する
ミンラット/マックスラット 入力画像ファイルの緯度範囲として。 [デフォルトは
-72.006/72.006]。 程度を判断できるため、通常は必要ありません。
ファイルサイズの検査。

-E こんなときに使えます -M 最終的なグリッドにまったく同じものを強制するように設定されていません
でリクエストされた地域 -R。 デフォルトでは、最終領域は次の直接投影です。
元のメルカトル領域からわずかに超えて広がるのが通常です。
要求された緯度範囲、さらに緯度のグリッド増分は行われません
経度の増分と一致します。 ただし、追加のリサンプリングにより小さな問題が発生します。
内挿エラーが発生するため、出力グリッドが内挿エラーと一致する必要がある場合にのみ使用してください。
要求された領域であり、x_inc = y_inc になります。 この場合、によって設定された領域は、 -R は
増分の倍数で与えられます (例: -R0/45/45/72)。

-I 示す 分 経度分単位の入力画像ピクセルの幅として。
[デフォルトは 2.0]。 ピクセル サイズは次から決定できるため、通常は必要ありません。
サイズの検査。

-M 球面メルカトル グリッドを出力します [デフォルトは地理的な経度/緯度グリッドです]。 の
img ファイルの球面メルカトル図法は保持されるため、領域 -R
ユーザーが設定した値はわずかに変更されます。 変更された領域はエッジに対応します
ピクセルのグループ[またはピクセルのグループ ナビゲーション ピクセル]。 グリッド ファイルのヘッダーは、x と
y 軸の長さは画像の西端と南端からの距離を表します。
ユーザーのデフォルト単位で測定され、 -Jm1と調整された -R。 を設定することにより
デフォルト プロジェクト_ 楕円 = Sphere、ユーザーは調整されたオーバーレイを作成できます -R
一致するように。 見る 例 下。 調整された -R グリッドにも書かれています
ヘッダーのコメントにあるので、後で見つけることができます。 見る -C に相対的な座標を設定するには
プロジェクションセンター。

-Nナビゲーション 入力画像ピクセルの値を平均します。 ナビゲーション by ナビゲーション 正方形を作成し、
このような正方形ごとに XNUMX つの出力ピクセル。 と一緒に使用する場合 -T3 平均を報告します
0 と 1 の間の制約。 -T2 出力は平均データ値になります
平均制約が 0.5 より大きいかどうかに応じて、NaN または NaN になります。 ナビゲーション 均等に分けなければなりません
imgfile のサイズをピクセル単位で変換します。 [デフォルト 1 平均化は行いません]。

-NS[階段]
img ファイルの値を乗算します。 階段 グリッドファイルに保存する前に。 [デフォルトは
1.0]。 最近の img ファイルの場合: img topo ファイルは (修正された) メートルに保存されます [-S1];
mGal*10 の自由空気重力ファイル [-SmGal を取得するには 0.1]; 垂直偏向ファイル
マイクロラジアン*10 [-Sマイクロラジアンを取得するには 0.1]、垂直重力勾配ファイル
エトヴォス*50 [-SEotvos を取得するには 0.02、または -SmGal/km を取得するには 0.002])。 いいえの場合 階段 与えられる
ファイル名を手がかりに調べてスケールを決定しようとします。

-Ttype type 制約情報のエンコードを処理します。 type = 0 は、いいえを示します
このような情報は img ファイル内にエンコードされます (1995 年より前のバージョンの
重力データ)、すべてのデータを取得します。 type > 0 は、制約情報が
エンコードされたもの (1995 年以降 (現在の) バージョンの img ファイル)
次のようにグリッド ファイルを作成します。 -T1 すべてのポイントでデータ値を取得し、 -T2 データを取得します
制約された点の値と補間点の NaN。 -T3 で 1 を取得します
拘束されたポイントでは 0、補間ポイントでは 1 [デフォルトは XNUMX]。

-V [レベル] （もっと ...）
詳細レベル [c] を選択します。 ここは特に参考になるのでおすすめです
座標がどのように調整されるか。

-Wマックスロン
示す マックスロン 入力画像ファイルの最大経度範囲として。 バージョン
1995年以来 マックスロン = 360.0、以前の一部のファイルは マックスロン = 390.0。
[デフォルトは360.0]。

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

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

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

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

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

地理的 例

　 -M 出力グリッドを地理的に表示する必要がある場合は、このオプションを除外する必要があります。
コーディネート。 リージョン内のデータを抽出するには -R-40/40/-70/-30から world_grav.img.7.2 &
再投影して地理座標を取得するには、試してみてください

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -V

img ファイル内の緯度の間隔はメルカトル単位で等距離であるため、
結果のグリッドは指定されたものと一致しません -R 緯度の間隔は正確に一致しません。
経度の間隔と等しくなります。 正確に一致する必要がある場合は、 -R そして同じ間隔で
経度と緯度では、 -E オプション：

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -E -V

メルカトル 例

img ファイルはメルカトル図法であるため、地理グリッドを抽出しないでください。
メルカトル図法を作成する計画がある場合。 それを行うと、投影が終了し、
グリッドを再投影し、短波長のディテールを失います。 使った方が良い -M そしてグリッドをプロットします
目的のメルカトル図法と同じスケールの線形投影法を使用します (GMT を参照)
例29）。 リージョン内のデータを抽出するには -R-40/40/-70/-30から world_grav.img.7.2走る

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav.nc -R-40/40/-70/-30 -V

なお、 -V オプションは、範囲が次のように調整されたことを示します。
-R-40/40/-70.0004681551/-29.9945810754。 私たちも使うことができます グリッド情報 グリッドファイルを見つけるには
ヘッダーはその領域が次であることを示します -R0/80/0/67.9666667 これは取得する x,y の範囲です。
球面メルカトル図法から -R-40/40/-70.0004681551/-29.9945810754 および
-Jm1. したがって、ship.lonlatgrav を取得し、それを使用して merc_grav.nc をサンプリングするには、次のようにします。

gmt セット PROJ_ELLIPSOID 球

gmt マッププロジェクト -R-40/40/-70.0004681551/-29.9945810754 -Jm1i ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | GMT マッププロジェクト \
-R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I > ship.lonlatgravsat

このような場合には、データを投影および投影解除する上記の方法を使用することをお勧めします。
アプリケーションの場合、補間ステップは XNUMX つしかないため ( グリッドトラック）。 XNUMX つあれば
まずグリッド ファイルを経度、緯度に変換してからサンプリングします。次の XNUMX つがあります。
補間ステップ (変換時およびサンプリング時)。

上記のグリッドから経度、緯度のグリッドを作成するには、使用できます

gmt grdproject merc_grav.nc -R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I -D2m -Ggrav.nc

場合によっては、これは簡単ではありません。 -R XNUMX つの座標系は一致しない可能性があります
良い。 このような場合は、 を使用することもできます (実際には、常に使用した方がよい場合があります)。

gmt grd2xyz merc_grav.nc | GMT マッププロジェクト \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -I | \
gmt 表面 -R-40/40/-70/70 -I2m -Ggrav.nc

上記の地域のメルカトル図を作成するには、次の gmt.conf 値を想定します。
PROJ_LENGTH_UNIT はインチです。 すると、上記の merc_grav.nc ファイルが投影されるので、 -Jm1i
幅は80インチです。 次を使用すると、幅 8 インチの地図を作成できます。 -Jxどの地図プログラムでも 0.1i
このグリッドに適用されます (例: グリッドコンター, 画像, グリッドビュー)、次にオーバーレイの場合
経度、緯度で作業します (例: psxy, コースト）上記の調整済みを使用できます -R & -Jm0.1 で取得
XNUMX つのシステムを一致させます。

しかし、私たちはこれよりも賢くすることができます。 入力画像ファイルのピクセル数が 2.0 であることに気づきました。
分幅（または grdinfo merc_grav.nc で nx と ny をチェック）すると、次のことがわかります。
merc_grav.nc は img ファイルのフル解像度を使用し、2400 x 2039 ピクセルを持っています。
幅 8 インチでは、300 インチあたり XNUMX ピクセルになります。 私たちはそれほど多くは必要ないと判断し、
100 インチあたり 3 ピクセルで満足できるため、データを 3 x XNUMX に平均化したいと考えます。
正方形。 (等高線プロットが必要な場合は、おそらくデータをより平均化することを選択するでしょう。
(例: 6 x 6)、滑らかな輪郭を取得します。) 2039 は 3 で割り切れないため、次の結果が得られます。
異なる調整済み -R この時：

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40/-70/-30 -N3 -V

今回は、調整された領域が次であることがわかります。 -R-40/40/-70.023256525/-29.9368261101 および
出力は 800 x 601 ピクセルで、私たちにとっては適切なサイズです。 これで人工物を作成できるようになりました
このために使用するイルミネーション ファイル 勾配勾配:

gmt grdgradient merc_grav_2.nc -Gillum.nc -A0/270 -Ne0.6

「grav.cpt」という CPT ファイルもあれば、カラー陰影起伏図を作成できます。
このような：

gmt grdimage merc_grav_2.nc -Iillum.nc -Cgrav.cpt -XXS0.1i -K > map.ps
gmt psbasemap -R-40/40/-70.023256525/-29.9368261101 -Jm0.1i -Ba10 -O >>map.ps

img ファイルから制約されたデータ値のみを緯度/経度で取得したいとします。
コーディネート。 それから実行します img2grd -T2 オプション、使用 grd2xyz 値をダンプするには、パイプします
grep -v NaN を使用して NaN を削除し、パイプスルーします。 地図プロジェクト 逆に
上記のような投影です。

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