地図上で対象領域を示す - rffbl気象学的データ解析メモ

地図上に対象領域を図示する

GMTで地図を描いて研究対象領域を点線で示す。
残念ながらGMTには四角の枠を描くコマンドはない。 psxyで線を4本描くことによって四角の枠を表現する。

Bashスクリプト例

#!/bin/bash
RANGE='135/140/32/37'
PROJ='M15'
PEN='6/red'

gmtdefaults -D > .gmtdefaults
gmtset PAPER_MEDIA A4+

pscoast -R$RANGE -J$PROJ -A5 -Df -W2 -Ba1f0.5WeSn -P -K > waku.eps
psxy -R$RANGE -J$PROJ -W$PEN -O -K <<EOF >> waku.eps
136 33
138 33
138 35
136 35
136 33
EOF
convert -quality 100 waku.eps waku.png

この例では136E 33N, 138E 33N, 138E 35N, 136E 35N
で領域を囲っている。