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

---

  使うコマンド　　makecpt, grdimage, grdcontour,basemap



bash スクリプト例
#
#   superposition of two grid files  2D
#
cntgrdfile=izu.altbathy.grd    #grid data for contouring
imggrdfile=izu.mag.grd    # grid data for grdimage (color)
region=140/145/27/30      # map region
proj=M10                    # projection and scale
ticks=f1a1g1            # tick info
frame=WSne              # frame info
scale_loc=5/-1/8/0.3h   # location of color scale
cptfile=mag_grad.cpt    # color table file
climit=-400/400/100     # color table limit
cint=500                # contour interval
# plot
gmt begin izu.mag+topo/2D
    gmt makecpt -Crainbow -T$climit -Z
    gmt grdimage $imggrdfile -R$region -J$proj
    gmt grdcontour $cntgrdfile -C$cint -A$aint -L$limit -Wthinnest,black
    gmt colorbar -Dx1/-1+w8/0.3+h
    gmt basemap -B$ticks -B$frame
gmt end show
#



Tips

• 任意のデータ(上の例は磁気異常）をプロットする際に、等深線で地形の概略を同時に示す。この例では磁気異常（カラー）の contourは描いていないが、grdcontourを２回使って違う色のペン（-W指定）でcontourも2種類描くという手も使える。その場合 は、片方のペン色を淡いグレーくらいにするといい。等深線は細いペンで描くこと。
  
• 地形のグリッドからgrdgradientで陰影ファイルを作成しておいて、grdimageで色は磁気異常（など）、影は地形にして重ねることもできるが、多くの場合これをやると地形の印象が強く出過ぎるのであまりお勧めできない。