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    # boundary info
boundaryz=f1000a4000    # 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
#
# plot
gmt begin izu.mag+topo.3D
    gmt grdgradient $reliefgrdfile -A$light -G$reliefintfile -Ne0.5
    gmt makecpt -Crainbow -T$climit -Z
    gmt grdview $reliefgrdfile -G$drapegrdfile -J$proj -Jz$vscale -I$intfile -p$view -B$boundary -Bz$boundaryz -BwSnE -Qi300
gmt end show
#



Tips

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