color filled contour map 1 - 彩色した図と等深線を重ねる(グラデーション)
使うコマンド makecpt, grdimage, grdcontour, colorbar, coast basemap
bash スクリプト例
# plot bathymetry color fill map
# 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+tcolor_fill
# boundary frame info
climit=-8000/0/1000
# color table min/max/interval
cint=500
# contour interval
limit=-10000/-500
# contour min/max
aint=2000
# annotation contour interval
grdfile=JTOPO1_30.grd
# input bathymetry grid file
#
# plot
gmt begin colfilcnt
gmt makecpt -Chaxby -T$climit -Z
gmt grdimage $grdfile -R$region -J$proj
gmt grdcontour $grdfile -C$cint -A$aint -L$limit -Wthinner
gmt colorbar -DJMR -B1000a2000
gmt coast -Di -Ggray -Wthin,black
gmt basemap -B$ticks -B$frame
gmt end show
Tips
- psscaleのスケールバーを横置きにした場合(hをつける)。スケールバーのアノテーションはcptファイルの行に対応して煩雑すぎることが多いので-Bオプションで明示的につける
- makecptコマンドで-Zオプションをつけると、カラーテーブルの色は滑らかに変化する