color image map - 水深に応じて彩色した図をつくる　

---

  使うコマンド　 makecpt, grdimage, psscale, pscoast,  psbasemap



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
grdfile=JTOPO1_30.grd                   # input bathymetry grid file
cptfile=haxby_grad.cpt                  # color table
psfile=colfil.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 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

• カラーテーブルのセンスが問われる。とりあえずGMTに装備されている標準カラーパレットからmakecptコマンドでつくって、不満があれば.cptファイルをエディタで開けてRGB値や境界値を微修正する。