GRDCONTOUR

名前
書式
説明
オプション
アスキー形式の精度

バグ
関連項目

名前

grdcontour − 2 次元のグリッドデータセットの等値線を作成する

書式

grdcontour grdfile −Ccont_int −Jparameters [ −A[-|annot_int][labelinfo] ] [ −B[p|s]parameters ] [ −Ddumpfile ] [ −Eazimuth/elevation ] [ −G[d|f|n|l|L|x|X]params ] [ −K ] [ −Llow/high ] [ −M[flag] ] [ −O ] [ −P ] [ −Qcut ] [ −Rwest/east/south/north[r] ] [ −Ssmoothfactor ] [ −T[+|-][gap/length][:LH] ] [ −U[/dx/dy/][label] ] [ −V ] [ −W[+][type]pen ] [ −X[a|c|r][x-shift[u]] ] [ −Y[a|c|r][y-shift[u]] ] [ −Z[factor[/shift]][p] ] [ −ccopies ] [ −bo[s|S|d|D][ncol] ]

説明

grdcontour は 2 次元のグリッドファイルを読んでグリッド全体の各等値線を たどることにより等値線マップを作成する。オプションとして 、 等 値 線 の x/y/z の位置を 1 つの複数セグメントのファイル又は複数に分割されたファイ ルにダンプすることができる。 PostScript コードが作成されて標準出力に 送 られる。描画に影響を及ぼす様々なオプションが利用できる。
オ プションフラグとそれに付属する引き数の間にスペースを入れないこと。オ プションフラグには大文字を使い、その付属コードには小文字を使うこと。

grdfile

等値線を作成する 2 次元のグリッドデータセット

−C

等値線の描き方は 3 通りの方法の中から 1 つを指定する。

cont_int に ".cpt" という拡張子が含まれている場合は、ファイルと し て開かれカラーパレットテーブルであると見なされる。その色の境 界が等値線レベルとして使われる。 cpt ファイルの最後の列に目盛フ ラ グが入っているとき、その等値線には目盛が入れられる。デフォル トでは全ての等値線にラベルが付けられる。全ての目盛を無効にす る には −A− を使うこと。
cont_int
が cpt ファイルではないファイルの場合は、 1 列目に等値 線レベル、 2 列目に C(ontour) 又は A(nnotate) が入っていると 見 な される。 C (又は c )でマークされたレベルには等値線が引かれ、 A (又は a )とマークされたレベルには目盛が入れられる。オプション として、 3 列目にその等値線レベルに対して固定された目盛の角度を 入れることもできる。
ファイルが見つからないときは、 cont_int は一定の等値線間隔で あ る と見なされる。 −A が設定されていて −C が設定されていないとき は、等値線間隔は目盛間隔と同じ値に設定される。
ファイルが与えられていて −T が設定されているときは、大文字 C 又 は A でマークされた等値線にだけケバが付けられる。全ての場合にお いて等値線はグリッドと同じ単位である。

−J

地図の図法を選択する。尺度は UNIT/degree 、1:xxxxx 、 又 は UNIT 単 位 の 幅 (オプションを大文字にする) である。 UNIT は .gmtdefaults4 で設定された MEASURE_UNIT によって cm 、インチ 、 又は m になるが、 c, i, m を scale/width の値に付けることによっ て、コマンドライン上で上書き変更することができる。地図の高さ 、 最 大 の大きさ、最小の大きさを指定するには、それぞれ h, +, - を width に付けること。

詳細は psbasemap の man ページを参照すること。

円筒図法:

−Jclon0/lat0/scale (カッシーニ図法)
−Jj
lon0/scale (ミラー図法)
−Jm
scale (メルカトル図法 - 投影中心はグリニッジ子午線と赤道)
−Jm
lon0/lat0/scale (メルカトル図法 - 中心子午線と標準緯度 を 設 定)
−Joa
lon0/lat0/azimuth/scale (斜めメルカトル図法 - 投影中心と方 位角を設定)
−Job
lon0/lat0/lon1/lat1/scale (斜めメルカトル図法 - 2 点を設定)
−Joc
lon0/lat0/lonp/latp/scale (斜めメルカトル図法 - 中心点と極 を設定)
−Jq
lon0/scale (等距円筒図法 (Plate Carree))
−Jt
lon0/scale (TM - 横メルカトル図法、赤道は y = 0)
−Jt
lon0/lat0/scale (TM - 横メルカトル図法、原点を設定)
−Ju
zone/scale (UTM - ユニバーサル横メルカトル図法)
−Jy
lon0/lats/scale (基本円筒図法)

