3 - D view 鯨観図(鳥瞰図)をつくる 

  使うコマンド  makecpt, grdgradient, grdimage, psbasemap



bash スクリプト例
#
#   make 3D view  
#
grdfile=izu.altbathy.grd    #input grid data
intfile=izu.altbathy.int    #intensity (light) grid
proj=M10                    # projection and scale
boundary=f1a1g1/f1a1g1/f1000a4000wSnEZ    # boundary info
cptfile=topo_grad.cpt    # color table file
light=90            # light angle
vscale=0.0005        # vertical scale
view=160/25           # view angle azimth/inclination
zlevel=-10000+ggray # zlevel , +g draw facade
psfile=izu.3d.ps        # output postscript file
#
gmt grdgradient $grdfile -A$light -G$intfile -Ne0.5 -V
gmt makecpt -Ctopo -T-10000/10000/1000 -Z  > $cptfile
#
gmt grdview $grdfile -J$proj -Jz$vscale -C$cptfile -I$intfile -p$view -B$boundary -Qi300 -N$zlevel -V  > $psfile
#


Tips