superposition of two grid files 3D  - 2種類のデータを重ねて表示（3次元）

---

  使うコマンド　　makecpt, grdgradient, grdview



bash スクリプト例
#
#   superposition of two grid files  2D
#
reliefgrdfile=izu.altbathy.grd    #grid data for relief
reliefintfile=izu.altbathy.int    #intensity (light) grid
drapegrdfile=izu.mag.grd    # grid data for grdimage (color)
region=140/145/27/30      # map region
proj=M10                    # projection and scale
vscale=0.0003        # vertical scale
boundary=f1a1g1/f1a1g1/f1000a4000wSnEZ    # boundary info
cptfile=mag_grad.cpt    # color table file
climit=-400/400/100     # color table limit
light=90            # light angle
view=160/35           # view angle azimth/inclination
psfile=izu.mag+topo.3D.ps        # output postscript file
#
# making color table
gmt makecpt -Crainbow -T$climit -V -Z > $cptfile
gmt grdgradient $reliefgrdfile -A$light -G$reliefintfile -Ne0.5 -V
#
# plot
gmt grdview $reliefgrdfile -G$drapegrdfile -J$proj -Jz$vscale -C$cptfile -I$reliefintfile -p$view -B$boundary -Qi300  -V  > $psfile
#



Tips

• 2種類のグリッドが同じスペック（範囲やグリッドサイズ）である必要はないが、より精度のよいほうにあわせて粗いデータを補間しておいたほうが見た目はよいものができる。
  
• 一般的には垂直方向の強調はあまり強くしないほうがよい。