方位図法:

−Jalon0/lat0/scale (ランベルト正積方位図法)
−Je
lon0/lat0/scale (正距方位図法)
−Jf
lon0/lat0/horizon/scale (心射図法)
−Jg
lon0/lat0/scale (正射図法)
−Js
lon0/lat0/[slat/]scale (一般平射図法)

円錐図法:

−Jblon0/lat0/lat1/lat2/scale (アルベルス正積円錐図法)
−Jd
lon0/lat0/lat1/lat2/scale (正距円錐図法)
−Jl
lon0/lat0/lat1/lat2/scale (ランベルト正角円錐図法)

様々な図法:

−Jhlon0/scale (ハンメル図法)
−Ji
lon0/scale (正弦曲線図法=サンソン図法)
−Jk
[f|s]lon0/scale (エッケルト IV (f) 及び VI (s) 図法)
−Jn
lon0/scale (ロビンソン図法)
−Jr
lon0/scale (ヴィンケル第 3 式図法)
−Jv
lon0/scale (ファン・デア・グリンテン図法)
−Jw
lon0/scale (モルワイデ図法)

非地理的図法:

−Jp[a]scale[/origin][r|z] (極座標 (theta,r))
−Jx
x-scale[d|l|ppow|t|T][/y-scale[d|l|ppow|t|T]] (線形、対数 、 指数スケーリング)

オプション

オプションフラグとそれに付属する引き数の間にスペースを入れないこと。

−A

annot_int はデータの単位の目盛間隔である。ファイルで等値線 レベルが与えられている場合は無視される [デフォルトでは目盛は 付 か ない]。 - を付けると −C で指定された全ての目盛が無効になる。 オプションの labelinfo はラベルの書式を制御するためのもので、制 御のための記述を結合した以下の文字列で構成される。

+aangle

目 盛の角度を固定する。 +an とすると等値線に垂直に、 +ap とすると等値線に平行[デフォルト]になる。 u 又は d を付け ると目盛の上端が隣接する上又は下の目盛の方向をそれぞれ向 く。

+cdx[/dy]

ラベルとオプションのテキストボックスの間の間隔を設定する 。 c|i|m|p を付けて単位を指定するか、 % を付けてラベルの フォントサイズ [15%] のパーセンテージを指定することが で きる。

+ffont

希 望 す る フ ォ ントを設定する[デフォルトは ANNOT_FONT_PRIMARY]。

+g[color]

不透明なテキストボックスを選ぶ[デフォルトは透明]。随意で 色を指定できる[デフォルトは PAGE_COLOR]。

+jjust

ラベルの位置を設定する[デフォルトは CM]。

+kcolor

テ キ ス ト ラ ベ ル の 色 を 設 定 す る[ デ フォルトは COLOR_BACKGROUND ]。

+o

角の丸い長方形のテキストボックスを選ぶ[デフォルト は 普 通の長方形]。曲がった文字列 (+v) 適用できないし、不透 明なテキストボックスでしか意味が無い。

+p[pen]

thin[ner|nest], テキストボックスの枠を描く[デフォルト は 枠無し]。随意で枠のペン属性を指定できる[デフォルトは幅 = 0.25p 、色 = 黒、テクスチャ = 実線]。 pen (ペン属性)はコ ン マで区切られた width, color, texture のリストで、それ ぞれはオプションである。 width (幅)は長さの単位(ポイント 、 セ ン チ メートル、インチ)又は faint, thin[ner|nest], thick[er|est], fat[ter|test], obese のいずれかで指定する ことができる。 color (色)はグレースケール (0−255) 又はカ ラー (r/g/b, それぞれ範囲は 0−255; h-s-v, 範囲は 0−360, 0−1, 0−1; c/m/y/k, それぞれ範囲は 0−100%; 有効なカラー 名) により指定できる。 texture (テクスチャ)はダ ッ シ ュ ‘-’ と点 ‘.’ の組合せである。

