3 - D view 鯨観図（鳥瞰図）をつくる　

---

  使うコマンド　　makecpt, grdgradient, grdview, 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

• このスクリプトでは、-Rオプションを指定していないので、入力グリッドをすべて書かせている。もちろん領域指定してもよい。
  
• 水平方向の縮尺図法とは独立に-Jz$vscaleで縦方向の縮尺を指定する。論文の図の場合は、縦横比（vertical exaggeration )をキャプションに書くことが多いので、比がきりのよい数になるように調整するほうがよいかも。
  
• 上の例はそうなっていないが、一般的には視点（$viewパラメタ）の斜めうしろから陰（$lightパラメタ）をつけると見栄えがよい場合が多い。
  
• -Qi300でできあがり300dpiでイメージを作成している。-Qsのsurfaceプロットのほうがきれいな場合もあるがケースバイケースか。