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

プログラム：

NAME

gmtselect - 複数の空間基準に基づいてデータ テーブルのサブセットを選択します

SYNOPSIS

gmtselect [ テーブル ] [ 最小面積[/min_level/最大レベル][+ ag|i|s|S][+r|l][pパーセント]] [
DIST[単位]/ptファイル ] [ 分解能[+]] [[[fn]] [ ポリゴンファイル ] [ -I[cflrsz]] [ パラメータ
] [[p]DIST[単位]/ラインファイル ] [ マスク値 ] [ 地域 ] [ 分[/マックス][+cコル]] [[[レベル]] [
-b] [ -d] [ -f] [ -g] [ -h] [ -i] [
-o] [ -:[i|o]]

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

DESCRIPTION

gmtselect 最初の 2 列から位置 (経度、緯度) を読み取るフィルターです。
of インファイル [または標準入力] を使用し、1 ～ 6 の基準を組み合わせて合格または拒否します。
記録。 レコードは、1) 内部にあるかどうかに基づいて選択できます。
長方形領域 (-R [そして -J])、2) 以内 DIST の任意の地点からの km ptファイル、3）以内 DIST
任意の路線のキロメートル ラインファイル, 4) のポリゴンの XNUMX つの内側 ポリゴンファイル、5）内部
地理的フィーチャ (海岸線に基づく)、または 6) 指定された範囲内の Z 値を持っています。 の
これら 6 つの基準ごとに、テストの意味を逆転させることができます。 -I オプションを選択します。
オプションを参照 -: (緯度、経度) ファイルの読み取り方法について。

REQUIRED 議論

なし

オプション 議論

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

-A最小面積[/min_level/最大レベル] [+ ag | i | s | S] [+ r | l] [+ pパーセント]
面積が小さいフィーチャー 最小面積 km ^ 2または階層レベルで
より低い min_level 以上 最大レベル プロットされません[デフォルトは
0/0/4（すべての機能）]。 レベル2（湖）には通常の湖と広い川が含まれています
私たちが通常湖として含む体; 追加 +r ただ川の湖を手に入れるために +l
通常の湖を手に入れるだけです。 デフォルトでは（+ ai）棚氷の境界を次のように選択します
南極大陸の海岸線。 追加 + ag 代わりに氷の接地線を選択します
海岸線として。 独自の南極海岸線を印刷したいエキスパートユーザー向け
と島々 psxy あなたは使うことができます + as 60S未満のすべてのGSHHG機能をスキップするまたは + aS 〜へ
代わりに、60Sより北のすべての機能をスキップしてください。 最後に、追加します +pパーセント 除外する
対応するフル解像度機能の面積の割合が少ないポリゴン
より パーセント。 詳細については、以下の GSHHG 情報を参照してください。 そうでない場合は無視されます -N is
設定します。

-CDIST[単位]/ptファイル
位置が以下の範囲内にあるすべてのレコードを渡します DIST ASCII 内の任意のポイントの
file ptファイル。 場合 DIST がゼロの場合、 の 3 番目の列 ptファイル 各ポイントが必要です
個々の影響範囲。 距離はデカルト座標でユーザー単位で表されます。 特定
-fg 球面の距離を示し、距離単位を追加します (単位を参照)。
または、 -R & -J が使用され、地理座標が投影されます。
前のマップ座標 (PROJ_LENGTH_UNIT によって決定される、cm、インチ、またはポイント単位)
デカルト距離は次のように比較されます。 DIST.