+rmin_rad

等値線の曲率半径が min_rad [デフォルトは 0 ]より小さな場 所にはラベルを配置しない。

+ssize

希望するフォントサイズをポイント単位で設定す る[ デフォルトは 9 ]。

+uunit

全 ての等値線のラベルに unit (単位)を付ける。 unit の頭にハイフン (-) が付けるとラベルの値と単位の間の 隙間がなくなる。 unit が付いていないときはグリッドファイ ルの中にある単位が用いられる [デフォルトでは単位は付かな い]。

+v

等値線に沿って曲がったラベルにする[デフォルトではま っすぐなラベルになる]。

+w

ラベルの角度を推定するのに使われる (x, y) の点数を指 定する[デフォルトは 10 ]。

+=prefix

prefix を全ての等値線のラベルの前に付ける。 prefix の頭 にハイフン (-) が付けるとラベルの値と単位の間の隙間が な くなる [デフォルトでは付かない]。

−B

地 図 の 枠 の 目 盛とケバの間隔を設定する。全ての詳細は psbasemap の man ページを参照すること。

−D

各等値線の (x,y,z) 座標を等値線のセグメントごとに別々のファ イ ルにダンプする。ファイルには dumpfile_cont_segment[_i].xyz ( 又は -b が選択されたときは .b )という名前が付く。ここで cont は 等値線の値で segment は各等値線に対するセグメント番号 (閉じた等 値線には _i が付く) である。プレフィックスが ’-’ で与えられた場 合 、ファイル名は C#_i (内部) 又は C#_e (外部) プラス拡張子にな り、 # はランニングナンバーになる。これにより DOS で動く GNU の ユ ーティリティでも利用できる短いファイル名にすることができる。 しかし −M−D と一緒に使うと代わりに 1 つのマルチセグメントの ファイルが作られる。

−E

遠近図法に使う視点の方位角と仰角を設定する[180/90]。

−G

等値線に沿ったラベルの配置を制御する。以下の 5 通りの制御ア ルゴリズムの中から選ぶこと。

−Gddist[c|i|m|p] 又は −GDdist[d|e|k|m|n]

小文字の d は描画上のラベルの間隔を c (cm), i (イン チ), m (メートル), p (ポイント)のいずれかの単位によって示すも のである。一方、大文字の De (m), k (km), m (マイル), n (海里), d (度)のいずれかの地図上の単位によって示すもの である[デフォルトは 10c 又は 4i ]。

−Gfffile.d

アスキーファイル ffile.d を読んで等値線に沿った位置に 合 うファイルの位置にラベルを配置する。不正確な位置や領域外 の点は飛ばされる。

−Gl|Lline1[,line2,...]

1 つ以上のコンマで区切られた直線セグメントの start ( 開 始) と stop (終了)の座標を与える。これらの線が等値線と交 差する位置にラベルは配置される。それぞれの line のフォー マ ッ ト は start/stop であり、 startstop は特定の点 lon/lat (経度/緯度)又は pstext で地図上の点を示すのに 使 わ れる [LCR][BMT] で与えられる 2 文字の XY のキーのいず れかである。さらに、 Z+ と Z- によりグリッドのそれぞれ最 大 、最小の位置を与えることもできる。 −GL は点の組を大円 として解釈する[デフォルトは直線]。

−Gnn_label

等値線に対する等距離なラベルの数を指定する[ 1 ]。大文 字 の −GN は線の開始点ちょうどからラベルの配置を開始する[デ フォルトでは線の中央に配置する]。 −GN-1 は等値線の開始点 に ラベルを配置し、一方 −GN+1 は等値線の終了点にラベルを 配置する。随意で /min_dist[c|i|m|p] を付けると連続するラ ベルの最小間隔が保たれる。

−Gx|Xxfile.d

複 数セグメントファイル xfile.d を読んで等値線と xfile.d の線の交差点にラベルを配置する。 −GX にすると最初に大 円 に沿って線をサンプルし直す。

−K

