plot survey points and lines - 観測点（線）の記入　

---

  使うコマンド　　makecpt, grdimage, psxy, pstext, pscoast,  psbasemap



bash スクリプト例
#  plot points and draw survey lines

# parameter setting
region=130/145/25/35                    # map region east/west/south/north
proj=M10                                # map projection and scale
ticks=a5f1g5                            # boundary tick info
frame=WSne+tindex                        # boundary frame info
climit=-8000/4000/1000                  # color table min/max/interval
textfont="+f8p,Helvetica,-=0.5p,blue+jBL"   # text size,font,color, justification
grdfile=JTOPO1_30.grd                   # input bathymetry grid file
pointlist=DSDP.list                     # input station location data
linelist=line.list                      # input line location data
cptfile=haxby_grad.cpt                  # color table
psfile=colfil_index.ps                  # output postscript file name
#
# making color table
gmt makecpt -Chaxby -T$climit -Z > $cptfile
#
# plot
gmt grdimage $grdfile -R$region -J$proj -C$cptfile -K -V > $psfile
gmt psscale -D12/3/9/0.5 -C$cptfile -Bf1000a2000 -K -O >> $psfile
gmt psxy $pointlist -R$region -J$proj -Sc0.15 -Gblue -K -O -V >> $psfile
gmt pstext $pointlist -R$region -J$proj -F$textfont -K -O -V >> $psfile
gmt psxy -R$region -J$proj -Sa0.30 -Gred -K -O -V << END >> $psfile
137.5134 29.49045
137.8729 28.25310
END
gmt psxy $linelist -R$region -J$proj -Wthin,green -K -O -V >> $psfile
gmt pscoast -R$region -J$proj -Di -Ggray -Wthin,black -K -V -O >> $psfile
gmt psbasemap -R$region -J$proj -B$ticks -B$frame -O  -V >> $psfile
#




Tips

• 観測点などを描く場合、数点であれば、図中星印を描いている場合のようにスクリプト直接記入でもよい。別にインプットファイルを用意する場合は、(経度、緯度、値）のテキストを用意する
• psxyでシンボルを指定しないと直線をつないで描いていく。複数の測線位置をひとつのファイルに入れる場合は、測線ごとにレコード区切り">"の行を入れる
• 上記の例では、測点番号は、BL(bottom left)で文字の左下が座標位置になるよう指定。フォントや位置指定を入力ファイルに書いておく方法もあって、点ごとに指定を細かく変えたい場合はその ほうがよい（でもちょっとの重複ならお絵かきソフトで後で動かすほうが簡単）