-D分解能[+]
そうでない場合は無視されます -N が設定されています。 使用する海岸線データセットの解像度を選択します
((f）ull、（h）igh、（i）中級、（l)わあ、または(c）失礼）。 解像度が落ちる
データセット間で最大 80%。 [デフォルトは l]。 (+) を追加すると、自動的に下位が選択されます。
要求された解決策が利用できない場合の解決策 [見つからない場合は中止]。 注記
海岸線は細部が異なるため、ポイントが一致することは保証されません。
別の解像度が選択された場合、内側[または外側]に残ります。

-E[fn] ポリゴン境界上の点をどのように正確に考慮するかを指定します。 デフォルトでは、
そのような点は多角形の内側にあるとみなされます。 追加 n および f 変更するには
この動作は、 -F & -N 境界点が
外側にあるものとみなされます。

-Fポリゴンファイル
位置が閉じた多角形のいずれか内にあるすべてのレコードを渡します。
複数セグメントのファイル ポリゴンファイル。 球状ポリゴン (経度、緯度) の場合は、
連続する点は経度 180 度以上離れています。 ご了承ください
ポリゴンファイル かどうかに関係なく、ASCII で指定する必要があります。 -bi 使用されている。

-I[cflrsz]
指定された各基準のテストの意味を逆にします。

c どの点の影響範囲内にもないレコードを選択してください。

f どのポリゴンの内側にもないレコードを選択してください。

l どの行からも指定された距離内にないレコードを選択します。

r 指定された長方形領域内にないレコードを選択します。

s によって指定された内部とはみなされないレコードを選択します -N （と -A, -D).

z で指定された範囲内にないレコードを選択します -Z.

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

-L[p]DIST[単位]/ラインファイル
位置が以下の範囲内にあるすべてのレコードを渡します DIST のいずれかの線分の
ASCII 複数セグメント ファイル ラインファイル。 場合 DIST がゼロの場合、それぞれをスキャンします
のサブヘッダー ptファイル 埋め込み型の場合 -DDIST 各行を設定する設定
個々の距離値。 距離はデカルト座標でユーザー単位で表されます。 特定 -fg
球面距離を示すには、距離単位を追加します (「単位」を参照)。 あるいは、
if -R & -J 地理座標が使用され、地図座標に投影されます。
デカルト以前の (PROJ_LENGTH_UNIT によって決定される、cm、インチ、m、またはポイント単位)
距離を比較する DIST。 使用 -LP 直交する点のみを保証するため
最も近い線分への投影は、セグメントの端点内に収まります
[デフォルトでは、線の終点を「超えた」点が考慮されます。

-Nマスク値
指定された地理的フィーチャ内に位置するすべてのレコードを渡します。 特定
1 つの形式のうち 2 つを使用してレコードをスキップするか (s)、保持するか (k) するかどうか:

-N乾湿.

-N海/陸/湖/島/池.

[デフォルトは s/k/s/k/s (すなわち s/k) で、乾燥した土地上のすべてのポイントを通過します]。

-NS[単位]xmin/xmax/イミン/ワイマックス[r] （もっと ...）
関心のある領域を指定します。 地図投影法が指定されていない場合は、暗黙的に設定されます。
-Jx1.

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

-Z分[/マックス][+cコル]
3 列目が (z; コル = 2) 指定された範囲内にあるか、または
NaN (使用 -s NaN レコードをスキップします)。 もし マックス が省略されている場合は、次のことをテストします。 z 等しい 分
その代わり。 入力ファイルには少なくとも XNUMX つの列が必要です。 最小値に制限がないことを示すには
または max の場合はハイフン (-) を指定します。 3 番目の列が絶対時間の場合は、次のことを忘れないでください。
供給 -f2T。 別の列を指定するには、追加します +cコル、および複数のテストを指定するには
ただ繰り返すだけです Z テストする列がある場合は、このオプションを何度でも使用できます。 注: いつ
複数の Z オプションが与えられると、 Iz オプションは使用できません。

-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は最初の列です）。

-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設定を使用してより多くの小数を指定します。

この注意は、バイナリまたは netCDF 入力と組み合わせた ASCII 出力にのみ適用されます。
-: オプション。 以下の注記も参照してください。

注意 ON 加工 ASCII 入力 記録

を使用していない限り、 -: オプションを使用すると、選択した ASCII 入力レコードがそのままコピーされます。
出力。 つまり、次のようなオプションが -FOT FORMAT_FLOAT_OUT や
FORMAT_GEO_OUT は出力には影響しません。 一方、それは許可します
引用符の有無にかかわらず、文字列を含む多様な内容を含むレコードを選択する
コメントやその他の数値以外のコンテンツ。

注意 ON 距離

オプションの場合 -C or -L が選択されている場合、距離はデカルトおよびユーザー単位で表示されます。 使用 -fg
km 単位の球面距離と地理 (経度、緯度) 座標を意味します。 あるいは、
指定する -R & -J 投影されたデカルト距離をマップ単位 (cm、インチ、または
PROJ_LENGTH_UNIT によって決定されるポイント)。

このプログラムは長年にわたって進化してきました。 もともと、 -R & -J 順番的には必須でした
地理データを処理するためのものでしたが、球面計算が完全にサポートされるようになりました。 したがって、
-J テストを投影データに適用したい場合にのみ使用してください。
オリジナルコーディネート。 もし -J で与えられた距離が使用されます -C & -L 投影される
距離。

注意 ON セグメント

入力ファイルのセグメント ヘッダーは、ファイルからの XNUMX つ以上のレコードが出力にコピーされます。
セグメントはテストに合格しました。 選択はセグメントごとではなく、常にポイントごとに行われます。 それ
これは、テストに合格したセグメントのポイントのみが出力に含まれることを意味します。 もし
線をクリップし、セグメントの端に新しい境界点を含めたい場合は、
代わりに gmtspatial を使用する必要があります。

例

pts.d 内のいずれかのポイントから 300 km 以内にあるデータセットのサブセットを抽出するには、
line.d の線から 100 km 以上離れた場所で、実行します。

gmt select lonlatfile -fg -C300k/pts.d -L100/lines.d -Il > サブセット

ここで指定する必要があります -fg これにより、プログラムは地理データを処理していることを認識します。

data.d 内のすべてのポイントを、陸地上のポイントを除き、指定された領域内に保持するには (
高解像度の海岸線によって決定されます)、使用します

gmt select data.d -R120/121/22/24 -Dh -Nk/s > サブセット

球状多角形 lonlatpath.d の内部または上にある quakes.d 内のすべての点を返すには、次のようにします。
試します

gmt select quakes.d -Flonlatpath.d -fg > subset1

原点 d の点から 5 cm 以内にある station.d 内のすべての点を返します。
特定の投影、試してみてください

gmt select station.d -C5/origin.d -R20/50/-10/20 -JM20c \
--PROJ_LENGTH_UNIT=cm > サブセット 2

GSHHS 情報

海岸線データベースはGSHHG（以前のGSHHS）であり、次のXNUMXつのソースからコンパイルされています。
World Vector Shorelines（WVS）、CIA World Data Bank II（WDBII）、および雪氷圏のアトラス
（AC、南極のみ）。 南極大陸を除いて、すべてのレベル1ポリゴン（海域）
境界）は、より正確なWVSから派生しますが、すべての高レベルのポリゴン（レベル
2-4、土地/湖、湖/湖の島、および
湖の島/湖の湖の境界）は、WDBIIから取得されます。 南極大陸
海岸線にはXNUMXつのフレーバーがあります。アイスフロントまたはグラウンディングラインで、 -A オプションを選択します。
WVS、WDBII、およびACデータを使用可能な形式に変換するために多くの処理が行われました。
GMT：線分から閉じたポリゴンを組み立て、重複をチェックし、
ポリゴン間の交差を修正します。 各ポリゴンの面積が決定されました
ユーザーが最小領域よりも小さいフィーチャを描画しないことを選択できるようにします（を参照）。 -A）; XNUMX
含まれるポリゴンの最上位の階層レベルを制限する場合もあります（4は
最大）。 4つの低解像度データベースは、フル解像度データベースから派生しました
Douglas-Peuckerライン簡略化アルゴリズムを使用します。 河川の分類と
境界線はWDBIIの境界線に従います。 GMTクックブックおよびテクニカルリファレンス付録Kを参照してください
詳細については、。

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