さらに PostScript コードが後に続くときに使用する [デフォル トでは描画システムを終了する]。

−L

範囲を制限する。 low より下や high より上のデータ値に対する 等値線は描かれない。

−M

−D と一緒に使うと 1 つのマルチセグメントのファイルが作られ る。第 1 列が flag でその後に等値線のレベルが書かれたヘッダ行の 後に各等値線のセクションが続く。

−O

上書き描画モードを選択する [デフォルトでは新しい描画システ ムに初期化する]。

−P

縦置き描画モードを選択する [ GMT のデフォルトは横置きである 。これを変更するには gmtdefaults を参照すること]。

−Q

cut の数より少ない点については等値線を描かない [全ての等値 線を描く]。

−R

xmin, xmax, ymin, ymax は対象とする領域を指定する。地理的領 域 に対しては、これらの境界は西、東、南、北に対応し、小数値を含 む度又は [+-]dd:mm[:ss.xxx][W|E|S|N] のフォーマットで指定できる 。東西南北の代わりに地図座標の左下と右上を使うときは r を付ける こと。 2 つの略号 −Rg −Rd は全球領域(それぞれ経度 0/360 又 は -180/+180 と緯度 -90/+90 )を表す。暦時間座標系に対しては相対時 間(選択された TIME_EPOCH に対して相対的で、選択された TIME_UNIT を単位とする。 t−JX|x に付けること)又は [date]T[clock] 形式 の絶対時刻(T−JX|x に付けること)のどちらかを与えることができ る。少なくとも dateclock のどちらかが存在していなければなら ない。 T は常に必要である。 date の文字列は [-]yyyy[-mm[-dd]] ( グレゴリオ暦)又は yyyy[-Www[-d]] (ISO 週暦)の形式でなければなら ない。 clock の文字列は hh:mm:ss[.xxx] の形式でなければならない 。 区切り文字を利用する場合にはその型と位置が指定どおりでなけれ ばならない(ただし入出力と描画フォーマットは変更可能である) [ デ フォルトはグリッドファイルで定義された領域]。

−S

等値線の線を大まかに (gridbox_size/smoothfactor) の間隔ごと にサンプルし直す。

−T

下向きのケバを最も内側の閉じた等値線に沿って各 gap ごとに描 く 。 gap とケバの長さを付けるか、デフォルト[ 0.5c/0.1c 又は 0.2i/0.04i ]を使うこと。 −T+ 又は −T- を指定するとそれぞれ局 所 的 な 高部又は低部にケバをつけることができる。 :LH を付けると L や H の文字が最も内側の閉じた等値線(局所的な高部や低部)の中心に 記 さ れる。 L と H は単一の文字ならば何にでもできる(例:LH、-+ 等)。ファイルが −C で与えられ −T が設定されている場合、大文字 C 又は A がマークされた等値線だけにケバ[及び目盛]が付けられる。

−U

図に Unix システムのタイムスタンプを入れる。図の左下隅に対 するスタンプの左下隅の相対的な位置を決めることができる。随意 で ラ ベ ル又は c (コマンドの文字列を入れる)を付けることができる。 GMT のパラメータ UNIX_TIMEUNIX_TIME_POS は外見に影響を及 ぼ す。詳細は gmtdefaults の man ページを参照すること。

−V

冗長モードを選択する。標準エラー出力に経過報告を送る [デフ ォルトでは「黙って」実行する]。

−W

type はもし存在すれば目盛付き等値線については a 、普通の 等 値 線については c である[デフォルト]。 pen にはその特定の線の属 性を設定する。目盛付き等値線のデフォルトは width = 3, color = black, texture = solid である。普通の等値線のデフォルトは width = 1 である。 pen はコンマで区切られた width, color, texture の リ ストで、それぞれはオプションである。 width は長さの単位(ポイ ント、センチメートル、 イ ン チ) 又 は faint, thin[ner|nest], thick[er|est], fat[ter|test], obese のいずれかで指定することが できる。 color はグレースケール (0−255) 又はカラー (r/g/b, それ ぞ れ範囲は 0−255; h-s-v, 範囲は 0−360, 0−1, 0−1; c/m/y/k, それ ぞれ範囲は 0−100%; 有効なカラー名) により指定できる 。 texture はダッシュ ‘-’ と点 ‘.’ の組合せである。 + フラグがあるときは等 値線の色は cpt ファイルから与えられる( −C 参照)。

