psvelo メモ

ベクトル場を矢印で描いたり,楕円を描いたりできる.詳細は

man psvelo

で確認できる.

ベクトル場,矢印

gmt psvelo -J -R -B -Se[velscale]/[confidence]/[fontsize] -A[params] -G[fill] inputfile.txt > output.ps

# input file
input="gradient.txt"

# params and options
proj="X12/12" # -J
region="-3/3/-3/3" # -R
# -B
axes="-Ba1f1 -Bx+lX-axis -By+lY-axis -BWSne+glightblue"
# -Se<velscale>/<confidence>/<fontsize>
sc=0.4
conf=0.0

# Output 1
psname="psvelo1.ps"
optA="0.03/0.25/0.10" # LineWidth/HeadLength/HeadSize
# GMT
gmt psvelo -J$proj -R$region $axes -A$optA -Gblack -Se$sc/$conf/0 -P -K -V $input > $psname
gmt psvelo -J -R -A$optA -Gblack -Se$sc/$conf/12 -P -K -O -V <<EOF >> $psname
# lon   lat    u1    u2  sig1  sig2  cor  (option) name 
-2.00  2.80  1.00  0.00  1.00  1.00  1.0     1.0 [unit]
EOF

## Output 2
psname="psvelo2.ps"
optA="+a25+e+gred" # vector attributes
# GMT
gmt psvelo -J$proj -R$region $axes -A$optA -Gblack -Se$sc/$conf/0 -P -K -V $input > $psname
gmt psvelo -J -R -A$optA -Gblack -Se$sc/$conf/12 -P -K -O -V <<EOF >> $psname
# lon   lat    u1    u2  sig1  sig2  cor  (option) name 
-2.00  2.80  1.00  0.00  1.00  1.00  1.0     1.0 [unit]
EOF

psvelo1.png psvelo2.png

-J, -R, -B

いつもの共通オプション.省略.

-Se[velscale]/[confidence]/[fontsize]

  • velscale: 矢印全体の大きさを調節する.
  • confidence: 楕円を描く際に使う信頼限界.ここでは関係ない.
  • fontsize: captionのフォントサイズ.スケールを表示する時に使う.表示させない場合は0でOK.

-A[params]

必須ではないがほとんどの場合で必要.矢印の形状について指定する.
指定方法は一つではなく

-A0.03/0.25/0.10 # LineWidth/HeadLength/HeadSize

という方法と

-A+a25+e # vector attributes

という方法がある.vector attributesの指定子は

  • +a: Arrow Headの傘(?)の開く角度(デフォルト30deg)
  • +e: 線の末端にArrow Headを置く (+s, +mは始点,中間点)
  • +r: Arrow Headの半分右側だけ表示する(+lは左側だけ)
  • +g[fill]: Arrow Headの内部の着色 (fill attributes)
  • +p[pen]: 線のスペック(pen attributes) など.

-G[fill]

fill attributesで矢印の着色.-Aからvector attributesを指定した場合はそれ優先. Arrow Headも線もどちらも同じ色になる.

input file

例えば filegradient.txt のようなファイル. それぞれの点で,各列に以下の情報が格納されているファイルを入力する.

  • 1,2列: lon(x), lat(y)
  • 3,4列: east(x)方向, north(y)方向成分
  • 5,6列: east(x)方向, north(y)方向の不確実性
  • 7列: 2方向成分の相関
  • 8列: (空白でも可) 名前 これは-Seを指定した場合で,-Srなどの場合は入れるべき変数が異なる.
    また上にあるような図を描く場合には,5〜7列目の値を変えても何も変わらない.

Front page   Edit Diff Attach Copy Rename Reload   New List of pages Search Recent changes   Help   RSS of recent changes
Last-modified: 2018-09-25 (Tue) 12:50:49 (28d)
Powered by PHP 5.6.21. HTML convert time: 0.049 sec.