mapprojectgmt-クラウドでのオンライン

プログラム：

NAME

mapproject-順方向および逆方向のマップ変換、データム変換、および測地学を実行します

SYNOPSIS

地図プロジェクト [ テーブル ] パラメータ 地域 [ b|B|f|F|o|O[ロン0/lat0]] [[[dx/dy]] [ c|i|p ] [
[データム]] [[[単位]] [[[x0/y0/] [[+|-]単位][+|-]] [] [ ライン.xy[/ [+|-]単位][+]] [
[a|c|g|m]] [[[d|e ] [] [[[h]から[/〜へ]] [[[レベル]] [[[w|h]] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o]]

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

DESCRIPTION

地図プロジェクト から（経度、緯度）の位置を読み取ります テーブル [または標準入力]および
指定された地図投影法と縮尺を使用して（x、y）座標を計算します。 オプションで、
（x、y）位置を読み取り、逆を実行して（経度、緯度）値を計算できます
変身。 これは、デジタル化によって得られた線形（x、y）ポイントを変換するために使用できます
地理座標への既知の投影のマップ。 に沿った距離も計算できます
トラック、固定点、またはラインへの最も近いアプローチ。 最後に、実行するために使用することができます
さまざまなデータム変換。 最初の2列の後に追加のデータフィールドが許可されます
（経度、緯度）または（x、y）が必要です。 オプションを参照してください -: 読み方について
（緯度、経度）ファイル。

REQUIRED 議論

-Jパラメータ （もっと ...）
地図投影法を選択します。

-NS[単位]xmin/xmax/イミン/ワイマックス[r] （もっと ...）
関心領域を指定します。 UTMプロジェクションの特殊なケース： -C 使用されている
& -R が指定されていない場合、領域は指定されたUTMゾーンと一致するように設定されます。
完全な楕円体解を保持するために（詳細については、制限を参照してください）
情報）。

オプション 議論

テーブル XNUMXつ以上のASCII（またはバイナリ。を参照） -bi[ncols][type]）を保持するデータテーブルファイル
データ列の数。 テーブルが指定されていない場合は、標準入力から読み取ります。

-Ab | B | f | F | o | O [ロン0/lat0]
-Af 固定小数点から（前方）方位角を計算します lon / lat 各データポイントに。
　 -アブ データポイントから固定小数点に戻る方位角を取得します。 使用する -青 取得するため
方位角（90/90）ではなく方向（-0/360）。 大文字 F, B or O 意志
測地線から地球中心の緯度に変換し、測地線の方位角を推定します
（現在の楕円体が球ではないと仮定します）。 固定小数点が指定されていない場合は、
前のポイントから方位角（または後方方位角）を計算します。

-NS[dx/dy]
投影座標の中心を地図投影中心に設定します[デフォルトは
左下隅]。 必要に応じて、追加する投影単位にオフセットを追加します（または
減算されたとき -I falseなどの投影座標に（から）に設定されます
特定の投影ゾーン[0/0]の東向きと北向き。 に使用される単位
offsetsは、有効なプロット距離の単位です（PROJ_LENGTH_UNITを参照）。 -F is
使用されます。この場合、オフセットはメートル単位です。

-Dc | i | p
PROJ_LENGTH_UNITを一時的にオーバーライドして、 c （cm）、 i （インチ）、または p （ポイント）
代わりは。 と一緒に使用することはできません -F.

-E [データム]
測地（lon、lat、height）からEarth Centered Earth Fixed（ECEF）に変換する
（x、y、z）座標（追加 -I 逆変換の場合）。 データムIDを追加します（を参照） -Qd)
または与える 楕円:dx,dy,dz コラボレー 楕円 楕円体IDの場合があります（を参照） -量的緩和）または
として与えられる a[、* inv_f *]、ここで a は準主軸であり、 inv_f 逆です
平坦化（省略した場合は0）。 もしも データム -または与えられていない場合は、WGS-84を想定しています。

-NS[単位]
1：1スケーリングを強制します。つまり、出力（または入力。を参照） -I）データは実際に予測されています
メートル。 他の単位を指定するには、目的の単位を追加します（UNITSを参照）。 それなし -F,
出力（または入力、を参照） -I）はPROJ_LENGTH_UNITで指定された単位にあります（ただし
-D).

-NS[x0/y0/] [[+ |-]単位] [+ |-]
トラックに沿った距離を計算する or オプションのポイントセットに -Gx0 / y0。 追加する
距離単位（UNITSを参照）。 c （入力を使用したデカルト距離
座標）または C （投影座標を使用したデカルト距離）。 The C 単位
必要 -R & -J 設定します。 不動点が与えられていない場合、累積を計算します
トラックに沿った距離。 追加する - 連続する間の増分距離を取得するには
ポイント。 追加する + 入力ファイルの2つの追加列を介してXNUMX番目のポイントを指定します。

-I 逆変換を実行します。つまり、（x、y）データから（経度、緯度）を取得します。

-Lライン.xy[/ [+ |-]単位] [+]
入力データポイントからで与えられた線までの最短距離を決定します
ASCIIマルチセグメントファイル ライン.xy。 の距離と座標
最も近いポイントは、XNUMXつの新しい列として出力に追加されます。 追加する
距離単位（UNITSを参照）、 c （入力座標を使用したデカルト距離）
or C （投影座標を使用したデカルト距離）。 The C ユニットが必要です -R &
-J 設定します。 最後に、追加します + 線分IDと小数部をレポートします
最も近いポイントのlon/latの代わりにポイント番号。

-N [a | c | g | m]
測地緯度から変換します（現在の楕円体を使用します。PROJ_ELLIPSOIDを参照してください）。
XNUMXつの異なる補助緯度のXNUMXつに（経度は影響を受けません）。 選ぶ
から authalic、 c非公式、 gエオセントリック、そして meridional緯度[天動説]。 使用する -I
補助緯度から測地緯度に変換します。

-Q [d | e すべての投影パラメータを一覧表示します。 データムのみを一覧表示するには、 -Qd。 リストするだけ
楕円体、使用 -量的緩和.

-S 地域外のポイントを抑制します。

-T [h]から[/〜へ]
データム間の座標変換 から & 〜へ 標準のMolodenskyを使用
変身。 使用する -NS 3番目の入力列の高さが楕円体より高い場合[デフォルト
高さ=0、つまり楕円体上にあると仮定します]。 データムIDを使用してデータムを指定します（を参照）
-Qd）または与える 楕円:dx,dy,dz コラボレー 楕円 楕円体IDの場合があります（を参照） -量的緩和）または
として与えられる a[、* inv_f *]、ここで a は準主軸であり、 inv_f 逆です
平坦化（省略した場合は0）。 もしも データム -または与えられていない場合は、WGS-84を想定しています。 -T 多分
と組み合わせて使用 -R -J 座標投影の前にデータムを変更するには
（追加 -I 逆投影の後にデータム変換を適用します）。 それを確認してください
PROJ_ELLIPSOID設定はあなたのケースに適しています。

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

-W [w | h]
マップの幅と高さを標準出力に出力します。 入力ファイルは読み込まれません。 のみ
幅または高さを出力し、追加します w or h、 それぞれ。 の単位
寸法は次の方法で変更できます -D.

-bi [ncols][NS] （もっと ...）
ネイティブバイナリ入力を選択します。 [デフォルトは2入力列]。

-bo [ncols][type] （もっと ...）
ネイティブバイナリ出力を選択します。 [デフォルトは入力と同じ]。

-d [i | o]データなし （もっと ...）
等しい入力列を置き換えます データなし NaNを使用して、出力で逆を実行します。

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

-g [a] x | y | d | X | Y | D | [コル] z [+ |-]ギャップ[う] （もっと ...）
データのギャップと改行を特定します。

-h [i | o] [n] [+ c] [+ d] [+ r発言] [+ rタイトル] （もっと ...）
ヘッダーレコードをスキップまたは生成します。

-iコルズ[l] [s階段] [oオフセット] [、...] （もっと ...）
入力列を選択します（0は最初の列です）。

-oコルズ[、...] （もっと ...）
出力列を選択します（0は最初の列です）。

-p [x | y | z]アジム/標高[/zレベル] [+ wロン0/lat0[/z0]] [+ vx0/y0] （もっと ...）
パースビューを選択します。

-s [コルズ] [a | r] （もっと ...）
NaNレコードの処理を設定します。

-：[i | o] （もっと ...）
入力および/または出力の1番目と2番目の列を交換します。

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

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

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

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

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

単位

マップ距離の単位については、 単位 d アーク度の場合、 m アーク分、および s アーク用
XNUMX番目、または e メーターの場合[デフォルト]、 f 足用、 k kmの場合、 M 法定マイルの場合、 n 航海用
マイル、そして u 米国の測量フィート用。 デフォルトでは、球形を使用してそのような距離を計算します
大円での近似。 プリペンド - 距離まで（または単位が距離ではない
与えられた）「地球平面説」の計算を実行する（より速いが正確ではない）または前に追加する + 〜へ
正確な測地線計算を実行します（低速ですが、より正確です）。

ASCII FORMAT 精度

数値データのASCII出力形式は、 gmt.conf
ファイル。 経度と緯度はFORMAT_GEO_OUTに従ってフォーマットされますが、その他は
値はFORMAT_FLOAT_OUTに従ってフォーマットされます。 有効なフォーマットは次のことができることに注意してください
出力の精度が低下し、下流でさまざまな問題が発生する可能性があります。 もしも
出力が十分な精度で書き込まれていないことがわかった場合は、バイナリへの切り替えを検討してください
出力（-ぼ 可能な場合）またはFORMAT_FLOAT_OUT設定を使用してより多くの小数を指定します。

例

メートル単位のUTM座標を地理的な場所に変換するには、ファイルutm.txtと
UTMゾーン（およびゾーンまたは半球）を知っている場合は、

gmt mapproject utm.txt -Ju + 11/1：1 -C -I -F

メルカトル図法で（経度、緯度）のファイルをcm単位の（x、y）位置に変換するには
0.5度あたりXNUMXcmの所定のスケールのグリッド、実行

gmt mapproject lonlatfile -R20 / 50/12/25 -Jm0.5c> xyfile

いくつかの2列、バイナリ、倍精度ファイルを（緯度、経度）で変換するには
横メルカトルグリッド（中央経度75W）のインチ単位の（x、y）位置に
縮尺=1：500000で、マップ領域の外に出るポイントを抑制し、実行します

gmtmapprojecttracks。*-R-80/ -70 / 20/40 -Jt-75 / 1：500000-：-S -Di -bo -bi2> tmfile.b

NAD27のファイルold.datの測地座標（lon、lat、height）を変換するには
CONUSデータム（Clarke-131楕円体を使用するデータムID 1866）からWGS 84まで、実行

gmt mapproject old.dat -Th131> new.dat

入力ファイルquakes.datの各ポイント間の最も近い距離（km）を計算するには
マルチセグメントASCIIファイルcoastline.xyで指定された線分は、次のように実行されます。

gmt mapproject quakes.dat -Lcoastline.xy / k> quake_dist.dat

制限事項

で設定された長方形の入力領域 -R 一般に、非長方形にマッピングされます
グリッド。 そうでもなければ -C が設定されている場合、このグリッドの左端の点はxvalue = 0.0であり、
最下部のポイントはyvalue=0.0になります。 したがって、地図をデジタル化する前に、極端な
を介して座標をマップ 地図プロジェクト 適切なスケールを使用して、（x、y）値を確認します
それらはにマップされます。 これらの値をデジタル化の設定時に使用して、
逆変換は正しく機能するか、または代わりに使用します awk スケーリングおよびシフトするには
変換前の（x、y）値。

一部の投影では、ユーザーが選択したにもかかわらず、球形のソリューションを使用できます。
楕円。 これは、ユーザーが -R 設定は、ドメインを超える領域を意味します
楕円体の級数展開が有効です。 条件は次のとおりです。（1）
ランベルト正角円錐図法（-JL）およびAlbers Equal-Area（-J.B.）球面解を使用します
マップの縮尺が1.0E7を超えたとき。 （2）横メルカトル図法（-JT）およびUTM（-JU）
で与えられた西または東の境界のいずれかで球形の解を使用します -R 以上です
中央子午線から10度、（3）カッシーニについても同じ（XNUMX）-JC）ただし、制限はあります
わずか4度。

楕円体 そして スフェロイド

GMTは、楕円体の数式が実装されていて、ユーザーが楕円体の数式を選択した場合、それらを使用します。
参照形状としての楕円体（PROJ_ELLIPSOIDを参照）。 ユーザーはいくつかのことに注意する必要があります
潜在的な落とし穴：（1）横メルカトル図法、アルバース図法、
ランベルト正角円錐図法では、マップされた領域が次の場合に楕円体式を使用します。
小さく、球面式に切り替えます（そして適切な補助を置き換えます
緯度）より大きな地図の場合。 楕円体の公式は次のように使用されます。（a）横方向
メルカトル図法：すべての点が中央子午線から10度以内にある場合、（b）円錐図法
縦方向の範囲が90度未満の場合、（c）すべての点でのカッシーニ図法
中央子午線から4度以内です。 （2）いくつかの歴史的なものと一致させようとしているとき
データ（たとえば、特定の投影法と特定の参照で取得された座標
楕円体）GMTではわずかに異なる結果が得られる場合があります。 可能性が高いXNUMXつ
この不一致の原因は、古い計算では有効数字が少ないことが多いことです。 為に
たとえば、Snyderの例では、Clarke 1866楕円体（彼によって次のように定義されています）を使用することがよくあります。
扁平率f=1 / 294.98）。 fから、離心率の0.00676862818乗はXNUMXになります。
（これはGMTが使用するものです）一方、Snyderは四捨五入して0.00676866を使用します。 この違いは
数十cmの不一致を与えます。 投影された座標を再現する必要がある場合
このわずかに異なる離心率では、独自の楕円体を指定する必要があります
Clarke 1866と同じパラメーターですが、f = 1/294.97861076です。 また、古いことに注意してください
データは異なるデータを参照する場合があります。ただし、どのデータが使用されたかがわからない場合は、
すべてのデータを共通のデータに変換すると、数十から数百の不一致が発生する可能性があります
メートル。 （3）最後に、PROJ_SCALE_FACTORには、一部のデフォルト値があることに注意してください。
投影されるため、生成された結果と一致させるために設定をオーバーライドする必要がある場合があります
他の設定で。

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