−X −Y

図の原点を現在の原点に対して (x-shift,y-shift) だけ動 か す 。随意で長さの単位 (c, i, m, p) を付けることができる。値の前 に a を付けると描画後に原点を元の位置に戻し、値の前に r [デフォ ルト]を付けると現在の原点を新しい位置として再設定する。 −O を使 うとデフォルトの (x-shift,y-shift) が (0,0) になり、そうでな い ときは (r1i, r1i) 又は (r2.5c, r2.5c) になる。代わりに、 c は現 在のページサイズに基づくページの中心を使って図を中央の座 標 (x 又は y) に動かす。

−Z

等値線作成を始める前にデータから shift を引いてその結果に factor を掛ける[1/0]。 ( −A, −C, −L はこのスケーリングの後の 値 を 適用する。) このグリッドファイルに 360 度周期の z 値が含まれ ているときは p を付けること (例えば位相データ、角度分布)。そ う すれば 0 等値線を決定するときに考慮される。

−c

図のコピー枚数を指定する[デフォルトは 1 ]。

−f

入 出 力列の特殊フォーマットを使う(時間又は地理データ)。 i(nput) 又は o(utput) を指定すること [デフォルトでは入力と出 力 の 両方]。コンマ区切りで 1 列以上(又は列の範囲)を与えること。各 列又は列の範囲の項目に対して T (絶対暦時刻)、 t (選 択 さ れ た TIME_EPOCH に対する相対時間)、 x (経度)、 y (緯度)、 f (小数) のいずれかを付けること。 −f[i|o]g−f[i|o]0x,1y (地理座標) の 省略形である。

−bo

バイナリ出力を選択する。単精度に対しては s を付けること [ デフォルトは d (倍精度)]。大文字の S (又は D) にするとバイト ス ワ ップを行う。随意で、バイナリ出力ファイルの必要な列数 ncol を 付けること。随意で、バイナリファイルの列数がこのプログラムで 必 要 な列数を越えるときは、バイナリファイルの列数 ncol を付けるこ と。

アスキー形式の精度

数値データのアスキー出力形式は .gmtdefaults4 ファイルのパラメータによっ て 決まる。緯度経度は OUTPUT_DEGREE_FORMAT に従ってフォーマットされ、一 方、他の値は D_FORMAT に従ってフォーマットされる。フォーマットの結果 と し て出力の精度が落ちて、その先の処理に影響する可能性があることに注意す ること。十分な精度で出力されていないと思ったら、バイナリ出力に変更する( も し使えるなら −bo とする)か、 D_FORMAT の設定を用いてより多くの桁数を 指定すること。

ファイル hawaii_grav.grd のデータを用いてメルカトル地図上に 0.5 インチ/ 度のスケールで 25 mGal 間隔の等値線を引き、 50 mGal 間隔で目盛を付け(フ ォントサイズは 10 を使用)、 1 度間隔のケバと 30 分間隔の格子線を描く。

grdcontour hawaii_grav.grd −Jm0.5i −C25 −A50+f10 −B1g30m > hawaii_grav.ps

ファイル image.grd の等値線図を、ファイル cont.d のレベルを用いて、線形 投影図上に、 0.1 cm/x-unit 及び 50 cm/y-unit で 20 (x) 及び 0.1 (y) 間 隔 のケバを用いて、等値線を多少滑らかにし、 "RMS Misfit" を図のタイトル とし、目盛付き等値線を太い赤線に、それ以外の等値線を細い青の破線で描 く 。

grdcontour image.grd −Jx0.1c/50.0c −Ccont.d −S4 −B20/0.1:."RMS Misfit": −Wathick,red −Wcthinnest,blue,- | lp

バグ

ラベルの位置を指定するのに等値線の位置の平均値しか用いていないため、 局 所 的な高部や低部のラベルは最も内側の等値線の外側に配置されることがある 。

関連項目

GMT(l), gmtdefaults(l), psbasemap(l), grdimage(l), grdview(l), pscontour(l)