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

---

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



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
aint=2000               # annotated contour interval
limit=-100000/-500      # min/max contours
psfile=izu.mag+topo.2D.ps        # output postscript file
#
# making color table
gmt makecpt -Crainbow -T$climit -V -Z > $cptfile
#
# plot
gmt grdimage $imggrdfile -R$region -J$proj -C$cptfile -K -V > $psfile
gmt grdcontour $cntgrdfile -R$region -J$proj -C$cint -A$aint -L$limit -K -O -Wthinnest,black >> $psfile
gmt psscale -D$scale_loc -C$cptfile -K -O  >> $psfile
gmt pscoast -R$region -J$proj -Di -Ggray -Wthin,black -K -V -O >> $psfile
gmt psbasemap -R$region -J$proj -B$ticks -B$frame -O  -V >> $psfile
#



Tips

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