PSMASK

名前
書式
説明
オプション

関連事項

名前

psmask − 地図上のデータの無い範囲をクリップあるいはマスクする

書式

psmask [xyzfile] −Ixinc[unit][=|+][/yinc[unit][=|+]] −Jparameters −Rwest/east/south/north[r] [ −B[p|s]parameters ] [ −Ddumpfile ] [ −Eazimuth/elevation ] [ −F ] [ −Gfill ] [ −H[i][nrec] ] [ −K ] [ −M[flag] ] [ −N ] [ −O ] [ −P ] [ −Ssearch_radius[m|c|k|K] ] [ −T ] [ −U[/dx/dy/][label] ] [ −V ] [ −X[a|c|r][x-shift[u]] ] [ −Y[a|c|r][y-shift[u]] ] [ −ccopies ] [ −:[i|o] ] [ −bi[s|S|d|D][ncol] ]

psmask −C [ −K ] [ −O ]

説明

psmask は (x,y,z) ファイル[または標準入力]を読み取り、この情報を用いて 有効なグリッドセルを探す。 1 つ以上のデータを含むグリッドセルだけが有効 で あると見なされる。随意で、影響の及ぶ半径を指定することができる。その 場合、データ点から半径 radius 以内にある全てのグリッドセルが有効であ る と 見なされる。さらに、オプションによりこれとは逆のセンスにすることがで きる。有効または無効な点を探しながら、 psmask は( −T を用いて) これらの ノ ードをマスクするタイルを描くか、興味の対象でない領域をクリップするよ うな多角形の輪郭を作成する。クリッピングを開始したときは 、 も う 一 度 psmask−C オプション付きで呼び出すまでその効果は継続する。

xyzfile

(x,y,z) の値を持つファイル (例えば surface を実行するために使わ れたもの)。ファイルが与えられていないときは、標準入力を読み取る 。バイナリファイルの場合は −b を参照すること。

−I

x_inc [及び随意で y_inc ] は格子間隔である。分を示すには m を、秒を示すには c を付けること。もし代わりに e, k, i, n のいず れ かの単位が付け加えられれば、それぞれ刻み幅はメートル、キロメ ートル、マイル、海里のいずれかと見なされ、その領域の中央の緯 度 における経度の等価な度単位に変換される(変換方法は ELLIPSOID に 依存する)。もし /y_inc が与えられて 0 に設定されて い れ ば 、 /x_inc と同じ値に再設定され、そうでない場合は緯度の度に変換され る。もし = が付加されていたら、対応する最大の x (東) 又は y ( 北) が与えられた刻み幅にちょうど一致するように微調整される [デ フォルトでは刻み幅は与えられた領域に合うように微調整される]。最 後に、刻み幅を与える代わりに、与えられた刻み幅に + を付加するこ とにより、必要な格子点の数を指定することができる。刻み幅は格 子 点 の数と領域から再計算される。再計算される刻み幅は格子保存グリ ッドか画素保存グリッドのどちらを選ぶかに依存する。詳細は付録 B を参照すること。

−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]] (線形、対数、 指数スケーリング)

−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] の形式でなければならない 。 区切り文字を利用する場合にはその型と位置が指定どおりでなけれ ばならない(ただし入出力と描画フォーマットは変更可能である)。

オプション

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

−B

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

−C

存在するクリップパスを終了する。入力ファイルや投影情報は不 要である。必ず −O を設定しなければならない。

−D

クリップのために作成された多角形をディスクにダンプする。 −T が 設定されているときは無視される。ダンプファイルの接頭辞がない ときは mask が用いられる (ファイルは mask_*.d という名前になる) 。

−E

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

−F

画素保存にする[デフォルトは格子保存]。

−G

クリップ用の多角形[またはタイル]を選択した色で塗りつぶす [ デフォルトでは塗りつぶさない]。グレースケール (0−255) 、 カ ラ ー(r/g/b, それぞれ範囲は 0−255; h-s-v, 範囲は 0−360, 0−1, 0−1; c/m/y/k, それぞれ範囲は 0−100%; 有効なカラー名)を指定すること。 代 わ りに、 −Gpdpi/pattern 、ただし pattern には用意された模様 (1-90) 又は Sun の 1, 8, 24 ビットのラスタファイルの名前を与 え る こともできる。 dpi には画像の解像度を設定すること。 1 ビット のラスタについては −GP を使うと反転画像になる。前景や背景の色を 指 定 す る には :Fcolor[B[color]] を付けること (透明にするには color = - とする)。個々の模様に関する情報については GMT テク ニ カルリファレンス & クックブックの付録 E を参照すること。

−H

入 力 フ ァ イルにヘッダ行がある時に使う。ヘッダ行の数は .gmtdefaults4 ファイルを編集することにより変更できる。このオ プ シ ョンが使われた場合、 GMT のデフォルトではヘッダ行は 1 行であ る。入力データだけにヘッダ行が必要な場合は −Hi を使うこと [デフ ォルトでは入力データにヘッダ行があればそれも出力される]。バイナ リデータには使わないこと。

−K

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

−M

複数セグメントファイル。セグメントは特別な行で区切られる。 アスキーファイルでは最初の文字は flag [デフォルトは ’>’ ]でなけ れ ばならない。バイナリファイルでは全てのフィールドが NaN で −b に出力する列数を明示しなければならない。デフォルトでは −M の 設 定は入力と出力の両方に適用される。別々に設定するには −Mi−Mo を使うこと。

−N

テストの意味を逆にする。すなわち、データが存在する範囲を ク リップする。

−O

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

−P

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

−S

影響半径を設定する。データ点から半径 radius 以内の格子点が 有効であると見なされる [デフォルトでは 0 であり、これはデータを 含 む格子セルのみが有効であることを意味する]。 m を付けると分、 c を付けると秒の意味になる。 k を付けると km の意味にな る( −R−I の単位が度であり、距離計算に速い平面地球を使用することを 示す)。より正確にするには、 K を使うと距離が地球面に沿って計 算 される。しかし、現在の ELLIPSOID が球体に設定されているときは球 面上での大円計算が用いられる。

−T

多角形のクリップパスの代わりにタイルを使用する。タイルの 色 又は模様を設定するには −G を使うこと。

−U

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

−V

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

−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) に動かす。

−:

入出力の(経度、緯度)と(緯度、経度)を入れ替える [デフォルト は(経度、緯度)]。 i を付けると入力のみが、 o を付けると出力のみ が入れ替わる [デフォルトでは両方入れ替わる]。

−bi

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

−c

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

多 角形のクリップパスを用いて等値線地図の支配データがない領域をマスクす る上書き用 PostScript ファイルを作成する。

psmask africa_grav.xyg −R20/40/20/40 −I5m −JM10i −O −K > mask.ps

同じ例で今度はタイルを使用する。

psmask africa_grav.xyg −R20/40/20/40 −I5m −JM10i −T −O −K −Gwhite > mask.ps

関連事項

GMT(l), grdmask(l), surface(l), psbasemap(l), psclip(l)